On the Metric $f(R)$ gravity Viability in Accounting for the Binned Supernovae Data

This paper investigates two metric $f(R)$ gravity models using binned Type Ia Supernovae data to explain the effective running of the Hubble constant, demonstrating that while a general phenomenological fit is possible, it yields unphysical scalar field masses unless an additional condition is imposed to restore a consistent Cauchy problem and ensure a finite, positive scalar field mass.

The Gist

Imagine the universe as a giant, expanding balloon. For decades, scientists have been trying to figure out exactly how fast this balloon is inflating. They use "Standard Candles" (like Type Ia Supernovae, which are exploding stars that all shine with the same brightness) to measure distances. By looking at how far away these stars are and how fast they are moving away, they calculate the Hubble Constant ($H_0$), which is essentially the speed limit of the universe's expansion.

Here is the problem: When scientists look at very distant stars (high redshift), the calculated speed limit seems to be slower than when they look at nearby stars. It's as if the universe's expansion speed is "running" or changing over time, which contradicts our standard model of cosmology (called $\Lambda$CDM), where the expansion rate should be more consistent.

This paper is a detective story about two different theories trying to explain this "running" speed, using a modified version of Einstein's gravity called $f(R)$ gravity.

The Cast of Characters

1. The Standard Model ($\Lambda$CDM): The current "gold standard" theory. It assumes the universe is made of normal matter, dark matter, and dark energy (a mysterious force pushing the universe apart). It's like a well-oiled machine that works great, but it's starting to show cracks when looking at the new data.
2. The "Running" Hubble Constant: The observation that the expansion speed seems to change depending on how far back in time (or how far away) we look.
3. The Scalar Field: In these modified gravity theories, gravity isn't just a force; it's like a fluid or a field that can wiggle. This "wiggle" has a mass. If the mass is negative or infinite, the theory breaks (like a car engine exploding).

The Investigation: Two Theories Tested

The authors tested two different ways to fix the universe's expansion speed using $f(R)$ gravity.

Theory 1: The "Rigid Blueprint" (The Failed Attempt)

The first approach was like trying to fit a square peg into a round hole by forcing the peg to change shape.

• The Idea: They tried to build a mathematical "blueprint" (a function called $f(z)$) that directly described how the expansion speed changes. They made this blueprint flexible enough to match the data perfectly.
• The Result: It fit the data! The math matched the observations of the exploding stars.
• The Catch: When they checked the physics behind the blueprint, it fell apart. The "scalar field" (the wiggle in gravity) developed a negative mass (which causes instability, like a ball rolling up a hill forever) or an infinite mass (which breaks the laws of physics).
• The Analogy: Imagine you build a bridge that perfectly matches the shape of a river. It looks great on paper. But when you try to build it, you realize the steel beams are made of jelly. The bridge fits the river, but it collapses under its own weight. This theory was mathematically pretty but physically impossible.

Theory 2: The "Flexible Mechanic" (The Successful Attempt)

The second approach was smarter. Instead of forcing the blueprint to fit the data, they let the "engine" of the universe dictate the rules.

• The Idea: They added a specific rule (an "additional condition") to their equations. Think of this as giving the universe's engine a new governor that prevents it from revving too high or stalling. This rule allows the "scalar field" to have a healthy, positive mass.
• The Result:
  • Physical Viability: The scalar field is stable. It has a positive mass and doesn't explode. It's a "healthy" theory.
  • Data Fit: It still fits the supernova data very well. In fact, for one of the datasets (Pantheon), it fits better than the standard model. For the larger dataset (Master), it's just as good as the standard model.
• The Analogy: Instead of forcing a jelly bridge, they built a suspension bridge. It's flexible. It sways with the wind (the changing expansion rate) but stays standing because the cables (the new rule) are strong and properly tensioned.

The Big Discovery: Why the "Extra Rule" Matters

The most important part of this paper isn't just that they found a working model; it's that they explained why the extra rule was necessary.

In the first failed theory, the scientists accidentally tied the hands of the universe. They forced the "scalar field" to start with zero speed and stay at a specific value. This was like telling a runner, "You must start at the starting line, but you cannot move your legs until you reach the finish line." It's a contradiction that leads to a crash.

The second theory realized that to make the math work, the universe needs the freedom to start moving its "legs" (the scalar field derivative) immediately. The "extra condition" they added essentially says: "Okay, you can start moving right away, as long as you follow this specific path."

This condition isn't just a random trick; it's a safety mechanism that ensures the theory doesn't break the laws of physics while still explaining the weird data.

The Conclusion

The paper concludes that:

1. The Standard Model is struggling: The data suggests the universe's expansion rate is changing in a way the standard model can't easily explain.
2. Modified Gravity is a contender: A specific type of modified gravity ($f(R)$) can explain this data.
3. But you have to do it right: You can't just throw random math at the problem. You need to ensure the "engine" (the scalar field) is healthy. The authors found a specific way to set up the equations so the theory is both mathematically accurate (fits the stars) and physically safe (doesn't break reality).

In short: They found a new way to describe gravity that explains why the universe seems to be expanding at a "running" speed, but only if you add a specific safety valve to keep the theory from collapsing. It's a step toward understanding the dark energy that is pushing our universe apart.

Technical Summary

1. Problem Statement

Recent analyses of binned Type Ia Supernovae (SNe Ia) data (specifically the Pantheon and Master samples) have revealed a decreasing trend in the inferred Hubble constant ($H_0$) as redshift ($z$) increases. This "effective running" of $H_0$ challenges the standard $\Lambda$CDM cosmological model and is linked to the Hubble Tension (the discrepancy between early-universe CMB measurements and late-universe distance ladder measurements).

The authors investigate whether metric $f(R)$ gravity in the Jordan frame can explain this phenomenon. A critical issue in previous studies is that while some $f(R)$ models fit the data, they often suffer from physical instabilities, specifically:

1. Tachyonic instabilities: The scalar field mass squared ($m^2$) becomes negative.
2. Divergences: The scalar field mass diverges at low redshift ($z \to 0$).
3. Over-constrained Cauchy problems: The reconstruction algorithms often implicitly fix both the scalar field value and its time derivative at the present epoch, violating the degrees of freedom required for a second-order differential equation.

The paper aims to determine if a physically viable $f(R)$ model exists that fits the binned SNe data without these pathologies and to clarify the dynamical origin of the "additional conditions" often imposed on Friedmann equations in such analyses.

2. Methodology

Data Sets

The study utilizes two binned datasets to reconstruct the effective Hubble parameter $H(z)$:

• Pantheon Sample: 1048 SNe Ia, binned into 40 equipopulated redshift bins ($z \lesssim 1.2$).
• Master Sample: A homogeneous compilation of 3714 unique SNe Ia (merging Pantheon, Pantheon+, JLA, and DES), binned into 20 equipopulated redshift bins ($0.0012 \leq z \leq 2.3$).
• Statistical Approach: The analysis employs non-linear fitting for the Pantheon sample and Markov Chain Monte Carlo (MCMC) methods (using Cobaya) for the Master Sample. Crucially, the likelihood function is tested within each bin to avoid assuming a Gaussian distribution where the data suggests heavy-tailed distributions (e.g., Student-t).

Theoretical Framework

The authors analyze metric $f(R)$ gravity in the Jordan frame, which is dynamically equivalent to a scalar-tensor theory with a scalar field $\phi = f_R$ and a potential $U(\phi)$.

• Diagnostic Tool: The "effective running Hubble constant" is defined as $H(z) = H_0 \frac{H_{NEW}(z)}{H_{\Lambda CDM}(z)}$.
• Dark Energy Parametrization: The Chevallier-Polarski-Linder (CPL) parametrization is used for the dark energy equation of state: $\omega(z) = \omega_0 + \omega_a \frac{z}{1+z}$.

Two Models Investigated

1. Model A (Phenomenological $f(z)$): The gravitational Lagrangian is encoded via a redshift-dependent function $f(z)$ introduced directly into the Friedmann equation. The function is approximated by a second-order Taylor expansion: $f(z) = 1 + az + bz^2$. This model attempts to reconstruct the scalar potential $U(\phi)$ and mass $m^2$ a posteriori.
2. Model B (Dynamical Constraint): A reformulation where an additional condition is imposed on the modified Friedmann equation: $6H\dot{\phi} = U(\phi)$. This condition treats the potential as a dynamical quantity determined by the evolution of the scalar field, rather than fixing the initial derivative of the scalar field to zero.

3. Key Contributions

1. Identification of a "No-Go" Theorem for Model A: The authors demonstrate that the standard phenomenological reconstruction (Model A), which fixes $\phi(0)=1$ and implicitly forces $\phi'(0)=0$, inevitably leads to a divergent scalar mass at $z=0$ and tachyonic instabilities ($m^2 < 0$) at $z > 0$. This proves that simply fitting a function $f(z)$ to data without respecting the Cauchy problem constraints yields unphysical theories.
2. Physical Justification for the Additional Condition: The paper provides a rigorous dynamical derivation showing that the "additional condition" ($6H\dot{\phi} = U(\phi)$) used in previous literature is not merely an ad-hoc phenomenological fix. It is necessary to:
   • Restore a well-posed Cauchy problem by allowing $\phi'(0) \neq 0$.
   • Ensure the scalar field mass remains finite and positive.
   • Avoid the over-constraining of the scalar field dynamics.
3. Reconstruction of a Viable $f(R)$ Model: By adopting Model B, the authors successfully reconstruct an $f(R)$ theory that fits the observational data while maintaining physical stability. They provide an analytical low-redshift reconstruction of the $f(R)$ function in the form of a quadratic expansion: $f(R) \approx m^2 B_0 + B_1 R + B_2 \frac{R^2}{m^2}$.

4. Results

Model A (Phenomenological $f(z)$)

• Statistical Fit: Provides a good fit to the Pantheon and Master data (better than $\Lambda$CDM for Pantheon in terms of $\chi^2_{red}$).
• Physical Viability: Failed.
  • The scalar mass squared $u^2(z)$ diverges at $z=0$.
  • $u^2(z)$ becomes negative for $z > 0$, indicating a tachyonic instability.
  • This failure is intrinsic to the formulation and independent of the dataset used.

Model B (Dynamical Constraint)

• Statistical Fit:
  • Pantheon Sample: The model is strongly preferred over $\Lambda$CDM ($\Delta AIC \approx 7$, $\Delta BIC \approx 5.3$). It also compares favorably with the Power Law model.
  • Master Sample: The model is statistically comparable to $\Lambda$CDM. While the Bayesian Information Criterion (BIC) slightly favors $\Lambda$CDM due to the penalty of extra parameters, the model remains competitive and physically consistent.
• Physical Viability: Successful.
  • The scalar field mass remains positive and finite across the entire redshift range ($0 \leq z \leq 2.3$).
  • The model avoids tachyonic instabilities.
  • The reconstructed $f(R)$ function shows a mild deviation from General Relativity ($B_1 \approx 0.92$) and a positive quadratic term ($B_2 > 0$), ensuring stability.
  • The deviation from GR at $z=0$ is interpreted as an effective description valid on cosmological scales, consistent with the chameleon mechanism which suppresses deviations in high-density local environments (like the Solar System).

5. Significance

• Resolution of the Hubble Tension: The study demonstrates that metric $f(R)$ gravity can account for the observed "running" of the Hubble constant in binned SNe data without requiring exotic dark energy or violating physical laws.
• Theoretical Rigor: It clarifies a long-standing ambiguity in modified gravity literature regarding the "additional condition" on the Friedmann equation. The authors prove this condition is a dynamical necessity to ensure a consistent Cauchy problem and a stable scalar sector, rather than a phenomenological trick.
• Viability of Modified Gravity: The paper establishes that viable $f(R)$ models exist which are consistent with late-time cosmological data. It highlights that the failure of previous attempts was due to over-constraining the initial conditions of the scalar field, not an inherent flaw in $f(R)$ gravity itself.
• Future Directions: The work suggests that future analyses of cosmological tensions must prioritize the physical consistency of the scalar degree of freedom (mass positivity and finiteness) alongside statistical goodness-of-fit.

In conclusion, the paper successfully bridges the gap between phenomenological data fitting and theoretical consistency in modified gravity, offering a robust $f(R)$ framework that explains the effective running of the Hubble constant while adhering to the fundamental requirements of stability and well-posed dynamics.