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

Nuclear Physics B

Statistical thermodynamics of the quantum spin-s XXZ Heisenberg model in the Weiss limit

J. Ricardo G. Mendonça

Published Oct 8, 2026• 0 citations
Publisher

Abstract

The statistical thermodynamics of the quantum spin- s XXZ Heisenberg model in the Weiss (mean-field) limit rests on the outer multiplicities of the irreducible representations in the Clebsch–Gordan decomposition of the n -fold spin- s tensor product, which determine the energy-level degeneracies and hence the partition function. Using known closed forms for these multiplicities in terms of generalized binomial coefficients of order 2 s + 1 , we write the exact finite-size partition function of the effective collective Hamiltonian and evaluate its thermodynamic limit by a large-deviation (Laplace) analysis validated against direct summations for systems of up to n = 10 4 spins. The free energy then follows from a two-parameter variational principle in the total spin and the magnetization per site. The ferromagnetic ordering is driven by the rotationally invariant part of the Hamiltonian, so that the Curie temperature T c = 1 3 J s ( s + 1 ) is independent of the anisotropy Δ, and the transition is of second order for every spin s , with classical mean-field exponents. The anisotropy reshapes the phase diagram: for Δ < 1 the field-driven first-order line of the isotropic ( Δ = 1 ) model opens into a dome of easy-plane ferromagnetic order, bounded by a continuous transition that restores the residual U(1) symmetry, and no first-order transitions or metastability occur throughout − 1 ≤ Δ < 1 ; the isotropic Curie-Weiss ferromagnet is recovered as Δ → 1. In two appendices we detail the asymptotics of the outer multiplicities down to the Gaussian-fluctuation prefactor and discuss the q -number and Jacobi-theta structures underlying the model, including the partial-theta-type structure that obstructs a standard theta-function reduction of the residual total-spin sum.

DOI
10.1016/j.nuclphysb.2026.117688

Related papers

Quantum spin-1 anisotropic ferromagnetic Heisenberg model in a crystal field: A variational approach

D. C. Carvalho, J. A. Plascak, L. M. Castro

A variational approach based on Bogoliubov inequality for the free energy is employed in order to treat the quantum spin-1 anisotropic ferromagnetic Heisenberg model in the presence of a crystal field. Within the Bogoliubov scheme an improved pair approximation has been used. The temperature-dependent thermodynamic functions have been obtained and provide much better results than the previous simple mean-field scheme. In one dimension, which is still nonintegrable for quantum spin-1, we get the exact results in the classical limit, or near-exact results in the quantum case, for the free energy, magnetization, and quadrupole moment, as well for the transition temperature. In two and three dimensions the corresponding global phase diagrams have been obtained as a function of the parameters of the Hamiltonian. First-order transition lines, second-order transition lines, tricritical and tetracritical points, and critical endpoints have been located through the analysis of the minimum of the Helmholtz free energy and a Landau-like expansion in the approximated free energy. Only first-order quantum transitions have been found at zero temperature. Limiting cases, such as isotropic Heisenberg, Blume-Capel, and Ising models, have been analyzed and compared to previous results obtained from other analytical approaches as well as from Monte Carlo simulations.

arXiv:1311.6851•Physical Review E•Sep 6, 2013•10 citations

Statistical Mechanical Theory of Ferromagnetism. High Density Behavior

R. Brout

The partition function of the Ising model of ferromagnetism is examined in the limit of high density in the anticipation that in the limit of infinite density one recovers the Weiss molecular field. The formal parameter of expansion is $\frac{1}{z}$ where $z$ is the number of spins in the range of the exchange potential (not restricted to nearest neighbor interactions). In the absence of long-range order, only ring diagrams in the cluster expansion contribute. These give a divergence in the specific heat at $k{T}_{c}=\ensuremath{\Sigma}{j\ensuremath{\ne}i}^{}{v}_{\mathrm{ij}}$ where ${v}_{\mathrm{ij}}$ is the exchange potential. This is the molecular field value for the Curie point ${T}_{c}$. In the presence of a magnetic field the partition function is evaluated for fixed magnetic moment $M$ in the same approximation, $M$ being determined by minimization. This results in a susceptibility differing from the molecular field theory and hence an inconsistency in the theory.The inconsistency is traced back to the observation that the acceptance of ring diagrams is equivalent to the gaussian model of Kac and Berlin which violates the sum rule ${{\ensuremath{\Sigma}}_{i=1}}^{N}{{\ensuremath{\mu}}_{i}}^{2}=N$. Here ${\ensuremath{\mu}}_{i}$ is the "spin" per particle and $N$ is the total number of particles. This condition is reinstated by insuring the sum rule. The result leads to the spherical model. Thus, a consistent high density approximation to the Ising model is the spherical model. Below the Curie point or for fixed magnetic field, $M$ is again held fixed and only the Fourier components of the spin density with nonvanishing wave vector are "sphericalized." The result leads to a physically acceptable model which becomes the molecular field theory at low temperatures or high fields and deviates in $O(\frac{1}{z})$ in general. Formally, the results are simply expressed in terms of a temperature dependent Weiss field. These results differ from the ordinary spherical model which is physically unacceptable below the Curie point. However, a molecular field modification of the spherical model due to M. Lax yields the same result when properly interpreted.It is shown that the above results are also valid (to the same approximation) in the quantum mechanical Heisenberg model, for temperatures above the Curie point.

Physical Review•May 15, 1960•181 citations

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

Critical behavior of the two-dimensional easy-plane ferromagnet

Alessandro Cuccoli, V. Tognetti, Paola Verrucchi +1 more

The critical behavior of 2D magnetic easy-plane systems has mainly been studied by the classical XY model. However, the z components of the spins have to be considered in order to describe real systems, and their fluctuations cannot positively be neglected when quantum effects are to be included, quantum spins being intrinsically three-component objects. Therefore, Monte Carlo simulations are performed for the Heisenberg ferromagnet with easy-plane anisotropy (XXZ model) on a two-dimensional square lattice with a twofold aim: first, to obtain accurate quantitative results about the critical behavior of the classical model, showing the relevant role played by the out-of-plane fluctuations; second, to open the way for approaching the quantum thermodynamics by means of the effective Hamiltonian method that reduces the quantum thermodynamics of the XXZ ferromagnet to the investigation of an effective classical model with temperature-dependent renormalized interaction parameters. Specific heat, magnetic susceptibility, and correlation length are calculated in the critical region for lattice sizes up to 128×128. These quantities preserve the Kosterlitz–Thouless behavior of the XY model. For the transition temperature of the classical XX0 model we obtain the estimate kBTc/(JS̃2)=0.69±0.01.

Journal of Applied Physics•Nov 15, 1994•10 citations

Lee-Yang theory of the Curie-Weiss model and its rare fluctuations

Aydin Deger, Christian Flindt

Phase transitions are typically accompanied by nonanalytic behaviors of the free energy, which can be explained by considering the zeros of the partition function in the complex plane of the control parameter and their approach to the critical value on the real axis as the system size is increased.Recent experiments have shown that partition function zeros are not just a theoretical concept.They can also be determined experimentally by measuring fluctuations of thermodynamic observables in systems of finite size.Motivated by this progress, we investigate here the partition function zeros for the Curie-Weiss model of spontaneous magnetization using our recently established cumulant method.Specifically, we extract the leading Fisher and Lee-Yang zeros of the Curie-Weiss model from the fluctuations of the energy and the magnetization in systems of finite size.We develop a finite-size scaling analysis of the partition function zeros, which is valid for mean-field models and which allows us to extract both the critical values of the control parameters and the critical exponents, even for small systems that are away from criticality.We also show that the Lee-Yang zeros carry important information about the rare magnetic fluctuations as they allow us to predict many essential features of the large-deviation statistics of the magnetization.This finding may constitute a profound connection between Lee-Yang theory and large-deviation statistics.

arXiv:2002.01269•Physical Review Research•Jul 2, 2020•28 citations

Heisenberg models and a particular isotropic model

Arjan J. van der Sijs

The Heisenberg model, a quantum-mechanical analogue of the Ising model, has a large ground-state degeneracy, due to the symmetry generated by the total spin. This symmetry is also responsible for degeneracies in the rest of the spectrum. We discuss the global structure of the spectrum of Heisenberg models with arbitrary couplings, using group-theoretical methods. The Hilbert space breaks up in blocks characterized by the quantum numbers of the total spin, S and M, and each block is shown to constitute the representation space of an explicitly given irreducible representation of the symmetric group ${\mathit{S}}_{\mathit{N}}$, consisting of permutations of the N spins in the system. In the second part of the paper we consider, as a concrete application, the model where each spin is coupled to all the other spins with equal strength. Its partition function is written as a single integral, elucidating its N dependence. This provides a useful framework for studying finite-size effects. We give explicit results for the heat capacity, revealing interesting behavior just around the phase transition.

arXiv:cond-mat/9306045•Physical review. B, Condensed matter•Sep 1, 1993•10 citations
The statistical thermodynamics of the quantum spin- s XXZ Heisenberg model in the Weiss (mean-field) limit rests on the outer multiplicities of the irreducible representations in the Clebsch–Gordan decomposition of the n -fold spin- s tensor product, which determine the energy-level degeneracies and hence the partition function. Using known closed forms for these multiplicities in terms of generalized binomial coefficients of order 2 s + 1 , we write the exact finite-size partition function of the effective collective Hamiltonian and evaluate its thermodynamic limit by a large-deviation (Laplace) analysis validated against direct summations for systems of up to n = 10 4 spins. The free energy then follows from a two-parameter variational principle in the total spin and the magnetization per site. The ferromagnetic ordering is driven by the rotationally invariant part of the Hamiltonian, so that the Curie temperature T c = 1 3 J s ( s + 1 ) is independent of the anisotropy Δ, and the transition is of second order for every spin s , with classical mean-field exponents. The anisotropy reshapes the phase diagram: for Δ < 1 the field-driven first-order line of the isotropic ( Δ = 1 ) model opens into a dome of easy-plane ferromagnetic order, bounded by a continuous transition that restores the residual U(1) symmetry, and no first-order transitions or metastability occur throughout − 1 ≤ Δ < 1 ; the isotropic Curie-Weiss ferromagnet is recovered as Δ → 1. In two appendices we detail the asymptotics of the outer multiplicities down to the Gaussian-fluctuation prefactor and discuss the q -number and Jacobi-theta structures underlying the model, including the partial-theta-type structure that obstructs a standard theta-function reduction of the residual total-spin sum.
A variational approach based on Bogoliubov inequality for the free energy is employed in order to treat the quantum spin-1 anisotropic ferromagnetic Heisenberg model in the presence of a crystal field. Within the Bogoliubov scheme an improved pair approximation has been used. The temperature-dependent thermodynamic functions have been obtained and provide much better results than the previous simple mean-field scheme. In one dimension, which is still nonintegrable for quantum spin-1, we get the exact results in the classical limit, or near-exact results in the quantum case, for the free energy, magnetization, and quadrupole moment, as well for the transition temperature. In two and three dimensions the corresponding global phase diagrams have been obtained as a function of the parameters of the Hamiltonian. First-order transition lines, second-order transition lines, tricritical and tetracritical points, and critical endpoints have been located through the analysis of the minimum of the Helmholtz free energy and a Landau-like expansion in the approximated free energy. Only first-order quantum transitions have been found at zero temperature. Limiting cases, such as isotropic Heisenberg, Blume-Capel, and Ising models, have been analyzed and compared to previous results obtained from other analytical approaches as well as from Monte Carlo simulations.
The partition function of the Ising model of ferromagnetism is examined in the limit of high density in the anticipation that in the limit of infinite density one recovers the Weiss molecular field. The formal parameter of expansion is 1z\frac{1}{z}z1​ where zzz is the number of spins in the range of the exchange potential (not restricted to nearest neighbor interactions). In the absence of long-range order, only ring diagrams in the cluster expansion contribute. These give a divergence in the specific heat at $k{T}_{c}=\ensuremath{\Sigma}{j\ensuremath{\ne}i}^{}{v}_{\mathrm{ij}}$ where vij{v}_{\mathrm{ij}}vij​ is the exchange potential. This is the molecular field value for the Curie point Tc{T}_{c}Tc​. In the presence of a magnetic field the partition function is evaluated for fixed magnetic moment MMM in the same approximation, MMM being determined by minimization. This results in a susceptibility differing from the molecular field theory and hence an inconsistency in the theory.The inconsistency is traced back to the observation that the acceptance of ring diagrams is equivalent to the gaussian model of Kac and Berlin which violates the sum rule ${{\ensuremath{\Sigma}}_{i=1}}^{N}{{\ensuremath{\mu}}_{i}}^{2}=N$. Here ${\ensuremath{\mu}}_{i}$ is the "spin" per particle and NNN is the total number of particles. This condition is reinstated by insuring the sum rule. The result leads to the spherical model. Thus, a consistent high density approximation to the Ising model is the spherical model. Below the Curie point or for fixed magnetic field, MMM is again held fixed and only the Fourier components of the spin density with nonvanishing wave vector are "sphericalized." The result leads to a physically acceptable model which becomes the molecular field theory at low temperatures or high fields and deviates in O(1z)O(\frac{1}{z})O(z1​) in general. Formally, the results are simply expressed in terms of a temperature dependent Weiss field. These results differ from the ordinary spherical model which is physically unacceptable below the Curie point. However, a molecular field modification of the spherical model due to M. Lax yields the same result when properly interpreted.It is shown that the above results are also valid (to the same approximation) in the quantum mechanical Heisenberg model, for temperatures above the Curie point.
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.
The critical behavior of 2D magnetic easy-plane systems has mainly been studied by the classical XY model. However, the z components of the spins have to be considered in order to describe real systems, and their fluctuations cannot positively be neglected when quantum effects are to be included, quantum spins being intrinsically three-component objects. Therefore, Monte Carlo simulations are performed for the Heisenberg ferromagnet with easy-plane anisotropy (XXZ model) on a two-dimensional square lattice with a twofold aim: first, to obtain accurate quantitative results about the critical behavior of the classical model, showing the relevant role played by the out-of-plane fluctuations; second, to open the way for approaching the quantum thermodynamics by means of the effective Hamiltonian method that reduces the quantum thermodynamics of the XXZ ferromagnet to the investigation of an effective classical model with temperature-dependent renormalized interaction parameters. Specific heat, magnetic susceptibility, and correlation length are calculated in the critical region for lattice sizes up to 128×128. These quantities preserve the Kosterlitz–Thouless behavior of the XY model. For the transition temperature of the classical XX0 model we obtain the estimate kBTc/(JS̃2)=0.69±0.01.
Phase transitions are typically accompanied by nonanalytic behaviors of the free energy, which can be explained by considering the zeros of the partition function in the complex plane of the control parameter and their approach to the critical value on the real axis as the system size is increased.Recent experiments have shown that partition function zeros are not just a theoretical concept.They can also be determined experimentally by measuring fluctuations of thermodynamic observables in systems of finite size.Motivated by this progress, we investigate here the partition function zeros for the Curie-Weiss model of spontaneous magnetization using our recently established cumulant method.Specifically, we extract the leading Fisher and Lee-Yang zeros of the Curie-Weiss model from the fluctuations of the energy and the magnetization in systems of finite size.We develop a finite-size scaling analysis of the partition function zeros, which is valid for mean-field models and which allows us to extract both the critical values of the control parameters and the critical exponents, even for small systems that are away from criticality.We also show that the Lee-Yang zeros carry important information about the rare magnetic fluctuations as they allow us to predict many essential features of the large-deviation statistics of the magnetization.This finding may constitute a profound connection between Lee-Yang theory and large-deviation statistics.
The Heisenberg model, a quantum-mechanical analogue of the Ising model, has a large ground-state degeneracy, due to the symmetry generated by the total spin. This symmetry is also responsible for degeneracies in the rest of the spectrum. We discuss the global structure of the spectrum of Heisenberg models with arbitrary couplings, using group-theoretical methods. The Hilbert space breaks up in blocks characterized by the quantum numbers of the total spin, S and M, and each block is shown to constitute the representation space of an explicitly given irreducible representation of the symmetric group SN{\mathit{S}}_{\mathit{N}}SN​, consisting of permutations of the N spins in the system. In the second part of the paper we consider, as a concrete application, the model where each spin is coupled to all the other spins with equal strength. Its partition function is written as a single integral, elucidating its N dependence. This provides a useful framework for studying finite-size effects. We give explicit results for the heat capacity, revealing interesting behavior just around the phase transition.