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:1706.05472 · Physical review. E

Dynamic scaling in the two-dimensional Ising spin glass with normal-distributed couplings

Na Xu, Kai-Hsin Wu, Shanon J. Rubin, Ying-Jer Kao, Anders W. Sandvik

Published Nov 2, 2017• 11 citations
arXivHTMLPDF

Abstract

We carry out simulated annealing and employ a generalized Kibble-Zurek scaling hypothesis to study the two-dimensional Ising spin glass with normal-distributed couplings. The system has an equilibrium glass transition at temperature T=0. From a scaling analysis when T→0 at different annealing velocities v, we find power-law scaling in the system size for the velocity required in order to relax toward the ground state, v∼L^{-(z+1/ν)}, the Kibble-Zurek ansatz where z is the dynamic critical exponent and ν the previously known correlation-length exponent, ν≈3.6. We find z≈13.6 for both the Edwards-Anderson spin-glass order parameter and the excess energy. This is different from a previous study of the system with bimodal couplings [Rubin et al., Phys. Rev. E 95, 052133 (2017)2470-004510.1103/PhysRevE.95.052133] where the dynamics is faster (z is smaller) and the above two quantities relax with different dynamic exponents (with that of the energy being larger). We argue that the different behaviors arise as a consequence of the different low-energy landscapes: for normal-distributed couplings the ground state is unique (up to a spin reflection), while the system with bimodal couplings is massively degenerate. Our results reinforce the conclusion of anomalous entropy-driven relaxation behavior in the bimodal Ising glass. In the case of a continuous coupling distribution, our results presented here also indicate that, although Kibble-Zurek scaling holds, the perturbative behavior normally applying in the slow limit breaks down, likely due to quasidegenerate states, and the scaling function takes a different form.

DOI
10.1103/physreve.96.052102

Related papers

Dual time scales in simulated annealing of a two-dimensional Ising spin glass

Shanon J. Rubin, Na Xu, Anders W. Sandvik

We apply a generalized Kibble-Zurek out-of-equilibrium scaling ansatz to simulated annealing when approaching the spin-glass transition at temperature T=0 of the two-dimensional Ising model with random J=±1 couplings. Analyzing the spin-glass order parameter and the excess energy as functions of the system size and the annealing velocity in Monte Carlo simulations with Metropolis dynamics, we find scaling where the energy relaxes slower than the spin-glass order parameter, i.e., there are two different dynamic exponents. The values of the exponents relating the relaxation time scales to the system length, τ∼L^{z}, are z=8.28±0.03 for the relaxation of the order parameter and z=10.31±0.04 for the energy relaxation. We argue that the behavior with dual time scales arises as a consequence of the entropy-driven ordering mechanism within droplet theory. We point out that the dynamic exponents found here for T→0 simulated annealing are different from the temperature-dependent equilibrium dynamic exponent z_{eq}(T), for which previous studies have found a divergent behavior: z_{eq}(T→0)→∞. Thus, our study shows that, within Metropolis dynamics, it is easier to relax the system to one of its degenerate ground states than to migrate at low temperatures between regions of the configuration space surrounding different ground states. In a more general context of optimization, our study provides an example of robust dense-region solutions for which the excess energy (the conventional cost function) may not be the best measure of success.

arXiv:1609.09024•Physical review. E•May 22, 2017•13 citations

Annealing Dynamics of the Two-Dimensional Fully Frustrated Ising Model

Kai-Hsin Wu, Anders W. Sandvik, Ying-Jer Kao

We study the dynamical scaling properties of the two-dimensional fully frustrated Ising model by simulated annealing. The model has a highly degenerate, critical ground state with well known equilibrium properties. To study its out-of-equilibrium behavior, we develop an annealing protocol that is convenient for extending the Kibble-Zurek mechanism to the case of an exponentially divergent correlation length in the limit of the temperature $T \to T_c = 0$. We find excellent agreement with the scaling forms derived, including a dependence on the initial temperature that is not present in the conventional Kibble-Zurek mechanism. The dynamic exponent $z$ is consistent among all quantities studied; fits of its dependence on the initial temperature give $z = 1.99(2)$ (susceptibility), $z = 2.072(8)$ (excess energy), and $z = 2.09(3)$ (Edwards-Anderson order parameter). These values are close to $z=2$, which is consistent with previous studies of relaxation after instantaneous quenches and of time dependent correlations in equilibrium, including the estimate $z = 2.004(9)$ from equilibrium autocorrelation functions.

arXiv:2610.25109•Oct 18, 2026•0 citations

Nontrivial critical behavior of the free energy in the two-dimensional Ising spin glass with bimodal interactions

Helmut G. Katzgraber, L. W. Lee, Ian A. Campbell

A detailed analysis of Monte Carlo data on the two-dimensional Ising spin glass with bimodal interactions shows that the free energy of the model has a nontrivial scaling. In particular, we show by studying the correlation length that much larger system sizes and lower temperatures are required to see the true critical behavior of the model in the thermodynamic limit. Our results agree with data by Lukic et al. in that the degenerate ground state is separated from the first excited state by an energy gap of 2J.

arXiv:cond-mat/0510668•Oct 25, 2005•1 citation

Bimodal and Gaussian Ising spin glasses in dimension two

P. H. Lundow, Ian A. Campbell

An analysis is given of numerical simulation data to size L=128 on the archetype square lattice Ising spin glasses (ISGs) with bimodal (±J) and Gaussian interaction distributions. It is well established that the ordering temperature of both models is zero. The Gaussian model has a nondegenerate ground state and thus a critical exponent η≡0, and a continuous distribution of energy levels. For the bimodal model, above a size-dependent crossover temperature T(*)(L) there is a regime of effectively continuous energy levels; below T(*)(L) there is a distinct regime dominated by the highly degenerate ground state plus an energy gap to the excited states. T(*)(L) tends to zero at very large L, leaving only the effectively continuous regime in the thermodynamic limit. The simulation data on both models are analyzed with the conventional scaling variable t=T and with a scaling variable τ(b)=T(2)/(1+T(2)) suitable for zero-temperature transition ISGs, together with appropriate scaling expressions. The data for the temperature dependence of the reduced susceptibility χ(τ(b),L) and second moment correlation length ξ(τ(b),L) in the thermodynamic limit regime are extrapolated to the τ(b)=0 critical limit. The Gaussian critical exponent estimates from the simulations, η=0 and ν=3.55(5), are in full agreement with the well-established values in the literature. The bimodal critical exponents, estimated from the thermodynamic limit regime analyses using the same extrapolation protocols as for the Gaussian model, are η=0.20(2) and ν=4.8(3), distinctly different from the Gaussian critical exponents.

arXiv:1506.07141•Physical review. E•Feb 11, 2016•11 citations

Universal dynamic scaling in three-dimensional Ising spin glasses

Cheng-Wei Liu, Anatoli Polkovnikov, Anders W. Sandvik +1 more

We use a nonequilibrium Monte Carlo simulation method and dynamical scaling to study the phase transition in three-dimensional Ising spin glasses. The transition point is repeatedly approached at finite velocity v (temperature change versus time) in Monte Carlo simulations starting at a high temperature. This approach has the advantage that the equilibrium limit does not have to be strictly reached for a scaling analysis to yield critical exponents. For the dynamic exponent we obtain z=5.85(9) for bimodal couplings distribution and z=6.00(10) for the Gaussian case. Assuming universal dynamic scaling, we combine the two results and obtain z=5.93±0.07 for generic 3D Ising spin glasses.

arXiv:1411.6745•Physical Review E•Aug 19, 2015•25 citations

Effective critical behavior of the two-dimensional Ising spin glass with bimodal interactions

Helmut G. Katzgraber, L. W. Lee, Ian A. Campbell

Monte Carlo data of the two-dimensional Ising spin glass with bimodal interactions are presented with the aim of understanding the low-temperature physics of the model. An analysis of the specific heat, spin-glass susceptibility, finite-size correlation length, and the Binder ratio is performed to try to verify a recent proposal in which for large system sizes and finite but low temperatures the effective critical exponents are identical to the critical exponents of the two-dimensional Ising spin glass with Gaussian interactions. Our results show that with present system sizes the recently proposed scenario, in which the two-dimensional Ising spin glass with bimodally distributed interactions is in the same universality class as the model with Gaussian-distributed disorder at low but finite temperatures, cannot be reliably proven.

arXiv:cond-mat/0609377•Physical Review B•Jan 10, 2007•32 citations
We carry out simulated annealing and employ a generalized Kibble-Zurek scaling hypothesis to study the two-dimensional Ising spin glass with normal-distributed couplings. The system has an equilibrium glass transition at temperature T=0. From a scaling analysis when T→0 at different annealing velocities v, we find power-law scaling in the system size for the velocity required in order to relax toward the ground state, v∼L^{-(z+1/ν)}, the Kibble-Zurek ansatz where z is the dynamic critical exponent and ν the previously known correlation-length exponent, ν≈3.6. We find z≈13.6 for both the Edwards-Anderson spin-glass order parameter and the excess energy. This is different from a previous study of the system with bimodal couplings [Rubin et al., Phys. Rev. E 95, 052133 (2017)2470-004510.1103/PhysRevE.95.052133] where the dynamics is faster (z is smaller) and the above two quantities relax with different dynamic exponents (with that of the energy being larger). We argue that the different behaviors arise as a consequence of the different low-energy landscapes: for normal-distributed couplings the ground state is unique (up to a spin reflection), while the system with bimodal couplings is massively degenerate. Our results reinforce the conclusion of anomalous entropy-driven relaxation behavior in the bimodal Ising glass. In the case of a continuous coupling distribution, our results presented here also indicate that, although Kibble-Zurek scaling holds, the perturbative behavior normally applying in the slow limit breaks down, likely due to quasidegenerate states, and the scaling function takes a different form.
We apply a generalized Kibble-Zurek out-of-equilibrium scaling ansatz to simulated annealing when approaching the spin-glass transition at temperature T=0 of the two-dimensional Ising model with random J=±1 couplings. Analyzing the spin-glass order parameter and the excess energy as functions of the system size and the annealing velocity in Monte Carlo simulations with Metropolis dynamics, we find scaling where the energy relaxes slower than the spin-glass order parameter, i.e., there are two different dynamic exponents. The values of the exponents relating the relaxation time scales to the system length, τ∼L^{z}, are z=8.28±0.03 for the relaxation of the order parameter and z=10.31±0.04 for the energy relaxation. We argue that the behavior with dual time scales arises as a consequence of the entropy-driven ordering mechanism within droplet theory. We point out that the dynamic exponents found here for T→0 simulated annealing are different from the temperature-dependent equilibrium dynamic exponent z_{eq}(T), for which previous studies have found a divergent behavior: z_{eq}(T→0)→∞. Thus, our study shows that, within Metropolis dynamics, it is easier to relax the system to one of its degenerate ground states than to migrate at low temperatures between regions of the configuration space surrounding different ground states. In a more general context of optimization, our study provides an example of robust dense-region solutions for which the excess energy (the conventional cost function) may not be the best measure of success.
We study the dynamical scaling properties of the two-dimensional fully frustrated Ising model by simulated annealing. The model has a highly degenerate, critical ground state with well known equilibrium properties. To study its out-of-equilibrium behavior, we develop an annealing protocol that is convenient for extending the Kibble-Zurek mechanism to the case of an exponentially divergent correlation length in the limit of the temperature T→Tc=0T \to T_c = 0T→Tc​=0. We find excellent agreement with the scaling forms derived, including a dependence on the initial temperature that is not present in the conventional Kibble-Zurek mechanism. The dynamic exponent zzz is consistent among all quantities studied; fits of its dependence on the initial temperature give z=1.99(2)z = 1.99(2)z=1.99(2) (susceptibility), z=2.072(8)z = 2.072(8)z=2.072(8) (excess energy), and z=2.09(3)z = 2.09(3)z=2.09(3) (Edwards-Anderson order parameter). These values are close to z=2z=2z=2, which is consistent with previous studies of relaxation after instantaneous quenches and of time dependent correlations in equilibrium, including the estimate z=2.004(9)z = 2.004(9)z=2.004(9) from equilibrium autocorrelation functions.
A detailed analysis of Monte Carlo data on the two-dimensional Ising spin glass with bimodal interactions shows that the free energy of the model has a nontrivial scaling. In particular, we show by studying the correlation length that much larger system sizes and lower temperatures are required to see the true critical behavior of the model in the thermodynamic limit. Our results agree with data by Lukic et al. in that the degenerate ground state is separated from the first excited state by an energy gap of 2J.
An analysis is given of numerical simulation data to size L=128 on the archetype square lattice Ising spin glasses (ISGs) with bimodal (±J) and Gaussian interaction distributions. It is well established that the ordering temperature of both models is zero. The Gaussian model has a nondegenerate ground state and thus a critical exponent η≡0, and a continuous distribution of energy levels. For the bimodal model, above a size-dependent crossover temperature T(*)(L) there is a regime of effectively continuous energy levels; below T(*)(L) there is a distinct regime dominated by the highly degenerate ground state plus an energy gap to the excited states. T(*)(L) tends to zero at very large L, leaving only the effectively continuous regime in the thermodynamic limit. The simulation data on both models are analyzed with the conventional scaling variable t=T and with a scaling variable τ(b)=T(2)/(1+T(2)) suitable for zero-temperature transition ISGs, together with appropriate scaling expressions. The data for the temperature dependence of the reduced susceptibility χ(τ(b),L) and second moment correlation length ξ(τ(b),L) in the thermodynamic limit regime are extrapolated to the τ(b)=0 critical limit. The Gaussian critical exponent estimates from the simulations, η=0 and ν=3.55(5), are in full agreement with the well-established values in the literature. The bimodal critical exponents, estimated from the thermodynamic limit regime analyses using the same extrapolation protocols as for the Gaussian model, are η=0.20(2) and ν=4.8(3), distinctly different from the Gaussian critical exponents.
We use a nonequilibrium Monte Carlo simulation method and dynamical scaling to study the phase transition in three-dimensional Ising spin glasses. The transition point is repeatedly approached at finite velocity v (temperature change versus time) in Monte Carlo simulations starting at a high temperature. This approach has the advantage that the equilibrium limit does not have to be strictly reached for a scaling analysis to yield critical exponents. For the dynamic exponent we obtain z=5.85(9) for bimodal couplings distribution and z=6.00(10) for the Gaussian case. Assuming universal dynamic scaling, we combine the two results and obtain z=5.93±0.07 for generic 3D Ising spin glasses.
Monte Carlo data of the two-dimensional Ising spin glass with bimodal interactions are presented with the aim of understanding the low-temperature physics of the model. An analysis of the specific heat, spin-glass susceptibility, finite-size correlation length, and the Binder ratio is performed to try to verify a recent proposal in which for large system sizes and finite but low temperatures the effective critical exponents are identical to the critical exponents of the two-dimensional Ising spin glass with Gaussian interactions. Our results show that with present system sizes the recently proposed scenario, in which the two-dimensional Ising spin glass with bimodally distributed interactions is in the same universality class as the model with Gaussian-distributed disorder at low but finite temperatures, cannot be reliably proven.