Part 8: Mathematical Relativity

Tensor Calculus

General Relativity requires a mathematical language that works in any coordinate system — one where the laws of physics look the same regardless of whether we use Cartesian, spherical, or any arbitrary coordinates. This language is tensor calculus. A tensor is a geometric object that transforms in a specific, predictable way under coordinate transformations.

The simplest tensors are familiar objects:

• Scalar (rank-0 tensor): A single number invariant under coordinate changes (e.g., temperature, rest mass)
• Vector (rank-1 tensor): Has one index — V^\mu (contravariant) or V_\mu (covariant). Transforms with one factor of the coordinate transformation matrix
• Rank-2 tensor: Has two indices — T^{\mu\nu}, T_{\mu\nu}, or T^\mu{}_\nu. The metric tensor g_{\mu\nu} is the most important example

The transformation law for a rank-2 covariant tensor under a coordinate change x^\mu \to x^{\mu'} is:

Einstein's summation convention is crucial: any index that appears once up and once down in a term is summed over all spacetime dimensions (0,1,2,3). So A^\mu B_\mu means \sum_{\mu=0}^{3} A^\mu B_\mu = A^0 B_0 + A^1 B_1 + A^2 B_2 + A^3 B_3.

The Metric Tensor

The metric tensor g_{\mu\nu} is the fundamental object of GR. It defines the geometry of spacetime by specifying how to measure distances and angles. The infinitesimal spacetime interval (proper distance) between two nearby events is:

In flat Minkowski spacetime (special relativity), the metric is simply the diagonal matrix:

This gives the familiar interval ds^2 = -c^2 dt^2 + dx^2 + dy^2 + dz^2 (using the "mostly plus" sign convention). In curved spacetime, g_{\mu\nu} becomes a function of position — it varies from point to point, encoding the gravitational field.

Key properties of the metric tensor:

• Symmetric: g_{\mu\nu} = g_{\nu\mu} (so it has 10 independent components in 4D, not 16)
• Non-degenerate: \det(g_{\mu\nu}) \neq 0 (invertible everywhere except at singularities)
• Index raising/lowering: V_\mu = g_{\mu\nu} V^\nu and V^\mu = g^{\mu\nu} V_\nu
• Inverse: g^{\mu\alpha} g_{\alpha\nu} = \delta^\mu_\nu (Kronecker delta)

For example, the metric of a 2-sphere of radius R in coordinates (\theta, \phi) is:

So g_{\theta\theta} = R^2, g_{\phi\phi} = R^2 \sin^2\theta, and all off-diagonal components are zero. This simple metric already contains curvature information — the sphere is not flat, and that fact is encoded in how g_{\phi\phi} depends on \theta.

Christoffel Symbols & the Covariant Derivative

In curved spacetime, ordinary partial derivatives of tensors do not produce tensors — they pick up extra terms due to the coordinate system changing from point to point. We need the covariant derivative \nabla_\mu, which accounts for curvature and produces proper tensors.

The Christoffel symbols (also called the connection coefficients) \Gamma^\lambda_{\mu\nu} are the "correction terms" that make the covariant derivative work. They are defined entirely in terms of the metric:

This formula, while intimidating at first glance, has a clear structure: each Christoffel symbol is built from first derivatives of the metric tensor. In flat spacetime with Cartesian coordinates, all \Gamma^\lambda_{\mu\nu} = 0. In curved spacetime or curvilinear coordinates, they encode how basis vectors change from point to point.

The covariant derivative of a vector is:

The geodesic equation — the equation of motion for a free particle in curved spacetime — is expressed using Christoffel symbols:

This says: a freely falling particle follows the straightest possible path (geodesic) through curved spacetime. The Christoffel symbols act like a "gravitational force" term — they are the mathematical manifestation of gravity in GR.

import numpy as np

# ============================================================
# COMPUTING CHRISTOFFEL SYMBOLS FOR A 2-SPHERE
# Demonstrates how the metric determines the connection
# ============================================================

# For a 2-sphere of radius R, coordinates are (theta, phi)
# Metric: g_ab = diag(R^2, R^2 * sin^2(theta))
# We compute Gamma^a_bc from the metric and its derivatives

R = 1.0  # unit sphere for simplicity

def metric_sphere(theta):
    """Return the metric tensor g_ab for a 2-sphere at given theta."""
    g = np.zeros((2, 2))
    g[0, 0] = R**2                      # g_theta_theta
    g[1, 1] = R**2 * np.sin(theta)**2   # g_phi_phi
    return g

def metric_inverse_sphere(theta):
    """Return the inverse metric g^ab for a 2-sphere."""
    g_inv = np.zeros((2, 2))
    g_inv[0, 0] = 1.0 / R**2
    g_inv[1, 1] = 1.0 / (R**2 * np.sin(theta)**2)
    return g_inv

def christoffel_symbols_sphere(theta, dtheta=1e-8):
    """
    Compute Christoffel symbols Gamma^a_bc numerically.
    
    Gamma^a_bc = (1/2) g^{ad} (dg_{db}/dx^c + dg_{dc}/dx^b - dg_{bc}/dx^d)
    
    For 2-sphere, only theta-derivatives are non-trivial (metric doesn't depend on phi).
    """
    dim = 2  # theta=0, phi=1
    
    # Numerical derivatives of metric w.r.t. theta (index 0)
    # dg/dphi = 0 by symmetry
    g_plus = metric_sphere(theta + dtheta)
    g_minus = metric_sphere(theta - dtheta)
    dg_dtheta = (g_plus - g_minus) / (2 * dtheta)
    
    # dg/dphi = 0 (metric is phi-independent)
    dg = np.zeros((2, 2, 2))  # dg[a][b][c] = dg_ab / dx^c
    dg[:, :, 0] = dg_dtheta   # derivative w.r.t. theta
    # dg[:, :, 1] = 0         # derivative w.r.t. phi (already zero)
    
    # Compute Christoffel symbols
    g_inv = metric_inverse_sphere(theta)
    gamma = np.zeros((2, 2, 2))  # Gamma^a_bc
    
    for a in range(dim):
        for b in range(dim):
            for c in range(dim):
                total = 0.0
                for d in range(dim):
                    # (1/2) g^{ad} * (dg_{db}/dx^c + dg_{dc}/dx^b - dg_{bc}/dx^d)
                    total += 0.5 * g_inv[a, d] * (
                        dg[d, b, c] + dg[d, c, b] - dg[b, c, d]
                    )
                gamma[a, b, c] = total
    
    return gamma

# Compute at theta = pi/4 (45 degrees)
theta_test = np.pi / 4
gamma = christoffel_symbols_sphere(theta_test)

# Known analytical results for comparison
gamma_analytical = {
    'Gamma^theta_phi_phi': -np.sin(theta_test) * np.cos(theta_test),
    'Gamma^phi_theta_phi': np.cos(theta_test) / np.sin(theta_test),
    'Gamma^phi_phi_theta': np.cos(theta_test) / np.sin(theta_test),
}

print("=" * 60)
print("CHRISTOFFEL SYMBOLS FOR A 2-SPHERE (R = 1)")
print(f"Evaluated at theta = pi/4 = {theta_test:.4f} rad")
print("=" * 60)

labels = [['theta', 'phi'], ['theta', 'phi']]
coord_names = ['theta', 'phi']

print(f"\n{'Symbol':<30} {'Numerical':<15} {'Analytical':<15}")
print("-" * 60)
for a in range(2):
    for b in range(2):
        for c in range(b, 2):  # symmetric in lower indices
            val = gamma[a, b, c]
            if abs(val) > 1e-10:
                name = f"Gamma^{coord_names[a]}_{coord_names[b]}_{coord_names[c]}"
                # Get analytical value
                analytical = gamma_analytical.get(name, 0.0)
                print(f"  {name:<28} {val:<15.8f} {analytical:<15.8f}")

print(f"\n--- Physical Interpretation ---")
print(f"  Gamma^theta_phi_phi = -sin(theta)cos(theta) = {gamma[0,1,1]:.6f}")
print(f"    -> Lines of latitude are NOT geodesics (they 'accelerate' in theta)")
print(f"  Gamma^phi_theta_phi = cot(theta) = {gamma[1,0,1]:.6f}")
print(f"    -> The phi-basis vector rotates as you move in theta")
print(f"\n  At the equator (theta=pi/2):")
gamma_eq = christoffel_symbols_sphere(np.pi / 2)
print(f"    Gamma^theta_phi_phi = {gamma_eq[0,1,1]:.6f} (= -sin(90)cos(90) = 0)")
print(f"    Gamma^phi_theta_phi = {gamma_eq[1,0,1]:.6f} (= cot(90) = 0)")
print(f"    -> At the equator, both corrections vanish (locally 'flat' directions)")

Differential Geometry

The proper mathematical setting for GR is a differentiable manifold — a space that looks locally like \mathbb{R}^n (flat Euclidean space) but may have non-trivial global structure and curvature. Spacetime is a 4-dimensional manifold equipped with a Lorentzian metric.

Manifolds & Tangent Spaces

A manifold is a topological space covered by overlapping coordinate patches (charts), with smooth transition functions between overlapping regions. The key idea: you cannot cover a sphere with a single coordinate system without singularities (the poles), but you can cover it with two or more overlapping charts.

At each point p on the manifold, there exists a tangent space T_p M — a vector space of all possible "velocities" at that point. Vectors live in tangent spaces, not on the manifold itself. A tangent vector is formally defined as a directional derivative operator acting on smooth functions.

Parallel Transport

In flat space, comparing vectors at different points is trivial — just move them. On a curved manifold, there is no unique way to compare vectors at different points. Parallel transport provides one prescription: move a vector along a curve while keeping it "as parallel as possible" according to the connection.

The condition for parallel transport of a vector V^\mu along a curve with tangent u^\alpha = dx^\alpha/d\lambda is:

On a curved surface, parallel-transporting a vector around a closed loop brings it back rotated — the rotation angle is directly related to the enclosed curvature. This is the geometric meaning of the Riemann tensor.

flowchart TD
    A["Metric Tensor
g_μν
10 components"] --> B["Christoffel Symbols
Γ^λ_μν
40 independent"]
    B --> C["Riemann Tensor
R^ρ_σμν
20 independent"]
    C --> D["Ricci Tensor
R_μν = R^λ_μλν
10 components"]
    C --> E["Ricci Scalar
R = g^μν R_μν
1 component"]
    D --> F["Einstein Tensor
G_μν = R_μν - ½g_μν R
10 components"]
    E --> F
    F --> G["Einstein Field Equations
G_μν = 8πG/c⁴ T_μν"]

    style A fill:#3B9797,color:#fff
    style B fill:#16476A,color:#fff
    style C fill:#132440,color:#fff
    style D fill:#16476A,color:#fff
    style E fill:#16476A,color:#fff
    style F fill:#BF092F,color:#fff
    style G fill:#BF092F,color:#fff
                            

Curvature Tensors

Curvature is what distinguishes GR from special relativity. The Riemann curvature tensor R^\rho{}_{\sigma\mu\nu} is the master curvature object — it contains all information about how spacetime is curved. It measures how much a vector rotates when parallel-transported around an infinitesimal closed loop.

The Riemann Tensor

Defined in terms of Christoffel symbols:

In 4D spacetime, the Riemann tensor has 4^4 = 256 components, but symmetries reduce the independent components to just 20:

• Antisymmetry in last two indices: R^\rho{}_{\sigma\mu\nu} = -R^\rho{}_{\sigma\nu\mu}
• Antisymmetry in first pair (with lowered index): R_{\rho\sigma\mu\nu} = -R_{\sigma\rho\mu\nu}
• Pair symmetry: R_{\rho\sigma\mu\nu} = R_{\mu\nu\rho\sigma}
• First Bianchi identity: R^\rho{}_{[\sigma\mu\nu]} = 0

Ricci Tensor & Ricci Scalar

The Ricci tensor is obtained by contracting (summing over) one index of the Riemann tensor:

The Ricci scalar (scalar curvature) is the trace of the Ricci tensor:

The Ricci scalar is a single number at each point that measures the overall "amount" of curvature. For a 2-sphere of radius R, the Ricci scalar is \mathcal{R} = 2/R^2 — smaller spheres are more curved.

The Einstein Tensor

Einstein's great insight was finding the right combination of curvature quantities to equate with the energy-momentum tensor. The Einstein tensor is:

This tensor has a crucial property: it is divergence-free (\nabla^\mu G_{\mu\nu} = 0), which automatically ensures local energy-momentum conservation (\nabla^\mu T_{\mu\nu} = 0). Einstein's field equations are:

where \Lambda is the cosmological constant and T_{\mu\nu} is the stress-energy tensor describing matter and energy content.

The Schwarzschild Solution

Just weeks after Einstein published his field equations in November 1915, Karl Schwarzschild found the first exact solution while serving on the Russian front in World War I. The Schwarzschild metric describes the spacetime geometry outside a spherically symmetric, non-rotating, uncharged mass M:

where the Schwarzschild radius is:

Key features of this metric:

• At r \to \infty, it reduces to the flat Minkowski metric (as it should — far from the mass, spacetime is approximately flat)
• At r = r_s, the event horizon appears: g_{tt} = 0 and g_{rr} \to \infty. This is a coordinate singularity (removable by choosing different coordinates like Kruskal-Szekeres), but it marks the boundary from which nothing can escape
• At r = 0, there is a genuine physical singularity where R^{\mu\nu\rho\sigma}R_{\mu\nu\rho\sigma} \to \infty
• For the Sun: r_s \approx 3 \text{ km} (well inside the Sun's radius, so the Schwarzschild exterior solution applies outside the Sun's surface)

The effective potential for orbital motion in the Schwarzschild geometry is:

The last term is the GR correction absent in Newtonian gravity — it causes Mercury's perihelion precession and allows the existence of an innermost stable circular orbit (ISCO) at r = 6GM/c^2 = 3r_s.

import numpy as np

# ============================================================
# SCHWARZSCHILD ORBITAL EFFECTIVE POTENTIAL
# Compares Newtonian gravity with GR corrections
# ============================================================

# Physical constants
G = 6.674e-11        # gravitational constant (m^3 kg^-1 s^-2)
c = 2.998e8          # speed of light (m/s)
M_sun = 1.989e30     # solar mass (kg)

# Black hole mass (10 solar masses)
M = 10 * M_sun
r_s = 2 * G * M / c**2  # Schwarzschild radius

print("=" * 60)
print("SCHWARZSCHILD EFFECTIVE POTENTIAL")
print("=" * 60)
print(f"\nBlack hole mass: M = {M/M_sun:.0f} M_sun")
print(f"Schwarzschild radius: r_s = {r_s:.2f} m = {r_s/1000:.4f} km")

# Angular momentum (normalized: L per unit mass)
# For circular orbit at r = 10*r_s in GR:
# L^2 = GMr^2 / (r - 3GM/c^2)  for circular orbits
r_orbit = 10 * r_s
L_squared = G * M * r_orbit**2 / (r_orbit - 3*G*M/c**2)
L = np.sqrt(L_squared)

print(f"\nOrbital radius: r = {r_orbit/r_s:.0f} r_s = {r_orbit:.2f} m")
print(f"Specific angular momentum: L = {L:.4e} m^2/s")

# Radial coordinate array (from 3*r_s to 50*r_s)
r = np.linspace(3 * r_s, 50 * r_s, 500)

# Newtonian effective potential (per unit mass, units of c^2)
V_newton = -G * M / r + L**2 / (2 * r**2)

# GR effective potential (Schwarzschild)
V_gr = -G * M / r + L**2 / (2 * r**2) - G * M * L**2 / (c**2 * r**3)

# Normalize to c^2 for cleaner numbers
V_newton_norm = V_newton / c**2
V_gr_norm = V_gr / c**2

# Find ISCO (innermost stable circular orbit)
r_isco = 6 * G * M / c**2  # = 3 * r_s
print(f"\nISCO radius: r_ISCO = 3 r_s = {r_isco:.2f} m")
print(f"ISCO radius: r_ISCO = {r_isco/r_s:.1f} r_s")

# Find potential minimum (stable circular orbit location)
# dV_gr/dr = 0 gives the circular orbit condition
# GM/r^2 - L^2/r^3 + 3GML^2/(c^2 r^4) = 0
# Numerically find the minimum
dV_gr = np.gradient(V_gr_norm, r)
min_idx = np.argmin(np.abs(dV_gr[50:]))  # skip inner region
r_min = r[50 + min_idx]

print(f"Potential minimum at: r = {r_min/r_s:.2f} r_s")

# Key comparison: GR correction magnitude at various radii
print(f"\n--- GR Correction Term vs Newtonian ---")
print(f"{'r/r_s':<10} {'V_Newton/c²':<18} {'V_GR/c²':<18} {'GR correction':<15}")
print("-" * 60)
for r_val in [4*r_s, 6*r_s, 10*r_s, 20*r_s, 50*r_s]:
    vn = (-G*M/r_val + L**2/(2*r_val**2)) / c**2
    vgr = vn - G*M*L**2/(c**2 * r_val**3) / c**2
    correction = -G*M*L**2/(c**2 * r_val**3) / c**2
    print(f"  {r_val/r_s:<8.0f} {vn:<18.8e} {vgr:<18.8e} {correction:<15.8e}")

# Compute orbital energy at ISCO
E_isco = (1 - 2*G*M/(c**2 * r_isco)) / np.sqrt(1 - 3*G*M/(c**2 * r_isco))
print(f"\nBinding energy at ISCO:")
print(f"  E/mc² = {E_isco:.6f}")
print(f"  Binding energy = {(1-E_isco)*100:.2f}% of rest mass")
print(f"  (This is the maximum energy extractable from accretion)")
print(f"  For comparison, nuclear fusion extracts only ~0.7% of rest mass")

The Kerr Metric

Real astrophysical black holes rotate. The Kerr metric, discovered by Roy Kerr in 1963, describes the geometry of a rotating, uncharged black hole. It is arguably the most important exact solution in astrophysics, as all observed black holes are expected to be well-described by the Kerr geometry (the "no-hair theorem" states that a black hole is completely characterized by mass M, angular momentum J, and charge Q).

In Boyer-Lindquist coordinates (t, r, \theta, \phi), the Kerr metric is:

where:

• a = J/(Mc) — the spin parameter (angular momentum per unit mass, divided by c)
• \Sigma = r^2 + a^2\cos^2\theta
• \Delta = r^2 - r_s r + a^2
• r_s = 2GM/c^2 — Schwarzschild radius

When a = 0 (no rotation), the Kerr metric reduces to Schwarzschild. The maximum possible spin is a = r_s/2 = GM/c^2 (an extremal Kerr black hole).

Key Features of Kerr Geometry

Two horizons: The event horizons occur where \Delta = 0:

Ergosphere: Outside the outer horizon, there exists an ergosphere — a region where spacetime is dragged so violently that nothing can remain stationary. Its outer boundary (the static limit) is:

In the ergosphere, the g_{tt} component becomes positive, meaning the time Killing vector becomes spacelike — even light must co-rotate with the black hole. This enables the Penrose process: energy can be extracted from a rotating black hole by exploiting the ergosphere geometry.

Frame dragging: The off-diagonal term g_{t\phi} (the dt \, d\phi cross-term) means that spacetime itself is dragged in the direction of rotation. An observer falling radially inward from infinity will acquire angular velocity — they are swept along by the rotating spacetime.

flowchart TD
    A["Einstein Field Equations
G_μν + Λg_μν = 8πG/c⁴ T_μν"] --> B["Vacuum Solutions
T_μν = 0"]
    A --> C["Matter Solutions
T_μν ≠ 0"]
    
    B --> D["Schwarzschild
Static, spherical
M only"]
    B --> E["Kerr
Stationary, axial
M, J"]
    B --> F["Kerr-Newman
Stationary, axial
M, J, Q"]
    B --> G["Gravitational Waves
Linearized solutions"]
    
    C --> H["Friedmann (FLRW)
Homogeneous cosmology"]
    C --> I["Oppenheimer-Snyder
Dust collapse"]
    C --> J["Reissner-Nordström
Charged, static"]
    
    D --> K["No Hair Theorem
All BH → Kerr-Newman"]
    E --> K
    F --> K

    style A fill:#BF092F,color:#fff
    style D fill:#3B9797,color:#fff
    style E fill:#3B9797,color:#fff
    style F fill:#3B9797,color:#fff
    style H fill:#16476A,color:#fff
    style K fill:#132440,color:#fff
                            

import numpy as np

# ============================================================
# KERR BLACK HOLE: ERGOSPHERE AND HORIZON STRUCTURE
# Visualizing the anatomy of a rotating black hole
# ============================================================

# Physical constants
G = 6.674e-11        # gravitational constant (m^3 kg^-1 s^-2)
c = 2.998e8          # speed of light (m/s)
M_sun = 1.989e30     # solar mass (kg)

# Black hole parameters
M = 10 * M_sun       # 10 solar mass black hole
r_s = 2 * G * M / c**2  # Schwarzschild radius

# Spin parameter: a = J/(Mc), must satisfy a <= r_s/2 = GM/c^2
a_max = G * M / c**2  # maximum spin (extremal Kerr)

print("=" * 60)
print("KERR BLACK HOLE GEOMETRY")
print("=" * 60)
print(f"\nBlack hole mass: M = {M/M_sun:.0f} M_sun")
print(f"Schwarzschild radius: r_s = {r_s:.2f} m")
print(f"Maximum spin parameter: a_max = GM/c² = {a_max:.2f} m")

# Explore different spin values
spin_fractions = [0.0, 0.3, 0.5, 0.7, 0.9, 0.95, 0.998]

print(f"\n{'a/a_max':<10} {'a (m)':<12} {'r+ (m)':<12} {'r- (m)':<12} {'r_ergo(eq)':<14} {'r_ergo(pole)':<14}")
print("-" * 74)

for spin_frac in spin_fractions:
    a = spin_frac * a_max
    
    # Event horizons: r± = r_s/2 ± sqrt((r_s/2)² - a²)
    discriminant = (r_s / 2)**2 - a**2
    if discriminant >= 0:
        r_plus = r_s / 2 + np.sqrt(discriminant)   # outer horizon
        r_minus = r_s / 2 - np.sqrt(discriminant)  # inner (Cauchy) horizon
    else:
        r_plus = np.nan  # naked singularity (unphysical)
        r_minus = np.nan
    
    # Ergosphere outer boundary: r_ergo(theta) = r_s/2 + sqrt((r_s/2)² - a²cos²θ)
    # At equator (theta = pi/2): cos(theta) = 0
    r_ergo_equator = r_s / 2 + np.sqrt((r_s / 2)**2 - a**2 * 0**2)  # = r_s
    # At pole (theta = 0): cos(theta) = 1
    if (r_s / 2)**2 - a**2 >= 0:
        r_ergo_pole = r_s / 2 + np.sqrt((r_s / 2)**2 - a**2 * 1**2)
    else:
        r_ergo_pole = np.nan
    
    print(f"  {spin_frac:<8.3f} {a:<12.2f} {r_plus:<12.2f} {r_minus:<12.2f} {r_ergo_equator:<14.2f} {r_ergo_pole:<14.2f}")

# Detailed analysis for a moderately spinning black hole (a = 0.9 a_max)
print(f"\n{'='*60}")
print(f"DETAILED ANALYSIS: a = 0.9 a_max (rapidly spinning)")
print(f"{'='*60}")
a = 0.9 * a_max
r_plus = r_s / 2 + np.sqrt((r_s/2)**2 - a**2)
r_minus = r_s / 2 - np.sqrt((r_s/2)**2 - a**2)

print(f"\n  Spin parameter: a = {a:.4f} m")
print(f"  Outer horizon: r+ = {r_plus:.4f} m = {r_plus/r_s:.4f} r_s")
print(f"  Inner horizon: r- = {r_minus:.4f} m = {r_minus/r_s:.4f} r_s")

# Ergosphere shape at various latitudes
print(f"\n  Ergosphere boundary r_ergo(θ):")
thetas = np.array([0, 15, 30, 45, 60, 75, 90])
for theta_deg in thetas:
    theta_rad = np.radians(theta_deg)
    r_ergo = r_s/2 + np.sqrt((r_s/2)**2 - a**2 * np.cos(theta_rad)**2)
    thickness = r_ergo - r_plus
    print(f"    θ = {theta_deg:>3}°: r_ergo = {r_ergo/r_s:.4f} r_s, "
          f"ergosphere thickness = {thickness/r_s:.4f} r_s")

# ISCO for Kerr (prograde and retrograde)
# Prograde ISCO shrinks as spin increases; retrograde grows
# For a/M -> 1: r_ISCO(prograde) -> GM/c^2, r_ISCO(retrograde) -> 9GM/c^2
r_isco_schwarz = 6 * G * M / c**2  # = 3 r_s
print(f"\n  ISCO comparisons:")
print(f"    Schwarzschild (a=0): r_ISCO = 6 GM/c² = 3 r_s = {r_isco_schwarz:.2f} m")
print(f"    Extremal prograde:   r_ISCO → GM/c² = 0.5 r_s (radiative efficiency ~42%!)")
print(f"    Extremal retrograde: r_ISCO → 9 GM/c² = 4.5 r_s")

# Energy extraction via Penrose process
# Maximum energy extractable = (1 - 1/sqrt(2)) * Mc^2 ≈ 29% for extremal Kerr
E_max_fraction = 1 - 1/np.sqrt(2)
print(f"\n  Penrose process (extremal Kerr):")
print(f"    Max extractable energy: {E_max_fraction*100:.1f}% of Mc²")
print(f"    For M = 10 M_sun: E_max = {E_max_fraction * M * c**2:.3e} J")

Practice Exercises

Test your understanding of the mathematical formalism:

Exercise 1: Metric & Christoffel Symbols

Consider the 2D metric ds^2 = dr^2 + r^2 d\phi^2 (flat plane in polar coordinates).

1. Write down all components of g_{\mu\nu} and g^{\mu\nu}.
2. Compute all non-zero Christoffel symbols \Gamma^a_{bc}.
3. Write down the geodesic equations. Show they give straight lines.
4. What happens if you change to ds^2 = dr^2 + \sinh^2(r) \, d\phi^2? Is this a flat space? (Hint: compute the Riemann tensor.)

Exercise 2: Schwarzschild Geometry

For a 5 M_\odot black hole:

1. Calculate the Schwarzschild radius in km.
2. What is the proper time for an observer to fall from the event horizon to the singularity? (Answer: \tau = \pi GM/c^3)
3. Calculate the tidal acceleration (stretching force) at r = 2r_s for a 2-meter tall astronaut. Is this survivable?
4. What is the radius of the photon sphere (unstable circular light orbit)?

Exercise 3: Kerr Black Hole

A stellar black hole of mass M = 10 M_\odot has spin parameter a = 0.95 \, GM/c^2.

1. Calculate the radii of both event horizons (r_+ and r_-).
2. Find the ergosphere radius at the equator and at the pole.
3. What is the maximum energy extractable by the Penrose process from this black hole?
4. If the ISCO for prograde orbits is at r \approx 2.32 \, GM/c^2, what is the radiative efficiency of an accretion disk? (Compare with Schwarzschild's ~5.7%.)

Exercise 4: Tensor Manipulation

Given a rank-2 tensor T^{\mu\nu} and the Minkowski metric \eta_{\mu\nu} = \text{diag}(-1,1,1,1):

1. Show that T^\mu{}_\mu = \eta_{\mu\nu}T^{\mu\nu} (the trace).
2. If T^{\mu\nu} = \text{diag}(\rho, p, p, p) (perfect fluid), compute T = T^\mu{}_\mu.
3. Compute T_{\mu\nu}T^{\mu\nu} for the same perfect fluid tensor.

Conclusion & Next Steps

We have now developed the full mathematical machinery of General Relativity. The journey from the metric tensor through Christoffel symbols, the Riemann curvature tensor, and finally to Einstein's field equations reveals an extraordinary mathematical structure: geometry tells matter how to move; matter tells geometry how to curve.

The key mathematical objects form a clear hierarchy:

• Metric tensor g_{\mu\nu} — defines spacetime geometry (distances, angles, causal structure)
• Christoffel symbols \Gamma^\lambda_{\mu\nu} — encode the "gravitational force" (connection, parallel transport)
• Riemann tensor R^\rho{}_{\sigma\mu\nu} — measures intrinsic curvature (tidal forces)
• Ricci tensor & scalar — contracted curvature, sourced by matter
• Einstein tensor G_{\mu\nu} — the divergence-free combination equated with stress-energy

The two exact solutions we studied — Schwarzschild (1916) and Kerr (1963) — remain the most important in astrophysics. Every black hole observed by LIGO, the EHT, or X-ray telescopes is ultimately described by the Kerr metric with just two parameters: mass and spin. The mathematical elegance of this "no-hair" simplicity is one of the deepest results in theoretical physics.

In the next part, we explore the frontiers where this mathematical framework begins to break down — quantum gravity, singularity theorems, wormholes, and the quest to unify GR with quantum mechanics.