# Heisenberg model (quantum)

The 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. In the prototypical Ising model, defined on a d-dimensional lattice, at each lattice site, a spin ${\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

For quantum mechanical reasons (see exchange interaction or the subchapter "quantum-mechanical origin of magnetism" in the article on 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) the Hamiltonian can be written in the form

${\displaystyle {\hat {H}}=-J\sum _{j=1}^{N}\sigma _{j}\sigma _{j+1}-h\sum _{j=1}^{N}\sigma _{j}}$

where ${\displaystyle J}$  is the coupling constant for a 1-dimensional model consisting of N dipoles, represented by classical vectors (or "spins") σj, subject to the periodic boundary condition ${\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 ${\displaystyle (\mathbb {C} ^{2})^{\otimes N}}$ , of dimension ${\displaystyle 2^{N}}$ . To define it, recall the Pauli spin-1/2 matrices

${\displaystyle \sigma ^{x}={\begin{pmatrix}0&1\\1&0\end{pmatrix}}}$
${\displaystyle \sigma ^{y}={\begin{pmatrix}0&-i\\i&0\end{pmatrix}}}$
${\displaystyle \sigma ^{z}={\begin{pmatrix}1&0\\0&-1\end{pmatrix}}}$

and for ${\displaystyle 1\leq j\leq N}$  and ${\displaystyle a\in \{x,y,z\}}$  denote ${\displaystyle \sigma _{j}^{a}=I^{\otimes j-1}\otimes \sigma ^{a}\otimes I^{\otimes N-j}}$ , where ${\displaystyle I}$  is the ${\displaystyle 2\times 2}$  identity matrix. In 3-dimensions, given a choice of real-valued coupling constants ${\displaystyle J_{x},J_{y},}$  and ${\displaystyle J_{z}}$ , the Hamiltonian is given by

${\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 ${\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, from which the thermodynamics of the system can be studied.

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:

${\displaystyle {\hat {H}}=-J_{z}\sum _{j=1}^{N}\sigma _{j}^{z}\sigma _{j+1}^{z}-gJ_{z}\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.[1] The duality transition of the Pauli matrices is ${\displaystyle \sigma _{i}^{z}=\prod _{j\leq i}S_{j}^{x}}$  and ${\displaystyle \sigma _{i}^{x}=S_{i}^{z}S_{i+1}^{z}}$ , where ${\displaystyle S^{x}}$  and ${\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:

${\displaystyle {\hat {H}}=-gJ_{z}\sum _{j=1}^{N}S_{j}^{z}S_{j+1}^{z}-gJ_{z}\sum _{j=1}^{N}S_{j}^{x}}$

but for the ${\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 ${\displaystyle g=1}$ .

It is common to name the model depending on the values of ${\displaystyle J_{x}}$ , ${\displaystyle J_{y}}$  and ${\displaystyle J_{z}}$ : if ${\displaystyle J_{x}\neq J_{y}\neq J_{z}}$ , the model is called the Heisenberg XYZ model; in the case of ${\displaystyle J=J_{x}=J_{y}\neq J_{z}=\Delta }$ , it is the Heisenberg XXZ model; if ${\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.[2] In the algebraic formulation, these are related to particular Quantum affine algebras and Elliptic Quantum Group in the XXZ and XYZ cases respectively.[3] Other approaches do so without Bethe ansatz.[4]

The physics of the Heisenberg model strongly depends on the sign of the coupling constant ${\displaystyle J}$  and the dimension of the space. For positive ${\displaystyle J}$  the ground state is always ferromagnetic. At negative ${\displaystyle J}$  the ground state is antiferromagnetic in two and three dimensions.[5] 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.

## Applications

• 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. For large temperatures linear dependence follows from the second law of thermodynamics.
• The six-vertex model can be solved using the Algebraic Bethe Ansatz for the Heisenberg Spin Chain (see Baxter, "Exactly Solved Models in Statistical Mechanics").
• The half-filled Hubbard model in the limit of strong repulsive interactions can be mapped onto a Heisenberg model with ${\displaystyle J<0}$  representing the strength of the superexchange interaction.

## References

• R.J. Baxter, Exactly solved models in statistical mechanics, London, Academic Press, 1982
• H. Bethe, Zur Theorie der Metalle, Zeitschrift für Physik A, 1931 doi:10.1007/BF01341708

## Notes

1. ^ Fisher, Matthew P. A. "Duality in low dimensional quantum field theories". Physics and Chemistry of Materials with Low-Dimens (PDF). Dordrecht: Springer Netherlands. pp. 419–438. doi:10.1007/978-1-4020-3463-3_13. ISBN 978-1-4020-1798-8.
2. ^
3. ^ Faddeev, L. D. (26 May 1996). "How Algebraic Bethe Ansatz works for integrable model". arXiv:hep-th/9605187v1.
4. ^ Rojas, Onofre; Souza, S.M. de; Corrêa Silva, E.V.; Thomaz, M.T. (2001). "Thermodynamics of the limiting cases of the XXZ model without Bethe ansatz". Brazilian Journal of Physics. FapUNIFESP (SciELO). 31 (4): 577–582. doi:10.1590/s0103-97332001000400008. ISSN 0103-9733.
5. ^ Tom Kennedy; Bruno Nachtergaele. "The Heisenberg Model - a Bibliography". Retrieved 6 Jun 2019.