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

Physical review. B, Condensed matter

Monte Carlo studies of Ising spin-glass systems: Investigations of the overlap length in spin glasses

J.-O. Andersson, Johan Mattsson, Peter Svedlindh

Published Jan 1, 1994• 13 citations
Publisher

Abstract

In order to determine if the effects of an overlap length are visible in an Ising Edwards-Anderson spin-glass model, we have performed Monte Carlo simulations that mimic experimental procedures. The spin system was quenched from a random initial condition to a temperature ${\mathit{T}}_{1}$, equilibrated in zero field for a time ${\mathit{t}}_{\mathit{w}1}$, whereupon the system was subjected to a perturbation in field \ensuremath{\Delta}H, temperature \ensuremath{\Delta}T, or bond \ensuremath{\Delta}J with a duration ${\mathit{t}}_{\mathit{w}2}$. After the perturbation of the original field, the temperature or bonds were restored. After waiting for a time ${\mathit{t}}_{\mathit{w}3}$, a weak magnetic field H was applied. The susceptibility \ensuremath{\chi}(t)=[${\mathit{M}}_{+}$(t)-${\mathit{M}}_{\mathrm{\ensuremath{-}}}$(t)]/2H, where ${\mathit{M}}_{\ifmmode\pm\else\textpm\fi{}}$(t) is the magnetization due to positive or negative field, and the spin autocorrelation function q(t) were studied as functions of time. For small perturbations, \ensuremath{\Delta}H${\mathit{H}}_{0}$, \ensuremath{\Delta}T${\mathit{T}}_{0}$, or \ensuremath{\Delta}J${\mathit{J}}_{0}$, respectively, the system is undisturbed. For larger perturbations, the correlations built up during ${\mathit{t}}_{\mathit{w}1}$ are gradually destroyed. Similar results are found in all three types of simulations and the qualitative features are in agreement with experimental results. These features are found for both two-dimensional and three-dimensional spin-glass systems. The results can be interpreted in the domain models for spin glasses. \ensuremath{\Delta}${\mathit{H}}_{0}$ is found to be strongly dependent on the duration ${\mathit{t}}_{\mathit{w}2}$ of the perturbation. If ${\mathit{t}}_{\mathit{w}1}$ is constant and ${\mathit{t}}_{\mathit{w}2}$ increases, \ensuremath{\Delta}${\mathit{H}}_{0}$ decreases. If ${\mathit{t}}_{\mathit{w}1}$ is set equal to ${\mathit{t}}_{\mathit{w}2}$ and ${\mathit{t}}_{\mathit{w}1}$=${\mathit{t}}_{\mathit{w}2}$ is increased, \ensuremath{\Delta}${\mathit{H}}_{0}$ decreases, in agreement with the expected behavior. \ensuremath{\Delta}${\mathit{H}}_{0}$ was found to increase with increasing temperature T, which is in agreement with experimental results, but contrary to the theoretical expectations.

DOI
10.1103/physrevb.49.1120

Related papers

Monte Carlo studies of Ising spin-glass systems: Aging behavior and crossover between equilibrium and nonequilibrium dynamics

J.-O. Andersson, Johan Mattsson, Peter Svedlindh

The crossover between quasiequilibrium and nonequilibrium spin-glass dynamics has been studied in Monte Carlo simulations of two- and three-dimensional short-range Ising spin-glass systems. The spin system was quenched in zero field to a low temperature T. After equilibrating for a time ${\mathit{t}}_{\mathit{w}}$, a weak magnetic field was applied and the time dependence of the magnetization M(t) and the spin autocorrelation function q(t) was studied. Our results show that the relaxation rate S(t)=\ensuremath{\partial}M(t)/\ensuremath{\partial} ${\mathrm{log}}_{10}$(t) exhibits a maximum at t\ensuremath{\approxeq}${\mathit{t}}_{\mathit{w}}$. If the temperature is lowered or raised, immediately prior to the field application, the apparent age (${\mathit{t}}_{\mathit{a}}$) of the system, as defined by the maximum in relaxation rate, is shifted to longer or shorter time scales, respectively. No significant differences between the aging behavior in two and three dimensions are found. The relation \ensuremath{\chi}(t)=[1-q(t)]/T holds only for short time scales (t\ensuremath{\ll}${\mathit{t}}_{\mathit{w}}$), clear deviations are found at longer time scales (t>${\mathit{t}}_{\mathit{w}}$). Since this relation should according to the fluctuation-dissipation theorem hold for a system in equilibrium, this result is interpreted as an indication of a crossover between equilibrium and nonequilibrium dynamics. Our results are qualitatively in agreement with experiments as well as with theoretical models for the nonequilibrium spin-glass dynamics.

Physical review. B, Condensed matter•Oct 1, 1992•65 citations

Replica Exchange Monte Carlo Simulations of the Ising Spin Glass: Static and Dynamic Properties

Burcu Yücesoy

Spin glasses have been the subject of intense study and considerable controversy for decades, and the low-temperature phase of short-range spin glasses is still poorly understood. Our main goal is to improve our understanding in this area and find an answer to the following question: Are there only a single pair or a countable infinity of pure states in the low temperature phase of the EA spin glass? To that aim we first start by introducing spin glasses and provide a brief history of their research, then proceed to describe our method of simulation, the parallel tempering Monte Carlo algorithm. Next, we present the results of a large-scale numerical study of the equilibrium three-dimensional Edwards-Anderson Ising spin glass with Gaussian disorder. In order to understand how the parallel tempering algorithm works, we measure various static, as well as dynamical quantities, such as the autocorrelation times and round-trip times for the parallel tempering Monte Carlo method. We examine the correlation between static and dynamic observables for _ 5000 disorder realizations and up to 1000 spins down to temperatures at 20% of the critical temperature, and our results show that autocorrelation times are directly correlated with the roughness of the free-energy landscape. In the following chapters, the three- and four-dimensional Edwards-Anderson and mean-field Sherrington-Kirkpatrick Ising spin glasses are studied again via large-scale Monte Carlo simulations at low temperatures, deep within the spinglass phase. Performing a careful statistical analysis of several thousand independent disorder realizations and using an observable that detects peaks in the overlap distribution, we show that the Sherrington-Kirkpatrick and Edwards-Anderson models have a distinctly different low-temperature behavior. We arrive to the following conclusion: The structure of the spin-glass overlap distribution for the Edwards-Anderson model suggests that its low-temperature phase has only a single pair of pure states. Finally we present results for several new observables, along with a few preliminary studies and suggestions for future research.

•Apr 23, 2019•0 citations

Numerical Study on Aging Dynamics in Ising Spin-Glass Models: Temperature-Change Protocols

達雄 小森, Hajime Yoshino, Hajime Takayama

By means of Monte Carlo simulations on the three-dimensional Ising spin-glass model, we have studied aging phenomena with various temperature($T$)-change protocols. Particularly, a $T$-shift protocol, in which a system is first quenched to and aged for a period $t_{\rm w1}$ at a temperature $T_1$, and subsequently aged at a new temperature $T_2$ is closely investigated. Most importantly, the mean size of domains, which is extracted from the replica-overlap function, was found to grow monotonically without any appreciable decrease by the $T$-change. We also found the relaxation of energy density and spin auto-correlation function during the $T$-shift process can be explained within the following picture: the dynamics finally crossovers to an isothermal aging of $T_2$ at around $t_2 \simeq t_{\rm w1}^{\rm (eff)}$ after the $T$-change, where $t_{\rm w1}^{\rm (eff)}$ is an effective waiting time related with $t_{\rm w1}$, but in the transient regime $t_2 \lt t_{\rm w1}^{\rm (eff)}$, adjustment of the population of thermally active droplets from that at $T_1$ to $T_2$ takes place very slowly. Implications of the results on aging phenomena in other $T$-change protocols such as $T$-cycling and continuous $T$-change with an intermittent stop observed by simulations as well as experiments are discussed.

arXiv:cond-mat/0001395•Jan 27, 2000•7 citations

Aging in two-dimensional Ising spin glasses

A. G. Schins, E. M. Dons, A. F. M. Arts +4 more

Aging of two-dimensional (2D) Ising spin glasses with short-range interactions ${\mathrm{Rb}}_{2}$${\mathrm{Cu}}_{1\mathrm{\ensuremath{-}}\mathit{x}}$${\mathrm{Co}}_{\mathit{x}}$${\mathrm{F}}_{4}$ (${\mathit{T}}_{\mathit{c}}$=0) is studied in the critical regime via the low-frequency ac susceptibility and via the dc magnetization for x=0.22 and 0.33. After quenching to low temperatures, slow logarithmic relaxation is observed for both the in-phase and out-of-phase ac susceptibilities. Using the droplet scaling model, the spin-glass correlation length R is concluded to increase after a temperature quench as R\ensuremath{\propto}(lnt${)}^{1/\mathrm{\ensuremath{\psi}}}$, with t in units of some microscopic time. The exponent \ensuremath{\psi}=1.0\ifmmode\pm\else\textpm\fi{}0.1 compares with recent theoretical results and numerical simulations for 2D short-range Ising spin glasses. The waiting-time dependence found for the field-cooled and zero-field-cooled dc magnetizations points to a markedly slower crossover from the quasiequilibrium to the nonequilibrium regime than in 3D. From experiments with different cooling histories, finally, the equilibrium spin-glass state is inferred to be strongly dependent on the temperature.

Physical review. B, Condensed matter•Dec 1, 1993•18 citations

Universality in three-dimensional Ising spin glasses: Nonequilibrium dynamics from Monte Carlo simulations

FEDERICO J. ROMÁ

The nonequilibrium dynamics of the three-dimensional Edwards-Anderson spin-glass model with different bond distributions is investigated by means of Monte Carlo simulation. A numerical method is used to determine the critical temperature and the scaling exponents of the correlation and the integrated response functions. The results obtained agree with those calculated in equilibrium simulations and suggest that the universality class does not depend on the exact form of the bond distribution.

arXiv:1007.4780•Physical Review B•Dec 13, 2010•8 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
In order to determine if the effects of an overlap length are visible in an Ising Edwards-Anderson spin-glass model, we have performed Monte Carlo simulations that mimic experimental procedures. The spin system was quenched from a random initial condition to a temperature T1{\mathit{T}}_{1}T1​, equilibrated in zero field for a time tw1{\mathit{t}}_{\mathit{w}1}tw1​, whereupon the system was subjected to a perturbation in field \ensuremath{\Delta}H, temperature \ensuremath{\Delta}T, or bond \ensuremath{\Delta}J with a duration tw2{\mathit{t}}_{\mathit{w}2}tw2​. After the perturbation of the original field, the temperature or bonds were restored. After waiting for a time tw3{\mathit{t}}_{\mathit{w}3}tw3​, a weak magnetic field H was applied. The susceptibility \ensuremath{\chi}(t)=[M+{\mathit{M}}_{+}M+​(t)-M\ensuremath−{\mathit{M}}_{\mathrm{\ensuremath{-}}}M\ensuremath−​(t)]/2H, where M\ifmmode±\else\textpm\fi{\mathit{M}}_{\ifmmode\pm\else\textpm\fi{}}M\ifmmode±\else\textpm\fi​(t) is the magnetization due to positive or negative field, and the spin autocorrelation function q(t) were studied as functions of time. For small perturbations, \ensuremath{\Delta}HH0{\mathit{H}}_{0}H0​, \ensuremath{\Delta}TT0{\mathit{T}}_{0}T0​, or \ensuremath{\Delta}JJ0{\mathit{J}}_{0}J0​, respectively, the system is undisturbed. For larger perturbations, the correlations built up during tw1{\mathit{t}}_{\mathit{w}1}tw1​ are gradually destroyed. Similar results are found in all three types of simulations and the qualitative features are in agreement with experimental results. These features are found for both two-dimensional and three-dimensional spin-glass systems. The results can be interpreted in the domain models for spin glasses. \ensuremath{\Delta}H0{\mathit{H}}_{0}H0​ is found to be strongly dependent on the duration tw2{\mathit{t}}_{\mathit{w}2}tw2​ of the perturbation. If tw1{\mathit{t}}_{\mathit{w}1}tw1​ is constant and tw2{\mathit{t}}_{\mathit{w}2}tw2​ increases, \ensuremath{\Delta}H0{\mathit{H}}_{0}H0​ decreases. If tw1{\mathit{t}}_{\mathit{w}1}tw1​ is set equal to tw2{\mathit{t}}_{\mathit{w}2}tw2​ and tw1{\mathit{t}}_{\mathit{w}1}tw1​=tw2{\mathit{t}}_{\mathit{w}2}tw2​ is increased, \ensuremath{\Delta}H0{\mathit{H}}_{0}H0​ decreases, in agreement with the expected behavior. \ensuremath{\Delta}H0{\mathit{H}}_{0}H0​ was found to increase with increasing temperature T, which is in agreement with experimental results, but contrary to the theoretical expectations.
The crossover between quasiequilibrium and nonequilibrium spin-glass dynamics has been studied in Monte Carlo simulations of two- and three-dimensional short-range Ising spin-glass systems. The spin system was quenched in zero field to a low temperature T. After equilibrating for a time tw{\mathit{t}}_{\mathit{w}}tw​, a weak magnetic field was applied and the time dependence of the magnetization M(t) and the spin autocorrelation function q(t) was studied. Our results show that the relaxation rate S(t)=\ensuremath{\partial}M(t)/\ensuremath{\partial} log10{\mathrm{log}}_{10}log10​(t) exhibits a maximum at t\ensuremath{\approxeq}tw{\mathit{t}}_{\mathit{w}}tw​. If the temperature is lowered or raised, immediately prior to the field application, the apparent age (ta{\mathit{t}}_{\mathit{a}}ta​) of the system, as defined by the maximum in relaxation rate, is shifted to longer or shorter time scales, respectively. No significant differences between the aging behavior in two and three dimensions are found. The relation \ensuremath{\chi}(t)=[1-q(t)]/T holds only for short time scales (t\ensuremath{\ll}tw{\mathit{t}}_{\mathit{w}}tw​), clear deviations are found at longer time scales (t>tw{\mathit{t}}_{\mathit{w}}tw​). Since this relation should according to the fluctuation-dissipation theorem hold for a system in equilibrium, this result is interpreted as an indication of a crossover between equilibrium and nonequilibrium dynamics. Our results are qualitatively in agreement with experiments as well as with theoretical models for the nonequilibrium spin-glass dynamics.
Spin glasses have been the subject of intense study and considerable controversy for decades, and the low-temperature phase of short-range spin glasses is still poorly understood. Our main goal is to improve our understanding in this area and find an answer to the following question: Are there only a single pair or a countable infinity of pure states in the low temperature phase of the EA spin glass? To that aim we first start by introducing spin glasses and provide a brief history of their research, then proceed to describe our method of simulation, the parallel tempering Monte Carlo algorithm. Next, we present the results of a large-scale numerical study of the equilibrium three-dimensional Edwards-Anderson Ising spin glass with Gaussian disorder. In order to understand how the parallel tempering algorithm works, we measure various static, as well as dynamical quantities, such as the autocorrelation times and round-trip times for the parallel tempering Monte Carlo method. We examine the correlation between static and dynamic observables for _ 5000 disorder realizations and up to 1000 spins down to temperatures at 20% of the critical temperature, and our results show that autocorrelation times are directly correlated with the roughness of the free-energy landscape. In the following chapters, the three- and four-dimensional Edwards-Anderson and mean-field Sherrington-Kirkpatrick Ising spin glasses are studied again via large-scale Monte Carlo simulations at low temperatures, deep within the spinglass phase. Performing a careful statistical analysis of several thousand independent disorder realizations and using an observable that detects peaks in the overlap distribution, we show that the Sherrington-Kirkpatrick and Edwards-Anderson models have a distinctly different low-temperature behavior. We arrive to the following conclusion: The structure of the spin-glass overlap distribution for the Edwards-Anderson model suggests that its low-temperature phase has only a single pair of pure states. Finally we present results for several new observables, along with a few preliminary studies and suggestions for future research.
By means of Monte Carlo simulations on the three-dimensional Ising spin-glass model, we have studied aging phenomena with various temperature(TTT)-change protocols. Particularly, a TTT-shift protocol, in which a system is first quenched to and aged for a period tw1t_{\rm w1}tw1​ at a temperature T1T_1T1​, and subsequently aged at a new temperature T2T_2T2​ is closely investigated. Most importantly, the mean size of domains, which is extracted from the replica-overlap function, was found to grow monotonically without any appreciable decrease by the TTT-change. We also found the relaxation of energy density and spin auto-correlation function during the TTT-shift process can be explained within the following picture: the dynamics finally crossovers to an isothermal aging of T2T_2T2​ at around t2≃tw1(eff)t_2 \simeq t_{\rm w1}^{\rm (eff)}t2​≃tw1(eff)​ after the TTT-change, where tw1(eff)t_{\rm w1}^{\rm (eff)}tw1(eff)​ is an effective waiting time related with tw1t_{\rm w1}tw1​, but in the transient regime t2<tw1(eff)t_2 \lt t_{\rm w1}^{\rm (eff)}t2​<tw1(eff)​, adjustment of the population of thermally active droplets from that at T1T_1T1​ to T2T_2T2​ takes place very slowly. Implications of the results on aging phenomena in other TTT-change protocols such as TTT-cycling and continuous TTT-change with an intermittent stop observed by simulations as well as experiments are discussed.
Aging of two-dimensional (2D) Ising spin glasses with short-range interactions Rb2{\mathrm{Rb}}_{2}Rb2​${\mathrm{Cu}}_{1\mathrm{\ensuremath{-}}\mathit{x}}$Cox{\mathrm{Co}}_{\mathit{x}}Cox​F4{\mathrm{F}}_{4}F4​ (Tc{\mathit{T}}_{\mathit{c}}Tc​=0) is studied in the critical regime via the low-frequency ac susceptibility and via the dc magnetization for x=0.22 and 0.33. After quenching to low temperatures, slow logarithmic relaxation is observed for both the in-phase and out-of-phase ac susceptibilities. Using the droplet scaling model, the spin-glass correlation length R is concluded to increase after a temperature quench as R\ensuremath{\propto}(lnt${)}^{1/\mathrm{\ensuremath{\psi}}}$, with t in units of some microscopic time. The exponent \ensuremath{\psi}=1.0\ifmmode\pm\else\textpm\fi{}0.1 compares with recent theoretical results and numerical simulations for 2D short-range Ising spin glasses. The waiting-time dependence found for the field-cooled and zero-field-cooled dc magnetizations points to a markedly slower crossover from the quasiequilibrium to the nonequilibrium regime than in 3D. From experiments with different cooling histories, finally, the equilibrium spin-glass state is inferred to be strongly dependent on the temperature.
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.