Centre-of-Momentum Frame

In 1954 the Berkeley Radiation Laboratory finished a machine built around a single number. An antiproton must be made in pairs — baryon number forbids making one alone — so the cheapest reaction on hydrogen is p+p→p+p+p+pˉp + p \to p + p + p + \bar p, and the four particles coming out weigh 4mpc2=3753 MeV4m_pc^2 = 3753\ \mathrm{MeV}. The naive budget is obvious: supply 3753 MeV3753\ \mathrm{MeV} of beam and you can afford them. The Bevatron was designed for 6.2 GeV of kinetic energy, nearly four times that. It was not over-engineering. It was arithmetic. Fire a 6.2 GeV6.2\ \mathrm{GeV} proton at a proton sitting still and the lab ledger reads: beam 7.138 GeV7.138\ \mathrm{GeV}, target rest energy 0.938 GeV0.938\ \mathrm{GeV}, total 8.077 GeV8.077\ \mathrm{GeV}. Yet only 3.893 GeV3.893\ \mathrm{GeV} of it can become new matter. The other 4.183 GeV4.183\ \mathrm{GeV} is unspendable — not lost, not radiated, but locked into the forward motion the products must inherit, because the beam brought momentum and momentum does not evaporate. The frame that accounts for this honestly is the one where the total three-momentum vanishes: the centre-of-momentum frame. Its total energy is s\sqrt{s}, and s\sqrt{s} — not the beam energy, not the sum of lab energies — is what a collision has to spend. For the Bevatron, s=3.893 GeV\sqrt{s} = 3.893\ \mathrm{GeV} against a requirement of 3.753 GeV3.753\ \mathrm{GeV}: a margin of 140 MeV140\ \mathrm{MeV}. It cleared the bar by less than a sixth of a proton, and the antiproton appeared in 1955.

Key takeaways
  • Signature (+,−,−,−)(+,-,-,-): p⋅p=(E/c)2−∣p⃗ ∣2=+ m2c2p\cdot p = (E/c)^2 - |\vec p\,|^2 = +\,m^2c^2 on the mass shell. Timelike norms are positive.
  • Convention lock: ss carries units of ENERGY², so s\sqrt{s} is an ENERGY. s≡(∑iEi)2−∣∑ip⃗i∣2c2s \equiv (\sum_i E_i)^2 - |\sum_i \vec p_i|^2c^2, and s=13 TeV\sqrt{s} = 13\ \mathrm{TeV} at the LHC. The four-momentum norm P⋅PP\cdot P is instead in momentum² units, with s=(P⋅P)c2s = (P\cdot P)c^2. Both are standard; mixing them inside one calculation is a factor-of-c2c^2 disaster, and the AllQFT base for this article does exactly that.
  • The CM frame is the frame in which ∑ip⃗i=0⃗\sum_i \vec p_i = \vec 0. It exists — uniquely, up to rotations — for any system whose total four-momentum is timelike, and it moves at V⃗CM=c2P⃗tot/Etot\vec V_{\mathrm{CM}} = c^2\vec P_{\mathrm{tot}}/E_{\mathrm{tot}} relative to the lab.
  • There Pμ=(Mc,0⃗ )P^\mu = (Mc,\vec 0\,), so s=Mc2\sqrt{s} = Mc^2: the system's invariant mass, times c2c^2, IS its centre-of-momentum energy. Invariant Mass owns the concept; this article supplies the frame that makes it a reading.
  • Owned result — the available energy. Only s\sqrt{s} can become rest mass, because in every other frame the products must still carry the system's momentum. Fixed target on a stationary m2m_2: s=m12c4+m22c4+2m2c2E1lab\sqrt{s} = \sqrt{m_1^2c^4 + m_2^2c^4 + 2m_2c^2E_1^{\mathrm{lab}}}, growing only as E1\sqrt{E_1}. A head-on collider gives s=2Ebeam\sqrt{s} = 2E_{\mathrm{beam}}.
  • Key idea: the lab says how much energy is PRESENT; the CM frame says how much is USABLE, and the ratio is exactly 1/γCM1/\gamma_{\mathrm{CM}}. Matching the LHC's 13 TeV with a fixed target would need a 9.0 × 10⁴ TeV beam — ninety thousand TeV to buy thirteen.
The core insight

Momentum conservation is not a constraint applied at the end of a collision calculation; it is a tax levied at the start. A system with P⃗tot≠0⃗\vec P_{\mathrm{tot}} \neq \vec 0 must end with P⃗tot\vec P_{\mathrm{tot}}, so its products can never all be at rest, so some of their energy is permanently spoken for. Boost to the frame where P⃗tot=0⃗\vec P_{\mathrm{tot}} = \vec 0 and the tax vanishes — everything there is spendable. That frame is the CM frame, its energy is s\sqrt{s}, and since s\sqrt{s} is a Lorentz scalar you may compute it wherever is convenient.

Where this lives today

Every event at a modern collider is processed through this arithmetic before anyone looks at it. Detectors measure four-momenta in the lab; physics questions are asked of invariants. A candidate particle's mass is the s\sqrt{s} of its decay products. The "hard scatter" energy of an LHC event is the s\sqrt{s} of the two colliding partons — not 13 TeV, because each parton carries only a fraction xx of its proton's momentum: s^=x1x2 s\sqrt{\hat s} = \sqrt{x_1x_2}\,\sqrt{s}, with the hard collision's CM frame boosted along the beam by an unknown amount. That is exactly why hadron-collider physicists measure transverse momentum and use rapidity differences — those survive an unknown longitudinal boost. The same logic sets the design of every machine ever built. Fixed target: cheap, dense, high event rate — the right choice for statistics on physics you can already reach, and still the workhorse of neutrino production and nuclear physics. Collider: expensive, thin target, low rate — and the only choice when what you need is s\sqrt{s}.

The physical picture

Take two particles. Each carries a four-momentum piμ=(Ei/c, p⃗i )p_i^\mu = (E_i/c,\ \vec p_i\,), transforming with Λ\Lambda and with its norm fixed by the particle's mass, pi⋅pi=+mi2c2p_i\cdot p_i = +m_i^2c^2 (Relativistic Four-Momentum owns this; Conservation of Four-Momentum owns the statement that the sum is conserved). Add them:

Pμ=p1μ+p2μ=(E1+E2c, p⃗1+p⃗2)P^\mu = p_1^\mu + p_2^\mu = \left(\frac{E_1 + E_2}{c},\ \vec p_1 + \vec p_2\right)

This is a four-vector, because four-vectors add. But it is emphatically not on any particle's mass shell, and treating it as though it were is the source of half the errors in this subject. Its norm, however, is still an invariant — and it is the single most useful number in collision physics. Give it a name: s≡(P⋅P)c2s \equiv (P\cdot P)c^2. Now ask a question that sounds trivial and is not: is there a frame in which the total three-momentum is zero? For a single massive particle, yes — that is the rest frame. For a system, yes precisely when PμP^\mu is timelike, which is to say when s>0s > 0. And where it exists, the physics is simplest: two-body states become back-to-back, decays symmetric, momentum bookkeeping nil.

Figure
PUNCHLINE — the energy ledger of the Bevatron shot. In the lab, 8.077 GeV of energy is present; only 3.893 GeV of it is SPENDABLE. The remaining 4.183 GeV is locked into the forward motion the products are obliged to inherit from the beam's momentum. The antiproton needs 4·m_p·c² = 3.753 GeV — the machine cleared it by 140 MeV. Available fraction = 1/γ_CM = 48.2 percent
p + p collision, beam kinetic energy T = 6.2 GeV, proton target at rest LAB frame — total energy present E_tot = E_beam + m_p c² = 7.138 + 0.938 = 8.077 GeV but momentum conservation takes its cut … the same collision, honestly accounted √s = 3.893 GeV SPENDABLE — the CM energy 4.183 GeV LOCKED — kinetic energy of the CM's own motion (V = 0.876c) threshold to make an antiproton: 4 m_p c² = 3.753 GeV margin 140 MeV ✓ √s / E_tot = 1/γ_CM = 1/2.0746 = 0.482 — barely half the energy in the room is money you can spend. The naive budget (3.753 GeV of beam) would have missed the antiproton by a factor of nearly four.

That figure is the whole subject in one picture. Green is s\sqrt{s}; red is everything the lab has but cannot use. As the beam energy climbs, red grows faster, and the available fraction s/Etot=1/γCM\sqrt{s}/E_{\mathrm{tot}} = 1/\gamma_{\mathrm{CM}} collapses.

Figure
The existence proof, drawn. In the E–pc plane the two particles' four-momenta each sit on their own mass hyperbola (p·p = +m²c²). Their SUM P does not — it sits strictly INSIDE the light cone, at s = 15.155 GeV². A boost slides P along the hyperbola √s = 3.893 GeV until its momentum component is zero: that landing point, P = (√s, 0), IS the CM frame. The one system with no CM frame is a pair of collinear photons, whose sum lies ON the cone (s = 0) and can never be stood upright
E (GeV) p c (GeV) 02468 2468 light cone: E = pc, s = 0 proton mass shell, m_p c² = 0.938 GeV the invariant hyperbola √s = 3.893 GeV target p₂ = (0.938, 0) beam p₁ = (7.138, 7.076) P = p₁ + p₂ = (8.077, 7.076) — LAB not on any mass shell; strictly inside the cone P = (3.893, 0) — CM boost by V = 0.876c Existence test: is P timelike? s = E_tot² − |P|²c² = 8.077² − 7.076² = 15.155 GeV² s ≥ 0 ⇒ a CM frame exists, at V = c²P/E_tot ≤ c ✓ s = 0 (collinear photons) ⇒ NO CM frame exists ✗

The boost that stands PP upright is unique up to a rotation, which changes nothing. So "the" CM frame is a well-defined object, and every statement made in it is a Lorentz-invariant statement about the collision.

The core equations

Start from the total four-momentum of an isolated system, which is a four-vector because it is a sum of four-vectors, and conserved because each collision conserves it:

Pμ  =  ∑ipiμ  =  (Etotc, P⃗tot),P⃗tot=∑ip⃗i(1)P^{\mu} \;=\; \sum_i p_i^{\mu} \;=\; \left(\frac{E_{\mathrm{tot}}}{c},\ \vec P_{\mathrm{tot}}\right), \qquad \vec P_{\mathrm{tot}} = \sum_i \vec p_i \qquad (1)

where: EtotE_{\mathrm{tot}}, P⃗tot\vec P_{\mathrm{tot}} — the totals in whatever frame you stand in. Physical reading: PμP^\mu is the only object in a collision that both transforms cleanly and survives it. The CM frame kills its spatial part, and the boost that gets you there is read straight off PμP^\mu:

P⃗tot  CM=0⃗⟺V⃗CM=c2 P⃗totEtot,γCM=Etots(2)\vec P_{\mathrm{tot}}^{\;\mathrm{CM}} = \vec 0 \qquad\Longleftrightarrow\qquad \vec V_{\mathrm{CM}} = \frac{c^2\,\vec P_{\mathrm{tot}}}{E_{\mathrm{tot}}}, \qquad \gamma_{\mathrm{CM}} = \frac{E_{\mathrm{tot}}}{\sqrt{s}} \qquad (2)

where: V⃗CM\vec V_{\mathrm{CM}} — the CM frame's velocity seen from the lab. Physical reading: not a mass-weighted average of particle velocities (that object is not even Lorentz-covariant), but the velocity that makes the momentum vanish. And ∣V⃗CM∣<c|\vec V_{\mathrm{CM}}| < c requires ∣P⃗tot∣c<Etot|\vec P_{\mathrm{tot}}|c < E_{\mathrm{tot}} — exactly s>0s > 0. Now the convention this article — and the four after it — lives or dies by:

  s  ≡  (∑iEi)2−∣∑ip⃗i∣2c2  =  (P⋅P) c2,s  =  ECM  =  Mc2  (3)\boxed{\; s \;\equiv\; \Big(\textstyle\sum_i E_i\Big)^2 - \Big|\textstyle\sum_i \vec p_i\Big|^2 c^2 \;=\; (P\cdot P)\,c^2, \qquad \sqrt{s} \;=\; E_{\mathrm{CM}} \;=\; Mc^2 \;} \qquad (3)

where: MM — the invariant mass of the system (Invariant Mass owns the concept; this article supplies the frame in which it is just energy over c2c^2). Physical reading: ss has units of energy², so s\sqrt{s} is an energy — not a momentum. Write s≥∑fmfc2\sqrt{s} \geq \sum_f m_fc^2 with an ss defined as (E/c)2−∣p⃗ ∣2(E/c)^2 - |\vec p\,|^2 and you have compared a momentum to an energy. Both conventions exist (P⋅PP\cdot P is in momentum² units, s=(P⋅P)c2s = (P\cdot P)c^2 converts). Pick one, state it, never drift. Evaluate (3) in the lab for the standard arrangement — beam m1m_1 at lab energy E1E_1, target m2m_2 at rest:

s  =   m12c4+m22c4+2 m2c2 E1lab     →  E1≫mc2      2 m2c2 E1lab(4)\sqrt{s} \;=\; \sqrt{\,m_1^2c^4 + m_2^2c^4 + 2\,m_2c^2\,E_1^{\mathrm{lab}}\,} \;\;\xrightarrow[\;E_1 \gg mc^2\;]{}\;\; \sqrt{2\,m_2c^2\,E_1^{\mathrm{lab}}} \qquad (4)

where: E1labE_1^{\mathrm{lab}} — the beam's total energy, not its kinetic energy. Physical reading: the useful energy grows only as the square root of the beam energy. Quadruple the accelerator, double the reach. A head-on collider of equal beams already has P⃗tot=0⃗\vec P_{\mathrm{tot}} = \vec 0, so s=2Ebeam\sqrt{s} = 2E_{\mathrm{beam}} — linear. Threshold Energy owns that comparison and plots it. Finally, in the CM frame a two-body system's energies and momenta are fixed entirely by ss and the two masses:

E1∗=s+(m12−m22)c42s,p∗  =  [s−(m1+m2)2c4][s−(m1−m2)2c4]2cs(5)E_1^{*} = \frac{s + (m_1^2 - m_2^2)c^4}{2\sqrt{s}}, \qquad p^{*} \;=\; \frac{\sqrt{\big[s - (m_1+m_2)^2c^4\big]\big[s - (m_1-m_2)^2c^4\big]}}{2c\sqrt{s}} \qquad (5)

where: p∗≡∣p⃗1  ∗∣=∣p⃗2  ∗∣p^{*} \equiv |\vec p_1^{\;*}| = |\vec p_2^{\;*}| — the common momentum magnitude, equal and opposite by construction. Physical reading: the cc belongs in the denominator — the bracket product is energy⁴, its root energy², and dividing by csc\sqrt{s} gives energy/cc, a momentum. Particle Decay Kinematics sets s=Mc2\sqrt{s} = Mc^2 here and reads off the monoenergetic decay momentum.

Deriving the frame

Derivation AThe centre-of-momentum frame exists, is unique, and moves at c squared P over E

Everything below rests on a frame we have so far only asserted. Three things need proving: that a frame with ∑ip⃗i=0⃗\sum_i\vec p_i = \vec 0 exists at all; that its velocity relative to the lab is c2P⃗tot/Etotc^2\vec P_{\mathrm{tot}}/E_{\mathrm{tot}}; and that this velocity is subluminal, since a "frame" moving at cc is not a frame. The three turn out to be one statement, decided entirely by the sign of the invariant ss — which is the pattern of the whole subject: a question about frames, answered by a scalar.

Step 1. Boost along the momentum. Let the lab see EtotE_{\mathrm{tot}} and P⃗tot\vec P_{\mathrm{tot}}, and boost with velocity V⃗=Vn^\vec V = V\hat n along n^=P⃗tot/∣P⃗tot∣\hat n = \vec P_{\mathrm{tot}}/|\vec P_{\mathrm{tot}}|. Since PμP^\mu is a four-vector, its longitudinal momentum component transforms as P∥′=γV(P∥−VEtot/c2)P'_{\parallel} = \gamma_V\big(P_{\parallel} - VE_{\mathrm{tot}}/c^2\big). Why: a boost mixes the time component only with the momentum component along the boost direction, so boosting along P⃗tot\vec P_{\mathrm{tot}} leaves nothing transverse to kill.

Step 2. Demand the momentum vanish. Setting P∥′=0P'_{\parallel} = 0 gives P∥=VEtot/c2P_{\parallel} = VE_{\mathrm{tot}}/c^2, hence V⃗CM=c2P⃗tot/Etot\boxed{\vec V_{\mathrm{CM}} = c^2\vec P_{\mathrm{tot}}/E_{\mathrm{tot}}}. Why: γV\gamma_V is never zero, so the bracket must vanish. Note what this velocity is NOT: a mass-weighted average of constituent velocities. Momentum, not mass, is what gets cancelled.

Step 3. Check it is a legal velocity. ∣V⃗CM∣/c=∣P⃗tot∣c/Etot|\vec V_{\mathrm{CM}}|/c = |\vec P_{\mathrm{tot}}|c/E_{\mathrm{tot}}, which is below 11 if and only if ∣P⃗tot∣2c2<Etot2|\vec P_{\mathrm{tot}}|^2c^2 < E_{\mathrm{tot}}^2 — precisely s>0s > 0. Why: existence of the CM frame is EXACTLY the statement that PμP^\mu is timelike. Not an extra assumption but a property of the system, and an invariant one, so every observer agrees whether a CM frame exists.

Step 4. Find the exception. The sum of two future-pointing timelike-or-null four-vectors is timelike unless both are null and parallel. Two photons flying the same way have PμP^\mu null: s=0s = 0, no CM frame, no mass. Turn one around and s=4Eγ2>0s = 4E_\gamma^2 > 0 — a CM frame appears out of nowhere. Why: the sharpest illustration that invariant mass belongs to the SYSTEM and is not a sum over its parts. Two massless particles can make a massive system; two pointed the same way cannot.

Step 5. Uniqueness. Two frames both seeing P⃗=0⃗\vec P = \vec 0 differ by a Lorentz transformation preserving (Mc,0⃗ )(Mc,\vec 0\,) — and the stabiliser of a timelike vector is exactly the rotation group. So the CM frame is unique up to orientation; any boost away from it immediately reintroduces momentum. Why: this licenses the phrase "the CM frame", and lets E1∗E_1^{*} and p∗p^{*} count as invariants of the collision rather than accidents of a choice. *Verified in verify_aq-sr-centre-of-momentum-frame.py: boosting (Etot/c,P⃗tot)(E_{\mathrm{tot}}/c, \vec P_{\mathrm{tot}}) with V⃗=c2P⃗/E\vec V = c^2\vec P/E returns spatial part exactly zero and time component s/c\sqrt{s}/c; ∣V∣<c⇔s>0|V| < c \Leftrightarrow s > 0; γCM=Etot/s\gamma_{\mathrm{CM}} = E_{\mathrm{tot}}/\sqrt{s}; and the two-photon cases give s=0s = 0 (collinear) and s=4Eγ2s = 4E_\gamma^2 (anti-collinear).*

Go deeper · Derivation B · 2 minThe CM energy is the invariant mass, and the fixed-target formula that follows
Derivation BThe CM energy is the invariant mass, and the fixed-target formula that follows

The invariant ss is built from energies and momenta, both of which change under a boost; calling it invariant says that this particular combination does not. That turns ss from a definition into a tool: evaluate it wherever the algebra collapses, then read the answer wherever you work. We evaluate it twice — in the CM frame, which supplies the meaning of s\sqrt{s}, and in the lab, which supplies the number an experimenter needs.

Step 1. Evaluate ss in the CM frame. There ∑ip⃗i=0⃗\sum_i\vec p_i = \vec 0 by definition, so s=(∑iEiCM)2−0=ECM2s = \big(\sum_i E_i^{\mathrm{CM}}\big)^2 - 0 = E_{\mathrm{CM}}^2, giving s=ECM\sqrt{s} = E_{\mathrm{CM}}. Why: the momentum term dies. This is the only frame in which ss is the perfect square of something physically obvious — which is why it is the frame that gives ss its meaning.

Step 2. Identify the invariant mass. In the CM frame Pμ=(ECM/c, 0⃗ )≡(Mc, 0⃗ )P^\mu = (E_{\mathrm{CM}}/c,\ \vec 0\,) \equiv (Mc,\ \vec 0\,), which defines MM. So s=Mc2\boxed{\sqrt{s} = Mc^2}. Why: a four-vector with no spatial part looks exactly like the four-momentum of a single particle at rest, and MM is what that fictitious particle would weigh — a genuine, gravitating mass. But it is NOT ∑imi\sum_i m_i, and Relativistic Inelastic Collisions is built on the difference.

Step 3. Evaluate the same ss in the lab: beam (m1m_1, energy E1E_1, momentum p⃗1\vec p_1) on a stationary target (m2m_2). Expanding, s=(E1+m2c2)2−∣p⃗1∣2c2s = \big(E_1 + m_2c^2\big)^2 - |\vec p_1|^2c^2. Why: nothing is assumed about the collision — this is arithmetic on the initial state, legitimate because ss is conserved (the four-momentum is) and invariant (the norm is).

Step 4. Kill ∣p⃗1∣2c2|\vec p_1|^2c^2 with the mass shell, ∣p⃗1∣2c2=E12−m12c4|\vec p_1|^2c^2 = E_1^2 - m_1^2c^4. Substituting: s=E12+2m2c2E1+m22c4−E12+m12c4s = E_1^2 + 2m_2c^2E_1 + m_2^2c^4 - E_1^2 + m_1^2c^4, so s=m12c4+m22c4+2m2c2E1lab\boxed{s = m_1^2c^4 + m_2^2c^4 + 2m_2c^2E_1^{\mathrm{lab}}}. Why: the E12E_1^2 terms cancel identically. That cancellation is why ss is LINEAR in the beam energy — hence why s\sqrt{s} grows only as E1\sqrt{E_1}. The entire economics of accelerator physics sits in it.

Step 5. Check both ends. At E1=m1c2E_1 = m_1c^2 the formula collapses to s=(m1+m2)2c4s = (m_1+m_2)^2c^4, so s=(m1+m2)c2\sqrt{s} = (m_1+m_2)c^2 — two particles at rest side by side have invariant mass m1+m2m_1 + m_2, as they must. At E1≫mc2E_1 \gg mc^2 the mass terms drop and s→2m2c2E1\sqrt{s} \to \sqrt{2m_2c^2E_1}. Why: the low-energy limit is the ONLY place where invariant mass is additive; the high-energy limit is where the subject gets expensive.

Step 6. Contrast the collider. Two equal beams of energy EE head-on have P⃗tot=0⃗\vec P_{\mathrm{tot}} = \vec 0 in the lab — the lab already IS the CM frame — so s=(2E)2s = (2E)^2 and s=2E\sqrt{s} = 2E, with no square root eating the gain. Why: a collider does nothing clever with the energy. It simply refuses to pay the momentum tax, by arranging for there to be no momentum. *Verified in verify_aq-sr-centre-of-momentum-frame.py: ss is recomputed symbolically in both frames and shown equal; the fixed-target formula is derived from the mass shell; both limits are taken; and the Bevatron numbers printed in this article (E1=7.138E_1 = 7.138, s=3.893\sqrt{s} = 3.893, VCM=0.876cV_{\mathrm{CM}} = 0.876c, γCM=2.075\gamma_{\mathrm{CM}} = 2.075) are reproduced to the digits shown.*

Go deeper · Derivation C · 2 minSplitting the CM energy between two bodies
Derivation CSplitting the CM energy between two bodies

Once the CM frame is fixed, a two-body state has no freedom left but its orientation: the momenta are equal and opposite, and their common magnitude — along with both energies — is fixed entirely by ss and the two masses. Worth doing carefully, because the result is the workhorse of the next three articles, and because the standard form of it is where factors of cc go to die.

Step 1. Write the constraints. Energy: E1∗+E2∗=sE_1^{*} + E_2^{*} = \sqrt{s}. Momentum: p⃗1  ∗=−p⃗2  ∗\vec p_1^{\;*} = -\vec p_2^{\;*}, so ∣p⃗1  ∗∣=∣p⃗2  ∗∣≡p∗|\vec p_1^{\;*}| = |\vec p_2^{\;*}| \equiv p^{*}. Why: two equations, two unknowns. The problem is already solved in principle; the rest is algebra.

Step 2. Subtract the mass shells. E1∗2−E2∗2=(p∗2c2+m12c4)−(p∗2c2+m22c4)=(m12−m22)c4E_1^{*2} - E_2^{*2} = \big(p^{*2}c^2 + m_1^2c^4\big) - \big(p^{*2}c^2 + m_2^2c^4\big) = (m_1^2 - m_2^2)c^4, since p∗p^{*} is common and cancels. Why: this is the trick — subtracting the mass shells eliminates the momentum entirely and leaves a pure statement about energies.

Step 3. Solve the linear system. Factor: (E1∗+E2∗)(E1∗−E2∗)=s (E1∗−E2∗)(E_1^{*} + E_2^{*})(E_1^{*} - E_2^{*}) = \sqrt{s}\,(E_1^{*} - E_2^{*}), so E1∗−E2∗=(m12−m22)c4/sE_1^{*} - E_2^{*} = (m_1^2 - m_2^2)c^4/\sqrt{s}. Add E1∗+E2∗=sE_1^{*} + E_2^{*} = \sqrt{s} and halve: E1∗=s+(m12−m22)c42s\boxed{E_1^{*} = \dfrac{s + (m_1^2 - m_2^2)c^4}{2\sqrt{s}}}, with E2∗E_2^{*} following by swapping labels. Why: the difference of squares turns a quadratic system linear. Read the structure — each body gets half of s\sqrt{s}, plus or minus a term measuring which is heavier.

Step 4. Extract the common momentum. Insert E1∗E_1^{*} into p∗2c2=E1∗2−m12c4p^{*2}c^2 = E_1^{*2} - m_1^2c^4 and factor the resulting difference of squares twice. The Källén (triangle) function λ(x,y,z)=x2+y2+z2−2xy−2yz−2zx\lambda(x,y,z) = x^2+y^2+z^2-2xy-2yz-2zx appears, and it factorises: λ(s,m12c4,m22c4)=[s−(m1+m2)2c4][s−(m1−m2)2c4]\lambda\big(s, m_1^2c^4, m_2^2c^4\big) = \big[s - (m_1+m_2)^2c^4\big]\big[s - (m_1-m_2)^2c^4\big], giving p∗=λ1/2(s, m12c4, m22c4)2cs\boxed{p^{*} = \dfrac{\lambda^{1/2}\big(s,\,m_1^2c^4,\,m_2^2c^4\big)}{2c\sqrt{s}}}. Why: the factorised form beats the symmetric one because it makes the physics visible — p∗p^{*} vanishes exactly when s=(m1+m2)c2\sqrt{s} = (m_1+m_2)c^2, the threshold, and goes imaginary below it, which is how the algebra tells you the state is forbidden.

Step 5. Audit the units. λ\lambda takes energy² arguments, so λ\lambda is energy⁴ and λ1/2\lambda^{1/2} energy²; dividing by 2cs2c\sqrt{s} (a velocity times an energy) leaves energy/velocity — a momentum. Put the cc in the numerator instead and p∗p^{*} comes out as velocity × energy, wrong by c2c^2. Why: the AllQFT base writes exactly that, having defined ss in momentum² units in one card and used it in energy² units in the next four. Units are the cheapest error-detector you own. *Verified in verify_aq-sr-centre-of-momentum-frame.py: E1∗+E2∗=sE_1^{*} + E_2^{*} = \sqrt{s} and E1∗2−p∗2c2=m12c4E_1^{*2} - p^{*2}c^2 = m_1^2c^4 hold symbolically; the Källén factorisation is proved by expansion; p∗p^{*} is dimensionally asserted to be an energy over cc; and the equal-mass case gives E∗=s/2E^{*} = \sqrt{s}/2 and the printed Bevatron value p∗=1.705p^{*} = 1.705 GeV/cc.*

Limiting cases

CaseResultMeaning
P⃗tot=0⃗\vec P_{\mathrm{tot}} = \vec 0 alreadyV⃗CM=0⃗\vec V_{\mathrm{CM}} = \vec 0, s=Etot\sqrt{s} = E_{\mathrm{tot}}the lab is the CM frame — what a collider buys
all particles at rests=(∑imi)c2\sqrt{s} = \big(\sum_i m_i\big)c^2the only case where invariant mass is additive
E1≫mc2E_1 \gg mc^2 (fixed target)s→2m2c2E1\sqrt{s} \to \sqrt{2m_2c^2E_1}available energy ∝Ebeam\propto \sqrt{E_{\mathrm{beam}}} — the root colliders exist to defeat
head-on equal beamss=2Ebeam\sqrt{s} = 2E_{\mathrm{beam}}available energy ∝Ebeam\propto E_{\mathrm{beam}} — linear
v≪cv \ll c throughoutV⃗CM→∑miv⃗i/∑mi\vec V_{\mathrm{CM}} \to \sum m_i\vec v_i/\sum m_iNewton's centre-of-mass velocity, recovered because Ei→mic2E_i \to m_ic^2
s→(m1+m2)2c4s \to (m_1+m_2)^2c^4p∗→0p^{*} \to 0threshold: the two bodies are at rest relative to each other in the CM frame
two photons, back-to-back, EγE_\gamma eachs=4Eγ2>0s = 4E_\gamma^2 > 0, M=2Eγ/c2M = 2E_\gamma/c^2a system of two massless particles has mass
two photons, collinears=0s = 0, M=0M = 0no CM frame exists — PμP^\mu is null and no boost stands it upright

Worked examples

1 · The Bevatron: was 6.2 GeV enough?

A proton beam of kinetic energy T=6.2 GeVT = 6.2\ \mathrm{GeV} strikes protons at rest in a liquid-hydrogen target. Find the CM energy, the CM frame's velocity, and decide whether the reaction p+p→p+p+p+pˉp+p\to p+p+p+\bar p (which requires four protons' worth of rest mass in the final state) is possible. Take mpc2=938.272 MeVm_pc^2 = 938.272\ \mathrm{MeV}.

Step 1 — identify. "Is this reaction possible" is a question about s\sqrt{s} and nothing else. The lab energies matter only as inputs to ss.

Step 2 — the beam's four-momentum. E1=T+mpc2=6.200+0.938=7.138 GeVE_1 = T + m_pc^2 = 6.200 + 0.938 = \mathbf{7.138\ GeV}; from the mass shell, p1c=E12−mp2c4=50.955−0.880=7.076 GeVp_1c = \sqrt{E_1^2 - m_p^2c^4} = \sqrt{50.955 - 0.880} = \mathbf{7.076\ GeV}.

Step 3 — apply the invariant. s=2mp2c4+2mpc2E1=0.880+0.880+2(0.938)(7.138)=15.155 GeV2s = 2m_p^2c^4 + 2m_pc^2E_1 = 0.880 + 0.880 + 2(0.938)(7.138) = \mathbf{15.155\ GeV^2}, so s=3.893 GeV\sqrt{s} = \mathbf{3.893\ GeV}. Cross-check from the totals: Etot=8.077 GeVE_{\mathrm{tot}} = 8.077\ \mathrm{GeV}, ∣P⃗∣c=7.076 GeV|\vec P|c = 7.076\ \mathrm{GeV}, and 8.0772−7.0762=65.230−50.075=15.155 GeV28.077^2 - 7.076^2 = 65.230 - 50.075 = 15.155\ \mathrm{GeV^2}. Same number, two routes — that is what "invariant" buys.

Step 4 — the frame itself. VCM/c=7.076/8.077=0.876V_{\mathrm{CM}}/c = 7.076/8.077 = \mathbf{0.876} and γCM=Etot/s=8.077/3.893=2.075\gamma_{\mathrm{CM}} = E_{\mathrm{tot}}/\sqrt{s} = 8.077/3.893 = \mathbf{2.075}, agreeing with 1/1−0.87621/\sqrt{1-0.876^2} as it must. There each proton carries E∗=s/2=1.947 GeVE^{*} = \sqrt{s}/2 = 1.947\ \mathrm{GeV} and p∗c=1.9472−0.880=1.705 GeVp^{*}c = \sqrt{1.947^2 - 0.880} = 1.705\ \mathrm{GeV}, back to back.

Step 5 — the verdict. The final state needs ∑fmfc2=4mpc2=3.753 GeV\sum_f m_fc^2 = 4m_pc^2 = \mathbf{3.753\ GeV}, and s=3.893>3.753\sqrt{s} = 3.893 > 3.753. Yes — by 140 MeV. The 6.2 GeV was not a round number; it was chosen to sit above this threshold with room for Fermi motion in nuclear targets.

Step 6 — the boundary number. Of the 6.200 GeV6.200\ \mathrm{GeV} of kinetic energy poured in, the part that could become new rest mass is s−2mpc2=3.893−1.877=2.017 GeV\sqrt{s} - 2m_pc^2 = 3.893 - 1.877 = \mathbf{2.017\ GeV}, or 32.5 %. Nearly two-thirds is spent dragging the collision forward at 0.876c0.876c.

2 · What a fixed target would cost at LHC energies

The LHC collides protons head-on at 6.5 TeV6.5\ \mathrm{TeV} per beam. What is s\sqrt{s}? And what beam energy would a fixed-target machine need to match it?

Step 1 — the collider is easy. Head-on equal beams: P⃗tot=0⃗\vec P_{\mathrm{tot}} = \vec 0, the lab already IS the CM frame, and s=2Ebeam=13 TeV\sqrt{s} = 2E_{\mathrm{beam}} = \mathbf{13\ TeV}. Every joule is spendable.

Step 2 — invert the fixed-target formula. From s=m12c4+m22c4+2m2c2E1s = m_1^2c^4 + m_2^2c^4 + 2m_2c^2E_1 with m1=m2=mpm_1 = m_2 = m_p: E1=(s−2mp2c4)/(2mpc2)E_1 = \big(s - 2m_p^2c^4\big)/(2m_pc^2).

Step 3 — apply. s=(13 TeV)2=1.69×108 GeV2s = (13\ \mathrm{TeV})^2 = 1.69\times10^8\ \mathrm{GeV^2} and 2mpc2=1.8765 GeV2m_pc^2 = 1.8765\ \mathrm{GeV}; the 2mp2c4=1.76 GeV22m_p^2c^4 = 1.76\ \mathrm{GeV^2} correction is negligible against 1.69×1081.69\times10^8. So E1=1.69×108/1.8765=9.006×107 GeVE_1 = 1.69\times10^8/1.8765 = \mathbf{9.006\times10^7\ GeV}.

Step 4 — the boundary number. That is 9.0×104 TeV=90 PeV\mathbf{9.0\times10^4\ TeV} = 90\ \mathrm{PeV} — about 6,900 times the LHC's per-beam energy, for the same physics. At LHC magnet technology (27 km for 6.5 TeV) such a ring would be of order 105 km10^5\ \mathrm{km} around: several Earth diameters. This is why the collider was invented; Threshold Energy puts the whole curve on a plot.

Figure
Available energy versus beam energy for a proton on a stationary proton. The curve is √s = √(2 m_p c² E_beam + 2 m_p²c⁴) — a SQUARE ROOT. Quadrupling the beam doubles the reach. Marked: the Bevatron shot (7.14 GeV of beam buys 3.89 GeV of √s, just clearing the 3.753 GeV antiproton line) and the brutal high-energy points. A head-on collider instead gives √s = 2E_beam, a straight line off the top of this frame
√s (GeV) beam energy E₁ (GeV) 05101520 050100150200250 fixed target: √s = √(2 m_p c² E₁ + 2 m_p²c⁴) grows as √E₁ — a flattening curve collider: √s = 2E₁ — LINEAR leaves the frame at E₁ = 8.5 GeV antiproton line: √s = 4 m_p c² = 3.753 GeV BEVATRON E₁ = 7.14 GeV ⇒ √s = 3.89 GeV 50 GeV ⇒ 9.78 GeV 100 GeV ⇒ 13.76 GeV 200 GeV ⇒ 19.42 GeV Quadruple the beam, double the reach. 50 → 200 GeV of beam (4×) takes √s from 9.78 to 19.42 GeV (2×). To match the LHC's 13 TeV you would need 9.0 × 10⁴ TeV of beam.

Physics of the portrait

The frame has no discoverer, and that is the interesting thing about it. It arrives twice — in 1687 as a theorem about the centre of mass, and around 1908 as a statement about a four-vector — and the second arrival quietly destroys the first while keeping its answer.

Origins · 2 minCentre-of-momentum frame and how this was discovered
CF
Centre-of-momentum frame
1687–1955 · Newton's Corollary IV · Minkowski's four-vector · the invariant s · the Bevatron's design energy

Newton's Principia (1687) already contains the essential move: Corollary IV to the Laws states that the common centre of gravity of a system of bodies does not alter its state of motion or rest by the mutual actions of those bodies. That is momentum conservation, dressed as a statement about a point, and it licenses the trick of solving a collision in the frame where that point stands still. It works because in Newtonian mechanics a point and a frame are interchangeable. Relativity severs them: once simultaneity is frame-dependent, the average position ∑mir⃗i/∑mi\sum m_i\vec r_i/\sum m_i stops being observer-independent, because different frames evaluate the average at different sets of events and get genuinely different answers. What survives is not the centre of mass but the condition ∑ip⃗i=0⃗\sum_i\vec p_i = \vec 0, which mentions no positions and needs no simultaneity convention. Minkowski's 1908 four-vector formalism made that natural: a system's total four-momentum is one geometric object, and asking for the frame where its spatial part vanishes is just asking to look along it. Planck, Lewis and Tolman, and later Møller worked out what a relativistic "centre" could even mean; the answer — that the frame is well defined while the point is not — is why CM is best read as centre of momentum. The idea earned its keep in 1954, when the Bevatron was built to a beam energy chosen by exactly this calculation, and in 1955 the antiproton appeared where s\sqrt{s} said it would.

History
1687
— Newton's Principia, Corollary IV: the common centre of gravity of an isolated system does not change its state of motion — momentum conservation, and the licence to work in the frame where it is at rest.
1908
— Minkowski's "Raum und Zeit" gives the four-vector formalism in which a system's total energy and momentum are one geometric object, and the CM condition becomes a statement about direction.
1943
— Rolf Widerøe patents the idea of colliding two beams head-on rather than firing one at a stationary target, precisely to avoid paying the momentum tax.
1954
— the Bevatron at Berkeley is commissioned with a design proton energy of 6.2 GeV kinetic, set by the requirement that √s clear 4·m_p·c² = 3.753 GeV.
1955
— Chamberlain, Segrè, Wiegand and Ypsilantis observe the antiproton at the Bevatron (Physical Review 100, 947); Chamberlain and Segrè share the 1959 Nobel Prize.
1971
— CERN's Intersecting Storage Rings collide proton on proton for the first time on 27 January, reaching √s = 62 GeV — a fixed-target machine would have needed a 2000 GeV beam.

Going further

Texts. Rindler, Introduction to Special Relativity, 2nd ed., §31 · Griffiths, Introduction to Elementary Particles, 2nd ed., §3.3–3.5 · Hagedorn, Relativistic Kinematics (1963), ch. 4 · Byckling & Kajantie, Particle Kinematics (1973), ch. 1–4 · Particle Data Group, Review of Particle Physics, "Kinematics" · Chamberlain, Segrè, Wiegand & Ypsilantis, "Observation of Antiprotons", Phys. Rev. 100, 947 (1955).

Provenance. Signature (+,−,−,−)(+,-,-,-): p⋅p=+m2c2p\cdot p = +m^2c^2, s>0s > 0 for a timelike total four-momentum, and ss in energy² units so that s\sqrt{s} is an energy (equation 3). The AllQFT base is internally inconsistent by a factor of c2c^2 and its p∗p^{*} card puts cc in the numerator; both are corrected above. verify_aq-sr-centre-of-momentum-frame.py parses this article and re-derives every number printed in it — E1=7.138E_1 = 7.138, p1c=7.076p_1c = 7.076, Etot=8.077E_{\mathrm{tot}} = 8.077, s=15.155s = 15.155, s=3.893\sqrt{s} = 3.893, VCM=0.876cV_{\mathrm{CM}} = 0.876c, γCM=2.075\gamma_{\mathrm{CM}} = 2.075, E∗=1.947E^{*} = 1.947, p∗c=1.705p^{*}c = 1.705, 4mpc2=3.7534m_pc^2 = 3.753, margin 140 MeV140\ \mathrm{MeV}, 32.5%32.5\%, LHC equivalent 9.006×107 GeV9.006\times10^7\ \mathrm{GeV}.

Scope. Thresholds are stated here and owned by Threshold Energy, which derives s≥∑fmfc2\sqrt{s} \geq \sum_f m_fc^2 and plots the asymmetry. Mandelstam tt and uu are not introduced; this article uses ss only. Invariant mass is owned by Invariant Mass, four-momentum conservation by Conservation of Four-Momentum.

Builds on. Relativistic Four-Momentum · Conservation of Four-Momentum · Invariant Mass · Relativistic Energy. Leads to. Relativistic Elastic Collisions · Relativistic Inelastic Collisions · Particle Decay Kinematics · Threshold Energy · Mandelstam Invariants.

Is the CM frame the frame in which the centre of mass is at rest?

Non-relativistically yes, and that is where the name comes from. Relativistically the two come apart and only one survives. The Newtonian centre of mass is a point, R⃗=∑mir⃗i/∑mi\vec R = \sum m_i\vec r_i/\sum m_i, and computing it means evaluating all the positions at the same time — so it depends on a simultaneity convention, so different frames locate it at genuinely different places. The CM frame is defined by a condition on the total four-momentum alone, ∑ip⃗i=0⃗\sum_i\vec p_i = \vec 0: no positions, no simultaneity. "CM" survives as an abbreviation but should be read as centre of momentum; when a particle-physics paper says "centre-of-mass frame", the momentum condition is what it means.

Why is √s and not the total lab energy the available energy?

Because momentum is conserved and therefore cannot be spent. If the incoming system carries P⃗tot\vec P_{\mathrm{tot}}, whatever comes out must carry P⃗tot\vec P_{\mathrm{tot}} too, so the products cannot all be at rest — collectively they must retain enough kinetic energy to supply that momentum, and that part cannot become rest mass. In the CM frame the constraint vanishes: the products may all be at rest, and the whole of s\sqrt{s} is convertible. Since s\sqrt{s} is invariant you may compute it in the lab and still trust it as the budget. The unavailable part is exactly Etot−s=s (γCM−1)E_{\mathrm{tot}} - \sqrt{s} = \sqrt{s}\,(\gamma_{\mathrm{CM}} - 1) — the kinetic energy of the CM frame's own motion.

Can the invariant mass of a system be less than the sum of its parts' masses?

Yes — every bound system does it. A helium nucleus weighs less than two protons plus two neutrons, and the deficit is the binding energy. Nothing above forbids that: s=Mc2\sqrt{s} = Mc^2 is the total CM energy, and a bound system's negative potential energy is part of that total. What is forbidden is a system of free particles with M<∑imiM < \sum_i m_i, since there s=∑iEiCM≥∑imic2\sqrt{s} = \sum_i E_i^{\mathrm{CM}} \geq \sum_i m_ic^2, with equality only if every one is at rest. Binding Energy and Mass Defect develops the bound case; Relativistic Inelastic Collisions the free one, where MM strictly exceeds the sum.

The formula V = c²P/E looks like it could exceed c. Can it?

Not for any physical system. ∣V⃗CM∣/c=∣P⃗tot∣c/Etot|\vec V_{\mathrm{CM}}|/c = |\vec P_{\mathrm{tot}}|c/E_{\mathrm{tot}}, and since every particle obeys Ei≥∣p⃗i∣cE_i \geq |\vec p_i|c (equality only if massless), the triangle inequality gives Etot=∑Ei≥∑∣p⃗i∣c≥∣∑p⃗i∣c=∣P⃗tot∣cE_{\mathrm{tot}} = \sum E_i \geq \sum|\vec p_i|c \geq |\sum\vec p_i|c = |\vec P_{\mathrm{tot}}|c. So ∣V⃗CM∣≤c|\vec V_{\mathrm{CM}}| \leq c always, with equality only in the degenerate case where every particle is massless and they all point the same way — which is exactly the one configuration with no CM frame. The formula is safe precisely because nothing with mass in it can saturate the bound.

If s is conserved AND invariant, are those the same statement?

No, and keeping them apart matters. Invariance means all observers of the same state agree on ss — it is a Lorentz scalar. Conservation means the initial and final states have the same ss — it follows from ∑pμ\sum p^\mu being conserved. The two are independent: a quantity can be invariant but not conserved (a particle's rest mass, in an inelastic collision), or conserved but not invariant (the total energy, conserved in every frame but with a different value in each). ss is both — which is what lets you evaluate it in the lab before the collision and spend it as the CM budget after.

Why is the fixed-target √s formula only linear in the beam energy under the root?

Because of a cancellation worth seeing. Expanding s=(E1+m2c2)2−∣p⃗1∣2c2s = (E_1 + m_2c^2)^2 - |\vec p_1|^2c^2 produces an E12E_1^2 from the square — and the mass shell, ∣p⃗1∣2c2=E12−m12c4|\vec p_1|^2c^2 = E_1^2 - m_1^2c^4, produces a −E12-E_1^2. They cancel exactly, leaving only the cross term 2m2c2E12m_2c^2E_1. The cancellation is not an accident: at high energy the beam four-momentum is nearly null, and a null vector contributes nothing to a norm by itself. All of ss comes from the cross term between beam and target — which is why ss is proportional to the target mass, and why a heavier target buys more s\sqrt{s} only as m2\sqrt{m_2}.

Before this

Where this leads