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.

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.

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.

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.