Linearized Einstein–Bianchi System

• The linearized Einstein–Bianchi system is a gauge-theoretic framework that linearizes Einstein’s equations and Bianchi identities around fixed backgrounds, yielding Maxwell-like formulations.
• Different formulations across flat, Einstein, and 3+1 settings highlight key aspects of wave propagation, gauge invariance, and the preservation of physical constraints.
• Discrete and computational realizations employ finite element, lattice, and conformal-Hessian methods to accurately simulate gravitational phenomena while preserving constraint integrity.

Searching arXiv for recent and foundational papers on the linearized Einstein–Bianchi system. The linearized Einstein–Bianchi system is the gauge-theoretic and constraint-propagating structure obtained by linearizing Einstein’s equations and their Bianchi identities around a fixed background, most commonly Minkowski spacetime or an Einstein background. In the weak-field flat-background setting, it appears as the linearized Einstein equations together with the divergence-free structure inherited from the contracted Bianchi identities, and it yields a Maxwell-like gravitoelectromagnetic system for the gravitoelectric and gravitomagnetic fields (Álvarez-Samaniego et al., 2017). On general Einstein backgrounds, the same structure can be formulated invariantly as a differential complex in which infinitesimal diffeomorphisms, metric perturbations, and linearized field equations are linked by natural bundle maps; in that setting the vanishing of the composition of the gauge and field operators encodes the linearized Bianchi identity (Eastwood, 2022). In $3+1$ and related decompositions, the same phrase also refers to the subsidiary linear first-order system governing the propagation of Einstein constraints, derived from the contracted Bianchi identities and shown to be symmetric hyperbolic under appropriate hypotheses (Rácz, 2014). A complementary tradition treats the electric and magnetic parts of the Weyl tensor as the primary variables and studies a first-order curl-div system for symmetric trace-free tensors, together with finite element or lattice discretizations designed to preserve the underlying complex and constraints (Hu et al., 2021, Brewin, 2011).

1. Linearized geometric framework

The standard weak-field construction starts from a metric decomposition

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$

with $\eta_{ik}$ the Minkowski metric and $h_{ik}$ the perturbation (Álvarez-Samaniego et al., 2017). In that setting only terms linear in $h_{ik}$ are retained, source velocities are taken to satisfy $|\vec v|\ll c$, and only first-order expressions for the curvature are used. The trace $h = h^i{}_i$ is introduced, and the trace-reversed perturbation

$\Psi^{ik} := h^{ik} - \eta^{ik} h$

is used together with the harmonic gauge condition

$\partial_i \Psi^{ik} = 0.$

Within these assumptions the linearized Riemann tensor is

$$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$

while the Ricci tensor and scalar curvature reduce to

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$0

with

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$1

(Álvarez-Samaniego et al., 2017).

On a general Einstein background $$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$2 satisfying

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$3

the perturbation is written

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$4

and the linearized geometry is encoded using the Calabi operator $$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$5 from projective differential geometry (Eastwood, 2022). In that formulation the linearized Riemann tensor is

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$6

and the linearized Einstein vacuum equation with cosmological constant is expressed through the Ricci contraction

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$7

as

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$8

(Eastwood, 2022). In explicit form,

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$9

and on an Einstein background this becomes

$\eta_{ik}$0

(Eastwood, 2022).

These two formulations are the flat-background and background-covariant versions of the same structure. The first emphasizes harmonic gauge and explicit wave operators; the second packages the system as a differential complex and makes gauge invariance intrinsic. A plausible implication is that the term “linearized Einstein–Bianchi system” has no single canonical realization, but rather a family of equivalent realizations adapted to flat-space perturbation theory, Einstein-background deformation theory, $\eta_{ik}$1 PDE analysis, or curvature-based formulations.

2. Einstein equations, Bianchi identities, and gauge structure

In the weak-field flat-background setting, Einstein’s equations

$\eta_{ik}$2

reduce, after using the harmonic gauge, to

$\eta_{ik}$3

where

$\eta_{ik}$4

up to the paper’s normalization (Álvarez-Samaniego et al., 2017). This is the “Einstein” part of the system. The “Bianchi” part comes from the full identities

$\eta_{ik}$5

and their contraction

$\eta_{ik}$6

which in linearized form becomes

$\eta_{ik}$7

Einstein’s equations then imply

$\eta_{ik}$8

so source conservation is not an independent assumption but a compatibility condition enforced by the Bianchi identities (Álvarez-Samaniego et al., 2017).

On Einstein backgrounds, the same compatibility is encoded by the complex

$\eta_{ik}$9

where

$h_{ik}$0

is the Killing operator and $h_{ik}$1 is the Einstein deformation operator (Eastwood, 2022). The key structural statement is

$h_{ik}$2

This means that pure-gauge perturbations $h_{ik}$3 automatically solve the linearized equations, and that the field equations satisfy differential identities expressing linearized diffeomorphism invariance (Eastwood, 2022). In this language, gauge-equivalent perturbations differ by elements of $h_{ik}$4, and physical perturbations modulo gauge correspond to

$h_{ik}$5

A recurrent misconception is that the Bianchi identity merely supplies a divergence condition after the field equations have been imposed. The invariant complex formulation shows a stronger statement: the linearized Bianchi identity is built into the operator structure itself through the annihilation of pure gauge modes by the field operator (Eastwood, 2022). In the flat weak-field picture, the same point appears as the preservation of the harmonic gauge and conservation of the linearized sources (Álvarez-Samaniego et al., 2017).

3. Maxwell-type and Weyl-curvature formulations

One of the most widely used realizations of the linearized Einstein–Bianchi system is the gravitoelectromagnetic or Weyl-electric/Weyl-magnetic formulation. In the weak-field approximation around Minkowski spacetime, the perturbation components define a gravitoelectric field $h_{ik}$6, a gravitomagnetic field $h_{ik}$7, a vector potential $h_{ik}$8, a space–time–mass density $h_{ik}$9, and a space–time–mass current density $h_{ik}$0, leading to

$h_{ik}$1

$h_{ik}$2

These equations, together with the continuity law

$h_{ik}$3

are identified as the GEM manifestation of the linearized Einstein–Bianchi system (Álvarez-Samaniego et al., 2017).

In the vacuum case, this reduces to the wave system

$h_{ik}$4

and hence

$h_{ik}$5

(Álvarez-Samaniego et al., 2017). This is the most direct analogue of Maxwell theory, with the sign difference in Gauss’s law reflecting the attractive nature of gravity.

A more invariant three-dimensional curvature formulation uses the electric and magnetic parts of the Weyl tensor, denoted $h_{ik}$6 and $h_{ik}$7, both symmetric and trace-free. In the finite-element literature the first-order system is written

$h_{ik}$8

with curl and divergence acting row-wise on matrix fields (Hu et al., 2021). Introducing

$h_{ik}$9

the system is recast as

$|\vec v|\ll c$0

(Hu et al., 2021). This version is especially important because it identifies the linearized Einstein–Bianchi system with a Hodge wave equation on a tensor complex.

The 2025 conformal-Hessian formulation sharpens this viewpoint by taking $|\vec v|\ll c$1 and $|\vec v|\ll c$2 to be explicitly symmetric and traceless and by using

$|\vec v|\ll c$3

with $|\vec v|\ll c$4 (Guo et al., 6 Aug 2025). There the differential constraints are encoded by the exact conformal Hessian complex

$|\vec v|\ll c$5

(Guo et al., 6 Aug 2025). This suggests that the Weyl-electric/Weyl-magnetic formulation is not merely an analogy with electromagnetism, but a precise instance of a Hodge-type evolution system on a geometric complex.

4. Constraint propagation and hyperbolic subsidiary systems

In $|\vec v|\ll c$6 formulations, the linearized Einstein–Bianchi system frequently means the subsystem governing the propagation of the Hamiltonian and momentum constraints. Let

$|\vec v|\ll c$7

and define

$|\vec v|\ll c$8

These are the Hamiltonian and momentum expressions associated with a foliation by hypersurfaces $|\vec v|\ll c$9 (Rácz, 2014). Choosing the reduced evolution equations

$h = h^i{}_i$0

and setting $h = h^i{}_i$1, one uses the contracted Bianchi identity together with $h = h^i{}_i$2 to derive a first-order linear homogeneous system for $h = h^i{}_i$3 (Rácz, 2014).

In adapted coordinates, the subsidiary system has the form

$h = h^i{}_i$4

with

$h = h^i{}_i$5

When the foliating hypersurfaces are Riemannian, $h = h^i{}_i$6 is positive definite, so $h = h^i{}_i$7 is positive definite and the system is linear first-order symmetric hyperbolic (Rácz, 2014). The consequence is the standard constraint-propagation result: if the constraints vanish on one leaf and the reduced evolution equations hold, then the constraints vanish throughout the domain of dependence (Rácz, 2014).

This subsidiary hyperbolic viewpoint is analytically distinct from the GEM or Weyl-curl formulations, but it addresses the same compatibility problem. In both cases the Bianchi identities provide transport equations for constraints that are linear and homogeneous in the constraint variables. A plausible implication is that “linearized Einstein–Bianchi system” is best understood as a structural label for the compatibility layer of linearized gravity rather than as a single fixed PDE system.

Related hyperbolic formulations arise in orthonormal-frame and Einstein–Friedrich systems. For a minimally coupled nonlinear scalar field, a frame representation yields a quasi-linear first-order symmetric hyperbolic system for $h = h^i{}_i$8, where $h = h^i{}_i$9 and $\Psi^{ik} := h^{ik} - \eta^{ik} h$0 are the electric and magnetic parts of the Weyl tensor and the Bianchi equations provide their evolution (Alho et al., 2010). Linearization around expanding FLRW backgrounds then leads to a dissipative linearized Einstein–Bianchi system whose variables decay exponentially under suitable assumptions on the scalar potential (Alho et al., 2010). In the $\Psi^{ik} := h^{ik} - \eta^{ik} h$1 tetrad formalism, a different structural result shows that for generic timelike or null congruences the Einstein equations themselves can be regarded as integrability conditions for the Jacobi, Ricci, and Bianchi equations (Bergh, 2013). That result concerns the nonlinear theory, but it clarifies why linearized Einstein–Bianchi formulations can be closed on curvature and connection variables.

5. Discrete and computational realizations

Several computational approaches are built directly around the linearized Einstein–Bianchi structure. One line of work develops conforming finite element complexes. In the divdiv-based formulation, the unknowns are a scalar $\Psi^{ik} := h^{ik} - \eta^{ik} h$2, a symmetric tensor $\Psi^{ik} := h^{ik} - \eta^{ik} h$3, and a trace-free tensor $\Psi^{ik} := h^{ik} - \eta^{ik} h$4, and the semidiscrete weak system is

$\Psi^{ik} := h^{ik} - \eta^{ik} h$5

$\Psi^{ik} := h^{ik} - \eta^{ik} h$6

$\Psi^{ik} := h^{ik} - \eta^{ik} h$7

for all admissible test functions (Hu et al., 2021). The crucial structure is the exact discrete complex

$\Psi^{ik} := h^{ik} - \eta^{ik} h$8

which ensures discrete Bianchi identities, compatibility of operators, and inf-sup stability (Hu et al., 2021). With Crank–Nicolson time discretization, the fully discrete scheme achieves $\Psi^{ik} := h^{ik} - \eta^{ik} h$9 error under the stated regularity assumptions (Hu et al., 2021).

The conformal-Hessian approach advances this by enforcing both symmetry and tracelessness strongly at the level of finite element spaces. Its discrete complex is

$\partial_i \Psi^{ik} = 0.$0

exact for $\partial_i \Psi^{ik} = 0.$1 on contractible domains (Guo et al., 6 Aug 2025). The corresponding semidiscrete and fully discrete Hodge-wave systems preserve symmetry and tracelessness exactly and yield convergence

$\partial_i \Psi^{ik} = 0.$2

(Guo et al., 6 Aug 2025). This is explicitly presented as a discretization of the linearized Einstein–Bianchi system near Minkowski space.

A very different discretization is the smooth-lattice Einstein–Bianchi system. In the Schwarzschild reduction, the independent curvature variables are

$\partial_i \Psi^{ik} = 0.$3

subject to the vacuum constraint

$\partial_i \Psi^{ik} = 0.$4

(Brewin, 2011). Together with lattice leg lengths $\partial_i \Psi^{ik} = 0.$5 and extrinsic curvatures $\partial_i \Psi^{ik} = 0.$6, the evolution equations include

$\partial_i \Psi^{ik} = 0.$7

$\partial_i \Psi^{ik} = 0.$8

where the curvature equations are obtained from the second Bianchi identity (Brewin, 2011). Constraint preservation is explicit: $\partial_i \Psi^{ik} = 0.$9 so the Hamiltonian-like curvature constraint propagates (Brewin, 2011). The $$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$0 vacuum generalization evolves all 20 independent Riemann components, uses the full uncontracted Bianchi identities for 14 evolution equations, the vacuum Einstein equations for 6 algebraic relations, and shows that each curvature component satisfies a wave equation, establishing hyperbolic structure and preservation of the associated 10 constraints (Brewin, 2011).

In anisotropic cosmology, the linearized Einstein–Bianchi structure is also realized in Hamiltonian perturbation theory. For Bianchi I backgrounds, Fourier-expanded ADM perturbations are reorganized so that the linear scalar and vector constraints become canonical momenta

$$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$1

while three configuration variables $$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$2 are gauge invariant (Agullo et al., 2020). The reduced Hamiltonian is

$$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$3

with $$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$4, and the equations of motion show scalar–tensor mixing through the anisotropic effective potential matrix $$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$5 (Agullo et al., 2020). There the linearized Bianchi identities are reflected in the first-class constraint algebra and the automatic propagation of the linear constraints.

6. Background dependence, variants, and scope

The linearized Einstein–Bianchi system depends strongly on the chosen background and analytical setting. Around Minkowski space, the weak-field approximation yields the most transparent wave or Maxwell-type expressions (Álvarez-Samaniego et al., 2017). Around Einstein backgrounds with $$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$6, the operator $$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$7 and the complex

$$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$8

provide the natural invariant formulation (Eastwood, 2022). On compact Riemannian manifolds with boundary and arbitrary interior geometry, the relevant object is the linearized Einstein equation with sources in linear Bianchi gauge,

$$R_{iklm} \simeq \frac{1}{2}\left( \partial_k\partial_m h_{il} + \partial_i\partial_l h_{km} - \partial_k\partial_l h_{im} - \partial_i\partial_m h_{kl} \right),$$9

supplemented by extended Cauchy data on the boundary (Leder, 14 Oct 2025). In that setting global solvability is characterized by a lifted divergence condition

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$00

together with $$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$01, and uniqueness follows from generalized Hodge theory plus unique continuation (Leder, 14 Oct 2025). This is a genuinely elliptic Einstein–Bianchi theory rather than a hyperbolic one.

Cosmological applications supply further variants. In nearly isotropic homogeneous Bianchi cosmologies, the metric perturbation is encoded by a homogeneous trace-free tensor $$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$02, and the linearized Einstein–Bianchi system becomes a finite system of coupled or decoupled ODEs for mode amplitudes, depending on the Bianchi type (Pontzen et al., 2010). There the momentum constraint

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$03

plays the role of the linearized Einstein–Bianchi coupling between metric anisotropy and matter tilt (Pontzen et al., 2010). In maximal–isothermal gauge for axisymmetric vacuum perturbations of Minkowski space, the reduced variables $$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$04 and $$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$05 satisfy a singular hyperbolic–elliptic system,

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$06

$$g_{ik} \simeq \eta_{ik} + h_{ik}, \qquad |h_{ik}|\ll |\eta_{ik}|,$$07

which is solved explicitly by Fourier and Hankel transforms (Dain et al., 2010). The paper does not name it an Einstein–Bianchi system explicitly, but the wave variable and elliptic gauge variable together encode a linear hyperbolic–elliptic realization of the same constraint-preserving structure (Dain et al., 2010).

A common misconception is that all linearized Einstein–Bianchi systems are first-order symmetric hyperbolic systems in the same sense. The record is more varied. Some are hyperbolic curvature systems (Brewin, 2011), some are subsidiary symmetric hyperbolic systems for constraints (Rácz, 2014), some are Hodge-wave systems on tensor complexes (Hu et al., 2021, Guo et al., 6 Aug 2025), and some are elliptic boundary value problems in Bianchi gauge (Leder, 14 Oct 2025). What unifies them is not a unique PDE normal form, but the coexistence of linearized Einstein equations, gauge structure from infinitesimal diffeomorphisms, and compatibility identities inherited from the Bianchi identities.

The subject therefore occupies a central position between linearized general relativity, Weyl-curvature evolution, gauge complexes, constraint propagation, and structure-preserving discretization. In each of these settings, the linearized Einstein–Bianchi system serves as the mechanism by which evolution equations, gauge freedom, and consistency conditions are organized into a closed mathematical structure (Álvarez-Samaniego et al., 2017, Eastwood, 2022, Rácz, 2014, Hu et al., 2021, Guo et al., 6 Aug 2025).