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

Classical Heisenberg Model

G S Joyce

Published Mar 10, 1967• 196 citations
Publisher

Abstract

Exact expressions for the partition function, spin pair correlation function, and susceptibility of the onedimensional isotropic classical Heisenberg model are obtained in zero external field with cyclic boundary conditions. It is shown that the methods used to derive these results enable the partition functions and susceptibilities of finite clusters of interacting classical spins to be evaluated in terms of the $3n\ensuremath{-}j$ symbols of Wigner. Exact results in one dimension are also obtained for the partition function and susceptibility of a "planar" classical Heisenberg model. In this model the spin vectors interact via a Heisenberg coupling but each spin vector is restricted to lie in a plane.The anisotropic classical Heisenberg model described by the Hamiltonian $\mathcal{H}=\ensuremath{-}\ensuremath{\Sigma}\stackrel{}{(\mathrm{ij})}2({{J}_{\mathrm{ij}}}^{x}{{s}_{i}}^{x}{{s}_{j}}^{x}+{{J}_{\mathrm{ij}}}^{y}{{s}_{i}}^{y}{{s}_{j}}^{y}+{{J}_{\mathrm{ij}}}^{z}{{s}_{i}}^{z}{{s}_{j}}^{z})\ensuremath{-}mH\ensuremath{\Sigma}\stackrel{N}{j=1}{{s}_{i}}^{z},$ where ${{s}_{i}}^{x}$, ${{s}_{i}}^{y}$, and ${{s}_{i}}^{z}$ are components of the unit vector ${\mathbf{s}}_{i}$, is also considered. A perturbation series for the zero-field free energy of the anisotropic model in one dimension with nearest-neighbor interactions ${{J}_{\mathrm{ij}}}^{x}={{J}_{\mathrm{ij}}}^{y}=J$ and ${{J}_{\mathrm{ij}}}^{z}=\ensuremath{\gamma}J$ is developed in powers of $\ensuremath{\gamma}\ensuremath{-}1$ using the isotropic model as the unperturbed system. Detailed calculations are performed to third order in $\ensuremath{\gamma}\ensuremath{-}1$. It is found that the perturbation series for the energy per spin breaks down as $T\ensuremath{\rightarrow}0$. A high-temperature series expansion for the anisotropic model, which is valid for a general interaction potential and lattice, is derived by generalizing the methods developed by Horwitz and Callen for the Ising model. This series is rearranged to give a simplified diagram expansion. Finally, a practical technique for calculating the high-temperature series expansions of the zero-field free energy and susceptibility of the isotropic classical Heisenberg model is presented.

DOI
10.1103/physrev.155.478

Related papers

Diagrammatic representation of the two-spin correlation function for the generalized heisenberg model

Eugene Stanley, M. Howard Lee

We consider a “generalized Heisenberg model”, with a Hamiltonian given by ℋ︁(D) = –JΣij S(D)j, where the quantities Si are isotropically-interacting D-dimensional classical spins and the summation is restricted to nearest-neighbor pairs of sites <i,j>. Thus when D = 1, 2, and 3, this model reduces to the S = 1/2 Ising, plane rotator, and classical Heisenberg models which are useful for describing phase transitions in, respectively, liquid-gas systems, Bose fluids, and magnetic materials. Unfortunately this Hamiltonian is exactly soluble only in the limit D → ∞, in which case the solution corresponds to the spherical model of Berlin and Kac. For finite D, we must consider approximation procedures in order to obtain any useful information. In this paper we develop a diagrammatic representation, valid for arbitrary D, for successive terms in the high-temperature expansions of the free energy and the spin-spin correlation function—from which we can obtain series expansions for all the thermodynamic properties of the system.

International Journal of Quantum Chemistry•Jun 18, 2009•3 citations

High-Temperature Critical Indices for the Classical Anisotropic Heisenberg Model

David Jasnow, Michael Wortis

Using the vertex renormalized form of the technique developed by Englert,1 high-temperature series expansions for the spin-spin correlation function of the classical anisotropic Heisenberg model are calculated for various lattices and anisotropies through order T−8 (close-packed lattices) and T−9 (loose-packed lattices). These series are combined and then extrapolated to give the high-temperature critical indices γ (zero-field susceptibility), ν (correlation range), and α (specific heat) as functions of anisotropy. The results are consistent with the hypothesis that the critical indices only change when there is a change in the symmetry of the system, e.g., in interpolating between the Ising and isotropic Heisenberg models, indices remain Ising-like until the system becomes isotropic, at which point they appear to change discontinuously. Previous results for the limiting cases are confirmed and extended. Series for the limiting cases (spin-infinity Ising, XY and isotropic Heisenberg models) as well as for the spin-½ Ising model2 have been derived for the B-site spinel lattice. Analysis of these series shows that the anomalies in the critical indices, reported on the basis of six-term expansions,3 disappear with the use of longer series.

Journal of Applied Physics•Mar 1, 1969•2 citations

Exact Results for the One-Dimensional, Anisotropic Classical Heisenberg Model

G S Joyce

We obtain the zero-field partition function and correlation functions of the one-dimensional, anisotropic classical Heisenberg model exactly and find similar results for an anisotropic "planar" spin model. The isotropic and the anisotropic cases have greatly different low-temperature behavior.

Physical Review Letters•Sep 4, 1967•61 citations

Statistical Mechanics of the Finite Heisenberg Model. I

Chikao Kawabata, Masuo Suzuki

The eigen-functions and eigen-values of finite three-dimensional Heisenberg models with anisotropy constant as a parameter (spin=1/2) were obtained exactly with the aid of group theoretical technique and by using a high speed computer. The calculated eigenvalues are identical with those obtained by Serber and Dresselhaus for the isotropic simple cubic 2×2×2 lattice. However, it was found that some terms in their full expressions were to be corrected. The zeros of the partition function of the above 2×2×2 lattice have been investigated as functions of the anisotropy parameter both in the complex temperature plane and in the complex magnetic field plane. Thermodynamic functions such as energy, specific heat, magnetization, and susceptibility have been numerically computed for these finite system.

Journal of the Physical Society of Japan•Jan 1, 1970•14 citations

Critical Indices for a System of Spins of Arbitrary Dimensionality Situated on a Lattice of Arbitrary Dimensionality

Eugene Stanley

A model Hamiltonian H=−JΣ(ij)Si·Sj is considered, where Si are isotropically interacting D-dimensional classical spins. Thus the S=½ Ising, classical planar, and classical Heisenberg models are obtained if D=1, 2, and 3, respectively; moreover, as D→∞, the free energy corresponding to H approaches the free energy of the exactly-soluble spherical model of Berlin and Kac. H can also be solved exactly providing the lattice dimensionality d=1. For lattices of higher dimensionality, we must resort to approximation techniques; hence we calculate high-temperature expansions for arbitrary D and d of the susceptibility, internal energy (specific heat), and the second moment. (With a few exceptions, such expansions have been available in the past only for D=1, 3 and d=1–3.) Somewhat smoother series to facilitate more reliable extrapolation are then obtained by re-expanding all of the series in terms of the nearest-neighbor spin correlation function for a one-dimensional lattice. We then study the dependence on spin and lattice dimensionality of the critical temperatures Tc(D, d) and the critical exponents γ(D, d) [(χ∼(T-Tc)−γ], α(D, d) [C∼(T-Tc)−α], ν(D, d) [inverse correlation range κ∼(T-Tc)−v]. These ``critical properties'' are all found to vary monotonically (and in most cases, smoothly) with D and d. The assumptions behind making critical phenomena predictions from a dozen or so terms (of an infinite series!) are given increased plausibility by the fact that we have extended the various expansions for the spherical model to 80–100 terms and find that none of the predictions based upon ∼10 terms is changed [e.g., Tc(∞, 2)=0, γ(∞, 3)=2, etc.].

Journal of Applied Physics•Mar 1, 1969•34 citations

High-Temperature Critical Indices for the Classical Anisotropic Heisenberg Model

David Jasnow, Michael Wortis

High-temperature series expansions for the spin-spin correlation function of the classical anisotropic Heinsenberg model are calculated for various lattices and anisotropies through order ${T}^{\ensuremath{-}8}$ (close-packed lattices) and ${T}^{\ensuremath{-}9}$ (loose-packed lattices). These series are combined and then extrapolated to give the high-temperature critical indices $\ensuremath{\gamma}$ (susceptibility), $\ensuremath{\nu}$ (correlation range), and $\ensuremath{\alpha}$ (specific heat) as functions of anisotropy. Our results are consistent with the hypothesis that the critical indices change only when there is a change in the symmetry of the system, e.g., in interpolating between the Ising and isotropic Heisenberg models, indices remain Ising-like until the system becomes isotropic, at which point they appear to change discontinuously. Previous results for the limiting cases are confirmed and extended.

Physical Review•Dec 10, 1968•196 citations
Exact expressions for the partition function, spin pair correlation function, and susceptibility of the onedimensional isotropic classical Heisenberg model are obtained in zero external field with cyclic boundary conditions. It is shown that the methods used to derive these results enable the partition functions and susceptibilities of finite clusters of interacting classical spins to be evaluated in terms of the 3n\ensuremath−j3n\ensuremath{-}j3n\ensuremath−j symbols of Wigner. Exact results in one dimension are also obtained for the partition function and susceptibility of a "planar" classical Heisenberg model. In this model the spin vectors interact via a Heisenberg coupling but each spin vector is restricted to lie in a plane.The anisotropic classical Heisenberg model described by the Hamiltonian H=\ensuremath−\ensuremathΣ(ij)2(Jijxsixsjx+Jijysiysjy+Jijzsizsjz)\ensuremath−mH\ensuremathΣj=1Nsiz,\mathcal{H}=\ensuremath{-}\ensuremath{\Sigma}\stackrel{}{(\mathrm{ij})}2({{J}_{\mathrm{ij}}}^{x}{{s}_{i}}^{x}{{s}_{j}}^{x}+{{J}_{\mathrm{ij}}}^{y}{{s}_{i}}^{y}{{s}_{j}}^{y}+{{J}_{\mathrm{ij}}}^{z}{{s}_{i}}^{z}{{s}_{j}}^{z})\ensuremath{-}mH\ensuremath{\Sigma}\stackrel{N}{j=1}{{s}_{i}}^{z},H=\ensuremath−\ensuremathΣ(ij)​2(Jij​xsi​xsj​x+Jij​ysi​ysj​y+Jij​zsi​zsj​z)\ensuremath−mH\ensuremathΣj=1N​si​z, where six{{s}_{i}}^{x}si​x, siy{{s}_{i}}^{y}si​y, and siz{{s}_{i}}^{z}si​z are components of the unit vector si{\mathbf{s}}_{i}si​, is also considered. A perturbation series for the zero-field free energy of the anisotropic model in one dimension with nearest-neighbor interactions Jijx=Jijy=J{{J}_{\mathrm{ij}}}^{x}={{J}_{\mathrm{ij}}}^{y}=JJij​x=Jij​y=J and Jijz=\ensuremathγJ{{J}_{\mathrm{ij}}}^{z}=\ensuremath{\gamma}JJij​z=\ensuremathγJ is developed in powers of \ensuremathγ\ensuremath−1\ensuremath{\gamma}\ensuremath{-}1\ensuremathγ\ensuremath−1 using the isotropic model as the unperturbed system. Detailed calculations are performed to third order in \ensuremathγ\ensuremath−1\ensuremath{\gamma}\ensuremath{-}1\ensuremathγ\ensuremath−1. It is found that the perturbation series for the energy per spin breaks down as T\ensuremath→0T\ensuremath{\rightarrow}0T\ensuremath→0. A high-temperature series expansion for the anisotropic model, which is valid for a general interaction potential and lattice, is derived by generalizing the methods developed by Horwitz and Callen for the Ising model. This series is rearranged to give a simplified diagram expansion. Finally, a practical technique for calculating the high-temperature series expansions of the zero-field free energy and susceptibility of the isotropic classical Heisenberg model is presented.
We consider a “generalized Heisenberg model”, with a Hamiltonian given by ℋ︁(D) = –JΣij S(D)j, where the quantities Si are isotropically-interacting D-dimensional classical spins and the summation is restricted to nearest-neighbor pairs of sites <i,j>. Thus when D = 1, 2, and 3, this model reduces to the S = 1/2 Ising, plane rotator, and classical Heisenberg models which are useful for describing phase transitions in, respectively, liquid-gas systems, Bose fluids, and magnetic materials. Unfortunately this Hamiltonian is exactly soluble only in the limit D → ∞, in which case the solution corresponds to the spherical model of Berlin and Kac. For finite D, we must consider approximation procedures in order to obtain any useful information. In this paper we develop a diagrammatic representation, valid for arbitrary D, for successive terms in the high-temperature expansions of the free energy and the spin-spin correlation function—from which we can obtain series expansions for all the thermodynamic properties of the system.
Using the vertex renormalized form of the technique developed by Englert,1 high-temperature series expansions for the spin-spin correlation function of the classical anisotropic Heisenberg model are calculated for various lattices and anisotropies through order T−8 (close-packed lattices) and T−9 (loose-packed lattices). These series are combined and then extrapolated to give the high-temperature critical indices γ (zero-field susceptibility), ν (correlation range), and α (specific heat) as functions of anisotropy. The results are consistent with the hypothesis that the critical indices only change when there is a change in the symmetry of the system, e.g., in interpolating between the Ising and isotropic Heisenberg models, indices remain Ising-like until the system becomes isotropic, at which point they appear to change discontinuously. Previous results for the limiting cases are confirmed and extended. Series for the limiting cases (spin-infinity Ising, XY and isotropic Heisenberg models) as well as for the spin-½ Ising model2 have been derived for the B-site spinel lattice. Analysis of these series shows that the anomalies in the critical indices, reported on the basis of six-term expansions,3 disappear with the use of longer series.
The eigen-functions and eigen-values of finite three-dimensional Heisenberg models with anisotropy constant as a parameter (spin=1/2) were obtained exactly with the aid of group theoretical technique and by using a high speed computer. The calculated eigenvalues are identical with those obtained by Serber and Dresselhaus for the isotropic simple cubic 2×2×2 lattice. However, it was found that some terms in their full expressions were to be corrected. The zeros of the partition function of the above 2×2×2 lattice have been investigated as functions of the anisotropy parameter both in the complex temperature plane and in the complex magnetic field plane. Thermodynamic functions such as energy, specific heat, magnetization, and susceptibility have been numerically computed for these finite system.
A model Hamiltonian H=−JΣ(ij)Si·Sj is considered, where Si are isotropically interacting D-dimensional classical spins. Thus the S=½ Ising, classical planar, and classical Heisenberg models are obtained if D=1, 2, and 3, respectively; moreover, as D→∞, the free energy corresponding to H approaches the free energy of the exactly-soluble spherical model of Berlin and Kac. H can also be solved exactly providing the lattice dimensionality d=1. For lattices of higher dimensionality, we must resort to approximation techniques; hence we calculate high-temperature expansions for arbitrary D and d of the susceptibility, internal energy (specific heat), and the second moment. (With a few exceptions, such expansions have been available in the past only for D=1, 3 and d=1–3.) Somewhat smoother series to facilitate more reliable extrapolation are then obtained by re-expanding all of the series in terms of the nearest-neighbor spin correlation function for a one-dimensional lattice. We then study the dependence on spin and lattice dimensionality of the critical temperatures Tc(D, d) and the critical exponents γ(D, d) [(χ∼(T-Tc)−γ], α(D, d) [C∼(T-Tc)−α], ν(D, d) [inverse correlation range κ∼(T-Tc)−v]. These ``critical properties'' are all found to vary monotonically (and in most cases, smoothly) with D and d. The assumptions behind making critical phenomena predictions from a dozen or so terms (of an infinite series!) are given increased plausibility by the fact that we have extended the various expansions for the spherical model to 80–100 terms and find that none of the predictions based upon ∼10 terms is changed [e.g., Tc(∞, 2)=0, γ(∞, 3)=2, etc.].
High-temperature series expansions for the spin-spin correlation function of the classical anisotropic Heinsenberg model are calculated for various lattices and anisotropies through order ${T}^{\ensuremath{-}8}$ (close-packed lattices) and ${T}^{\ensuremath{-}9}$ (loose-packed lattices). These series are combined and then extrapolated to give the high-temperature critical indices $\ensuremath{\gamma}$ (susceptibility), $\ensuremath{\nu}$ (correlation range), and $\ensuremath{\alpha}$ (specific heat) as functions of anisotropy. Our results are consistent with the hypothesis that the critical indices change only when there is a change in the symmetry of the system, e.g., in interpolating between the Ising and isotropic Heisenberg models, indices remain Ising-like until the system becomes isotropic, at which point they appear to change discontinuously. Previous results for the limiting cases are confirmed and extended.