Testing Gauss-Bonnet Gravity with DESI BAO Data

Using Pantheon Plus supernovae, cosmic chronometer, and DESI BAO data, this study employs MCMC simulations to show that both power-law and exponential f(G) gravity models are statistically favored over the standard ΛCDM model, with the exponential variant uniquely predicting a future transition back to a decelerating phase.

The Gist

Imagine the universe as a giant, expanding balloon. For decades, scientists have been trying to figure out exactly how that balloon is inflating.

The standard story, known as the $\Lambda$CDM model, is like a recipe that says: "The balloon is expanding because of an invisible, mysterious force called 'Dark Energy' that pushes it outward." It works well, but it has some annoying glitches, like the fact that we have no idea what this "Dark Energy" actually is.

This paper proposes a different idea. Instead of adding a mysterious new ingredient (Dark Energy), what if we just tweak the recipe for gravity itself?

The Big Idea: Rewriting the Rules of Gravity

The authors are testing a theory called $f(G)$ Gravity.

• The Old Rule (General Relativity): Think of gravity as a rigid, unchangeable law of physics, like the rules of chess. It works perfectly for moving pieces on a board, but when you try to explain the whole universe, the rules seem to break down.
• The New Rule ($f(G)$ Gravity): This theory suggests that gravity isn't a rigid rulebook but more like a smart, shape-shifting material. It can stretch and change its behavior depending on how much "stuff" (matter) is around it.

The authors tested two specific "shapes" this gravity material could take:

1. The Power-Law Model: Imagine gravity acting like a rubber band that gets stronger or weaker in a predictable, mathematical way as the universe stretches.
2. The Exponential Model: Imagine gravity acting like a smart thermostat. It reacts very quickly to changes in the universe's size, adjusting its "push" in a complex, exponential curve.

The Experiment: Checking the Receipts

To see which theory is right, the authors didn't just sit in a lab; they went to the cosmic "supermarket" to check the receipts. They used three massive datasets to see which gravity recipe fits the data best:

1. Supernovae (Pantheon Plus): These are exploding stars that act like "standard candles." They tell us how far away things are, helping us measure the speed of the universe's expansion.
2. Cosmic Chronometers: These are like the universe's own clocks (old galaxies). By looking at how fast time is passing for them, we can measure the expansion rate directly.
3. DESI BAO Data: This is the newest and most detailed map of the universe's structure. It's like a giant 3D scan of where galaxies are sitting, showing the "fossilized" sound waves from the Big Bang.

The Results: The New Recipes Win

When the authors ran their numbers (using a super-computer method called MCMC, which is like running millions of simulations to find the perfect fit), they found something surprising:

The tweaked gravity recipes ($f(G)$) fit the data better than the standard "Dark Energy" recipe.

• The Verdict: Both new models were statistically "favored" over the standard model. It's as if you tried three different brands of coffee, and the two new, experimental blends tasted better to the judges than the famous brand everyone has been drinking for years.
• The Winner: The Power-Law model (the rubber band version) was the slight favorite, but the Exponential model (the thermostat version) was also a strong contender.

The Twist: A Surprising Future

Here is the most fascinating part. The standard model predicts the universe will keep speeding up forever, like a car on cruise control that never slows down.

However, the Exponential Model predicted something weird:

• It suggests that after speeding up for a while, the universe might slow down again in the distant future (around a time we call "redshift -0.1").
• The Analogy: Imagine a car that accelerates hard, but then the driver suddenly takes their foot off the gas and gently presses the brake. The universe might be heading for a "deceleration phase" instead of an eternal speed-up. This is a unique prediction that the standard model doesn't make.

Why Does This Matter?

This paper is a big deal because it suggests we might not need "Dark Energy" at all. Instead of inventing a mysterious invisible force to explain why the universe is expanding, we might just need to admit that our understanding of gravity is a little too simple.

In short: The universe isn't being pushed by a ghost (Dark Energy); it's just that gravity is a bit more flexible and creative than we thought. And if this new theory is right, the universe's party might eventually slow down and cool off, rather than speeding up forever.

Technical Summary

1. Problem Statement

The standard $\Lambda$CDM model, while successful in fitting many cosmological observations, faces significant theoretical and observational challenges, including the cosmological constant problem, the $H_0$ tension (discrepancy in Hubble constant measurements), and the $\sigma_8$ tension (discrepancy in matter clustering). These issues have motivated the exploration of Modified Gravity theories as alternatives to General Relativity (GR) to explain cosmic acceleration without invoking a fine-tuned cosmological constant.

Specifically, this paper focuses on $f(G)$ gravity, a modification where the action depends on an arbitrary function of the Gauss-Bonnet invariant ($G$), rather than just the Ricci scalar ($R$). While $f(G)$ theories offer potential to unify early inflation and late-time acceleration, their viability depends on identifying specific functional forms that are consistent with observational data. The authors aim to observationally constrain two specific $f(G)$ models using the latest high-precision cosmological data, particularly the newly released DESI (Dark Energy Spectroscopic Instrument) Baryon Acoustic Oscillation (BAO) measurements.

2. Methodology

Theoretical Framework

The authors work within a flat Friedmann-Robertson-Walker (FRW) universe. They modify the gravitational action by adding an $f(G)$ term:
$$S = \int d^4x \sqrt{-g} \left( \frac{R}{2\kappa^2} + f(G) + \mathcal{L}_m \right)$$
From this, they derive modified Friedmann equations where the Gauss-Bonnet fluid contributes an effective energy density ($\rho_G$) and pressure ($p_G$). They define a dimensionless Hubble parameter $E(z) = H(z)/H_0$ and derive normalized Friedmann equations for two specific models:

1. Model I (Power-law): $f(G) = \alpha G^\beta$
2. Model II (Exponential): $f(G) = \alpha G_0 \left( 1 - e^{-p G/G_0} \right)$

In both cases, the parameters are constrained such that the models recover the standard $\Lambda$CDM limit under specific conditions (e.g., $\beta \to 0$ or $p \to \infty$).

Datasets

The analysis utilizes three primary observational datasets:

1. Pantheon Plus (PP): 1,701 Type Ia Supernovae (SNe Ia) covering $0.001 \leq z \leq 2.3$.
2. Cosmic Chronometers (CC): 36 direct measurements of the Hubble parameter $H(z)$ derived from differential ages of passively evolving galaxies.
3. DESI BAO: The first-year data release from DESI, providing 12 measurements of comoving distances ($D_M, D_H, D_V$) across various tracers (BGS, LRGs, ELGs, Quasars, Ly$\alpha$ forest) in the redshift range $0.1 < z < 4.2$.

The authors analyze two dataset combinations:

• Dataset I: PP + CC
• Dataset II: PP + CC + DESI BAO

Statistical Analysis

• Parameter Estimation: The authors employ Markov Chain Monte Carlo (MCMC) simulations to constrain the free parameters ($\Omega_m, h, \beta$ for Model I; $\Omega_m, h, p$ for Model II) by numerically solving the modified Friedmann equations.
• Model Selection: To compare the performance of the $f(G)$ models against the standard $\Lambda$CDM model, they use information criteria:
  • Akaike Information Criterion (AICc): Corrected for small sample sizes.
  • Bayesian Information Criterion (BIC): Penalizes model complexity more heavily.
  • The model with the lowest $\Delta$AICc and $\Delta$BIC (relative to the reference $\Lambda$CDM) is considered statistically favored.

3. Key Contributions

1. First Application of DESI BAO to $f(G)$ Gravity: This is one of the first studies to utilize the high-precision DESI Year-1 BAO data to constrain Gauss-Bonnet modified gravity models, breaking parameter degeneracies that exist when using only SNe Ia and CC data.
2. Comparative Analysis of Functional Forms: The paper rigorously compares a power-law form against an exponential form within the $f(G)$ framework, determining which functional form is more consistent with current cosmological observations.
3. Prediction of Future Cosmic Evolution: The analysis reveals a unique dynamical behavior in the exponential model, predicting a future transition from acceleration back to deceleration, a feature absent in both the power-law model and $\Lambda$CDM.

4. Key Results

Parameter Constraints

Using the combined Dataset II (PP + CC + DESI BAO):

• $\Lambda$CDM: $\Omega_m = 0.280 \pm 0.011$, $h = 0.7139 \pm 0.0082$.
• Model I (Power-law): $\Omega_m = 0.230^{+0.027}_{-0.024}$, $h = 0.7115 \pm 0.0083$, $\beta = 0.192^{+0.068}_{-0.062}$.
• Model II (Exponential): $\Omega_m = 0.275 \pm 0.011$, $h = 0.7085 \pm 0.0084$, $p = 4.18^{+0.35}_{-0.40}$.

The $f(G)$ models generally yield slightly lower values for $\Omega_m$ and $h$ compared to $\Lambda$CDM, though within $1\sigma$ or $2\sigma$ confidence levels.

Statistical Preference

• Model I vs. $\Lambda$CDM: Model I is strongly favored.
  • $\Delta$AICc = $-12.29$
  • $\Delta$BIC = $-6.91$
• Model II vs. $\Lambda$CDM: Model II is also favored, but less so than Model I.
  • $\Delta$AICc = $-10.72$
  • $\Delta$BIC = $-5.27$
• Conclusion: Both $f(G)$ models provide a better statistical fit to the data than the standard $\Lambda$CDM model, with the Power-law (Model I) being the preferred candidate.

Cosmic Dynamics (Deceleration Parameter $q(z)$)

• Transition Redshift ($z_t$): The redshift at which the universe transitions from deceleration to acceleration:
  • Model I: $z_t \approx 0.793$
  • Model II: $z_t \approx 0.788$
  • $\Lambda$CDM: $z_t \approx 0.720$
  • Note: The $f(G)$ models predict an earlier transition than $\Lambda$CDM, consistent with some observational constraints.
• Present-day Deceleration ($q_0$):
  • Model I: $q_0 = -0.51$
  • Model II: $q_0 = -0.287$
  • $\Lambda$CDM: $q_0 = -0.577$
• Future Evolution: A distinctive result for Model II is the prediction of a future transition at $z \approx -0.1$, where the universe would return to a decelerating phase. In contrast, Model I and $\Lambda$CDM predict continued acceleration into the infinite future.

5. Significance and Conclusion

The paper demonstrates that $f(G)$ gravity is a viable alternative to $\Lambda$CDM for describing the universe's expansion history. The inclusion of DESI BAO data significantly tightens constraints on the model parameters.

• Statistical Robustness: The negative $\Delta$AICc and $\Delta$BIC values indicate that the modified gravity models are not only consistent with data but statistically preferred over the standard model, suggesting that the "dark energy" component might be an artifact of modified gravity rather than a new fluid.
• Model Selection: The Power-law model ($f(G) \propto G^\beta$) outperforms the Exponential model, suggesting that simple power-law deviations from the Gauss-Bonnet invariant are more compatible with current data.
• Theoretical Implications: The prediction of a future deceleration phase in the exponential model offers a testable hypothesis for future surveys. If future data confirms a return to deceleration, it would strongly support specific non-linear $f(G)$ forms over the standard $\Lambda$CDM paradigm.

The authors conclude that $f(G)$ gravity can successfully describe both the matter-dominated and dark-energy-dominated eras without invoking a cosmological constant, though further work is needed to test these models against perturbation data (large-scale structure growth) and to resolve current cosmological tensions ($H_0$ and $S_8$).