⚛️ general relativity

Comparison of $f(R,T)$ Gravity with Multiple Datasets

This paper analyzes $f(R,T)=R+\lambda T^\epsilon$ gravity models using multiple cosmological datasets, finding that the best-fit parameter $\epsilon=0.010^{+0.013}_{-0.021}$ is consistent with the standard cosmological model where $\epsilon=0$ .

Original authors: Vincent R. Siggia, Eric D. Carlson

Published 2026-07-20

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Vincent R. Siggia, Eric D. Carlson

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, invisible ocean. For decades, scientists have been trying to figure out why this ocean is not just sitting still, but actually speeding up, expanding faster and faster every second. The most popular map we have for this journey is called the "standard model," which relies on a mysterious, invisible push known as "dark energy" and a heavy, invisible anchor called "dark matter." But this map has some cracks in it. The numbers we get from looking at the oldest light in the universe don't quite match the numbers we get from looking at exploding stars nearby. It's like two GPS devices giving you different routes to the same destination.

To fix this, some scientists have wondered if the rules of the road themselves—the laws of gravity—might need a little tweak. Instead of gravity being a simple, unchangeable law, maybe it changes slightly depending on how much "stuff" (matter and energy) is around it. This idea is called modified gravity. It's like asking if a car's engine behaves differently when it's carrying a heavy load versus when it's empty. If gravity works differently, it might explain why the universe is accelerating without needing to invent new, invisible substances. This paper dives into one specific version of this idea to see if it can fix the map better than the old one.

The authors of this paper, Vincent R. Siggia and Eric D. Carlson, decided to test a specific theory called f(R, T) gravity. Think of standard gravity as a recipe that uses only one ingredient: the curvature of space (which they call R). This new theory adds a second ingredient: the "trace of the stress-energy" (which they call T). In plain English, T is a measure of how much matter and energy are packed into a specific spot. The theory suggests that gravity isn't just about the shape of space, but also about how much "stuff" is in it, and that these two things mix together in a specific way controlled by a number called ϵ (epsilon).

In their previous work, the authors looked only at data from exploding stars (Type Ia supernovae) and found that almost any value for this mixing number ϵ seemed to work, even values that were very different from zero. It was like trying to bake a cake and finding that you could use almost any amount of sugar and it still tasted okay. However, they suspected that if they looked at the whole picture, the answer might be more specific.

In this new study, they didn't just look at exploding stars. They gathered a massive "tasting panel" of data from four different sources to see which version of the recipe fits best:

The Cosmic Microwave Background (CMB): The faint afterglow of the Big Bang, like the static on an old TV that tells us what the universe looked like as a baby.
Baryon Acoustic Oscillations (BAO): Fossilized sound waves from the early universe that left a specific pattern in how galaxies are spaced out, acting like a cosmic ruler.
Cosmic Chronometers: Pairs of aging galaxies that act like clocks, helping scientists measure how fast the universe is expanding at different times.
Type Ia Supernovae: Updated data from the "Pantheon+SH0ES" project, which are the exploding stars used as standard candles to measure distance.

When they crunched the numbers with all this new, high-quality data, the results were much stricter than before. The "free-for-all" of possible values for ϵ disappeared. The data strongly suggested that the mixing number ϵ must be very close to zero. Specifically, they found the best fit is ϵ = 0.010 with a margin of error of +0.013 and -0.021.

What does this mean? A value of ϵ = 0 is exactly the same as the standard model (the one with just dark energy and dark matter). The authors found that the new, fancy gravity theory doesn't actually need to be "fancy" at all to explain the data. The best fit is essentially the standard recipe. While the new theory could work with a tiny tweak (that small 0.010 number), the data doesn't force us to use it. The standard model, where gravity behaves exactly as Einstein originally described without this extra "stuff" ingredient, still fits the universe's behavior perfectly well.

So, the story ends with a bit of a twist: the scientists went looking for a new, more complex way to explain the universe's expansion, but the evidence they gathered suggests that the old, simple explanation is still the champion. The universe is accelerating, but it seems to be doing so without needing to rewrite the fundamental laws of gravity. The "tweak" they tested is likely just a tiny fluctuation in the data, and the most probable answer is that the standard model is doing just fine.

Technical Summary: Comparison of $f(R, T)$ Gravity with Multiple Datasets

Problem Statement
The standard cosmological model ( $\Lambda$ CDM), comprising a cosmological constant and cold dark matter, successfully explains the accelerated expansion of the Universe observed via Type Ia supernovae (SNe Ia) and Cosmic Microwave Background (CMB) measurements. However, persistent tensions remain, specifically the Hubble tension between CMB and SNe Ia distance measurements, and discrepancies between CMB observations and matter distribution inferred from Baryon Acoustic Oscillations (BAO). To address these issues, alternative theories of gravity have been proposed, including $f(R)$ gravity and its extension, $f(R, T)$ gravity, where the action depends on the Ricci scalar $R$ and the trace of the stress-energy tensor $T$ . Previous work by the authors, focusing exclusively on SNe Ia data, suggested that a wide range of parameters in the model $f(R, T) = R + \lambda T^\epsilon$ were consistent with observations. This paper investigates whether incorporating a broader suite of cosmological datasets constrains the parameter $\epsilon$ more tightly, potentially ruling out significant deviations from the standard $\Lambda$ CDM limit ( $\epsilon = 0$ ).

Methodology
The authors analyze the specific model $f(R, T) = R + \lambda T^\epsilon$ within a flat Friedmann-Lemaître-Robertson-Walker (FLRW) metric. Unlike their previous studies which ignored radiation to focus on late-time SNe Ia, this analysis incorporates radiation effects to accommodate CMB and early-universe data.

Theoretical Framework:
- The action is defined as $S = \int d^4x\sqrt{-g} [L_m + L_r - \frac{1}{2\kappa^2}(R + \lambda T^\epsilon)]$ , where $L_m$ and $L_r$ represent matter and radiation Lagrangians.
- The field equations are derived, noting that in $f(R, T)$ gravity, the stress-energy tensor is not conserved ( $\nabla_\mu T^{\mu\nu} \neq 0$ ). Consequently, an effective conserved current density $J'^\mu$ is defined.
- The authors derive the equations of motion for the scale factor $a(t)$ and the Hubble parameter $H(z)$ , expressing them in terms of dimensionless variables $\tau$ and $x$ to handle the non-linear dependence on $T$ .
- The model predicts a modified expansion history where the matter density term acquires a dependence on $\epsilon$ , and an effective dark energy component emerges from the $T^\epsilon$ term.
Datasets and Statistical Analysis:
- The study utilizes five free parameters: $\epsilon$ (the model parameter), $\xi$ (a dimensionless coupling related to $\lambda$ ), $\omega_b$ (baryon density), $H_0$ (Hubble constant), and $M$ (absolute magnitude of SNe Ia).
- The analysis minimizes a total $\chi^2$ $χ^{2}$ statistic combining four distinct datasets:
  - CMB: Distance priors ( $R$ , $\ell_A$ ) from Planck 2018 results.
  - BAO: Model-independent ratios $d_z(z)$ and $A(z)$ from the Dark Energy Spectroscopic Instrument (DESI DR2) and previous surveys (SDSS, 6dF, WiggleZ, etc.), including covariance matrices for correlated data points.
  - Cosmic Chronometers (CC): Direct measurements of $H(z)$ from passively evolving galaxies.
  - Type Ia Supernovae: The Pantheon+SH0ES dataset (1701 light curves of 1550 SNe), which includes correlations and the local Hubble constant measurement from SH0ES.
- The parameter $\epsilon$ is treated as the primary parameter of interest, while the others are marginalized over as nuisance parameters. The probability distribution $\rho(\epsilon)$ is calculated by integrating the likelihood function $\exp(-\chi^2/2)$ over the nuisance parameters, subject to theoretical upper bounds on $\xi$ derived from the model's asymptotic behavior.

Key Contributions

Incorporation of Radiation: The authors extend the $f(R, T)$ formalism to include radiation contributions, a necessary step for analyzing CMB and BAO data which probe epochs where radiation is non-negligible.
Unified Analysis: This work represents the first comprehensive constraint on the $f(R, T) = R + \lambda T^\epsilon$ model using a simultaneous fit to CMB, BAO, CC, and the updated Pantheon+SH0ES SNe Ia data.
Theoretical Constraints: The authors identify and enforce theoretical upper bounds on the parameter $\xi$ based on the requirement that the scale factor remains real and well-defined, preventing unphysical breakdowns of the theory in the future evolution of the Universe.

Results
The statistical analysis yields a highly constrained probability distribution for the parameter $\epsilon$ .

The best-fit value is found to be $\epsilon = 0.010^{+0.013}_{-0.021}$ at the 68% confidence level.
The probability distribution $\rho(\epsilon)$ shows a precipitous drop for $\epsilon > 0.027$ . This suppression is driven by the theoretical upper bound on $\xi$ ( $\xi_{lim}$ ); as $\epsilon$ increases, the best-fit $\xi_{min}$ approaches and eventually exceeds $\xi_{lim}$ , rendering those models physically invalid or highly disfavored.
The standard $\Lambda$ CDM model, corresponding to $\epsilon = 0$ , lies well within the tight error bars of the best-fit result.
While previous work suggested that negative values of $\epsilon$ (down to $-\infty$ ) were competitive with $\Lambda$ CDM, the inclusion of the expanded supernova data and DESI DR2 BAO data in this study significantly tightens the constraints, rendering the model indistinguishable from the standard model within current observational precision.

Significance and Claims
The paper concludes that current cosmological data do not provide a preference for the specific $f(R, T) = R + \lambda T^\epsilon$ modification of gravity over the standard $\Lambda$ CDM model. The authors state that the inclusion of the latest supernova data (Pantheon+SH0ES) and the DESI DR2 BAO measurements strongly constrains $\epsilon$ to be very small. Consequently, the results are "easily compatible with the $\Lambda$ CDM value of $\epsilon = 0$ ," suggesting that, at present, there is no observational evidence necessitating this specific modification to General Relativity. The study serves to narrow the viable parameter space for this class of modified gravity theories, reinforcing the robustness of the standard cosmological model against this particular alternative.

More like this