Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2505.07038

Temperature chaos may emerge many thermodynamic states in spin glasses

Wenlong Wang

Published May 11, 2025• 0 citations
arXivHTMLPDF

Abstract

We present a large-scale simulation of the three-dimensional and mean-field spin glasses down to a very low but finite temperature. We extrapolate pertinent observables, e.g., the disorder-averaged central weight to zero temperature, finding that many thermodynamic states at a finite temperature and two ground states at zero temperature are fully compatible. While the disorder-averaged central weight monotonically decreases with decreasing temperature, this is far from true for individual samples. This motivates us to link this behaviour with the well-known temperature chaos. At an observing temperature, a sample may or may not have pure state coexistence depending on whether it is undergoing temperature chaos, which is a random process. Therefore, temperature chaos is likely responsible for the emergence of many pure states, providing a natural and intuitive explanation for the coexistence of expensive domain-wall excitations and many pure states at the disorder-averaged level.

DOI
10.48550/arxiv.2505.07038

Related papers

Evidence against a mean-field description of short-range spin glasses revealed through thermal boundary conditions

Wenlong Wang, Jonathan L. Machta, Helmut G. Katzgraber

A theoretical description of the low-temperature phase of short-range spin glasses has remained elusive for decades. In particular, it is unclear if theories that assert a single pair of pure states, or theories that are based on infinitely many pure states---such as replica symmetry breaking---best describe realistic short-range systems. To resolve this controversy, the three-dimensional Edwards-Anderson Ising spin glass in thermal boundary conditions is studied numerically using population annealing Monte Carlo. In thermal boundary conditions all eight combinations of periodic vs antiperiodic boundary conditions in the three spatial directions appear in the ensemble with their respective Boltzmann weights, thus minimizing finite-size corrections due to domain walls. From the relative weighting of the eight boundary conditions for each disorder instance a sample stiffness is defined, and its typical value is shown to grow with system size according to a stiffness exponent. An extrapolation to the large-system-size limit is in agreement with a description that supports the droplet picture and other theories that assert a single pair of pure states. The results are, however, incompatible with the mean-field replica symmetry breaking picture, thus highlighting the need to go beyond mean-field descriptions to accurately describe short-range spin-glass systems.

arXiv:1408.0438•Physical Review B•Nov 12, 2014•32 citations

Evidence of many thermodynamic states of the three-dimensional Ising spin glass

Wenlong Wang, Mats Wallin, Jack Lidmar

We present a large-scale simulation of the three-dimensional Ising spin glass with Gaussian disorder to low temperatures and large sizes using optimized population annealing Monte Carlo.Our primary focus is investigating the number of pure states regarding a controversial statistic, characterizing the fraction of centrally peaked disorder instances, of the overlap function order parameter.We observe that this statistic is subtly and sensitively influenced by the slight fluctuations of the integrated central weight of the disorder-averaged overlap function, making the asymptotic growth behavior very difficult to identify.Modified statistics effectively reducing this correlation are studied, and essentially monotonic growth trends are obtained.The effect of temperature is also studied, finding a larger growth rate at a higher temperature.Our state-of-the-art simulation and variance reduction data analysis suggest that the many pure states picture is most likely and coherent.

arXiv:2005.05805•Physical Review Research•Nov 16, 2020•6 citations

Fragility of the Equilibrium State in the Four-Dimensional ±JIsing Spin Glass

M. Sasaki, Koji Hukushima, Hajime Yoshino +1 more

Domain-wall free-energy, the difference in the free-energy between the system with the periodic boundary condition (BC) and that with the anti-periodic BC, is measured in the 4 dimensional ±J Ising spin glass. Temperature dependence of the domain-wall free-energy for each sample shows that the favorable BC changes with temperature, suggesting the existence of temperature chaos. We also investigate fragility of the equilibrium state against small changes of couplings { Jij } (bond chaos) by measuring the domain-wall free-energy. As a result, we find scaling behavior which suggests that the chaos exponent and the scaling function for temperature perturbation are the same as those for bond perturbation, meaning that temperature chaos and bond chaos belong to the same universarity class on fragility.

Progress of Theoretical Physics Supplement•Jan 1, 2005•0 citations

Large Scale Excitations in Disordered Systems

Marta Sales i Pardo

Disorder is present in many systems in nature and in many different versions. For instance, the dislocations of a crystal lattice, or the randomness of the interaction between magnetic moments. One of the most studied examples is that of spin glasses because they are simple to model but keep most of the very complex features that many disordered systems have. The frustration of the ground state configuration is responsible for the existence of a gap less spectrum of excitations and a rugged and complex free-energy landscape which bring about a very slow relaxation towards the equilibrium state. The main concern of the thesis has been to study what the properties of the typical excitation, i.e. those excitations that are large and contribute dominantly to the physics in the frozen phase. The existence of these large excitations brings about large fluctuations of the order parameter, and we have shown in these theses that this feature can be exploited to study the transition of any spin glass model. Moreover, we have shown that the information about these excitations can be extracted from the statistics of the lowest lying excitations. This is because due to the random nature of spin glasses, the physics obtained from averaging over the whole spectrum of excitations of an infinite sample is equivalent to averaging over many finite systems where only the ground state and the first excitation are considered. The novelty of this approach is that we do not need to make any assumption on what are typical excitations like because we can compute them exactly using numerical methods. Finally, we have investigated the dynamics and more specifically the link between the problem of chaos and the rejuvenation phenomena observed experimentally. Rejuvenation means that when lowering the temperature the aging process restarts again from scratch. This is potentially linked with the chaos assumption which states that equilibrium configurations at two different properties are not correlated. Chaos is a large scale phenomenon possible if entropy fluctuations are large. However, in this thesis we have shown that the response to temperature changes can be large in the absence of chaos close to a localization transition where the Boltzmann weight condenses in a few states. This has been observed in simulation of the Sinai model in which this localization is realized dynamically. In this model, since at low temperatures the system gets trapped in the very deep states, the dynamics is only local, so that only small excitations contribute to the rejuvenation signal that we have been able to observe. Thus, in agreement with the hierarchical picture, rejuvenation is possible even in the absence of chaos and reflects the start of the aging process of small length scales.

•Nov 22, 2002•0 citations

Static chaos and scaling behavior in the spin-glass phase

Félix Ritort

We discuss the problem of static chaos in spin glasses. In the case of magnetic-field perturbations, we propose a scaling theory for the spin-glass phase. Using the mean-field approach we argue that some pure states are suppressed by the magnetic field and their free-energy cost is determined by the finite-temperature fixed-point exponents. In this framework, numerical results suggest that mean-field chaos exponents are probably exact in finite dimensions. If we use the droplet approach, numerical results suggest that the zero-temperature fixed-point exponent \ensuremath{\theta} is very close to (d-3)/2. In both approaches d=3 is the lower critical dimension in agreement with recent numerical simulations.

arXiv:cond-mat/9404020•Physical review. B, Condensed matter•Sep 1, 1994•53 citations

Low-temperature behavior of two-dimensional Gaussian Ising spin glasses

J. Houdayer, Alexander K. Hartmann

We perform Monte Carlo simulations of large two-dimensional Gaussian Ising spin glasses down to very low temperatures $\ensuremath{\beta}=1∕T=50$. Equilibration is ensured by using a cluster algorithm including Monte Carlo moves consisting of flipping fundamental excitations. We study the thermodynamic behavior using the Binder cumulant, the spin-glass susceptibility, the distribution of overlaps, the overlap with the ground state, and the specific heat. We confirm that ${T}_{c}=0$. All results are compatible with an algebraic divergence of the correlation length with an exponent $\ensuremath{\nu}$. We find $\ensuremath{-}1∕\ensuremath{\nu}=\ensuremath{-}0.295(40)$, which is compatible with the value for the domain-wall and droplet exponent $\ensuremath{\theta}\ensuremath{\approx}\ensuremath{-}0.29$ found previously in ground-state studies. Hence the thermodynamic behavior of this model seems to be governed by one single exponent.

arXiv:cond-mat/0402036•Physical Review B•Jul 1, 2004•133 citations
We present a large-scale simulation of the three-dimensional and mean-field spin glasses down to a very low but finite temperature. We extrapolate pertinent observables, e.g., the disorder-averaged central weight to zero temperature, finding that many thermodynamic states at a finite temperature and two ground states at zero temperature are fully compatible. While the disorder-averaged central weight monotonically decreases with decreasing temperature, this is far from true for individual samples. This motivates us to link this behaviour with the well-known temperature chaos. At an observing temperature, a sample may or may not have pure state coexistence depending on whether it is undergoing temperature chaos, which is a random process. Therefore, temperature chaos is likely responsible for the emergence of many pure states, providing a natural and intuitive explanation for the coexistence of expensive domain-wall excitations and many pure states at the disorder-averaged level.
A theoretical description of the low-temperature phase of short-range spin glasses has remained elusive for decades. In particular, it is unclear if theories that assert a single pair of pure states, or theories that are based on infinitely many pure states---such as replica symmetry breaking---best describe realistic short-range systems. To resolve this controversy, the three-dimensional Edwards-Anderson Ising spin glass in thermal boundary conditions is studied numerically using population annealing Monte Carlo. In thermal boundary conditions all eight combinations of periodic vs antiperiodic boundary conditions in the three spatial directions appear in the ensemble with their respective Boltzmann weights, thus minimizing finite-size corrections due to domain walls. From the relative weighting of the eight boundary conditions for each disorder instance a sample stiffness is defined, and its typical value is shown to grow with system size according to a stiffness exponent. An extrapolation to the large-system-size limit is in agreement with a description that supports the droplet picture and other theories that assert a single pair of pure states. The results are, however, incompatible with the mean-field replica symmetry breaking picture, thus highlighting the need to go beyond mean-field descriptions to accurately describe short-range spin-glass systems.
We present a large-scale simulation of the three-dimensional Ising spin glass with Gaussian disorder to low temperatures and large sizes using optimized population annealing Monte Carlo.Our primary focus is investigating the number of pure states regarding a controversial statistic, characterizing the fraction of centrally peaked disorder instances, of the overlap function order parameter.We observe that this statistic is subtly and sensitively influenced by the slight fluctuations of the integrated central weight of the disorder-averaged overlap function, making the asymptotic growth behavior very difficult to identify.Modified statistics effectively reducing this correlation are studied, and essentially monotonic growth trends are obtained.The effect of temperature is also studied, finding a larger growth rate at a higher temperature.Our state-of-the-art simulation and variance reduction data analysis suggest that the many pure states picture is most likely and coherent.
Domain-wall free-energy, the difference in the free-energy between the system with the periodic boundary condition (BC) and that with the anti-periodic BC, is measured in the 4 dimensional ±J Ising spin glass. Temperature dependence of the domain-wall free-energy for each sample shows that the favorable BC changes with temperature, suggesting the existence of temperature chaos. We also investigate fragility of the equilibrium state against small changes of couplings { Jij } (bond chaos) by measuring the domain-wall free-energy. As a result, we find scaling behavior which suggests that the chaos exponent and the scaling function for temperature perturbation are the same as those for bond perturbation, meaning that temperature chaos and bond chaos belong to the same universarity class on fragility.
Disorder is present in many systems in nature and in many different versions. For instance, the dislocations of a crystal lattice, or the randomness of the interaction between magnetic moments. One of the most studied examples is that of spin glasses because they are simple to model but keep most of the very complex features that many disordered systems have. The frustration of the ground state configuration is responsible for the existence of a gap less spectrum of excitations and a rugged and complex free-energy landscape which bring about a very slow relaxation towards the equilibrium state. The main concern of the thesis has been to study what the properties of the typical excitation, i.e. those excitations that are large and contribute dominantly to the physics in the frozen phase. The existence of these large excitations brings about large fluctuations of the order parameter, and we have shown in these theses that this feature can be exploited to study the transition of any spin glass model. Moreover, we have shown that the information about these excitations can be extracted from the statistics of the lowest lying excitations. This is because due to the random nature of spin glasses, the physics obtained from averaging over the whole spectrum of excitations of an infinite sample is equivalent to averaging over many finite systems where only the ground state and the first excitation are considered. The novelty of this approach is that we do not need to make any assumption on what are typical excitations like because we can compute them exactly using numerical methods. Finally, we have investigated the dynamics and more specifically the link between the problem of chaos and the rejuvenation phenomena observed experimentally. Rejuvenation means that when lowering the temperature the aging process restarts again from scratch. This is potentially linked with the chaos assumption which states that equilibrium configurations at two different properties are not correlated. Chaos is a large scale phenomenon possible if entropy fluctuations are large. However, in this thesis we have shown that the response to temperature changes can be large in the absence of chaos close to a localization transition where the Boltzmann weight condenses in a few states. This has been observed in simulation of the Sinai model in which this localization is realized dynamically. In this model, since at low temperatures the system gets trapped in the very deep states, the dynamics is only local, so that only small excitations contribute to the rejuvenation signal that we have been able to observe. Thus, in agreement with the hierarchical picture, rejuvenation is possible even in the absence of chaos and reflects the start of the aging process of small length scales.
We discuss the problem of static chaos in spin glasses. In the case of magnetic-field perturbations, we propose a scaling theory for the spin-glass phase. Using the mean-field approach we argue that some pure states are suppressed by the magnetic field and their free-energy cost is determined by the finite-temperature fixed-point exponents. In this framework, numerical results suggest that mean-field chaos exponents are probably exact in finite dimensions. If we use the droplet approach, numerical results suggest that the zero-temperature fixed-point exponent \ensuremath{\theta} is very close to (d-3)/2. In both approaches d=3 is the lower critical dimension in agreement with recent numerical simulations.
We perform Monte Carlo simulations of large two-dimensional Gaussian Ising spin glasses down to very low temperatures $\ensuremath{\beta}=1∕T=50$. Equilibration is ensured by using a cluster algorithm including Monte Carlo moves consisting of flipping fundamental excitations. We study the thermodynamic behavior using the Binder cumulant, the spin-glass susceptibility, the distribution of overlaps, the overlap with the ground state, and the specific heat. We confirm that Tc=0{T}_{c}=0Tc​=0. All results are compatible with an algebraic divergence of the correlation length with an exponent $\ensuremath{\nu}$. We find $\ensuremath{-}1∕\ensuremath{\nu}=\ensuremath{-}0.295(40)$, which is compatible with the value for the domain-wall and droplet exponent $\ensuremath{\theta}\ensuremath{\approx}\ensuremath{-}0.29$ found previously in ground-state studies. Hence the thermodynamic behavior of this model seems to be governed by one single exponent.