© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:cond-mat/9507044 · Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics

Monte Carlo study of the random-field Ising model

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

Published Jan 1, 1996• 81 citations

Abstract

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.

DOI: 10.1103/physreve.53.393

Related papers

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

Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model

Alan M. Ferrenberg, Jiahao Xu, David Paul Landau

While the three-dimensional Ising model has defied analytic solution, various numerical methods like Monte Carlo, Monte Carlo renormalization group, and series expansion have provided precise information about the phase transition. Using Monte Carlo simulation that employs the Wolff cluster flipping algorithm with both 32-bit and 53-bit random number generators and data analysis with histogram reweighting and quadruple precision arithmetic, we have investigated the critical behavior of the simple cubic Ising Model, with lattice sizes ranging from 16^{3} to 1024^{3}. By analyzing data with cross correlations between various thermodynamic quantities obtained from the same data pool, e.g., logarithmic derivatives of magnetization and derivatives of magnetization cumulants, we have obtained the critical inverse temperature K_{c}=0.221654626(5) and the critical exponent of the correlation length ν=0.629912(86) with precision that exceeds all previous Monte Carlo estimates.

arXiv:1806.03558•Physical review. E•Apr 4, 2018•191 citations

92 Years of the Ising Model: A High Resolution Monte Carlo Study

Jiahao Xu, Alan M. Ferrenberg, David Paul Landau

Using extensive Monte Carlo simulations that employ the Wolff cluster flipping and data analysis with histogram reweighting and quadruple precision arithmetic, we have investigated the critical behavior of the simple cubic Ising model with lattice sizes ranging from 16 3 to 1024 3 . By analyzing data with cross correlations between various thermodynamic quantities obtained from the same data pool, we obtained the critical inverse temperature K c = 0.221 654 626(5) and the critical exponent of the correlation length ν = 0.629 912(86) with precision that improves upon previous Monte Carlo estimates.

Journal of Physics Conference Series•Apr 1, 2018•16 citations

Static critical behavior of three-dimensional classical Heisenberg models: A high-resolution Monte Carlo study

Kun Chen, Alan M. Ferrenberg, David Paul Landau

Using both recently developed cluster-algorithm and histogram methods, we have carried out a high-resolution Monte Carlo study of static critical properties of classical ferromagnetic Heisenberg models. Extensive Monte Carlo simulations were performed at several temperatures in the critical region, using an improved cluster-updating scheme, on L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L simple-cubic and body-centered-cubic systems with L\ensuremath{\le}40. Thermodynamic quantities as a function of temperature in the vicinity of the critical point were obtained by an optimized multiple-histogram method, and the critical temperature and static critical exponents were extracted using finite-size scaling. Our best estimates for the inverse critical temperatures are 0.693 035(37) for the simple-cubic system and 0.486 798(12) for the body-centered-cubic system. Estimated static critical exponents for both systems agree with each other within their respective error bars, and the mean estimates \ensuremath{\nu}=0.7048(30) and \ensuremath{\gamma}=1.3873(85) are also in excellent agreement with field-theoretic predictions 0.705(3) and 1.386(4).

Physical review. B, Condensed matter•Aug 1, 1993•331 citations

Ising model on three-dimensional random lattices: A Monte Carlo study

Wolfhard Janke E, Ramon Villanova

We report single-cluster Monte Carlo simulations of the Ising model on three-dimensional Poissonian random lattices with up to $128000\ensuremath{\approx}{50}^{3}$ sites which are linked together according to the Voronoi-Delaunay prescription. For each lattice size quenched averages are performed over 96 realizations. By using reweighting techniques and finite-size scaling analyses we investigate the critical properties of the model in the close vicinity of the phase transition point. Our random lattice data provide strong evidence that, for the available system sizes, the resulting effective critical exponents are indistinguishable from recent high-precision estimates obtained in Monte Carlo studies of the Ising model and ${\ensuremath{\varphi}}^{4}$ field theory on three-dimensional regular cubic lattices.

arXiv:cond-mat/0211446•Physical review. B, Condensed matter•Oct 31, 2002•31 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

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.

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.

While the three-dimensional Ising model has defied analytic solution, various numerical methods like Monte Carlo, Monte Carlo renormalization group, and series expansion have provided precise information about the phase transition. Using Monte Carlo simulation that employs the Wolff cluster flipping algorithm with both 32-bit and 53-bit random number generators and data analysis with histogram reweighting and quadruple precision arithmetic, we have investigated the critical behavior of the simple cubic Ising Model, with lattice sizes ranging from 16^{3} to 1024^{3}. By analyzing data with cross correlations between various thermodynamic quantities obtained from the same data pool, e.g., logarithmic derivatives of magnetization and derivatives of magnetization cumulants, we have obtained the critical inverse temperature K_{c}=0.221654626(5) and the critical exponent of the correlation length ν=0.629912(86) with precision that exceeds all previous Monte Carlo estimates.

Using extensive Monte Carlo simulations that employ the Wolff cluster flipping and data analysis with histogram reweighting and quadruple precision arithmetic, we have investigated the critical behavior of the simple cubic Ising model with lattice sizes ranging from 16 3 to 1024 3 . By analyzing data with cross correlations between various thermodynamic quantities obtained from the same data pool, we obtained the critical inverse temperature K c = 0.221 654 626(5) and the critical exponent of the correlation length ν = 0.629 912(86) with precision that improves upon previous Monte Carlo estimates.

Using both recently developed cluster-algorithm and histogram methods, we have carried out a high-resolution Monte Carlo study of static critical properties of classical ferromagnetic Heisenberg models. Extensive Monte Carlo simulations were performed at several temperatures in the critical region, using an improved cluster-updating scheme, on L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L simple-cubic and body-centered-cubic systems with L\ensuremath{\le}40. Thermodynamic quantities as a function of temperature in the vicinity of the critical point were obtained by an optimized multiple-histogram method, and the critical temperature and static critical exponents were extracted using finite-size scaling. Our best estimates for the inverse critical temperatures are 0.693 035(37) for the simple-cubic system and 0.486 798(12) for the body-centered-cubic system. Estimated static critical exponents for both systems agree with each other within their respective error bars, and the mean estimates \ensuremath{\nu}=0.7048(30) and \ensuremath{\gamma}=1.3873(85) are also in excellent agreement with field-theoretic predictions 0.705(3) and 1.386(4).

We report single-cluster Monte Carlo simulations of the Ising model on three-dimensional Poissonian random lattices with up to $128000\ensuremath{\approx}{50}^{3}$ sites which are linked together according to the Voronoi-Delaunay prescription. For each lattice size quenched averages are performed over 96 realizations. By using reweighting techniques and finite-size scaling analyses we investigate the critical properties of the model in the close vicinity of the phase transition point. Our random lattice data provide strong evidence that, for the available system sizes, the resulting effective critical exponents are indistinguishable from recent high-precision estimates obtained in Monte Carlo studies of the Ising model and ${\ensuremath{\varphi}}^{4}$ field theory on three-dimensional regular cubic lattices.

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.