FLRW Cosmological Model

• The FLRW cosmological model is a solution to Einstein’s equations assuming spatial homogeneity and isotropy with a time-dependent scale factor.
• It employs the Friedmann equations to characterize cosmic dynamics, differentiating universes by curvature, matter content, and dark energy effects.
• Extensions of the model include modified gravity theories and complex topologies, with its predictions tested via CMB, BAO, and numerical simulations.

The Friedmann–Lemaître–Robertson–Walker (FLRW) cosmological model is the foundational solution class in relativistic cosmology describing universes that are spatially homogeneous and isotropic at large scales. Introduced through successive developments by Friedmann, Lemaître, Robertson, and Walker, the FLRW spacetimes provide the geometric and dynamical framework for the concordance cosmological model and most explicit models of cosmic evolution. The metric admits a constant-curvature spatial foliation with a time-dependent scale factor, governed by Einstein’s field equations coupled to perfect-fluid or more general matter sectors. The scope of FLRW includes a rich space of models differentiated by curvature, topology, matter content, and modifications of the gravitational action.

1. Mathematical Structure and Symmetry

The FLRW metric is a unique (up to isometry) Lorentzian metric admitting maximally symmetric spatial hypersurfaces and a geodesic, shear-free, irrotational cosmological flow. In comoving coordinates $(t, \chi, \theta, \phi)$, the metric reads: $$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$ where $a(t)$ is the scale factor and $S_k(\chi)$ is determined by the normalized spatial curvature $k\in\{+1,0,-1\}$: $$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$ The metric can be recast into other coordinate systems (e.g., with proper radial or stereographic coordinates), and alternative canonical forms adapted to arbitrary observers and frames have been constructed. A fully general geometric characterization establishes that, given a foliation by constant-curvature spatial leaves with geodesic, shear-free velocity field, imposing the vanishing of the spatial gradient and Laplacian of the expansion $\theta$ along a single timelike curve forces the spacetime to be FLRW. In general relativity, these conditions correspond to vanishing energy flux and its divergence measured by the cosmological flow, yielding both uniqueness and adaptability to arbitrary observer worldlines (Mars et al., 2023).

2. Dynamical Evolution: Friedmann Equations

Einstein’s equations coupled to a perfect fluid with energy density $\rho(t)$ and isotropic pressure $p(t)$, possibly augmented by a cosmological constant $\Lambda$, reduce under FLRW symmetry to two ordinary differential equations:

• Energy constraint (“first Friedmann equation”):

$$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$0

• Acceleration (“second Friedmann equation”):

$$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$1

With an equation of state $$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$2, the continuity equation

$$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$3

yields $$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$4. The parameter $$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$5 defines the sign of spatial curvature: $$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$6 (spherical), $$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$7 (Euclidean), $$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$8 (hyperbolic). Energy density contributions from radiation, baryons, cold dark matter, and dark energy can each be tracked with separate $$ds^2 = -dt^2 + a(t)^2 \left[ d\chi^2 + S_k^2(\chi)\, (d\theta^2 + \sin^2\theta\,d\phi^2) \right],$$9 and $a(t)$0, and the normalised constraint is expressed in terms of the density parameters,

$a(t)$1

(Öztaş et al., 2021, Buchert et al., 2015).

The dynamical system can also be formulated as Hamiltonian mechanics on a minisuperspace, where $a(t)$2 evolves as a "particle" in an effective potential, enabling both direct integration and qualitative phase-space analysis. Modern methods—including Hamilton–Jacobi group invariants—yield general solutions for $a(t)$3 for any barotropic $a(t)$4 and cosmological constant, and connect conserved quantities with Noether symmetries of the cosmic evolution (Ahmadi-Azar et al., 4 Jan 2026).

3. Spatial Topologies and Classifications

While the FLRW metric is usually written with simply connected spatial slices ($a(t)$5), each case admits multiply connected compactifications by quotienting with a discrete, fixed-point-free group of isometries $a(t)$6:

• Spherical forms: $a(t)$7 (e.g., lens spaces, projective space), countably infinite possibilities (Fagundes, 2011).
• Flat forms: Only six compact orientable Euclidean manifolds exist (including the 3-torus $a(t)$8), classified by Wolf.
• Hyperbolic forms: Infinitely many distinct closed hyperbolic 3-manifolds, catalogued by their fundamental group; the Weeks manifold provides the minimal-volume example.

The physical volume $a(t)$9 depends on both the scale factor and the topological manifold. For $S_k(\chi)$0, Weeks’s SnapPea catalog contains thousands of explicit closed hyperbolic 3-manifolds. Each compactification yields a distinct FLRW cosmological model with identical local geometry and expansion history but differing global topological signatures (Fagundes, 2011).

Observational implications include "multiple image" effects—light from a single object can traverse different topological paths—and "circles-in-the-sky" in the Cosmic Microwave Background, where matched circles correspond to the intersection of the last-scattering surface with topological images. CMB searches for such signatures currently limit the minimal fundamental domain scale but cannot yet exclude large classes of compact topologies.

4. Qualitative Cosmological Behavior and Parameter Space

The space of FLRW models is classified according to curvature, matter content, and cosmological constant. Qualitative mechanical analogy assigns the Friedmann equation the structure of a unit-mass particle evolving in a potential $S_k(\chi)$1, with total energy set by the curvature parameter. The system supports a taxonomy of behaviors:

• $S_k(\chi)$2, $S_k(\chi)$3: Closed, recollapsing universes (big bang $S_k(\chi)$4 maximum $S_k(\chi)$5 $S_k(\chi)$6 big crunch).
• $S_k(\chi)$7, $S_k(\chi)$8: Einstein–de Sitter, eternally expanding, $S_k(\chi)$9 for dust.
• $k\in\{+1,0,-1\}$0, $k\in\{+1,0,-1\}$1: Open, forever-expanding models (Milne, $k\in\{+1,0,-1\}$2).
• $k\in\{+1,0,-1\}$3: Asymptotic de Sitter acceleration; late-time behavior transitions to $k\in\{+1,0,-1\}$4.
• Phantom energy ($k\in\{+1,0,-1\}$5): Super-accelerated expansion, big-rip singularities.
• Stiff matter and mixed fluids: Early time $k\in\{+1,0,-1\}$6, $k\in\{+1,0,-1\}$7, transitioning to dust ($k\in\{+1,0,-1\}$8) and/or vacuum domination ($k\in\{+1,0,-1\}$9) (Dariescu et al., 2017, Sonego et al., 2011).

Dynamical systems analysis reveals late-time attractors: for open models with standard fluids, the Milne (curvature-dominated) universe is a generic asymptote unless the equation-of-state is sufficiently soft ($$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$0) (Leon et al., 2021, Kohli et al., 2016). Near Minkowski spacetime, the full Einstein-perfect-fluid FLRW system reduces to a degenerate Bogdanov–Takens normal form, permitting rigorous bifurcation analysis and clarifying stability boundaries between expanding and contracting regimes (Kohli et al., 2016).

5. Observational Foundations, Limitations, and Extensions

The FLRW model forms the background for extraction of cosmological parameters from high-precision data:

• CMB angular power spectrum: Constraints on $$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$1, $$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$2, $$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$3 to percent level (Buchert et al., 2015).
• Baryon Acoustic Oscillations and SNe Ia: Standard rulers and candles, evidence for acceleration ($$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$4).
• Hubble tension and anomalies: Persistent discrepancy between local and CMB-derived $$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$5, as well as hints for quasi-Euclidean universes with very low matter content ($$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$6, $$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$7), challenging the conventional $$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$8 paradigm (Öztaş et al., 2021).

Averaging inhomogeneities (the "backreaction problem") is addressed by the Buchert formalism, which adds effective terms to the Friedmann equations representing averaged variance and shear but finds corrections at $$S_{k}(\chi)= \begin{cases} \sin \chi, & k=+1 \ \chi, & k=0 \ \sinh\chi, & k=-1 \end{cases}$$9 insufficient to mimic dark energy (Fleury, 1 Sep 2025). Light propagation effects—weak lensing, "empty-beam" corrections, inhomogeneity-induced biases—are subdominant in the mean for SNe or BAO distance-redshift relations, though individual lines of sight may experience percent-level deviations.

Not all departures from the strict FLRW framework are excluded by current observations; anomalies in the CMB (e.g., power asymmetry), large-scale structure, and environment-dependent BAO features point to the necessity of broader tests of homogeneity and isotropy (Buchert et al., 2015, Fleury, 1 Sep 2025).

6. Extensions: Modified Gravity and Mathematical Correspondences

Beyond general relativity, the FLRW metric and its dynamical reduction accommodate broad classes of modified gravity:

• $\theta$0 gravity: Modified Friedmann equations with higher-curvature terms, reinterpretation as effective fluids. $\theta$1 models may unify inflation and late-time acceleration, or replicate $\theta$2CDM evolution (Sáez-Gómez, 2011).
• Gauss–Bonnet and Hořava–Lifshitz gravities: The former introduces higher-order curvature invariants; the latter breaks Lorentz symmetry for improved UV behavior, alters the cosmic evolution equations, and can accommodate hybrid expansion histories, subject to experimental constraints.
• Scalar-tensor reconstruction: Any desirable $\theta$3 can be engineered by reconstructing suitable scalar potentials or $\theta$4 functions, mapping the problem to a linear ODE and promoting a geometric-driven approach to cosmic acceleration.

Additionally, mathematical analogies exist between FLRW cosmology (in arbitrary dimension) and dynamical systems such as the time-dependent, harmonically trapped Bose–Einstein condensate, where the width of the BEC evolution under the Gross–Pitaevskii equation corresponds in detail to curvature-normalized FLRW expansion (D'Ambroise et al., 2010). The mapping hinges on the Ermakov–Milne–Pinney equation, with curvature and cosmological constant encoded in conserved quantities and additive shifts.

7. Numerical Methods and Computational Implementation

Numerically evolving the FLRW system is standard for both pedagogical and research applications. Implementations with Planck 2018 parameters and explicit Euler integration schemes reproduce all major expansion regimes (radiation-, matter-, dark-energy-dominated), permit exploration with varying equations of state, and are modular for including further physics such as dynamical dark energy or topological compactification (Marasca et al., 27 Nov 2025). The discretized acceleration equation (with $\theta$5 and $\theta$6 as state variables) captures qualitative transitions and can be upgraded to higher-order methods or generalized to include curvature and more exotic components.

• "An Infinite Number of Closed FLRW Universes for Any Value of the Spatial Curvature" (Fagundes, 2011)
• "Physical Basis for the Symmetries in the Friedmann-Robertson-Walker Metric" (Melia, 2016)
• "Exact solutions of the FLRW cosmological model via invariants of the Hamilton-Jacobi method" (Ahmadi-Azar et al., 4 Jan 2026)
• "Averaging Generalized Scalar Field Cosmologies I: Locally Rotationally Symmetric Bianchi III and open FLRW models" (Leon et al., 2021)
• "A Degenerate Bogdanov-Takens Normal Form for FLRW Cosmologies" (Kohli et al., 2016)
• "Numerical implementation of flat FLRW models of cosmic expansion with Planck 2018 cosmological parameters" (Marasca et al., 27 Nov 2025)
• "Elements of cosmology beyond FLRW" (Fleury, 1 Sep 2025)
• "New characterization of Robertson-Walker geometries involving a single timelike curve" (Mars et al., 2023)
• "Re-evaluation of $\theta$7 of the normalised Friedmann-Lemaître-Robertson-Walker model" (Öztaş et al., 2021)
• "Observational Challenges for the Standard FLRW Model" (Buchert et al., 2015)
• "Qualitative study of perfect-fluid Friedmann-Lemaître-Robertson-Walker models with a cosmological constant" (Sonego et al., 2011)
• "A dynamic correspondence between Bose-Einstein condensates and FLRW and Bianchi I cosmology" (D'Ambroise et al., 2010)
• "Spatially-Flat Robertson-Walker models with combined CDM and stiff matter sources and the corresponding thermodynamics" (Dariescu et al., 2017)
• "On Friedmann-Lemaître-Robertson-Walker cosmologies in non-standard gravity" (Sáez-Gómez, 2011)