Skip to content

The Yang–Baxter Equation

For matrix-valued factorized scattering, the Yang–Baxter equation states that two admissible orders of three pairwise collisions define the same map on the three-particle internal space. It is a consistency equation for the tensor structure of a two-body amplitude. It does not by itself impose unitarity, crossing, analyticity, the correct particle spectrum, completeness, or existence of a local QFT.

Required background. Elasticity, factorization, and their hypotheses supplies the reduction of many-body scattering to ordered two-body factors. Helpful background. Direct sums, tensor products, and index structure supplies tensor-factor notation, while representations, intertwiners, and invariants supplies the symmetry decomposition of matrix amplitudes.

Let VaV_a denote the internal one-particle space. In the equal-mass example below, all particles transform in one space VV. Define the two-body operator S(θ):VVVVS(\theta):V\otimes V\to V\otimes V by

S(θ)a,b=c,dSabcd(θ)c,d.S(\theta)\lvert a,b\rangle = \sum_{c,d}S_{ab}^{cd}(\theta)\lvert c,d\rangle .

Lower indices are inputs, upper indices are outputs, and the first index always labels the left tensor factor. This declaration prevents the common transpose error in component formulas.

On V1V2V3V_1\otimes V_2\otimes V_3, S12S_{12} acts on factors 1 and 2 and as the identity on factor 3. The definition of S13S_{13} includes the tensor-factor permutation needed to act on factors 1 and 3; it is not obtained by merely renaming adjacent component indices.

For on-shell rapidities, Lorentz invariance makes each amplitude depend on θij=θiθj\theta_{ij}=\theta_i-\theta_j. We write operator products in the usual mathematical order: the rightmost factor acts first.

Take three wave packets with θ1>θ2>θ3\theta_1>\theta_2>\theta_3. Factorization reduces their scattering to three two-body encounters. Moving the encounters apart with a higher conserved charge changes their order without changing the total amplitude. The two sequences are

S12(θ12)S13(θ13)S23(θ23)S_{12}(\theta_{12}) S_{13}(\theta_{13}) S_{23}(\theta_{23})

and

S23(θ23)S13(θ13)S12(θ12).S_{23}(\theta_{23}) S_{13}(\theta_{13}) S_{12}(\theta_{12}).

Their equality is the Yang–Baxter equation:

S12(θ12)S13(θ13)S23(θ23)=S23(θ23)S13(θ13)S12(θ12).S_{12}(\theta_{12}) S_{13}(\theta_{13}) S_{23}(\theta_{23}) = S_{23}(\theta_{23}) S_{13}(\theta_{13}) S_{12}(\theta_{12}).

The physical derivation in factorized relativistic scattering appears in Zamolodchikov and Zamolodchikov 1979, § 2, pp. 257–260. The same consistency relation arose in Yang’s coordinate Bethe-ansatz solution of a one-dimensional many-body problem Yang 1967, pp. 1312–1315.

With the declared index convention, the coefficient of i,j,k\lvert i,j,k\rangle on the left is

α,β,γSbcβγ(θ23)Saγαk(θ13)Sαβij(θ12),\sum_{\alpha,\beta,\gamma} S_{bc}^{\beta\gamma}(\theta_{23}) S_{a\gamma}^{\alpha k}(\theta_{13}) S_{\alpha\beta}^{ij}(\theta_{12}),

whereas the coefficient on the right is

α,β,γSabαβ(θ12)Sαciγ(θ13)Sβγjk(θ23).\sum_{\alpha,\beta,\gamma} S_{ab}^{\alpha\beta}(\theta_{12}) S_{\alpha c}^{i\gamma}(\theta_{13}) S_{\beta\gamma}^{jk}(\theta_{23}).

The order can be checked directly: on the left, the first collision is bcβγb\,c\to\beta\,\gamma, then aγαka\,\gamma\to\alpha\,k, then αβij\alpha\,\beta\to i\,j. On the right the first collision is abαβa\,b\to\alpha\,\beta. This verbal trace is a useful defense against silently reversing the operator product.

Let P(xy)=yxP(x\otimes y)=y\otimes x be the permutation operator on VVV\otimes V, and define

R(u)=uI+iP.R(u)=u\,I+iP.

The additive-parameter Yang–Baxter equation is

R12(u)R13(u+v)R23(v)=R23(v)R13(u+v)R12(u).\begin{aligned} R_{12}(u)R_{13}(u+v)R_{23}(v) ={}& R_{23}(v)R_{13}(u+v)R_{12}(u). \end{aligned}

To verify it, expand each side in powers of uu, vv, and ii. The identity, linear, and quadratic-permutation terms match trivially after using Pij2=IP_{ij}^2=I. The cubic and mixed terms reduce to the symmetric-group relations

P12P13=P23P12,P13P23=P12P13,P12P13P23=P23P13P12.P_{12}P_{13}=P_{23}P_{12}, \qquad P_{13}P_{23}=P_{12}P_{13}, \qquad P_{12}P_{13}P_{23}=P_{23}P_{13}P_{12}.

Each relation can also be checked on a basis tensor xyzx\otimes y\otimes z. For example, both sides of the last relation send it to zyxz\otimes y\otimes x when the rightmost factor acts first.

This RR matrix is an algebraic solution, not yet a relativistic two-body S matrix. A physical amplitude needs a scalar normalization and analytic structure compatible with unitarity, crossing, spectrum, and statistics.

The matrix equation constrains how internal labels are rearranged. It is logically independent of the following conditions.

  • Braiding unitarity: in a simple self-conjugate channel, S12(θ)S21(θ)=IS_{12}(\theta)S_{21}(-\theta)=I. Hermitian analyticity is a related but distinct reality statement.
  • Crossing: an amplitude with an incoming antiparticle is related to a crossed outgoing-particle amplitude through the charge-conjugation matrices and the shift θiπθ\theta\mapsto i\pi-\theta.
  • Analyticity: the physical sheet and physical strip must be declared; allowed poles and zeros have physical interpretations.
  • Spectrum and completeness: the representation content and stable particles entering the tensor spaces must be supplied. Yang–Baxter cannot reveal a species that was omitted from the ansatz.
  • Local existence: a consistent meromorphic matrix need not by itself construct local observables or a Hilbert-space QFT.

A symmetry-invariant ansatz often decomposes as

S(θ)=RSR(θ)PR,S(\theta)=\sum_R S_R(\theta)\,\mathcal P_R,

where PR\mathcal P_R projects VVV\otimes V onto irreducible representation RR. Yang–Baxter relates the channel functions SRS_R; a common scalar factor can remain undetermined. The exact S-matrix bootstrap then imposes the other analytic conditions and exposes the remaining CDD freedom.

The integrability exact-data chain places Yang–Baxter alongside, rather than above, these independent tests. The exact and rigorous status comparison keeps algebraic consistency distinct from construction.

For a doublet with no particle production, the simplest invariant ansatz is

Sabcd(θ)=A(θ)δacδbd+B(θ)δadδbc=A(θ)I+B(θ)P.S_{ab}^{cd}(\theta) =A(\theta)\delta_a^c\delta_b^d +B(\theta)\delta_a^d\delta_b^c =A(\theta)I+B(\theta)P.

Substituting into the Yang–Baxter equation and using the permutation algebra reduces the tensor equation to functional relations for AA and BB. The rational solution has B/A=i/uB/A=i/u up to an overall scalar factor when the spectral parameter uu is additive. The matrix ratio is constrained; the overall factor remains available for unitarity, crossing, pole data, and CDD choices.

Reversing the component order. State whether lower or upper indices are inputs and remember that the rightmost operator acts first. A correct-looking formula with the opposite convention can represent the transposed process.

Using Yang–Baxter as a substitute for crossing. The equations constrain different continuations and compositions. Both must be checked.

Calling every R matrix an S matrix. An R matrix can solve the algebraic equation without physical unitarity, relativistic crossing, a declared sheet, or a local-QFT realization.

  1. Starting from a,b,c\lvert a,b,c\rangle, reproduce the two displayed component expressions by applying the rightmost operator first.
Solution

For the left product,

a,b,cS23β,γSbcβγa,β,γS13α,β,γSbcβγSaγαkα,β,k,\lvert a,b,c\rangle \xrightarrow{S_{23}} \sum_{\beta,\gamma}S_{bc}^{\beta\gamma} \lvert a,\beta,\gamma\rangle \xrightarrow{S_{13}} \sum_{\alpha,\beta,\gamma} S_{bc}^{\beta\gamma}S_{a\gamma}^{\alpha k} \lvert\alpha,\beta,k\rangle ,

and S12S_{12} supplies SαβijS_{\alpha\beta}^{ij}. Applying S12S_{12}, then S13S_{13}, then S23S_{23} gives the right expression.

  1. Verify P12P13P23=P23P13P12P_{12}P_{13}P_{23}=P_{23}P_{13}P_{12} on xyzx\otimes y\otimes z.
Solution

On the left, P23P_{23} gives xzyx\otimes z\otimes y, P13P_{13} gives yzxy\otimes z\otimes x, and P12P_{12} gives zyxz\otimes y\otimes x. On the right, P12P_{12} gives yxzy\otimes x\otimes z, P13P_{13} gives zxyz\otimes x\otimes y, and P23P_{23} gives the same zyxz\otimes y\otimes x.

  • Yang, C. N. “Some Exact Results for the Many-Body Problem in One Dimension with Repulsive Delta-Function Interaction.” Physical Review Letters 19 (1967): 1312–1315. DOI.
  • Zamolodchikov, Alexander B., and Alexey B. Zamolodchikov. “Factorized S-Matrices in Two Dimensions as the Exact Solutions of Certain Relativistic Quantum Field Models.” Annals of Physics 120 (1979): 253–291. DOI.