Big bang stability and isotropisation for the Einstein-scalar field equations in the ekpyrotic regime

Abstract: It has been shown that, in spacetime dimensions $n\geq 3$, that the Kasner-scalar field solutions to the Einstein-scalar fields equations with potential $V_0 e^{-s φ}$, where $s<s_c=\sqrt{\frac{8(n-1)}{n-2}}$ and $V_0\in \mathbb{R}$, are nonlinearly stable to the past and terminate at a quiescent big bang singularity over the full range of sub-critical Kasner exponents. In particular, the spatially homogeneous and isotropic solutions, the Friedman-Lemaitre-Robertson-Walker (FLRW) spacetimes, to the Einstein-scalar field equations are stable in this sense for $s<s_c$ and $V_0\in\mathbb{R}$. While perturbations of the sub-critical Kasner-scalar field family of solutions, including the FLRW solutions, are asymptotically velocity term dominated (AVTD) near the big bang, they do not in general isotropise near the big bang singularity. Rather, they remain highly anisotropic, even for small perturbations of the isotropic FLRW solutions. In this article, we establish the stability of FLRW solutions to the Einstein-scalar field equations for the potential parameter values $s>s_c$ and $V_0<0$. Such scalar field potentials are known in the literature as \textit{ekpyrotic}. In particular, we prove that the FLRW solutions to the Einstein-scalar field equations are nonlinearly stable to the past and terminate at a quiescent, crushing AVTD big bang singularity. A distinguishing property of these perturbed spacetimes is that they isotropise towards the big bang.

• The paper demonstrates that in the ekpyrotic regime, FLRW solutions are nonlinearly stable and isotropise near the big bang.
• It employs Fuchsian analytic methods to derive precise asymptotic expansions and energy estimates for curvature decay and scalar field dynamics.
• Results reveal a Ricci-dominated singularity with suppressed anisotropies, challenging traditional views on Kasner behavior.

Nonlinear Stability and Isotropisation at the Big Bang in the Ekpyrotic Einstein-Scalar Field Regime

---

Background: Big Bang Singularities and Scalar Field Dynamics

This paper rigorously investigates the asymptotic behavior and stability of spatially flat FLRW solutions to the Einstein-scalar field equations with exponential potentials in spacetime dimensions $n \geq 3$, focusing on the so-called ekpyrotic regime ($s > s_c$ and $V_0 < 0$). The analysis addresses longstanding questions about the nature of cosmological singularities: specifically, whether generic cosmological solutions can remain stable and isotropic near the big bang under nontrivial matter content, and how the scalar field potential affects this behavior.

The context is framed in terms of the BKL conjecture, which posits that generic spacelike singularities are local and oscillatory. However, prior work shows that the presence of a massless scalar field or stiff fluids tends to suppress such oscillations, leading instead to asymptotically velocity term dominated (AVTD) asymptotics near the big bang. Historically, much of the rigorous work has been limited to the "Kasner regime" ($s < s_c$) or free massless scalar field cases ($V_0 = 0$), where pointwise anisotropy persists even for perturbations of isotropic solutions.

---

Regime Classification and Main Results

The Einstein-scalar field equations with an exponential potential $V(\phi) = V_0 e^{-s \phi}$ admit two qualitatively distinct regimes:

• Kasner regime ($s < s_c$): Solutions are AVTD, but perturbations generically do not isotropise—anisotropies persist near the singularity. The effective equation-of-state parameter asymptotes to $w = 1$ (stiff fluid).
• Ekpyrotic regime ($s > s_c$, $V_0 < 0$): The steepness of the potential fundamentally alters the behavior. The FLRW solution is not only nonlinearly stable under perturbations, but all nearby solutions isotropise as $s > s_c$0 (the big bang), and terminate at a quiescent AVTD singularity. The effective equation-of-state parameter $s > s_c$1 approaches a constant $s > s_c$2, corresponding to "ultrastiff" behavior.

Strong numerical results:

• For $s > s_c$3 (ekpyrotic), the FLRW solution is a past attractor for the spatially flat class.
• The physical shear, spatial curvature, and Weyl components decay relative to the mean (Hubble-normalized) expansion as $s > s_c$4. The Ricci curvature blows up at the same rate as the expansion squared, establishing Ricci-dominated singularity structure.

Bold/contradictory theoretical claims:

• The ekpyrotic regime entirely suppresses anisotropies even for generic small perturbations (contradicting the typical persistence of Kasner-like behavior).
• The scalar field is not dominated by its kinetic energy near the big bang in this regime, unlike all prior stiff-fluid models.

---

Technical Approach and Fuchsian Analysis

To prove full nonlinear stability, the authors use a Fuchsian analytic framework based on a detailed conformal and orthonormal-frame reformulation of the Einstein-scalar equations, building upon techniques from prior work in the massless case. The analysis includes:

• Dynamical-system reduction for FLRW models, with the key autonomous ODE for Hubble-normalized kinetic energy $s > s_c$5.
• Identification and characterization of fixed points: $s > s_c$6 (ekpyrotic), $s > s_c$7 (Kasner).
• Gauge choices and Lagrangian variable synchronization to track the scalar-field time to the big bang.
• Construction of high-order, renormalized, differentiated variables to overcome block-triangular and non-symmetric Fuchsian singular matrices: only through control of high-order spatial derivatives does the system admit energy estimates for the quiescent singularity and isotropisation.
• Precise asymptotic expansions and decay rates for all geometric and curvature quantities, including Ricci and Weyl tensors, showing that isotropisation and Ricci dominance are achieved.

---

Implications and Future Directions

Practical and theoretical significance:

• Establishes that in ekpyrotic cosmologies, the universe will dynamically evolve toward isotropy at the big bang, even under nonlinear, spatially inhomogeneous perturbations.
• Resolves several questions concerning the cosmic no-hair property for contracting universes: the ekpyrotic mechanism provides a mathematically well-defined route to suppressing anisotropic relics, with possible implications for pre-inflationary or cyclic cosmological models.
• The results indicate that AVTD behavior can be decoupled from anisotropy: ultrastiff matter (ekpyrotic scalar fields) achieves AVTD but with isotropisation, contrasting previous cases where AVTD was compatible with persistent anisotropies.

Speculations for future AI and mathematical cosmology:

• The Fuchsian method adopted here, including spatially localized and differentiated techniques to handle block-triangular singular structures, could be adapted for stability analysis of more general Einstein-matter systems.
• Understanding the precise mechanisms for isotropisation in dynamical systems with exponential potentials may be relevant for deep-learning analogs of PDEs, where quiescent boundary behavior and suppression of nonlinear instabilities are desired.
• These results may inform generative modeling of early-universe scenarios, enabling probabilistic models to enforce isotropy constraints dynamically.

---

Conclusion

This work demonstrates the nonlinear past stability and isotropisation of spatially flat FLRW solutions in the ekpyrotic regime for the Einstein-scalar field system. The spectral and energy estimates rigorously quantify the approach to an isotropic and Ricci-dominated big bang singularity, with asymptotic invariants for the scalar field converging to spatially constant values. The mathematical structure developed overcomes previous analytic obstacles and sets a robust template for future stability studies in general relativity and hyperbolic PDEs.

These results redefine the landscape of possible generic singularity behaviors in cosmological models with scalar fields, clarifying the distinction between Kasner and ekpyrotic regimes and providing precise controls over curvature blow-up and isotropisation. The technical strategies and physical conclusions will inform both mathematical cosmology and future analytic approaches for nonlinear wave systems in geometric contexts.