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, Condensed matter

Critical behavior of the three-dimensional Ising spin glass

Alan J. Bray, M A Moore

Published Jan 1, 1985• 141 citations
Publisher

Abstract

The three-dimensional nearest-neighbor Ising spin-glass model with Gaussian bonds (〈${J}^{2}$〉=1) is studied numerically using the ``large-cell'' renormalization group. Estimates for the critical coupling and the thermal exponent are ${K}_{c}$=1/${T}_{c}$=1.2\ifmmode\pm\else\textpm\fi{}0.1 and ${y}_{t}$=1/\ensuremath{\nu}=0.30\ifmmode\pm\else\textpm\fi{}0.05, respectively.

DOI
10.1103/physrevb.31.631

Related papers

Critical behavior of three-dimensional Ising spin glass models

Martin Hasenbusch B, Andrea Pelissetto C, Ettore Vicari

We perform high-statistics Monte Carlo simulations of three-dimensional Ising spin glass models on cubic lattices of size $L$: the $\ifmmode\pm\else\textpm\fi{}J$ (Edwards-Anderson) Ising model for two values of the disorder parameter $p$, $p=0.5$ and $p=0.7$ (up to $L=28$ and $L=20$, respectively), and the bond-diluted bimodal model for bond-occupation probability ${p}_{b}=0.45$ (up to $L=16$). The finite-size behavior of the quartic cumulants at the critical point allows us to check very accurately that these models belong to the same universality class. Moreover, it allows us to estimate the scaling-correction exponent $\ensuremath{\omega}$ related to the leading irrelevant operator: $\ensuremath{\omega}=1.0(1)$. Shorter Monte Carlo simulations of the bond-diluted bimodal models at ${p}_{b}=0.7$ and ${p}_{b}=0.35$ (up to $L=10$) and of the Ising spin glass model with Gaussian bond distribution (up to $L=8$) also support the existence of a unique Ising spin glass universality class. A careful finite-size analysis of the Monte Carlo data which takes into account the analytic and the nonanalytic corrections to scaling allows us to obtain precise and reliable estimates of the critical exponents. We obtain $\ensuremath{\nu}=2.45(15)$ and $\ensuremath{\eta}=\ensuremath{-}0.375(10)$.

arXiv:0809.3329•Physical Review B•Dec 17, 2008•119 citations

The critical behavior of 3D Ising spin glass models: universality and scaling corrections

Martin Hasenbusch, Andrea Pelissetto, Ettore Vicari

We perform high-statistics Monte Carlo simulations of three three-dimensional Ising spin glass models: the ± J Ising model for two values of the disorder parameter p , p = 1/2 and 0.7, and the bond-diluted ± J model for bond-occupation probability p b = 0.45. A finite-size scaling analysis of the quartic cumulants at the critical point shows conclusively that these models belong to the same universality class and allows us to estimate the scaling-correction exponent ω related to the leading irrelevant operator, ω = 1.0(1). We also determine the critical exponents ν and η. Taking into account the scaling corrections, we obtain ν = 2.53(8) and η = −0.384(9).

arXiv:0710.1980•Journal of Statistical Mechanics Theory and Experiment•Feb 1, 2008•50 citations

Domain-wall renormalization-group study of the three-dimensional random Ising model

W. L. McMillan

The three-dimensional random Ising model with a Gaussian distribution of nearest-neighbor interactions is studied for the pure spin-glass case where the average interaction vanishes. The distribution of domain-wall energies at zero temperature is calculated using a Monte Carlo-quench algorithm to find the ground-state energy for finite lattices. A renormalization-group transformation is set up which preserves the domain-wall energy distribution when the lattice parameter is changed. In the strong-coupling regime (zero temperature) the model iterates toward strong coupling and therefore exhibits a spin-glass phase transition at nonzero temperature. The thermal exponent is $\ensuremath{\nu}=3.0\ifmmode\pm\else\textpm\fi{}1.0$ and the heat capacity exponent is $\ensuremath{\alpha}=\ensuremath{-}7\ifmmode\pm\else\textpm\fi{}3$.

Physical review. B, Condensed matter•Jul 1, 1984•185 citations

Domain-wall renormalization-group study of the three-dimensional random Ising model at finite temperature

W. L. McMillan

The three-dimensional random Ising model with a Gaussian distribution of nearest-neighbor interactions is studied for the pure spin-glass case where the average interaction vanishes. The distribution of domain-wall free energies at finite temperature is calculated with the use of the Metropolis Monte Carlo algorithm for finite lattices. A renormalization-group transformation is set up which preserves the domain-wall free-energy distribution when the lattice parameter is changed. The spin-glass transition temperature is found to be ${T}_{g}$=(1.0\ifmmode\pm\else\textpm\fi{}0.2)J\ifmmode \tilde{}\else \~{}\fi{} where J\ifmmode \tilde{}\else \~{}\fi{} is the variance of the Gaussian interaction distribution. The thermal exponent is \ensuremath{\nu}=1.8\ifmmode\pm\else\textpm\fi{}0.5 and the heat-capacity exponent is \ensuremath{\alpha}=-3.4\ifmmode\pm\else\textpm\fi{}1.5. The heat capacity exhibits a rounded peak at higher temperatures.

Physical review. B, Condensed matter•Jan 1, 1985•109 citations

Search for a transition in the three-dimensional J Ising spin-glass

R. N. Bhatt, A. P. Young

The three-dimensional Ising spin-glass in zero field with nearest-neighbor interactions having \ifmmode\pm\else\textpm\fi{}J distribution is studied by Monte Carlo simulations for samples of linear dimension L with 3\ensuremath{\le}L\ensuremath{\le}20. Results for the probability distribution, ${\mathrm{P}}_{\mathrm{L}}$(q), of the spin-glass order parameter are analyzed by finite-size scaling. Data for T\ensuremath{\ge}1.2 are consistent with a conventional phase transition at ${\mathrm{T}}_{\mathrm{c}}$${=1.2}_{\mathrm{\ensuremath{-}}0.2}^{+0.1}$, with exponents \ensuremath{\nu}=1.3\ifmmode\pm\else\textpm\fi{}0.3 and \ensuremath{\eta}=-0.3\ifmmode\pm\else\textpm\fi{}0.2. However, results at lower temperature indicate marginal behavior and suggest that the lower critical dimensional is close to three.

Physical Review Letters•Mar 4, 1985•392 citations

Numerical studies of Ising spin glasses in two, three, and four dimensions

R. N. Bhatt, A. P. Young

We present the results of numerical simulations on Ising spin glasses in zero magnetic field with nearest-neighbor interactions on hyper-cubic lattices in two, three, and four dimensions with both Gaussian and \ifmmode\pm\else\textpm\fi{}J bond distributions. Finite-size scaling is used to analyze the results. In two dimensions (d=2) we agree with earlier work that the transition temperature is at ${T}_{c}$=0, and obtain the correlation-length exponent \ensuremath{\nu}, and the exponent \ensuremath{\eta}, at the zero-temperature transition for the \ifmmode\pm\else\textpm\fi{}J model. In d=3 dimensions we concentrate on results for the Gaussian distribution, since our results for the \ifmmode\pm\else\textpm\fi{}J distribution have been presented earlier. As expected, we find similar results for the two distributions, namely a nonzero ${T}_{c}$ but evidence that d=3 is close to the lower critical dimension. In a four-dimensional spin glass with Gaussian bonds we find that only a modest amount of computer time is required to show that ${T}_{c}$ is nonzero with a long-range-ordered phase below ${T}_{c}$. Our estimates for critical exponents in d=4 dimensions agree well with results from recent high-temperature-series expansions.

Physical review. B, Condensed matter•Apr 1, 1988•301 citations
The three-dimensional nearest-neighbor Ising spin-glass model with Gaussian bonds (〈J2{J}^{2}J2〉=1) is studied numerically using the ``large-cell'' renormalization group. Estimates for the critical coupling and the thermal exponent are Kc{K}_{c}Kc​=1/Tc{T}_{c}Tc​=1.2\ifmmode\pm\else\textpm\fi{}0.1 and yt{y}_{t}yt​=1/\ensuremath{\nu}=0.30\ifmmode\pm\else\textpm\fi{}0.05, respectively.
We perform high-statistics Monte Carlo simulations of three-dimensional Ising spin glass models on cubic lattices of size LLL: the $\ifmmode\pm\else\textpm\fi{}J$ (Edwards-Anderson) Ising model for two values of the disorder parameter ppp, p=0.5p=0.5p=0.5 and p=0.7p=0.7p=0.7 (up to L=28L=28L=28 and L=20L=20L=20, respectively), and the bond-diluted bimodal model for bond-occupation probability pb=0.45{p}_{b}=0.45pb​=0.45 (up to L=16L=16L=16). The finite-size behavior of the quartic cumulants at the critical point allows us to check very accurately that these models belong to the same universality class. Moreover, it allows us to estimate the scaling-correction exponent $\ensuremath{\omega}$ related to the leading irrelevant operator: $\ensuremath{\omega}=1.0(1)$. Shorter Monte Carlo simulations of the bond-diluted bimodal models at pb=0.7{p}_{b}=0.7pb​=0.7 and pb=0.35{p}_{b}=0.35pb​=0.35 (up to L=10L=10L=10) and of the Ising spin glass model with Gaussian bond distribution (up to L=8L=8L=8) also support the existence of a unique Ising spin glass universality class. A careful finite-size analysis of the Monte Carlo data which takes into account the analytic and the nonanalytic corrections to scaling allows us to obtain precise and reliable estimates of the critical exponents. We obtain $\ensuremath{\nu}=2.45(15)$ and $\ensuremath{\eta}=\ensuremath{-}0.375(10)$.
We perform high-statistics Monte Carlo simulations of three three-dimensional Ising spin glass models: the ± J Ising model for two values of the disorder parameter p , p = 1/2 and 0.7, and the bond-diluted ± J model for bond-occupation probability p b = 0.45. A finite-size scaling analysis of the quartic cumulants at the critical point shows conclusively that these models belong to the same universality class and allows us to estimate the scaling-correction exponent ω related to the leading irrelevant operator, ω = 1.0(1). We also determine the critical exponents ν and η. Taking into account the scaling corrections, we obtain ν = 2.53(8) and η = −0.384(9).
The three-dimensional random Ising model with a Gaussian distribution of nearest-neighbor interactions is studied for the pure spin-glass case where the average interaction vanishes. The distribution of domain-wall energies at zero temperature is calculated using a Monte Carlo-quench algorithm to find the ground-state energy for finite lattices. A renormalization-group transformation is set up which preserves the domain-wall energy distribution when the lattice parameter is changed. In the strong-coupling regime (zero temperature) the model iterates toward strong coupling and therefore exhibits a spin-glass phase transition at nonzero temperature. The thermal exponent is $\ensuremath{\nu}=3.0\ifmmode\pm\else\textpm\fi{}1.0$ and the heat capacity exponent is $\ensuremath{\alpha}=\ensuremath{-}7\ifmmode\pm\else\textpm\fi{}3$.
The three-dimensional random Ising model with a Gaussian distribution of nearest-neighbor interactions is studied for the pure spin-glass case where the average interaction vanishes. The distribution of domain-wall free energies at finite temperature is calculated with the use of the Metropolis Monte Carlo algorithm for finite lattices. A renormalization-group transformation is set up which preserves the domain-wall free-energy distribution when the lattice parameter is changed. The spin-glass transition temperature is found to be Tg{T}_{g}Tg​=(1.0\ifmmode\pm\else\textpm\fi{}0.2)J\ifmmode \tilde{}\else \~{}\fi{} where J\ifmmode \tilde{}\else \~{}\fi{} is the variance of the Gaussian interaction distribution. The thermal exponent is \ensuremath{\nu}=1.8\ifmmode\pm\else\textpm\fi{}0.5 and the heat-capacity exponent is \ensuremath{\alpha}=-3.4\ifmmode\pm\else\textpm\fi{}1.5. The heat capacity exhibits a rounded peak at higher temperatures.
The three-dimensional Ising spin-glass in zero field with nearest-neighbor interactions having \ifmmode\pm\else\textpm\fi{}J distribution is studied by Monte Carlo simulations for samples of linear dimension L with 3\ensuremath{\le}L\ensuremath{\le}20. Results for the probability distribution, PL{\mathrm{P}}_{\mathrm{L}}PL​(q), of the spin-glass order parameter are analyzed by finite-size scaling. Data for T\ensuremath{\ge}1.2 are consistent with a conventional phase transition at Tc{\mathrm{T}}_{\mathrm{c}}Tc​${=1.2}_{\mathrm{\ensuremath{-}}0.2}^{+0.1}$, with exponents \ensuremath{\nu}=1.3\ifmmode\pm\else\textpm\fi{}0.3 and \ensuremath{\eta}=-0.3\ifmmode\pm\else\textpm\fi{}0.2. However, results at lower temperature indicate marginal behavior and suggest that the lower critical dimensional is close to three.
We present the results of numerical simulations on Ising spin glasses in zero magnetic field with nearest-neighbor interactions on hyper-cubic lattices in two, three, and four dimensions with both Gaussian and \ifmmode\pm\else\textpm\fi{}J bond distributions. Finite-size scaling is used to analyze the results. In two dimensions (d=2) we agree with earlier work that the transition temperature is at Tc{T}_{c}Tc​=0, and obtain the correlation-length exponent \ensuremath{\nu}, and the exponent \ensuremath{\eta}, at the zero-temperature transition for the \ifmmode\pm\else\textpm\fi{}J model. In d=3 dimensions we concentrate on results for the Gaussian distribution, since our results for the \ifmmode\pm\else\textpm\fi{}J distribution have been presented earlier. As expected, we find similar results for the two distributions, namely a nonzero Tc{T}_{c}Tc​ but evidence that d=3 is close to the lower critical dimension. In a four-dimensional spin glass with Gaussian bonds we find that only a modest amount of computer time is required to show that Tc{T}_{c}Tc​ is nonzero with a long-range-ordered phase below Tc{T}_{c}Tc​. Our estimates for critical exponents in d=4 dimensions agree well with results from recent high-temperature-series expansions.