Dark Matter Mass Function Research Articles

Inelastic Scattering of Dark Matter with Heavy Cosmic Rays

We investigate the impact of inelastic collisions between dark matter (DM) and heavy cosmic ray (CR) nuclei on CR propagation. We approximate the fragmentation cross-sections for DM-CR collisions using collider-measured proton-nuclei scattering cross-sections, allowing us to assess how these collisions affect the spectra of CR boron and carbon. We derive new CR spectra from DM-CR collisions by incorporating their cross-sections into the source terms and solving the diffusion equation for the complete network of reactions involved in generating secondary species. In a specific example with a coupling strength of b χ = 0.1 and a DM mass of m χ = 0.1 GeV, considering a simplified scenario where DM interacts exclusively with oxygen, a notable modification in the boron-to-carbon spectrum due to the DM-CR interaction is observed. Particularly, the peak within the spectrum, spanning from 0.1 to 10 GeV, experiences an enhancement of approximately 1.5 times. However, in a more realistic scenario where DM particles interact with all CRs, this peak can be amplified to twice its original value. Utilizing the latest data from AMS-02 and DAMPE on the boron-to-carbon ratio, we estimate a 95% upper limit for the effective inelastic cross-section of DM-proton as a function of DM mass. Our findings reveal that at m χ ≃ 2 MeV, the effective inelastic cross-section between DM and protons must be less than O(10−32)cm2 .

Blazar Constraints on Neutrino-Dark Matter Scattering.

Neutrino emission in coincidence with gamma rays has been observed from the blazar TXS 0506+056 by the IceCube telescope. Neutrinos from the blazar had to pass through a dense spike of dark matter (DM) surrounding the central black hole. The observation of such a neutrino implies new upper bounds on the neutrino-DM scattering cross section as a function of DM mass. The constraint is stronger than existing ones for a range of DM masses, if the cross section rises linearly with energy. For constant cross sections, competitive bounds are also possible, depending on details of the DM spike.

A global analysis of dark matter signals from 27 dwarf spheroidal galaxies using 11 years of Fermi-LAT observations

We search for a dark matter signal in 11 years of Fermi-LAT gamma-ray data from 27 Milky Way dwarf spheroidal galaxies with spectroscopically measured J-factors. Our analysis includes uncertainties in J-factors and background normalisations and compares results from a Bayesian and a frequentist perspective. We revisit the dwarf spheroidal galaxy Reticulum II, confirming that the purported gamma-ray excess seen in Pass 7 data is much weaker in Pass 8, independently of the statistical approach adopted. We introduce for the first time posterior predictive distributions to quantify the probability of a dark matter detection from another dwarf galaxy given a tentative excess. A global analysis including all 27 dwarfs shows no indication for a signal in nine annihilation channels. We present stringent new Bayesian and frequentist upper limits on the dark matter cross section as a function of dark matter mass. The best-fit dark matter parameters associated with the Galactic Centre excess are excluded by at least 95% confidence level/posterior probability in the frequentist/Bayesian framework in all cases. However, from a Bayesian model comparison perspective, dark matter annihilation within the dwarfs is not strongly disfavoured compared to a background-only model. These results constitute the highest exposure analysis on the most complete sample of dwarfs to date. Posterior samples and likelihood maps from this study are publicly available.

Gamma-Ray Sensitivity to Dark Matter Subhalo Modelling at High Latitudes

Searches for “dark” subhaloes in gamma-ray point-like source catalogues are among promising strategies for indirect dark matter detection. Such a search is nevertheless affected by uncertainties related, on the one hand, to the modelling of the dark matter subhalo distribution in Milky-Way-like galaxies, and, on the other hand, to the sensitivity of gamma-ray instruments to the dark matter subhalo signals. In the present work, we assess the detectability of dark matter subhaloes in Fermi-LAT catalogues, taking into accounts uncertainties associated with the modelling of the galactic subhalo population. We use four different halo models bracketing a large set of uncertainties. For each model, adopting an accurate detection threshold of the LAT to dark matter subhalo signals and comparing model predictions with the number of unassociated point-sources in Fermi-LAT catalogues, we derive upper limits on the annihilation cross section as a function of dark matter mass. Our results show that, even in the best-case scenario (i.e., DMonly subhalo model), which does not include tidal disruption from baryons, the limits on the dark matter parameter space are less stringent than current gamma-ray limits from dwarf spheroidal galaxies. Comparing the results obtained with the different subhalo models, we find that baryonic effects on the subhalo population are significant and lead to dark matter constraints that are less stringent by a factor of ∼2 to ∼5. This uncertainty comes from the unknown resilience of dark matter subhaloes to tidal disruption.

Cold dark energy constraints from the abundance of galaxy clusters

We constrain cold dark energy of negligible sound speed using galaxy cluster abundance observations. In contrast to standard quasi-homogeneous dark energy, negligible sound speed implies clustering of the dark energy fluid at all scales, allowing us to measure the effects of dark energy perturbations at cluster scales. We compare those models and set the stage for using non-linear information from semi-analytical modelling in cluster growth data analyses. For this, we recalibrate the halo mass function with non-linear characteristic quantities, the spherical collapse threshold and virial overdensity, that account for model and redshift dependent behaviours, as well as an additional mass contribution for cold dark energy. We present the first constraints from this cold dark matter plus cold dark energy mass function using our cluster abundance likelihood, which self-consistently accounts for selection effects, covariances and systematic uncertainties. We combine cluster growth data with CMB, SNe Ia and BAO data, and find a shift between cold versus quasi-homogeneous dark energy of up to $1\sigma$. We make a Fisher matrix forecast of constraints attainable with cluster growth data from the on-going Dark Energy Survey (DES). For DES, we predict $\sim$50$\%$ tighter constraints on $\left(\Omega_\mathrm{m},w \right)$ for cold dark energy versus $w$CDM models, with the same free parameters. Overall, we show that cluster abundance analyses are sensitive to cold dark energy, an alternative, viable model that should be routinely investigated alongside the standard dark energy scenario.

Indirect dark matter signatures in the cosmic dark ages. I. Generalizing the bound ons-wave dark matter annihilation from Planck results

Recent measurements of the cosmic microwave background (CMB) anisotropies by Planck provide a sensitive probe of dark matter annihilation during the cosmic dark ages, and specifically constrain the annihilation parameter $f_\mathrm{eff} \langle \sigma v \rangle/m_\chi$. Using new results (Paper II) for the ionization produced by particles injected at arbitrary energies, we calculate and provide $f_\mathrm{eff}$ values for photons and $e^+e^-$ pairs injected at keV-TeV energies; the $f_\mathrm{eff}$ value for any dark matter model can be obtained straightforwardly by weighting these results by the spectrum of annihilation products. This result allows the sensitive and robust constraints on dark matter annihilation presented by the Planck Collaboration to be applied to arbitrary dark matter models with $s$-wave annihilation. We demonstrate the validity of this approach using principal component analysis. As an example, we integrate over the spectrum of annihilation products for a range of Standard Model final states to determine the CMB bounds on these models as a function of dark matter mass, and demonstrate that the new limits generically exclude models proposed to explain the observed high-energy rise in the cosmic ray positron fraction. We make our results publicly available at http://nebel.rc.fas.harvard.edu/epsilon.

THE COSMIC EVOLUTION OF FAINT SATELLITE GALAXIES AS A TEST OF GALAXY FORMATION AND THE NATURE OF DARK MATTER

The standard cosmological model based on cold dark matter (CDM) predicts a large number of subhalos for each galaxy-size halo. Matching the subhalos to the observed properties of luminous satellites of galaxies in the local universe poses a significant challenge to our understanding of the astrophysics of galaxy formation. We show that the cosmic evolution and host mass dependence of the luminosity function of satellites provide a powerful new diagnostic to disentangle astrophysical effects from variations in the underlying dark matter mass function. We illustrate this by comparing recent observations of satellites between redshifts 0.1 < z < 0.8 based on Hubble Space Telescope images, with predictions from three different state-of-the-art semi-analytic models applied to CDM power spectra, with one model also applied to a warm dark matter (WDM) spectrum. We find that even though CDM models provide a reasonable fit to the local luminosity function of satellites around galaxies comparable to the Milky Way, they do not reproduce the data as well for different redshifts and host galaxy stellar masses, indicating that further improvements in the description of star formation are likely needed. The WDM model matches the observed mass dependence and redshift evolution of satellite galaxies more closely, indicating that a modification of the underlying power spectrum may offer an alternative solution to this tension. We conclude by presenting predictions for the color distribution of satellite galaxies to demonstrate how future observations will be able to further distinguish between these models and to help constrain baryonic and non-baryonic physics.

THE ARECIBO LEGACY FAST ALFA SURVEY. X. THE H I MASS FUNCTION AND $\Omega _{\rm H\,{\mathsc{i}}}$ FROM THE 40% ALFALFA SURVEY

The Arecibo Legacy Fast ALFA (ALFALFA) survey has completed source extraction for 40% of its total sky area, resulting in the largest sample of HI-selected galaxies to date. We measure the HI mass function from a sample of 10,119 galaxies with 6.2 < log (M_HI/M_Sun) < 11.0 and with well-described mass errors that accurately reflect our knowledge of low-mass systems. We characterize the survey sensitivity and its dependence on profile velocity width, the effect of large-scale structure, and the impact of radio frequency interference in order to calculate the HIMF with both the 1/Vmax and 2DSWML methods. We also assess a flux-limited sample to test the robustness of the methods applied to the full sample. These measurements are in excellent agreement with one another; the derived Schechter function parameters are phi* = 4.8 (+/- 0.3) * 10^-3, log (M*/M_Sun) + 2 log(h_70) = 9.96 (+/- 0.2), and alpha = -1.33 (+/- 0.02). We find Omega_HI = 4.3 (+/- 0.3) * 10^-4, 16% larger than the 2005 HIPASS result, and our Schechter function fit extrapolated to log (M_HI/M_Sun) = 11.0 predicts an order of magnitude more galaxies than HIPASS. The larger values of Omega_HI and of M* imply an upward adjustment for estimates of the detection rate of future large-scale HI line surveys with, e.g., the Square Kilometer Array. A comparison with simulated galaxies from the Millennium Run and a treatment of photoheating as a method of baryon removal from HI-selected halos indicates that the disagreement between dark matter mass functions and baryonic mass functions may soon be resolved.

Spherical collapse in quintessence models with zero speed of sound

We study the spherical collapse model in the presence of quintessence with negligible speed of sound. This case is particularly motivated for w < −1 as it is required by stability. As pressure gradients are negligible, quintessence follows dark matter during the collapse. The spherical overdensity behaves as a separate closed FLRW universe, so that its evolution can be studied exactly. We derive the critical overdensity for collapse and we use the extended Press-Schechter theory to study how the clustering of quintessence affects the dark matter mass function. The effect is dominated by the modification of the linear dark matter growth function. A larger effect occurs on the total mass function, which includes the quintessence overdensities. Indeed, here quintessence constitutes a third component of virialized objects, together with baryons and dark matter, and contributes to the total halo mass by a fraction ∼ (1+w)ΩQ/Ωm. This gives a distinctive modification of the total mass function at low redshift.

The Evolving Faint End of the Luminosity Function

We investigate the evolution of the faint-end slope of the luminosity function, $\alpha$, using semi-analytical modeling of galaxy formation. In agreement with observations, we find that the slope can be fitted well by $\alpha (z) =a+b z$, with a=-1.13 and b=-0.1. The main driver for the evolution in $\alpha$ is the evolution in the underlying dark matter mass function. Sub-L_* galaxies reside in dark matter halos that occupy a different part of the mass function. At high redshifts, this part of the mass function is steeper than at low redshifts and hence $\alpha$ is steeper. Supernova feedback in general causes the same relative flattening with respect to the dark matter mass function. The faint-end slope at low redshifts is dominated by field galaxies and at high redshifts by cluster galaxies. The evolution of $\alpha(z)$ in each of these environments is different, with field galaxies having a slope b=-0.14 and cluster galaxies b=-0.05. The transition from cluster-dominated to field-dominated faint-end slope occurs roughly at a redshift $z_* \sim 2$, and suggests that a single linear fit to the overall evolution of $\alpha(z)$ might not be appropriate. Furthermore, this result indicates that tidal disruption of dwarf galaxies in clusters cannot play a significant role in explaining the evolution of $\alpha(z)$ at z< z_*. In addition we find that different star formation efficiencies a_* in the Schmidt-Kennicutt-law and supernovae-feedback efficiencies $\epsilon$ generally do not strongly influence the evolution of $\alpha(z)$.

The faint end of the galaxy luminosity function

We present and discuss optical measurements of the faint end of the galaxy luminosity function down to M R = -10 in five different local environments of varying galaxy density and morphological content. The environments we studied, in order of decreasing galaxy density, are the Virgo Cluster, the NGC 1407 Group, the Coma I Group, the Leo Group and the NGC 1023 Group. Our results come from a deep wide-angle survey with the National Astronomical Observatories of Japan Subaru 8-m Telescope on Mauna Kea and are sensitive down to very faint surface-brightness levels. Galaxies were identified as group or cluster members on the basis of their surface brightness and morphology. The faintest galaxies in our sample have R ∼ 22.5. There were thousands of fainter galaxies but we cannot distinguish cluster members from background galaxies at these faint limits so do not attempt to determine a luminosity function fainter than M R = -10. In all cases, there are far fewer dwarfs than the numbers of low-mass haloes anticipated by cold dark matter theory. The mean logarithmic slope of the luminosity function between M R = -18 and M R = -10 is α ≃ -1.2, far shallower than the cold dark matter mass function slope of α ≃ -1.8. We would therefore need to be missing about 90 per cent of the dwarfs at the faint end of our sample in all the environments we study to achieve consistency with CDM theory. It is unlikely that such large numbers of dwarfs are missed because (i) the data are deep enough that we are sensitive to very low surface brightness galaxies, and (ii) the seeing is good enough that we can have some confidence in our ability to distinguish high surface brightness dwarfs from background galaxies brighter than R = 22.5. One caveat is that we miss compact members taken to be background galaxies, but such objects (like M32) are thought to be rare.

MASS ESTIMATES FOR GALACTIC DARK MATTER OBJECTS AS A TEST OF A FRACTAL COSMOLOGICAL PARADIGM

Empirical mass estimates for the enigmatic galactic dark matter objects (GDMOs), which account for most of the mass of the universe, but remain unidentified, are compared with the definitive predictions of a discrete fractal model of the cosmos. Trends in the reported mass estimates are identified and uncertainties are discussed. Intriguing similarities among various stellar and dark matter mass functions are noted, and a possible identification of the galactic dark matter objects is suggested.

The Expected Mass Function for Low‐Mass Galaxies in a Cold Dark Matter Cosmology: Is There a Problem?

It is well known that the mass function for halos in cold dark matter (CDM) cosmology is a relatively steep power law [n(M) ∝ M-2] for low masses, possibly too steep to be consistent with observations. But how steep is the galaxy mass function? We have analyzed the stellar and gas mass functions of the first luminous objects formed in a ΛCDM universe, as calculated in the numerical simulation described by Gnedin in 2000. We have found that suppression of the low-mass end of the galaxy mass function emerges naturally, implying that perhaps no new physics beyond the standard model is needed. Although previous pure semianalytic and combined N-body/semianalytic treatments of the galaxy formation problem have suggested this result, the calculation found here is derived from the first fully hydrodynamic simulation with (1) sufficient mass and spatial resolution to individually resolve low-mass galaxies, (2) a physically plausible star formation and UV photon generation algorithm, and (3) allowance for inhomogeneity in radiation emission and absorption. In particular, we find that the stellar mass function is consistently flat at the low-mass end. Moreover, while the gas mass function follows the dark matter mass function until reionization at z ~ 7, between z = 7 and z = 4 the gas mass function also flattens considerably at the low-mass end. At z = 4, the gas and stellar mass functions are fitted by a Schechter function with α ≈ -1.2 ± 0.1, significantly shallower than the dark matter halo mass function and consistent with some recent observations. In addition, the low-mass objects found in the simulation are consistent with the observed relationship between luminosity and velocity dispersion for dwarf spheroidals. The flattening of the mass function is a result of two factors: (1) the dark matter halo mass (differential) function is consistent with the Press-Schechter formulation at low masses n(M) ∝ M-2 and (2) heating/cooling and ionization processes seem to cause baryons to collect in halos with the relationship Mb ∝ M at low masses. Combining these two factors yields n(Mb) ∝ M, comparable to the simulation results. Finally, although the simulation ends at z = 4, processes such as mergers and additional gas ejection can only further suppress the mass function at low masses from z = 4 to the present epoch. If these processes are significant, the fact that the simulation at z = 4 produces a turnover in the gas mass function corresponding approximately to that observed in dwarf galaxies today suggests that ΛCDM could end up underproducing low-mass galaxies!

Probing the Galactic Dark matter mass function using microlensing and direct searches

If compact baryonic objects contribute significantly to the dark matter in our Galaxy, their mass function will present vital clues for galaxy formation theories and star formation processes in the early Universe. Here we discuss what one might expect to learn about the mass function of Galactic dark matter from microlensing and from direct searches in the infrared and optical wavebands. Current microlensing results from the EROS collaboration already constrain halo mass functions which extend below 10 −4 M bd , whilst recent HST observations place strong constraints on disc and halo dark matter mass functions extending above 0.1 M ⊙ . Infrared observations should either detect or constrain objects larger than 0.01 M ⊙ in the near future. Objects below 0.01 M ⊙ should be detectable through microlensing, although the prospects of determining their mass function depend critically on a number of factors.