Covariant 1+3 Decomposition in GR

• The covariant 1+3 framework is a geometric, observer-adapted method that splits spacetime into temporal and spatial parts using a unit timelike vector field.
• It decomposes kinematic and dynamical equations—such as expansion, shear, and vorticity—into evolution and constraint components for precise analysis.
• The formalism extends to applications like cosmological perturbations, gravitoelectromagnetism, and backreaction studies, offering a robust tool for GR research.

The covariant 1+3 description of General Relativity ("threading" formalism) is a geometric, coordinate-free framework for splitting spacetime into temporal and spatial components relative to a choice of observer congruence. This method introduces a unit timelike vector field $u^a$—interpreted as the 4-velocity of a class of observers—which underpins the decomposition of all tensor fields. The 1+3 covariant approach establishes a closed first-order system of kinematical and dynamical equations, connects geometric properties to observer measurements, and naturally accommodates generic symmetry reductions, fluid models, and perturbative analyses. Unlike the traditional 3+1 (ADM) foliation, 1+3 is adapted to observer worldlines and enables manifestly covariant treatment of quantities including vorticity and non-geodesic acceleration (Tsamparlis et al., 2018, Park, 2018, &&&2&&&, Danehkar, 2022, Magni, 2012).

1. Observer Congruence and Projection Tensor

Let $(M,g_{ab})$ be a Lorentzian 4-manifold with metric signature $(-,+,+,+)$. The observer congruence $u^a$ satisfies $u^a u_a = -1$ and defines at each point a local rest space—spanned by the projection tensor $h_{ab} = g_{ab} + u_a u_b$. The following identities hold:

• $h_{ab}u^b=0$ (spatiality),
• $h_{ab}=h_{ba}$ (symmetry),
• $h_a{}^{c}h_{c}{}^{b}=h_a{}^{b}$ (idempotence),
• $h_a{}^{a}=3$.

For any tensor $T$, its purely spatial part is obtained by full projection of all indices with $h_{ab}$. This splitting is the foundation for decomposing the geometry and dynamics of spacetime into observer-adapted temporal and spatial domains, without a choice of coordinates or basis (Park, 2018, Roy, 2014, Tsamparlis et al., 2018).

2. Kinematic Decomposition of $\nabla_b u_a$

The covariant derivative $\nabla_b u_a$ splits uniquely as:

$$\nabla_b u_a = -u_a \dot u_b + \omega_{ab} + \sigma_{ab} + \tfrac{1}{3}\theta h_{ab}$$

• 4-acceleration: $\dot u_a := u^b\nabla_b u_a$ ($\dot u_a u^a = 0$),
• Expansion scalar: $\theta := \nabla_a u^a = h^{ab}\nabla_a u_b$, measuring volume change rate,
• Shear: $$\sigma_{ab} := h_a{}^c h_b{}^d (u_{(c;d)} - \tfrac{1}{3}\theta h_{cd})$$ (symmetric, trace-free, spatial),
• Vorticity: $\omega_{ab} := h_a{}^c h_b{}^d u_{[c;d]}$ (antisymmetric, spatial), encapsulates local rotation about $u^a$.

Associated invariants:

• $\sigma^2 = \frac{1}{2}\sigma_{ab}\sigma^{ab}$,
• $\omega^2 = \frac{1}{2}\omega_{ab}\omega^{ab}$.

The geometrical and physical interpretation of these quantities is standard: $\theta$ controls expansion, $\sigma_{ab}$ shape distortion, $\omega_{ab}$ local rotation, and $\dot u_a$ non-geodesic acceleration (Tsamparlis et al., 2018, Park, 2018, Roy, 2014, Magni, 2012, Danehkar, 2022).

3. Covariant 1+3 Decomposition of Field Equations

Energy-Momentum Tensor

Any symmetric tensor $T_{ab}$ (the energy-momentum content) splits as:

$$T_{ab} = \mu u_a u_b + 2q_{(a}u_{b)} + p h_{ab} + \pi_{ab}$$

• $\mu := T_{ab}u^a u^b$: energy density,
• $p := \frac{1}{3}h^{ab}T_{ab}$: isotropic pressure,
• $q_a := -h_a{}^b T_{bc}u^c$, ($q_a u^a = 0$): momentum flux/heat flux,
• $\pi_{ab} := h_a{}^c h_b{}^d(T_{cd} - p h_{cd})$: anisotropic stress (trace-free, symmetric, spatial).

These quantities become the dynamical variables determining the gravitational evolution in the frame of the defined observers (Tsamparlis et al., 2018, Roy, 2014, Danehkar, 2022).

Propagation and Constraint Equations

The Ricci identity ($2\nabla_{[a}\nabla_{b]}u_c = R_{abc}{}^{d}u_d$) and Einstein's equations $G_{ab}+\Lambda g_{ab}=T_{ab}$ decompose as follows:

• Raychaudhuri equation (volume evolution):

$$\dot\theta = -\tfrac13\theta^2 - 2(\sigma^2 - \omega^2) + \nabla^a\dot u_a - \tfrac12(\mu+3p) + \Lambda$$

• Shear propagation:

$$h_a{}^c h_b{}^d \dot\sigma_{cd} = -\tfrac23\theta\sigma_{ab} - \sigma_{c\langle a}\sigma_{b\rangle}{}^c - \omega_{c\langle a}\omega_{b\rangle}{}^c + \nabla_{\langle a}\dot u_{b\rangle} + E_{ab} - \tfrac12\pi_{ab}$$

• Vorticity propagation:

$$h_a{}^c h_b{}^d \dot\omega_{cd} = -\tfrac23\theta\omega_{ab} + 2\sigma_{[a}{}^c\omega_{b]c} - \tfrac12 h_a{}^c h_b{}^d \nabla_{[c}\dot u_{d]}$$

Constraint relations include vorticity divergence, shear–vorticity constraint, momentum-constraint, and energy-constraint equations. These collectively govern initial data and ensure compatibility with the full system (Tsamparlis et al., 2018, Park, 2018, Roy, 2014, Magni, 2012).

4. Covariant Split of Curvature and Weyl Gravitoelectromagnetism

Riemann and Weyl Tensor Decomposition

The Riemann curvature tensor $R_{abcd}$ is split into its Ricci and Weyl (conformal) components:

$R_{abcd} = C_{abcd} + \text{matter terms (Ricci)}$

The Weyl tensor $C_{abcd}$ is decomposed via the observer congruence into:

• Electric part: $E_{ab} := C_{acbd}u^c u^d$ (symmetric, trace-free, spatial),
• Magnetic part: $$H_{ab} := \frac{1}{2}\epsilon_{acd}C^{cd}{}_{be}u^e$$ (symmetric, trace-free, spatial).

These gravitoelectric and gravitomagnetic fields represent tidal, frame-dragging, and gravitational wave phenomena. Their evolution and differential constraints—derived from Bianchi identities—form a system formally analogous to Maxwell's equations, but with coupling to fluid and kinematic variables (Danehkar, 2022, Tsamparlis et al., 2018, Magni, 2012).

Maxwell-like System

The propagation and constraint equations are:

Equation	Main Dependents	Purpose
Electric propagation	$E_{ab},\,H_{ab},\,\sigma_{ab},\,\omega_{ab},\,a_a$	Evolution of gravitoelectric
Magnetic propagation	$E_{ab},\,H_{ab},\,\sigma_{ab},\,\omega_{ab},\,a_a$	Evolution of gravitomagnetic
Electric constraint	$E_{ab},\,\pi_{ab},\,\mu$	Compatibility in divergence
Magnetic constraint	$H_{ab},\,q_a,\,\omega_a$	Compatibility, divergence

Analogous to Gauss's, Ampère's, and Faraday's laws in electromagnetism, these equations manifest the free gravitational degrees of freedom and their interaction with matter and spacetime kinematics (Danehkar, 2022, Tsamparlis et al., 2018, Magni, 2012).

5. Geometric Relations: Gauss–Codazzi–Ricci Splitting

By projection of $R_{abcd}$, spatial frames induce three fundamental geometric relations (Park, 2018, Roy, 2014):

• Gauss relation: The spatial curvature involves both intrinsic 3-space curvature and second fundamental form,
• Codazzi relation: Mixed time-space derivatives encode spatial variations of extrinsic curvature,
• Ricci relation: Temporal evolution of extrinsic curvature encodes interaction between geometric and physical variables.

$$\begin{gathered} h_a{}^p h_b{}^q h_c{}^r h_d{}^s R_{pqrs} = {}^{(3)}R_{abcd} + B_{ac}B_{bd} - B_{ad}B_{bc} \ u^p h_a{}^q h_b{}^r h_c{}^s R_{pqrs} = D_b B_{ac} - D_c B_{ab} + 2a_{a}\omega_{bc} \ u^p u^q h_a{}^r h_b{}^s R_{pqrs} = D_b a_a + a_a a_b + B_{ac}B^c{}_b - \dot B_{ab} \end{gathered}$$

These underpin the geometric source terms in the evolution and constraint equations (Park, 2018, Roy, 2014).

6. Extensions: Tetrad, 1+1+2 Splitting, Backreaction

Tetrad and 1+1+2 Semi-Covariant Methods

Generalizing further, a local orthonormal tetrad $\{e_a\}$ can be introduced with $e_0=u$. Ricci identities yield a first-order symmetric hyperbolic PDE system for all relevant variables. The 1+1+2 split refines the spatial decomposition by selecting a preferred spacelike direction $n_a$, facilitating analysis of problems with further symmetry, such as axisymmetry or radial perturbations. All derivative operators and tensor decompositions split accordingly, yielding a finer-grained evolution and constraint system (Danehkar, 2022, Tsamparlis et al., 2018).

Backreaction in Cosmology

Application of the 1+3 covariant formalism in cosmological averaging yields the Buchert backreaction framework, where noncommutativity of time and spatial averaging renders a correction to the expansion law:

$${\cal Q} = \tfrac{2}{3}(\langle\Theta^2\rangle - \langle\Theta\rangle^2) - 2\langle\sigma^2\rangle$$

The averaged Friedmann equations are modified accordingly, and in irrotational dust models this correction may drive late-time acceleration in cosmology without resort to exotic matter (Magni, 2012).

7. Comparison with 3+1 Formalism and Physical Applications

When the observer congruence $u^a$ is hypersurface-orthogonal ($\omega_{ab}=0$), the 1+3 covariant equations reduce to the standard 3+1 (ADM) evolution and constraint equations, with the extrinsic curvature becoming symmetric and the spatial connection torsion-free (Park, 2018, Roy, 2014). The computational distinction is:

• 3+1: Eulerian (fields on hypersurfaces), coordinate- and frame-dependent,
• 1+3: Lagrangian (fields along worldlines), manifestly covariant.

Physical applications span relativistic fluid models, cosmological perturbation theory, analysis of gravitational waves, and studies of exact solutions. The physical meaning of each variable is manifest, facilitating their use in both analytic and numerical studies. The extensions to multi-fluid models, inclusion of tilt effects, and detailed gravitoelectromagnetic analysis enable comprehensive treatment of observational signatures—including gravitational lensing, redshift, and wave propagation—directly from the geometric structure (Danehkar, 2022, Tsamparlis et al., 2018, Magni, 2012).

---

This covariant 1+3 scheme provides a rigorous and transparent foundation for the analysis of general relativistic systems, maintaining manifest covariance, clarifying observer-dependent physics, and unifying geometric, kinematical, and dynamical perspectives (Tsamparlis et al., 2018, Park, 2018, Roy, 2014, Danehkar, 2022, Magni, 2012).