Beecroft Institute for Particle Astrophysics and Cosmology

Abstract:

The concentrations of dark matter haloes provide crucial information about their internal structure and how it depends on mass and redshift – the so-called concentration–mass–redshift relation, denoted c(M, z). We present here an extensive study of the cosmology-dependence of c(M, z) that is based on a suite of 72 gravity-only, full N-body simulations in which the following cosmological parameters were varied: σ8Mb, ns, h, Mν, w0, and wa. We characterize the impact of these parameters on concentrations for different halo masses and redshifts. In agreement with previous works, and for all cosmologies studied, we find that there exists a tight correlation between the characteristic densities of dark matter haloes within their scale radii, r−2, and the critical density of the universe at a suitably defined formation time. This finding, when combined with excursion set modelling of halo formation histories, allows us to accurately predict the concentrations of dark matter haloes as a function of mass, redshift, and cosmology. We use our simulations to test the reliability of a number of published models for predicting halo concentration and highlight when they succeed or fail to reproduce the cosmological c(M, z) relation.

Abstract:

Decaying or annihilating dark matter particles could be detected through gamma-ray emission from the species they decay or annihilate into. This is usually done by modeling the flux from specific dark matter-rich objects such as the Milky Way halo, Local Group dwarfs, and nearby groups. However, these objects are expected to have significant emission from baryonic processes as well, and the analyses discard gamma-ray data over most of the sky. Here we construct full-sky templates for gamma-ray flux from the large-scale structure within ∼200 Mpc by means of a suite of constrained N-body simulations (csiborg) produced using the Bayesian Origin Reconstruction from Galaxies algorithm. Marginalizing over uncertainties in this reconstruction, small-scale structure, and parameters describing astrophysical contributions to the observed gamma-ray sky, we compare to observations from the Fermi Large Area Telescope to constrain dark matter annihilation cross sections and decay rates through a Markov chain Monte Carlo analysis. We rule out the thermal relic cross section for s-wave annihilation for all mχ7 GeV/c2 at 95% confidence if the annihilation produces gluons or quarks less massive than the bottom quark. We infer a contribution to the gamma-ray sky with the same spatial distribution as dark matter decay at 3.3σ. Although this could be due to dark matter decay via these channels with a decay rate Γ≈6×10-28 s-1, we find that a power-law spectrum of index p=-2.75-0.46+0.71, likely of baryonic origin, is preferred by the data.