Study of Galaxy Distributions with SDSS DR14 Data and Measurement of Neutrino Masses

Show more

1. Introduction

We measure neutrino masses by comparing the predictions of the ΛCDM model with measurements of the power spectrum of linear density perturbations $P\left(k\right)$ . We consider three measurements of $P\left(k\right)$ : 1) the Sachs-Wolfe effect of fluctuations of the Cosmic Microwave Background (CMB) which is a direct measurement of density fluctuations [1] [2] ; 2) the relative mass fluctuations ${\sigma}_{8}$ in randomly placed spheres of radius ${r}_{s}=8/h$ Mpc with gravitational lensing and studies of rich galaxy clusters [2] [3] ; and 3) measurements of $P\left(k\right)$ inferred from galaxy clustering with the Sloan Digital Sky Survey [4] [5] [6] . Baryon Acoustic Oscillations (BAO) were considered separately [7] [8] and are not included in the present study, except for the final combinations.

To be specific, we consider three active neutrino eigenstates with nearly the same mass, so $\sum}\text{\hspace{0.05em}}{m}_{\nu}=3{m}_{\nu$ . This is a useful scenario to consider because the current limits on ${m}_{\nu}^{2}$ are much larger than the mass-squared-differences $\Delta {m}^{2}$ and $\Delta {m}_{21}^{2}$ obtained from neutrino oscillations [3] .

Figures 1-4 illustrate the problem at hand. Figures 1-3 present measurements of the “reconstructed” galaxy power spectrum ${P}_{\text{gal}}\left(k\right)$ obtained from the SDSS-III BOSS survey [4] , while Figure 4 presents the corresponding “standard” ${P}_{\text{gal}}\left(k\right)$ . The “reconstructed” ${P}_{\text{gal}}\left(k\right)$ is obtained from the directly measured “standard” ${P}_{\text{gal}}\left(k\right)$ by subtracting peculiar motions to obtain the power spectrum prior to non-linear clustering. Also shown are various fits to this data (with floating normalization), and to measurements of the Sachs-Wolfe effect, and ${\sigma}_{8}$ . The Sachs-Wolfe effect normalizes $P\left(k\right)$ , within its uncertainty, in the approximate range of $lo{g}_{10}\left(k/\left(h\text{\hspace{0.17em}}{\text{Mpc}}^{-1}\right)\right)$ from −3.1 to −2.7, while ${\sigma}_{8}$ is most sensitive to the range −1.3 to −0.6. Full details will be given in the main body of this article.

The fit in Figure 1 corresponds to the function

Figure 1. Comparison of ${P}_{\text{gal}}\left(k\right)$ obtained from the SDSS-III BOSS survey [4] (“reconstructed”) with ${b}^{2}P\left(k\right)$ obtained from a fit of Equation (5) with $\sum}\text{\hspace{0.05em}}{m}_{\nu}=0$ eV to the Sachs-Wolfe effect, ${\sigma}_{8}$ , and ${P}_{\text{gal}}\left(k\right)$ . The fit obtains $A=8738\text{\hspace{0.17em}}{\text{Mpc}}^{3}$ , ${k}_{\text{eq}}=0.068h\text{\hspace{0.17em}}{\text{Mpc}}^{-1}$ , $\eta =4.5$ , and ${b}^{2}=1.8$ , with ${\chi}^{2}=24.7$ for 19 degrees of freedom. Also shown for comparison is the curve with the same parameters, except $\sum}\text{\hspace{0.05em}}{m}_{\nu}=0.6$ eV.

Figure 2. Same as Figure 1, except that the curve “ $\sum}\text{\hspace{0.05em}}{m}_{\nu}=0.6$ eV” is fit. We obtain $A=9312\text{\hspace{0.17em}}{\text{Mpc}}^{3}$ , ${k}_{\text{eq}}=0.080h\text{\hspace{0.17em}}{\text{Mpc}}^{-1}$ , $\eta =4.2$ , and $\kappa =1.8$ , with ${\chi}^{2}=21.8$ for 19 degrees of freedom. Note that $\sum}\text{\hspace{0.05em}}{m}_{\nu$ is largely degenerate with the remaining parameters in Equation (5), unless we are able to constrain ${k}_{\text{eq}}$ .

Figure 3. Comparison of ${P}_{\text{gal}}\left(k\right)$ obtained from the SDSS-III BOSS survey [4] (“reconstructed”) with ${b}^{2}P\left(k\right)$ obtained from a fit of Equation (6) to the Sachs-Wolfe effect, ${\sigma}_{8}$ , and ${P}_{\text{gal}}\left(k\right)$ . The fit obtains $\sum}\text{\hspace{0.05em}}{m}_{\nu}=0.014\pm 0.079$ eV, ${N}^{2}=\left(1.41\pm 0.12\right)\times {10}^{-10}$ , and ${b}^{2}=1.7\pm 0.1$ , with ${\chi}^{2}=47$ for 20 degrees of freedom (so the statistical uncertainties shown need to be multiplied by $\sqrt{47/20}$ ). Also shown is the fit with $\sum}\text{\hspace{0.05em}}{m}_{\nu}=0.6$ eV fixed for comparison.

Figure 4. Comparison of ${P}_{\text{gal}}\left(k\right)$ obtained from the SDSS-III BOSS survey [4] (“standard”) with ${b}^{2}P\left(k\right)$ obtained from a fit of Equation (6) to the Sachs-Wolfe effect, ${\sigma}_{8}$ , and ${P}_{\text{gal}}\left(k\right)$ . The fit obtains $\sum}\text{\hspace{0.05em}}{m}_{\nu}=0.163\pm 0.061$ eV, ${N}^{2}=\left(1.56\pm 0.12\right)\times {10}^{-10}$ , and ${b}^{2}=2.2\pm 0.2$ , with ${\chi}^{2}=33.9$ for 20 degrees of freedom (so the statistical uncertainties shown need to be multiplied by $\sqrt{33.9/20}$ ). Also shown is the fit with $\sum}\text{\hspace{0.05em}}{m}_{\nu}=0.6$ eV fixed for comparison.

${P}^{\prime}\left(k\right)\equiv \frac{{\left(\text{2\pi}\right)}^{3}{a}_{20}^{2}A{w}^{n}}{{\left(1+\eta w+{w}^{2}\right)}^{2}}\mathrm{,}$ (1)

where
$w\equiv k/{k}_{\text{eq}}$ . Unless otherwise noted, we take the Harrison-Zel’dovich index n = 1 which is close to observations. The parameters A,
$\eta $ , and
${k}_{\text{eq}}$ , as well as the normalization factor b^{2}, are free in the fit. The uncertainties of two data points that fall on BAO peaks are multiplied by three (since BAO is not included in
${P}^{\prime}\left(k\right)$ ).

Also shown in Figure 1 is the suppression of $P\left(k\right)$ for k greater than

${k}_{\text{nr}}=0.018\cdot {\Omega}_{m}^{1/2}{\left(\frac{{\displaystyle \sum}\text{\hspace{0.05em}}{m}_{\nu}}{1\text{\hspace{0.17em}}\text{eV}}\right)}^{1/2}h\text{\hspace{0.17em}}{\text{Mpc}}^{-1}$ (2)

due to free-streaming of massive neutrinos that can not cluster on these small scales, and, more importantly, to the slower growth of structure with massive neutrinos [9] . The suppression factor for $k\gg {k}_{\text{nr}}$ for one massive neutrino, or three degenerate massive neutrinos, is

$f\left(k,{\displaystyle \sum}\text{\hspace{0.05em}}{m}_{\nu}\right)\equiv \frac{P{\left(k\right)}^{{f}_{\nu}}}{P{\left(k\right)}^{{f}_{\nu}=0}}=1-8{f}_{\nu},$ (3)

where ${f}_{\nu}={\Omega}_{\nu}/{\Omega}_{m}$ [9] . ${\Omega}_{m}$ is the total (dark plus baryonic plus neutrino) matter density today relative to the critical density, and includes the contribution ${\Omega}_{\nu}$ of neutrinos that are non-relativistic today. ${\Omega}_{\nu}={h}^{-2}{\displaystyle \sum}\text{\hspace{0.05em}}{m}_{\nu}/93.04$ eV for three left-handed plus right-handed Majorana neutrino eigenstates, or three eigenstates of left-handed Dirac neutrinos plus three right-handed Dirac anti-neutrinos, that are non-relativistic today (right-handed Dirac neutrinos and left-handed Dirac anti-neutrinos are assumed to not have reached thermal and chemical equilibrium with the Standard Model particles). We take $f\left(k,{\displaystyle \sum}\text{\hspace{0.05em}}{m}_{\nu}\right)=1$ for $k<{k}_{\text{nr}}/0.604$ , and

$f\left(k,{\displaystyle \sum}\text{\hspace{0.05em}}{m}_{\nu}\right)=1-0.407\frac{{\displaystyle \sum}\text{\hspace{0.05em}}{m}_{\nu}}{0.6\text{\hspace{0.17em}}\text{eV}}\left[1-{\left(\frac{{k}_{\text{nr}}}{0.604k}\right)}^{0.494}\right]$ (4)

for $k>{k}_{\text{nr}}/0.604$ and $\sum}\text{\hspace{0.05em}}{m}_{\nu}<1.1$ eV, for galaxy formation at a redshift $z=0.5$ [9] .

Figure 2 is the same as Figure 1 except that the function

$P\left(k\right)={P}^{\prime}\left(k\right)f\left(k\mathrm{,}{\displaystyle \sum}\text{\hspace{0.05em}}{m}_{\nu}\right)$ (5)

with $\sum}\text{\hspace{0.05em}}{m}_{\nu}=0.6$ eV is fit. We see that the parameter $\sum}\text{\hspace{0.05em}}{m}_{\nu$ is largely degenerate with the parameters A, $\eta $ , and ${k}_{\text{eq}}$ , so that only a weak sensitivity to $\sum}\text{\hspace{0.05em}}{m}_{\nu$ is obtained unless we are able to constrain ${k}_{\text{eq}}$ . The power spectrum ${P}^{\prime}\left(k\right)$ of Equation (1) neglects the growth of structure inside the horizon while radiation dominates.

The fits in Figure 3 and Figure 4 make full use of the ΛCDM theory. The fitted function is

(6)

where is given by [2] :

(7)

with

(8)

C is a function of, and we take [2] . is a function given in Reference [2] . The pivot point is chosen to not upset Equation (41) below. The fit depends on h, , and the spectral index n, so we define [3] , [7] , and [3] , and obtain, tentatively,

(9)

by minimizing the with respect to, and. The statistical uncertainty has been multiplied by the square root of the per degree of freedom. This result corresponds to the “reconstructed” data in Figure 3. The systematic uncertainties included are from the top-hat window function instead of the gaussian window function, and an alternative value of (details will be given in Section 3). Not included is the systematic uncertainty due to the possible scale dependence of the galaxy bias b.

To obtain, we would like to measure the density at redshift z, but we only have information on the peaks of that have gone non-linear collapsing into visible galaxies. How accurate is the measurement of with galaxies? The measurement of in Ref. [4] is based on a procedure described in [10] based on “the usual assumption that the galaxies form a Poisson sample [11] of the density field”. In other words, the assumption is that the number density of point galaxies is equal to its expected mean (which depends on the position dependent galaxy selection criteria), modulated by the perturbation of the density field:

(10)

Both sides of this equation are measured or calculated at the same length scale, and at the same time. The “galaxy bias” b is explicitly assumed to be scale invariant. If we choose a region of space such that is constant, then the galaxy power spectrum (derived from) should be proportional, under the above assumption, to the power spectrum of linear density perturbations (derived from) up to corrections:

(11)

It is due to this bias b that we have freed the normalization of the measured in the fits corresponding to Figures 1-4.

In the following Sections we study galaxy distributions with SDSS DR14 data and with simulations, in order to understand their connection with the underlying power spectrum of linear density fluctuations. In the end we return to the measurement of neutrino masses.

2. The Hierarchical Formation of Galaxies

This Section allows a precise definition of, and an understanding of the connection between and. We generate galaxies as follows (see [12] for full details). The evolution of the Universe in the homogeneous approximation is described by the Friedmann equation

(12)

The expansion parameter has been normalized to 1 at the present time:. has been normalized so. Therefore is the present Hubble expansion rate. With these normalizations we have. The matter density is, where is the critical density of the Universe. We are interested in the period after the density of matter exceeds the density of radiation. For our simulations we assume flat space, i.e., we neglect the radiation density, take constant [7] , and the present Hubble expansion rate with [3] . The solution to Equation (12) with these parameters is shown by the curve “” in Figure 5. The present age of the universe with these parameters is Gyr.

Setting we obtain the critical universe with expansion parameter

(13)

also shown in Figure 5. We note that. Let us now add density fluctuations to this critical universe and consider a density peak. The growing mode for this density peak is obtained by adding a negative to the critical Universe. This prescription is exact if the density peak is spherically symmetric. An example with “expansion parameter” is presented in Figure 5. Note that grows to maximum expansion and then collapses to zero at time, and, in our model [12] , a galaxy forms. In the example of Figure 5 the galaxy forms at redshift. is the linear approximation to. In the linear approximation for growing modes the density fluctuations relative to grow in proportion to:

(14)

while. At the time, when the galaxy forms, in the linear approximation (which has already broken down).

Figure 5. Expansion parameters as a function of time t of four solutions of the Friedmann equation. From top to bottom, a_{1} corresponds to, a_{2} corresponds to, a_{3} is the linear approximation to, while a_{4} is the exact solution to. a_{4} is the exact solution for the growing mode of a spherically symmetric density peak that collapses to a galaxy at t_{4}. In all cases.

In the linear approximation the density due to Fourier components of wavevector is

(15)

where

(16)

are random phases. The sum of the Fourier series is over comoving wavevectors that satisfy periodic boundary conditions in a rectangular box of volume:

(17)

where, , , , and

(18)

where, and.

Inverting Equation (16) obtains

(19)

where is the comoving coordinate in the linear approximation. The power spectrum of density fluctuations

(20)

is defined in the linear approximation corresponding to, and is approximately independent of V for large V. Averaging over in a bin of obtains. Note that

(21)

Each term in this equation is approximately independent of V. The Fourier transform of the power spectrum is the correlation function:

(22)

The generation of galaxies at time t proceeds as follows. We start with, calculate, and search for local maximums of inside a comoving volume. If a maximum exceeds 1.69 we generate a galaxy of radius

(23)

and dark plus baryonic plus neutrino mass

(24)

if it “fits”, i.e. if it does not overlap previously generated galaxies. I is increased by 1 unit to generate galaxies of a smaller generation, until is reached. See Figure 6.

The peculiar velocity of the generated galaxies is

(25)

and their peculiar displacement is

(26)

is the proper coordinate of a galaxy at the time t of its generation, and. The comoving coordinate of this galaxy, i.e. its position extrapolated to the present time, is the corresponding. causes the difference between the data points in Figure 3 and Figure 4 at large k.

Note in Figure 6 that the formation of galaxies is hierarchical: small galaxies form first, and, as time goes on, density perturbations grow, and groups of galaxies coalesce into larger galaxies in an ongoing process until dark energy dominates and the hierarchical formation of galaxies comes to an end. The distribution of galaxies of generation I depend only on for. Also,

Figure 6. The hierarchical formation of galaxies [12] . Three Fourier components of the density in the linear approximation are shown. Note that in the linear approximation. When reaches 1.69 in the linear approximation the exact solution diverges and a galaxy forms. As time goes on, density perturbations grow, and groups of galaxies of one generation coalesce into larger galaxies of a new generation as shown on the right.

luminous galaxies occupy a total volume (luminous plus dark) less than 1/2.69 of space.

Neutrinos with eV become non-relativistic after the densities of radiation and matter become equal, as illustrated in Figure 7.

3. Fluctuation Amplitude σ_{8}

is the root-mean-square fluctuation of total mass relative to the mean in randomly placed volumes of radius Mpc. We use a “gaussian window function”

(27)

which smoothly defines a volume

(28)

Note that

(29)

The Fourier transform of is

(30)

Then

Figure 7. Example with eV for each of 3 active neutrino eigenstates. Neutrinos become non-relativistic at. The matter density relative to the critical density is for, and for. The densities of matter and radiation become equal at.

(31)

An alternative window function is the “top hat” function for, and for. Then

(32)

Direct measurements obtain [3]

(33)

80% of is due to k/h in the range 0.05 to 0.25 Mpc^{−1}. For comparison, from the 6-parameter ΛCDM fit [3] ,.

4. The Sachs-Wolfe Effect

The spherical harmonic expansion of the CMB temperature fluctuation is

(34)

Averaging over m obtains. The variable that is measured is [2]

(35)

For the dominant contribution to is from the Sachs-Wolfe effect [1] [2] [3] . This range corresponds to. The Sachs-Wolfe effect relates temperature fluctuations of the CMB to perturbations of the gravitational potential [2] :

(36)

When expressed as a function of comoving coordinates, is independent of time when matter dominates. The primordial power spectrum of gravitational potential fluctuations is assumed to have the form [2]

(37)

The relation between and is [2] . In the present analysis, unless otherwise stated, we assume the Harrison-Zel’dovich power spectrum with n = 1, which is close to observations [3] . For, [2]

(38)

where the “quadrupole moment” Q is measured to be

(39)

from the 1996 COBE results (see list of references in [2] ). Then, for,

(40)

and for,

(41)

independently of. Detailed integration obtains results within 10% for.

5. Data and Simulations

The data are obtained from the publicly available SDSS DR14 catalog [5] [6] , see acknowledgement. We consider objects classified as GALAXY, with redshift z in the range 0.4 to 0.6, with redshift error, passing quality selection flags. We further select galaxies in the northern galactic cap, in a “rectangular” volume with Mpc along the line of sight (corresponding to redshift), Mpc (corresponding to an angle 86^{0} across the sky), and Mpc (corresponding to an angle 32^{0}). In total 222470 galaxies pass these selections. The distributions of these galaxies are shown in Figure 8.

Unless otherwise specified, the simulations have Mpc, , , , , , and the input power spectrum of density fluctuations is (5) with Mpc^{3}, Mpc^{−1}, , and eV. We generate galaxies at redshift, corresponding to Gyr, and. This reference simulation has 34,444 galaxies, which is near the limit we can generate with available computing resources.

Figure 8. Distributions of 222470 SDSS DR14 galaxies in a “rectangular” box of dimensions Mpc along the line of sight (corresponding to redshift), Mpc (corresponding to an angle 86^{0} across the sky), and Mpc (corresponding to an angle 32^{0}).

Some definitions are in order. For data we define the absolute red magnitude of a galaxy MAGr at redshift z as the SDSS DR14 variable -modelMag_r corrected to the reference redshift 0.35. Similarly, we define the absolute green magnitude of a galaxy MAGg at redshift z as the SDSS DR14 variable -modelMag_g corrected to the reference redshift 0.35. For a simulated galaxy we define the absolute magnitude MAG, where M is defined by Equation (24). Note that MAGr and MAGg are derived from observed luminosities, while MAG is derived from the total (baryonic plus dark plus neutrino) mass of the simulation. These quantities can only be compared if the luminosity-to-mass ratio is known.

The number of galaxies per unit volume depends on the limiting magnitude of the survey, or on of the simulation.

6. Distributions of Galaxies in SDSS DR14 Data and in Simulation

We would like to obtain from Equation (19) and Equation (20). Unfortunately we do not have access to the relative density fluctuation. Instead we have access to the positions of galaxies and their luminosities. The relation between luminosity and mass of galaxies depends on many variables and is largely unknown, so we focus on the information contained in the positions of galaxies.

Let

(42)

be the number density of point galaxies at redshift z as a function of the comoving coordinate. We have applied periodic boundary conditions in a comoving volume, so has the discrete values of Equation (17). is real, so. The number of galaxies in V is. To invert Equation (42), we multiply it by, integrate over V, and obtain a sum over galaxies j:

(43)

The first term on the right hand side of Equation (43) is the result of a coherent sum of terms corresponding to mode. The second term is the result of an incoherent sum which we have approximated to, where the phase is arbitrary. We define the “galaxy power spectrum”

(44)

and obtain

(45)

The transition between signal and noise occurs at for our data sample, and ≈3.49 for our reference simulation. To test these ideas we can select a narrow range of MAGr, MAGg, or MAG to shift the noise upwards, compare Figures 9-11 (which plot the first term on the right hand side of Equation (45) and include the noise at large k).

Averaging over in a bin of obtains. The factor V is inserted so that becomes independent of the arbitrary choice of V for large V. The function defines statistically the distribution of galaxies. The variables in Equation (16) and Equation (45) should not be confused: there is not necessarily a one-to-one relation between them.

Results for data are presented in Figures 9-11. We note that the galaxy bias b depends on MAGr and MAGg. Even tho, at small k. For this reason in Figure 9 extends to higher k than in Figure 10 and Figure 11 before saturating with noise. Figure 12 presents the noise subtracted galaxy power spectrum, obtained from Figure 9, compared with calculated with the indicated parameters. Their ratio is the bias b^{2}.

Results for the simulations are presented in Figures 13-15. In Figure 15 we compare the reference simulation with, with simulations with (“steeper slope”), or (“less slope”). Note that the function varies between ≈1.3 to ≈0.5 in the region of interest. We observe, qualitatively, that the slope of has a larger effect on than the amplitude A. A comparison of the

Figure 9. Galaxy power spectrum (plus noise visible at large k) from SDSS DR14 data in a volume, at redshift. The fit is with for 29 degrees of freedom, where x and y are the axis in this figure.

Figure 10. Galaxy power spectrum (plus noise visible at large k) in bins of MAGr from SDSS DR14 galaxies with redshift.

Figure 11. Galaxy power spectrum (plus noise visible at large k) in bins of MAGg from SDSS DR14 galaxies with redshift.

Figure 12. Noise subtracted galaxy power spectrum, obtained from Figure 9, compared with calculated with the indicated parameters. Their ratio is b^{2}.

Figure 13. Galaxy power spectrum (plus noise visible at large k) from simulations with three amplitudes A. All other parameters of the simulation are given in Section 5.

Figure 14. Galaxy power spectrum (plus noise visible at large k) from simulations with eV, with three amplitudes A. Other parameters are, and.

simulations in Figure 15 with from data in Figure 9 favors a power spectrum “steeper” than in the reference simulation. The reference simulation has parameters of similar to the ones obtained from the fit in Figure 1 which assumes scale invariant b, and eV. The reference simulation is also similar to the fit “eV” in Figure 12 (taken from Figure 3 which assumes scale invariant b). A steeper implies as shown in Figure 12 by the curve “eV”, and corresponds to a bias b with positive slope as in Equation (55) below.

Figure 15. Galaxy power spectrum (plus noise visible at large k) of the reference simulation, a simulation with with steeper slope (of Equation (1) is multiplied by), and a simulation with less slope (of Equation (1) is divided by).

7. Luminosity and Mass Distributions of Galaxies

Distributions of MAGr and MAGg from data, and MAG from several simulations are presented in Figure 16 and Figure 17. From these figures it is possible to obtain the “mean” luminosity-to-mass ratios. We note that these figures do not show useful sensitivity to.

8. Test of Scale Invariance of the Galaxy Bias b

In this Section we test the scale invariance of the bias b defined in Equation (11). To do so, we count galaxies in an array of spheres of radii, and obtain their mean, and their root-mean-square (rms). All spheres have their center at redshift to ensure the homogeneity of the galaxy selections. The results for, and 256/h Mpc are presented in Table 1. The (rms)^{2} has a contribution from, and a contribution from statistical fluctuations:

(46)

We compare obtained from galaxy counts, with the relative mass fluctuations obtained from Equation (6) and Equation (31). The ratio of these two quantities divided by a correction factor

[2] is the bias b.

Figure 16. Distributions of MAGr and MAGg of SDSS DR14 data, and distributions of MAG of several simulations (see definitions in Section 5). The difference between the MAGr or MAGg of data and MAG of simulations determines the “mean” galaxy L/M ratio.

Figure 17. Same as Figure 16 with additional simulations.

The measured bias b is a function of, , h and the spectral index n. Results for and are presented in Table 1. The last column is the of the five b’s of spheres with, assuming these b’s are scale invariant with respect to their weighted average. Additional measurements of are presented in Figure 18. Assuming that b is scale invariant we obtain

(47)

with minimum for four degrees of freedom. We have defined, and.

Table 1. Mean galaxy counts in spheres of radius. All spheres have their center at redshift. The number of spheres is. Note that the observed root-mean-square (rms) fluctuation relative to the mean is larger than the corresponding statistical fluctuation, i.e.. is calculated with Equation (6) and Equation (31) with N^{2} chosen so to set the scale for b (e.g. for eV, or for eV). Both the galaxy counts and are obtained with the top-hat window function. The true standard deviation is obtained from. The measured “bias” is defined as. The last column is the of the five b’s of spheres with, assuming these b’s are scale invariant., , and.

In conclusion, the galaxy bias b is scale invariant within the statistical uncertainties of b presented in Table 1, provided satisfies Equation (47), else scale invariance is broken. Note in Table 1 that the variation of b with scale depends on.

9. Measurement of Neutrino Masses with the Sachs-Wolfe Effect and σ_{8}

We return to the measurement of neutrino masses. Since the galaxy bias b may be scale dependent, in this Section we exclude measurements of with galaxies.

Figure 18. of five measurements of bias b assumed to be scale invariant with respect to their weighted mean as a function of, for several values of the Hubble parameter h, and the spectral index n. The five measurements of b correspond to scales, and 256/h Mpc.

The ΛCDM model is described by Equation (6) that has three free parameters: N^{2}, n, and. We keep n fixed. We vary the two parameters N^{2} and to minimize a with two terms corresponding to two observables: the Sachs-Wolfe effect (N^{2} from Equation (41)), and given by Equation (33). We therefore have zero degrees of freedom. The result is a function of h, , and the spectral index n, so we define [3] , [7] , and [3] , and obtain

(48)

Note that in the “6 parameter ΛCDM fit” [3] , which assumes eV,

. Here, and below, the systematic uncertainties are obtained by repeating the fits with the top-hat window function instead of the gaussian window function for (and for if applicable), and also with obtained with the “6 parameter ΛCDM fit” [3] , instead of from direct measurements, Equation (33).

The fit of Equation (48) is compared with measurements of obtained from the SDSS-III BOSS survey [4] in Figure 19. It is interesting to note that the discrepancy, i.e. the drop of in the range, is also observed in Figure 12.

For comparison, reference [8] obtains and

(49)

Figure 19. Comparison of obtained from the SDSS-III BOSS survey [4] (“reconstructed”) with obtained from a fit of Equation (6) to the Sachs-Wolfe effect and only. The fit obtains eV with zero degrees of freedom. and are fixed.

where, from a study of BAO with SDSS DR13 galaxies. We allow to vary by one standard deviation, i.e. [7] . To combine the independent measurements (48) and (49) we add one more term to the corresponding to the measurement (49), so we now have one degree of freedom. We obtain

(50)

with for one degree of freedom, so the two independent measurements of, Equation (48) and Equation (49), are consistent. Note that the uncertainty of h dominates the uncertainty of in Equation (50).

We now free h and add one term to the corresponding to [3] , and obtain

(51)

(52)

with for one degree of freedom. The systematic uncertainties in Equation (51) now include. The 1σ, 2σ, and 3σ contours are presented in Figure 20.

If instead we set from the direct measurement of the Hubble expansion rate [3] , we obtain

(53)

Figure 20. Contours corresponding to 1, 2, and 3 standard deviations in the plane, from Sachs-Wolfe, , , and BAO measurements. Points on the contours have, and 9, respectively, where has been minimized with respect to N^{2}. The total uncertainty of is dominated by the uncertainty of h. In this figure n = 1, and the systematic uncertainties, presented in Equation (51), are not included.

(54)

with for 1 degree of freedom. The corresponding 1σ, 2σ, and 3σ contours are presented in Figure 21. Note that the fitted h does not change significantly.

10. Measurement of Neutrino Masses with the Sachs-Wolfe Effect, σ_{8}, and P_{gal}(k)

We repeat the fit of Figure 3, which includes the “reconstructed” SDSS-III BOSS measurements [4] , but this time we allow the galaxy bias b to depend on scale:. Minimizing the with respect to, N^{2}, n, , , and, we obtain

(55)

with for 18 degrees of freedom. The uncertainties have been multiplied by. Confidence contours are presented in Figure 22. Fixing

Figure 21. Same as Figure 20 but.

Figure 22. Contours corresponding to 1, 2, and 3 standard deviations in the plane, from Sachs-Wolfe, , , and measurements. Points on the contours have, and 9, respectively, where has been minimized with respect to N^{2}, n, , and.

obtains, so including the scale dependence of b is necessary.

11. Measurement of Neutrino Masses with the Sachs-Wolfe Effect, σ_{8}, and Galaxy Fluctuations

We repeat the measurements of of Section 9 but add 4 more experimental constraints: of galaxy counts in spheres of radius , and 128/h Mpc, which are listed in Table 1. Spheres of radius 8/h Mpc were not considered because they have. Spheres of radius 256/h Mpc were excluded because there are only 4 spheres of this radius, and the difference between the rms for the top-hat and gaussian window functions turns out to be large (while consistent results are obtained for the other radii). We add two more parameters to be fit: and which define the bias, with for, and 128/h Mpc, respectively. Note that we do not obtain a good fit with fixed bias, and so have introduced a “bias slope”.

From the Sachs-Wolfe effect, , and the measurements we obtain

(56)

with for 2 degrees of freedom. The variables that minimize the are, , , and. This result may be compared with (48).

Freeing, and keeping fixed, we obtain

(57)

with for 2 degrees of freedom.

Combining with the BAO measurement (49) we obtain

(58)

with for 3 degrees of freedom. The variables that minimize the are, , , and. Freeing, and keeping fixed, we obtain

(59)

with for 3 degrees of freedom.

Finally, freeing n, and minimizing the with respect to, N^{2}, n, , , and, we obtain

(60)

with for 2 degrees of freedom. The parameter correlation coefficients, defined in [3] , are

Note that we have measured the amplitude N^{2} and spectral index n of, and the bias including its slope for the SDSS DR14 galaxy selections at redshift z = 0.5. 1, 2, 3, and 4 standard deviation contours are presented in Figure 23.

Figure 23. Contours corresponding to 1, 2, 3, and 4 standard deviations in the plane, from Sachs-Wolfe, , , BAO, and measurements. Points on the contours have, and 16, respectively, where has been minimized with respect to N^{2}, n, , and. The total uncertainty of is dominated by the uncertainty of h. In this figure the systematic uncertainties, presented in Equation (60), are not included.

Figure 23 and Equation (60) are our final results.

12. Conclusions

We have studied galaxy distributions with Sloan Digital Sky Survey SDSS DR14 data and with simulations searching for variables that can constrain neutrino masses. Fitting the predictions of the ΛCDM model to the Sachs-Wolfe effect, , , fluctuations of galaxy counts in spheres of radii ranging from 16/h to 128/h Mpc, BAO measurements, and, in various combinations, with free spectral index n, and free galaxy bias and galaxy bias slope, we obtain consistent measurements of. The uncertainty of is dominated by the uncertainty of h, so we have presented confidence contours in the plane.

Fitting the predictions of the ΛCDM model to the Sachs-Wolfe effect and we obtain (48). Fitting the predictions of the ΛCDM model to the Sachs-Wolfe effect, , and galaxy number fluctuations in spheres of radius, and 128/h, we obtain (56). These results are consistent with the measurement (49) with BAO. Combining these last two independent measurements we obtain

(61)

Note that the uncertainty of is dominated by the uncertainty of h. A global fit with obtains eV, , and the amplitude and spectral index of:, and. The fit also returns the galaxy bias b including its scale dependence.

Figure 23 and Equation (60) are our final results. These results follow from the data analyzed and the assumptions of the validity of the ΛCDM model and. The measured is anticorrelated with h. All steps in this analysis have been fully described.

Note added in proof: Let us comment on Equations (49) and (56). Equation (49) is mainly determined by the precise measurement of the sound horizon angle by the Planck experiment, and by the assumption that the BAO wave stalls at redshift. Equation (49) tells us that lies on the diagonal shown in Figure 23 (with some uncertainty from). Equation (56) is a constraint mainly between and n with large uncertainties. To determine we need as input a value for h (or a value for n). In this article we have taken from [3] . If we obtain eV, and. If however we obtain eV, and. And if, we obtain eV. Alternatively, if we fix, then and eV. Or if we fix [3] , then and eV. At the Guadeloupe 2018 Conference, Adam Riess, representing the SH_{0}ES Team, presented the latest direct measurement of the expansion parameter:, which corresponds to negative! The solution may come from an unexpected direction: gravitational waves from merging black holes are a “standard siren”. The single black hole merger GW170817 already obtains, see the talk by Archil Kobakhidze!

Acknowledgements

Funding for the Sloan Digital Sky Survey (SDSS) has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Aeronautics and Space Administration, the National Science Foundation, the U.S. Department of Energy, the Japanese Monbukagakusho, and the Max Planck Society. The SDSS Web site is http://www.sdss.org/.

The SDSS is managed by the Astrophysical Research Consortium (ARC) for the Participating Institutions. The Participating Institutions are The University of Chicago, Fermilab, the Institute for Advanced Study, the Japan Participation Group, The Johns Hopkins University, Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, University of Pittsburgh, Princeton University, the United States Naval Observatory, and the University of Washington.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Sachs, R.K. and Wolfe, A.M. (1967) Perturbations of a Cosmological Model and Angular Variations of the Microwave Background. The Astrophysical Journal, 147, 73.

https://doi.org/10.1086/148982

[2] Steven Weinberg (2008) Cosmology. Oxford University Press, New York.

[3] Patrignani, C., et al. (2016) Review of Particle Physics. Chinese Physics C, 40, 100001.

https://doi.org/10.1088/1674-1137/40/10/100001

[4] Anderson, L., et al. (2012) The Clustering of Galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: Baryon Acoustic Oscillations in the Data Release 9 Spectroscopic Galaxy Sample. MNRAS, 427, 3435.

https://doi.org/10.1111/j.1365-2966.2012.22066.x

[5] Blanton, M.R., et al. (2017) Sloan Digital Sky Survey IV: Mapping the Milky Way, Nearby Galaxies, and the Distant Universe. The Astronomical Journal, 154, Article ID: 28.

[6] Dawson, K.S., et al. (2013) The Baryon Oscillation Spectroscopic Survey of SDSS-III. The Astronomical Journal, 145, Article ID: 10.

[7] 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.

https://doi.org/10.4236/ijaa.2017.71002

[8] Hoeneisen, B. (2018) Constraints on Neutrino Masses from Baryon Acoustic Oscillation Measurements. International Journal of Astronomy and Astrophysics, 8, 1-5.

https://doi.org/10.4236/ijaa.2018.81001

[9] Lesgourgues J. and Pastor S. (2006) Massive Neutrinos and Cosmology. Physics Reports, 429, 307.

https://doi.org/10.1016/j.physrep.2006.04.001

[10] Feldman H.A., Kaiser N. and Peacock J.A. (1994) Power-Spectrum Analysis of Three-Dimensional Redshift Surveys. ApJ, 426, 23.

https://doi.org/10.1086/174036

[11] Peebles, P.J.E. (1980) The Large-Scale Structure of the Universe. Princeton University Press.

[12] Hoeneisen, B. (2000) A Simple Model of the Hierarchical Formation of Galaxies. arXiv:astro-ph/0009071