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

Numerical simulations on the 4d Heisenberg spin glass

Barbara Coluzzi B, Preprint Roma N

Published Jan 1, 1994• 19 citations
OpenAlex

Abstract

We study the 4d Heisenberg spin glass model with Gaussian nearestneighbor interactions. We use finite size scaling to analyze the data. We find a behavior consistent with a finite temperature spin glass transition. Our estimates for the critical exponents agree with the results from "-expansion. PACS Numbers 7510N-0270 Preprint Roma1 n o 1054 The lower critical dimension d l of the short range models remains one of the most controversial questions in the spin glass [1]-[3] theory. In the last ten years a consistent number of works has been devoted to studying the Ising model, which seems to be [4, 5] very close to d l in d = 3. In the case of the short range isotropic Heisenberg spin glass, the conclusions of various computer simulations in d = 3 [6, 7] agree that the system is below d l . Using domain wall renormalization group techniques, it was argued [2, 3] that d l = 4 for this model. To our knowledge, there are no previous numerical simulations on the 4d Heisenberg spin glass...

Related papers

Numerical simulations on the 4D Heisenberg spin glass

Barbara Coluzzi B

We study the 4D Heisenberg spin-glass model with Gaussian nearest-neighbour interactions. We use finite-size scaling to analyse the data. We find a behaviour consistent with a finite temperature spin-glass transition. Our estimates for the critical exponents agree with the results from the epsilon -expansion.

arXiv:cond-mat/9410053•Journal of Physics A Mathematical and General•Feb 7, 1995•21 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 Phase Transitions in a ppp spin Glass Model: A Numerical Study

Prerak Gupta, Auditya Sharma, Bharadwaj Vedula +2 more

We investigate the balanced $M=4$, $p=4$ spin-glass model for a one-dimensional long-range proxy for the finite dimensional short-range $p$-spin glass model to examine the nature of the glass transition beyond mean-field theory. We perform large-scale Monte Carlo equilibrated simulations for both fully connected and power-law diluted versions of the model. The critical temperatures extracted from the finite-size scaling (FSS) analysis of spin-glass susceptibility are in good agreement with theoretical predictions for $σ= 0, 0.25$, and 0.55. For these values of the long-range exponent $σ$ (which is the power of the decrease of the interactions between the spins with their separation), one might have expected that mean-field theory would provide a good description of the system. However, the spin-overlap distribution and the value of the $λ$-parameter do not provide numerical evidence for a one-step replica symmetry breaking (1RSB) phase transition. Instead, our results indicate a direct transition from the paramagnetic state to a full replica symmetry broken phase, with a renormalized value of $λ\equiv ω_2/ω_1 < 1$ suggesting a continuous FRSB transition, despite this ratio being equal to 2 at mean-field level. A value of $λ> 1$ is required for the discontinuous 1RSB transition. We argue that strong finite-size effects and closely spaced transition temperatures remove the expected 1RSB transition for the system sizes which we can study. For values of the exponent $σ= 0.85$, which roughly corresponds to a three dimensional system, we find that the renormalized value of $λ$ is again less than 1, with no signs of either the 1RSB transition or the continuous FRSB transition, suggesting that the Kauzmann temperature $T_K$ in three dimensions might be zero and the complete absence of phase transitions in structural glasses.

arXiv:2604.15433•Apr 16, 2026•0 citations

The Ising Spin Glass in dimension four; non-universality

P. H. Lundow, Ian A. Campbell

Extensive simulations are made on Ising Spin Glasses (ISG) with Gaussian, Laplacian and bimodal interaction distributions in dimension four. Standard finite size scaling analyses near and at criticality provide estimates of the critical inverse temperatures $β_c$, critical exponents, and critical values of a number of dimensionless parameters. Independent estimates are obtained for $β_c$ and the exponent $ν$ from thermodynamic derivative peak data. A detailed explanation is given of scaling in the thermodynamic limit with the ISG scaling variable $τ= 1-β^2/β_c^2$ and the appropriate scaling expressions. Data over the entire paramagnetic range of temperatures are analysed in order to obtain further estimates of the critical exponents together with correction to scaling terms. The Privman-Fisher ansatz then leads to compact scaling expressions for the whole paramagnetic regime and for all sample sizes L. Comparisons between the 4d ISG models show that the critical dimensionless parameters characteristic of a universality class, and the susceptibility and correlation length critical exponents $γ$ and $ν$, depend on the form of the interaction distribution. From these observations it can be deduced that critical exponents are not universal in ISGs, at least in dimension four.

arXiv:1402.1991•Feb 9, 2014•1 citation

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

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.

arXiv:cond-mat/9903391•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Oct 1, 1999•94 citations
We study the 4d Heisenberg spin glass model with Gaussian nearestneighbor interactions. We use finite size scaling to analyze the data. We find a behavior consistent with a finite temperature spin glass transition. Our estimates for the critical exponents agree with the results from "-expansion. PACS Numbers 7510N-0270 Preprint Roma1 n o 1054 The lower critical dimension d l of the short range models remains one of the most controversial questions in the spin glass [1]-[3] theory. In the last ten years a consistent number of works has been devoted to studying the Ising model, which seems to be [4, 5] very close to d l in d = 3. In the case of the short range isotropic Heisenberg spin glass, the conclusions of various computer simulations in d = 3 [6, 7] agree that the system is below d l . Using domain wall renormalization group techniques, it was argued [2, 3] that d l = 4 for this model. To our knowledge, there are no previous numerical simulations on the 4d Heisenberg spin glass...
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.
We investigate the balanced M=4M=4M=4, p=4p=4p=4 spin-glass model for a one-dimensional long-range proxy for the finite dimensional short-range ppp-spin glass model to examine the nature of the glass transition beyond mean-field theory. We perform large-scale Monte Carlo equilibrated simulations for both fully connected and power-law diluted versions of the model. The critical temperatures extracted from the finite-size scaling (FSS) analysis of spin-glass susceptibility are in good agreement with theoretical predictions for σ=0,0.25σ= 0, 0.25σ=0,0.25, and 0.55. For these values of the long-range exponent σσσ (which is the power of the decrease of the interactions between the spins with their separation), one might have expected that mean-field theory would provide a good description of the system. However, the spin-overlap distribution and the value of the λλλ-parameter do not provide numerical evidence for a one-step replica symmetry breaking (1RSB) phase transition. Instead, our results indicate a direct transition from the paramagnetic state to a full replica symmetry broken phase, with a renormalized value of λ≡ω2/ω1<1λ\equiv ω_2/ω_1 < 1λ≡ω2​/ω1​<1 suggesting a continuous FRSB transition, despite this ratio being equal to 2 at mean-field level. A value of λ>1λ> 1λ>1 is required for the discontinuous 1RSB transition. We argue that strong finite-size effects and closely spaced transition temperatures remove the expected 1RSB transition for the system sizes which we can study. For values of the exponent σ=0.85σ= 0.85σ=0.85, which roughly corresponds to a three dimensional system, we find that the renormalized value of λλλ is again less than 1, with no signs of either the 1RSB transition or the continuous FRSB transition, suggesting that the Kauzmann temperature TKT_KTK​ in three dimensions might be zero and the complete absence of phase transitions in structural glasses.
Extensive simulations are made on Ising Spin Glasses (ISG) with Gaussian, Laplacian and bimodal interaction distributions in dimension four. Standard finite size scaling analyses near and at criticality provide estimates of the critical inverse temperatures βcβ_cβc​, critical exponents, and critical values of a number of dimensionless parameters. Independent estimates are obtained for βcβ_cβc​ and the exponent ννν from thermodynamic derivative peak data. A detailed explanation is given of scaling in the thermodynamic limit with the ISG scaling variable τ=1−β2/βc2τ= 1-β^2/β_c^2τ=1−β2/βc2​ and the appropriate scaling expressions. Data over the entire paramagnetic range of temperatures are analysed in order to obtain further estimates of the critical exponents together with correction to scaling terms. The Privman-Fisher ansatz then leads to compact scaling expressions for the whole paramagnetic regime and for all sample sizes L. Comparisons between the 4d ISG models show that the critical dimensionless parameters characteristic of a universality class, and the susceptibility and correlation length critical exponents γγγ and ννν, depend on the form of the interaction distribution. From these observations it can be deduced that critical exponents are not universal in ISGs, at least in dimension four.
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.
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.