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.
Update 2018-03-26: I’m making public a companion xAct notebook which contains most of the calculations on this web page. You can also get this notebook from the xAct-contrib/examples repository.
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. Wald1). 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 Shapiro3). 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 Shapiro3 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 (see the companion notebook). 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 (see the companion notebook). The result is simply