Investigating Late-Time Dark Energy and Massive Neutrinos in Light of DESI Y1 BAO

Abstract: Baryonic Acoustic Oscillation (BAO) data from the Dark Energy Spectroscopic Instrument (DESI), in combination with Cosmic Microwave Background (CMB) data and Type Ia Supernovae (SN) luminosity distances, suggests a dynamical evolution of the dark energy equation of state with a phantom phase ($w < -1$) in the past when the so-called $w_0w_a$ parametrization $w(a) = w_0 + w_a(1-a)$ is assumed. In this work, we investigate more general dark energy models that also allow a phantom equation of state. We consider three cases: an equation of state with a transition feature, a model-agnostic equation of state with constant values in chosen redshift bins, and a k-essence model. Since the dark energy equation of state is correlated with neutrino masses, we reassess constraints on the neutrino mass sum focusing on the model-agnostic equation of state. We find that the combination of DESI BAO with Planck 2018 CMB data and SN data from Pantheon, Pantheon+, or Union3 is consistent with an oscillatory dark energy equation of state, while a monotonic behavior is preferred by the DESY5 SN data. Performing model comparison techniques, we find that the $w_0w_a$ parametrization remains the simplest dark energy model that can provide a better fit to DESI BAO, CMB, and all SN datasets than $\Lambda$CDM. Constraints on the neutrino mass sum assuming dynamical dark energy are relaxed compared to $\Lambda$CDM and we show that these constraints are tighter in the model-agnostic case relative to $w_0w_a$ model by $70\%-90\%$.