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:cond-mat/0309238 · Physical Review Letters

Critical Thermodynamics of the Two-Dimensional ±J Ising Spin Glass

Jovanka Lukic, Anna Galluccio, Enzo Marinari, Olivier Martin, Giovanni Rinaldi

Published Mar 18, 2004• 43 citations
arXivHTMLPDF

Abstract

We compute the exact partition function of 2d Ising spin glasses with binary couplings. In these systems, the ground state is highly degenerate and is separated from the first excited state by a gap of size 4J. Nevertheless, we find that the low temperature specific heat density scales as exp(-2J/T), corresponding to an "effective" gap of size 2J; in addition, an associated crossover length scale grows as exp(J/T). We justify these scalings via the degeneracy of the low lying excitations and by the way low energy domain walls proliferate in this model.

DOI
10.1103/physrevlett.92.117202

Related papers

Energy gap of the bimodal two-dimensional Ising spin glass

W Atisattapong, J Poulter

An exact algorithm is used to compute the degeneracies of the excited states of the bimodal Ising spin glass in two dimensions. It is found that the specific heat at arbitrary low temperature is not a self-averaging quantity and has a distribution that is neither normal nor lognormal. Nevertheless, it is possible to estimate the most likely value and this is found to scale as L 3 T -2 exp (−4 J / kT ), for a L × L lattice. Our analysis also explains, for the first time, why a correlation length ξ∼exp (2 J / kT ) is consistent with an energy gap of 2 J . Our method allows us to obtain results for up to 10 5 disorder realizations with L ⩽96. Distributions of second and third excitations are also shown.

arXiv:0804.1857•New Journal of Physics•Sep 11, 2008•7 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

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

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

Excitations of the bimodal Ising spin glass on the brickwork lattice

W Atisattapong, J Poulter

An exact algorithm is used to investigate the distributions of the degeneracies of low-energy excited states for the bimodal Ising spin glass on the brickwork lattice. Since the distributions are extreme and do not self-average, we base our conclusions on the most likely values of the degeneracies. Our main result is that the degeneracy of the first excited state per ground state and per spin is finite in the thermodynamic limit. This is very different from the same model on a square lattice where a divergence proportional to the linear lattice size is expected. The energy gap for the brickwork lattice is obviously 2 J on finite systems and predicted to be the same in the thermodynamic limit. Our results suggest that a 2 J gap is universal for planar bimodal Ising spin glasses. The distribution of the second contribution to the internal energy has a mode close to zero and we predict that the low-temperature specific heat is dominated by the leading term proportional to T −2 exp (−2 J / kT ).

arXiv:0901.1538•New Journal of Physics•Jun 22, 2009•3 citations

Finite-size scaling analysis of exact ground states for ± J spin glass models in two dimensions

Naoki Kawashima, Heiko Rieger

With the help of exact ground states obtained by a polynomial algorithm we compute the domain wall energy Δ at zero temperature for the bond-random and the site-random Ising spin glass model in two dimensions. We find that in both models the stability of the ferromagnetic and the spin glass order ceases to exist at a unique concentration p c for the ferromagnetic bonds. In the vicinity of this critical point, the size and concentration dependence of the first and second moment of the domain wall energy are, for both models, described by a common finite-size scaling form. Moreover, below this concentration the stiffness exponent turns out to be slightly negative, θ S = − 0.056(6), indicating the absence of any intermediate spin glass phase at non-zero temperature.

arXiv:cond-mat/9612116•Europhysics Letters (EPL)•Jul 1, 1997•81 citations
We compute the exact partition function of 2d Ising spin glasses with binary couplings. In these systems, the ground state is highly degenerate and is separated from the first excited state by a gap of size 4J. Nevertheless, we find that the low temperature specific heat density scales as exp(-2J/T), corresponding to an "effective" gap of size 2J; in addition, an associated crossover length scale grows as exp(J/T). We justify these scalings via the degeneracy of the low lying excitations and by the way low energy domain walls proliferate in this model.
An exact algorithm is used to compute the degeneracies of the excited states of the bimodal Ising spin glass in two dimensions. It is found that the specific heat at arbitrary low temperature is not a self-averaging quantity and has a distribution that is neither normal nor lognormal. Nevertheless, it is possible to estimate the most likely value and this is found to scale as L 3 T -2 exp (−4 J / kT ), for a L × L lattice. Our analysis also explains, for the first time, why a correlation length ξ∼exp (2 J / kT ) is consistent with an energy gap of 2 J . Our method allows us to obtain results for up to 10 5 disorder realizations with L ⩽96. Distributions of second and third excitations are also shown.
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.
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.
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.
An exact algorithm is used to investigate the distributions of the degeneracies of low-energy excited states for the bimodal Ising spin glass on the brickwork lattice. Since the distributions are extreme and do not self-average, we base our conclusions on the most likely values of the degeneracies. Our main result is that the degeneracy of the first excited state per ground state and per spin is finite in the thermodynamic limit. This is very different from the same model on a square lattice where a divergence proportional to the linear lattice size is expected. The energy gap for the brickwork lattice is obviously 2 J on finite systems and predicted to be the same in the thermodynamic limit. Our results suggest that a 2 J gap is universal for planar bimodal Ising spin glasses. The distribution of the second contribution to the internal energy has a mode close to zero and we predict that the low-temperature specific heat is dominated by the leading term proportional to T −2 exp (−2 J / kT ).
With the help of exact ground states obtained by a polynomial algorithm we compute the domain wall energy Δ at zero temperature for the bond-random and the site-random Ising spin glass model in two dimensions. We find that in both models the stability of the ferromagnetic and the spin glass order ceases to exist at a unique concentration p c for the ferromagnetic bonds. In the vicinity of this critical point, the size and concentration dependence of the first and second moment of the domain wall energy are, for both models, described by a common finite-size scaling form. Moreover, below this concentration the stiffness exponent turns out to be slightly negative, θ S = − 0.056(6), indicating the absence of any intermediate spin glass phase at non-zero temperature.