Fundamental Cosmic Anisotropy and its Ramifications II: Perturbations in Bianchi spacetimes, and fixed in the Newtonian gauge

Abstract: The standard cosmological model is challenged by an ever-growing collection of observations, which invites (and stimulates) inquiry into possible additions and/or alterations. One such alteration comes from letting cosmic isotropy -- as demanded by the cosmological principle -- go, whilst maintaining only homogeneity. This study concerns Bianchi models, a class of anisotropic, homogeneous spacetimes, and in particular their perturbations. Knowledge of their properties under perturbations (such as allowed wavemodes) aids in understanding cosmological signatures of such universes, e.g. CMBs, and thus allows for comparsion to observation and the theory of the standard model. This study develops linear perturbation theory of general Bianchi models, by working in a frame such that metric components depend solely on (cosmic) time. Perturbation equations in the Newtonian gauge, but for arbitrary metric, are derived for energy density $ρ$, (relativistic) pressure $p$, momentum density $q$, and anisotropic stress $π$, for the case of scalar and pure tensor perturbations. For the former, the equations for density and pressure are combined to yield the equivalent of the Mukhanov-Sasaki equation for Bianchi models. For a specific choice of metric and fluid flow $u$, the Friedmann equations for Bianchi models are also formulated, as this knowledge is necessary to fully formulate the perturbation equations. Finally, the obtained results are applied to the formulations of density contrasts in an Einstein-de Sitter universe and a Bianchi I universe.

• This study extends linear cosmological perturbation theory to a wide range of general Bianchi spaces.
• The research establishes an isentropic equation for perturbations - HAIPE - that applies to all spatially homogeneous spacetimes, and identifying features that impact hydrodynamic and scalar factors.
• This work validates findings by matching Einstein–de Sitter solutions with Newtonian predictions, emphasizing shear's impact on density contrast.
• follow up questions

Motivation and scope

The standard cosmological model rests on the Cosmological Principle, which postulates both spatial homogeneity and isotropy at the largest scales. A growing body of observational results—potential anisotropy in the Hubble parameter, anisotropy in cosmic acceleration, alignments in quasar and supernova distributions, and local bulk flows—has motivated a systematic re-examination of the isotropy assumption while retaining homogeneity. The natural theoretical framework for such an inquiry is the class of Bianchi spacetimes: homogeneous, generally anisotropic cosmologies whose Killing vector fields generate three-dimensional Lie algebras. This paper, the second in a series by Scholtens, Seri, Waalkens, and van de Weygaert (2604.19706), develops linear cosmological perturbation theory for general Bianchi models, without restricting to Bianchi I or to near-FLRW configurations.

Prior work on perturbations of anisotropic universes has concentrated almost exclusively on Bianchi I, whose metric admits distinct scale factors $a(t)$, $b(t)$, $c(t)$ along orthogonal directions and remains analytically tractable with FLRW-like techniques. Analyses of other Bianchi types typically assume closeness to isotropy or collapse to a single scale factor. The present work removes these restrictions by working in a non-coordinate frame adapted to the homogeneity group.

Non-coordinate frames as the technical foundation

The central methodological device is a frame $\{\mathbf{X}_\mu\}$ satisfying three properties: $\mathbf{X}_0$ is timelike everywhere; $[\mathbf{X}_0,\mathbf{X}_\mu]=0$; and the spatial frame vectors close into a Lie algebra $[\mathbf{X}_i,\mathbf{X}_j]=C^a{}_{ij}\mathbf{X}_a$ with structure constants $C^a{}_{ij}$. In this frame any spatially homogeneous metric takes the separated form

$$\mathbf{g}(x)=g_{\alpha\beta}(t)\,\mathbf{e}^\alpha\mathbf{e}^\beta,$$

with all metric components depending solely on cosmic time. This reduces the Einstein field equations from partial to ordinary differential equations—the same simplification exploited in the orthonormal-frame formalism of Taub, Jantzen, and Ellis–MacCallum. The cost is geometric: the Levi-Civita connection acquires explicit structure-constant terms,

$$\Gamma^\sigma{}_{\mu\nu}=\tfrac{1}{2}g^{\alpha\sigma}\left(2g_{\alpha(\mu,\nu)}-g_{\mu\nu,\alpha}+2C^\beta{}_{\alpha(\mu}g_{\nu)\beta}\right)-\tfrac{1}{2}C^\sigma{}_{\mu\nu},$$

and is symmetric if and only if the frame is coordinated. The Riemann tensor gains an additional term proportional to $b(t)$0. Importantly, the authors verify directly that the Ricci identity—and hence the Raychaudhuri equation—remains valid in non-coordinate frames, which underpins the entire subsequent perturbative analysis.

First-order perturbation framework

Following Bruni et al., perturbations are defined via pullback along a diffeomorphism generated by a vector field transverse to the leaves of a one-parameter family of spacetimes, yielding Lie-derivative-based perturbation tensors on the background. The paper then computes first-order variations of (i) the metric (with careful treatment of index raising, noting that index movement does not commute with perturbation), (ii) the energy-momentum tensor decomposed in full generality as

$b(t)$1

including momentum density $b(t)$2 and anisotropic stress $b(t)$3—both essential once isotropy is dropped—and (iii) the Einstein tensor, computed with the xPand/xAct package and then simplified by hand. Notably, the fluid flow is only assumed normalized ($b(t)$4); geodesy is not assumed, so accelerations $b(t)$5 are carried through all expressions. As a consistency check, the general EMT perturbation equations reduce to the standard Mukhanov results upon substitution of conformally flat FLRW data.

Scalar perturbations and the HAIPE equation

Fixing the Newtonian gauge ($b(t)$6), the metric perturbation is

$b(t)$7

and the authors derive the four projected perturbed Einstein equations for $b(t)$8, $b(t)$9, $c(t)$0, and $c(t)$1 for arbitrary background metric. Under three assumptions—isentropic perturbations ($c(t)$2), no scalar anisotropic stress ($c(t)$3, equivalently the weak-field limit), and a barotropic equation of state $c(t)$4—the density and pressure equations combine into a single master equation, termed the homogeneous-and-anisotropic, isentropic perturbation equation (HAIPE):

$c(t)$5

Here $c(t)$6 is the projected wave operator (distinct from the hypersurface Laplacian $c(t)$7). This is the paper's headline result: it generalizes the Mukhanov–Sasaki equation—which applies only to conformally flat FLRW—to arbitrary spatially homogeneous spacetimes with arbitrary normalized fluid flow. Two structural features deserve emphasis. First, shear enters with coefficient $c(t)$8 and vorticity with $c(t)$9: vorticity is strongly suppressed for radiation-like equations of state ($\{\mathbf{X}_\mu\}$0) but enhanced for stiff matter ($\{\mathbf{X}_\mu\}$1). Second, when the momentum density vanishes—as it does for the diagonal-metric, stationary-flow case treated below—the right-hand side vanishes identically and the HAIPE becomes a damped but undriven wave equation in $\{\mathbf{X}_\mu\}$2.

Tensor perturbations

For pure tensor modes $\{\mathbf{X}_\mu\}$3 with $\{\mathbf{X}_\mu\}$4 and $\{\mathbf{X}_\mu\}$5, the authors derive the corresponding projected perturbation equations and their combined isentropic form. These constitute the gravitational-wave propagation equations in generic Bianchi backgrounds; they simplify considerably when $\{\mathbf{X}_\mu\}$6, recovering the standard FLRW case. The combined tensor equation constrains only one combination of the degrees of freedom in $\{\mathbf{X}_\mu\}$7, so it cannot by itself determine the evolution given initial conditions—a limitation stated explicitly in the paper.

Friedmann equations with structure constants

Specializing to a diagonal metric $\{\mathbf{X}_\mu\}$8 with stationary flow $\{\mathbf{X}_\mu\}$9, the kinematical quantities follow directly: $\mathbf{X}_0$0, $\mathbf{X}_0$1, and vanishing vorticity; the flow is geodesic regardless of structure constants. Using Wald's expressions for the Einstein tensor of such metrics, the Friedmann equations take the form

$\mathbf{X}_0$2

with analogous second and constraint equations. The structure-constant contributions scale as $\mathbf{X}_0$3, playing precisely the role of spatial curvature in open/closed FLRW models—an interpretation confirmed via the Gauss equation, where these terms equal half the intrinsic 3-scalar curvature of the homogeneous hypersurfaces. A notable structural result: nonzero momentum density and anisotropic stress require both nonzero shear (distinct scale factors) and nonzero structure constants (a genuinely non-coordinate frame).

Density contrasts: validation and shear amplification

Two applications test and illustrate the formalism. For the Einstein–de Sitter case ($\mathbf{X}_0$4, $\mathbf{X}_0$5, $\mathbf{X}_0$6, $\mathbf{X}_0$7), the HAIPE reduces to $\mathbf{X}_0$8 with solutions $\mathbf{X}_0$9 and $[\mathbf{X}_0,\mathbf{X}_\mu]=0$0. Naively inserting these into the density contrast expression yields leading behavior $[\mathbf{X}_0,\mathbf{X}_\mu]=0$1 and constancy, whereas Newtonian theory gives $[\mathbf{X}_0,\mathbf{X}_\mu]=0$2. The discrepancy is resolved by observing that the constant-$[\mathbf{X}_0,\mathbf{X}_\mu]=0$3 term enters the density contrast multiplied by the small amplitude of $[\mathbf{X}_0,\mathbf{X}_\mu]=0$4 itself; discarding that subleading term recovers exactly the Newtonian growth rates. This agreement provides the paper's main validation of the derived equations.

For a pressureless Bianchi I universe, expanding the density contrast in powers of $[\mathbf{X}_0,\mathbf{X}_\mu]=0$5 gives

$[\mathbf{X}_0,\mathbf{X}_\mu]=0$6

showing that shear systematically amplifies existing density contrasts of either sign. The authors offer a physical argument: sheared volume elements have larger surface area per unit volume, enhancing energy exchange with the environment, so overdensities accumulate faster. Since $[\mathbf{X}_0,\mathbf{X}_\mu]=0$7 asymptotically for Bianchi I (and more broadly, per Wald's asymptotic results), this shear-driven enhancement is transient—but its transient magnitude could leave observable imprints on structure formation statistics.

Limitations and open questions

Several restrictions bound the applicability of the results. The gauge choice is fixed to Newtonian from the outset because the general theory of gauge transformations in non-coordinate frames "would cloud any useful conclusions"—so gauge-invariant formulations of the HAIPE remain unaddressed. The diagonal-metric specialization assumes stationary, geodesic flow; the authors note that non-stationary yet geodesic flows exist for certain structure constants and would encode peculiar velocities, at the price of losing the identification of fluid-flow derivatives with temporal derivatives. The tensor-sector master equation is a single constraint insufficient to determine the two polarizations. Finally, the EdS validation required neglecting a term proportional to the small amplitude of $[\mathbf{X}_0,\mathbf{X}_\mu]=0$8; a fully quantitative reconciliation with Newtonian growth theory, rather than the qualitative one given here, is not carried out.

Conclusion

This work supplies the missing perturbative machinery for general Bianchi cosmologies: gauge-fixed linear perturbation equations valid for arbitrary homogeneous metrics and normalized flows, a Mukhanov–Sasaki-type master equation (HAIPE) governing isentropic scalar perturbations, gravitational-wave propagation equations, and explicit Friedmann equations exhibiting the curvature-like role of Lie-algebra structure constants. Validation against Einstein–de Sitter growth rates, together with the prediction that shear amplifies density contrasts, establishes the framework's internal consistency and physical content. The announced application to CMB simulations in Bianchi V spacetimes via the Sachs–Wolfe effect will determine whether these formal results translate into discriminating observational signatures of fundamental cosmic anisotropy.