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:0704.2748

Droplets in the two-dimensional +-J spin glass: evidence for (non-) universality

Alexander K. Hartmann

Published Apr 20, 2007• 0 citations
arXivHTMLPDF

Abstract

Using mappings to computer-science problems and by applying sophisticated algorithms, one can study numerically many problems much better compared to applying standard approaches like Monte Carlo simulations. Here, using calculations of ground states of suitable perturbed systems, droplets are obtained in two-dimensional +-J spin glasses, which are in the focus of a currently very lifely debate. Since a sophisticated matching algorithm is applied here, exact ground states of large systems up to L^2=256^2 spins can be generated. Furthermore, no equilibration or extrapolation to T=0 is necessary. Three different +-J models are studied here: a) with open boundary conditions, b) with fixed boundary conditions and c) a diluted system where a fraction p=0.125 of all bonds is zero. For large systems, the droplet energy shows for all three models a power-law behavior E_D L^θ'_D with θ'_D<0. This is different from previous studies of domain walls, where a convergence to a constant non-zero value (θ_dw=0) has been found for such models. After correcting for the non-compactness of the droplets, the results are likely to be compatible with θ_D= -0.29 for all three models. This is in accordance with the Gaussian system where θ_D=-0.287(4) (ν=3.5 via ν=-1/θ_D). Nevertheless, the disorder-averaged spin-spin correlation exponent ηis determined here via the probability to have a non-zero-energy droplet, and η~0.22$ is found for all three models, this being in contrast to the behavior of the model with Gaussian interactions, where exactly η=0.

DOI
10.48550/arxiv.0704.2748

Related papers

Droplets in the two-dimensional ±J Ising spin glass

Alexander K. Hartmann

By using problem-dependent mappings to computer-science problems and by applying sophisticated algorithms, one can numerically study some important problems much better compared to applying standard approaches, such as Monte Carlo simulations. Here, by using calculations of ground states of suitable perturbed systems, droplets are obtained in two-dimensional $\ifmmode\pm\else\textpm\fi{}J$ spin glasses, which are in the focus of a currently very lively debate. By using a sophisticated matching algorithm, exact ground states of large systems up to ${L}^{2}={256}^{2}$ spins can be generated. Furthermore, no equilibration or extrapolation to $T=0$ is necessary. Three different $\ifmmode\pm\else\textpm\fi{}J$ models are studied here: (a) with free boundary conditions, (b) with fixed boundary conditions, and (c) a diluted system wherein a fraction $p=0.125$ of all bonds is zero. For large systems, the droplet energy shows for all three models a power-law behavior ${E}_{\text{dp}}\ensuremath{\sim}{L}^{{\ensuremath{\theta}}_{\text{dp}}^{\ensuremath{'}}}$ with ${\ensuremath{\theta}}_{\text{dp}}^{\ensuremath{'}}<0$. This is different from the previous studies of domain walls, wherein a convergence to a constant nonzero value $({\ensuremath{\theta}}_{\text{DW}}=0)$ was found. After correcting for the noncompactness of the droplets, the results are likely to be compatible with ${\ensuremath{\theta}}_{\text{dp}}\ensuremath{\approx}\ensuremath{-}0.29$ for all three models. This is in accordance with the Gaussian system wherein ${\ensuremath{\theta}}_{\text{dp}}=\ensuremath{-}0.287(4)$ ($\ensuremath{\nu}\ensuremath{\approx}3.5$ via $\ensuremath{\nu}=\ensuremath{-}1/{\ensuremath{\theta}}_{\text{dp}}$). Nevertheless, the disorder-averaged spin-spin correlation exponent $\ensuremath{\eta}$ is determined here via the probability to have a nonzero-energy droplet, and $\ensuremath{\eta}\ensuremath{\approx}0.22$ is found for all three models, in contrast to the model with the Gaussian interactions, where $\ensuremath{\eta}=0$ is exact.

Physical Review B•Apr 18, 2008•20 citations

Large-scale low-energy excitations in the two-dimensional Ising spin glass

Alexander K. Hartmann, A. P. Young

We study large-scale, low-energy excitations in the Ising spin glass with Gaussian interactions in two dimensions at zero temperature, using an optimization algorithm to determine exact ground states. Periodic boundary conditions are applied. Our results for the fractal dimension of the surface, ${d}_{s},$ and stiffness exponent, ${\ensuremath{\theta}}^{\ensuremath{'}},$ for ``droplet'' excitations, are in reasonable agreement with estimates from ``domain wall'' calculations, and so support the predictions of the ``droplet theory.'' Restricting our analysis to small lattices, we do not find an effective value of ${\ensuremath{\theta}}^{\ensuremath{'}}$ close to $\ensuremath{-}0.47$ as has been recently proposed. The effects of averaging over droplets of different sizes are studied and are also found to be too small to give ${\ensuremath{\theta}}^{\ensuremath{'}}\ensuremath{\approx}\ensuremath{-}0.47$ for smaller sizes. Larger corrections to finite-size scaling would be needed in three and four dimensions in order for the numerical data used to support the ``TNT'' scenario to be compatible with the droplet theory prediction that the stiffness exponent is the same for droplets and domain walls.

arXiv:cond-mat/0205659•Physical review. B, Condensed matter•Sep 16, 2002•45 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

Overcoming System-Size Limitations in Spin Glasses

Helmut G. Katzgraber, Mathias Körner, Frauke Liers +1 more

In order to overcome the limitations of small system sizes in spin-glass simulations, we investigate the one-dimensional Ising spin chain with power-law interactions.The model has the advantage over traditional higher-dimensional Hamiltonians in that a large range of system sizes can be studied.In addition, the universality class of the model can be changed by tuning the power law exponent, thus allowing us to scan from the mean-field to long-range and short-range universality classes.We illustrate the advantages of this model by studying the nature of the spin glass state where our results hint towards a replica symmetry breaking scenario.We also compute ground-state energy distributions and show that mean-field and non-mean-field models are intrinsically different. * )Recently Krzakala and Martin, as well as Palassini and Young, suggest an intermediate picture in which droplets have a fractal surface, and their energy is finite in the thermodynamic limit.See Ref. 5) and references therein for details.* * ) Note that we explicitly distinguish between θ , the stiffness exponent for droplet excitations, and θ, the exponent for domain walls, because DP predicts θ = θ , which we want to test here.

arXiv:cond-mat/0509469•Progress of Theoretical Physics Supplement•Jan 1, 2005•11 citations

Lower critical dimension of Ising spin glasses

Alexander K. Hartmann, A. P. Young

Exact ground states of two-dimensional Ising spin glasses with Gaussian and bimodal $(\ifmmode\pm\else\textpm\fi{}J)$ distributions of the disorder are calculated using a ``matching''algorithm, which allows large system sizes of up to ${N=480}^{2}$ spins to be investigated. We study domain walls induced by two rather different types of boundary condition changes, and, in each case, analyze the system-size dependence of an appropriately defined ``defect energy,'' which we denote by $\ensuremath{\Delta}E.$ For Gaussian disorder, we find a power-law behavior $\ensuremath{\Delta}E\ensuremath{\sim}{L}^{\ensuremath{\theta}},$ with $\ensuremath{\theta}=\ensuremath{-}0.266(2)$ and $\ensuremath{\theta}=\ensuremath{-}0.282(2)$ for the two types of boundary condition changes. These results are in reasonable agreement with each other, allowing for small systematic effects. They also agree well with earlier work on smaller sizes. The negative value indicates that two dimensions is below the lower critical dimension ${d}_{c}.$ For the $\ifmmode\pm\else\textpm\fi{}J$ model, we obtain a different result, namely that the domain-wall energy saturates at a nonzero value for $\stackrel{\ensuremath{\rightarrow}}{L}\ensuremath{\infty},$ so $\ensuremath{\theta}=0,$ indicating that the lower critical dimension for the $\ifmmode\pm\else\textpm\fi{}J$ model is exactly ${d}_{c}=2.$

arXiv:cond-mat/0107308•Physical review. B, Condensed matter•Oct 5, 2001•131 citations

One-dimensional Ising spin-glass with power-law interaction: real-space renormalization at zero temperature

Cécile Monthus

For the one-dimensional long-ranged Ising spin-glass with random couplings decaying with the distance r as J ( r ) ∼ r − σ and distributed with the Lévy symmetric stable distribution of index 1 < μ ≤ 2 (including the usual Gaussian case μ = 2), we consider the region σ > 1/ μ where the energy is extensive. We study two real space renormalization procedures at zero temperature, namely a simple box decimation that leads to explicit calculations, and a strong disorder decimation that can be studied numerically on large sizes. The droplet exponent governing the scaling of the renormalized couplings J L ∝ L θ μ ( σ ) is found to be whenever the long-ranged couplings are relevant θ μ ( σ ) ≥ −1. For the statistics of the ground state energy over disordered samples, we obtain that the droplet exponent θ μ ( σ ) governs the leading correction to extensivity of the averaged value . The characteristic scale of the fluctuations around this average is of order , and the rescaled variable is Gaussian distributed for μ = 2, or displays the negative power-law tail in 1/(− u ) 1+ μ for u → − ∞ in the Lévy case 1 < μ < 2. Finally we apply the zero-temperature renormalization procedure to the related Dyson hierarchical spin-glass model where the same droplet exponent appears.

arXiv:1403.1098•Journal of Statistical Mechanics Theory and Experiment•Jun 30, 2014•9 citations
Using mappings to computer-science problems and by applying sophisticated algorithms, one can study numerically many problems much better compared to applying standard approaches like Monte Carlo simulations. Here, using calculations of ground states of suitable perturbed systems, droplets are obtained in two-dimensional +-J spin glasses, which are in the focus of a currently very lifely debate. Since a sophisticated matching algorithm is applied here, exact ground states of large systems up to L^2=256^2 spins can be generated. Furthermore, no equilibration or extrapolation to T=0 is necessary. Three different +-J models are studied here: a) with open boundary conditions, b) with fixed boundary conditions and c) a diluted system where a fraction p=0.125 of all bonds is zero. For large systems, the droplet energy shows for all three models a power-law behavior E_D L^θ'_D with θ'_D<0. This is different from previous studies of domain walls, where a convergence to a constant non-zero value (θ_dw=0) has been found for such models. After correcting for the non-compactness of the droplets, the results are likely to be compatible with θ_D= -0.29 for all three models. This is in accordance with the Gaussian system where θ_D=-0.287(4) (ν=3.5 via ν=-1/θ_D). Nevertheless, the disorder-averaged spin-spin correlation exponent ηis determined here via the probability to have a non-zero-energy droplet, and η~0.22$ is found for all three models, this being in contrast to the behavior of the model with Gaussian interactions, where exactly η=0.
By using problem-dependent mappings to computer-science problems and by applying sophisticated algorithms, one can numerically study some important problems much better compared to applying standard approaches, such as Monte Carlo simulations. Here, by using calculations of ground states of suitable perturbed systems, droplets are obtained in two-dimensional $\ifmmode\pm\else\textpm\fi{}J$ spin glasses, which are in the focus of a currently very lively debate. By using a sophisticated matching algorithm, exact ground states of large systems up to L2=2562{L}^{2}={256}^{2}L2=2562 spins can be generated. Furthermore, no equilibration or extrapolation to T=0T=0T=0 is necessary. Three different $\ifmmode\pm\else\textpm\fi{}J$ models are studied here: (a) with free boundary conditions, (b) with fixed boundary conditions, and (c) a diluted system wherein a fraction p=0.125p=0.125p=0.125 of all bonds is zero. For large systems, the droplet energy shows for all three models a power-law behavior ${E}_{\text{dp}}\ensuremath{\sim}{L}^{{\ensuremath{\theta}}_{\text{dp}}^{\ensuremath{'}}}$ with ${\ensuremath{\theta}}_{\text{dp}}^{\ensuremath{'}}<0$. This is different from the previous studies of domain walls, wherein a convergence to a constant nonzero value $({\ensuremath{\theta}}_{\text{DW}}=0)$ was found. After correcting for the noncompactness of the droplets, the results are likely to be compatible with ${\ensuremath{\theta}}_{\text{dp}}\ensuremath{\approx}\ensuremath{-}0.29$ for all three models. This is in accordance with the Gaussian system wherein ${\ensuremath{\theta}}_{\text{dp}}=\ensuremath{-}0.287(4)$ ($\ensuremath{\nu}\ensuremath{\approx}3.5$ via $\ensuremath{\nu}=\ensuremath{-}1/{\ensuremath{\theta}}_{\text{dp}}$). Nevertheless, the disorder-averaged spin-spin correlation exponent $\ensuremath{\eta}$ is determined here via the probability to have a nonzero-energy droplet, and $\ensuremath{\eta}\ensuremath{\approx}0.22$ is found for all three models, in contrast to the model with the Gaussian interactions, where $\ensuremath{\eta}=0$ is exact.
We study large-scale, low-energy excitations in the Ising spin glass with Gaussian interactions in two dimensions at zero temperature, using an optimization algorithm to determine exact ground states. Periodic boundary conditions are applied. Our results for the fractal dimension of the surface, ds,{d}_{s},ds​, and stiffness exponent, ${\ensuremath{\theta}}^{\ensuremath{'}},$ for ``droplet'' excitations, are in reasonable agreement with estimates from ``domain wall'' calculations, and so support the predictions of the ``droplet theory.'' Restricting our analysis to small lattices, we do not find an effective value of ${\ensuremath{\theta}}^{\ensuremath{'}}$ close to $\ensuremath{-}0.47$ as has been recently proposed. The effects of averaging over droplets of different sizes are studied and are also found to be too small to give ${\ensuremath{\theta}}^{\ensuremath{'}}\ensuremath{\approx}\ensuremath{-}0.47$ for smaller sizes. Larger corrections to finite-size scaling would be needed in three and four dimensions in order for the numerical data used to support the ``TNT'' scenario to be compatible with the droplet theory prediction that the stiffness exponent is the same for droplets and domain walls.
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.
In order to overcome the limitations of small system sizes in spin-glass simulations, we investigate the one-dimensional Ising spin chain with power-law interactions.The model has the advantage over traditional higher-dimensional Hamiltonians in that a large range of system sizes can be studied.In addition, the universality class of the model can be changed by tuning the power law exponent, thus allowing us to scan from the mean-field to long-range and short-range universality classes.We illustrate the advantages of this model by studying the nature of the spin glass state where our results hint towards a replica symmetry breaking scenario.We also compute ground-state energy distributions and show that mean-field and non-mean-field models are intrinsically different. * )Recently Krzakala and Martin, as well as Palassini and Young, suggest an intermediate picture in which droplets have a fractal surface, and their energy is finite in the thermodynamic limit.See Ref. 5) and references therein for details.* * ) Note that we explicitly distinguish between θ , the stiffness exponent for droplet excitations, and θ, the exponent for domain walls, because DP predicts θ = θ , which we want to test here.
Exact ground states of two-dimensional Ising spin glasses with Gaussian and bimodal $(\ifmmode\pm\else\textpm\fi{}J)$ distributions of the disorder are calculated using a ``matching''algorithm, which allows large system sizes of up to N=4802{N=480}^{2}N=4802 spins to be investigated. We study domain walls induced by two rather different types of boundary condition changes, and, in each case, analyze the system-size dependence of an appropriately defined ``defect energy,'' which we denote by $\ensuremath{\Delta}E.$ For Gaussian disorder, we find a power-law behavior $\ensuremath{\Delta}E\ensuremath{\sim}{L}^{\ensuremath{\theta}},$ with $\ensuremath{\theta}=\ensuremath{-}0.266(2)$ and $\ensuremath{\theta}=\ensuremath{-}0.282(2)$ for the two types of boundary condition changes. These results are in reasonable agreement with each other, allowing for small systematic effects. They also agree well with earlier work on smaller sizes. The negative value indicates that two dimensions is below the lower critical dimension dc.{d}_{c}.dc​. For the $\ifmmode\pm\else\textpm\fi{}J$ model, we obtain a different result, namely that the domain-wall energy saturates at a nonzero value for $\stackrel{\ensuremath{\rightarrow}}{L}\ensuremath{\infty},$ so $\ensuremath{\theta}=0,$ indicating that the lower critical dimension for the $\ifmmode\pm\else\textpm\fi{}J$ model is exactly dc=2.{d}_{c}=2.dc​=2.
For the one-dimensional long-ranged Ising spin-glass with random couplings decaying with the distance r as J ( r ) ∼ r − σ and distributed with the Lévy symmetric stable distribution of index 1 < μ ≤ 2 (including the usual Gaussian case μ = 2), we consider the region σ > 1/ μ where the energy is extensive. We study two real space renormalization procedures at zero temperature, namely a simple box decimation that leads to explicit calculations, and a strong disorder decimation that can be studied numerically on large sizes. The droplet exponent governing the scaling of the renormalized couplings J L ∝ L θ μ ( σ ) is found to be whenever the long-ranged couplings are relevant θ μ ( σ ) ≥ −1. For the statistics of the ground state energy over disordered samples, we obtain that the droplet exponent θ μ ( σ ) governs the leading correction to extensivity of the averaged value . The characteristic scale of the fluctuations around this average is of order , and the rescaled variable is Gaussian distributed for μ = 2, or displays the negative power-law tail in 1/(− u ) 1+ μ for u → − ∞ in the Lévy case 1 < μ < 2. Finally we apply the zero-temperature renormalization procedure to the related Dyson hierarchical spin-glass model where the same droplet exponent appears.