Skip to content

Canonical Constraints and Total Hamiltonian

Scrollable formula. Use Left and Right arrow keys to reveal hidden content.

Scrollable table. Use Left and Right arrow keys to reveal additional columns.

Scrollable figure. Use Left and Right arrow keys to inspect hidden content.

The Legendre transform of the ADM action does not produce an unconstrained Hamiltonian system. Lapse and shift multiply, respectively, a scalar and a spatial-vector constraint, while their own conjugate momenta vanish. The resulting total Hamiltonian separates bulk gauge evolution from boundary generators and makes the two local gravitational degrees of freedom visible.

Required background. Presymplectic Potentials and Currents supplies the canonical meaning of the action’s boundary variation, while The Gibbons-Hawking-York Boundary Term supplies the differentiable spacetime action from which the split begins.

Helpful background. The Einstein Equation supplies the covariant constraints being recast, Hypersurface Normals and Tangential Projectors supplies normal and tangential projections, and Formulations of General Relativity locates the Hamiltonian system among equivalent formulations.

Legendre transform of the gravitational action

Section titled “Legendre transform of the gravitational action”

The velocity of the spatial metric is

γ˙ij=−2αKij+Diβj+Djβi,\dot\gamma_{ij}=-2\alpha K_{ij}+D_i\beta_j+D_j\beta_i,

and the momentum appropriate to this sign of KijK_{ij} is

πij=γ16πG(Kγij−Kij).\pi^{ij}=\frac{\sqrt\gamma}{16\pi G} \left(K\gamma^{ij}-K^{ij}\right).

Consequently,

−2απijKij=αγ8πG(KijKij−K2).-2\alpha\pi^{ij}K_{ij} =\frac{\alpha\sqrt\gamma}{8\pi G} \left(K_{ij}K^{ij}-K^2\right).

The shift term can be integrated by parts:

2πijDiβj=−2βiDjπji+2Di(πijβj).2\pi^{ij}D_i\beta_j =-2\beta^iD_j\pi^j{}_i +2D_i(\pi^{ij}\beta_j).

Subtracting the ADM Lagrangian and collecting the divergence into a boundary term gives

Hgrav=∫Σd3x(αC⊥g+βiCig)+B,H_{\mathrm{grav}} =\int_\Sigma d^3x \left(\alpha\mathcal C_\perp^{\rm g} +\beta^i\mathcal C_i^{\rm g}\right)+B,

where

C⊥g=16πGγ(πijπij−12π2)−γ16πG((3)R−2Λ),\mathcal C_\perp^{\rm g} =\frac{16\pi G}{\sqrt\gamma} \left(\pi_{ij}\pi^{ij}-\frac12\pi^2\right) -\frac{\sqrt\gamma}{16\pi G} \left({}^{(3)}R-2\Lambda\right),

and

Cig=−2Djπji.\mathcal C_i^{\rm g}=-2D_j\pi^j{}_i.

This Legendre transform and its surface contribution are developed in Arnowitt, Deser, and Misner 1962/2004, §§3–5, PDF; the same sign-consistent construction is given in Gourgoulhon 2007, §4.5.2, pp. 67–69, PDF.

Define the projections of stress energy by

ρ=Tμνnμnν,ji=−γiμnνTμν.\rho=T_{\mu\nu}n^\mu n^\nu, \qquad j_i=-\gamma_i{}^\mu n^\nu T_{\mu\nu}.

For minimally coupled matter with no metric-velocity coupling, its Legendre transform adds the corresponding matter canonical terms and changes the gravitational constraints to

C⊥=C⊥g+γ ρ,Ci=−2Djπji−γ ji.\mathcal C_\perp =\mathcal C_\perp^{\rm g}+\sqrt\gamma\,\rho, \qquad \mathcal C_i =-2D_j\pi^j{}_i-\sqrt\gamma\,j_i.

Varying lapse and shift gives

C⊥=0,Ci=0.\mathcal C_\perp=0, \qquad \mathcal C_i=0.

In geometric variables these are exactly

(3)R+K2−KijKij−2Λ=16πGρ,{}^{(3)}R+K^2-K_{ij}K^{ij}-2\Lambda=16\pi G\rho, DjKji−DiK=8πGji.D_jK^j{}_i-D_iK=8\pi G j_i.

The signs in the vector constraint are linked: with ji=−γiμnνTμνj_i=-\gamma_i{}^\mu n^\nu T_{\mu\nu}, the identity −2Djπji=γji-2D_j\pi^j{}_i=\sqrt\gamma j_i holds on the momentum constraint.

A real scalar field supplies a useful independent sign check. For Lm=−−g[(∇ϕ)2/2+V(ϕ)]\mathcal L_{\rm m}=-\sqrt{-g}[(\nabla\phi)^2/2+V(\phi)], define v=nμ∇μϕ=(ϕ˙−βiDiϕ)/αv=n^\mu\nabla_\mu\phi=(\dot\phi-\beta^iD_i\phi)/\alpha. Then

pϕ=γ v,ρ=12v2+12DiϕDiϕ+V,ji=−vDiϕ.p_\phi=\sqrt\gamma\,v, \qquad \rho=\frac12v^2+\frac12D_i\phi D^i\phi+V, \qquad j_i=-vD_i\phi.

The matter Legendre transform is therefore

pϕϕ˙−Lm=αγρ+βipϕDiϕ=αγρ−βiγji.p_\phi\dot\phi-\mathcal L_{\rm m} =\alpha\sqrt\gamma\rho+\beta^i p_\phi D_i\phi =\alpha\sqrt\gamma\rho-\beta^i\sqrt\gamma j_i.

The minus sign multiplying jij_i is thus required by the scalar’s ordinary canonical momentum, not a discretionary choice. Gauge fields can also have their own Gauss constraints; the full total Hamiltonian then includes those additional multipliers.

Primary constraints and the total Hamiltonian

Section titled “Primary constraints and the total Hamiltonian”

Because neither α˙\dot\alpha nor β˙i\dot\beta^i occurs in the action,

pα=0,pi=0p_\alpha=0, \qquad p_i=0

are primary constraints. The Dirac total Hamiltonian is

HT=∫Σd3x(αC⊥+βiCi+upα+uipi)+B.H_T=\int_\Sigma d^3x \left( \alpha\mathcal C_\perp+\beta^i\mathcal C_i +u p_\alpha+u^ip_i \right)+B.

The arbitrary multipliers uu and uiu^i determine the velocities of lapse and shift. Preserving their primary constraints in time produces the scalar and vector constraints. The metric Hamilton equation begins

γ˙ij=32πGαγ(πij−12πγij)+Lβγij,\dot\gamma_{ij} =\frac{32\pi G\alpha}{\sqrt\gamma} \left(\pi_{ij}-\frac12\pi\gamma_{ij}\right) +\mathcal L_\beta\gamma_{ij},

which returns γ˙ij=−2αKij+Lβγij\dot\gamma_{ij}=-2\alpha K_{ij}+\mathcal L_\beta\gamma_{ij}. The momentum equation is equivalent to the spatial projection of the Einstein equation.

At a regular point of the pure-gravity constraint surface on a closed slice, the twelve canonical fields (γij,πij)(\gamma_{ij},\pi^{ij}) are reduced by four first-class constraints and four gauge directions, leaving four phase-space functions, or two configuration degrees of freedom, per point. On a region with boundary, BB is not optional bookkeeping: it is chosen so that HTH_T has well-defined functional derivatives under the stated boundary conditions.

Inertial Minkowski data have γij=δij\gamma_{ij}=\delta_{ij}, πij=0\pi^{ij}=0, and (3)R=0{}^{(3)}R=0, so both vacuum constraints vanish. For Milne slices,

(3)R=−6τ2,πij=−γ8πGτγij.{}^{(3)}R=-\frac6{\tau^2}, \qquad \pi^{ij}=-\frac{\sqrt\gamma}{8\pi G\tau}\gamma^{ij}.

The momentum constraint vanishes by spatial homogeneity. The curvature and momentum terms in C⊥\mathcal C_\perp cancel, confirming that a nonzero canonical momentum can describe flat spacetime.

Exercise: verify the Milne scalar constraint

Section titled “Exercise: verify the Milne scalar constraint”

Using the Milne data above with Λ=0\Lambda=0 and no matter, compute πijπij−12π2\pi_{ij}\pi^{ij}-\tfrac12\pi^2 and verify C⊥=0\mathcal C_\perp=0.

Solution

Let A=γ/(8πGτ)A=\sqrt\gamma/(8\pi G\tau), so πij=−Aγij\pi^{ij}=-A\gamma^{ij} and π=−3A\pi=-3A. Then

πijπij−12π2=3A2−92A2=−32A2.\pi_{ij}\pi^{ij}-\frac12\pi^2 =3A^2-\frac92A^2=-\frac32A^2.

The kinetic contribution is therefore

16πGγ(−32A2)=−3γ8πGτ2.\frac{16\pi G}{\sqrt\gamma} \left(-\frac32A^2\right) =-\frac{3\sqrt\gamma}{8\pi G\tau^2}.

The curvature contribution is

−γ16πG(−6τ2)=3γ8πGτ2.-\frac{\sqrt\gamma}{16\pi G} \left(-\frac6{\tau^2}\right) =\frac{3\sqrt\gamma}{8\pi G\tau^2}.

Their sum vanishes. Milne data satisfy the vacuum scalar constraint even though neither term vanishes separately.

  • Arnowitt, Richard, Stanley Deser, and Charles W. Misner. “The Dynamics of General Relativity.” In Gravitation: An Introduction to Current Research, edited by Louis Witten, 227–265. Wiley, 1962; reprinted as arXiv:gr-qc/0405109 (2004). arXiv.
  • Gourgoulhon, Éric. 3+1 Formalism and Bases of Numerical Relativity. Lecture notes, arXiv:gr-qc/0703035v1 (2007). arXiv.