We present a catalog of 20,977 extended low-surface-brightness galaxies (LSBGs) identified in ~ 5000 deg$^2$ from the first three years of imaging data from the Dark Energy Survey (DES). Based on a single-component S\'ersic model fit, we define extended LSBGs as galaxies with $g$-band effective radii $R_{\scriptsize{eff}} > 2.5"$ and mean surface brightness $\bar{\mu}_{\scriptsize{eff}}(g) > 24.3$ mag arcsec$^{-2}$. We find that the distribution of LSBGs is strongly bimodal in $(g-r)$ vs. $(g-i$) color space. We divide our sample into red ($g-i \geq 0.59$) and blue ($g-i<0.59$) galaxies and study the properties of the two populations. Redder LSBGs are more clustered than their blue counterparts, and are correlated with the distribution of nearby ($z < 0.10$) bright galaxies. Red LSBGs constitute $\sim 35\%$ of our LSBG sample, and $\sim 30\%$ of these are located within 1 deg of low-redshift galaxy groups and clusters (compared to $\sim 8\%$ of the blue LSBGs). For nine of the most prominent galaxy groups and clusters, we calculate the physical properties of associated LSBGs assuming a redshift derived from the host system. In these systems, we identify 108 objects that can be classified as ultra-diffuse galaxies, defined as LSBGs with projected physical effective radii $R_{\scriptsize{eff}} > 1.5$ kpc. The wide-area sample of LSBGs in DES can be used to test the role of environment on models of LSBG formation and evolution.

The radial spatial distribution of low-mass satellites around a Milky Way (MW)-like host is an important benchmark for simulations of small-scale structure. The distribution is sensitive to the disruption of subhalos by the central disk and can indicate whether the disruption observed in simulations of MW analogs is artificial (i.e., numeric) or physical in origin. We consider a sample of 12 well-surveyed satellite systems of MW-like hosts in the Local Volume that are complete to $M_V<-9$ and within 150 projected kpc. We investigate the radial distribution of satellites and compare with $\Lambda$CDM cosmological simulations, including big-box cosmological simulations and high resolution zoom in simulations of MW sized halos. We find that the observed satellites are significantly more centrally concentrated than the simulated systems. Several of the observed hosts, including the MW, are $\sim2\sigma$ outliers relative to the simulated hosts in being too concentrated, while none of the observed hosts are less centrally concentrated than the simulations. This result is robust to different ways of measuring the radial concentration. We find that this discrepancy is more significant for bright, $M_V<-12$ satellites, suggestive that this is not the result of observational incompleteness. We argue that the discrepancy is possibly due to artificial disruption in the simulations, but, if so, this has important ramifications for what stellar to halo mass relation is allowed in the low-mass regime by the observed abundance of satellites.

Low-mass satellites around Milky Way (MW)-like galaxies are important probes of small scale structure and galaxy formation. However, confirmation of satellite candidates with distance measurements remains a key barrier to fast progress in the Local Volume (LV). We measure the surface brightness fluctuation (SBF) distances to recently cataloged candidate dwarf satellites around 10 massive hosts within $D<12$ Mpc to confirm association. The satellite systems of these hosts are complete and mostly cleaned of contaminants down to $M_g{\sim}-9$ to $-10$, within the area of the search footprints. Joining this sample with hosts surveyed to comparable or better completeness in the literature, we explore how well cosmological simulations combined with common stellar to halo mass relations (SHMR) match observed satellite luminosity functions in the classical satellite luminosity regime. Adopting a SHMR that matches hydrodynamic simulations, the predicted overall satellite abundance agrees well with the observations. The MW is remarkably typical in its luminosity function amongst LV hosts. Contrary to recent results, we find that the host-to-host scatter predicted by the model is in close agreement with the scatter between the observed systems, once the different masses of the observed systems are taken into account. However, we find significant evidence that the observed systems have more bright and fewer faint satellites than the SHMR model predicts, necessitating a higher normalization of the SHMR around halo masses of $10^{11}$ \msun\ than present in common SHMRs. These results demonstrate the utility of nearby satellite systems in inferring the galaxy-subhalo connection in the low-mass regime.

Understanding the non-linear dynamics of satellite halos (a.k.a. ``sub-halos'') is important for predicting the abundance and distribution of dark matter substructures and satellite galaxies, and for distinguishing among microphysical dark matter models using observations. Typically, modeling these dynamics requires large N-body simulations with high resolution. Semi-analytic models can provide a more efficient way to describe the key physical processes such as dynamical friction, tidal mass loss, and tidal heating, with only a few free parameters. In this work, we present a fast Monte Carlo Markov Chain fitting approach to explore the parameter space of such a sub-halo non-linear evolution model. We use the dynamical models described in an earlier work and calibrate the models to two sets of high-resolution cold dark matter N-body simulations, ELVIS and Caterpillar. Compared to previous calibrations that used manual parameter tuning, our approach provides a more robust way to determine the best-fit parameters and their posterior probabilities. We find that jointly fitting for the sub-halo mass and maximum velocity functions can break the degeneracy between tidal stripping and tidal heating parameters, as well as providing better constraints on the strength of dynamical friction. We show that our semi-analytic simulation can accurately reproduce N-body simulations statistics, and that the calibration results for the two sets of N-body simulations agree at 95% confidence level. Dynamical models calibrated in this work will be important for future dark matter substructure studies.

We present the first satellite system of the Large Binocular Telescope Satellites Of Nearby Galaxies Survey (LBT-SONG), a survey to characterize the close satellite populations of Large Magellanic Cloud to Milky Way-mass, star-forming galaxies in the Local Volume. In this paper, we describe our satellite finding and completeness measurement methodology and apply this framework to NGC 628, an isolated galaxy with ~1/4 the stellar mass of the Milky Way. We present two new dwarf satellite galaxy candidates: NGC 628 dwA, and dwB with $\text{M}_{\text{V}}$ = $-12.2$ and $-7.7$, respectively. NGC 628 dwA is a classical dwarf while NGC 628 dwB is a low-luminosity galaxy that appears to have been quenched after reionization. Completeness corrections indicate that the presence of these two satellites is consistent with CDM predictions. The satellite colors indicate that the galaxies are neither actively star-forming nor do they have the purely ancient stellar populations characteristic of ultrafaint dwarfs. Instead, and consistent with our previous work on the NGC 4214 system, they show signs of recent quenching, further indicating that environmental quenching can play a role in modifying satellite populations even for hosts smaller than the Milky Way.

The magnifications of compact-source lenses are extremely sensitive to the presence of low mass dark matter halos along the entire sight line from the source to the observer. Traditionally, the study of dark matter structure in compact-source strong gravitational lenses has been limited to radio-loud systems, as the radio emission is extended and thus unaffected by microlensing which can mimic the signal of dark matter structure. An alternate approach is to measure quasar nuclear-narrow line emission, which is free from microlensing and present in virtually all quasar lenses. In this paper, we double the number of systems which can be used for gravitational lensing analyses by presenting measurements of narrow-line emission from a sample of 8 quadruply imaged quasar lens systems, WGD J0405-3308, HS 0810+2554, RX J0911+0551, SDSS J1330+1810, PS J1606-2333, WFI 2026-4536, WFI 2033-4723 and WGD J2038-4008. We describe our updated grism spectral modelling pipeline, which we use to measure narrow-line fluxes with uncertainties of 2-10\%, presented here. We fit the lensed image positions with smooth mass models and demonstrate that these models fail to produce the observed distribution of image fluxes over the entire sample of lenses. Furthermore, typical deviations are larger than those expected from macromodel uncertainties. This discrepancy indicates the presence of perturbations caused by small-scale dark matter structure. The interpretation of this result in terms of dark matter models is presented in a companion paper.

Our multi-view metric learning framework enables robust characterization of star categories by directly learning to discriminate in a multi-faceted feature space, thus, eliminating the need to combine feature representations prior to fitting the machine learning model. We also demonstrate how to extend standard multi-view learning, which employs multiple vectorized views, to the matrix-variate case which allows very novel variable star signature representations. The performance of our proposed methods is evaluated on the UCR Starlight and LINEAR datasets. Both the vector and matrix-variate versions of our multi-view learning framework perform favorably --- demonstrating the ability to discriminate variable star categories.

We present the construction of a novel time-domain signature extraction methodology and the development of a supporting supervised pattern detection algorithm. We focus on the targeted identification of eclipsing binaries that demonstrate a feature known as the O'Connell effect. Our proposed methodology maps stellar variable observations to a new representation known as distribution fields (DFs). Given this novel representation, we develop a metric learning technique directly on the DF space that is capable of specifically identifying our stars of interest. The metric is tuned on a set of labeled eclipsing binary data from the Kepler survey, targeting particular systems exhibiting the O'Connell effect. The result is a conservative selection of 124 potential targets of interest out of the Villanova Eclipsing Binary Catalog. Our framework demonstrates favorable performance on Kepler eclipsing binary data, taking a crucial step in preparing the way for large-scale data volumes from next-generation telescopes such as LSST and SKA.

We report the discovery of a $>1^\circ$ ($\sim50$ kpc) long stellar tidal stream emanating from the dwarf galaxy DDO 44, a likely satellite of Local Volume galaxy NGC 2403 located $\sim70$ kpc in projection from its companion. NGC 2403 is a roughly Large Magellanic Cloud stellar-mass galaxy 3 Mpc away, residing at the outer limits of the M 81 group. We are mapping a large region around NGC 2403 as part of our MADCASH (Magellanic Analogs' Dwarf Companions and Stellar Halos) survey, reaching point source depths (90% completeness) of ($g, i$) = (26.5, 26.2). Density maps of old, metal-poor RGB stars reveal tidal streams extending on two sides of DDO 44, with the streams directed toward NGC 2403. We estimate total luminosities of the original DDO 44 system (dwarf and streams combined) to be $M_{i, \rm{tot}} = -13.4$ and $M_{g, \rm{tot}} = -12.6$, with $\sim25-30\%$ of the luminosity in the streams. Analogs of $\sim$LMC-mass hosts with massive tidally disrupting satellites are rare in the Illustris simulations, especially at large separations such as that of DDO 44. The few analogs that are present in the models suggest that even low-mass hosts can efficiently quench their massive satellites.

Using N-body simulations of the Large Magellanic Cloud (LMC's) passage through the Milky Way (MW), tailored to reproduce observed kinematic properties of both galaxies, we show that the high-speed tail of the Solar Neighborhood dark matter distribution is overwhelmingly of LMC origin. Two populations contribute at high speeds: 1) Particles that were once bound to the LMC, and 2) MW halo particles that have been accelerated owing to the response of the halo to the recent passage of the LMC. These particles reach speeds of 700-900 km/s with respect to the Earth, above the local escape speed of the MW. The high-speed particles follow trajectories similar to the Solar reflex motion, with peak velocities reached in June. For low-mass dark matter, these high-speed particles can dominate the signal in direct-detection experiments, extending the reach of the experiments to lower mass and elastic scattering cross sections even with existing data sets. Our study shows that even non-disrupted MW satellite galaxies can leave a significant dark-matter footprint in the Solar Neighborhood.

We investigate the case for environmental quenching of the Fornax-mass satellite DDO 113, which lies only 9 kpc in projection from its host, the Large-Magellanic-Cloud-mass galaxy NGC 4214. DDO 113 was quenched about 1 Gyr ago and is virtually gas-free, while analogs in the field are predominantly star-forming and gas-rich. We use deep imaging obtained with the Large Binocular Telescope to show that DDO 113 exhibits no evidence of tidal disruption to a surface brightness of $\mu_V\sim29$ mag $\text{arcsec}^{-2}$, based on both unresolved emission and resolved stars. Mass-analogs of DDO 113 in Illustris-1 with similar hosts, small projected separations, and no significant tidal stripping first fell into their host halo 2--6 Gyr ago, showing that tidal features (or lack thereof) can be used to constrain infall times in systems where there are few other constraints on the orbit of the satellite. With the infall time setting the clock for environmental quenching mechanisms, we investigate the plausibility of several such mechanisms. We find that strangulation, the cessation of cold gas inflows, is likely the dominant quenching mechanism for DDO 113, requiring a time-averaged mass-loading factor of $\eta=6-11$ for star-formation-driven outflows that is consistent with theoretical and observational constraints. Motivated by recent numerical work, we connect DDO 113's strangulation to the presence of a cool circumgalactic medium (CGM) around NGC 4214. This discovery shows that the CGM of low-mass galaxies can affect their satellites significantly and motivates further work on understanding the baryon cycle in low-mass galaxies.

Stellar halos offer fossil evidence for hierarchical structure formation. Since halo assembly is predicted to be scale-free, stellar halos around low-mass galaxies constrain properties such as star formation in the accreted subhalos and the formation of dwarf galaxies. However, few observational searches for stellar halos in dwarfs exist. Here we present gi photometry of resolved stars in isolated Local Group dwarf irregular galaxy IC 1613 ($M_{\star} \sim 10^8 M_{\odot})$. These Subaru/Hyper Suprime-Cam observations are the widest and deepest of IC 1613 to date. We measure surface density profiles of young main-sequence, intermediate to old red giant branch, and ancient horizontal branch stars outside of 12' ($\sim 2.6$ kpc; 2.5 half-light radii) from the IC 1613 center. All of the populations extend to ~24' ($\sim 5.2$ kpc; 5 half-light radii), with the older populations best fit by a broken exponential in these outer regions. Comparison with earlier studies sensitive to IC 1613's inner regions shows that the density of old stellar populations steepens substantially with distance from the center; we trace the $g$-band effective surface brightness to an extremely faint limit of $\sim 33.7$ mag arcsec$^{-2}$. Conversely, the distribution of younger stars follows a single, shallow exponential profile in the outer regions, demonstrating different formation channels for the younger and older components of IC 1613. The outermost, intermediate-age and old stars have properties consistent with those expected for accreted stellar halos, though future observational and theoretical work is needed to definitively distinguish this scenario from other possibilities.

Critical probes of dark matter come from tests of its elastic scattering with nuclei. The results are typically assumed to be model-independent, meaning that the form of the potential need not be specified and that the cross sections on different nuclear targets can be simply related to the cross section on nucleons. For point-like spin-independent scattering, the assumed scaling relation is $\sigma_{\chi A} \propto A^2 \mu_A^2 \sigma_{\chi N}\propto A^4 \sigma_{\chi N}$, where the $A^2$ comes from coherence and the $\mu_A^2\simeq A^2 m_N^2$ from kinematics for $m_\chi\gg m_A$. Here we calculate where model independence ends, i.e., where the cross section becomes so large that it violates its defining assumptions. We show that the assumed scaling relations generically fail for dark matter-nucleus cross sections $\sigma_{\chi A} \sim 10^{-32}-10^{-27}\;\text{cm}^2$, significantly below the geometric sizes of nuclei, and well within the regime probed by underground detectors. Last, we show on theoretical grounds, and in light of existing limits on light mediators, that point-like dark matter cannot have $\sigma_{\chi N}\gtrsim10^{-25}\;\text{cm}^2$, above which many claimed constraints originate from cosmology and astrophysics. The most viable way to have such large cross sections is composite dark matter, which introduces significant additional model dependence through the choice of form factor. All prior limits on dark matter with cross sections $\sigma_{\chi N}>10^{-32}\;\text{cm}^2$ with $m_\chi\gtrsim 1\;\text{GeV}$ must therefore be re-evaluated and reinterpreted.

A census of the satellite population around dwarf galaxy primary hosts in environments outside the Local Group is essential to understanding $\Lambda$CDM galaxy formation and evolution on the smallest scales. We present deep optical Hubble Space Telescope imaging of the gas-rich, faint dwarf galaxy Antlia B ($M_V = -9.4$) -- a likely satellite of NGC 3109 ($D = 1.3$ Mpc) -- discovered as part of our ongoing survey of primary host galaxies similar to the Magellanic Clouds. We derive a new tip of the red giant branch (TRGB) distance of $D = 1.35 \pm 0.06$ Mpc ($m-M = 25.65 \pm 0.10$), consistent with membership in the nearby NGC 3109 dwarf association. The color-magnitude diagram shows both a prominent old, metal-poor stellar component and confirms a small population of young, blue stars with ages $\lesssim 1$ Gyr. We use the color-magnitude diagram fitting algorithm MATCH to derive the star formation history and find that it is consistent with the typical dwarf irregular or transitional dwarf galaxy (dTrans) in the Local Group. Antlia B shows relatively constant stellar mass growth for the first $\sim 10-11$ Gyr and almost no growth in the last $\sim 2-3$ Gyr. Despite being gas-rich, Antlia B shows no evidence of active star formation (i.e., no H$\alpha$ emission) and should therefore be classified as a dTrans dwarf. Both Antlia B and the Antlia dwarf (dTrans) are likely satellites of NGC 3109 suggesting that the cessation of ongoing star formation in these galaxies may be environmentally driven. Future work studying the gas kinematics and distribution in Antlia B will explore this scenario in greater detail. Our work highlights the fact that detailed studies of nearby dwarf galaxies in a variety of environments may continue to shed light on the processes that drive the star formation history and evolution of dwarf galaxies more generally.

The current disparity in computational knowledge is a critical hindrance to the diversity and success of the field. Recommendations are outlined for policies and funding models to enable the growth and retention of a new generation of computational researchers that reflect the demographics of the undergraduate population in Astronomy and Physics.

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

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We discuss how astrophysical observations with the Maunakea Spectroscopic Explorer (MSE), a high-multiplexity (about 4300 fibers), wide field-of-view (1.5 square degree), large telescope aperture (11.25 m) facility, can probe the particle nature of dark matter. MSE will conduct a suite of surveys that will provide critical input for determinations of the mass function, phase-space distribution, and internal density profiles of dark matter halos across all mass scales. N-body and hydrodynamical simulations of cold, warm, fuzzy and self-interacting dark matter suggest that non-trivial dynamics in the dark sector could have left an imprint on structure formation. Analysed within these frameworks, the extensive and unprecedented datasets produced by MSE will be used to search for deviations away from cold and collisionless dark matter model. MSE will provide an improved estimate of the local density of dark matter, critical for direct detection experiments, and will improve estimates of the J-factor for indirect searches through self-annihilation or decay into Standard Model particles. MSE will determine the impact of low mass substructures on the dynamics of Milky Way stellar streams in velocity space, and will allow for estimates of the density profiles of the dark matter halos of Milky Way dwarf galaxies using more than an order of magnitude more tracers. In the low redshift Universe, MSE will provide critical redshifts to pin down the luminosity functions of vast numbers of satellite systems, and MSE will be an essential component of future strong lensing measurements to constrain the halo mass function. Across nearly all mass scales, the improvements offered by MSE, in comparison to other facilities, are such that the relevant analyses are limited by systematics rather than statistics.

The Sun is an excellent laboratory for astroparticle physics but remains poorly understood at GeV--TeV energies. Despite the immense relevance for both cosmic-ray propagation and dark matter searches, only in recent years has the Sun become a target for precision gamma-ray astronomy with the Fermi-LAT instrument. Among the most surprising results from the observations is a hard excess of GeV gamma-ray flux that strongly anti-correlates with solar activity, especially at the highest energies accessible to Fermi-LAT. Most of the observed properties of the gamma-ray emission cannot be explained by existing models of cosmic-ray interactions with the solar atmosphere. GeV--TeV gamma-ray observations of the Sun spanning an entire solar cycle would provide key insights into the origin of these gamma rays, and consequently improve our understanding of the Sun's environment as well as the foregrounds for new physics searches, such as dark matter. These can be complemented with new observations with neutrinos and cosmic rays. Together these observations make the Sun a new testing ground for particle physics in dynamic environments.

The 2020s are poised to continue the past two decades of significant advances based on observations of dwarf galaxies in the nearby universe. Upcoming wide-field photometric surveys will probe substantially deeper than previous data sets, pushing the discovery frontier for new dwarf galaxies to fainter magnitudes, lower surface brightnesses, and larger distances. These dwarfs will be compelling targets for testing models of galaxy formation and cosmology, including the properties of dark matter and possible modifications to gravity. However, most of the science that can be extracted from nearby dwarf galaxies relies on spectroscopy with large telescopes. We suggest that maximizing the scientific impact of near-future imaging surveys will require both major spectroscopic surveys on 6-10m telescopes and multiplexed spectroscopy with even larger apertures.

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.

We analyze the Sun as a source for the indirect detection of dark matter through a search for gamma rays from the solar disk. Capture of dark matter by elastic interactions with the solar nuclei followed by annihilation to long-lived mediators can produce a detectable gamma-ray flux. We search three years of data from the High Altitude Water Cherenkov Observatory and find no statistically significant detection of TeV gamma-ray emission from the Sun. Using this, we constrain the spin-dependent elastic scattering cross section of dark matter with protons for dark matter masses above 1 TeV, assuming an unstable mediator with a favorable lifetime. The results complement constraints obtained from Fermi-LAT observations of the Sun and together cover WIMP masses between 4 GeV and $10^6$ GeV. The cross section constraints for mediator decays to gamma rays can be as strong as $\sim10^{-45}$ cm$^{-2}$, which is more than four orders of magnitude stronger than current direct-detection experiments for 1 TeV dark matter mass. The cross-section constraints at higher masses are even better, nearly 7 orders of magnitude better than the current direct-detection constraints for 100 TeV dark matter mass. This demonstration of sensitivity encourages detailed development of theoretical models in light of these powerful new constraints.

Steady gamma-ray emission up to at least 200 GeV has been detected from the solar disk in the Fermi-LAT data, with the brightest, hardest emission occurring during solar minimum. The likely cause is hadronic cosmic rays undergoing collisions in the Sun's atmosphere after being redirected from ingoing to outgoing in magnetic fields, though the exact mechanism is not understood. An important new test of the gamma-ray production mechanism will follow from observations at higher energies. Only the High Altitude Water Cherenkov (HAWC) Observatory has the required sensitivity to effectively probe the Sun in the TeV range. Using three years of HAWC data from November 2014 to December 2017, just prior to the solar minimum, we search for 1--100 TeV gamma rays from the solar disk. No evidence of a signal is observed, and we set strong upper limits on the flux at a few $10^{-12}$ TeV$^{-1}$ cm$^{-2}$ s$^{-1}$ at 1 TeV. Our limit, which is the most constraining result on TeV gamma rays from the Sun, is $\sim10\%$ of the theoretical maximum flux (based on a model where all incoming cosmic rays produce outgoing photons), which in turn is comparable to the Fermi-LAT data near 100 GeV. The prospects for a first TeV detection of the Sun by HAWC are especially high during solar minimum, which began in early 2018.

The solar disk is a bright source of multi-GeV gamma rays, due to the interactions of hadronic cosmic rays with the solar atmosphere. However, the underlying production mechanism is not understood, except that its efficiency must be greatly enhanced by magnetic fields that redirect some cosmic rays from ingoing to outgoing before they interact. To elucidate the nature of this emission, we perform a new analysis of solar atmospheric gamma rays with 9 years of Fermi-LAT data, which spans nearly the full 11-year solar cycle. We detect significant gamma-ray emission from the solar disk from 1 GeV up to $\gtrsim200$ GeV. The overall gamma-ray spectrum is much harder ($\sim E_{\gamma}^{-2.2}$) than the cosmic-ray spectrum ($\sim E_{\rm CR}^{-2.7}$). We find a clear anticorrelation between the solar cycle phase and the gamma-ray flux between 1-10 GeV. Surprisingly, we observe a spectral dip between $\sim$30-50 GeV in an otherwise power-law spectrum. This was not predicted, is not understood, and may provide crucial clues to the gamma-ray emission mechanism. The flux above 100 GeV, which is brightest during the solar minimum, poses exciting opportunities for HAWC, LHAASO, IceCube, and KM3NeT.

The nature of dark matter is one of the most pressing questions in particle physics. Yet all our present knowledge of the dark sector to date comes from its gravitational interactions with astrophysical systems. Moreover, astronomical results still have immense potential to constrain the particle properties of dark matter. We introduce a simple 2D parameter space which classifies models in terms of a particle physics interaction strength and a characteristic astrophysical scale on which new physics appears, in order to facilitate communication between the fields of particle physics and astronomy. We survey the known astrophysical anomalies that are suggestive of non-trivial dark matter particle physics, and present a theoretical and observational program for future astrophysical measurements that will shed light on the nature of dark matter.

Golovich et al. 2017b presents an optical imaging and spectroscopic survey of 29 radio relic merging galaxy clusters. In this paper, we study this survey to identify substructure and quantify the dynamics of the mergers. Using a combined photometric and spectroscopic approach, we identify the minimum number of substructures in each system to describe the galaxy populations and estimate the line of sight velocity difference between likely merging subclusters. We find that the line-of-sight velocity components of the mergers are typically small compared with the maximum three dimensional relative velocity (usually $<1000$ km s$^{-1}$ and often consistent with zero). This suggests that the merger axes of these systems are generally in or near the plane of the sky matching findings in magneto-hydrodynamical simulations. In 28 of the 29 systems we identify substructures in the galaxy population aligned with the radio relic(s) and presumed associated merger induced shock. From this ensemble, we identify eight systems to include in a `gold' sample that is prime for further observation, modeling, and simulation study. Additional papers will present weak lensing mass maps and dynamical modeling for each merging system, ultimately leading to new insight into a wide range of astrophysical phenomena at some of the largest scales in the universe.

A critical challenge to the cold dark matter (CDM) paradigm is that there are fewer satellites observed around the Milky Way than found in simulations of dark matter substructure. We show that there is a match between the observed satellite counts corrected by the detection efficiency of the Sloan Digital Sky Survey (for luminosities $L \gtrsim$ 340 L$_\odot$) and the number of luminous satellites predicted by CDM, assuming an empirical relation between stellar mass and halo mass. The "missing satellites problem", cast in terms of number counts, is thus solved. We also show that warm dark matter models with a thermal relic mass smaller than 4 keV are in tension with satellite counts, putting pressure on the sterile neutrino interpretation of recent X-ray observations. Importantly, the total number of Milky Way satellites depends sensitively on the spatial distribution of satellites, possibly leading to a "too many satellites" problem. Measurements of completely dark halos below $10^8$ M$_\odot$, achievable with substructure lensing and stellar stream perturbations, are the next frontier for tests of CDM.

The observed multi-GeV gamma-ray emission from the solar disk --- sourced by hadronic cosmic rays interacting with gas, and affected by complex magnetic fields --- is not understood. Utilizing an improved analysis of the Fermi-LAT data that includes the first resolved imaging of the disk, we find strong evidence that this emission is produced by two separate mechanisms. Between 2010-2017 (the rise to and fall from solar maximum), the gamma-ray emission is dominated by a polar component. Between 2008-2009 (solar minimum) this component remains present, but the total emission is instead dominated by a new equatorial component with a brighter flux and harder spectrum. Most strikingly, although 6 gamma rays above 100 GeV are observed during the 1.4 years of solar minimum, none are observed during the next 7.8 years. These features, along with a 30-50 GeV spectral dip which will be discussed in a companion paper, were not anticipated by theory. To understand the underlying physics, Fermi and HAWC observations of the imminent Cycle 25 solar minimum are crucial.

As is well known, dark matter direct detection experiments will ultimately be limited by a "neutrino floor," due to the scattering of nuclei by MeV neutrinos from, e.g., nuclear fusion in the Sun. Here we point out the existence of a new "neutrino floor" that will similarly limit indirect detection with the Sun, due to high-energy neutrinos from cosmic-ray interactions with the solar atmosphere. We have two key findings. First, solar atmospheric neutrinos $\lesssim 1$ TeV cause a sensitivity floor for standard WIMP scenarios, for which higher-energy neutrinos are absorbed in the Sun. This floor will be reached once the present sensitivity is improved by just one order of magnitude. Second, for neutrinos $\gtrsim 1$ TeV, which can be isolated by muon energy loss rate, solar atmospheric neutrinos should soon be detectable in IceCube. Discovery will help probe the complicated effects of solar magnetic fields on cosmic rays. These events will be backgrounds to WIMP scenarios with long-lived mediators, for which higher-energy neutrinos can escape from the Sun.

Multi-band photometric and multi-object spectroscopic surveys of merging galaxy clusters allow for the characterization of the distributions of constituent dark matter and galaxy populations, constraints on the dynamics of the merging subclusters, and an understanding of galaxy evolution of member galaxies. We present deep photometric observations from Subaru/SuprimeCam and a catalog of $\sim$5400 spectroscopic cluster members from Keck/DEIMOS across 29 merging galaxy clusters ranging in redshift from $z=0.07$ to $0.55$. The ensemble is compiled based on the presence of radio relics, which highlight cluster scale collisionless shocks in the intra-cluster medium. Together with the spectroscopic and photometric information, the velocities, timescales, and geometries of the respective merging events may be tightly constrained. In this preliminary analysis, the velocity distributions of 28 of the 29 clusters are shown to be well fit by single Gaussians. This indicates that radio relic mergers largely occur transverse to the line of sight and/or near apocenter. In this paper, we present our optical and spectroscopic surveys, preliminary results, and a discussion of the value of radio relic mergers for developing accurate dynamical models of each system.

We present deep, wide-field Subaru Hyper Suprime-Cam photometry of two recently discovered satellites of the Milky Way (MW): Columba I and Triangulum II. The color magnitude diagrams of both objects point to exclusively old and metal-poor stellar populations. We re-derive structural parameters and luminosities of these satellites, and find $M_{\rm V, Col~I} = -4.2\pm0.2$ for Col I and $M_{\rm V, Tri~II} = -1.2\pm0.4$ for Tri II, with corresponding half-light radii of $r_{\rm h, Col~I} = 117\pm17$ pc and $r_{\rm h, Tri~II} = 21\pm4$ pc. The properties of both systems are consistent with observed scaling relations for MW dwarf galaxies. Based on archival data, we derive upper limits on the neutral gas content of these dwarfs, and find that they lack HI, as do the majority of observed satellites within the MW virial radius. Neither satellite shows evidence of tidal stripping in the form of extensions or distortions in matched-filter stellar density maps or surface density profiles. However, the smaller Tri II system is relatively metal-rich for its luminosity (compared to other MW satellites), possibly because it has been tidally stripped. Through a suite of orbit simulations, we show that Tri II is approaching pericenter of its eccentric orbit, a stage at which tidal debris is unlikely to be seen. In addition, we find that Tri II may be on its first infall into the MW, which helps explain its unique properties among MW dwarfs. Further evidence that Tri II is likely an ultra-faint dwarf comes from its stellar mass function, which is similar to those of other MW dwarfs.

Recent discovery of many dwarf satellite galaxies in the direction of the Small and Large Magellanic Clouds (SMC and LMC) provokes questions of their origins, and what they can reveal about galaxy evolution theory. Here, we predict the satellite stellar mass function of Magellanic Cloud-mass host galaxies using abundance matching and reionization models applied to the \textit{Caterpillar} simulations. Specifically focusing on the volume within $50$~kpc of the LMC, we predict a mean of 4-8 satellites with stellar mass $M_* > 10^4\, \mathrm{M_\odot}$, and 3-4 satellites with $80 < M_* \leq 3000\, \mathrm{M_\odot}$. Surprisingly, all $12$ currently known satellite candidates have stellar masses of $80 < M_* \leq 3000\, \mathrm{M_\odot}$. Reconciling the dearth of large satellites and profusion of small satellites is challenging and may require a combination of a major modification of the \mstarmhalo relationship (steep, but with an abrupt flattening at $10^3\, \mathrm{M_\odot}$), late reionization for the Local Group ($z_{\rm{reion}} \lesssim 9$ preferred), and/or strong tidal stripping. We can more robustly predict that $\sim 53\%$ of satellites within this volume were accreted together with the LMC and SMC, and $\sim 47\%$ were only ever Milky Way satellites. Observing satellites of isolated LMC-sized field galaxies is essential to placing the LMC in context, and to better constrain the \mstarmhalo relationship. Modeling known LMC-sized galaxies within $8$ Mpc, we predict 1-6 (2-12) satellites with $M_* > 10^5\, \mathrm{M_\odot}$ ($M_* > 10^4\, \mathrm{M_\odot}$) within the virial volume of each, and 1-3 (1-7) within a single $1.5^{\circ}$ diameter field of view, making their discovery likely.

A recent surge in the discovery of new ultrafaint dwarf satellites of the Milky Way has inspired the idea of searching for faint satellites, $10^3\, \mathrm{M_{\odot}} < M_* < 10^6\, \mathrm{M_{\odot}}$, around less massive field galaxies in the Local Group. Such satellites would be subject to weaker environmental influences than Milky Way satellites, and could lead to new insights on low mass galaxy formation. In this paper, we predict the number of luminous satellites expected around field dwarf galaxies by applying several abundance matching models and a reionization model to the dark-matter only \textit{Caterpillar} simulation suite. For three of the four abundance matching models used, we find a $>99\%$ chance that at least one satellite with stellar mass $M_*> 10^5\, \mathrm{M_{\odot}}$ exists around the combined five Local Group field dwarf galaxies with the largest stellar mass. When considering satellites with $M_*> 10^4\, \mathrm{M_{\odot}}$, we predict a combined $5-25$ satellites for the five largest field dwarfs, and $10-50$ for the whole Local Group field dwarf population. Because of the relatively small number of predicted dwarfs, and their extended spatial distribution, a large fraction each Local Group dwarf's virial volume will need to be surveyed to guarantee discoveries. We compute the predicted number of satellites in a given field of view of specific Local Group galaxies, as a function of minimum satellite luminosity, and explicitly obtain such values for the Solitary Local dwarfs survey. Uncertainties in abundance matching and reionization models are large, implying that comprehensive searches could lead to refinements of both models.

Dwarf galaxies are known to have remarkably low star formation efficiency due to strong feedback. Adopting the dwarf galaxies of the Milky Way as a laboratory, we explore a flexible semi-analytic galaxy formation model to understand how the feedback processes shape the satellite galaxies of the Milky Way. Using Markov-Chain Monte-Carlo, we exhaustively search a large parameter space of the model and rigorously show that the general wisdom of strong outflows as the primary feedback mechanism cannot simultaneously explain the stellar mass function and the mass--metallicity relation of the Milky Way satellites. An extended model that assumes that a fraction of baryons is prevented from collapsing into low-mass halos in the first place can be accurately constrained to simultaneously reproduce those observations. The inference suggests that two different physical mechanisms are needed to explain the two different data sets. In particular, moderate outflows with weak halo mass dependence are needed to explain the mass--metallicity relation, and prevention of baryons falling into shallow gravitational potentials of low-mass halos (e.g. "pre-heating") is needed to explain the low stellar mass fraction for a given subhalo mass.

The Sun is a bright source of GeV gamma rays, due to cosmic rays interacting with solar matter and photons. Key aspects of the underlying processes remain mysterious. The emission in the TeV range, for which there are neither observational nor theoretical studies, could provide crucial clues. The new experiments HAWC (running) and LHAASO (planned) can look at the Sun with unprecedented sensitivity. In this paper, we predict the very high-energy (up to 1000 TeV) gamma-ray flux from the solar disk and halo, due to cosmic-ray hadrons and electrons ($e^++e^-$), respectively. We neglect solar magnetic effects, which is valid at TeV energies; at lower energies, this gives a theoretical lower bound on the disk flux and a theoretical upper bound on the halo flux. We show that the solar-halo gamma-ray flux allows the first test of the $\sim 5$--70 TeV cosmic-ray electron spectrum. Further, we show that HAWC can immediately make an even stronger test with nondirectional observations of cosmic-ray electrons. Together, these gamma-ray and electron studies will provide new insights about the local density of cosmic rays and their interactions with the Sun and its magnetic environment. These studies will also be an important input to tests of new physics, including 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.

Strong gravitational lensing provides a powerful test of Cold Dark Matter (CDM) as it enables the detection and mass measurement of low mass haloes even if they do not contain baryons. Compact lensed sources such as Active Galactic Nuclei (AGN) are particularly sensitive to perturbing subhalos, but their use as a test of CDM has been limited by the small number of systems which have significant radio emission which is extended enough avoid significant lensing by stars in the plane of the lens galaxy, and red enough to be minimally affected by differential dust extinction. Narrow-line emission is a promising alternative as it is also extended and, unlike radio, detectable in virtually all optically selected AGN lenses. We present first results from a WFC3 grism narrow-line survey of lensed quasars, for the quadruply lensed AGN HE0435-1223. Using a forward modelling pipeline which enables us to robustly account for spatial blending, we measure the [OIII] 5007 \AA~ flux ratios of the four images. We find that the [OIII] fluxes and positions are well fit by a simple smooth mass model for the main lens. Our data rule out a $M_{600}>10^{8} (10^{7.2}) M_\odot$ NFW perturber projected within $\sim$1\farcs0 (0\farcs1) arcseconds of each of the lensed images, where $M_{600}$ is the perturber mass within its central 600 pc. The non-detection is broadly consistent with the expectations of $\Lambda$CDM for a single system. The sensitivity achieved demonstrates that powerful limits on the nature of dark matter can be obtained with the analysis of $\sim20$ narrow-line lenses.

We study the atomic physics and the astrophysical implications of a model in which the dark matter is the analog of hydrogen in a secluded sector. The self interactions between dark matter particles include both elastic scatterings as well as inelastic processes due to a hyperfine transition. The self-interaction cross sections are computed by numerically solving the coupled Schr\"{o}dinger equations for this system. We show that these self interactions exhibit the right velocity dependence to explain the low dark matter density cores seen in small galaxies while being consistent with all constraints from observations of clusters of galaxies. For a viable solution, the dark hydrogen mass has to be in 10--100 GeV range and the dark fine-structure constant has to be larger than 0.02. Precisely for this range of parameters, we show that significant cooling losses may occur due to inelastic excitations to the hyperfine state and subsequent decays, with implications for the evolution of low-mass halos and the early growth of supermassive black holes. Cooling from excitations to higher $n$ levels of dark hydrogen and subsequent decays is possible at the cluster scale, with a strong dependence on halo mass. Finally, we show that the minimum halo mass is in the range of $10^{3.5}$ to $10^7 M_\odot$ for the viable regions of parameter space, significantly larger than the typical predictions for weakly-interacting dark matter models. This pattern of observables in cosmological structure formation is unique to this model, making it possible to rule in or rule out hidden sector hydrogen as a viable dark matter model.

We investigate the effects of self-interacting dark matter (SIDM) on the tidal stripping and evaporation of satellite galaxies in a Milky Way-like host. We use a suite of five zoom-in, dark-matter-only simulations, two with velocity-independent SIDM cross sections, two with velocity-dependent SIDM cross sections, and one cold dark matter simulation for comparison. After carefully assigning stellar mass to satellites at infall, we find that stars are stripped at a higher rate in SIDM than in CDM. In contrast, the total bound dark matter mass loss rate is minimally affected, with subhalo evaporation having negligible effects on satellites for viable SIDM models. Centrally located stars in SIDM haloes disperse out to larger radii as cores grow. Consequently, the half-light radius of satellites increases, stars become more vulnerable to tidal stripping, and the stellar mass function is suppressed. We find that the ratio of core radius to tidal radius accurately predicts the relative strength of enhanced SIDM stellar stripping. Velocity-independent SIDM models show a modest increase in the stellar stripping effect with satellite mass, whereas velocity-dependent SIDM models show a large increase in this effect towards lower masses, making observations of ultra-faint dwarfs prime targets for distinguishing between and constraining SIDM models. Due to small cores in the largest satellites of velocity-dependent SIDM, no identifiable imprint is left on the all-sky properties of the stellar halo. While our results focus on SIDM, the main physical mechanism of enhanced tidal stripping of stars apply similarly to satellites with cores formed via other means.

Merging galaxy clusters have been touted as one of the best probes for constraining self-interacting dark matter, but few simulations exist to back up this claim. We simulate equal mass mergers of 10$^{15}$ M$_\odot$ halos, like the El Gordo and Sausage clusters, with cosmologically-motivated halo and merger parameters, and with velocity-independent dark-matter self-interactions. Although the standard lore for merging clusters is that self-interactions lead to large separations between the galaxy and dark-matter distributions, we find that maximal galaxy-dark matter offsets of $\lesssim~20$~kpc form for a self-interaction cross section of $\sigma_\text{SI}/m_\chi$ = 1 cm$^2$/g. This is an order of magnitude smaller than those measured in observed equal mass and near equal mass mergers, and is likely to be even smaller for lower-mass systems. While competitive cross-section constraints are thus unlikely to emerge from offsets, we find other signatures of self-interactions which are more promising. Intriguingly, we find that after dark matter halos coalesce, the collisionless galaxies (and especially the Brightest Cluster Galaxy [BGC]) oscillate around the center of the merger remnant on stable orbits of 100 kpc for $\sigma_\text{SI}/m_\chi = 1$~cm$^2$/g for at least several Gyr, well after the clusters have relaxed. If BCG miscentering in relaxed clusters remains a robust prediction of SIDM under the addition of gas physics, substructure, merger mass ratios (e.g., 10:1 like the Bullet Cluster), and complex cosmological merger histories, the observed BCG offsets may constrain $\sigma_\text{SI}/m_\chi \lesssim$ 0.1 cm$^2$/g---the tightest constraint yet.

We report the discovery of the faintest known dwarf galaxy satellite of an LMC stellar-mass host beyond the Local Group, based on deep imaging with Subaru/Hyper Suprime-Cam. MADCASH J074238+652501-dw lies $\sim$35 kpc in projection from NGC 2403, a dwarf spiral galaxy at $D$$\approx$3.2 Mpc. This new dwarf has $M_{g} = -7.4\pm0.4$ and a half-light radius of $168\pm70$ pc, at the calculated distance of $3.39\pm0.41$ Mpc. The color-magnitude diagram reveals no evidence of young stellar populations, suggesting that MADCASH J074238+652501-dw is an old, metal-poor dwarf similar to low luminosity dwarfs in the Local Group. The lack of either detected HI gas ($M_{\rm HI}/L_{V} < 0.69 M_\odot/L_\odot$, based on Green Bank Telescope observations) or $GALEX$ NUV/FUV flux enhancement is consistent with a lack of young stars. This is the first result from the MADCASH (Magellanic Analog Dwarf Companions And Stellar Halos) survey, which is conducting a census of the stellar substructure and faint satellites in the halos of Local Volume LMC analogs via resolved stellar populations. Models predict a total of $\sim$4-10 satellites at least as massive as MADCASH J074238+652501-dw around a host with the mass of NGC 2403, with 2-3 within our field of view, slightly more than the one such satellite observed in our footprint.

Many properties of the Milky Way's dark matter halo, including its mass assembly history, concentration, and subhalo population, remain poorly constrained. We explore the connection between these properties of the Milky Way and its satellite galaxy population, especially the implication of the presence of the Magellanic Clouds for the properties of the Milky Way halo. Using a suite of high-resolution $N$-body simulations of Milky Way-mass halos with a fixed final Mvir ~ 10^{12.1}Msun, we find that the presence of Magellanic Cloud-like satellites strongly correlates with the assembly history, concentration, and subhalo population of the host halo, such that Milky Way-mass systems with Magellanic Clouds have lower concentration, more rapid recent accretion, and more massive subhalos than typical halos of the same mass. Using a flexible semi-analytic galaxy formation model that is tuned to reproduce the stellar mass function of the classical dwarf galaxies of the Milky Way with Markov-Chain Monte-Carlo, we show that adopting host halos with different mass-assembly histories and concentrations can lead to different best-fit models for galaxy-formation physics, especially for the strength of feedback. These biases arise because the presence of the Magellanic Clouds boosts the overall population of high-mass subhalos, thus requiring a different stellar-mass-to-halo-mass ratio to match the data. These biases also lead to significant differences in the mass--metallicity relation, the kinematics of low-mass satellites, the number counts of small satellites associated with the Magellanic Clouds, and the stellar mass of Milky Way itself. Observations of these galaxy properties can thus provide useful constraints on the properties of the Milky Way halo.

The solar disk is a bright gamma-ray source. Surprisingly, its flux is about one order of magnitude higher than predicted. As a first step toward understanding the physical origin of this discrepancy, we perform a new analysis in 1-100 GeV using 6 years of public Fermi-LAT data. Compared to the previous analysis by the Fermi Collaboration, who analyzed 1.5 years of data and detected the solar disk in 0.1-10 GeV, we find two new and significant results: 1. In the 1-10 GeV flux (detected at $>5\sigma$), we discover a significant time variation that anticorrelates with solar activity. 2. We detect gamma rays in 10-30 GeV at $>5\sigma$, and in 30-100 GeV at $> 2\sigma$. The time variation strongly indicates that solar-disk gamma rays are induced by cosmic rays and that solar atmospheric magnetic fields play an important role. Our results provide essential clues for understanding the underlying gamma-ray production processes, which may allow new probes of solar atmospheric magnetic fields, cosmic rays in the solar system, and possible new physics. Finally, we show that the Sun is a promising new target for ground-based TeV gamma-ray telescopes such as HAWC and LHAASO.

Cosmological observations indicate that most of the matter in the Universe is Dark Matter. Dark Matter in the form of Weakly Interacting Massive Particles (WIMPs) can be detected directly, via its elastic scattering off target nuclei. Most current direct detection experiments only measure the energy of the recoiling nuclei. However, directional detection experiments are sensitive to the direction of the nuclear recoil as well. Due to the Sun's motion with respect to the Galactic rest frame, the directional recoil rate has a dipole feature, peaking around the direction of the Solar motion. This provides a powerful tool for demonstrating the Galactic origin of nuclear recoils and hence unambiguously detecting Dark Matter. Furthermore, the directional recoil distribution depends on the WIMP mass, scattering cross section and local velocity distribution. Therefore, with a large number of recoil events it will be possible to study the physics of Dark Matter in terms of particle and astrophysical properties. We review the potential of directional detectors for detecting and characterizing WIMPs.

With the advent of digital astronomy, new benefits and new challenges have been presented to the modern day astronomer. No longer can the astronomer rely on manual processing, instead the profession as a whole has begun to adopt more advanced computational means. This paper focuses on the construction and application of a novel time-domain signature extraction methodology and the development of a supporting supervised pattern classification algorithm for the identification of variable stars. A methodology for the reduction of stellar variable observations (time-domain data) into a novel feature space representation is introduced. The methodology presented will be referred to as Slotted Symbolic Markov Modeling (SSMM) and has a number of advantages which will be demonstrated to be beneficial; specifically to the supervised classification of stellar variables. It will be shown that the methodology outperformed a baseline standard methodology on a standardized set of stellar light curve data. The performance on a set of data derived from the LINEAR dataset will also be shown

Identifying the true theory of dark matter depends crucially on accurately characterizing interactions of dark matter (DM) with other species. In the context of DM direct detection, we present a study of the prospects for correctly identifying the low-energy effective DM-nucleus scattering operators connected to UV-complete models of DM-quark interactions. We take a census of plausible UV-complete interaction models with different low-energy leading-order DM-nuclear responses. For each model (corresponding to different spin-, momentum-, and velocity-dependent responses), we create a large number of realizations of recoil-energy spectra, and use Bayesian methods to investigate the probability that experiments will be able to select the correct scattering model within a broad set of competing scattering hypotheses. We conclude that agnostic analysis of a strong signal (such as Generation-2 would see if cross sections are just below the current limits) seen on xenon and germanium experiments is likely to correctly identify momentum dependence of the dominant response, ruling out models with either "heavy" or "light" mediators, and enabling downselection of allowed models. However, a unique determination of the correct UV completion will critically depend on the availability of measurements from a wider variety of nuclear targets, including iodine or fluorine. We investigate how model-selection prospects depend on the energy window available for the analysis. In addition, we discuss accuracy of the DM particle mass determination under a wide variety of scattering models, and investigate impact of the specific types of particle-physics uncertainties on prospects for model selection.

The Rosetta lander Philae successfully landed on the nucleus of comet 67P/Churyumov-Gerasimenko on 12 November 2014. Philae carries the Dust Impact Monitor (DIM) on board, which is part of the Surface Electric Sounding and Acoustic Monitoring Experiment (SESAME). DIM employs piezoelectric PZT sensors to detect impacts by sub-millimeter and millimeter-sized ice and dust particles that are emitted from the nucleus and transported into the cometary coma. The DIM sensor measures dynamical data like flux and the directionality of the impacting particles. Mass and speed of the particles can be constrained assuming density and elastic particle properties. DIM was operated during three mission phases of Philae at the comet: (1) Before Philae's separation from Rosetta at distances of about 9.6 km, 11.8 km, and 25.3 km from the nucleus barycenter. In this mission phase particles released from the nucleus on radial trajectories remained undetectable because of significant obscuration by the structures of Rosetta, and no dust particles were indeed detected. (2) During Philae's descent to its nominal landing site Agilkia, DIM detected one approximately millimeter-sized particle at a distance of 5.0 km from the nucleus' barycenter, corresponding to an altitude of 2.4 km from the surface. This is the closest ever dust detection at a cometary nucleus by a dedicated in-situ dust detector. (3) At Philae's final landing site, Abydos, DIM detected no dust impact which may be due to low cometary activity in the vicinity of Philae, or due to shading by obstacles close to Philae, or both. Laboratory calibration experiments showed that the material properties of the detected particle are compatible with a porous particle having a bulk density of approximately $250\, \mathrm{kg\,m^{-3}}$. The particle could have been lifted off from the comet's surface by sublimating water ice.

We investigate the effect of self-interacting dark matter (SIDM) on the density profiles of $V_{\rm max} \simeq 40~km~s^{-1}$ isolated dwarf dark matter halos -- the scale of relevance for the too big to fail problem (TBTF) -- using very high-resolution cosmological zoom simulations. Each halo has millions of particles within its virial radius. We find that SIDM models with cross sections per unit mass spanning the range \sigma/m = $0.5 - 50$ $cm^2~g^{-1}$ alleviate TBTF and produce constant density cores of size 300-1000 pc, comparable to the half-light radii of $M_\star$ ~ $10^{5-7}$ $M_\odot$ dwarfs. The largest, lowest density cores develop for cross sections in the middle of this range, \sigma/m ~ $5-10~cm^2~g^{-1}$. Our largest SIDM cross section run (\sigma/m = $50~cm^2~g^{-1}$) develops a slightly denser core owing to mild core-collapse behavior, but it remains less dense than the CDM case and retains a constant density core profile. Our work suggests that SIDM cross sections as large or larger than $50~cm^2~g^{-1}$ remain viable on velocity scales of dwarf galaxies ($v_{\rm rms}$ ~ $40~km~s^{-1}$). The range of SIDM cross sections that alleviate TBTF and the cusp/core problem spans at least two orders of magnitude and therefore need not be particularly fine-tuned.

We predict the spatial distribution and number of Milky Way dwarf galaxies to be discovered in the DES and LSST surveys, by completeness correcting the observed SDSS dwarf population. We apply most massive in the past, earliest forming, and earliest infall toy models to a set of dark matter-only simulated Milky Way/M31 halo pairs from Exploring the Local Volume In Simulations (ELVIS). The observed spatial distribution of Milky Way dwarfs in the LSST-era will discriminate between the earliest infall and other simplified models for how dwarf galaxies populate dark matter subhalos. Inclusive of all toy models and simulations, at 90% confidence we predict a total of 37-114 L $\gtrsim 10^3$L$_{\odot}$ dwarfs and 131-782 L $\lesssim 10^3$L$_{\odot}$ dwarfs within 300 kpc. These numbers of L $\gtrsim 10^3$L$_{\odot}$ dwarfs are dramatically lower than previous predictions, owing primarily to our use of updated detection limits and the decreasing number of SDSS dwarfs discovered per sky area. For an effective $r_{\rm limit}$ of 25.8 mag, we predict: 3-13 L $\gtrsim 10^3$L$_{\odot}$ and 9-99 L $\lesssim 10^3$L$_{\odot}$ dwarfs for DES, and 18-53 L $\gtrsim 10^3$L$_{\odot}$ and 53-307 L $\lesssim 10^3$L$_{\odot}$ dwarfs for LSST. These enormous predicted ranges ensure a coming decade of near-field excitement with these next generation surveys.

We study prospects of dark-matter direct-detection searches for probing non-relativistic effective theory for WIMP-baryon scattering. We simulate a large set of noisy recoil-energy spectra for different scattering scenarios (beyond the standard momentum-independent contact interaction), for Generation 2 and futuristic experiments. We analyze these simulations and quantify the probability of successfully identifying the operator governing the scattering, if a WIMP signal is observed. We find that the success rate depends on a combination of factors: the WIMP mass, the mediator mass, the type of interaction, and the experimental energy window. For example, for a 20 GeV WIMP, Generation 2 is only likely to identify the right operator if the interaction is Coulomb-like, and is unlikely to do so in any other case. For a WIMP with a mass of 200 GeV or higher, success is almost guaranteed. We also find that, regardless of the scattering model and the WIMP parameters, a single Generation 2 experiment is unlikely to successfully discern the momentum dependence of the underlying operator on its own, but prospects improve drastically when experiments with different target materials and energy windows are analyzed jointly. Furthermore, we examine the quality of parameter estimation and degeneracies in the multi-dimensional parameter space of the effective theory. We find in particular that the resulting WIMP mass estimates can be severely biased if data are analyzed assuming the standard (momentum-independent) operator while the actual operator has momentum-dependence. Finally, we evaluate the ultimate reach of direct detection, finding that the prospects for successful operator selection prior to reaching the irreducible backgrounds are excellent, if the signal is just below the current limits, but slim if Generation 2 does not report WIMP detection.

We present a set of N-body simulations of a class of models in which an unstable dark matter particle decays into a stable non-interacting dark matter particle, with decay lifetime comparable to the Hubble time. We study the effects of the kinematic recoil velocity received by the stable dark matter on the structures of dark matter halos ranging from galaxy-cluster to Milky Way mass scales. For Milky Way-mass halos, we use high-resolution, zoom-in simulations to explore the effects of decays on Galactic substructure. In general, halos with circular velocities comparable to the magnitude of kick velocity are most strongly affected by decays. We show that decaying dark matter models with lifetimes comparable to Hubble time and recoil speeds about 20-40 km/s can significantly reduce both the abundance of Galactic subhalos and the internal densities of the subhalos. We also compare subhalo circular velocity profiles with observational constraints on the Milky Way dwarf satellite galaxies. Interestingly, we find that decaying dark matter models that do not violate current astrophysical constraints, can significantly mitigate both the well-documented "missing satellites problem" and the more recent "too big to fail problem" associated with the abundances and densities of Local Group dwarf satellite galaxies. A relatively unique feature of late decaying dark matter models is that they predict significant evolution of halos as a function of time. This is an important consideration because at high redshifts, prior to decays, decaying models exhibit the same sequence of structure formation as cold dark matter. We conclude that models of decaying dark matter make predictions that are relevant for the interpretation of observations of small galaxies in the Local Group and can be tested or constrained by the kinematics of Local Group dwarf galaxies as well as by forthcoming large-scale surveys.

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.

We present 3D kinematic observations of stars within the central 0.5 pc of the Milky Way nuclear star cluster using adaptive optics imaging and spectroscopy from the Keck telescopes. Recent observations have shown that the cluster has a shallower surface density profile than expected for a dynamically relaxed cusp, leading to important implications for its formation and evolution. However, the true three dimensional profile of the cluster is unknown due to the difficulty in de-projecting the stellar number counts. Here, we use spherical Jeans modeling of individual proper motions and radial velocities to constrain for the first time, the de-projected spatial density profile, cluster velocity anisotropy, black hole mass ($M_\mathrm{BH}$), and distance to the Galactic center ($R_0$) simultaneously. We find that the inner stellar density profile of the late-type stars, $\rho(r)\propto r^{-\gamma}$ to have a power law slope $\gamma=0.05_{-0.60}^{+0.29}$, much more shallow than the frequently assumed Bahcall $\&$ Wolf slope of $\gamma=7/4$. The measured slope will significantly affect dynamical predictions involving the cluster, such as the dynamical friction time scale. The cluster core must be larger than 0.5 pc, which disfavors some scenarios for its origin. Our measurement of $M_\mathrm{BH}=5.76_{-1.26}^{+1.76}\times10^6$ $M_\odot$ and $R_0=8.92_{-0.55}^{+0.58}$ kpc is consistent with that derived from stellar orbits within 1$^{\prime\prime}$ of Sgr A*. When combined with the orbit of S0-2, the uncertainty on $R_0$ is reduced by 30% ($8.46_{-0.38}^{+0.42}$ kpc). We suggest that the MW NSC can be used in the future in combination with stellar orbits to significantly improve constraints on $R_0$.

The search for weakly-interacting massive particle (WIMP) dark matter is multi-pronged. Ultimately, the WIMP-dark-matter picture will only be confirmed if different classes of experiments see consistent signals and infer the same WIMP properties. In this work, we review the ideas, methods, and status of direct-detection searches. We focus in particular on extracting WIMP physics (WIMP interactions and phase-space distribution) from direct-detection data in the early discovery days when multiple experiments see of order dozens to hundreds of events. To demonstrate the essential complementarity of different direct-detection experiments in this context, we create mock data intended to represent the data from the near-future Generation 2 experiments. We consider both conventional supersymmetry-inspired benchmark points (with spin-independent and -dependent elastic cross sections just below current limits), as well as benchmark points for other classes of models (inelastic and effective-operator paradigms). We also investigate the effect on parameter estimation of loosening or dropping the assumptions about the local WIMP phase-space distribution. We arrive at two main conclusions. Firstly, teasing out WIMP physics with experiments depends critically on having a wide set of detector target materials, spanning a large range of target nuclear masses and spin-dependent sensitivity. It is also highly desirable to obtain data from low-threshold experiments. Secondly, a general reconstruction of the local WIMP velocity distribution, which will only be achieved if there are multiple experiments using different target materials, is critical to obtaining a robust and unbiased estimate of the WIMP mass.

As part of the Snowmass process, the Cosmic Frontier Indirect-Detection subgroup (CF2) has drawn on input from the Cosmic Frontier and the broader Particle Physics community to produce this document. The purposes of this report are to identify opportunities for dark matter science through indirect detection, to give an overview of the primary scientific drivers for indirect searches for dark matter, and to survey current and planned experiments that have, as a large part of their scientific program, the goal of searching for indirect (or astrophysical) signatures of dark matter. We primarily address existing experiments with a large U.S. role, or future experiments where a U.S. contribution is sought. We also address the limitations of this technique, and answer the tough questions relevant to this subgroup posed by the HEP community through the Snowmass process.

We present an analysis of high-resolution N-body simulations of decaying dark matter cosmologies focusing on the statistical properties of the transmitted Lyman-alpha forest flux in the high-redshift intergalactic medium. In this type of model a dark matter particle decays into a slightly less massive stable dark matter daughter particle and a comparably light particle. The small mass splitting provides a non-relativistic kick velocity V_k to the daughter particle resulting in free-streaming and subsequent damping of small-scale density fluctuations. Current Lyman-alpha forest power spectrum measurements probe comoving scales up to ~ 2-3 h^-1 Mpc at redshifts z ~ 2-4, providing one of the most robust ways to probe cosmological density fluctuations on relatively small scales. The suppression of structure growth due to the free-streaming of dark matter daughter particles also has a significant impact on the neutral hydrogen cloud distribution, which traces the underlying dark matter distribution well at high redshift. We exploit Lyman-alpha forest power spectrum measurements to constrain the amount of free-streaming of dark matter in such models and thereby place limits on decaying dark matter based only on the dynamics of cosmological perturbations without any assumptions about the interactions of the decay products. We find that SDSS 1D Lyman-alpha forest power spectrum data place a lifetime-dependent upper limit V_k < 30-70 km/s for decay lifetimes < 10 Gyr. This is the most stringent model-independent bound on invisible dark matter decays with small mass splittings. For large mass splittings (large V_k), Lyman-alpha forest data restrict the dark matter lifetime to Gamma^-1 > 40 Gyr. Forthcoming BOSS data should be able to provide more stringent constraints on exotic dark matter, mainly because the larger BOSS quasar spectrum sample will significantly reduce statistical errors.

The cold dark matter (CDM) cosmological model has been remarkably successful in explaining cosmic structure over an enormous span of redshift, but it has faced persistent challenges from observations that probe the innermost regions of dark matter halos and the properties of the Milky Way's dwarf galaxy satellites. We review the current observational and theoretical status of these "small scale controversies." Cosmological simulations that incorporate only gravity and collisionless CDM predict halos with abundant substructure and central densities that are too high to match constraints from galaxy dynamics. The solution could lie in baryonic physics: recent numerical simulations and analytic models suggest that gravitational potential fluctuations tied to efficient supernova feedback can flatten the central cusps of halos in massive galaxies, and a combination of feedback and low star-formation efficiency could explain why most of the dark matter subhalos orbiting the Milky Way do not host visible galaxies. However, it is not clear that this solution can work in the lowest mass galaxies where discrepancies are observed. Alternatively, the small-scale conflicts could be evidence of more complex physics in the dark sector itself. For example, elastic scattering from strong dark matter self-interactions can alter predicted halo mass profiles, leading to good agreement with observations across a wide range of galaxy mass. Gravitational lensing and dynamical perturbations of tidal streams in the stellar halo provide evidence for an abundant population of low mass subhalos in accord with CDM predictions. These observational approaches will get more powerful over the next few years.

If dark matter has a large self-interaction scattering cross section, then interactions among dark-matter particles will drive galaxy and cluster halos to become spherical in their centers. Work in the past has used this effect to rule out velocity-independent, elastic cross sections larger than sigma/m ~ 0.02 cm^2/g based on comparisons to the shapes of galaxy cluster lensing potentials and X-ray isophotes. In this paper, we use cosmological simulations to show that these constraints were off by more than an order of magnitude because (a) they did not properly account for the fact that the observed ellipticity gets contributions from the triaxial mass distribution outside the core set by scatterings, (b) the scatter in axis ratios is large and (c) the core region retains more of its triaxial nature than estimated before. Including these effects properly shows that the same observations now allow dark matter self-interaction cross sections at least as large as sigma/m = 0.1 cm^2/g. We show that constraints on self-interacting dark matter from strong-lensing clusters are likely to improve significantly in the near future, but possibly more via central densities and core sizes than halo shapes.

We use cosmological simulations to study the effects of self-interacting dark matter (SIDM) on the density profiles and substructure counts of dark matter halos from the scales of spiral galaxies to galaxy clusters, focusing explicitly on models with cross sections over dark matter particle mass \sigma/m = 1 and 0.1 cm^2/g. Our simulations rely on a new SIDM N-body algorithm that is derived self-consistently from the Boltzmann equation and that reproduces analytic expectations in controlled numerical experiments. We find that well-resolved SIDM halos have constant-density cores, with significantly lower central densities than their CDM counterparts. In contrast, the subhalo content of SIDM halos is only modestly reduced compared to CDM, with the suppression greatest for large hosts and small halo-centric distances. Moreover, the large-scale clustering and halo circular velocity functions in SIDM are effectively identical to CDM, meaning that all of the large-scale successes of CDM are equally well matched by SIDM. From our largest cross section runs we are able to extract scaling relations for core sizes and central densities over a range of halo sizes and find a strong correlation between the core radius of an SIDM halo and the NFW scale radius of its CDM counterpart. We construct a simple analytic model, based on CDM scaling relations, that captures all aspects of the scaling relations for SIDM halos. Our results show that halo core densities in \sigma/m = 1 cm^2/g models are too low to match observations of galaxy clusters, low surface brightness spirals (LSBs), and dwarf spheroidal galaxies. However, SIDM with \sigma/m ~ 0.1 cm^2/g appears capable of reproducing reported core sizes and central densities of dwarfs, LSBs, and galaxy clusters without the need for velocity dependence. (abridged)

We explore the ability of directional nuclear-recoil detectors to constrain the local velocity distribution of weakly interacting massive particle (WIMP) dark matter by performing Bayesian parameter estimation on simulated recoil-event data sets. We discuss in detail how directional information, when combined with measurements of the recoil-energy spectrum, helps break degeneracies in the velocity-distribution parameters. We also consider the possibility that velocity structures such as cold tidal streams or a dark disk may also be present in addition to the Galactic halo. Assuming a carbon-tetrafluoride detector with a 30-kg-yr exposure, a 50-GeV WIMP mass, and a WIMP-nucleon spin-dependent cross-section of 0.001 pb, we show that the properties of a cold tidal stream may be well constrained. However, measurement of the parameters of a dark-disk component with a low lag speed of ~50 km/s may be challenging unless energy thresholds are improved.

We discuss irreducible statistical limitations of future ton-scale dark matter direct detection experiments. We focus in particular on the coverage of confidence intervals, which quantifies the reliability of the statistical method used to reconstruct the dark matter parameters, and the bias of the reconstructed parameters. We study 36 benchmark dark matter models within the reach of upcoming ton-scale experiments. We find that approximate confidence intervals from a profile-likelihood analysis exactly cover or over-cover the true values of the WIMP parameters, and are hence conservative. We evaluate the probability that unavoidable statistical fluctuations in the data might lead to a biased reconstruction of the dark matter parameters, or large uncertainties on the reconstructed parameter values. We show that this probability can be surprisingly large, even for benchmark models leading to a large event rate of order a hundred counts. We find that combining data sets from two different targets leads to improved coverage properties, as well as a substantial reduction of statistical bias and uncertainty on the dark matter parameters.

From astronomical observations, we know that dark matter exists, makes up 23% of the mass budget of the Universe, clusters strongly to form the load-bearing frame of structure for galaxy formation, and hardly interacts with ordinary matter except gravitationally. However, this information is not enough to identify the particle specie(s) that make up dark matter. As such, the problem of determining the identity of dark matter has largely shifted to the fields of astroparticle and particle physics. In this talk, I will review the current status of the search for the nature of dark matter. I will provide an introduction to possible particle candidates for dark matter and highlight recent experimental astroparticle- and particle-physics results that constrain the properties of those candidates. Given the absence of detections in those experiments, I will advocate a return of the problem of dark-matter identification to astronomy, and show what kinds of theoretical and observational work might be used to pin down the nature of dark matter once and for all. This talk is intended for a broad astronomy audience.

Once weakly-interacting massive particles (WIMPs) are unambiguously detected in direct-detection experiments, the challenge will be to determine what one may infer from the data. Here, I examine the prospects for reconstructing the local speed distribution of WIMPs in addition to WIMP particle-physics properties (mass, cross sections) from next-generation cryogenic and liquid-noble direct-detection experiments. I find that the common method of fixing the form of the velocity distribution when estimating constraints on WIMP mass and cross sections means losing out on the information on the speed distribution contained in the data and may lead to biases in the inferred values of the particle-physics parameters. I show that using a more general, empirical form of the speed distribution can lead to good constraints on the speed distribution. Moreover, one can use Bayesian model-selection criteria to determine if a theoretically-inspired functional form for the speed distribution (such as a Maxwell-Boltzmann distribution) fits better than an empirical model. The shape of the degeneracy between WIMP mass and cross sections and their offset from the true values of those parameters depends on the hypothesis for the speed distribution, which has significant implications for consistency checks between direct-detection and collider data. In addition, I find that the uncertainties on theoretical parameters depends sensitively on the upper end of the energy range used for WIMP searches. Better constraints on the WIMP particle-physics parameters and speed distribution are obtained if the WIMP search is extended to higher energy (~ 1 MeV).

We analyze subhalos in the Via Lactea II (VL2) cosmological simulation to look for correlations among their infall times and z = 0 dynamical properties. We find that the present day orbital energy is tightly correlated with the time at which subhalos last crossed into the virial radius. This energy-infall correlation provides a means to infer infall times for Milky Way satellite galaxies. Assuming that the Milky Way's assembly can be modeled by VL2, we show that the infall times of some satellites are well constrained given only their Galactocentric positions and line-of-sight velocities. The constraints sharpen for satellites with proper motion measurements. We find that Carina, Ursa Minor, and Sculptor were all accreted early, more than 8 Gyr ago. Five other dwarfs, including Sextans and Segue 1, are also probable early accreters, though with larger uncertainties. On the other extreme, Leo T is just falling into the Milky Way for the first time while Leo I fell in \sim 2 Gyr ago and is now climbing out of the Milky Way's potential after its first perigalacticon. The energies of several other dwarfs, including Fornax and Hercules, point to intermediate infall times, 2 - 8 Gyr ago. We compare our infall time estimates to published star formation histories and find hints of a dichotomy between ultrafaint and classical dwarfs. The classical dwarfs appear to have quenched star formation after infall but the ultrafaint dwarfs tend to be quenched long before infall, at least for the cases in which our uncertainties allow us to discern differences. Our analysis suggests that the Large Magellanic Cloud crossed inside the Milky Way virial radius recently, within the last \sim 4 billion years.

We have analyzed high resolution N-body simulations of dark matter halos, focusing specifically on the evolution of angular momentum. We find that not only is individual particle angular momentum not conserved, but the angular momentum of radial shells also varies over the age of the Universe by up to factors of a few. We find that torques from external structure are the most likely cause for this distribution shift. Since the model of adiabatic contraction that is often applied to model the effects of galaxy evolution on the dark-matter density profile in a halo assumes angular momentum conservation, this variation implies that there is a fundamental limit on the possible accuracy of the adiabatic contraction model in modeling the response of DM halos to the growth of galaxies.

We explore a dark-matter model in which there are two dark-matter species nearly degenerate in mass, with epsilon = Delta M/M << 1. The heavier particle undergoes two-body decay with a half-life tau, to the lighter dark-matter particle and a noninteracting massless particle. Unlike previous work on decaying dark matter, we explore the regime tau > 100 Myr and non-relativistic kick speeds vk / c = epsilon. Using a set of N-body simulations of isolated dark-matter halos, we show how halos change as a function of tau and vk. We find that tau < 40 Gyr is ruled out for vk > 20 km s^{-1} (epsilon > 10^{-4}) when we compare the simulations to observations of dwarf-galaxy- to cluster-mass dark matter halos. We highlight which set of observations should provide better future constraints for decays and other types of dark-matter physics.

We consider constraints on a phenomenological dark-matter model consisting of two nearly degenerate particle species using observed properties of the Milky Way satellite galaxy population. The two parameters of this model, assuming the particle masses are >~ GeV, are v_k, the recoil speed of the daughter particle, and tau, the lifetime of the parent particle. The satellite constraint that spans the widest range of v_k is the number of satellites that have a mass within 300 pc M300 > 5 x 10^6 solar masses, although constraints based on M300 in the classical dwarfs and the overall velocity function are competitive for v_k >~ 50 km/s. In general, we find that tau <~ 30 Gyr is ruled out for 20 km/s <~ v_k <~ 200 km/s, although we find that the limits on tau for fixed v_k can change constraints by a factor of ~3 depending on the star-formation histories of the satellites. We advocate using the distribution of M300 in Milky Way satellites determined by next-generation all-sky surveys and follow-up spectroscopy as a probe of dark-matter properties.

We consider models in which a dark-matter particle decays to a slightly less massive daughter particle and a noninteracting massless particle. The decay gives the daughter particle a small velocity kick. Self-gravitating dark-matter halos that have a virial velocity smaller than this velocity kick may be disrupted by these particle decays, while those with larger virial velocities will be heated. We use numerical simulations to follow the detailed evolution of the total mass and density profile of self-gravitating systems composed of particles that undergo such velocity kicks as a function of the kick speed (relative to the virial velocity) and the decay time (relative to the dynamical time). We show how these decays will affect the halo mass-concentration relation and mass function. Using measurements of the halo mass-concentration relation and galaxy-cluster mass function to constrain the lifetime--kick-velocity parameter space for decaying dark matter, we find roughly that the observations rule out the combination of kick velocities greater than 100 km/s and decay times less than a few times the age of the Universe.

Recently, several papers have appeared that examine the process of capturing dark-matter particles from the Galactic halo to orbits bound to the Solar System. The authors of these papers predict large enhancements to the local dark-matter density via gravitational three-body interactions with planets. However, these conclusions are wrong; these papers do not include the inverse process to capture, namely the ejection of dark-matter particles by three-body gravitational encounters. We emphasize previous work that shows that by including both capture and ejection of dark matter from the Solar System, the density of dark matter bound to the Solar System is small compared to the local Galactic dark-matter density.

I consider constraints on a phenomenological decaying-dark-matter model, in which two weakly-interacting massive particle (WIMP) species have a small mass splitting, and in which the heavier particle decays to the lighter particle and a massless particle on cosmological timescales. The decay parameter space is parameterized by $v_k$, the speed of the lighter particle in the center-of-mass frame of the heavier particle prior to decay, and the decay time $\tau$. Since I consider the case in which dark-matter halos have formed before there has been significant decay, I focus on the effects of decay in already-formed halos. I show that the $v_k-\tau$ parameter space may be constrained by observed properties of dark-matter halos. I highlight which set of observations is likely to yield the cleanest constraints on $v_k-\tau$ parameter space, and calculate the constraints in those cases in which the effect of decay on the observables can be calculated without N-body simulations of decaying dark matter. I show that for $v_k \gtrsim 5\times 10^3$ km s$^{-1}$, the z=0 galaxy cluster mass function and halo mass-concentration relation constrain $\tau \gtrsim$ 40 Gyr, and that precise constraints on $\tau$ for smaller $v_k$ will require N-body simulations.

The next decade will bring massive new data sets from experiments of the direct detection of weakly interacting massive particle (WIMP) dark matter. The primary goal of these experiments is to identify and characterize the dark-matter particle species. However, mapping the data sets to the particle-physics properties of dark matter is complicated not only by the considerable uncertainties in the dark-matter model, but by its poorly constrained local distribution function (the "astrophysics" of dark matter). In this Letter, I propose a shift in how to do direct-detection data analysis. I show that by treating the astrophysical and particle physics uncertainties of dark matter on equal footing, and by incorporating a combination of data sets into the analysis, one may recover both the particle physics and astrophysics of dark matter. Not only does such an approach yield more accurate estimates of dark-matter properties, but may illuminate how dark matter coevolves with galaxies.

One method to search for particle dark matter is to hunt down its annihilation products. In the Solar System, three potential types of signals of annihilation have been identified: neutrinos and gamma-rays from the Sun, and neutrinos from the Earth. Each of these signals depends sensitively on the orbital evolution of dark matter once it becomes bound to the Solar System. I will review progress on characterizing these signals based on recent improvements in the determination of the properties of the bound dark matter population.

As disc galaxies form in a hierarchical cosmology, massive merging satellites are preferentially dragged towards the disc plane. The material accreted from these satellites forms a dark matter disc that contributes 0.25 - 1.5 times the non-rotating halo density at the solar position. Here, we show the importance of the dark disc for indirect dark matter detection in neutrino telescopes. Previous predictions of the neutrino flux from WIMP annihilation in the Earth and the Sun have assumed that Galactic dark matter is spherically distributed with a Gaussian velocity distribution, the standard halo model. Although the dark disc has a local density comparable to the dark halo, its higher phase space density at low velocities greatly enhances capture rates in the Sun and Earth. For typical dark disc properties, the resulting muon flux from the Earth is increased by three orders of magnitude over the SHM, while for the Sun the increase is an order of magnitude. This significantly increases the sensitivity of neutrino telescopes to fix or constrain parameters in WIMP models. The flux from the Earth is extremely sensitive to the detailed properties of the dark disc, while the flux from the Sun is more robust. The enhancement of the muon flux from the dark disc puts the search for WIMP annihilation in the Earth on the same level as the Sun for WIMP masses < 100 GeV.

In this last paper in a series of three on weakly interacting massive particle (WIMP) dark matter in the solar system, we focus on WIMPs bound to the system by gravitationally scattering off of planets. We present simulations of WIMP orbits in a toy solar system consisting of only the Sun and Jupiter. As previous work suggested, we find that the density of gravitationally captured WIMPs at the Earth is small and largely insensitive to the details of elastic scattering in the Sun. However, we find that the density of gravitationally captured WIMPs may be affected by external Galactic gravitational fields. If such fields are unimportant, the density of gravitationally captured WIMPs at the Earth should be similar to the maximum density of WIMPs captured in the solar system by elastic scattering in the Sun. Using standard assumptions about the halo WIMP distribution function, we find that the gravitationally captured WIMPs contribute negligibly to direct detection event rates. While these WIMPs do dominate the annihilation rate of WIMPs in the Earth, the resulting event rate in neutrino telescopes is too low to be observed in next-generation neutrino telescopes.

We calculate the annihilation rate of weakly interacting massive particles (WIMPs) in the Sun as a function of their mass and elastic scattering cross section. One byproduct of the annihilation, muon neutrinos, may be observed by the next generation of neutrino telescopes. Previous estimates of the annihilation rate assumed that any WIMPs from the Galactic dark halo that are captured in the Sun by elastic scattering off solar nuclei quickly reach thermal equilibrium in the Sun. We show that the optical depth of the Sun to WIMPs and the gravitational forces from planets both serve to decrease the annihilation rate below these estimates. While we find that the sensitivity of upcoming km^3-scale neutrino telescopes to ~100 GeV WIMPs is virtually unchanged from previous estimates, the sensitivity of these experiments to ~10 TeV WIMPs may be an order of magnitude less than the standard calculations would suggest. The new estimates of the annihilation rates should guide future experiment design and improve the mapping from neutrino event rates to WIMP parameter space.

The next generation of dark matter (DM) direct detection experiments and neutrino telescopes will probe large swaths of dark matter parameter space. In order to interpret the signals in these experiments, it is necessary to have good models of both the halo DM streaming through the solar system and the population of DM bound to the solar system. In this paper, the first in a series of three on DM in the solar system, we present simulations of orbits of DM bound to the solar system by solar capture in a toy solar system consisting of only the Sun and Jupiter, assuming that DM consists of a single species of weakly interacting massive particle (WIMP). We describe how the size of the bound WIMP population depends on the WIMP mass, spin-independent cross section, and spin-dependent cross section. Using a standard description of the Galactic DM halo, we find that the maximum enhancement to the direct detection event rate, consistent with current experimental constraints on the WIMP-nucleon cross section, is < 1% relative to the event rate from halo WIMPs, while the event rate from neutrinos from WIMP annihilation in the center of the Earth is unlikely to meet the threshold of next-generation, km^3-sized (IceCube, KM3NeT) neutrino telescopes.

Semi-analytic treatments of the evolution of orbits of weakly interacting massive particles (WIMPs) in the solar system suggest that the WIMPs bound to the solar system may enhance the direct detection rate relative to that of the unbound population by up to a factor of order unity, and boost the flux of neutrinos from WIMP annihilation in the Earth by up to two orders of magnitude. To test these important but uncertain results, we perform a suite of numerical orbit integrations to explore the properties of the bound WIMP population as a function of the WIMP mass and the scattering cross section with baryonic matter. For regions of WIMP parameter space presently allowed by experiments, we find that (i) the bound WIMP population enhances the direct detection rate by at most ~1% relative to the rate from unbound halo WIMPs; (ii) it is unlikely that planned km^3-scale neutrino telescopes will detect neutrinos from WIMP annihilation in the Earth; (iii) the event rate from neutrinos produced by WIMP annihilation in the Sun may be much smaller than implied by the usual calculations, which assume that WIMPs scattered onto bound orbits are rapidly thermalized in the Sun.

We present results from the first robust investigation of galaxy morphology as a function of environment at z>1.5. Our study is motivated by the fact that star-forming galaxies contained within a protocluster at z=2.3 in the HS1700+64 field have significantly older ages and larger stellar masses on average than those at similar redshifts but more typical environmental densities. In the analysis of HST/ACS images, we apply non-parametric statistics to characterize the rest-frame UV morphologies of a sample of 85 UV-selected star-forming galaxies at z=1.7-2.9, 22 of which are contained in the protocluster. The remaining 63 control-sample galaxies are not in the protocluster but have a similar mean redshift of <z>~2.3. We find no environmental dependence for the distributions of morphological properties. Combining the measured morphologies with the results of population synthesis modeling, we find only weak correlations, if any, between morphological and stellar population properties such as stellar mass, age, extinction and star-formation rate. Given the incomplete census of the protocluster galaxy population, and the lack of correlation between rest-frame UV morphology and star-formation history at z~2 within our sample, the absence of environmental trends in the distribution of morphological properties is not surprising. Additionally, using a larger sample of photometric candidates, we compare morphological distributions for 282 UV-selected and 43 near-IR-selected galaxies. While the difference in the degree of nebulosity between the two samples appears to be a byproduct of the fainter average rest-frame UV surface brightness of the near-IR-selected galaxies, we find that, among the lowest surface brightness galaxies, the near-IR-selected objects have significantly smaller angular sizes (abridged).