Lecture 13

Three-body consistency: the Yang–Baxter equation

Course map: P3, the quantum consistency condition. Lecture 11 solved the delta gas with a two-body phase, a number, and Lecture 12 watched the classical Toda chain scatter with numbers at its crossings, which commute. This lecture gives the particles a label. The two-body S-matrix becomes a matrix, the order of the collisions starts to matter, and asking that it not matter gives the Yang–Baxter equation. Yang’s matrix satisfies it, and it is the R-matrix from which Lecture 14 builds the Heisenberg chain’s commuting transfer matrices.

13.1 The question

Numbers commute, so for Lecture 11’s gas and Lecture 12’s chain the two routes of Figure 12.6 could not disagree. Give each particle a spin, and a collision can exchange the spins as well as the momenta. So the question is

when particles carry labels, the S-matrix is a matrix: what makes factorised scattering consistent?

The answer: the Yang–Baxter equation (13.6), which the three-body contact conditions force (Theorem 13.3), and which Yang’s matrix satisfies (Proposition 13.4). The worked example is two labelled particles (Proposition 13.1).

13.2 Particles with a label

Let each particle carry a label in Cd\C^{d} — a spin component, a hyperfine state — and let the delta be blind to it. A two-particle wavefunction takes values in Cd⊗Cd\C^{d}\otimes\C^{d}, the first factor for particle 11’s label, and is symmetric (bosons) or antisymmetric (fermions) under exchanging positions and labels together: ψ(x1,x2)=±P ψ(x2,x1), \psi(x_1,x_2)=\pm P\,\psi(x_2,x_1), (13.1) with PP the permutation of the two factors (Definition 9.4). On x1<x2x_1<x_2 write ψ(x1,x2)=A ei(k1x1+k2x2)+A′ ei(k2x1+k1x2),A,A′∈Cd⊗Cd,\psi(x_1,x_2)=A\,e^{\ii(k_1x_1+k_2x_2)}+A'\,e^{\ii(k_2x_1+k_1x_2)},\qquad A,A'\in\C^{d}\otimes\C^{d}, and let (13.1) fix ψ\psi on x2<x1x_2<x_1.

Proposition 13.1 (Yang’s matrix). The delta interaction forces A′=Y(k1−k2)AA'=Y(k_1-k_2)A, with YB(k)=kP−ick+ic  (bosons),YF(k)=−kP+ick+ic  (fermions). Y_{\mathrm B}(k)=\frac{kP-\ii c}{k+\ii c}\ \ \text{(bosons)},\qquad Y_{\mathrm F}(k)=-\frac{kP+\ii c}{k+\ii c}\ \ \text{(fermions)}. (13.2)

Proof. PP commutes with everything here, so split A,A′A,A' by Π±=12(1±P)\Pi_\pm=\tfrac12(\Id\pm P) into parts symmetric and antisymmetric in the labels. Where (13.1) makes ψ\psi symmetric in the positions — Π+\Pi_+ for bosons, Π−\Pi_- for fermions — the problem is Proposition 11.1’s, and A′=S(k)AA'=S(k)A. In the other part ψ\psi is antisymmetric in the positions and vanishes at x1=x2x_1=x_2, so the delta does not act, and continuity forces A′=−AA'=-A. Hence YB=S Π+−Π−Y_{\mathrm B}=S\,\Pi_+-\Pi_- and YF=S Π−−Π+Y_{\mathrm F}=S\,\Pi_--\Pi_+. With S−1=−2ic/(k+ic)S-1=-2\ii c/(k+\ii c) and S+1=2k/(k+ic)S+1=2k/(k+\ii c) these are (13.2). 

For spin-12\tfrac12 fermions, then, the interacting phase lives in the spin singlet, and the triplet scatters as free fermions do: two fermions with parallel spins never meet.

Reading

In terms of the tt and rr of §11.3’s reading, YB=t(k) P+r(k)Y_{\mathrm B}=t(k)\,P+r(k): transmission carries the labels across the collision, so the particle now on the left has the label that was on the right — the PP — and reflection leaves them where they were. For fermions YF=−t(k) P+r(k)Y_{\mathrm F}=-t(k)\,P+r(k).

Following the momentum. YY is written in the order of positions (slot 11 is the particle on the left). Figure 11.5 follows the momenta, and there the natural object is X(u)=PY,u=k1−k2c:XB(u)=R(u)u+i,XF(u)=R(−u)u+i=−L(u)u+i, X(u)=PY,\qquad u=\frac{k_1-k_2}c:\qquad X_{\mathrm B}(u)=\frac{\Rmat(u)}{u+\ii},\quad X_{\mathrm F}(u)=\frac{\Rmat(-u)}{u+\ii}=-\frac{L(u)}{u+\ii}, (13.3) where R(u)=u−iP\Rmat(u)=u-\ii P is the course’s rational R-matrix and L(u)=u+iPL(u)=u+\ii P the Lax operator of the Heisenberg chain (Lecture 14). Both obey unitarity, X(u)X(−u)=1X(u)X(-u)=\Id. XBX_{\mathrm B} is the R-matrix of Lecture 14; XFX_{\mathrm F}, which Yang found for spin-12\tfrac12 fermions, is the u+iPu+\ii P of the reading at the end of this lecture up to a scalar, and it is why the algebraic Bethe ansatz of Lecture 15 solves the labelled gases of optional lecture C.

13.3 Three-body consistency: the Yang–Baxter equation

Three labelled particles, on x1<x2<x3x_1<x_2<x_3, with the factorised ansatz ψ=∑σ∈S3Aσ ei(kσ(1)x1+kσ(2)x2+kσ(3)x3),Aσ∈(Cd)⊗3,k1>k2>k3. \psi=\sum_{\sigma\in S_3}A_\sigma\,e^{\ii(k_{\sigma(1)}x_1+k_{\sigma(2)}x_2+k_{\sigma(3)}x_3)}, \qquad A_\sigma\in(\C^{d})^{\otimes3},\quad k_1>k_2>k_3 . (13.4) The term σ=id\sigma=\mathrm{id}, fastest on the left, is the incoming wave; σ=w0\sigma=w_0, the order reversed, is the outgoing one; the other four come after one or two collisions. The ansatz is non-diffraction written down; the question is whether the contact conditions allow it.

The walls. On the wall xj+1=xjx_{j+1}=x_j the terms pair up as before, σ\sigma with σsj\sigma s_j, sj=(j  j+1)s_j=(j\;j{+}1), and the third particle is a spectator. So Proposition 13.1, applied to the labels in slots j,j+1j,j+1, says Aσsj=Yj,j+1(kσ(j)−kσ(j+1))Aσ A_{\sigma s_j}=Y_{j,j+1}\big(k_{\sigma(j)}-k_{\sigma(j+1)}\big)A_\sigma (13.5) for every σ\sigma and jj: all the conditions. Each is invertible, and read from σsj\sigma s_j it says the same thing: it is one contact condition.

Following the momenta. Reorder the slots of AσA_\sigma so that slot aa holds the label of the particle that carries kak_a: Bσ=ΠσAσB_\sigma=\Pi_\sigma A_\sigma, where Πσ\Pi_\sigma moves the factor in slot jj to slot σ(j)\sigma(j).

Lemma 13.2 (A crossing is an XX). In terms of BB, (13.5) reads Bσsj=Xab(uab) Bσ,a=σ(j), b=σ(j+1),uab=ka−kbc,B_{\sigma s_j}=X_{ab}(u_{ab})\,B_\sigma,\qquad a=\sigma(j),\ b=\sigma(j+1),\quad u_{ab}=\frac{k_a-k_b}c, with XabX_{ab} acting on slots aa and bb (first leg on aa) (Lecture 9’s leg notation).

Proof. Πσsj=ΠσPj,j+1\Pi_{\sigma s_j}=\Pi_\sigma P_{j,j+1}: both send slot jj to slot σ(j+1)\sigma(j+1) and slot j+1j+1 to σ(j)\sigma(j). So Bσsj=ΠσPj,j+1Yj,j+1Aσ=ΠσXj,j+1Πσ−1BσB_{\sigma s_j}=\Pi_\sigma P_{j,j+1}Y_{j,j+1}A_\sigma=\Pi_\sigma X_{j,j+1}\Pi_\sigma^{-1}B_\sigma, and conjugating by Πσ\Pi_\sigma moves the legs of XX from slots j,j+1j,j+1 to σ(j),σ(j+1)\sigma(j),\sigma(j+1). 

So each crossing of momentum lines aa and bb multiplies by Xab(uab)X_{ab}(u_{ab}): Figure 11.5 with matrices. From the incoming wave to the outgoing one there are two routes. Along id→s1→s1s2→w0\mathrm{id}\to s_1\to s_1s_2\to w_0 the crossings are 1212, 1313, 2323; along id→s2→s2s1→w0\mathrm{id}\to s_2\to s_2s_1\to w_0 they are 2323, 1313, 1212 (Figure 13.1).

Figure 13.1. The Yang–Baxter equation, read off three labelled particles. Each line carries a momentum and a label; a crossing of lines aa and bb acts on their labels with Xab(uab)X_{ab}(u_{ab}), uab=(ka−kb)/cu_{ab}=(k_a-k_b)/c (Lemma 13.2). Time runs up, so the first crossing is the rightmost factor. The two orders are the two routes from the incoming wave to the outgoing one; the contact conditions allow a wavefunction only if they agree.

Theorem 13.3 (Three-body consistency). Let the two-body problem give a relation of the form (13.5) with X=PYX=PY invertible and unitary, X21(−u)X12(u)=1X_{21}(-u)X_{12}(u)=\Id — which for XX commuting with PP, as Yang’s do, is X(u)X(−u)=1X(u)X(-u)=\Id. The contact conditions on the three-body ansatz (13.4) have a solution for every incoming amplitude AidA_{\mathrm{id}} if and only if X12(u12) X13(u13) X23(u23)=X23(u23) X13(u13) X12(u12). X_{12}(u_{12})\,X_{13}(u_{13})\,X_{23}(u_{23}) =X_{23}(u_{23})\,X_{13}(u_{13})\,X_{12}(u_{12}) . (13.6) The outgoing amplitude is then Bw0=X12X13X23 BidB_{w_0}=X_{12}X_{13}X_{23}\,B_{\mathrm{id}}.

Proof. The six amplitudes sit at the vertices of a hexagon, the permutations of S3S_3, and the conditions sit on its six edges, each an invertible map from one vertex to the next (Lemma 13.2); unitarity says that the map read from the other end is its inverse, so each edge carries one relation. Given BidB_{\mathrm{id}}, going round either side fixes everything up to Bw0B_{w_0}; the two values agree for every BidB_{\mathrm{id}} if and only if the two products of Figure 13.1 are equal, which is (13.6). If they are equal, the amplitudes so defined satisfy every edge, and Aσ=Πσ−1BσA_\sigma=\Pi_\sigma^{-1}B_\sigma. 

This is “derived, not postulated”: nothing was assumed but that the three-body wave is a sum of rearranged plane waves, and the equation is what that costs. The logic runs one way: the equation is necessary for a consistent factorised ansatz; it does not by itself prove that an interaction is non-diffractive.

Proposition 13.4 (Yang’s matrix satisfies it). XBX_{\mathrm B} and XFX_{\mathrm F} of (13.3) satisfy (13.6) for all real k1,k2,k3k_1,k_2,k_3, for every dd.

Proof. The scalar factors 1/(uab+i)1/(u_{ab}+\ii) are the same three on both sides. For XBX_{\mathrm B} what is left is (13.6) for R(u)=u−iP\Rmat(u)=u-\ii P. Expand both sides; the transpositions obey P12P13=P13P23=P23P12=:CP_{12}P_{13}=P_{13}P_{23}=P_{23}P_{12}=:C and P13P12=P12P23=P23P13=C−1P_{13}P_{12}=P_{12}P_{23}=P_{23}P_{13}=C^{-1}, and P12P13P23=P13=P23P13P12P_{12}P_{13}P_{23}=P_{13}=P_{23}P_{13}P_{12}. The terms with no PP, one PP and three PP’s agree term by term. The terms with two are −(u23+u12)C−u13C−1-(u_{23}+u_{12})C-u_{13}C^{-1} on the left and −u13C−(u12+u23)C−1-u_{13}C-(u_{12}+u_{23})C^{-1} on the right, and they agree because u13=u12+u23:u_{13}=u_{12}+u_{23} : the spectral parameters are differences of momenta, so they add. For XFX_{\mathrm F}, replace every uu by −u-u. 

The proof used only the form u−iPu-\ii P and the additivity of uu; the variant u−i(1+ϵu)Pu-\ii(1+\epsilon u)P, just as harmless-looking, fails (13.6) for every ϵ≠0\epsilon\ne0 (the check script; the reading at the end).

Scalar SS: automatic. Without labels, d=1d=1, P=1P=1 and X=(u−i)/(u+i)=SX=(u-\ii)/(u+\ii)=S: both sides are S12S13S23S_{12}S_{13}S_{23}, and the order of the crossings does not matter. The Toda chain is in the same position: no labels, numbers (displacements) at the crossings, and its braid relation (Figure 12.6) is the scalar, automatic case. The equation earns its keep only with labels.

NN particles sketched. The amplitudes sit on the permutations of SNS_N, the conditions on the edges σ→σsj\sigma\to\sigma s_j, and they are consistent if every closed walk multiplies to 1\Id. By Coxeter’s presentation of SNS_N (cited from Gaudin §10.1) the closed walks are generated by an edge and back (unitarity, X21(−u)X12(u)=1X_{21}(-u)X_{12}(u)=\Id), the squares sjsl=slsjs_js_l=s_ls_j, ∣j−l∣≥2\abs{j-l}\ge2 (automatic: four different lines), and the hexagons sjsj+1sj=sj+1sjsj+1s_js_{j+1}s_j=s_{j+1}s_js_{j+1} (the Yang–Baxter equation). So the Yang–Baxter equation and unitarity make the NN-body factorised ansatz consistent; the check script confirms it for N=4N=4.

13.4 Factorised scattering

The NN-body S-matrix is the ordered product of two-body matrices, one per pair of momentum lines, in the order the lines cross, SN=∏crossings, in time orderXab(uab), S_N=\prod_{\text{crossings, in time order}}X_{ab}(u_{ab}), (13.7) and by the Yang–Baxter equation every order of crossings gives the same product: the many-body problem is determined by the two-body one. On a ring, carrying one particle round gives, as in Theorem 11.3, a condition — now on a vector of labels: an eigenvalue problem for a product of XX’s, a transfer matrix, solved by the nested algebraic Bethe ansatz (optional lecture C).

13.5 What it buys

  • The consistency condition of Part II, derived rather than postulated: a factorised three-body wave exists for every incoming amplitude exactly when the two orders of Figure 13.1 agree.

  • Yang’s matrix, R(u)=u−iP\Rmat(u)=u-\ii P up to a scalar, the R-matrix from which Lecture 14 builds the Heisenberg chain’s commuting transfer matrices.

  • Factorised scattering for NN particles, sketched: with unitarity, the equation makes every order of crossings give the same product (13.7).

13.6 Caveats, in one place

  • Status. Proposition 13.1, Theorem 13.3 and Proposition 13.4 are proved. The NN-body statement is sketched, with Coxeter’s presentation of SNS_N cited, and checked for N=4N=4. The Zamolodchikov–Faddeev algebra of the reading is cited.

  • Necessary, not sufficient. The equation is what a consistent factorised ansatz costs. It does not by itself prove that an interaction is non-diffractive.

  • Unitarity is a second condition, X(u)X(−u)=1X(u)X(-u)=\Id for matrices that commute with PP; Yang’s satisfy it.

  • Without labels the equation is automatic, so Lecture 11’s gas and Lecture 12’s chain never tested it.

Reading

13.7 Three lines crossing, and the equation written down

Suppose three-body scattering is a sequence of two-body events. There are two ways for three lines to cross pairwise. Draw each particle’s path in space and time (a world line), with time running up. Three world lines with three different speeds meet pairwise in three crossings, and the middle line can pass either below or above the crossing of the other two (Figure 13.2). The two orders are different sequences of two-body events, and they describe the same three-body collision. They must agree.

Figure 13.2. The picture behind the Yang–Baxter equation. Three particles, the same three speeds in both panels, time running up; only the starting offsets differ. The middle world line passes on one side of the crossing of the other two or on the other, so the three two-body events (dots, labelled by the pair) happen in opposite orders. If three-body scattering is nothing but its two-body events, both orders must give the same result. Figure 12.6 is its classical form for the Toda chain, and Figure 13.1 its quantum form.

This requirement is the Yang–Baxter equation, and it is a condition on the two-body scattering alone. For particles with nothing but a momentum it is automatic: each two-body event then contributes a number (classically, a shift in position; quantum mechanically, a phase), and numbers do not care about order. It bites when the particles carry an internal label — a spin, a colour — that a two-body event can exchange. Then each event is a matrix, matrices do not commute, and the equation is a strong constraint: a system whose two-body scattering fails it cannot scatter factorisably. Yang met it in 1967 as the consistency condition for exactly this kind of many-body scattering; Baxter met it in 1972 in the statistical mechanics of two-dimensional lattices.

The equation does not organise everything. The geometry of Part I — phase space, the Liouville–Arnold theorem, action-angle variables, the failure of integrability — owes it nothing. Nor does all of the algebra: the Lax pairs of Lecture 8 and the reduction behind Calogero–Moser (optional lecture B) do not come from it. What rests on it is one construction: commuting transfer matrices built from an r\rmat-matrix. There the classical equation is the model-independent condition that makes the bracket consistent — sufficient, not forced, as Lecture 14’s reading explains — and it is realised by the periodic chain of Lecture 9; the quantum equation plays the same part from Lecture 14 on. (The open chain of Lecture 9’s readings has an r\rmat-matrix bracket too, but of another kind, and its involution does not need the equation.) That construction is the course’s spine from Part I into Part II.

For readers who have met spin: let two particles each carry a spin-12\frac12, and let them be fermions with a contact interaction. Their two-body scattering acts on the four-dimensional space C2⊗C2\C^{2}\otimes\C^{2} of the pair by XF(u)=−L(u)/(u+i)X_{\mathrm F}(u)=-L(u)/(u+\ii) of (13.3), with numerator L(u)=u 1+i P,L(u)=u\,\Id+\ii\,\perm, where P\perm swaps the two spins and uu is the difference of the two particles’ rapidities, a parametrisation of their momenta. The scalar −1/(u+i)-1/(u+\ii) makes XFX_{\mathrm F} unitary, and it is the same three factors on both sides of the equation below, so the equation can be written with LL. For three particles, write L12(u)L_{12}(u) for L(u)L(u) acting on the spins of particles 11 and 22 and as the identity on the third, and so on. The two orders of Figure 13.2 are the two sides of L12(u−v) L13(u) L23(v)=L23(v) L13(u) L12(u−v),L_{12}(u-v)\,L_{13}(u)\,L_{23}(v) =L_{23}(v)\,L_{13}(u)\,L_{12}(u-v), an identity between 8×88\times8 matrices. For this LL it holds for all uu and vv; for the apparently harmless variant u 1+i(1+cu)Pu\,\Id+\ii(1+cu)\perm it fails for every c≠0c\ne0 (scripts/check-l0-algebra.py, group 6). This is Yang’s solution (1967), the one behind the Heisenberg chain of Lecture 10. This lecture derives the equation from three-body consistency, and Lecture 14 builds the rest of Part II on it.

Reading

13.8 The two coordinate Bethe ansätze, side by side

Table 13.1. The two coordinate Bethe ansätze side by side.
Heisenberg chain (Lecture 10) delta gas (Lecture 11)
particles magnons on a lattice bosons on a line
variable rapidity λ\lambda, eik=(λ+i/2)/(λ−i/2)e^{\ii k}=(\lambda+\ii/2)/(\lambda-\ii/2) momentum kk
two-body SS (λ−i)/(λ+i)(\lambda-\ii)/(\lambda+\ii) (k−ic)/(k+ic)(k-\ii c)/(k+\ii c), the same at k=cλk=c\lambda
Bethe equations eikjN∏l≠jS(λj−λl)=1e^{\ii k_jN}\prod_{l\ne j}S(\lambda_j-\lambda_l)=1 eikjL∏l≠jS(kj−kl)=1e^{\ii k_jL}\prod_{l\ne j}S(k_j-k_l)=1
energy −∑2J/(4λj2+1)-\sum2J/(4\lambda_j^{2}+1) ∑kj2\sum k_j^{2}
complex roots strings: bound states none for c>0c>0
eigenfunction sketched proved
completeness numerics; a reformulation (cited) proved (cited)

Reading

13.9 The Zamolodchikov–Faddeev algebra

Factorised scattering has an algebraic shorthand. Give each asymptotic particle an operator Z(k)Z(k) that creates it, with one component per label, and write Za(k)Z_a(k) for Z(k)Z(k) with its label in tensor slot aa (Lecture 9’s leg notation). Postulate the exchange relation Z1(k1) Z2(k2)=X12(u12) Z2(k2) Z1(k1):Z_1(k_1)\,Z_2(k_2)=X_{12}(u_{12})\,Z_2(k_2)\,Z_1(k_1): reordering two neighbours costs one two-body matrix. A product of three can be reversed in two ways, giving X12X13X23 Z3Z2Z1X_{12}X_{13}X_{23}\,Z_3Z_2Z_1 and X23X13X12 Z3Z2Z1X_{23}X_{13}X_{12}\,Z_3Z_2Z_1: the two sides of (13.6). Exchanging twice returns the product to itself only if X(u)X(−u)=1X(u)X(-u)=\Id (for matrices that commute with PP, as Yang’s do): unitarity. This is the Zamolodchikov–Faddeev algebra (Zamolodchikov and Zamolodchikov 1979; Faddeev 1980) cited. In relativistic field theory, where it was born, the two-body matrix obeys a third condition, crossing symmetry, which relates a particle going forward in time to an antiparticle going backwards. It has no counterpart for the non-relativistic particles of this course, which is one reason field theory is outside it (Lecture 0).

Summary

  • Labels (Proposition 13.1): two labelled particles scatter with Yang’s matrix YY; along the momentum lines X=PYX=PY, which for bosons is R(u)/(u+i)\Rmat(u)/(u+\ii) with R(u)=u−iP\Rmat(u)=u-\ii P and u=(k1−k2)/cu=(k_1-k_2)/c.
  • Three-body consistency (Theorem 13.3): given unitarity, the factorised ansatz has a solution for every incoming amplitude if and only if X12X13X23=X23X13X12X_{12}X_{13}X_{23}=X_{23}X_{13}X_{12}, the Yang–Baxter equation.
  • Yang’s matrix satisfies it (Proposition 13.4), because the spectral parameters are differences of momenta and add; a harmless-looking variant fails.
  • NN particles (sketched): unitarity and the equation make the ansatz consistent, and the NN-body S-matrix is an ordered product of two-body matrices (13.7).
  • Next: Lecture 14 meets the same equation as the RLL relation, and builds commuting transfer matrices from it.

References

C. N. Yang, Phys. Rev. Lett. 19 (1967) 1312 (labelled particles, and the equation); M. Gaudin, The Bethe Wavefunction, trans. J.-S. Caux, Cambridge University Press (2014), §10.1 (Coxeter’s presentation and the NN-body ansatz); B. Sutherland, Beautiful Models, World Scientific (2004), Chapter 2 (non-diffraction); R. J. Baxter, Exactly Solved Models in Statistical Mechanics, Academic Press (1982) (the equation on lattices); A. B. Zamolodchikov and Al. B. Zamolodchikov, Ann. Phys. 120 (1979) 253; L. D. Faddeev, Sov. Sci. Rev. C1 (1980) 107 (the exchange algebra).

End of Lecture 13

Three-body consistency: the Yang–Baxter equation