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

Droplets in the two-dimensional ±J Ising spin glass

Alexander K. Hartmann

Published Apr 18, 2008• 20 citations
Publisher

Abstract

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.

DOI
10.1103/physrevb.77.144418

Related papers

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

Alexander K. Hartmann

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.

arXiv:0704.2748•Apr 20, 2007•0 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

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

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

Energy size effects of two-dimensional Ising spin glasses

Ian A. Campbell, Alexander K. Hartmann, Helmut G. Katzgraber

We analyze exact ground-state energies of two-dimensional Ising spin glasses with either Gaussian or bimodal nearest-neighbor interactions for large system sizes and for three types of boundary conditions: free on both axes, periodic on both axes, and free on one axis and periodic on the other. We find accurate values for bulk-, edge-, and corner-site energies. Fits for the system with Gaussian bonds are excellent for all types of boundary conditions over the whole range of system sizes. In particular, the leading behavior for nonfree boundary conditions is governed by the stiffness exponent $\ensuremath{\theta}\ensuremath{\approx}\ensuremath{-}0.282$ describing the scaling of domain-wall and droplet excitations. For the system with a bimodal distribution of bonds the fit is good for free boundary conditions but worse for other geometries, particularly for periodic-free boundary conditions where there appear to be unorthodox corrections to scaling up to large sizes. Finally, by introducing hard bonds we test explicitly for the Gaussian case the relationship between domain walls and the standard scaling behavior.

arXiv:cond-mat/0402627•Physical Review B•Aug 31, 2004•35 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
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±\else\textpm\fiJ\ifmmode\pm\else\textpm\fi{}J\ifmmode±\else\textpm\fiJ 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±\else\textpm\fiJ\ifmmode\pm\else\textpm\fi{}J\ifmmode±\else\textpm\fiJ 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 Edp\ensuremath∼L\ensuremathθdp\ensuremath′{E}_{\text{dp}}\ensuremath{\sim}{L}^{{\ensuremath{\theta}}_{\text{dp}}^{\ensuremath{'}}}Edp​\ensuremath∼L\ensuremathθdp\ensuremath′​ with \ensuremathθdp\ensuremath′<0{\ensuremath{\theta}}_{\text{dp}}^{\ensuremath{'}}<0\ensuremathθdp\ensuremath′​<0. This is different from the previous studies of domain walls, wherein a convergence to a constant nonzero value (\ensuremathθDW=0)({\ensuremath{\theta}}_{\text{DW}}=0)(\ensuremathθDW​=0) was found. After correcting for the noncompactness of the droplets, the results are likely to be compatible with \ensuremathθdp\ensuremath≈\ensuremath−0.29{\ensuremath{\theta}}_{\text{dp}}\ensuremath{\approx}\ensuremath{-}0.29\ensuremathθdp​\ensuremath≈\ensuremath−0.29 for all three models. This is in accordance with the Gaussian system wherein \ensuremathθdp=\ensuremath−0.287(4){\ensuremath{\theta}}_{\text{dp}}=\ensuremath{-}0.287(4)\ensuremathθdp​=\ensuremath−0.287(4) (\ensuremathν\ensuremath≈3.5\ensuremath{\nu}\ensuremath{\approx}3.5\ensuremathν\ensuremath≈3.5 via \ensuremathν=\ensuremath−1/\ensuremathθdp\ensuremath{\nu}=\ensuremath{-}1/{\ensuremath{\theta}}_{\text{dp}}\ensuremathν=\ensuremath−1/\ensuremathθdp​). Nevertheless, the disorder-averaged spin-spin correlation exponent \ensuremathη\ensuremath{\eta}\ensuremathη is determined here via the probability to have a nonzero-energy droplet, and \ensuremathη\ensuremath≈0.22\ensuremath{\eta}\ensuremath{\approx}0.22\ensuremathη\ensuremath≈0.22 is found for all three models, in contrast to the model with the Gaussian interactions, where \ensuremathη=0\ensuremath{\eta}=0\ensuremathη=0 is exact.
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.
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.
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.
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.
We analyze exact ground-state energies of two-dimensional Ising spin glasses with either Gaussian or bimodal nearest-neighbor interactions for large system sizes and for three types of boundary conditions: free on both axes, periodic on both axes, and free on one axis and periodic on the other. We find accurate values for bulk-, edge-, and corner-site energies. Fits for the system with Gaussian bonds are excellent for all types of boundary conditions over the whole range of system sizes. In particular, the leading behavior for nonfree boundary conditions is governed by the stiffness exponent $\ensuremath{\theta}\ensuremath{\approx}\ensuremath{-}0.282$ describing the scaling of domain-wall and droplet excitations. For the system with a bimodal distribution of bonds the fit is good for free boundary conditions but worse for other geometries, particularly for periodic-free boundary conditions where there appear to be unorthodox corrections to scaling up to large sizes. Finally, by introducing hard bonds we test explicitly for the Gaussian case the relationship between domain walls and the standard scaling behavior.
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.