# Notes on the Weyl, E/B, and 3+1 decompositions of Riemann, and curvature invariants

More notes I had that needed to live somewhere. I hope you find them useful.

# Weyl decomposition

The most basic ingredient here is splitting the Riemann tensor into
trace and tracefree pieces. This should be covered in every standard
GR textbook (see e.g. Wald^{1}). The decomposition is, in
dimensions,

where is the tracefree Weyl tensor, and as usual, the Ricci tensor is and the Ricci scalar is .

Now, this is not a very pleasant expression to remember. It can be
brought into a slightly nicer form by defining an ‘impose index
symmetry’ operation with some group , which is a subgroup of the
signed permutation group
(signed symmetric group / hyperoctahedral group), which is just the
symmetric group where each
permutation is further given a sign, . The ‘Riemann
symmetric’ group is generated by the elements . That is, the Riemann tensor is
separately antisymmetric on , separately antisymmetric on
, and symmetric under the *simultaneous* interchange of
. You can check that is a group
with 8 elements.

The ‘impose index symmetry’ operation on a tensor expression is given by

where is the order of the group, is the sign
(+1 or -1) associated with the element , and
means to permute indices as specified by the permutation
. This is a generalization of the symmetrization and
antisymmetrization operations that are commonly used,
corresponding to respectively the whole symmetric group with all plus
signs, or the whole symmetric group where signs are
determined by the parity of the permutation . This general
operation is implemented in xAct as
`ImposeSymmetry[]`

.

Having defined this operation, we can (slightly) improve Eq. \eqref{eq:R-decomp-1} via

# Curvature invariants 1

Each of the Weyl/Ricci tensor/Ricci scalar pieces of the decomposition
are orthogonal to each other in the following way. Construct the
*Kretschmann invariant*, . In terms of the
Weyl decomposition, it is given by

Note: the coefficients in Eq. \eqref{eq:K-Weyl-decomp} are the same as those appearing in Eq. \eqref{eq:R-decomp-proj}. This has something to do with the impose-symmetry operation being a properly normalized idempotent projection operation. In this case, it is actually a Young projector for the standard Young tableau filled with {{1,3},{2,4}}.

From here forward we focus on dimension . Now we can also define the (left) dual Riemann tensor, via

With this dual, we can define the *Pontryagin invariant*,

where is the (left) dual Weyl tensor, defined in the same way as the dual to Riemann. Interestingly, the Pontryagin is completely insensitive to the Ricci tensor.

For future purposes, we will also define a complex tensor where the real part is Weyl and the imaginary part is the (left) dual of Weyl,

The square of this tensor contains both Kretschmann (minus the Ricci parts) and Pontryagin:

This quantity is related to the so-called invariant which is
usually defined in spinor language (see e.g. Stephani et al.^{2}). The completely symmetric ‘Weyl
spinor’ related to the tensor is

where tensor indices are related to spinor index pairs and so on. In other words, the soldering form (or van der Waerden symbols) are to be inferred on the LHS of Eq. \eqref{eq:tilde-C-to-spinor}.

The aforementioned invariant is given by

where the Weyl scalars in some spinor dyad basis are defined as usual (the subscript counts the number of ’s contracted into the spinor),

Similarly, for the invariant,

A mnemonic for remembering in terms of the Weyl scalars is given by the determinant expression

# 3+1 and electric/magnetic decomposition

Now let us suppose we have a 3+1 decomposition of the spacetime,
i.e. that (a region of) the spacetime is foliated by spacelike
hypersurfaces (see standard references on this topic, e.g. Baumgarte
and Shapiro^{3}). We have a timelike unit normal , with
. This induces a splitting of the tangent spaces into
time components (parallel to ) and ‘spatial’ components
(orthogonal to ). We have the spatial metric induced on the
spatial hypersurface , with ,
given by

The metric-compatible derivative of is , the metric-compatible derivative of is , the extrinsic curvature is .

To avoid confusion, we will now label the 4-Riemann (and 4-Ricci) with , and the 3-Riemann (and 3-Ricci, a.k.a. “thricci”) with .

We can now decompose all 4-dimensional objects using this splitting of the tangent bundle. We will need to split the 4-dimensional volume form, via (Warning: conventions may vary!)

There is no need for 3/4 superscripts because just reading the number of indices tells you whether this is the 3-volume or 4-volume form.

Now we proceed to splitting the Weyl tensor into space- and time pieces. This proceeds by analogy to the Maxwell tensor. The “electric” part of Weyl is defined via

while the “magnetic” part of Weyl is defined via (Warning: conventions may vary!)

From these two it is obvious that

Each of and are symmetric and trace-free spatial tensors. This is as far as you can decompose Weyl without imposing additional structure (like some preferred direction).

We can re-express Weyl in terms of , , the timelike unit vector and spatial metric, and 3-epsilon tensor. I don’t know if there is a clever way to do this other than to take the below expression and verify that it’s correct. Let me know if you have some better insight. Regardless, the most compact way of writing this is again using the above ‘impose index symmetry’ operation in Eq. \eqref{eq:impose-sym-def} with , the index symmetries of Riemann. The expression is (thanks to Alfonso García-Parrado Gómez-Lobo)

# Curvature invariants 2

Obviously we would like to express curvature invariants in terms of the electric and magnetic Weyl tensors. This is a snap if you use xTensor (otherwise there is a lot of algebra). The relationship is

Now it is easy to see that Weyl squared and Pontryagin are simply

Kretshmann of course also adds the 4-Ricci terms to Eq. \eqref{eq:Weyl-squared-in-E-B}, via the relationship in Eq. \eqref{eq:K-Weyl-decomp}.

Similarly, the invariant turns out to be

# Computing from Gauss-Codazzi relations

Let’s suppose you have at your disposal a numerical relativity codebase (such as SpEC) and thus you have access to the 3-metric , extrinsic curvature , the 3-Ricci and 4-Ricci ; and you would like to compute and . Certainly you should not attempt to construct the full Riemann or Weyl tensors, because they have 256 numbers per grid point but there are only 10 algebraically independent components in and !

Therefore we would like an expression for how to compute
and in terms of the available quantities, without all the
overhead of computing Riemann or Weyl. To do so requires using the
Gauss-Codazzi equations (see
e.g. Baumgarte and Shapiro^{3} for a didactic exposition). The Gauss
and Codazzi equations relate the orthogonal decomposition of the
4-curvature to the intrinsic and extrinsic curvature of the
spatial submanifold. These are implemented
in xAct with the command `GaussCodazzi[]`

.

The Gauss equation relates the completely spatial projection,

The recipe for finding a formula for is as follows:

- Use the 4-metric to trace over two indices
- Use the Weyl decomposition \eqref{eq:R-decomp-1} for on the left hand side
- Insert the E/B decomposition \eqref{eq:weyl-in-E-B} for the Weyl tensor on the left hand side
- You should now have with free indices and no prefactor

I highly recommend using xAct for this calculation, to avoid errors. The result is

The three terms on the second line all vanish when in a 4-Ricci-flat manifold, e.g. in vacuum in GR.

To find an expression for the magnetic part involves the Codazzi equation. This is the decomposition of 4-Riemann when one index is contracted with and the remaining three are projected onto spatial directions. The Codazzi equation is

The recipe for finding a formula for is as follows:

- Contract with the 3-epsilon
- Use the Weyl decomposition \eqref{eq:R-decomp-1} for on the left hand side
- Insert the E/B decomposition \eqref{eq:weyl-in-E-B} for the Weyl tensor on the left hand side
- Symmetrize on indices
- You should now have with free indices and no prefactor

Once again you should use xAct. The result is simply