Skip to content

Heisenberg Chain

The spin-1/21/2 Heisenberg chain is a one-dimensional array of two-level quantum spins coupled by rotationally invariant nearest-neighbor exchange, with an exact ferromagnetic ground multiplet for J<0J<0 and a Bethe-ansatz-integrable, gapless spinon liquid for J>0J>0.

This dossier fixes the chain’s model data, normalization, exact-status claims, limiting cases, observables, and validation targets. The Heisenberg Model teaching article owns the broader treatment of arbitrary spin, graphs, microscopic exchange, spin waves, dimensionality, and physical interpretation. Exact Solutions Preview owns the general Bethe-ansatz and integrability framework.

Unless a variant is stated explicitly, this dossier uses the following chain.

FieldBaseline choice
sitesLL spin-1/21/2 degrees of freedom
local spaceHj≃C2\mathcal H_j\simeq\mathbb C^2
spin conventiondimensionless sj=σj/2\mathbf s_j=\boldsymbol\sigma_j/2
geometrya uniform one-dimensional nearest-neighbor chain
exchangeisotropic xx, yy, and zz coupling
couplingreal JJ with units of energy
signJ>0J>0 antiferromagnetic; J<0J<0 ferromagnetic
boundariesopen or periodic, stated explicitly
fieldabsent
energy offsetnone
anisotropy, disorder, and longer rangeabsent
thermodynamic antiferromagnetic referenceeven periodic rings followed by L→∞L\to\infty

This baseline is narrower than “the Heisenberg model” as a family. Changing the onsite spin, graph, dimension, bond signs, unit cell, anisotropy, or field can change the low-energy phase and exact solution status. Those changes are variants, not implicit parts of one numerical contract.

The many-body Hilbert space is

HL=⨂j=1LC2,dim⁡HL=2L.\mathcal H_L = \bigotimes_{j=1}^{L} \mathbb C^2, \qquad \dim\mathcal H_L = 2^L.

Dimensionless local spins obey

[siα,sjβ]=iδijϵαβγsiγ,si2=34,\begin{aligned} [s_i^\alpha,s_j^\beta] &= i\delta_{ij} \epsilon_{\alpha\beta\gamma} s_i^\gamma, \\ \mathbf s_i^2 &= \frac34, \end{aligned}

with

si=12σi.\mathbf s_i = \frac12 \boldsymbol\sigma_i.

The computational basis diagonalizes every sjzs_j^z:

sz∣↑⟩=12∣↑⟩,sz∣↓⟩=−12∣↓⟩.s^z\lvert\uparrow\rangle = \frac12\lvert\uparrow\rangle, \qquad s^z\lvert\downarrow\rangle = -\frac12\lvert\downarrow\rangle.

A basis state can be encoded by a bit string, but the Hamiltonian is not diagonal in that basis. On one bond,

si⋅sj=sizsjz+12(si+sj−+si−sj+).\mathbf s_i\cdot\mathbf s_j = s_i^zs_j^z + \frac12 \left( s_i^+s_j^- + s_i^-s_j^+ \right).

The first term measures parallel or antiparallel zz projections. The remaining terms exchange an up-down pair and are responsible for quantum fluctuations in the antiferromagnet.

For an open chain,

HO=J∑j=1L−1sj⋅sj+1.H_{\mathrm O} = J \sum_{j=1}^{L-1} \mathbf s_j\cdot\mathbf s_{j+1}.

For a periodic ring,

HP=J∑j=1Lsj⋅sj+1,sL+1=s1.H_{\mathrm P} = J \sum_{j=1}^{L} \mathbf s_j\cdot\mathbf s_{j+1}, \qquad \mathbf s_{L+1} = \mathbf s_1.

The open chain has L−1L-1 undirected bonds; an ordinary ring has LL. Each physical bond appears once. The L=2L=2 periodic shorthand is exceptional because the terms labeled j=1j=1 and j=2j=2 refer to the same undirected pair. A two-site ring contract must say whether it contains one bond or a deliberately doubled multigraph edge.

Because sj=σj/2\mathbf s_j=\boldsymbol\sigma_j/2,

H=J4∑⟨j,k⟩σj⋅σk.H = \frac J4 \sum_{\langle j,k\rangle} \boldsymbol\sigma_j \mathbin{\cdot} \boldsymbol\sigma_k.

A spectrum four times too large usually comes from inserting Pauli matrices directly into the spin-operator Hamiltonian without the factor 1/41/4.

Physical angular momenta are

Sj=ℏsj.\mathbf S_j = \hbar\mathbf s_j.

In that convention the same exchange is

Jsi⋅sj=Jℏ2Si⋅Sj.J\mathbf s_i\cdot\mathbf s_j = \frac{J}{\hbar^2} \mathbf S_i\cdot\mathbf S_j.

Writing JSi⋅SjJ\mathbf S_i\cdot\mathbf S_j instead assigns different units to the symbol JJ.

For two spins,

sij=si+sj,\mathbf s_{ij} = \mathbf s_i+\mathbf s_j,

and

si⋅sj=12(sij2−si2−sj2).\mathbf s_i\cdot\mathbf s_j = \frac12 \left( \mathbf s_{ij}^2 - \mathbf s_i^2 - \mathbf s_j^2 \right).

Two spin-1/21/2 spaces decompose into a singlet and a triplet:

12⊗12=0⊕1.\frac12\otimes\frac12 = 0\oplus1.

The exchange eigenvalues are

bond spinsi⋅sjE/J0−3/4−3/41+1/4+1/4\begin{array}{c|c|c} \text{bond spin} & \mathbf s_i\cdot\mathbf s_j & E/J \\ \hline 0 & -3/4 & -3/4 \\ 1 & +1/4 & +1/4 \end{array}

so J>0J>0 favors the singlet and J<0J<0 favors the triplet. Equivalently, the swap operator on two spin-1/21/2 sites satisfies

Pij=2si⋅sj+12I.\mathsf P_{ij} = 2\mathbf s_i\cdot\mathbf s_j + \frac12\mathbb I.

The bond projectors are

ΠS=14−si⋅sj,ΠT=34+si⋅sj.\begin{aligned} \Pi_{S} &= \frac14 - \mathbf s_i\cdot\mathbf s_j, \\ \Pi_{T} &= \frac34 + \mathbf s_i\cdot\mathbf s_j. \end{aligned}

These formulas are exact local checks and do not imply that every bond in an extended antiferromagnet can simultaneously occupy a singlet.

Define total spin

Stot=∑j=1Lsj.\mathbf S_{\mathrm{tot}} = \sum_{j=1}^{L} \mathbf s_j.

Every exchange bond is a scalar under a common spin rotation, hence

[H,Stotα]=0,[H,Stot2]=0.\begin{aligned} [H,S_{\mathrm{tot}}^\alpha] &= 0, \\ [H,\mathbf S_{\mathrm{tot}}^2] &= 0. \end{aligned}

The internal continuous symmetry is SU(2)SU(2) at the Hilbert-space level, with physical rotations represented projectively through the corresponding SO(3)SO(3) action. Energy eigenstates can be organized into total-spin multiplets labeled by

Stot=0,12,1,…,L2.S_{\mathrm{tot}} = 0,\frac12,1,\ldots,\frac L2.

A spin-SS multiplet contains 2S+12S+1 states with

m=−S,−S+1,…,S.m = -S,-S+1,\ldots,S.

Numerically equal energies across all these mm values are a stringent symmetry check.

The conserved component

Stotz=∑jsjzS_{\mathrm{tot}}^z = \sum_j s_j^z

allows a product-basis block decomposition. If MM sites are down,

Stotz=L2−M,S_{\mathrm{tot}}^z = \frac L2-M,

and the block dimension is

dim⁡HM=(LM).\dim\mathcal H_M = \binom LM.

An StotzS_{\mathrm{tot}}^z block does not isolate one total-spin irreducible representation; it generally contains one state from every compatible multiplet. Recovering the complete spectrum requires all magnetization blocks or an explicit multiplet reconstruction.

The uniform periodic chain has translation and reflection symmetry. The open chain retains reflection about its midpoint but not one-site translation symmetry.

With no magnetic field, physical time reversal sends

sj↦−sj\mathbf s_j \mapsto -\mathbf s_j

and leaves every exchange bond invariant. For LL spin-1/21/2 sites, the many-body time-reversal operator obeys

Θ2=(−1)L.\Theta^2 = (-1)^L.

Odd chains therefore carry Kramers degeneracy in the absence of time-reversal-breaking terms. This antiunitary statement is distinct from the unitary SU(2)SU(2) multiplet degeneracy, although both may act on the same levels.

RegimeExact or controlled pictureMain caution
J=0J=0all 2L2^L states are degenerateno exchange scale remains
L=2L=2 with one bondexact singlet–triplet dimera periodic shorthand may double the bond
J<0J<0maximal-spin ground multipletlow-energy magnons are quadratic
J>0J>0, even ringsinglet ground state and gapless spinon continuum as L→∞L\to\inftyfinite-size gaps vanish with logarithmic corrections
J>0J>0, odd ringunavoidable antiferromagnetic frustrationground state cannot be a singlet for odd LL
open chainL−1L-1 bonds and boundary correctionsmomentum is not a good quantum number
strong uniform fieldpolarized state above the saturation fieldfield reduces SU(2)SU(2) to U(1)U(1)
higher onsite spininteger and half-integer chains can differ qualitativelynot the spin-1/21/2 baseline

For J<0J<0 on a connected chain, every bond can attain its triplet value simultaneously. The ground states form the maximal-spin multiplet

Stot=L2.S_{\mathrm{tot}} = \frac L2.

With NbN_b bonds, its energy is

EF=JNb4.E_F = \frac{J N_b}{4}.

The fully polarized product state is one member; global rotations generate the remaining LL states in the (L+1)(L+1)-dimensional multiplet.

For J>0J>0, a singlet minimizes one isolated bond, but adjacent bonds share sites. An extended chain cannot be a product of singlets on every nearest-neighbor bond while preserving one-site translation symmetry. The classical Néel product state also fails to be an eigenstate because the transverse exchange terms move spin deviations.

For an even bipartite chain, the ground state is a total-spin singlet. It has strong staggered correlations but no nonzero staggered magnetization in the translation-invariant one-dimensional ground state. The thermodynamic chain is critical rather than Néel ordered.

An odd antiferromagnetic ring cannot alternate up and down consistently around the boundary. Its graph is nonbipartite, its minimum total spin is half-integer, and its low-energy quantum numbers differ from an even ring. Odd and even sequences should not be mixed casually in finite-size extrapolations.

The chain contains several exact structures, but they have different scopes.

Quantity or regimeStatus
two-site spectrumexact by angular-momentum addition
ferromagnetic ground manifoldexact for uniform J<0J<0
ferromagnetic one-magnon sectorexact plane-wave diagonalization
uniform spin-1/21/2 periodic chainBethe-ansatz integrable
thermodynamic antiferromagnetic energy densityexact
antiferromagnetic low-energy velocity and continuum edgesexact
finite-temperature bulk thermodynamicsexact integral-equation framework
generic spin correlationsexact representations exist, but evaluation and asymptotics are nontrivial
arbitrary higher-spin chainnot automatically solved by the spin-1/21/2 Bethe equations
higher dimensions or frustrated graphsnot generically Bethe-ansatz solvable
generic anisotropy, disorder, or longer-range exchangeintegrability must be reassessed

One rapidity convention for the periodic spin-1/21/2 chain uses roots λj\lambda_j satisfying

(λj+i/2λj−i/2)L=∏k=1k≠jMλj−λk+iλj−λk−i.\left( \frac{ \lambda_j+i/2 }{ \lambda_j-i/2 } \right)^L = \prod_{\substack{k=1\\k\ne j}}^{M} \frac{ \lambda_j-\lambda_k+i }{ \lambda_j-\lambda_k-i }.

For the unshifted Hamiltonian in this dossier, the energy is

E=JL4−J2∑j=1M1λj2+1/4.E = \frac{JL}{4} - \frac J2 \sum_{j=1}^{M} \frac{1}{ \lambda_j^2+1/4 }.

These equations pin the convention and the meaning of exactness; they do not replace the root-classification, completeness, thermodynamic-limit, or correlation-function machinery. The XXZ Spin Chain and Exact Solutions Preview own those developments.

Let J<0J<0 and ∣F⟩\lvert F\rangle be fully polarized. A normalized localized spin deviation is

∣j⟩=sj−∣F⟩\lvert j\rangle = s_j^-\lvert F\rangle

for spin 1/21/2. Within the one-deviation subspace,

(H−EF)∣j⟩=J2(∣j−1⟩+∣j+1⟩−2∣j⟩).\left( H-E_F \right) \lvert j\rangle = \frac J2 \left( \lvert j-1\rangle + \lvert j+1\rangle - 2\lvert j\rangle \right).

Momentum states have exact excitation energy

εF(k)=∣J∣[1−cos⁡(ka)].\varepsilon_F(k) = \lvert J\rvert \left[ 1-\cos(ka) \right].

Near k=0k=0,

εF(k)≃∣J∣2(ka)2.\varepsilon_F(k) \simeq \frac{\lvert J\rvert}{2} (ka)^2.

The k=0k=0 state belongs to the degenerate ground multiplet. The first dispersing mode of a periodic ring has k=2π/(La)k=2\pi/(La) and an energy of order L−2L^{-2}.

For the even periodic spin-1/21/2 chain with J>0J>0,

lim⁡L→∞E0(L)LJ=14−ln⁡2≃−0.443147180559945.\lim_{L\to\infty} \frac{E_0(L)}{LJ} = \frac14-\ln2 \simeq -0.443147180559945.

The Néel product state has exchange expectation −J/4-J/4 per bond. The exact value is lower because transverse exchange produces correlated quantum fluctuations.

The thermodynamic antiferromagnetic chain is gapless. Its elementary fractional excitations are spinons with spin 1/21/2, while a local spin operator creates states with integer total spin and therefore accesses multi-spinon continua.

For physical wave number qq, the exact lower boundary of the two-spinon continuum is

εL(q)=πJ2∣sin⁡(qa)∣.\varepsilon_{\mathrm L}(q) = \frac{\pi J}{2} \left| \sin(qa) \right|.

The upper boundary of the two-spinon contribution is

εU(q)=πJ∣sin⁡qa2∣.\varepsilon_{\mathrm U}(q) = \pi J \left| \sin\frac{qa}{2} \right|.

The lower boundary is linear near a gapless wavevector, giving spin velocity

v=πJa2ℏ.v = \frac{\pi Ja}{2\hbar}.

The low-energy continuum theory has central charge

c=1.c = 1.

Marginally irrelevant interactions generate logarithmic corrections to simple scaling forms. Finite-size gaps, correlations, and entanglement should therefore not be fit blindly to pure power laws over short size ranges.

The most basic labels are

Stot2andStotz.\mathbf S_{\mathrm{tot}}^2 \quad\text{and}\quad S_{\mathrm{tot}}^z.

A uniform magnetization per site is

m=1L∑jsj.\mathbf m = \frac1L \sum_j\mathbf s_j.

On a bipartite even chain, a staggered operator is

mst=1L∑j(−1)jsj.\mathbf m_{\mathrm{st}} = \frac1L \sum_j (-1)^j \mathbf s_j.

In an SU(2)SU(2)-invariant finite-system state, ⟨mst⟩\langle\mathbf m_{\mathrm{st}}\rangle vanishes. Rotationally invariant diagnostics include ⟨mst2⟩\langle\mathbf m_{\mathrm{st}}^2\rangle and long-distance spin correlations.

Define

Cαβ(r)=⟨sjαsj+rβ⟩.C^{\alpha\beta}(r) = \left\langle s_j^\alpha s_{j+r}^\beta \right\rangle.

In an SU(2)SU(2)-invariant state,

⟨siαsjβ⟩=δαβ3⟨si⋅sj⟩.\left\langle s_i^\alpha s_j^\beta \right\rangle = \frac{ \delta_{\alpha\beta} }{3} \left\langle \mathbf s_i\cdot\mathbf s_j \right\rangle.

This tensor relation is a useful numerical symmetry test. It need not hold after choosing a symmetry-broken state, applying a field, or introducing anisotropy.

For a translation-invariant ring, one common static convention is

S(q)=1L∑j,ke−iqa(j−k)⟨sj⋅sk⟩.S(q) = \frac1L \sum_{j,k} e^{-iq a(j-k)} \left\langle \mathbf s_j\cdot\mathbf s_k \right\rangle.

Ferromagnetic correlations accumulate near q=0q=0; antiferromagnetic correlations accumulate near q=π/aq=\pi/a. The peak must be scaled with LL before it is interpreted as long-range order.

A dynamic structure factor additionally resolves energy transfer:

Sαβ(q,ω)=12πL∫−∞∞dt eiωt×∑j,ke−iqa(j−k)⟨sjα(t)skβ(0)⟩.\begin{aligned} S^{\alpha\beta}(q,\omega) ={}& \frac{1}{2\pi L} \int_{-\infty}^{\infty} \mathrm dt\, e^{i\omega t} \\ &\times \sum_{j,k} e^{-iqa(j-k)} \left\langle s_j^\alpha(t) s_k^\beta(0) \right\rangle. \end{aligned}

It distinguishes a sharp ferromagnetic magnon from the antiferromagnetic spinon continuum. Structure Factors owns Fourier, normalization, detailed-balance, and experimental conventions.

A finite-size gap must state the compared sectors:

ΔL=Etarget(L)−E0(L).\Delta_L = E_{\mathrm{target}}(L) - E_0(L).

The antiferromagnetic bulk gap vanishes, but finite rings have nonzero level spacings. The ferromagnetic ground multiplet is already degenerate, so the gap to the first state outside that multiplet is the meaningful dispersing scale.

For a contiguous interval, the ground-state entanglement entropy is a sensitive distinction between the gapped variants and the critical spin-1/21/2 antiferromagnet. At criticality its leading logarithm is governed by c=1c=1, with boundary-dependent coefficients and subleading corrections.

Rotational invariance organizes levels by total spin. The sign of JJ selects whether a bond prefers low or high combined spin, but that local preference does not by itself establish a thermodynamic phase.

All ferromagnetic bonds can minimize their energy in the same maximal-spin state. The resulting product-state representative is exact, and the quadratic magnon follows from the broken nonrelativistic spin-rotation generators. At any nonzero temperature, however, a short-range one-dimensional isotropic chain has no true ferromagnetic long-range order.

Antiferromagnetic bonds compete because each site is shared. The ground state lowers its energy through a coherent superposition of spin configurations rather than a classical alternating product. In one dimension the outcome is a critical singlet with algebraic correlations and fractional spinons, not a Néel-ordered state with a single sharp magnon branch.

SU(2)SU(2) symmetry and Bethe integrability are distinct structures. Many higher-dimensional or perturbed Heisenberg Hamiltonians remain rotationally invariant but lose the commuting transfer-matrix hierarchy that makes the uniform chain solvable.

The spin-1/21/2 antiferromagnetic chain is gapless, while the uniform spin-1 chain has a nonzero Haldane gap and boundary spin-1/21/2 modes under open boundaries. Higher-dimensional bipartite antiferromagnets can order at zero temperature. These are variants of the Heisenberg family, not contradictions of the baseline record.

For two spin-1/21/2 sites with one bond,

H2=Js1⋅s2.H_2 = J \mathbf s_1\cdot\mathbf s_2.

Let

S=s1+s2.\mathbf S = \mathbf s_1+\mathbf s_2.

Since

s1⋅s2=12(S2−32),\mathbf s_1\cdot\mathbf s_2 = \frac12 \left( \mathbf S^2 - \frac32 \right),

the total-spin singlet has

ES=−3J4,E_S = -\frac{3J}{4},

and the triplet has

ET=+J4.E_T = +\frac{J}{4}.

The singlet is

∣S⟩=∣↑↓⟩−∣↓↑⟩2,\lvert S\rangle = \frac{ \lvert\uparrow\downarrow\rangle - \lvert\downarrow\uparrow\rangle }{\sqrt2},

while the triplet is spanned by

∣T+⟩=∣↑↑⟩,∣T0⟩=∣↑↓⟩+∣↓↑⟩2,∣T−⟩=∣↓↓⟩.\begin{aligned} \lvert T_+\rangle &= \lvert\uparrow\uparrow\rangle, \\ \lvert T_0\rangle &= \frac{ \lvert\uparrow\downarrow\rangle + \lvert\downarrow\uparrow\rangle }{\sqrt2}, \\ \lvert T_-\rangle &= \lvert\downarrow\downarrow\rangle. \end{aligned}

The singlet–triplet separation is

ET−ES=J.E_T-E_S = J.

Thus antiferromagnetic exchange selects an entangled singlet, while ferromagnetic exchange selects the triplet multiplet. The trace vanishes:

ES+3ET=0,E_S+3E_T = 0,

as expected because every nonidentity Pauli string has zero trace.

MB-B002 fixes

H=J∑j=14sj⋅sj+1,s5=s1,H = J \sum_{j=1}^{4} \mathbf s_j\cdot\mathbf s_{j+1}, \qquad \mathbf s_5 = \mathbf s_1,

with J=1J=1, four undirected bonds, and no offset. The complete spectrum is

E/JE/JdegeneracySU(2)SU(2) content
−2-21one singlet
−1-13one triplet
007one singlet and two triplets
+1+15one spin-two multiplet

The degeneracies sum to

1+3+7+5=16=24.1+3+7+5 = 16 = 2^4.

The fixed-MM product-basis block dimensions are

1, 4, 6, 4, 1.1,\ 4,\ 6,\ 4,\ 1.

Basis-independent checks are

Tr⁡H=0,116Tr⁡H2=34J2,[H,Stotz]=0,[H,Stot2]=0.\begin{aligned} \operatorname{Tr}H &= 0, \\ \frac1{16} \operatorname{Tr}H^2 &= \frac34J^2, \\ [H,S_{\mathrm{tot}}^z] &= 0, \\ [H,\mathbf S_{\mathrm{tot}}^2] &= 0. \end{aligned}

The fully polarized state has one-quarter unit of exchange on each bond, so

E↑↑↑↑=J.E_{\uparrow\uparrow\uparrow\uparrow} = J.

A result larger by a factor of four indicates Pauli matrices were used as the dimensionless spin operators. A wrong polarized energy with otherwise plausible low levels indicates a bond-list error. Correct eigenvalues with broken (2S+1)(2S+1) degeneracies indicate a faulty total-spin operator or an inconsistent tensor-factor ordering.

For a simple spin-1/21/2 graph with NbN_b distinct bonds and uniform JJ, orthogonality of Pauli strings gives

12LTr⁡H2=3Nb16J2.\frac1{2^L} \operatorname{Tr}H^2 = \frac{3N_b}{16} J^2.

Cross terms between distinct bonds vanish under the full trace, including bonds that share one site. For the four-site ring Nb=4N_b=4, reproducing 3J2/43J^2/4 checks the complete matrix rather than one selected eigenpair.

After the finite-ring test passes, an antiferromagnetic sequence can be compared with

e0J=14−ln⁡2.\frac{e_0}{J} = \frac14-\ln2.

This is a later scaling test, not part of the tiny-ring acceptance condition. Even and odd rings carry different frustration and total-spin structure, so a controlled extrapolation should normally use one parity of LL.

The artifact catalog reserves

notebooks/many-body-statistical/heisenberg_chain_ed.ipynb

for basis construction, symmetry labels, and the complete small-ring spectrum. It is planned rather than published. Until it passes the release gates on Reproducible Notebooks, MB-B002 is the authoritative numerical contract.

VariantWhat changesCanonical route
general spin-ss Heisenberg modellocal dimension becomes 2s+12s+1Heisenberg Model
XXZ chainaxial exchange anisotropy Δ\Delta reduces SU(2)SU(2) to U(1)U(1)XXZ Chain dossier
XYZ chainall three exchange couplings differSpin-1/21/2 Chain dossier
chain in a fieldadd −hStotz-hS_{\mathrm{tot}}^z and track magnetization sectorsHeisenberg Model
spin-1 chainHaldane gap and open-chain edge degrees of freedomHeisenberg Model
frustrated J1J_1–J2J_2 chaincompeting next-neighbor exchangegeneric Bethe integrability is lost
ladder or higher dimensionchange graph and coordinationspin waves, tensor networks, Monte Carlo, or other methods may apply
Hubbard strong-coupling limitexchange emerges as J=4t2/UJ=4t^2/U under stated assumptionsHubbard Chain

The Heisenberg Chain model card and Hamiltonian card are compact lookup surfaces. This dossier is the convention-complete record; the teaching article remains the canonical derivation and interpretation route.

  • Comparing exchange constants before checking whether the operators are s\mathbf s, σ\boldsymbol\sigma, or dimensionful S\mathbf S.
  • Reversing the sign convention and calling J>0J>0 ferromagnetic for the Hamiltonian used here.
  • Treating the Néel product state as the exact antiferromagnetic ground state.
  • Forgetting the spin-exchange terms si+sj−+si−sj+s_i^+s_j^-+s_i^-s_j^+.
  • Counting L−1L-1 bonds on a periodic ring or LL bonds on an open chain.
  • Double-counting the single undirected bond in a two-site periodic shorthand.
  • Ignoring frustration and Kramers structure on odd antiferromagnetic rings.
  • Assuming an StotzS_{\mathrm{tot}}^z block contains only one total-spin multiplet.
  • Calling the antiferromagnetic spinon continuum a single magnon dispersion.
  • Assuming SU(2)SU(2) symmetry is sufficient for Bethe integrability.
  • Applying the spin-1/21/2 gapless result to the integer-spin chain.
  • Inferring long-range order or a bulk gap from one small cluster.

1. Convert between Pauli and spin conventions

Section titled “1. Convert between Pauli and spin conventions”

Suppose a code assembles

H~=J~∑⟨j,k⟩σj⋅σk.\widetilde H = \widetilde J \sum_{\langle j,k\rangle} \boldsymbol\sigma_j \mathbin{\cdot} \boldsymbol\sigma_k.

What value of J~\widetilde J reproduces the baseline Hamiltonian with exchange JJ? What factor error appears if one sets J~=J\widetilde J=J?

Solution

Since

sj⋅sk=14σj⋅σk,\mathbf s_j\cdot\mathbf s_k = \frac14 \boldsymbol\sigma_j \mathbin{\cdot} \boldsymbol\sigma_k,

the matching coefficient is

J~=J4.\widetilde J = \frac J4.

Setting J~=J\widetilde J=J multiplies every energy and every gap by four. Degeneracies and dimensionless eigenvectors may still look correct, which is why the dimer energies and trace moment are essential normalization checks.

Use the two exchange eigenvalues to show that

ΠS=14−s1⋅s2,ΠT=34+s1⋅s2.\Pi_S = \frac14-\mathbf s_1\cdot\mathbf s_2, \qquad \Pi_T = \frac34+\mathbf s_1\cdot\mathbf s_2.

Verify projector completeness and orthogonality.

Solution

On the singlet, s1⋅s2=−3/4\mathbf s_1\cdot\mathbf s_2=-3/4, so ΠS\Pi_S has eigenvalue 11 and ΠT\Pi_T has eigenvalue 00. On the triplet, the exchange eigenvalue is 1/41/4, so the roles reverse. The singlet and triplet span the full four-dimensional two-spin space, hence

ΠS+ΠT=I.\Pi_S+\Pi_T = \mathbb I.

Because both are spectral projectors of the Hermitian exchange operator, they obey

ΠS2=ΠS,ΠT2=ΠT,ΠSΠT=0.\Pi_S^2=\Pi_S, \qquad \Pi_T^2=\Pi_T, \qquad \Pi_S\Pi_T=0.

Let ∣E,S,m⟩\lvert E,S,m\rangle be an eigenstate of HH, Stot2\mathbf S_{\mathrm{tot}}^2, and StotzS_{\mathrm{tot}}^z. Show why all nonzero states obtained by applying Stot±S_{\mathrm{tot}}^\pm have the same energy.

Solution

Rotational invariance gives

[H,Stot±]=0.[H,S_{\mathrm{tot}}^\pm] = 0.

Therefore

HStot±∣E,S,m⟩=Stot±H∣E,S,m⟩=EStot±∣E,S,m⟩.\begin{aligned} H S_{\mathrm{tot}}^\pm \lvert E,S,m\rangle &= S_{\mathrm{tot}}^\pm H \lvert E,S,m\rangle \\ &= E S_{\mathrm{tot}}^\pm \lvert E,S,m\rangle. \end{aligned}

Whenever the ladder operation is nonzero, it produces another state at energy EE with mm shifted by one. Repeating the operation generates all 2S+12S+1 members of the multiplet. A numerical spectrum that splits these levels violates the baseline SU(2)SU(2) symmetry.

For a simple graph with NbN_b distinct spin-1/21/2 bonds, prove

12LTr⁡H2=3Nb16J2.\frac1{2^L} \operatorname{Tr}H^2 = \frac{3N_b}{16}J^2.
Solution

Write each bond as

hjk=J4∑α=x,y,zσjασkα.h_{jk} = \frac J4 \sum_{\alpha=x,y,z} \sigma_j^\alpha\sigma_k^\alpha.

The 3Nb3N_b Pauli strings are distinct and nonidentity. Pauli-string orthogonality gives

12LTr⁡(QaQb)=δab.\frac1{2^L} \operatorname{Tr} \left( Q_aQ_b \right) = \delta_{ab}.

Thus only the squares of identical strings survive in the full trace:

12LTr⁡H2=3Nb(J4)2=3Nb16J2.\frac1{2^L} \operatorname{Tr}H^2 = 3N_b \left( \frac J4 \right)^2 = \frac{3N_b}{16}J^2.

The argument also shows why cross terms from adjacent bonds vanish despite sharing one site.

5. Obtain the ferromagnetic finite-size scale

Section titled “5. Obtain the ferromagnetic finite-size scale”

Use

εF(k)=∣J∣[1−cos⁡(ka)]\varepsilon_F(k) = \lvert J\rvert \left[ 1-\cos(ka) \right]

to find the first nonzero one-magnon energy on a periodic ring and its large-LL behavior.

Solution

Allowed momenta are

kn=2πnLa.k_n = \frac{2\pi n}{La}.

The n=0n=0 state is a member of the ground multiplet. The smallest nonzero magnitude is 2π/(La)2\pi/(La), so

ΔLmag=∣J∣[1−cos⁡(2πL)].\Delta_L^{\mathrm{mag}} = \lvert J\rvert \left[ 1-\cos \left( \frac{2\pi}{L} \right) \right].

Using 1−cos⁡x=x2/2+O(x4)1-\cos x=x^2/2+O(x^4),

ΔLmag=2π2∣J∣L2+O(L−4).\Delta_L^{\mathrm{mag}} = \frac{2\pi^2\lvert J\rvert}{L^2} + O(L^{-4}).

The quadratic scaling is the finite-size signature of the ferromagnetic magnon near k=0k=0.

6. Compare the Néel estimate with the exact energy

Section titled “6. Compare the Néel estimate with the exact energy”

For the antiferromagnetic even chain, compute the exchange expectation per bond in the alternating product state and compare it with the exact thermodynamic energy density.

Solution

On an antiparallel product bond,

⟨sizsi+1z⟩=−14,\left\langle s_i^zs_{i+1}^z \right\rangle = -\frac14,

while the transverse terms have zero expectation. Thus

⟨H⟩NeˊelLJ=−14\frac{ \langle H\rangle_{\mathrm{N\acute eel}} }{ LJ } = -\frac14

for a periodic even ring. The exact thermodynamic value is

e0J=14−ln⁡2≃−0.443147.\frac{e_0}{J} = \frac14-\ln2 \simeq -0.443147.

It is lower, as the variational principle requires. The difference quantifies energy gained from quantum correlations that the classical product state omits.

  1. W. Heisenberg, “Zur Theorie des Ferromagnetismus,” Zeitschrift für Physik 49, 619–636 (1928), doi:10.1007/BF01328601.
  2. H. Bethe, “Zur Theorie der Metalle. I. Eigenwerte und Eigenfunktionen der linearen Atomkette,” Zeitschrift für Physik 71, 205–226 (1931), doi:10.1007/BF01341708.
  3. T. Holstein and H. Primakoff, “Field dependence of the intrinsic domain magnetization of a ferromagnet,” Physical Review 58, 1098–1113 (1940), doi:10.1103/PhysRev.58.1098.
  4. J. des Cloizeaux and J. J. Pearson, “Spin-wave spectrum of the antiferromagnetic linear chain,” Physical Review 128, 2131–2135 (1962), doi:10.1103/PhysRev.128.2131.
  5. N. D. Mermin and H. Wagner, “Absence of ferromagnetism or antiferromagnetism in one- or two-dimensional isotropic Heisenberg models,” Physical Review Letters 17, 1133–1136 (1966), doi:10.1103/PhysRevLett.17.1133.
  6. F. D. M. Haldane, “Nonlinear field theory of large-spin Heisenberg antiferromagnets,” Physical Review Letters 50, 1153–1156 (1983), doi:10.1103/PhysRevLett.50.1153.
  7. S. R. White and D. A. Huse, “Numerical renormalization-group study of low-lying eigenstates of the antiferromagnetic S=1S=1 Heisenberg chain,” Physical Review B 48, 3844–3852 (1993), doi:10.1103/PhysRevB.48.3844.
  8. A. Auerbach, Interacting Electrons and Quantum Magnetism, Springer (1994).
  9. M. Takahashi, Thermodynamics of One-Dimensional Solvable Models, Cambridge University Press (1999).
  10. T. Giamarchi, Quantum Physics in One Dimension, Oxford University Press (2004).
  11. S. Sachdev, Quantum Phase Transitions, 2nd ed., Cambridge University Press (2011), doi:10.1017/CBO9780511973765.