Constraints on Neutrino Physics from DESI DR2 BAO and DR1 Full Shape

The Dark Energy Spectroscopic Instrument (DESI) Collaboration has obtained robust measurements of baryon acoustic oscillations (BAO) in the redshift range, $0.1 < z < 4.2$, based on the Lyman-$\alpha$ forest and galaxies from Data Release 2 (DR2). We combine these measurements with external cosmic microwave background (CMB) data from Planck and ACT to place our tightest constraints yet on the sum of neutrino masses. Assuming the cosmological $\Lambda$CDM model and three degenerate neutrino states, we find $\sum m_\nu<0.0642$ eV (95%) with a marginalized error of $\sigma(\sum m_\nu)=0.020$ eV. We also constrain the effective number of neutrino species, finding N_\rm{eff} = 3.23^{+0.35}_{-0.34} (95%), in line with the Standard Model prediction. When accounting for neutrino oscillation constraints, we find a preference for the normal mass ordering and an upper limit on the lightest neutrino mass of $m_l < 0.023$ eV (95%). However, we determine using frequentist and Bayesian methods that our constraints are in tension with the lower limits derived from neutrino oscillations. Correcting for the physical boundary at zero mass, we report a 95% Feldman-Cousins upper limit of $\sum m_\nu<0.053$ eV, breaching the lower limit from neutrino oscillations. Considering a more general Bayesian analysis with an effective cosmological neutrino mass parameter, $\sum m_{\nu,\rm{eff}}$, that allows for negative energy densities and removes unsatisfactory prior weight effects, we derive constraints that are in $3\sigma$ tension with the same oscillation limit. In the absence of unknown systematics, this finding could be interpreted as a hint of new physics not necessarily related to neutrinos. The preference of DESI and CMB data for an evolving dark energy model offers one possible solution. In the $w_0w_a$CDM model, we find $\sum m_\nu<0.163$ eV (95%), relaxing the neutrino tension. [Abridged]