Local Dark Matter Velocity Distribution Research Articles

Are there any extragalactic high speed dark matter particles in the Solar neighborhood?

We use the APOSTLE suite of cosmological hydrodynamical simulations of the Local Group to examine the high speed tail of the local dark matter velocity distribution in simulated Milky Way analogues. The velocity distribution in the Solar neighborhood is well approximated by a generalized Maxwellian distribution sharply truncated at a well-defined maximum “escape” speed. The truncated generalized Maxwellian distribution accurately models the local dark matter velocity distribution of all our Milky Way analogues, with no evidence for any separate extragalactic high-speed components. The local maximum speed is well approximated by the terminal velocity expected for particles able to reach the Solar neighborhood in a Hubble time from the farthest confines of the Local Group. This timing constraint means that the local dark matter velocity distribution is unlikely to contain any high-speed particles contributed by the Virgo Supercluster “envelope”, as argued in recent work. Particles in the Solar neighborhood with speeds close to the local maximum speed can reach well outside the virial radius of the Galaxy, and, in that sense, belong to the Local Group envelope posited in earlier work. The local manifestation of such envelope is thus not a distinct high-speed component, but rather simply the high-speed tail of the truncated Maxwellian distribution.

The local dark matter distribution in self-interacting dark matter halos

We study the effects of dark matter self-interactions on the local dark matter distribution in selected Milky Way-like galaxies in the EAGLE hydrodynamical simulations. The simulations were run with two different self-interacting dark matter models, a constant and velocity-dependent self-interaction cross-section. We find that the local dark matter velocity distribution of the Milky Way-like halos in the simulations with dark matter self-interactions and baryons are generally similar to those extracted from cold collisionless dark matter simulations with baryons. In both cases, the local dark matter speed distributions agree well with their best fit Maxwellian distributions. Including baryons in the simulations with or without dark matter self-interactions increases the local dark matter density and shifts the dark matter speed distributions to higher speeds. To study the implications for direct detection, we compute the dark matter halo integrals obtained directly from the simulations and compare them to those obtained from the best fit Maxwellian velocity distribution. We find that a Maxwellian distribution provides a good fit to the halo integrals of most halos, without any significant difference between the results of different dark matter self-interaction models.

Equilibrium axisymmetric halo model for the Milky Way and its implications for direct and indirect dark matter searches

We for the first time provide self-consistent axisymmetric phase-space distribution models for the Milky Way's dark matter (DM) halo which are carefully matched against the latest kinematic measurements through Bayesian analysis. By using broad priors on the individual galactic components, we derive conservative estimates for the astrophysical factors entering the interpretation of direct and indirect DM searches. While the resulting DM density profiles are in good agreement with previous studies, implying $\rho_\odot \approx 10^{-2} \, M_\odot / \mathrm{pc}^3$, the presence of baryonic disc leads to significant differences in the local DM velocity distribution in comparison with the standard halo model. For direct detection, this implies roughly 30% stronger cross-section limits at DM masses near detectors maximum sensitivity and up to an order of magnitude weaker limits at the lower end of the mass range. Furthermore, by performing Monte-Carlo simulations for the upcoming DARWIN and DarkSide-20k experiments, we demonstrate that upon successful detection of heavy DM with coupling just below the current limits, the carefully constructed axisymmetric models can eliminate bias and reduce uncertainties by more then 50% in the reconstructed DM coupling and mass, but also help in a more reliable determination of the scattering operator. Furthermore, the velocity anisotropies induced by the baryonic disc can lead to significantly larger annual modulation amplitude and sizable differences in the directional distribution of the expected DM-induced events. For indirect searches, we provide the differential $J$-factors and compute several moments of the relative velocity distribution that are needed for predicting the rate of velocity-dependent annihilations. However, we find that accurate predictions are still hindered by large uncertainties regarding the DM distribution near the galactic center.

Under the FIRElight: Stellar Tracers of the Local Dark Matter Velocity Distribution in the Milky Way

The Gaia era opens new possibilities for discovering the remnants of disrupted satellite galaxies in the Solar neighborhood. If the population of local accreted stars is correlated with the dark matter sourced by the same mergers, one can then map the dark matter distribution directly. Using two cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies from the Latte suite of Fire-2 simulations, we find a strong correlation between the velocity distribution of stars and dark matter at the solar circle that were accreted from luminous satellites. This correspondence holds for dark matter that is either relaxed or in kinematic substructure called debris flow, and is consistent between two simulated hosts with different merger histories. The correspondence is more problematic for streams because of possible spatial offsets between the dark matter and stars. We demonstrate how to reconstruct the dark matter velocity distribution from the observed properties of the accreted stellar population by properly accounting for the ratio of stars to dark matter contributed by individual mergers. After demonstrating this method using the Fire-2 simulations, we apply it to the Milky Way and use it to recover the dark matter velocity distribution associated with the recently discovered stellar debris field in the Solar neighborhood. Based on results from Gaia, we estimate that $42 ^{+26}_{-22}\%$ of the local dark matter that is accreted from luminous mergers is in debris flow.

On the correlation between the local dark matter and stellar velocities

The dark matter velocity distribution in the Solar neighbourhood is an important astrophysical input which enters in the predicted event rate of dark matter direct detection experiments. It has been recently suggested that the local dark matter velocity distribution can be inferred from that of old or metal-poor stars in the Milky Way. We investigate this potential relation using six high resolution magneto-hydrodynamical simulations of Milky Way-like galaxies of the Auriga project. We do not find any correlation between the velocity distributions of dark matter and old stars in the Solar neighbourhood. Likewise, there are no strong correlations between the local velocity distributions of dark matter and metal-poor stars selected by applying reasonable cuts on metallicity. In some simulated galaxies, extremely metal-poor stars have a velocity distribution that is statistically consistent with that of the dark matter, but the sample of such stars is so small that we cannot draw any strong conclusions.

Directional axion detection

We develop a formalism to describe extensions of existing axion haloscope designs to those that possess directional sensitivity to incoming dark matter axion velocities. The effects are measurable if experiments are designed to have dimensions that approach the typical coherence length for the local axion field. With directional sensitivity, axion detection experiments would have a greatly enhanced potential to probe the local dark matter velocity distribution. We develop our formalism generally, but apply it to specific experimental designs, namely resonant cavities and dielectric disk haloscopes. We demonstrate that these experiments are capable of measuring the daily modulation of the dark matter signal and using it to reconstruct the three-dimensional velocity distribution. This allows one to measure the Solar peculiar velocity, probe the anisotropy of the dark matter velocity ellipsoid and identify cold substructures such as the recently discovered streams near to Earth. Directional experiments can also identify features over much shorter timescales, potentially facilitating the mapping of debris from axion miniclusters.

Optimized velocity distributions for direct dark matter detection

We present a method to calculate, without making assumptions about the local dark matter velocity distribution, the maximal and minimal number of signal events in a direct detection experiment given a set of constraints from other direct detection experiments and/or neutrino telescopes. The method also allows to determine the velocity distribution that optimizes the signal rates. We illustrate our method with three concrete applications: i) to derive a halo-independent upper limit on the cross section from a set of null results, ii) to confront in a halo-independent way a detection claim to a set of null results and iii) to assess, in a halo-independent manner, the prospects for detection in a future experiment given a set of current null results.

Reconstructing the three-dimensional local dark matter velocity distribution

Directionally sensitive dark matter (DM) direct detection experiments present the only way to observe the full three-dimensional velocity distribution of the Milky Way halo local to Earth. In this work we compare methods for extracting information about the local DM velocity distribution from a set of recoil directions and energies in a range of hypothetical directional and non-directional experiments. We compare a model independent empirical parameterisation of the velocity distribution based on an angular discretisation with a model dependent approach which assumes knowledge of the functional form of the distribution. The methods are tested under three distinct halo models which cover a range of possible phase space structures for the local velocity distribution: a smooth Maxwellian halo, a tidal stream and a debris flow. In each case we use simulated directional data to attempt to reconstruct the shape and parameters describing each model as well as the DM particle properties. We find that the empirical parametrisation is able to make accurate unbiased reconstructions of the DM mass and cross section as well as capture features in the underlying velocity distribution in certain directions without any assumptions about its true functional form. We also find that by extracting directionally averaged velocity parameters with this method one can discriminate between halo models with different classes of substructure.

ASSESSING ASTROPHYSICAL UNCERTAINTIES IN DIRECT DETECTION WITH GALAXY SIMULATIONS

ABSTRACT We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence of a baryonic disk increases the amount of high velocity dark matter, resulting in velocity distributions that are more similar to the Maxwellian Standard Halo Model than predicted from dark matter-only simulations. Furthermore, the velocity structures present in baryonic simulations possess a greater diversity than expected from dark matter-only simulations. We show that the impact on the direct detection experiments LUX, DAMA/Libra, and CoGeNT using our simulated velocity distributions, and explore how resolution and halo mass within the Milky Way’s estimated mass range impact the results. A Maxwellian fit to the velocity distribution tends to overpredict the amount of dark matter in the high velocity tail, even with baryons, and thus leads to overly optimistic direct detection bounds on models that are dependent on this region of phase space for an experimental signal. Our work further demonstrates that it is critical to transform simulated velocity distributions to the lab frame of reference, due to the fact that velocity structure in the solar neighborhood appears when baryons are included. There is more velocity structure present when baryons are included than in dark matter-only simulations. Even when baryons are included, the importance of the velocity structure is not as apparent in the Galactic frame of reference as in the Earth frame.

A new halo-independent approach to dark matter direct detection analysis

Uncertainty in the local dark matter velocity distribution is a key difficulty in the analysis of data from direct detection experiments. Here we propose a new approach for dealing with this uncertainty, which does not involve any assumptions about the structure of the dark matter halo. Given a dark matter model, our method yields the velocity distribution which best describes a set of direct detection data as a finite sum of streams with optimised speeds and densities. The method is conceptually simple and numerically very efficient. We give an explicit example in which the method is applied to determining the ratio of proton to neutron couplings of dark matter from a hypothetical set of future data.

The impact of the phase-space density on the indirect detection of dark matter

We study the indirect detection of dark matter when the local dark matter velocity distribution depends upon position, as expected for the Milky Way and its dwarf spheroidal satellites, and the annihilation cross-section is not purely s-wave. Using a phase-space distribution consistent with the dark matter density profile, we present estimates of cosmic and gamma-ray fluxes from dark matter annihilations. The expectations for the indirect detection of dark matter can differ significantly from the usual calculation that assumes that the velocity of the dark matter particles follows a Maxwell-Boltzmann distribution.

Inelastic dark matter, non-standard halos and the DAMA/LIBRA results

The DAMA collaboration have claimed to detect particle dark matter (DM) via an annual modulation in their observed recoil event rate. This appears to be in strong disagreement with the null results of other experiments if interpreted in terms of elastic DM scattering, while agreement for a small region of parameter space is possible for inelastic DM (iDM) due to the altered kinematics of the collision. To date most analyses assume a simple galactic halo DM velocity distribution, the Standard Halo Model, but direct experimental support for the SHM is severely lacking and theoretical studies indicate possible significant differences. We investigate the dependence of DAMA and the other direct detection experiments on the local DM velocity distribution, utilizing the results of the Via Lactea and Dark Disc numerical simulations. We also investigate effects of varying the solar circular velocity, the DM escape velocity, and the DAMA quenching factor within experimental limits. Our data set includes the latest ZEPLIN-III results, as well as full publicly available data sets. Due to the more sensitive dependence of the inelastic cross section on the velocity distribution, we find that with Via Lactea the DAMA results can be consistent with all other experiments over an enlarged region of iDM parameter space, with higher mass particles being preferred, while Dark Disc does not lead to an improvement. A definitive test of DAMA for iDM requires heavy element detectors.

Velocity distributions and annual-modulation signatures of weakly-interacting massive particles

An annual modulation in the event rate of the NaI detector of the DAMA collaboration has been used to infer the existence of particle dark matter in the Galactic halo. Bounds on the WIMP mass and WIMP-nucleon cross section have been derived. These analyses have assumed that the local dark-matter velocity distribution is either isotropic or has some bulk rotation. Here we consider the effects of possible structure in the WIMP velocity distribution on the annual-modulation amplitude. We show that if we allow for a locally anisotropic velocity dispersion tensor, the interpretation of direct detection experiments could be altered significantly. We also show that uncertainties in the velocity distribution that arise from uncertainties in the radial density profile are less important if the velocity dispersion is assumed to be isotropic.