Cold or Warm Dark Matter?: A Study of Galaxy Stellar Mass Distributions ()
1. Introduction
Most current cosmological observations are well described by the cold dark matter ΛCDM model with only six independent parameters, and a few assumptions that are consistent with present observations: flat space, a cosmological constant, and scale invariant adiabatic primordial density perturbations [1]. This economical description of the universe is apparently in agreement with all observations on large scales, but seems to have tensions with some small scale phenomena: the “cusp vs core” problem of spiral galaxies, i.e. simulations obtain a cusp while observations find a core, and the “missing satellite” problem [2]. The ΛCDM model assumes that dark matter has a negligible free-streaming length. However, fits to spiral galaxy rotation curves obtain a non-negligible dark matter free-streaming length [3]. This free-streaming cuts off the power spectrum of linear density perturbations at a comoving wavenumber kfs. Adding this parameter to the ΛCDM model obtains the warm dark matter model (ΛWDM).
We compare the observed galaxy stellar mass distributions in the redshift range
with expectations of the cold and warm dark matter models, and obtain the cut-off wavenumber kfs. The notation and cosmological parameters are as in Reference [1].
The outline of this article is as follows. In Section 2 we obtain predictions, based on the Press-Schechter formalism, of the stellar mass distributions for the cold and warm dark matter models. This formalism is valid only at redshifts
as discussed in Section 3. In Section 4 we present measurements of kfs by comparing predictions with data in the redshift range
. Section 5 verifies the compatibility between predictions and the galaxy with largest observed spectroscopic redshift to date. We close with conclusions.
2. Predictions of the Stellar Mass Distributions
Let
be the power spectrum of linear density perturbations in the cold dark matter ΛCDM model as defined in Reference [4], Equation (8.1.42). k is the comoving wavenumber. If dark matter is warm,
becomes replaced by
, where
is a cut-off factor. The cut-off is due to free-streaming of the warm dark matter particles.
In Reference [3] we consider a step-function cut-off factor. In that approximation, the first galaxies to form have the transition mass
(1)
where
. Galaxies with larger masses form bottom up by hierarchical clustering. Once saturation is reached, galaxies that would have formed with mass
may “not fit”, loose mass to neighboring galaxies, and collapse with mass less than
. These are stripped down galaxies, they populate all masses, and are the only galaxies that form with mass less than
in the step function approximation [3].
In the present article we take
(2)
This smooth cut-off is approximately the Born approximation of the calculation presented in Reference [5]. The true cut-off factor has a longer tail at large k than the Born approximation [5]. To study the effect of the tail, we also consider the cut-off factor
(3)
All figures, except Figure 13, include the tail: its effect is relatively small.
As we shall see in the following, the smooth cut-off results in bottom up hierarchical clustering, as in the ΛCDM model, up to saturation at redshift
, and thereafter seems to become dominated by the generation of stripped down galaxies. Irregular “clumpy galaxies”, that resemble beads on filaments or sheets [6], that are dynamically unstable and break up, may also contribute to the galaxy stellar mass function [6] [7].
The mean of the square of the fractional mass fluctuation in a sphere of comoving radius
(smoothed by a gaussian window function), and mass
, at redshift z, is [4]
(4)
while density perturbations are still linear. For simplicity, we have assumed the cut-off factor (2).
is a correction due to the cosmological constant;
for
, respectively [4]. For
Mpc,
[1] is becoming non-linear at the present time.
fixes the normalization of (4).
The Press-Schechter stellar mass function [8] is obtained from (4) as follows. The mass fraction locked up in halos with mass greater than
at redshift z is identified with the probability that the relative fluctuation of mass
exceeds 1.686:
(5)
where
. Then
is identified with the mass fraction in halos with masses between
and
. This identification is valid so long as the galaxies do not break up, or loose mass to neighboring galaxies, and have time to cluster. The Press-Schechter stellar mass function is then obtained after some algebra, and the inclusion of a “fudge factor” 2 [8], justified in [9]:
(6)
where
(7)
and
. Equation (7) is valid in the spherical collapse approximation. A calculation that takes into account the average ellipticity and prolateness of perturbations, is the ellipsoidal collapse approximation, pioneered by R.K. Sheth and G. Tormen [10] [11], that replaces
by
:
(8)
with
. Good fits to simulations are obtained with
[11]. The factor 0.84 depends on the algorithm used to identify the collapsed halos, e.g. on the “link length” of the “friends-of-friends” algorithm, and also on the simulation volume. We note that Equations (6), (7) and (8), have no free parameters, except kfs.
Figures 1-3 present galaxy stellar mass function calculations for the ΛCDM model, and for ΛWDM with
and
, respectively . We have converted from the halo mass
to the stellar mass
as follows:
[3]. This uncertainty should be kept in mind when comparing the figures with observations.
Figure 1. Calculated stellar mass functions with the Press-Schechter, Ellipsoidal Collapse with
, and Ellipsoidal Collapse with
, approximations, for the ΛCDM model, at several redshifts. These distributions are valid before saturation, i.e. for
.
Figure 2. Calculated stellar mass functions with the Press-Schechter, Ellipsoidal Collapse with
, and Ellipsoidal Collapse with
, approximations, for the ΛWDM model with
, at several redshifts. These distributions are valid before saturation, i.e. for
.
Figure 3. Calculated stellar mass functions with the Press-Schechter, Ellipsoidal Collapse with
, and Ellipsoidal Collapse with
, approximations, for the ΛWDM model with
, at several redshifts. These distributions are valid before saturation, i.e. for
.
3. The Stellar Mass Distribution from SDSS Data
We analyze Sloan Digital Sky Survey (SDSS) data release DR16 [12]. We include all data in the right ascension range 145˚ to 230˚, and declination range 0˚ to 50˚. By eye inspection of each redshift bin of this sky patch, we see only mild extraneous features such as zones with different exposure. The galaxy properties, including stellar mass, stellar age, star formation rate (SFR), and metallicity, are obtained from the photon spectra in filters u, g, r, and i, by several stellar population synthesis (SPS) models. The results that we analyze are placed in the following SDSS DR16 classes: stellarMassFSPSGranWideDust [13], stellarMassStarformingPort [14], stellarMassPCAWiscBC03 [15], and stellarMassPCAWiscM11 [15]. The SPS of these classes are described in the cited references. The galaxy stellar mass distributions for these SPS are presented in Figures 4-7, for several redshift bins. The units are counts per unit
(dex) and unit comoving volume (Mpc3).
is the galaxy stellar mass returned by the SPS. The reduction of the distributions at low mass are due to the relative luminosity threshold of the observations. To obtain the galaxy stellar mass functions it is still necessary to divide by the stellar mass completeness factor (which is over 80% at
, and decreases at higher z [16]).
In Figure 4 we observe mass distributions that increase with redshift z at the high mass end. This top down evolution is also observed by the Dark Energy Survey (DES), see Figure 7 of Reference [17]. If we assume that the mass corresponding to a threshold factor 1/2 scales as the square of the luminosity distance, then the shift of the distributions to the right for
should be even larger.
Figure 4. Galaxy counts per dex per comoving volume
[dex−1Mpc−3] as a function of
from SDSS DR16 data in class stellarMassFSPSGranWideDust [13], in bins of redshift z of width ±0.5 (bin
has width ±0.025).
Figure 5. Galaxy counts per dex per comoving volume
[dex−1Mpc−3] as a function of
from SDSS DR16 data in class stellarMassStarformingPort [14], in bins of redshift z of width ±0.5 (bin
has width ±0.025).
Figure 6. Galaxy counts per dex per comoving volume
[dex−1Mpc−3] as a function of
from SDSS DR16 data in class stellarMassPCAWiscM11 [15], in bins of redshift z of width ±0.5 (bin
has width ±0.025).
Figure 7. Galaxy counts per dex per comoving volume
[dex−1∙Mpc−3] as a function of
from SDSS DR16 data in class stellarMassPCAWiscBC03 [15], in bins of redshift z of width ±0.5 (bin
has width ±0.025).
The top down evolution is observed even when the expected mass is replaced by the median mass minus one standard deviation, so the excess at high mass is not due to a statistical fluctuation. However, Figure 5 presents galaxy stellar mass distributions that do not change significantly with redshift. In Figure 6 and Figure 7 the evolution is slightly top down. In summary, at our current level of understanding, in the redshift range
the galaxy stellar mass function either evolves top down, or is stationary within observational uncertainties.
Let us compare the observed stellar mass function at
, e.g. Figure 4, with the calculations in Figures 1-3. We find that at
the calculations at
already matches the observation at
. This “saturation” at the high mass end is not understood. At
we obtain
, 4 and 2 for
, 1.6 Mpc−1 and 0.8 Mpc−1, respectively. At these
for
the probability
is of order 0.01, stripped down galaxies form, and the Press-Schechter formalism breaks down. Galaxy merging requires dissipation. The “saturation” observed at
may be due to the long time required for “dry” mergers of galaxies with little gas content. In conclusion, to measure
, we need to compare observations with calculations at
, before the saturation sets in.
Note that the predictions become insensitive to
for
. Therefore, to measure
, we verify that prediction and data are in agreement for
. For future convenience,
for
, respectively.
4. Measurements of kfs from Stellar Mass Distributions with z ≈ 5.5 to z ≈ 8.5
Reference [18] presents a compilation of measured stellar mass functions for redshifts
to
, and estimates the systematic uncertainties imposing continuity equation constraints. The measurements with
to
[19] [20] [21] are compared with calculations in Figures 8-10. From these figures we obtain the measurements of kfs summarized in Table 1. Note that the bin centered at
already shows signs of saturation at the high mass end, see Figure 11.
Taking the Ellipsoidal Collapse model with
as the preferred prediction with an uncertainty
(see Table 1), the contribution of correlated systematic uncertainties of the data obtained in Reference [18], ±0.15 Mpc−1,
Table 1. Measurements of the warm dark matter cut-off wavenumber kfs obtained from Figures 8-10, assuming the validity of the Press-Schechter, Ellipsoidal Collapse with
, and Ellipsoidal Collapse with
, approximations. The total uncertainties shown include statistical uncertainties, and systematic uncertainties estimated in [18].
Figure 8. Calculated stellar mass functions with the Press-Schechter, Ellipsoidal Collapse with
, and Ellipsoidal Collapse with
, approximations, for ΛCDM, and ΛWDM with
and 0.8 Mpc−1, at redshift
, compared with observations [18] [19].
Figure 9. Calculated stellar mass functions with the Press-Schechter, Ellipsoidal Collapse with
, and Ellipsoidal Collapse with
, approximations, for ΛCDM, and ΛWDM with
and 0.8 Mpc−1, at redshift
, compared with observations [18] [19] [20].
Figure 10. Calculated stellar mass functions with the Press-Schechter, Ellipsoidal Collapse with
, and Ellipsoidal Collapse with
, approximations, for ΛCDM, and ΛWDM with
and 0.8 Mpc−1, at redshift
, compared with observations [18] [19] [20].
Figure 11. Calculated stellar mass functions with the Press-Schechter, Ellipsoidal Collapse with
, and Ellipsoidal Collapse with
approximations, for ΛCDM, and ΛWDM with
and 0.8 Mpc−1, at redshift
, compared with observations at
[18] [19] [20] [21]. Note the onset of “saturation” at the high mass end (that is not understood).
an uncertainty due to
, ±0.2, and statistical uncertainties, we obtain our final measurement:
. This result is insensitive to the “tail” in (3).
(Note: The present measurement of kfs superceeds the estimate in Reference [3] that was based on data in SDSS DR15, class stellarMassFSPSGranWideDust that shows strong top down galaxy evolution, see Figure 4.)
5. Estimate of kfs from Galaxy GN-z11
The galaxy with largest spectroscopically confirmed redshift to date is GN-z11 with
[22]. Its stellar mass is estimated to be
. One such galaxy was found in a comoving search volume
, for
. Figure 12 compares this single galaxy with expectations corresponding to the cut-off factor (3). To illustrate the effect of the cut-off factor tail, Figure 13 presents the expectations corresponding to the gaussian cut-off factor (2). From this single galaxy we obtain
.
6. Conclusion
Comparing measurements of stellar mass distributions of galaxies in the redshift range
with expectations, we obtain the warm dark matter cut-off wavenumber
. This result is in agreement with the independent measurements obtained by fitting spiral galaxy rotation curves (demonstrating that the cut-off kfs is due to warm dark matter free-streaming), and is consistent with the scenario of dark matter with no freeze-in and no freeze-out, see Table 2 [3] [23] [24] [25]. Detailed properties of warm dark matter can be derived from kfs [3]. The observed stellar mass functions disfavor the ΛCDM model.
Figure 12. Calculated stellar mass functions with the Press-Schechter, Ellipsoidal Collapse with
, and Ellipsoidal Collapse with
, approximations, for ΛCDM, and ΛWDM with
and 0.8 Mpc−1, at redshift
, compared with one observed galaxy GN-z11 (assuming one similar galaxy per dex) [22]. The cut-off factor is given in Equation (3). This graph obtains kfs of the order of 1.1 Mpc−1.
Figure 13. Same as Figure 12, but with the gaussian cut-off factor (2).
Table 2. Update of Table 2 of Reference [3]. Summary of three independent measurements of the adiabatic invariant
[3], the expansion parameter at which dark matter particles become non-relativistic
, the cut-off wavenumber of warm dark matter kfs, the transition galaxy mass
and the mass mh of dark matter particles (for the case of zero chemical potential). The top (bottom) table is for fermions with
(bosons with
).
Acknowledgements
Funding for the Sloan Digital Sky Survey IV has been provided by the Alfred P. Sloan Foundation, the U.S. Department of Energy Office of Science, and the Participating Institutions. SDSS-IV acknowledges support and resources from the Center for High-Performance Computing at the University of Utah. The SDSS web site is http://www.sdss.org.
SDSS-IV is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS Collaboration including the Brazilian Participation Group, the Carnegie Institution for Science, Carnegie Mellon University, the Chilean Participation Group, the French Participation Group, Harvard-Smithsonian Center for Astrophysics, Instituto de Astrofísica de Canarias, The Johns Hopkins University, Kavli Institute for the Physics and Mathematics of the Universe (IPMU)/University of Tokyo, the Korean Participation Group, Lawrence Berkeley National Laboratory, Leibniz Institut für Astrophysik Potsdam (AIP), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Max-Planck-Institut für Astrophysik (MPA Garching), Max-Planck-Institut für Extraterrestrische Physik (MPE), National Astronomical Observatories of China, New Mexico State University, New York University, University of Notre Dame, Observatário Nacional/MCTI, The Ohio State University, Pennsylvania State University, Shanghai Astronomical Observatory, United Kingdom Participation Group, Universidad Nacional Autónoma de México, University of Arizona, University of Colorado Boulder, University of Oxford, University of Portsmouth, University of Utah, University of Virginia, University of Washington, University of Wisconsin, Vanderbilt University, and Yale University.