The cosmology of $f(R, L_m)$ gravity: constraining the background and perturbed dynamics

This paper investigates the cosmological viability of a specific $f(R, L_m)$ gravity model by constraining its parameters using MCMC simulations on combined observational Hubble, supernova, and large-scale structure datasets, ultimately comparing its ability to explain late-time accelerated expansion and cosmic structure growth against the standard $\Lambda$CDM model.

The Gist

Imagine the universe as a giant, expanding balloon. For decades, scientists have been trying to figure out why this balloon is not just expanding, but speeding up its expansion. The standard explanation (the "Old Guard" theory, called ΛCDM) says there's an invisible, mysterious force called Dark Energy pushing the balloon apart, like a hidden hand inflating it faster and faster.

But what if the "hand" isn't actually there? What if the rules of how the balloon stretches are just different than we thought?

This paper is about testing a new set of rules for the universe, called $f(R, L_m)$ gravity. Here is a simple breakdown of what the authors did, using everyday analogies.

1. The New Rulebook: $f(R, L_m)$ Gravity

In standard physics (Einstein's General Relativity), gravity is like a trampoline. If you put a heavy bowling ball (a star) on it, the fabric curves, and marbles (planets) roll toward it.

The authors propose a new rulebook where the trampoline fabric doesn't just react to the weight (matter), but also to how the fabric itself is stretching.

• $R$ represents the curvature of the fabric (geometry).
• $L_m$ represents the matter sitting on it.
• The new theory says these two things are "glued" together in a complex way. It's like saying the trampoline fabric has a mind of its own and changes its tension based on how many people are jumping on it.

The authors created a specific formula for this new gravity: $f(R, L_m) = \lambda R + \beta L_m^\alpha + \eta$.
Think of this formula as a recipe.

• $\lambda, \beta, \alpha, \eta$ are the ingredients (numbers) we need to find.
• If you pick the "standard" ingredients, the recipe tastes exactly like the Old Guard theory (ΛCDM).
• If you tweak the ingredients, you get a new flavor of gravity that might explain the universe's acceleration without needing mysterious Dark Energy.

2. The Taste Test: Checking the Recipe

You can't just guess the right ingredients; you have to taste the dish and see if it matches reality. The authors went to the cosmic kitchen and gathered the latest "taste tests" (data) from the universe:

• The Expansion History (OHD & SNIa): They looked at how fast the universe was expanding at different times in the past.

  • Analogy: Imagine looking at a time-lapse video of the balloon inflating. They checked if the new recipe predicts the balloon's size at every frame correctly.
  • Data used: 57 measurements of expansion speed and 1,048 measurements of exploding stars (Supernovae) which act as "standard candles" to measure distance.

• The Structure Growth (f & fσ8): They looked at how galaxies clump together to form clusters, like how dust motes clump in a sunbeam.

  • Analogy: If the universe is expanding too fast (due to Dark Energy), the "clumps" of galaxies should struggle to form. If the new gravity recipe is right, the clumps should form at a specific rate.
  • Data used: 14 and 30 measurements of how fast these galaxy clusters are growing.

3. The Results: Does the New Recipe Work?

The authors used a super-computer simulation (MCMC) to find the perfect mix of ingredients ($\alpha, \beta, \lambda$, etc.) that fits the data best.

The Good News:

• When they looked at the expansion speed (how fast the balloon is inflating), the new recipe worked surprisingly well! It fit the data almost as well as the standard theory.
• When they looked at galaxy growth (how the clumps form), the new recipe also showed "substantial support." It predicted the growth of the universe's structure very accurately.

The Bad News (The "But..."):

• The new recipe didn't work perfectly with every type of data.
• Specifically, when they used only the exploding star data (SNIa), the new recipe was rejected. It didn't taste right.
• When they combined the data, the new recipe was "okay" but not "amazing" according to some strict statistical judges (called BIC). It's like a dish that is delicious to some people but too spicy for others.

4. The Verdict

The paper concludes that this new $f(R, L_m)$ gravity is a viable contender.

• The Standard Theory (ΛCDM) is like a classic, reliable dish (e.g., a plain cheese pizza) that everyone knows and loves.
• The New Theory ($f(R, L_m)$) is like a gourmet fusion pizza.
  • It tastes great with certain toppings (expansion data and galaxy growth data).
  • It fails with others (some specific distance measurements).
  • It offers a fascinating alternative: maybe we don't need "Dark Energy" as a separate ingredient; maybe the "glue" between space and matter just works differently than we thought.

In Summary:
The universe is a complex puzzle. The authors tried a new piece (a new gravity formula) to see if it fits better than the old one. The piece fits well in some spots (explaining why the universe is speeding up and how galaxies form) but doesn't fit perfectly everywhere yet. It's a promising new direction that needs more testing with even more data before we can say, "Yes, this is the new rulebook for the universe."

Technical Summary

Here is a detailed technical summary of the paper "The cosmology of f(R,Lm) gravity: constraining the background and perturbed dynamics."

1. Problem Statement

The paper addresses the ongoing cosmological challenge of explaining the late-time accelerated expansion of the universe and the growth of cosmic structures. While the standard $\Lambda$CDM model (General Relativity with a Cosmological Constant) is successful, it relies on the enigmatic nature of Dark Energy ($\Lambda$) and Dark Matter.

The authors investigate an alternative Modified Theory of Gravity (MTG) known as $f(R, L_m)$ gravity, where the gravitational action depends on both the Ricci scalar ($R$) and the matter Lagrangian density ($L_m$). Specifically, they focus on a model that has not been comprehensively constrained by both background expansion data and large-scale structure (LSS) data simultaneously. The goal is to determine if this specific $f(R, L_m)$ model can viably replace or compete with $\Lambda$CDM when tested against modern observational datasets.

2. Methodology

Theoretical Framework

• Action: The study utilizes the generalized action $S = \int f(R, L_m)\sqrt{-g} d^4x$.
• Model Form: The authors propose a specific functional form:
  $$f(R, L_m) = \lambda R + \beta L_m^\alpha + \eta$$
  where $\lambda, \beta, \alpha,$ and $\eta$ are constants.
  • This model is designed to recover General Relativity (GR) when $\alpha=1, \beta=1, \lambda=1/2$.
  • It recovers the $\Lambda$CDM limit when $\eta = -\Lambda$.
• Field Equations: The modified field equations are derived by varying the action with respect to the metric. The authors assume a spatially flat Friedmann-Lemaître-Robertson-Walker (FLRW) metric and a perfect fluid matter source.
• Matter Lagrangian: Following standard conventions, the matter Lagrangian is chosen as $L_m = \rho$ (energy density), corresponding to dust-like matter.

Observational Data

The study employs a comprehensive set of datasets to constrain the model parameters:

1. Background Expansion Data:
   • OHD: 57 Observational Hubble Parameter data points (from cosmic chronometers and BAO).
   • SNIa: 1048 Type Ia Supernova distance moduli (Pantheon sample).
   • Combined: OHD + SNIa.
2. Large Scale Structure (LSS) Data (Perturbations):
   • $f$-data: 14 growth rate data points.
   • $f\sigma_8$-data: 30 redshift-space distortion data points.

Statistical Analysis & Tools

• Parameter Estimation: Markov Chain Monte Carlo (MCMC) simulations using the emcee (MCMC hammer) and GetDist packages were used to constrain the free parameters ($\Omega_m, H_0, \alpha, \beta, \lambda, \gamma, \sigma_8$).
• Model Comparison: The viability of the $f(R, L_m)$ model is compared against the standard $\Lambda$CDM model using:
  • Akaike Information Criterion (AIC): $\Delta AIC \leq 2$ (substantial support), $4 \leq \Delta AIC \leq 7$(less support),$\Delta AIC \geq 10$ (no support).
  • Bayesian Information Criterion (BIC): Similar thresholds applied to $\Delta BIC$.

3. Key Contributions

1. Comprehensive Constraint: Unlike previous works that focused primarily on background evolution, this paper provides a dual analysis of both the background expansion history and the linear perturbation dynamics (structure growth) for the specific $f(R, L_m) = \lambda R + \beta L_m^\alpha + \eta$ model.
2. Perturbation Formalism: The authors derive the linear cosmological perturbation equations within the $f(R, L_m)$ framework using the 1+3 covariant and gauge-invariant formalism. They derive the second-order evolution equation for the density contrast $\delta_m$ and the growth rate $f(z)$.
3. Growth Index ($\gamma$): The study calculates the growth index $\gamma$ for this modified gravity model, comparing it to the standard $\Lambda$CDM value ($\gamma \approx 6/11 \approx 0.545$).
4. Statistical Rigor: A detailed statistical evaluation using AIC and BIC criteria is performed across multiple distinct datasets to determine the relative "acceptance" or "rejection" of the model compared to $\Lambda$CDM.

4. Key Results

Background Cosmology (OHD, SNIa, OHD+SNIa)

• Parameter Constraints: Using the joint OHD+SNIa dataset, the best-fit parameters were found to be:
  • $\Omega_m \approx 0.287$
  • $H_0 \approx 71.72$ km/s/Mpc
  • $\alpha \approx 1.091$
  • $\beta \approx 1.237$
  • $\lambda \approx 0.630$
• Statistical Verdict:
  • OHD Data: The model shows substantial observational support ($\Delta AIC \leq 2$).
  • OHD+SNIa Data: The model shows less observational support ($4 \leq \Delta AIC \leq 7$).
  • SNIa Data Alone: The model is rejected ($\Delta AIC > 7$) compared to $\Lambda$CDM.
  • BIC Analysis: Based on $\Delta BIC$, the model generally lacks strong support for SNIa and combined datasets, suggesting the added complexity of the model is not fully justified by the background data alone.

Perturbed Dynamics (LSS Data: $f$ and $f\sigma_8$)

• Growth Index ($\gamma$):
  • For $\Lambda$CDM: $\gamma \approx 0.549 \pm 0.029$ (using $f$-data).
  • For $f(R, L_m)$: $\gamma \approx 0.555 \pm 0.014$ (using $f$-data).
  • The values are very close, indicating the model mimics standard growth behavior well.
• Parameter Constraints (using $f$-data):
  • $\Omega_m \approx 0.242$, $\alpha \approx 1.15$, $\beta \approx 1.12$, $\lambda \approx 0.72$.
• Parameter Constraints (using $f\sigma_8$-data):
  • $\Omega_m \approx 0.284$, $\sigma_8 \approx 0.799$, $\alpha \approx 0.766$, $\beta \approx 1.08$, $\lambda \approx 0.279$.
• Statistical Verdict (LSS):
  • $\Delta AIC$: The model shows substantial observational support for both $f$ and $f\sigma_8$ datasets ($\Delta AIC \approx 3.77$ and $2.40$ respectively, which falls in the "substantial" to "less support" range, but closer to substantial than the SNIa rejection).
  • $\Delta BIC$: The model shows less observational support ($\Delta BIC \approx 5.25$ and $5.40$), indicating that while the fit is good, the penalty for extra parameters is significant.

5. Significance and Conclusion

The paper concludes that the $f(R, L_m)$ gravity model is a viable alternative to $\Lambda$CDM, but its acceptance depends heavily on the dataset used:

• Strengths: The model is strongly supported by Hubble parameter (OHD) data and shows substantial support when analyzing structure growth ($f$ and $f\sigma_8$). It successfully reproduces the late-time acceleration and structure growth rates observed in the universe.
• Weaknesses: The model struggles when constrained solely by Supernova (SNIa) data, where it is statistically rejected in favor of $\Lambda$CDM. The Bayesian Information Criterion (BIC) generally penalizes the model more than the AIC, suggesting that the extra degrees of freedom ($\alpha, \beta, \lambda$) may not be strictly necessary given current data precision.

Final Outlook: The authors assert that while the model is promising, a definitive conclusion on its viability requires more extensive data, including larger sample sizes and future datasets (e.g., from Euclid or LSST), to distinguish between the subtle differences in the growth of structure predicted by $f(R, L_m)$ and standard $\Lambda$CDM. The study successfully bridges the gap between theoretical modification of gravity and observational constraints on both background and perturbed cosmology.