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:cond-mat/0207438 · Physical review. B, Condensed matter

Ground-state numerical study of the three-dimensional random-field Ising model

Ilija Dukovski, Jonathan L. Machta

Published Jan 21, 2003• 40 citations
arXivHTMLPDF

Abstract

The random field Ising model in three dimensions with Gaussian random fields is studied at zero temperature for system sizes up to ${60}^{3}.$ For each realization of the normalized random fields, the strength of the random field, $\ensuremath{\Delta}$ and a uniform external, H is adjusted to find the finite-size critical point. The finite-size critical point is identified as the point in the $H\ensuremath{-}\ensuremath{\Delta}$ plane where three degenerate ground states have the largest discontinuities in the magnetization. The discontinuities in the magnetization and bond energy between these ground states are used to calculate the magnetization and specific heat critical exponents and both exponents are found to be near zero.

DOI
10.1103/physrevb.67.014413

Related papers

Numerical study of the three-dimensional random-field Ising model at zero and positive temperature

Yong Wu, Jonathan L. Machta

In this paper the three-dimensional random-field Ising model is studied at both zero temperature and positive temperature. Critical exponents are extracted at zero temperature by finite size scaling analysis of large discontinuities in the bond energy. The heat capacity exponent $\ensuremath{\alpha}$ is found to be near zero. The ground states are determined for a range of external field and disorder strength near the zero temperature critical point and the scaling of ground state tilings of the field-disorder plane is discussed. At positive temperature the specific heat and the susceptibility are obtained using the Wang-Landau algorithm. It is found that sharp peaks are present in these physical quantities for some realizations of systems sized ${16}^{3}$ and larger. These sharp peaks result from flipping large domains and correspond to large discontinuities in ground state bond energies. Finally, zero temperature and positive temperature spin configurations near the critical line are found to be highly correlated suggesting a strong version of the zero temperature fixed point hypothesis.

arXiv:cond-mat/0602314•Physical Review B•Aug 21, 2006•40 citations

Critical behavior of the three-dimensional random-field Ising model: Two-exponent scaling and discontinuous transition

Heiko Rieger

In extensive Monte Carlo simulations the phase transition of the random-field Ising model in three dimensions is investigated. The values of the critical exponents are determined via finite-size scaling. For a Gaussian distribution of the random fields it is found that the correlation length \ensuremath{\xi} diverges with an exponent \ensuremath{\nu}=1.1\ifmmode\pm\else\textpm\fi{}0.2 at the critical temperature and that \ensuremath{\chi}\ensuremath{\sim}${\ensuremath{\xi}}^{2\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}}$ with \ensuremath{\eta}=0.50\ifmmode\pm\else\textpm\fi{}0.05 for the connected susceptibility and ${\mathrm{\ensuremath{\chi}}}_{\mathrm{dis}}$\ensuremath{\sim}${\ensuremath{\xi}}^{4\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}$ with \ensuremath{\eta}\ifmmode\bar\else\textasciimacron\fi{}=1.03\ifmmode\pm\else\textpm\fi{}0.05 for the disconnected susceptibility. Together with the amplitude ratio A=${\mathrm{lim}}_{\mathit{T}\ensuremath{\rightarrow}\mathit{T}\mathit{c}}$${\mathrm{\ensuremath{\chi}}}_{\mathrm{dis}}$/${\mathrm{\ensuremath{\chi}}}^{2}$(${\mathit{h}}_{\mathit{r}}$/T${)}^{2}$ being close to one this gives further support for a two-exponent scaling scenario implying \ensuremath{\eta}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\eta}. The magnetization behaves discontinuously at the transition, i.e., \ensuremath{\beta}=0. However, no divergence for the specific heat and in particular no latent heat is found. Also the probability distribution of the magnetization does not show a multipeak structure that would be characteristic for the phase-coexistence at first-order phase-transition points.

arXiv:cond-mat/9503041•Physical review. B, Condensed matter•Sep 1, 1995•138 citations

Specific-heat exponent of random-field systems via ground-state calculations

Alexander K. Hartmann, A. P. Young

Exact ground states of three-dimensional random field Ising magnets with Gaussian distribution of the disorder are calculated using graph-theoretical algorithms. Systems for different strengths h of the random fields and sizes up to ${N=96}^{3}$ are considered. By numerically differentiating the bond-energy with respect to h a specific-heat-like quantity is obtained, which does not appear to diverge at the critical point but rather exhibits a cusp. We also consider the effect of a small uniform magnetic field, which allows us to calculate the $T=0$ susceptibility. From a finite-size scaling analysis, we obtain the critical exponents $\ensuremath{\nu}=1.32(7),$ $\ensuremath{\alpha}=\ensuremath{-}0.63(7),$ $\ensuremath{\eta}=0.50(3)$ and find that the critical strength of the random field is ${h}_{c}=2.28(1).$ We discuss the significance of the result that $\ensuremath{\alpha}$ appears to be strongly negative.

arXiv:cond-mat/0105310•Physical review. B, Condensed matter•Nov 13, 2001•81 citations

Random-field Ising model: Insight from zero-temperature simulations

Panagiotis E. Theodorakis, Nikolaos G. Fytas

We enlighten some critical aspects of the three-dimensional (d = 3) random-field Ising model from simulations performed at zero temperature.We consider two different, in terms of the field distribution, versions of model, namely a Gaussian random-field Ising model and an equal-weight trimodal random-field Ising model.By implementing a computational approach that maps the ground-state of the system to the maximum-flow optimization problem of a network, we employ the most up-to-date version of the push-relabel algorithm and simulate large ensembles of disorder realizations of both models for a broad range of random-field values and systems sizes V = L × L × L, where L denotes linear lattice size and L max = 156.Using as finite-size measures the sampleto-sample fluctuations of various quantities of physical and technical origin, and the primitive operations of the push-relabel algorithm, we propose, for both types of distributions, estimates of the critical field h c and the critical exponent ν of the correlation length, the latter clearly suggesting that both models share the same universality class.Additional simulations of the Gaussian random-field Ising model at the best-known value of the critical field provide the magnetic exponent ratio β/ν with high accuracy and clear out the controversial issue of the critical exponent α of the specific heat.Finally, we discuss the infinite-limit size extrapolation of energyand order-parameter-based noise to signal ratios related to the self-averaging properties of the model, as well as the critical slowing down aspects of the algorithm.

arXiv:1501.02338•Condensed Matter Physics•Dec 1, 2014•9 citations

Fluctuations and criticality in the random-field Ising model

Panagiotis E. Theodorakis, Ioannis Georgiou, Nikolaos G. Fytas

We investigate the critical properties of the $d=3$ random-field Ising model with a Gaussian field distribution at zero temperature. By implementing suitable graph-theoretical algorithms, we perform a large-scale numerical simulation of the model for a vast range of values of the disorder strength $h$ and system sizes $\mathcal{V}=L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L$, with $L\ensuremath{\le}156$. Using the sample-to-sample fluctuations of various quantities and proper finite-size scaling techniques we estimate with high accuracy the critical disorder strength ${h}_{\mathrm{c}}$ and the correlation length exponent $\ensuremath{\nu}$. Additional simulations in the area of the estimated critical-field strength and relevant scaling analysis of the bond energy suggest bounds for the specific heat critical exponent $\ensuremath{\alpha}$ and the violation of the hyperscaling exponent $\ensuremath{\theta}$. Finally, a data collapse analysis of the order parameter and disconnected susceptibility provides accurate estimates for the critical exponent ratios $\ensuremath{\beta}/\ensuremath{\nu}$ and $\overline{\ensuremath{\gamma}}/\ensuremath{\nu}$, respectively.

Physical Review E•Mar 8, 2013•11 citations

Critical exponents of the three-dimensional random field Ising model

Heiko Rieger, A. P. Young

The phase transition of the three-dimensional random field Ising model with a discrete (+or-h) field distribution is investigated by extensive Monte Carlo simulations. Values of the critical exponents for the correlation length, specific heat, susceptibility, disconnected susceptibility and magnetization are determined simultaneously via finite size scaling. While the magnetization appears to be discontinuous, the specific heat appears to saturate, indicating no latent heat. Sample-to-sample fluctuations of the susceptibility are consistent with the droplet picture for the transition.

arXiv:cond-mat/9303045•Journal of Physics A Mathematical and General•Oct 21, 1993•109 citations
The random field Ising model in three dimensions with Gaussian random fields is studied at zero temperature for system sizes up to 603.{60}^{3}.603. For each realization of the normalized random fields, the strength of the random field, \ensuremathΔ\ensuremath{\Delta}\ensuremathΔ and a uniform external, H is adjusted to find the finite-size critical point. The finite-size critical point is identified as the point in the H\ensuremath−\ensuremathΔH\ensuremath{-}\ensuremath{\Delta}H\ensuremath−\ensuremathΔ plane where three degenerate ground states have the largest discontinuities in the magnetization. The discontinuities in the magnetization and bond energy between these ground states are used to calculate the magnetization and specific heat critical exponents and both exponents are found to be near zero.
In this paper the three-dimensional random-field Ising model is studied at both zero temperature and positive temperature. Critical exponents are extracted at zero temperature by finite size scaling analysis of large discontinuities in the bond energy. The heat capacity exponent $\ensuremath{\alpha}$ is found to be near zero. The ground states are determined for a range of external field and disorder strength near the zero temperature critical point and the scaling of ground state tilings of the field-disorder plane is discussed. At positive temperature the specific heat and the susceptibility are obtained using the Wang-Landau algorithm. It is found that sharp peaks are present in these physical quantities for some realizations of systems sized 163{16}^{3}163 and larger. These sharp peaks result from flipping large domains and correspond to large discontinuities in ground state bond energies. Finally, zero temperature and positive temperature spin configurations near the critical line are found to be highly correlated suggesting a strong version of the zero temperature fixed point hypothesis.
In extensive Monte Carlo simulations the phase transition of the random-field Ising model in three dimensions is investigated. The values of the critical exponents are determined via finite-size scaling. For a Gaussian distribution of the random fields it is found that the correlation length \ensuremath{\xi} diverges with an exponent \ensuremath{\nu}=1.1\ifmmode\pm\else\textpm\fi{}0.2 at the critical temperature and that \ensuremath{\chi}\ensuremath{\sim}${\ensuremath{\xi}}^{2\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}}$ with \ensuremath{\eta}=0.50\ifmmode\pm\else\textpm\fi{}0.05 for the connected susceptibility and ${\mathrm{\ensuremath{\chi}}}_{\mathrm{dis}}$\ensuremath{\sim}${\ensuremath{\xi}}^{4\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}$ with \ensuremath{\eta}\ifmmode\bar\else\textasciimacron\fi{}=1.03\ifmmode\pm\else\textpm\fi{}0.05 for the disconnected susceptibility. Together with the amplitude ratio A=${\mathrm{lim}}_{\mathit{T}\ensuremath{\rightarrow}\mathit{T}\mathit{c}}$${\mathrm{\ensuremath{\chi}}}_{\mathrm{dis}}$/${\mathrm{\ensuremath{\chi}}}^{2}$(hr{\mathit{h}}_{\mathit{r}}hr​/T)2{)}^{2})2 being close to one this gives further support for a two-exponent scaling scenario implying \ensuremath{\eta}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\eta}. The magnetization behaves discontinuously at the transition, i.e., \ensuremath{\beta}=0. However, no divergence for the specific heat and in particular no latent heat is found. Also the probability distribution of the magnetization does not show a multipeak structure that would be characteristic for the phase-coexistence at first-order phase-transition points.
Exact ground states of three-dimensional random field Ising magnets with Gaussian distribution of the disorder are calculated using graph-theoretical algorithms. Systems for different strengths h of the random fields and sizes up to N=963{N=96}^{3}N=963 are considered. By numerically differentiating the bond-energy with respect to h a specific-heat-like quantity is obtained, which does not appear to diverge at the critical point but rather exhibits a cusp. We also consider the effect of a small uniform magnetic field, which allows us to calculate the T=0T=0T=0 susceptibility. From a finite-size scaling analysis, we obtain the critical exponents $\ensuremath{\nu}=1.32(7),$ $\ensuremath{\alpha}=\ensuremath{-}0.63(7),$ $\ensuremath{\eta}=0.50(3)$ and find that the critical strength of the random field is hc=2.28(1).{h}_{c}=2.28(1).hc​=2.28(1). We discuss the significance of the result that $\ensuremath{\alpha}$ appears to be strongly negative.
We enlighten some critical aspects of the three-dimensional (d = 3) random-field Ising model from simulations performed at zero temperature.We consider two different, in terms of the field distribution, versions of model, namely a Gaussian random-field Ising model and an equal-weight trimodal random-field Ising model.By implementing a computational approach that maps the ground-state of the system to the maximum-flow optimization problem of a network, we employ the most up-to-date version of the push-relabel algorithm and simulate large ensembles of disorder realizations of both models for a broad range of random-field values and systems sizes V = L × L × L, where L denotes linear lattice size and L max = 156.Using as finite-size measures the sampleto-sample fluctuations of various quantities of physical and technical origin, and the primitive operations of the push-relabel algorithm, we propose, for both types of distributions, estimates of the critical field h c and the critical exponent ν of the correlation length, the latter clearly suggesting that both models share the same universality class.Additional simulations of the Gaussian random-field Ising model at the best-known value of the critical field provide the magnetic exponent ratio β/ν with high accuracy and clear out the controversial issue of the critical exponent α of the specific heat.Finally, we discuss the infinite-limit size extrapolation of energyand order-parameter-based noise to signal ratios related to the self-averaging properties of the model, as well as the critical slowing down aspects of the algorithm.
We investigate the critical properties of the d=3d=3d=3 random-field Ising model with a Gaussian field distribution at zero temperature. By implementing suitable graph-theoretical algorithms, we perform a large-scale numerical simulation of the model for a vast range of values of the disorder strength hhh and system sizes $\mathcal{V}=L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L$, with $L\ensuremath{\le}156$. Using the sample-to-sample fluctuations of various quantities and proper finite-size scaling techniques we estimate with high accuracy the critical disorder strength hc{h}_{\mathrm{c}}hc​ and the correlation length exponent $\ensuremath{\nu}$. Additional simulations in the area of the estimated critical-field strength and relevant scaling analysis of the bond energy suggest bounds for the specific heat critical exponent $\ensuremath{\alpha}$ and the violation of the hyperscaling exponent $\ensuremath{\theta}$. Finally, a data collapse analysis of the order parameter and disconnected susceptibility provides accurate estimates for the critical exponent ratios $\ensuremath{\beta}/\ensuremath{\nu}$ and $\overline{\ensuremath{\gamma}}/\ensuremath{\nu}$, respectively.
The phase transition of the three-dimensional random field Ising model with a discrete (+or-h) field distribution is investigated by extensive Monte Carlo simulations. Values of the critical exponents for the correlation length, specific heat, susceptibility, disconnected susceptibility and magnetization are determined simultaneously via finite size scaling. While the magnetization appears to be discontinuous, the specific heat appears to saturate, indicating no latent heat. Sample-to-sample fluctuations of the susceptibility are consistent with the droplet picture for the transition.