Skip to content

Diffeomorphism Constraints and Gauge Orbits

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.

In canonical gravity, a spacetime diffeomorphism is represented by constraints only after its behavior at the boundary has been fixed. Deformations that die off at the boundary generate gauge orbits; deformations that act on admissible boundary data can instead carry energy, momentum, or angular momentum. The same local formula can therefore represent either redundancy or physical symmetry, and the differentiability of its generator tells the two cases apart.

Required background. Variational and Symplectic Geometry supplies Hamiltonian generators and surface terms. Equivalence Principles and Gravitational Measurement supplies the operational meaning of coordinates and reference frames.

Helpful background. Flat-Spacetime Quantum Field Theory supplies constrained quantization, Einstein Equations supplies the covariant dynamics, Quantum Information supplies subsystem language, and Quantum Gravity sets this volume’s regimes and evidence standards.

Take a globally hyperbolic four-dimensional region foliated by spacelike hypersurfaces Σt\Sigma_t, with induced metric hijh_{ij}, canonical momentum πij\pi^{ij}, lapse NN, and shift NiN^i. Matter terms may be added without changing the logic. Up to the boundary action, the Hamiltonian form is

S=∫dt[∫Σtd3x(πijh˙ij−NH⊥−NiHi)−H∂Σ].S=\int dt\left[\int_{\Sigma_t}d^3x \left(\pi^{ij}\dot h_{ij}-N\mathcal H_\perp-N^i\mathcal H_i\right)-H_{\partial\Sigma}\right].

The lapse and shift are Lagrange multipliers, giving the Hamiltonian and momentum constraints H⊥≈0\mathcal H_\perp\approx0 and Hi≈0\mathcal H_i\approx0. Choose Kij=(h˙ij−DiNj−DjNi)/(2N)K_{ij}=(\dot h_{ij}-D_iN_j-D_jN_i)/(2N) and the momentum density πij=h (Kij−hijK)/(16πG)\pi^{ij}=\sqrt h\,(K^{ij}-h^{ij}K)/(16\pi G), where h=det⁡hijh=\det h_{ij}. Then

Hi=−2Djπji+Himatter,\mathcal H_i=-2D_j\pi^j{}_i+\mathcal H_i^{\mathrm{matter}},

while H⊥\mathcal H_\perp combines the spatial scalar curvature and a quadratic expression in πij\pi^{ij}. The complete canonical decomposition and its relation to the Einstein equation are derived in Arnowitt, Deser, and Misner 1962/2008, §§3–4.

On a closed hypersurface, the smeared momentum constraint

D[ξ]=∫Σd3x ξiHiD[\xi]=\int_\Sigma d^3x\,\xi^i\mathcal H_i

generates spatial diffeomorphisms:

{hij,D[ξ]}=Lξhij,{πij,D[ξ]}=Lξπij,\{h_{ij},D[\xi]\}=\mathcal L_\xi h_{ij}, \qquad \{\pi^{ij},D[\xi]\}=\mathcal L_\xi\pi^{ij},

where πij\pi^{ij} is a tensor density. The constraints are first class; their brackets close with structure functions. This closure expresses hypersurface-deformation consistency, not an ordinary Lie algebra with constant structure coefficients.

The surface term decides the transformation

Section titled “The surface term decides the transformation”

For clarity, omit matter and write the ordinary tensor pij=πij/hp^{ij}=\pi^{ij}/\sqrt h. Integration by parts gives

∫Σd3x πijLξhij=∫Σd3x ξiHi+2∫∂Σd2x σ nipijξj.\int_\Sigma d^3x\,\pi^{ij}\mathcal L_\xi h_{ij} =\int_\Sigma d^3x\,\xi^i\mathcal H_i +2\int_{\partial\Sigma}d^2x\,\sqrt\sigma\,n_ip^{ij}\xi_j.

Here nin_i points outward within Σ\Sigma and σ\sigma is the determinant of the induced two-metric. The surface measure already carries its density, so the integrand uses pijp^{ij}, not the weight-one bulk density πij\pi^{ij}. This is the covariant divergence theorem applied to h pij\sqrt h\,p^{ij}. If ξi\xi^i vanishes at the boundary, the surface term vanishes and D[ξ]D[\xi] is a sum of constraints. On the constraint surface its Hamiltonian flow moves along a gauge orbit.

If ξi\xi^i approaches an allowed nonzero boundary transformation, the bare bulk constraint is generally not functionally differentiable. One must add a surface functional Q[ξ]Q[\xi] so that

G[ξ]=∫Σd3x ξiHi+Q[ξ]G[\xi]=\int_\Sigma d^3x\,\xi^i\mathcal H_i+Q[\xi]

has a well-defined variation. On shell, G[ξ]=Q[ξ]G[\xi]=Q[\xi] need not vanish. In asymptotically flat gravity, the corresponding lapse and shift limits yield ADM energy, momentum, and angular momentum. Regge and Teitelboim give the systematic differentiability analysis and asymptotic generators Regge and Teitelboim 1974, §§3–6.

Thus “diffeomorphisms are gauge” is too coarse. Proper diffeomorphisms preserve the boundary data and have vanishing generators; allowed boundary symmetries can act nontrivially on the physical boundary sector. A charge vanishing on one particular state does not make its generator gauge on all states. Nor is “large” synonymous with “charged”: large diffeomorphisms can be disconnected from the identity, a separate global issue. The classification depends on falloff conditions, topology, corners, and the chosen ensemble.

Let ΓC\Gamma_C be the constraint surface in phase space. The reduced phase space is schematically

Γred=ΓC/G0,\Gamma_{\mathrm{red}}=\Gamma_C/\mathcal G_0,

where G0\mathcal G_0 is the group generated by proper constraints, not necessarily every transformation preserving the form of the boundary conditions. A Dirac observable OO satisfies {O,CA}≈0\{O,C_A\}\approx0 for all proper first-class constraints CAC_A. A boundary charge may act nontrivially on OO and organize it into a physical symmetry representation.

There are two useful checks. On compact Σ\Sigma without boundary, Q[ξ]Q[\xi] is absent and the spatial generator reduces to the constraint. In the weak-field asymptotically flat limit, a constant asymptotic time translation instead gives the nonzero ADM Hamiltonian. Treating both transformations as gauge would incorrectly set the total energy to zero.

The quotient can also be singular where gauge orbits change dimension, and global gauge fixing can fail. The local first-class analysis therefore does not guarantee a globally smooth reduced phase space. Quantum constructions may use Dirac quantization instead, imposing constraints on states while representing boundary charges as operators.

Starting from Hi=−2Djπji\mathcal H_i=-2D_j\pi^j{}_i, integration by parts gives

∫ΣξiHi=2∫ΣπijDiξj−2∫∂Σσ nipijξj.\int_\Sigma \xi^i\mathcal H_i =2\int_\Sigma \pi^{ij}D_i\xi_j -2\int_{\partial\Sigma}\sqrt\sigma\,n_ip^{ij}\xi_j.

Because πij\pi^{ij} is symmetric, the volume term is ∫ΣπijLξhij\int_\Sigma\pi^{ij}\mathcal L_\xi h_{ij}. Moving the surface term to the other side reproduces the earlier identity and checks the sign of Q[ξ]Q[\xi]. At a fixed finite boundary, take ξ\xi field independent, tangent to the boundary, and preserving the allowed boundary data; these conditions eliminate the residual normal boundary variation and make the spatial canonical generator differentiable. A normal deformation requires a separately specified boundary variational problem; integration by parts alone does not establish integrability of its charge. Reversing the boundary orientation reverses the charge integral. Field-dependent ξ\xi also requires additional terms in the charge algebra.

Constraint does not mean zero Hamiltonian in every setting. The bulk Hamiltonian is constrained, but a boundary Hamiltonian can generate physical time translations and be nonzero.

Coordinate covariance is not the full gauge test. An expression may be a scalar under passive coordinate changes yet still depend on where a dynamical point is placed along a gauge orbit. The relevant test is its Poisson bracket with proper constraints.

Exercise. Let Σ\Sigma have a fixed boundary, ignore matter, and take ξ\xi field independent and tangent to the boundary with boundary-preserving data. Starting from D0[ξ]=∫ΣξiHiD_0[\xi]=\int_\Sigma\xi^i\mathcal H_i, derive the boundary term needed so that the generator equals ∫ΣπijLξhij\int_\Sigma\pi^{ij}\mathcal L_\xi h_{ij}. State when the resulting generator is pure gauge.

Solution

Using Hi=−2Djπji\mathcal H_i=-2D_j\pi^j{}_i and integrating by parts gives

D0[ξ]=∫Σd3x πijLξhij−2∫∂Σd2x σ nipijξj.D_0[\xi]=\int_\Sigma d^3x\,\pi^{ij}\mathcal L_\xi h_{ij} -2\int_{\partial\Sigma}d^2x\,\sqrt\sigma\,n_ip^{ij}\xi_j.

Therefore the differentiable generator with the desired canonical action is

D[ξ]=D0[ξ]+2∫∂Σd2x σ nipijξj.D[\xi]=D_0[\xi]+2\int_{\partial\Sigma}d^2x\,\sqrt\sigma\,n_ip^{ij}\xi_j.

If ξi\xi^i obeys the proper-gauge boundary conditions so the surface functional vanishes throughout the allowed phase space, D[ξ]≈0D[\xi]\approx0 and generates a gauge orbit. If the allowed boundary symmetry instead has a nontrivial charge functional, it must remain in the physical symmetry algebra. Its numerical value may vanish on a special state without changing that classification.

  • 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. New York: Wiley, 1962. Republication, arXiv:gr-qc/0405109.
  • Regge, Tullio, and Claudio Teitelboim. “Role of Surface Integrals in the Hamiltonian Formulation of General Relativity.” Annals of Physics 88 (1974): 286–318. DOI.