Nonlinear Matter Power Spectrum from relativistic $N$-body Simulations: $\Lambda_{\rm s}$CDM versus $\Lambda$CDM

This paper presents relativistic $N$-body simulations demonstrating that the sign-switching $\Lambda_{\rm s}$CDM model produces a distinct, redshift-dependent enhancement in the nonlinear matter power spectrum at group and poor-cluster scales, offering a unique, falsifiable signature that differentiates it from standard $\Lambda$CDM and potentially aligns with the observed peak in cosmic star formation.

The Gist

The Big Picture: A Cosmic Plot Twist

Imagine the history of our Universe as a long movie. For decades, the standard script (called $\Lambda$CDM) has been the box-office hit. It says the Universe started with a Big Bang, expanded, and is currently being pushed apart by a mysterious force called "Dark Energy" (represented by the Greek letter Lambda, $\Lambda$). This force acts like a constant, unchanging engine that never speeds up or slows down.

However, recent scenes in the movie have started to look a bit glitchy. When astronomers look at the early Universe (the "opening credits"), they get one set of numbers for how fast the Universe is expanding. But when they look at the modern Universe (the "current scene"), they get different numbers. These "plot holes" are known as the Hubble Tension and the $S_8$ Tension.

This paper proposes a radical rewrite of the script. Instead of a constant engine, what if the Dark Energy engine has a gear shifter?

The New Script: The "Sign-Switching" Cosmological Constant ($\Lambda_s$CDM)

The authors, led by Özgür Akarsu and colleagues, suggest that about 10 billion years ago (around redshift $z \approx 2$), the Universe didn't just keep coasting. Instead, the Dark Energy underwent a sudden mirror flip.

• Before the flip (The AdS Phase): The Dark Energy was actually negative. Imagine this as a cosmic "brake" or a gentle gravity that was trying to pull things together, but not enough to stop the expansion.
• The Flip: Suddenly, the sign switched.
• After the flip (The dS Phase): The Dark Energy became positive (like our current standard model), acting as a repulsive force pushing the Universe apart.

Think of it like driving a car up a hill.

1. Standard Model ($\Lambda$CDM): You are driving up a hill with a constant, gentle slope.
2. New Model ($\Lambda_s$CDM): You are driving up a hill, but halfway up, the road suddenly flattens out and then starts sloping downhill (the negative phase), making it easier to go faster. Then, just as you reach the top, the road flips back to an uphill slope (the positive phase), but you've already built up so much speed that you zoom past the finish line faster than the standard model predicted.

The Experiment: Simulating the Universe

To test if this "gear shift" story makes sense, the scientists didn't just do math on paper. They built a virtual universe inside a supercomputer.

They used a tool called gevolution, which is like a high-tech video game engine that follows the strict rules of Einstein's General Relativity. They ran two simulations side-by-side:

1. The Standard Movie: A universe with the normal, constant Dark Energy.
2. The Twist Movie: A universe where Dark Energy flips from negative to positive.

They filled these virtual universes with billions of particles representing matter (stars, gas, dark matter) and watched how they clumped together to form galaxies and clusters over time.

The Discovery: The "Crest" of Power

Here is the most exciting part of their finding.

In the "Twist Movie," the Universe didn't just expand differently; the structure of the Universe changed in a very specific way.

• The "Braking" Effect: During the early "negative" phase, the expansion of the Universe slowed down slightly. In physics, when expansion slows, gravity gets a better chance to pull matter together. It's like a runner slowing down on a track; they have more time to build up momentum. This caused matter to clump together faster than in the standard model.
• The "Boost": When the switch flipped to positive, the Universe started expanding faster again. This tried to pull the clumps apart, but the matter had already built up too much momentum. The extra clumping couldn't be erased.

The Result: The scientists found a distinct "bump" or crest in the data.

• Imagine the Universe is a giant ocean. The standard model predicts waves of a certain size.
• The new model predicts that in a specific size range (the size of galaxy groups and small clusters), the waves are 15% to 25% taller than expected.
• This "bump" appeared right around the time the Dark Energy flipped its sign, and it has drifted to larger scales as the Universe aged.

Why This Matters: Solving the Mystery

Why is this "bump" so important?

1. It Solves the Tensions: The standard model struggles to explain why the Universe is expanding faster today than the early Universe suggests it should. The "gear shift" model naturally explains this by saying, "We had a period of slower expansion that let things speed up later."
2. It's a Fingerprint: You can't fake this result. If you just turned up the volume on the whole universe (making everything slightly denser), you would get a uniform change. But this model creates a specific shape—a localized bump at a specific size. It's like a unique fingerprint left at the scene of the crime.
3. The "Cosmic Noon" Connection: The time when this bump is strongest corresponds to "Cosmic Noon"—the era in the Universe's history when stars were being born at their fastest rate. The authors suggest this extra gravitational "push" might be the reason why star formation peaked so strongly at that specific time.

The Bottom Line

This paper says: "The Universe might have a secret gear shifter."

Instead of Dark Energy being a boring, constant force, it might have flipped from a "pulling" force to a "pushing" force billions of years ago. This flip left a permanent scar on the Universe: a specific pattern of galaxy clusters that are slightly larger and more numerous than our current theories predict.

What's next?
The authors are telling astronomers: "Look at the galaxy clusters and the way light bends around them (weak lensing) at these specific sizes. If you find this 'bump,' we've found the gear shifter. If you don't, the standard script might be right after all."

It's a testable, falsifiable idea that turns a cosmic mystery into a concrete target for future telescopes like Euclid and the James Webb Space Telescope.

Technical Summary

1. Problem Statement

The standard $\Lambda$CDM cosmological model, while successful at describing early-universe physics (CMB, BBN), faces persistent tensions with late-time observations. The most prominent are:

• The $H_0$ Tension: A $>6\sigma$ discrepancy between the Hubble constant inferred from CMB data and local distance ladder measurements.
• The $S_8$ Tension: A discrepancy in the clustering amplitude of matter ($S_8 \equiv \sigma_8 \sqrt{\Omega_m/0.3}$) between CMB predictions and weak lensing surveys.
• The Growth-Index ($\gamma$) Tension: Observations suggest a lower present-day growth rate ($f_0$) than predicted by standard General Relativity (GR) within $\Lambda$CDM.

The paper investigates the $\Lambda_s$CDM model, a phenomenological extension where the cosmological constant ($\Lambda$) undergoes a rapid "sign-switching" transition from negative (Anti-de Sitter, AdS-like) to positive (de Sitter, dS-like) at a specific redshift $z^\dagger \approx 2$. While linear perturbation theory suggests this model can alleviate these tensions by modifying the expansion history, the nonlinear regime (where structure formation becomes complex) had not been rigorously tested. The authors aim to determine if the unique dynamical history of $\Lambda_s$CDM leaves a distinct, falsifiable signature in the nonlinear matter power spectrum that cannot be mimicked by simple parameter rescaling (e.g., changing $\sigma_8$).

2. Methodology

The authors employed relativistic N-body simulations to model structure formation in the nonlinear regime, avoiding the approximations of standard Newtonian codes which may fail in scenarios with rapid background evolution changes.

• Simulation Code: They used gevolution, a particle-mesh code that solves the Einstein equations in the Poisson gauge, ensuring a consistent General Relativity framework.
• Model Setup:
  • $\Lambda_s$CDM: Implemented an "abrupt" limit where the cosmological constant switches sign instantaneously at $z^\dagger$ (modeled as a step function).
  • $\Lambda$CDM: Standard baseline model.
  • Parameters: Two independent best-fit parameter sets were used based on observational data:
    1. Planck-only: CMB data from Planck.
    2. Full Dataset: Combined Planck + BAO + Supernovae (Pantheon+SH0ES) + Weak Lensing (KiDS-1000).
  • Physics: Simulations included baryons and Cold Dark Matter (CDM) as a single particle ensemble ($4160^3$ particles in a $(2080 \text{ Mpc}/h)^3$ box). Massive neutrinos ($\sum m_\nu = 0.06$ eV) were included and evolved linearly.
• Analysis:
  • Calculated the total matter power spectrum $P(k)$ at redshifts $z = 15, 2, 1, 0$.
  • Compared results against linear theory predictions from CLASS (Cosmic Linear Anisotropy Solving System).
  • Analyzed the ratio $P_{\Lambda_s\text{CDM}} / P_{\Lambda\text{CDM}}$ to isolate the specific effects of the sign-switching mechanism, minimizing numerical resolution errors.

3. Key Contributions

• First Nonlinear Test: This is the first study to use full relativistic N-body simulations to test the $\Lambda_s$CDM scenario in the nonlinear regime.
• Identification of a Unique Signature: The authors identified a specific, redshift-dependent "crest" in the power spectrum ratio that serves as a "smoking gun" for the sign-switching mechanism.
• Distinction from $\sigma_8$ Rescaling: They demonstrated that the observed features cannot be reproduced by a scale-independent change in the amplitude of fluctuations (like tuning $\sigma_8$ or $A_s$), proving the effect is dynamical and structural.
• Connection to Cosmic Noon: The paper links the timing of the power spectrum enhancement to the epoch of peak cosmic star formation ("cosmic noon"), suggesting a gravitational prior for observed astrophysical phenomena.

4. Key Results

The simulations revealed a distinct evolutionary pattern in the matter power spectrum ratio ($P_{\Lambda_s\text{CDM}} / P_{\Lambda\text{CDM}}$):

1. Mechanism of Enhancement:

   • AdS Phase ($z > z^\dagger$): The negative cosmological constant reduces the Hubble expansion rate ($H(z)$), thereby reducing Hubble friction ($2H\dot{\delta}$). This dynamically enhances the growth of matter perturbations compared to $\Lambda$CDM.
   • dS Phase ($z < z^\dagger$): After the switch, the expansion rate is higher than in $\Lambda$CDM, increasing friction and suppressing subsequent growth. However, the earlier amplification is preserved, leading to a net excess in power.

2. The "Crest" Feature:

   • A localized peak (crest) appears in the ratio $P_{\Lambda_s\text{CDM}} / P_{\Lambda\text{CDM}}$.
   • Amplitude: The crest reaches 20–25% enhancement near the transition redshift.
   • Location & Drift:
     • At $z \approx 2$ (Planck-only) or $z \approx 1$ (Full dataset), the peak is at wavenumbers $k \simeq 1\text{--}3 \, h \text{ Mpc}^{-1}$.
     • As the universe evolves to $z=0$, the crest migrates to larger physical scales (lower $k$), settling at $k \simeq 0.6\text{--}1.0 \, h \text{ Mpc}^{-1}$.
   • Persistence: By $z=0$, a robust 15–20% uplift remains at $k \simeq 0.6\text{--}1.0 \, h \text{ Mpc}^{-1}$.

3. Physical Scales:

   • The affected wavenumbers correspond to group and poor-cluster environments (masses $\sim 10^{13} \text{--} 10^{14} M_\odot$).
   • This scale is highly relevant for weak lensing, galaxy-galaxy lensing, cluster counts, and the thermal Sunyaev-Zel'dovich (tSZ) effect.

4. Validation:

   • On large scales ($k \lesssim 0.1 \, h \text{ Mpc}^{-1}$), the simulation results perfectly match linear theory (CLASS), validating the code.
   • The nonlinear enhancement is confirmed to be a genuine response to the background dynamics, not a numerical artifact.

5. Significance and Implications

• Falsifiability: The paper provides a concrete, falsifiable target for upcoming surveys (Euclid, LSST, CMB-S4). Detecting a redshift-evolving, localized excess in power at $k \sim 1 \, h \text{ Mpc}^{-1}$ would strongly support the $\Lambda_s$CDM model, whereas a null result would rule out this specific sign-switching scenario.
• Solving Tensions: The results support the hypothesis that $\Lambda_s$CDM can simultaneously resolve the $H_0$, $S_8$, and $\gamma$ tensions without invoking modified gravity, by naturally suppressing the present-day growth rate while maintaining early-universe consistency.
• Astrophysical Connection: The timing of the maximum power enhancement aligns with the "cosmic noon" ($z \sim 1\text{--}2$), the epoch of peak star formation. The model suggests that the enhanced gravitational backdrop during the AdS phase may have modulated the efficiency of halo collapse and star formation during this critical era.
• Future Directions: The authors call for hydrodynamical simulations to include baryonic feedback and further exploration of smooth (non-abrupt) transitions to refine the constraints on the transition redshift $z^\dagger$ and the transition sharpness.

In conclusion, this work establishes that the sign-switching cosmological constant leaves a unique, non-trivial imprint on the nonlinear matter power spectrum, offering a new pathway to test fundamental physics beyond the standard $\Lambda$CDM paradigm.