Interacting bosonic dark energy and fermionic dark matter in Einstein scalar Gauss-Bonnet gravity

This paper investigates a cosmological model where a Gauss-Bonnet-coupled scalar field (dark energy) interacts with fermionic dark matter, demonstrating that both exponential and power-law potentials yield expansion histories consistent with current observational data and the standard $\Lambda$CDM model while allowing for testable deviations in gravitational wave propagation.

The Gist

The Big Picture: A Cosmic Dance Between Ghosts

Imagine the universe is a giant, expanding stage. For a long time, scientists thought the show was being run by two invisible actors: Dark Matter (the heavy anchor holding galaxies together) and Dark Energy (the mysterious force pushing the universe apart).

For decades, we assumed these two actors never spoke to each other. They just did their own thing. But this paper asks a bold question: What if they are actually dancing together?

The authors propose a new story where Dark Matter and Dark Energy aren't just floating fluids; they are specific types of "fields" (like invisible waves) that are constantly exchanging energy. They are using a special rulebook called Einstein-Scalar-Gauss-Bonnet (EsGB) gravity to write this story.

The Special Rulebook: The "Gauss-Bonnet" Twist

To understand the rulebook, think of gravity not just as a simple pull, but as a fabric.

• Standard Gravity (Einstein): Imagine a trampoline. If you put a bowling ball on it, it curves. That's gravity.
• The New Rulebook (EsGB): This is like a trampoline made of a special, stretchy material that has a hidden "memory" or a "glitch" in its weave. This glitch is called the Gauss-Bonnet term. It's a mathematical feature that comes from advanced theories like String Theory (the idea that everything is made of tiny vibrating strings).

In this paper, the "Dark Energy" actor is a scalar field (a wave) that is tangled with this special "glitch" in gravity. This connection is crucial because it changes how gravity behaves, especially regarding Gravitational Waves (ripples in spacetime, like sound waves in a pond).

The Plot: Two Scenarios

The scientists tested two versions of this story to see which one fits the real universe best:

1. Scenario A (The "Speed Bump"): In this version, the interaction between the fields makes Gravitational Waves travel at a slightly different speed than light.
   • The Catch: We know from a famous event in 2017 (when two neutron stars crashed) that Gravitational Waves and light arrive at almost the exact same time. So, this version has to be very careful not to break that rule.
2. Scenario B (The "Perfect Sync"): In this version, the math is tweaked so that Gravitational Waves must travel at the exact speed of light. This is the "safe" version that definitely matches our current observations.

The Investigation: Watching the Movie Play Out

The authors didn't just write equations; they ran a simulation of the universe's entire history, from the Big Bang to today. They treated the universe like a video game where they could watch the "energy bars" of different components rise and fall.

They looked for a specific sequence of events:

1. The Radiation Era: The universe is hot and full of light particles.
2. The Matter Era: Gravity takes over, clumping matter together to form stars and galaxies.
3. The Dark Energy Era: The universe starts expanding faster and faster.

The Result: Both versions of their model successfully recreated this history. They showed that the universe could start with radiation, move to matter, and end up with Dark Energy dominating, just like our real universe does.

The Detective Work: Checking Against the Evidence

To see if their story is true, they compared their simulation against real-world data:

• The Cosmic Microwave Background (CMB): The "baby picture" of the universe.
• Supernovae: Exploding stars used as "standard candles" to measure distance.
• DESI & Roman Telescope Data: New, high-tech surveys mapping the positions of millions of galaxies.

They found that their "dancing" models fit the data almost as well as the standard model (called $\Lambda$CDM), which assumes Dark Matter and Dark Energy don't talk to each other. In fact, when they included future data from the Roman Space Telescope (a telescope that hasn't launched yet but is simulated in the paper), their model looked even better than the standard one.

Why This Matters

1. Solving the "Hubble Tension": There is a current crisis in physics where different methods of measuring the universe's expansion rate give different answers. This model suggests that because Dark Matter and Dark Energy are exchanging energy, it might slightly change the expansion rate, potentially helping to solve this puzzle.
2. No "Phantom" Energy: Some theories predict "phantom energy," which would cause the universe to rip itself apart in a "Big Rip." This model naturally avoids that disaster, keeping the universe stable.
3. Future Proofing: The model predicts subtle differences that future telescopes (like the Roman Space Telescope) will be able to spot. If we see those differences, we'll know that Dark Matter and Dark Energy are indeed dancing together.

The Bottom Line

This paper suggests that the invisible forces running our universe are more connected than we thought. Instead of two strangers ignoring each other, Dark Matter and Dark Energy might be partners in a cosmic dance, influenced by a deep, stringy structure of gravity. While we can't prove it yet, the math works, it fits the data we have, and it promises to be testable with the next generation of space telescopes.

Technical Summary

1. Problem Statement

The standard $\Lambda$CDM cosmological model, while successful, faces significant theoretical and observational challenges:

• Theoretical Issues: The cosmological constant ($\Lambda$) suffers from the fine-tuning problem and the vacuum energy discrepancy with quantum field theory predictions.
• Observational Tensions: There are persistent discrepancies, most notably the Hubble tension (disagreement between $H_0$ inferred from CMB vs. local distance ladders) and the $S_8$ tension (clustering of matter).
• Gravitational Wave Constraints: The joint detection of GW170817 (gravitational waves) and GRB 170817A (gamma rays) imposed a stringent constraint that the speed of gravitational waves ($c_T$) must equal the speed of light ($c$) to within $10^{-15}$. This rules out many modified gravity theories that predict $c_T \neq c$.
• Dark Sector Interaction: The nature of the interaction between Dark Matter (DM) and Dark Energy (DE) remains unknown. Phenomenological fluid models often introduce arbitrary coupling terms.

Objective: The authors propose a theoretically motivated framework where DE is a bosonic scalar field coupled to the Gauss-Bonnet (GB) term, and DM is a fermionic field. They investigate how a non-gravitational interaction between these fields affects cosmic evolution while satisfying the GW speed constraint.

2. Methodology

Theoretical Framework

The study is based on Einstein-Scalar-Gauss-Bonnet (EsGB) gravity, a low-energy limit of string/M-theory. The action includes:

• A scalar field $\phi$ (DE) non-minimally coupled to the GB invariant $\mathcal{G}$ via a function $f(\phi)$.
• A fermionic field $\psi$ (DM) described by a Dirac spinor.
• An interaction Lagrangian $L_{int} = -M(\phi)\bar{\psi}\psi$, representing a Yukawa-type coupling where the fermion mass depends on the scalar field. This provides a natural, particle-physics-based mechanism for energy exchange between DM and DE, avoiding ad-hoc source terms used in fluid models.

Dynamical System Analysis

The authors reformulate the field equations in a spatially flat FLRW background into an autonomous dynamical system.

• Variables: They introduce dimensionless variables ($x, y, z, u, v, \Omega_r, \Omega_b$) representing kinetic energy, potential energy, fermion density, and coupling terms.
• Scenarios: Two scenarios are analyzed regarding the GW speed constraint ($|c_T^2 - 1| < 5 \times 10^{-16}$):
  1. Model I: $c_T \neq 1$ (but constrained to be within observational bounds). Uses power-law potentials ($V \propto \phi^2$) and linear coupling ($M \propto \phi$).
  2. Model II: $c_T = 1$ exactly. Uses exponential potentials ($V \propto e^{-\gamma \phi}$) and exponential coupling ($M \propto e^{\mu \phi}$).
• Numerical Evolution: Instead of a full analytical fixed-point classification (which is complex due to high phase-space dimensionality), they numerically evolve the system from the deep radiation era ($N \approx -20$) to the future ($N > 0$) to check for stability and physical viability (bounded energy densities, correct sequence of epochs).

Observational Constraints

The model parameters are constrained using a comprehensive dataset via Markov Chain Monte Carlo (MCMC) methods:

• CMB: Planck 2018 distance priors.
• Hubble Parameter: Cosmic Chronometers (CC) and DESI BAO (DR2).
• Supernovae: Pantheon+ (PP), DESY5, and mock data from the upcoming Roman Space Telescope (simulated high-redshift SNe Ia).
• Statistical Tools: Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) are used to compare the models against standard $\Lambda$CDM.

3. Key Contributions

1. Particle-Based Interaction: The paper moves beyond fluid approximations by modeling DM as a fermionic field and DE as a scalar field, deriving the interaction naturally from the action via a field-dependent mass term. This resolves the arbitrariness often found in phenomenological coupling functions.
2. GW-Consistent EsGB Models: The authors explicitly construct models that satisfy the GW170817 constraint. They demonstrate that in Model II, the condition $c_T=1$ can be satisfied identically by specific choices of the coupling function, while Model I remains viable if the coupling is sufficiently suppressed at late times.
3. Robust Cosmological History: The dynamical analysis confirms that the models naturally reproduce the standard sequence of cosmic epochs (Radiation $\to$ Matter $\to$ Dark Energy) without fine-tuning initial conditions, provided the interaction strength is perturbative in the early universe.
4. Future Projections: The study incorporates mock Roman Space Telescope data, providing a forecast for how future high-precision observations could distinguish these interacting models from $\Lambda$CDM.

4. Results

• Dynamical Evolution: Both Model I and Model II successfully evolve from a radiation-dominated era to a matter-dominated era and finally to a late-time accelerated expansion phase. The effective equation of state ($w_{eff}$) transitions from $\approx 1/3$ to $\approx 0$ and finally to $w_{eff} < -1/3$, approaching $-1$ at late times.
• GW Speed Behavior:
  • In Model I, the GB coupling function dynamically evolves to very small values at low redshifts ($v \lesssim 10^{-23}$), effectively suppressing the GB term today and recovering General Relativity (GR) while allowing modifications at high redshifts.
  • In Model II, the constraint $c_T=1$ is built into the model structure via the specific form of the coupling function.
• Parameter Constraints:
  • The models yield a Hubble constant $H_0 \approx 68\text{--}70$ km/s/Mpc, showing a moderate alleviation of the Hubble tension compared to standard $\Lambda$CDM, particularly when Roman data is included.
  • The dark matter density parameter $\Omega_\psi$ is found to be consistent with current estimates ($\approx 0.27\text{--}0.29$).
• Statistical Comparison:
  • For current datasets (Pantheon+, DESY5), the interacting models show a weak preference over $\Lambda$CDM according to AIC, though BIC penalizes them slightly due to extra parameters.
  • Crucially, for the mock Roman dataset, both interacting models show strong statistical support (significant reduction in $\Delta$AIC and $\Delta$BIC) over $\Lambda$CDM, suggesting that future high-redshift supernova data could decisively favor interacting dark sector models.
• Stability: The models avoid the "phantom regime" ($w < -1$) and ghost instabilities, maintaining theoretical consistency throughout the evolution.

5. Significance

This work provides a theoretically robust alternative to the standard $\Lambda$CDM model that addresses the Hubble tension and the nature of dark sector interactions without violating gravitational wave constraints.

• Theoretical Viability: By grounding the interaction in particle physics (Yukawa coupling) and string theory (EsGB), it offers a more fundamental explanation for dark sector dynamics than phenomenological fluid models.
• Observational Relevance: The study highlights the critical role of future high-redshift observations (specifically from the Roman Space Telescope) in testing modified gravity and interacting dark energy scenarios. It suggests that the standard $\Lambda$CDM model may be an approximation of a more complex, interacting dark sector that will be revealed by next-generation data.
• Resolution of Tensions: The models offer a mechanism to slightly shift the inferred $H_0$ and matter density, potentially reconciling local and early-universe measurements without introducing pathological instabilities.

In conclusion, the paper demonstrates that Einstein-Scalar-Gauss-Bonnet gravity with interacting fermionic dark matter and bosonic dark energy is a viable, stable, and observationally competitive framework that warrants further investigation with upcoming cosmological surveys.