Tanisha Jhaveri, Tanvi Karwal, Thomas Crawford, Wayne Hu, Ali Rida Khalife, Lennart Balkenhol, Fei Ge

Recent cosmological data reveal tension between parameters inferred from measurements of the cosmic microwave background (CMB), baryon acoustic oscillations (BAO), and supernovae (SN) under $Λ$CDM. Typical dynamical dark energy parameterizations (such as $w_0w_a$) that seek to jointly resolve these tensions have an equation of state parameter that crosses into the phantom regime, leading to potential instabilities for physical models. We show that the BAO (early-time) and SN (late-time) sides of the tension can instead be treated independently. Early dark energy (EDE) can reduce the tension between CMB-BAO data by changing the calibration of the sound horizon at the drag epoch $r_d$, with a $Δχ^2 = -{9.4}$ relative to $Λ$CDM, raising $H_0$ to 70.87 $\rm km s^{-1}Mpc^{-1}$. EDE alone cannot bring consistency between CMB, BAO, and SN data, but combining with a thawing-quintessence component of dark energy reduces tensions between the three datasets, with $Δχ^2=-12.6$ relative to $Λ$CDM without a phantom component, vs. $Δχ^2=-15.8$ for $w_0 w_a$ with one. We consider different SN datasets, using the most recent DES Dovekie catalog as our default while assessing differences with the original DESY5 and Pantheon+ catalogs. While the significance of adding thawing quintessence changes, the EDE solution to the CMB-BAO tension remains nearly unaffected. Moreover, though we do not include direct Hubble constant measurements in these $Δχ^2$ values, the EDE solution reduces the Hubble tension with the Local Distance Network value from $7σ$ in $Λ$CDM to $2-3σ$ depending on the SN dataset, nominally the equivalent of an extra $Δχ^2 \sim -40$ or more.