⚛️ general relativity

Sub-Horizon and Quasi-Static Approximations in Interacting Dark Sector Models Involving Scalar Field Dark Energy

This paper demonstrates that while the sub-horizon and quasi-static approximations are generally accurate for non-interacting dark energy models, their validity significantly deteriorates in interacting scalar-field scenarios—particularly on large scales and with strong coupling—necessitating full relativistic treatments for precision cosmological simulations and observations.

Original authors: Priyanka Saha, Sharvari Nadkarni-Ghosh, Kaushik bhattacharya

Published 2026-09-28

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Priyanka Saha, Sharvari Nadkarni-Ghosh, Kaushik bhattacharya

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The universe is not a static stage but a dynamic, expanding fabric where invisible forces shape the distribution of matter. For decades, the prevailing model of cosmology has relied on a simple picture: a universe filled with ordinary matter, invisible dark matter, and a mysterious dark energy that pushes the cosmos apart. To understand how galaxies form and cluster, scientists must track how tiny ripples in this cosmic fluid grow over time. In the standard model, these ripples behave predictably, allowing researchers to use simplified rules that treat gravity much like it works on Earth, ignoring the complex, relativistic effects of Einstein's theory. These rules work well when looking at small patches of the sky, far from the edge of the observable universe. However, as new telescopes prepare to map the entire sky with unprecedented precision, scientists are asking whether these simplified rules still hold true when the universe behaves in more complicated ways, specifically when dark matter and dark energy might be talking to each other.

A team of researchers at the Indian Institute of Technology, Kanpur, has investigated this question by studying models where dark matter and dark energy interact directly. In the standard view, these two invisible components evolve independently, but in these alternative theories, they exchange energy and momentum, creating a more complex dance of forces. The researchers wanted to know if the simplified "sub-horizon" approximation—a method that assumes we are looking at scales much smaller than the entire universe and that changes happen slowly enough to ignore certain time-dependent effects—could still accurately describe this interaction. They built a complete, rigorous mathematical description of how these interacting fields evolve and then compared it against the results produced by the simplified methods. Their goal was to see if the shortcuts used in massive computer simulations, which model the formation of cosmic structures, would lead to significant errors when applied to these interacting models.

The study focused on two specific types of simplified approaches. The first, known as the sub-horizon approximation, assumes that the ripples in the universe are so small compared to the size of the cosmos that the expansion of space can be largely ignored for the calculation of gravity. The second, the quasi-static approximation, assumes that the gravitational potential changes so slowly that it can be treated as fixed in time. In the standard model of the universe, where dark energy is a constant, these shortcuts work perfectly well on small scales. However, the researchers found that when dark matter and dark energy interact, the situation changes dramatically. They discovered that the simplified methods fail to capture the full physics of the interaction. Even when looking at scales well within the horizon, the simplified equations miss crucial terms that arise from the exchange of energy between the dark sectors.

The results showed that the accuracy of these approximations depends heavily on the strength of the interaction. For weak interactions, the simplified methods still provide a reasonably good picture, with errors remaining at the level of a few percent. But as the interaction between the dark components grows stronger, the errors in the simplified models increase rapidly. The researchers found that for stronger interactions, the simplified approach could deviate from the true physical evolution by several percent, even on scales that are typically considered safe for these approximations. This is a significant finding because a difference of a few percent is not negligible in modern cosmology; it is comparable to the precision of the data coming from next-generation surveys. The study also revealed that the "quasi-static" approach, which keeps more of the relativistic details, generally performed better than the sub-horizon approximation, staying closer to the true evolution across the range of scales tested.

Perhaps the most critical insight from this work is that the validity of these shortcuts cannot be assumed just because a scale is small. In the past, scientists believed that if they were looking at a region much smaller than the observable universe, the simplified Newtonian rules would always apply. This paper demonstrates that in interacting dark energy models, the presence of the interaction itself introduces new time scales and dynamics that break this assumption. The researchers found that the error introduced by using the simplified equations grows as the interaction strength increases and as one looks at larger scales within the universe. This means that for models where dark matter and dark energy talk to each other, the standard "Newtonian" equations used in many large-scale computer simulations are not fully recovering the correct physics.

The implications for future research are substantial. Upcoming surveys and massive hydrodynamical simulations, which aim to map billions of galaxies and test the nature of dark energy, often rely on these simplified equations to save computational time. If the universe follows an interacting dark energy model, these simulations could be introducing systematic errors that mimic or hide the very signals scientists are trying to detect. The researchers concluded that before employing these models in large-volume simulations or making precision comparisons with observational data, scientists must benchmark the approximations against the full, complex equations for the specific interaction strength they are studying. There is no universal rule that guarantees the simplified methods will work; instead, their validity must be checked for each specific scenario. This work serves as a necessary warning that as our telescopes become more powerful and our simulations more ambitious, we cannot rely on old shortcuts without verifying that they still hold true in a universe where the dark sector might be far more connected than we previously thought.

Technical Summary: Sub-Horizon and Quasi-Static Approximations in Interacting Dark Sector Models

Problem Statement
Upcoming Large Scale Structure (LSS) surveys (e.g., Euclid, DESI, Rubin Observatory) aim to map matter distributions on scales approaching one gigaparsec to constrain cosmological parameters. These analyses rely on theoretical predictions derived from $N$ -body codes or perturbative techniques. In standard $\Lambda$ CDM and non-interacting scalar-field models, the Sub-Horizon Approximation (SHA) and Quasi-Static Approximation (QSA) are widely used to simplify the relativistic Einstein equations into Newtonian fluid equations (continuity, Euler, and Poisson) on scales well within the horizon ( $k \gg H$ ).

However, in interacting dark sector models where dark matter couples to a dynamical scalar field (quintessence or phantom), additional characteristic timescales and dynamical degrees of freedom are introduced. The validity of the SHA and QSA in these interacting scenarios is not guaranteed, as the coupling may introduce rapid oscillations or time-dependent terms that violate the assumptions underlying these approximations. While these approximations are computationally efficient, their accuracy in interacting models has received comparatively little attention, raising concerns about their reliability for precision cosmology and large-volume hydrodynamical simulations.

Methodology
The authors investigate the validity of the SHA and QSA by deriving and comparing three sets of linear perturbation equations for interacting scalar-field dark energy models:

Full Relativistic System: The complete set of linearized Einstein and conservation equations without approximations.
Sub-Horizon Approximation (SHA): Derived by taking the limit $k \gg H$ (neglecting terms suppressed by $H/k$ ) while retaining time derivatives.
Quasi-Static Approximation (QSA): Derived by neglecting time derivatives of the metric potentials ( $\Phi' \ll H\Phi$ ) while retaining spatial gradients.

The study considers three cosmological scenarios:

Interacting Scalar Field: Dark matter couples to a scalar field via an interaction term $\rho_{int} = \gamma \rho_m e^{-\beta \phi}$ .
Non-Interacting Scalar Field: The same scalar field model but with the interaction parameter $\gamma = 0$ .
$\Lambda$ CDM: The standard model with a cosmological constant.

The authors utilize quintessence ( $\epsilon_d = 1$ ) and phantom ( $\epsilon_d = -1$ ) scalar fields. Background evolution is constrained to match $\Lambda$ CDM within 1% to isolate perturbative effects. Initial conditions for perturbations are extracted from the Boltzmann code CLASS. The authors numerically integrate the full, SHA, and QSA equations using an adaptive Runge-Kutta method (RK45) over a range of comoving wavenumbers ( $k \in [10^{-4}, 10^{-1}] \text{ Mpc}^{-1}$ ) and interaction strengths ( $\gamma \in [0, 0.09]$ ).

Accuracy is quantified by comparing the matter growth factor ( $D$ ), growth rate ( $f$ ), and linear matter power spectrum ( $P_m$ ) against the full relativistic solution. Additionally, the authors evaluate whether the SHA can capture the physical differences between interacting models and $\Lambda$ CDM by comparing "full-full" deviations against "SHA-SHA" deviations.

Key Contributions and Results

Breakdown of Newtonian Recovery in Interacting Models:
In $\Lambda$ CDM and non-interacting scalar-field models, the SHA successfully reduces the relativistic equations to the standard Newtonian fluid equations. However, the authors find that in interacting models, applying the SHA does not fully recover the standard Newtonian equations. Additional interaction-dependent terms remain in the evolution equations, meaning the standard Newtonian treatment is insufficient even on sub-horizon scales for these models.
Accuracy Degradation with Interaction and Scale:
- Non-Interacting/Weakly Interacting: The SHA reproduces the full evolution with percent-level accuracy on sufficiently sub-horizon scales.
- Strong Interactions: The accuracy of the SHA progressively deteriorates as the interaction strength ( $\gamma$ ) increases.
- Scale Dependence: The error increases towards larger scales (smaller $k$ ). For strong interactions, the SHA error exceeds 5% even for modes well inside the horizon ( $k \gg H$ ).
- QSA Performance: The QSA generally provides a closer description of the full relativistic evolution than the SHA across the scales considered, as it retains more of the relativistic structure (specifically, it does not expand in powers of $H/k$ as aggressively).
Impact on Model Discrimination:
The authors demonstrate that the difference between comparing a model to $\Lambda$ CDM using full equations versus using SHA equations becomes increasingly significant for stronger interactions and smaller wavenumbers. This indicates that the SHA fails to capture the full extent of the physical differences introduced by the interaction. The "mismatch" ( $D_{mismatch}$ ) between full and SHA predictions grows with $\gamma$ , suggesting that using SHA in interacting models may lead to incorrect conclusions about the viability of alternative dark energy models.
Relevance to Simulations and Surveys:
The errors introduced by the SHA can reach several percent on scales accessible to large-volume hydrodynamical simulations (e.g., the FLAMINGO project, with fundamental modes $k \sim 2.24 \times 10^{-3} \text{ Mpc}^{-1}$ ). Since modern surveys (like DESI) utilize full-shape clustering information including large-scale modes, a scale-dependent theoretical error of a few percent can become comparable to the statistical precision of the data, potentially biasing parameter estimation.

Significance and Claims
The paper argues that the validity of the SHA cannot be assumed solely based on the sub-horizon condition ( $k \gg H$ ) when dealing with interacting dark energy. The presence of interaction introduces a new scale and time dependence that invalidates the standard Newtonian reduction.

The authors conclude that:

The SHA is reliable only for sufficiently weak interactions and large wavenumbers.
For stronger interactions or larger simulation volumes, the approximation must be benchmarked against the full relativistic equations before being employed in $N$ -body or hydrodynamical simulations.
The QSA offers a potentially more accurate intermediate description than SHA for these models but also requires model-specific validation.
As simulations and surveys move to larger volumes and higher precision, approximation errors that were previously negligible become dominant theoretical systematics. Therefore, the range of validity for these approximations must be explicitly established relative to the specific interaction parameters and the scales of the intended application.

The work does not propose new observational strategies or future experimental setups but provides a necessary theoretical benchmark to ensure the consistency of current and future large-scale structure analyses involving interacting dark sector models.

More like this