The dark matter explanation of the 3.5 keV line is strongly disfavored by our work in Dessert et al. 2020. Boyarsky et al. 2020 questions that conclusion: modeling additional background lines is claimed to weaken the limit sufficiently to re-allow a dark matter interpretation. We respond as follows. 1) A more conservative limit is obtained by modeling additional lines; this point appeared in its entirety in our work in Dessert et al., though we also showed that the inclusion of such lines is not necessary. 2) Despite suggestions in Boyarsky et al., even the more conservative limits strongly disfavor a decaying dark matter origin of the 3.5 keV line.

Observations of nearby galaxies and galaxy clusters have reported an unexpected X-ray emission line around 3.5 kilo-electron volts (keV). Proposals to explain this line include decaying dark matter$-$in particular, that the decay of sterile neutrinos with a mass around 7 keV could match the available data. If this interpretation is correct, the 3.5 keV line should also be emitted by dark matter in the halo of the Milky Way. We used more than 30 megaseconds of XMM-Newton (X-ray Multi-Mirror Mission) blank-sky observations to test this hypothesis, finding no evidence of the 3.5-keV line emission from the Milky Way halo. We set an upper limit on the decay rate of dark matter in this mass range, which is inconsistent with the possibility that the 3.5-keV line originates from dark matter decay.

Axion dark matter (DM) may convert to radio-frequency electromagnetic radiation in the strong magnetic fields around neutron stars. The radio signature of such a process would be an ultra-narrow spectral peak at a frequency determined by the mass of the axion particle. We analyze data we collected from the Robert C. Byrd Green Bank Telescope in the L-band and the Effelsberg 100-m Telescope in the L-Band and S-band from a number of sources expected to produce bright signals of axion-photon conversion, including the Galactic Center of the Milky Way and the nearby isolated neutron stars RX J0720.4-3125 and RX J0806.4-4123. We find no evidence for axion DM and are able to set some of the strongest constraints to-date on the existence of axion DM in the highly-motivated mass range between ~5-11 $\mu$eV.

The Galactic Center Excess (GCE) of GeV gamma rays can be explained as a signal of annihilating dark matter or of emission from unresolved astrophysical sources, such as millisecond pulsars. Evidence for the latter is provided by a statistical procedure---referred to as Non-Poissonian Template Fitting (NPTF)---that distinguishes the smooth distribution of photons expected for dark matter annihilation from a "clumpy" photon distribution expected for point sources. In this paper, we perform an extensive study of the NPTF on simulated data, exploring its ability to recover the flux and luminosity function of unresolved sources at the Galactic Center. When astrophysical background emission is perfectly modeled, we find that the NPTF successfully distinguishes between the dark matter and point source hypotheses when either component makes up the entirety of the GCE. When the GCE is a mixture of dark matter and point sources, the NPTF may fail to reconstruct the correct contribution of each component. We further study the impact of mismodeling the Galactic diffuse backgrounds, finding that while a dark matter signal could be attributed to point sources in some outlying cases for the scenarios we consider, the significance of a true point source signal remains robust. Our work enables us to comment on a recent study by Leane and Slatyer (2019) that questions prior NPTF conclusions because the method does not recover an artificial dark matter signal injected on actual Fermi data. We demonstrate that the failure of the NPTF to extract an artificial dark matter signal can be natural when point sources are present in the data---with the effect further exacerbated by the presence of diffuse mismodeling---and does not on its own invalidate the conclusions of the NPTF analysis in the Inner Galaxy.

The Fermi Large Area Telescope has observed an excess of ~GeV energy gamma rays from the center of the Milky Way, which may arise from near-thermal dark matter annihilation. Firmly establishing the dark matter origin for this excess is however complicated by challenges in modeling diffuse cosmic-ray foregrounds as well as unresolved astrophysical sources, such as millisecond pulsars. Non-Poissonian Template Fitting (NPTF) is one statistical technique that has previously been used to show that at least some fraction of the GeV excess is likely due to a population of dim point sources. These results were recently called into question by Leane and Slatyer (2019), who showed that a synthetic dark matter annihilation signal injected on top of the real Fermi data is not recovered by the NPTF procedure. In this work, we perform a dedicated study of the Fermi data and explicitly show that the central result of Leane and Slatyer (2019) is likely driven by the fact that their choice of model for the Galactic foreground emission does not provide a sufficiently good description of the data. We repeat the NPTF analyses using a state-of-the-art model for diffuse gamma-ray emission in the Milky Way and introduce a novel statistical procedure, based on spherical-harmonic marginalization, to provide an improved description of the Galactic diffuse emission in a data-driven fashion. With these improvements, we find that the NPTF results continue to robustly favor the interpretation that the Galactic Center excess is due, in part, to unresolved astrophysical point sources across the analysis variations that we have explored.

Axions may be produced thermally inside the cores of neutron stars (NSs), escape the stars due to their weak interactions with matter, and subsequently convert into X-rays in the magnetic fields surrounding the NSs. We perform a hard X-ray search from 2 - 8 keV arising from this emission mechanism from the nearby Magnificent Seven isolated NSs using archival XMM-Newton and Chandra data. These NSs are unique in that they have only been detected in the UV and soft X-ray bands; their observed fluxes are thought to arise from thermal emission at temperatures ~100 eV from the NS surfaces. We use the X-ray datasets to set the strongest limits to-date on the product of the axion-photon and axion-nucleon couplings for axion masses below ~$10^{-4}$ eV. Moreover, we interpret an observed excess of hard X-rays from the Magnificent Seven in the context of the axion model, and we discuss possible future measurements to help rule out or confirm this possibility.

We report significant hard X-ray excesses in the energy range 2-8 keV for two nearby isolated neutron stars RX J1856.6-3754 and RX J0420.0-5022. These neutron stars have previously been observed in soft X-rays to have nearly thermal spectra at temperatures ~100 eV, which are thought to arise from the warm neutron star surfaces. We find non-trivial hard X-ray spectra well above the thermal surface predictions with archival data from the XMM-Newton and Chandra X-ray telescopes. We analyze possible systematic effects that could generate such spurious signals, such as nearby X-ray point sources and pileup of soft X-rays, but we find that the hard X-ray excesses are robust to these systematics. We also investigate possible sources of hard X-ray emission from the neutron stars and find no satisfactory explanation with known mechanisms, suggesting that a novel source of X-ray emission is at play. We do not find high-significance hard X-ray excesses from the other five Magnificent Seven isolated neutron stars.

The presence of a population of point sources in a dataset modifies the underlying neutrino-count statistics from the Poisson distribution. This deviation can be exactly quantified using the non-Poissonian template fitting technique, and in this work we present the first application this approach to the IceCube high-energy neutrino dataset. Using this method, we search in 7 years of IceCube data for point-source populations correlated with the disk of the Milky Way, the Fermi bubbles, the Schlegel, Finkbeiner, and Davis dust map, or with the isotropic extragalactic sky. No evidence for such a population is found in the data using this technique, and in the absence of a signal we establish constraints on population models with source count distribution functions that can be described by a power-law with a single break. The derived limits can be interpreted in the context of many possible source classes. In order to enhance the flexibility of the results, we publish the full posterior from our analysis, which can be used to establish limits on specific population models that would contribute to the observed IceCube neutrino flux.

Dark matter comprises the bulk of the matter in the universe but its particle nature and cosmological origin remain mysterious. Knowledge of the dark matter density distribution in the Milky Way Galaxy is crucial to both our understanding of the standard cosmological model and for grounding direct and indirect searches for the particles comprising dark matter. Current measurements of Galactic dark matter content rely on model assumptions to infer the forces acting upon stars from the distribution of observed velocities. Here, we propose to apply the precision radial velocity method, optimized in recent years for exoplanet astronomy, to measure the change in the velocity of stars over time, thereby providing a direct probe of the local gravitational potential in the Galaxy. Using numerical simulations, we develop a realistic strategy to observe the differential accelerations of stars in our Galactic neighborhood with next-generation telescopes, at the level of $10^{-8}$ cm/s$^{2}$. Our simulations show that detecting accelerations at this level with an ensemble of $10^{3}$ stars requires the effect of stellar noise on radial velocity measurements to be reduced to $<10$ cm/s. The measured stellar accelerations may then be used to extract the local dark matter density and morphological parameters of the density profile.

Utilizing the Fermi measurement of the gamma-ray spectrum toward the Galactic Center, we derive some of the strongest constraints to date on the dark matter (DM) lifetime in the mass range from hundreds of MeV to above an EeV. Our profile-likelihood based analysis relies on 413 weeks of Fermi Pass 8 data from 200 MeV to 2 TeV, along with up-to-date models for diffuse gamma-ray emission within the Milky Way. We model Galactic and extragalactic DM decay and include contributions to the DM-induced gamma-ray flux resulting from both primary emission and inverse-Compton scattering of primary electrons and positrons. For the extragalactic flux, we also calculate the spectrum associated with cascades of high-energy gamma-rays scattering off of the cosmic background radiation. We argue that a decaying DM interpretation for the 10 TeV-1 PeV neutrino flux observed by IceCube is disfavored by our constraints. Our results also challenge a decaying DM explanation of the AMS-02 positron flux. We interpret the results in terms of individual final states and in the context of simplified scenarios such as a hidden-sector glueball model.

Ultracompact dark matter (DM) minihalos at masses at and below $10^{-12}$ $M_\odot$ arise in axion DM models where the Peccei-Quinn (PQ) symmetry is broken after inflation. The minihalos arise from density perturbations that are generated from the non-trivial axion self interactions during and shortly after the collapse of the axion-string and domain-wall network. We perform high-resolution simulations of this scenario starting at the epoch before the PQ phase transition and ending at matter-radiation equality. We characterize the spectrum of primordial perturbations that are generated and comment on implications for efforts to detect axion DM. We also measure the DM density at different simulated masses and argue that the correct DM density is obtained for $m_a = 25.2 \pm 11.0$ $\mu\mathrm{eV}$.

Direct measurement of acceleration is a key scientific goal for both cosmology and exoplanets. For cosmology, the concept of redshift drift (more than 60 years old by the 2020s) could directly establish the Friedmann-Lema{\^\i}tre-Robertson-Walker model. It would increase the dark energy figure of merit by a factor of 3 beyond Stage 4 experiments, in combination with cosmic microwave background measurements. For exoplanets, the same technology required provides unprecedented radial velocity accuracy, enabling detection of Earth mass planets in the habitable zone. Other science cases include mapping the Milky Way gravitational potential and testing its dark matter distribution.

White dwarf (WD) stars may radiate keV-energy axions produced in their stellar cores. This has been extensively studied as an extra channel by which WDs may cool, with some analyses even suggesting that axions can help explain the observed WD luminosity function. We show that the radiated axions may convert into X-rays in the strong magnetic fields surrounding the WDs, leading to observable X-ray signatures. We use Suzaku observations of the WD RE J0317-853 to set the strongest constraints to-date on the combination of the axion-electron ($g_{aee}$) times axion-photon ($g_{a\gamma\gamma}$) couplings, and we show that dedicated observations of magnetic WDs by telescopes such as Chandra, XMM-Newton, and NuSTAR could increase the sensitivity to $|g_{aee} g_{a\gamma\gamma}|$ by over an order of magnitude, allowing for a definitive test of the axion-like-particle explanation of the stellar cooling anomalies.

It has been suggested that radio telescopes may be sensitive to axion dark matter that resonantly converts to radio photons in the magnetospheres surrounding neutron stars (NSs). In this work, we closely examine this possibility by calculating the radiated power from and projected sensitivity to axion dark matter conversion in ensembles of NSs within astrophysical systems like galaxies and globular clusters. We use population synthesis and evolution models to describe the spatial distributions of NSs within these systems and the distributions of NS properties. Focusing on three specific targets for illustration, the Galactic Center of the Milky Way, the globular cluster M54 in the Sagittarius dwarf galaxy, and the Andromeda galaxy, we show that narrow-band radio observations with telescopes such as the Green Bank Telescope and the future Square Kilometer Array may be able to probe the quantum chromodynamics axion over roughly two orders of magnitude in mass, starting at a fraction of a $\mu$eV.

If a component of the dark matter has dissipative interactions, it could collapse to form a thin dark disk in our Galaxy that is coplanar with the baryonic disk. It has been suggested that dark disks could explain a variety of observed phenomena, including periodic comet impacts. Using the first data release from the $Gaia$ space observatory, we search for a dark disk via its effect on stellar kinematics in the Milky Way. Our new limits disfavor the presence of a thin dark matter disk, and we present updated measurements on the total matter density in the Solar neighborhood.

The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta.

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 axion like 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$\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 $\sim$2$\sigma$ significance. A co-rotating dark disk, with lag speed $\sim$50 km$/$s, that is $\sim$20$\%$ of the local DM 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.

We use 413 weeks of publicly-available $\textit{Fermi}$ Pass 8 gamma-ray data, combined with recently-developed galaxy group catalogs, to search for evidence of dark matter annihilation in extragalactic halos. In our study, we use luminosity-based mass estimates and mass-to-concentration relations to infer the $J$-factors and associated uncertainties for hundreds of galaxy groups within a redshift range $z \lesssim 0.03$. We employ a conservative substructure boost-factor model, which only enhances the sensitivity by an $\mathcal{O}(1)$ factor. No significant evidence for dark matter annihilation is found and we exclude thermal relic cross sections for dark matter masses below $\sim$30 GeV to 95% confidence in the $b\bar{b}$ annihilation channel. These bounds are comparable to those from Milky Way dwarf spheroidal satellite galaxies. The results of our analysis increase the tension, but do not rule out, the dark matter interpretation of the Galactic Center excess. We provide a catalog of the galaxy groups used in this study and their inferred properties, which can be broadly applied to searches for extragalactic dark matter.

Dark matter in the halos surrounding galaxy groups and clusters can annihilate to high-energy photons. Recent advancements in the construction of galaxy group catalogs provide many thousands of potential extragalactic targets for dark matter. In this paper, we outline a procedure to infer the dark matter signal associated with a given galaxy group. Applying this procedure to a catalog of sources, one can create a full-sky map of the brightest extragalactic dark matter targets in the nearby Universe ($z\lesssim 0.03$), supplementing sources of dark matter annihilation from within the Local Group. As with searches for dark matter in dwarf galaxies, these extragalactic targets can be stacked together to enhance the signals associated with dark matter. We validate this procedure on mock $\textit{Fermi}$ gamma-ray data sets using a galaxy catalog constructed from the $\texttt{DarkSky}$ $N$-body cosmological simulation and demonstrate that the limits are robust, at $\mathcal{O}(1)$ levels, to systematic uncertainties on halo mass and concentration. We also quantify other sources of systematic uncertainty arising from the analysis and modeling assumptions. Our results suggest that a stacking analysis using galaxy group catalogs provides a powerful opportunity to discover extragalactic dark matter and complements existing studies of Milky Way dwarf galaxies.

We show that axion dark matter (DM) may be detectable through narrow radio lines emitted from neutron stars. The neutron star magnetosphere hosts a strong magnetic field and a plasma frequency that increases towards the neutron star surface. As the axions pass through the magnetosphere, they can resonantly convert into radio photons in a narrow region around the radius at which the plasma frequency equals the axion mass. The bandwidth of the signal is set by the small DM velocity dispersion far away from the neutron star. We solve the axion-photon mixing equations, including a full treatment of the magnetized plasma and associated anisotropic dielectric tensor, to obtain the conversion probability. We discuss possible neutron-star targets and how they may probe the QCD axion parameter space in the mass range of ~0.2-40 $\mu$eV.

We propose a novel method utilizing stellar kinematic data to detect low-mass substructure in the Milky Way's dark matter halo. By probing characteristic wakes that a passing dark matter subhalo leaves in the phase space distribution of ambient halo stars, we estimate sensitivities down to subhalo masses $\sim 10^7\,M_\odot$ or below. The detection of such subhalos would have implications for dark-matter and cosmological models that predict modifications to the halo-mass function at low halo masses. We develop an analytic formalism for describing the perturbed stellar phase-space distributions, and we demonstrate through simulations the ability to detect subhalos using the phase-space model and a likelihood framework. Our method complements existing methods for low-mass subhalo searches, such as searches for gaps in stellar streams, in that we can localize the positions and velocities of the subhalos today.

The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a wide field-of-view observatory sensitive to 500 GeV - 100 TeV gamma rays and cosmic rays. With its observations over 2/3 of the sky every day, the HAWC observatory is sensitive to a wide variety of astrophysical sources, including possible gamma rays from dark matter. Dark matter annihilation and decay in the Milky Way Galaxy should produce gamma-ray signals across many degrees on the sky. The HAWC instantaneous field-of-view of 2 sr enables observations of extended regions on the sky, such as those from dark matter in the Galactic halo. Here we show limits on the dark matter annihilation cross-section and decay lifetime from HAWC observations of the Galactic halo with 15 months of data. These are some of the most robust limits on TeV and PeV dark matter, largely insensitive to the dark matter morphology. These limits begin to constrain models in which PeV IceCube neutrinos are explained by dark matter which primarily decays into hadrons.

The $\textit{Fermi}$-LAT Collaboration recently presented a new catalog of gamma-ray sources located within the $40^{\circ} \times 40^{\circ}$ region around the Galactic Center~(Ajello et al. 2017) -- the Second Fermi Inner Galaxy (2FIG) catalog. Utilizing this catalog, they analyzed models for the spatial distribution and luminosity function of sources with a pulsar-like gamma-ray spectrum. Ajello et al. 2017 v1 also claimed to detect, in addition to a disk-like population of pulsar-like sources, an approximately 7$\sigma$ preference for an additional centrally concentrated population of pulsar-like sources, which they referred to as a "Galactic Bulge" population. Such a population would be of great interest, as it would support a pulsar interpretation of the gamma-ray excess that has long been observed in this region. In an effort to further explore the implications of this new source catalog, we attempted to reproduce the results presented by the $\textit{Fermi}$-LAT Collaboration, but failed to do so. Mimicking as closely as possible the analysis techniques undertaken in Ajello et al. 2017, we instead find that our likelihood analysis favors a very different spatial distribution and luminosity function for these sources. Most notably, our results do not exhibit a strong preference for a "Galactic Bulge" population of pulsars. Furthermore, we find that masking the regions immediately surrounding each of the 2FIG pulsar candidates does $\textit{not}$ significantly impact the spectrum or intensity of the Galactic Center gamma-ray excess. Although these results refute the claim of strong evidence for a centrally concentrated pulsar population presented in Ajello et al. 2017, they neither rule out nor provide support for the possibility that the Galactic Center excess is generated by a population of low-luminosity and currently largely unobserved pulsars.

The Milky Way's Galactic Center harbors a gamma-ray excess that is a candidate signal of annihilating dark matter. Dwarf galaxies remain predominantly dark in their expected commensurate emission. In this work we quantify the degree of consistency between these two observations through a joint likelihood analysis. In doing so we incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center Extended gamma-ray excess (GCE) detected by the Fermi Gamma-Ray Space Telescope. The preferred range of annihilation rates and masses expands when including these unknowns. Even so, using two recent determinations of the Milky Way halo's local density leave the GCE preferred region of single-channel dark matter annihilation models to be in strong tension with annihilation searches in combined dwarf galaxy analyses. A third, higher Milky Way density determination, alleviates this tension. Our joint likelihood analysis allows us to quantify this inconsistency. We provide a set of tools for testing dark matter annihilation models' consistency within this combined dataset. As an example, we test a representative inverse Compton sourced self-interacting dark matter model, which is consistent with both the GCE and dwarfs.

One of the next frontiers in dark-matter direct-detection experiments is to explore the MeV to GeV mass regime. Such light dark matter does not carry enough kinetic energy to produce an observable nuclear recoil, but it can scatter off electrons, leading to a measurable signal. We introduce a semi-analytic approach to characterize the resulting electron-scattering events in atomic and semiconductor targets, improving on previous analytic proposals that underestimate the signal at high recoil energies. We then use this procedure to study the time-dependent properties of the electron-scattering signal, including the modulation fraction, higher-harmonic modes and modulation phase. The time dependence can be distinct in a non-trivial way from the nuclear scattering case. Additionally, we show that dark-matter interactions inside the Earth can significantly distort the lab-frame phase-space distribution of sub-GeV dark matter.

This white paper summarizes the workshop "U.S. Cosmic Visions: New Ideas in Dark Matter" held at University of Maryland on March 23-25, 2017.

We present NPTFit, an open-source code package, written in python and cython, for performing non-Poissonian template fits (NPTFs). The NPTF is a recently-developed statistical procedure for characterizing the contribution of unresolved point sources (PSs) to astrophysical data sets. The NPTF was first applied to Fermi gamma-ray data to give evidence that the excess of ~GeV gamma-rays observed in the inner regions of the Milky Way likely arises from a population of sub-threshold point sources, and the NPTF has since found additional applications studying sub-threshold extragalactic sources at high Galactic latitudes. The NPTF generalizes traditional astrophysical template fits to allow for the ability to search for populations of unresolved PSs that may follow a given spatial distribution. NPTFit builds upon the framework of the fluctuation analyses developed in X-ray astronomy, and thus likely has applications beyond those demonstrated with gamma-ray data. The NPTFit package utilizes novel computational methods to perform the NPTF efficiently. The code is available at https://github.com/bsafdi/NPTFit and up-to-date and extensive documentation may be found at http://nptfit.readthedocs.io

Observations by the Fermi-LAT have uncovered a bright, spherically symmetric excess surrounding the center of the Milky Way galaxy. The spectrum of the gamma-ray excess peaks sharply at an energy ~2 GeV, exhibiting a hard spectrum at lower energies, and falls off quickly above an energy ~5 GeV. The spectrum of the excess above ~10 GeV is potentially an important discriminator between different physical models for its origin. We focus our study on observations of the gamma-ray excess at energies exceeding 10 GeV, finding: (1) a statistically significant excess remains in the energy range 9.5-47.5 GeV, which is not degenerate with known diffuse emission templates such as the Fermi Bubbles, (2) the radial profile of the excess at high energies remains relatively consistent with data near the spectral peak, (3) the data above ~5 GeV prefer a slightly greater ellipticity with a major axis oriented perpendicular to the Galactic plane. Using the recently developed non-Poissonian template fit, we find mild evidence for a point-source origin for the high-energy excess, although given the statistical and systematic uncertainties we show that a smooth origin of the high-energy emission cannot be ruled out. We discuss the implication of these findings for pulsar and dark matter models of the gamma-ray excess. Finally we provide a number of updated measurements of the gamma-ray excess, utilizing novel diffuse templates and the Pass 8 dataset.

Astrophysical sources outside the Milky Way, such as active galactic nuclei and star-forming galaxies, leave their imprint on the gamma-ray sky as nearly isotropic emission referred to as the Extragalactic Gamma-Ray Background (EGB). While the brightest of these sources may be individually resolved, their fainter counterparts contribute diffusely. In this work, we use a recently-developed analysis method, called the Non-Poissonian Template Fit, on up to 93 months of publicly-available data from the Fermi Large Area Telescope to determine the properties of the point sources that comprise the EGB. This analysis takes advantage of photon-count statistics to probe the aggregate properties of these source populations below the sensitivity threshold of published catalogs. We measure the source-count distributions and point-source intensities, as a function of energy, from 2 GeV to 2 TeV. We find that the EGB is dominated by point sources, likely blazars, in all seven energy sub-bins considered. These results have implications for the interpretation of IceCube's PeV neutrinos, which may originate from sources that contribute to the non-blazar component of the EGB. Additionally, we comment on implications for future TeV observatories such as the Cherenkov Telescope Array. We provide sky maps showing locations most likely to contain these new sources at both low (< 50 GeV) and high (> 50 GeV) energies for use in future observations and cross-correlation studies.

When ultralight axion dark matter encounters a static magnetic field, it sources an effective electric current that follows the magnetic field lines and oscillates at the axion Compton frequency. We propose a new experiment to detect this axion effective current. In the presence of axion dark matter, a large toroidal magnet will act like an oscillating current ring, whose induced magnetic flux can be measured by an external pickup loop inductively coupled to a SQUID magnetometer. We consider both resonant and broadband readout circuits and show that a broadband approach has advantages at small axion masses. We estimate the reach of this design, taking into account the irreducible sources of noise, and demonstrate potential sensitivity to axion-like dark matter with masses in the range of 10^{-14}-10^{-6} eV. In particular, both the broadband and resonant strategies can probe the QCD axion with a GUT-scale decay constant.

We present a new method to characterize unresolved point sources (PSs), generalizing traditional template fits to account for non-Poissonian photon statistics. We apply this method to Fermi Large Area Telescope gamma-ray data to characterize PS populations at high latitudes and in the Inner Galaxy. We find that PSs (resolved and unresolved) account for ~50% of the total extragalactic gamma-ray background in the energy range ~1.9 to 11.9 GeV. Within 10$^\circ$ of the Galactic Center with $|b| \geq 2^\circ$, we find that ~5-10% of the flux can be accounted for by a population of unresolved PSs, distributed consistently with the observed ~GeV gamma-ray excess in this region. The excess is fully absorbed by such a population, in preference to dark-matter annihilation. The inferred source population is dominated by near-threshold sources, which may be detectable in future searches.

Data from the Fermi Large Area Telescope suggests that there is an extended excess of GeV gamma-ray photons in the Inner Galaxy. Identifying potential astrophysical sources that contribute to this excess is an important step in verifying whether the signal originates from annihilating dark matter. In this paper, we focus on the potential contribution of unresolved point sources, such as millisecond pulsars (MSPs). We propose that the statistics of the photons---in particular, the flux probability density function (PDF) of the photon counts below the point-source detection threshold---can potentially distinguish between the dark-matter and point-source interpretations. We calculate the flux PDF via the method of generating functions for these two models of the excess. Working in the framework of Bayesian model comparison, we then demonstrate that the flux PDF can potentially provide evidence for an unresolved MSP-like point-source population.

Neutrino capture on tritium has emerged as a promising method for detecting the cosmic neutrino background (CvB). We show that relic neutrinos are captured most readily when their spin vectors are anti-aligned with the polarization axis of the tritium nuclei and when they approach along the direction of polarization. As a result, CvB observatories may measure anisotropies in the cosmic neutrino velocity and spin distributions by polarizing the tritium targets. A small dipole anisotropy in the CvB is expected due to the peculiar velocity of the lab frame with respect to the cosmic frame and due to late-time gravitational effects. The PTOLEMY experiment, a tritium observatory currently under construction, should observe a nearly isotropic background. This would serve as a strong test of the cosmological origin of a potential signal. The polarized-target measurements may also constrain non-standard neutrino interactions that would induce larger anisotropies and help discriminate between Majorana versus Dirac neutrinos.

The cosmic neutrino background (CvB), produced about one second after the Big Bang, permeates the Universe today. New technological advancements make neutrino capture on beta-decaying nuclei (NCB) a clear path forward towards the detection of the CvB. We show that gravitational focusing by the Sun causes the expected neutrino capture rate to modulate annually. The amplitude and phase of the modulation depend on the phase-space distribution of the local neutrino background, which is perturbed by structure formation. These results also apply to searches for sterile neutrinos at NCB experiments. Gravitational focusing is the only source of modulation for neutrino capture experiments, in contrast to dark-matter direct-detection searches where the Earth's time-dependent velocity relative to the Sun also plays a role.

We show that Big Bang Nucleosynthesis (BBN) significantly constrains axion-like dark matter. The axion acts like an oscillating QCD $\theta$ angle that redshifts in the early universe, increasing the neutron-proton mass difference at neutron freeze-out. An axion-like particle that couples too strongly to QCD results in the underproduction of 4He during BBN and is thus excluded. The BBN bound overlaps with much of the parameter space that would be covered by proposed searches for time-varying neutron EDMs. The QCD axion does not couple strongly enough to affect BBN.

The scattering rate at dark-matter direct-detection experiments should modulate annually due to the Earth's orbit around the Sun. The rate is typically thought to be extremized around June 1, when the relative velocity of the Earth with respect to the dark-matter wind is maximal. We point out that gravitational focusing can alter this modulation phase. Unbound dark-matter particles are focused by the Sun's gravitational potential, affecting their phase-space density in the lab frame. Gravitational focusing can result in a significant overall shift in the annual-modulation phase, which is most relevant for dark matter with low scattering speeds. The induced phase shift for light O(10) GeV dark matter may also be significant, depending on the threshold energy of the experiment.

The count rate at dark-matter direct-detection experiments should modulate annually due to the motion of the Earth around the Sun. We show that higher-frequency modulations, including daily modulation, are also present and in some cases are nearly as strong as the annual modulation. These higher-order modes are particularly relevant if (i) the dark matter is light, O(10) GeV, (ii) the scattering is inelastic, or (iii) velocity substructure is present; for these cases, the higher-frequency modes are potentially observable at current and ton-scale detectors. We derive simple expressions for the harmonic modes as functions of the astrophysical and geophysical parameters describing the Earth's orbit, using an updated expression for the Earth's velocity that corrects a common error in the literature. For an isotropic halo velocity distribution, certain ratios of the modes are approximately constant as a function of nuclear recoil energy. Anisotropic distributions can also leave observable features in the harmonic spectrum. Consequently, the higher-order harmonic modes are a powerful tool for identifying a potential signal from interactions with the Galactic dark-matter halo.