Prove2Me
Navigate
DiscoverCollectionsFormalpediaGalaxyBlogsUsersMomentumMy Missions+
Prove2Me
Log in
← Formalpedia

No-hair: spacetimes, causal structure, stationary black holes

Definition
NoHair_spacetime

by Lucas · Oct 4, 2026 · Mathlib 0df444a (Lean v4.33.1)

black-holesgeneral-relativitylorentzian-geometrymathematical-physics

Spacetimes, causal structure and stationary black holes.

A spacetime is a C∞C^\inftyC∞ manifold MMM modelled on R4\mathbb R^4R4, carrying a field ggg of bilinear forms on the tangent spaces (the metric), a field FFF of bilinear forms (the electromagnetic field) and a vector field TTT. Components of a bilinear field BBB along a map ψ:R4→M\psi:\mathbb R^4\to Mψ:R4→M are (ψ∗B)μν(y)=Bψ(y)(dψy∂μ,dψy∂ν)(\psi^*B)_{\mu\nu}(y)=B_{\psi(y)}(d\psi_y\partial_\mu,d\psi_y\partial_\nu)(ψ∗B)μν​(y)=Bψ(y)​(dψy​∂μ​,dψy​∂ν​); along the inverse of a chart these are the usual coordinate components. The notions defined are:

  1. ggg is Lorentzian if at each point some basis of the tangent space has ggg-Gram matrix diag(−1,1,1,1)\mathrm{diag}(-1,1,1,1)diag(−1,1,1,1); BBB is smooth (resp. analytic) if its components are C∞C^\inftyC∞ (resp. CωC^\omegaCω) in every chart; FFF is a 2-form if antisymmetric; TTT is a time orientation if continuous and everywhere timelike.
  2. A tangent vector vvv is future timelike if g(v,v)<0g(v,v)<0g(v,v)<0, g(v,T)<0g(v,T)<0g(v,T)<0, and future causal if v≠0v\neq0v=0, g(v,v)≤0g(v,v)\le0g(v,v)≤0, g(v,T)<0g(v,T)<0g(v,T)<0. Chronological and causal futures/pasts IU±(p)I^\pm_U(p)IU±​(p), JU±(p)J^\pm_U(p)JU±​(p) inside U⊆MU\subseteq MU⊆M are defined by C1C^1C1 curves on [0,1][0,1][0,1] lying in UUU with future timelike (resp. causal) velocity; JU±(p)J^\pm_U(p)JU±​(p) contain ppp. For a set SSS, I±(S)=⋃p∈SIM±(p)I^\pm(S)=\bigcup_{p\in S}I^\pm_M(p)I±(S)=⋃p∈S​IM±​(p).
  3. UUU is globally hyperbolic if it contains no closed future causal curve and every causal diamond JU+(p)∩JU−(q)J^+_U(p)\cap J^-_U(q)JU+​(p)∩JU−​(q), p,q∈Up,q\in Up,q∈U, is compact.
  4. (g,F)(g,F)(g,F) is electrovacuum if the Einstein–Maxwell system holds at every point of every chart.
  5. A symmetry flow is a smooth R\mathbb RR-action φt\varphi_tφt​ on MMM preserving ggg and FFF; its generator is X(x)=ddt∣t=0φt(x)X(x)=\frac{d}{dt}\big|_{t=0}\varphi_t(x)X(x)=dtd​​t=0​φt​(x).
  6. A stationary asymptotically flat end consists of R>0R>0R>0, C∈RC\in\mathbb RC∈R and a diffeomorphism ψ\psiψ from {(t,x⃗):∣x⃗∣>R}\{(t,\vec x):|\vec x|>R\}{(t,x):∣x∣>R} onto an open set Mext⊆MM_{\rm ext}\subseteq MMext​⊆M with ψ(y+s∂t)=φs(ψ(y))\psi(y+s\partial_t)=\varphi_s(\psi(y))ψ(y+s∂t​)=φs​(ψ(y)), ∣gμν−ημν∣≤C/∣x⃗∣|g_{\mu\nu}-\eta_{\mu\nu}|\le C/|\vec x|∣gμν​−ημν​∣≤C/∣x∣, ∣∂σgμν∣≤C/∣x⃗∣2|\partial_\sigma g_{\mu\nu}|\le C/|\vec x|^2∣∂σ​gμν​∣≤C/∣x∣2 and ∣Fμν∣≤C/∣x⃗∣2|F_{\mu\nu}|\le C/|\vec x|^2∣Fμν​∣≤C/∣x∣2 in these coordinates.
  7. The domain of outer communications is ⟨⟨Mext⟩⟩=I+(Mext)∩I−(Mext)\langle\langle M_{\rm ext}\rangle\rangle=I^+(M_{\rm ext})\cap I^-(M_{\rm ext})⟨⟨Mext​⟩⟩=I+(Mext​)∩I−(Mext​), the black-hole region is M∖I−(Mext)M\setminus I^-(M_{\rm ext})M∖I−(Mext​) and the future event horizon is ∂I−(Mext)\partial I^-(M_{\rm ext})∂I−(Mext​).
  8. Stationary black hole (the standing hypotheses): ggg Lorentzian and smooth, FFF a smooth 2-form, TTT a time orientation, (g,F)(g,F)(g,F) electrovacuum, φ\varphiφ a symmetry flow with an asymptotically flat end as in 6, ⟨⟨Mext⟩⟩\langle\langle M_{\rm ext}\rangle\rangle⟨⟨Mext​⟩⟩ globally hyperbolic, nonempty black-hole region and connected event horizon.
  9. A set D⊆MD\subseteq MD⊆M is a Kerr–Newman domain of outer communications with parameters (m,a,e)(m,a,e)(m,a,e) if there are β∈R\beta\in\mathbb Rβ∈R and a diffeomorphism Ψ\PsiΨ from {r>r+}\{r>r_+\}{r>r+​} onto DDD with Ψ∗g=gKN\Psi^*g=g^{KN}Ψ∗g=gKN and Ψ∗F=cos⁡β FKN+sin⁡β ⋆FKN\Psi^*F=\cos\beta\,F^{KN}+\sin\beta\,\star F^{KN}Ψ∗F=cosβFKN+sinβ⋆FKN (constant duality rotation).
  10. Auxiliary hypotheses used by the milestones: static on DDD (the stationary generator KKK satisfies K∧dK=0K\wedge dK=0K∧dK=0 there), axial symmetry (a 2π2\pi2π-periodic symmetry flow commuting with φ\varphiφ, with a fixed point and not the identity), and non-degenerate horizon (a symmetry flow whose generator XXX is null on H∩D‾H\cap\overline DH∩D, nonzero somewhere there, and satisfies ∂μ(g(X,X))=−2κXμ\partial_\mu(g(X,X))=-2\kappa X_\mu∂μ​(g(X,X))=−2κXμ​ there with a constant κ≠0\kappa\neq0κ=0).

These notions make the informal phrases of the no-hair theorem ("stationary", "black hole", "event horizon", "described by the Kerr–Newman metric") precise.

Formalization Note Tangent spaces are Mathlib's TangentSpace 𝓘(ℝ, ℝ⁴) x (identified with R4\mathbb R^4R4); derivatives of maps are mfderiv, which returns 000 at points of non-differentiability, so all statements that use it also impose smoothness.

Definition code
import Mathlib
import Definitions.Def_NoHair_coordinates

/-!
# No-hair mission: spacetimes, causal structure, stationary black holes

A spacetime is a smooth manifold `M` modelled on `ℝ⁴`. Tangent spaces `T_xM` are identified
with `ℝ⁴` through Mathlib's convention (components in the preferred chart at `x`); every notion
below is chart-independent because it is expressed through `mfderiv` or through the components
in *every* chart of the atlas.
-/

noncomputable section

open scoped ContDiff Manifold BigOperators
open Set

namespace NoHair

variable {M : Type*} [TopologicalSpace M] [ChartedSpace E4 M] [IsManifold 𝓘(ℝ, E4) ∞ M]

/-- A field of bilinear forms on the tangent spaces of `M` (used for the metric `g` and for the
electromagnetic field `F`). -/
abbrev BilinField (M : Type*) [TopologicalSpace M] [ChartedSpace E4 M] :=
  ∀ x : M, TangentSpace 𝓘(ℝ, E4) x →L[ℝ] TangentSpace 𝓘(ℝ, E4) x →L[ℝ] ℝ

/-- A vector field on `M`. -/
abbrev VecField (M : Type*) [TopologicalSpace M] [ChartedSpace E4 M] :=
  ∀ x : M, TangentSpace 𝓘(ℝ, E4) x

/-- Components of the pull-back of a bilinear field `B` along a map `ψ : ℝ⁴ → M` at `y`:
`(ψ^* B)_{μν}(y) = B_{ψ y}(dψ_y ∂_μ, dψ_y ∂_ν)`. For `ψ = e.symm` with `e` a chart this is the
coordinate matrix of `B` in that chart. -/
def pullbackCoeffs (B : BilinField M) (ψ : E4 → M) (y : E4) : Matrix (Fin 4) (Fin 4) ℝ :=
  fun μ ν => B (ψ y) (mfderiv 𝓘(ℝ, E4) 𝓘(ℝ, E4) ψ y (coordBasis μ))
    (mfderiv 𝓘(ℝ, E4) 𝓘(ℝ, E4) ψ y (coordBasis ν))

/-- Components `X_μ = g(X, ∂_μ)` of the covector `g(X, ·)` along a map `ψ : ℝ⁴ → M`. -/
def loweredCoeffs (g : BilinField M) (X : VecField M) (ψ : E4 → M) (y : E4) : Fin 4 → ℝ :=
  fun μ => g (ψ y) (X (ψ y)) (mfderiv 𝓘(ℝ, E4) 𝓘(ℝ, E4) ψ y (coordBasis μ))

/-- `B` has smooth components in every chart of the atlas. -/
def IsSmoothBilinField (B : BilinField M) : Prop :=
  ∀ e ∈ atlas E4 M, ∀ μ ν : Fin 4,
    ContDiffOn ℝ ∞ (fun y => pullbackCoeffs B e.symm y μ ν) e.target

/-- `B` has real-analytic components in every chart of the atlas. -/
def IsAnalyticBilinField (B : BilinField M) : Prop :=
  ∀ e ∈ atlas E4 M, ∀ μ ν : Fin 4,
    ContDiffOn ℝ ω (fun y => pullbackCoeffs B e.symm y μ ν) e.target

/-- `g` is a Lorentzian metric of signature `(-,+,+,+)`: at every point there is a basis of the
tangent space in which `g` has components `diag(-1,1,1,1)`. -/
def IsLorentzian (g : BilinField M) : Prop :=
  ∀ x : M, ∃ b : Module.Basis (Fin 4) ℝ (TangentSpace 𝓘(ℝ, E4) x),
    ∀ μ ν, g x (b μ) (b ν) = minkowski μ ν

/-- `F` is a 2-form (antisymmetric at every point). -/
def IsTwoForm (F : BilinField M) : Prop := ∀ x v w, F x v w = - F x w v

/-- `T` is a time orientation for `g`: a continuous, everywhere timelike vector field. -/
def IsTimeOrientation (g : BilinField M) (T : VecField M) : Prop :=
  Continuous (fun x => (⟨x, T x⟩ : TangentBundle 𝓘(ℝ, E4) M)) ∧ ∀ x, g x (T x) (T x) < 0

/-- A tangent vector is future-directed timelike. -/
def IsFutureTimelike (g : BilinField M) (T : VecField M) (x : M)
    (v : TangentSpace 𝓘(ℝ, E4) x) : Prop :=
  g x v v < 0 ∧ g x v (T x) < 0

/-- A tangent vector is future-directed causal (nonzero, timelike or null, future-pointing). -/
def IsFutureCausal (g : BilinField M) (T : VecField M) (x : M)
    (v : TangentSpace 𝓘(ℝ, E4) x) : Prop :=
  v ≠ 0 ∧ g x v v ≤ 0 ∧ g x v (T x) < 0

/-- Velocity `γ'(t)` of a curve `γ : ℝ → M`. -/
def velocity (γ : ℝ → M) (t : ℝ) : TangentSpace 𝓘(ℝ, E4) (γ t) :=
  mfderiv 𝓘(ℝ, ℝ) 𝓘(ℝ, E4) γ t (1 : ℝ)

/-- `γ` restricted to `[0,1]` is a `C¹` future-directed timelike curve lying in `U`. -/
def IsFutureTimelikeCurveIn (g : BilinField M) (T : VecField M) (U : Set M) (γ : ℝ → M) :
    Prop :=
  ContMDiff 𝓘(ℝ, ℝ) 𝓘(ℝ, E4) 1 γ ∧
    ∀ t ∈ Icc (0 : ℝ) 1, γ t ∈ U ∧ IsFutureTimelike g T (γ t) (velocity γ t)

/-- `γ` restricted to `[0,1]` is a `C¹` future-directed causal curve lying in `U`. -/
def IsFutureCausalCurveIn (g : BilinField M) (T : VecField M) (U : Set M) (γ : ℝ → M) : Prop :=
  ContMDiff 𝓘(ℝ, ℝ) 𝓘(ℝ, E4) 1 γ ∧
    ∀ t ∈ Icc (0 : ℝ) 1, γ t ∈ U ∧ IsFutureCausal g T (γ t) (velocity γ t)

/-- Chronological future `I⁺_U(p)` of `p` inside `U`. -/
def chronoFuture (g : BilinField M) (T : VecField M) (U : Set M) (p : M) : Set M :=
  {q | ∃ γ : ℝ → M, IsFutureTimelikeCurveIn g T U γ ∧ γ 0 = p ∧ γ 1 = q}

/-- Chronological past `I⁻_U(p)` of `p` inside `U`. -/
def chronoPast (g : BilinField M) (T : VecField M) (U : Set M) (p : M) : Set M :=
  {q | p ∈ chronoFuture g T U q}

/-- Causal future `J⁺_U(p)` of `p` inside `U` (contains `p`). -/
def causalFuture (g : BilinField M) (T : VecField M) (U : Set M) (p : M) : Set M :=
  {p} ∪ {q | ∃ γ : ℝ → M, IsFutureCausalCurveIn g T U γ ∧ γ 0 = p ∧ γ 1 = q}

/-- Causal past `J⁻_U(p)` of `p` inside `U` (contains `p`). -/
def causalPast (g : BilinField M) (T : VecField M) (U : Set M) (p : M) : Set M :=
  {q | p ∈ causalFuture g T U q}

/-- Chronological future of a set, `I⁺(S) = ⋃_{p ∈ S} I⁺(p)` (curves in all of `M`). -/
def chronoFutureSet (g : BilinField M) (T : VecField M) (S : Set M) : Set M :=
  ⋃ p ∈ S, chronoFuture g T univ p

/-- Chronological past of a set, `I⁻(S) = ⋃_{p ∈ S} I⁻(p)` (curves in all of `M`). -/
def chronoPastSet (g : BilinField M) (T : VecField M) (S : Set M) : Set M :=
  ⋃ p ∈ S, chronoPast g T univ p

/-- The open set `U` with the restricted metric is globally hyperbolic: it is causal (no closed
future-directed causal curve in `U`) and all causal diamonds `J⁺_U(p) ∩ J⁻_U(q)` are compact. -/
def IsGloballyHyperbolicOn (g : BilinField M) (T : VecField M) (U : Set M) : Prop :=
  (∀ p ∈ U, ¬ ∃ γ : ℝ → M, IsFutureCausalCurveIn g T U γ ∧ γ 0 = p ∧ γ 1 = p) ∧
  ∀ p ∈ U, ∀ q ∈ U, IsCompact (causalFuture g T U p ∩ causalPast g T U q)

/-- `(g, F)` solves the Einstein–Maxwell equations (electrovacuum) in every chart of the atlas. -/
def IsElectrovac (g F : BilinField M) : Prop :=
  ∀ e ∈ atlas E4 M, ∀ y ∈ e.target,
    EinsteinMaxwellAt (pullbackCoeffs g e.symm) (pullbackCoeffs F e.symm) y

/-- `φ` is a smooth one-parameter group of diffeomorphisms (an `ℝ`-action) of `M` preserving
the metric `g` and the field `F`. -/
def IsSymmetryFlow (g F : BilinField M) (φ : ℝ → M → M) : Prop :=
  (∀ x, φ 0 x = x) ∧ (∀ s t x, φ (s + t) x = φ s (φ t x)) ∧
  ContMDiff (𝓘(ℝ, ℝ).prod 𝓘(ℝ, E4)) 𝓘(ℝ, E4) ∞ (fun p : ℝ × M => φ p.1 p.2) ∧
  (∀ t x v w, g (φ t x) (mfderiv 𝓘(ℝ, E4) 𝓘(ℝ, E4) (φ t) x v)
      (mfderiv 𝓘(ℝ, E4) 𝓘(ℝ, E4) (φ t) x w) = g x v w) ∧
  (∀ t x v w, F (φ t x) (mfderiv 𝓘(ℝ, E4) 𝓘(ℝ, E4) (φ t) x v)
      (mfderiv 𝓘(ℝ, E4) 𝓘(ℝ, E4) (φ t) x w) = F x v w)

/-- The infinitesimal generator `X(x) = d/dt|_{t=0} φ_t(x)` of a flow (a Killing field when `φ`
is a symmetry flow). -/
def flowGenerator (φ : ℝ → M → M) : VecField M :=
  fun x => mfderiv 𝓘(ℝ, ℝ) 𝓘(ℝ, E4) (fun t => φ t x) 0 (1 : ℝ)

/-- The asymptotic region `{(t, x⃗) : |x⃗| > R}` of coordinate space. -/
def exteriorRegion (R : ℝ) : Set E4 := {y | R < spatialRadius y}

/-- `ψ` maps the open set `S ⊆ ℝ⁴` diffeomorphically onto `D ⊆ M`. -/
def IsDiffeoOnto (ψ : E4 → M) (S : Set E4) (D : Set M) : Prop :=
  ∃ χ : M → E4, MapsTo ψ S D ∧ MapsTo χ D S ∧ LeftInvOn χ ψ S ∧ RightInvOn χ ψ D ∧
    ContMDiffOn 𝓘(ℝ, E4) 𝓘(ℝ, E4) ∞ ψ S ∧ ContMDiffOn 𝓘(ℝ, E4) 𝓘(ℝ, E4) ∞ χ D

/-- A stationary asymptotically flat end for `(g, F)` adapted to the stationary flow `φ`:
coordinates `ψ` identify `{|x⃗| > R}` (with `R > 0`) with an open subset of `M` such that
1. the stationary flow is time translation: `ψ(y + s ∂_t) = φ_s(ψ y)`;
2. `|g_{μν} - η_{μν}| ≤ C/|x⃗|` and `|∂_σ g_{μν}| ≤ C/|x⃗|²`;
3. `|F_{μν}| ≤ C/|x⃗|²`. -/
structure AsymptoticallyFlatEnd (g F : BilinField M) (φ : ℝ → M → M) where
  /-- Radius of the asymptotic region. -/
  R : ℝ
  /-- Asymptotic coordinates. -/
  ψ : E4 → M
  /-- Decay constant. -/
  C : ℝ
  R_pos : 0 < R
  diffeo : IsDiffeoOnto ψ (exteriorRegion R) (ψ '' exteriorRegion R)
  isOpen_image : IsOpen (ψ '' exteriorRegion R)
  equivariant : ∀ s : ℝ, ∀ y ∈ exteriorRegion R, ψ (y + s • coordBasis 0) = φ s (ψ y)
  metric_decay : ∀ y ∈ exteriorRegion R, ∀ μ ν,
    |pullbackCoeffs g ψ y μ ν - minkowski μ ν| ≤ C / spatialRadius y
  metric_deriv_decay : ∀ y ∈ exteriorRegion R, ∀ μ ν σ,
    |coordPartial (fun z => pullbackCoeffs g ψ z μ ν) σ y| ≤ C / spatialRadius y ^ 2
  field_decay : ∀ y ∈ exteriorRegion R, ∀ μ ν,
    |pullbackCoeffs F ψ y μ ν| ≤ C / spatialRadius y ^ 2

/-- The asymptotic region `M_ext = ψ({|x⃗| > R})` of an end. -/
def AsymptoticallyFlatEnd.region {g F : BilinField M} {φ : ℝ → M → M}
    (E : AsymptoticallyFlatEnd g F φ) : Set M :=
  E.ψ '' exteriorRegion E.R

/-- Domain of outer communications `⟨⟨M_ext⟩⟩ = I⁺(M_ext) ∩ I⁻(M_ext)`. -/
def domainOfOuterCommunications (g : BilinField M) (T : VecField M) (Mext : Set M) : Set M :=
  chronoFutureSet g T Mext ∩ chronoPastSet g T Mext

/-- Black-hole region `M \ I⁻(M_ext)`. -/
def blackHoleRegion (g : BilinField M) (T : VecField M) (Mext : Set M) : Set M :=
  (chronoPastSet g T Mext)ᶜ

/-- Future event horizon `∂ I⁻(M_ext)`. -/
def futureEventHorizon (g : BilinField M) (T : VecField M) (Mext : Set M) : Set M :=
  frontier (chronoPastSet g T Mext)

/-- The standing hypotheses of the no-hair (black-hole uniqueness) problem:
`(M, g, F, T)` is a smooth, time-oriented, 4-dimensional Lorentzian electrovacuum spacetime,
stationary under the symmetry flow `φ`, with a stationary asymptotically flat end `E`, whose
domain of outer communications is globally hyperbolic, which contains a black hole
(nonempty black-hole region) with a connected future event horizon. -/
structure IsStationaryBlackHole (g F : BilinField M) (T : VecField M) (φ : ℝ → M → M)
    (E : AsymptoticallyFlatEnd g F φ) : Prop where
  lorentzian : IsLorentzian g
  smooth_metric : IsSmoothBilinField g
  smooth_field : IsSmoothBilinField F
  twoForm : IsTwoForm F
  timeOrientation : IsTimeOrientation g T
  electrovac : IsElectrovac g F
  stationary : IsSymmetryFlow g F φ
  globallyHyperbolic :
    IsGloballyHyperbolicOn g T (domainOfOuterCommunications g T E.region)
  blackHole_nonempty : (blackHoleRegion g T E.region).Nonempty
  horizon_connected : IsConnected (futureEventHorizon g T E.region)

/-- The open set `D ⊆ M`, with the restrictions of `(g, F)`, is isometric to the Kerr–Newman
domain of outer communications with parameters `(m, a, e)`, the electromagnetic field
matching the Kerr–Newman field up to a constant duality rotation by angle `β`:
there is a diffeomorphism `Ψ` from `{r > r₊}` (Kerr–Schild coordinates) onto `D` with
`Ψ^* g = g_KN` and `Ψ^* F = cos β · F_KN + sin β · ⋆F_KN`. -/
def IsKerrNewmanDOC (g F : BilinField M) (D : Set M) (m a e : ℝ) : Prop :=
  ∃ β : ℝ, ∃ Ψ : E4 → M, IsDiffeoOnto Ψ (kerrNewmanDOC m a e) D ∧
    ∀ y ∈ kerrNewmanDOC m a e,
      pullbackCoeffs g Ψ y = kerrNewmanMetric m a e y ∧
      pullbackCoeffs F Ψ y = Real.cos β • kerrNewmanField a e y
        + Real.sin β • hodgeDual (kerrNewmanMetric m a e y) (kerrNewmanField a e y)

/-- The stationary flow is *static* on `D`: its generator `K` is hypersurface-orthogonal there,
`K ∧ dK = 0`, i.e. `K_[λ ∂_μ K_ν] = 0` in every chart. -/
def IsStaticOn (g : BilinField M) (φ : ℝ → M → M) (D : Set M) : Prop :=
  ∀ e ∈ atlas E4 M, ∀ y ∈ e.target, e.symm y ∈ D → ∀ l μ ν : Fin 4,
    let k := fun z => loweredCoeffs g (flowGenerator φ) e.symm z
    k y l * (coordPartial (fun z => k z ν) μ y - coordPartial (fun z => k z μ) ν y)
    + k y μ * (coordPartial (fun z => k z l) ν y - coordPartial (fun z => k z ν) l y)
    + k y ν * (coordPartial (fun z => k z μ) l y - coordPartial (fun z => k z l) μ y) = 0

/-- `ρ` is an axial symmetry commuting with the stationary flow `φ`: a `2π`-periodic symmetry
flow of `(g, F)` commuting with `φ`, with a nonempty fixed-point set (the axis) and not the
identity. -/
def IsAxialSymmetry (g F : BilinField M) (φ ρ : ℝ → M → M) : Prop :=
  IsSymmetryFlow g F ρ ∧ (∀ x, ρ (2 * Real.pi) x = x) ∧
  (∀ s t x, ρ s (φ t x) = φ t (ρ s x)) ∧
  (∃ x, ∀ s, ρ s x = x) ∧ (∃ s x, ρ s x ≠ x)

/-- The event horizon `H` is a non-degenerate Killing horizon: there is a symmetry flow `χ` of
`(g, F)` whose generator `X` is null on `H ∩ cl(D)`, not identically zero there, and satisfies
`∂_μ (g(X,X)) = -2κ X_μ` on `H ∩ cl(D)` for a constant surface gravity `κ ≠ 0`. -/
def IsNondegenerateHorizon (g F : BilinField M) (H D : Set M) : Prop :=
  ∃ χ : ℝ → M → M, ∃ κ : ℝ, IsSymmetryFlow g F χ ∧ κ ≠ 0 ∧
    (∃ x ∈ H ∩ closure D, flowGenerator χ x ≠ 0) ∧
    ∀ x ∈ H ∩ closure D,
      g x (flowGenerator χ x) (flowGenerator χ x) = 0 ∧
      ∀ e ∈ atlas E4 M, ∀ y ∈ e.target, e.symm y = x → ∀ μ : Fin 4,
        coordPartial (fun z => g (e.symm z) (flowGenerator χ (e.symm z))
          (flowGenerator χ (e.symm z))) μ y
          = -2 * κ * loweredCoeffs g (flowGenerator χ) e.symm y μ

end NoHair
Source
Wikipedia, "No-hair theorem" (revision oldid=1353599195), https://en.wikipedia.org/w/index.php?title=No-hair_theorem&oldid=1353599195
Read-back

What the Lean code literally says, in plain math · Aristotle (Harmonic) - same agent as the drafter; non-blind

Non-blind read-back. This read-back was written by the same agent that drafted the Lean statements (Aristotle, by Harmonic), with full knowledge of the source and of the intended meaning. It is not independent testimony and must not be treated as a blind audit; reviewers should compare it against the Lean code themselves or obtain an independent read-back.

Throughout, MMM is a topological space with a charted-space structure modelled on R4\mathbb R^4R4 that is a C∞C^\inftyC∞ manifold in Mathlib's sense; no Hausdorff, second-countability, connectedness or orientability assumption is made. A "bilinear field" assigns to each x∈Mx\in Mx∈M a continuous bilinear form on TxM≅R4T_xM\cong\mathbb R^4Tx​M≅R4, a "vector field" assigns a vector of TxMT_xMTx​M (no regularity unless stated). dψyd\psi_ydψy​ denotes Mathlib's manifold derivative (zero where ψ\psiψ is not differentiable).

  1. Pull-back components (ψ∗B)μν(y)=Bψy(dψyeμ,dψyeν)(\psi^*B)_{\mu\nu}(y)=B_{\psi y}(d\psi_y e_\mu,d\psi_y e_\nu)(ψ∗B)μν​(y)=Bψy​(dψy​eμ​,dψy​eν​); lowered components Xμ(y)=gψy(X(ψy),dψyeμ)X_\mu(y)=g_{\psi y}(X(\psi y),d\psi_ye_\mu)Xμ​(y)=gψy​(X(ψy),dψy​eμ​).
  2. Smooth / analytic bilinear field: for every chart eee in the atlas and all μ,ν\mu,\nuμ,ν, y↦(e−1)∗Bμν(y)y\mapsto(e^{-1})^*B_{\mu\nu}(y)y↦(e−1)∗Bμν​(y) is C∞C^\inftyC∞ / CωC^\omegaCω on the chart's target.
  3. Lorentzian: for every xxx there is a basis bbb of TxMT_xMTx​M with gx(bμ,bν)=ημνg_x(b_\mu,b_\nu)=\eta_{\mu\nu}gx​(bμ​,bν​)=ημν​, η=diag(−1,1,1,1)\eta=\mathrm{diag}(-1,1,1,1)η=diag(−1,1,1,1). 2-form: Fx(v,w)=−Fx(w,v)F_x(v,w)=-F_x(w,v)Fx​(v,w)=−Fx​(w,v) for all x,v,wx,v,wx,v,w. Time orientation: x↦(x,T(x))x\mapsto(x,T(x))x↦(x,T(x)) is continuous into the tangent bundle and gx(T,T)<0g_x(T,T)<0gx​(T,T)<0 everywhere.
  4. Future timelike: g(v,v)<0g(v,v)<0g(v,v)<0 and g(v,T)<0g(v,T)<0g(v,T)<0. Future causal: v≠0v\ne0v=0, g(v,v)≤0g(v,v)\le0g(v,v)≤0, g(v,T)<0g(v,T)<0g(v,T)<0. A future timelike (causal) curve in UUU is a C1C^1C1 map γ:R→M\gamma:\mathbb R\to Mγ:R→M such that for every t∈[0,1]t\in[0,1]t∈[0,1], γ(t)∈U\gamma(t)\in Uγ(t)∈U and γ′(t)\gamma'(t)γ′(t) is future timelike (causal). IU+(p)I^+_U(p)IU+​(p) = endpoints γ(1)\gamma(1)γ(1) of such timelike curves with γ(0)=p\gamma(0)=pγ(0)=p; IU−(p)={q:p∈IU+(q)}I^-_U(p)=\{q:p\in I^+_U(q)\}IU−​(p)={q:p∈IU+​(q)}; JU+(p)={p}∪{γ(1)}J^+_U(p)=\{p\}\cup\{\gamma(1)\}JU+​(p)={p}∪{γ(1)} over causal curves from ppp; JU−J^-_UJU−​ analogously. For a set SSS, I±(S)I^\pm(S)I±(S) is the union over p∈Sp\in Sp∈S of IM±(p)I^\pm_M(p)IM±​(p) (curves anywhere in MMM).
  5. Globally hyperbolic on UUU: no p∈Up\in Up∈U admits a future causal curve in UUU from ppp to ppp, and for all p,q∈Up,q\in Up,q∈U the set JU+(p)∩JU−(q)J^+_U(p)\cap J^-_U(q)JU+​(p)∩JU−​(q) is compact.
  6. Electrovacuum: for every atlas chart eee and every yyy in its target, the coordinate Einstein–Maxwell predicate holds at yyy for the chart components of ggg and FFF.
  7. Symmetry flow φ:R×M→M\varphi:\mathbb R\times M\to Mφ:R×M→M: φ0=id\varphi_0=\mathrm{id}φ0​=id, φs+t=φs∘φt\varphi_{s+t}=\varphi_s\circ\varphi_tφs+t​=φs​∘φt​, jointly C∞C^\inftyC∞, and ggg and FFF are invariant: gφtx(dφtv,dφtw)=gx(v,w)g_{\varphi_tx}(d\varphi_tv,d\varphi_tw)=g_x(v,w)gφt​x​(dφt​v,dφt​w)=gx​(v,w) and likewise for FFF. Generator: X(x)=X(x)=X(x)= derivative of t↦φt(x)t\mapsto\varphi_t(x)t↦φt​(x) at t=0t=0t=0 applied to 111.
  8. Asymptotically flat end for (g,F,φ)(g,F,\varphi)(g,F,φ): data R,ψ,CR,\psi,CR,ψ,C with R>0R>0R>0; there is χ:M→R4\chi:M\to\mathbb R^4χ:M→R4 such that ψ\psiψ maps ΩR={y:(y1)2+(y2)2+(y3)2>R}\Omega_R=\{y:\sqrt{(y^1)^2+(y^2)^2+(y^3)^2}>R\}ΩR​={y:(y1)2+(y2)2+(y3)2​>R} into ψ(ΩR)\psi(\Omega_R)ψ(ΩR​), χ\chiχ maps ψ(ΩR)\psi(\Omega_R)ψ(ΩR​) into ΩR\Omega_RΩR​, χ∘ψ=id\chi\circ\psi=\mathrm{id}χ∘ψ=id on ΩR\Omega_RΩR​, ψ∘χ=id\psi\circ\chi=\mathrm{id}ψ∘χ=id on ψ(ΩR)\psi(\Omega_R)ψ(ΩR​), both are C∞C^\inftyC∞ on these sets; ψ(ΩR)\psi(\Omega_R)ψ(ΩR​) is open; for all sss and y∈ΩRy\in\Omega_Ry∈ΩR​, ψ(y+se0)=φs(ψy)\psi(y+se_0)=\varphi_s(\psi y)ψ(y+se0​)=φs​(ψy); and for y∈ΩRy\in\Omega_Ry∈ΩR​, writing ∣y⃗∣|\vec y|∣y​∣ for the spatial radius: ∣(ψ∗g)μν(y)−ημν∣≤C/∣y⃗∣|(\psi^*g)_{\mu\nu}(y)-\eta_{\mu\nu}|\le C/|\vec y|∣(ψ∗g)μν​(y)−ημν​∣≤C/∣y​∣, ∣∂σ(ψ∗g)μν(y)∣≤C/∣y⃗∣2|\partial_\sigma(\psi^*g)_{\mu\nu}(y)|\le C/|\vec y|^2∣∂σ​(ψ∗g)μν​(y)∣≤C/∣y​∣2, ∣(ψ∗F)μν(y)∣≤C/∣y⃗∣2|(\psi^*F)_{\mu\nu}(y)|\le C/|\vec y|^2∣(ψ∗F)μν​(y)∣≤C/∣y​∣2. The region of the end is ψ(ΩR)\psi(\Omega_R)ψ(ΩR​).
  9. DOC =I+(region)∩I−(region)=I^+(\text{region})\cap I^-(\text{region})=I+(region)∩I−(region); black-hole region = complement of I−(region)I^-(\text{region})I−(region); future event horizon = topological frontier of I−(region)I^-(\text{region})I−(region).
  10. Stationary black hole (g,F,T,φ,E)(g,F,T,\varphi,E)(g,F,T,φ,E): all of — ggg Lorentzian; ggg and FFF smooth; FFF a 2-form; TTT a time orientation for ggg; electrovacuum; φ\varphiφ a symmetry flow; DOC globally hyperbolic (curves restricted to DOC); black-hole region nonempty; future event horizon connected (in particular nonempty).
  11. Kerr–Newman DOC with parameters (m,a,e)(m,a,e)(m,a,e) for a set D⊆MD\subseteq MD⊆M: there exist a real β\betaβ and Ψ:R4→M\Psi:\mathbb R^4\to MΨ:R4→M that is a diffeomorphism (in the sense of item 8) from {r>r+}\{r>r_+\}{r>r+​} onto DDD, such that for every yyy with r(y)>r+r(y)>r_+r(y)>r+​: the matrix (Ψ∗g)(y)(\Psi^*g)(y)(Ψ∗g)(y) equals the Kerr–Newman metric matrix at yyy, and (Ψ∗F)(y)=cos⁡β FKN(y)+sin⁡β ⋆FKN(y)(\Psi^*F)(y)=\cos\beta\,F^{KN}(y)+\sin\beta\,\star F^{KN}(y)(Ψ∗F)(y)=cosβFKN(y)+sinβ⋆FKN(y), the Hodge dual taken with the Kerr–Newman metric matrix at yyy. No condition on m,a,em,a,em,a,e is included here.
  12. Static on DDD: for every atlas chart eee, every yyy in its target with e−1(y)∈De^{-1}(y)\in De−1(y)∈D, and all λ,μ,ν\lambda,\mu,\nuλ,μ,ν, with kμk_\mukμ​ the lowered components of the generator of φ\varphiφ along e−1e^{-1}e−1: kλ(∂μkν−∂νkμ)+kμ(∂νkλ−∂λkν)+kν(∂λkμ−∂μkλ)=0k_\lambda(\partial_\mu k_\nu-\partial_\nu k_\mu)+k_\mu(\partial_\nu k_\lambda-\partial_\lambda k_\nu)+k_\nu(\partial_\lambda k_\mu-\partial_\mu k_\lambda)=0kλ​(∂μ​kν​−∂ν​kμ​)+kμ​(∂ν​kλ​−∂λ​kν​)+kν​(∂λ​kμ​−∂μ​kλ​)=0 at yyy.
  13. Axial symmetry ρ\rhoρ (relative to φ\varphiφ): ρ\rhoρ is a symmetry flow of (g,F)(g,F)(g,F), ρ2π=id\rho_{2\pi}=\mathrm{id}ρ2π​=id, ρs∘φt=φt∘ρs\rho_s\circ\varphi_t=\varphi_t\circ\rho_sρs​∘φt​=φt​∘ρs​, some point is fixed by every ρs\rho_sρs​, and some ρs\rho_sρs​ moves some point.
  14. Non-degenerate horizon for (H,D)(H,D)(H,D): there are a symmetry flow χ\chiχ of (g,F)(g,F)(g,F) with generator XXX and a real κ≠0\kappa\neq0κ=0 such that X(x)≠0X(x)\neq0X(x)=0 for some x∈H∩D‾x\in H\cap\overline Dx∈H∩D, and for every x∈H∩D‾x\in H\cap\overline Dx∈H∩D: gx(X,X)=0g_x(X,X)=0gx​(X,X)=0 and, for every atlas chart eee and yyy with e−1(y)=xe^{-1}(y)=xe−1(y)=x and every μ\muμ, ∂μ(z↦g(X,X)(e−1z))(y)=−2κXμ(y)\partial_\mu\big(z\mapsto g(X,X)(e^{-1}z)\big)(y)=-2\kappa X_\mu(y)∂μ​(z↦g(X,X)(e−1z))(y)=−2κXμ​(y) (XμX_\muXμ​ lowered along e−1e^{-1}e−1). Note the chart clause only refers to yyy with e−1(y)=xe^{-1}(y)=xe−1(y)=x; yyy need not lie in the chart target.
Human review
  • Endorsed by Shuze Chen · Oct 8, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Lucas · Oct 8, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paperProve2Me Team
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me