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

Europhysics Letters (EPL)

Analytical Solution of the 2 d Classical Heisenberg Model

Jacques Curély

Published Nov 20, 1995• 60 citations
Publisher

Abstract

For the first time we exactly establish the term of higher degree in the polynomial expansion of the zero-field partition function. In the infinite-layer limit, we derive analytical expressions for the spin-spin correlations, the susceptibility and the specific heat, and we examine their respective behaviours near the critical temperature T c = 0 K. We calculate the analytical values of critical exponents α, γ, η and ν.

DOI
10.1209/0295-5075/32/6/012

Related papers

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

Study of the magnetic properties of two-dimensional (2D) classical square heisenberg antiferromagnets II- spin correlations and susceptibility

Jacques Curély

In this part labelled II we examine the spin correlations and the susceptibility. We use a similar method which has allowed to derive a closed-form expression of the zero-field partition function ZN(0), for 2D square lattices composed of (2N+1)2classical spins isotropically coupled [1]. We rigorously show that the spin correlation vanishes in the zero-field limit, except at T=0 K. Thus, the critical temperature is TC=0 K, in agreement with Mermin-Wagner's theorem. For calculating the spin-spin correlation, we show that it is necessary to distinguish a correlation domain in which the correlation path is confined and a wing domain (Theorem 1). In the thermodynamic limit (N→+∞), we give a general closed-form expression for the spin-spin correlation between any two lattice sites. We prove that all the possible paths have the same analytic expression and correspond to the shortest ones in agreement with the classical principle of least action and its quantum version (Theorem 2). As a result and for the first time, we derive the closed-form expression for the susceptibility, without any approximation. We finally test previous experimental fits and we show that the use of a truncated expansion for the susceptibility was totally justified.

International Conference on Optimization of Electrical and Electronic Equipment•May 1, 2010•2 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

Critical Behavior of the Spin System with Anisotropic Exchange Interaction. II. Two-Dimensional Lattice

Takuma Ishikawa, Takehiko Oguchi

The phase transition for the spin systems described by the Hamiltonian H =-2 J ∑ < i , j > (ξ S i x S j x +η S i y S j y + S i z S j z )- g µ H ∑ i S i z is studied by means of the high-temperature series expansion of the zero field susceptibility. The first five terms in the susceptibility series are explicitly obtained for the two-dimensional square lattice in the cases of (i) 0≤ξ≤1, η=0 and (ii) 0≤ξ=η≤1. The critical temperature T c and the critical exponent γ (susceptibility) are obtained as functions of the spin S and the anisotropy parameter ξ. The critical temperature has a finite value in the XY model (ξ=1, η=0) as well as the Heisenberg model (ξ=η=1). The simple formulae expressing T c and γ valid for arbitrary S and ξ are proposed.

Journal of the Physical Society of Japan•Oct 1, 1971•29 citations

Classical Heisenberg Model

G S Joyce

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.

Physical Review•Mar 10, 1967•196 citations

Exponent Inequalities for Critical Point Spin Correlation Functions

Michael Ellis Fisher

Simple but rigorous derivations are presented for the inequalities2−η≤dγ′/(2β+γ′)≤dγ′/(2−α′),2−η≤d(δ−1)/(δ+1),γ≤(2−η)ν.(Proofs of the first two inequalities have been given by Buckingham and Gunton.1) Here d is the dimensionality and α′, β, γ, γ′ and δ are the exponents characterizing a ferromagnet near its critical point,2 while η and ν describe the decay of the corresponding spin-spin correlation functions2 Γ2(r)=〈S0zSrz〉, according to ∑ |r|<RΓ2(r)∼R2−η(T=Tc,H=0,R→∞),κ1∼(T−Tc)ν,(T→Tc,H=0),where the inverse range of correlation κ1 may be defined by κ1−2=[Σrr2Γ2(r)]/[ΣrΓ2(r)]. The corresponding exponent ηE for the energy-energy correlation function (essentially Γ4(r, a)−[Γ(a)]2, where Γ4(r, a)=〈S0zSazSrzSr+az〉) is proved to satisfy 2−ηE≤dα′,2−ηE≤dαc/(1+ζ+αc),where the specific heat CM at T=Tc diverges with the magnetization as M−αc while the energy derivative | ∂U/∂M |Tc varies as Mζ. (Mean Field or classical values are αc=0, ζ=1.) The proofs are based on the ``intuitively obvious'' properties: (a) positivity, namely, Γ1=〈S0z〉≥0,Γ(r)−Γ12≥0,Γ4(r, a)−[Γ2(a)]2≥0, and (b) monotonic increase of Γ1, Γ2 and Γ4 under increase of magnetic field and decrease of temperature. These properties are known3 to be rigorously valid for Ising models of arbitrary spin, lattice structure and ferromagnetic coupling (Jij≥0). Their validity for real magnets together with the experimental observation δ≤4.7 leads to the significant conclusion η≥0.05.

Journal of Applied Physics•Mar 1, 1969•8 citations
For the first time we exactly establish the term of higher degree in the polynomial expansion of the zero-field partition function. In the infinite-layer limit, we derive analytical expressions for the spin-spin correlations, the susceptibility and the specific heat, and we examine their respective behaviours near the critical temperature T c = 0 K. We calculate the analytical values of critical exponents α, γ, η and ν.
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.
In this part labelled II we examine the spin correlations and the susceptibility. We use a similar method which has allowed to derive a closed-form expression of the zero-field partition function ZN(0), for 2D square lattices composed of (2N+1)2classical spins isotropically coupled [1]. We rigorously show that the spin correlation vanishes in the zero-field limit, except at T=0 K. Thus, the critical temperature is TC=0 K, in agreement with Mermin-Wagner's theorem. For calculating the spin-spin correlation, we show that it is necessary to distinguish a correlation domain in which the correlation path is confined and a wing domain (Theorem 1). In the thermodynamic limit (N→+∞), we give a general closed-form expression for the spin-spin correlation between any two lattice sites. We prove that all the possible paths have the same analytic expression and correspond to the shortest ones in agreement with the classical principle of least action and its quantum version (Theorem 2). As a result and for the first time, we derive the closed-form expression for the susceptibility, without any approximation. We finally test previous experimental fits and we show that the use of a truncated expansion for the susceptibility was totally justified.
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.].
The phase transition for the spin systems described by the Hamiltonian H =-2 J ∑ < i , j > (ξ S i x S j x +η S i y S j y + S i z S j z )- g µ H ∑ i S i z is studied by means of the high-temperature series expansion of the zero field susceptibility. The first five terms in the susceptibility series are explicitly obtained for the two-dimensional square lattice in the cases of (i) 0≤ξ≤1, η=0 and (ii) 0≤ξ=η≤1. The critical temperature T c and the critical exponent γ (susceptibility) are obtained as functions of the spin S and the anisotropy parameter ξ. The critical temperature has a finite value in the XY model (ξ=1, η=0) as well as the Heisenberg model (ξ=η=1). The simple formulae expressing T c and γ valid for arbitrary S and ξ are proposed.
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 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 ${{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.
Simple but rigorous derivations are presented for the inequalities2−η≤dγ′/(2β+γ′)≤dγ′/(2−α′),2−η≤d(δ−1)/(δ+1),γ≤(2−η)ν.(Proofs of the first two inequalities have been given by Buckingham and Gunton.1) Here d is the dimensionality and α′, β, γ, γ′ and δ are the exponents characterizing a ferromagnet near its critical point,2 while η and ν describe the decay of the corresponding spin-spin correlation functions2 Γ2(r)=〈S0zSrz〉, according to ∑ |r|<RΓ2(r)∼R2−η(T=Tc,H=0,R→∞),κ1∼(T−Tc)ν,(T→Tc,H=0),where the inverse range of correlation κ1 may be defined by κ1−2=[Σrr2Γ2(r)]/[ΣrΓ2(r)]. The corresponding exponent ηE for the energy-energy correlation function (essentially Γ4(r, a)−[Γ(a)]2, where Γ4(r, a)=〈S0zSazSrzSr+az〉) is proved to satisfy 2−ηE≤dα′,2−ηE≤dαc/(1+ζ+αc),where the specific heat CM at T=Tc diverges with the magnetization as M−αc while the energy derivative | ∂U/∂M |Tc varies as Mζ. (Mean Field or classical values are αc=0, ζ=1.) The proofs are based on the ``intuitively obvious'' properties: (a) positivity, namely, Γ1=〈S0z〉≥0,Γ(r)−Γ12≥0,Γ4(r, a)−[Γ2(a)]2≥0, and (b) monotonic increase of Γ1, Γ2 and Γ4 under increase of magnetic field and decrease of temperature. These properties are known3 to be rigorously valid for Ising models of arbitrary spin, lattice structure and ferromagnetic coupling (Jij≥0). Their validity for real magnets together with the experimental observation δ≤4.7 leads to the significant conclusion η≥0.05.