The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significantly reduces the number of classifiable stars. The approach advocated here is applicable to a much larger portion of Gaia DR2. A method known as "transfer learning" is shown to be effective through extensive testing on a set of mock Gaia catalogs that are based on the FIRE cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies. The machine is first trained on simulated data using only 5D kinematics as inputs and is then further trained on a cross-matched Gaia/RAVE data set, which improves sensitivity to properties of the real Milky Way. The result is a catalog that identifies around 767,000 accreted stars within Gaia DR2. This catalog can yield empirical insights into the merger history of the Milky Way and could be used to infer properties of the dark matter distribution.

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.

We present an in-depth exploration of the phenomenon of dynamical friction in a universe where the dark matter is composed entirely of so-called Fuzzy Dark Matter (FDM), ultralight bosons of mass $m\sim\mathcal{O}(10^{-22})\,$eV. We review the classical treatment of dynamical friction before presenting analytic results in the case of FDM for point masses, extended mass distributions, and FDM backgrounds with finite velocity dispersion. We then test these results against a large suite of fully non-linear simulations that allow us to assess the regime of applicability of the analytic results. We apply these results to a variety of astrophysical problems of interest, including infalling satellites in a galactic dark matter background, and determine that \emph{(1)}~for FDM masses $m\gtrsim 10^{-21}\, {\rm eV}\, c^{-2}$, the timing problem of the Fornax dwarf spheroidal's globular clusters is no longer solved and \emph{(2)}~the effects of FDM on the process of dynamical friction for satellites of total mass $M$ and relative velocity $v_{\rm rel}$ should require detailed numerical simulations for $\left(M/10^9~M_{\odot}\right) \left(m/10^{-22}~{\rm eV}\right)\left(100~{\rm km}~{\rm s}^{-1}/v_{\rm rel}\right) \sim 1$, parameters which would lie outside the validated range of applicability of any currently developed analytic theory, due to transient wave structures in the time-dependent regime.

There are many well-known correlations between dark matter and baryons that exist on galactic scales. These correlations can essentially be encompassed by a simple scaling relation between observed and baryonic accelerations, historically known as the Mass Discrepancy Acceleration Relation (MDAR). The existence of such a relation has prompted many theories that attempt to explain the correlations by invoking additional fundamental forces on baryons. The standard lore has been that a theory that reduces to the MDAR on galaxy scales but behaves like cold dark matter (CDM) on larger scales provides an excellent fit to data, since CDM is desirable on scales of clusters and above. However, this statement should be revised in light of recent results showing that a fundamental force that reproduces the MDAR is challenged by Milky Way dynamics. In this study, we test this claim on the example of Superfluid Dark Matter. We find that a standard CDM model is strongly preferred over a static superfluid profile. This is due to the fact that the superfluid model over-predicts vertical accelerations, even while reproducing galactic rotation curves. Our results establish an important criterion that any dark matter model must satisfy within the Milky Way.

The PTOLEMY project aims to develop a scalable design for a Cosmic Neutrino Background (CNB) detector, the first of its kind and the only one conceived that can look directly at the image of the Universe encoded in neutrino background produced in the first second after the Big Bang. The scope of the work for the next three years is to complete the conceptual design of this detector and to validate with direct measurements that the non-neutrino backgrounds are below the expected cosmological signal. In this paper we discuss in details the theoretical aspects of the experiment and its physics goals. In particular, we mainly address three issues. First we discuss the sensitivity of PTOLEMY to the standard neutrino mass scale. We then study the perspectives of the experiment to detect the CNB via neutrino capture on tritium as a function of the neutrino mass scale and the energy resolution of the apparatus. Finally, we consider an extra sterile neutrino with mass in the eV range, coupled to the active states via oscillations, which has been advocated in view of neutrino oscillation anomalies. This extra state would contribute to the tritium decay spectrum, and its properties, mass and mixing angle, could be studied by analyzing the features in the beta decay electron spectrum.

Galactic rotation curves are often considered the first robust evidence for the existence of dark matter. However, even in the presence of a dark matter halo, other galactic-scale observations, such as the Baryonic Tully-Fisher Relation and the Radial Acceleration Relation, remain challenging to explain. This has motivated long-distance, infrared modifications to gravity as an alternative to the dark matter hypothesis as well as various DM theories with similar phenomenology. In general, the standard lore has been that any model that reduces to the phenomenology of MOdified Newtonian Dynamics (MOND) on galactic scales explains essentially all galaxy-scale observables. We present a framework to test precisely this statement using local Milky Way observables, including the vertical acceleration field, the rotation curve, the baryonic surface density, and the stellar disk profile. We focus on models that predict scalar amplifications of gravity, i.e., models that increase the magnitude but do not change the direction of the gravitational acceleration. We find that models of this type are disfavored relative to a simple dark matter halo model because the Milky Way data requires a substantial amplification of the radial acceleration with little amplification of the vertical acceleration. We conclude that models which result in a MOND-like force struggle to simultaneously explain both the rotational velocity and vertical motion of nearby stars in the Milky Way.

We use the distribution of accreted stars in SDSS-Gaia DR2 to demonstrate that a non-trivial fraction of the dark matter halo within Galactocentric radii of 7.5-10 kpc and $|z| > 2.5$ kpc is in substructure, and thus may not be in equilibrium. Using a mixture likelihood analysis, we separate the contributions of an old, isotropic stellar halo and a younger anisotropic population. The latter dominates and is uniform within the region studied. It can be explained as the tidal debris of a disrupted massive satellite on a highly radial orbit, and is consistent with mounting evidence from recent studies. Simulations that track the tidal debris from such mergers find that the dark matter traces the kinematics of its stellar counterpart. If so, our results indicate that a component of the nearby dark matter halo that is sourced by luminous satellites is in kinematic substructure referred to as debris flow. These results challenge the Standard Halo Model, which is discrepant with the distribution recovered from the stellar data, and have important ramifications for the interpretation of direct detection experiments.

In Ostdiek et al. (2019), we developed a deep neural network classifier that only relies on phase-space information to obtain a catalog of accreted stars based on the second data release of Gaia (DR2). In this paper, we apply two clustering algorithms to identify velocity substructure within this catalog. We focus on the subset of stars with line-of-sight velocity measurements that fall in the range of Galactocentric radii $r \in [6.5, 9.5]$ kpc and vertical distances $|z| < 3$ kpc. Known structures such as Gaia Enceladus and the Helmi stream are identified. The largest previously-unknown structure, Nyx, first introduced in Necib et al. (2019a), is a vast stream consisting of at least 500 stars in the region of interest. This study displays the power of the machine learning approach by not only successfully identifying known features, but also discovering new kinematic structures that may shed light on the merger history of the Milky Way.

Massive dwarf galaxies that merge with the Milky Way on prograde orbits can be dragged into the disk plane before being completely disrupted. Such mergers can contribute to an accreted stellar disk and a dark matter disk. We present evidence for Nyx, a vast new stellar stream in the vicinity of the Sun, that may provide the first indication that such an event occurred in the Milky Way. We identify about 500 stars that have coherent radial and prograde motion in this stream using a catalog of accreted stars built by applying deep learning methods to the second Gaia data release. Nyx is concentrated within $\pm 2$ kpc of the Galactic midplane and spans the full radial range studied (6.5-9.5 kpc). The kinematics of Nyx stars are distinct from those of both the thin and thick disk. In particular, its rotational speed lags the disk by $\sim 80$ km/s and its stars follow more eccentric orbits. A small number of Nyx stars have chemical abundances or inferred ages; from these, we deduce that Nyx stars have a peak metallicity of [Fe/H] $\sim -0.5$ and ages $\sim $10-13 Gyr. Taken together with the kinematic observations, these results strongly favor the interpretation that Nyx is the remnant of a disrupted dwarf galaxy. To further justify this interpretation, we explicitly demonstrate that metal-rich, prograde streams like Nyx can be found in the disk plane of Milky Way-like galaxies using the FIRE hydrodynamic simulations. Future spectroscopic studies will be able to validate whether Nyx stars originate from a single progenitor.

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 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 \textit{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 $\sim 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 $\sim6\sigma$ and predicts fewer high-speed particles. As a result, the Standard Halo Model overpredicts the nuclear scattering rate for dark matter masses below $\sim 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.

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 briefly comment on the paper "On the correlation between the local dark matter and stellar velocities" (Bozorgnia et al. 2018), which purports to test our previous claims regarding the correlation of stellar and dark matter kinematics (Herzog-Arbeitman et al. 2017; Necib et al. 2018). The conclusions of these papers are not accurately presented in Bozorgnia et al. (2018), and we thus clarify them here.

The Gaia era opens new possibilities for discovering the remnants of disrupted satellite galaxies in the Solar neighborhood. If the population of local accreted stars is correlated with the dark matter sourced by the same mergers, one can then map the dark matter distribution directly. Using two cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies from the Latte suite of Fire-2 simulations, we find a strong correlation between the velocity distribution of stars and dark matter at the solar circle that were accreted from luminous satellites. This correspondence holds for dark matter that is either relaxed or in kinematic substructure called debris flow, and is consistent between two simulated hosts with different merger histories. The correspondence is more problematic for streams because of possible spatial offsets between the dark matter and stars. We demonstrate how to reconstruct the dark matter velocity distribution from the observed properties of the accreted stellar population by properly accounting for the ratio of stars to dark matter contributed by individual mergers. After demonstrating this method using the Fire-2 simulations, we apply it to the Milky Way and use it to recover the dark matter velocity distribution associated with the recently discovered stellar debris field in the Solar neighborhood. Based on results from Gaia, we estimate that $42 ^{+26}_{-22}\%$ of the local dark matter that is accreted from luminous mergers is in debris flow.

We present a detailed description of the electromagnetic filter for the PTOLEMY project to directly detect the Cosmic Neutrino Background (CNB). Starting with an initial estimate for the orbital magnetic moment, the higher-order drift process of ExB is configured to balance the gradient-B drift motion of the electron in such a way as to guide the trajectory into the standing voltage potential along the mid-plane of the filter. As a function of drift distance along the length of the filter, the filter zooms in with exponentially increasing precision on the transverse velocity component of the electron kinetic energy. This yields a linear dimension for the total filter length that is exceptionally compact compared to previous techniques for electromagnetic filtering. The parallel velocity component of the electron kinetic energy oscillates in an electrostatic harmonic trap as the electron drifts along the length of the filter. An analysis of the phase-space volume conservation validates the expected behavior of the filter from the adiabatic invariance of the orbital magnetic moment and energy conservation following Liouville's theorem for Hamiltonian systems.

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

We propose to achieve the proof-of-principle of the PTOLEMY project to directly detect the Cosmic Neutrino Background (CNB). Each of the technological challenges described in [1,2] will be targeted and hopefully solved by the use of the latest experimental developments and profiting from the low background environment provided by the LNGS underground site. The first phase will focus on the graphene technology for a tritium target and the demonstration of TES microcalorimetry with an energy resolution of better than 0.05 eV for low energy electrons. These technologies will be evaluated using the PTOLEMY prototype, proposed for underground installation, using precision HV controls to step down the kinematic energy of endpoint electrons to match the calorimeter dynamic range and rate capabilities. The second phase will produce a novel implementation of the EM filter that is scalable to the full target size and which demonstrates intrinsic triggering capability for selecting endpoint electrons. Concurrent with the CNB program, we plan to exploit and develop the unique properties of graphene to implement an intermediate program for direct directional detection of MeV dark matter [3,4]. This program will evaluate the radio-purity and scalability of the graphene fabrication process with the goal of using recently identified ultra-high radio-purity CO2 sources. The direct detection of the CNB is a snapshot of early universe dynamics recorded by the thermal relic neutrino yield taken at a time that predates the epochs of Big Bang Nucleosynthesis, the Cosmic Microwave Background and the recession of galaxies (Hubble Expansion). Big Bang neutrinos are believed to have a central role in the evolution of the Universe and a direct measurement with PTOLEMY will unequivocally establish the extent to which these predictions match present-day neutrino densities.

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 Milky Way dark matter halo is formed from the accretion of smaller subhalos. These sub-units also harbor stars---typically old and metal-poor---that are deposited in the Galactic inner regions by disruption events. In this Letter, we show that the dark matter and metal-poor stars in the Solar neighborhood share similar kinematics due to their common origin. Using the high-resolution Eris simulation, which traces the evolution of both the dark matter and baryons in a realistic Milky-Way analog galaxy, we demonstrate that metal-poor stars are indeed effective tracers for the local, virialized dark matter velocity distribution. The local dark matter velocities can therefore be inferred from observations of the stellar halo made by the Sloan Digital Sky Survey within 4 kpc of the Sun. This empirical distribution differs from the Standard Halo Model in important ways and suggests that the bounds on the spin-independent scattering cross section may be weakened for dark matter masses below $\sim$10 GeV. Data from Gaia will allow us to further refine the expected distribution for the smooth dark matter component, and to test for the presence of local substructure.

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 propose two-dimensional materials as targets for direct detection of dark matter. Using graphene as an example, we focus on the case where dark matter scattering deposits sufficient energy on a valence-band electron to eject it from the target. We show that the sensitivity of graphene to dark matter of MeV to GeV mass can be comparable, for similar exposure and background levels, to that of semiconductor targets such as silicon and germanium. Moreover, a two-dimensional target is an excellent directional detector, as the ejected electron retains information about the angular dependence of the incident dark matter particle. This proposal can be implemented by the PTOLEMY experiment, presenting for the first time an opportunity for directional detection of sub-GeV dark matter.

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.

This report, based on the Dark Sectors workshop at SLAC in April 2016, summarizes the scientific importance of searches for dark sector dark matter and forces at masses beneath the weak-scale, the status of this broad international field, the important milestones motivating future exploration, and promising experimental opportunities to reach these milestones over the next 5-10 years.

Rotation curve measurements provided the first strong indication that a significant fraction of matter in the Universe is non-baryonic. Since then, a tremendous amount of progress has been made on both the theoretical and experimental fronts in the search for this missing matter, which we now know constitutes nearly 85% of the Universe's matter density. These series of lectures, first given at the TASI 2015 summer school, provide an introduction to the basics of dark matter physics. They are geared for the advanced undergraduate or graduate student interested in pursuing research in high-energy physics. The primary goal is to build an understanding of how observations constrain the assumptions that can be made about the astro- and particle physics properties of dark matter. The lectures begin by delineating the basic assumptions that can be inferred about dark matter from rotation curves. A detailed discussion of thermal dark matter follows, motivating Weakly Interacting Massive Particles, as well as lighter-mass alternatives. As an application of these concepts, the phenomenology of direct and indirect detection experiments is discussed in detail.

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.

Tidal debris from infalling satellites can leave observable structure in the phase-space distribution of the Galactic halo. Such substructure can be manifest in the spatial and/or velocity distributions of the stars in the halo. This paper focuses on a class of substructure that is purely kinematic in nature, with no accompanying spatial features. To study its properties, we use a simulated stellar halo created by dynamically populating the Via Lactea II high-resolution N-body simulation with stars. A significant fraction of the stars in the inner halo of Via Lactea share a common speed and metallicity, despite the fact that they are spatially diffuse. We argue that this kinematic substructure is a generic feature of tidal debris from older mergers and may explain the detection of radial-velocity substructure in the inner halo made by the Sloan Extension for Galactic Understanding and Exploration. The GAIA satellite, which will provide the proper motions of an unprecedented number of stars, should further characterize the kinematic substructure in the inner halo. Our study of the Via Lactea simulation suggests that the stellar halo can be used to map the speed distribution of the local dark-matter halo, which has important consequences for dark-matter direct-detection experiments.

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.

A fermion triplet of SU(2)_L - a wino - is a well-motivated dark matter candidate. This work shows that present-day wino annihilations are constrained by indirect detection experiments, with the strongest limits coming from H.E.S.S. and Fermi. The bounds on wino dark matter are presented as a function of mass for two scenarios: thermal (winos constitute a subdominant component of the dark matter for masses less than 3.1 TeV) and non-thermal (winos comprise all the dark matter). Assuming the NFW halo model, the H.E.S.S. search for gamma-ray lines excludes the 3.1 TeV thermal wino; the combined H.E.S.S. and Fermi results completely exclude the non-thermal scenario. Uncertainties in the exclusions are explored. Indirect detection may provide the only probe for models of anomaly plus gravity mediation where the wino is the lightest superpartner and scalars reside at the 100 TeV scale.

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.

Direct detection experiments, which are designed to detect the scattering of dark matter off nuclei in detectors, are a critical component in the search for the Universe's missing matter. The count rate in these experiments should experience an annual modulation due to the relative motion of the Earth around the Sun. This modulation, not present for most known background sources, is critical for solidifying the origin of a potential signal as dark matter. In this article, we review the physics of annual modulation, discussing the practical formulae needed to interpret a modulating signal. We focus on how the modulation spectrum changes depending on the particle and astrophysics models for the dark matter. For standard assumptions, the count rate has a cosine dependence with time, with a maximum in June and a minimum in December. Well-motivated generalizations of these models, however, can affect both the phase and amplitude of the modulation. We show how a measurement of an annually modulating signal could teach us about the presence of substructure in the Galactic halo or about the interactions between dark and baryonic matter. Although primarily a theoretical review, we briefly discuss the current experimental situation for annual modulation and future experimental directions.

We show that subhalos falling into the Milky Way create a flow of tidally-stripped debris particles near the galactic center with characteristic velocity behavior. In the Via Lactea-II N-body simulation, this unvirialized component constitutes a few percent of the local density and has velocities peaked at 340 km/s in the solar neighborhood. Such velocity substructure has important implications for surveys of low-metallicity stars, as well as direct detection experiments sensitive to dark matter with large scattering thresholds.

There is evidence for a 130 GeV gamma-ray line at the Galactic Center in the Fermi Large Area Telescope data. Dark matter candidates that explain this feature should also annihilate to Standard Model particles, resulting in a continuous spectrum of photons. To study this continuum, we analyze the Fermi data down to 5 GeV, restricted to the inner 3 degrees of the Galaxy. We place a strong bound on the ratio of continuum photons to monochromatic line photons that is independent of uncertainties in the dark matter density profile. The derived constraints exclude neutralino dark matter as an explanation for the line.

Tidal stripping of dark matter from subhalos falling into the Milky Way produces narrow, cold tidal streams as well as more spatially extended "debris flows" in the form of shells, sheets, and plumes. Here we focus on the debris flow in the Via Lactea II simulation, and show that this incompletely phase-mixed material exhibits distinctive high velocity behavior. Unlike tidal streams, which may not necessarily intersect the Earth's location, debris flow is spatially uniform at 8 kpc and thus guaranteed to be present in the dark matter flux incident on direct detection experiments. At Earth-frame speeds greater than 450 km/s, debris flow comprises more than half of the dark matter at the Sun's location, and up to 80% at even higher speeds. Therefore, debris flow is most important for experiments that are particularly sensitive to the high speed tail of the dark matter distribution, such as searches for light or inelastic dark matter or experiments with directional sensitivity. We show that debris flow yields a distinctive recoil energy spectrum and a broadening of the distribution of incidence direction.

We analyze the recently released CoGeNT data with a focus on their time-dependent properties. Using various statistical techniques, we confirm the presence of modulation in the data, and find a significant component at high (E_{ee} > 1.5$ keVee) energies. We find that standard elastic WIMPs in a Maxwellian halo do not provide a good description of the modulation. We consider the possibility of non-standard halos, using halo independent techniques, and find a good agreement with the DAMA modulation for Q_{Na} \approx 0.3, but disfavoring interpretations with Q_{Na} = 0.5. The same techniques indicate that CDMS-Ge should see an O(1) modulation, and XENON100 should have seen 10-30 events (based upon the modulation in the 1.5-3.1 keVee range), unless L_{eff} is smaller than recent measurements. Models such as inelastic dark matter provide a good fit to the modulation, but not the spectrum. We note that tensions with XENON could be alleviated in such models if the peak is dominantly in April, when XENON data are not available due to noise.

Dark matter density profiles based upon Lambda-CDM cosmology motivate an ansatz velocity distribution function with fewer high velocity particles than the Maxwell-Boltzmann distribution or proposed variants. The high velocity tail of the distribution is determined by the outer slope of the dark matter halo, the large radius behavior of the Galactic dark matter density. N-body simulations of Galactic halos reproduce the high velocity behavior of this ansatz. Predictions for direct detection rates are dramatically affected for models where the threshold scattering velocity is within 30% of the escape velocity.

The XENON100 and CRESST experiments will directly test the inelastic dark matter explanation for DAMA's 8.9? sigma anomaly. This article discusses how predictions for direct detection experiments depend on uncertainties in quenching factor measurements, the dark matter interaction with the Standard Model and the halo velocity distribution. When these uncertainties are accounted for, an order of magnitude variation is found in the number of expected events at CRESST and XENON100.

Inelastic dark matter reconciles the DAMA anomaly with other null direct detection experiments and points to a non-minimal structure in the dark matter sector. In addition to the dominant inelastic interaction, dark matter scattering may have a subdominant elastic component. If these elastic interactions are suppressed at low momentum transfer, they will have similar nuclear recoil spectra to inelastic scattering events. While upcoming direct detection experiments will see strong signals from such models, they may not be able to unambiguously determine the presence of the subdominant elastic scattering from the recoil spectra alone. We show that directional detection experiments can separate elastic and inelastic scattering events and discover the underlying dynamics of dark matter models.

We present Submillimeter Array observations of the Herbig Ae star HD169142 in 1.3 millimeter continuum emission and 12CO J=2-1 line emission at 1.5 arcsecond resolution that reveal a circumstellar disk. The continuum emission is centered on the star position and resolved, and provides a mass estimate of about 0.02 solar masses for the disk. The CO images show patterns in position and velocity that are well matched by a disk in Keplerian rotation with low inclination to the line-of-sight. We use radiative transfer calculations based on a flared, passive disk model to constrain the disk parameters by comparison to the spectral line emission. The derived disk radius is 235 AU, and the inclination is 13 degrees. The model also necessitates modest depletion of the CO molecules, similar to that found in Keplerian disks around T Tauri stars.