Skip to content

The spin-1/21/2 XXZ chain is an axially anisotropic nearest-neighbor Heisenberg chain whose single parameter Δ\Delta connects a polarized ferromagnet, a gapless c=1c=1 Luttinger liquid, the isotropic antiferromagnet, and a gapped Néel phase while retaining Bethe-ansatz integrability.

The phase diagram below is at zero temperature and zero longitudinal field, with positive transverse exchange. The same parameter changes the interaction between Jordan–Wigner fermions, so only one point is a free-fermion model.

Helpful background. Heisenberg Chain explains the isotropic spin normalization and multiplets; XXZ Spin Chain develops the Bethe and Luttinger descriptions.

The reference chain has these physical choices.

Field Baseline choice
sites LL spin-1/21/2 degrees of freedom
local space Hj≃C2\mathcal H_j\simeq\mathbb C^2
operators dimensionless sj=σj/2\mathbf s_j=\boldsymbol\sigma_j/2
geometry uniform one-dimensional nearest-neighbor chain
transverse exchange J>0J>0
longitudinal exchange JΔJ\Delta with real dimensionless Δ\Delta
field h=0h=0
boundaries open or periodic, stated explicitly
energy offset none
disorder, dimerization, and longer range absent
thermodynamic reference even periodic rings followed by L→∞L\to\infty

The zero-field phase labels below assume J>0J>0. A longitudinal field is a useful integrable deformation, but its value must be stated because it changes magnetization, filling, and parts of the phase diagram. Negative JJ, transverse fields, bond alternation, and frustrated exchange are not silently absorbed into the baseline.

The 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 spin operators satisfy

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

and

sjα=12σjα.s_j^\alpha = \frac12 \sigma_j^\alpha.

The product basis uses

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

Writing

sj±=sjx±isjys_j^\pm = s_j^x \pm i s_j^y

gives the bond identity

sjxsj+1x+sjysj+1y=12(sj+sj+1−+sj−sj+1+).s_j^xs_{j+1}^x + s_j^ys_{j+1}^y = \frac12 \left( s_j^+s_{j+1}^- + s_j^-s_{j+1}^+ \right).

The longitudinal term is diagonal in the product basis. The transverse term exchanges an adjacent up-down pair and conserves the total number of down spins.

For an open chain,

HO=J∑j=1L−1(sjxsj+1x+sjysj+1y+Δsjzsj+1z).H_{\mathrm O} = J \sum_{j=1}^{L-1} \big( s_j^xs_{j+1}^x + s_j^ys_{j+1}^y + \Delta s_j^zs_{j+1}^z \big).

For a periodic chain,

HP=J∑j=1L(sjxsj+1x+sjysj+1y+Δsjzsj+1z),H_{\mathrm P} = J \sum_{j=1}^{L} \big( s_j^xs_{j+1}^x + s_j^ys_{j+1}^y + \Delta s_j^zs_{j+1}^z \big),

with

sL+1α=s1α.s_{L+1}^{\alpha} = s_1^\alpha.

The open chain has L−1L-1 undirected bonds; an ordinary ring has LL. As in other nearest-neighbor models, the L=2L=2 periodic shorthand can count the same undirected pair twice. A two-site check therefore uses one explicitly open bond.

A longitudinal field measured in energy units adds

Hh=−h∑j=1Lsjz.H_h = -h \sum_{j=1}^{L} s_j^z.

In terms of Pauli matrices,

H=J4∑⟨j,k⟩(σjxσkx+σjyσky+Δσjzσkz)−h2∑jσjz.H = \frac J4 \sum_{\langle j,k\rangle} \big( \sigma_j^x\sigma_k^x + \sigma_j^y\sigma_k^y + \Delta \sigma_j^z\sigma_k^z \big) - \frac h2 \sum_j\sigma_j^z.

A result four times too large in the exchange spectrum or twice too large in the field response indicates that Pauli and spin conventions were mixed.

The exchange constants are

Jx=Jy=J,Jz=JΔ.J_x = J_y = J, \qquad J_z = J\Delta.

Δ\Delta is dimensionless. The labels easy plane and easy axis refer to this tensor, not to a classical orientation already chosen by a finite quantum ground state.

For every Δ\Delta and longitudinal field hh,

Stotz=∑j=1LsjzS_{\mathrm{tot}}^z = \sum_{j=1}^{L} s_j^z

is conserved:

[H,Stotz]=0.[H,S_{\mathrm{tot}}^z] = 0.

The internal continuous symmetry is generically U(1)U(1). If MM sites are down relative to the all-up state,

Stotz=L2−M,dim⁡HM=(LM).S_{\mathrm{tot}}^z = \frac L2-M, \qquad \dim\mathcal H_M = \binom LM.

This is both a physical charge decomposition and the natural exact-diagonalization basis.

At

Δ=1,h=0,\Delta=1, \qquad h=0,

the Hamiltonian becomes isotropic:

H=J∑⟨j,k⟩sj⋅sk.H = J \sum_{\langle j,k\rangle} \mathbf s_j\cdot\mathbf s_k.

All three components of total spin commute with HH, enhancing U(1)U(1) to SU(2)SU(2). Heisenberg Chain describes the resulting multiplets and exact finite spectra.

At h=0h=0, a global π\pi rotation about xx maps

sjz↦−sjzs_j^z \mapsto -s_j^z

and exchanges magnetization sectors M↔L−MM\leftrightarrow L-M. It leaves the XXZ exchange invariant.

Physical time reversal sends every spin component to its negative and also preserves the zero-field Hamiltonian. A nonzero longitudinal field breaks both magnetization reversal and time reversal, while retaining the axial U(1)U(1) symmetry.

The uniform periodic chain has one-site translation and reflection symmetry. The open chain retains reflection about its midpoint but not translation. Bond alternation, disorder, twists, and boundary fields can change the spatial symmetry without necessarily breaking magnetization conservation.

For J>0J>0, h=0h=0, and the thermodynamic ground state:

Anisotropy Ground-state regime Bulk gap Long-distance character
Δ<−1\Delta<-1 zz-polarized ferromagnet nonzero away from −1-1 saturated magnetization
Δ=−1\Delta=-1 first-order endpoint quadratic soft mode enlarged ground space on an even bipartite chain
−1<Δ<1-1<\Delta<1 easy-plane Luttinger liquid zero algebraic correlations and c=1c=1
Δ=1\Delta=1 isotropic antiferromagnet zero c=1c=1 with logarithmic corrections
Δ>1\Delta>1 easy-axis antiferromagnet nonzero thermodynamic Néel order

The critical interval is an entire line of fixed points, not one isolated critical point. Its correlation exponents vary continuously with Δ\Delta. The endpoint at Δ=−1\Delta=-1 is first order and has a vanishing linear-mode velocity; the transition at Δ=1\Delta=1 is Berezinskii–Kosterlitz–Thouless and opens a gap through an essential singularity.

A finite chain does not literally display all thermodynamic order parameters. Phase assignments require size sequences, long-distance correlations, sector information, or a symmetry-breaking source taken to zero after the thermodynamic limit.

At

Δ=0,\Delta=0,

Jordan–Wigner maps the chain to free spinless fermions. At zero field and half filling,

e0J=−1π,v=Jaℏ.\frac{e_0}{J} = -\frac1\pi, \qquad v = \frac{Ja}{\hbar}.

The fermions are free, but transverse spin operators contain nonlocal strings. Their correlations are not reduced to a one-body spin observable.

At

Δ=1,\Delta=1,

the model is the antiferromagnetic spin-1/21/2 Heisenberg chain. Exact bulk fingerprints include

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

and

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

The point remains gapless, but marginally irrelevant interactions produce slow logarithmic corrections.

On an even bipartite chain, a staggered π\pi rotation about zz maps the Δ=−1\Delta=-1 Hamiltonian to an isotropic ferromagnetic Heisenberg chain. Its ground-space dimension is therefore

L+1L+1

for spin 1/21/2. The low-energy dispersion is quadratic, so the endpoint is not an ordinary finite-velocity conformal point.

For Δ≫1\Delta\gg1, the longitudinal exchange dominates and the two alternating zz configurations organize the low-energy antiferromagnetic sector. Transverse exchange dresses them with quantum fluctuations. For Δ≪−1\Delta\ll-1, the two fully polarized zz product states are exact zero-field ground states.

The uniform nearest-neighbor spin-1/21/2 chain is Bethe-ansatz integrable for the standard anisotropies and a longitudinal field. Exactness remains quantity and boundary dependent.

Quantity or regime Status
finite-volume spectrum encoded by Bethe roots after sector and boundary conventions are fixed
Δ=0\Delta=0 exactly quadratic after Jordan–Wigner
one-magnon sector exact plane-wave result
zero-field phase boundaries exact
bulk energy and dressed thermodynamics exact integral-equation framework
zero-field KK and vv in the critical interval exact
finite-size conformal spectrum universal leading form with model-specific corrections
correlation exponents in the critical interval exact once the bosonization convention is fixed
correlation amplitudes and dynamical functions substantially harder; exact representations need further evaluation
generic next-neighbor, disordered, or transversely driven chain not generically integrable

Bethe solvability does not make the model free. For Δ≠0\Delta\ne0, quasiparticle momenta or rapidities are coupled by nontrivial scattering phases, and spin correlations require matrix elements as well as energies.

Fermions, Magnons, and Luttinger Parameters

Section titled “Fermions, Magnons, and Luttinger Parameters”

With

sjz=nj−12,s_j^z = n_j-\frac12,

the open chain maps to

HO=J2∑j=1L−1(cj†cj+1+cj+1†cj)+JΔ∑j=1L−1(nj−12)(nj+1−12)−h∑j=1L(nj−12).H_{\mathrm O} ={} \frac J2 \sum_{j=1}^{L-1} \left( c_j^\dagger c_{j+1} + c_{j+1}^\dagger c_j \right) + J\Delta \sum_{j=1}^{L-1} \left( n_j-\frac12 \right) \left( n_{j+1}-\frac12 \right) - h \sum_{j=1}^{L} \left( n_j-\frac12 \right).

Thus the spinless-fermion hopping and nearest-neighbor interaction are

t=J2,V=JΔ,t = \frac J2, \qquad V = J\Delta,

up to a removable sign of tt from a staggered phase convention on this open chain. On a periodic odd ring, the same transformation changes the boundary twist and cannot silently remove the sign. Δ=0\Delta=0 is free; Δ≠0\Delta\ne0 retains a density interaction.

For a periodic chain in a longitudinal field, the all-up energy is

EF=JΔL4−hL2.E_F = \frac{J\Delta L}{4} - \frac{hL}{2}.

One down spin with physical wave number qq has excitation energy

ε(q)=h−JΔ+Jcos⁡(qa).\varepsilon(q) = h - J\Delta + J\cos(qa).

For an even periodic ring, or in the thermodynamic limit, the minimum lies at q=π/aq=\pi/a. For Δ≥−1\Delta\ge-1 the saturation field is

h≥hsat,hsat=J(1+Δ).h \ge h_{\mathrm{sat}}, \qquad h_{\mathrm{sat}} = J(1+\Delta).

This is both a necessary one-magnon condition and a sufficient ground-state condition in that regime. To see sufficiency, write each bond with half of the field on each end. Its four energies are

JΔ/4−h/2,JΔ/4+h/2,−JΔ/4±J/2.J\Delta/4-h/2,\quad J\Delta/4+h/2,\quad -J\Delta/4\pm J/2.

For h≥J(1+Δ)≥0h\ge J(1+\Delta)\ge0, the all-up state minimizes every bond, hence their sum. At equality other states can become degenerate. A finite odd ring need not contain q=π/aq=\pi/a and has a different finite-size threshold.

For Δ<−1\Delta<-1, the thermodynamic equilibrium positive saturation threshold is instead h=0h=0: any h>0h>0 selects all-up over all-down. The negative value J(1+Δ)J(1+\Delta) is merely a one-magnon stability bound for the all-up branch, not its global ground-state boundary. At zero field in that regime, the magnon gap is

Δmag=J(−Δ−1).\Delta_{\mathrm{mag}} = J(-\Delta-1).

It closes at Δ=−1\Delta=-1.

In the critical interval, write

Δ=cos⁡γ,0<γ<π.\Delta = \cos\gamma, \qquad 0<\gamma<\pi.

At zero magnetization, one common bosonization convention gives

K=π2(π−γ)K = \frac{\pi}{ 2(\pi-\gamma) }

and

v=Jaℏπsin⁡γ2γ.v = \frac{Ja}{\hbar} \frac{ \pi\sin\gamma }{ 2\gamma }.

The limiting checks are

ΔKv01Ja/ℏ1−1/2πJa/(2ℏ)−1+∞0\begin{array}{c|c|c} \Delta & K & v \\ \hline 0 & 1 & Ja/\hbar \\ 1^- & 1/2 & \pi Ja/(2\hbar) \\ -1^+ & \infty & 0 \end{array}

and

c=1c = 1

throughout −1<Δ≤1-1<\Delta\le1, with endpoint and logarithmic qualifications as stated above.

With the same KK convention, leading zero-magnetization correlations behave schematically as

⟨s0zsrz⟩∼−K2π2r2+Az(−1)rr2K,⟨s0+sr−⟩∼A⊥(−1)rr1/(2K).\begin{aligned} \langle s_0^zs_r^z\rangle \sim {}& -\frac{K}{2\pi^2r^2} + A_z \frac{(-1)^r}{r^{2K}}, \\ \langle s_0^+s_r^-\rangle \sim {}& A_\perp \frac{(-1)^r}{r^{1/(2K)}}. \end{aligned}

Here rr is an integer site separation. The amplitudes are nonuniversal. At Δ=1\Delta=1, multiplicative logarithms modify the pure powers. At finite magnetization, both KK and the oscillation wavevectors change.

Boundary Conditions and Finite-Size Effects

Section titled “Boundary Conditions and Finite-Size Effects”

Three boundary statements must be kept distinct:

  1. open and periodic spin chains have different bond counts;
  2. periodic Jordan–Wigner fermions acquire a boundary sign fixed by total fermion parity;
  3. Bethe equations depend on periodicity, twists, magnetization sector, and rapidity convention.

These effects are subextensive in a bulk energy density but decisive for exact finite spectra, momenta, and gaps.

An even ring is the clean zero-field thermodynamic sequence. In the easy-axis antiferromagnet, an odd periodic ring frustrates the alternating order. Near the BKT point, very large correlation lengths and logarithmic corrections make modest system sizes especially deceptive.

A twist can be imposed through

sL+1±=e±iΦs1±.s_{L+1}^{\pm} = e^{\pm i\Phi} s_1^\pm.

The curvature of the ground energy with Φ\Phi probes spin stiffness, but conventions differ by factors of LL, π\pi, and ℏ\hbar.

The uniform magnetization per site is

mz=1L∑j⟨sjz⟩.m_z = \frac1L \sum_j \langle s_j^z\rangle.

For an even chain, a staggered operator is

Mstz=∑j(−1)jsjz.M_{\mathrm{st}}^z = \sum_j (-1)^j s_j^z.

In a finite spin-reversal eigenstate, ⟨Mstz⟩\langle M_{\mathrm{st}}^z\rangle can vanish even in the easy-axis regime. Useful diagnostics include ⟨(Mstz)2⟩\langle(M_{\mathrm{st}}^z)^2\rangle, the staggered structure factor, long-distance zzzz correlations, and a controlled symmetry-breaking source.

The anisotropy makes longitudinal and transverse correlators inequivalent:

Czz(r)=⟨sjzsj+rz⟩,C_{zz}(r) = \langle s_j^z s_{j+r}^z\rangle,

and

C+−(r)=⟨sj+sj+r−⟩.C_{+-}(r) = \langle s_j^+s_{j+r}^-\rangle.

Their Fourier transforms distinguish zz-ferromagnetic, critical easy-plane, and Néel regimes. Structure Factors explains normalization and scattering conventions.

The finite gap

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

must identify magnetization, momentum, parity, and boundary sectors. In the critical phase, low gaps scale as 1/L1/L with corrections. In either gapped phase they approach nonzero bulk scales, aside from symmetry-partner or boundary splittings.

At criticality, interval entanglement has the leading c=1c=1 logarithm. In the gapped phases it saturates beyond the correlation length, with possible finite-cat or boundary contributions.

The continuity equation for sjzs_j^z defines one current convention:

Jjz=Jℏ(sjxsj+1y−sjysj+1x).\mathcal J_j^z = \frac J\hbar \left( s_j^x s_{j+1}^y - s_j^y s_{j+1}^x \right).

Transport claims require temperature, magnetization, frequency, system-size, and order-of-limits conventions. Integrability can protect ballistic contributions in selected regimes but does not make every transport coefficient trivial.

For −1<Δ<1-1<\Delta<1 away from the free point Δ=0\Delta=0, the chain is both interacting and exactly integrable. Its low-energy physics is collective: the Luttinger parameter continuously changes correlation exponents even though the central charge remains 11.

Fractionalization and variable descriptions

Section titled “Fractionalization and variable descriptions”

The same microscopic system can be described as local spins, interacting spinless fermions, Bethe quasiparticles, or a compact boson at low energy. These descriptions are complementary rather than interchangeable. A local spin operator may be nonlocal in fermions, and the continuum theory omits lattice-scale observables.

At Δ=−1\Delta=-1, a quadratic soft mode and enlarged ground space accompany a first-order boundary. At Δ=1\Delta=1, a marginal interaction drives a BKT transition whose gap opens essentially rather than through a simple power. Treating both as ordinary critical endpoints loses the central physics.

The easy-axis phase has two symmetry-related Néel states in the thermodynamic description, while a finite ring can retain translation or spin-reversal symmetry. Conversely, large finite staggered correlations inside the critical phase do not imply a nonzero order parameter.

A short-range one-dimensional XXZ chain has no nonzero-temperature transition into true ferromagnetic or Néel long-range order. Low-temperature quantum-critical windows and long correlation lengths are crossovers, not additional equilibrium phase transitions.

Minimal Worked Example: One Anisotropic Bond

Section titled “Minimal Worked Example: One Anisotropic Bond”

For one open bond at zero field,

H2=J(s1xs2x+s1ys2y+Δs1zs2z).H_2 = J \left( s_1^xs_2^x + s_1^ys_2^y + \Delta s_1^zs_2^z \right).

The parallel states have

E↑↑=E↓↓=JΔ4.E_{\uparrow\uparrow} = E_{\downarrow\downarrow} = \frac{J\Delta}{4}.

In the zero-magnetization sector, define

∣t0⟩=∣↑↓⟩+∣↓↑⟩2,∣s⟩=∣↑↓⟩−∣↓↑⟩2.\begin{aligned} \lvert t_0\rangle &= \frac{ \lvert\uparrow\downarrow\rangle + \lvert\downarrow\uparrow\rangle }{\sqrt2}, \\ \lvert s\rangle &= \frac{ \lvert\uparrow\downarrow\rangle - \lvert\downarrow\uparrow\rangle }{\sqrt2}. \end{aligned}

Their energies are

Et0=J(2−Δ)4,Es=−J(2+Δ)4.\begin{aligned} E_{t_0} &= \frac{J(2-\Delta)}{4}, \\ E_s &= -\frac{J(2+\Delta)}{4}. \end{aligned}

At Δ=1\Delta=1, the three triplet states become degenerate at J/4J/4 and the singlet lies at −3J/4-3J/4. At Δ=−1\Delta=-1, the singlet crosses the two parallel states. The trace is zero for every Δ\Delta:

2E↑↑+Et0+Es=0.2E_{\uparrow\uparrow} + E_{t_0} + E_s = 0.

The normalized second moment is

14Tr⁡H22=J216(2+Δ2).\frac14 \operatorname{Tr}H_2^2 = \frac{J^2}{16} \left( 2+\Delta^2 \right).

The dimer catches normalization and the Δ=−1\Delta=-1 crossing, but it does not reproduce the thermodynamic Luttinger phase or BKT transition.

Use L=4L=4, periodic boundaries, four distinct undirected bonds,

J=1,Δ=12,h=0,J=1, \qquad \Delta=\frac12, \qquad h=0,

and no constant offset. The anisotropy is interacting after Jordan–Wigner and is not an SU(2)SU(2) point.

The complete spectrum is

E/JE/J degeneracy
(−1−33)/4(-1-\sqrt{33})/4 1
−1-1 2
−1/2-1/2 1
00 7
+1/2+1/2 2
+1+1 2
(−1+33)/4(-1+\sqrt{33})/4 1

The ground energy is

E0J=−1.686140661634507.\frac{E_0}{J} = -1.686140661634507.

The degeneracies sum to 1616, and the fixed-MM block dimensions are

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

Useful independent comparisons are

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

At h=0h=0, spin reversal also maps the MM block to the L−ML-M block with the same spectrum.

For a simple graph with NbN_b distinct XXZ bonds,

12LTr⁡H2=NbJ216(2+Δ2)\frac1{2^L} \operatorname{Tr}H^2 = \frac{ N_bJ^2 }{16} \left( 2+\Delta^2 \right)

at zero field. Distinct Pauli strings are orthogonal under the full trace, so cross-bond terms vanish even for adjacent bonds. With Nb=4N_b=4 and Δ=1/2\Delta=1/2, this gives 9J2/169J^2/16.

After the finite matrix passes, separate thermodynamic calculations can test

e0J∣Δ=0=−1π\left. \frac{e_0}{J} \right|_{\Delta=0} = -\frac1\pi

and

e0J∣Δ=1=14−ln⁡2.\left. \frac{e_0}{J} \right|_{\Delta=1} = \frac14-\ln2.

These values use the unshifted spin-operator Hamiltonian. Adding −JΔ/4-J\Delta/4 per bond changes energy densities but not eigenstates or gaps.

Variant What changes Further calculation
isotropic point set Δ=1\Delta=1 and recover SU(2)SU(2) Heisenberg Chain dossier
XX model set Δ=0\Delta=0 and obtain free fermions Spinless Fermion Chains
longitudinal field change magnetization and fermion filling XXZ Spin Chain
twist probe stiffness through boundary phase Φ\Phi Boundary Conditions on Lattices
dimerized chain alternate bond strengths translation and phase structure change
J1J_1–J2J_2 chain add frustrating next-neighbor exchange generic Bethe integrability is lost
random-field chain break translation and study disorder dynamics Many-Body Localization Preview
higher-spin XXZ chain enlarge the onsite representation phase and integrability claims must be re-established

For the isotropic limit, the Heisenberg Chain Hamiltonian card gives a compact operator formula.

Show directly that the exchange Hamiltonian commutes with StotzS_{\mathrm{tot}}^z.

Solution

The longitudinal term is built from sjzs_j^z operators and therefore commutes with StotzS_{\mathrm{tot}}^z. For one transverse term,

[sjz+sj+1z,sj+sj+1−]=sj+sj+1−−sj+sj+1−=0.\begin{aligned} \left[ s_j^z+s_{j+1}^z, s_j^+s_{j+1}^- \right] ={}& s_j^+s_{j+1}^- - s_j^+s_{j+1}^- \\ = {}& 0. \end{aligned}

The Hermitian-conjugate term also commutes, so

[H,Stotz]=0.[H,S_{\mathrm{tot}}^z] = 0.

Each transverse exchange moves a down spin but does not create or destroy one.

Diagonalize the zero-magnetization block of the one-bond Hamiltonian in the basis {∣↑↓⟩,∣↓↑⟩}\{\lvert\uparrow\downarrow\rangle,\lvert\downarrow\uparrow\rangle\}.

Solution

The block is

HM=1=J(−Δ/41/21/2−Δ/4).H_{M=1} = J \begin{pmatrix} -\Delta/4 & 1/2 \\ 1/2 & -\Delta/4 \end{pmatrix}.

Its symmetric and antisymmetric eigenvectors are ∣t0⟩\lvert t_0\rangle and ∣s⟩\lvert s\rangle, with eigenvalues

Et0=J(2−Δ)4E_{t_0} = \frac{J(2-\Delta)}4

and

Es=−J(2+Δ)4.E_s = -\frac{J(2+\Delta)}4.

The parallel states are separate one-dimensional magnetization sectors with energy JΔ/4J\Delta/4.

Write the one-bond Hamiltonian entirely in Pauli matrices. If a code omits the factor 1/41/4, how do the dimer energies and the trace moment change?

Solution

The correct Pauli form is

H2=J4(σ1xσ2x+σ1yσ2y+Δσ1zσ2z).H_2 = \frac J4 \left( \sigma_1^x\sigma_2^x + \sigma_1^y\sigma_2^y + \Delta\sigma_1^z\sigma_2^z \right).

Omitting 1/41/4 multiplies the Hamiltonian and every energy by four. Since the second trace moment is quadratic in HH, it becomes sixteen times too large:

14Tr⁡H~22=16[J216(2+Δ2)].\frac14\operatorname{Tr}\widetilde H_2^2 = 16 \left[ \frac{J^2}{16} \left( 2+\Delta^2 \right) \right].

Use the one-magnon dispersion to derive the positive saturation field for J>0J>0 and Δ>−1\Delta>-1.

Solution

The excitation energy above the all-up state is

ε(q)=h−JΔ+Jcos⁡(qa).\varepsilon(q) = h-J\Delta+J\cos(qa).

For J>0J>0, its minimum occurs at qa=πqa=\pi, where cos⁡(qa)=−1\cos(qa)=-1. Stability requires

εmin⁡=h−J(Δ+1)≥0.\varepsilon_{\min} = h-J(\Delta+1) \ge 0.

Therefore

hsat=J(1+Δ).h_{\mathrm{sat}} = J(1+\Delta).

For an even periodic chain, the bond-energy argument above also proves sufficiency. On a finite odd ring, minimize over the actually allowed wave numbers. The result depends on the spin-operator normalization and the field term being −h∑jsjz-h\sum_js_j^z.

Derive the zero-field normalized second moment for a simple graph of NbN_b XXZ bonds.

Solution

Each bond contributes three Pauli strings with coefficients

J4,J4,JΔ4.\frac J4, \qquad \frac J4, \qquad \frac{J\Delta}{4}.

Every string is traceless, and distinct full-system Pauli strings are orthogonal. Therefore only identical-string squares survive:

12LTr⁡H2=Nb[(J4)2+(J4)2+(JΔ4)2]=NbJ216(2+Δ2).\begin{aligned} \frac1{2^L} \operatorname{Tr}H^2 &= N_b \left[ \left( \frac J4 \right)^2 + \left( \frac J4 \right)^2 + \left( \frac{J\Delta}{4} \right)^2 \right] \\ &= \frac{N_bJ^2}{16} \left( 2+\Delta^2 \right). \end{aligned}

Evaluate KK and vv at Δ=0\Delta=0, and take their limits as Δ→1−\Delta\to1^- and Δ→−1+\Delta\to-1^+.

Solution

With Δ=cos⁡γ\Delta=\cos\gamma,

K=π2(π−γ),v=Jaℏπsin⁡γ2γ.K = \frac{\pi}{2(\pi-\gamma)}, \qquad v = \frac{Ja}{\hbar} \frac{\pi\sin\gamma}{2\gamma}.

At Δ=0\Delta=0, γ=π/2\gamma=\pi/2, so

K=1,v=Jaℏ.K=1, \qquad v=\frac{Ja}{\hbar}.

As Δ→1−\Delta\to1^-, γ→0\gamma\to0 and sin⁡γ/γ→1\sin\gamma/\gamma\to1, giving

K→12,v→πJa2ℏ.K\to\frac12, \qquad v\to\frac{\pi Ja}{2\hbar}.

As Δ→−1+\Delta\to-1^+, γ→π\gamma\to\pi, so

K→∞,v→0.K\to\infty, \qquad v\to0.

The last limit signals the singular ferromagnetic endpoint rather than an ordinary finite-velocity conformal point.

  • 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.
  • F. Franchini, An Introduction to Integrable Techniques for One-Dimensional Quantum Systems, Springer (2017), doi:10.1007/978-3-319-48487-7.
  • M. Gaudin, The Bethe Wavefunction, translated by J.-S. Caux, Cambridge University Press (2014), doi:10.1017/CBO9781107053885.
  • T. Giamarchi, Quantum Physics in One Dimension, Oxford University Press (2004), doi:10.1093/acprof:oso/9780198525004.001.0001.
  • F. D. M. Haldane, “‘Luttinger liquid theory’ of one-dimensional quantum fluids. I,” Journal of Physics C 14, 2585–2609 (1981), doi:10.1088/0022-3719/14/19/010.
  • A. Luther and I. Peschel, “Calculation of critical exponents in two dimensions from quantum field theory in one dimension,” Physical Review B 12, 3908–3917 (1975), doi:10.1103/PhysRevB.12.3908.
  • R. Orbach, “Linear antiferromagnetic chain with anisotropic coupling,” Physical Review 112, 309–316 (1958), doi:10.1103/PhysRev.112.309.
  • M. Takahashi, Thermodynamics of One-Dimensional Solvable Models, Cambridge University Press (1999), doi:10.1017/CBO9780511524332.
  • C. N. Yang and C. P. Yang, “One-dimensional chain of anisotropic spin-spin interactions. I. Proof of Bethe’s hypothesis for ground state in a finite system,” Physical Review 150, 321–327 (1966), doi:10.1103/PhysRev.150.321.
  • C. N. Yang and C. P. Yang, “One-dimensional chain of anisotropic spin-spin interactions. II. Properties of the ground-state energy per lattice site for an infinite system,” Physical Review 150, 327–339 (1966), doi:10.1103/PhysRev.150.327.
  • C. N. Yang and C. P. Yang, “One-dimensional chain of anisotropic spin-spin interactions. III. Applications,” Physical Review 151, 258–264 (1966), doi:10.1103/PhysRev.151.258.