Interacting Scalar Fields as Dark Energy and Dark Matter in Einstein scalar Gauss Bonnet Gravity

This paper investigates two models of interacting scalar fields in Einstein scalar Gauss-Bonnet gravity as unified dark energy and dark matter candidates, demonstrating that while they align with current supernova data and the $\Lambda$CDM model, they show a statistically significant preference over $\Lambda$CDM when incorporating future Roman Space Telescope mock data at higher redshifts.

The Gist

Imagine the universe as a giant, expanding balloon. For decades, scientists have been trying to figure out what's inside the balloon and what's pushing it to inflate faster and faster.

We know about the "stuff" we can see (stars, planets, you, me), but that's only about 5% of the total. The rest is a mysterious 95% split into two invisible giants:

1. Dark Matter: The invisible glue that holds galaxies together.
2. Dark Energy: The mysterious force pushing the universe apart.

This paper is like a detective story where the authors try to solve the mystery of these two giants by treating them not as invisible fluids, but as vibrating strings of energy (scalar fields) that are talking to each other.

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

1. The Setup: Two Dancers on a Stage

The authors propose a new dance floor (a theory of gravity called Einstein-scalar-Gauss-Bonnet gravity).

• Dancer A (Dark Energy): A slow-moving, energetic dancer who wants to push the balloon apart.
• Dancer B (Dark Matter): A heavy, fast-vibrating dancer who acts like a heavy weight, holding things together.
• The Interaction: Usually, scientists think these two dancers ignore each other. This paper suggests they are holding hands and exchanging energy. When one gets tired, the other gets a boost.

2. The "Speed Limit" Rule (The Gravitational Wave Constraint)

In this new dance theory, there's a rule about how fast "ripples" in the fabric of space (gravitational waves) can travel.

• The Problem: Some versions of this theory predicted that these ripples would travel at a different speed than light.
• The Reality Check: We just observed a neutron star collision, and we know for a fact that gravitational waves travel at the exact same speed as light.
• The Fix: The authors had to tweak their theory so that the "dance moves" (the coupling function) don't change the speed of these ripples. This forced them to simplify the rules, making the model cleaner and more universal.

3. The Two Models: Different Dance Styles

The authors tested two different ways the dancers could hold hands (two different interaction formulas):

• Model I (The Exponential Hug): The strength of their connection grows or shrinks very quickly, like a hug that gets tighter or looser exponentially.
• Model II (The Power-Law Grip): The connection strength follows a steady, mathematical curve, like a spring that stretches at a predictable rate.

They ran computer simulations to see if these dances would lead to a stable universe or if the dancers would trip and fall (creating a universe that collapses or explodes).

4. The Results: A Perfect Mimicry

Here is the surprising part: Both models work perfectly well.

• The "Chameleon" Effect: When they looked at the universe today (low redshift), these complex, vibrating-string models looked almost exactly like the standard, boring model everyone uses (called $\Lambda$CDM). It's like a high-tech robot wearing a t-shirt and jeans; from a distance, it looks just like a human.
• Why? The "extra" gravity effects (from the Gauss-Bonnet term) get so small today that they are practically invisible. The universe behaves normally.

5. The Twist: The Future Telescope (Roman Space Telescope)

If the models look the same today, how do we tell them apart?

• The Time Machine: The authors looked at the "past" (high redshift) and the "future" (using data from the upcoming Roman Space Telescope).
• The Reveal: At very high distances (looking far back in time), the two models start to drift away from the standard model.
  • Think of it like two cars driving side-by-side on a highway. They look identical at 60 mph. But if you check their GPS logs from 100 miles ago, you realize one took a slightly different route.
• The Verdict: The standard model ($\Lambda$CDM) is the "safe" choice. But the new models are statistically preferred when you include the future Roman telescope data. This suggests that if we look far enough back, we might finally see the "vibrating strings" doing their thing.

6. Why This Matters

• Solving the "Hubble Tension": There is a big argument in physics about how fast the universe is expanding (the Hubble constant). Different measurements give different answers. These models don't fully fix that argument yet, but they offer a new, physically consistent way to think about the problem.
• Stability: The authors proved that these models are stable. The universe won't collapse or tear apart because of these interactions.
• Particle Physics Connection: Instead of just guessing how dark matter and energy interact, they built the theory from the ground up using principles from particle physics (like how particles get mass). This makes the theory feel more "real" and less like a mathematical trick.

The Bottom Line

The authors have built a sophisticated, mathematically sound theory where Dark Matter and Dark Energy are two interacting fields of energy.

• Today: They look exactly like the standard model we already know.
• Tomorrow (with new telescopes): They might reveal a hidden layer of complexity that the standard model misses.

It's a bit like discovering that the "ghost" haunting your house is actually just a very complex, vibrating piece of furniture that only makes noise when the wind blows from a specific direction. We can't hear it right now, but with a better microphone (the Roman Telescope), we might finally hear the music.

Technical Summary

1. Problem Statement

The standard cosmological model ($\Lambda$CDM) faces significant theoretical challenges, including the nature of dark energy (DE) and dark matter (DM), the Hubble tension ($H_0$ discrepancy), and the theoretical limitations of a cosmological constant. While interacting dark sector models have been proposed to alleviate these tensions, many rely on fluid approximations where interactions are introduced phenomenologically via continuity equations without a fundamental Lagrangian. This often leads to ambiguities in the interaction term and potential stability issues (ghosts/gradient instabilities).

This paper addresses these issues by:

1. Proposing a fundamental particle physics framework where both DE and DM are described by scalar fields ($\phi$ and $\psi$) interacting via a Lagrangian potential.
2. Embedding this interaction within Einstein-scalar-Gauss-Bonnet (EsGB) gravity, a modification of General Relativity motivated by string/M-theory.
3. Investigating whether such a system can reproduce the observed cosmic acceleration and structure formation while satisfying strict observational constraints, particularly the speed of gravitational waves ($c_T \approx c$).

2. Methodology

Theoretical Framework

• Action: The authors consider an action where two scalar fields ($\phi$ for DE, $\psi$ for DM) are coupled to the Gauss-Bonnet (GB) invariant $\mathcal{G}$ and interact via a potential $W(\phi, \psi)$.
  • DE Field ($\phi$): Modeled as a quintessence field with an exponential potential $V(\phi)$. It couples non-minimally to the GB term via a function $f(\phi)$.
  • DM Field ($\psi$): Modeled as a coherent scalar field with a quadratic potential $U(\psi) = \frac{1}{2}m^2\psi^2$. In the limit $m \gg H$, this mimics pressureless cold dark matter (CDM).
  • Interaction: A non-gravitational interaction term $W(\phi, \psi) = \frac{1}{2}\xi(\phi)\psi^2$ is introduced, effectively modifying the DM mass to $m_{\text{eff}}^2 = m^2 + \xi(\phi)$.
• Gravitational Wave Constraint: The GB coupling $f(\phi)$ affects the speed of gravitational waves ($c_T$). To satisfy the tight observational bound $|c_T^2 - 1| < 10^{-16}$ (from GW170817), the authors impose $c_T = 1$. This constraint forces the time derivative of the coupling function to scale with the scale factor ($\dot{f} \propto a$), rendering the analysis model-independent regarding the specific form of $f$.

Dynamical System Analysis

• The authors formulate the background equations of motion as an autonomous system of first-order differential equations using dimensionless variables ($x, y, z, \Omega_{\psi}, \Omega_{GB}$, etc.).
• Fluid Approximation: To avoid computationally expensive numerical integration of rapid oscillations of the DM field, they employ a WKB approximation to derive a fluid-equivalent continuity equation for $\psi$.
• Two Interaction Models:
  • Model I: Exponential interaction $\xi(\phi) = M \exp\left(\frac{\lambda \kappa^2 \phi^2}{2}\right)^{n_c}$.
  • Model II: Power-law interaction $\xi(\phi) = M (\kappa \phi)^{n_c}$.
• Stability: The stability of the system is analyzed numerically by evolving the system from the deep radiation era ($N \approx -20$) to the future ($N > 0$) to ensure the existence of a stable de Sitter attractor and the absence of ghost/gradient instabilities.

Data Analysis

• Datasets: The models are constrained using a combination of:
  • CC: Cosmic Chronometers (Hubble parameter).
  • PP: Pantheon+ Type Ia Supernovae.
  • DES: Improved Dark Energy Survey (DES-Dovekie) Type Ia Supernovae.
  • DESI: Baryon Acoustic Oscillations (BAO) from DESI Release II.
  • Planck: Compressed CMB data (distance priors).
  • Roman: Mock data from the upcoming Roman Space Telescope (high-redshift SNe).
• Statistical Tools: Markov Chain Monte Carlo (MCMC) sampling using emcee and GetDist. Model comparison is performed using Akaike (AIC) and Bayesian (BIC) Information Criteria.

3. Key Contributions

1. Fundamental Interaction Lagrangian: Unlike fluid-based interacting models, this work derives the interaction directly from an action principle, ensuring consistent equations of motion and avoiding ad-hoc stability fixes.
2. Model-Independent GB Analysis: By enforcing $c_T=1$, the authors simplify the GB coupling function, allowing the study of the interaction dynamics without being bogged down by the specific functional form of the GB coupling, while remaining consistent with GW observations.
3. Initial Condition Strategy: The paper provides a robust methodology for setting initial conditions in dynamical systems. Instead of using a computationally expensive "shooting method" to match early-universe conditions, they sample variables at the present epoch ($N=0$) and demonstrate that for attractor solutions, this reproduces the correct past evolution.
4. Future Survey Projections: The inclusion of Roman Space Telescope mock data allows the authors to project how future high-redshift observations can distinguish between these interacting models and $\Lambda$CDM.

4. Results

• Dynamical Stability: Both Model I and Model II exhibit a stable de Sitter attractor in the late-time universe for specific ranges of parameters ($\lambda, n_c$). The systems successfully transition from radiation domination to matter domination and finally to acceleration.
• Cosmological Evolution:
  • The interacting models closely mimic the $\Lambda$CDM expansion history ($H(z)$) at low redshifts.
  • The GB contribution ($\Omega_{GB}$) is tightly constrained to be negligible ($\lesssim 10^{-25}$) at low redshifts, effectively recovering General Relativity today, but may become significant at high redshifts.
  • The equation of state ($w_{\text{eff}}$) remains in the accelerating regime ($-1 < w < -0.3$) and does not cross the phantom divide ($w < -1$), ensuring stability.
• Parameter Constraints:
  • Hubble Constant: The models yield $H_0 \approx 68.2 - 68.3$ km/s/Mpc for Pantheon+ and DES data, slightly lower than SH0ES but consistent with Planck. For Roman mock data, $H_0 \approx 70.1$ km/s/Mpc.
  • Unconstrained Parameters: For current data (Pantheon+, DES), the interaction parameters ($\lambda, n_c$) and GB coupling remain largely unconstrained (flat posteriors), indicating the models are degenerate with $\Lambda$CDM at current precision.
  • Roman Data: The Roman mock data breaks this degeneracy. At high redshifts ($z > 1$), the models show a statistically significant departure from $\Lambda$CDM, with tight constraints on the interaction parameters.
• Statistical Preference:
  • Current Data: AIC shows a weak preference for the interacting models over $\Lambda$CDM, while BIC penalizes them due to extra parameters (no strong support).
  • Roman Data: Both AIC and BIC show strong statistical support for the interacting models over $\Lambda$CDM when Roman data is included, suggesting future surveys will favor dynamical dark energy.

5. Significance

• Viability of Scalar-Field DM: The study confirms that a coherent scalar field with a quadratic potential is a viable candidate for dark matter that can interact with dark energy without violating observational constraints.
• Resolution of Tensions: While the models do not fully resolve the $H_0$ tension (yielding values closer to Planck than SH0ES), they provide a physically consistent framework that avoids the theoretical pathologies often found in phenomenological fluid models.
• Future Observability: The paper highlights that the distinction between standard $\Lambda$CDM and these interacting EsGB models is subtle at low redshifts but becomes pronounced at high redshifts ($z > 1$). This underscores the critical importance of upcoming high-redshift surveys like the Roman Space Telescope in testing the nature of dark energy and modified gravity.
• Theoretical Consistency: The work demonstrates that EsGB gravity can be compatible with gravitational wave speed constraints while still offering rich cosmological dynamics through scalar-field interactions, bridging the gap between string-inspired gravity and observational cosmology.