Mandelstam Variables

Three numbers describe a two-body collision, one equation ties them together, and so only two are free. `aq-sr-mandelstam-invariants` builds those numbers, proves the constraint s+t+u=∑imi2c4s + t + u = \sum_i m_i^2c^4, and shows that in elastic scattering tt is never positive. This article picks up where that one deliberately stops. Two questions are left open there, and both have sharp answers. First: given ss, exactly which values of tt can occur? Not a vague "some range" — a closed interval whose width is

tmax⁡−tmin⁡  =  λ(s,m12c4,m22c4) λ(s,m32c4,m42c4)s,t_{\max} - t_{\min} \;=\; \frac{\sqrt{\lambda\big(s, m_1^2c^4, m_2^2c^4\big)}\ \sqrt{\lambda\big(s, m_3^2c^4, m_4^2c^4\big)}}{s},

a Källén function for the incoming pair times one for the outgoing pair. Second: what is the shape of the physical region? The usual picture is three straight-sided wedges. That picture is exactly right for equal masses and wrong in general, and the correct boundary is a cubic curve — the vanishing of the 3×33\times3 Gram determinant of the momenta, the same determinant that turned out to be the Källén function in `aq-sr-relativistic-collision-kinematics`, now one dimension up. And then the reason the apparatus exists. The same three invariants, read in three different regions of the same plane, describe three different reactions. That is crossing symmetry, and this article owns the kinematic half of it: the substitution rule, the three channels, and a proof that their physical regions never overlap.

Key takeaways
  • Convention lock. (+,−,−,−)(+,-,-,-); Pμ=(E, p⃗ c)P^\mu = (E,\ \vec p\,c); P⋅P=+m2c4P\cdot P = +m^2c^4. ss, tt, uu are in ENERGY², so s\sqrt{s} is an ENERGY and the sum rule carries c4c^4: s+t+u=∑imi2c4s+t+u = \sum_i m_i^2 c^4. `aq-sr-mandelstam-invariants` boxes it; this article recalls it and does not re-derive it.
  • Owned result 1 — the tt-range at fixed ss. tt is an affine function of cos⁡θ∗\cos\theta^{*}, so it sweeps a closed interval as cos⁡θ∗\cos\theta^{*} runs over [−1,+1][-1,+1], of width λ12λ34/s\sqrt{\lambda_{12}}\sqrt{\lambda_{34}}/s with λ12≡λ(s,m12c4,m22c4)\lambda_{12} \equiv \lambda(s, m_1^2c^4, m_2^2c^4).
  • Owned result 2 — the physical region is a Gram determinant (Kibble, 1960). With Gij=Pi⋅PjG_{ij} = P_i\cdot P_j for i,j∈{1,2,3}i,j \in \{1,2,3\}, the boundary is det⁡G=0\det G = 0 and the interior is det⁡G>0\det G > 0. In invariants, det⁡G\det G is a cubic in (s,t,u)(s,t,u).
  • And for equal masses it collapses. det⁡G=14 s t u\det G = \tfrac14\,s\,t\,u exactly. The boundary is stu=0stu = 0, the ss-channel region is {t≤0, u≤0}\{t \leq 0,\ u \leq 0\} — and s≥4m2c4s \geq 4m^2c^4 is then a consequence, not an extra condition.
  • The wedges are not straight in general. Machine-checked counterexample: π−p→K0Λ\pi^- p \to K^0\Lambda at s=1.700\sqrt{s} = 1.700 GeV has u=+0.189u = +0.189 GeV² at backward angles. u≤0u \leq 0 is an equal-mass accident.
  • Crossing symmetry, kinematically. One amplitude, three reactions: 1+2→3+41+2\to3+4, 1+3ˉ→2ˉ+41+\bar 3 \to \bar 2 + 4, 1+4ˉ→2ˉ+31+\bar 4 \to \bar 2 + 3. The rule is Pi→−PiP_i \to -P_i per crossed leg, and the three physical regions are provably disjoint.
  • Key idea: the sum rule makes (s,t,u)(s,t,u) a plane, the Gram determinant carves the physical region out of it, and crossing says the pieces it carves are three different experiments described by one function. The whole kinematics of two-body scattering is one cubic curve.
The core insight

The sum rule s+t+u=∑imi2c4s+t+u = \sum_i m_i^2c^4 says the three invariants are a point on a plane. That is one line, and it is where most treatments stop. But not every point on that plane is a collision: momenta must be on-shell, angles must be real, and the four four-vectors of a 2→22 \to 2 process must fit inside a three-dimensional subspace of spacetime, because conservation leaves only three independent and the scattering is planar. The condition for three four-vectors to span a subspace of the right causal type is exactly the sign of their Gram determinant — and expanded in ss, tt and uu, that determinant is a cubic. The physical region is where the cubic is positive; its boundary is where the momenta go collinear. Everything else — the tt-range, the three channels, their disjointness — reads off the same curve.

Where this lives today

Every measured cross-section at a hadron collider is reported as a function of ss and tt, for the reason this article proves: those are the invariants, and the physical region is a definite, computable set. When an experiment quotes dσ/dt\mathrm{d}\sigma/\mathrm{d}t it is quoting something that means the same thing in every frame. The Gram-determinant boundary is not a curiosity either. Kinematic fitting codes test candidate events against det⁡G≥0\det G \geq 0 to reject combinations no real collision could produce, and the same determinant — one row larger — draws the boundary of a Dalitz plot for three-body final states, which is how resonances in three-body decays are found at all. Crossing survives most robustly of the three. The dynamical statement — that a single analytic function M(s,t)\mathcal{M}(s,t) serves all three channels — is a theorem of quantum field theory, proved from locality and the spectral conditions, and it is why a calculation of Compton scattering doubles as a calculation of pair annihilation. The kinematics stops here. We have shown that the three reactions live in three disjoint regions of one plane, that the plane's physical part is bounded by a determinant, and that nothing in the invariants remembers which experiment you were doing. Why one function covers all three is a question for the field theory, and `aq-sr-mandelstam-invariants` is right to send it there.

The physical picture

A two-body reaction 1+2→3+41+2 \to 3+4 has four four-momenta, and conservation kills one: P4=P1+P2−P3P_4 = P_1 + P_2 - P_3. Three are left, each on its mass shell. What remains free is the relative geometry of three four-vectors, and after the sum rule that is exactly two numbers — take them to be ss and tt. So a collision is a point in a plane. The question this article is really about is: which points? Not all of them. Pick ss and tt at random and you will usually be describing something impossible — a scattering angle whose cosine is 1.41.4. The set of (s,t)(s,t) that can occur is the physical region, and its boundary is where the geometry degenerates: where the outgoing particle goes exactly forward or exactly backward, and the three independent four-momenta stop spanning a three-dimensional space. That is a Gram determinant condition, and it is the same idea `aq-sr-relativistic-collision-kinematics` used to explain the Källén function — there, two momenta and a 2×22\times2 determinant; here, three momenta and a 3×33\times3 one. Kinematic boundaries are always determinants going to zero, because a boundary is always a degeneracy. The last piece is what Mandelstam was actually after. The invariants do not know which reaction produced them. Feed them a point in a different region of the same plane and they describe 1+3ˉ→2ˉ+41 + \bar 3 \to \bar 2 + 4 — a different experiment, with different particles in the beam. Crossing symmetry is the statement that one function of (s,t,u)(s,t,u) serves all three, and the kinematics of it — which region is which, and why they never overlap — we can prove here with no quantum field theory at all.

The core equations

The invariants, recalled from the owner (`aq-sr-mandelstam-invariants`), with the metric re-derived from η\eta, not transcribed. For 1+2→3+41+2\to3+4, with Pμ=(E, p⃗ c)P^\mu = (E,\ \vec p\,c) so that P⋅P=ημνPμPν=E2−∣p⃗ ∣2c2=m2c4P\cdot P = \eta_{\mu\nu}P^\mu P^\nu = E^2 - |\vec p\,|^2c^2 = m^2c^4:

s=(P1+P2)2,t=(P1−P3)2,u=(P1−P4)2(1)s = (P_1+P_2)^2, \qquad t = (P_1-P_3)^2, \qquad u = (P_1-P_4)^2 \qquad (1)

where: all three are in energy². Physical reading: s\sqrt{s} is the total CM energy, and tt is the invariant square of the momentum transfer. The constraint, recalled and not re-derived — `aq-sr-mandelstam-invariants` proves it:

  s+t+u  =  ∑i=14mi2c4  (2)\boxed{\;s + t + u \;=\; \sum_{i=1}^{4} m_i^2c^4\;} \qquad (2)

where: note the c4c^4. Physical reading: three invariants, one equation, two free numbers — geometrically, any collision is a point on a plane. The general tt, at fixed ss, as a function of the CM scattering angle:

  t  =  m12c4+m32c4−2E1∗E3∗+2 p1∗p3∗c2cos⁡θ∗  (3)\boxed{\;t \;=\; m_1^2c^4 + m_3^2c^4 - 2E_1^{*}E_3^{*} + 2\,p_1^{*}p_3^{*}c^2\cos\theta^{*}\;} \qquad (3)

where: Ei∗E_i^{*} and pi∗p_i^{*} come from `aq-sr-centre-of-momentum-frame`: E1∗=(s+m12c4−m22c4)/(2s)E_1^{*} = (s + m_1^2c^4 - m_2^2c^4)/(2\sqrt{s}), p1∗=λ121/2/(2cs)p_1^{*} = \lambda^{1/2}_{12}/(2c\sqrt{s}), with λ12≡λ(s,m12c4,m22c4)\lambda_{12} \equiv \lambda(s, m_1^2c^4, m_2^2c^4), and 3→43\to4 for the outgoing pair. Physical reading: tt is affine in cos⁡θ∗\cos\theta^{*} — the general, unequal-mass statement, of which the elastic t=−2p2c2(1−cos⁡θ)t = -2p^2c^2(1-\cos\theta) boxed by `aq-sr-mandelstam-invariants` is the case m1=m3m_1 = m_3, m2=m4m_2 = m_4. Hence the first owned result — the width of the physical tt-interval:

  tmax⁡−tmin⁡  =  4 p1∗p3∗c2  =  λ12 λ34s  (4)\boxed{\;t_{\max} - t_{\min} \;=\; 4\,p_1^{*}p_3^{*}c^2 \;=\; \frac{\sqrt{\lambda_{12}}\ \sqrt{\lambda_{34}}}{s}\;} \qquad (4)

where: tmax⁡t_{\max} is at cos⁡θ∗=+1\cos\theta^{*} = +1 (forward), tmin⁡t_{\min} at −1-1 (backward). Physical reading: the entire allowed range of momentum transfer, from the four masses and ss alone. It vanishes when either pair is at its threshold — no phase space, nothing to transfer. And the second — the boundary of the physical region (Kibble, 1960). With Gij=Pi⋅PjG_{ij} = P_i\cdot P_j for the three independent momenta:

  det⁡G  ≥  0in the physical region,det⁡G=0  on its boundary  (5)\boxed{\;\det G \;\geq\; 0 \quad\text{in the physical region}, \qquad \det G = 0 \ \text{ on its boundary}\;} \qquad (5)

where: det⁡G\det G, in invariants, is a cubic in (s,t,u)(s,t,u). Physical reading: the boundary is where the three momenta stop spanning three dimensions — exactly forward or exactly backward scattering. For equal masses the cubic factorises completely:

  det⁡G  =  14 s t u(all mi=m)  (6)\boxed{\;\det G \;=\; \tfrac{1}{4}\,s\,t\,u \qquad (\text{all } m_i = m)\;} \qquad (6)

where: the sum rule has been used. Physical reading: the three straight-sided wedges of the textbook picture are the three components of stu≥0stu \geq 0 that run to infinity. In the ss-channel wedge s>0s>0, t≤0t\leq0, u≤0u\leq0, so stu≥0stu \geq 0 ✓ — and s≥4m2c4s \geq 4m^2c^4 follows from t,u≤0t,u \leq 0 and the sum rule rather than being imposed.

Figure
THE PHYSICAL REGION IS A GRAM DETERMINANT. The Mandelstam plane, drawn in the symmetric 120-degree coordinates in which s, t and u are the three signed distances to the three axes (the sum rule s + t + u = sum of m-squared c-to-the-fourth makes that possible). For EQUAL MASSES the Gram determinant collapses to det G = s t u / 4, so its zero set is the three straight lines s = 0, t = 0, u = 0, and the physical regions are the three shaded wedges running to infinity — one per channel, each with one invariant above threshold and the other two negative. The small central triangle also has stu greater than zero, but NO channel is above threshold there: that is the MANDELSTAM TRIANGLE, unphysical, and the domain in which the amplitude is real. For UNEQUAL masses the straight lines bend: the boundary is the full cubic curve det G = 0, and the region can poke into u greater than zero, as pi-minus p to K-zero Lambda actually does
u = 0 t = 0 s = 0 s-CHANNEL 1 + 2 → 3 + 4 s ≥ 4m²c⁴ · t ≤ 0 · u ≤ 0 u-CHANNEL 1 + 4̄ → 2̄ + 3 t-CHANNEL 1 + 3̄ → 2̄ + 4 Mandelstam triangle stu POSITIVE but NO channel physical det G = ¼ s t u (equal masses). Boundary = det G = 0 = the momenta going collinear.

The boundary, and the three channels

Derivation AThe t-range at fixed s, for arbitrary masses, and why its width is a product of two Kaellen functions

`aq-sr-mandelstam-invariants` boxes t=−2p2c2(1−cos⁡θ)t = -2p^2c^2(1-\cos\theta), which is exact and which is elastic: it assumes m1=m3m_1 = m_3 and m2=m4m_2 = m_4, so p1∗=p3∗≡pp_1^{*} = p_3^{*} \equiv p. Real experiments are usually not elastic. This derivation does the general case, and the payoff is a closed formula for how much momentum transfer a given s\sqrt{s} makes available.

Step 1. Work in the CM frame, where p⃗1 ∗=−p⃗2 ∗\vec p_1^{\,*} = -\vec p_2^{\,*}, p⃗3 ∗=−p⃗4 ∗\vec p_3^{\,*} = -\vec p_4^{\,*}, and θ∗\theta^{*} is the angle between p⃗1 ∗\vec p_1^{\,*} and p⃗3 ∗\vec p_3^{\,*}. Expand t=(P1−P3)⋅(P1−P3)t = (P_1 - P_3)\cdot(P_1-P_3):

t=m12c4+m32c4−2 P1 ⁣⋅ ⁣P3=m12c4+m32c4−2E1∗E3∗+2p1∗p3∗c2cos⁡θ∗.t = m_1^2c^4 + m_3^2c^4 - 2\,P_1\!\cdot\! P_3 = m_1^2c^4 + m_3^2c^4 - 2E_1^{*}E_3^{*} + 2p_1^{*}p_3^{*}c^2\cos\theta^{*} .

Why: <strong>only the last term contains the angle, and it enters linearly.</strong> tt is an affine function of cos⁡θ∗\cos\theta^{*}, so as cos⁡θ∗\cos\theta^{*} sweeps its physical range [−1,+1][-1,+1], tt sweeps a closed interval and nothing else.

Step 2. Read off the endpoints. Forward scattering (cos⁡θ∗=+1\cos\theta^{*} = +1) gives the largest tt; backward (cos⁡θ∗=−1\cos\theta^{*} = -1) the smallest:

tmax⁡,min⁡=m12c4+m32c4−2E1∗E3∗±2p1∗p3∗c2.t_{\max,\min} = m_1^2c^4 + m_3^2c^4 - 2E_1^{*}E_3^{*} \pm 2p_1^{*}p_3^{*}c^2 .

Why: no dynamics has entered. The amplitude decides how often each tt occurs; kinematics decides which are possible at all, and that is settled by the masses and ss.

Step 3. Take the difference and substitute the CM momenta from `aq-sr-centre-of-momentum-frame`, pi∗=λ1/2/(2cs)p_i^{*} = \lambda^{1/2}/(2c\sqrt{s}):

tmax⁡−tmin⁡=4 p1∗p3∗c2=4c2⋅λ122cs⋅λ342cs= λ12 λ34s t_{\max} - t_{\min} = 4\,p_1^{*}p_3^{*}c^2 = 4c^2\cdot\frac{\sqrt{\lambda_{12}}}{2c\sqrt{s}}\cdot\frac{\sqrt{\lambda_{34}}}{2c\sqrt{s}} = \boxed{\ \frac{\sqrt{\lambda_{12}}\,\sqrt{\lambda_{34}}}{s}\ }

with λ12=λ(s,m12c4,m22c4)\lambda_{12} = \lambda(s, m_1^2c^4, m_2^2c^4) and λ34=λ(s,m32c4,m42c4)\lambda_{34} = \lambda(s, m_3^2c^4, m_4^2c^4). Why: <strong>the incoming pair and the outgoing pair each contribute one Källén function, and the width is their geometric product over ss.</strong> Every factor of cc cancels, as it must. And by the Heron reading of `aq-sr-relativistic-collision-kinematics`, each λ\lambda is minus sixteen times a squared triangle area — so the width of the tt-range is set by two failed triangles.

Step 4. Sanity-check against the owner. For elastic scattering, λ12=λ34\lambda_{12} = \lambda_{34}, so the width is λ12/s=4p∗2c2\lambda_{12}/s = 4p^{*2}c^2; and tmax⁡=0t_{\max} = 0 because E1∗=E3∗E_1^{*} = E_3^{*} and p1∗=p3∗p_1^{*} = p_3^{*} make the constant term 2(m12c4−E1∗2+p1∗2c2)=02(m_1^2c^4 - E_1^{*2} + p_1^{*2}c^2) = 0. So tt runs over [−4p∗2c2, 0][-4p^{*2}c^2,\ 0] — exactly `aq-sr-mandelstam-invariants`' t=−2p2c2(1−cos⁡θ)t = -2p^2c^2(1-\cos\theta). Why: the general result must contain the special one, and it does. <strong>Equation (3) is the elastic formula with the assumption removed.</strong> *Verified in verify_aq-sr-mandelstam-variables.py: equation (3) is derived symbolically from t=(P1−P3)2t = (P_1-P_3)^2 in the CM frame; the width identity 4p1∗p3∗c2=λ12λ34/s4p_1^{*}p_3^{*}c^2 = \sqrt{\lambda_{12}\lambda_{34}}/s is proved in sympy; the elastic limit is checked to reproduce the owner's boxed formula exactly; and the π−p→K0Λ\pi^-p \to K^0\Lambda numbers printed below are computed at 50 digits.*

Go deeper · Derivation B · 2 minThe physical region is where a Gram determinant is positive, and for equal masses that determinant is s t u over four
Derivation BThe physical region is where a Gram determinant is positive, and for equal masses that determinant is s t u over four

The Mandelstam plane is usually drawn with three straight-sided wedges and a note that the picture is schematic. It is worth knowing what the exact statement is, because the exact statement is short, it is provable, and it is not three straight lines unless the masses are equal. The boundary is a cubic, and the cubic is a determinant. This is Kibble's theorem (Phys. Rev. 117, 1159, 1960), whose abstract puts it in one line: the physical regions "are bounded by a cubic curve in the plane of the two independent invariants". We derive it.

Step 1 — count. After conservation, three momenta are independent: take P1P_1, P2P_2, P3P_3. In the CM frame all four lie in the (ct,x,y)(ct, x, y) hyperplane — the collision is planar — so the three span at most a three-dimensional subspace of Minkowski space, of signature (+,−,−)(+,-,-). Why: this is the geometric content of "two-body scattering".

Step 2 — the determinant. Build Gij=Pi⋅PjG_{ij} = P_i\cdot P_j, i,j∈{1,2,3}i,j \in \{1,2,3\}. Every entry is ss, tt, uu or a mass:

G=(m12c412(s−m12c4−m22c4)12(m12c4+m32c4−t)⋅m22c412(m22c4+m32c4−u)⋅⋅m32c4)G = \begin{pmatrix} m_1^2c^4 & \tfrac12\big(s - m_1^2c^4 - m_2^2c^4\big) & \tfrac12\big(m_1^2c^4 + m_3^2c^4 - t\big) \\[3pt] \cdot & m_2^2c^4 & \tfrac12\big(m_2^2c^4 + m_3^2c^4 - u\big) \\[3pt] \cdot & \cdot & m_3^2c^4 \end{pmatrix}

Why: <strong>the Gram matrix of the process is a table of Mandelstam variables.</strong> P2⋅P3P_2\cdot P_3 is fixed by uu because u=(P2−P3)2u = (P_2-P_3)^2 as well as (P1−P4)2(P_1-P_4)^2 — the two forms agree by conservation.

Step 3 — the sign. A 3-plane of signature (+,−,−)(+,-,-) has det⁡(induced metric)=(+1)(−1)(−1)=+1>0\det(\text{induced metric}) = (+1)(-1)(-1) = +1 > 0. So for any real, physical configuration, det⁡G>0\det G > 0; and det⁡G=0\det G = 0 exactly when the three momenta become linearly dependent — which happens exactly when the scattering is collinear, cos⁡θ∗=±1\cos\theta^{*} = \pm1. Hence

  physical region={det⁡G≥0},boundary={det⁡G=0}={cos⁡θ∗=±1}  \boxed{\;\text{physical region} = \{\det G \geq 0\}, \qquad \text{boundary} = \{\det G = 0\} = \{\cos\theta^{*} = \pm 1\}\;}

Why: <strong>a boundary is a degeneracy, and a degeneracy is a determinant going to zero.</strong> The same principle made the Källén function a 2×22\times2 Gram determinant in `aq-sr-relativistic-collision-kinematics`. One more particle, one more row.

Step 4 — expand, and watch the textbook picture fall out. Written in invariants, det⁡G\det G is a cubic in (s,t,u)(s,t,u). Set all four masses equal to mm, impose the sum rule s+t+u=4m2c4s+t+u = 4m^2c^4, and it factorises completely:

  det⁡G=14 s t u(all masses equal)  \boxed{\;\det G = \tfrac14\,s\,t\,u \qquad (\text{all masses equal})\;}

So the boundary is stu=0stu = 0: the three straight lines of the textbook figure. In the ss-channel, t≤0t \leq 0 and u≤0u \leq 0, hence stu≥0stu \geq 0 ✓ — and, by the sum rule, s=4m2c4−t−u≥4m2c4s = 4m^2c^4 - t - u \geq 4m^2c^4. The threshold is not an extra condition. It is a theorem. Why: <strong>the equal-mass wedge is exact for a reason</strong> — the cubic degenerates into three lines. For unequal masses it does not, the boundary curves, and the region escapes the wedge: π−p→K0Λ\pi^-p \to K^0\Lambda at s=1.700\sqrt{s} = 1.700 GeV has u=+0.189 GeV2>0u = +0.189\ \mathrm{GeV}^2 > 0 at backward angles, which the straight wedge forbids. `aq-sr-mandelstam-invariants` draws that picture explicitly as a schematic and proves t≤0t \leq 0 for the elastic case; both statements survive intact. <strong>What does not survive is reading the schematic as the general boundary.</strong> *Verified in verify_aq-sr-mandelstam-variables.py: det⁡G\det G is built from the invariant entries and shown symbolically to be a cubic in (s,t,u)(s,t,u); the equal-mass factorisation det⁡G=stu/4\det G = stu/4 is proved after imposing the sum rule; det⁡G>0\det G > 0 is confirmed at physical interior points and det⁡G=0\det G = 0 at cos⁡θ∗=±1\cos\theta^{*} = \pm1; and the u=+0.189u = +0.189 GeV² counterexample is exhibited with its four PDG masses.*

Figure
t IS AFFINE IN cos-theta-star, SO THE PHYSICAL RANGE IS A CLOSED INTERVAL. The straight line is equation (3) for pi-minus p to K-zero Lambda at root-s = 1.700 GeV: t plotted against cos theta-star. At cos = +1 (forward) t reaches its maximum, minus 0.1077 GeV-squared; at cos = minus 1 (backward) its minimum, minus 0.6871 GeV-squared. The full width is 0.5794 GeV-squared, and that number equals root of lambda-12 times root of lambda-34, divided by s — a Kaellen function for the incoming pion-proton pair times one for the outgoing kaon-Lambda pair. NOTICE the dashed line: u changes sign inside the physical region, at cos theta-star = minus 0.3465, and is POSITIVE for everything more backward than 110.3 degrees. The straight-sided wedge picture says u cannot be positive here. The straight-sided wedge picture is the equal-mass special case
cos θ* = −1 cos θ* = +1 t t_max = −0.1077 GeV² (forward) t_min = −0.6871 GeV² (backward) width 0.5794 GeV² = √λ₁₂ √λ₃₄ / s u = 0 here cos θ* = −0.3465 · θ* = 110.3° u POSITIVE in here up to u = +0.189 GeV² π⁻p → K⁰Λ at √s = 1.700 GeV — and u is positive over a third of the angular range
Go deeper · Derivation C · 2 minCrossing symmetry as a substitution rule, and a proof that the three channels never overlap
Derivation CCrossing symmetry as a substitution rule, and a proof that the three channels never overlap

Here is the observation the invariants were invented for. Nothing in the definitions of ss, tt and uu says which particles were in the beam. Give the same three numbers to a different reaction and they still make sense — they just mean something else. That is crossing symmetry, and its kinematic half can be proved outright. Its dynamical half — that one analytic function M(s,t)\mathcal{M}(s,t) really does serve all three reactions — is quantum field theory, and `aq-sr-mandelstam-invariants` correctly defers it.

Step 1 — the rule. To cross a leg is to move a particle between the initial and final states, replacing it by its antiparticle and reversing its four-momentum:

particle i in ⟷ antiparticle ıˉ out,Pi ⟶ −Pi.\text{particle } i \text{ in} \ \longleftrightarrow\ \text{antiparticle } \bar\imath \text{ out}, \qquad P_i \ \longrightarrow\ -P_i .

Why: this is the Feynman–Stückelberg move of `aq-sr-reinterpretation-principle-causal-paradoxes`, applied to an external leg. For a <strong>timelike</strong> four-momentum that relabelling is legitimate in every frame — which is precisely why crossing is allowed here and the tachyonic version is not.

Step 2 — the three channels. Starting from 1+2→3+41+2\to3+4 and crossing one leg at a time:

channelreactionthe invariant that plays the role of CM energy²
ss1+2→3+41 + 2 \to 3 + 4s=(P1+P2)2s = (P_1+P_2)^2
tt1+3ˉ→2ˉ+41 + \bar 3 \to \bar 2 + 4t=(P1−P3)2→(P1+P3ˉ)2t = (P_1 - P_3)^2 \to (P_1 + P_{\bar 3})^2
uu1+4ˉ→2ˉ+31 + \bar 4 \to \bar 2 + 3u=(P1−P4)2→(P1+P4ˉ)2u = (P_1 - P_4)^2 \to (P_1 + P_{\bar 4})^2

Why: <strong>the minus signs in the definitions of tt and uu were the crossing all along.</strong> t=(P1−P3)2t = (P_1 - P_3)^2 becomes (P1+P3ˉ)2(P_1 + P_{\bar 3})^2 the moment you flip P3→−P3ˉP_3 \to -P_{\bar 3} — and then tt is manifestly the squared total energy of an incoming pair: positive, and able to sit on a resonance.

Step 3 — physicality. Each channel needs its own invariant above its own threshold: s≥(m1+m2)2c4s \geq (m_1+m_2)^2c^4, or t≥(m1+m3)2c4t \geq (m_1+m_3)^2c^4, or u≥(m1+m4)2c4u \geq (m_1+m_4)^2c^4. Why: `aq-sr-threshold-energy`'s condition, applied to whichever pair is incoming.

Step 4 — disjointness, proved. Take equal masses mm, so all three thresholds are 4m2c44m^2c^4. Suppose a point lay in the ss- and tt-channel regions at once. Then s≥4m2c4s \geq 4m^2c^4 and t≥4m2c4t \geq 4m^2c^4, so u=4m2c4−s−t≤−4m2c4<0u = 4m^2c^4 - s - t \leq -4m^2c^4 < 0. But the ss-channel region requires t≤0t \leq 0 (equation 6, and `aq-sr-mandelstam-invariants`), contradicting t≥4m2c4>0t \geq 4m^2c^4 > 0. The overlap is empty.

  the s-, t- and u-channel physical regions are pairwise disjoint  \boxed{\;\text{the } s\text{-, } t\text{- and } u\text{-channel physical regions are pairwise disjoint}\;}

Why: <strong>the sum rule does the work.</strong> Three quantities summing to a fixed small number cannot all be large, and each channel needs its invariant large and the other two negative. One collision is never two experiments — but one <strong>function</strong> is. *Verified in verify_aq-sr-mandelstam-variables.py: the crossing substitution P3→−P3P_3 \to -P_3 is applied symbolically and shown to turn tt into the squared CM energy of the crossed pair; the three thresholds are stated and the pairwise-disjointness argument is checked by exhaustive search over the equal-mass Mandelstam plane; and det⁡G\det G is confirmed non-negative in all three wedges and negative outside them.*

Limiting cases

CaseResultMeaning
cos⁡θ∗=±1\cos\theta^{*} = \pm1t=tmax⁡t = t_{\max} or tmin⁡t_{\min}, det⁡G=0\det G = 0collinear scattering: the Gram determinant degenerates. This is the boundary
elastic (m1=m3m_1{=}m_3, m2=m4m_2{=}m_4)tmax⁡=0t_{\max} = 0, width =4p∗2c2= 4p^{*2}c^2recovers `aq-sr-mandelstam-invariants`' t=−2p2c2(1−cos⁡θ)t = -2p^2c^2(1-\cos\theta) exactly
all masses equaldet⁡G=14stu\det G = \tfrac14 stuthe boundary degenerates into three straight lines; the wedges are exact
unequal massesdet⁡G=0\det G = 0 is a cubic curvethe wedges bend, and u>0u > 0 becomes possible inside the ss-channel region
s→(m3+m4)2c4s \to (m_3+m_4)^2c^4λ34→0\lambda_{34} \to 0, width →0\to 0at the production threshold there is no angle left to scatter through
all mi→0m_i \to 0s+t+u→0s+t+u \to 0, det⁡G→14stu\det G \to \tfrac14 stuthe massless limit — the form most often quoted

Worked examples

1 · Elastic pppp at the LHC: the forward peak occupies a billionth of the plane

`aq-sr-mandelstam-invariants` works elastic pppp scattering at s=13\sqrt{s} = 13 TeV and finds that t=−0.1t = -0.1 GeV² is a deflection of 48.7 μ48.7\ \murad. Recall that, then ask the question it does not: how much of the available tt-range is that?

Step 1 — recall. At s=13\sqrt{s} = 13 TeV, p∗c=6500.0p^{*}c = \mathbf{6500.0} GeV (exactly 6499.99993236499.9999323), and t=−0.1t = -0.1 GeV² is a deflection of 48.7 μrad\mathbf{48.7\ \mu rad} with a 53.353.3 MeV recoil. Owned by `aq-sr-mandelstam-invariants`; not re-derived.

Step 2 — the full range, from equation (4). Elastic, so the width is 4p∗2c2=s−4mp2c44p^{*2}c^2 = s - 4m_p^2c^4:

t∈[−(s−4mp2c4), 0]=[−1.6900×108, 0] GeV2.t \in \big[-(s - 4m_p^2c^4),\ 0\big] = \big[\mathbf{-1.6900\times10^{8}},\ 0\big]\ \mathrm{GeV^2}.

Step 3 — the boundary number.

∣t∣∣t∣max⁡=0.11.6900×108=5.917×10−10.\frac{|t|}{|t|_{\max}} = \frac{0.1}{1.6900\times10^{8}} = \mathbf{5.917\times10^{-10}} .

The entire elastic peak the LHC's forward detectors can see occupies about one billionth of the physically available tt-range. Everything the machine measures elastically is crushed against the t=0t = 0 wall.

Step 4 — check the geometry. The sum rule gives u=4mp2c4−s−t=−1.6900×108u = 4m_p^2c^4 - s - t = \mathbf{-1.6900\times10^{8}} GeV², and equation (6) gives

det⁡G=14 s t u=7.140×1014 GeV6  >  0 ✓\det G = \tfrac14\,s\,t\,u = \mathbf{7.140\times10^{14}}\ \mathrm{GeV^6} \;>\; 0 \ \checkmark

Positive, so the point is inside the physical region — and the smallness of ∣t∣|t| is exactly how close to the boundary it sits.

Step 5 — the moral. The elastic peak is not a small feature of a big region; it is a sliver against a wall. That is why the invariant cross-section is written dσ/dt\mathrm{d}\sigma/\mathrm{d}t and plotted logarithmically, and why forward physics needs detectors hundreds of metres down the beamline.

2 · π−p→K0Λ\pi^- p \to K^0\Lambda: an inelastic reaction, and a positive uu

Strangeness production at s=1.700\sqrt{s} = 1.700 GeV. Four different masses. Find the tt-range from equation (4), then find where uu changes sign.

Step 1 — identify. PDG masses, in MeV: mπ=139.570m_\pi = 139.570, mp=938.272m_p = 938.272, mK0=497.611m_{K^0} = 497.611, mΛ=1115.683m_\Lambda = 1115.683. The incoming pair needs s≥1077.84\sqrt{s} \geq 1077.84 MeV, the outgoing pair s≥1613.29\sqrt{s} \geq \mathbf{1613.29} MeV. At s=1700\sqrt{s} = 1700 MeV both clear, the second by only 8787 MeV.

Step 2 — the Källén functions and the CM momenta.

λ12=3.8922×1012,λ34=7.2050×1011  MeV4 ⟹ p1∗=580.25,  p3∗=249.65 MeV/c.\lambda_{12} = \mathbf{3.8922\times10^{12}},\qquad \lambda_{34} = \mathbf{7.2050\times10^{11}}\ \ \mathrm{MeV^4} \ \Longrightarrow\ p_1^{*} = \mathbf{580.25},\ \ p_3^{*} = \mathbf{249.65}\ \mathrm{MeV}/c .

The outgoing momentum is much the smaller — that is what "8787 MeV above threshold" buys you.

Step 3 — the tt-range, two ways. From equation (3) with E1∗=596.80E_1^{*} = 596.80 MeV and E3∗=556.73E_3^{*} = 556.73 MeV:

tmax⁡=−0.10769 GeV2 (forward),tmin⁡=−0.68714 GeV2 (backward).t_{\max} = \mathbf{-0.10769\ GeV^2}\ (\text{forward}), \qquad t_{\min} = \mathbf{-0.68714\ GeV^2}\ (\text{backward}).

Width =0.57945= \mathbf{0.57945} GeV². Check against equation (4): λ12λ34/s=0.57945\sqrt{\lambda_{12}}\sqrt{\lambda_{34}}/s = \mathbf{0.57945} GeV² ✓. And note tmax⁡≠0t_{\max} \neq 0: inelastic scattering cannot reach zero momentum transfer, because the incoming and outgoing momenta differ.

Step 4 — the boundary number. The sum rule gives u=∑imi2c4−s−tu = \sum_i m_i^2c^4 - s - t, with ∑imi2c4=2.3922\sum_i m_i^2c^4 = \mathbf{2.3922} GeV². At backward scattering:

u=2.3922−2.8900+0.68714=+0.18934 GeV2  >  0.u = 2.3922 - 2.8900 + 0.68714 = \mathbf{+0.18934\ GeV^2} \;>\; 0 .

uu is positive, and stays positive for every cos⁡θ∗<−0.3465\cos\theta^{*} < \mathbf{-0.3465} — everything more backward than θ∗=110.3∘\theta^{*} = \mathbf{110.3^\circ}, a third of the angular range.

Step 5 — the moral. The straight-sided wedge, which requires u≤0u \leq 0, has just been violated by an ordinary strong-interaction reaction. Nothing is wrong: the wedge is the equal-mass degeneration of the cubic det⁡G=0\det G = 0, and this process has four different masses. When the masses differ, use the determinant.

Physics of the portrait

Origins · 2 minCrossing symmetry and how this was discovered
CS
Crossing symmetry
1954–1958 · Gell-Mann, Goldberger and the analytic amplitude · Mandelstam's plane · and the variable that was never a variable

Crossing symmetry is one of the few ideas in physics that arrived as a consistency requirement before anyone knew what it was a requirement of. By the early 1950s it was clear that the amplitude for γ+p→γ+p\gamma + p \to \gamma + p and the amplitude for γ+γ→p+pˉ\gamma + \gamma \to p + \bar p were not independent functions — they were the same function, evaluated at different arguments. Murray Gell-Mann, Marvin Goldberger and Walter Thirring put the causality-and-analyticity machinery in place in Physical Review 95, 1612 (1954), and Gell-Mann and Goldberger made the crossing relation explicit for low-energy photon scattering in Physical Review 96, 1433 (1954). The content is startling on first meeting and inevitable on the second: an antiparticle in the final state is a particle in the initial state with reversed four-momentum, so the algebra cannot tell the two reactions apart — and neither, it turns out, can the dynamics. Stanley Mandelstam — born in Johannesburg on 12 December 1928, an engineering graduate of the Witwatersrand who took his physics doctorate at Birmingham in 1956 — supplied the coordinates in Physical Review 112, 1344 (1958). His purpose was not bookkeeping. He wanted to write the amplitude as an analytic function of two complex variables and reconstruct it from its singularities alone, and for that he needed a plane to be analytic on. The invariants gave him one, the sum rule made it two-dimensional, and the three physical regions became three boundary pieces of one analytic function's domain. The Mandelstam representation — a double dispersion relation over that plane — was the high-water mark of the S-matrix programme: particle physics from analyticity, unitarity and crossing, with no Lagrangian at all. The programme did not win; quantum field theory did. But the variables outlived it completely, and so did the picture: every collider paper still quotes s\sqrt{s}, every cross-section is still differential in tt, and the amplitude is still, when convenient, the same function in three different experiments.

History
1928
— Stanley Mandelstam is born in Johannesburg.
1954
— Gell-Mann, Goldberger and Thirring build the causality-to-analyticity machinery (Physical Review 95, 1612); Gell-Mann and Goldberger then make the crossing relation explicit for low-energy photon scattering (Physical Review 96, 1433). Crossing symmetry has a name.
1958
— Mandelstam introduces s, t and u, and the double dispersion relation over the plane they span, in Physical Review 112, 1344.
1960
— Kibble, at Caltech, shows the boundary of the physical region is a cubic curve in the plane of the two independent invariants (Physical Review 117, 1159). It is the vanishing of a Gram determinant of the momenta.
1973
— Byckling and Kajantie's Particle Kinematics collects the physical-region theory in the form still used, with Gram determinants throughout.
2016
— Mandelstam dies in Berkeley. The S-matrix programme he helped build is long superseded; the plane named after him is on the first page of every collider paper.

Going further

Texts. Byckling & Kajantie, Particle Kinematics (Wiley, 1973) — the definitive treatment of physical regions and Gram determinants · Eden, Landshoff, Olive & Polkinghorne, The Analytic S-Matrix (1966) · Halzen & Martin, Quarks and Leptons, ch. 4. Primary: M. Gell-Mann and M. L. Goldberger, Phys. Rev. 96, 1433 (1954); S. Mandelstam, Phys. Rev. 112, 1344 (1958); T. W. B. Kibble, Phys. Rev. 117, 1159 (1960).

Provenance. Signature (+,−,−,−)(+,-,-,-); Pμ=(E,p⃗ c)P^\mu = (E,\vec p\,c); ss, tt, uu in energy², so s+t+u=∑imi2c4s+t+u = \sum_i m_i^2 \mathbf{c^4} — the convention shipped by the owner, `aq-sr-mandelstam-invariants`, which owns the sum rule and derives it; this article recalls it and does not re-derive it, and likewise recalls its t≤0t \leq 0 result for the elastic case. Base defects corrected: eq[3] prints the constraint as ∑imi2c2\sum_i m_i^2 c^2 (momentum² units) — it should be c4c^4; and eq[5] writes the fixed-target invariant as s=(m1c)2+(m2c)2+2m2E1lab/cs = (m_1c)^2 + (m_2c)^2 + 2m_2E_1^{\rm lab}/c, whose final term carries a spurious 1/c1/c and does not match the units of the other two. Both are fixed here. The base's DERIV[0] (the Källén function) duplicates `aq-sr-centre-of-momentum-frame`'s owned payload and is recalled, not re-derived. All numbers above are asserted in verify_aq-sr-mandelstam-variables.py.

Builds on. `aq-sr-mandelstam-invariants` · `aq-sr-relativistic-collision-kinematics` · `aq-sr-center-of-momentum-invariant-mass` · `aq-sr-centre-of-momentum-frame`. Leads to. `aq-sr-compton-scattering`.

Why is the physical region not just the three wedges I have always seen drawn?

Because those wedges are the equal-mass case. The exact boundary is det⁡G=0\det G = 0, a cubic curve, and a cubic degenerates into three straight lines only when it factorises — which, after the sum rule, is exactly what happens when all four masses coincide, giving det⁡G=stu/4\det G = stu/4. With four different masses the cubic does not factorise, the sides bend, and consequences that looked structural turn out to be accidents. The sharpest: u≤0u \leq 0 in the ss-channel is an equal-mass accident. π−p→K0Λ\pi^-p \to K^0\Lambda at s=1.700\sqrt{s} = 1.700 GeV reaches u=+0.189u = +0.189 GeV², and nothing is wrong with it.

Does that contradict the claim that t is never positive?

No, and the distinction is worth being careful about. `aq-sr-mandelstam-invariants` proves t≤0t \leq 0 using two elastic identities — t=−2p2c2(1−cos⁡θ)t = -2p^2c^2(1-\cos\theta) and t=−2m2c2Trecoilt = -2m_2c^2T_{\rm recoil} — and both are exactly right in the setting they are proved in. This article adds the general formula, equation (3), which contains the elastic one as the case m1=m3m_1 = m_3, m2=m4m_2 = m_4. In our inelastic example tt is still negative everywhere (tmax⁡=−0.1077t_{\max} = -0.1077 GeV²) — it simply cannot reach zero. It is uu, not tt, whose sign is the equal-mass accident.

What is a Gram determinant doing in a scattering calculation?

Deciding whether the momenta fit. A 2→22 \to 2 process has three independent four-momenta, and they must span a subspace with the causal structure of a real collision — signature (+,−,−)(+,-,-), hence positive Gram determinant. When det⁡G\det G reaches zero the three vectors have become linearly dependent, which physically means the scattering has gone collinear: forward or backward, no angle left. The boundary of the physical region is a degeneracy, and degeneracies are determinants. The same idea one dimension down is why the Källén function of `aq-sr-relativistic-collision-kinematics` is a 2×22\times2 Gram determinant; one dimension up, it draws the Dalitz plot.

Is crossing symmetry a theorem or an observation?

Both, depending which half you mean. The kinematic half is a theorem, proved in this article: the substitution Pi→−PiP_i \to -P_i turns the definition of tt into the squared CM energy of a crossed pair, and the three channels are physical in three disjoint regions. The dynamical half — that one analytic function M(s,t)\mathcal{M}(s,t) takes the correct value in all three — is a theorem of quantum field theory, proved from locality and analyticity, and not something relativity alone can give you. `aq-sr-mandelstam-invariants` defers it explicitly, and so does this article.

If only two of s, t, u are independent, why carry the third?

For symmetry, and for exchange. The sum rule makes uu redundant as information but not as notation: with identical particles in the final state there is no fact of the matter about which one was "deflected", so the amplitude must be built symmetrically from tt and uu, and dropping uu gives the wrong answer. And in the crossing picture uu labels a genuine third physical process with its own wedge. The redundancy is the price of a description in which the three channels look alike.

What is the Mandelstam triangle in the middle of the plane?

The region where stu>0stu > 0 but no channel is above threshold — all three invariants positive and all three below 4m2c44m^2c^4. No real collision lives there, and that is what makes it valuable: with no physical process to generate an imaginary part, the amplitude is real and analytic throughout it. Mandelstam's double dispersion relation is built by starting in that triangle and continuing outwards to the three physical wedges. It is the one part of the plane where nothing happens, and it is the part the whole programme was anchored to.

Before this

Where this leads