physics2u
PhysicsMathsChemistryBiology
Tier
⌕ Search ⌘K
Derivation

Newtonian Limit Fixes the Coupling Constant

Statement

In the weak-field (\(g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu}\), \(|h_{\mu\nu}|\ll 1\)), slow-motion (\(v\ll c\)), and static (\(\partial_0 h_{\mu\nu}\approx 0\)) limit, the Einstein field equations \(G_{\mu\nu}=\kappa\,T_{\mu\nu}\) reduce to the Newtonian Poisson equation \(\nabla^2\Phi=4\pi G\rho\) provided the coupling is fixed to \(\kappa=\dfrac{8\pi G}{c^{4}}\), with the identification \(g_{00}=-\left(1+\dfrac{2\Phi}{c^{2}}\right)\).

Why it matters

General relativity is postulated, not deduced: the Einstein–Hilbert action fixes the tensor structure of the field equations but leaves an overall multiplicative constant \(\kappa\) undetermined. Only correspondence with an already-tested theory can pin down that number. The Newtonian limit is that correspondence: it forces the abstract geometry to reproduce two centuries of celestial mechanics and, in doing so, welds \(G\) and \(c\) into the single relativistic coupling \(8\pi G/c^{4}\).

The smallness of \(8\pi G/c^{4}\approx 2\times10^{-43}\ \text{s}^2\,\text{kg}^{-1}\,\text{m}^{-1}\) is why spacetime is so stiff: enormous stress-energy produces minute curvature. Every quantitative test of GR — perihelion precession, light bending, the Hulse–Taylor binary, LIGO waveforms — inherits its absolute scale from this one matching calculation.

Assumptions
Weak field: the metric is a small perturbation of flat spacetime, \(g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu}\) with \(|h_{\mu\nu}|\ll 1\), so we may linearise and raise/lower indices with \(\eta\). Drop it and products \(h\partial h\) survive, the equations stay nonlinear, and no closed Poisson form emerges.
Slow motion: source and test particles move with \(v\ll c\), so \(u^{i}\ll u^{0}\approx c\) and \(T_{00}=\rho c^{2}\) dominates every other component. Drop it and momentum flux and pressure terms contribute, changing the effective source.
Static / quasi-static fields: time derivatives are negligible, \(\partial_0 h_{\mu\nu}\approx 0\), so \(\Box\to\nabla^2\). Drop it and radiative terms \(-\frac{1}{c^2}\partial_t^2 h\) reinstate the wave operator, giving linearised gravitational waves rather than a static potential.
Pressureless dust source: matter is modelled as \(T_{\mu\nu}=\rho\,u_\mu u_\nu\) with negligible internal pressure \(p\ll\rho c^2\). Drop it and the Newtonian source becomes \(\rho+3p/c^2\) (the Tolman "active gravitational mass"), correcting Poisson at relativistic pressures.
Harmonic-free identification is gauge-consistent: the physical identification of \(h_{00}\) with \(\Phi\) presumes a coordinate (gauge) choice in which the linearised Ricci tensor collapses to \(-\tfrac12\nabla^2 h_{00}\). A different gauge redistributes \(h_{\mu\nu}\); the observable \(\nabla^2\Phi\) is gauge-invariant but intermediate lines are not.
Derivation
1
\[ \frac{d^2 x^\mu}{d\tau^2}+\Gamma^{\mu}{}_{\alpha\beta}\frac{dx^\alpha}{d\tau}\frac{dx^\beta}{d\tau}=0 \]
Prior result: the geodesic equation from extremal proper time, the equation of motion for a free test particle. A
2
\[ \frac{d^2 x^{i}}{d\tau^2}+\Gamma^{i}{}_{00}\left(\frac{dx^0}{d\tau}\right)^2\approx 0 \]
Slow motion: \(dx^i/d\tau\ll dx^0/d\tau\), so only the \(\alpha=\beta=0\) term in the double sum survives. A
3
\[ \Gamma^{i}{}_{00}=\tfrac12 g^{i\lambda}\big(2\,\partial_0 g_{\lambda 0}-\partial_\lambda g_{00}\big)\;\xrightarrow{\ \text{static}\ }\;-\tfrac12\,\eta^{ij}\partial_j h_{00}=-\tfrac12\,\partial_i h_{00} \]
Christoffel definition; static kills \(\partial_0\) terms; weak field replaces \(g^{i\lambda}\to\eta^{ij}\) and \(g_{00}\to\eta_{00}+h_{00}\) with \(\partial_j\eta=0\). B
4
\[ \frac{dx^0}{d\tau}=c\frac{dt}{d\tau}\approx c,\qquad d\tau\approx dt \]
Slow motion and weak field give \(d\tau/dt\to 1\) to leading order; higher corrections are \(O(v^2/c^2, h)\). A
5
\[ \frac{d^2 x^{i}}{dt^2}=-\Gamma^{i}{}_{00}\,c^2=\tfrac12 c^2\,\partial_i h_{00} \]
Substitute steps 3–4 into step 2 and use \(d^2x^i/d\tau^2\approx d^2x^i/dt^2\). A
6
\[ \frac{d^2 x^{i}}{dt^2}=-\partial_i\Phi \quad\Longrightarrow\quad \tfrac12 c^2\,\partial_i h_{00}=-\partial_i\Phi \quad\Longrightarrow\quad h_{00}=-\frac{2\Phi}{c^2} \]
Match to Newton's second law with a potential; equality of gradients for all \(i\) fixes \(h_{00}\) up to a constant absorbed into the boundary condition \(\Phi\to 0\). B
7
\[ g_{00}=-1+h_{00}=-\left(1+\frac{2\Phi}{c^2}\right) \]
Assemble the \(00\) metric component; this is the geometric encoding of the Newtonian potential. A
8
\[ G_{\mu\nu}=\kappa\,T_{\mu\nu}\quad\Longleftrightarrow\quad R_{\mu\nu}=\kappa\left(T_{\mu\nu}-\tfrac12 g_{\mu\nu}T\right) \]
Prior result (Einstein field equations from the Einstein–Hilbert action); trace-reverse using \(G_{\mu\nu}=R_{\mu\nu}-\tfrac12 g_{\mu\nu}R\) and \(R=-\kappa T\). B
9
\[ T_{00}=\rho c^2,\qquad T=g^{\mu\nu}T_{\mu\nu}\approx \eta^{00}T_{00}=-\rho c^2 \]
Dust source \(T_{\mu\nu}=\rho u_\mu u_\nu\) with \(u^0\approx c\), \(u^i\approx0\), \(u_0\approx -c\); trace taken with \(\eta^{00}=-1\). C
10
\[ T_{00}-\tfrac12 g_{00}T\approx \rho c^2-\tfrac12(-1)(-\rho c^2)=\tfrac12\rho c^2 \]
Evaluate the trace-reversed source for the \(00\) component, using \(g_{00}\approx\eta_{00}=-1\) to leading order. B
11
\[ R_{00}\;\xrightarrow{\ \text{static, linear}\ }\;\partial_i\Gamma^{i}{}_{00}=\partial_i\!\left(-\tfrac12\partial_i h_{00}\right)=-\tfrac12\nabla^2 h_{00} \]
Linearised Ricci tensor \(R_{00}=\partial_\lambda\Gamma^\lambda{}_{00}-\partial_0\Gamma^\lambda{}_{\lambda 0}\); the second term and all \(\Gamma\Gamma\) products vanish in the static, first-order limit. C
12
\[ R_{00}=-\tfrac12\nabla^2\!\left(-\frac{2\Phi}{c^2}\right)=\frac{1}{c^2}\nabla^2\Phi \]
Substitute the geodesic identification \(h_{00}=-2\Phi/c^2\) from step 6. A
13
\[ \frac{1}{c^2}\nabla^2\Phi=\kappa\cdot\tfrac12\rho c^2 \quad\Longrightarrow\quad \nabla^2\Phi=\tfrac12\kappa c^4\,\rho \]
Set the geometric \(R_{00}\) (step 12) equal to the source \(R_{00}=\kappa(T_{00}-\tfrac12 g_{00}T)\) (steps 8, 10). B
14
\[ \tfrac12\kappa c^4=4\pi G \quad\Longrightarrow\quad \boxed{\;\kappa=\frac{8\pi G}{c^4}\;} \]
Demand agreement with the empirically established Poisson equation \(\nabla^2\Phi=4\pi G\rho\); the coupling is now fixed. A
Result
\[ G_{\mu\nu}=\frac{8\pi G}{c^4}\,T_{\mu\nu},\qquad g_{00}=-\left(1+\frac{2\Phi}{c^2}\right),\qquad \nabla^2\Phi=4\pi G\rho \]

Reading. The undetermined constant in Einstein's equations is not free: forcing gravity to reproduce Newton in the weak, slow, static corner of parameter space fixes it to \(8\pi G/c^4\). Geometrically, the Newtonian potential lives in the time–time part of the metric, and its Laplacian is (up to \(1/c^2\)) the \(00\) Ricci curvature; the \(8\pi\) is the fingerprint of the trace-reversal, while the \(1/c^4\) sets the extreme rigidity of spacetime.

Units check. \([\kappa]=\dfrac{[G]}{[c^4]}=\dfrac{\text{m}^3\,\text{kg}^{-1}\,\text{s}^{-2}}{\text{m}^4\,\text{s}^{-4}}=\text{m}^{-1}\,\text{kg}^{-1}\,\text{s}^{2}\). With \([T_{00}]=\text{J}\,\text{m}^{-3}=\text{kg}\,\text{m}^{-1}\,\text{s}^{-2}\), the product \([\kappa T_{00}]=\text{m}^{-2}=[G_{\mu\nu}]\), the dimension of curvature. Consistent.

Limiting cases
  • \(c\to\infty\): \(\kappa\to 0\) and \(h_{00}=-2\Phi/c^2\to 0\) — spacetime becomes infinitely stiff and flat; only the Poisson content survives, recovering pure Newtonian gravity.
  • \(\Phi\to 0\) (far field / no source): \(\nabla^2\Phi=0\), the Laplace equation; the metric relaxes to \(\eta_{\mu\nu}\) and geodesics are straight lines (special relativity).
  • \(p\to\rho c^2\) (relativistic pressure): the neglected pressure terms enter and the source becomes \(\rho+3p/c^2\); Poisson is superseded by the full Tolman–Oppenheimer–Volkoff structure.
  • Vacuum with \(\Lambda\): restoring a cosmological constant sends \(\nabla^2\Phi=4\pi G\rho-\Lambda c^2\), a constant background curvature added to the Newtonian source.
Breaks when
  • Strong field (\(|h_{\mu\nu}|\sim 1\), e.g. near a black-hole horizon where \(2GM/rc^2\to 1\)): linearisation fails, the neglected \(h\partial h\) terms are comparable to the kept ones, and no single scalar \(\Phi\) captures the geometry.
  • Fast motion / radiation (\(v\sim c\) or \(\partial_t h\ne 0\)): dropping \(u^i\) and \(\partial_0 h\) is illegitimate; the wave operator \(\Box\) reappears and one obtains gravitational waves, not a static potential.
  • Non-negligible pressure or anisotropic stress (relativistic gases, radiation-dominated fluids, neutron-star cores): \(T_{00}\) no longer dominates, the trace changes, and the effective Newtonian source acquires \(3p/c^2\) corrections.
Failure modes
  • Skipping the trace reversal: equating \(R_{00}=\kappa T_{00}\) directly instead of \(R_{00}=\kappa(T_{00}-\tfrac12 g_{00}T)\) yields \(\kappa=4\pi G/c^4\), wrong by a factor of two.
  • Sign of \(h_{00}\): forgetting that \(\eta_{00}=-1\) in the \((-{+}{+}{+})\) signature flips the sign, giving \(g_{00}=-(1-2\Phi/c^2)\) and a repulsive Newtonian force.
  • Using mass density \(\rho\) where energy density \(\rho c^2\) is required: setting \(T_{00}=\rho\) drops a factor \(c^2\) and mis-scales \(\kappa\).
  • Keeping spatial 4-velocity components: retaining \(u^i\) in \(T_{\mu\nu}\) contaminates the source with momentum flux that has no Newtonian counterpart.
  • Confusing \(\Box\) with \(\nabla^2\): not invoking the static assumption leaves \(-\tfrac1{c^2}\partial_t^2\Phi\) in the equation, so a "static" derivation silently smuggles in radiative terms.
Discussion

The calculation is a two-front matching. The geodesic equation supplies the kinematic half: it tells us how a test particle responds to geometry and thereby forces \(h_{00}=-2\Phi/c^2\). The field equation supplies the dynamic half: it tells us how matter sources geometry and thereby relates \(\nabla^2 h_{00}\) to \(\rho\). Only when both halves are demanded to agree with Newton — one for the force law, one for the source law — is the single number \(\kappa\) over-determined and pinned. It is a genuine consistency miracle that the same \(\kappa\) works for both.

The factor \(8\pi\) is often mistaken for a fundamental constant of nature. It is not: it is a bookkeeping consequence of the trace-reversal and of Newton's own \(4\pi\) (which itself comes from the surface area of a sphere in Gauss's law). Had we defined \(\Phi\) with a different normalisation, or written Gauss's law in Heaviside–Lorentz form, the numerical prefactor would shift. What is physical is the combination that appears in any observable, e.g. the Schwarzschild radius \(r_s=2GM/c^2\), which the same limit reproduces from \(g_{00}=-(1+2\Phi/c^2)=-(1-2GM/rc^2)\).

At the level of rigour, the identification \(R_{00}=-\tfrac12\nabla^2 h_{00}\) is gauge-dependent: the linearised Ricci tensor contains extra terms \(\tfrac12(\partial_i\partial_0 h_{0}{}^{i}+\dots)\) that vanish only under the static assumption and an implicit harmonic-type gauge. The gauge-invariant content is that the trace of the tidal tensor \(\partial_i\partial_j\Phi\) equals \(4\pi G\rho\); individual metric components are coordinate artefacts. This is why the "potential" \(\Phi\) is best regarded as a shorthand for a curvature invariant, not a field with independent ontological status — a point that becomes essential when one passes to the parametrised post-Newtonian (PPN) formalism, where multiple potentials with independent coefficients are needed at the next order.

Common misconceptions. (i) That GR "derives" Newtonian gravity — it does not; it recovers it, and only after \(\kappa\) is fed in from Newton by hand. (ii) That the weak-field metric requires only \(g_{00}\); in fact bending of light needs \(g_{ij}=(1-2\Phi/c^2)\delta_{ij}\) too, and using \(g_{00}\) alone underpredicts light deflection by exactly a factor of two. (iii) That \(\Phi<0\) implies \(g_{00}>-1\) means "less time" — the sign chase must be done in a fixed signature convention.

Worked examples

Example 1 — Metric perturbation at the Sun's surface.

1
\[ h_{00}=-\frac{2\Phi}{c^2},\qquad \Phi(R)=-\frac{GM_\odot}{R_\odot}\;\Longrightarrow\; h_{00}=\frac{2GM_\odot}{R_\odot c^2} \]
Newtonian surface potential of a spherical mass; substitute into the identification from step 6. A
2
\[ h_{00}=\frac{2(6.674\times10^{-11})(1.989\times10^{30})}{(6.96\times10^{8})(2.998\times10^{8})^2} \]
Insert \(G,\ M_\odot,\ R_\odot,\ c\) in SI. A
3
\[ h_{00}=\frac{2(6.674\times10^{-11})(1.989\times10^{30})}{(6.96\times10^{8})(8.988\times10^{16})}=\frac{2.655\times10^{20}}{6.256\times10^{25}}\approx 4.24\times10^{-6} \]
Arithmetic; units of \(GM/(Rc^2)\) are dimensionless as required. A
\[ h_{00}\approx 4.2\times10^{-6} \]

Reading. The Sun's gravity distorts the metric at only the part-per-million level, comfortably validating \(|h_{\mu\nu}|\ll1\) and hence the entire weak-field limit for the solar system.

Units check. \([GM/Rc^2]=\dfrac{(\text{m}^3\text{kg}^{-1}\text{s}^{-2})(\text{kg})}{(\text{m})(\text{m}^2\text{s}^{-2})}=1\), dimensionless. Consistent.

Example 2 — Tidal-tensor trace of a uniform star via the fixed coupling.

1
\[ \nabla^2\Phi=\tfrac12\kappa c^4\rho=\frac{8\pi G}{c^4}\cdot\tfrac12 c^4\rho=4\pi G\rho \]
Use the fixed coupling to write the interior source in terms of density alone. B
2
\[ \rho=\bar\rho_\odot=\frac{M_\odot}{\tfrac43\pi R_\odot^3}=\frac{1.989\times10^{30}}{\tfrac43\pi(6.96\times10^{8})^3}\approx 1.41\times10^{3}\ \text{kg}\,\text{m}^{-3} \]
Mean solar density from mass and volume of a sphere. A
3
\[ \nabla^2\Phi=4\pi(6.674\times10^{-11})(1.41\times10^{3})\approx 1.18\times10^{-6}\ \text{s}^{-2} \]
Insert numbers; \(\nabla^2\Phi=\mathrm{tr}(\partial_i\partial_j\Phi)\) is the trace of the tidal tensor. B
\[ \nabla^2\Phi\approx 1.2\times10^{-6}\ \text{s}^{-2}=\frac{c^2}{2}\,R_{00} \]

Reading. The same physical number can be read as a Newtonian tidal divergence (\(4\pi G\bar\rho\)) or as the \(00\)-Ricci curvature (\(R_{00}=\nabla^2\Phi/c^2\)); the coupling \(8\pi G/c^4\) is exactly what makes these two readings numerically identical.

Units check. \([4\pi G\rho]=(\text{m}^3\text{kg}^{-1}\text{s}^{-2})(\text{kg}\,\text{m}^{-3})=\text{s}^{-2}\), the dimension of \(\nabla^2\Phi\). Consistent.

Problems
  1. (A) Recover the redshift factor. A clock sits at radius \(r\) in the weak field of mass \(M\). Using \(g_{00}=-(1+2\Phi/c^2)\) and \(d\tau=\sqrt{-g_{00}}\,dt\), show the fractional rate difference between \(r\) and infinity to first order.
    Solution \(d\tau/dt=\sqrt{1+2\Phi/c^2}\approx 1+\Phi/c^2\) with \(\Phi=-GM/r\). Hence \(\dfrac{d\tau}{dt}\approx 1-\dfrac{GM}{rc^2}\). A lower clock (more negative \(\Phi\)) runs slow; for Earth's surface, \(GM/rc^2\approx (6.674\times10^{-11})(5.972\times10^{24})/[(6.371\times10^{6})(8.988\times10^{16})]\approx 6.95\times10^{-10}\), the gravitational redshift GPS must correct for.
  2. (A) Perturbation at Earth's surface. Compute \(h_{00}=2GM_\oplus/(R_\oplus c^2)\) with \(M_\oplus=5.972\times10^{24}\ \text{kg}\), \(R_\oplus=6.371\times10^{6}\ \text{m}\).
    Solution Numerator \(2GM_\oplus=2(6.674\times10^{-11})(5.972\times10^{24})=7.97\times10^{14}\). Denominator \(R_\oplus c^2=(6.371\times10^{6})(8.988\times10^{16})=5.726\times10^{23}\). Thus \(h_{00}=7.97\times10^{14}/5.726\times10^{23}\approx 1.39\times10^{-9}\) — a thousand times weaker than the Sun's surface value, deep in the weak-field regime.
  3. (B) The factor-of-two trap. Suppose a student writes \(R_{\mu\nu}=\kappa'\,T_{\mu\nu}\) (no trace reversal) and demands the Newtonian limit. What \(\kappa'\) results, and why is the theory nonetheless inconsistent?
    Solution With \(R_{00}=\kappa' T_{00}=\kappa'\rho c^2\) and \(R_{00}=\nabla^2\Phi/c^2\), one gets \(\nabla^2\Phi=\kappa'\rho c^4\). Matching \(4\pi G\rho\) gives \(\kappa'=4\pi G/c^4\), half the correct value. The theory fails because \(\nabla^\mu R_{\mu\nu}\ne0\) in general (only \(\nabla^\mu G_{\mu\nu}=0\) by the Bianchi identity), so \(R_{\mu\nu}=\kappa'T_{\mu\nu}\) is incompatible with local energy–momentum conservation \(\nabla^\mu T_{\mu\nu}=0\). Trace reversal is mandatory, and it is exactly what restores the missing factor of two.
  4. (B) Trace of the field equations. Starting from \(G_{\mu\nu}=\kappa T_{\mu\nu}\), take the trace with \(g^{\mu\nu}\) (in 4 dimensions) to show \(R=-\kappa T\), and hence write the trace-reversed form.
    Solution \(g^{\mu\nu}G_{\mu\nu}=g^{\mu\nu}(R_{\mu\nu}-\tfrac12 g_{\mu\nu}R)=R-\tfrac12(4)R=R-2R=-R\). Setting equal to \(\kappa g^{\mu\nu}T_{\mu\nu}=\kappa T\) gives \(-R=\kappa T\), i.e. \(R=-\kappa T\). Substituting back, \(R_{\mu\nu}=G_{\mu\nu}+\tfrac12 g_{\mu\nu}R=\kappa T_{\mu\nu}+\tfrac12 g_{\mu\nu}(-\kappa T)=\kappa(T_{\mu\nu}-\tfrac12 g_{\mu\nu}T)\), the form used in step 8.
  5. (C) Pressure correction to the source. For a perfect fluid \(T_{\mu\nu}=(\rho c^2+p)u_\mu u_\nu/c^2+p\,g_{\mu\nu}\), redo the trace-reversed \(00\) source in the slow-motion limit and show the Newtonian source becomes \(\rho+3p/c^2\).
    Solution Slow motion: \(u_\mu u_\nu\to(u_0)^2\delta^0_\mu\delta^0_\nu\) with \((u_0)^2=c^2\). Then \(T_{00}=(\rho c^2+p)+p\,g_{00}\approx(\rho c^2+p)-p=\rho c^2\), and the trace \(T=g^{\mu\nu}T_{\mu\nu}=(\rho c^2+p)(-1)+4p=-\rho c^2+3p\). Hence \(T_{00}-\tfrac12 g_{00}T=\rho c^2-\tfrac12(-1)(-\rho c^2+3p)=\rho c^2-\tfrac12\rho c^2+\tfrac32 p=\tfrac12(\rho c^2+3p)\). Matching \(R_{00}=\nabla^2\Phi/c^2=\kappa\cdot\tfrac12(\rho c^2+3p)\) with \(\kappa=8\pi G/c^4\) gives \(\nabla^2\Phi=4\pi G\!\left(\rho+\dfrac{3p}{c^2}\right)\). Pressure gravitates — the origin of the Tolman active mass and a driver of relativistic stellar collapse.