DeWitt supermetric, metric traces and the ADM conjugate momentum (Isham §3.3.3)
DefinitionIshamProblemOfTime_DeWittPointwise algebra of the ADM canonical formalism of general relativity, as in Isham, §3.3.3. Fix a point of the spatial slice . The spatial metric at is a real matrix (assumed positive definite in every theorem that uses it), its inverse is the matrix inverse , and . Contravariant two-tensors such as the momentum and the extrinsic curvature are real matrices.
- DeWitt supermetric (eq. (3.3.21)): the bilinear form
It is the coefficient of the kinetic term of the super-Hamiltonian (eq. (3.3.20)).
- Metric trace (Isham's ) and metric square .
- Conjugate momentum (eq. (3.3.16)): , where .
Conventions of the Lean file: is written ; for a matrix with this evaluates to (so the form vanishes identically), which is why every theorem assumes positive definite.
import Mathlib
/-!
# DeWitt supermetric (Isham, *Canonical Quantum Gravity and the Problem of Time*, §3.3.3)
Pointwise (fixed `x ∈ Σ`) algebra of the ADM canonical variables. A spatial three-metric
at a point is a real `3 × 3` matrix `g` (components `g_{ab}`, assumed positive definite in
the theorems), and contravariant symmetric tensors such as `p^{ab}` and `K^{ab}` are real
`3 × 3` matrices. The inverse metric `g^{ab}` is the matrix inverse `g⁻¹`, and `|g|` is
`det g`.
-/
namespace IshamProblemOfTime
open Matrix
/-- Isham eq. (3.3.21): the DeWitt supermetric, as a bilinear form on contravariant
two-tensors,
`𝒢_{abcd} p^{ab} q^{cd}` with
`𝒢_{abcd} = ½ |g|^{-1/2} (g_{ac} g_{bd} + g_{bc} g_{ad} - g_{ab} g_{cd})`. -/
noncomputable def deWittSupermetric (g p q : Matrix (Fin 3) (Fin 3) ℝ) : ℝ :=
(1 / 2) * (Real.sqrt g.det)⁻¹ *
∑ a : Fin 3, ∑ b : Fin 3, ∑ c : Fin 3, ∑ d : Fin 3,
(g a c * g b d + g b c * g a d - g a b * g c d) * p a b * q c d
/-- The trace `g_{ab} T^{ab}` (written `T_c{}^c` in Isham) of a contravariant two-tensor
`T` with respect to the metric `g`. -/
def metricTrace (g T : Matrix (Fin 3) (Fin 3) ℝ) : ℝ :=
∑ a : Fin 3, ∑ b : Fin 3, g a b * T a b
/-- The full contraction `K_{ab} K^{ab} = g_{ac} g_{bd} K^{ab} K^{cd}` of a contravariant
two-tensor `K` with itself, indices lowered with `g`. -/
def metricSquare (g K : Matrix (Fin 3) (Fin 3) ℝ) : ℝ :=
∑ a : Fin 3, ∑ b : Fin 3, ∑ c : Fin 3, ∑ d : Fin 3, g a c * g b d * K a b * K c d
/-- Isham eq. (3.3.16): the momentum conjugate to `g_{ab}` expressed through the extrinsic
curvature,
`p^{ab} = -(|g|^{1/2} / κ²) (K^{ab} - g^{ab} K_c{}^c)`. -/
noncomputable def conjugateMomentum (κ : ℝ) (g K : Matrix (Fin 3) (Fin 3) ℝ) :
Matrix (Fin 3) (Fin 3) ℝ :=
(-(Real.sqrt g.det / κ ^ 2)) • (K - metricTrace g K • g⁻¹)
end IshamProblemOfTime
Read-back
What the Lean code literally says, in plain math · Aristotle (Harmonic) - NON-BLIND: same agent that drafted the statements
Non-blind read-back — not independent testimony. This read-back was written by the same agent (Aristotle, by Harmonic) that drafted the Lean statements, with full knowledge of the source paper and of the intended meaning. It was not produced blind by an independent auditor, so it must not be mistaken for independent testimony; reviewers should treat it as the drafter's own rendering of the code and check the Lean text directly.
The file introduces four objects, all built from real matrices (indices range over three values). No hypothesis (symmetry, positivity, invertibility) is built into any of them.
1. deWittSupermetric. For arbitrary real matrices it returns the real number
Edge cases: the real square root returns for a negative argument and the inverse of is taken to be , so whenever the value is for all . The matrices need not be symmetric; the formula is applied verbatim to all nine entries.
2. metricTrace. For matrices : (no symmetry assumed).
3. metricSquare. For matrices : .
4. conjugateMomentum. For a real number and matrices it returns the matrix
where is the matrix inverse (which is the zero matrix when ). Division by returns when , in which case is the zero matrix; when the factor is and again .