A dark photon may kinetically mix with the Standard Model photon, leading to observable cosmological signatures. The mixing is resonantly enhanced when the dark photon mass matches the primordial plasma frequency, which depends sensitively on the underlying spatial distribution of electrons. Crucially, inhomogeneities in this distribution can have a significant impact on the nature of resonant conversions. We develop and describe, for the first time, a general analytic formalism to treat resonant oscillations in the presence of inhomogeneities. Our formalism follows from the theory of level crossings of random fields and only requires knowledge of the one-point probability distribution function (PDF) of the underlying electron number density fluctuations. We validate our formalism using simulations and illustrate the photon-to-dark photon conversion probability for several different choices of PDFs that are used to characterize the low-redshift Universe.

A dark photon may kinetically mix with the ordinary photon, inducing oscillations with observable imprints on cosmology. Oscillations are resonantly enhanced if the dark photon mass equals the ordinary photon plasma mass, which tracks the free electron number density. Previous studies have assumed a homogeneous Universe; in this Letter, we introduce for the first time an analytic formalism for treating resonant oscillations in the presence of inhomogeneities of the photon plasma mass. We apply our formalism to determine constraints from Cosmic Microwave Background photons oscillating into dark photons, and from heating of the primordial plasma due to dark photon dark matter converting into low-energy photons. Including the effect of inhomogeneities demonstrates that prior homogeneous constraints are not conservative, and simultaneously extends current experimental limits into a vast new parameter space.

Astrometric weak gravitational lensing is a powerful probe of the distribution of matter on sub-Galactic scales, which harbor important information about the fundamental nature of dark matter. We propose a novel method that utilizes angular power spectra to search for the correlated pattern of apparent motions of celestial objects induced from time-dependent lensing by a population of Galactic subhalos. Application of this method to upcoming astrometric datasets will allow for the direct measurement of the properties of Galactic substructure, with implications for the underlying particle physics. We show that, with near-future astrometric observations, it may be possible to statistically detect populations of cold dark matter subhalos, compact objects, as well as density fluctuations sourced by scalar field dark matter. Currently-unconstrained parameter space will already be accessible using upcoming data from the ongoing Gaia mission.

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.

The subtle and unique imprint of dark matter substructure on extended arcs in strong lensing systems contains a wealth of information about the properties and distribution of dark matter on small scales and, consequently, about the underlying particle physics. However, teasing out this effect poses a significant challenge since the likelihood function for realistic simulations of population-level parameters is intractable. We apply recently-developed simulation-based inference techniques to the problem of substructure inference in galaxy-galaxy strong lenses. By leveraging additional information extracted from the simulator, neural networks are efficiently trained to estimate likelihood ratios associated with population-level parameters characterizing substructure. Through proof-of-principle application to simulated data, we show that these methods can provide an efficient and principled way to simultaneously analyze an ensemble of strong lenses, and can be used to mine the large sample of lensing images deliverable by near-future surveys for signatures of dark matter substructure.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

The Milky Way halo is the brightest source of dark matter annihilation on the sky. Indeed, the potential strength of the Galactic dark matter signal can supersede that expected from dwarf galaxies and galaxy groups even in regions away from the Inner Galaxy. In this paper, we present the results of a search for dark matter annihilation in the smooth Milky Way halo for $|b| > 20^\circ$ and $r < 50^\circ$ using 413 weeks of Fermi Pass 8 data within the energy range of $\sim$0.8-50 GeV. We exclude thermal dark matter with mass below $\sim$70 GeV that annihilates to $b\bar{b}$ at the 95% confidence level using the p6v11 cosmic-ray foreground model, providing the strongest limits on the annihilation cross section in this mass range. These results exclude the region of dark matter parameter space that is consistent with the excess of $\sim$GeV photons observed at the Galactic Center for the $b\bar{b}$ annihilation channel and, for the first time, start probing the $\tau^+\tau^-$ explanation. We explore how these results depend on uncertainties in the foregrounds by varying over a set of reasonable models.

We present new constraints on sub-GeV dark matter particles scattering off electrons in argon based on an analysis of ionization signal data from the DarkSide-50 detector.

We are at the dawn of a data-driven era in astrophysics and cosmology. A large number of ongoing and forthcoming experiments combined with an increasingly open approach to data availability offer great potential in unlocking some of the deepest mysteries of the Universe. Among these is understanding the nature of dark matter (DM)---one of the major unsolved problems in particle physics. Characterizing DM through its astrophysical signatures will require a robust understanding of its distribution in the sky and the use of novel statistical methods. The first part of this thesis describes the implementation of a novel statistical technique which leverages the "clumpiness" of photons originating from point sources (PSs) to derive the properties of PS populations hidden in astrophysical datasets. This is applied to data from the Fermi satellite at high latitudes ($|b| > 30$\deg) to characterize the contribution of PSs of extragalactic origin. We find that the majority of extragalactic gamma-ray emission can be ascribed to unresolved PSs having properties consistent with known sources such as active galactic nuclei. This leaves considerably less room for significant dark matter contribution. The second part of this thesis poses the question: "what is the best way to look for annihilating dark matter in extragalactic sources?" and attempts to answer it by constructing a pipeline to robustly map out the distribution of dark matter outside the Milky Way using galaxy group catalogs. This framework is then applied to Fermi data and existing group catalogs to search for annihilating dark matter in extragalactic galaxies and clusters.

We present the results of a search for dark matter WIMPs in the mass range below 20 GeV/c^2 using a target of low-radioactivity argon. The data were obtained using the DarkSide-50 apparatus at Laboratori Nazionali del Gran Sasso (LNGS). The analysis is based on the ionization signal, for which the DarkSide-50 time projection chamber is fully efficient at 0.1 keVee. The observed rate in the detector at 0.5 keVee is about 1.5 events/keVee/kg/day and is almost entirely accounted for by known background sources. We obtain a 90% C.L. exclusion limit above 1.8 GeV/c^2 for the spin-independent cross section of dark matter WIMPs on nucleons, extending the exclusion region for dark matter below previous limits in the range 1.8-6 GeV/c^2.

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.

Cosmological measurements over the next decade will enable us to shed light on the content and evolution of the Universe. Complementary measurements of the Cosmic Microwave Background (CMB) and Baryon Acoustic Oscillations are expected to allow an indirect determination of the sum of neutrino masses, within the framework of the flat $\Lambda$CDM model. However, possible deviations from $\Lambda$CDM such as a non-zero cosmological curvature or a dark energy equation of state with $w\neq -1$ would leave similar imprints on the expansion rate of the Universe and clustering of matter. We show how future CMB measurements can be combined with late-time measurements of galaxy clustering and cosmic shear from the Large Synoptic Survey Telescope to alleviate this degeneracy. Together, they are projected to reduce the uncertainty on the neutrino mass sum to 30 meV within this more general cosmological model. Achieving a 3$\sigma$ measurement of the minimal 60 meV mass (or 4$\sigma$ assuming $w=-1$) will require a five-fold improved measurement of the optical depth to reionization, obtainable through a large-scale CMB polarization measurement.

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 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.

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.

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.

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

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.