← ArticlesGeneral Relativity: Geometry, Fields and Applications: /20 SummaryEngineering · PhysicsLesson 12/13← PrevNext →
GuidePublished 14 Aug 202626 min readBy Kevin JoginPhysicsApplied Classical PhysicsGeneral Relativity: GeometryFields and Applications

Engineering · Physics · Applied Classical Physics

General Relativity: Geometry, Fields and Applications: /20 Summary

Engineering handbook for general relativity: geometry, fields and applications, covering context and scope, /20 summary, from special to general relativity.

Executive summary

This handbook section converts the supplied engineering material into a practical, source-controlled reference. It concentrates on the following learning outcomes.

Context and scope
/20 Summary
From Special to General Relativity
The Equivalence Principle — Three Versions
Three Predictions Before the Theory
Gravitational Redshift

Context and scope

the supplied physics reference, Applications of Classical Physics — Chapters 24–28 From SR to GR · Curved Spacetime · Stars & Black Holes · Gravitational Waves · Cosmology



/20 Summary

  • General Relativity = geometric theory of gravity. Gravitation is not a force — it's the curvature of 4-D spacetime sourced by energy-momentum.
  • Two-line summary (Wheeler):

    Matter tells spacetime how to curve; spacetime tells matter how to move.

  • Einstein's field equations: GαβRαβ12gαβR=8πGc4TαβG_{\alpha\beta} \equiv R_{\alpha\beta} - \tfrac{1}{2}g_{\alpha\beta}R = \frac{8\pi G}{c^4} T_{\alpha\beta} 10 nonlinear PDEs for the metric gαβg_{\alpha\beta}. The single equation that determines stars, black holes, gravitational waves, and the universe.
  • The Equivalence Principle (EP) is the conceptual foundation: locally, free-fall is indistinguishable from inertial motion in SR. Gravity = local consequence of curvature.
  • Geodesics are the GR analog of straight lines — paths that parallel-transport their own tangent vector. Free particles follow timelike geodesics; light follows null geodesics.
  • Four physical regimes, each a chapter:
    1. Spherical stars — Schwarzschild metric, Tolman-Oppenheimer-Volkoff equation, neutron stars + Chandrasekhar limit.
    2. Black holes — Schwarzschild (static) + Kerr (rotating); event horizons; ergospheres; ISCO; thermodynamics + Hawking radiation.
    3. Gravitational waves — ripples in gαβg_{\alpha\beta} from accelerating masses; LIGO/Virgo/KAGRA/LISA detections from 2015 onward.
    4. Cosmology — FRW metric for homogeneous-isotropic universe; Friedmann equations; ΛCDM (Ωm0.32\Omega_m \approx 0.32, ΩΛ0.68\Omega_\Lambda \approx 0.68).
  • Tested to extreme precision: Mercury perihelion, light bending, Shapiro delay, gravitational redshift, frame-dragging (Gravity Probe B), pulsar timing (Hulse-Taylor 0.2% on Ṗ\dot P), GW150914 inspiral-merger-ringdown.


Master Map

mindmap
  root((GR))
    Ch.24 SR→GR
      Equivalence principle
        Weak EP
        Einstein EP
        Strong EP
      Free fall = inertial
      Gravitational redshift
      Light bending
      Spacetime curvature intro
      Newtonian limit
    Ch.25 Fundamentals
      Manifold + tangent space
      Metric g_αβ
      Christoffel Γ
      Covariant deriv ∇
      Parallel transport
      Geodesic
      Riemann R^α_βγδ
      Ricci R_αβ, R
      Einstein G_αβ
      Field equations
      Bianchi identity
      Killing vectors
      Conservation
    Ch.26 Stars & BH
      Schwarzschild metric
      TOV equation
      White dwarfs
      Neutron stars
      Schwarzschild BH
        Event horizon
        Singularity
        Kruskal
      Kerr BH
        Ergosphere
        ISCO
        Penrose process
      BH thermodynamics
        T_Hawking
        S_BH = A/4
        Four laws
    Ch.27 GW & Tests
      Linearized GR
      TT gauge
      Plus / cross polarizations
      Quadrupole formula
      Binary inspirals
      LIGO/Virgo/KAGRA/LISA
      Tests of GR
        Mercury perihelion
        Light bending
        Shapiro delay
        Frame dragging
    Ch.28 Cosmology
      Cosmological principle
      FRW metric
      Friedmann eqs
      Hubble H₀
      Critical density
      ΛCDM
      Distance measures
      Redshift
      CMB
      BBN
      Inflation
      Structure formation
      Dark matter / energy


From Special to General Relativity


The Equivalence Principle — Three Versions

Version Statement Implication
Weak EP (WEP) Inertial mass = gravitational mass All bodies fall identically in vacuum
Einstein EP (EEP) WEP + local SR + no preferred frame Local physics in free-fall = SR
Strong EP (SEP) EEP + gravitational binding energy also obeys WEP Self-gravitating bodies also fall identically

Experimental tests of WEP: Eötvös experiments (10⁻⁹), MicroSCOPE (10⁻¹⁵), Lunar Laser Ranging.

Key consequence (EEP): every local experiment in a freely-falling lab reproduces special-relativistic physics. Gravity is not a force at a point — it's an effect that appears over finite extents through tidal effects (geodesic deviation).


Three Predictions Before the Theory

Using only EEP + SR, Einstein derived three observable consequences. Each has now been verified to high precision.


Gravitational Redshift

Photon climbing out of potential Φ\Phi loses energy. Frequency shift:

Δνν=ΔΦc2\frac{\Delta\nu}{\nu} = -\frac{\Delta\Phi}{c^2}

(positive ΔΦ\Delta\Phi = higher potential = redshift). Pound-Rebka 1959: 22.6 m tower, Δν/ν2.5×1015\Delta\nu/\nu \sim 2.5\times 10^{-15}, agreement to 1%. Modern: atomic clocks at different elevations, 10⁻¹⁸ precision.


Light Deflection

Newtonian (mass with v=cv = c): Δθ=2GM/(bc2)\Delta\theta = 2GM/(b c^2). GR doubles this:

Δθ=4GMbc2\Delta\theta = \frac{4GM}{bc^2}

For Sun grazing incidence (b=Rb = R_\odot): Δθ=1.75\Delta\theta = 1.75''. Eddington 1919 eclipse expedition — confirmed GR, made Einstein world-famous. Modern radio interferometry: 10410^{-4} precision.


Time Dilation in Gravitational Potential

A clock at potential Φ\Phi ticks slower than one at potential 0:

dτ=dt1+2Φ/c2d\tau = dt\sqrt{1 + 2\Phi/c^2}

GPS satellites must correct for this: 45μ\sim 45\,\mus/day "gain" (gravity) vs 7μ7\,\mus/day "loss" (motion, SR). Without correction, position errors 10\sim 10 km/day.


Spacetime as Curved Manifold

Einstein's insight (1907–1915): the EEP demands gravity be encoded in the geometry of spacetime itself, not added as a force.

Tidal effects between nearby free-falling observers reveal local curvature (geodesic deviation, §25.6). A small-enough region looks flat (EEP); but extended structure cannot be transformed away.


Newtonian Limit

For weak fields and slow motion, write gαβ=ηαβ+hαβg_{\alpha\beta} = \eta_{\alpha\beta} + h_{\alpha\beta} with |h|1|h|\ll 1. The 00-component dominates:

g00=(1+2Φ/c2),Φ=GM/r (point mass)g_{00} = -(1 + 2\Phi/c^2), \quad \Phi = -GM/r \text{ (point mass)}

The geodesic equation reduces to Newton's: ẍi=iΦ\ddot x^i = -\partial^i \Phi

Confirming GR contains Newtonian gravity in the appropriate limit.



Fundamentals of General Relativity


Spacetime Manifold

A 4-D manifold is a topological space locally homeomorphic to 4\mathbb{R}^4, with smooth coordinate charts. Tangent space TpT_p\mathcal{M} at point pp is the linear space of vectors at pp.

Basis vectors: 𝐞α=/xα\mathbf{e}_\alpha = \partial/\partial x^\alpha in chart {xα}\{x^\alpha\}. Dual basis: 𝐞α=dxα\mathbf{e}^\alpha = dx^\alpha, with 𝐞α,𝐞β=δβα\langle \mathbf{e}^\alpha, \mathbf{e}_\beta\rangle = \delta^\alpha_\beta.


Metric Tensor

The metric gαβg_{\alpha\beta} defines:

ds2=gαβdxαdxβds^2 = g_{\alpha\beta}\, dx^\alpha dx^\beta

  • Signature (,+,+,+)(-,+,+,+) (B&T convention).
  • Symmetric: gαβ=gβαg_{\alpha\beta} = g_{\beta\alpha} ⇒ 10 independent components in 4-D.
  • Has inverse gαβg^{\alpha\beta} with gαγgγβ=δβαg^{\alpha\gamma}g_{\gamma\beta} = \delta^\alpha_\beta.
  • Lowers/raises indices: Vα=gαβVβV_\alpha = g_{\alpha\beta} V^\beta.

Examples:

Metric gαβg_{\alpha\beta} What it describes
Minkowski ηαβ=diag(1,1,1,1)\eta_{\alpha\beta} = \text{diag}(-1,1,1,1) Flat (SR)
Schwarzschild gtt=(12M/r)g_{tt} = -(1-2M/r), grr=(12M/r)1g_{rr} = (1-2M/r)^{-1}, gθθ=r2g_{\theta\theta} = r^2, gϕϕ=r2sin2θg_{\phi\phi} = r^2\sin^2\theta Outside spherical mass
Kerr (Boyer-Lindquist) (complicated; see Ch. 26) Rotating BH
FRW c2dt2+a(t)2γijdxidxj-c^2 dt^2 + a(t)^2 \gamma_{ij}\,dx^i dx^j Homog./isotropic universe
pp-wave (gravitational wave) 2dudv+dx2+dy2+hij(u)xixjdu2-2du\,dv + dx^2 + dy^2 + h_{ij}(u)\,x^i x^j\, du^2 Plane GW

Christoffel Symbols (Connection Coefficients)

The covariant derivative requires a way to compare vectors at different points — the connection. For the Levi-Civita connection (metric-compatible, torsion-free):

Γβγα=12gαδ(βgδγ+γgδβδgβγ)\boxed{\Gamma^\alpha_{\beta\gamma} = \tfrac{1}{2}g^{\alpha\delta}(\partial_\beta g_{\delta\gamma} + \partial_\gamma g_{\delta\beta} - \partial_\delta g_{\beta\gamma})}

Properties:

  • Symmetric in lower indices: Γβγα=Γγβα\Gamma^\alpha_{\beta\gamma} = \Gamma^\alpha_{\gamma\beta}.
  • Not a tensor — transforms inhomogeneously under coordinate changes (extra 2x\partial^2 x term).
  • Vanishes at any point in a freely-falling frame (Riemann normal coordinates).

Covariant Derivative

For a vector field VαV^\alpha:

βVα=βVα+ΓβγαVγ\nabla_\beta V^\alpha = \partial_\beta V^\alpha + \Gamma^\alpha_{\beta\gamma}V^\gamma

For a covector VαV_\alpha:

βVα=βVαΓβαγVγ\nabla_\beta V_\alpha = \partial_\beta V_\alpha - \Gamma^\gamma_{\beta\alpha}V_\gamma

For higher-rank tensors: one +Γ+\Gamma for each upper index, Γ-\Gamma for each lower.

Key facts:

  • γgαβ=0\nabla_\gamma g_{\alpha\beta} = 0 (metric-compatibility).
  • \nabla is a derivation: Leibniz rule.
  • Produces tensor of one higher rank from tensor of given rank.

Notation: sometimes write β=;β\nabla_\beta = ;_\beta and β=,β\partial_\beta = ,_\beta: Vα;β=Vα,β+ΓβγαVγV^\alpha{}_{;\beta} = V^\alpha{}_{,\beta} + \Gamma^\alpha_{\beta\gamma}V^\gamma.


Parallel Transport & Geodesics

Parallel transport of VαV^\alpha along curve xα(λ)x^\alpha(\lambda):

DVαdλdxβdλβVα=dVαdλ+ΓβγαdxβdλVγ=0\frac{D V^\alpha}{d\lambda} \equiv \frac{dx^\beta}{d\lambda}\nabla_\beta V^\alpha = \frac{dV^\alpha}{d\lambda} + \Gamma^\alpha_{\beta\gamma}\frac{dx^\beta}{d\lambda}V^\gamma = 0

Geodesic equation (curve that parallel-transports its own tangent):

d2xαdτ2+Γβγαdxβdτdxγdτ=0\boxed{\frac{d^2 x^\alpha}{d\tau^2} + \Gamma^\alpha_{\beta\gamma}\frac{dx^\beta}{d\tau}\frac{dx^\gamma}{d\tau} = 0}

— with τ\tau proper time (timelike) or affine parameter λ\lambda (null/spacelike).

Variational form: geodesics extremize dτ=gαβẋαẋβdλ\int d\tau = \int\sqrt{-g_{\alpha\beta}\dot x^\alpha\dot x^\beta}\,d\lambda.

Path type gαβẋαẋβg_{\alpha\beta}\dot x^\alpha\dot x^\beta Particle
Timelike <0<0 Massive (max τ\tau)
Null =0=0 Photon
Spacelike >0>0 (not a physical worldline)

Riemann Curvature Tensor

Measures failure of parallel transport around a closed loop, or equivalently the failure of αβ=βα\nabla_\alpha\nabla_\beta = \nabla_\beta\nabla_\alpha:

[γ,δ]Vα=RαβγδVβ[\nabla_\gamma, \nabla_\delta] V^\alpha = R^\alpha{}_{\beta\gamma\delta} V^\beta

Explicit formula:

Rαβγδ=γΓβδαδΓβγα+ΓμγαΓβδμΓμδαΓβγμ\boxed{R^\alpha{}_{\beta\gamma\delta} = \partial_\gamma\Gamma^\alpha_{\beta\delta} - \partial_\delta\Gamma^\alpha_{\beta\gamma} + \Gamma^\alpha_{\mu\gamma}\Gamma^\mu_{\beta\delta} - \Gamma^\alpha_{\mu\delta}\Gamma^\mu_{\beta\gamma}}


Symmetries (drastically reduce independent components)

  1. Antisymmetry first pair: Rαβγδ=RβαγδR_{\alpha\beta\gamma\delta} = -R_{\beta\alpha\gamma\delta}
  2. Antisymmetry last pair: Rαβγδ=RαβδγR_{\alpha\beta\gamma\delta} = -R_{\alpha\beta\delta\gamma}
  3. Symmetry under pair exchange: Rαβγδ=RγδαβR_{\alpha\beta\gamma\delta} = R_{\gamma\delta\alpha\beta}
  4. First Bianchi: Rα[βγδ]=0R_{\alpha[\beta\gamma\delta]} = 0
  5. Second Bianchi: [ϵRαβ]γδ=0\nabla_{[\epsilon} R_{\alpha\beta]\gamma\delta} = 0

In 4-D: 20 independent components.


Geodesic Deviation Equation

Tidal force = relative acceleration of neighboring geodesics:

D2ξαdτ2=Rαβγδuβuγξδ\frac{D^2\xi^\alpha}{d\tau^2} = -R^\alpha{}_{\beta\gamma\delta}\, u^\beta u^\gamma \xi^\delta

where ξα\xi^\alpha = separation vector, uαu^\alpha = 4-velocity. This is the precise sense in which curvature = gravity: tidal forces are observable, frame-independent.


Ricci Tensor, Ricci Scalar, Einstein Tensor

Ricci tensor (contraction of Riemann): Rαβ=RγαγβR_{\alpha\beta} = R^\gamma{}_{\alpha\gamma\beta}

10 independent components in 4-D; symmetric.

Ricci (scalar) curvature: R=gαβRαβR = g^{\alpha\beta} R_{\alpha\beta}

Einstein tensor: GαβRαβ12gαβR\boxed{G_{\alpha\beta} \equiv R_{\alpha\beta} - \tfrac{1}{2}g_{\alpha\beta}R}

— combination chosen so that αGαβ=0\nabla^\alpha G_{\alpha\beta} = 0 identically (contracted Bianchi), matching conservation αTαβ=0\nabla^\alpha T_{\alpha\beta} = 0.


Einstein's Field Equations

Gαβ+Λgαβ=8πGc4Tαβ\boxed{G_{\alpha\beta} + \Lambda g_{\alpha\beta} = \frac{8\pi G}{c^4} T_{\alpha\beta}}

with cosmological constant Λ\Lambda (originally added by Einstein for static universe; later discovered to be needed for accelerating expansion).

Properties:

  • 10 nonlinear coupled PDEs for gαβg_{\alpha\beta}.
  • 4 of them are constraints (initial-value): G0βG^0{}_\beta involves only first time derivatives.
  • Remaining 6 are dynamical; with diffeomorphism freedom (4 functions) ⇒ 2 physical DOF per spacetime point. Two physical DOF = two polarizations of gravitational waves.
  • Reduce to Newton: g00(1+2Φ/c2)g_{00} \approx -(1 + 2\Phi/c^2), G00/c22ΦG_{00}/c^2 \to \nabla^2\Phi, T00/c2ρc2T_{00}/c^2 \to \rho c^2, giving 2Φ=4πGρ\nabla^2\Phi = 4\pi G\rho.

Lagrangian formulation (Einstein-Hilbert action): S=c416πGRgd4x+SmatterS = \frac{c^4}{16\pi G}\int R\,\sqrt{-g}\,d^4x + S_{\text{matter}}

Variation w.r.t. gαβg_{\alpha\beta} yields the field equations.


Stress-Energy Tensor in Curved Spacetime

Same forms as flat-space (Part I, Ch. 2), with ηαβgαβ\eta_{\alpha\beta} \to g_{\alpha\beta} and \partial \to \nabla:

Matter TαβT^{\alpha\beta}
Perfect fluid (ρ+P/c2)uαuβ+Pgαβ(\rho + P/c^2)u^\alpha u^\beta + P g^{\alpha\beta}
Dust (pressureless) ρuαuβ\rho u^\alpha u^\beta
EM field (1/4π)(FαμFβμ14gαβFμνFμν)(1/4\pi)(F^{\alpha\mu}F^\beta{}_\mu - \tfrac{1}{4}g^{\alpha\beta}F^{\mu\nu}F_{\mu\nu})
Cosmological constant Λc4/(8πG)gαβ-\Lambda c^4/(8\pi G)\,g^{\alpha\beta}
Scalar field ϕ\phi αϕβϕgαβ[12(ϕ)2+V(ϕ)]\nabla^\alpha\phi\nabla^\beta\phi - g^{\alpha\beta}[\tfrac{1}{2}(\nabla\phi)^2 + V(\phi)]

Conservation: βTαβ=0\nabla_\beta T^{\alpha\beta} = 0 — guarantees local energy-momentum conservation (with all gravitational effects via Γ\Gamma).


Symmetries & Killing Vectors

A Killing vector ξα\xi^\alpha generates a spacetime symmetry; satisfies Killing's equation:

αξβ+βξα=0\nabla_\alpha\xi_\beta + \nabla_\beta\xi_\alpha = 0

Conserved quantity along geodesics:

ξαpα=const\xi^\alpha p_\alpha = \text{const}

where pα=muαp^\alpha = m u^\alpha (timelike) or photon 4-momentum (null).

Killing vector Conservation
Time-translation t\partial_t (Schwarzschild) Energy at infinity EE
Rotation ϕ\partial_\phi (axisymmetric) Angular momentum LL
Boosts (Minkowski) Center-of-mass position
FRW spatial Killing Comoving 3-momentum


Relativistic Stars and Black Holes


The Schwarzschild Solution

Unique spherically symmetric vacuum solution of Einstein's equations (Birkhoff's theorem — even time-dependent spherical solutions reduce to Schwarzschild outside):

ds2=(1rsr)c2dt2+dr21rs/r+r2(dθ2+sin2θdϕ2)\boxed{ds^2 = -\left(1 - \frac{r_s}{r}\right)c^2 dt^2 + \frac{dr^2}{1 - r_s/r} + r^2(d\theta^2 + \sin^2\theta\,d\phi^2)}

with Schwarzschild radius:

rs=2GMc2r_s = \frac{2GM}{c^2}

Numerical values:

  • Earth: rs=8.87r_s = 8.87 mm
  • Sun: rs=2.95r_s = 2.95 km
  • 1M1 M_\odot neutron star: rs=2.95r_s = 2.95 km, actual radius ~10 km
  • Sgr A* (4.1×106M4.1\times 10^6 M_\odot): rs=1.2×1010r_s = 1.2\times 10^{10} m ≈ 0.08 AU
  • M87 (6.5×109M6.5\times 10^9 M_\odot): rs=1.9×1013r_s = 1.9\times 10^{13} m ≈ 130 AU

Coordinate singularities at r=rsr = r_s (removable; just the event horizon) and at θ=0,π\theta = 0, \pi (familiar polar issue). True curvature singularity at r=0r = 0 (RαβγδRαβγδ1/r6R_{\alpha\beta\gamma\delta}R^{\alpha\beta\gamma\delta} \propto 1/r^6).


Tests in Schwarzschild Geometry


Perihelion Precession

For bound timelike orbits, periapse advances by

Δϕ=6πGMc2a(1e2) per orbit\Delta\phi = \frac{6\pi GM}{c^2 a(1-e^2)}\text{ per orbit}

Mercury: observed 43''/century unexplained by Newtonian (after extracting other planets' pulls); Einstein's prediction matched immediately — pivotal 1915 confirmation.


Light Bending

Null geodesic with impact parameter bb deflects by

Δϕ=4GMbc2+O((GM/bc2)2)\Delta\phi = \frac{4GM}{bc^2} + O((GM/bc^2)^2)


Shapiro Delay (1964)

Round-trip light travel through gravitational potential is delayed:

Δt=4GMc3ln4rErRb2\Delta t = \frac{4GM}{c^3}\ln\frac{4 r_E r_R}{b^2}

Tested: Cassini probe (2002), agreement 105\sim 10^{-5}.


Stable Circular Orbits

  • Schwarzschild ISCO at r=6GM/c2=3rsr = 6GM/c^2 = 3 r_s.
  • Photon sphere at r=3GM/c2=1.5rsr = 3GM/c^2 = 1.5 r_s (unstable circular photon orbit).
  • No stable circular orbits inside ISCO; matter spirals in.

Inside the Star — TOV Equation

For static, spherically symmetric perfect fluid:

dPdr=G(ρ+P/c2)(M(r)+4πr3P/c2)r2(12GM(r)/rc2)\boxed{\frac{dP}{dr} = -\frac{G\,(\rho + P/c^2)\,(M(r) + 4\pi r^3 P/c^2)}{r^2(1 - 2GM(r)/rc^2)}}

with dM/dr=4πr2ρdM/dr = 4\pi r^2 \rho. Tolman-Oppenheimer-Volkoff (TOV) equation is the GR generalization of Newtonian hydrostatic equilibrium dP/dr=GMρ/r2dP/dr = -G M\rho/r^2.

Three GR corrections (all >0> 0, all destabilizing — make stars less stable):

  1. P/c2P/c^2 adds to mass-energy density.
  2. Pressure itself gravitates (+P+P term in numerator).
  3. Spatial curvature in denominator.

Equilibrium Stellar Sequences

Compact object Pressure source Mass limit rs/Rr_s/R
White dwarf Electron degeneracy $M_{\rm Ch} \approx 1.4 M_\odot$ 104\sim 10^{-4}
Neutron star Nucleon degeneracy + strong int. 2\sim 23M3 M_\odot (TOV-limited) 0.3\sim 0.3
Black hole None (gravitational collapse) - 1

Chandrasekhar limit (1931): for fully degenerate, ultra-relativistic electron gas (Pρ4/3P\propto\rho^{4/3}), no stable hydrostatic solution above

$$M_{\rm Ch} = \frac{0.7732 (\hbar c)^{3/2}}{G^{3/2} m_p^2 \mu_e^2} \approx 1.456 \left(\frac{2}{\mu_e}\right)^2 M_\odot$$

For μe=2\mu_e = 2 (carbon/oxygen WD): $M_{\rm Ch} \approx 1.44 M_\odot$.

Type Ia supernovae detonate as a WD accretes to near $M_{\rm Ch}$ — standard candles, used to discover dark energy (1998).

TOV limit for neutron stars: depends on EOS; modern limit 2.3M\sim 2.3 M_\odot. Observations (PSR J0740+6620: 2.08±0.07M2.08\pm 0.07 M_\odot) constrain dense-matter EOS sharply.


Black Holes — Schwarzschild

Event horizon at r=rsr = r_s:

  • One-way membrane: timelike worldlines inside cannot exit.
  • gttg_{tt} and grrg_{rr} swap signs ⇒ inside, rr becomes timelike and tt spacelike. Falling toward r=0r=0 is as inevitable as moving forward in time.
  • Singularity at r=0r = 0 is a future endpoint in spacetime for all interior observers.

Eddington-Finkelstein coordinates remove horizon singularity: ds2=(12M/r)dv2+2dvdr+r2dΩ2ds^2 = -(1 - 2M/r)dv^2 + 2 dv\, dr + r^2 d\Omega^2 with advanced time v=t+r+2Mln|r/2M1|v = t + r + 2M\ln|r/2M - 1|. Smooth across horizon.

Kruskal extension reveals full maximal extension: includes the BH, an exterior region, a "white hole," and a second asymptotic exterior. Wormhole structure (not physically traversable for Schwarzschild).

Mass formula: the only parameter — Schwarzschild BH is fully characterized by mass MM alone. "No-hair" theorem: stationary BH in Einstein-Maxwell theory characterized by (M,J,Q)(M, J, Q) only — independent of how it formed.


Kerr Black Holes (Rotating)

The unique stationary axisymmetric vacuum BH solution. Two parameters: mass MM and angular momentum JJ. Spin parameter a=J/(Mc)a = J/(Mc).

Boyer-Lindquist coordinates: ds2=(12GMrΣc2)c2dt24GMarsin2θΣc2cdtdϕ+ΣΔdr2+Σdθ2+(r2+a2+2GMa2rsin2θΣc2)sin2θdϕ2ds^2 = -\left(1 - \frac{2GMr}{\Sigma c^2}\right)c^2 dt^2 - \frac{4 GMar\sin^2\theta}{\Sigma c^2}\,c\, dt\,d\phi + \frac{\Sigma}{\Delta}\,dr^2 + \Sigma\, d\theta^2 + \left(r^2 + a^2 + \frac{2GMa^2 r\sin^2\theta}{\Sigma c^2}\right)\sin^2\theta\, d\phi^2

with Σ=r2+a2cos2θ\Sigma = r^2 + a^2\cos^2\theta, Δ=r2rsr+a2\Delta = r^2 - r_s r + a^2.

Region Defined by Property
Outer horizon r+=M+M2a2r_+ = M + \sqrt{M^2 - a^2} (geom. units) Event horizon
Inner horizon r=MM2a2r_- = M - \sqrt{M^2 - a^2} Cauchy horizon (unstable)
Ergosphere r+<r<rE=M+M2a2cos2θr_+ < r < r_E = M + \sqrt{M^2 - a^2\cos^2\theta} Frame dragging forbids stationarity
Ring singularity r=0r = 0 in equatorial plane True singularity

Extremal limit: aMa \to M, r+=r=Mr_+ = r_- = M. For a>Ma > M: naked singularity (cosmic censorship hypothesis: physically excluded).


Ergosphere & Penrose Process

Within ergosphere, no static observer (frame-dragging). Particle splits into two; one falls in with E<0E < 0, other escapes with Eout>EinE_{\text{out}} > E_{\text{in}}. Net extracts rotational energy from BH (up to ~29% of mass for extremal Kerr).

Astrophysical relevance: Blandford-Znajek mechanism — electromagnetic Penrose-like extraction powering AGN jets, GRBs.


Innermost Stable Circular Orbit (ISCO)

Depends on spin and direction:

Spin a/Ma/M Prograde ISCO Retrograde ISCO
0 (Schwarzschild) 6GM/c26 GM/c^2 6GM/c26 GM/c^2
0.998 (near-extremal) 1.24GM/c21.24 GM/c^2 9GM/c29 GM/c^2
aMa \to M (extremal) GM/c2GM/c^2 9GM/c29 GM/c^2

Maximum binding energy at ISCO:

  • Schwarzschild: 5.7%\sim 5.7\% of rest mass.
  • Maximally rotating prograde: 42%\sim 42\%.
  • Compare nuclear burning H→He: 0.7%0.7\%. Accretion onto BHs is the most efficient energy source in nature.

Black Hole Thermodynamics

Bekenstein (1972) + Hawking (1974): BHs are thermodynamic objects.


Four Laws of BH Mechanics

Law Statement
0th Surface gravity κ\kappa constant over horizon
1st dM=(κ/8πG)dA+ΩHdJ+ΦHdQdM = (\kappa/8\pi G)\, dA + \Omega_H\, dJ + \Phi_H\, dQ
2nd dA0dA \ge 0 (classical) — area theorem
3rd Cannot reach κ=0\kappa = 0 (extremal) in finite steps

with horizon area AA, angular velocity ΩH\Omega_H, electric potential ΦH\Phi_H. Looks exactly like ordinary thermodynamics if TκT \sim \kappa, SAS \sim A.


Hawking Temperature

Quantum field theory in BH background gives thermal radiation at:

TH=κ2πckB=c38πGMkB\boxed{T_H = \frac{\hbar \kappa}{2\pi c k_B} = \frac{\hbar c^3}{8\pi G M k_B}}

For solar-mass BH: TH6×108T_H \approx 6\times 10^{-8} K — far below CMB; no net evaporation.

Evaporation time: tevapM3×1026t_{\text{evap}} \sim M^3 \times 10^{-26} yr/kg³. Solar BH: 106710^{67} yr. Primordial BHs 1015\sim 10^{15} g (small asteroid): evaporating now, end with gamma-ray burst.


Bekenstein-Hawking Entropy

SBH=kBc3A4G=kBA4P2\boxed{S_{BH} = \frac{k_B c^3 A}{4 G \hbar} = \frac{k_B A}{4\ell_P^2}}

with Planck length P=G/c31.6×1035\ell_P = \sqrt{G\hbar/c^3} \approx 1.6\times 10^{-35} m.

Holographic principle: entropy of any system bounded by area, not volume. Foundational hint at quantum gravity (AdS/CFT, etc.).


Generalized Second Law

d(Smatter+SBH)0d(S_{\text{matter}} + S_{BH}) \ge 0

Combined ordinary + BH entropy non-decreasing. Restores 2nd law: ordinary entropy can decrease (e.g., matter falls in) only if BH entropy increases more.



Gravitational Waves & Experimental Tests


Linearized GR

Background flat + small perturbation: gαβ=ηαβ+hαβg_{\alpha\beta} = \eta_{\alpha\beta} + h_{\alpha\beta}, |h|1|h| \ll 1.

Define trace-reversed: hαβ=hαβ12ηαβh\bar h_{\alpha\beta} = h_{\alpha\beta} - \tfrac{1}{2}\eta_{\alpha\beta} h.

In Lorenz gauge αhαβ=0\partial^\alpha \bar h_{\alpha\beta} = 0:

hαβ=16πGc4Tαβ\boxed{\Box \bar h_{\alpha\beta} = -\frac{16\pi G}{c^4} T_{\alpha\beta}}

— wave equation with source. In vacuum, plane-wave solutions exist.


Transverse-Traceless (TT) Gauge

Residual gauge freedom fixes: h0α=0\bar h^{0\alpha} = 0, hii=0h^i{}_i = 0. In TT gauge, hαβTT=hαβTT\bar h^{TT}_{\alpha\beta} = h^{TT}_{\alpha\beta} has only two independent components corresponding to the two polarizations.

For wave traveling in zz-direction:

hαβTT=(00000h+h×00h×h+00000)cos(ωtkz)h^{TT}_{\alpha\beta} = \begin{pmatrix} 0 & 0 & 0 & 0 \\ 0 & h_+ & h_\times & 0 \\ 0 & h_\times & -h_+ & 0 \\ 0 & 0 & 0 & 0 \end{pmatrix}\cos(\omega t - kz)

Polarization Effect on ring of test masses
h+h_+ (plus) Stretches along xx, squeezes along yy, then vice versa
h×h_\times (cross) Same pattern rotated 45°

Both polarizations propagate at cc (in GR; constrained to |cGW/c1|<1015|c_{GW}/c - 1| < 10^{-15} by GW170817 + GRB).


Quadrupole Formula — Generation

Far-field radiation from a source:

hijTT=2Grc4Q̈ijTT(tr/c)\boxed{h^{TT}_{ij} = \frac{2G}{rc^4}\,\ddot Q^{TT}_{ij}(t - r/c)}

where Qij=(xixj13δijr2)ρd3xQ_{ij} = \int(x_i x_j - \tfrac{1}{3}\delta_{ij}r^2)\rho\,d^3x is the traceless mass quadrupole moment.

No monopole (mass conservation forbids it). No dipole (momentum conservation: 1st time derivative of mass dipole = momentum, conserved → 2nd derivative vanishes). GWs are intrinsically quadrupolar.


Power Radiated

P=G5c5QijQijP = \frac{G}{5c^5}\,\langle\dddot Q_{ij}\dddot Q^{ij}\rangle

The factor G/c52.7×1053G/c^5 \approx 2.7\times 10^{-53} W⁻¹ explains why GWs are so weak. To radiate 1 W requires Q21053\dddot Q^2 \sim 10^{53} in SI units — an enormous accelerating quadrupole.

For binary with masses m1,m2m_1, m_2, separation aa, circular orbit: P=32G45c5(m1m2)2(m1+m2)a5P = \frac{32 G^4}{5 c^5}\,\frac{(m_1 m_2)^2(m_1 + m_2)}{a^5}

Order of magnitude for systems:

Source Power (W)
Earth orbiting Sun 200\sim 200
Hulse-Taylor binary pulsar 7×1024\sim 7\times 10^{24}
Final inspiral of 30M30 M_\odot BHs 1049\sim 10^{49} (briefly outshines all stars in observable universe combined!)

Binary Inspiral

Orbit decay drives the two bodies together. Chirp signal — increasing frequency and amplitude.


Hulse-Taylor Pulsar PSR B1913+16 (1974)

Ṗorbobs/ṖorbGR=1.0013±0.0021\dot P_{\text{orb}}^{\text{obs}}/\dot P_{\text{orb}}^{\text{GR}} = 1.0013 \pm 0.0021

Indirect confirmation, Nobel 1993 (Hulse, Taylor). Best test of GR in strong-field, dynamic regime — until LIGO.


Chirp Mass

For binary observable through GW, the combination

=(m1m2)3/5(m1+m2)1/5\mathcal{M} = \frac{(m_1 m_2)^{3/5}}{(m_1 + m_2)^{1/5}}

is the chirp mass — directly extracted from GW frequency evolution. Sets amplitude and frequency drift.


Three Stages of Compact Binary Coalescence

flowchart LR
    A[Inspiral<br/>Post-Newtonian] --> B[Merger<br/>Numerical Relativity]
    B --> C[Ringdown<br/>Quasinormal modes]
Phase Physics Method
Inspiral Slow PN expansion Analytic
Merger Strong-field, dynamical Numerical relativity
Ringdown Perturbed BH Quasinormal mode spectroscopy

Direct Detection


LIGO — Laser Interferometer Gravitational-wave Observatory

  • 4 km Michelson interferometers (Hanford, Livingston)
  • Sensitivity h1021h \sim 10^{-21} — fractional length change ΔL/L1022\Delta L / L \sim 10^{-22} over 4 km = 4×10194\times 10^{-19} m ≈ 1/10000 proton diameter!
  • First detection: GW150914 (Sept 14, 2015): merger of 30M\sim 30 M_\odot + 35M\sim 35 M_\odot BHs → 62M62 M_\odot + radiated 3Mc25×10473 M_\odot c^2 \approx 5\times 10^{47} J in GWs.
  • Nobel 2017 (Weiss, Barish, Thorne).

Multi-Messenger: GW170817

Binary neutron-star merger detected coincidently with gamma-ray burst (GRB 170817A) and kilonova (AT 2017gfo) across the EM spectrum. Established:

  • Speed of gravity = cc to 101510^{-15} precision.
  • Origin of heavy r-process elements (gold, platinum etc.).
  • Independent H0H_0 measurement (standard siren).

Current and Future Detectors

Detector Type Frequency band Sensitivity / status
LIGO (US, ×2) 4 km laser interferometer 10 Hz – kHz ~150 events (O4 ongoing)
Virgo (Italy) 3 km laser same Operating
KAGRA (Japan) Cryogenic, underground same Operating
LIGO India Planned, 2030\sim 2030 same
LISA (ESA/NASA) 2.5 Mkm space interferometer mHz Launch 2035; SMBH mergers, EMRIs
Einstein Telescope 10 km underground 1 Hz – kHz Planned, 2035\sim 2035
Cosmic Explorer 40 km surface same Concept
Pulsar timing arrays (NANOGrav, EPTA) Indirect via $\Delta t_{\rm arr}$ nHz SMBH binaries; first evidence 2023

Tests of GR (Beyond GW)

Test Prediction Verified to
Light bending (Sun) 1.751.75'' 10410^{-4} (VLBI)
Mercury perihelion 43.043.0''/century 10310^{-3}
Gravitational redshift Δν/ν=gΔz/c2\Delta\nu/\nu = g\Delta z/c^2 10410^{-4} (Pound-Rebka), 101810^{-18} (clocks)
Shapiro delay Cassini probe 10510^{-5}
Lunar laser ranging (EP, Strong EP) WEP & SEP 101310^{-13} (WEP); 10410^{-4} (SEP)
Geodetic precession Gravity Probe B 0.3%0.3\%
Frame dragging (Lense-Thirring) Gravity Probe B 20%\sim 20\%
Hulse-Taylor $\dot P_{\rm orb}$ GR quadrupole formula 0.2%0.2\%
GW propagation speed $ c_{GW} - c
Black hole shadow Photon sphere ratio 10%\sim 10\% (Event Horizon Telescope)

No verified deviations from GR to date. Despite decades of attempts at modified gravity (MOND, f(R)f(R), scalar-tensor, ...), Einstein wins.



Cosmology


The Cosmological Principle

On scales 100\gtrsim 100 Mpc, the universe is:

  1. Homogeneous (same density everywhere)
  2. Isotropic (same in every direction)

Verified by CMB temperature uniformity (ΔT/T105\Delta T/T \sim 10^{-5}) and large-scale galaxy distribution.


The FRW Metric

Most general metric consistent with cosmological principle:

ds2=c2dt2+a(t)2[dr21kr2+r2dΩ2]\boxed{ds^2 = -c^2 dt^2 + a(t)^2\left[\frac{dr^2}{1-kr^2} + r^2 d\Omega^2\right]}

with scale factor a(t)a(t) and curvature constant k{1,0,+1}k \in \{-1, 0, +1\} (open, flat, closed).

Comoving coordinates: galaxies (on average) at fixed (r,θ,ϕ)(r,\theta,\phi); expansion is encoded in a(t)a(t).

Hubble parameter:

$H(t)ȧa,H0=H(t0)67H(t) \equiv \frac{\dot a}{a}, \quad H_0 = H(t_0) \approx 677373 km/s/Mpc \approx 2.2\times 10^{-18},\text{s}^{-1}$$

(Persistent tension between "early universe" H067.4H_0 \approx 67.4 from Planck/CMB and "late" H073.0H_0 \approx 73.0 from local Cepheid+SNIa — known as the Hubble tension.)


Friedmann Equations

Plugging FRW into Einstein equations with perfect-fluid source:

H2=8πG3ρkc2a2+Λc23\boxed{H^2 = \frac{8\pi G}{3}\rho - \frac{k c^2}{a^2} + \frac{\Lambda c^2}{3}}

äa=4πG3(ρ+3Pc2)+Λc23\boxed{\frac{\ddot a}{a} = -\frac{4\pi G}{3}\left(\rho + \frac{3 P}{c^2}\right) + \frac{\Lambda c^2}{3}}

Conservation equation (continuity): ρ̇+3H(ρ+P/c2)=0\dot\rho + 3 H(\rho + P/c^2) = 0

(Any two of the three equations are independent.)


The Cosmic Energy Budget

Define density parameters Ωi=ρi/ρc\Omega_i = \rho_i/\rho_c where ρc=3H2/(8πG)\rho_c = 3 H^2/(8\pi G). Friedmann becomes:

Ωm+Ωr+ΩΛ+Ωk=1\Omega_m + \Omega_r + \Omega_\Lambda + \Omega_k = 1

Best-fit ΛCDM (Planck 2018):

Component Ω\Omega today Equation of state w=P/(ρc2)w = P/(\rho c^2) ρa?\rho \propto a^?
Matter (mm): baryon + dark 0.315±0.0070.315\pm0.007 0 a3a^{-3}
Radiation (rr): photons + neutrinos 9×1059\times 10^{-5} 1/31/3 a4a^{-4}
Dark energy (Λ\Lambda) 0.685±0.0070.685\pm 0.007 1-1 (so far) a0a^0 (constant)
Curvature (kk) 0.001±0.0020.001\pm 0.002 1/3-1/3 a2a^{-2}
Baryons (within Ωm\Omega_m) 0.0493±0.00020.0493\pm 0.0002 00 a3a^{-3}
Cold dark matter (within Ωm\Omega_m) 0.265±0.0070.265\pm 0.007 00 a3a^{-3}

Hubble distance: DH=c/H04.3D_H = c/H_0 \approx 4.3 Gpc.


Cosmic History

flowchart LR
    A[t=0<br/>Big Bang] --> B[Inflation<br/>t~10⁻³⁶ s]
    B --> C[Radiation-dom<br/>BBN at t~minutes]
    C --> D[Matter-dom<br/>z~3400]
    D --> E[Recombination<br/>z~1090, t~380 kyr<br/>CMB!]
    E --> F[Dark Ages]
    F --> G[Reionization<br/>z~7]
    G --> H[Structure formation]
    H --> I[Λ-dom<br/>z~0.4]
    I --> J[Today<br/>t=13.8 Gyr]
Epoch Redshift Time Temperature Physics
Planck 104310^{-43} s 103210^{32} K (Quantum gravity)
Inflation 103610^{-36} s Exponential expansion; δρ/ρ\delta\rho/\rho from quantum fluctuations
Electroweak 101210^{-12} s 101510^{15} K Higgs sym. breaking
Quark-hadron 10610^{-6} s 101210^{12} K Quark confinement
Big Bang Nucleosynthesis 4×108\sim 4\times 10^8 3 min 10910^9 K 4^4He, D, 3^3He, 7^7Li forge
Matter-radiation equality 3400 50 kyr 9000 K ρm=ρr\rho_m = \rho_r
Recombination 1090 380 kyr 3000 K e+pe^- + p \to H; CMB released
Dark ages 30–1090 10710^710810^8 yr 30–3000 K No starlight
Reionization ~7 10910^9 yr ~30 K First stars/quasars; H reionized
Λ takes over 0.4 9 Gyr 4 K ρΛ=ρm\rho_\Lambda = \rho_m
Now 0 13.8 Gyr 2.725 K CMB temperature

Engineering use and verification

State the model, coordinate system, assumptions, boundary conditions and validity range before using an equation. Track dimensions and sign conventions through each derivation, test limiting cases, and distinguish mathematical possibility from physical realisability. Where a model informs engineering design, compare it with measurement or a second method and quantify the effect of idealisations rather than hiding them inside numerical precision.

  • Confirm scope, assumptions, interfaces and required outcome.
  • Check dimensions, sign conventions, boundary conditions and limiting cases.
  • Identify current project, customer and regulatory requirements.
  • Separate source examples from mandatory acceptance criteria.
  • Check calculations, tables and selections by an independent method.
  • Verify safety, maintainability and credible failure modes.
  • Record evidence, revisions, approvals and unresolved limitations.
  • Validate the result under representative operating conditions.

Continue learning

Plasma Physics: Kinetics, Fields and Waves: Final TakeawaysGuide · PhysicsNEXT LESSON →General Relativity: Geometry, Fields and Applications: Distance MeasuresGuide · PhysicsPlasma Physics: Kinetics, Fields and Waves: /20 SummaryGuide · PhysicsFluid Dynamics: Conservation, Flow and Stability: Petschek ReconnectionGuide · Physics