We extend the analysis presented in “Study of baryon acoustic oscillations with SDSS DR13 data and measurements of and ”  to include neutrino masses. The present analysis has three steps: 1) we calculate the distance of propagation , in units of , referred to the present time, of sound waves in the photon-electron-baryon plasma until decoupling by numerical integration of Equation (16) and Equation (17) of Ref.  ; 2) we fit the Friedmann equation of evolution of the universe to 21 independent Baryon Acoustic Oscillation (BAO) distance measurements listed in  used as uncalibrated standard rulers and obtain the length d of these rulers, in units of , referred to the present time; and 3) we set
to constrain the sum of neutrino masses . c is the speed of light, and is the present day Hubble expansion parameter.
2. Constraints on Neutrino Masses
The main body of this article assumes: 1) flat space, i.e. , and 2) constant dark energy density relative to the critical density, i.e. independent of the expansion parameter a. These constraints are in agreement with all observations to date   . Results without these constraints are presented in Appendix 1. Results with partial data sets are presented in Appendix 2.
To be specific we consider three active neutrino flavors with three eigenstates with nearly the same mass , so . This is a useful scenario to consider since our current limits on are much larger than the mass-squared-differences and obtained from neutrino oscillations  . These neutrinos become non-relativistic at a neutrino temperature or a photon temperature . The corresponding expansion parameter is .
The matter density relative to the present critical density is for . includes the density of Dirac or Majorana neutrinos that are non-relativistic today. Note that for Dirac neutrinos we are considering the scenario in which right-handed neutrinos and left-handed anti-neutrinos are sterile and never achieved thermal equilibrium. Our results can be amended for other specific scenarios. For we take the matter density to be . The radiation density is for , where for three flavors of Dirac (mostly) left-handed neutrinos and right-handed anti-neutrinos. We also take for three active flavors of Majorana left-handed and right-handed neutrinos. For , we take the radiation density to be . The present density of photons relative to the critical density is  .
The data used to obtain d are 18 independent BAO distance measurements with Sloan Digital Sky Survey (SDSS) data release DR13 galaxies in the redshift range to 0.7    summarized in Table 3 of  , two BAO distance measurements in the Lyman-alpha forest (Lyα) at (cross-correlation  ) and (auto-correlation  ) summarized in Section 6 of  , and the Cosmic Microwave Background (CMB) correlation angle   , used as an uncalibrated standard ruler. Note that the correlation angle is also determined by BAO. These 21 independent BAO measurements and full details of the fitting method are presented in  .
As a reference we take
(at 68% confidence) from “Planck TT + low P + lensing” data (that does not contain BAO information)  . is the present density of baryons relative to the critical density.
Due to correlations and non-linearities we obtain our final result (Equation (9) below) with a global fit. The following equations are included to illustrate the dependence of and d on the cosmological parameters h, and in limited ranges of interest. Integrating the comoving sound speed of the photon-baryon-electron plasma until with  we obtain
To obtain d we minimize the with 21 terms, corresponding to the 21 BAO observables, with respect to and d, and obtain and
with per degree of freedom 19.8/19, and correlation coefficient 0.989 between and d (this high correlation coefficient is due to the high precision of ). Setting we obtain
A more precise result is obtained with a global fit by minimizing the with 21 terms varying and directly. We obtain and
with , and correlation coefficient 0.924 between and . This is our main result. Equation (9) is obtained from BAO measurements alone, and is written in a way that can be combined with independent constraints on the cosmological parameters , and , such as measurements of the power spectrum of density fluctuations , the CMB, and direct measurements of the Hubble parameter.
Setting and we obtain the following upper bound on the mass of active neutrinos :
Appendix 1. Removing constraints
Freeing and keeping constant we obtain , , and
Fixing and letting we obtain , , and
Freeing and letting we obtain , , , and
Appendix 2. Removing data.
In this Appendix we apply the constraints and constant. Removing the measurement of we obtain and
Removing the measurement of and the two Lya measurements we obtain and
Keeping only the measurement of we need to fix in order to get zero degrees of freedom and have a unique solution. The best way to fix is with BAO measurements, and that is the purpose of the present study.
 Hoeneisen, B. (2017) Study of Baryon Acoustic Oscillations with SDSS DR13 Data and Measurements of and . International Journal of Astronomy and Astrophysics, 7, 11-27.
 Albareti, F.D., et al. (2016) The Thirteenth Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data from the SDSS-IV Survey Mapping Nearby Galaxies at Apache Point Observatory. SDSS Collaboration, arXiv:1608.02013.
 Dawson, K.S., Kneib, J.-P., Percival, W.J., et al. (2016) The SDSS-IV Extended Baryon Oscillation Spectroscopic Survey: Overview and Early Data. Astronomical Journal, 151, 44.