Black Holes and General Relativity

Introduction to General Relativity

General relativity (GR)—Einstein's 1915 theory of gravitation—describes gravity not as a force but as the curvature of four-dimensional spacetime caused by mass and energy. The Einstein field equations G_mu_nu = 8*pi*G/c^4 * T_mu_nu relate the Einstein tensor (describing spacetime geometry via the metric g_mu_nu and its derivatives) to the stress-energy tensor T_mu_nu (describing the distribution of energy, momentum, and stress in matter and fields). "Matter tells spacetime how to curve; curved spacetime tells matter how to move" (John Wheeler's summary). GR reduces to Newtonian gravity in the weak-field, slow-motion limit, but produces qualitatively new predictions: time dilation in gravitational fields, perihelion precession of Mercury (43 arcseconds/century—unexplained by Newtonian gravity), gravitational deflection of light (Eddington 1919 solar eclipse observation confirming 1.75 arcsecond deflection), gravitational red shift, frame dragging, gravitational waves, black holes, and an expanding universe.

GR has passed every experimental test with extraordinary precision: gravitational wave detection by LIGO (2015) confirmed GR's prediction with waveforms matching black hole merger simulations to 1%; gravitational time dilation in GPS satellites (clocks run 45 microseconds/day faster in weaker field, 7 microseconds/day slower from special relativistic time dilation—net +38 microseconds/day) is corrected in GPS algorithms; strong-field tests from binary pulsar PSR B1913+16 (orbital energy loss from gravitational wave emission matching GR to 0.2%—Nobel 1993 Hulse, Taylor); PPN (parameterised post-Newtonian) parameter gamma = 1.000021±0.000023 from Cassini spacecraft radio signal Shapiro delay measurement—all confirming metric theory of gravity to very high precision.

Black Holes

Stellar Black Holes and Schwarzschild Metric

Black holes form when sufficient mass is concentrated within the Schwarzschild radius r_S = 2GM/c^2—for a solar mass, r_S = 2.95 km (below 10 km neutron star radius and only reached by stellar remnants above ~3 solar masses after supernova). The Schwarzschild metric (exterior vacuum solution): ds^2 = -(1-r_S/r)c^2*dt^2 + (1-r_S/r)^-1*dr^2 + r^2*(d theta^2 + sin^2(theta)*d phi^2)—describes the exact geometry outside a non-rotating spherically symmetric mass. Event horizon at r = r_S: a null surface beyond which no light or matter can escape to infinity—not a physical singularity (locally flat Minkowski spacetime for freely falling observer) but a coordinate singularity in Schwarzschild coordinates (removed in Kruskal-Szekeres extension). Gravitational time dilation at r: d tau/dt = sqrt(1 - r_S/r)—clocks run slower in stronger gravity, stopping at r_S. Photon sphere at r = 3GM/c^2 = 1.5 r_S: circular photon orbits—unstable, but relevant to black hole shadow image. Innermost stable circular orbit (ISCO) for Schwarzschild BH: r_ISCO = 6GM/c^2 = 3*r_S—matter spiralling inward releases ~6% of rest mass energy before reaching ISCO (vs ~42% for maximally spinning Kerr BH)—fundamental limit on X-ray binary accretion luminosity efficiency.

Rotating Black Holes and Kerr Geometry

Real astrophysical black holes rotate—the Kerr metric (1963) describes rotating black holes characterised by mass M and angular momentum J = a*M*c (spin parameter a, 0 ≤ a ≤ GM/c^2). The ergosphere—region between outer event horizon r+ and static limit r_ergosphere = GM/c^2 + sqrt((GM/c^2)^2 - a^2 cos^2 theta) where spacetime rotation drags all observers to co-rotate—enables the Penrose process (theoretical): particle splitting in ergosphere can extract rotational energy, with outgoing fragment having more energy than the original. Frame dragging (Lense-Thirring effect): rotating mass drags spacetime around it, causing precession of gyroscopes—measured by Gravity Probe B (Stanford, NASA, 2004-2011): geodetic precession 6606.1 mas/yr (vs GR prediction 6606.1), frame-dragging 37.2 mas/yr (vs GR 39.2)—confirmed to 19% accuracy for frame-dragging. Magnetically arrested disk and SANE accretion flows around maximally spinning Kerr black holes are thought to power relativistic jets—Blandford-Znajek mechanism extracting BH spin energy electromagnetically via force-free magnetosphere threading the event horizon.

Gravitational Waves

LIGO and Gravitational Wave Astronomy

Gravitational waves—ripples in spacetime metric propagating at c—are emitted by accelerating asymmetric masses: inspiral and merger of compact binaries (black holes, neutron stars) are the strongest sources. LIGO (Laser Interferometer Gravitational-Wave Observatory): two L-shaped detectors (each arm 4 km) in Washington State and Louisiana; laser Michelson interferometry measures differential arm length changes delta L/L ~ 10^-21 (strain h)—equivalent to measuring 1/1000 proton diameter over 4 km. GW150914 (September 14, 2015, first detection): two ~36 and ~29 solar mass black holes merging at ~410 Mpc, peak gravitational wave strain h ~10^-21, peak luminosity ~3.6×10^56 W (exceeding entire observable universe electromagnetic luminosity for ~0.2 seconds)—confirmed GR's 100-year-old prediction (Nobel 2017 Weiss, Barish, Thorne). LIGO-Virgo-KAGRA network (O3, O4 observing runs 2020-2023): ~200 confirmed compact binary merger events in GWTC catalogues—including neutron star mergers (GW170817: first multi-messenger event with gamma-ray burst, kilonova optical counterpart, X-ray, radio—revealing r-process heavy element nucleosynthesis in neutron star merger). Future: LISA (Laser Interferometer Space Antenna, 2.5 Mkm arm ESA mission, launch ~2034): detects millihertz GW sources (supermassive BH mergers at cosmological distances, extreme mass ratio inspirals, galactic compact binary population)—opening low-frequency GW astronomy window.

Cosmology and Expanding Universe

Friedmann-Lemaître-Robertson-Walker (FLRW) metric describes a homogeneous, isotropic expanding universe: ds^2 = -c^2*dt^2 + a(t)^2*(dr^2/(1-kr^2) + r^2*d Omega^2) where a(t) = scale factor. Friedmann equations from GR applied to FLRW metric: (dot a/a)^2 = H^2 = 8*pi*G/3*rho - kc^2/a^2 + Lambda*c^2/3 (H = Hubble parameter; rho = energy density; k = spatial curvature; Lambda = cosmological constant). Hubble constant H_0 = 67-73 km/s/Mpc (Hubble tension: CMB-based Planck value 67.4 vs. Cepheid/SN 73.0—~5 sigma discrepancy potentially indicating new physics). Cosmic inflation: exponential expansion in first ~10^-36 to ~10^-32 seconds smooths curvature, establishes causal flatness, and generates primordial density fluctuations from quantum vacuum to seeds of large-scale structure. Dark energy (68%): accelerating expansion (SNe Ia Nobel 1998) attributed to cosmological constant Lambda or dynamical scalar field—equation of state w = p/(rho*c^2) ≈ -1 consistent to date with all probes (DES, Planck, BAO). JWST (James Webb Space Telescope): discovering fully formed massive galaxies at z>10 (400 Myr after Big Bang)—challenging standard model of galaxy formation and dark matter halo assembly, potentially requiring modifications to early-universe physics.

Examples and Applications

Example 1: First Image of a Black Hole

Event Horizon Telescope (EHT) collaboration—global network of 8 millimetre/submillimetre radio telescope arrays combined by very-long-baseline interferometry (VLBI) to achieve Earth-diameter baseline resolution of ~20 microarcseconds—captured the first image of a black hole shadow in M87 (April 2019). M87*: mass ~6.5 billion solar masses at 55 million light-years, event horizon diameter ~23 billion km—angular resolution required ~40 microarcseconds at 230 GHz (1.3 mm wavelength). Black hole shadow: central dark region (photon capture cross-section ~27*sqrt(3)*GM/c^2 in Schwarzschild) surrounded by bright photon ring from lensed emission—size, shape, and brightness asymmetry consistent with Kerr GR predictions favouring retrograde jet relative to disk rotation. First image of Sagittarius A* (Milky Way center, 4 million solar masses, 2022): variability timescale ~8 minutes (comparable to image integration time) complicated reconstruction from time-averaged images—revealing ring morphology consistent with GR predictions for a Kerr BH. Next generation EHT (ngEHT): adding more stations (Atacama 345 GHz, Greenland, Africa) improving dynamic range and enabling time-resolved imaging of accretion flow dynamics and jet formation mechanism.

Example 2: Neutron Stars and Pulsars

Neutron stars—remnants of massive star supernovae with typically 1.2-2.2 solar masses in ~10 km radius—represent the densest observable matter (central density 3-8 times nuclear saturation density rho_0 = 2.3×10^17 kg/m³). Pulsars: rapidly rotating neutron stars emitting radio beams aligned with magnetic axis—Jodrell Bank, Parkes, FAST telescope timing of pulse arrival times to nanosecond precision tracks orbital and spin dynamics enabling precision tests of GR and dense matter EOS. Double pulsar PSR J0737-3039 (Burgay et al. 2003): first known double pulsar system—five GR effects measured including Shapiro delay, orbital decay from GW emission, and geodetic precession; system orbital decay matches GR to 1.3×10^-4 accuracy—the most stringent GR test in strong-field regime. Neutron star equation of state: NICER (NASA X-ray telescope on ISS) mass and radius measurements of PSR J0030 and J0437 via X-ray pulse profile modelling at ~15%/5% precision; combined with GW170817 tidal deformability constraint Lambda_1.4 < 580—constraining dense matter EOS in 2-4 rho_0 range. Magnetars: ultra-highly magnetised neutron stars (B ~ 10^11 T) emitting X-ray/gamma giant flares from sudden magnetic field reconfiguration—December 2004 giant flare (SGR 1806-20) briefly exceeded total solar luminosity by 10^6×, detected by every operating gamma-ray detector in the solar system.

Example 3: GPS and Relativistic Corrections

Global Positioning System (GPS) satellite clocks experience two opposing GR effects requiring precise corrections for accurate navigation. Special relativistic time dilation: GPS satellites orbit at v = 3.87 km/s—clocks run slow by 7 microseconds/day relative to ground clocks (moving clocks tick slower by sqrt(1-v^2/c^2) ≈ 1 - v^2/(2c^2)). General relativistic gravitational time dilation: GPS orbit at 20,200 km altitude has weaker gravitational field—clocks run fast by 45 microseconds/day relative to ground (clocks at higher gravitational potential tick faster—gravitational blueshift). Net effect: +38 microseconds/day (GPS satellite clocks run fast). Without correction, GPS position would drift ~11 km/day. Correction: GPS satellite clocks are pre-programmed to run 38 microseconds/day slower before launch. Additionally, the Sagnac effect (Earth rotation in non-rotating frame—receivers on rotating Earth are non-inertial)—up to 133 nanoseconds correction for equatorial paths—is incorporated. If GR were neglected in GPS design, position errors would accumulate within hours making GPS useless—a daily lived demonstration of applied relativistic physics in every smartphone navigation application globally.

Example 4: Gravitational Lensing

Light follows geodesics in curved spacetime—massive objects, galaxies, and galaxy clusters bend light rays from background sources (gravitational lensing). Weak lensing: statistical measurement of correlated shape distortions (shear) of background galaxies by foreground large-scale structure—maps dark matter distribution projected along line of sight; Dark Energy Survey (DES) and Rubin LSST will map billions of galaxy shapes. Strong lensing: point like sources (quasar, gamma-ray burst) lensed by galaxy into Einstein ring or multiple images when source-lens-observer alignment within Einstein radius theta_E = sqrt(4GM*D_LS/c^2/D_L/D_S)—Abell 2744 galaxy cluster forms 67 multiple images of 5 background galaxies enabling high-precision mass map to 1% accuracy. Microlensing: stellar-mass lensing of background stars (OGLE survey toward Galactic bulge): bright amplification events lasting days to weeks from planetary mass objects, stellar mass BHs, neutron stars—discovery of hundreds of exoplanets (OGLE method); confirmed isolated stellar mass black holes via astrometric microlensing (HST measurement of deflected source position in addition to brightness). Time delay between multiple images of lensed quasars constrains Hubble constant H0 (H0LiCOW/TDCOSMO programme—TDC-derived H0 = 73.3 +2.4/-2.3 km/s/Mpc consistent with distance ladder, ~2 sigma above CMB value).

Example 5: Hawking Radiation and Black Hole Thermodynamics

Stephen Hawking (1974) showed that quantum field theory in curved spacetime predicts black holes emit thermal radiation—Hawking radiation—from pair production near the event horizon. Virtual particle-antiparticle pairs created from quantum vacuum fluctuations near horizon: one particle falls in (negative energy in Killing frame), the other escapes as real particle—black hole loses mass at temperature T_H = hbar*c^3/(8*pi*G*M*k_B) = 6.17×10^-8 * (M_sun/M) K. For solar mass black hole T_H ~ 60 nanokelvin—undetectable against CMB (2.7 K); primordial microscopic black holes M ~ 10^12 kg would have T_H ~ 10^11 K and evaporate today. Hawking radiation spectra thermal with greybody factors from photon potential barrier near horizon. Information paradox: Hawking radiation appears perfectly thermal—encoding no information about infalling matter—violating unitarity of quantum mechanics (pure states evolving to mixed states). Resolutions: black hole complementarity (Susskind: no observer sees both in-fall and radiation); firewall proposal (AMPS 2012: infalling observer encounters Planck-energy firewall at horizon to preserve unitarity); and ER=EPR (Maldacena, Susskind 2013: wormhole geometry related to quantum entanglement)—one of the most actively debated problems at the intersection of GR and quantum mechanics as of 2025.

Example 6: Pulsar Timing Arrays

Pulsar Timing Arrays (PTAs)—collections of millisecond pulsars monitored as a Galactic-scale gravitational wave detector—detect nanohertz gravitational waves from supermassive black hole binary systems inspiral long before merger. Millisecond pulsars (MSPs): spun up via accretion in binary systems to ~300-700 Hz rotation rates; pulse timing stability rivalling atomic clocks (timing residuals 100 ns over years). GW background from SMBH binary population: a gravitational wave passes between the pulsar and Earth—changes effective pulsar distance—introducing correlated timing residuals across sky with Hellings-Downs angular correlation function characterising the isotropic GW background signature. PTA results (2023): NANOGrav (North American Nanohertz Observatory for Gravitational Waves) 15-year data set, combined with PPTA (Parkes), EPTA (European), InPTA (Indian), and CPTA (Chinese)—all independently found evidence of gravitational wave background signal with strain h_c ~ 2.4×10^-15 at nHz frequencies and Hellings-Downs angular correlation—marking the probable detection of the nanohertz GW background from SMBH binaries in galaxy merger remnants throughout the universe. IPTA (International PTA)—combining all PTA datasets—currently characterising spectral index and searching for individual resolvable SMBH binary sources in PTA data.

Example 7: Compact Object Mergers

Compact object mergers—binary black hole (BBH), binary neutron star (BNS), and neutron star-black hole (NSBH) coalescence—are the primary gravitational wave source population detected by LIGO-Virgo-KAGRA. GW waveform encodes: inspiral chirp (frequency and amplitude sweep from 10's Hz to kHz in minutes-seconds), merger, and ringdown (quasinormal modes of final Kerr black hole decaying exponentially—spectrum uniquely identifying final mass and spin consistent with no-hair theorem tests). BBH population: mass spectrum from ~5 to ~100 solar masses; lower mass gap (3-5 M_sun) and upper mass gap (~55-130 M_sun, pair-instability supernova) in black hole formation distribution. GW170817 multimessenger: LIGO BNS detection → GRB 170817A (fermi, INTEGRAL, 1.7 second delay probing Lorentz invariance violation to 10^-16 enhancement factor per unit energy) → AT2017gfo kilonova optical transient decaying from blue (lanthanide-poor) to red (lanthanide-rich) over days—confirming r-process production of ~0.05 solar masses of heavy elements (strontium spectroscopically confirmed in kilonova spectrum by absorption feature). Merging neutron star remnant (HMNS or direct BH) depends on total mass relative to maximum TOV mass—observationally constraining nuclear EOS above nuclear saturation density.

Example 8: Tests of General Relativity

GR has been tested across 17 orders of magnitude in gravitational potential and 9 orders of magnitude in curvature—from solar system (weak field) to double pulsars (moderate strong field) to BH mergers (extreme strong field). Parametrised post-Newtonian (PPN) framework: ten dimensionless parameters characterising any metric theory of gravity—GR uniquely predicts gamma = beta = 1, xi = alpha_1 = alpha_2 = alpha_3 = zeta_1 = zeta_2 = zeta_3 = zeta_4 = 0; all measured to be consistent with GR values at 10^-4 to 10^-6 precision through solar system experiments (Shapiro delay, lunar laser ranging, gyroscope precession). Strong equivalence principle: gravitational self-energy of compact object contributes equally to inertial and gravitational mass—tested to 10^-3 from pulsar timing in wide binary systems (Nordtvedt effect). Gravitational wave tests: no-hair theorem (inspiral GW spectrum plus ringdown quasinormal modes constrain BH mass and spin; consistency between inspiral and ringdown masses tests Kerr hypothesis); polarisation measurement (only two tensor modes in GR; constraints from multi-detector network on vector/scalar polarisation modes); and graviton mass upper limit m_g < 1.27×10^-23 eV/c^2 from GW dispersion—most stringent graviton mass bound, all consistent with GR massless spin-2 graviton prediction.

Try it live

Everything above runs in your browser — open Black Holes & Spacetime Curvature and change the parameters while it is running. Nothing is installed, nothing is uploaded, the whole model lives in one tab.