🔭 astrophysics

Non-linear Cosmological Perturbations for Coupled Dark Energy

This paper derives analytical and numerical one-loop perturbation kernels for a minimal modified gravity model with constant dark energy-dark matter coupling, providing fitting functions to test such theories in future large-scale surveys.

Original authors: Bilal Tüdes, Luca Amendola

Published 2026-04-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Bilal Tüdes, Luca Amendola

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

Imagine the universe as a giant, expanding ocean. For decades, scientists have been trying to understand the currents and waves in this ocean to figure out what's actually driving the water's movement.

In the standard story (called $\Lambda$ CDM), the ocean is pushed by two main forces: Dark Matter (which acts like heavy anchors pulling things together) and Dark Energy (which acts like a mysterious wind pushing everything apart). Usually, we assume these two forces ignore each other completely—they just do their own thing.

But this paper asks a "What if?" question: What if Dark Matter and Dark Energy are actually holding hands? What if they are constantly talking to each other, trading energy and influencing how the universe expands?

The authors, Bilal Tudes and Luca Amendola, are trying to figure out exactly how this "hand-holding" changes the shape of the waves in the cosmic ocean, especially when those waves get a little choppy and messy.

Here is the breakdown of their work, using some everyday analogies:

1. The Problem: The "Smooth" vs. The "Messy"

For a long time, scientists only looked at the universe on a very large scale, where everything looks smooth and calm. It's like looking at the ocean from a satellite; you just see a flat, blue surface. In this "smooth" world, the math is easy.

But as we get better telescopes, we can zoom in. Now we see the choppy waves, the whitecaps, and the turbulence. This is the non-linear regime. It's like looking at the ocean from a boat in a storm. The waves crash into each other, twist, and turn in complex ways.

To predict how these waves behave, scientists use a tool called Perturbation Theory. Think of this like a recipe.

Level 1 (Linear): You just mix flour and water. Easy.
Level 2 & 3 (Non-linear): Now you have to account for how the flour clumps, how the water splashes, and how the bowl shakes. You need a much more complex recipe to predict the final dough.

2. The "Kernels": The Secret Sauce

In this paper, the authors are calculating the "Kernels."

If you imagine the universe's structure as a giant soup, the "kernels" are the secret spices that tell you how the ingredients (matter) mix together when they collide.

In the standard model (where Dark Matter and Dark Energy don't talk), we know exactly how much salt and pepper to add.
In this new model (Coupled Dark Energy), the "spices" are different because the ingredients are interacting. The authors had to derive new recipes for these spices.

They calculated these recipes up to the "third order," which means they looked at very complex interactions where three different waves crash into each other simultaneously.

3. The Two Scenarios They Tested

The authors tested two different ways this "hand-holding" between Dark Matter and Dark Energy could happen:

Scenario A: The "Constant Chat" (Exponential Potential): Imagine Dark Energy and Dark Matter are having a constant, steady conversation. The strength of their connection never changes.
Scenario B: The "Gentle Nudge" (Linear Potential): Imagine the connection is very weak and changes slowly, almost like a gentle nudge. This is interesting because if the nudge is zero, it looks exactly like our standard universe, making it a "minimal" change to our current understanding.

4. The "ϕMDE": The Starting Line

To solve these complex equations, you need a starting point. The authors found a specific moment in the universe's history called the $\phi$ -Matter-Dominated Epoch (ϕMDE).

Think of this as the "starting line" of a race. Even if the runners (Dark Matter and Dark Energy) have different strategies later on, they all start from the same track. The authors calculated the exact position and speed of the runners at this starting line, which allows them to predict how the race will play out later.

5. The Result: A New Toolkit for Future Surveys

The most important part of this paper is that the authors didn't just do the math; they created user-friendly tools (called "fitting functions").

Imagine they built a GPS app for cosmologists.

Before: If you wanted to study a specific type of universe, you had to run a supercomputer simulation for weeks to get the answer.
Now: You can just plug in the numbers (how strong the coupling is, what time in the universe you are looking at), and the "app" gives you the answer instantly with 99% accuracy.

Why Does This Matter?

We are about to launch massive new telescopes (like Euclid and LSST) that will map millions of galaxies. These telescopes will see the "choppy waves" of the universe in incredible detail.

If we use the old "smooth ocean" recipes to analyze this new data, we might get the wrong answer. We might think we found a new law of physics when we just used the wrong recipe.

This paper provides the correct recipe for a universe where Dark Matter and Dark Energy are coupled. By using these new kernels, scientists can:

Test Gravity: Check if Einstein's theory of gravity holds up or if it needs a tweak.
Break Deadlocks: Solve puzzles where different theories look the same until you look at the "choppy waves."
Find the Truth: Determine if Dark Matter and Dark Energy are actually friends (coupled) or strangers (uncoupled).

In short: The authors have written the instruction manual for the next generation of cosmic maps, ensuring we don't get lost when we finally zoom in on the messy, beautiful turbulence of our universe.

1. Problem Statement

Testing gravity on cosmological scales is a primary goal of upcoming large-scale structure (LSS) surveys (e.g., Euclid, DESI, LSST). While the linear regime (scales $>50$ Mpc/h) has been extensively studied, the mildly non-linear regime ( $\sim 10$ Mpc/h) contains crucial information to break cosmological parameter degeneracies.

Most alternative gravity models introduce a non-minimal coupling between a scalar field (dark energy) and matter. The authors focus on the Coupled Dark Energy (CDE) model, the simplest modified gravity scenario where a scalar field $\phi$ is coupled to cold dark matter (CDM) via a constant coupling constant $\beta$ , while baryons remain uncoupled (satisfying local gravity constraints).

The Gap: Previous work on CDE has largely been restricted to linear perturbation theory. To utilize the full power of future surveys, one must extend analysis to the non-linear regime. This requires deriving the perturbation kernels (functions describing mode coupling) for CDE, which differ from the standard $\Lambda$ CDM or Einstein-de Sitter (EdS) models due to the interaction term.

2. Methodology

The authors employ Eulerian Standard Perturbation Theory (SPT) to derive the non-linear evolution of density ( $\delta$ ) and velocity divergence ( $\theta$ ) fields up to the one-loop order (third-order kernels).

Background Dynamics:
- They define a flat FLRW background with interacting dark matter and dark energy.
- The system is rewritten using dimensionless dynamical variables ( $x, y, z, v$ ) representing kinetic energy, potential energy, radiation, and baryons.
- They analyze the autonomous dynamical system to identify critical points, specifically the $\phi$ -Matter-Dominated Epoch ( $\phi$ MDE), which serves as the initial condition for late-time evolution.
- Two scalar field potentials are considered:
  1. Exponential Potential: $V(\phi) = V_0 e^{-\mu\phi}$ (theoretical motivation).
  2. Linear Potential: $V(\phi) \approx V_0 - V_0\mu\phi$ (analytical tractability, approximates a minimal modified gravity model with $\Lambda$ CDM + coupling).
Perturbation Theory Framework:
- The authors adopt the kernel structure derived in Ref. [25], utilizing extended Galilean invariance to constrain the form of the kernels.
- The density contrast $\delta$ and velocity divergence $\theta$ are expanded as convolutions of linear modes weighted by kernels $F_n$ and $G_n$ .
- The kernels are decomposed into basis functions ( $\alpha, \beta, \gamma$ ) with time-dependent coefficients ( $a^{(n)}_i, d^{(n)}_i$ ).
Derivation of Evolution Equations:
- By inserting the perturbative expansions into the continuity and Euler equations (modified by the coupling $\beta$ ), the authors derive a system of coupled differential equations for the time-dependent kernel coefficients.
- These equations (Eqs. 3.35–3.40) depend on the linear growth rate $f$ and the source term $S$ , which are specific to the CDE model.
Solution Strategy:
- Analytical: Exact solutions are derived for the EdS limit and the $\phi$ MDE epoch. Analytical solutions for the background variables ( $x, y$ ) are found for the linear potential case.
- Numerical: For general CDE evolution (exponential potential and late-time linear potential), the differential equations for the coefficients are solved numerically starting from $\phi$ MDE initial conditions.
- Fitting: High-precision fitting functions (polynomials in $\beta, \mu, \eta$ ) are constructed for the numerical results to facilitate their use in data analysis.

3. Key Contributions

Derivation of One-Loop Kernels for CDE: The paper provides the first complete derivation of the second- and third-order perturbation kernels ( $F_2, F_3, G_2, G_3$ ) specifically for a coupled dark energy model with a constant coupling $\beta$ .
Analytical and Numerical Solutions:
- Closed-form analytical solutions for the kernel coefficients in the $\phi$ MDE epoch.
- Analytical background solutions for the linear potential case.
- Numerical integration of the full time evolution for both linear and exponential potentials.
Fitting Functions: The authors provide accurate fitting functions for all kernel coefficients and the growth rate $f$ over the relevant parameter space ( $\beta \in [0, 0.2]$ , $\mu \in [0, 0.2]$ , $\eta \in [-1.4, 0]$ ). These are ready for immediate implementation in power spectrum and bispectrum forecasts.
Initial Conditions: They rigorously establish the $\phi$ MDE epoch as the correct initial condition for numerical integration, ensuring the system starts in a matter-dominated phase consistent with CDE dynamics.

4. Key Results

Kernel Coefficients: The time-dependent coefficients ( $a^{(n)}_i, d^{(n)}_i$ $a_{i}^{(n)}, d_{i}^{(n)}$ ) deviate significantly from the standard EdS values (e.g., $a^{(2)}_1 = 10/7$ $a_{1}^{(2)} = 10/7$ in EdS) as $\beta$ $β$ increases.
- For $\beta \in (0, 0.2)$ , deviations from EdS values can reach up to 2% for $a^{(2)}_1$ , 6% for $d^{(2)}_1$ , and up to 10% for $d^{(3)}_5$ .
Detectability: These deviations are large enough to be detected by forthcoming LSS surveys, implying that non-linear clustering data can place tighter constraints on the coupling parameter $\beta$ and the potential slope $\mu$ than linear data alone.
Linear vs. Exponential Potential: The study confirms that the linear potential is an excellent approximation for the exponential potential in the relevant regime, with kernel coefficient deviations between the two models being $< 3\%$ (mostly $< 1\%$ ).
Growth Rate: The linear growth rate $f$ is modified by the coupling, with analytical solutions showing $f = 1 + 2\beta^2$ in the $\phi$ MDE limit.

5. Significance

Bridge to Observables: This work bridges the gap between theoretical CDE models and observational data. By providing the necessary non-linear kernels, it enables the calculation of the power spectrum and bispectrum in redshift space for CDE models.
Model Discrimination: The derived kernels allow future surveys to distinguish between standard $\Lambda$ CDM, uncoupled dark energy, and coupled dark energy models by analyzing the non-linear clustering of galaxies.
Minimal Modified Gravity: The paper highlights the "minimal" nature of the CDE model (adding only one parameter $\beta$ to $\Lambda$ CDM when $\mu=0$ ), making it a prime target for testing the nature of dark energy and gravity.
Practical Utility: The provided fitting functions (Appendix A) allow cosmologists to immediately incorporate these effects into Monte Carlo Markov Chain (MCMC) analyses and forecast studies without needing to re-solve the differential equations for every parameter set.

In conclusion, Tüdes and Amendola have successfully extended perturbation theory to the non-linear regime for coupled dark energy, providing the essential theoretical tools required to test these models against the high-precision data expected from the next generation of cosmological surveys.