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

Some Critical Properties of the Heisenberg Model

Roger Bowers, Michael E. Woolf

Published Jan 10, 1969• 48 citations
Publisher

Abstract

A detailed series-extrapolation investigation of some critical properties ($T\ensuremath{\rightarrow}{T}_{c}+$) of the Heisenberg model is presented. Systematic use is made of ratio and Pad\'e-approximant techniques, some of the latter being new. Attention is confined to those nearest-neighbor and order-two equivalent-neighbor models which are based on the simple cubic, body-centered cubic, and face-centered cubic lattices. For the face-centered cubic nearest-neighbor classical Heisenberg model, the following estimates are obtained for familiar critical exponents: $\ensuremath{\alpha}=0, \ensuremath{-}\frac{1}{8}<~{\ensuremath{\alpha}}_{s}<~\ensuremath{-}\frac{1}{16}$, $\ensuremath{\gamma}=1.375\ifmmode\pm\else\textpm\fi{}0.002$, $0.6875<~\ensuremath{\nu}<~0.7125$, $0<~\ensuremath{\eta}<~0.07$. Estimates for $\ensuremath{\gamma}$ are also obtained for the other nearest-neighbor and equivalent-neighbor classical Heisenberg models, and it is conjectured that $\ensuremath{\gamma}=1\frac{3}{8}$ may be the exact result in all the cases considered. For the face-centered cubic order-two equivalent-neighbor spin-\textonehalf{} Heisenberg ferromagnet, it is estimated that $\ensuremath{\gamma}=1.3744\ifmmode\pm\else\textpm\fi{}0.0008$. A similar but less precise result is found for the corresponding simple cubic model, and the spin-\textonehalf{} nearest-neighbor susceptibility series are briefly examined. It is suggested that consideration should be given to the possibility that for ferromagnetic Heisenberg interactions $\ensuremath{\gamma}$ has the same value $1\frac{3}{8}$ irrespective of spin and the particular three-dimensional nearest-neighbor or finite-order equivalent-neighbor model considered. Other quantities estimated include various Curie points and singularity amplitudes and the critical values of certain thermodynamic functions. Some attention is given to the simple classical Heisenberg antiferromagnet.

DOI
10.1103/physrev.177.917

Related papers

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

Ising models with n>1 : A series-expansion approach

Michael Plischke, J. Oitmaa

We investigate the Ising model with nearest-neighbor ferromagnetic and next-nearest-neighbor antiferromagnetic interactions on the square and body-centered-cubic lattices. For $0<\frac{{J}_{1}}{{J}_{2}}<\frac{3}{2}$ on the bcc lattice the model has two nontrivially degenerate ground states and has been conjectured to have $\mathrm{XY}$ ($n=2$) critical exponents. For $\frac{{J}_{1}}{{J}_{2}}=\frac{3}{2}$ there are three degenerate ground states and Heisenberg critical exponents ($n=3$) are expected. We have obtained nine terms in the high-temperature expansion of the staggered susceptibility. For $0<\frac{{J}_{1}}{{J}_{2}}<~1$ we find $\ensuremath{\gamma}=1.25$ ($n=1$). At $\frac{{J}_{1}}{{J}_{2}}=\frac{3}{2}$ the results are not conclusive, but $\ensuremath{\gamma}=1.33$ ($n=2$) is not ruled out by the series analysis whereas $\ensuremath{\gamma}=1.4$ ($n=3$) is not probable. Results for the square lattice are inconclusive.

Physical review. B, Condensed matter•Jan 1, 1979•16 citations

N -vector spin models on the simple-cubic and the body-centered-cubic lattices: A study of the critical behavior of the susceptibility and of the correlation length by high-temperature series extended to order β21

P. Butera, M. Comi

High-temperature expansions for the free energy, the susceptibility, and the second correlation moment of the classical $N$-vector model [also known as the $O(N)$ symmetric classical spin-Heisenberg model or as the lattice $O(N)$ nonlinear $\ensuremath{\sigma}$ model] on the simple-cubic and the body-centered-cubic lattices are extended to order ${\ensuremath{\beta}}^{21}$ for arbitrary $N$. The series for the second field derivative of the susceptibility is extended to order ${\ensuremath{\beta}}^{17}$. We report here on the analysis of the computed series for the susceptibility and the (second moment) correlation length which yields updated estimates of the critical parameters for various values of the spin dimensionality $N$, including $N=0$ (the self-avoiding walk model), $N=1$ (the Ising spin-1/2 model), $N=2$ (the $\mathrm{XY}$ model), and $N=3$ (the classical Heisenberg model). For all values of $N$ we confirm a good agreement with the present renormalization-group estimates. A study of the series for the other observables will appear in a forthcoming paper.

arXiv:hep-lat/9703018•Physical review. B, Condensed matter•Oct 1, 1997•139 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

Extended high-temperature series for the classical Heisenberg model in three dimensions

Sati McKenzie, Cyril Domb, D L Hunter

Extended series expansions are derived for the high-temperature susceptibility of the classical Heisenberg model, on three-dimensional lattices. Series coefficients are presented to twelfth order for the simple cubic (SC) and face centred cubic (FCC) lattices and to eleventh order for the body centred cubic (BCC) lattice. The results are in agreement with earlier calculations apart from a small discrepancy at the tenth order on the FCC lattice. Extrapolation studies on the extended series are used to obtain revised estimates for the critical points (K c ) and the susceptibility exponent ( gamma ). On the FCC lattice, the authors also investigate the possibility of a confluent non-analytic correction to the dominant singularity. While the coefficients are consistent with the presence of such a correction term with an exponent ( Delta 1 ) of 0.55, as predicted by renormalisation group (RG) calculations, the amplitude of the correction term appears to be very small compared with that of the first analytic correction term.

Journal of Physics A Mathematical and General•Dec 1, 1982•30 citations

High-temperature series for the susceptibility of the spin- S Ising model

John P. VAN DYKE, William J. Camp

We have extended the series for the zero-field susceptibility of the spin-$S$ Ising model to eighth order in the reduced temperature $K$, on the triangular, simple cubic, body-centered-cubic, and face-centered-cubic lattices. The coefficients of these series ${h}_{n}(S)$ are expressed as simple polynomials in $X=S(S+1)$. For the face-centered-cubic lattice, an accurate polynomial fit to the critical point ${K}_{c}(S)$ is presented; and the apparent spin dependence of the critical exponent $\ensuremath{\gamma}$ is briefly discussed. The series are quite well behaved for all $S$. However, the large-$S$ series seems to exhibit more rapid apparent convergence.

Physical review. B, Solid state•Apr 1, 1974•42 citations
A detailed series-extrapolation investigation of some critical properties (T\ensuremath→Tc+T\ensuremath{\rightarrow}{T}_{c}+T\ensuremath→Tc​+) of the Heisenberg model is presented. Systematic use is made of ratio and Pad'e-approximant techniques, some of the latter being new. Attention is confined to those nearest-neighbor and order-two equivalent-neighbor models which are based on the simple cubic, body-centered cubic, and face-centered cubic lattices. For the face-centered cubic nearest-neighbor classical Heisenberg model, the following estimates are obtained for familiar critical exponents: \ensuremathα=0,\ensuremath−18< \ensuremathαs< \ensuremath−116\ensuremath{\alpha}=0, \ensuremath{-}\frac{1}{8}<~{\ensuremath{\alpha}}_{s}<~\ensuremath{-}\frac{1}{16}\ensuremathα=0,\ensuremath−81​< \ensuremathαs​< \ensuremath−161​, \ensuremathγ=1.375\ifmmode±\else\textpm\fi0.002\ensuremath{\gamma}=1.375\ifmmode\pm\else\textpm\fi{}0.002\ensuremathγ=1.375\ifmmode±\else\textpm\fi0.002, 0.6875< \ensuremathν< 0.71250.6875<~\ensuremath{\nu}<~0.71250.6875< \ensuremathν< 0.7125, 0< \ensuremathη< 0.070<~\ensuremath{\eta}<~0.070< \ensuremathη< 0.07. Estimates for \ensuremathγ\ensuremath{\gamma}\ensuremathγ are also obtained for the other nearest-neighbor and equivalent-neighbor classical Heisenberg models, and it is conjectured that \ensuremathγ=138\ensuremath{\gamma}=1\frac{3}{8}\ensuremathγ=183​ may be the exact result in all the cases considered. For the face-centered cubic order-two equivalent-neighbor spin-\textonehalf{} Heisenberg ferromagnet, it is estimated that \ensuremathγ=1.3744\ifmmode±\else\textpm\fi0.0008\ensuremath{\gamma}=1.3744\ifmmode\pm\else\textpm\fi{}0.0008\ensuremathγ=1.3744\ifmmode±\else\textpm\fi0.0008. A similar but less precise result is found for the corresponding simple cubic model, and the spin-\textonehalf{} nearest-neighbor susceptibility series are briefly examined. It is suggested that consideration should be given to the possibility that for ferromagnetic Heisenberg interactions \ensuremathγ\ensuremath{\gamma}\ensuremathγ has the same value 1381\frac{3}{8}183​ irrespective of spin and the particular three-dimensional nearest-neighbor or finite-order equivalent-neighbor model considered. Other quantities estimated include various Curie points and singularity amplitudes and the critical values of certain thermodynamic functions. Some attention is given to the simple classical Heisenberg antiferromagnet.
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.
We investigate the Ising model with nearest-neighbor ferromagnetic and next-nearest-neighbor antiferromagnetic interactions on the square and body-centered-cubic lattices. For 0<J1J2<320<\frac{{J}_{1}}{{J}_{2}}<\frac{3}{2}0<J2​J1​​<23​ on the bcc lattice the model has two nontrivially degenerate ground states and has been conjectured to have XY\mathrm{XY}XY (n=2n=2n=2) critical exponents. For J1J2=32\frac{{J}_{1}}{{J}_{2}}=\frac{3}{2}J2​J1​​=23​ there are three degenerate ground states and Heisenberg critical exponents (n=3n=3n=3) are expected. We have obtained nine terms in the high-temperature expansion of the staggered susceptibility. For 0<J1J2< 10<\frac{{J}_{1}}{{J}_{2}}<~10<J2​J1​​< 1 we find $\ensuremath{\gamma}=1.25$ (n=1n=1n=1). At J1J2=32\frac{{J}_{1}}{{J}_{2}}=\frac{3}{2}J2​J1​​=23​ the results are not conclusive, but $\ensuremath{\gamma}=1.33$ (n=2n=2n=2) is not ruled out by the series analysis whereas $\ensuremath{\gamma}=1.4$ (n=3n=3n=3) is not probable. Results for the square lattice are inconclusive.
High-temperature expansions for the free energy, the susceptibility, and the second correlation moment of the classical NNN-vector model [also known as the O(N)O(N)O(N) symmetric classical spin-Heisenberg model or as the lattice O(N)O(N)O(N) nonlinear $\ensuremath{\sigma}$ model] on the simple-cubic and the body-centered-cubic lattices are extended to order ${\ensuremath{\beta}}^{21}$ for arbitrary NNN. The series for the second field derivative of the susceptibility is extended to order ${\ensuremath{\beta}}^{17}$. We report here on the analysis of the computed series for the susceptibility and the (second moment) correlation length which yields updated estimates of the critical parameters for various values of the spin dimensionality NNN, including N=0N=0N=0 (the self-avoiding walk model), N=1N=1N=1 (the Ising spin-1/2 model), N=2N=2N=2 (the XY\mathrm{XY}XY model), and N=3N=3N=3 (the classical Heisenberg model). For all values of NNN we confirm a good agreement with the present renormalization-group estimates. A study of the series for the other observables will appear in a forthcoming paper.
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).
Extended series expansions are derived for the high-temperature susceptibility of the classical Heisenberg model, on three-dimensional lattices. Series coefficients are presented to twelfth order for the simple cubic (SC) and face centred cubic (FCC) lattices and to eleventh order for the body centred cubic (BCC) lattice. The results are in agreement with earlier calculations apart from a small discrepancy at the tenth order on the FCC lattice. Extrapolation studies on the extended series are used to obtain revised estimates for the critical points (K c ) and the susceptibility exponent ( gamma ). On the FCC lattice, the authors also investigate the possibility of a confluent non-analytic correction to the dominant singularity. While the coefficients are consistent with the presence of such a correction term with an exponent ( Delta 1 ) of 0.55, as predicted by renormalisation group (RG) calculations, the amplitude of the correction term appears to be very small compared with that of the first analytic correction term.
We have extended the series for the zero-field susceptibility of the spin-SSS Ising model to eighth order in the reduced temperature KKK, on the triangular, simple cubic, body-centered-cubic, and face-centered-cubic lattices. The coefficients of these series hn(S){h}_{n}(S)hn​(S) are expressed as simple polynomials in X=S(S+1)X=S(S+1)X=S(S+1). For the face-centered-cubic lattice, an accurate polynomial fit to the critical point Kc(S){K}_{c}(S)Kc​(S) is presented; and the apparent spin dependence of the critical exponent $\ensuremath{\gamma}$ is briefly discussed. The series are quite well behaved for all SSS. However, the large-SSS series seems to exhibit more rapid apparent convergence.