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/0501619 · Physical Review Letters

Ground States and Thermal States of the Random Field Ising Model

Yong Wu, Jonathan L. Machta

Published Sep 23, 2005• 48 citations
arXivHTMLPDF

Abstract

The random field Ising model is studied numerically at both zero and positive temperature. Ground states are mapped out in a region of random field and external field strength. Thermal states and thermodynamic properties are obtained for all temperatures using the Wang-Landau algorithm. The specific heat and susceptibility typically display sharp peaks in the critical region for large systems and strong disorder. These sharp peaks result from large domains flipping. For a given realization of disorder, ground states and thermal states near the critical line are found to be strongly correlated--a concrete manifestation of the zero temperature fixed point scenario.

DOI
10.1103/physrevlett.95.137208

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

Computational studies of classical disordered spin systems

Yong Wu

In this thesis classical disordered spin systems, in particular, the random field Ising model (RFIM), are studied using numerical methods. The previous experimental and theoretical research on the RFIM is reviewed. Critical exponents of the RFIM are extracted at zero temperature. The heat capacity exponent α is found to be near zero. The ground states within a certain region with varying external field and disorder are portrayed, and it is speculated that if related quantities are scaled properly, ground state portraits at different system sizes are statistically similar. The performance of the Wang-Landau (WL) algorithm, which is used in our positive temperature study of the RFIM, is measured by the mean first-passage time. It is proved that the mean first, passage time of the mean field Ising model scales as N 2 log N, where N is the number of spins. The logarithmic correction term to the unbiased Markovian random walk is due to the boundaries. When applied to the Ising model with nearest neighbor interactions, the WL algorithm suffers from power law slowing down. Carefully designed cluster updating, however, can accelerate the WL algorithm and eliminate the slowing down. The WL algorithm is used in our finite temperature study of the RFIM. Physical quantities, such as the specific heat and the susceptibility, are obtained by the WL algorithm over a broad range of temperature. It is discovered that for some realizations with system size larger than 8 3, sharp peaks are present in these physical quantities. These sharp peaks result from flipping a large domain, and furthermore, the critical domains are strongly correlated to the domains found at the zero temperature transition. This observation suggests a strong version of the zero temperature fixed point picture. The MKRG approximation is employed to test the correlation between the critical states and the ground states in much larger systems. It is found that the fraction of correlated realizations drops as system size grows, while there are still significant number of correlated realizations even for systems as large as 10243. Whether substantial correlations will remain in the thermodynamic limit is unclear from our present data.

•Jan 1, 2006•0 citations

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

Ilija Dukovski, Jonathan L. Machta

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.

arXiv:cond-mat/0207438•Physical review. B, Condensed matter•Jan 21, 2003•40 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

Thermal critical behavior and universality aspects of the three-dimensional random-field Ising model

Anastasios Malakis, Nikolaos G. Fytas

The three-dimensional bimodal random-field Ising model is investigated using the N-fold version of the Wang-Landau algorithm. The essential energy subspaces are determined by the recently developed critical minimum energy subspace technique, and two implementations of this scheme are utilized. The random fields are obtained from a bimodal discrete (±Δ) distribution, and we study the model for various values of the disorder strength Δ, Δ=0.5,1,1.5 and 2, on cubic lattices with linear sizes L=4–24. We extract information for the probability distributions of the specific heat peaks over samples of random fields. This permits us to obtain the phase diagram and present the finite-size behavior of the specific heat. The question of saturation of the specific heat is re-examined and it is shown that the open problem of universality for the random-field Ising model is strongly influenced by the lack of self-averaging of the model. This property appears to be substantially depended on the disorder strength.

•Jan 1, 2006•5 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
The random field Ising model is studied numerically at both zero and positive temperature. Ground states are mapped out in a region of random field and external field strength. Thermal states and thermodynamic properties are obtained for all temperatures using the Wang-Landau algorithm. The specific heat and susceptibility typically display sharp peaks in the critical region for large systems and strong disorder. These sharp peaks result from large domains flipping. For a given realization of disorder, ground states and thermal states near the critical line are found to be strongly correlated--a concrete manifestation of the zero temperature fixed point scenario.
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 this thesis classical disordered spin systems, in particular, the random field Ising model (RFIM), are studied using numerical methods. The previous experimental and theoretical research on the RFIM is reviewed. Critical exponents of the RFIM are extracted at zero temperature. The heat capacity exponent α is found to be near zero. The ground states within a certain region with varying external field and disorder are portrayed, and it is speculated that if related quantities are scaled properly, ground state portraits at different system sizes are statistically similar. The performance of the Wang-Landau (WL) algorithm, which is used in our positive temperature study of the RFIM, is measured by the mean first-passage time. It is proved that the mean first, passage time of the mean field Ising model scales as N 2 log N, where N is the number of spins. The logarithmic correction term to the unbiased Markovian random walk is due to the boundaries. When applied to the Ising model with nearest neighbor interactions, the WL algorithm suffers from power law slowing down. Carefully designed cluster updating, however, can accelerate the WL algorithm and eliminate the slowing down. The WL algorithm is used in our finite temperature study of the RFIM. Physical quantities, such as the specific heat and the susceptibility, are obtained by the WL algorithm over a broad range of temperature. It is discovered that for some realizations with system size larger than 8 3, sharp peaks are present in these physical quantities. These sharp peaks result from flipping a large domain, and furthermore, the critical domains are strongly correlated to the domains found at the zero temperature transition. This observation suggests a strong version of the zero temperature fixed point picture. The MKRG approximation is employed to test the correlation between the critical states and the ground states in much larger systems. It is found that the fraction of correlated realizations drops as system size grows, while there are still significant number of correlated realizations even for systems as large as 10243. Whether substantial correlations will remain in the thermodynamic limit is unclear from our present data.
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{\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.
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 three-dimensional bimodal random-field Ising model is investigated using the N-fold version of the Wang-Landau algorithm. The essential energy subspaces are determined by the recently developed critical minimum energy subspace technique, and two implementations of this scheme are utilized. The random fields are obtained from a bimodal discrete (±Δ) distribution, and we study the model for various values of the disorder strength Δ, Δ=0.5,1,1.5 and 2, on cubic lattices with linear sizes L=4–24. We extract information for the probability distributions of the specific heat peaks over samples of random fields. This permits us to obtain the phase diagram and present the finite-size behavior of the specific heat. The question of saturation of the specific heat is re-examined and it is shown that the open problem of universality for the random-field Ising model is strongly influenced by the lack of self-averaging of the model. This property appears to be substantially depended on the disorder strength.
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.