Skip to content

Covariant Phase Space and Horizon Hamiltonians

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.

Covariant phase space turns a spacetime symmetry into a charge variation without choosing a time coordinate or canonical decomposition. The key object is not a charge density by itself but the boundary form δQξ−ξ⋅θ\delta Q_\xi-\xi\mathbin{\cdot}\theta. Its integral is the variation of a Hamiltonian only when the boundary conditions make it finite, conserved, and integrable on the chosen solution space. At a bifurcation surface the horizon generator vanishes, exposing the gravitational area term that enters the stationary first law. For Maxwell fields, the boundary on which the electric work term appears depends on the gauge representative and its regularity.

Required background. Variational and Symplectic Geometry supplies field-space variations, symplectic currents, and gauge degeneracies. Helpful background. Exact Solutions of Einstein’s Equations supplies stationary test families; Causality and Global Structure supplies the boundaries of a hypersurface; Einstein Equations fixes the equations of motion; Gravitational Actions and Boundary Terms fixes the action variation; Black-Hole Geometry and Mechanics supplies the normalized horizon generator.

From a Lagrangian to the symplectic current

Section titled “From a Lagrangian to the symplectic current”

Let Ψ\Psi denote all dynamical fields and let L[Ψ]L[\Psi] be the covariant four-form Lagrangian. Its first variation can be written

δL=E(Ψ) δΨ+dθ(Ψ;δΨ),\delta L=E(\Psi)\,\delta\Psi+d\theta(\Psi;\delta\Psi),

where E=0E=0 denotes the field equations and θ\theta is the presymplectic potential three-form. For two commuting variations, define

ω(Ψ;δ1Ψ,δ2Ψ)=δ1θ(Ψ;δ2Ψ)−δ2θ(Ψ;δ1Ψ).\omega(\Psi;\delta_1\Psi,\delta_2\Psi) =\delta_1\theta(\Psi;\delta_2\Psi) -\delta_2\theta(\Psi;\delta_1\Psi).

Its integral over a Cauchy hypersurface is the presymplectic form. “Pre” is important: diffeomorphisms that vanish appropriately at the boundary are degenerate gauge directions, while transformations with nontrivial boundary action can carry charges.

For the Einstein–Hilbert Lagrangian, the metric potential can be represented in components by

θabc(g;δg)=116πGϵdabc(∇dδg−∇eδgde),\theta_{abc}(g;\delta g) =\frac{1}{16\pi G}\epsilon_{dabc} \left(\nabla^d\delta g-\nabla_e\delta g^{de}\right),

where δg=gabδgab\delta g=g_{ab}\delta g^{ab} because the independent variation here is δgab\delta g^{ab}. If instead hab=δgabh_{ab}=\delta g_{ab} is the independent variation, then δgab=−hab\delta g^{ab}=-h^{ab} and the same potential reads ϵdabc(∇ehde−∇dh)/(16πG)\epsilon_{dabc}(\nabla_eh^{de}-\nabla^dh)/(16\pi G). The two familiar formulas therefore agree after translating the varied metric variable.

For a vector field ξa\xi^a, define the diffeomorphism Noether current

Jξ=θ(Ψ;LξΨ)−ξ⋅L.J_\xi=\theta(\Psi;\mathcal L_\xi\Psi)-\xi\mathbin{\cdot}L.

Covariance implies that, on shell, JξJ_\xi is closed and can locally be written Jξ=dQξJ_\xi=dQ_\xi. For Einstein gravity the Noether charge two-form is

Qξ=−116πG⋆dξ,Q_\xi=-\frac{1}{16\pi G}\star d\xi,

which is one half of the time-translation Komar normalization. The missing half is not an error: the Hamiltonian variation also contains −ξ⋅θ-\xi\mathbin{\cdot}\theta, and the relation between Noether and ADM charges uses both terms.

When the background solves the equations, the perturbation solves the linearized equations, and ξ\xi is field independent during the variation, set

kξ(δΨ)=δQξ−ξ⋅θ(Ψ;δΨ).k_\xi(\delta\Psi) =\delta Q_\xi-\xi\mathbin{\cdot}\theta(\Psi;\delta\Psi).

The fundamental identity is

ΩΣ(δΨ,LξΨ)=∫∂Σkξ.\Omega_\Sigma(\delta\Psi,\mathcal L_\xi\Psi) =\int_{\partial\Sigma}k_\xi.

On a selected closed boundary component SS, the candidate charge one-form is

δHξ[S]=∫Skξ.\delta H_\xi[S]=\int_S k_\xi.

If ξ\xi is a symmetry of the background, the left side of the full-boundary identity vanishes. For ∂Σ=S∞∪(−B)\partial\Sigma=S_\infty\cup(-B), the orientations therefore give

0=∫∂Σkξ=∫S∞kξ−∫Bkξ,0=\int_{\partial\Sigma}k_\xi =\int_{S_\infty}k_\xi-\int_B k_\xi,

so the same charge variation can be evaluated at infinity or at the horizon. Iyer and Wald define the symplectic current and Noether charge in §§III–IV, derive the closed boundary form in Iyer and Wald 1994, §VI, equations (73)–(77), Open PDF, p. 14, and construct the Hamiltonian in Iyer and Wald 1994, §VI, equations (78)–(82), Open PDF, p. 15.

The displayed boundary integral is initially a one-form on solution space. It equals the variation of a function HξH_\xi only if its integral around closed loops vanishes. For field-independent ξ\xi, the local obstruction is controlled by the symplectic flux:

∫Sξ⋅ω(Ψ;δ1Ψ,δ2Ψ)=0\int_S \xi\mathbin{\cdot}\omega (\Psi;\delta_1\Psi,\delta_2\Psi)=0

for all allowed tangent variations on the selected boundary component SS, up to the boundary prescription used. Nonzero radiation flux, varying sources, or insufficient falloff can prevent integrability or conservation. One may then obtain a flux-balance law rather than a time-independent Hamiltonian charge.

The finite-inner-boundary version of this question is developed on Isolated-Horizon First Laws, where exactness of the horizon one-form selects the permitted evolution fields and energies.

The vector field may itself depend on the solution, as happens when χ=t+ΩHϕ\chi=t+\Omega_H\phi is compared across a family with changing ΩH\Omega_H. Then one must specify whether ξ\xi or its coefficients are held fixed in the variation and include the corresponding correction. Silently varying the generator double-counts work terms.

There are also standard ambiguities: L↦L+dμL\mapsto L+d\mu, θ↦θ+δμ+dY\theta\mapsto\theta+\delta\mu+dY, and a corresponding shift of QξQ_\xi. Boundary conditions, counterterms, reference solutions, and the desired Hamiltonian select a finite representative. Exact-form ambiguity does not license arbitrary physical charges; consistent changes must be made at every boundary and tested against the variational identity.

Take a stationary black hole with horizon generator χ=t+ΩHϕ\chi=t+\Omega_H\phi and a smooth bifurcation surface BB. The defining property χa∣B=0\chi^a|_B=0 eliminates χ⋅θ\chi\mathbin{\cdot}\theta there. Moreover,

∇aχb=κϵabon B,\nabla_a\chi_b=\kappa\epsilon_{ab} \quad\text{on }B,

so the Einstein Noether charge gives, under variations between nearby stationary solutions with the comparison generator handled consistently,

∫BδQχ=κ8πGδAH.\int_B\delta Q_\chi =\frac{\kappa}{8\pi G}\delta A_H.

At infinity, the same Hamiltonian variation is δM−ΩHδJ\delta M-\Omega_H\delta J in vacuum. Equality of the two boundary evaluations yields the vacuum stationary first law. The full smooth-field bifurcation-surface argument and its assumptions are given in Iyer and Wald 1994, §VI, Theorem 6.1 and equations (92)–(102), Open PDF, pp. 17–18.

For Einstein–Maxwell theory, define the potential difference, which is invariant under admissible symmetry-compatible gauge transformations,

ΦH=−χaAa∣H+χaAa∣∞.\Phi_H=-\chi^aA_a\big|_H+\chi^aA_a\big|_\infty.

There are two consistent ways to allocate its work term between the boundaries.

Maxwell representativePermitted inner boundaryInfinity boundary formInner boundary form
Aa,δAa=O(r−1)A_a,\delta A_a=O(r^{-1}) at infinity; only their pullbacks are smooth on the future horizonA suitable future-horizon cross-section SHS_H; when BB exists, choose SHS_H strictly to its futureδM−ΩHδJ\delta M-\Omega_H\delta JκδAH/(8πG)+ΦHδQ\kappa\delta A_H/(8\pi G)+\Phi_H\delta Q
AaA_a and δAa\delta A_a smooth through BB, with stationary gauge dataThe bifurcation surface BBδM−ΩHδJ−ΦHδQ\delta M-\Omega_H\delta J-\Phi_H\delta QκδAH/(8πG)\kappa\delta A_H/(8\pi G)

The first row and second row are two complete boundary prescriptions. Mixing the vanishing-at-infinity outer expression from the first with the smooth-BB inner expression from the second would incorrectly lose the electric work term. Gao contrasts the smooth-bifurcation allocation with the future-horizon route in Gao 2003, §I, equations (2)–(8), and §III, equations (27)–(50), Open PDF.

The detailed Reissner–Nordström gauge check, the full boundary evaluation, and the extremal alternative are carried out on First Law for Stationary Einstein–Maxwell Black Holes. The structural point here is that boundary regularity and gauge reference determine where a work term is displayed, while the complete Hamiltonian identity is unchanged.

A proposed horizon Hamiltonian should pass four tests. Its variation must be finite under the declared falloff. Symplectic flux through excluded boundary segments must vanish or be included as flux. Mixed second variations must commute on the permitted solution family. Finally, the infinity and horizon expressions must use the same tat^a, axial periodicity, electromagnetic gauge representative and reference, admissible inner cross-section, and orientation.

For Schwarzschild, δHt=δM\delta H_t=\delta M at infinity and δHt=κδAH/(8πG)\delta H_t=\kappa\delta A_H/(8\pi G) at the bifurcation sphere. Since AH=16πG2M2A_H=16\pi G^2M^2 and κ=1/(4GM)\kappa=1/(4GM), the horizon expression equals δM\delta M. This check tests the Noether prefactor and the absence of an unaccounted boundary term.

Suppose a proposed two-parameter boundary charge has variation

δH=(2x+y) δx+x δy.\delta H=(2x+y)\,\delta x+x\,\delta y.

Is it integrable? Find HH up to a constant. What would change if the coefficient of δy\delta y were 2x2x?

Solution

For the displayed one-form, ∂y(2x+y)=1\partial_y(2x+y)=1 and ∂x(x)=1\partial_x(x)=1, so mixed derivatives agree. An integral is

H=x2+xy+C.H=x^2+xy+C.

Indeed,

dH=(2x+y) dx+x dy,dH=(2x+y)\,dx+x\,dy,

which reproduces both coefficients of the proposed variation.

If the second coefficient were 2x2x, the mixed derivatives would be 11 and 22. The one-form would not be closed, so no path-independent Hamiltonian would exist on the full (x,y)(x,y) space. One would need to restrict the allowed variations, add justified boundary data, or interpret the mismatch as flux; calling the original expression δH\delta H would not make it integrable.

First Law for Stationary Einstein–Maxwell Black Holes now evaluates the two boundaries and checks the result directly on Kerr–Newman.

  • Gao, Sijie. “First Law of Black Hole Mechanics in Einstein–Maxwell and Einstein–Yang–Mills Theories.” Physical Review D 68 (2003), 044016. DOI. Open PDF, arXiv:gr-qc/0304094v2.
  • Iyer, Vivek, and Robert M. Wald. “Some Properties of Noether Charge and a Proposal for Dynamical Black Hole Entropy.” Physical Review D 50 (1994): 846–864. DOI.