⚛️ high-energy theory

Revisiting the Hubble tension with an inverse power-law early dark energy

This paper proposes a new early dark energy model with an inverse power-law potential that alleviates the Hubble tension by yielding a best-fit $H_0$ value of approximately 70.56 km s $^{-1}$ Mpc $^{-1}$ when fitted to CMB, BAO, SN, and H0DN datasets, while also discussing its testability with future observations and associated fine-tuning challenges.

Original authors: Quan Zhou, Sibo Zheng

Published 2026-08-28

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Quan Zhou, Sibo Zheng

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

The universe is expanding, a fact established nearly a century ago. But when astronomers try to measure exactly how fast this expansion is happening today, they run into a stubborn disagreement. One method looks at the afterglow of the Big Bang, a faint light that has traveled across the cosmos for billions of years. When scientists analyze this ancient light, it suggests the universe is expanding at a certain, slower pace. Another method looks at the stars and galaxies right here in our cosmic neighborhood, measuring their distances and speeds directly. These local measurements consistently point to a faster expansion rate. The difference between these two numbers is small in absolute terms, but statistically significant enough to suggest that our current understanding of the universe's history might be missing a piece of the puzzle. This discrepancy is known as the Hubble tension, and it has become one of the most pressing mysteries in modern cosmology.

To solve this, researchers have proposed that something unusual happened in the very early universe, just a fraction of a second after the Big Bang. They call this hypothetical component "early dark energy." Unlike the dark energy we know today, which is pushing the universe apart at an accelerating rate, this early version would have been a temporary burst of energy that existed for a brief moment before fading away. If such a burst occurred, it would have altered the conditions of the early universe just enough to change how we interpret the ancient light, potentially reconciling the slow expansion rate seen in the past with the fast rate measured today. However, building a model for this energy that fits the data without creating new problems has proven difficult.

In a recent study, physicists Quan Zhou and Sibo Zheng from Chongqing University in China proposed a new way to describe this early dark energy. Instead of the complex shapes used in previous theories, they suggested a potential that behaves like a gentle slope with a long, tapering tail. Imagine a hill that is flat at the top but then slopes down gradually, never quite reaching the bottom but getting infinitely close. In their model, a field of energy sits on this slope. For a long time, the field remains stuck near the top, acting like a constant energy source. But as the universe expands and cools, the field eventually begins to roll down the slope. This rolling motion releases energy, creating the temporary burst needed to tweak the early universe's history. Once the field rolls far enough down, its energy becomes negligible, effectively disappearing from the cosmic equation and leaving the universe to evolve as we see it today.

The researchers tested this idea by feeding their model into powerful computer simulations that track the evolution of the universe from the Big Bang to the present day. They compared the results of their simulations against a vast collection of real-world data. This data included the detailed map of the cosmic microwave background, measurements of how galaxies are spaced out across the universe, observations of distant exploding stars, and the latest direct measurements of the expansion rate from a collaboration known as H0DN. When they ran the numbers, they found that their new model could indeed bring the two conflicting measurements of the expansion rate closer together. Specifically, for a specific version of their model where the slope of the energy hill has a particular steepness, the predicted expansion rate settled at a value of approximately 70.20 kilometers per second per megaparsec, with a best-fit value of 70.56. This result sits comfortably between the lower value predicted by the ancient light and the higher value measured locally, effectively reducing the tension between the two.

However, the solution is not without its own costs. While the model successfully eases the Hubble tension, it introduces a different kind of friction. The same adjustments that raise the expansion rate also tend to increase the clumpiness of matter in the universe. This creates a new conflict with observations regarding how much matter is clumped together, a problem known as the S8 tension. The researchers acknowledge that their model makes this existing problem slightly worse, even as it solves the expansion rate puzzle. Furthermore, the model relies on specific initial conditions that require a degree of fine-tuning, meaning the starting settings of the energy field must be set with high precision to work correctly. This suggests that while the mathematical shape of their potential is a promising candidate, the deeper theoretical reasons for why the universe would start in such a precise configuration remain unclear.

The study also highlights that this new model leaves a distinct fingerprint on the universe that future telescopes could detect. The researchers calculated that the presence of this early energy burst would cause tiny, specific deviations in the patterns of the cosmic microwave background and in the way matter is distributed across the cosmos. These deviations are currently too small to be seen with existing data, but they are large enough that upcoming missions, which will map the universe with much greater precision, should be able to spot them. If future observations confirm these specific patterns, it would provide strong evidence that this type of early dark energy really did exist. If they do not, the model will likely be ruled out.

Ultimately, this work represents a step forward in the ongoing effort to understand the universe's expansion. It offers a concrete, testable alternative to the standard model of cosmology, one that is mathematically distinct from previous attempts. While it does not provide a perfect, all-encompassing solution, it demonstrates that an energy field with a simple, inverse power-law shape can successfully bridge the gap between conflicting measurements of the Hubble constant. The path forward involves gathering more precise data to see if the universe's history truly bears the mark of this fleeting, early energy, or if the mystery of the Hubble tension requires an even more radical explanation.

Technical Summary: Revisiting the Hubble tension with an inverse power-law early dark energy

Problem Statement
The paper addresses the "Hubble tension," a significant discrepancy between the Hubble constant ( $H_0$ ) inferred from Cosmic Microwave Background (CMB) data under the standard $\Lambda$ CDM model ( $H_0 \approx 67.36$ km s $^{-1}$ Mpc $^{-1}$ ) and direct local measurements (e.g., SH0ES, $H_0 \approx 73.17$ km s $^{-1}$ Mpc $^{-1}$ ). The authors note that while Early Dark Energy (EDE) scenarios are among the few candidates capable of reducing this tension to below $3\sigma$ , existing models (such as axion-like, Rock 'n' Roll, and $\alpha$ -attractor EDE) suffer from theoretical fine-tuning issues. Specifically, axion-like models with index $n \ge 3$ face challenges in ultraviolet (UV) completion, while "new EDE" models lack symmetry protection against radiative corrections.

Methodology
The authors propose a new EDE toy model featuring a scalar field $\phi$ with a potential defined by an inverse power-law asymptotic tail:
$V(\phi) = V_0 \left[ 1 + \left(\frac{\phi}{f}\right)^2 \right]^{-n}$
where $V_0$ sets the energy scale, $f$ controls the field range, and $n$ determines the steepness.

The study employs the following methodological framework:

Numerical Implementation: The background and linear perturbation equations for the scalar field are implemented into the Boltzmann solver CLASS (via CLASS EDE) and interfaced with Cobaya for Markov Chain Monte Carlo (MCMC) analysis.
Parameter Mapping: A "shooting algorithm" is utilized to map the fundamental model parameters ( $V_0, f$ ) to phenomenological parameters ( $f_{\text{EDE}}$ , the maximum fraction of EDE energy density, and $z_c$ , the redshift at which this maximum occurs) required for the MCMC chains.
Datasets: The model is constrained using:
- CMB: Planck 2018 low- $\ell$ and high- $\ell$ (NPIPE CamSpec) temperature/polarization spectra and lensing.
- BAO: DESI DR2 isotropic and anisotropic distance constraints.
- SN: Pantheon+ Type Ia supernova data.
- Local $H_0$ : A Gaussian prior from the H0DN Collaboration ( $H_0 = 73.50 \pm 0.81$ km s $^{-1}$ Mpc $^{-1}$ ).
Statistical Metrics: Goodness-of-fit is evaluated using the Akaike Information Criterion ( $\Delta$ AIC) and the difference in maximum a posteriori estimators ( $Q_{\text{DMAP}}$ ) to quantify the tension relief and model preference.

Key Contributions and Results
The authors analyze the model for two representative indices, $n=8$ and $n=10$ , finding the following results:

Without Local $H_0$ Prior: When fitting only CMB+BAO+SN, the EDE model is disfavored compared to $\Lambda$ CDM (positive $\Delta$ AIC values of +6.2 and +5.1). The EDE fraction is constrained to be very small ( $f_{\text{EDE}} < 0.018$ ), and the tension remains unresolved.
With Local $H_0$ Prior: Including the H0DN prior significantly alters the constraints. The model prefers a non-negligible EDE fraction ( $f_{\text{EDE}} \approx 0.067$ $f_{EDE} \approx 0.067$ for $n=8$ $n = 8$ and $0.079$ for $n=10$ $n = 10$ ).
- Hubble Constant: For $n=10$ , the model yields a best-fit $H_0 \approx 70.56$ km s $^{-1}$ Mpc $^{-1}$ and a 68% confidence level value of $H_0 \approx 70.20^{+0.60}_{-0.47}$ km s $^{-1}$ Mpc $^{-1}$ .
- Tension Relief: The inclusion of the EDE model reduces the statistical conflict between the datasets and the $\Lambda$ CDM prediction. The $Q_{\text{DMAP}}$ drops from $6.25\sigma$ (in $\Lambda$ CDM) to $4.79\sigma$ for $n=10$ .
- Model Preference: The $\Delta$ AIC becomes negative ($-11.1$ for $n=10$ ), indicating a statistical preference for the EDE model over $\Lambda$ CDM when the local $H_0$ prior is included.
S8 Tension: Consistent with other EDE models, the authors find that the inferred value of $S_8$ increases (to $\approx 0.828$ for $n=10$ ), thereby worsening the $S_8$ tension.
Observable Signals: The paper predicts that the model leaves detectable imprints on CMB power spectra (specifically shifts in TT and EE at large $\ell$ ) and matter power spectra (deviations up to $\sim 10\%$ near $k \sim 1$ h/Mpc). These deviations are currently below cosmic variance limits but may be testable with future CMB and large-scale structure data.

Significance and Claims
The paper claims to offer a new theoretical avenue for addressing the Hubble tension through a scalar field potential with an inverse power-law tail. The authors highlight two main aspects of their work:

Phenomenological Success: The model successfully alleviates the Hubble tension to a level comparable to other successful EDE scenarios when local $H_0$ data is included, yielding $H_0$ values around $70.2 $–$ 70.6$ km s $^{-1}$ Mpc $^{-1}$ .
Theoretical Distinction: The proposed potential differs from existing axion-like or "new EDE" potentials, potentially offering a different approach to the fine-tuning problems, although the authors explicitly acknowledge that the model still suffers from fine-tuning issues regarding initial conditions.

The authors modestly conclude that while the model alleviates the Hubble tension, it does not resolve the $S_8$ tension. They suggest that future work could explore steeper inverse power-law potentials (e.g., $V(\phi) = V_0 \text{sech}^n(\phi/f)$ ) to further reduce tension levels and investigate the embedding of such potentials into UV-complete theories. They also note that the inclusion of galaxy clustering data could further impact the numerical results.

More like this