Caustics, Cold Flows, and Annual Modulation

Earth's rotation about the Sun produces an annual modulation in the expected scattering rate at direct dark matter detection experiments. The annual modulation as a function of the recoil energy ER imparted by the dark matter particle to a target nucleus is expected to vary depending on the detector material. However, for most interactions a change of variables from ER to vmin, the minimum speed a dark matter particle must have to impart a fixed ER to a target nucleus, produces an annual modulation independent of the target element. We recently showed that if the dark matter-nucleus cross section contains a non-factorizable target and dark matter velocity dependence, the annual modulation as a function of vmin can be target dependent. Here we examine more extensively the necessary conditions for target-dependent modulation, its observability in present-day experiments, and the extent to which putative signals could identify a dark matter-nucleus differential cross section with a non-factorizable dependence on the dark matter velocity.

Light yield determination in large sodium iodide detectors applied in the search for dark matter

We present the first and so far the only simulations to follow the fine-grained phase-space structure of galaxy haloes formed from generic LCDM initial conditions. We integrate the geodesic deviation equation in tandem with the N-body equations of motion, demonstrating that this can produce numerically converged results for the properties of fine-grained phase-space streams and their associated caustics, even in the inner regions of haloes. Our effective resolution for such structures is many orders of magnitude better than achieved by conventional techniques on even the largest simulations. We apply these methods to the six Milky Way-mass haloes of the Aquarius Project. At 8 kpc from halo centre a typical point intersects about 10^14 streams with a very broad range of individual densities; the ~10^6 most massive streams contribute about half of the local dark matter density. As a result, the velocity distribution of dark matter particles should be very smooth with the most massive fine-grained stream contributing about 0.1% of the total signal. Dark matter particles at this radius have typically passed 200 caustics since the Big Bang. The peak densities on present-day caustics in the inner halo almost all lie well below the mean local dark matter density. As a result caustics provide a negligible boost (<0.1%) to the predicted local dark matter annihilation rate. The effective boost is larger in the outer halo but never exceeds about 10%. Thus fine-grained streams and their associated caustics have no effect on the detectability of dark matter, either directly in Earth-bound laboratories, or indirectly through annihilation radiation, with the exception that resonant cavity experiments searching for axions may see the most massive local fine-grained streams because of their extreme localisation in energy/momentum space. (abridged)

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.

Observations of galaxies and primordial radiation suggest that the Universe is made mostly of non-luminous dark matter1,2. Several new types of fundamental particle have been proposed as candidates for dark matter3, such as weakly interacting massive particles (WIMPs)4,5. These particles would be expected to interact with nuclei in suitable detector materials on Earth, for example, causing them to recoil. However, no definitive signal from such dark-matter interactions has been detected despite concerted efforts by many collaborations6. One exception is the much-debated claim by the DAMA collaboration of a statistically significant (more than nine standard deviations) annual modulation in the rate of nuclear interaction events. Annual modulation is expected because of the variation in Earth's velocity relative to the Galaxy's dark-matter halo that arises from Earth's orbital motion around the Sun. DAMA observed a modulation in the rate of interaction events in their detector7-9 with a period and phase consistent with that expected for WIMPs10-12. Several groups have been working to develop experiments with the aim of reproducing DAMA's results using the same target medium (sodium iodide)13-17. To determine whether there is evidence for an excess of events above the expected background in sodium iodide and to look for evidence of an annual modulation, the COSINE-100 experiment uses sodium iodide as the target medium to carry out a model-independent test of DAMA's claim. Here we report results from the initial operation of the COSINE-100 experiment related to the first task18,19. We observe no excess of signal-like events above the expected background in the first 59.5 days of data from COSINE-100. Assuming the so-called standard dark-matter halo model, this result rules out spin-independent [corrected] WIMP-nucleon interactions as the cause of the annual modulation observed by the DAMA collaboration.20-23. The exclusion limit on the WIMP-sodium interaction cross-section is 1.14 × 10-40 cm2 for 10-GeV c-2 WIMPs at a 90% confidence level. The COSINE-100 experiment will continue to collect data for two more years, enabling a model-independent test of the annual modulation observed by the DAMA collaboration.

The next generation of axion direct-detection experiments may rule out or confirm axions as the dominant source of dark matter. We develop a general likelihood-based framework for studying the time-series data at such experiments, with a focus on the role of dark matter astrophysics, to search for signatures of the QCD axion or axionlike particles. We illustrate how in the event of a detection the likelihood framework may be used to extract measures of the local dark matter phase-space distribution, accounting for effects such as annual modulation and gravitational focusing, which is the perturbation to the dark matter phase-space distribution by the gravitational field of the Sun. Moreover, we show how potential dark matter substructure, such as cold dark matter streams or a thick dark disk, could impact the signal. For example, we find that when the bulk dark matter halo is detected at $5\ensuremath{\sigma}$ global significance, the unique time-dependent features imprinted by the dark matter component of the Sagittarius stream, even if only a few percent of the local dark matter density, may be detectable at $\ensuremath{\sim}2\ensuremath{\sigma}$ significance. A corotating dark disk, with lag speed $\ensuremath{\sim}50\text{ }\text{ }\mathrm{km}/\mathrm{s}$, that is $\ensuremath{\sim}20%$ of the local dark matter density could dominate the signal, while colder but as-of-yet unknown substructure may be even more important. Our likelihood formalism, and the results derived with it, are generally applicable to any time-series-based approach to axion direct detection.

Laser shadowgraphy measurements of abrasive particle spatial, size and velocity distributions through micro-masks used in abrasive jet micro-machining

Aromatic organic compounds, because of their small excitation energies $\ensuremath{\sim}\mathcal{O}(\mathrm{few}\text{ }\mathrm{eV})$ and scintillating properties, are promising targets for detecting dark matter of mass $\ensuremath{\sim}\mathcal{O}(\mathrm{few}\text{ }\mathrm{MeV})$. Additionally, their planar molecular structures lead to large anisotropies in the electronic wave functions, yielding a significant daily modulation in the event rate expected to be observed in crystals of these molecules. We characterize the daily modulation rate of dark matter interacting with an anisotropic scintillating organic crystal such as trans-stilbene, and show that daily modulation is an $\ensuremath{\sim}\mathcal{O}(1)$ fraction of the total rate for small DM masses and comparable to, or larger than, the $\ensuremath{\sim}10%$ annual modulation fraction at large DM masses. As we discuss in detail, this modulation provides significant leverage for detecting or excluding dark matter scattering, even in the presence of a non-negligible background rate. Assuming a nonmodulating background rate of $1/\mathrm{min}/\mathrm{kg}$ that scales with total exposure, we find that a $100\text{ }\text{ }\mathrm{kg}\ifmmode\cdot\else\textperiodcentered\fi{}\mathrm{yr}$ experiment is sensitive to the cross section corresponding to the correct relic density for dark matter masses between 1.3--14 MeV (1.5--1000 MeV) if dark matter interacts via a heavy (light) mediator. This modulation can be understood using an effective velocity scale ${v}^{*}=\mathrm{\ensuremath{\Delta}}E/{q}^{*}$, where $\mathrm{\ensuremath{\Delta}}E$ is the electronic transition energy and ${q}^{*}$ is a characteristic momentum scale of the electronic orbitals. We also characterize promising future directions for the development of scintillating organic crystals as dark matter detectors.

The DAMA collaboration reported an annually modulated rate with a phase compatible with a Dark Matter induced signal. We point out that a slowly varying rate can bias or even simulate an annual modulation if data are analyzed in terms of residuals computed by subtracting approximately yearly averages starting from a fixed date, rather than a background continuous in time. In the most extreme case, the amplitude and phase of the annual modulation reported by DAMA could be alternatively interpreted as a decennial growth of the rate. This possibility appears mildly disfavoured by a detailed study of the available data, but cannot be safely excluded. In general, a decreasing or increasing rate could partially reduce or enhance a true annual modulation, respectively. The issue could be clarified by looking at the full time-dependence of the DAMA total rate, not explicitly published so far.

Both experiment and theory recently showed that the H + D2(v = 0, j = 0) → HD(v' = 4, j') + D reactions at a collision energy of 1.97 eV display a seemingly anomalous HD product angular distribution that moves in the backward direction as the value of j' increases and the corresponding energy available for product recoil decreases. This behavior was attributed to the presence of a centrifugal barrier along the reaction path. Here, we show, using fully quantum mechanical calculations, that for low recoil energies, the collision mechanism is nearly independent of the HD internal state and the HD product becomes aligned, with its rotational angular momentum j' pointing perpendicular to the recoil momentum k'. As the kinetic energy to overcome this barrier becomes limited, the three atoms adopt a nearly collinear configuration in the transition-state region to permit reaction, which strongly polarizes the resulting HD product. These results are expected to be general for any chemical reaction in the low recoil energy limit.

The lack of new physics at the LHC so far weakens the argument for TeV scale thermal dark matter. On the other hand, heavier, nonthermal dark matter is generally difficult to test experimentally. Here we consider the interesting and generic case of hypercharged dark matter, which can allow for heavy dark matter masses without spoiling testability. Planned direct detection experiments will be able to see a signal for masses up to an incredible 1010 GeV, and this can further serve to probe the reheating temperature up to about 109 GeV, as determined by the nonthermal dark matter relic abundance. The Z-mediated nature of the dark matter scattering may be determined in principle by comparing scattering rates on different detector nuclei, which in turn can reveal the dark matter mass. We will discuss the extent to which future experiments may be able to make such a determination.

We propose a novel method to significantly enhance the signal rate in qubit-based dark matter detection experiments with the help of quantum interference. Various quantum sensors possess ideal properties for detecting wavelike dark matter, and qubits, commonly employed in quantum computers, are excellent candidates for dark matter detectors. We demonstrate that, by designing an appropriate quantum circuit to manipulate the qubits, the signal rate scales proportionally to n_{q}^{2}, with n_{q} being the number of sensor qubits, rather than linearly with n_{q}. Consequently, in the dark matter detection with a substantial number of sensor qubits, a significant increase in the signal rate can be expected. We provide a specific example of a quantum circuit that achieves this enhancement by coherently combining the phase evolution in each individual qubit due to its interaction with dark matter. We also demonstrate that the circuit is fault tolerant to dephasing noises, a critical quantum noise source in quantum computers. The enhancement mechanism proposed here is applicable to various modalities for quantum computers, provided that the quantum operations relevant to enhancing the dark matter signal can be applied to these devices.

Measurements of solid concentration and particle velocity distributions near the wall of a cyclone

The substructures of the Galactic dark matter halo such as dark matter subhalos and dwarf galaxies have very low velocity dispersions, which makes them useful in constraining the scenario of Sommerfeld-enhanced dark matter annihilation. We calculate the velocity distribution of dark matter particles in the dark matter halo substructures using the Eddington's formula with NFW density profile. We parameterize the effect of Sommerfeld enhancement of s-wave dark matter annihilation on the gamma-ray flux as the Sommerfeld-enhanced J-factors, and explicitly calculate their values for 15 known dwarf spheroidal galaxies. Using the results from the analysis of the unassociated point-sources in the Fermi-LAT 3FGL and Fermi-LAT gamma-ray data observation on the 15 dwarf galaxies, we derive upper limits on the dark matter annihilation cross sections with Sommerfeld enhancement. We find that in a wide region of parameter space, the constraints can be a few orders of magnitude more stringent than that in the case without the Sommerfeld enhancement. For dark matter annihilation dominantly into e+e−, μ+μ−, τ+τ− and bb̄, those constraints can exclude thermal relic dark matter for the dark matter mass below about 1 TeV.

In regions of very high dark matter density such as the Galactic centre, the capture and annihilation of WIMP dark matter by stars has the potential to significantly alter their evolution. We describe the dark stellar evolution code DarkStars, and present a series of grids of WIMP-influenced stellar models for main sequence stars. We describe changes in which occur as a function of the rate of energy injection by WIMPs, for stars of 0.3-2.0 solar masses and metallicities Z = 0.0003-0.02. We show what rates of energy injection can be obtained using realistic orbital parameters for stars at the Galactic centre, including detailed consideration of the velocity and density profiles of dark matter. Capture and annihilation rates are strongly boosted when stars follow elliptical rather than circular orbits. If there is a spike of dark matter induced by the supermassive black hole at the Galactic centre, single solar-mass stars following orbits with periods as long as 50 years and eccentricities as low as 0.9 could be significantly affected. Binary systems with similar periods about the Galactic centre could be affected on even less eccentric orbits. The most striking evidence of this scenario would be the existence of a binary consisting of a low-mass protostar and a higher-mass evolved star. The observation of low-mass stars and/or binaries on such orbits would either provide a detection of WIMP dark matter, or place stringent limits on the combination of the WIMP mass, spin-dependent nuclear-scattering cross-section, halo density and velocity distribution near the Galactic centre. In some cases, the limits on the WIMP mass and spin-dependent nuclear-scattering cross-section would be of comparable sensitivity to current direct-detection experiments.