Local Dark Matter Distribution Research Articles

Slant, Fan, and Narrow: The Response of Stellar Streams to a Tilting Galactic Disk

Stellar streams are sensitive tracers of the gravitational potential, which is typically assumed to be static in the inner Galaxy. However, massive mergers like Gaia-Sausage-Enceladus can impart torques on the stellar disk of the Milky Way that result in the disk tilting at rates of up to 10°–20° Gyr−1. Here, we demonstrate the effects of disk tilting on the morphology and kinematics of stellar streams. Through a series of numerical experiments, we find that streams with nearby apocenters (r apo ≲ 20 kpc) are sensitive to disk tilting, with the primary effect being changes to the stream’s on-sky track and width. Interestingly, disk tilting can produce both more diffuse streams and more narrow streams, depending on the orbital inclination of the progenitor and the direction in which the disk is tilting. Our model of Pal 5's tidal tails for a tilting rate of 15° Gyr−1 is in excellent agreement with the observed stream’s track and width, and reproduces the extreme narrowing of the trailing tail. We also find that failure to account for a tilting disk can bias constraints on shape parameters of the Milky Way’s local dark matter distribution at the level of 5%–10%, with the direction of the bias changing for different streams. Disk tilting could therefore explain discrepancies in the Milky Way’s dark matter halo shape inferred using different streams.

The impact of the Large Magellanic Cloud on dark matter direct detection signals

We study the effect of the Large Magellanic Cloud (LMC) on the dark matter (DM) distribution in the Solar neighborhood, utilizing the Auriga magneto-hydrodynamical simulations of Milky Way (MW) analogues that have an LMC-like system. We extract the local DM velocity distribution at different times during the orbit of the LMC around the MW in the simulations. As found in previous idealized simulations of the MW-LMC system, we find that the DM particles in the Solar neighborhood originating from the LMC analogue dominate the high speed tail of the local DM speed distribution. Furthermore, the native DM particles of the MW in the Solar region are boosted to higher speeds as a result of a response to the LMC's motion.We simulate the signals expected in near future xenon, germanium, and silicon direct detection experiments, considering DM interactions with target nuclei or electrons. We find that the presence of the LMC causes a considerable shift in the expected direct detection exclusion limits towards smaller cross sections and DM masses, with the effect being more prominent for low mass DM. Hence, our study shows, for the first time, that the LMC's influence on the local DM distribution is significant even in fully cosmological MW analogues.

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.

Effects on the local dark matter distribution due to the large magellanic cloud

ABSTRACT We study the local dark matter distribution in two models for the Milky Way (MW)-Large Magellanic Cloud (LMC) interaction. The effect of the LMC on the local dark matter distribution is dependent on the evolution of the MW-LMC system, such that a static model is insufficient to accurately model the dark matter velocity distribution in the solar neighbourhood. An evolved model boosts local LMC dark matter particle velocities by nearly 50 per cent, to a median value of ≈750 km s−1. MW dark matter particles also experience a velocity boost caused by disc reflex motion owing to the infall of the LMC. We study the implications of LMC particles in the solar neighbourhood for dark matter detection experiments. The directionality of LMC particles is distinguishable from the MW particles, with a difference in the apparent origin centroid location between the MW and LMC particles of 26 ± 6°. This unique identifier along with their high velocities can be utilized by directional detectors to search for dark matter particles originating in the LMC.

Discrimination of anisotropy in dark matter velocity distribution with directional detectors

Directional detection of dark matter has sensitivity for both recoil energy and direction of nuclear recoil. It opens the way to measure local velocity distribution of dark matter. In this paper, we study possibility to discriminate isotropic distribution and anisotropic one suggested by a N-body simulation with directional detector. Numerical simulation is performed for two cases according to the detectors, one corresponds to angular histogram and the other is energy-angular distribution of the signals. We reveal that the anisotropy of velocity distribution can be discriminated at 90% C.L. with chi-squared test if O($10^4$) signals are obtained.

Global analyses of Higgs portal singlet dark matter models using GAMBIT

We present global analyses of effective Higgs portal dark matter models in the frequentist and Bayesian statistical frameworks. Complementing earlier studies of the scalar Higgs portal, we use GAMBIT to determine the preferred mass and coupling ranges for models with vector, Majorana and Dirac fermion dark matter. We also assess the relative plausibility of all four models using Bayesian model comparison. Our analysis includes up-to-date likelihood functions for the dark matter relic density, invisible Higgs decays, and direct and indirect searches for weakly-interacting dark matter including the latest XENON1T data. We also account for important uncertainties arising from the local density and velocity distribution of dark matter, nuclear matrix elements relevant to direct detection, and Standard Model masses and couplings. In all Higgs portal models, we find parameter regions that can explain all of dark matter and give a good fit to all data. The case of vector dark matter requires the most tuning and is therefore slightly disfavoured from a Bayesian point of view. In the case of fermionic dark matter, we find a strong preference for including a CP-violating phase that allows suppression of constraints from direct detection experiments, with odds in favour of CP violation of the order of 100:1. Finally, we present DDCalc2.0.0, a tool for calculating direct detection observables and likelihoods for arbitrary non-relativistic effective operators.

Vertical waves in the solar neighbourhood inGaiaDR2

The vertical structure and dynamics of stars in our local Galactic neighbourhood contains much information about the local distribution of visible and dark matter and of perturbations to the Milky Way disc. We use data on the positions and velocities of stars in the solar neighbourhood from \gaia\ DR2 and large spectroscopic surveys to investigate the vertical number counts and mean-velocity trend as a function of distance from the local Galactic mid-plane. We perform a detailed measurement of the wave-like North-South asymmetry in the vertical number counts, which reveals a number of deficits at heights $\approx 0.4\,\mathrm{kpc}$, $\approx 0.9\,\mathrm{kpc}$, and $\approx 1.5\,\mathrm{kpc}$, and peaks at $\approx 0.2\,\mathrm{kpc}$, $\approx 0.7\,\mathrm{kpc}$, and $\approx 1.1\,\mathrm{kpc}$. We find that the asymmetry pattern is independent of colour. The mean vertical velocity is almost constant to $<1\,\mathrm{km\,s}^{-1}$ within a few 100 pc from the mid-plane and then displays a North-South symmetric dip at $\approx0.5\,\mathrm{kpc}$ with an amplitude of $\approx 2\,\mathrm{km\,s}^{-1}$ that is a plausible velocity counterpart to the main number-count dip at a similar height. Thus, with \gaia\ DR2 we confirm at high fidelity that the local Galactic disc is undergoing a wave-like oscillation and a dynamically-consistent observational picture of the perturbed local vertical structure emerges for the first time. We also present the most precise and accurate determination of the Sun's height above the local Galactic mid-plane, correcting for any asymmetry in the vertical density: $z_\odot = 20.8 \pm 0.3\,\mathrm{pc}$.

Galactic dark matter population as the source of neutrino masses

We propose that neutrino masses can be zero $in ~vacuo$ and may be generated by the local distribution of dark matter through a feeble long range scalar force. We discuss potential phenomenological constraints and implications of this framework. Our model typically implies that the relic neutrino background left over from the Big Bang is mostly absent in our Galactic neighborhood. Hence, a positive detection signal from future proposed experiments, such as PTOLEMY, could in principle falsify our scenario.

The metal-poor stellar halo in RAVE-TGAS and its implications for the velocity distribution of dark matter

The local velocity distribution of dark matter plays an integral role in interpreting the results from direct detection experiments. We previously showed that metal-poor halo stars serve as excellent tracers of the virialized dark matter velocity distribution using a high-resolution hydrodynamic simulation of a Milky Way-like halo. In this paper, we take advantage of the first Gaia data release, coupled with spectroscopic measurements from the RAdial Velocity Experiment (RAVE), to study the kinematics of stars belonging to the metal-poor halo within an average distance of ∼5 kpc of the Sun. We study stars with iron abundances [Fe/H] < −1.5 and −1.8 that are located more than 1.5 kpc from the Galactic plane. Using a Gaussian mixture model analysis, we identify the stars that belong to the halo population, as well as some kinematic outliers. We find that both metallicity samples have similar velocity distributions for the halo component, within uncertainties. Assuming that the stellar halo velocities adequately trace the virialized dark matter, we study the implications for direct detection experiments. The Standard Halo Model, which is typically assumed for dark matter, is discrepant with the empirical distribution by ∼6σ, predicts fewer high-speed particles, and is anisotropic. As a result, the Standard Halo Model overpredicts the nuclear scattering rate for dark matter masses below ∼10 GeV. The kinematic outliers that we identify may potentially be correlated with dark matter substructure, though further study is needed to establish this correspondence.

Astrophysical uncertainties on the local dark matter distribution and direct detection experiments

The differential event rate in weakly interacting massive particle (WIMP) direct detection experiments depends on the local dark matter density and velocity distribution. Accurate modelling of the local dark matter distribution is therefore required to obtain reliable constraints on the WIMP particle physics properties. Data analyses typically use a simple standard halo model which might not be a good approximation to the real Milky Way (MW) halo. We review observational determinations of the local dark matter density, circular speed and escape speed and also studies of the local dark matter distribution in simulated MW-like galaxies. We discuss the effects of the uncertainties in these quantities on the energy spectrum and its time and direction dependence. Finally, we conclude with an overview of various methods for handling these astrophysical uncertainties.

Axion astronomy with microwave cavity experiments

Terrestrial searches for the conversion of dark matter axions or axion-like particles into photons inside magnetic fields are sensitive to the phase space structure of the local Milky Way halo. We simulate signals in a hypothetical future experiment based on the Axion Dark Matter eXperiment (ADMX) that could be performed once the axion has been detected and a frequency range containing the axion mass has been identified. We develop a statistical analysis to extract astrophysical parameters, such as the halo velocity dispersion and laboratory velocity, from such data and find that with only a few days integration time a level of precision can be reached matching that of astronomical observations. For longer experiments lasting up to a year in duration we find that exploiting the modulation of the power spectrum in time allows accurate measurements of the Solar peculiar velocity with an accuracy that would improve upon astronomical observations. We also simulate signals based on results from N-body simulations and find that finer substructure in the form of tidal streams would show up prominently in future data, even if only a subdominant contribution to the local dark matter distribution. In these cases it would be possible to reconstruct all the properties of a dark matter stream using the time and frequency dependence of the signal. Finally we consider the detection prospects for a network of streams from tidally disrupted axion miniclusters. These features appear much more prominently in the resolved spectrum than suggested by calculations based on a scan over a range of resonant frequencies, making the detection of axion minicluster streams more viable than previously thought. These results confirm that haloscope experiments in a post-discovery era are able to perform "axion astronomy".

Directional detection of dark matter streams

Directional detection of WIMPs, in which the energies and directions of the recoiling nuclei are measured, currently presents the only prospect for probing the local velocity distribution of Galactic dark matter. We investigate the extent to which future directional detectors would be capable of probing dark matter substructure in the form of streams. We analyse the signal expected from a Sagittarius-like stream and also explore the full parameter space of stream speed, direction, dispersion and density. Using a combination of non-parametric directional statistics, a profile likelihood ratio test and Bayesian parameter inference we find that within acceptable exposure times (O(10) kg yr for cross sections just below the current exclusion limits) future directional detectors will be sensitive to a wide range of stream velocities and densities. We also examine and discuss the importance of the energy window of the detector.

What is the probability that direct detection experiments have observed dark matter?

In Dark Matter direct detection we are facing the situation of some experiments reporting positive signals which are in conflict with limits from other experiments. Such conclusions are subject to large uncertainties introduced by the poorly known local Dark Matter distribution. We present a method to calculate an upper bound on the joint probability of obtaining the outcome of two potentially conflicting experiments under the assumption that the Dark Matter hypothesis is correct, but completely independent of assumptions about the Dark Matter distribution. In this way we can quantify the compatibility of two experiments in an astrophysics independent way. We illustrate our method by testing the compatibility of the hints reported by DAMA and CDMS-Si with the limits from the LUX and SuperCDMS experiments. The method does not require Monte Carlo simulations but is mostly based on using Poisson statistics. In order to deal with signals of few events we introduce the so-called ``signal length'' to take into account energy information. The signal length method provides a simple way to calculate the probability to obtain a given experimental outcome under a specified Dark Matter and background hypothesis.

THE DISTRIBUTION OF DARK MATTER IN THE MILKY WAY’S DISK

We present an analysis of the effects of dissipational baryonic physics on the local dark matter (DM) distribution at the location of the Sun, with an emphasis on the consequences for direct detection experiments. Our work is based on a comparative analysis of two cosmological simulations of a Milky Way halo: Eris, a full hydrodynamic simulation, and ErisDark, its DM-only realization. We find that two distinct processes lead in Eris to a 30% enhancement of DM in the disk plane at the location of the Sun: the accretion and disruption of satellites resulting in a DM component with net angular momentum and the contraction of baryons pulling DM into the disk plane without forcing it to co-rotate. Owing to its particularly quiescent merger history for dark halos of Milky Way mass, the co-rotating dark disk in Eris is less massive than what has been suggested by previous work, contributing only 9% of the local DM density. The speed distribution in Eris is broadened and shifted to higher speeds compared to its DM-only twin simulation ErisDark. At high speeds $f(v)$ falls more steeply in Eris than in ErisDark or the SHM, easing the tension between recent results from the CDMS-II and XENON100 experiments. The non-Maxwellian aspects of $f(v)$ are still present, but much less pronounced in Eris than in DM-only runs. The weak dark disk increases the time-averaged scattering rate by only a few percent at low recoil energies. On the high velocity tail, however, the increase in typical speeds due to baryonic contraction results in strongly enhanced mean scattering rates compared to ErisDark, although they are still suppressed compared to the SHM. Similar trends are seen regarding the amplitude of the annual modulation, while the modulated fraction is increased compared to the SHM and decreased compared to ErisDark.

Anisotropic dark matter distribution functions and impact on WIMP direct detection

Dark matter N-body simulations suggest that the velocity distribution of dark matter is anisotropic. In this work we employ a mass model for the Milky Way whose parameters are determined from a fit to kinematical data. Then we adopt an ansatz for the dark matter phase space distribution which allows to construct self-consistent halo models which feature a degree of anisotropy as a function of the radius such as suggested by the simulations. The resulting velocity distributions are then used for an analysis of current data from dark matter direct detection experiments. We find that velocity distributions which are radially biased at large galactocentric distances (up to the virial radius) lead to an increased high velocity tail of the local dark matter distribution. This affects the interpretation of data from direct detection experiments, especially for dark matter masses around 10 GeV, since in this region the high velocity tail is sampled. We find that the allowed regions in the dark matter mass-cross section plane as indicated by possible hints for a dark matter signal reported by several experiments as well as conflicting exclusion limits from other experiments shift in a similar way when the halo model is varied. Hence, it is not possible to improve the consistency of the data by referring to anisotropic halo models of the type considered in this work.

PROBING THE HALO FROM THE SOLAR VICINITY TO THE OUTER GALAXY: CONNECTING STARS IN LOCAL VELOCITY STRUCTURES TO LARGE-SCALE CLOUDS

(Abridged) This paper presents the first connections made between two local features in velocity-space found in a survey of M giant stars and stellar spatial inhomogeneities on global scales. Comparison to cosmological, chemodynamical stellar halo models confirm that the M giant population is particularly sensitive to rare, recent and massive accretion events. These events can give rise to local observed velocity sequences - a signature of a small fraction of debris from a common progenitor, passing at high velocity through the survey volume, near the pericenters of their eccentric orbits. The majority of the debris is found in much larger structures, whose morphologies are more cloud-like than stream-like and which lie at the orbital apocenters. Adopting this interpretation, the full-space motions represented by the observed velocity features are derived under the assumption that the members within each sequence share a common velocity. Orbit integrations are then used to trace the past and future trajectories of these stars across the sky revealing plausible associations with large, previously-discovered, cloud-like structures. The connections made between nearby velocity structures and these distant clouds represent preliminary steps towards developing coherent maps of such giant debris systems. These maps promise to provide new insights into the origin of debris clouds, new probes of Galactic history and structure, and new constraints on the high-velocity tails of the local dark matter distribution that are essential for interpreting direct detection experiments.

Analyzing direct dark matter detection data with unrejected background events by the AMIDAS website

In this talk I have presented the data analysis results of extracting properties of halo WIMPs: the mass and the (ratios between the) spin-independent and spin-dependent couplings/cross sections on nucleons by the AMIDAS website by taking into account possible unrejected background events in the analyzed data sets. Although non-standard astronomical setup has been used to generate pseudodata sets for our analyses, it has been found that, without prior information/assumption about the local density and velocity distribution of halo Dark Matter, these WIMP properties have been reconstructed with ∼ 2% to ≲ 30% deviations from the input values.

Probing the local velocity distribution of WIMP dark matter with directional detectors

We explore the ability of directional nuclear-recoil detectors to constrain the local velocity distribution of weakly interacting massive particle (WIMP) dark matter by performing Bayesian parameter estimation on simulated recoil-event data sets. We discuss in detail how directional information, when combined with measurements of the recoil-energy spectrum, helps break degeneracies in the velocity-distribution parameters. We also consider the possibility that velocity structures such as cold tidal streams or a dark disk may also be present in addition to the Galactic halo. Assuming a CF4 detector with a 30-kg-yr exposure, a 50-GeV WIMP mass, and a WIMP-nucleon spin-dependent cross-section of 10−3pb, we show that the properties of a cold tidal stream may be well constrained. However, measurement of the parameters of a dark-disk component with a low lag speed of ∼ 50km/s may be challenging unless energy thresholds are improved.

Toward a consistent picture for CRESST, CoGeNT, and DAMA

Three dark matter direct detection experiments (DAMA/LIBRA, CoGeNT, and CRESST-II) have each reported signals which are not consistent with known backgrounds, but resemble that predicted for a dark matter particle with a mass of roughly $\sim$10 GeV and an elastic scattering cross section with nucleons of $\sim$$10^{-41}$--$10^{-40}$ cm$^2$. In this article, we compare the signals of these experiments and discuss whether they can be explained by a single species of dark matter particle, without conflicting with the constraints of other experiments. We find that the spectrum of events reported by CoGeNT and CRESST-II are consistent with each other and with the constraints from CDMS-II, although some tension with xenon-based experiments remains. Similarly, the modulation signals reported by DAMA/LIBRA and CoGeNT appear to be compatible, although the corresponding amplitude of the observed modulations are a factor of at least a few higher than would be naively expected, based on the event spectra reported by CoGeNT and CRESST-II. This apparent discrepancy could potentially be resolved if tidal streams or other non-Maxwellian structures are present in the local distribution of dark matter.

Interpreting dark matter direct detection independently of the local velocity and density distribution

We demonstrate precisely what particle physics information can be extracted from a single direct detection observation of dark matter while making absolutely no assumptions about the local velocity distribution and local density of dark matter. Our central conclusions follow from a very simple observation: the velocity distribution of dark matter is positive definite, $f(v)\ensuremath{\ge}0$. We demonstrate the utility of this result in several ways. First, we show a falling deconvoluted recoil spectrum (deconvoluted of the nuclear form factor), such as from ordinary elastic scattering, can be ``mocked up'' by any mass of dark matter above a kinematic minimum. As an example, we show that dark matter much heavier than previously considered can explain the CoGeNT excess. Specifically, ${m}_{\ensuremath{\chi}}&lt;{m}_{\mathrm{Ge}}$ can be in just as good agreement as light dark matter, while ${m}_{\ensuremath{\chi}}&gt;{m}_{\mathrm{Ge}}$ depends on understanding the sensitivity of xenon to dark matter at very low recoil energies, ${E}_{R}\ensuremath{\lesssim}6\text{ }\text{ }\mathrm{keVnr}$. Second, we show that any rise in the deconvoluted recoil spectrum represents distinct particle physics information that cannot be faked by an arbitrary $f(v)$. As examples of resulting nontrivial particle physics, we show that inelastic dark matter and dark matter with a form factor can both yield such a rise.