Nonlinear stability of Einstein-de Sitter universes

Abstract: The Einstein-de Sitter universe is the prevailing model used in cosmology to describe the cold dark matter-dominated epoch of the universe. This model is a spatially homogeneous and isotropic spacetime undergoing decelerated expansion, and is linearly unstable under the Einstein-Euler equations with a pressureless fluid equation of state. We show that every initial data set for the Einstein-Euler equations on $\mathbb{T}^3$ with a near-flat metric and positive fluid energy density converges to a flat metric under the Einstein-Euler flow with a polytropic equation of state. This means the metric asymptotes to an Einstein-de Sitter spacetime. In particular, this settles the question of whether the Einstein-de Sitter model can be nonlinearly stable for an appropriate matter model.

• The paper establishes a strong nonlinear stability theorem for Einstein–de Sitter universes under the Einstein–Euler system with a polytropic equation of state (n > 3), proving convergence to EdS-like asymptotics.
• It utilizes a refined analytical approach with sourced wave gauges, conformal rescaling, and high-order L2 energy estimates to control metric and matter nonlinearities over cosmological timescales.
• The results validate pressure-stabilized, matter-dominated models in cosmology and provide a framework for stability analysis in decelerated, non-accelerating universes.

Nonlinear Stability of Einstein–de Sitter Universes: An Expert Analysis

Background and Motivation

The Einstein–de Sitter (EdS) universe occupies a central place in mathematical cosmology as the idealized model for the cold dark matter–dominated epoch. Structurally, the EdS spacetime is spatially flat, homogeneous, isotropic, and exhibits a decelerating cosmic expansion with $a(t) \propto t^{2/3}$, filled with a pressureless dust fluid. While extensively utilized in both theoretical and observational cosmology, EdS is known to be linearly unstable under the Einstein–Euler (dust) system: small scalar perturbations in the dust density can grow unboundedly, as presaged by early Jeans-type analysis and numerical results.

This paper establishes a strong nonlinear stability result for a continuous family of solutions close to EdS, but crucially with a polytropic rather than purely dust equation of state. The result addresses a long-standing gap: the nonlinear stability of matter-dominated cosmological spacetimes. Specifically, the authors prove that for initial data on $T^3$ with a near-flat metric and positive fluid energy density, the coupled Einstein–Euler flow with a polytropic equation of state (with index $n > 3$) is nonlinearly stable—every such initial data set converges to a spatially homogeneous (EdS-like) spacetime, with the metric asymptoting to the EdS form and the density to the expected scaling.

Main Results and Claims

The critical theorem asserts that the set of nonlinearly stable solutions to the Einstein–Euler system with polytropic index $n > 3$ is open, future geodesically complete, and exhibits the same leading-order asymptotics as EdS. Explicitly, all relevant metric and fluid quantities converge to time-dependent backgrounds with precisely computed rates, matching the expansion law and density decay of EdS. Symmetry breaking and localized perturbations dissipate in a homogenizing manner:

• As $t \rightarrow \infty$, the spatial inhomogeneities decay and all metric components and energy densities approach constants after rescaling by the EdS behaviors.
• The solution is future causally geodesically complete.
• The class of stable solutions includes all sufficiently small, near-EdS perturbations with positive (polytropic) pressure and energy density.

Numerically, the decay exponents for perturbations in the system match those previously observed for polytropic fluids in works such as [FajMar:polytropic:26].

Importantly, the proof is fully nonlinear: it controls both metric and matter nonlinearities, and establishes boundedness and specific decay of perturbations over cosmological timescales. Moreover, the authors show that for $n \le 3$, stability fails—fluid homogenization and hence spacetime stability are lost.

Technical Approach and Proof Outline

The analysis proceeds by:

• Formulating the Einstein–Euler system in sourced wave gauge with carefully designed gauge source functions. The damping structure in the gauge is crucial for long-term control.
• Employing conformal rescalings so that the background EdS solution becomes a (perturbed) Minkowski metric, making decay estimates tractable.
• Decomposing dynamics into spatial averages (ODEs for mode zero) and oscillatory components (PDEs for higher modes), leveraging Poincaré inequalities to control the latter.
• Constructing high-order $L^2$ energy functionals with gauge-adapted correction terms, allowing for coercive bulk damping as in previous global nonlinear stability analyses (cf. [ChrKlai:minkowskistab:93], [Ringström]).
• Utilizing a commutator scheme in which vectorfield (material) derivatives act as preferred commutators, exploiting the improved decay rates inherent to the matter sector (fluid equations).
• Identifying and controlling all critical or nonintegrable source terms via gauge choices (notably, choosing the damping parameter $\kappa$ sufficiently large in the gauge source).
• Employing detailed ODE analysis for the spatially averaged sector to prevent secular growth and guarantee decay to the flat background.

The proof handles the strong metric–matter coupling, the absence of spatial dispersion due to the $T^3$ topology, and the fact that the leading decay rates are barely integrable—necessitating optimized correction mechanisms both in the energies and the gauge source.

Connections to Prior Work and Context

Previous stability results in general relativity are restricted to either vacuum models (Minkowski [ChrKlai:minkowskistab:93], Kerr family) or rapidly expanding (accelerated) FLRW spacetimes with positive cosmological constant, for which exponential expansion is the key stabilizing mechanism. For decelerated models (like EdS, with $a(t) \sim t^{2/3}$), linear analyses predict instability for dust, yet with added pressure the situation changes. Recent linear and numerical results (e.g., [Fajetal:decelEuler2:25], [FajMar:polytropic:26]) indicate $T^3$0 polytropes can stabilize the fluid and, consequently, the geometry.

This work realizes and rigorously establishes that insight at the fully nonlinear level.

Theoretical and Practical Implications

Practically, this result retroactively justifies the use of EdS-like models for the large-scale structure and evolution of the actual universe during the matter-dominated epoch—provided one augments dust with a stiff enough equation of state (subgrid pressure, $T^3$1). It reinforces that homogenization (i.e., collapse of perturbations to spatial averages) is achievable in the nonlinear regime for a wide class of physically relevant initial data.

Theoretically, the proof strategy delivers a new blueprint for handling cosmological stability for strongly coupled matter–gravity systems without exponential expansion; it shows how precise gauge fixing, corrected energies, and careful ODE analysis for non-dispersive settings can deliver robust decay.

Implications extend toward:

• Stability of broader classes of matter-dominated cosmological models (potentially those motivated by kinetic theory or more general fluids).
• Pathways to addressing more realistic scenarios with coexisting matter fields, such as coupling dark matter, dark energy, or electromagnetic fields.
• Further analysis of fluid equations at or near critical indices ($T^3$2), where homogenization is marginally attainable.

Limitations and Open Questions

• The stability mechanism requires a polytropic index $T^3$3. For $T^3$4, fluid and metric instabilities prevail.
• The analysis is restricted to compact topologies ($T^3$5); questions remain for open (noncompact) universes.
• The approach assumes initial data close to homogeneous, time-symmetric (spatially flat) backgrounds; extension to genuinely large perturbations or other Bianchi types is nontrivial.

Open mathematical problems include the critical case $T^3$6, coupled matter systems with structure formation (nonlinear clustering), and rigorous matching of stability results across cosmological transitions (radiation/matter/dark energy epochs).

Conclusion

This work furnishes the first nonlinear stability theorem for matter-dominated cosmological models with physically reasonable matter content, under the Einstein–Euler system with stiff ($T^3$7) polytropic fluids. It demonstrates that, beyond the linearly unstable dust case, pressure-stabilized fluids can drive nonlinear convergence to spatial homogeneity and EdS asymptotics. The methodology, combining customized gauge-fixing, conformal normalization, corrected high-order energies, and precise ODE decay control, is likely to influence the broader program of nonlinear stability in general relativistic cosmology, particularly for decelerated (non-accelerating) spacetimes.

References:

• (2607.09326) Louie Bernhardt, David Fajman, Zoe Wyatt, "Nonlinear stability of Einstein–de Sitter universes"
• Supplementary: [Fajetal:decelEuler2:25], [FajMar:polytropic:26], [ChrKlai:minkowskistab:93]