1. Introduction
In the curved space-time of general relativity, gravitational waves propagate with the speed of light (the velocity limit), correcting the Newtonian description of gravity. Whenever a reference system is chosen, time needs to be singled out. In the 3 + 1 decomposition of space-time, known as the Arnowitt–Deser–Misner (ADM) formalism of gravity [
1], the constant time three-surfaces form a foliation. While the time parameter is constant on each hypersurface, it changes monotonically from one hypersurface to the other. The role of the four-dimensional metric
(with
independent components) is taken by the metric induced on the three-dimensional hypersurfaces (six variables) and their extrinsic curvature (six variables), which generate canonical pairs. Einstein equations are replaced by the Hamiltonian evolution of these canonical pairs. Beside these, there are constraint equations to be fulfilled in each instant (on each hypersurface). These are the Hamiltonian and diffeomorphism constraints. In the generic case of the 3 + 1 decompositions, there is no preferred time; the formalism has to be valid for any possible temporal choices (“many-fingered time” formalism [
2,
3,
4]). Preferred choices arise by either imposing coordinate conditions [
5,
6] or by filling space-time with an adequate reference fluid [
7,
8]. Although the 3 + 1 decomposition breaks space-time covariance, a manifestly covariant canonical formalism based on the hyperspace (defined by all space-like hypersurfaces) has been also proposed [
9,
10,
11,
12]. The ADM decomposition has been generalized for some of the modified gravity theories, as well, e.g., for the
gravitational theories [
13].
If in addition, a spatial direction plays a special role, the 2 + 1 + 1 decomposition of space-time may prove useful. This special role can be provided by a Killing-symmetry, for example the radial directions in either spherical or cylindrical symmetric space-times are such singled-out directions. We do not explore however simplifications arising from the imposition of symmetries, as for example applying mini-superspace or midi-superspace approaches [
14,
15]. The scenario we have in mind is to discuss generic perturbations of a background with certain symmetry. In the most generic case, both singled-out directions have expansion, shear and vorticity [
16]. The corresponding optical scalars were explored in the discussion of perturbations of spherically-symmetric space-times [
17], also for the discussion of gravitational waves in anisotropic Kantowski–Sachs space-times [
18].
In another, much simpler 2 + 1 + 1 decomposition formalism, the decomposition is made along a perpendicular double foliation [
19,
20]. This formalism has been employed in the framework of dark matter and dark energy-motivated scalar-tensor gravitational theories in the discussion of the odd sector perturbations of spherically-symmetric gravity in the effective field theory approach [
21]. The requirement of perpendicularity however consumes one gauge degree of freedom by fixing a metric function to vanish. This has posed no problem in the discussion of the odd sector; however for the even sector, it generates an arbitrary function in the solution, hampering the physical interpretation of perturbations. Therefore, a modified 2 + 1 + 1 decomposition formalism would be desirable, which keeps the relative simplicity of the formalism of [
19,
20] (as compared to the formalism exploring optical scalars [
16]), but employs 10 metric functions instead of nine, hence becoming suitable for the discussion of the even sector. Such a formalism could be worked out at the price of relaxing the perpendicularity requirement [
22].
In this conference report, we summarize the main feature of this new formalism and sketch the derivation of the Hamiltonian formalism, without insisting on the involved computational details and related proofs of the statements, which are given in [
22] together with additional details.
Latin indices denote four-dimensional space-time indices. Boldface lower-case (as ) or uppercase (as ) Latin letters count two-dimensional or four-dimensional basis vectors.
2. The Nonorthogonal Double Foliation
We generalize the orthogonal 2 + 1 + 1 decomposition of [
19,
20] such that the hypersurfaces
of constant
t and
of constant
with normal vectors
and
, respectively, are nonorthogonal, as presented in
Figure 1.
Their intersection is the surface
, with an adapted vector basis
. We introduce two orthonormal bases adapted to the two foliations, as follows:
and
. The four-dimensional metric can then be decomposed in both:
Here, is the metric induced on .
The temporal and selected spatial evolution vectors in the
basis are:
They define a coordinate-basis, the duality relations of which imply [
22]:
making manifest that
is tangent to
. The shift component
arises due to the non-orthogonality of the foliations and generates all new terms arising as compared to the formalism presented in References [
19,
20], where
was imposed. With the introduction of a nonvanishing
, full gauge freedom is re-established, with 10 metric components in the formalism (three for
, two for
and
each, one for each of the lapses
N,
M and shift component
). At times, it will be convenient to parametrize this 10th metric function as
and also employ the notations
,
. This is especially convenient in proving [
22] that the two bases are related by a Lorentz rotation:
also to derive the decomposition of the evolution vectors in the basis
:
Note that is manifestly tangent to the hypersurface .
Finally, in the basis
, it is straightforward to check:
reassuring (due to the Frobenius theorem) that
is hypersurface-orthogonal and:
implying that the vector
has vorticity. Similarly, in the basis
, we find that
is hypersurface-orthogonal and
has vorticity:
Hence the 10th metric function
bears a double interpretation: (1) it gives the angle of the Lorentz rotation between the two bases and, (2) generates the vorticity of the complementary basis vectors
and
. More details of these interpretations are presented in [
22].
3. The 2 + 1 + 1 Decomposition of Covariant Derivatives
The projected covariant derivative of any tensor
defined on
arises by projecting in all indices with
:
The
D-derivative obtained in this way is related to the connection compatible with the two-metric due to the property:
It will be of particular importance to 2 + 1 + 1 decompose the covariant derivatives of the basis vectors. We found:
where
,
,
and
are extrinsic curvatures of the surface
;
and
are normal fundamental forms;
,
,
and
are normal fundamental scalars [
23]. The quantities
and
are defined similarly to the normal fundamental forms, but they also contain the contributions of the vorticities of the corresponding vectors. Finally
,
,
and
are the projections onto
of the nongravitational accelerations of the respective observers (among which those moving along
and
are physical).
The set of above quantities is not independent. As shown in detail in [
22], it is enough to select the sets
,
,
and
in order to express all the others. In particular, for orthogonal foliations, all starry quantities reduce to non-starred ones. Beside, the set
is related to time derivatives of the metric variables:
while the set
is connected to their
-derivatives only:
Moreover, the accelerations can be expressed as
D-derivatives of the lapses:
4. Hamiltonian Dynamics
The Einstein–Hilbert action:
can be rewritten by employing the twice contracted Gauss identity [
22] and the decomposition
as:
which besides scalars contains only tensors and vectors defined on
. The total covariant divergence is not yet decomposed; however upon decomposition, it will generate only boundary terms. The set of variables
comprises the generalized coordinates,
the generalized velocities, while
can be perceived as shorthand notations for the
-derivatives of the generalized coordinates. Similarly to the 3 + 1 decomposition, time derivatives of
do not emerge in the action. The generalized momenta arise as derivatives with respect to the time derivatives of the generalized coordinates as:
Then, the action can be rewritten in an already Hamiltonian form as [
22]:
where
Q is a sum of boundary terms, given explicitly in [
22], while:
is the Hamiltonian constraint,
and:
are the diffeomorphism constraints. Note that as expected,
only appear as Lagrange-multipliers.
The evolution equations for the generalized coordinates and momenta then emerge as the Hamiltonian equations written for the gravitational Hamiltonian density:
They are explicitly worked out in [
22].
5. Conclusions
We generalized the formalism of [
19,
20] by allowing for nonorthogonal foliations. As the main benefit, this led to the reestablishment of the full gauge freedom, allowing a generic discussion of perturbations. We gave a two-fold geometrical interpretation the 10th metric variable as the angle of the Lorentz rotation of the basis vectors and the measure of the vorticity of the basis vectors.
In the ADM formalism, the induced metric and extrinsic curvature of the hypersurface play the role of Hamiltonian coordinates and momenta. In the new formalism, we identified those geometrical quantities characterizing the embedding, which bear a dynamical role (they contain time derivatives). Non-dynamical geometrical quantities appear only in the basis ; hence, we employed that for the 2 + 1 + 1 decomposition of the Einstein–Hilbert action. From among the geometric variables, we identified those that combine into canonical pairs and proceeded with performing the Hamiltonian analysis. We identified the 2 + 1 + 1 decomposed gravitational Hamiltonian, as well as the Hamiltonian and momentum constraints in terms of canonical coordinates and momenta.
We intend to apply this formalism both for the discussion of the even sector of perturbations of spherically-symmetric gravity in the effective field theories of gravity and for the Hamiltonian treatment of canonically-quantizable cylindrical gravitational waves. The first of these has the potential to address the stability of dark matter halo models in scalar-tensor gravity. Furthermore, for the discussion of gravitational waves in space-times with particular symmetries, the 2 + 1 + 1 decomposition of the Weyl-tensor would be an asset.