No Arabic abstract
We determine the dark matter pair-wise relative velocity distribution in a set of Milky Way-like halos in the Auriga and APOSTLE simulations. Focusing on the smooth halo component, the relative velocity distribution is well-described by a Maxwell-Boltzmann distribution over nearly all radii in the halo. We explore the implications for velocity-dependent dark matter annihilation, focusing on four models which scale as different powers of the relative velocity: Sommerfeld, s-wave, p-wave, and d-wave models. We show that the J-factors scale as the moments of the relative velocity distribution, and that the halo-to-halo scatter is largest for d-wave, and smallest for Sommerfeld models. The J-factor is strongly correlated with the dark matter density in the halo, and is very weakly correlated with the velocity dispersion. This implies that if the dark matter density in the Milky Way can be robustly determined, one can accurately predict the dark matter annihilation signal, without the need to identify the dark matter velocity distribution in the Galaxy.
We calculate the effective $J$-factors, which determine the strength of indirect detection signals from dark matter annihilation, for 25 dwarf spheroidal galaxies (dSphs). We consider several well-motivated assumptions for the relative velocity dependence of the dark matter annihilation cross section: $sigma_A v$: $s$-wave (velocity independent), $p$-wave ($sigma_A v propto v^2$), $d$-wave ($sigma_A v propto v^4$), and Sommerfeld-enhancement in the Coulomb limit ($sigma_A v propto 1/v$). As a result we provide the largest and most updated sample of J-factors for velocity-dependent annihilation models. For each scenario, we use Fermi-LAT gamma-ray data to constrain the annihilation cross section. Due to the assumptions made in our gamma-ray data analysis, our bounds are comparable to previous bounds on both the $p$-wave and Sommerfeld-enhanced cross sections using dSphs. Our bounds on the $d$-wave cross section are the first such bounds using indirect detection data.
We obtain predictions for the properties of cold dark matter annihilation radiation using high resolution hydrodynamic zoom-in cosmological simulations of Milky Way-like galaxies (APOSTLE project) carried out as part of the Evolution and Assembly of GaLaxies and their Environments (EAGLE) programme. Galactic halos in the simulation have significantly different properties from those assumed in the standard halo model often used in dark matter detection studies. The formation of the galaxy causes a contraction of the dark matter halo, whose density profile develops a steeper slope than the Navarro-Frenk-White (NFW) profile between $rapprox1.5$ kpc and $rapprox10$ kpc. At smaller radii, $rlesssim1.5$ kpc, the halos develop a flatter than NFW slope. This unexpected feature may be specific to our particular choice of subgrid physics model but nevertheless the dark matter density profiles agree within 30% as the mass resolution is increased by a factor 150. The inner regions of the halos are almost perfectly spherical (axis ratios $b/a > 0.97$ within $r=1$ kpc) and there is no offset larger than 45 pc between the centre of the stellar distribution and the centre of the dark halo. The morphology of the predicted dark matter annihilation radiation signal is in broad agreement with $gamma$-ray observations at large Galactic latitudes ($bgtrsim3^circ$). At smaller angles, the inferred signal in one of our four galaxies is similar to that which is observed but it is significantly weaker in the other three.
We present the effective $J$-factors for the Milky Way for scenarios in which dark matter annihilation is p-wave or d-wave suppressed. We find that the velocity suppression of dark matter annihilation can have a sizable effect on the morphology of a potential dark matter annihilation signal in the Galactic Center. The gamma-ray flux from the innermost region of the Galactic Center is in particular suppressed. We find that for dark matter density profiles with steep inner slopes, the morphology of the Inner Galaxy gamma-ray emission in p-wave models can be made similar to the morphology in standard s-wave models. This similarity may suggest that model discrimination between s-wave and p-wave is challenging, for example, when fitting the Galactic Center excess. However, we show that it is difficult to simultaneously match s- and p-wave morphologies at both large and small angular scales. The $J$-factors we calculate may be implemented with astrophysical foreground models to self-consistently determine the morphology of the excess with velocity-suppressed dark matter annihilation.
We identify Local Group (LG) analogs in the IllustrisTNG cosmological simulation, and use these to study two mass estimators for the LG: one based on the timing argument (TA) and one based on the virial theorem (VT). Including updated measurements of the Milky Way-M31 tangential velocity and the cosmological constant, we show that the TA mass estimator slightly overestimates the true median LG-mass, though the ratio of the TA to the true mass is consistent at the approximate 90% c.l. These are in broad agreement with previous results using dark matter-only simulations. We show that the VT estimator better estimates the true LG-mass, though there is a larger scatter in the virial mass to true mass ratio relative to the corresponding ratio for the TA. We attribute the broader scatter in the VT estimator to several factors, including the predominantly radial orbits for LG satellite galaxies, which differs from the VT assumption of isotropic orbits. With the systematic uncertainties we derive, the updated measurements of the LG mass at 90% c.l. are $4.75_{-2.41}^{+2.22} times 10^{12}$ M$_odot$ from the TA and $2.0_{-1.5}^{+2.1} times 10^{12}$ M$_odot$ from the VT.
Studies of cluster mass and velocity anisotropy profiles are useful tests of dark matter models, and of the assembly history of clusters of galaxies. These studies might be affected by unknown systematics caused by projection effects. We aim at testing observational methods for the determination of mass and velocity anisotropy profiles of clusters of galaxies. Particularly, we focus on the MAMPOSSt technique (Mamon et al. 2013). We use results from two semi-analytic models of galaxy formation coupled with high-resolution N-body cosmological simulations, the catalog of De Lucia & Blaizot (2007) and the FIRE catalog based on the new GAlaxy Evolution and Assembly model. We test the reliability of the Jeans equation in recovering the true mass profile when full projected phase-space information is available. We examine the reliability of the MAMPOSSt method in estimating the true mass and velocity anisotropy profiles of the simulated halos when only projected phase-space information is available, as in observations. The spherical Jeans equation provides a reliable tool for the determination of cluster mass profiles, also for subsamples of tracers separated by galaxy color. Results are equally good for prolate and oblate clusters. Using only projected phase-space information, MAMPOSSt provides estimates of the mass profile with a standard deviation of 35-69 %, and a negative bias of 7-17 %, nearly independent of radius, and that we attribute to the presence of interlopers in the projected samples. The bias changes sign, that is, the mass is over-estimated, for prolate clusters with their major axis aligned along the line-of-sight. MAMPOSSt measures the velocity anisotropy profiles accurately in the inner cluster regions, with a slight overestimate in the outer regions, both for the whole sample of observationally-identified cluster members and separately for red and blue galaxies.