Precision Constraints on New Dark Energy Parametrization from DESI BAO DR2

1. Problem Statement

The standard cosmological model, $\Lambda$CDM (Lambda Cold Dark Matter), successfully explains a wide range of cosmological observations but faces significant theoretical and observational challenges:

• Theoretical Issues: It suffers from the cosmological constant problem (fine-tuning) and the cosmic coincidence problem.
• Observational Tensions: There is a persistent discrepancy known as the Hubble Tension, where early-universe measurements (e.g., Planck CMB) yield $H_0 \approx 67.4$ km s$^{-1}$ Mpc$^{-1}$, while late-universe distance-ladder measurements (e.g., SH0ES) yield $H_0 \approx 73.0$ km s$^{-1}$ Mpc$^{-1}$. This tension exceeds $5\sigma$ to $7\sigma$ depending on the datasets used.
• Static Dark Energy: The standard model assumes Dark Energy (DE) is a cosmological constant ($\Lambda$) with a fixed equation of state (EoS) $w = -1$. However, dynamical DE models are proposed to resolve these tensions and provide a physical basis for acceleration.

The authors aim to investigate a new single-parameter parametrization of the DE equation of state to test if a dynamical DE model can better fit recent high-precision data, specifically the DESI BAO Data Release 2 (DR2), and potentially alleviate the Hubble tension.

2. Methodology

2.1. Theoretical Framework

The authors assume a spatially flat, homogeneous, and isotropic FLRW universe. They introduce a new parametrization for the DE equation of state, $w_{de}(z)$, proposed in a previous study [42], defined as:
$$w_{de}(z) = -1 + \frac{(1+z)^\alpha \exp\left(-\frac{2\alpha}{1+z}\right) \alpha [\arctan(1+z)]^\alpha}{3(1 + (1+z)^{2\alpha})}$$

• Parameter $\alpha$: A free parameter. If $\alpha = 0$, the model reduces to $\Lambda$CDM ($w=-1$).
• Physical Interpretation:
  • $\alpha > 0$: Quintessence-like behavior ($w > -1$).
  • $\alpha < 0$: Phantom-like behavior ($w < -1$).
  • The model allows for a dynamical evolution of DE density $\rho_{de}$ dependent on redshift.

2.2. Datasets

The study utilizes three primary observational datasets to constrain the model parameters ($\alpha$, $H_0$, $\Omega_{m0}$):

1. Observational Hubble Data (OHD): 66 data points derived from the Cosmic Chronometer (CC) method, covering redshifts $0.07 \leq z \leq 1.965$.
2. Pantheon+SH0ES (PPS): A combined dataset of 1701 Type Ia Supernovae (SNe Ia) light curves from the Pantheon+ sample, anchored by SH0ES Cepheid host distances, covering $z \in [0.01, 2.26]$.
3. DESI BAO DR2: The second data release from the Dark Energy Spectroscopic Instrument, providing 9 Baryon Acoustic Oscillation (BAO) measurements (including BGS, LRG, ELG, QSO, and Ly-$\alpha$ forest) across a wide redshift range.

2.3. Statistical Analysis

• Likelihood Analysis: The authors employ a Markov Chain Monte Carlo (MCMC) algorithm using the emcee Python package to explore the parameter space.
• Priors: Broad flat priors are used ($H_0 \in [40, 100]$, $\Omega_{m0} \in [0.01, 0.6]$, $\alpha \in [-1, 2]$) to avoid bias and allow the data to drive the constraints.
• Model Selection: The goodness of fit is evaluated using $\chi^2_{min}$, reduced $\chi^2$, Akaike Information Criterion (AIC), and Bayesian Information Criterion (BIC) to compare the new model against the standard $\Lambda$CDM.

3. Key Contributions

1. Application of DESI BAO DR2: This is one of the first studies to apply the latest DESI BAO DR2 data to constrain this specific new DE parametrization.
2. Broad Prior Exploration: Unlike previous studies on this parametrization, this work uses significantly broader prior ranges, allowing for a more robust exploration of degeneracies and convergence to high-likelihood regions.
3. Hubble Tension Analysis: The study explicitly investigates whether this dynamical DE model can shift the inferred $H_0$ value to an intermediate range between Planck and SH0ES, potentially offering a partial resolution to the tension.
4. Physical Diagnostics: The paper provides a comprehensive physical analysis of the model's evolution, including the deceleration parameter, universe age, and the $Om(z)$ diagnostic test.

4. Key Results

4.1. Parameter Constraints

• Equation of State Parameter ($\alpha$):
  • Joint analysis (OHD + PPS + DESI BAO DR2) yields $\alpha = 0.239 \pm 0.07$ (68% CL).
  • This value is non-null, indicating a statistically significant deviation from the standard $\Lambda$CDM model ($\alpha=0$).
  • The positive value suggests a quintessence-like behavior ($w > -1$).
• Hubble Constant ($H_0$):
  • The model accommodates varying $H_0$ values depending on the dataset.
  • Joint Analysis: $H_0 = 68.40 \pm 0.23$ km s$^{-1}$ Mpc$^{-1}$.
  • This value lies between the Planck CMB estimate ($\sim 67.4$) and the SH0ES estimate ($\sim 73.0$), suggesting the model shifts the expansion history to an intermediate value.
  • In contrast, the $\Lambda$CDM model constrained by the same joint data yields $H_0 = 73.38 \pm 0.15$ km s$^{-1}$ Mpc$^{-1}$, aligning closely with SH0ES but failing to bridge the gap with early-universe data.
• Matter Density ($\Omega_{m0}$): Consistent with Planck CMB measurements, ranging from $0.28$ to $0.36$ depending on the dataset combination.

4.2. Model Comparison (AIC/BIC)

• OHD Only: $\Lambda$CDM is mildly favored.
• PPS Only: The new model is preferred by AIC ($\Delta AIC \approx -2.34$), though BIC favors $\Lambda$CDM.
• OHD + DESI BAO DR2: Strongly favors $\Lambda$CDM.
• Full Joint Dataset (OHD + PPS + DESI BAO DR2): The new model is strongly preferred with $\Delta AIC \approx -16.95$ and $\Delta BIC \approx -11.47$. This indicates that when all datasets are combined, the dynamical DE model provides a significantly better fit than $\Lambda$CDM.

4.3. Physical Behavior

• Deceleration Parameter ($q$): The model confirms a smooth transition from a past deceleration phase to the current acceleration.
  • Future Evolution: A unique finding is a second transition from acceleration back to deceleration in the future ($z < 0$) when PPS data is included. This feature is absent in datasets without SNe Ia, suggesting it is driven by late-time luminosity distance constraints.
  • Current value: $q_0 \approx -0.42$ (Joint).
• Age of the Universe: The estimated age is $t_0 \approx 12.80 - 13.67$ Gyr, consistent with Planck 2018 predictions ($13.78 \pm 0.20$ Gyr).
• $Om(z)$ Diagnostic: The $Om(z)$ parameter is not constant (unlike in $\Lambda$CDM) and shows a positive, increasing trend with decreasing redshift. This confirms the presence of a dynamical DE component and a quintessence-like nature.

5. Significance and Conclusion

The paper demonstrates that a new single-parameter DE parametrization is statistically favored over the standard $\Lambda$CDM model when constrained by the latest OHD, Pantheon+, and DESI BAO DR2 datasets.

• Resolution of Tensions: While the model does not completely resolve the Hubble tension, it successfully shifts the inferred $H_0$ to an intermediate value ($68.40$) that is compatible with both early and late-universe observations, unlike the standard model which forces a choice between them.
• Dynamical Nature of DE: The non-zero constraint on $\alpha$ provides observational evidence for a dynamical dark energy sector rather than a static cosmological constant.
• Future Outlook: The authors conclude that while the model shows minimal departure from $\Lambda$CDM in some aspects, it offers more robust estimates for cosmological parameters. Future work should incorporate CMB, SPT, and ACT data to further constrain the model and verify the predicted future transition to deceleration.

In summary, this work provides a rigorous, data-driven validation of a dynamical dark energy model, highlighting its potential to address current cosmological anomalies through the integration of the latest DESI BAO data.