arXiv:cond-mat/0208546 · Physical review. B, Condensed matter

Stiffness exponent of two-dimensional Ising spin glasses for nonperiodic boundary conditions using aspect-ratio scaling

Alexander K. Hartmann, Alan J. Bray, Adam C. Carter, M A Moore, A. P. Young

Published Dec 4, 2002• 52 citations

Abstract

We study the scaling behavior of domain-wall energies in two-dimensional Ising spin glasses with Gaussian and bimodal distributions of the interactions and different types of boundary conditions. The domain walls are generated by changing the boundary conditions at $T=0$ in one direction. The ground states of the original and perturbed system are calculated numerically by applying an efficient matching algorithm. Systems of size $L\ifmmode\times\else\texttimes\fi{}M$ with different aspect-ratios $1/8<~L/M<~64$ are considered. For Gaussian interactions, using the aspect-ratio scaling approach, we find a stiffness exponent $\ensuremath{\theta}=\ensuremath{-}0.287(4),$ which is independent of the boundary conditions in contrast to earlier results. Furthermore, we find a scaling behavior of the domain-wall energy as predicted by the aspect-ratio approach. Finally, we show that this approach does not work for the bimodal distribution of interactions.

DOI: 10.1103/physrevb.66.224401

Related papers

Energy size effects of two-dimensional Ising spin glasses

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

arXiv:cond-mat/0402627Physical Review BAug 31, 200435 citations

Lower critical dimension of Ising spin glasses

Alexander K. Hartmann, A. P. Young

arXiv:cond-mat/0107308Physical review. B, Condensed matterOct 5, 2001131 citations

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

Alexander K. Hartmann, A. P. Young

arXiv:cond-mat/0205659Physical review. B, Condensed matterSep 16, 200245 citations

Ground-state and domain-wall energies in the spin-glass region of the two-dimensional ±J random-bond Ising model

Ronald Fisch, Alexander K. Hartmann

arXiv:cond-mat/0603729Physical Review BMay 10, 20076 citations

Domain-wall free energy of spin-glass models: Numerical method and boundary conditions

Koji Hukushima

arXiv:cond-mat/9903391Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topicsOct 1, 199994 citations

Zero-temperature spin-glass–ferromagnetic transition: Scaling analysis of the domain-wall energy

Cécile Monthus, Thomas Garel

arXiv:1401.6342Physical Review BMay 13, 201413 citations

arXiv:cond-mat/0208546 · Physical review. B, Condensed matter

Stiffness exponent of two-dimensional Ising spin glasses for nonperiodic boundary conditions using aspect-ratio scaling

Alexander K. Hartmann, Alan J. Bray, Adam C. Carter, M A Moore, A. P. Young

Published Dec 4, 2002• 52 citations

arXiv HTML PDF

Abstract

T=0

in one direction. The ground states of the original and perturbed system are calculated numerically by applying an efficient matching algorithm. Systems of size

L\ifmmode\times\else\texttimes\fi{}M

with different aspect-ratios

1/8<~L/M<~64

are considered. For Gaussian interactions, using the aspect-ratio scaling approach, we find a stiffness exponent

\ensuremath{\theta}=\ensuremath{-}0.287(4),

which is independent of the boundary conditions in contrast to earlier results. Furthermore, we find a scaling behavior of the domain-wall energy as predicted by the aspect-ratio approach. Finally, we show that this approach does not work for the bimodal distribution of interactions.

DOI: 10.1103/physrevb.66.224401

Related papers

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.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/0402627Physical Review BAug 31, 200435 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.$ 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/0107308Physical review. B, Condensed matterOct 5, 2001131 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.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/0205659Physical review. B, Condensed matterSep 16, 200245 citations

Ground-state and domain-wall energies in the spin-glass region of the two-dimensional ±J random-bond Ising model

Ronald Fisch, Alexander K. Hartmann

The statistics of the ground-state and domain-wall energies for the two-dimensional random-bond Ising model on square lattices with independent, identically distributed bonds having probability $p$ of ${J}_{ij}=\ensuremath{-}1$ and $(1\ensuremath{-}p)$ of ${J}_{ij}=+1$ are studied. We are able to consider large samples of up to ${320}^{2}$ spins by using sophisticated matching algorithms. We study $L\ifmmode\times\else\texttimes\fi{}L$ systems, but we also consider $L\ifmmode\times\else\texttimes\fi{}M$ samples, for different aspect ratios $R=L∕M$ . We find that the scaling behavior of the ground-state energy and its sample-to-sample fluctuations inside the spin-glass region $({p}_{c}\ensuremath{\leqslant}p\ensuremath{\leqslant}1\ensuremath{-}{p}_{c})$ are characterized by simple scaling functions. In particular, the fluctuations exhibit a cusplike singularity at ${p}_{c}$ . Inside the spin-glass region, the average domain-wall energy converges to a finite nonzero value as the sample size becomes infinite, holding $R$ fixed. Here, large finite-size effects are visible, which can be explained for all $p$ by a single exponent $\ensuremath{\omega}\ensuremath{\approx}2∕3$, provided higher-order corrections to scaling are included. Finally, we confirm the validity of aspect-ratio scaling for $R\ensuremath{\rightarrow}0$: the distribution of the domain-wall energies converges to a Gaussian for $R\ensuremath{\rightarrow}0$, although the domain walls of neighboring subsystems of size $L\ifmmode\times\else\texttimes\fi{}L$ are not independent.The statistics of the ground-state and domain-wall energies for the two-dimensional random-bond Ising model on square lattices with independent, identically distributed bonds having probability $p$ of ${J}_{ij}=\ensuremath{-}1$ and $(1\ensuremath{-}p)$ of ${J}_{ij}=+1$ are studied. We are able to consider large samples of up to ${320}^{2}$ spins by using sophisticated matching algorithms. We study $L\ifmmode\times\else\texttimes\fi{}L$ systems, but we also consider $L\ifmmode\times\else\texttimes\fi{}M$ samples, for different aspect ratios $R=L∕M$. We find that the scaling behavior of the ground-state energy and its sample-to-sample fluctuations inside the spin-glass region $({p}_{c}\ensuremath{\leqslant}p\ensuremath{\leqslant}1\ensuremath{-}{p}_{c})$ are characterized by simple scaling functions. In particular, the fluctuations exhibit a cusplike singularity at ${p}_{c}$. Inside the spin-glass region, the average domain-wall energy converges to a finite nonzero value as the sample size becomes infinite, holding $R$ fixed. Here, large finite-size effects are visible, which can be explained for all $p$ by a single exponent $\ensuremath{\omega}\ensuremath{\approx}2∕3$, provided higher-order corrections to scaling are included. Finally, we confirm the validity of aspect-ratio scaling for $R\ensuremath{\rightarrow}0$: the distribution of the domain-wall energies converges to a Gaussian for $R\ensuremath{\rightarrow}0$, although the domain walls of neighboring subsystems of size $L\ifmmode\times\else\texttimes\fi{}L$ are not independent.

arXiv:cond-mat/0603729Physical Review BMay 10, 20076 citations

Domain-wall free energy of spin-glass models: Numerical method and boundary conditions

Koji Hukushima

An efficient Monte Carlo method is extended to evaluate directly domain-wall free energy for randomly frustrated spin systems. Using the method, critical phenomena of spin-glass phase transition are investigated in the 4d+/-J Ising model under the replica boundary condition. Our values of the critical temperature and exponent, obtained by finite-size scaling, are in good agreement with those of the standard Monte Carlo and the series expansion studies. In addition, two exponents, the stiffness exponent and the fractal dimension of the domain wall, which characterize the ordered phase, are obtained. The latter value is larger than d-1, indicating that the domain wall is really rough in the 4d Ising spin-glass phase.An efficient Monte Carlo method is extended to evaluate directly domain-wall free energy for randomly frustrated spin systems. Using the method, critical phenomena of spin-glass phase transition are investigated in the 4d+/-J Ising model under the replica boundary condition. Our values of the critical temperature and exponent, obtained by finite-size scaling, are in good agreement with those of the standard Monte Carlo and the series expansion studies. In addition, two exponents, the stiffness exponent and the fractal dimension of the domain wall, which characterize the ordered phase, are obtained. The latter value is larger than d-1, indicating that the domain wall is really rough in the 4d Ising spin-glass phase.

arXiv:cond-mat/9903391Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topicsOct 1, 199994 citations

Zero-temperature spin-glass–ferromagnetic transition: Scaling analysis of the domain-wall energy

Cécile Monthus, Thomas Garel

For the Ising model with Gaussian random coupling of average ${J}_{0}$ and unit variance, the zero-temperature spin-glass--ferromagnetic transition as a function of the control parameter ${J}_{0}$ can be studied via the size- $L$ -dependent renormalized coupling defined as the domain-wall energy ${J}^{R}(L)\ensuremath{\equiv}{E}_{\text{GS}}^{(\text{AF})}(L)\ensuremath{-}{E}_{\text{GS}}^{(\text{F})}(L)$ (i.e., the difference between the ground-state energies corresponding to antiferromagnetic and ferromagnetic boundary conditions in one direction). We study numerically the critical exponents of this zero-temperature transition within the Migdal-Kadanoff approximation as a function of the dimension $d=2,3,4,5,6$ . We then compare with the mean-field spherical model. Our main conclusion is that in low dimensions, the critical stiffness exponent ${\ensuremath{\theta}}^{c}$ is clearly bigger than the spin-glass stiffness exponent ${\ensuremath{\theta}}^{\text{SG}}$, but that they turn out to coincide in high enough dimension and in the mean-field spherical model. We also discuss the finite-size scaling properties of the averaged value and of the width of the distribution of the renormalized couplings.For the Ising model with Gaussian random coupling of average ${J}_{0}$ and unit variance, the zero-temperature spin-glass--ferromagnetic transition as a function of the control parameter ${J}_{0}$ can be studied via the size-$L$-dependent renormalized coupling defined as the domain-wall energy ${J}^{R}(L)\ensuremath{\equiv}{E}_{\text{GS}}^{(\text{AF})}(L)\ensuremath{-}{E}_{\text{GS}}^{(\text{F})}(L)$ (i.e., the difference between the ground-state energies corresponding to antiferromagnetic and ferromagnetic boundary conditions in one direction). We study numerically the critical exponents of this zero-temperature transition within the Migdal-Kadanoff approximation as a function of the dimension $d=2,3,4,5,6$. We then compare with the mean-field spherical model. Our main conclusion is that in low dimensions, the critical stiffness exponent ${\ensuremath{\theta}}^{c}$ is clearly bigger than the spin-glass stiffness exponent ${\ensuremath{\theta}}^{\text{SG}}$, but that they turn out to coincide in high enough dimension and in the mean-field spherical model. We also discuss the finite-size scaling properties of the averaged value and of the width of the distribution of the renormalized couplings.

arXiv:1401.6342Physical Review BMay 13, 201413 citations

T=0

in one direction. The ground states of the original and perturbed system are calculated numerically by applying an efficient matching algorithm. Systems of size

L\ifmmode\times\else\texttimes\fi{}M

with different aspect-ratios

1/8<~L/M<~64

are considered. For Gaussian interactions, using the aspect-ratio scaling approach, we find a stiffness exponent

\ensuremath{\theta}=\ensuremath{-}0.287(4),

{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.

{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.

The statistics of the ground-state and domain-wall energies for the two-dimensional random-bond Ising model on square lattices with independent, identically distributed bonds having probability

p

of ${J}_{ij}=\ensuremath{-}1$ and $(1\ensuremath{-}p)$ of

{J}_{ij}=+1

are studied. We are able to consider large samples of up to

{320}^{2}

spins by using sophisticated matching algorithms. We study $L\ifmmode\times\else\texttimes\fi{}L$ systems, but we also consider $L\ifmmode\times\else\texttimes\fi{}M$ samples, for different aspect ratios

R=L∕M

. We find that the scaling behavior of the ground-state energy and its sample-to-sample fluctuations inside the spin-glass region $({p}_{c}\ensuremath{\leqslant}p\ensuremath{\leqslant}1\ensuremath{-}{p}_{c})$ are characterized by simple scaling functions. In particular, the fluctuations exhibit a cusplike singularity at

{p}_{c}

. Inside the spin-glass region, the average domain-wall energy converges to a finite nonzero value as the sample size becomes infinite, holding

R

fixed. Here, large finite-size effects are visible, which can be explained for all

p

by a single exponent $\ensuremath{\omega}\ensuremath{\approx}2∕3$, provided higher-order corrections to scaling are included. Finally, we confirm the validity of aspect-ratio scaling for $R\ensuremath{\rightarrow}0$: the distribution of the domain-wall energies converges to a Gaussian for $R\ensuremath{\rightarrow}0$, although the domain walls of neighboring subsystems of size $L\ifmmode\times\else\texttimes\fi{}L$ are not independent.

For the Ising model with Gaussian random coupling of average

{J}_{0}

and unit variance, the zero-temperature spin-glass--ferromagnetic transition as a function of the control parameter

{J}_{0}

can be studied via the size-

L

-dependent renormalized coupling defined as the domain-wall energy ${J}^{R}(L)\ensuremath{\equiv}{E}_{\text{GS}}^{(\text{AF})}(L)\ensuremath{-}{E}_{\text{GS}}^{(\text{F})}(L)$ (i.e., the difference between the ground-state energies corresponding to antiferromagnetic and ferromagnetic boundary conditions in one direction). We study numerically the critical exponents of this zero-temperature transition within the Migdal-Kadanoff approximation as a function of the dimension

d=2,3,4,5,6

. We then compare with the mean-field spherical model. Our main conclusion is that in low dimensions, the critical stiffness exponent ${\ensuremath{\theta}}^{c}$ is clearly bigger than the spin-glass stiffness exponent ${\ensuremath{\theta}}^{\text{SG}}$, but that they turn out to coincide in high enough dimension and in the mean-field spherical model. We also discuss the finite-size scaling properties of the averaged value and of the width of the distribution of the renormalized couplings.