Covariant Phase-Space Method

• Covariant Phase-Space Method is a geometric framework for constructing the symplectic structure of Lagrangian field theories while preserving manifest covariance.
• It systematically computes presymplectic potentials and currents to derive conserved charges and manage boundary contributions in gravity and gauge theories.
• The method employs dual exterior calculi on spacetime and field space to analyze symmetries, integrability, and flux in a unified, boundary-aware formulation.

The covariant phase-space method is a manifestly covariant construction of the symplectic structure and conserved charges of a Lagrangian field theory from the space of its solutions. In its standard form, one starts from a spacetime manifold $\mathcal M$, a Lagrangian $d$-form $L[\Phi]$, and the decomposition

$$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$

where $E_I=0$ are the equations of motion and $\Theta$ is the presymplectic potential. The presymplectic current is the field-space exterior derivative of $\Theta$,

$$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$

and its integral over a Cauchy slice $\Sigma$ defines the presymplectic form

$$\Omega[\delta_1\Phi,\delta_2\Phi;\Phi]=\int_\Sigma \omega.$$

On shell, $d$0, so $d$1 is independent of the choice of $d$2 after quotienting gauge directions. In recent treatments, the formalism is described as acting simultaneously on spacetime and on the infinite-dimensional solution phase space, with parallel geometric structures on both manifolds (Varo, 2023, Golshani et al., 15 Mar 2026).

1. Definition and geometric setting

The covariant phase space is the space of all solutions of the Euler–Lagrange equations. In the notation of the gravity-with-boundaries literature,

$d$3

while in the geometric solution-space formulation it is the Solution Phase Space (SPS) $d$4, whose points are solutions $d$5 modulo proper gauge transformations (Varo, 2023, Golshani et al., 15 Mar 2026).

A central structural point is that the method uses two exterior calculi. On spacetime one has $d$6, interior products $d$7, and Lie derivatives $d$8; on field space one has $d$9, interior products $L[\Phi]$0, and Lie derivatives $L[\Phi]$1. The “fundamental link” is that a spacetime symmetry $L[\Phi]$2 induces a field-space vector $L[\Phi]$3 via

$L[\Phi]$4

This parallelism is not merely notational. It underlies the construction of surface charges, the analysis of integrability, and later developments involving connections, torsion, and curvature on solution space (Golshani et al., 15 Mar 2026).

In the local jet-bundle picture used for theories with boundaries, all bulk and boundary forms live in a bicomplex graded by spacetime degree and field-variation degree. The familiar first-variation formula then becomes the starting point of a bicomplex descent. The covariant phase-space two-form is obtained by pulling the presymplectic current back to the solution space and integrating it over a Cauchy slice. This is the sense in which the method replaces a noncovariant $L[\Phi]$5 split by a construction directly on the solution space (Margalef-Bentabol et al., 2020).

2. Boundaries, corners, and the relative bicomplex

For manifolds with timelike boundaries, the boundary-free formalism is insufficient. One must add boundary terms to obtain a well-posed variational principle, and one must track boundary and corner contributions in the presymplectic structure. The relative bicomplex framework does this by replacing an ordinary Lagrangian form with a Lagrangian pair

$L[\Phi]$6

or equivalently $L[\Phi]$7, together with a relative differential $L[\Phi]$8 combining the bulk exterior derivative with pullbacks to the boundary (Margalef-Bentabol et al., 2020, Varo, 2023).

In this framework the generalized first-variation formula takes the form

$L[\Phi]$9

where $$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$0 are the bulk Euler–Lagrange forms, $$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$1 are the boundary “equations,” and $$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$2 is the relative presymplectic potential. The relative presymplectic current is then

$$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$3

and the presymplectic form on a Cauchy pair $$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$4 is

$$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$5

In the formulation of (Varo, 2023), the same structure appears as

$$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$6

Three general statements organize the boundary theory. First, any local Lagrangian pair admits a well-defined relative presymplectic potential. Second, the relative presymplectic current is $$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$7-closed on shell, so $$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$8 is independent of $$\delta L = E_I[\Phi]\,\delta\phi^I + d\,\Theta[\delta\Phi;\Phi],$$9. Third, infinitesimal gauge symmetries generate conserved charges that depend only on boundary data (Varo, 2023).

A distinctive claim of the relative bicomplex literature is that it gives a formal equivalence between the relative version of a theory with boundary and the non-relative version of the same theory with no boundary. This is used to transfer bicomplex identities, Cartan formulae, and exactness theorems from the boundaryless setting to the theory with boundary and corners. This also underlies the claim that the framework is “clean and ambiguity-free” for the solution space and its symplectic structure (Margalef-Bentabol et al., 2020, G. et al., 2021).

3. Symmetries, Noether currents, and Hamiltonian charges

For a vertical symmetry $E_I=0$0 satisfying

$E_I=0$1

the relative Noether current pair is

$E_I=0$2

On shell, $E_I=0$3, and the associated charge on a Cauchy pair is

$E_I=0$4

For diffeomorphisms generated by $E_I=0$5, one obtains on shell a Noether current $E_I=0$6, so the Hamiltonian generator is a codimension-two boundary integral,

$E_I=0$7

a statement emphasized in the relative-bicomplex treatment of gravity with timelike boundaries (Varo, 2023).

The boundary literature also gives an algorithmic derivation of the Hamiltonian in the presence of spatial boundaries. Starting from an action

$E_I=0$8

one demands that, at solutions of the equations of motion, the action is stationary modulo future/past boundary terms under arbitrary variations obeying the spatial boundary conditions. From this, both the symplectic structure and the Hamiltonian for any diffeomorphism that preserves the theory are constructed, and the term denoted $E_I=0$9 in the original literature emerges naturally. The same procedure produces an additional boundary term $\Theta$0, which can already be nonvanishing in general relativity with sufficiently permissive boundary conditions (Harlow et al., 2019).

A related issue is ambiguity in the presymplectic potential under shifts $\Theta$1. For subregions, this produces explicit boundary ambiguities in the naive prescription $\Theta$2. A path-integral derivation resolves this by replacing the hypersurface integral with a contour integral around a partial Cauchy surface,

$\Theta$3

where $\Theta$4. This makes the symplectic form unambiguous and implies that all gauge transformations, “large” or “small,” lie in the kernel of $\Theta$5 for the subregion construction (Kirklin, 2019).

The same boundary analysis is also the setting in which the Poisson bracket on covariant phase space is shown to coincide directly with the Peierls bracket, without non-covariant intermediate steps (Harlow et al., 2019).

4. Gravitational realizations: metric, tetrad, and Palatini formulations

In gravity, the method is often presented through the comparison of metric and tetrad formulations on manifolds with boundary. For metric Einstein–Hilbert plus Gibbons–Hawking–York,

$\Theta$6

The presymplectic potential is

$\Theta$7

For the tetrad formulation with fields $\Theta$8,

$\Theta$9

$\Theta$0

and one finds

$\Theta$1

A direct pull-back argument gives, off shell,

$\Theta$2

and hence on shell

$\Theta$3

The same pattern holds for Palatini versus Einstein–Hilbert variables, yielding on shell

$\Theta$4

In the tetrad description, the only extra degeneracies are those along internal Lorentz rotations in $\Theta$5 (Varo, 2023, G. et al., 2021).

This equivalence result addresses a long-standing comparison problem between metric and tetrad covariant phase spaces with boundaries. In the cohomological formulation, the two presymplectic structures differ only by relative exact terms, which vanish after integration over a Cauchy pair by the relative Stokes theorem (G. et al., 2021).

5. Integrability, slicing dependence, and flux

Surface-charge integrability is one of the most technically sensitive parts of the method. In the standard expression,

$\Theta$6

$\Theta$7 need not be exact. The Wald–Zoupas criterion requires $\Theta$8, or equivalently $\Theta$9 for a chosen slicing. A central result of the recent geometric SPS analysis is that this criterion is inherently slicing-dependent. The Frobenius theorem on SPS then yields a slicing-independent condition: $$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$0

necessary and sufficient for the existence of an integrable slicing (Golshani et al., 15 Mar 2026).

This geometric reformulation introduces an SPS connection, torsion, and curvature. In Cartan form,

$$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$1

Within this language, “fake” flux is encoded in the pure-gauge connection $$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$2, while “genuine” flux is measured by non-vanishing torsion $$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$3. The latter is slicing-independent and directly proportional to the physical gravitational News tensor or bulk symplectic flux at the boundary. The same framework gives a boundary Liouville theorem when $$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$4 and classifies theories according to whether torsion and curvature vanish (Golshani et al., 15 Mar 2026).

A different resolution of non-integrability appears when the boundary itself is allowed to fluctuate. Using a step-function reformulation on a boundary-less manifold, one obtains a total presymplectic potential

$$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$5

and a symplectic form with both bulk and corner pieces. In this setup, all surface charges become integrable once $$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$6, and the charge takes the form

$$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$7

The same formalism yields a charge algebra with a possible central term, and a flux law $$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$8, with nonzero $$\omega[\delta_1\Phi,\delta_2\Phi;\Phi] =\delta_1\Theta[\delta_2\Phi;\Phi]-\delta_2\Theta[\delta_1\Phi;\Phi],$$9 encoding radiation memory (Adami et al., 2024).

These developments suggest that some charge non-integrability previously attributed directly to radiation can instead arise from solution-space slicing choices, while other instances correspond to genuine boundary flux. The distinction is explicit only in the newer SPS geometry (Golshani et al., 15 Mar 2026).

6. Representative extensions and applications

The method has become a standard tool for boundary-sensitive problems in gravity and gauge theory. In asymptotically locally AdS spacetimes with non-smooth conformal boundaries, a covariant phase-space treatment with relaxed boundary conditions and codimension-two corner terms matches holographic Noether charges to Wald Hamiltonians. For accelerating, non-rotating dyonic AdS$\Sigma$0 black holes, this yields the first law with the “thermodynamic length” terms associated with cosmic-string tensions (Kim et al., 2023).

In BRST language, the formalism has been extended to a trigraded BRST covariant phase space. Starting from a BRST-invariant gauge-fixed Lagrangian, one obtains a BRST-invariant presymplectic form whose graded canonical brackets reduce to the classical covariant phase-space brackets when ghosts, antighosts, and antifields are set to zero. In that framework, the global BRST Noether charges are equivalent to the usual classical corner charges of large gauge transformations, and a unified Ward identity separates into the BRST–BV quantum master equation in the bulk and a boundary identity for large gauge transformations (Baulieu et al., 2024).

A semi-classical generalization replaces classical matter by a quantum state and defines a semi-classical symplectic two-form as the sum of the gravitational symplectic form and the Berry curvature of the matter state,

$\Sigma$1

This two-form is slice-independent and satisfies a quantum generalization of the Hollands–Iyer–Wald identity. For gauge-invariant subregions, the quantum contribution is replaced by the Berry curvature of special purifications involving the Connes cocycle (Bhattacharya et al., 22 Oct 2025).

The method has also been used to analyze boundary-supported phase spaces in string theory. For tensile open strings in a constant Kalb–Ramond background, the presymplectic current splits into a bulk kinetic term plus an exact boundary term, and imposing mixed boundary conditions yields a boundary symplectic form whose inverse reproduces the Seiberg–Witten noncommutativity parameter. In the tensionless limit with constant $\Sigma$2-field, the bulk contribution vanishes and the symplectic form localizes entirely on the endpoints, so the physical phase space becomes purely boundary-supported (Das et al., 14 Apr 2026).

Across these applications, the recurrent theme is that the covariant phase-space method is not merely a recipe for $\Sigma$3. In its contemporary form, it is a boundary-aware, cohomological, and increasingly geometric framework for presymplectic structures, Noether currents, Hamiltonians, fluxes, and charge algebras in diffeomorphism-invariant theories (Margalef-Bentabol et al., 2020).