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/9503041 · Physical review. B, Condensed matter

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

Heiko Rieger

Published Sep 1, 1995• 138 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevb.52.6659

Related papers

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

Critical exponents of four-dimensional random-field Ising systems

Alexander K. Hartmann

The ferromagnet-to-paramagnet transition of the four-dimensional random-field Ising model with Gaussian distribution of the random fields is studied. Exact ground states of systems with sizes up to ${32}^{4}$ are obtained using graph theoretical algorithms. The magnetization, the disconnected susceptibility, the susceptibility, and a specific-heat-like quantity are calculated. Using finite-size scaling techniques, the corresponding critical exponents are obtained: $\ensuremath{\beta}=0.13(5),$ $\overline{\ensuremath{\gamma}}=2.83(50),$ $\ensuremath{\gamma}=1.42(20).$ For the specific heat both $\ensuremath{\alpha}=0$ (logarithmic divergence) and an algebraic divergence with $\ensuremath{\alpha}=0.26(5)$ are compatible with the data. Furthermore, values for the critical randomness ${h}_{c}=4.18(1)$ and the correlation-length exponent $\ensuremath{\nu}=0.78(10)$ were found. These independently obtained exponents are compatible with all (hyper)scaling relations and support the two-exponent scenario $(\overline{\ensuremath{\gamma}}=2\ensuremath{\gamma}).$

arXiv:cond-mat/0201441•Physical review. B, Condensed matter•May 7, 2002•30 citations

Monte Carlo study of the random-field Ising model

Mark E. J. Newman, G. T. Barkema

Using a cluster-flipping Monte Carlo algorithm combined with a generalization of the histogram reweighting scheme of Ferrenberg and Swendsen [Phys. Rev. Lett. 61, 2635 (1988); 63, 1195 (1989)], we have studied the equilibrium properties of the thermal random-field Ising model on a cubic lattice in three dimensions. We have equilibrated systems of L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L spins, with values of L up to 32, and for these systems the cluster-flipping method appears to a large extent to overcome the slow equilibration seen in single-spin-flip methods. From the results of our simulations we have extracted values for the critical exponents and the critical temperature and randomness of the model by finite size scaling. For the exponents we find \ensuremath{\nu}=1.02\ifmmode\pm\else\textpm\fi{}0.06, \ensuremath{\beta}=0.06\ifmmode\pm\else\textpm\fi{}0.07, \ensuremath{\gamma}=1.9\ifmmode\pm\else\textpm\fi{}0.2, and \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2.9\ifmmode\pm\else\textpm\fi{}0.2, where \ensuremath{\nu}, \ensuremath{\beta}, \ensuremath{\gamma}, and \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{} govern the critical singularities in the correlation length, magnetization, and connected and disconnected susceptibilities, respectively. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9507044•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jan 1, 1996•81 citations

Critical behavior of the random-field Ising model

Misha Gofman, Joan Adler, Amnon Aharony A +2 more

We study the critical properties of the random field Ising model in general dimension d using high-temperature expansions for the susceptibility, \ensuremath{\chi}=${\ensuremath{\sum}}_{\mathit{j}}$ [〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$${\mathrm{〉}}_{\mathit{T}}$-〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{〉}}_{\mathit{T}}$〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$${\mathrm{〉}}_{\mathit{T}}$${]}_{\mathit{h}}$ and the structure factor, G=${\ensuremath{\sum}}_{\mathit{j}}$[〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$${\mathrm{〉}}_{\mathit{T}}$${]}_{\mathit{h}}$, where 〈${\mathrm{〉}}_{\mathit{T}}$ indicates a canonical average at temperature T for an arbitrary configuration of random fields and [${]}_{\mathit{h}}$ indicates an average over random fields. We treated two distributions of random fields, the bimodal in which each ${\mathit{h}}_{\mathit{i}}$=\ifmmode\pm\else\textpm\fi{}${\mathit{h}}_{0}$ and a Gaussian distribution in which each ${\mathit{h}}_{\mathit{i}}$ has variance ${\mathit{h}}_{0}^{2}$. We obtained series for \ensuremath{\chi} and G in the form ${\ensuremath{\sum}}_{\mathit{n}=1,15}$${\mathit{a}}_{\mathit{n}}$(g,d)(J/T${)}^{\mathit{n}}$, where J is the exchange constant and the coefficients ${\mathit{a}}_{\mathit{n}}$(g,d) are polynomials in g\ensuremath{\equiv}${\mathit{h}}_{0}^{2}$/${\mathit{J}}^{2}$ and in d. We assume that as T approaches its critical value, ${\mathit{T}}_{\mathit{c}}$, one has \ensuremath{\chi}\ensuremath{\sim}(T-${\mathit{T}}_{\mathit{c}}$${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}}$ and G\ensuremath{\sim}(T-${\mathit{T}}_{\mathit{c}}$${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}$. For dimensions above d=2 we find a range of values of g for which the critical exponents obtained from our series seem not to depend on g.For large values of g our results show a g dependence which is attributable to either a tricritical point or a first-order transition. All our results for critical exponents suggest that \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}, in agreement with the two-exponent scaling picture. In addition we have also constructed series for the amplitude ratio, A=(G/${\mathrm{\ensuremath{\chi}}}^{2}$)(${\mathit{T}}^{2}$)/(${\mathit{gJ}}^{2}$). We find that A approaches a constant value as T\ensuremath{\rightarrow}${\mathit{T}}_{\mathit{c}}$ (consistent with \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}) with A\ensuremath{\approxeq}1. It appears that A is somewhat larger for the bimodal than for the Gaussian model, in agreement with a recent analysis at high d. \textcopyright{} 1996 The American Physical Society.

Physical review. B, Condensed matter•Mar 1, 1996•59 citations

Dynamical Monte Carlo studies of the three-dimensional bimodal random-field Ising model

Wanjie Xiong, Yafei Dai

The nature of the three-dimensional random-field Ising model with a bimodal probability distribution is investigated using finite-time scaling combined with the Monte Carlo renormalization group method, in the presence of a linearly varying temperature. Our results support the existence of a first-order phase transition for this model, obtained through reducing the influence of finite-size effects. The critical exponents are estimated to be ν = 0.74(3), β = 0.27(1), α = − 0.023(5), and γ = 1.48(3) with corrections to the scaling. The Rushbrooke scaling law is satisfied with these exponents, and in turn identifies the asymptotic critical behavior.

Journal of Statistical Mechanics Theory and Experiment•May 16, 2012•6 citations

On the phase transition of the 3D random field Ising model

Marco Picco, Nicolas Sourlas

We present new numerical simulations of the random field Ising model in three dimensions at zero temperature. The critical exponents are found to agree with previous results. We study the magnetic susceptibility by applying a small magnetic field perturbation. We find the critical amplitude ratio of the magnetic susceptibilities to be very large, equal to 233.1 ± 1.5. We find strong sample to sample fluctuations which obey finite size scaling. The probability distribution of the size of small energy excitations is maximally non-self-averaging, obeying a double peak distribution, and is finite size scaling invariant. We also study the approach to the thermodynamic limit of the ground state magnetization at the phase transition.

arXiv:1310.2364•Journal of Statistical Mechanics Theory and Experiment•Mar 20, 2014•11 citations
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ξ2\ensuremath−\ensuremathη{\ensuremath{\xi}}^{2\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}}\ensuremathξ2\ensuremath−\ensuremathη with \ensuremath{\eta}=0.50\ifmmode\pm\else\textpm\fi{}0.05 for the connected susceptibility and \ensuremathχdis{\mathrm{\ensuremath{\chi}}}_{\mathrm{dis}}\ensuremathχdis​\ensuremath{\sim}\ensuremathξ4\ensuremath−\ensuremathη\ifmmode\elseˉ\textasciimacron\fi{\ensuremath{\xi}}^{4\mathrm{\ensuremath{-}}\mathrm{\ensuremath{\eta}}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}\ensuremathξ4\ensuremath−\ensuremathη\ifmmode\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}}/\ensuremathχ2{\mathrm{\ensuremath{\chi}}}^{2}\ensuremathχ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.
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.
The ferromagnet-to-paramagnet transition of the four-dimensional random-field Ising model with Gaussian distribution of the random fields is studied. Exact ground states of systems with sizes up to 324{32}^{4}324 are obtained using graph theoretical algorithms. The magnetization, the disconnected susceptibility, the susceptibility, and a specific-heat-like quantity are calculated. Using finite-size scaling techniques, the corresponding critical exponents are obtained: $\ensuremath{\beta}=0.13(5),$ $\overline{\ensuremath{\gamma}}=2.83(50),$ $\ensuremath{\gamma}=1.42(20).$ For the specific heat both $\ensuremath{\alpha}=0$ (logarithmic divergence) and an algebraic divergence with $\ensuremath{\alpha}=0.26(5)$ are compatible with the data. Furthermore, values for the critical randomness hc=4.18(1){h}_{c}=4.18(1)hc​=4.18(1) and the correlation-length exponent $\ensuremath{\nu}=0.78(10)$ were found. These independently obtained exponents are compatible with all (hyper)scaling relations and support the two-exponent scenario $(\overline{\ensuremath{\gamma}}=2\ensuremath{\gamma}).$
Using a cluster-flipping Monte Carlo algorithm combined with a generalization of the histogram reweighting scheme of Ferrenberg and Swendsen [Phys. Rev. Lett. 61, 2635 (1988); 63, 1195 (1989)], we have studied the equilibrium properties of the thermal random-field Ising model on a cubic lattice in three dimensions. We have equilibrated systems of L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L spins, with values of L up to 32, and for these systems the cluster-flipping method appears to a large extent to overcome the slow equilibration seen in single-spin-flip methods. From the results of our simulations we have extracted values for the critical exponents and the critical temperature and randomness of the model by finite size scaling. For the exponents we find \ensuremath{\nu}=1.02\ifmmode\pm\else\textpm\fi{}0.06, \ensuremath{\beta}=0.06\ifmmode\pm\else\textpm\fi{}0.07, \ensuremath{\gamma}=1.9\ifmmode\pm\else\textpm\fi{}0.2, and \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2.9\ifmmode\pm\else\textpm\fi{}0.2, where \ensuremath{\nu}, \ensuremath{\beta}, \ensuremath{\gamma}, and \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{} govern the critical singularities in the correlation length, magnetization, and connected and disconnected susceptibilities, respectively. \textcopyright{} 1996 The American Physical Society.
We study the critical properties of the random field Ising model in general dimension d using high-temperature expansions for the susceptibility, \ensuremath{\chi}=${\ensuremath{\sum}}_{\mathit{j}}$ [〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$〉T{\mathrm{〉}}_{\mathit{T}}〉T​-〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$〉T{\mathrm{〉}}_{\mathit{T}}〉T​〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$〉T{\mathrm{〉}}_{\mathit{T}}〉T​]h{]}_{\mathit{h}}]h​ and the structure factor, G=${\ensuremath{\sum}}_{\mathit{j}}$[〈${\mathrm{\ensuremath{\sigma}}}_{\mathit{i}}$${\mathrm{\ensuremath{\sigma}}}_{\mathit{j}}$〉T{\mathrm{〉}}_{\mathit{T}}〉T​]h{]}_{\mathit{h}}]h​, where 〈〉T{\mathrm{〉}}_{\mathit{T}}〉T​ indicates a canonical average at temperature T for an arbitrary configuration of random fields and []h{]}_{\mathit{h}}]h​ indicates an average over random fields. We treated two distributions of random fields, the bimodal in which each hi{\mathit{h}}_{\mathit{i}}hi​=\ifmmode\pm\else\textpm\fi{}h0{\mathit{h}}_{0}h0​ and a Gaussian distribution in which each hi{\mathit{h}}_{\mathit{i}}hi​ has variance h02{\mathit{h}}_{0}^{2}h02​. We obtained series for \ensuremath{\chi} and G in the form ${\ensuremath{\sum}}_{\mathit{n}=1,15}$an{\mathit{a}}_{\mathit{n}}an​(g,d)(J/T)n{)}^{\mathit{n}})n, where J is the exchange constant and the coefficients an{\mathit{a}}_{\mathit{n}}an​(g,d) are polynomials in g\ensuremath{\equiv}h02{\mathit{h}}_{0}^{2}h02​/J2{\mathit{J}}^{2}J2 and in d. We assume that as T approaches its critical value, Tc{\mathit{T}}_{\mathit{c}}Tc​, one has \ensuremath{\chi}\ensuremath{\sim}(T-Tc{\mathit{T}}_{\mathit{c}}Tc​${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}}$ and G\ensuremath{\sim}(T-Tc{\mathit{T}}_{\mathit{c}}Tc​${)}^{\mathrm{\ensuremath{-}}\ensuremath{\gamma}\mathrm{\ifmmode\bar\else\textasciimacron\fi{}}}$. For dimensions above d=2 we find a range of values of g for which the critical exponents obtained from our series seem not to depend on g.For large values of g our results show a g dependence which is attributable to either a tricritical point or a first-order transition. All our results for critical exponents suggest that \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}, in agreement with the two-exponent scaling picture. In addition we have also constructed series for the amplitude ratio, A=(G/${\mathrm{\ensuremath{\chi}}}^{2}$)(T2{\mathit{T}}^{2}T2)/(gJ2{\mathit{gJ}}^{2}gJ2). We find that A approaches a constant value as T\ensuremath{\rightarrow}Tc{\mathit{T}}_{\mathit{c}}Tc​ (consistent with \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=2\ensuremath{\gamma}) with A\ensuremath{\approxeq}1. It appears that A is somewhat larger for the bimodal than for the Gaussian model, in agreement with a recent analysis at high d. \textcopyright{} 1996 The American Physical Society.
The nature of the three-dimensional random-field Ising model with a bimodal probability distribution is investigated using finite-time scaling combined with the Monte Carlo renormalization group method, in the presence of a linearly varying temperature. Our results support the existence of a first-order phase transition for this model, obtained through reducing the influence of finite-size effects. The critical exponents are estimated to be ν = 0.74(3), β = 0.27(1), α = − 0.023(5), and γ = 1.48(3) with corrections to the scaling. The Rushbrooke scaling law is satisfied with these exponents, and in turn identifies the asymptotic critical behavior.
We present new numerical simulations of the random field Ising model in three dimensions at zero temperature. The critical exponents are found to agree with previous results. We study the magnetic susceptibility by applying a small magnetic field perturbation. We find the critical amplitude ratio of the magnetic susceptibilities to be very large, equal to 233.1 ± 1.5. We find strong sample to sample fluctuations which obey finite size scaling. The probability distribution of the size of small energy excitations is maximally non-self-averaging, obeying a double peak distribution, and is finite size scaling invariant. We also study the approach to the thermodynamic limit of the ground state magnetization at the phase transition.