Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Energy and Momentum from Symmetry: a Unified Derivation

This note derives the kinetic energy and momentum of a free particle in both Galilean and special relativity from a single common axiom set. The two theories are two specializations of one functional equation; the only difference is a small “package” — the boost (velocity vs rapidity) together with one companion axiom (mass conserved vs total energy conserved). Everything else is identical.

The detailed case-by-case derivations live in Kinetic Energy is Quadratic: a Galilean-Invariance Derivation; this note’s role is to exhibit their common root.

The common axioms

Both theories share, verbatim:

Axiom (A4) is the heart: it says the collision’s two bookkeeping quantities — the mass MM and the dissipation QQ — are each frame-invariant. What distinguishes the two theories is which of them is independently fixed.

The package: the only difference

Each theory selects a boost and fixes one of the two frame-invariant scalars in (A4):

Galilean packageRelativistic package
boostvelocities add: ξ=v\xi=v, uubu\mapsto u-brapidities add: ξ=ϕ=artanh(v/c)\xi=\phi=\operatorname{artanh}(v/c), ϕϕβ\phi\mapsto\phi-\beta
fixed scalarmass conserved: M=imiM=\sum_i m_itotal energy conserved: Q=0Q=0
free scalarQQ (the heat)MM (the blob’s rest mass — it grows)

The two rows are not independent choices: they come locked as a pair (see Why the packages are locked). Picking the Galilean boost forces mass conservation for a sensible theory; picking the Lorentz boost forces total-energy conservation.

The single functional equation

Derive (1). Take the symmetric collision: two masses mm at ±ξ\pm\xi, blob at rest (A3) of mass MM. The lab-frame balance (A4) is

2mE(ξ)=ME(0)+Q.2m\,\mathcal E(\xi)=M\,\mathcal E(0)+Q.

View it from a frame boosted by β\beta (A2). The particles move to ξβ\xi-\beta and ξβ-\xi-\beta; the blob, still of mass MM (frame-invariant, A4), moves to β-\beta. The balance there is

mE(ξβ)+mE(ξβ)=ME(β)+Q.m\,\mathcal E(\xi-\beta)+m\,\mathcal E(-\xi-\beta)=M\,\mathcal E(\beta)+Q.

(The blob’s energy is ME(β)M\mathcal E(\beta) by evenness; the dissipation QQ is the same number, frame-invariant by A4.) Subtract (3) from (4) and use evenness:

E(ξ+β)+E(ξβ)=2E(ξ)+Mm[E(β)E(0)].\boxed{\,\mathcal E(\xi+\beta)+\mathcal E(\xi-\beta)=2\,\mathcal E(\xi)+\frac{M}{m}\bigl[\mathcal E(\beta)-\mathcal E(0)\bigr]\,.}

This is the entire theory before the package is chosen. Now specialize.

Galilean specialization

M=2mM=2m (mass conserved), ξ=v\xi=v. Then M/m=2M/m=2 and (5) becomes

E(v+b)+E(vb)=2E(v)+2[E(b)E(0)].\mathcal E(v+b)+\mathcal E(v-b)=2\,\mathcal E(v)+2\bigl[\mathcal E(b)-\mathcal E(0)\bigr].

Writing the kinetic energy T(v):=E(v)E(0)T(v):=\mathcal E(v)-\mathcal E(0) (so T(0)=0T(0)=0),

T(v+b)+T(vb)=2T(v)+2T(b),T(v+b)+T(v-b)=2\,T(v)+2\,T(b),

the parallelogram law. With regularity (A5) its solutions are T(v)=cv2T(v)=c\,v^2 (the Hamel-basis pathologies excluded by measurability). Hence

Ekin=mT(v)=cmv2.E_{\text{kin}}=m\,T(v)=c\,m\,v^2.

The constant cc is fixed by energy units (Newton’s convention c=12c=\tfrac12, giving 12mv2\tfrac12mv^2); the rest energy E(0)\mathcal E(0) is a free additive constant — in Galilean mechanics energy is shiftable, so we conventionally set E(0)=0\mathcal E(0)=0.

Relativistic specialization

Q=0Q=0 (total energy conserved), so (3) fixes M=2mE(ξ)/E(0)M=2m\,\mathcal E(\xi)/\mathcal E(0), i.e. M/m=2E(ξ)/E(0)M/m=2\,\mathcal E(\xi)/\mathcal E(0). With ξ=ϕ\xi=\phi, (5) becomes

E(ϕ+β)+E(ϕβ)=2E(ϕ)+2E(ϕ)E(0)[E(β)E(0)]=2E(0)E(ϕ)E(β),\mathcal E(\phi+\beta)+\mathcal E(\phi-\beta)=2\,\mathcal E(\phi)+\frac{2\,\mathcal E(\phi)}{\mathcal E(0)}\bigl[\mathcal E(\beta)-\mathcal E(0)\bigr]=\frac{2}{\mathcal E(0)}\,\mathcal E(\phi)\,\mathcal E(\beta),

d’Alembert’s equation. With regularity and the physical selection (energy \ge rest, minimized at ϕ=0\phi=0), its solution is

E(ϕ)=E(0)coshϕ.\mathcal E(\phi)=\mathcal E(0)\cosh\phi.

The non-relativistic limit fixes the constant: Ekin=m[E(ϕ)E(0)]=mE(0)(coshϕ1)mE(0)v2/(2c2)E_{\text{kin}}=m[\mathcal E(\phi)-\mathcal E(0)]=m\,\mathcal E(0)(\cosh\phi-1)\approx m\,\mathcal E(0)\,v^2/(2c^2), and matching 12mv2\tfrac12mv^2 forces

E(0)=c2,henceE=γmc2,Ekin=(γ1)mc2.\mathcal E(0)=c^2, \qquad\text{hence}\qquad \mathcal E=\gamma\,mc^2,\quad E_{\text{kin}}=(\gamma-1)mc^2.

Note the contrast with the Galilean case: here E(0)=c2\mathcal E(0)=c^2 is forced, not a free constant — Lorentz covariance forbids shifting the rest energy. And the blob’s mass, M=2mcoshϕ=2γmM=2m\cosh\phi=2\gamma m, grows: the dissipated kinetic energy became rest mass, Qrel=ΔMc2Q_{\text{rel}}=\Delta M\,c^2 — itself a derived theorem, not an assumption.

Momentum, unified

Momentum is derived from the same collision by the same move (write momentum balance in the boosted frame), giving a second functional equation that is the exact companion of (5).

Assume momentum is extensive in mass (a body of mass mm at ξ\xi carries momentum mπ(ξ)m\,\pi(\xi), with π\pi odd) and conserved in every frame. The boosted-frame balance for the symmetric collision reads

mπ(ξβ)+mπ(ξβ)=Mπ(β).m\,\pi(\xi-\beta)+m\,\pi(-\xi-\beta)=M\,\pi(-\beta).

Oddness on both sides gives the unified momentum equation

π(ξ+β)π(ξβ)=Mmπ(β).\boxed{\,\pi(\xi+\beta)-\pi(\xi-\beta)=\frac{M}{m}\,\pi(\beta)\,.}

Galilean (M/m=2M/m=2, ξ=v\xi=v): π(v+b)π(vb)=2π(b)\pi(v+b)-\pi(v-b)=2\pi(b). Swapping vbv\leftrightarrow b and using oddness yields π(v+b)+π(vb)=2π(v)\pi(v+b)+\pi(v-b)=2\pi(v); adding gives π(v+b)=π(v)+π(b)\pi(v+b)=\pi(v)+\pi(b), Cauchy’s equation, whose regular solution is π(v)=v\pi(v)=v (the slope is a free unit convention, set to 1 so that p=mvp=mv). Thus

p=mπ(v)=mv.p=m\pi(v)=mv.

Relativistic (M/m=2coshϕM/m=2\cosh\phi, ξ=ϕ\xi=\phi): π(ϕ+β)π(ϕβ)=2coshϕπ(β)\pi(\phi+\beta)-\pi(\phi-\beta)=2\cosh\phi\,\pi(\beta). Swapping and adding gives π(ϕ+β)=coshϕπ(β)+coshβπ(ϕ)\pi(\phi+\beta)=\cosh\phi\,\pi(\beta)+\cosh\beta\,\pi(\phi), the sinh\sinh addition formula, whose regular solution is π(ϕ)=csinhϕ\pi(\phi)=c\sinh\phi — here the prefactor is forced to the invariant speed cc by the non-relativistic limit pmvp\to mv. Thus

p=mcsinhϕ=γmv.p=mc\sinh\phi=\gamma mv.

So the same two equations, (5) and (13), produce both {Ekin=12mv2, p=mv}\{E_{\text{kin}}=\tfrac12mv^2,\ p=mv\} and {Ekin=(γ1)mc2, p=γmv}\{E_{\text{kin}}=(\gamma-1)mc^2,\ p=\gamma mv\}, distinguished only by the package. The energy equation fixes the “even” function (cosh\cosh or quadratic) and the momentum equation its “odd” companion (sinh\sinh or linear); together they assemble the Galilean pair (p,Ekin)=(mv,12mv2)(p,E_{\text{kin}})=(mv,\tfrac12mv^2) or the relativistic four-momentum Pμ=(E/c,p)=(γmc,γmv)P^\mu=(E/c,p)=(\gamma mc,\gamma mv).

Four-momentum and the conservation law

The momentum derived above is not an isolated three-vector. In both spacetimes it is the spatial part of a four-momentum Pμ=muμP^\mu=m\,u^\mu (mass times four-velocity). The relativistic time slot is E/cE/c; the Galilean time slot is mm. Thus the Galilean object is more precisely a mass-momentum four-vector, not an energy-momentum four-vector.

This section builds that object in full for both theories and proves a key structural fact: once the appropriate components are conserved, they assemble into a single frame-independent conservation law; conversely, frame-independence couples the components, so momentum conservation cannot hold without the time component (mass, or energy) being conserved too.

What “four-vector” means here

There are two levels of covariance to keep distinct.

First, in flat spacetime with inertial Cartesian/adapted frames, the admissible inertial frame changes are the ten continuous symmetries: rotations, translations, boosts, and time translations (plus optional discrete reflections). These are the Galilei transformations in Galilean spacetime and the Poincare transformations in Minkowski spacetime. Checking how PμP^\mu transforms under these transformations proves that its components form the correct inertial-frame four-vector.

Second, in the differential-geometric sense, a tensor is a coordinate-independent object whose components transform tensorially under arbitrary smooth coordinate changes. That stronger statement is proved not by checking every coordinate system, but by giving a coordinate-free definition: uμu^\mu is the tangent to the worldline, and Pμ=muμP^\mu=m u^\mu with mm a scalar. Then the tensor transformation law follows by the chain rule. The Galilei and Poincare formulas below are the inertial-coordinate special cases of this general tensorial statement.

Also, tensorial means covariant, not invariant. Components may change from one frame to another; what matters is that they change by the correct transformation law.

Definition versus conservation

Two claims must be kept separate.

First, Pμ=muμP^\mu=m u^\mu is the definition of four-momentum. Since uμu^\mu is the worldline tangent and mm is a scalar, this also proves that PμP^\mu is a vector. This is a kinematic/geometric statement: it determines how PμP^\mu transforms between frames, but it does not by itself say that PμP^\mu is conserved.

Second, conservation is an additional dynamical statement:

inPμ=outPμ.\sum_{\rm in}P^\mu=\sum_{\rm out}P^\mu.

The earlier collision setup supplies or assumes the component conservation laws: in Galilean mechanics, mass conservation plus momentum conservation; in relativity, energy conservation plus momentum conservation. Once these components are conserved in one inertial frame, the vector character of PμP^\mu packages them into one law and makes that law frame-independent. Thus Pμ=muμP^\mu=m u^\mu proves tensoriality, while the collision/symmetry argument supplies conservation.

The four-velocity

Parametrize a particle’s worldline by the time natural to each spacetime.

Relativistic — Minkowski spacetime. The invariant parameter is the proper time τ\tau, related to coordinate time by dτ=dt/γd\tau=dt/\gamma. The four-velocity is

uμ:=dxμdτ=γdxμdt=(γc, γv),γ=11v2/c2.u^\mu:=\frac{dx^\mu}{d\tau}=\gamma\frac{dx^\mu}{dt}=(\gamma c,\ \gamma\mathbf v), \qquad \gamma=\frac{1}{\sqrt{1-v^2/c^2}}.

It is a four-vector; under a Lorentz transformation uμ=Λμνuνu'^\mu=\Lambda^\mu{}_\nu u^\nu, with constant Minkowski norm ημνuμuν=γ2(c2v2)=c2\eta_{\mu\nu}u^\mu u^\nu=\gamma^2(c^2-v^2)=c^2.

Galilean — absolute time. Time is absolute, so in adapted coordinates the natural parameter is coordinate time tt. Equivalently, the Newton-Cartan clock one-form τμ\tau_\mu normalizes the worldline tangent by τμuμ=1\tau_\mu u^\mu=1. In such adapted coordinates the four-velocity is

uμ:=dxμdt=(1, v).u^\mu:=\frac{dx^\mu}{dt}=(1,\ \mathbf v).

Under a Galilean boost, including a possible rotation (x=Rxbt+a\mathbf x'=R\mathbf x-\mathbf b\,t+\mathbf a, t=t+a0t'=t+a^0),

u0=u0,u=Rubu0,u'^0=u^0,\qquad \mathbf u'=R\mathbf u-\mathbf b\,u^0,

so for u0=1u^0=1 one has v=Rvb\mathbf v'=R\mathbf v-\mathbf b — precisely the transformation law of a Galilean four-vector.

The four-momentum Pμ=muμP^\mu=m\,u^\mu

Multiply by the rest mass mm (a frame-invariant scalar in both theories):

four-velocity uμu^\mufour-momentum Pμ=muμP^\mu=m\,u^\mu
relativistic(γc, γv)(\gamma c,\ \gamma\mathbf v)(γmc, γmv)=(E/c, p)(\gamma mc,\ \gamma m\mathbf v)=(E/c,\ \mathbf p)
Galilean(1, v)(1,\ \mathbf v)(m, mv)=(m, p)(m,\ m\mathbf v)=(m,\ \mathbf p)

The same formula Pμ=muμP^\mu=m\,u^\mu holds in both. Since mm is a scalar and uμu^\mu is the worldline tangent, PμP^\mu is a genuine vector in the generally covariant formulation. The only difference is what occupies the time slot in adapted inertial coordinates — E/cE/c (relativistic) or mm (Galilean).

For a Galilean boost this gives

P0=P0=m,P=RPbP0=Rpmb.P'^0=P^0=m,\qquad \mathbf P'=R\mathbf P-\mathbf b\,P^0 =R\mathbf p-m\mathbf b.

The term mb-m\mathbf b has the elementary meaning m(vb)=mvmbm(\mathbf v-\mathbf b)=m\mathbf v-m\mathbf b: changing to a frame moving with velocity b\mathbf b subtracts the frame’s velocity from every particle velocity, and therefore subtracts mbm\mathbf b from its momentum. This is why p\mathbf p alone is not a Galilean spacetime vector. It is a three-vector under rotations, but under boosts it is not closed by itself; it needs the mass component P0=mP^0=m.

The metric then reads off the “mass-shell”:

Conservation, and a one-line proof of frame-independence

For a collision with incoming momenta PiμP_i^\mu and outgoing PjμP_j^\mu, conservation of four-momentum is the single statement

inPμ=outPμ.\sum_{\rm in}P^\mu=\sum_{\rm out}P^\mu.

Because PμP^\mu is a four-vector it transforms linearly, Pμ=ΛμνPνP'^\mu=\Lambda^\mu{}_\nu P^\nu, so (23) holds in every inertial frame:

inPμ=ΛμνinPν=ΛμνoutPν=outPμ.\sum_{\rm in}P'^\mu=\Lambda^\mu{}_\nu\sum_{\rm in}P^\nu =\Lambda^\mu{}_\nu\sum_{\rm out}P^\nu=\sum_{\rm out}P'^\mu.

Linearity of Λ\Lambda is the whole proof. Reading off components gives the conservation laws:

The components are coupled (the heart of the matter)

The frame-independence proof conceals the key point: a boost mixes the time and space components, so conserving one component alone does not survive a change of frame. The components are locked together. In relativity the mixing is two-way: E/cE/c and p\mathbf p are not separate Lorentz tensors, only components of the four-vector PμP^\mu. In Galilean spacetime the mixing is triangular: the mass component is invariant, but the spatial momentum still receives a contribution from it. Thus p\mathbf p is a tensor only for the reduced rotation subgroup, not for the full Galilei group with boosts.

Galilean — momentum conservation forces mass conservation. Assume only that momentum is conserved in the lab, inp=outp\sum_{\rm in}\mathbf p=\sum_{\rm out}\mathbf p, and ask whether it holds in a frame boosted by b\mathbf b. Since p=pmb\mathbf p'=\mathbf p-m\mathbf b,

inp=inpb ⁣inm,outp=outpb ⁣outm.\sum_{\rm in}\mathbf p'=\sum_{\rm in}\mathbf p-\mathbf b\!\sum_{\rm in}m,\qquad \sum_{\rm out}\mathbf p'=\sum_{\rm out}\mathbf p-\mathbf b\!\sum_{\rm out}m.

Subtracting,

inpoutp=(inpoutp)=0 (assumption)    b(inmoutm).\sum_{\rm in}\mathbf p'-\sum_{\rm out}\mathbf p' =\underbrace{\bigl(\sum_{\rm in}\mathbf p-\sum_{\rm out}\mathbf p\bigr)}_{=\,0\text{ (assumption)}} \;-\;\mathbf b\bigl(\sum_{\rm in}m-\sum_{\rm out}m\bigr).

For this to vanish for every boost b\mathbf b,

inm=outm.\sum_{\rm in}m=\sum_{\rm out}m.

So frame-independent momentum conservation entails mass conservation. The two are inseparable: they are the time and space slots of one four-vector, and the boost couples them.

Relativistic — momentum conservation forces energy conservation. Take a Lorentz boost of rapidity β\beta along xx,

E=γβ(Eβcpx),px=γβ ⁣(pxβcE),E'=\gamma_\beta(E-\beta c\,p_x),\qquad p'_x=\gamma_\beta\!\left(p_x-\frac{\beta}{c}E\right),

and assume only inpx=outpx\sum_{\rm in}p_x=\sum_{\rm out}p_x. Then

inpxoutpx=γβ ⁣[(inpxoutpx)=0    βc(inEoutE)],\sum_{\rm in}p'_x-\sum_{\rm out}p'_x =\gamma_\beta\!\left[\underbrace{\bigl(\sum_{\rm in}p_x-\sum_{\rm out}p_x\bigr)}_{=\,0} \;-\;\frac{\beta}{c}\bigl(\sum_{\rm in}E-\sum_{\rm out}E\bigr)\right],

which vanishes for every β\beta only if inE=outE\sum_{\rm in}E=\sum_{\rm out}E. So in relativity frame-independent momentum conservation entails energy conservation, and conversely. The invariant E2p2c2=m2c4E^2-p^2c^2=m^2c^4 is the same coupling in static form.

What this explains, and the Galilean asymmetry

The coupling is exactly what the package demanded and what our derivations produced, now seen geometrically. In the Galilean package the companion axiom was mass conservation; the four-vector shows it is forced on us by frame-independent momentum conservation — the two come as one. In the relativistic package the companion axiom was total-energy conservation; the four-vector shows it is forced by frame-independent momentum conservation just the same.

The asymmetry, now sharp: the Galilean four-momentum (m,p)(m,\mathbf p) contains mass and momentum, not energy. Energy conservation in Galilean mechanics is therefore a separate law — outside the four-vector — and (as derived in Kinetic Energy is Quadratic: a Galilean-Invariance Derivation) it follows from the frame-invariance of heat, not from four-vector structure. Relativity needs no such extra law: its four-momentum already carries the energy in the time slot, so the single statement (23) covers everything. One conservation law (relativistic) versus one plus an extra (Galilean) — that is the precise cost of the degenerate metric.

Adding gravity: from inertial forces to curvature

The derivation so far was the flat-spacetime theory. If we write the same theory in arbitrary coordinates, including accelerating or rotating coordinates, the equations already contain connection terms. Those terms are fictitious forces: they come from a non-inertial coordinate choice, and the corresponding curvature is still zero. True gravity begins when the connection has nonzero curvature, so that no coordinate choice can remove the effect everywhere.

The bottom-up extension to gravity follows the same pattern in the Galilean and relativistic theories. The common steps are:

  1. Universal free fall. All test bodies fall the same way, independently of their mass and composition.

  2. Geometrize the universal force. Since the acceleration is universal, it should not be a force attached separately to each body. It is encoded in a common spacetime connection \nabla, and free particles follow geodesics.

  3. Only tidal effects are invariant. A uniform gravitational field can be removed locally by a freely falling frame. Thus the connection coefficients themselves are not the invariant gravitational field; curvature is.

  4. Choose the conserved source. Galilean gravity is sourced by mass density. Relativistic gravity is sourced by the full stress-energy tensor.

  5. Locality and second-order simplicity. The field equation should be local, covariant, and second order in the gravitational variables.

  6. Fix the normalization. The remaining constant is fixed by the observed Newtonian weak/static limit; symmetry fixes the form, not the numerical value of GG.

Universal free fall gives the geodesic equation

The equivalence principle is the shared starting point:

Galilean / Newton-CartanRelativistic / GR
clock/metric structureclock one-form τμ\tau_\mu and spatial metric hμνh^{\mu\nu}Lorentzian metric gμνg_{\mu\nu}
normalizationτμuμ=1\tau_\mu u^\mu=1gμνuμuν=c2g_{\mu\nu}u^\mu u^\nu=c^2
free-fall equationuννuμ=0u^\nu\nabla_\nu u^\mu=0uννuμ=0u^\nu\nabla_\nu u^\mu=0

So the equation of motion is literally the same tensor equation in both theories:

uννuμ=0.\boxed{\,u^\nu\nabla_\nu u^\mu=0\,}.

The difference lies in the background geometric structure: Newton-Cartan has absolute time and a degenerate spatial metric; GR has a non-degenerate Lorentzian metric.

In adapted Newton-Cartan coordinates, τ=dt\tau=dt and uμ=(1,x˙)u^\mu=(1,\dot{\mathbf x}). For an ordinary Newtonian potential Φ\Phi, the gravitational connection component is

Γ00i=iΦ,\Gamma^i_{00}=\partial^i\Phi,

so (30) becomes

d2xidt2+Γ00i=0d2xidt2=iΦ.\frac{d^2x^i}{dt^2}+\Gamma^i_{00}=0 \qquad\Longleftrightarrow\qquad \frac{d^2x^i}{dt^2}=-\partial^i\Phi.

In GR, the same geodesic equation in the weak, slowly moving, static limit with

g00=1+2Φc2g_{00}=1+\frac{2\Phi}{c^2}

also reduces to

d2xidt2=iΦ.\frac{d^2x^i}{dt^2}=-\partial^i\Phi.

Thus both theories agree on the motion law in the Newtonian limit.

Curvature is the gravitational field

A freely falling coordinate system can set the connection coefficients to zero at one event. Therefore Γνρμ\Gamma^\mu_{\nu\rho} itself is not the invariant gravitational field. The invariant obstruction to removing gravity throughout a region is curvature.

In Newton-Cartan theory, with Γ00i=iΦ\Gamma^i_{00}=\partial^i\Phi, the tidal tensor is the second derivative of the potential:

Ri0j0=jiΦ,R00=iiΦ=2Φ.R^i{}_{0j0}=\partial_j\partial^i\Phi, \qquad R_{00}=\partial_i\partial^i\Phi=\nabla^2\Phi.

A uniform field has iΦ=constant\partial_i\Phi=\text{constant} and hence zero curvature; it is locally removable by an accelerating frame. Tidal gravity is what remains.

In GR, the same statement is expressed by the Riemann tensor RμνρσR^\mu{}_{\nu\rho\sigma}. In a local inertial frame the connection can be made to vanish at a point, but the curvature tensor generally cannot.

The field equation: curvature sourced by the conserved quantity

The final step is to connect curvature to matter.

Galilean / Newton-Cartan. The source selected by the flat theory is mass. Let ρ\rho be the mass density. Locality, spatial rotational invariance, and second-order simplicity select the scalar second-derivative equation

2Φ=Cρ.\nabla^2\Phi=C\rho.

Writing C=4πGC=4\pi G gives Poisson’s equation,

2Φ=4πGρ.\boxed{\,\nabla^2\Phi=4\pi G\rho\,}.

In Newton-Cartan tensor form this is

Rμν=4πGρτμτν.\boxed{\,R_{\mu\nu}=4\pi G\,\rho\,\tau_\mu\tau_\nu\,}.

In adapted coordinates, the only nontrivial component is R00=4πGρR_{00}=4\pi G\rho, exactly (37).

For a point mass, ρ=Mδ(3)(x)\rho=M\delta^{(3)}(\mathbf x), the rotationally symmetric solution is

Φ(r)=GMr,a=Φ=GMr2r^.\Phi(r)=-\frac{GM}{r}, \qquad \mathbf a=-\nabla\Phi=-\frac{GM}{r^2}\,\hat{\mathbf r}.

Thus the inverse-square law follows from locality, three-dimensional rotational symmetry, and the second-order Poisson equation; GG is the empirical coupling constant.

Relativistic / GR. The source selected by the relativistic flat theory is the full stress-energy tensor TμνT_{\mu\nu}, not just mass density. The left-hand side must be a symmetric, generally covariant, divergence-free tensor built locally from the metric and at most two derivatives. In four dimensions, the Lovelock uniqueness theorem gives

Gμν+ΛgμνG_{\mu\nu}+\Lambda g_{\mu\nu}

as the unique possibility of this type, where Gμν:=Rμν12RgμνG_{\mu\nu}:=R_{\mu\nu}-\tfrac12 Rg_{\mu\nu} and μGμν=0\nabla^\mu G_{\mu\nu}=0 by the Bianchi identity. Therefore the field equation must be

Gμν+Λgμν=8πGc4Tμν.\boxed{\,G_{\mu\nu}+\Lambda g_{\mu\nu} =\frac{8\pi G}{c^4}\,T_{\mu\nu}\,}.

The coefficient is fixed by the weak/static limit: with T00ρc2T_{00}\approx\rho c^2 and g00=1+2Φ/c2g_{00}=1+2\Phi/c^2, the 00 component of (41) reduces to

2Φ=4πGρ.\nabla^2\Phi=4\pi G\rho.

The cosmological constant Λ\Lambda is allowed by the same covariance and second-order axioms; it is an additional empirical constant, set to zero if one demands exactly flat empty spacetime.

The parallel structure is therefore:

universal free fallconnection/geodesicscurvature as gravitycurvature sourced by mass or stress-energy.\text{universal free fall} \Rightarrow \text{connection/geodesics} \Rightarrow \text{curvature as gravity} \Rightarrow \text{curvature sourced by mass or stress-energy}.

With gravity present, the flat-spacetime global conservation law is also replaced by a local covariant statement. In GR this is μTμν=0\nabla_\mu T^{\mu\nu}=0 for matter (with the gravitational field included geometrically rather than as a separate local stress-energy tensor). In Newton-Cartan theory, mass conservation remains the local continuity equation for the mass current, while momentum balance is expressed using the Newton-Cartan connection. Global conserved energies or momenta require extra spacetime symmetries, just as in GR.

The Galilean theory is the cc\to\infty limit

The two packages are not parallel alternatives; one is the singular limit of the other. In the relativistic solution let cc\to\infty at fixed vv. Then rapidity ϕ=artanh(v/c)v/c0\phi=\operatorname{artanh}(v/c)\to v/c\to0, coshϕ1+v2/(2c2)\cosh\phi\to 1+v^2/(2c^2), and

E=c2coshϕc2+12v2,M=2mcoshϕ2m,p=mcsinhϕmv.\mathcal E=c^2\cosh\phi\to c^2+\tfrac12 v^2,\qquad M=2m\cosh\phi\to 2m,\qquad p=mc\sinh\phi\to mv.

The rest energy c2c^2 diverges and drops out of the kinetic energy Ekin=Ec212v2E_{\text{kin}}=\mathcal E-c^2\to\tfrac12 v^2; the blob’s mass M2mM\to2m (mass conservation re-emerges); the mass excess ΔM=2m(coshϕ1)mv2/c20\Delta M=2m(\cosh\phi-1)\to mv^2/c^2\to0 vanishes, while the heat Q=ΔMc2mv2Q=\Delta M\,c^2\to mv^2 stays finite. In words: as cc\to\infty the relativistic bookkeeping “QQ becomes rest mass” degenerates into the Galilean bookkeeping “mass is conserved, QQ is heat” — the heat is the cc\to\infty shadow of the relativistic rest-mass gain. This is why the two packages are locked: the Galilean one is the singular limit of the relativistic one, not an independent option.

Why the packages are locked

One might hope to mix and match — e.g. Lorentz boosts with mass conservation. Equation (5) shows this fails: with M=2mM=2m (mass conserved) but ξ=ϕ\xi=\phi (rapidity), it becomes the parallelogram law in the rapidity variable, T(ϕ+β)+T(ϕβ)=2T(ϕ)+2T(β)T(\phi+\beta)+T(\phi-\beta)=2T(\phi)+2T(\beta), whose solution T(ϕ)ϕ2T(\phi)\propto\phi^2 grows only quadratically with rapidity and does not match the observed γmc2\gamma mc^2. Conversely, Galilean boosts with Q=0Q=0 (total energy conserved) give d’Alembert’s equation in the velocity variable, T(v)cosh(v)T(v)\propto\cosh(v), exponential in speed. Both cross-combinations are mathematically consistent but describe worlds that are not ours.

So the boost and the bookkeeping are not independent dials: a physically correct theory requires them matched. Choosing the boost effectively chooses its companion axiom.

Summary

GalileanRelativistic
boostvvbv\mapsto v-bϕϕβ\phi\mapsto\phi-\beta
companion axiommass conserved (M=2mM=2m)total energy conserved (Q=0Q=0)
energy equationparallelogram lawd’Alembert’s equation
kinetic energy12mv2\tfrac12 mv^2(γ1)mc2(\gamma-1)mc^2
momentum equationCauchy’s equationsinh\sinh addition formula
momentummvmvγmv\gamma mv
four-momentum(m,p)(m,\mathbf p), mass-momentum(E/c,p)(E/c,\mathbf p), energy-momentum
energy in four-vector?no; $E=\mathbf p
free fall with gravityuννuμ=0u^\nu\nabla_\nu u^\mu=0, τμuμ=1\tau_\mu u^\mu=1uννuμ=0u^\nu\nabla_\nu u^\mu=0, gμνuμuν=c2g_{\mu\nu}u^\mu u^\nu=c^2
gravitational field equationRμν=4πGρτμτνR_{\mu\nu}=4\pi G\rho\,\tau_\mu\tau_\nuGμν+Λgμν=8πGTμν/c4G_{\mu\nu}+\Lambda g_{\mu\nu}=8\pi G T_{\mu\nu}/c^4
rest energy E(0)\mathcal E(0)0 (free, shiftable)c2c^2 (forced)
mass in collisionconserved (M=2mM=2m)grows (M=2γmM=2\gamma m)

The common axioms (A1)–(A5) — energy extensive and isotropic, relativity, parity, the frame-invariant collision balance, regularity — produce one energy equation (5) and one momentum equation (13). The sole difference is the package: {velocity boost, mass conserved} versus {rapidity boost, energy conserved}. From this single fork, both 12mv2, mv\tfrac12mv^2,\ mv and (γ1)mc2, γmv(\gamma-1)mc^2,\ \gamma mv unfold.