Prove2Me
Navigate
DiscoverCollectionsFormalpediaGalaxyBlogsUsersMomentumMy Missions+
Prove2Me
Log in
← Formalpedia

Lattice superspace and the Wheeler–DeWitt operator

Definition
WheelerDeWittSuperspace

by Mazecto · Oct 5, 2026 · Mathlib 0df444a (Lean v4.33.1)

general-relativitymathematical-physicsquantum-gravity

Lattice model of the configuration space of quantum geometrodynamics. Space is a finite set XXX of sites; a configuration hhh assigns to every site xxx the nine components hab(x)h_{ab}(x)hab​(x) of a spatial metric, all treated as independent real coordinates, and it is physical when every h(x)h(x)h(x) is symmetric positive definite.

The module defines:

  • the volume density det⁡h\sqrt{\det h}deth​;
  • the DeWitt supermetric Gabcd=(hachbd+hadhbc−habhcd)/(2det⁡h)G_{abcd}=(h_{ac}h_{bd}+h_{ad}h_{bc}-h_{ab}h_{cd})/(2\sqrt{\det h})Gabcd​=(hac​hbd​+had​hbc​−hab​hcd​)/(2deth​) and its inverse Gabcd=12det⁡h (hachbd+hadhbc−2habhcd)G^{abcd}=\tfrac12\sqrt{\det h}\,(h^{ac}h^{bd}+h^{ad}h^{bc}-2h^{ab}h^{cd})Gabcd=21​deth​(hachbd+hadhbc−2habhcd), the associated quadratic forms and index maps;
  • the coordinate partial derivatives ∂ab(x)=∂/∂hab(x)\partial_{ab}(x)=\partial/\partial h_{ab}(x)∂ab​(x)=∂/∂hab​(x) as genuine Fréchet derivatives, of complex wavefunctions and of real functions;
  • the kinetic term ∑Gabcd ∂ab∂cdΨ\sum G_{abcd}\,\partial_{ab}\partial_{cd}\Psi∑Gabcd​∂ab​∂cd​Ψ with the supermetric on the left;
  • the vacuum Wheeler–DeWitt operator −2κℏ2 G ∂∂Ψ−det⁡h2κ(R−2Λ)Ψ-2\kappa\hbar^2\,G\,\partial\partial\Psi-\tfrac{\sqrt{\det h}}{2\kappa}(R-2\Lambda)\Psi−2κℏ2G∂∂Ψ−2κdeth​​(R−2Λ)Ψ at each site, with a supplied potential RRR (a lattice scalar curvature, allowed to couple sites);
  • the Einstein–Hamilton–Jacobi functional;
  • the classical lattice constraints Hz=2κ Gabcd(h(z)) pab(z)pcd(z)+Vz(h)H_z=2\kappa\,G_{abcd}(h(z))\,p^{ab}(z)p^{cd}(z)+V_z(h)Hz​=2κGabcd​(h(z))pab(z)pcd(z)+Vz​(h) with Vz=−det⁡h(z)2κ(R−2Λ)V_z=-\tfrac{\sqrt{\det h(z)}}{2\kappa}(R-2\Lambda)Vz​=−2κdeth(z)​​(R−2Λ), and the canonical Poisson bracket on (h,p)(h,p)(h,p).

The cell volume is normalized to one. Coupling constants are named kappa, hbar, Lam.

Definition code
import Mathlib.LinearAlgebra.Matrix.PosDef
import Mathlib.LinearAlgebra.Matrix.NonsingularInverse
import Mathlib.LinearAlgebra.Matrix.Trace
import Mathlib.Analysis.SpecialFunctions.Sqrt
import Mathlib.Analysis.SpecialFunctions.Pow.Real
import Mathlib.Analysis.SpecialFunctions.ExpDeriv
import Mathlib.Analysis.Calculus.IteratedDeriv.Defs

set_option autoImplicit false

namespace WheelerDeWittSuperspace

noncomputable section

open Matrix


/-- Components `h_ab` of a spatial metric at one point, as coordinates on `ℝ⁹`.
All nine entries are independent coordinates; physical metrics are the symmetric
positive-definite ones. -/
abbrev Metric3 := Fin 3 → Fin 3 → ℝ

/-- A configuration: the metric components at each of finitely many spatial sites
(a lattice regularization of a spatial three-geometry). -/
abbrev Config (X : Type*) := X → Metric3

/-- The metric at a site, as a matrix. -/
def metricAt {X : Type*} (h : Config X) (x : X) : Matrix (Fin 3) (Fin 3) ℝ :=
  Matrix.of (h x)

/-- Physical configurations: the metric is symmetric positive definite at every site. -/
def IsPhysical {X : Type*} (h : Config X) : Prop :=
  ∀ x, (metricAt h x).PosDef

/-- Volume density `√h` of a matrix metric. -/
def volume (m : Matrix (Fin 3) (Fin 3) ℝ) : ℝ := Real.sqrt m.det

/-- The DeWitt supermetric `G_abcd = (h_ac h_bd + h_ad h_bc - h_ab h_cd) / (2√h)`
(Kiefer, eq. (5)); it pairs two momenta `p^ab`, `p^cd`. -/
def deWitt (m : Matrix (Fin 3) (Fin 3) ℝ) (a b c d : Fin 3) : ℝ :=
  (m a c * m b d + m a d * m b c - m a b * m c d) / (2 * volume m)

/-- The inverse supermetric `G^abcd = √h (h^ac h^bd + h^ad h^bc - 2 h^ab h^cd) / 2`;
it pairs two velocities `k_ab`, `k_cd`. -/
def invDeWitt (m : Matrix (Fin 3) (Fin 3) ℝ) (a b c d : Fin 3) : ℝ :=
  volume m * (m⁻¹ a c * m⁻¹ b d + m⁻¹ a d * m⁻¹ b c - 2 * m⁻¹ a b * m⁻¹ c d) / 2

/-- `G(p,p) = Σ G_abcd p^ab p^cd` on momenta. -/
def superQuad (m p : Matrix (Fin 3) (Fin 3) ℝ) : ℝ :=
  ∑ a, ∑ b, ∑ c, ∑ d, deWitt m a b c d * p a b * p c d

/-- Index lowering by the supermetric: `(G p)_ab = Σ G_abcd p^cd`. -/
def superMap (m p : Matrix (Fin 3) (Fin 3) ℝ) : Matrix (Fin 3) (Fin 3) ℝ :=
  Matrix.of fun a b => ∑ c, ∑ d, deWitt m a b c d * p c d

/-- Index raising by the inverse supermetric: `(G⁻¹ k)^ab = Σ G^abcd k_cd`. -/
def invSuperMap (m k : Matrix (Fin 3) (Fin 3) ℝ) : Matrix (Fin 3) (Fin 3) ℝ :=
  Matrix.of fun a b => ∑ c, ∑ d, invDeWitt m a b c d * k c d

/-- The supermetric on velocities, `√h (k_ab k^ab - k²)`, written with traces. -/
def velocityQuad (m k : Matrix (Fin 3) (Fin 3) ℝ) : ℝ :=
  volume m * ((m⁻¹ * k * m⁻¹ * k).trace - (m⁻¹ * k).trace ^ 2)

/-- Coordinate direction of `h_ab(x)` in configuration space. -/
def coordDir {X : Type*} [DecidableEq X] (x : X) (a b : Fin 3) : Config X :=
  fun y i j => if y = x ∧ i = a ∧ j = b then 1 else 0

/-- Partial derivative `∂/∂h_ab(x)` of a complex wavefunction (a genuine Fréchet
derivative in a coordinate direction). -/
def partialD {X : Type*} [Fintype X] [DecidableEq X] (Ψ : Config X → ℂ) (x : X) (a b : Fin 3) :
    Config X → ℂ :=
  fun h => fderiv ℝ Ψ h (coordDir x a b)

/-- Partial derivative `∂/∂h_ab(x)` of a real function on configuration space. -/
def partialR {X : Type*} [Fintype X] [DecidableEq X] (S : Config X → ℝ) (x : X) (a b : Fin 3) :
    Config X → ℝ :=
  fun h => fderiv ℝ S h (coordDir x a b)

/-- Kinetic term with the supermetric to the left of both derivatives:
`Σ G_abcd(h(x)) ∂_ab(x) ∂_cd(x) Ψ`. -/
def kinetic {X : Type*} [Fintype X] [DecidableEq X] (Ψ : Config X → ℂ) (h : Config X) (x : X) : ℂ :=
  ∑ a, ∑ b, ∑ c, ∑ d,
    (deWitt (metricAt h x) a b c d : ℂ) * partialD (partialD Ψ x c d) x a b h

/-- The vacuum Wheeler–DeWitt operator at site `x` (Kiefer, eq. (7) with `ρ = 0`):
`-2kappahbar² G ∂∂Ψ - (2kappa)⁻¹ √h (R - 2Lam) Ψ`. The potential `R` is supplied data
(a lattice scalar curvature), allowed to depend on the whole configuration. -/
def wdw {X : Type*} [Fintype X] [DecidableEq X] (kappa hbar Lam : ℝ) (R : Config X → X → ℝ)
    (Ψ : Config X → ℂ) (h : Config X) (x : X) : ℂ :=
  -2 * (kappa : ℂ) * (hbar : ℂ) ^ 2 * kinetic Ψ h x -
    ((2 * kappa)⁻¹ * volume (metricAt h x) * (R h x - 2 * Lam) : ℝ) * Ψ h

/-- The Einstein–Hamilton–Jacobi functional: the classical Hamiltonian constraint
with momenta `p^ab = ∂S/∂h_ab`. -/
def hamiltonJacobi {X : Type*} [Fintype X] [DecidableEq X] (kappa Lam : ℝ) (R : Config X → X → ℝ)
    (S : Config X → ℝ) (h : Config X) (x : X) : ℝ :=
  2 * kappa * (∑ a, ∑ b, ∑ c, ∑ d,
      deWitt (metricAt h x) a b c d * partialR S x a b h * partialR S x c d h) -
    (2 * kappa)⁻¹ * volume (metricAt h x) * (R h x - 2 * Lam)


/-! ### Constraint algebra on the lattice -/

/-- The potential part of the constraint at site `z`: `-(2κ)⁻¹ √h(z) (R(h,z) - 2Λ)`.
It depends on the whole configuration as soon as `R` couples neighbouring sites. -/
def potentialAt {X : Type*} (kappa Lam : ℝ) (R : Config X → X → ℝ) (h : Config X) (z : X) : ℝ :=
  -(2 * kappa)⁻¹ * volume (metricAt h z) * (R h z - 2 * Lam)

/-- Classical phase space: metric components and conjugate momenta at every site. -/
abbrev Phase (X : Type*) := Config X × Config X

/-- The classical lattice Hamiltonian constraint at site `z`:
`2κ G_abcd(h(z)) p^ab(z) p^cd(z) + V_z(h)`. -/
def classicalConstraint {X : Type*} (kappa Lam : ℝ) (R : Config X → X → ℝ) (z : X)
    (w : Phase X) : ℝ :=
  2 * kappa * (∑ a, ∑ b, ∑ c, ∑ d,
      deWitt (metricAt w.1 z) a b c d * w.2 z a b * w.2 z c d) +
    potentialAt kappa Lam R w.1 z

/-- Canonical Poisson bracket on the lattice phase space,
`{F,G} = Σ_{z,a,b} (∂F/∂h_ab(z) ∂G/∂p^ab(z) - ∂F/∂p^ab(z) ∂G/∂h_ab(z))`. -/
def poisson {X : Type*} [Fintype X] [DecidableEq X] (F G : Phase X → ℝ) (w : Phase X) : ℝ :=
  ∑ z, ∑ a, ∑ b,
    (fderiv ℝ F w (coordDir z a b, 0) * fderiv ℝ G w (0, coordDir z a b) -
      fderiv ℝ F w (0, coordDir z a b) * fderiv ℝ G w (coordDir z a b, 0))

end

end WheelerDeWittSuperspace
Source
C. Kiefer, Quantum Geometrodynamics: whence, whither?, Gen. Relativ. Gravit. 41 (2009) 877-901, https://arxiv.org/abs/0812.0295; B. S. DeWitt, Quantum Theory of Gravity I. The Canonical Theory, Phys. Rev. 160 (1967) 1113-1148, https://doi.org/10.1103/PhysRev.160.1113
Read-back

What the Lean code literally says, in plain math · claude-opus-5-5

Read-back of the definition module WheelerDeWittSuperspace. All definitions below live in one namespace. Throughout, indices a,b,c,d,i,ja,b,c,d,i,ja,b,c,d,i,j range over {0,1,2}\{0,1,2\}{0,1,2}, and XXX is an arbitrary type of "sites" (any universe). Unless stated otherwise, no finiteness, nonemptiness or decidability is assumed on XXX.

Global conventions (junk values). The definitions use total functions, so the following conventions apply silently everywhere:

  • t\sqrt{t}t​ is the real square root with t=0\sqrt{t}=0t​=0 for every t<0t<0t<0.
  • Real (and complex) division by 000 and the inverse 0−10^{-1}0−1 both give 000; in particular (2κ)−1=0(2\kappa)^{-1}=0(2κ)−1=0 when κ=0\kappa=0κ=0.
  • For a real 3×33\times 33×3 matrix mmm, m−1m^{-1}m−1 means (det⁡m)−1 adj(m)(\det m)^{-1}\,\mathrm{adj}(m)(detm)−1adj(m). If det⁡m=0\det m = 0detm=0 this is the zero matrix. No symmetry or positivity of mmm is needed for m−1m^{-1}m−1 to be defined.
  • The Fréchet derivative Df(h)Df(h)Df(h) of a function at a point where it is not differentiable is the zero linear map, so every directional derivative built from it is 000 there.

1. Metric components Metric3. This is the type of all functions {0,1,2}×{0,1,2}→R\{0,1,2\}\times\{0,1,2\}\to\mathbb{R}{0,1,2}×{0,1,2}→R, i.e. arrays (hab)(h_{ab})(hab​) of nine real numbers, identified with R9\mathbb{R}^9R9. All nine entries are independent: there is no symmetry condition (habh_{ab}hab​ and hbah_{ba}hba​ are distinct coordinates) and no positivity or nondegeneracy condition. As a normed space it carries the product structure (sup norm).

2. Configurations Config X. A configuration is any function h:X→h : X \toh:X→ Metric3, i.e. a family h(x)=(hab(x))h(x) = (h_{ab}(x))h(x)=(hab​(x)) of nine arbitrary real numbers per site x∈Xx\in Xx∈X. If XXX is empty there is exactly one configuration (the empty family).

3. metricAt. For a configuration hhh and a site xxx, metricAt⁡(h,x)\operatorname{metricAt}(h,x)metricAt(h,x) is the 3×33\times33×3 real matrix whose (a,b)(a,b)(a,b) entry is hab(x)h_{ab}(x)hab​(x). Below I write h(x)h(x)h(x) for this matrix.

4. IsPhysical. For a configuration hhh on XXX, IsPhysical⁡(h)\operatorname{IsPhysical}(h)IsPhysical(h) is the proposition that for every site x∈Xx\in Xx∈X the matrix h(x)h(x)h(x) is positive definite in Mathlib's sense: it is symmetric (hab(x)=hba(x)h_{ab}(x)=h_{ba}(x)hab​(x)=hba​(x) for all a,ba,ba,b, real Hermitian) and ∑i,jvi hij(x) vj>0\sum_{i,j} v_i\,h_{ij}(x)\,v_j > 0∑i,j​vi​hij​(x)vj​>0 for every nonzero v∈R3v\in\mathbb{R}^3v∈R3. If XXX is empty this holds vacuously. This predicate is only defined here; none of the other definitions assumes it.

5. volume. For any real 3×33\times33×3 matrix mmm,

vol⁡(m)=det⁡m,\operatorname{vol}(m) = \sqrt{\det m},vol(m)=detm​,

which equals 000 whenever det⁡m≤0\det m \le 0detm≤0 (including all non-symmetric or indefinite matrices with negative determinant).

6. deWitt. For any real 3×33\times33×3 matrix mmm and indices a,b,c,da,b,c,da,b,c,d,

Gabcd(m)=macmbd+madmbc−mabmcd2vol⁡(m).G_{abcd}(m) = \frac{m_{ac}m_{bd} + m_{ad}m_{bc} - m_{ab}m_{cd}}{2\operatorname{vol}(m)} .Gabcd​(m)=2vol(m)mac​mbd​+mad​mbc​−mab​mcd​​.

The entries are used exactly as given (no symmetrisation of mmm). If vol⁡(m)=0\operatorname{vol}(m)=0vol(m)=0 (i.e. det⁡m≤0\det m\le 0detm≤0), then Gabcd(m)=0G_{abcd}(m)=0Gabcd​(m)=0 for all indices.

7. invDeWitt. For any real 3×33\times33×3 matrix mmm, writing n=m−1n = m^{-1}n=m−1 (the matrix inverse with the zero-matrix convention above),

Gabcd(m)=vol⁡(m) (nacnbd+nadnbc−2 nabncd)2.G^{abcd}(m) = \frac{\operatorname{vol}(m)\,\bigl(n_{ac}n_{bd} + n_{ad}n_{bc} - 2\,n_{ab}n_{cd}\bigr)}{2}.Gabcd(m)=2vol(m)(nac​nbd​+nad​nbc​−2nab​ncd​)​.

If det⁡m≤0\det m \le 0detm≤0 then vol⁡(m)=0\operatorname{vol}(m)=0vol(m)=0, so Gabcd(m)=0G^{abcd}(m)=0Gabcd(m)=0 (and if det⁡m=0\det m=0detm=0 also n=0n=0n=0). Note the coefficient of the last term is 222 here, versus 111 in GabcdG_{abcd}Gabcd​.

8. superQuad. For real 3×33\times33×3 matrices mmm and ppp (with ppp arbitrary, not required symmetric),

superQuad⁡(m,p)=∑a,b,c,dGabcd(m) pab pcd.\operatorname{superQuad}(m,p) = \sum_{a,b,c,d} G_{abcd}(m)\,p_{ab}\,p_{cd}.superQuad(m,p)=a,b,c,d∑​Gabcd​(m)pab​pcd​.

9. superMap. For real 3×33\times33×3 matrices m,pm,pm,p, the matrix with (a,b)(a,b)(a,b) entry

(superMap⁡(m,p))ab=∑c,dGabcd(m) pcd.\bigl(\operatorname{superMap}(m,p)\bigr)_{ab} = \sum_{c,d} G_{abcd}(m)\,p_{cd}.(superMap(m,p))ab​=c,d∑​Gabcd​(m)pcd​.

10. invSuperMap. For real 3×33\times33×3 matrices m,km,km,k, the matrix with (a,b)(a,b)(a,b) entry

(invSuperMap⁡(m,k))ab=∑c,dGabcd(m) kcd.\bigl(\operatorname{invSuperMap}(m,k)\bigr)_{ab} = \sum_{c,d} G^{abcd}(m)\,k_{cd}.(invSuperMap(m,k))ab​=c,d∑​Gabcd(m)kcd​.

11. velocityQuad. For real 3×33\times33×3 matrices m,km,km,k, with m−1m^{-1}m−1 as above (zero if det⁡m=0\det m=0detm=0),

velocityQuad⁡(m,k)=vol⁡(m)(tr⁡(m−1k m−1k)−(tr⁡(m−1k))2),\operatorname{velocityQuad}(m,k) = \operatorname{vol}(m)\Bigl(\operatorname{tr}\bigl(m^{-1}k\,m^{-1}k\bigr) - \bigl(\operatorname{tr}(m^{-1}k)\bigr)^2\Bigr),velocityQuad(m,k)=vol(m)(tr(m−1km−1k)−(tr(m−1k))2),

using ordinary matrix products in the order written; kkk is not required to be symmetric. It is 000 whenever det⁡m≤0\det m \le 0detm≤0.

12. coordDir. Assuming decidable equality on XXX, for a site xxx and indices a,ba,ba,b, coordDir⁡(x,a,b)\operatorname{coordDir}(x,a,b)coordDir(x,a,b) is the configuration eab(x)e^{(x)}_{ab}eab(x)​ given by

(eab(x))ij(y)={1if y=x, i=a, j=b,0otherwise.\bigl(e^{(x)}_{ab}\bigr)_{ij}(y) = \begin{cases}1 & \text{if } y=x,\ i=a,\ j=b,\\ 0 & \text{otherwise.}\end{cases}(eab(x)​)ij​(y)={10​if y=x, i=a, j=b,otherwise.​

It is the unit vector of the single coordinate hab(x)h_{ab}(x)hab​(x); in particular for a≠ba\ne ba=b it does not also move hba(x)h_{ba}(x)hba​(x).

13. partialD. Assume XXX is finite with decidable equality (so Config X ≅R9∣X∣\cong \mathbb{R}^{9|X|}≅R9∣X∣ is a finite-dimensional real normed space). For a function Ψ:\Psi :Ψ: Config X →C\to\mathbb{C}→C (no regularity assumed), a site xxx and indices a,ba,ba,b, ∂ab(x)Ψ\partial_{ab}(x)\Psi∂ab​(x)Ψ is the function on configurations

h↦DΨ(h)[eab(x)],h \mapsto D\Psi(h)\bigl[e^{(x)}_{ab}\bigr],h↦DΨ(h)[eab(x)​],

where DΨ(h)D\Psi(h)DΨ(h) is the real Fréchet derivative of Ψ\PsiΨ at hhh (with C\mathbb{C}C viewed as a real normed space). At any hhh where Ψ\PsiΨ is not Fréchet differentiable, this value is 000 (even if the one-variable partial derivative in that direction happens to exist).

14. partialR. Same as partialD, but for a real-valued S:S :S: Config X →R\to\mathbb{R}→R: ∂ab(x)S(h)=DS(h)[eab(x)]\partial_{ab}(x)S(h) = DS(h)[e^{(x)}_{ab}]∂ab​(x)S(h)=DS(h)[eab(x)​], equal to 000 wherever SSS is not Fréchet differentiable.

15. kinetic. With XXX finite and decidable equality, for Ψ:\Psi :Ψ: Config X →C\to\mathbb{C}→C, a configuration hhh and a site xxx,

kinetic⁡(Ψ,h,x)=∑a,b,c,dGabcd(h(x))⋅∂ab(x)(∂cd(x)Ψ)(h) ∈C,\operatorname{kinetic}(\Psi,h,x) = \sum_{a,b,c,d} G_{abcd}\bigl(h(x)\bigr)\cdot \partial_{ab}(x)\bigl(\partial_{cd}(x)\Psi\bigr)(h)\ \in\mathbb{C},kinetic(Ψ,h,x)=a,b,c,d∑​Gabcd​(h(x))⋅∂ab​(x)(∂cd​(x)Ψ)(h) ∈C,

where the real coefficient Gabcd(h(x))G_{abcd}(h(x))Gabcd​(h(x)) is cast into C\mathbb{C}C. The inner derivative is in direction (c,d)(c,d)(c,d) and the outer one in direction (a,b)(a,b)(a,b); the outer derivative is the Fréchet derivative of the function h′↦∂cd(x)Ψ(h′)h'\mapsto \partial_{cd}(x)\Psi(h')h′↦∂cd​(x)Ψ(h′), which is 000 wherever that function is not differentiable. The coefficient is evaluated at h(x)h(x)h(x) and placed outside both derivatives. All coefficients vanish if det⁡h(x)≤0\det h(x)\le 0deth(x)≤0.

16. wdw. With XXX finite and decidable equality, for real parameters κ\kappaκ (kappa), ℏ\hbarℏ (hbar), Λ\LambdaΛ (Lam), an arbitrary function R:R :R: Config X ×X→R\times X\to\mathbb{R}×X→R (any "potential" depending on the whole configuration and on the site; no relation to curvature is imposed), Ψ:\Psi :Ψ: Config X →C\to\mathbb{C}→C, a configuration hhh and a site xxx,

wdw⁡(Ψ,h,x)=−2 κ ℏ2 kinetic⁡(Ψ,h,x)  −  [(2κ)−1vol⁡(h(x))(R(h,x)−2Λ)] Ψ(h) ∈C,\operatorname{wdw}(\Psi,h,x) = -2\,\kappa\,\hbar^2\,\operatorname{kinetic}(\Psi,h,x) \;-\; \Bigl[(2\kappa)^{-1}\operatorname{vol}\bigl(h(x)\bigr)\bigl(R(h,x)-2\Lambda\bigr)\Bigr]\,\Psi(h)\ \in\mathbb{C},wdw(Ψ,h,x)=−2κℏ2kinetic(Ψ,h,x)−[(2κ)−1vol(h(x))(R(h,x)−2Λ)]Ψ(h) ∈C,

where the bracket is computed in R\mathbb{R}R and cast to C\mathbb{C}C. No sign or nonzero assumptions on κ,ℏ,Λ\kappa,\hbar,\Lambdaκ,ℏ,Λ: if κ=0\kappa=0κ=0 both terms vanish (since (2⋅0)−1=0(2\cdot 0)^{-1}=0(2⋅0)−1=0) and wdw⁡=0\operatorname{wdw}=0wdw=0; if ℏ=0\hbar=0ℏ=0 only the potential term remains.

17. hamiltonJacobi. With XXX finite and decidable equality, for real κ,Λ\kappa,\Lambdaκ,Λ, arbitrary RRR as above, S:S :S: Config X →R\to\mathbb{R}→R, a configuration hhh and a site xxx,

HJ⁡(S,h,x)=2κ∑a,b,c,dGabcd(h(x)) ∂ab(x)S(h) ∂cd(x)S(h)  −  (2κ)−1vol⁡(h(x))(R(h,x)−2Λ) ∈R.\operatorname{HJ}(S,h,x) = 2\kappa\sum_{a,b,c,d} G_{abcd}\bigl(h(x)\bigr)\,\partial_{ab}(x)S(h)\,\partial_{cd}(x)S(h) \;-\; (2\kappa)^{-1}\operatorname{vol}\bigl(h(x)\bigr)\bigl(R(h,x)-2\Lambda\bigr)\ \in\mathbb{R}.HJ(S,h,x)=2κa,b,c,d∑​Gabcd​(h(x))∂ab​(x)S(h)∂cd​(x)S(h)−(2κ)−1vol(h(x))(R(h,x)−2Λ) ∈R.

Here there is no ℏ\hbarℏ factor; the first-derivative partials are those of item 14 (zero where SSS is not differentiable). If κ=0\kappa=0κ=0 the value is 000.

18. potentialAt. For any XXX (no finiteness needed), real κ,Λ\kappa,\Lambdaκ,Λ, arbitrary RRR, a configuration hhh and a site zzz,

Vz(h)=−(2κ)−1vol⁡(h(z))(R(h,z)−2Λ),V_z(h) = -(2\kappa)^{-1}\operatorname{vol}\bigl(h(z)\bigr)\bigl(R(h,z)-2\Lambda\bigr),Vz​(h)=−(2κ)−1vol(h(z))(R(h,z)−2Λ),

which is 000 if κ=0\kappa=0κ=0 or det⁡h(z)≤0\det h(z)\le 0deth(z)≤0.

19. Phase space Phase X. The type of pairs w=(h,p)w=(h,p)w=(h,p) of two configurations on XXX: hhh ("metric") and ppp ("momenta"), each assigning nine arbitrary real numbers to every site, with no symmetry condition. As a normed space it is the product of two copies of Config X.

20. classicalConstraint. For any XXX, real κ,Λ\kappa,\Lambdaκ,Λ, arbitrary RRR, a site zzz and a phase-space point w=(h,p)w=(h,p)w=(h,p),

Cz(w)=2κ∑a,b,c,dGabcd(h(z)) pab(z) pcd(z)  +  Vz(h),C_z(w) = 2\kappa\sum_{a,b,c,d} G_{abcd}\bigl(h(z)\bigr)\,p_{ab}(z)\,p_{cd}(z) \;+\; V_z(h),Cz​(w)=2κa,b,c,d∑​Gabcd​(h(z))pab​(z)pcd​(z)+Vz​(h),

with VzV_zVz​ from item 18 (so the potential enters with the minus sign built into VzV_zVz​). It is 000 if κ=0\kappa=0κ=0; its kinetic part is 000 if det⁡h(z)≤0\det h(z)\le0deth(z)≤0.

21. poisson. With XXX finite and decidable equality, for arbitrary functions F,G:F,G :F,G: Phase X →R\to\mathbb{R}→R (no regularity assumed) and a point www,

{F,G}(w)=∑z∈X∑a,b(DF(w)[(eab(z),0)]  DG(w)[(0,eab(z))]  −  DF(w)[(0,eab(z))]  DG(w)[(eab(z),0)]),\{F,G\}(w) = \sum_{z\in X}\sum_{a,b}\Bigl( DF(w)\bigl[(e^{(z)}_{ab},0)\bigr]\; DG(w)\bigl[(0,e^{(z)}_{ab})\bigr] \;-\; DF(w)\bigl[(0,e^{(z)}_{ab})\bigr]\; DG(w)\bigl[(e^{(z)}_{ab},0)\bigr]\Bigr),{F,G}(w)=z∈X∑​a,b∑​(DF(w)[(eab(z)​,0)]DG(w)[(0,eab(z)​)]−DF(w)[(0,eab(z)​)]DG(w)[(eab(z)​,0)]),

where DDD is the Fréchet derivative on the product space, (e,0)(e,0)(e,0) moves only the metric coordinate hab(z)h_{ab}(z)hab​(z) and (0,e)(0,e)(0,e) only the momentum coordinate pab(z)p_{ab}(z)pab​(z). All 9 index pairs (a,b)(a,b)(a,b) are summed separately (the (a,b)(a,b)(a,b) and (b,a)(b,a)(b,a) coordinates are independent). If FFF (resp. GGG) is not differentiable at www, its derivative there is 000 and the corresponding products vanish. If XXX is empty, {F,G}=0\{F,G\}=0{F,G}=0 identically.


Edge cases summarised. (i) Nothing outside IsPhysical restricts metrics: every definition accepts non-symmetric, singular or indefinite arrays, where vol⁡=0\operatorname{vol}=0vol=0 forces Gabcd=Gabcd=velocityQuad⁡=0G_{abcd}=G^{abcd}=\operatorname{velocityQuad}=0Gabcd​=Gabcd=velocityQuad=0, and the kinetic and potential terms vanish. (ii) All nine entries habh_{ab}hab​ (and pabp_{ab}pab​) are separate coordinates, and partial derivatives are taken along single entries. (iii) Derivative-based quantities silently return 000 at non-differentiable points. (iv) κ=0\kappa=0κ=0 kills every term carrying (2κ)−1(2\kappa)^{-1}(2κ)−1 or the factor 2κ2\kappa2κ. (v) For empty XXX, there is one configuration, IsPhysical holds vacuously, site-indexed quantities have no site to be evaluated at, and the Poisson bracket is identically 000.

Human review
  • Endorsed by Shuze Chen · Oct 8, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Mazecto · 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