We investigate the impact of prior models on the upper bound of the sum of neutrino masses, $\sum m_{\nu }$. Using data from the large scale structure of galaxies, cosmic microwave background, type Ia supernovae, and big bang nucleosynthesis, we argue that cosmological neutrino mass and hierarchy determination should be pursued using exact models, since approximations might lead to incorrect and nonphysical bounds. We compare constraints from physically motivated neutrino mass models (i.e., ones respecting oscillation experiments) to those from models using standard cosmological approximations. The former give a consistent upper bound of $\sum m_{\nu }\lesssim 0.26\mathrm{eV}$ (95% CI) and yield the first approximation-independent upper bound for the lightest neutrino mass species, $m_0^{\nu }<0.086\mathrm{eV}$ (95% CI). By contrast, one of the approximations, which is inconsistent with the known lower bounds from oscillation experiments, yields an upper bound of $\sum m_{\nu }\lesssim 0.15\mathrm{eV}$ (95% CI); this differs substantially from the physically motivated upper bound.

• Received 8 November 2018
• Revised 19 April 2019

DOI:https://doi.org/10.1103/PhysRevLett.123.081301