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

Physical Review Letters

Monte Carlo calculation of free energy, critical point, and surface critical behavior of three-dimensional Heisenberg ferromagnets

M. P. Nightingale, Henk W. J. Blöte

Published Apr 11, 1988• 81 citations
Publisher

Abstract

A transfer-matrix Monte Carlo technique is developed to compute the free energy of three-dimensional, classical Heisenberg ferromagnets. From the free energy of systems with periodic and antiperiodic boundary conditions, helicity moduli are calculated. From these the critical couplings for simple-cubic (sc) and face-centered-cubic lattices are estimated by use of finite-size scaling. For the simple-cubic lattice, the critical dimension of the surface magnetization is estimated with standard Monte Carlo methods, yielding a result in excellent agreement with the \ensuremath{\epsilon}-expansion work of Diehl and N\"usser.

DOI
10.1103/physrevlett.60.1562

Related papers

Critical behavior of the helicity modulus for the classical Heisenberg model

Robert Goodell Brown, Mikael Ciftan

The critical scaling of the helicity modulus of the classical O(3) $3d$ Heisenberg ferromagnet is studied directly. Monte Carlo methods that impose either an antiperiodic boundary condition or a finite twist of definite handedness across otherwise periodic boundaries in one lattice direction are used to measure scale-dependent enthalpy variations in a simple cubic lattice at the ferromagnetic critical temperature. Finite-size scaling is then used to determine the critical exponents ${v}_{E}$ and ${v}_{F}$ for helicity and, by evaluating three independent hyperscaling-linked pairs of $\ensuremath{\nu}$ and $\ensuremath{\alpha}$, to test hyperscaling for this model. It is observed that antiperiodic boundary conditions in particular constrain the lattice to have a nonzero topological charge, establishing a connection between topological charge and helicity in the model.

Physical Review B•Dec 14, 2006•12 citations

High-accuracy Monte Carlo study of the three-dimensional classical Heisenberg ferromagnet

Pawel Peczak, Alan M. Ferrenberg, D. P. Landau

Using extensive Monte Carlo simulations, we study the equilibrium properties of the simple-cubic, classical Heisenberg ferromagnet. We employ very long runs for L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L lattices to obtain high-precision data for the magnetization probability distribution. Using finite-size scaling for L\ensuremath{\le}24 and an optimized multiple-histogram data analysis, we obtain an accurate value of the inverse critical temperature J/${\mathit{k}}_{\mathit{B}}$${\mathit{T}}_{\mathit{c}}$=0.6929\ifmmode\pm\else\textpm\fi{}0.0001, which is higher than previously accepted estimates. Calculated values of various static exponents are in excellent agreement with renormalization-group and \ensuremath{\epsilon}-expansion predictions.

Physical review. B, Condensed matter•Mar 1, 1991•286 citations

Transfer-matrix Monte Carlo estimates of critical points in the simple-cubic Ising, planar, and Heisenberg models

M. P. Nightingale, Henk W. J. Blöte

The principle and the efficiency of the Monte Carlo transfer-matrix algorithm are discussed. Enhancements of this algorithm are illustrated by applications to several phase transitions in lattice spin models. We demonstrate how the statistical noise can be reduced considerably by a similarity transformation of the transfer matrix using a variational estimate of its leading eigenvector, in analogy with a common practice in various quantum Monte Carlo techniques. Here we take the two-dimensional coupled XY-Ising model as an example. Furthermore, we calculate interface free energies of finite three-dimensional O(n) models, for the three cases n=1, 2, and 3. Application of finite-size scaling to the numerical results yields estimates of the critical points of these three models. The statistical precision of the estimates is satisfactory for the modest amount of computer time spent. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9602089•Physical review. B, Condensed matter•Jul 1, 1996•21 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

Monte Carlo Calculation of the Scaling Equation of State for the Classical Heisenberg Ferromagnet

Kurt Binder, Heiner Muller-Krumbhaar

Using a Monte Carlo procedure with a self-consistent-field boundary condition, the magnetization of a simple cubic classical Heisenberg ferromagnet with nearest-neighbor interactions only is calculated. For reduced temperatures and reduced fields in the range $0.02<~|\frac{T}{{T}_{c}}\ensuremath{-}1|$, $\frac{H}{J}<~0.5$ the data of this computer experiment can be described in terms of the "effective" critical exponents $0.31<~\ensuremath{\beta}<~0.35$, $1.33<~\ensuremath{\gamma}<~1.39$, and $4.9<~\ensuremath{\delta}<~5.3$. Although the present data are not very close to the critical point they do obey the homogeneity requirement and determine the scaling function rather precisely. This function agrees very well with the scaling function for the face-centered-cubic classical Heisenberg magnet derived recently by Milo\ifmmode \check{s}\else \v{s}\fi{}evi\ifmmode \acute{c}\else \'{c}\fi{} and Stanley using high-temperature---series-expansion techniques. This agreement supports their hypothesis that neither critical exponents nor the scaling function depend on the lattice structure in the Heisenberg model.

Physical review. B, Solid state•Apr 1, 1973•67 citations

Monte Carlo study of critical behavior in the three-dimensional classical Heisenberg ferromagnet

Kun Chen, Alan M. Ferrenberg, David Paul Landau

The static critical properties of classical, isotropic Heisenberg ferromagnets have been studied numerically. Using an improved cluster algorithm we performed extensive Monte Carlo simulations on L×L×L simple cubic and body-centered cubic systems with L≤40 over a temperature range covering the critical region. Temperature dependences of thermodynamic quantities were determined with an optimized multiple-histogram method, and both critical temperatures and critical exponents were extracted from finite-size scaling of several different thermodynamic quantities. It was found that within their respective error bars the critical exponents for both ferromagnets agree with each other and with field theoretical predictions.

Journal of Applied Physics•May 15, 1993•22 citations
A transfer-matrix Monte Carlo technique is developed to compute the free energy of three-dimensional, classical Heisenberg ferromagnets. From the free energy of systems with periodic and antiperiodic boundary conditions, helicity moduli are calculated. From these the critical couplings for simple-cubic (sc) and face-centered-cubic lattices are estimated by use of finite-size scaling. For the simple-cubic lattice, the critical dimension of the surface magnetization is estimated with standard Monte Carlo methods, yielding a result in excellent agreement with the \ensuremath{\epsilon}-expansion work of Diehl and N"usser.
The critical scaling of the helicity modulus of the classical O(3) 3d3d3d Heisenberg ferromagnet is studied directly. Monte Carlo methods that impose either an antiperiodic boundary condition or a finite twist of definite handedness across otherwise periodic boundaries in one lattice direction are used to measure scale-dependent enthalpy variations in a simple cubic lattice at the ferromagnetic critical temperature. Finite-size scaling is then used to determine the critical exponents vE{v}_{E}vE​ and vF{v}_{F}vF​ for helicity and, by evaluating three independent hyperscaling-linked pairs of $\ensuremath{\nu}$ and $\ensuremath{\alpha}$, to test hyperscaling for this model. It is observed that antiperiodic boundary conditions in particular constrain the lattice to have a nonzero topological charge, establishing a connection between topological charge and helicity in the model.
Using extensive Monte Carlo simulations, we study the equilibrium properties of the simple-cubic, classical Heisenberg ferromagnet. We employ very long runs for L\ifmmode\times\else\texttimes\fi{}L\ifmmode\times\else\texttimes\fi{}L lattices to obtain high-precision data for the magnetization probability distribution. Using finite-size scaling for L\ensuremath{\le}24 and an optimized multiple-histogram data analysis, we obtain an accurate value of the inverse critical temperature J/kB{\mathit{k}}_{\mathit{B}}kB​Tc{\mathit{T}}_{\mathit{c}}Tc​=0.6929\ifmmode\pm\else\textpm\fi{}0.0001, which is higher than previously accepted estimates. Calculated values of various static exponents are in excellent agreement with renormalization-group and \ensuremath{\epsilon}-expansion predictions.
The principle and the efficiency of the Monte Carlo transfer-matrix algorithm are discussed. Enhancements of this algorithm are illustrated by applications to several phase transitions in lattice spin models. We demonstrate how the statistical noise can be reduced considerably by a similarity transformation of the transfer matrix using a variational estimate of its leading eigenvector, in analogy with a common practice in various quantum Monte Carlo techniques. Here we take the two-dimensional coupled XY-Ising model as an example. Furthermore, we calculate interface free energies of finite three-dimensional O(n) models, for the three cases n=1, 2, and 3. Application of finite-size scaling to the numerical results yields estimates of the critical points of these three models. The statistical precision of the estimates is satisfactory for the modest amount of computer time spent. \textcopyright{} 1996 The American Physical Society.
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).
Using a Monte Carlo procedure with a self-consistent-field boundary condition, the magnetization of a simple cubic classical Heisenberg ferromagnet with nearest-neighbor interactions only is calculated. For reduced temperatures and reduced fields in the range $0.02<~|\frac{T}{{T}_{c}}\ensuremath{-}1|$, HJ< 0.5\frac{H}{J}<~0.5JH​< 0.5 the data of this computer experiment can be described in terms of the "effective" critical exponents $0.31<~\ensuremath{\beta}<~0.35$, $1.33<~\ensuremath{\gamma}<~1.39$, and $4.9<~\ensuremath{\delta}<~5.3$. Although the present data are not very close to the critical point they do obey the homogeneity requirement and determine the scaling function rather precisely. This function agrees very well with the scaling function for the face-centered-cubic classical Heisenberg magnet derived recently by Milo\ifmmode \check{s}\else \v{s}\fi{}evi\ifmmode \acute{c}\else \'{c}\fi{} and Stanley using high-temperature---series-expansion techniques. This agreement supports their hypothesis that neither critical exponents nor the scaling function depend on the lattice structure in the Heisenberg model.
The static critical properties of classical, isotropic Heisenberg ferromagnets have been studied numerically. Using an improved cluster algorithm we performed extensive Monte Carlo simulations on L×L×L simple cubic and body-centered cubic systems with L≤40 over a temperature range covering the critical region. Temperature dependences of thermodynamic quantities were determined with an optimized multiple-histogram method, and both critical temperatures and critical exponents were extracted from finite-size scaling of several different thermodynamic quantities. It was found that within their respective error bars the critical exponents for both ferromagnets agree with each other and with field theoretical predictions.