Signatures of Extended Dark Energy Parametrisations in Structure Formation under Background Constraints

This study demonstrates that while extended dark energy models remain consistent with $\Lambda$CDM at the background level, their distinct expansion histories and matter densities induce significant, coherent non-linear signatures in structure formation—such as enhanced small-scale power and altered halo formation epochs—that can be effectively constrained by large-scale structure observables.

The Gist

Imagine the universe as a giant, expanding balloon. For decades, scientists have been trying to figure out exactly how fast this balloon is inflating and what's inside it. The standard theory, called ΛCDM (Lambda Cold Dark Matter), suggests the balloon is being pushed apart by a mysterious, unchanging force called "Dark Energy" (represented by the Greek letter Lambda, Λ).

However, recent measurements have created a bit of a headache. When we look at the early universe (like a baby photo), the expansion speed seems different than when we look at the universe today (like a current photo). This has led scientists to ask: Is Dark Energy actually constant, or is it changing over time?

This paper is a detective story where the authors test four different theories about how Dark Energy behaves, using a mix of real-world data and computer simulations.

The Four Suspects (The Models)

The authors compared the standard theory against three "alternative suspects":

1. The Standard Model (ΛCDM): Dark Energy is a constant, unchanging force. (The "Status Quo").
2. The Constant Wobbler (wCDM): Dark Energy is still constant, but its strength is slightly different from the standard guess.
3. The Shapeshifter (CPL): Dark Energy changes its strength slowly and smoothly over time, like a dimmer switch being turned up or down.
4. The Flexible Sculptor (Chebyshev): This is the most complex one. It treats Dark Energy like a flexible piece of clay that can be molded into almost any shape over time, not just a simple curve.

The Investigation: Two Steps

The authors didn't just guess; they followed a strict two-step process.

Step 1: The Background Check (The "Speedometer")

First, they looked at the "background" of the universe. They used real data from telescopes (like the DESI survey and Planck satellite) to measure the expansion history of the universe. Think of this as checking the speedometer of the universe's car.

• The Result: All four models fit the speedometer data reasonably well. They all look like they could be driving the same car. The "Shapeshifter" and "Flexible Sculptor" models showed some interesting deviations, suggesting Dark Energy might be a bit more active than we thought, but they weren't breaking the rules.

Step 2: The Crash Test (The "N-Body Simulations")

Here is where the paper gets really cool. Just because the cars have the same speedometer doesn't mean they handle the same way on a bumpy road.

The authors took the best-fit numbers from Step 1 and plugged them into a massive computer simulation. They created a virtual universe for each model and watched how gravity pulled matter together to form galaxies and clusters (the "cosmic web").

The Analogy: Imagine four different chefs making a cake. They all use the same recipe for the batter (the early universe), but they use slightly different ovens (the expansion history).

• The Standard Chef makes a perfect, standard cake.
• The Flexible Sculptor Chef makes a cake that rises faster and ends up with a much denser, more massive center.

The Big Discoveries

When they looked at the "cakes" (the simulated universes) at the end of the simulation (today), they found some surprising differences:

1. The "Power" of the Clumps:
   In the "Flexible Sculptor" (Chebyshev) and "Shapeshifter" (CPL) models, the universe formed clumps of matter (galaxies and clusters) faster and earlier than in the standard model. It's as if gravity got a head start. By the time the universe reached its current age, these models had built significantly more massive galaxy clusters.

2. The "Migration" of Mass:
   In the early days (high redshift), the alternative models had more small galaxies. But as time went on, those small galaxies merged together. By the time we reached "today" (redshift 0), the alternative models had fewer small galaxies but huge, massive giants that were twice as common as in the standard model. It's like a city where small houses are being demolished to build massive skyscrapers much faster than usual.

3. The "Universal Shape" of the Clumps:
   Despite the differences in how many galaxies formed and when, the actual shape of the galaxies (how dense they are from the center to the edge) remained surprisingly similar across all models.
   The Metaphor: Imagine four different groups of people building sandcastles. One group builds them faster, and another builds bigger ones. But if you look at the shape of a single sandcastle (the ratio of the wet sand in the center to the dry sand on the outside), they all look almost identical. The "internal structure" of the galaxies is governed by gravity, which works the same way in all these models.

The Verdict

The paper concludes that even if Dark Energy only changes a little bit in the background (the speedometer), those small changes get amplified over billions of years. They act like a butterfly effect, leading to a universe with more massive galaxy clusters and earlier structure formation than we currently expect.

The "Flexible Sculptor" model (Chebyshev) showed the most dramatic effects, suggesting that if Dark Energy is indeed flexible, we should see a lot more giant galaxy clusters in the sky than the standard model predicts.

Why does this matter?
It tells us that looking at the "speed" of the universe isn't enough. We also need to count the galaxies and look at how they are clustered. If future telescopes (like the Vera C. Rubin Observatory) find more massive galaxy clusters than expected, it could be the smoking gun that proves Dark Energy is not a constant, but a dynamic, changing force.

Technical Summary

Here is a detailed technical summary of the paper "Signatures of Extended Dark Energy Parametrisations in Structure Formation under Background Constraints".

1. Problem Statement

The standard $\Lambda$CDM cosmological model faces increasing observational tensions, most notably the discrepancy in the Hubble constant ($H_0$) between local and early-Universe measurements and inconsistencies in the growth rate of large-scale structure (LSS) inferred from weak lensing. While the nature of dark energy (DE) remains uncertain, it could be a cosmological constant ($\Lambda$) or a dynamical field.

The central problem addressed is whether modest deviations in the dark energy equation of state (EoS), $w(z)$, constrained only by background observations (expansion history), propagate into distinct, observable signatures in the non-linear regime of structure formation. Previous studies often focused on linear perturbation theory or specific models; this work aims to connect background constraints directly to non-linear clustering observables (power spectrum, halo mass function, density profiles) using N-body simulations for a broader class of DE models.

2. Methodology

A. Cosmological Models

The study compares four models, all assuming a flat Friedmann–Robertson–Walker geometry:

1. $\Lambda$CDM: Standard model with $w = -1$.
2. $w$CDM: Constant EoS parameter $w_0 \neq -1$.
3. CPL (Chevallier–Polarski–Linder): Dynamical EoS parameterized as $w(z) = w_0 + w_1 \frac{z}{1+z}$.
4. Chebyshev Expansion: A flexible, non-parametric expansion of $w(z)$ using Chebyshev polynomials up to fourth order, matched to $w=-1$ at high redshift ($z > 3.5$) to satisfy CMB constraints.

B. Background Constraints

The models were constrained using a joint Bayesian analysis (MCMC) of four independent observational probes:

• BAO: DESI DR1 data (BGS, LRGs, ELGs, QSOs, Lyman-$\alpha$).
• CMB: Planck 2018 compressed likelihood (acoustic scale $\ell_A$, shift parameter $R$, and $\Omega_b h^2$).
• Cosmic Chronometers (CC): Direct $H(z)$ measurements from passive galaxies.
• Strong Lensing Systems (SLS): 143 systems with elliptical galaxy lenses.

The total chi-squared ($\chi^2_{total}$) was minimized to obtain best-fitting parameters for each model.

C. N-body Simulations

Using the best-fitting background parameters, the authors generated initial conditions and ran dark-matter-only N-body simulations using the RAMSES code.

• Setup: Box size $100 , h^{-1}\text{Mpc}$,$256^3$particles, Adaptive Mesh Refinement (AMR) up to level 18 (resolving$\sim 3.8 , h^{-1}\text{kpc}$).
• Initial Conditions: Generated using the MUSIC code. Crucially, while the primordial power spectrum ($A_s, n_s$) was fixed to Planck values for all models, the linear transfer function and growth factor were model-specific, derived from the best-fitting $\Omega_m h^2$ and $H(z)$.
• Reference: A "Planck-$\Lambda$CDM" simulation was run using standard Planck 2018 parameters for comparison.

3. Key Contributions

• Unified Framework: The paper establishes a physically consistent pipeline linking background observational constraints directly to non-linear structure formation via N-body simulations for multiple DE parametrisations.
• Chebyshev Model Implementation: It implements and tests a flexible Chebyshev expansion for $w(z)$ in a full cosmological simulation context, moving beyond simple constant or linear $w(z)$ approximations.
• Decoupling Effects: It disentangles the effects of the physical matter density ($\Omega_m h^2$) on the transfer function from the effects of the expansion history $H(z)$ on the growth factor.

4. Key Results

A. Background Constraints

• All models remain broadly consistent with $\Lambda$CDM at the background level.
• $w$CDM: Favors a phantom-like EoS ($w_0 \approx -1.22$) and a higher Hubble parameter ($h \approx 0.73$).
• CPL: Shows large uncertainties but is consistent with dynamical evolution ($w_0 \approx -0.86, w_1 \approx -1.47$).
• Chebyshev: The zeroth coefficient ($C_0$) deviates significantly from the $\Lambda$CDM prior, driving the model away from $w=-1$.

B. Matter Power Spectrum ($P(k)$)

• Hierarchy: A clear hierarchy in power spectrum amplitude emerges: Chebyshev > CPL > $w$CDM > $\Lambda$CDM.
• Mechanism: The differences are driven by two factors:
  1. $\Omega_m h^2$: Variations shift the matter-radiation equality epoch. Models with higher $\Omega_m h^2$ (like $w$CDM) suppress small-scale power less, enhancing the linear amplitude.
  2. $H(z)$: Modified expansion histories alter the linear growth factor $D(z)$, amplifying differences at late times.
• Non-linear Evolution: The Chebyshev and CPL models exhibit enhanced small-scale power and earlier transition to the non-linear regime compared to $\Lambda$CDM.

C. Halo Mass Function (HMF)

• High Redshift ($z \gtrsim 2$): CPL and Chebyshev models show an enhanced abundance of low-to-intermediate mass haloes, indicating earlier structure formation.
• Low Redshift ($z \lesssim 1$): The excess mass shifts toward higher masses. The Chebyshev model shows up to twice the abundance of the most massive haloes compared to Planck-$\Lambda$CDM.
• $w$CDM Anomaly: Despite having enhanced power, $w$CDM shows a suppression of high-mass haloes. This is attributed to its higher expansion rate ($H(z)$) and lower $\Omega_m$, which increase the critical density $\rho_c$, reducing the enclosed mass for a fixed virial radius ($R_{200c}$).

D. Halo Density Profiles

• Universality: When expressed in scaled units ($r/R_{200c}$), the internal density profiles of haloes show a high degree of universality across all cosmologies.
• Convergence: Deviations are minimal (few percent level) and confined to the innermost regions ($r/R_{200c} \lesssim 0.2-0.4$).
• Implication: This suggests that internal halo structure is governed primarily by non-linear gravitational dynamics (General Relativity) rather than the specific background DE model, provided the models are not extreme.

5. Significance and Conclusion

• Constraining Power: The study demonstrates that large-scale structure observables (specifically the matter power spectrum and halo mass function) are highly sensitive to background-level variations in dark energy, even when those variations are statistically consistent with current background data.
• Non-linear Signatures: Modest changes in $w(z)$ translate into coherent, amplified signatures in the non-linear regime. The Chebyshev model, due to its flexibility, produces the most distinct deviations, highlighting the potential of future surveys (DESI, Euclid, Rubin) to rule out or constrain complex DE parametrisations.
• Limitations: The authors note that the small simulation volume ($100 , h^{-1}\text{Mpc}$) and single realization limit statistical robustness for the most massive clusters. Future work requires larger volumes and baryonic physics to fully characterize cosmic variance.

In summary, the paper provides a robust numerical framework showing that while dark energy models may look similar in the background expansion, they leave distinct "fingerprints" on the cosmic web, particularly in the abundance and clustering of dark matter haloes.