Quantum Heisenberg model

Statistical model in quantum mechanics of magnetic materials

The quantum Heisenberg model, developed by Werner Heisenberg, is a statistical mechanical model used in the study of critical points and phase transitions of magnetic systems, in which the spins of the magnetic systems are treated quantum mechanically. It is related to the prototypical Ising model, where at each site of a lattice, a spin  σ i ∈ { ± 1 } {\displaystyle \sigma _{i}\in \\{\pm 1\\}} represents a microscopic magnetic dipole to which the magnetic moment is either up or down. Except the coupling between magnetic dipole moments, there is also a multipolar version of Heisenberg model called the multipolar exchange interaction.

## Overview

[edit]

For quantum mechanical reasons (see exchange interaction or Magnetism § Quantum-mechanical origin of magnetism), the dominant coupling between two dipoles may cause nearest-neighbors to have lowest energy when they are aligned. Under this assumption (so that magnetic interactions only occur between adjacent dipoles) and on a 1-dimensional periodic lattice, the Hamiltonian can be written in the form

H ^ = − J ∑ j = 1 N σ j σ j + 1 − h ∑ j = 1 N σ j {\displaystyle {\hat {H}}=-J\sum _{j=1}^{N}\sigma _{j}\sigma _{j+1}-h\sum _{j=1}^{N}\sigma _{j}} ,

where  J {\displaystyle J} is the coupling constant and dipoles are represented by classical vectors (or "spins") σj, subject to the periodic boundary condition  σ N + 1 = σ 1 {\displaystyle \sigma _{N+1}=\sigma _{1}} . The Heisenberg model is a more realistic model in that it treats the spins quantum-mechanically, by replacing the spin by a quantum operator acting upon the tensor product ( C 2 ) ⊗ N {\displaystyle (\mathbb {C} ^{2})^{\otimes N}} , of dimension  2 N {\displaystyle 2^{N}} . To define it, recall the Pauli spin-1/2 matrices

σ x = ( 0 1 1 0 ) {\displaystyle \sigma ^{x}={\begin{pmatrix}0&1\\\1&0\end{pmatrix}}} ,

σ y = ( 0 − i i 0 ) {\displaystyle \sigma ^{y}={\begin{pmatrix}0&-i\\\i&0\end{pmatrix}}} ,

σ z = ( 1 0 0 − 1 ) {\displaystyle \sigma ^{z}={\begin{pmatrix}1&0\\\0&-1\end{pmatrix}}} ,

and for  1 ≤ j ≤ N {\displaystyle 1\leq j\leq N} and  a ∈ { x , y , z } {\displaystyle a\in \\{x,y,z\\}} denote  σ j a = I ⊗ j − 1 ⊗ σ a ⊗ I ⊗ N − j {\displaystyle \sigma _{j}^{a}=I^{\otimes j-1}\otimes \sigma ^{a}\otimes I^{\otimes N-j}} , where  I {\displaystyle I} is the  2 × 2 {\displaystyle 2\times 2} identity matrix. Given a choice of real-valued coupling constants  J x , J y , {\displaystyle J_{x},J_{y},} and  J z {\displaystyle J_{z}} , the Hamiltonian is given by

H ^ = − 1 2 ∑ j = 1 N ( J x σ j x σ j + 1 x + J y σ j y σ j + 1 y + J z σ j z σ j + 1 z + h σ j z ) {\displaystyle {\hat {H}}=-{\frac {1}{2}}\sum _{j=1}^{N}(J_{x}\sigma _{j}^{x}\sigma _{j+1}^{x}+J_{y}\sigma _{j}^{y}\sigma _{j+1}^{y}+J_{z}\sigma _{j}^{z}\sigma _{j+1}^{z}+h\sigma _{j}^{z})}

where the  h {\displaystyle h} on the right-hand side indicates the external magnetic field, with periodic boundary conditions. The objective is to determine the spectrum of the Hamiltonian, from which the partition function can be calculated and the thermodynamics of the system can be studied.

It is common to name the model depending on the values of  J x {\displaystyle J_{x}} ,  J y {\displaystyle J_{y}} and  J z {\displaystyle J_{z}} : if  J x ≠ J y ≠ J z {\displaystyle J_{x}\neq J_{y}\neq J_{z}} , the model is called the Heisenberg XYZ model; in the case of  J = J x = J y ≠ J z = Δ {\displaystyle J=J_{x}=J_{y}\neq J_{z}=\Delta } , it is the Heisenberg XXZ model; if  J x = J y = J z = J {\displaystyle J_{x}=J_{y}=J_{z}=J} , it is the Heisenberg XXX model. The spin 1/2 Heisenberg model in one dimension may be solved exactly using the Bethe ansatz.[1] In the algebraic formulation, these are related to particular quantum affine algebras and elliptic quantum groups in the XXZ and XYZ cases respectively.[2] Other approaches do so without Bethe ansatz.[3]

### XXX model

[edit]

The physics of the Heisenberg XXX model strongly depends on the sign of the coupling constant  J {\displaystyle J} and the dimension of the space. For positive  J {\displaystyle J} the ground state is always ferromagnetic. At negative  J {\displaystyle J} the ground state is antiferromagnetic in two and three dimensions.[4] In one dimension the nature of correlations in the antiferromagnetic Heisenberg model depends on the spin of the magnetic dipoles. If the spin is integer then only short-range order is present. A system of half-integer spins exhibits quasi-long range order.

A simplified version of Heisenberg model is the one-dimensional Ising model, where the transverse magnetic field is in the x-direction, and the interaction is only in the z-direction:

H ^ = − J ∑ j = 1 N σ j z σ j + 1 z − g J ∑ j = 1 N σ j x {\displaystyle {\hat {H}}=-J\sum _{j=1}^{N}\sigma _{j}^{z}\sigma _{j+1}^{z}-gJ\sum _{j=1}^{N}\sigma _{j}^{x}} .

At small g and large g, the ground state degeneracy is different, which implies that there must be a quantum phase transition in between. It can be solved exactly for the critical point using the duality analysis.[5] The duality transition of the Pauli matrices is  σ i z = ∏ j ≤ i S j x {\textstyle \sigma _{i}^{z}=\prod _{j\leq i}S_{j}^{x}} and  σ i x = S i z S i + 1 z {\displaystyle \sigma _{i}^{x}=S_{i}^{z}S_{i+1}^{z}} , where  S x {\displaystyle S^{x}} and  S z {\displaystyle S^{z}} are also Pauli matrices which obey the Pauli matrix algebra. Under periodic boundary conditions, the transformed Hamiltonian can be shown is of a very similar form:

H ^ = − g J ∑ j = 1 N S j z S j + 1 z − J ∑ j = 1 N S j x {\displaystyle {\hat {H}}=-gJ\sum _{j=1}^{N}S_{j}^{z}S_{j+1}^{z}-J\sum _{j=1}^{N}S_{j}^{x}}

but for the  g {\displaystyle g} attached to the spin interaction term. Assuming that there's only one critical point, we can conclude that the phase transition happens at  g = 1 {\displaystyle g=1} .

## Solution by Bethe ansatz

[edit]

Main article: Bethe ansatz

### XXX1/2 model

[edit]

Following the approach of Ludwig Faddeev (1996), the spectrum of the Hamiltonian for the XXX model  H = 1 4 ∑ α , n ( σ n α σ n + 1 α − 1 ) {\displaystyle H={\frac {1}{4}}\sum _{\alpha ,n}(\sigma _{n}^{\alpha }\sigma _{n+1}^{\alpha }-1)} can be determined by the Bethe ansatz. In this context, for an appropriately defined family of operators  B ( λ ) {\displaystyle B(\lambda )} dependent on a spectral parameter  λ ∈ C {\displaystyle \lambda \in \mathbb {C} } acting on the total Hilbert space  H = ⨂ n = 1 N h n {\displaystyle {\mathcal {H}}=\bigotimes _{n=1}^{N}h_{n}} with each  h n ≅ C 2 {\displaystyle h_{n}\cong \mathbb {C} ^{2}} , a Bethe vector is a vector of the form  Φ ( λ 1 , ⋯ , λ m ) = B ( λ 1 ) ⋯ B ( λ m ) v 0 {\displaystyle \Phi (\lambda _{1},\cdots ,\lambda _{m})=B(\lambda _{1})\cdots B(\lambda _{m})v_{0}} where  v 0 = ⨂ n = 1 N | ↑ ⟩ {\displaystyle v_{0}=\bigotimes _{n=1}^{N}|\uparrow \,\rangle } . If the  λ k {\displaystyle \lambda _{k}} satisfy the Bethe equation ( λ k + i / 2 λ k − i / 2 ) N = ∏ j ≠ k λ k − λ j + i λ k − λ j − i , {\displaystyle \left({\frac {\lambda _{k}+i/2}{\lambda _{k}-i/2}}\right)^{N}=\prod _{j\neq k}{\frac {\lambda _{k}-\lambda _{j}+i}{\lambda _{k}-\lambda _{j}-i}},} then the Bethe vector is an eigenvector of  H {\displaystyle H} with eigenvalue  − ∑ k 1 2 1 λ k 2 + 1 / 4 {\displaystyle -\sum _{k}{\frac {1}{2}}{\frac {1}{\lambda _{k}^{2}+1/4}}} .

The family  B ( λ ) {\displaystyle B(\lambda )} as well as three other families come from a transfer matrix T ( λ ) {\displaystyle T(\lambda )} (in turn defined using a Lax matrix), which acts on  H {\displaystyle {\mathcal {H}}} along with an auxiliary space  h a ≅ C 2 {\displaystyle h_{a}\cong \mathbb {C} ^{2}} , and can be written as a  2 × 2 {\displaystyle 2\times 2} block matrix with entries in  E n d ( H ) {\displaystyle \mathrm {End} ({\mathcal {H}})} ,  T ( λ ) = ( A ( λ ) B ( λ ) C ( λ ) D ( λ ) ) , {\displaystyle T(\lambda )={\begin{pmatrix}A(\lambda )&B(\lambda )\\\C(\lambda )&D(\lambda )\end{pmatrix}},} which satisfies fundamental commutation relations (FCRs) similar in form to the Yang–Baxter equation used to derive the Bethe equations. The FCRs also show there is a large commuting subalgebra given by the generating function F ( λ ) = t r a ( T ( λ ) ) = A ( λ ) + D ( λ ) {\displaystyle F(\lambda )=\mathrm {tr} _{a}(T(\lambda ))=A(\lambda )+D(\lambda )} , as  [ F ( λ ) , F ( μ ) ] = 0 {\displaystyle [F(\lambda ),F(\mu )]=0} , so when  F ( λ ) {\displaystyle F(\lambda )} is written as a polynomial in  λ {\displaystyle \lambda } , the coefficients all commute, spanning a commutative subalgebra which  H {\displaystyle H} is an element of. The Bethe vectors are in fact simultaneous eigenvectors for the whole subalgebra.

### XXXs model

[edit]

For higher spins, say spin  s {\displaystyle s} , replace  σ α {\displaystyle \sigma ^{\alpha }} with  S α {\displaystyle S^{\alpha }} coming from the Lie algebra representation of the Lie algebra s l ( 2 , C ) {\displaystyle {\mathfrak {sl}}(2,\mathbb {C} )} , of dimension  2 s + 1 {\displaystyle 2s+1} . The XXXs Hamiltonian  H = ∑ α , n ( S n α S n + 1 α − ( S n α S n + 1 α ) 2 ) {\displaystyle H=\sum _{\alpha ,n}(S_{n}^{\alpha }S_{n+1}^{\alpha }-(S_{n}^{\alpha }S_{n+1}^{\alpha })^{2})} is solvable by Bethe ansatz with Bethe equations  ( λ k + i s λ k − i s ) N = ∏ j ≠ k λ k − λ j + i λ k − λ j − i . {\displaystyle \left({\frac {\lambda _{k}+is}{\lambda _{k}-is}}\right)^{N}=\prod _{j\neq k}{\frac {\lambda _{k}-\lambda _{j}+i}{\lambda _{k}-\lambda _{j}-i}}.}

### XXZs model

[edit]

For spin  s {\displaystyle s} and a parameter  γ {\displaystyle \gamma } for the deformation from the XXX model, the BAE (Bethe ansatz equation) is  ( sinh ⁡ ( λ k + i s γ ) sinh ⁡ ( λ k − i s γ ) ) N = ∏ j ≠ k sinh ⁡ ( λ k − λ j + i γ ) sinh ⁡ ( λ k − λ j − i γ ) . {\displaystyle \left({\frac {\sinh(\lambda _{k}+is\gamma )}{\sinh(\lambda _{k}-is\gamma )}}\right)^{N}=\prod _{j\neq k}{\frac {\sinh(\lambda _{k}-\lambda _{j}+i\gamma )}{\sinh(\lambda _{k}-\lambda _{j}-i\gamma )}}.} Notably, for  s = 1 2 {\displaystyle s={\frac {1}{2}}} these are precisely the BAEs for the six-vertex model, after identifying  γ = 2 η {\displaystyle \gamma =2\eta } , where  η {\displaystyle \eta } is the anisotropy parameter of the six-vertex model.[6][7] This was originally thought to be coincidental until Baxter showed the XXZ Hamiltonian was contained in the algebra generated by the transfer matrix  T ( ν ) {\displaystyle T(\nu )} ,[8] given exactly by  H X X Z 1 / 2 = − i sin ⁡ 2 η d d ν log ⁡ T ( ν ) | ν = − i η − 1 2 cos ⁡ 2 η 1 ⊗ N . {\displaystyle H_{XXZ_{1/2}}=-i\sin 2\eta {\frac {d}{d\nu }}\log T(\nu ){\Big |}_{\nu =-i\eta }-{\frac {1}{2}}\cos 2\eta 1^{\otimes N}.}

## Applications

[edit]

* Another important object is entanglement entropy. One way to describe it is to subdivide the unique ground state into a block (several sequential spins) and the environment (the rest of the ground state). The entropy of the block can be considered as entanglement entropy. At zero temperature in the critical region (thermodynamic limit) it scales logarithmically with the size of the block. As the temperature increases the logarithmic dependence changes into a linear function.[9] For large temperatures linear dependence follows from the second law of thermodynamics.
* The Heisenberg model provides an important and tractable theoretical example for applying density matrix renormalisation.
* The six-vertex model can be solved using the algebraic Bethe ansatz for the Heisenberg spin chain (Baxter 1982).
* The half-filled Hubbard model in the limit of strong repulsive interactions can be mapped onto a Heisenberg model with  J < 0 {\displaystyle J<0} representing the strength of the superexchange interaction.
* Limits of the model as the lattice spacing is sent to zero (and various limits are taken for variables appearing in the theory) describes integrable field theories, both non-relativistic such as the nonlinear Schrödinger equation, and relativistic, such as the  S 2 {\displaystyle S^{2}} sigma model, the  S 3 {\displaystyle S^{3}} sigma model (which is also a principal chiral model) and the sine-Gordon model.
* Calculating certain correlation functions in the planar or large  N {\displaystyle N} limit of N = 4 supersymmetric Yang–Mills theory[10]

## Extended symmetry

[edit]

The integrability is underpinned by the existence of large symmetry algebras for the different models. For the XXX case this is the Yangian Y ( s l 2 ) {\displaystyle Y({\mathfrak {sl}}_{2})} , while in the XXZ case this is the quantum group s l q ( 2 ) ^ {\displaystyle {\hat {{\mathfrak {sl}}_{q}(2)}}} , the q-deformation of the affine Lie algebra of  s l 2 ^ {\displaystyle {\hat {{\mathfrak {sl}}_{2}}}} , as explained in the notes by Faddeev (1996).

These appear through the transfer matrix, and the condition that the Bethe vectors are generated from a state  Ω {\displaystyle \Omega } satisfying  C ( λ ) ⋅ Ω = 0 {\displaystyle C(\lambda )\cdot \Omega =0} corresponds to the solutions being part of a highest-weight representation of the extended symmetry algebras.

## See also

[edit]

* Classical Heisenberg model
* DMRG of the Heisenberg model
* Quantum rotor model
* t-J model
* J1 J2 model
* Majumdar–Ghosh model
* AKLT model
* Multipolar exchange interaction

## References

[edit]

* R.J. Baxter, Exactly solved models in statistical mechanics, London, Academic Press, 1982
* Heisenberg, W. (1 September 1928). "Zur Theorie des Ferromagnetismus" [On the theory of ferromagnetism]. Zeitschrift für Physik (in German). 49 (9): 619–636\. Bibcode:1928ZPhy...49..619H. doi:10.1007/BF01328601. S2CID 122524239.
* Bethe, H. (1 March 1931). "Zur Theorie der Metalle" [On the theory of metals]. Zeitschrift für Physik (in German). 71 (3): 205–226\. Bibcode:1931ZPhy...71..205B. doi:10.1007/BF01341708. S2CID 124225487.

## Notes

[edit]

1. ^ Bonechi, F; Celeghini, E; Giachetti, R; Sorace, E; Tarlini, M (7 August 1992). "Heisenberg XXZ model and quantum Galilei group". Journal of Physics A: Mathematical and General. 25 (15): L939 – L943. arXiv:hep-th/9204054. Bibcode:1992JPhA...25L.939B. doi:10.1088/0305-4470/25/15/007. S2CID 119046025.
2. ^ Faddeev, L. D. (26 May 1996). "How Algebraic Bethe Ansatz works for integrable model". arXiv:hep-th/9605187v1.
3. ^ Rojas, Onofre; Souza, S.M. de; Corrêa Silva, E.V.; Thomaz, M.T. (December 2001). "Thermodynamics of the limiting cases of the XXZ model without Bethe ansatz". Brazilian Journal of Physics. 31 (4): 577–582\. Bibcode:2001BrJPh..31..577R. doi:10.1590/s0103-97332001000400008.
4. ^ Tom Kennedy; Bruno Nachtergaele. "The Heisenberg Model - a Bibliography". Retrieved 6 Jun 2019.
5. ^ Fisher, Matthew P. A. (2004). "Duality in low dimensional quantum field theories". Strong interactions in low dimensions. Physics and Chemistry of Materials with Low-Dimens. Vol. 25\. pp. 419–438\. doi:10.1007/978-1-4020-3463-3_13. ISBN 978-1-4020-1798-8.
6. ^ Lieb, Elliott H. (24 April 1967). "Exact Solution of the Problem of the Entropy of Two-Dimensional Ice". Physical Review Letters. 18 (17): 692–694\. Bibcode:1967PhRvL..18..692L. doi:10.1103/PhysRevLett.18.692.
7. ^ Dorey, Patrick; Dunning, Clare; Tateo, Roberto (10 August 2007). "The ODE/IM correspondence". Journal of Physics A: Mathematical and Theoretical. 40 (32): R205 – R283. arXiv:hep-th/0703066. doi:10.1088/1751-8113/40/32/R01. ISSN 1751-8113. S2CID 14281617.
8. ^ Baxter, Rodney J (1 April 1972). "One-dimensional anisotropic Heisenberg chain". Annals of Physics. 70 (2): 323–337\. Bibcode:1972AnPhy..70..323B. doi:10.1016/0003-4916(72)90270-9. ISSN 0003-4916.
9. ^ Korepin, V. E. (5 March 2004). "Universality of Entropy Scaling in One Dimensional Gapless Models". Physical Review Letters. 92 (9): 096402. arXiv:cond-mat/0311056. Bibcode:2004PhRvL..92i6402K. doi:10.1103/PhysRevLett.92.096402. PMID 15089496. S2CID 20620724.
10. ^ Beisert, Niklas (1 December 2004). "The dilatation operator of N=4 super Yang–Mills theory and integrability". Physics Reports. 405 (1): 1–202\. arXiv:hep-th/0407277. Bibcode:2004PhR...405....1B. doi:10.1016/j.physrep.2004.09.007. S2CID 118949332.

* v
* t
* e

Statistical mechanics

Theory|

* Principle of maximum entropy
* ergodic theory

Statistical thermodynamics|

* Ensembles
* partition functions
* equations of state
* thermodynamic potential:
* U
* H
* F
* G
* Maxwell relations

Models|

* Ferromagnetism models
* Ising
* Potts
* Heisenberg
* percolation
* Particles with force field
* depletion force
* Lennard-Jones potential

Mathematical approaches|

* Boltzmann equation
* H-theorem
* Vlasov equation
* BBGKY hierarchy
* stochastic process
* mean-field theory and conformal field theory

Critical phenomena|

* Phase transition
* Critical exponents
* correlation length
* size scaling

Entropy|

* Boltzmann
* Shannon
* Tsallis
* Rényi
* von Neumann

Applications|

* Statistical field theory
* elementary particle
* superfluidity
* Condensed matter physics
* Complex system
* chaos
* information theory
* Boltzmann machine

* v
* t
* e

Integrable systems

Geometric integrability|

* Frobenius integrability
* Foliation
* Liouville integrability

In classical mechanics|

| Examples|

* Harmonic oscillator
* Central force systems
* Kepler system
* Two body problem
* Integrable tops
* Euler
* Kovalevskaya
* Lagrange
* Garnier integrable system
* Hitchin system

|

Theory|

* Liouville–Arnold theorem
* Action-angle variables
* Superintegrable Hamiltonian system

In quantum mechanics|

* Quantum harmonic oscillator
* Hydrogen atom
* Pöschl–Teller potential

Integrable PDEs/Classical integrable field theories|

| Examples|

* KdV equation
* KdV hierarchy
* Sine-Gordon equation
* Nonlinear Schrödinger equation
* Gross–Neveu model
* Thirring model
* Kadomtsev–Petviashvili equation

|

Theory|

* Bäcklund transformation
* Lax pairs
* Infinitely many integrals of motion
* Soliton solutions
* Topological soliton
* Inverse scattering transform

ASDYM as a master theory|

* Anti-self-dual Yang–Mills equations
* Twistor correspondence
* Ward conjecture

Integrable Quantum Field theories|

* Quantum Sine-Gordon
* Quantum KdV
* Quantum Liouville
* Thirring model
* Toda field theory
* Principal chiral model

| Master theories|

* Four-dimensional Chern–Simons theory (Lagrangian)
* Affine Gaudin models (Hamiltonian)
* Six-dimensional holomorphic Chern–Simons theory

|

Exactly solvable statistical lattice models|

* Ising model in one- and two-dimensions
* Square ice model
* Eight-vertex model
* Hard hexagon model
* Chiral Potts model

Exactly solvable quantum spin chains|

| Examples|

* Quantum Heisenberg model
* Gaudin model

|

Theory|

* Algebraic Bethe ansatz
* Quantum inverse scattering method
* Yang–Baxter equation

Contributors|

| Classical mechanics and geometry|

* Vladimir Arnold
* Leonhard Euler
* Ferdinand Georg Frobenius
* Nigel Hitchin
* Sofia Kovalevskaya
* Joseph-Louis Lagrange
* Joseph Liouville
* Siméon Denis Poisson

|

PDEs|

* Clifford S. Gardner
* John M. Greene
* Martin David Kruskal
* Peter Lax
* Robert Miura

IQFTs|

* Alexander Zamolodchikov
* Alexei Zamolodchikov

Classical and quantum statistical lattices|

* Rodney Baxter
* Ludvig Faddeev
* Elliott H. Lieb
* Yang Chen-Ning

*[v]: View this template
*[t]: Discuss this template
*[e]: Edit this template