Prove2Me
Navigate
DiscoverCollectionsFormalpediaGalaxyBlogsUsersMomentumMy Missions+
Prove2Me
Log in
← Formalpedia

Commutator of two lattice Wheeler–DeWitt constraints

Proved
WheelerDeWittSuperspace.wdw_commutator

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

general-relativitymathematical-physicsquantum-gravity

Let Vz(h)=−det⁡h(z)2κ(R(h,z)−2Λ)V_z(h)=-\tfrac{\sqrt{\det h(z)}}{2\kappa}(R(h,z)-2\Lambda)Vz​(h)=−2κdeth(z)​​(R(h,z)−2Λ) with each R(⋅,z)R(\cdot,z)R(⋅,z) of class C2C^2C2, let Ψ\PsiΨ be C4C^4C4, let hhh be physical and x≠yx\neq yx=y. Then

[H^x,H^y]Ψ=−2κℏ2(Gabcd(h(x))(2 ∂abxVy ∂cdxΨ+∂abx∂cdxVy Ψ)−(x↔y)),[\widehat{\mathcal H}_x,\widehat{\mathcal H}_y]\Psi=-2\kappa\hbar^2\Big(G_{abcd}(h(x))\big(2\,\partial^{x}_{ab}V_y\,\partial^{x}_{cd}\Psi+\partial^{x}_{ab}\partial^{x}_{cd}V_y\,\Psi\big)-(x\leftrightarrow y)\Big),[Hx​,Hy​]Ψ=−2κℏ2(Gabcd​(h(x))(2∂abx​Vy​∂cdx​Ψ+∂abx​∂cdx​Vy​Ψ)−(x↔y)),

with ∂abx=∂/∂hab(x)\partial^{x}_{ab}=\partial/\partial h_{ab}(x)∂abx​=∂/∂hab​(x). The commutator is a first-order operator whose coefficients measure how the potential at one site depends on the metric at the other.

Preamble
import Definitions.Def_WheelerDeWittSuperspace

set_option autoImplicit false

open Matrix
Formal statement
namespace WheelerDeWittSuperspace

/-- Milestone 6 (quantum lattice commutator): two lattice Wheeler–DeWitt constraints
fail to commute by an explicit first-order operator built from the cross-derivatives
of the potentials. -/
theorem wdw_commutator {X : Type*} [Fintype X] [DecidableEq X] (kappa hbar Lam : ℝ)
    (hkappa : kappa ≠ 0) (hhbar : hbar ≠ 0)
    (R : Config X → X → ℝ) (hR : ∀ z, ContDiff ℝ 2 (fun h => R h z))
    (Ψ : Config X → ℂ) (hΨ : ContDiff ℝ 4 Ψ) (h : Config X) (hh : IsPhysical h)
    (x y : X) (hxy : x ≠ y) :
    wdw kappa hbar Lam R (fun h' => wdw kappa hbar Lam R Ψ h' y) h x -
      wdw kappa hbar Lam R (fun h' => wdw kappa hbar Lam R Ψ h' x) h y =
    -2 * (kappa : ℂ) * (hbar : ℂ) ^ 2 *
      ((∑ a, ∑ b, ∑ c, ∑ d, (deWitt (metricAt h x) a b c d : ℂ) *
          (2 * (partialR (fun h' => potentialAt kappa Lam R h' y) x a b h : ℂ) *
              partialD Ψ x c d h +
            (partialR (partialR (fun h' => potentialAt kappa Lam R h' y) x c d) x a b h : ℂ) *
              Ψ h)) -
       (∑ a, ∑ b, ∑ c, ∑ d, (deWitt (metricAt h y) a b c d : ℂ) *
          (2 * (partialR (fun h' => potentialAt kappa Lam R h' x) y a b h : ℂ) *
              partialD Ψ y c d h +
            (partialR (partialR (fun h' => potentialAt kappa Lam R h' x) y c d) y a b h : ℂ) *
              Ψ h))) := by sorry

end WheelerDeWittSuperspace
Source
R. Loll, Discrete approaches to quantum gravity in four dimensions, Living Rev. Relativ. 1 (1998), Section 2.12, https://arxiv.org/abs/gr-qc/9805049; B. Dittrich, https://arxiv.org/abs/0810.3594, p. 6; B. Bahr, B. Dittrich, https://arxiv.org/abs/0905.1670; M. Bander, Phys. Rev. D 36 (1987) 2297
Read-back

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

Read-back of wdw_commutator. This theorem gives an exact formula, at one physical configuration, for the commutator of two lattice Wheeler–DeWitt operators at two distinct sites.

Setting and binders. Let XXX be any finite type with decidable equality (the "sites"). A configuration is a map h:X→R3×3h : X \to \mathbb{R}^{3\times 3}h:X→R3×3. All nine entries hab(x)h_{ab}(x)hab​(x), a,b∈{0,1,2}a,b\in\{0,1,2\}a,b∈{0,1,2}, count as independent real coordinates, so symmetry is not built in. The configuration space C=RX×3×3\mathcal{C}=\mathbb{R}^{X\times 3\times 3}C=RX×3×3 is a finite-dimensional real normed space with the sup norm. All derivatives below are real Fréchet derivatives on C\mathcal{C}C. The statement quantifies over:

  • real numbers κ,ℏ,Λ\kappa,\hbar,\Lambdaκ,ℏ,Λ with κ≠0\kappa\neq 0κ=0 and ℏ≠0\hbar\neq 0ℏ=0. Λ\LambdaΛ is unrestricted;
  • a real "potential datum" R:C×X→RR:\mathcal{C}\times X\to\mathbb{R}R:C×X→R, (h,z)↦R(h,z)(h,z)\mapsto R(h,z)(h,z)↦R(h,z), which may depend on the whole configuration. The only assumption is that for every site zzz the map h↦R(h,z)h\mapsto R(h,z)h↦R(h,z) is C2C^2C2 on all of C\mathcal{C}C;
  • a complex-valued function Ψ:C→C\Psi:\mathcal{C}\to\mathbb{C}Ψ:C→C that is C4C^4C4 as a map between real vector spaces;
  • a configuration hhh that is physical: for every site xxx, the matrix h(x)h(x)h(x) is symmetric positive definite (symmetric, and v⊤h(x)v>0v^\top h(x) v>0v⊤h(x)v>0 for all v≠0v\neq 0v=0). In particular det⁡h(x)>0\det h(x)>0deth(x)>0 at every site;
  • two sites x,y∈Xx,y\in Xx,y∈X with x≠yx\neq yx=y.

If XXX has fewer than two elements, no such x,yx,yx,y exist and the statement is vacuous.

Definitions unfolded.

  • Coordinate derivatives. eab(x)∈Ce^{(x)}_{ab}\in\mathcal{C}eab(x)​∈C is the configuration that is 111 in entry (a,b)(a,b)(a,b) at site xxx and 000 everywhere else. For a function FFF on C\mathcal{C}C (real or complex valued), set
∂ab(x)F(h):=DF(h)[eab(x)],\partial^{(x)}_{ab}F(h) := DF(h)\big[e^{(x)}_{ab}\big],∂ab(x)​F(h):=DF(h)[eab(x)​],

the full Fréchet derivative at hhh applied to that direction. By the Fréchet-derivative convention, DF(h)DF(h)DF(h) is taken to be 000 wherever FFF is not differentiable. Iterated derivatives are nested: ∂ab(x)∂cd(x)F\partial^{(x)}_{ab}\partial^{(x)}_{cd}F∂ab(x)​∂cd(x)​F means first form h′↦∂cd(x)F(h′)h'\mapsto\partial^{(x)}_{cd}F(h')h′↦∂cd(x)​F(h′), then take ∂ab(x)\partial^{(x)}_{ab}∂ab(x)​ of that function.

  • Volume density. h(z):=det⁡h(z)\sqrt{h(z)} := \sqrt{\det h(z)}h(z)​:=deth(z)​, with the real square root, which is 000 for negative arguments.
  • DeWitt supermetric at site zzz, for the matrix m=h(z)m=h(z)m=h(z):
Gabcd(m)=macmbd+madmbc−mabmcd2det⁡m,G_{abcd}(m)=\frac{m_{ac}m_{bd}+m_{ad}m_{bc}-m_{ab}m_{cd}}{2\sqrt{\det m}},Gabcd​(m)=2detm​mac​mbd​+mad​mbc​−mab​mcd​​,

with the convention that division by 000 gives 000.

  • Potential at site zzz:
Vz(h):=−12κ det⁡h(z) (R(h,z)−2Λ).V_z(h) := -\frac{1}{2\kappa}\,\sqrt{\det h(z)}\,\big(R(h,z)-2\Lambda\big).Vz​(h):=−2κ1​deth(z)​(R(h,z)−2Λ).
  • Lattice Wheeler–DeWitt operator at site zzz, acting on any Φ:C→C\Phi:\mathcal{C}\to\mathbb{C}Φ:C→C:
(HzΦ)(h):=−2κℏ2∑a,b,c,d=02Gabcd(h(z))  ∂ab(z)∂cd(z)Φ(h)  −  12κdet⁡h(z) (R(h,z)−2Λ) Φ(h).(\mathcal{H}_z\Phi)(h) := -2\kappa\hbar^2\sum_{a,b,c,d=0}^{2} G_{abcd}(h(z))\;\partial^{(z)}_{ab}\partial^{(z)}_{cd}\Phi(h)\;-\;\frac{1}{2\kappa}\sqrt{\det h(z)}\,\big(R(h,z)-2\Lambda\big)\,\Phi(h).(Hz​Φ)(h):=−2κℏ2a,b,c,d=0∑2​Gabcd​(h(z))∂ab(z)​∂cd(z)​Φ(h)−2κ1​deth(z)​(R(h,z)−2Λ)Φ(h).

The supermetric stands to the left of both derivatives, and the second term is exactly Vz(h) Φ(h)V_z(h)\,\Phi(h)Vz​(h)Φ(h). Real quantities are embedded in C\mathbb{C}C.

Assertion. At the given physical hhh and the given sites x≠yx\neq yx=y:

(Hx(HyΨ))(h)−(Hy(HxΨ))(h)=−2κℏ2(Tx,y−Ty,x),\big(\mathcal{H}_x(\mathcal{H}_y\Psi)\big)(h)-\big(\mathcal{H}_y(\mathcal{H}_x\Psi)\big)(h) = -2\kappa\hbar^2\Big(T_{x,y}-T_{y,x}\Big),(Hx​(Hy​Ψ))(h)−(Hy​(Hx​Ψ))(h)=−2κℏ2(Tx,y​−Ty,x​),

where

Tx,y=∑a,b,c,d=02Gabcd(h(x))(2 ∂ab(x)Vy(h)  ∂cd(x)Ψ(h)  +  ∂ab(x)(∂cd(x)Vy)(h)  Ψ(h))T_{x,y}=\sum_{a,b,c,d=0}^{2}G_{abcd}(h(x))\Big(2\,\partial^{(x)}_{ab}V_y(h)\;\partial^{(x)}_{cd}\Psi(h)\;+\;\partial^{(x)}_{ab}\big(\partial^{(x)}_{cd}V_y\big)(h)\;\Psi(h)\Big)Tx,y​=a,b,c,d=0∑2​Gabcd​(h(x))(2∂ab(x)​Vy​(h)∂cd(x)​Ψ(h)+∂ab(x)​(∂cd(x)​Vy​)(h)Ψ(h))

and Ty,xT_{y,x}Ty,x​ is the same expression with xxx and yyy swapped everywhere: the supermetric is G(h(y))G(h(y))G(h(y)), the derivatives are taken at site yyy, the potential is VxV_xVx​, and Ψ\PsiΨ is the same function.

The function HyΨ\mathcal{H}_y\PsiHy​Ψ is the whole map h′↦(HyΨ)(h′)h'\mapsto(\mathcal{H}_y\Psi)(h')h′↦(Hy​Ψ)(h′) on C\mathcal{C}C, and Hx\mathcal{H}_xHx​ is applied to it. On the left this needs Fréchet derivatives of HyΨ\mathcal{H}_y\PsiHy​Ψ at points near hhh, not only at hhh. The identity is claimed at this single point hhh only. In the first-order term the index pair (a,b)(a,b)(a,b) goes on the derivative of VyV_yVy​ and (c,d)(c,d)(c,d) on the derivative of Ψ\PsiΨ, and the sum is not symmetrized over swapping the pairs. The left side contains Λ\LambdaΛ only through the operators, and the right side only through Vx,VyV_x,V_yVx​,Vy​.

Edge cases and points to note.

  • R(⋅,z)R(\cdot,z)R(⋅,z) is assumed C2C^2C2 for every site, globally on C\mathcal{C}C. No locality is assumed: R(h,y)R(h,y)R(h,y) may depend on h(x)h(x)h(x). This is why ∂(x)Vy\partial^{(x)}V_y∂(x)Vy​ can be nonzero.
  • Ψ\PsiΨ is assumed C4C^4C4 on all of C\mathcal{C}C. Nothing restricts Ψ\PsiΨ or RRR to physical configurations. Only the evaluation point hhh is physical, and nearby non-symmetric configurations still lie in the domain of the derivatives.
  • det⁡\sqrt{\det}det​ and GabcdG_{abcd}Gabcd​ are defined by junk conventions where det⁡≤0\det\le 0det≤0 (square root 000, division by 000 gives 000). Because det⁡h(z)>0\det h(z)>0deth(z)>0 at the physical hhh, they are given by the displayed formulas on a neighbourhood of hhh.
  • The non-vanishing hypotheses κ≠0\kappa\neq 0κ=0 and ℏ≠0\hbar\neq 0ℏ=0 exclude the junk value (2⋅0)−1=0(2\cdot 0)^{-1}=0(2⋅0)−1=0 in VzV_zVz​. If ℏ=0\hbar=0ℏ=0, both sides would plainly be 000, so that case is excluded but would not be problematic.
  • The hypotheses are jointly satisfiable, for example with XXX of two points, h(x)=h(y)=Ih(x)=h(y)=Ih(x)=h(y)=I, R≡0R\equiv 0R≡0 and Ψ≡0\Psi\equiv 0Ψ≡0, so the theorem is not vacuous whenever ∣X∣≥2|X|\ge 2∣X∣≥2.
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