Relativistic Beam Dynamics in Accelerators

A proton at the top of the LHC's ramp is travelling at 0.999 999 990 50.999\,999\,990\,5 of the speed of light. Give it 1%1\% more momentum — a colossal 6868 GeV, more energy than the entire Tevatron beam — and its speed increases by 1.91.9 parts in 101010^{10}. The velocity is frozen. And once you accept that, almost every design decision in a modern accelerator follows, because a machine cannot steer a beam by its speed if the speed will not move. It must steer it by geometry — by how long a path the particle takes. That substitution is the whole of relativistic beam dynamics. It is why a particle in a synchrotron has two different masses — γm\gamma m when you push it sideways and γ3m\gamma^3 m when you push it forwards. It is why an accelerator's bunches hold themselves together instead of smearing out. And it is why every proton machine has one energy at which it briefly goes blind — the transition energy, where the frozen-velocity term and the path-length term cancel exactly, longitudinal focusing vanishes, and the radio-frequency system must jump its phase by 180°180° in a few milliseconds or lose the beam.

Key takeaways
  • Signature (+,−,−,−)(+,-,-,-): ds2=c2dt2−dx2\mathrm{d}s^2 = c^2\mathrm{d}t^2 - \mathrm{d}x^2, p ⁣⋅ ⁣p=+m2c2p\!\cdot\!p = +m^2c^2, timelike positive.
  • Relativistic inertia is anisotropic: F=γm a+γ3m (β ⁣⋅ ⁣a)β\mathbf{F} = \gamma m\,\mathbf{a} + \gamma^3 m\,(\boldsymbol\beta\!\cdot\!\mathbf{a})\boldsymbol\beta. Push across the motion and the particle's inertia is γm\gamma m; push along it and the inertia is γ3m\gamma^3 m. At high γ\gamma it is γ2\gamma^2 times harder to change a particle's speed than its direction.
  • Consequence, and it is the article's spine: the velocity freezes. dββ=1γ2dpp\dfrac{\mathrm{d}\beta}{\beta} = \dfrac{1}{\gamma^2}\dfrac{\mathrm{d}p}{p}. At the LHC's γ=7247\gamma = 7247, a 1%1\% momentum kick moves the speed by 1.9×10−101.9\times10^{-10}.
  • Steering is therefore all rigidity: p[GeV/c]=0.299792 B[T] ρ[m]p[\mathrm{GeV}/c] = 0.299792\,B[\mathrm{T}]\,\rho[\mathrm{m}]. Focusing is all geometry: alternating quadrupoles, net-focusing by 1/f∗=L/f2>01/f^* = L/f^2 > 0, betatron oscillations obeying Hill's equation.
  • Owned result: the transition energy is where relativity hands over to geometry. The phase-slip factor η=αc−1/γ2\eta = \alpha_c - 1/\gamma^2 is a race between a geometric term (αc\alpha_c: more momentum ⟹ longer orbit ⟹ arrives late) and a purely relativistic term (1/γ21/\gamma^2: more momentum ⟹ faster ⟹ arrives early). They cancel at γt=1/αc\boxed{\gamma_t = 1/\sqrt{\alpha_c}}, where longitudinal focusing vanishes and the stable RF phase flips from ϕs\phi_s to π−ϕs\pi - \phi_s. The CERN PS (γt=6.1\gamma_t = 6.1) must cross it every cycle; the LHC (γt=53.8\gamma_t = 53.8, injecting at γ=480\gamma = 480) was designed never to.
  • Acceleration shrinks the beam for free: normalized emittance εN=βγ ε\varepsilon_N = \beta\gamma\,\varepsilon is invariant, so ε∝1/βγ\varepsilon \propto 1/\beta\gamma. The LHC's beam shrinks by γ\gamma-ratio 1515 from injection to collision with no cooling at all, reaching a 10 μm\mathbf{10\ \mu m} spot.
  • Key idea: at β→1\beta \to 1 you cannot control a beam through its velocity, because its velocity will not move. Everything an accelerator does, it does through path length and rigidity.
The core insight

Newton's second law has one mass in it. Relativity's has two — and they differ by γ2\gamma^2, which at the LHC is 5×1075\times10^{7}. A modern accelerator is best understood as a machine built entirely around that ratio: it gives up on controlling the beam's speed (hopeless), and instead controls where the beam is and how long its path is. Rigidity, betatron oscillation, phase stability and transition energy are four names for the same surrender.

Where this lives today

Transition crossing is still an operational scar. The CERN PS, injecting at γ=2.49\gamma = 2.49 and extracting at γ=28.7\gamma = 28.7, must pass γt=6.09\gamma_t = 6.09 every cycle, and the beam loss and emittance blow-up around that instant have been fought for six decades with γt\gamma_t-jump quadrupoles and RF gymnastics. Modern designs increasingly do what the LHC did: choose a lattice that puts the problem outside the machine's range.

And the spin gives you a thermometer. In a bending field a particle's spin precesses faster than its momentum by the anomalous factor, so the spin tune — extra precessions per turn — is ν=aγ\nu = a\gamma, with a=(g−2)/2a = (g-2)/2 a known constant. Since ν\nu is a pure number and aea_e is known to parts in 10910^{9}, measuring ν\nu measures γ\gamma, and hence the beam energy, absolutely. One unit of spin tune is exactly mec2/ae=440.65 MeVm_ec^2/a_e = \mathbf{440.65\ MeV} of electron beam energy. At LEP this resonant depolarization calibrated the beam energy to better than 11 MeV — precise enough that the measured energy visibly responded to the tides of the solid Earth and to leakage currents from the Geneva–Paris TGV line. It is this chapter's sharpest case of relativistic kinematics used as a metrological standard. The covariant spin-precession equation underwriting it is `aq-sr-bmt-equation`, and the kinematic rotation beneath that is `aq-sr-thomas-precession`; the same aγa\gamma physics, turned around to measure aa rather than γ\gamma, is `aq-sr-muon-g2-time-dilation`.

The physical picture

Take p=γmv\mathbf{p} = \gamma m\mathbf{v} and differentiate. The Lorentz factor depends on the speed, so the product rule produces two terms:

F=dpdt=γm a  +  γ3m (β ⁣⋅ ⁣a) β.\mathbf{F} = \frac{\mathrm{d}\mathbf{p}}{\mathrm{d}t} = \gamma m\,\mathbf{a} \;+\; \gamma^3 m\,(\boldsymbol\beta\!\cdot\!\mathbf{a})\,\boldsymbol\beta .

The second term lives entirely along the direction of motion. So if you push sideways — β ⁣⋅ ⁣a=0\boldsymbol\beta\!\cdot\!\mathbf{a} = 0 — it vanishes, and the particle responds with an effective mass of γm\gamma m. If you push forwards, both terms act along β\boldsymbol\beta and the effective mass is γ3m\gamma^3 m. Force and acceleration are not even parallel in general. This is not a curiosity. At the LHC's top energy the ratio γ2\gamma^2 is 5.3×1075.3\times10^{7}. Bending a 6.86.8 TeV proton through the full 2727 km ring is routine engineering with 8.338.33 T dipoles; speeding it up by any noticeable amount is essentially impossible. The machine spends eight years of civil engineering and 12321232 superconducting dipoles to change the beam's direction, and adds only a few MeV per turn to its speed — and the speed barely notices. The cleanest statement of the freeze is a logarithmic derivative. From β=pc/E\beta = pc/E and E2=p2c2+m2c4E^2 = p^2c^2 + m^2c^4:

dEE=β2 dpp⟹  dββ=1γ2 dpp  \frac{\mathrm{d}E}{E} = \beta^2\,\frac{\mathrm{d}p}{p} \qquad\Longrightarrow\qquad \boxed{\;\frac{\mathrm{d}\beta}{\beta} = \frac{1}{\gamma^2}\,\frac{\mathrm{d}p}{p}\;}

A 1%1\% momentum spread in an LHC bunch corresponds to a velocity spread of 1.9×10−101.9\times10^{-10}. Every proton in the beam, fast ones and slow ones alike, is travelling at the same speed to ten significant figures. The only thing that distinguishes them is the length of the path they take.

Figure
The whole article in one plot. The particle's TWO effective masses, as a function of γ. Transverse inertia (blue) is γm — push sideways and the particle responds like a body of mass γm. Longitudinal inertia (red) is γ³m — push forwards and it responds like a body γ² times heavier still. The gap between the two curves IS the design constraint of every accelerator: at the LHC's γ = 7247 the ratio is 5.3×10⁷, so bending the beam is easy and speeding it up is hopeless. Below, the consequence: the fractional velocity response dβ/β per unit dp/p collapses as 1/γ². At LHC injection (γ = 480) a 1% momentum kick moves the speed by 4.3×10⁻⁸; at top energy by 1.9×10⁻¹⁰; at LEP2's γ = 204,501 by 2.4×10⁻¹³. The velocity is frozen, and everything else in this article is the consequence
effective mass / m γ 110³10⁶10⁹ 11010²10³10⁴10⁵ TRANSVERSE: γm push sideways — bending is easy LONGITUDINAL: γ³m push forwards — γ² times harder LHC top: γ = 7247 ratio γ² = 5.3×10⁷ THE CONSEQUENCE — the velocity is frozen: dβ/β = (1/γ²) · dp/p a 1% momentum kick changes the SPEED by: LHC injection, 450 GeV p γ = 480 4.3 × 10⁻⁸ LHC collision, 6.8 TeV p γ = 7 247 1.9 × 10⁻¹⁰ LEP2, 104.5 GeV e⁻ γ = 204 501 2.4 × 10⁻¹³ ⟹ steer by PATH LENGTH, not by speed

The core equations

Relativistic inertia, resolved into its two components:

  F=dpdt=γm a+γ3m (β ⁣⋅ ⁣a) β  (1)\boxed{\;\mathbf{F} = \frac{\mathrm{d}\mathbf{p}}{\mathrm{d}t} = \gamma m\,\mathbf{a} + \gamma^3 m\,(\boldsymbol\beta\!\cdot\!\mathbf{a})\,\boldsymbol\beta\;} \qquad (1)

where: a=dv/dt\mathbf{a} = \mathrm{d}\mathbf{v}/\mathrm{d}t is the coordinate acceleration. Physical reading: transverse inertia γm\gamma m, longitudinal inertia γ3m\gamma^3 m. Force and acceleration are not parallel unless the push is purely along or purely across the motion. The frozen velocity, which everything below rests on:

  dββ=1γ2 dpp  ,dEE=β2 dpp(2)\boxed{\;\frac{\mathrm{d}\beta}{\beta} = \frac{1}{\gamma^{2}}\,\frac{\mathrm{d}p}{p}\;}, \qquad \frac{\mathrm{d}E}{E} = \beta^{2}\,\frac{\mathrm{d}p}{p} \qquad (2)

where: exact, at all speeds. Physical reading: this is equation (1)'s longitudinal statement, in the currency an accelerator actually uses. As γ→∞\gamma \to \infty, momentum and energy track each other perfectly while the velocity refuses to move at all. Magnetic rigidity — how you steer something whose speed you cannot change:

  Bρ=pq  ⟹p [GeV/c]=0.299792  B [T]  ρ [m](3)\boxed{\;B\rho = \frac{p}{q}\;} \qquad\Longrightarrow\qquad p\,[\mathrm{GeV}/c] = 0.299792\;B\,[\mathrm{T}]\;\rho\,[\mathrm{m}] \qquad (3)

where: the coefficient is c×10−9c \times 10^{-9}, not a fitted constant. Physical reading: rigidity depends on momentum alone — not on energy, not on mass, not on speed. It is validated against a real machine in worked example 1. Transverse motion: Hill's equation, and its Courant–Snyder solution.

x′′(s)+K(s) x(s)=0,Kx=1ρ2−k,Ky=+k(4)x''(s) + K(s)\,x(s) = 0,\qquad K_x = \frac{1}{\rho^2} - k,\quad K_y = +k \qquad (4)
x(s)=ε β(s)  cos⁡ ⁣(ψ(s)+ϕ0),ψ(s)=∫ ⁣dsβ(s)(5)x(s) = \sqrt{\varepsilon\,\beta(s)}\;\cos\!\big(\psi(s) + \phi_0\big),\qquad \psi(s) = \int\!\frac{\mathrm{d}s}{\beta(s)} \qquad (5)

where: ss is arc length (the independent variable — not time), k=(q/p) ∂By/∂xk = (q/p)\,\partial B_y/\partial x is the normalized quadrupole gradient, and β(s)\beta(s) is the Twiss beta function, an envelope, entirely unrelated to v/cv/c. Physical reading: the beam does not travel down the centre of the pipe; it oscillates about it, and εβ(s)\sqrt{\varepsilon\beta(s)} is how far out it swings. The relativity is hidden inside k∝1/pk \propto 1/p: a stiffer beam is harder to focus, in exactly the way it is harder to bend. Emittance and adiabatic damping:

  εN=βγ ε=invariant under acceleration  ⟹ε∝1βγ(6)\boxed{\;\varepsilon_N = \beta\gamma\,\varepsilon = \text{invariant under acceleration}\;} \qquad\Longrightarrow\qquad \varepsilon \propto \frac{1}{\beta\gamma} \qquad (6)

where: ε\varepsilon is the geometric emittance (phase-space area over π\pi), εN\varepsilon_N the normalized emittance. Physical reading: acceleration cools the beam for free. Adding longitudinal momentum without adding transverse momentum tilts every trajectory closer to the axis, so the beam's angular spread shrinks as 1/βγ1/\beta\gamma while the invariant εN\varepsilon_N is untouched. And the longitudinal master equation — the phase-slip factor and the transition energy:

  η  =  αc−1γ2  ,  γt=1αc  (7)\boxed{\;\eta \;=\; \alpha_c - \frac{1}{\gamma^{2}}\;},\qquad \boxed{\;\gamma_t = \frac{1}{\sqrt{\alpha_c}}\;} \qquad (7)

where: αc≡(dC/C)/(dp/p)\alpha_c \equiv (\mathrm{d}C/C)/(\mathrm{d}p/p) is the momentum compaction — a pure number fixed by the lattice geometry, typically small and positive. Physical reading: η\eta decides whether a higher-momentum particle arrives early or late. The two terms fight, and the article's owned result is what happens when they draw.

Why a bunch holds together — and the one energy at which it does not

Derivation AThe phase-slip factor and the transition energy where relativity takes over from geometry

A synchrotron does not accelerate a continuous stream; it accelerates bunches, and it holds each bunch together with a feedback loop so stable that a proton can survive 10810^{8} turns. The loop works by making a particle's arrival time depend on its momentum — so that a particle with too much energy is punished with a badly-timed kick. This card derives the quantity that decides the sign of that punishment, shows it is a race between a geometric effect and a relativistic one, and locates the unique energy at which the race is a draw and the whole feedback loop switches off.

Step 1. Ask how the revolution time depends on momentum. With CC the orbit circumference and βc\beta c the speed, T=C/(βc)T = C/(\beta c), so

dTT=dCC−dββ.\frac{\mathrm{d}T}{T} = \frac{\mathrm{d}C}{C} - \frac{\mathrm{d}\beta}{\beta} .

Why: two, and only two, things can make a particle late — a longer path, or a slower speed. The whole of longitudinal dynamics is the competition between those two terms, and nothing else enters.

Step 2. Evaluate the geometric term. A higher-momentum particle is stiffer, bends less, and swings out to a longer orbit. Define the momentum compaction by exactly that:

dCC≡αc dpp,αc>0.\frac{\mathrm{d}C}{C} \equiv \alpha_c\,\frac{\mathrm{d}p}{p},\qquad \alpha_c > 0 .

Why: αc\alpha_c is pure lattice geometry — where the dipoles are, how the dispersion runs. It contains no relativity whatever, and for a typical ring it is small: 3.45×10−43.45\times10^{-4} for the LHC, 0.0270.027 for the CERN PS.

Step 3. Evaluate the relativistic term — and this is the one that carries the physics. From equation (2),

dββ=1γ2 dpp.\frac{\mathrm{d}\beta}{\beta} = \frac{1}{\gamma^2}\,\frac{\mathrm{d}p}{p} .

Why: <strong>this is the frozen velocity, and it is the entire relativistic content of the argument.</strong> Non-relativistically, γ→1\gamma \to 1 and this term is dp/p\mathrm{d}p/p at full strength — more momentum means proportionally more speed. Relativistically it is throttled by 1/γ21/\gamma^2 and dies. The particle's speed stops responding while its orbit keeps growing.

Step 4. Subtract, and name the result.

dTT=(αc−1γ2)dpp  ≡  η dpp,  η=αc−1γ2  \frac{\mathrm{d}T}{T} = \left(\alpha_c - \frac{1}{\gamma^2}\right)\frac{\mathrm{d}p}{p} \;\equiv\; \eta\,\frac{\mathrm{d}p}{p},\qquad \boxed{\;\eta = \alpha_c - \frac{1}{\gamma^2}\;}

Why: read the sign. <strong>Below transition</strong> (γ\gamma small, 1/γ2>αc1/\gamma^2 > \alpha_c, η<0\eta < 0): the speed-up wins, and a higher-momentum particle arrives <strong>early</strong>. <strong>Above transition</strong> (γ\gamma large, 1/γ2→01/\gamma^2 \to 0, η→αc>0\eta \to \alpha_c > 0): the longer path wins, and a higher-momentum particle arrives <strong>late</strong>. The particle's response to being given energy <strong>reverses sign</strong>, purely because of relativity.

Step 5. Find the draw, and see what breaks. Setting η=0\eta = 0:

  γt=1αc  \boxed{\;\gamma_t = \frac{1}{\sqrt{\alpha_c}}\;}

At this energy, and only here, a momentum error produces no timing error at all. The RF cavity can no longer tell an off-energy particle from an on-energy one, longitudinal focusing vanishes, the bunch stops being held together, and the synchrotron oscillation frequency Ωs∝η\Omega_s \propto \sqrt{\eta} goes to zero. Crossing it requires the stable phase to jump from ϕs\phi_s to π−ϕs\pi - \phi_s — a 180°180° flip of the entire RF system, executed in milliseconds while the beam is briefly unfocused. Why: <strong>this is the article&#39;s owned result and it is a genuinely relativistic accident.</strong> αc\alpha_c is geometry; 1/γ21/\gamma^2 is relativity; and γt\gamma_t is the energy at which relativity finishes taking over from geometry. In a Newtonian world 1/γ2≡1>αc1/\gamma^2 \equiv 1 > \alpha_c always, η\eta would be negative for every machine ever built, and <strong>transition energy would not exist</strong>.

Step 6. Check it against real machines. The CERN PS has αc=0.027\alpha_c = 0.027, hence γt=6.09\gamma_t = 6.09; it injects at γ=2.49\gamma = 2.49 and extracts at γ=28.7\gamma = 28.7, so it crosses transition on every single cycle, and has done since 1959 — the phase jump is routine operations, and the beam loss around it is a permanent nuisance. The LHC has αc=3.45×10−4\alpha_c = 3.45\times10^{-4}, hence γt=53.8\gamma_t = 53.8; it injects at γ=480\gamma = 480, already far above. Why: <strong>that was a design choice, not luck.</strong> The LHC never crosses transition because a lattice was chosen to put γt\gamma_t below the injection energy. Equation (7) is not a curiosity to be admired — it is a constraint that shapes the machine. *Verified in verify_aq-sr-relativistic-beam-dynamics.py: sympy derives (dβ/β)/(dp/p)=1/γ2(\mathrm{d}\beta/\beta)/(\mathrm{d}p/p) = 1/\gamma^2 and (dE/E)/(dp/p)=β2(\mathrm{d}E/E)/(\mathrm{d}p/p) = \beta^2 directly from E2=p2c2+m2c4E^2 = p^2c^2 + m^2c^4; η=αc−1/γ2\eta = \alpha_c - 1/\gamma^2 and γt=αc−1/2\gamma_t = \alpha_c^{-1/2} are solved symbolically; the sign of η\eta is asserted negative below γt\gamma_t and positive above; and the PS/LHC transition arithmetic is recomputed at 50 dps and grepped against this article.*

Go deeper · Derivation B · 2 minWhy alternating a focusing lens with a defocusing one focuses
Derivation BWhy alternating a focusing lens with a defocusing one focuses

A quadrupole magnet cannot focus a beam. Maxwell forbids it: a field that converges a beam in xx necessarily diverges it in yy, and the sum of the two focusing strengths is exactly zero. For twenty years this looked like a hard limit, and it forced early synchrotrons to use weak, gentle bending fields and beam pipes the size of a person. The escape — found by Christofilos in 1950 and published by Courant, Livingston and Snyder in 1952 — is that a sequence of alternating lenses can focus even though every individual lens is neutral. The goal here is to prove it in three lines of matrix algebra, and to see exactly where the free lunch comes from.

Step 1. Write the thin-lens transfer matrices, acting on the vector (x,x′)(x, x'):

MF=(10−1/f1),MD=(10+1/f1),ML=(1L01).M_F = \begin{pmatrix} 1 & 0 \\ -1/f & 1\end{pmatrix},\qquad M_D = \begin{pmatrix} 1 & 0 \\ +1/f & 1\end{pmatrix},\qquad M_L = \begin{pmatrix} 1 & L \\ 0 & 1\end{pmatrix}.

Why: MFM_F kicks a ray back toward the axis in proportion to its displacement; MDM_D kicks it away, by exactly as much. They are inverses. Multiplied together with nothing between them they give the identity — no focusing, as Maxwell demands.

Step 2. Put a drift between them and multiply, defocusing last:

M=MD ML MF=(1−L/fL−L/f21+L/f).M = M_D\,M_L\,M_F = \begin{pmatrix} 1 - L/f & L \\ -L/f^2 & 1 + L/f\end{pmatrix}.

Why: the order matters, and the drift is the whole trick. It is what stops the two lenses from cancelling.

Step 3. Read the net focal length off the lower-left element, which is by definition −1/f∗-1/f^*:

  1f∗=Lf2  >  0  \boxed{\;\frac{1}{f^{*}} = \frac{L}{f^{2}} \;>\; 0\;}

Net focusing — and, remarkably, whichever lens comes first, since L/f2L/f^2 is even in ff. Why: here is the mechanism, and it is pure geometry. The focusing lens bends the ray <strong>toward the axis</strong>. It therefore arrives at the defocusing lens <strong>closer in</strong>, where a quadrupole's kick — proportional to displacement — is <strong>weaker</strong>. The defocusing lens pushes back, but from a shorter lever arm. It cannot undo what was done. The net is convergent, second-order in L/fL/f, and it costs nothing but drift space.

Step 4. Note what it bought. A strong-focusing lattice confines the beam to millimetres instead of tens of centimetres, so the magnet aperture — and with it the iron, the power and the cost — collapses. Every synchrotron built since 1952 uses it. Why: this is why the transverse half of the subject is written in β(s)\beta(s) and ε\varepsilon rather than in force laws. Once the lattice alternates, the motion is a driven oscillation and the useful variables are amplitude and phase. *Verified in verify_aq-sr-relativistic-beam-dynamics.py: sympy multiplies the three matrices symbolically, confirms M21=−L/f2M_{21} = -L/f^2 exactly, confirms det⁡M=1\det M = 1 (symplectic, as Liouville requires), confirms the result is unchanged when MFM_F and MDM_D are swapped, and confirms that with no drift the two lenses cancel to the identity — the drift is the entire mechanism.*

Figure
THE OWNED RESULT. The phase-slip factor η = α_c − 1/γ², plotted against γ, is a RACE between two curves. The blue line is the geometric term α_c — pure lattice, no relativity, a constant: more momentum means a longer orbit, so the particle arrives LATE. The red curve is the relativistic term 1/γ², which starts at 1 (Newton) and dies as the velocity freezes: more momentum means more speed, so the particle arrives EARLY. Below transition the red curve wins and η is negative. Above it the blue line wins and η is positive. They cross at γ_t = 1/√α_c, where η = 0, the RF can no longer tell an off-momentum particle from an on-momentum one, longitudinal focusing VANISHES, and the stable phase must jump by 180°. The CERN PS (α_c = 0.027, γ_t = 6.09) injects at γ = 2.49 and extracts at γ = 28.7 — it crosses transition on every cycle. The LHC (α_c = 3.45×10⁻⁴, γ_t = 53.8) injects at γ = 480, already far above, and never crosses at all. And note the Newtonian limit at the far left: at γ = 1 the red curve sits at 1, above α_c for ANY real lattice, so η is negative always — TRANSITION ENERGY DOES NOT EXIST IN NEWTONIAN PHYSICS. It is made entirely of the frozen velocity
γ 110⁻¹10⁻²10⁻³ 11010²10³10⁴10⁵ size of each term GEOMETRY: α_c = 3.45×10⁻⁴ (LHC) constant. more momentum ⟹ longer orbit ⟹ arrives LATE RELATIVITY: 1/γ² dies as the velocity freezes more momentum ⟹ faster ⟹ arrives EARLY η = 0 γ_t = 1/√α_c = 53.8 NO FOCUSING. RF must jump 180° γ = 1 (Newton): 1/γ² = 1 > α_c always ⟹ η < 0 ALWAYS. no transition energy exists η = α_c − 1/γ² below γ_t : η < 0 — high momentum arrives EARLY above γ_t : η > 0 — high momentum arrives LATE CERN PS: γ_t = 6.09, runs γ = 2.49 → 28.7 — CROSSES LHC: γ_t = 53.8, injects at γ = 480 — NEVER CROSSES LHC inject γ = 480

Limiting cases

CaseResultMeaning
γ→1\gamma \to 1F→ma\mathbf{F} \to m\mathbf{a}, both inertias →m\to mNewton. Longitudinal and transverse mass coincide; there is no anisotropy to exploit
γ→1\gamma \to 1 in (7)η→αc−1\eta \to \alpha_c - 1for any real lattice αc≪1\alpha_c \ll 1, so η<0\eta < 0 always. Transition energy does not exist in a Newtonian world
γ≫1\gamma \gg 1η→αc\eta \to \alpha_cthe velocity is fully frozen; timing is pure geometry. All high-energy rings live here
γ=γt\gamma = \gamma_tη=0\eta = 0, Ωs=0\Omega_s = 0no longitudinal focusing. The bunch is unheld; the RF must jump phase by 180°180°
γ→∞\gamma \to \inftyε=εN/βγ→0\varepsilon = \varepsilon_N/\beta\gamma \to 0adiabatic damping: the beam shrinks for free, without cooling
c→∞c \to \inftyγ→1\gamma \to 1, dβ/β→dp/p\mathrm{d}\beta/\beta \to \mathrm{d}p/pspeed tracks momentum perfectly. Every effect in this article vanishes

Worked examples

1 · Steer a 6.8 TeV proton — and check the formula against the real machine first

The LHC's dipoles produce 8.338.33 T. What bending radius do they need for a 77 TeV proton, and does the ring fit in the 26.726.7 km LEP tunnel it inherited?

Step 1 — identify. This is rigidity, equation (3): p=0.299792 Bρp = 0.299792\,B\rho, with the coefficient c×10−9c \times 10^{-9} and no γ\gamma anywhere. Momentum is all a magnet sees.

Step 2 — validate the formula before trusting it. The LHC's actual bending radius is ρ=2804\rho = 2804 m. Then

p=0.299792×8.33×2804=7.00 TeV/c,p = 0.299792 \times 8.33 \times 2804 = \mathbf{7.00\ TeV}/c ,

which is the design energy to three figures. The formula reproduces a real machine, so we may use it.

Step 3 — check the geometry. The dipoles occupy only the arcs, so 2πρ=2π(2804)=17.62\pi\rho = 2\pi(2804) = 17.6 km of bending, inside a 26.726.7 km ring — the remaining 9.19.1 km is straight sections for the RF, the focusing lattice and the four experiments. The bending radius is not the ring radius, and confusing them is a standard error.

Step 4 — the boundary number. Run 3 operates at 6.86.8 TeV, not 7.07.0, because the dipoles were never fully trained; from (3), 6.86.8 TeV needs B=8.09B = 8.09 T. The energy of the LHC is set by one number: how much field its magnets will hold without quenching. Everything else — the tunnel, the optics, the physics reach of equation (6) in `aq-sr-collider-fixed-target-energetics` — followed from that.

2 · Free cooling: why the LHC beam is 23 times narrower at collision than at injection

The LHC injects protons at 450450 GeV with a normalized emittance εN=2.5 μ\varepsilon_N = 2.5\ \mum, and collides them at 6.86.8 TeV. No cooling system is applied. How much does the beam shrink, and what is the spot size at the interaction point?

Step 1 — identify. Equation (6): εN=βγ ε\varepsilon_N = \beta\gamma\,\varepsilon is invariant, so the geometric emittance falls as 1/βγ1/\beta\gamma. With β≈1\beta \approx 1 throughout, ε∝1/γ\varepsilon \propto 1/\gamma.

Step 2 — the two Lorentz factors. γinj=450/0.9383=480\gamma_{\text{inj}} = 450/0.9383 = 480 and γcoll=6800/0.9383=7247\gamma_{\text{coll}} = 6800/0.9383 = \mathbf{7247}. The ratio is 15.1\mathbf{15.1}: the geometric emittance shrinks by a factor of 1515 during the ramp, and nothing was done to it.

Step 3 — the emittances. εinj=2.5×10−6/480=5.21×10−9\varepsilon_{\text{inj}} = 2.5\times10^{-6}/480 = 5.21\times10^{-9} m·rad; εcoll=2.5×10−6/7247=3.45×10−10\varepsilon_{\text{coll}} = 2.5\times10^{-6}/7247 = \mathbf{3.45\times10^{-10}} m·rad.

Step 4 — the spot. The beam radius is σ=ε β(s)\sigma = \sqrt{\varepsilon\,\beta(s)} from (5). At injection the optics are relaxed, β≈11\beta \approx 11 m, giving σ=5.21×10−9×11=240 μm\sigma = \sqrt{5.21\times10^{-9} \times 11} = \mathbf{240\ \mu m}. At the collision point the optics are squeezed to β∗=0.30\beta^* = 0.30 m, giving

σ∗=3.45×10−10×0.30=10.2 μm.\sigma^{*} = \sqrt{3.45\times10^{-10} \times 0.30} = \mathbf{10.2\ \mu m} .

Step 5 — the boundary number. A 10 μ10\ \mum spot — a fifth the width of a human hair — carrying 6.86.8 TeV per proton, held for hours. And the total narrowing of 23.5\mathbf{23.5} in radius came from two sources only: acceleration (free, relativistic, factor 15.1=3.89\sqrt{15.1} = \mathbf{3.89}) and optics (paid for, factor 11/0.3=6.06\sqrt{11/0.3} = \mathbf{6.06}). Relativity did more than a third of the work and sent no bill.

Physics of the portrait

Origins · 2 minEmittance and how this was discovered
E
Emittance
1950–1958 · Christofilos invents strong focusing and does not publish · Courant, Livingston and Snyder publish it · Courant and Snyder give the theory its language

The quantity that named this article is the one nobody set out to invent. Emittance is the area a beam occupies in phase space — position against angle — and it matters because Liouville's theorem says you cannot shrink it without dissipation. It is the beam's irreducible untidiness, and it is what a machine designer spends their life fighting. The fight was transformed in 1952. Ernest Courant, M. Stanley Livingston and Hartland Snyder published The Strong-Focusing Synchrotron (Phys. Rev. 88, 1190 (1952)), showing that alternating quadrupole gradients confine a beam far more tightly than the weak, gentle focusing then in use. The magnet apertures collapsed, and with them the cost of a high-energy machine. The idea had in fact been found two years earlier by Nicholas Christofilos, a Greek-American elevator engineer working alone, who filed a patent application on 10 March 1950 and never published a paper; the Brookhaven group acknowledged his priority when they learned of it. Then Courant and Snyder gave it a language. Theory of the Alternating-Gradient Synchrotron (Annals of Physics 3, 1 (1958)) is one of the most-used papers in accelerator physics: it introduced the Twiss parameters, the betatron amplitude function β(s)\beta(s), the phase advance, the tune, and the Courant–Snyder invariant that is the emittance. Every accelerator built since has been designed in their variables — which is why the beam envelope in a 2727 km superconducting proton ring is still described by a function called β\beta, defined in a paper about a machine a thousand times weaker.

History
1928
— Rolf Widerøe's linear resonance accelerator; the multi-pass principle that all later machines share.
1944
— Vladimir Veksler publishes phase stability (Compt. Rend. Acad. Sci. U.R.S.S. 43, 444); the RF bucket becomes possible.
1945
— Edwin McMillan independently proposes the synchrotron on the same principle (Phys. Rev. 68, 143).
1950
— Nicholas Christofilos conceives alternating-gradient focusing and files a patent application; he does not publish.
1952
— Courant, Livingston and Snyder publish the strong-focusing synchrotron (Phys. Rev. 88, 1190).
1958
— Courant and Snyder's Theory of the Alternating-Gradient Synchrotron (Annals of Physics 3, 1) gives the field its permanent language: Twiss parameters, beta function, tune, emittance.
1959
— The CERN Proton Synchrotron begins operating at 26 GeV; with a transition gamma of 6.1 it must cross transition on every cycle, and has done ever since.
1989
— LEP begins operation; resonant depolarization later calibrates its beam energy through the spin tune, and proves sensitive enough to detect Earth tides and a nearby railway.
2008
— The LHC is commissioned in LEP's tunnel; its lattice puts the transition gamma at 53.8, below the injection gamma of 480, so it never crosses transition at all.

Going further

Texts. Wiedemann, Particle Accelerator Physics, 4th ed. · S.Y. Lee, Accelerator Physics, 4th ed. · Sands, The Physics of Electron Storage Rings, SLAC-121 (1970) · Courant & Snyder, Ann. Phys. 3, 1 (1958) · Courant, Livingston & Snyder, Phys. Rev. 88, 1190 (1952) · Chao & Tigner (eds.), Handbook of Accelerator Physics and Engineering.

Provenance. Signature (+,−,−,−)(+,-,-,-). Every number is recomputed at 50 dps in verify_aq-sr-relativistic-beam-dynamics.py, which reads this file and asserts each printed value against what it computed. Base status: the base's eq[6] (η=αc−1/γ2\eta = \alpha_c - 1/\gamma^2, γt=1/αc\gamma_t = 1/\sqrt{\alpha_c}), its FODO result 1/f∗=L/f21/f^* = L/f^2, its rigidity coefficient 0.2997920.299792, and its spin tune ν=aγ\nu = a\gamma were each re-derived symbolically and are correct. Its eq[8], the practical synchrotron-loss formula U0[keV]=88.5 E4[GeV]/ρ[m]U_0[\mathrm{keV}] = 88.5\,E^4[\mathrm{GeV}]/\rho[\mathrm{m}], was independently recomputed from first principles as 88.46\mathbf{88.46} and is correct — and it is the anchor that convicts the neighbouring base: `aq-sr-synchrotron-radiation`'s per-turn loss formula disagrees with it by a factor of cc. That defect is repaired in that article, not this one.

Scope. This article owns the relativistic dynamics of a beam: the inertia anisotropy, the frozen velocity, rigidity, the phase-slip factor and the transition energy. The radiation the beam emits while being bent — its γ4\gamma^4 scaling, its per-turn cost, and why it killed LEP — is `aq-sr-synchrotron-radiation`. The energetics of what you do with the beam once you have it is `aq-sr-collider-fixed-target-energetics`. The covariant spin-precession equation behind the spin tune is `aq-sr-bmt-equation`, and the kinematic rotation underneath it is `aq-sr-thomas-precession`; both are used here only through the single result ν=aγ\nu = a\gamma, which is quoted, not derived.

Builds on. Relativistic Momentum · Relativistic Force · The Lorentz Force. Leads to. Synchrotron Radiation · Muon Storage-Ring Time Dilation · Collider vs Fixed-Target Energetics.

Why is longitudinal inertia γ³m but transverse inertia only γm?

Because changing a particle's speed changes γ\gamma, and changing its direction does not. Differentiate p=γmv\mathbf{p} = \gamma m \mathbf{v}: the product rule gives γ˙mv+γmv˙\dot{\gamma} m\mathbf{v} + \gamma m \dot{\mathbf{v}}, and γ˙=γ3(β ⁣⋅ ⁣β˙)\dot\gamma = \gamma^3(\boldsymbol\beta\!\cdot\!\dot{\boldsymbol\beta}). For a purely transverse push, β ⁣⋅ ⁣a=0\boldsymbol\beta\!\cdot\!\mathbf{a} = 0, so γ˙=0\dot\gamma = 0, the first term dies, and only γma\gamma m\mathbf{a} survives. For a purely longitudinal push, both terms act, and they sum to γ3ma\gamma^3 m\mathbf{a}. The extra γ2\gamma^2 is precisely the cost of the energy that a speed change requires and a direction change does not. It is not two masses in any deep sense — it is one dp/dt\mathrm{d}\mathbf{p}/\mathrm{d}t, resolved into components.

If the velocity is frozen, how does an accelerator accelerate anything?

It does not accelerate the speed; it accelerates the momentum and the energy, which are unbounded even though the speed is not. An RF cavity adds energy every turn; EE and pp climb without limit; β\beta creeps toward 11 and stops mattering. Everything the machine subsequently does is a response to the momentum: the dipoles must be ramped in step with pp to keep the orbit fixed (equation 3), the quadrupole gradients must be ramped too because k∝1/pk \propto 1/p, and the RF must track a revolution frequency that changes almost not at all. The word "accelerator" is a historical accident. It is a momentum-raiser.

Is the beta function in Hill's equation the same beta as v/c?

No, and this is the single worst notation collision in accelerator physics. The Twiss beta function β(s)\beta(s) is an envelope with units of length — the LHC squeezes it to β∗=0.30\beta^* = 0.30 m at the collision point — and it describes how far the beam swings off-axis. The relativistic beta β=v/c\beta = v/c is dimensionless, close to 11, and effectively constant. They appear in the same equation in (6): εN=βγ ε\varepsilon_N = \beta\gamma\,\varepsilon uses the velocity. Nothing but context distinguishes them, and both are universal. Read carefully, and check the units.

What actually goes wrong at transition, physically?

The bunch stops being able to correct itself. Normally a particle with too much energy is nudged back by the RF, because its momentum error translates into a timing error, and the RF punishes bad timing. At γ=γt\gamma = \gamma_t the translation coefficient η\eta is zero: a particle can have a large momentum error and still arrive exactly on time. The RF is therefore blind to it, the restoring force vanishes (Ωs∝∣η∣→0\Omega_s \propto \sqrt{|\eta|} \to 0), and the bunch — held together by nothing — begins to lengthen and to spread in momentum. Worse, η\eta then reverses sign, so the phase that was stable becomes unstable: the RF must jump from ϕs\phi_s to π−ϕs\pi - \phi_s within milliseconds. The CERN PS has done this on every cycle since 1959, and it still costs beam.

Why does the LHC not have this problem?

Because its designers made γt\gamma_t irrelevant by construction. The transition gamma depends only on the lattice — γt=1/αc\gamma_t = 1/\sqrt{\alpha_c} — and the LHC's momentum compaction αc=3.45×10−4\alpha_c = 3.45\times10^{-4} puts γt=53.8\gamma_t = 53.8. But the LHC does not build its own beam from rest: it inherits one from the SPS at 450450 GeV, where γ=480\gamma = 480, already nine times above transition. The LHC therefore spends its entire operating life above transition and never crosses it. The problem was not solved. It was placed outside the machine's range — which is usually the better engineering.

Does the beam radiate while all this is going on?

Yes, and for electrons it dominates everything. A bent relativistic charge radiates with a power scaling as γ4/ρ2\gamma^4/\rho^2, so a per-turn energy loss of U0[keV]=88.5 E4[GeV]/ρ[m]U_0[\mathrm{keV}] = 88.5\,E^4[\mathrm{GeV}]/\rho[\mathrm{m}] — which at LEP2's 104.5104.5 GeV and ρ=3026\rho = 3026 m came to 3.493.49 GeV per turn, or 3.3%3.3\% of the entire beam energy, thrown away every 8989 microseconds and replaced by a vast RF plant. That E4/m4E^4/m^4 law is what ended LEP and what put protons in its tunnel instead. It is the subject of `aq-sr-synchrotron-radiation`, and it is deliberately not derived here: this article is about how you steer a beam, not about what it emits while you do.

Before this

Where this leads