Extended dark energy analysis using DESI DR2 BAO measurements

This paper asks whether the late-time acceleration of the Universe is consistent with a cosmological constant (ΛCDM) or instead requires dynamical dark energy with a time-varying equation of state. The question matters because ΛCDM is extremely successful but conceptually incomplete: it treats dark energy as a constant vacuum term, while many theoretical ideas motivate evolving dark sectors. The authors focus on new, high-precision baryon acoustic oscillation (BAO) measurements from DESI Data Release 2 (DR2), and they test whether the inferred dark-energy behavior is robust to modeling choices—both parametric (fixed functional forms for the equation of state) and non-parametric (reconstructions with more flexibility).

The study is a combined-probe Bayesian inference analysis. The primary dataset is DESI DR2 BAO, spanning redshifts from approximately 0.1 to 4.2 and organized into multiple tracer samples: BGS (using isotropic DV/rd constraints at 0.1 < z < 0.4), LRG bins (0.4 < z < 0.6 and 0.6 < z < 0.8), LRG+ELG (0.8 < z < 1.1), ELG (1.1 < z < 1.6), and QSO (0.8 < z < 2.1), plus Lyα BAO (1.8 < z < 4.2). To anchor early-universe physics and break degeneracies, they include Planck CMB information using a full likelihood setup (high-1 TTTEEE plus low-1 TT and low-1 EE, plus lensing from Planck NPIPE PR4 and ACT DR6) implemented via Cobaya. For some non-parametric reconstructions that can struggle with negative dark-energy density, they also use a compressed CMB prior described by $({\theta }_{\ast },{\omega }_b,{\omega }_{bc}{)}_{\mathrm{CMB}}$, which captures key geometric and sound-horizon information while marginalizing over late-time effects.

For late-time distance calibration, the analysis combines DESI with Type Ia supernovae (SNe Ia) from three compilations: PantheonPlus (1550 spectroscopically confirmed SNe, 0.001 < z < 2.26), Union3 (2087 SNe, 0.01 < z < 2.26), and DESY5 (1635 photometric SNe, 0.1 < z < 1.13, plus 194 historical low-z SNe). The paper emphasizes Union3 in figures as a conservative choice due to larger uncertainties, but it reports how conclusions vary across the three SNe datasets.

Methodologically, the authors run Markov Chain Monte Carlo sampling (Metropolis-Hastings) in Cobaya, using CAMB with modifications to allow generalized dark-energy equations of state. They assume spatial flatness (motivated by prior DESI evidence that curvature is not significant). For perturbations and the possibility of crossing the phantom divide, they use the parametrized post-Friedmann (PPF) framework, which permits transitions across $w=-1$. They adopt priors summarized in their Table I, including a baseline $w_0$ and $w_a$ prior structure and additional priors for non-parametric expansions.

The key baseline parameterization is the Chevallier-Polarski-Linder (CPL) form, written as $w(z)=w_0+w_a\frac{z}{1+z}$ (equivalently $w(a)=w_0+w_a(1-a)$). In the DESI DR2 BAO + CMB + Union3 combination, the authors report a preference for a dynamical equation of state away from ΛCDM: the data favor the quadrant $w_0>-1$ and $w_a<0$, implying that $w(z)$ was phantom-like ($w(z)<-1$) in the past and evolves to $w(z)>-1$ today. They quantify the deviation from ΛCDM as an improvement in fit corresponding to roughly $2.8$ to $4.2\sigma$ depending on the dataset combination; they also state that DESI DR2 increases the statistical significance relative to DR1, with improvements in fit ranging from $-21.0\le \Delta {\chi }^2\le -10.7$ (where $\Delta {\chi }^2$ is defined relative to the ΛCDM best fit). They further note that DESI+CMB alone already suggests a deviation at about $\sim 3\sigma$, independent of any SNe compilation.

To test robustness, the paper explores multiple alternative two-parameter $w(z)$ forms (BA, EXP, LOG, and JBP) and finds that, except for JBP, these parameterizations yield similar low-redshift phantom-crossing behavior and comparable fit quality. They report explicit $\Delta {\chi }_{\mathrm{MAP}}^2$ values relative to ΛCDM for DESI+CMB+Union3: BA $\Delta {\chi }^2=-17.3$, EXP $=-17.5$, LOG $=-17.6$, JBP $=-13.6$, and CPL $=-17.4$. This indicates that the data do not strongly prefer one functional form over another; rather, they constrain a general trend.

A central theme is whether the apparent “phantom crossing” is genuine or an artifact of parameterization. The authors introduce crossing statistics using Chebyshev polynomial expansions of $w(z)$ and of the normalized dark-energy density $f_{\mathrm{DE}}(z)={\rho }_{\mathrm{DE}}(z)/{\rho }_{\mathrm{DE},0}$. They emphasize that expanding $f_{\mathrm{DE}}(z)$ allows the effective energy density to change sign, broadening the model space (and capturing behaviors that $w(z)$-only expansions may miss). Using a low-order Chebyshev expansion (they show results for $N=3$), they find reconstructed behaviors that agree well with the CPL $w_0{w}_a$ trends, including a smooth evolution that is consistent with a turnover and a crossing near $z\sim 0.5$ (with the exact crossing redshift depending on the chosen parameterization).

They also perform non-parametric reconstructions via two techniques. First, they use binning with smooth transitions controlled by a hyperbolic tangent smoothing scale (they set $s=0.02$ for bin edges). They test multiple binning schemes and find the tightest constraints at low redshift, where the data prefer $w(z)$ values more than $3\sigma$ away from $-1$ in the lowest redshift bin. At higher redshift bins, deviations are within $2\sigma$ of ΛCDM. For $f_{\mathrm{DE}}(z)$ binning (using compressed CMB priors to avoid computational issues with $f_{\mathrm{DE}}<0$), they observe a turnover in the range $0.5<z<1.0$ with $f_{\mathrm{DE}}(z)>0$ at around $2\sigma$ for most bins. They further apply PCA to decorrelate binned parameters, finding that the most informative principal components are localized at low redshift (the first component is localized in $0.1<z\le 0.3$ and has an uncertainty at least $20\times$ smaller than later components).

Second, they use Gaussian process (GP) regression to reconstruct $w(z)$ and $f_{\mathrm{DE}}(z)$ as smooth functions with minimal assumptions. They impose $w(z\ge z_{\max })=-1$ with $z_{\max }=10$ and use a squared-exponential kernel centered on $w=-1$ with hyperparameters controlling smoothness. The GP reconstructions align closely with the CPL $w_0{w}_a$ posterior predictions and indicate a phantom-like deviation at low redshift with hints of crossing near $z\approx 0.5$. They also show that including CMB tightens constraints (especially on ${\Omega }_m$), while DESI+SNe alone allows a broader range of $w(z)$ shapes.

To interpret the deviations physically, the authors examine three dark-energy dynamical classes: (i) thawing quintessence (minimally coupled scalar fields with $w\ge -1$), (ii) emergent dark energy (dark energy negligible for most of cosmic history and emerging recently), and (iii) mirage dark energy (a phenomenological class aligned with a specific degeneracy direction in the $w_0{w}_a$ plane). They report that thawing and emergent classes are not strongly favored, while mirage performs remarkably well, capturing the data with one additional degree of freedom.

They quantify model preference using both $\Delta {\chi }_{\mathrm{MAP}}^2$ and deviance information criterion (DIC). For DESI+CMB+Union3, the w0wa model yields $\Delta {\chi }_{\mathrm{MAP}}^2=-17.4$ and $\Delta \mathrm{DIC}=-17.2$ relative to ΛCDM. The calibrated thawing class gives only mild improvements (e.g., $\Delta {\chi }^2$ around $-2.5$ for PantheonPlus and $-7.1$ for DESY5, with DIC values less favorable than w0wa), while algebraic thawing improves more ($\Delta {\chi }^2=-2.9$ for Union3 and $-13.2$ for DESY5). The emergent class shows essentially no improvement for Union3 ($\Delta {\chi }^2\approx -0.1$) and small negative DIC changes. In contrast, the mirage class achieves strong improvements comparable to w0wa: for Union3, $\Delta {\chi }^2=-16.2$ and $\Delta \mathrm{DIC}=-18.7$.

A key limitation acknowledged by the authors is that “phantom crossing” inferred from $w_0{w}_a$ or other flexible reconstructions may be spurious: the parameterization may mimic observables without reflecting a true underlying $w(z)$ crossing. They address this by comparing to a phantom-restricted thawing model (algebraic thawing) that enforces $w(z)\ge -1$. Although the thawing model can fit somewhat better than ΛCDM, it is substantially less favored than w0wa, suggesting that the data prefer a sharp evolution feature (a rapid increase and then decrease in dark-energy density) that non-crossing models struggle to reproduce without fine-tuning.

The paper also includes validation on mock datasets. In two mocks (one generated from ΛCDM and one from a w0waCDM model), both binning and GP reconstructions recover the true $w(z)$ within $1\sigma$ in most cases. They also test an extreme thawing-like mock and show that their non-parametric implementation may fail to recover some extreme behaviors, indicating that priors and reconstruction choices can limit performance in unusual regions of model space.

Practically, the results imply that if the observed deviations are not due to unmodeled systematics, then dark energy likely evolves at low redshift ($z\lesssim 0.3$) and may exhibit an effective phantom-like behavior around $z\sim 0.5$. This matters for theorists building dark-sector models and for survey analysts planning next-generation cross-checks. The authors emphasize that decisive tests will require complementary probes beyond BAO and background distances, including redshift-space distortions and peculiar velocities (growth information), improved low-redshift supernova measurements, and future CMB experiments to tighten early-universe constraints and break degeneracies.

Overall, the paper’s core contribution is a robustness study: across multiple parametric forms and two non-parametric reconstruction methods, the inferred dark-energy evolution is stable and consistent with a dynamical deviation from ΛCDM, with the strongest and most localized evidence at low redshift and an apparent crossing/turnover feature near $z\approx 0.5$. While the phantom-crossing interpretation remains theoretically challenging and potentially model-dependent, the authors conclude that the evidence for dynamical dark energy is robust under modeling choices and that ΛCDM is disfavored at the several-sigma level with DESI DR2 BAO combined with CMB and SNe.