Rapidity

Board a rocket and boost to half the speed of light. From that rocket, launch a second rocket at half the speed of light relative to the first. Then a third, from the second. Do it ten times, each boost identical, each one a perfectly ordinary half-light-speed shove relative to the ship you just left. Arithmetic says you are now doing five times the speed of light. Arithmetic is wrong. You are doing 0.9999661c0.9999661c — and the exact fraction, if you want it, is 29524/2952529524/29525. Ten enormous boosts have bought you a hundredth of a percent, and no number of further boosts will ever buy you the last sliver. Something has clearly failed to add. But something did add, perfectly and linearly, all ten times: a quantity called the rapidity, φ=artanh⁡β\varphi = \operatorname{artanh}\beta. Each half-light-speed boost adds exactly artanh⁡(0.5)=0.549306\operatorname{artanh}(0.5) = 0.549306 to it, so after ten you have φ=5.493061\varphi = 5.493061, on the nose, with no correction terms and no relativistic funny business. Velocity is what a speedometer reads. Rapidity is what a boost actually does. And the reason relativistic velocities behave so strangely is simply that we have been keeping score with the wrong variable.

Key takeaways
  • Signature (+,−,−,−)(+,-,-,-): Δs2=c2Δt2−Δx2\Delta s^2 = c^2\Delta t^2 - \Delta x^2, timelike >0> 0. The identity γ2−(γβ)2=1\gamma^2 - (\gamma\beta)^2 = 1 is exactly cosh⁡2φ−sinh⁡2φ=1\cosh^2\varphi - \sinh^2\varphi = 1, so a boost is a hyperbolic rotation.
  • Rapidity φ=artanh⁡β=12ln⁡1+β1−β\varphi = \operatorname{artanh}\beta = \tfrac12\ln\frac{1+\beta}{1-\beta}, and then β=tanh⁡φ\beta = \tanh\varphi, γ=cosh⁡φ\gamma = \cosh\varphi, γβ=sinh⁡φ\gamma\beta = \sinh\varphi.
  • Collinear boosts add their rapidities: Λ(φ1)Λ(φ2)=Λ(φ1+φ2)\Lambda(\varphi_1)\Lambda(\varphi_2) = \Lambda(\varphi_1+\varphi_2), exactly. Einstein's velocity-composition law is the tanh⁡\tanh addition formula, in disguise.
  • Owned result: velocity is not a group parameter; rapidity is. The 1+11{+}1D boost group is (R,+)(\mathbb{R},+) — and artanh⁡\operatorname{artanh} is, up to an overall scale, the only function of β\beta that adds. cc is unreachable because artanh⁡\operatorname{artanh} has a pole at β=1\beta = 1: reaching light speed costs infinite rapidity.
  • Exact anchor: β=0.6  ⟺  φ=ln⁡2\beta = 0.6 \iff \varphi = \ln 2, giving γ=5/4\gamma = 5/4 and γβ=3/4\gamma\beta = 3/4.
  • Key idea: boosts compose by addition, always — just not addition of velocity. Change variables from β\beta to φ\varphi and every awkward relativistic formula collapses into ordinary arithmetic.
The core insight

Relativity did not break addition; we were adding the wrong thing. A boost is a rotation through an imaginary angle — a hyperbolic rotation — and rotations compose by adding angles. The hyperbolic angle is the rapidity. Once you parametrise a boost by φ\varphi instead of β\beta, the velocity-addition law, the γ\gamma factor, the impossibility of reaching cc, and the whole apparatus of hyperbolic identities become one fact restated four ways.

Where this lives today

Rapidity is not a pedagogical convenience. It is the working variable of experimental particle physics. At a hadron collider the interesting collisions happen between partons carrying unknown, unequal fractions of the beam momentum, so the centre of momentum of any given event is moving along the beam axis by an unknown amount — and every measured quantity that depends on velocity along that axis is therefore contaminated. Rapidity solves this exactly: because a boost along the beam simply adds a constant to the rapidity of every particle, differences of rapidity are invariant under the unknown longitudinal boost, and so are shapes of rapidity distributions. Detectors are built with their granularity uniform in rapidity for precisely this reason. The numbers are worth seeing. A 6.5 TeV6.5\ \mathrm{TeV} proton in the LHC has γ=E/mc2=6500/0.938272=6927.6\gamma = E/mc^2 = 6500/0.938272 = 6927.6, hence a rapidity of φ=arcosh⁡γ=9.536\varphi = \operatorname{arcosh}\gamma = 9.536. Its β\beta is 0.999999990.99999999 — a number that conveys nothing and cannot usefully be plotted. Its rapidity, 9.549.54, is a perfectly ordinary quantity that sits comfortably on a linear axis alongside the rapidities of the particles it produces. And notice the scaling: doubling the beam energy adds only ln⁡2=0.69\ln 2 = 0.69 to the rapidity. The entire history of accelerator physics, in one number. Rapidity also gives the cleanest statement of what constant acceleration means. A rocket whose crew feels a constant gg underfoot is not one whose velocity increases linearly — that would take it past cc within a year. It is one whose rapidity increases linearly: φ=aτ/c\varphi = a\tau/c, where τ\tau is the proper time aboard and aa the proper acceleration. One Julian year of 1g1g gives φ=(9.80665)(3.15576×107)/(2.99792458×108)=1.0323\varphi = (9.80665)(3.15576\times10^7)/(2.99792458\times10^8) = 1.0323, hence β=tanh⁡(1.0323)=0.7748\beta = \tanh(1.0323) = 0.7748 and γ=1.5818\gamma = 1.5818. The Galilean bookkeeping would have read gτ=3.09×108 m/sg\tau = 3.09\times10^8\ \mathrm{m/s}, which exceeds cc — the failure is not subtle. (The resulting worldline is a hyperbola, and its geometry and its horizon belong to the articles on Relativistic Acceleration and on Rindler Coordinates; they are not treated here.)

The physical picture

Start from the one equation everyone already knows and nobody reads carefully:

γ=11−β2⟹γ2−(γβ)2=1.\gamma = \frac{1}{\sqrt{1-\beta^2}} \qquad\Longrightarrow\qquad \gamma^2 - (\gamma\beta)^2 = 1.

Look at what that says. The pair (γ, γβ)(\gamma,\ \gamma\beta) — which is exactly the pair of numbers appearing in the boost matrix — lies on the unit hyperbola X2−Y2=1X^2 - Y^2 = 1, for every β\beta whatsoever. It cannot lie anywhere else. And a point on the unit hyperbola has a standard parametrisation, one that has been sitting in the textbooks since the 1750s: X=cosh⁡φX = \cosh\varphi, Y=sinh⁡φY = \sinh\varphi. So there is a number φ\varphi such that γ=cosh⁡φ\gamma = \cosh\varphi and γβ=sinh⁡φ\gamma\beta = \sinh\varphi. Divide the second by the first and you get β=tanh⁡φ\beta = \tanh\varphi. Invert, and φ=artanh⁡β\varphi = \operatorname{artanh}\beta. That is the entire definition of rapidity, and it was forced on us: nothing was chosen, we simply named the parameter that the geometry already had. The payoff is immediate and enormous. In Euclidean space, the transformation that mixes xx and yy is a rotation, and two rotations compose by adding angles — nobody would dream of asking how two rotation angles "combine" beyond adding them. In Minkowski space, the transformation that mixes ctct and xx is a boost, and two boosts compose by adding rapidities. Same statement, one sign changed. What made relativistic velocity addition look bizarre is that β=tanh⁡φ\beta = \tanh\varphi is a nonlinear function of the thing that actually adds — like insisting on adding the sines of two angles and then being astonished that sin⁡(θ1+θ2)≠sin⁡θ1+sin⁡θ2\sin(\theta_1 + \theta_2) \ne \sin\theta_1 + \sin\theta_2.

Interactive Simulation
Minkowski Boost ↔ Euclidean Rotation — a boost is the hyperbolic analogue of a rotation; the rapidity φ = atanh β is the hyperbolic angle. Compare the two side by side.
Figure
PUNCHLINE — N successive boosts of 0.5c each, relative to the ship you just left. GALILEO says the speeds add: 5c after ten boosts (the dashed line, straight off the top of the frame). RELATIVITY says β saturates: β_N = (3^N − 1)/(3^N + 1) exactly, reaching 0.9999661c at N = 10 and never reaching 1. RAPIDITY says φ adds, exactly and forever: φ_N = N × 0.549306, a straight line with no ceiling. Galileo was right about the shape of the law — he was wrong about which variable obeys it
speed (units of c) · rapidity φ N boosts 0123456 246810 β = 1 (light speed) GALILEO: 0.5N c → 5c at N = 10 (absurd) 0.5 0.8 RELATIVITY: β saturates 0.9999661 at N = 10 RAPIDITY: φ adds φ = 5.4931 at N = 10 0.5493 EXACTLY, at each N N = 1 · β = 2/4 = 0.5 N = 2 · β = 8/10 = 0.8 N = 3 · β = 26/28 = 0.9286 N = 10 · β = 29524/29525 = 0.9999661 φ = N × 0.549306 always Galileo would say 5c γ = 121.502 at N = 10 βN = tanh(N · artanh ½) = (3N − 1)/(3N + 1) — an EXACT rational, for every N the boosts always add. It is the VELOCITY that refuses to — because β = tanh φ is nonlinear in φ

The core equations

Rapidity is defined by inverting the hyperbolic tangent, and everything else follows.

φ  ≡  artanh⁡β  =  12ln⁡ ⁣1+β1−β,β=tanh⁡φ(1)\varphi \;\equiv\; \operatorname{artanh}\beta \;=\; \tfrac{1}{2}\ln\!\frac{1+\beta}{1-\beta}, \qquad\qquad \beta = \tanh\varphi \qquad (1)

where: β=v/c\beta = v/c. Physical reading: φ\varphi is the hyperbolic angle of the boost. It runs over the whole real line, −∞<φ<∞-\infty < \varphi < \infty, while β\beta is confined to (−1,1)(-1, 1) — and the map between them is a bijection. Every attainable velocity has a rapidity; no rapidity, however large, corresponds to β=1\beta = 1. The two factors that dominate every relativistic formula are now just cosh⁡\cosh and sinh⁡\sinh:

γ=cosh⁡φ,γβ=sinh⁡φ,γ2−(γβ)2=cosh⁡2φ−sinh⁡2φ=1(2)\gamma = \cosh\varphi, \qquad \gamma\beta = \sinh\varphi, \qquad \gamma^2 - (\gamma\beta)^2 = \cosh^2\varphi - \sinh^2\varphi = 1 \qquad (2)

where: γ=(1−β2)−1/2\gamma = (1-\beta^2)^{-1/2}. Physical reading: the famous identity γ2(1−β2)=1\gamma^2(1-\beta^2) = 1 is the hyperbolic Pythagoras. It was never anything else. This is why γ\gamma and γβ\gamma\beta are always the two legs of a hyperbola and never of a circle. The boost matrix, written in φ\varphi, stops looking like a formula and starts looking like a rotation:

Λ(φ)=(cosh⁡φ−sinh⁡φ−sinh⁡φcosh⁡φ)acting on(ctx)(3)\Lambda(\varphi) = \begin{pmatrix}\cosh\varphi & -\sinh\varphi \\ -\sinh\varphi & \cosh\varphi\end{pmatrix} \quad\text{acting on}\quad \begin{pmatrix}ct \\ x\end{pmatrix} \qquad (3)

where: Λ\Lambda is the standard boost along +x+x. Physical reading: compare the Euclidean rotation (cos⁡θ−sin⁡θsin⁡θcos⁡θ)\begin{pmatrix}\cos\theta & -\sin\theta \\ \sin\theta & \cos\theta\end{pmatrix}. The trigonometric functions have become hyperbolic and one sign has moved; that is the whole difference between turning your axes and changing your speed. It also makes det⁡Λ=cosh⁡2−sinh⁡2=1\det\Lambda = \cosh^2 - \sinh^2 = 1 obvious at a glance. And the composition law is now the simplest equation in relativity:

Λ(φ1) Λ(φ2)=Λ(φ1+φ2),tanh⁡(φ1+φ2)=tanh⁡φ1+tanh⁡φ21+tanh⁡φ1tanh⁡φ2(4)\Lambda(\varphi_1)\,\Lambda(\varphi_2) = \Lambda(\varphi_1 + \varphi_2), \qquad\qquad \tanh(\varphi_1+\varphi_2) = \frac{\tanh\varphi_1 + \tanh\varphi_2}{1 + \tanh\varphi_1\tanh\varphi_2} \qquad (4)

where: the second identity is the standard tanh⁡\tanh addition formula. Physical reading: read it right to left. Written in β\beta, it is βtot=(β1+β2)/(1+β1β2)\beta_{\text{tot}} = (\beta_1+\beta_2)/(1+\beta_1\beta_2) — Einstein's velocity-composition law, the one that looks like an ad-hoc relativistic correction. Written in φ\varphi, it is φtot=φ1+φ2\varphi_{\text{tot}} = \varphi_1 + \varphi_2. They are the same equation. The "correction" was an artefact of the coordinate we chose for speed.

Where rapidity comes from

Derivation AThe boost is a rotation through a hyperbolic angle

The claim that a Lorentz boost "is" a hyperbolic rotation is usually presented as an elegant re-writing, an optional flourish for the mathematically inclined. It is not optional. The hyperbolic form is forced by a single algebraic fact about γ\gamma, and once you see the forcing there is no other parametrisation available. What follows is that derivation, done in the (+,−,−,−)(+,-,-,-) signature, with no appeal to any prior knowledge of hyperbolic functions beyond their definitions.

Step 1. Write the boost and extract its two entries. The standard boost along +x+x sends ct↦γ(ct−βx)ct \mapsto \gamma(ct - \beta x) and x↦γ(x−β ct)x \mapsto \gamma(x - \beta\, ct), so its matrix is built entirely out of the two numbers γ\gamma and γβ\gamma\beta. Why: whatever a boost is, it is determined by this pair. Understanding boosts means understanding the pair.

Step 2. Discover the constraint the pair obeys. From γ=(1−β2)−1/2\gamma = (1-\beta^2)^{-1/2} we get γ2(1−β2)=1\gamma^2(1-\beta^2) = 1, that is, γ2−(γβ)2=1\gamma^2 - (\gamma\beta)^2 = 1. So the point (X,Y)=(γ, γβ)(X, Y) = (\gamma,\ \gamma\beta) lies on the unit hyperbola X2−Y2=1X^2 - Y^2 = 1 — for every β\beta, without exception. Why: this is the geometric content of the Lorentz factor, and it is exactly the metric's minus sign showing through. In a Euclidean world the constraint would read X2+Y2=1X^2 + Y^2 = 1 and we would be on a circle.

Step 3. Parametrise the hyperbola. Any point on the branch X≥1X \ge 1 of X2−Y2=1X^2 - Y^2 = 1 can be written X=cosh⁡φX = \cosh\varphi, Y=sinh⁡φY = \sinh\varphi for exactly one real φ\varphi, because cosh⁡2−sinh⁡2=1\cosh^2 - \sinh^2 = 1 identically and sinh⁡\sinh is a bijection from R\mathbb{R} onto R\mathbb{R}. Therefore γ=cosh⁡φ\gamma = \cosh\varphi and γβ=sinh⁡φ\gamma\beta = \sinh\varphi. Why: this is not an ansatz; it is the standard parametrisation of the curve we have just been forced onto, and it is unique.

Step 4. Read off β\beta and invert. Dividing, β=sinh⁡φ/cosh⁡φ=tanh⁡φ\beta = \sinh\varphi/\cosh\varphi = \tanh\varphi, so φ=artanh⁡β=12ln⁡1+β1−β\boxed{\varphi = \operatorname{artanh}\beta = \tfrac12\ln\frac{1+\beta}{1-\beta}}. Substituting into the matrix gives Λ(φ)\Lambda(\varphi) with entries cosh⁡φ\cosh\varphi and −sinh⁡φ-\sinh\varphi. Why: the logarithmic form is worth memorising — it is what makes φ\varphi computable, and it is where the pole at β=1\beta = 1 will come from.

Step 5. Check the anchor. Put β=3/5\beta = 3/5. Then φ=12ln⁡8/52/5=12ln⁡4=ln⁡2\varphi = \tfrac12\ln\frac{8/5}{2/5} = \tfrac12\ln 4 = \ln 2, exactly. And indeed cosh⁡(ln⁡2)=12(2+12)=5/4=γ\cosh(\ln 2) = \tfrac12(2 + \tfrac12) = 5/4 = \gamma and sinh⁡(ln⁡2)=12(2−12)=3/4=γβ\sinh(\ln 2) = \tfrac12(2 - \tfrac12) = 3/4 = \gamma\beta, with (5/4)2−(3/4)2=1(5/4)^2 - (3/4)^2 = 1. Why: β=0.6\beta = 0.6 is the one velocity at which everything in relativity comes out as a clean fraction, and the reason is precisely that its rapidity is ln⁡2\ln 2, so eφ=2e^{\varphi} = 2 and e−φ=12e^{-\varphi} = \tfrac12. *Verified symbolically in verify_aq-sr-rapidity.py: ΛTηΛ=η\Lambda^{\mathsf T}\eta\Lambda = \eta and det⁡Λ=1\det\Lambda = 1 in both the β\beta-form and the φ\varphi-form; the two forms are identical at β=tanh⁡φ\beta = \tanh\varphi; and artanh⁡(3/5)=ln⁡2\operatorname{artanh}(3/5) = \ln 2, cosh⁡(ln⁡2)=5/4\cosh(\ln 2) = 5/4, sinh⁡(ln⁡2)=3/4\sinh(\ln 2) = 3/4 exactly.*

Go deeper · Derivation B · 2 minVelocity is not a group parameter. Rapidity is.
Derivation BVelocity is not a group parameter. Rapidity is.

Here is the article's owned result, and it is a statement about group theory, not about kinematics. The collinear boosts form a group; a one-dimensional continuous group has a parameter in which composition is simply addition; and the question is which function of the velocity that parameter is. The answer is artanh⁡β\operatorname{artanh}\beta — and, crucially, nothing else. Rapidity is not merely a convenient change of variable. It is the unique additive coordinate on the boost group, and velocity is provably not one.

Step 1. Show the rapidities add. Multiply two boosts in the form (3). The entries of the product are cosh⁡φ1cosh⁡φ2+sinh⁡φ1sinh⁡φ2=cosh⁡(φ1+φ2)\cosh\varphi_1\cosh\varphi_2 + \sinh\varphi_1\sinh\varphi_2 = \cosh(\varphi_1+\varphi_2) and −(sinh⁡φ1cosh⁡φ2+cosh⁡φ1sinh⁡φ2)=−sinh⁡(φ1+φ2)-(\sinh\varphi_1\cosh\varphi_2 + \cosh\varphi_1\sinh\varphi_2) = -\sinh(\varphi_1+\varphi_2), by the hyperbolic addition formulas. So Λ(φ1)Λ(φ2)=Λ(φ1+φ2)\Lambda(\varphi_1)\Lambda(\varphi_2) = \Lambda(\varphi_1+\varphi_2), with no approximation and no restriction. Why: this is a one-parameter group, and φ\varphi is its parameter. Composition of boosts is addition of reals — the simplest group there is.

Step 2. Collect the group axioms and name the group. Λ(0)=1\Lambda(0) = \mathbb{1} (identity), Λ(φ)−1=Λ(−φ)\Lambda(\varphi)^{-1} = \Lambda(-\varphi) (inverses), and composition is ++ on R\mathbb{R}, which is associative and commutative. Hence the 1+1D boost group is (R, +)\boxed{\text{the }1{+}1\text{D boost group is }(\mathbb{R},\,+)}. Why: boosts along a line are, as an abstract group, nothing more exotic than the real numbers under addition. All the apparent strangeness lives in the map φ↦β\varphi \mapsto \beta, not in the group.

Step 3. Recover Einstein's law as a corollary. Apply tanh⁡\tanh to Step 1: βtot=tanh⁡(φ1+φ2)=tanh⁡φ1+tanh⁡φ21+tanh⁡φ1tanh⁡φ2=β1+β21+β1β2\beta_{\text{tot}} = \tanh(\varphi_1 + \varphi_2) = \dfrac{\tanh\varphi_1 + \tanh\varphi_2}{1+\tanh\varphi_1\tanh\varphi_2} = \dfrac{\beta_1 + \beta_2}{1 + \beta_1\beta_2}. Why: the velocity-composition law is not an extra postulate, a correction, or a coincidence. It is the tanh⁡\tanh addition formula wearing a physicist's hat, and it has been a theorem of hyperbolic trigonometry since Riccati and Saladini wrote it down in the 1760s — a century and a half before anyone needed it.

Step 4. Show velocity fails. Compose β=0.6\beta = 0.6 with itself: (0.6+0.6)/(1+0.36)=1.2/1.36=15/17=0.882353…(0.6+0.6)/(1+0.36) = 1.2/1.36 = 15/17 = 0.882353\ldots, not 1.21.2. And γ\gamma fails too: γ(15/17)=17/8=2.125\gamma(15/17) = 17/8 = 2.125, whereas γ(0.6)+γ(0.6)=2.5\gamma(0.6) + \gamma(0.6) = 2.5. But the rapidities: ln⁡2+ln⁡2=ln⁡4\ln 2 + \ln 2 = \ln 4, and tanh⁡(ln⁡4)=(4−14)/(4+14)=15/17\tanh(\ln 4) = (4 - \tfrac14)/(4 + \tfrac14) = 15/17 — dead on, with cosh⁡(ln⁡4)=17/8\cosh(\ln 4) = 17/8 and sinh⁡(ln⁡4)=15/8\sinh(\ln 4) = 15/8. Why: neither β\beta nor γ\gamma is additive, and no amount of wishing makes them so. Only φ\varphi is.

Step 5. Prove the uniqueness — the part that makes this a theorem rather than a nice observation. Suppose some function f(β)f(\beta) is additive: f(β1⊕β2)=f(β1)+f(β2)f(\beta_1 \oplus \beta_2) = f(\beta_1) + f(\beta_2), where ⊕\oplus is the composition law. Compose an infinitesimal ϵ\epsilon onto β\beta: expanding (β+ϵ)/(1+βϵ)(\beta+\epsilon)/(1+\beta\epsilon) to first order gives β⊕ϵ=β+ϵ(1−β2)+O(ϵ2)\beta \oplus \epsilon = \beta + \epsilon(1-\beta^2) + O(\epsilon^2). Additivity then requires f(β)+f′(β) ϵ(1−β2)=f(β)+f′(0) ϵf(\beta) + f'(\beta)\,\epsilon(1-\beta^2) = f(\beta) + f'(0)\,\epsilon, so f′(β)=f′(0)/(1−β2)f'(\beta) = f'(0)/(1-\beta^2). Integrating, and using f(0)=0f(0) = 0 (the identity boost has zero parameter), f(β)=f′(0)⋅artanh⁡β\boxed{f(\beta) = f'(0)\cdot\operatorname{artanh}\beta}. Why: this closes the case. Any additive parametrisation of the boost group is rapidity, up to a choice of units. There is no second candidate hiding somewhere — the differential equation has one solution.

Step 6. Read off the impossibility of cc. Since φ=artanh⁡β\varphi = \operatorname{artanh}\beta and artanh⁡\operatorname{artanh} has a pole at β=1\beta = 1, reaching light speed requires infinite rapidity: φ→∞\varphi \to \infty as β→1−\beta \to 1^-. Equivalently, tanh⁡φ<1\tanh\varphi < 1 for every finite φ\varphi, however large. Why: this is the sharpest statement of why cc is a barrier. It is not that a massive body needs infinite energy as a brute fact — it is that in the variable in which boosting is additive, light speed sits at infinity. You cannot get there by adding finitely many finite amounts. *Verified symbolically in verify_aq-sr-rapidity.py: Λ(φ1)Λ(φ2)=Λ(φ1+φ2)\Lambda(\varphi_1)\Lambda(\varphi_2) = \Lambda(\varphi_1+\varphi_2) exactly; Λ(0)=1\Lambda(0) = \mathbb{1}; Λ(φ)−1=Λ(−φ)\Lambda(\varphi)^{-1} = \Lambda(-\varphi); tanh⁡(φ1+φ2)\tanh(\varphi_1+\varphi_2) equals the Einstein composition; 0.6⊕0.6=15/170.6 \oplus 0.6 = 15/17 with γ=17/8\gamma = 17/8 (and γ\gamma is confirmed non-additive); β⊕ϵ=β+ϵ(1−β2)+O(ϵ2)\beta \oplus \epsilon = \beta + \epsilon(1-\beta^2) + O(\epsilon^2); the ODE f′=k/(1−β2)f' = k/(1-\beta^2) with f(0)=0f(0) = 0 has the unique solution kartanh⁡βk\operatorname{artanh}\beta; and artanh⁡β→∞\operatorname{artanh}\beta \to \infty as β→1−\beta \to 1^-.*

Go deeper · Derivation C · 2 minThe ladder of ten boosts, in closed form
Derivation CThe ladder of ten boosts, in closed form

The rocket in the overview is not a thought experiment that needs numerical integration. Because rapidities add, NN identical boosts is just NN times one rapidity, and the whole ladder collapses to a single exact formula that you can evaluate by hand for any NN. This is what a good change of variables buys you: a problem that is intractable in β\beta (compose the Einstein law with itself ten times, and try to keep the fractions straight) becomes a one-line multiplication in φ\varphi.

Step 1. Compute the rung. One boost of β=1/2\beta = 1/2 has rapidity φ1=artanh⁡12=12ln⁡1+1/21−1/2=12ln⁡3=0.549306…\varphi_1 = \operatorname{artanh}\tfrac12 = \tfrac12\ln\frac{1 + 1/2}{1 - 1/2} = \tfrac12\ln 3 = 0.549306\ldots Why: this is the only number in the whole problem. Everything else is arithmetic on it.

Step 2. Climb. After NN boosts, φN=Nφ1=N2ln⁡3\varphi_N = N\varphi_1 = \tfrac{N}{2}\ln 3, so eφN=3N/2e^{\varphi_N} = 3^{N/2}. Then

βN=tanh⁡φN=eφN−e−φNeφN+e−φN=3N/2−3−N/23N/2+3−N/2=3N−13N+1\beta_N = \tanh\varphi_N = \frac{e^{\varphi_N} - e^{-\varphi_N}}{e^{\varphi_N} + e^{-\varphi_N}} = \frac{3^{N/2} - 3^{-N/2}}{3^{N/2} + 3^{-N/2}} = \boxed{\frac{3^N - 1}{3^N + 1}}

after multiplying top and bottom by 3N/23^{N/2}. Why: an exact rational, for every NN, with no square roots and no floating point anywhere. The awkwardness of the Einstein law has completely evaporated.

Step 3. Check the low rungs against the composition law directly. N=1N=1: (3−1)/(3+1)=1/2(3-1)/(3+1) = 1/2. ✓ N=2N=2: (9−1)/(9+1)=4/5=0.8(9-1)/(9+1) = 4/5 = 0.8 — and Einstein's law gives (12+12)/(1+14)=1/1.25=0.8(\tfrac12+\tfrac12)/(1+\tfrac14) = 1/1.25 = 0.8. ✓ N=3N=3: (27−1)/(27+1)=13/14=0.928571…(27-1)/(27+1) = 13/14 = 0.928571\ldots Why: the closed form is not a trick; it agrees with composing the velocities one at a time, as it must.

Step 4. Reach the top. N=10N = 10: φ=5ln⁡3=ln⁡243=5.493061…\varphi = 5\ln 3 = \ln 243 = 5.493061\ldots, and

β10=310−1310+1=5904859050=2952429525=0.9999661304…\beta_{10} = \frac{3^{10}-1}{3^{10}+1} = \frac{59048}{59050} = \frac{29524}{29525} = 0.9999661304\ldots

with γ10=cosh⁡(5ln⁡3)=12 ⁣(243+1243)=59050486=121.502058…\gamma_{10} = \cosh(5\ln 3) = \tfrac12\!\left(243 + \tfrac{1}{243}\right) = \frac{59050}{486} = 121.502058\ldots Why: ten half-light-speed boosts, and you are still 0.0000339c0.0000339c short of light. The Galilean bookkeeping would have said 5c5c — off not by a correction but by a category error.

Step 5. Contemplate the last sliver. To halve the gap 1−β1 - \beta you must add roughly 12ln⁡2=0.347\tfrac12\ln 2 = 0.347 to φ\varphi, since 1−tanh⁡φ≈2e−2φ1 - \tanh\varphi \approx 2e^{-2\varphi} for large φ\varphi. Every further boost of 0.5c0.5c adds 0.5490.549 and so shrinks the gap by a factor of 33 — forever, without ever closing it. Why: this is the geometry of the pole. Light speed is not a wall you can reach and be stopped at; it is a point at infinity in the only variable that composes. *Verified symbolically in verify_aq-sr-rapidity.py: artanh⁡(1/2)=12ln⁡3\operatorname{artanh}(1/2) = \tfrac12\ln 3; the closed form βN=(3N−1)/(3N+1)\beta_N = (3^N-1)/(3^N+1) is confirmed against tanh⁡(Nartanh⁡12)\tanh(N\operatorname{artanh}\tfrac12) for N=1…12N = 1\ldots12; β10=29524/29525\beta_{10} = 29524/29525 and γ10=59050/486\gamma_{10} = 59050/486 exactly, with γ102−(γβ)102=1\gamma_{10}^2 - (\gamma\beta)_{10}^2 = 1; and the companion β=0.6\beta = 0.6 ladder βN=(4N−1)/(4N+1)\beta_N = (4^N-1)/(4^N+1) is confirmed for N=1…8N = 1\ldots8.*

Figure
The rapidity ruler. Equal boosts are equal steps in φ — the top scale is uniform, forever. The SAME boosts, read in β on the bottom scale, pile up against the wall at β = 1 and never touch it. Marked: the exact anchor β = 0.6 at φ = ln 2 = 0.693, the two-boost point 15/17 at φ = ln 4, the three-boost point 63/65 at φ = ln 8 — and, far off to the right, a 6.5 TeV LHC proton at φ = 9.54, which in β is 0.99999999
RAPIDITY φ — uniform, unbounded 012345678910 → ∞ VELOCITY β = tanh φ — crowded, capped at 1 β = 1 · UNREACHABLE 00.76160.96400.995 φ = 1 → β = 0.762 φ = 2 → β = 0.964 φ = 3 → β = 0.995 φ = 4 → β = 0.9993 ln 2 β = 0.6 ln 4 · β = 15/17 ln 8 · β = 63/65 φ = 9.54 · LHC proton Three identical 0.6c boosts: φ = ln 2 + ln 2 + ln 2 = ln 8, so β = 63/65 = 0.9692, γ = 65/16 = 4.0625 Galileo would have said 1.8c. The rapidities added; the velocities never do.

Limiting cases

CaseResultMeaning
φ→0\varphi \to 0β≈φ\beta \approx \varphi, γ≈1+φ2/2\gamma \approx 1 + \varphi^2/2rapidity is v/cv/c at low speed — which is why Galilean velocity addition works so well
β\beta smallφ=β+β3/3+⋯\varphi = \beta + \beta^3/3 + \cdotsthe first correction is third order: at β=0.1\beta = 0.1, φ\varphi and β\beta differ by 0.03%0.03\%
φ\varphi largeγ≈eφ/2\gamma \approx e^{\varphi}/2γ\gamma grows exponentially in the additive variable — one more boost multiplies it
φ\varphi large1−β≈2e−2φ1 - \beta \approx 2e^{-2\varphi}the gap to light speed shrinks geometrically and never closes
β→1−\beta \to 1^-φ→+∞\varphi \to +\inftythe pole. Light speed costs infinite rapidity
Photonφ\varphi undefined (infinite)light has no rest frame, and correspondingly no finite rapidity
φ<0\varphi < 0boost in −x-xrapidity is signed, and Λ(−φ)=Λ(φ)−1\Lambda(-\varphi) = \Lambda(\varphi)^{-1}
Non-collinear boostsrapidities do not simply addadditivity is a collinear statement; see Hyperbolic Geometry of Boosts

Worked examples

1 · Ten boosts of half light speed

A rocket boosts to 0.5c0.5c. From it, a second rocket boosts to 0.5c0.5c relative to the first. Ten stages. What is the final speed relative to the original frame?

Step 1 — identify. Each stage is a collinear boost, so the rapidities add and nothing else does. The velocity-composition law is the awkward route; the rapidity route is one multiplication. Compute the rapidity of one stage: φ1=artanh⁡(0.5)=12ln⁡3=0.5493061\varphi_1 = \operatorname{artanh}(0.5) = \tfrac12\ln 3 = 0.5493061.

Step 2 — apply. Ten identical stages: φ10=10×0.5493061=5.4930614\varphi_{10} = 10 \times 0.5493061 = 5.4930614. Exactly, with no correction terms. (In closed form, φ10=5ln⁡3=ln⁡243\varphi_{10} = 5\ln 3 = \ln 243.)

Step 3 — convert back. β10=tanh⁡(5ln⁡3)\beta_{10} = \tanh(5\ln 3). Since eφ10=243e^{\varphi_{10}} = 243, \beta_{10} = \frac{243 - 1/243}{243 + 1/243} = \frac{59048}{59050} = \frac{29524}{29525} = 0.9999661304.

Step 4 — cross-check on the low rungs. Two stages by Einstein's law: (0.5+0.5)/(1+0.25)=0.8(0.5+0.5)/(1+0.25) = 0.8; by the closed form, (32−1)/(32+1)=8/10=0.8(3^2-1)/(3^2+1) = 8/10 = 0.8. ✓ Three stages: 13/14=0.92857113/14 = 0.928571. ✓ The two routes agree, as they must — they are the same law.

Step 5 — the boundary number. The Galilean answer is 10×0.5c=5c10 \times 0.5c = 5c. The true answer is 0.99997c0.99997c. The discrepancy is not a small correction to be added on; it is the difference between adding the right variable and the wrong one. Meanwhile γ10=cosh⁡(5ln⁡3)=12(243+1243)=121.502\gamma_{10} = \cosh(5\ln 3) = \tfrac12(243 + \tfrac{1}{243}) = 121.502 — so although the speed has barely moved since stage five, the ship's energy and momentum have gone up by a factor of over a hundred. That is the physical content of the saturation: the boosts are still doing real work; they are simply no longer buying speed.

Step 6 — read the structure. Each further boost multiplies eφe^{\varphi} by 3\sqrt{3}, hence multiplies γ\gamma by (asymptotically) 3\sqrt{3} and divides the gap 1−β1 - \beta by 33. You can do this forever. You will never arrive.

2 · The exact anchor, and three boosts of 0.6c

Show that β=0.6\beta = 0.6 has rapidity exactly ln⁡2\ln 2, and use it to compose three such boosts without ever touching the velocity-addition law.

Step 1 — identify. Use the logarithmic form, φ=12ln⁡1+β1−β\varphi = \tfrac12\ln\frac{1+\beta}{1-\beta}, which turns a rapidity into a logarithm of a rational number.

Step 2 — the anchor. φ=12ln⁡1+3/51−3/5=12ln⁡8/52/5=12ln⁡4=ln⁡2=0.693147\varphi = \tfrac12\ln\dfrac{1 + 3/5}{1 - 3/5} = \tfrac12\ln\dfrac{8/5}{2/5} = \tfrac12\ln 4 = \ln 2 = 0.693147. So eφ=2e^{\varphi} = 2 and e−φ=12e^{-\varphi} = \tfrac12, exactly, and hence γ=cosh⁡φ=12(2+12)=54\gamma = \cosh\varphi = \tfrac12(2 + \tfrac12) = \tfrac54 and γβ=sinh⁡φ=12(2−12)=34\gamma\beta = \sinh\varphi = \tfrac12(2 - \tfrac12) = \tfrac34. Check: (5/4)2−(3/4)2=25/16−9/16=1(5/4)^2 - (3/4)^2 = 25/16 - 9/16 = 1. ✓

Step 3 — apply, three times. Three collinear boosts of 0.6c0.6c: φtot=3ln⁡2=ln⁡8\varphi_{\text{tot}} = 3\ln 2 = \ln 8. So eφtot=8e^{\varphi_{\text{tot}}} = 8, and immediately \beta_{\text{tot}} = \frac{8 - 1/8}{8 + 1/8} = \frac{63}{65} = 0.969231, \qquad \gamma_{\text{tot}} = \frac{8 + 1/8}{2} = \frac{65}{16} = 4.0625, \qquad \gamma\beta = \frac{63}{16}. Check: (65/16)2−(63/16)2=(652−632)/256=256/256=1(65/16)^2 - (63/16)^2 = (65^2 - 63^2)/256 = 256/256 = 1. ✓

Step 4 — cross-check the hard way. Compose two: (0.6+0.6)/(1+0.36)=15/17(0.6+0.6)/(1+0.36) = 15/17. Compose the third: (15/17+3/5)/(1+9/17)=63/65(15/17 + 3/5)/(1 + 9/17) = 63/65. ✓ Identical — but note how much more work it was, and that the intermediate 15/1715/17 had to be carried exactly to avoid rounding.

Step 5 — the boundary number. Galileo would report 1.8c1.8c. Relativity reports 0.969c0.969c, and it reports it as the exact rational 63/6563/65, with γ=65/16\gamma = 65/16 exactly. The general rule for this ladder is βN=(4N−1)/(4N+1)\beta_N = (4^N - 1)/(4^N + 1), since eφN=2Ne^{\varphi_N} = 2^N — which is why β=0.6\beta = 0.6 is the velocity that makes every relativity problem come out in clean fractions. Its rapidity is a logarithm of 22.

Physics of the portrait

Rapidity has no discoverer to put in a portrait, because the mathematics was finished long before the physics needed it. The portrait therefore belongs to the functions themselves.

Origins · 1 minHyperbolic functions and how this was discovered
HF
Hyperbolic functions
1757–1911 · Bologna · Mulhouse · Zagreb · Cambridge

The functions were introduced by Vincenzo Riccati (1707–1775) in Bologna, who defined sinh⁡\sinh and cosh⁡\cosh in 1757 in the first volume of his Opusculorum ad res physicas et mathematicas pertinentium by way of the geometry of the unit hyperbola X2−Y2=1X^2 - Y^2 = 1 — deliberately in parallel with the circular functions on X2+Y2=1X^2 + Y^2 = 1. In the Institutiones Analyticae (1765–1767), written with his student Girolamo Saladini, he set down the addition formulas, including the one that would become Einstein's velocity-composition law. Johann Heinrich Lambert (1728–1777) systematised them in the following decade, credited Riccati for the names, and fixed the abbreviations still in use. For a century and a half they were a curiosity of analysis. Then Minkowski's 1908 geometry made cosh⁡\cosh and sinh⁡\sinh the natural coordinates of a boost; Vladimir Varičak (1865–1942), a Croatian mathematician who had spent 1903–1908 working on Bolyai–Lobachevsky geometry, published in 1910 the reading of special relativity as hyperbolic geometry, in which velocities compose by the hyperbolic triangle rule; and in 1911 Alfred Arthur Robb (1873–1936), in Optical Geometry of Motion, gave the hyperbolic angle its physical name: rapidity. A formula written for the pleasure of it in 1765 turned out, in 1911, to be a law of nature.

History
1757
— Vincenzo Riccati defines the hyperbolic sine and cosine from the geometry of the unit hyperbola.
1765
— Riccati and Girolamo Saladini publish the hyperbolic addition formulas in the Institutiones Analyticae.
1770
— Johann Heinrich Lambert systematises the functions and fixes the modern notation.
1908
— Minkowski's Cologne lecture makes spacetime a geometry with an indefinite metric, in which boosts are hyperbolic rotations.
1910
— Vladimir Varičak reads special relativity as Bolyai-Lobachevsky geometry: velocities compose by the hyperbolic triangle rule.
1911
— Alfred A. Robb, in Optical Geometry of Motion, names the hyperbolic angle of a boost the rapidity.
1912
— Varičak's non-Euclidean interpretation of relativity is published in full.

Going further

Texts. Rindler, Introduction to Special Relativity, 2nd ed., §§2, 9 · Taylor & Wheeler, Spacetime Physics, 2nd ed., ch. 8 ("velocity parameter") · Landau & Lifshitz, The Classical Theory of Fields, §4 · Robb, Optical Geometry of Motion (Heffer, Cambridge, 1911) · Varičak, "Die Relativtheorie und die Lobatschefskijsche Geometrie", Physikalische Zeitschrift 11, 93 (1910).

Provenance. Signature (+,−,−,−)(+,-,-,-). All numbers recomputed in verify_aq-sr-rapidity.py (92 checks): the anchor artanh⁡(3/5)=ln⁡2\operatorname{artanh}(3/5) = \ln 2 with γ=5/4\gamma = 5/4, γβ=3/4\gamma\beta = 3/4; the group law Λ(φ1)Λ(φ2)=Λ(φ1+φ2)\Lambda(\varphi_1)\Lambda(\varphi_2) = \Lambda(\varphi_1+\varphi_2); the uniqueness of artanh⁡\operatorname{artanh} as the additive parametrisation; the ladders βN=(3N−1)/(3N+1)\beta_N = (3^N-1)/(3^N+1) and (4N−1)/(4N+1)(4^N-1)/(4^N+1); β10=29524/29525=0.9999661304\beta_{10} = 29524/29525 = 0.9999661304 with γ10=59050/486\gamma_{10} = 59050/486; the LHC proton at φ=9.536\varphi = 9.536; and 1g1g for one Julian year giving φ=1.0323\varphi = 1.0323, β=0.7748\beta = 0.7748.

Scope. Non-collinear composition, the rotation it induces, and the Rindler treatment of constant proper acceleration are not derived here — see Hyperbolic Geometry of Boosts, Thomas Precession, and Relativistic Acceleration. See also: aq-sr-rapidity-hyperbolic-boosts.

Builds on. Hyperbolic Geometry of Boosts · The Spacetime Interval · Four-Vectors. Leads to. Lorentz Invariance · Causal Structure of Spacetime · Spacetime Diagrams.

Is rapidity a real physical quantity, or just a change of variable?

It is a change of variable, and that is exactly what makes it physical. Nothing you can measure is "the rapidity" in the way a speedometer measures speed — but the group of boosts is a physical structure, and rapidity is its canonical coordinate, the one in which composition is addition. Compare angle: no instrument reads "angle" directly either, yet nobody doubts that rotating by 30°30° and then 40°40° is rotating by 70°70°. Rapidity is the angle of a boost, and it earns its reality the same way.

If rapidity is so much better, why is every textbook formula written in terms of γ and β?

Historical accident, mostly, plus the fact that β\beta is what you measure. But the moment a problem involves composing boosts — chained frames, successive accelerations, collider kinematics, the relativistic rocket — the β\beta formulation becomes a swamp and the φ\varphi formulation becomes trivial. A useful working rule: if you are transforming between two frames, use γ\gamma and β\beta; if you are composing transformations, switch to φ\varphi first and switch back at the end.

Does rapidity add for boosts in different directions?

No, and this is the one trap in the subject. Rapidity is additive along a single axis; that is a statement about a one-parameter subgroup. Compose a boost along xx with a boost along yy and the result is not a pure boost at all — it is a boost together with a rotation, whose angle is not zero and whose existence is the source of Thomas precession. Treating rapidity as a vector that adds componentwise is wrong, and it is wrong at second order, which is exactly where the interesting physics is. See Hyperbolic Geometry of Boosts.

What is the rapidity of a photon?

Infinite, which is another way of saying that it does not have one. φ=artanh⁡(1)\varphi = \operatorname{artanh}(1) diverges. This is the same fact as "light has no rest frame": there is no finite boost that brings a photon to rest, because there is no finite rapidity you can subtract from infinity. Physically, the light cone is invariant — a boost slides null directions along themselves, rescaling them by e±φe^{\pm\varphi} but never tilting them off the cone — so no amount of rapidity gets you there.

Why do particle physicists use pseudorapidity instead?

Because true rapidity y=12ln⁡E+pzcE−pzcy = \tfrac12\ln\frac{E + p_z c}{E - p_z c} requires knowing the particle's energy, which requires knowing its mass, which a tracking detector often does not. Pseudorapidity η=−ln⁡tan⁡(θ/2)\eta = -\ln\tan(\theta/2), where θ\theta is the polar angle from the beam, depends only on the direction of the track — and in the ultra-relativistic limit E≫mc2E \gg mc^2 it converges to the true rapidity. It is a purely geometric stand-in, adopted because it is measurable, and it inherits rapidity's crucial property that longitudinal boosts shift it by a constant.

The formula φ = aτ/c says a 1g rocket gains rapidity forever. Is there no limit at all?

None on the rapidity, and that is the point. After one year of proper time φ=1.03\varphi = 1.03; after ten, φ=10.3\varphi = 10.3 and γ=15 000\gamma = 15\,000; after fifty, φ=51.6\varphi = 51.6 and γ\gamma is around 102210^{22}. The velocity saturates immediately and thereafter changes by amounts too small to write down, but the rapidity — and with it the energy, the momentum, and the crew's ability to cross the galaxy within a lifetime — grows without bound. Nothing forbids arbitrarily large rapidity. What is forbidden is arbitrarily large velocity, and those two statements only look like the same statement if you are using the wrong variable.

Before this

Where this leads