We present deep $g$- and $r$-band Magellan/Megacam photometry of two dwarf galaxy candidates discovered in the Dark Energy Survey (DES), Grus I and Indus II (DES J2038-4609). For the case of Grus I, we resolved the main sequence turn-off (MSTO) and $\sim 2$ mags below it. The MSTO can be seen at $g_0\sim 24$ with a photometric uncertainty of $0.03$ mag. We show Grus I to be consistent with an old, metal-poor ($\sim 13.3$ Gyr, [Fe/H]$\sim-1.9$) dwarf galaxy. We derive updated distance and structural parameters for Grus I using this deep, uniform, wide-field data set. We find an azimuthally averaged half-light radius more than two times larger ($\sim 151^{+21}_{-31}$ pc; $\sim 4.^{\prime} 16^{+0.54}_{-0.74}$) and an absolute $V$-band magnitude $\sim-4.1$ that is $\sim 1$ magnitude brighter than previous studies. We obtain updated distance, ellipticity, and centroid parameters which are in agreement with other studies within uncertainties. Although our photometry of Indus II is $\sim 2-3$ magnitudes deeper than the DES Y1 Public release, we find no coherent stellar population at its reported location. The original detection was located in an incomplete region of sky in the DES Y2Q1 data set and was flagged due to potential blue horizontal branch member stars. The best fit isochrone parameters are physically inconsistent with both dwarf galaxies and globular clusters. We conclude that Indus II is likely a false-positive, flagged due to a chance alignment of stars along the line of sight.

There is increasing evidence that a substantial fraction of Milky Way satellite galaxies align in a rotationally-supported plane of satellites, a rare configuration in cosmological simulations of galaxy formation. It has been suggested that other Milky Way substructures (namely young halo globular clusters and stellar/gaseous streams) similarly tend to align with this plane, accordingly dubbed the Vast Polar Structure (VPOS). Using systemic proper motions inferred from Gaia data, we find that globular cluster orbital poles are not clustered in the VPOS direction, though the population with the highest VPOS membership fraction is the young halo clusters (~30%). We additionally provide a current census of stellar streams, including new streams discovered using the Dark Energy Survey and Gaia datasets, and find that stellar stream normals are also not clustered in the direction of the VPOS normal. We also find that, based on orbit modeling, there is a likely association between NGC 3201 and the Gj\"{o}ll stellar stream and that, based on its orbital pole, NGC 4147 is likely not a Sagittarius globular cluster. That the Milky Way's accreted globular clusters and streams do not align in the same planar configuration as its satellites suggests that the plane of satellites is either a particularly stable orbital configuration or a population of recently accreted satellites. Neither of these explanations is particularly likely in light of other recent studies, leaving the plane of satellites problem as one of the more consequential open problems in galaxy formation and cosmology.

We study the sensitivity of future xenon- and argon-based dark matter and neutrino detection experiments to low-energy atmospheric neutrinos. Not accounting for experimental backgrounds, the primary obstacle for identifying nuclear recoils induced by atmospheric neutrinos is the tail of the electron recoil distribution due to $pp$ solar neutrinos. We use the NEST code to model the solar and atmospheric neutrino signals in a xenon detector and find that an exposure of 700 tonne-years will produce a $5\sigma$ detection of atmospheric neutrinos. We explore the effect of different detector properties and find that a sufficiently long electron lifetime is essential to the success of such a measurement.

The Magnificent CE$\nu$NS Workshop (2018) was held November 2 & 3 of 2018 on the University of Chicago campus and brought together theorists, phenomenologists, and experimentalists working in numerous areas but sharing a common interest in the process of coherent elastic neutrino-nucleus scattering (CE$\nu$NS). This is a collection of abstract-like summaries of the talks given at the meeting, including links to the slides presented. This document and the slides from the meeting provide an overview of the field and a snapshot of the robust CE$\nu$NS-related efforts both planned and underway.

We calculate the effective $J$-factors, which determine the strength of indirect detection signals from dark matter annihilation, for 25 dwarf spheroidal galaxies (dSphs). We consider several well-motivated assumptions for the relative velocity dependence of the dark matter annihilation cross section: $\sigma_A v$: $s$-wave (velocity independent), $p$-wave ($\sigma_A v \propto v^2$), $d$-wave ($\sigma_A v \propto v^4$), and Sommerfeld-enhancement in the Coulomb limit ($\sigma_A v \propto 1/v$). As a result we provide the largest and most updated sample of J-factors for velocity-dependent annihilation models. For each scenario, we use Fermi-LAT gamma-ray data to constrain the annihilation cross section. Due to the assumptions made in our gamma-ray data analysis, our bounds are comparable to previous bounds on both the $p$-wave and Sommerfeld-enhanced cross sections using dSphs. Our bounds on the $d$-wave cross section are the first such bounds using indirect detection data.

We present the effective $J$-factors for the Milky Way for scenarios in which dark matter annihilation is p-wave or d-wave suppressed. We find that the velocity suppression of dark matter annihilation can have a sizable effect on the morphology of a potential dark matter annihilation signal in the Galactic Center. The gamma-ray flux from the innermost region of the Galactic Center is in particular suppressed. We find that for dark matter density profiles with steep inner slopes, the morphology of the Inner Galaxy gamma-ray emission in p-wave models can be made similar to the morphology in standard s-wave models. This similarity may suggest that model discrimination between s-wave and p-wave is challenging, for example, when fitting the Galactic Center excess. However, we show that it is difficult to simultaneously match s- and p-wave morphologies at both large and small angular scales. The $J$-factors we calculate may be implemented with astrophysical foreground models to self-consistently determine the morphology of the excess with velocity-suppressed dark matter annihilation.

We consider the effect of velocity-dependent dark matter annihilation on the angular distribution of gamma rays produced in dark matter subhalos. We assume that the dark matter potential is spherically symmetric, characterized by a scale radius and scale density, and the velocity distribution is isotropic. We find that the effect of velocity-dependent dark matter annihilation is largely determined by dimensional analysis; the angular size of gamma-ray emission from an individual subhalo is rescaled by a factor which depends on the form of the dark matter distribution, but not on the halo parameters, while the relative normalization of the gamma-ray flux from different mass subhalos is rescaled by a factor which depends on the halo parameters, but not on the form of the dark matter distribution. We apply our results to a Navarro-Frenk-White profile for the case of an individual subhalo and comment on the application of these results to a distribution of subhalos.

We obtain distances to a sample of RR Lyrae in the central core of the Sagittarius dwarf spheroidal galaxy from OGLE data. We use these distances, along with RR Lyrae from \emph{Gaia} DR2, to measure the shape of the stellar distribution within the central $\sim$ 2 kpc. The best-fit stellar distribution is triaxial, with axis ratios 1 : 0.76 : 0.43. A prolate spheroid model is ruled out at high statistical significance relative to the triaxial model. The major axis is aligned nearly parallel to the sky plane as seen by an Earth-based observer and is nearly perpendicular to the direction of the Galactic center. This result may be compared to cosmological simulations which generally predict that the major axis of the dark matter distribution of subhalos is aligned with the Galactic center. The triaxial structure that we obtain can provide important constraints on the Sagittarius progenitor, as well as the central dark matter distribution under the assumption of dynamical equilibrium.

Indirect detection strategies of particle Dark Matter (DM) in Dwarf spheroidal satellite galaxies (dSphs) typically entail searching for annihilation signals above the astrophysical background. To robustly compare model predictions with the observed fluxes of product particles, most analyses of astrophysical data -- which are generally frequentist -- rely on estimating the abundance of DM by calculating the so-called $\textit{J-factor}$. This quantity is usually inferred from the kinematic properties of the stellar population of a dSph using Jeans equation, commonly by means of Bayesian techniques which entail the presence (and additional systematic uncertainty) of prior choice. Here, extending earlier work, we develop a scheme to derive the profile likelihood for $J$-factors of dwarf spheroidals for models with five or more free parameters. We validate our method on a publicly available simulation suite, released by the Gaia Challenge, finding satisfactory statistical properties for coverage and bias. We present the profile likelihood function and maximum likelihood estimates for the $J$-factor of ten dSphs . As an illustration, we apply these profiles likelihood to recently published analyses of gamma-ray data with the Fermi Large Area Telescope to derive new, consistent upper limits on the DM annihilation cross-section. We do this for a subset of systems, generally referred to as $\textit{classical dwarfs}$. The implications of these findings for DM searches are discussed, together with future improvements and extensions of this technique.

We analyse the orbital kinematics of the Milky Way (MW) satellite system utilizing the latest systemic proper motions for 38 satellites based on data from Gaia Data Release 2. Combining these data with distance and line-of-sight velocity measurements from the literature, we use a likelihood method to model the velocity anisotropy, $\beta$, as a function of Galactocentric distance and compare the MW satellite system with those of simulated MW-mass haloes from the APOSTLE and Auriga simulation suites. The anisotropy profile for the MW satellite system increases from $\beta\sim -2$ at $r\sim20$ kpc to $\beta\sim 0.5$ at $r\sim200$ kpc, indicating that satellites closer to the Galactic centre have tangentially-biased motions while those farther out have radially-biased motions. The motions of satellites around APOSTLE host galaxies are nearly isotropic at all radii, while the $\beta(r)$ profiles for satellite systems in the Auriga suite, whose host galaxies are substantially more massive in baryons than those in APOSTLE, are more consistent with that of the MW satellite system. This shape of the $\beta(r)$ profile may be attributed to the central stellar disc preferentially destroying satellites on radial orbits, or intrinsic processes from the formation of the Milky Way system.

We use the Fermi-LAT gamma-ray observatory to search for gamma-ray emission from four nearby, debris disk-hosting main sequence stars: $\tau$ Ceti, $\epsilon$ Eridani, Fomalhaut, and Vega. For three stars ($\tau$ Ceti, Fomalhaut, and Vega), we establish upper limits that are consistent with theoretical expectations. For $\epsilon$ Eridani, we find a possible spatially coincident source with a soft energy spectrum of $dN/dE \sim E^{-3.6}$. However, at this stage we are unable to rule out that this emission is due to a more extended feature in the diffuse background. In the interpretation that the emission is due to $\epsilon$ Eridani, the $> 100$ MeV gamma-ray luminosity is $\sim 10^{27}$ erg/s $\simeq 3\times 10^{-7}$ L$_\odot$, which is $\sim 10^{10}$ times the gamma-ray luminosity from the disk of the quiet Sun. We find $\lesssim 2 \sigma$ evidence of source variability over a $\sim 7$ year timescale. In the interpretation that the gamma-ray emission from $\epsilon$ Eridani itself, we consider two possible models: 1) cosmic-ray collisions with solid bodies in the debris disk which extends out $\sim$60 AU from the host star, and 2) emission from the stellar activity. For the former model, assuming a total disk mass consistent with infrared measurements, we find that the size distribution of bodies is steeper than expected for a collisional cascade state. If confirmed as being associated with $\epsilon$ Eridani, this would be the first indication of gamma-ray emission from the vicinity of a main sequence star other than the Sun.

Dark matter (DM) is the most abundant material in the Universe, but has so far been detected only via its gravitational effects. Several theories suggest that pairs of DM particles can annihilate into a flash of light at gamma-ray wavelengths. While gamma-ray emission has been observed from environments where DM is expected to accumulate, such as the centre of our Galaxy, other high energy sources can create a contaminating astrophysical gamma-ray background, thus making DM detection difficult. In principle, dwarf galaxies around the Milky Way are a better place to look -- they contain a greater fraction of DM with no astrophysical gamma-ray background -- but they are too distant for gamma-rays to have been seen. A range of observational evidence suggests that Omega Centauri (omega Cen or NGC 5139), usually classified as the Milky Way's largest globular cluster, is really the core of a captured and stripped dwarf galaxy. Importantly, Omega Cen is ten times closer to us than known dwarfs. Here we show that not only does Omega Cen contain DM with density as high as compact dwarf galaxies, but also that it emits gamma-rays with an energy spectrum matching that expected from the annihilation of DM particles with mass 31$\pm$4 GeV (68\% confidence limit). No astrophysical sources have been found that would otherwise explain Omega Cen's gamma-ray emission, despite deep multi-wavelength searches. We anticipate our results to be the starting point for even deeper radio observations of Omega Cen. If multi-wavelength searches continue to find no astrophysical explanations, this pristine, nearby clump of DM will become the best place to study DM interactions through forces other than gravity.

We present a chemical abundance analysis of four additional confirmed member stars of Tucana III, a Milky Way satellite galaxy candidate in the process of being tidally disrupted as it is accreted by the Galaxy. Two of these stars are centrally located in the core of the galaxy while the other two stars are located in the eastern and western tidal tails. The four stars have chemical abundance patterns consistent with the one previously studied star in Tucana III: they are moderately enhanced in $r$-process elements, i.e. they have $<$[Eu/Fe]$> \approx +$0.4 dex. The non-neutron-capture elements generally follow trends seen in other dwarf galaxies, including a metallicity range of 0.44 dex and the expected trend in $\alpha$-elements, i.e., the lower metallicity stars have higher Ca and Ti abundance. Overall, the chemical abundance patterns of these stars suggest that Tucana III was an ultra-faint dwarf galaxy, and not a globular cluster, before being tidally disturbed. As is the case for the one other galaxy dominated by $r$-process enhanced stars, Reticulum II, Tucana III's stellar chemical abundances are consistent with pollution from ejecta produced by a binary neutron star merger, although a different $r$-process element or dilution gas mass is required to explain the abundances in these two galaxies if a neutron star merger is the sole source of $r$-process enhancement.

Many studies have shown that RR Lyrae variable stars (RRL) are powerful stellar tracers of Galactic halo structure and satellite galaxies. The Dark Energy Survey (DES), with its deep and wide coverage (g ~ 23.5 mag) in a single exposure; over 5000 deg$^{2}$) provides a rich opportunity to search for substructures out to the edge of the Milky Way halo. However, the sparse and unevenly sampled multiband light curves from the DES wide-field survey (median 4 observations in each of grizY over the first three years) pose a challenge for traditional techniques used to detect RRL. We present an empirically motivated and computationally efficient template fitting method to identify these variable stars using three years of DES data. When tested on DES light curves of previously classified objects in SDSS stripe 82, our algorithm recovers 89% of RRL periods to within 1% of their true value with 85% purity and 76% completeness. Using this method, we identify 5783 RRL candidates, ~31% of which are previously undiscovered. This method will be useful for identifying RRL in other sparse multiband data sets.

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

For nearly 40 years, dark matter has been widely assumed to be cold and collisionless. Cold dark matter models make fundamental predictions for the behavior of dark matter on small (<10 kpc) scales. These predictions include cuspy density profiles at the centers of dark matter halos and a halo mass function that increases as dN/dM ~ M^-1.9 down to very small masses. We suggest two observational programs relying on extremely large telescopes to critically test these predictions, and thus shed new light on the nature of dark matter. (1) Combining adaptive optics-enabled imaging with deep spectroscopy to measure the three-dimensional motions of stars within a sample of Local Group dwarf galaxies that are the cleanest dark matter laboratories known in the nearby universe. From these observations the inner slope of the dark matter density profile can be determined with an accuracy of 0.20 dex, enabling a central cusp to be distinguished from a core at 5 sigma significance. (2) Diffraction-limited AO imaging and integral field spectroscopy of gravitationally lensed galaxies and quasars to quantify the abundance of dark substructures in the halos of the lens galaxies and along the line of sight. Observations of 50 lensed arcs and 50 multiply-imaged quasars will be sufficient to measure the halo mass function over the range 10^7 < M < 10^10 Msun at cosmological scales, independent of the baryonic and stellar composition of those structures. These two observational probes provide complementary information about the small scale structure, with a joint self-consistent analysis mitigating limitations of either probe. This program will produce the strongest existing constraints on the properties of dark matter on small scales, allowing conclusive tests of alternative warm, fuzzy, and self-interacting dark matter models.

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 study the prospects for measuring the low-energy components of the solar neutrino flux in future direct dark matter detection experiments. We show that for a depletion of $^{136}$Xe by a factor of 100 relative to its natural abundance, and an extension to electron recoil energies of $\sim$ MeV, future xenon experiments with exposure $\sim 200$ ton-yr can detect the CNO component of the solar neutrino flux at $\sim 3 \sigma$ significance. A CNO detection will provide important insight into metallicity of the solar interior. Precise measurement of low-energy solar neutrinos, including as $pp$, $^7$Be, and $pep$ components, will further improve constraints on the "neutrino luminosity" of the Sun, thereby providing constraints on alternative sources of energy production. We find that a measurement of $L_{\nu}/L_{\odot}$ of order one percent is possible with the above exposure, improving on current bounds from a global analysis of solar neutrino data by a factor of about seven.

Direct dark matter detection experiments will soon be sensitive to neutrinos from astrophysical sources, including the Sun, the atmosphere, and supernova. This sets an important benchmark for these experiments, and opens up a new window in neutrino physics and astrophysics. The detection of these neutrinos will be complementary to accelerator and reactor-based experiments which study neutrinos over the same energy range. Here we review the physics and astrophysics that can be extracted from the detection of these neutrinos, highlighting the potential for identifying new physics in the form of light mediators that arise from kinetic mixing and hidden sectors, and $\sim$ eV-scale sterile neutrinos. We discuss how the physics reach of these experiments will complement searches for new physics at the LHC and dedicated neutrino experiments.

For many dark matter models, the annihilation cross section to two-body final states is difficult to probe with current experiments because the dominant annihilation channel is velocity or helicity suppressed. The inclusion of gauge boson radiation for three-body final states can lift the helicity suppression, allowing a velocity-independent cross section to dominate the annihilation process, and providing an avenue to constrain these models. Here we examine experimental constraints on dark matter that annihilates to two leptons plus a bremsstrahlung boson, $\bar{\ell}+\ell+\gamma/W/Z$. We consider experimental constraints on photon final states from Fermi-LAT using both diffuse photon data and data from dwarf spheroidal galaxies, and compare to the implied constraints from 21 cm measurements. Diffuse photon line searches are generally the strongest over the entire mass regime. We in particular highlight the model in which dark matter annihilates to two neutrinos and a photon, and show that these models are more strongly constrained through photon measurements than through existing neutrino bounds.

We apply a Bayesian method to model multi-epoch radial velocity measurements in the ultra-faint dwarf galaxy Reticulum II, fully accounting for the effects of binary orbital motion and systematic offsets between different spectroscopic datasets. We find that the binary fraction of Ret II is higher than 0.5 at the 90% confidence level, if the mean orbital period is assumed to be 30 years or longer. Despite this high binary fraction, we infer a best-fit intrinsic dispersion of 2.8$_{-1.2}^{+0.7}$ km/s, which is smaller than previous estimates, but still indicates Ret II is a dark-matter dominated galaxy. We likewise infer a $\lesssim$ 1% probability that Ret II's dispersion is due to binaries rather than dark matter, corresponding to the regime $M_{\odot}/L_{\odot} \lesssim$ 2. Our inference of a high close binary fraction in Ret II echoes previous results for the Segue 1 ultra-faint dwarf and is consistent with studies of Milky Way halo stars that indicate a high close binary fraction tends to exist in metal-poor environments.

Measuring the dark matter distribution in dwarf spheroidal galaxies (dSphs) from stellar kinematics is crucial for indirect dark matter searches, as these distributions set the fluxes for both dark matter annihilation (J-Factor) and decay (D-Factor). Here we produce a compilation of J and D-Factors for dSphs, including new calculations for several newly-discovered Milky Way (MW) satellites, for dSphs outside of the MW virial radius, and for M31 satellites. From this compilation we test for scaling relations between the J and D-factors and physical properties of the dSphs such as the velocity dispersion ($\sigma_{\mathrm{los}}$), the distance ($d$), and the stellar half-light radius ($r_{1/2}$). We find that the following scaling relation minimizes the residuals as compared to different functional dependencies on the observed dSphs properties $J(0.5 {\rm deg}) = 10^{17.72} \left(\sigma_{\mathrm{los}}/5\,{\rm km \, s^{-1}}\right)^4 \left(d / 100\,{\rm kpc}\right)^{-2}\left( r_{1/2}/100 \,{\rm pc} \right)^{-1}$. We find this relation has considerably smaller scatter as compared to the simpler relations that scale only as $1/d^2$. We further explore scalings with luminosity ($L_V$), and find that the data do not strongly prefer a scaling including $L_V$ as compared to a pure $1/d^2$ scaling. The scaling relations we derive can be used to estimate the J-Factor without the full dynamical analysis, and will be useful for estimating limits on particle dark matter properties from new systems that do not have high-quality stellar kinematics.

Using deep wide-field photometry three-year data (Y3) from the Dark Energy Survey (DES), we present a panoramic study of the Fornax dwarf spheroidal galaxy. The data presented here -- a small subset of the full survey -- uniformly covers a region of 25 square degrees centered on the galaxy to a depth of g ~ 23.5. We use this data to study the structural properties of Fornax, overall stellar population, and its member stars in different evolutionary phases. We also search for possible signs of tidal disturbance. Fornax is found to be significantly more spatially extended than what early studies suggested. No statistically significant distortions or signs of tidal disturbances were found down to a surface brightness limit of ~32.1 mag/arcsec^2. However, there are hints of shell-like features located ~30' - 40' from the center of Fornax that may be stellar debris from past merger events. We also find that intermediate age and young main-sequence populations show different orientation at the galaxy center and have many substructures. The deep DES Y3 data allows us to characterize the age of those substructures with great accuracy, both those previously known and those newly discovered in this work, on the basis of their color-magnitude diagram morphology. We find that the youngest overdensities are all found on the Eastern side of Fornax, where the Fornax field population itself is slightly younger than in the West. The high quality DES Y3 data reveals that Fornax has many rich structures, and provides insights into its complex formation history.

Indirect dark matter searches targeting dwarf spheroidal galaxies (dSphs) have matured rapidly during the past decade. This has been because of the substantial increase in kinematic data sets from the dSphs, the new dSphs that have been discovered, and the operation of the Fermi-LAT and many ground-based gamma-ray experiments. Here we review the analysis methods that have been used to determine the dSph dark matter distributions, in particular the "J-factors," comparing and contrasting them, and detailing the underlying systematics that still affect the analysis. We discuss prospects for improving measurements of dark matter distributions, and how these interplay with future indirect dark matter searches.

Recently the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) reported the detection of a 21cm absorption signal stronger than astrophysical expectations. In this paper we study the impact of radiation from dark matter (DM) decay and primordial black holes (PBH) on the 21cm radiation temperature in the reionization epoch, and impose a constraint on the decaying dark matter and PBH energy injection in the intergalactic medium, which can heat up neutral hydrogen gas and weaken the 21cm absorption signal. We consider decay channels DM$\rightarrow e^+e^-, \gamma\gamma$, $\mu^+\mu^-$, $b\bar{b}$ and the $10^{15-17}$g mass range for primordial black holes, and require the heating of the neutral hydrogen does not negate the 21cm absorption signal. For $e^+e^-$, $\gamma\gamma$ final states and PBH cases we find strong 21cm bounds that can be more stringent than the current extragalactic diffuse photon bounds. For the DM$\rightarrow e^+e^-$ channel, the lifetime bound is $\tau_{\rm DM}> 10^{27}$s for sub-GeV dark matter. The bound is $\tau_{\rm DM}\ge 10^{26}$s for sub-GeV DM$\rightarrow \gamma\gamma$ channel and reaches $10^{27}$s at MeV DM mass. For $b\bar{b}$ and $\mu^+\mu^-$ cases, the 21 cm constraint is better than all the existing constraints for $m_{\rm DM}<20$ GeV where the bound on $\tau_{\rm DM}\ge10^{26}$s. For both DM decay and primordial black hole cases, the 21cm bounds significantly improve over the CMB damping limits from Planck data.

We compare the transverse velocity dispersions recently measured within the Sculptor dwarf spheroidal galaxy to the predictions of previously published dynamical models. These provide good fits to the observed number count and velocity dispersion profiles of metal-rich and metal-poor stars both in cored and in cusped potentials. At the projected radius where the proper motions were measured, these models predict transverse dispersions in the range 6 to 9.5 km/s, with the tangential dispersion about 1 km/s larger than the (projected) radial dispersion. Both dispersions are predicted to be about 1 km/s larger for metal-poor than for metal-rich stars. At this projected radius, cored and cusped potentials predict almost identical transverse dispersions. The measured tangential dispersion ($8.5 \pm 3.2$ km/s) agrees remarkably well with these predictions, while the measured radial dispersion ($11.5 \pm 4.3$ km/s) differs only at about the $1\sigma$ level. Thus, the proper motion data are in excellent agreement with previous data but do not help to distinguish between cored and cusped potentials. This will require velocity dispersion data (either from proper motions or from radial velocities) with uncertainties well below 1 km/s over a range of projected radii.

We perform a search for stellar streams around the Milky Way using the first three years of multi-band optical imaging data from the Dark Energy Survey (DES). We use DES data covering $\sim 5000$ sq. deg. to a depth of $g > 23.5$ with a relative photometric calibration uncertainty of $< 1 \%$. This data set yields unprecedented sensitivity to the stellar density field in the southern celestial hemisphere, enabling the detection of faint stellar streams to a heliocentric distance of $\sim 50$ kpc. We search for stellar streams using a matched-filter in color-magnitude space derived from a synthetic isochrone of an old, metal-poor stellar population. Our detection technique recovers four previously known thin stellar streams: Phoenix, ATLAS, Tucana III, and a possible extension of Molonglo. In addition, we report the discovery of eleven new stellar streams. In general, the new streams detected by DES are fainter, more distant, and lower surface brightness than streams detected by similar techniques in previous photometric surveys. As a by-product of our stellar stream search, we find evidence for extra-tidal stellar structure associated with four globular clusters: NGC 288, NGC 1261, NGC 1851, and NGC 1904. The ever-growing sample of stellar streams will provide insight into the formation of the Galactic stellar halo, the Milky Way gravitational potential, as well as the large- and small-scale distribution of dark matter around the Milky Way.

We examine the bounds on resonantly-produced sterile neutrino dark matter from phase-space densities of Milky Way dwarf spheroidal galaxies (dSphs). The bounds result from a derivation of the dark matter coarse-grained phase-space density from the stellar kinematics, which allows us to explore bounds from some of the most compact dSphs without suffering the resolution limitation from N-body simulations that conventional methods have. We find that the strongest constraints come from very compact dSphs, such as Draco II and Segue 1. We additionally forecast the constraining power of a few dSph candidates that do not yet have associated stellar kinematic data, and show that they can improve the bounds if they are confirmed to be highly dark matter dominated systems. Our results demonstrate that compact dSphs provide important constraints on sterile neutrino dark matter that are comparable to other methods using as Milky Way satellite counts. In particular, if more compact systems are discovered from current or future surveys such as LSST or HSC, it should be possible to test models that explain the 3.5 keV X-ray line signal with a 7.1 keV sterile neutrino particle decay.

We present new proper motion (PM) measurements of the dwarf spheroidal galaxies (dSphs) Draco and Sculptor using multi-epoch images obtained with the Hubble Space Telescope ACS/WFC. Our PM results have uncertainties far lower than previous measurements, even made with the same instrument. The PM results for Draco and Sculptor are (mu_W,mu_N)_Dra = (-0.0562+/-0.0099,-0.1765+/-0.0100) mas/yr and (mu_W,mu_N)_Scl = (-0.0296+/-0.0209,-0.1358 +/-0.0214) mas/yr. The implied Galactocentric velocity vectors for Draco and Sculptor have radial and tangential components: (V_rad,V_tan)_Dra = (-88.6,161.4) +/- (4.4,5.6) km/s; and (V_rad,V_tan)_Scl = (72.6,200.2) +/- (1.3,10.8) km/s. We study the detailed orbital history of both Draco and Sculptor via numerical orbit integrations. Orbital periods of Draco and Sculptor are found to be 1-2 and 2-5 Gyrs, respectively, accounting for uncertainties in the MW mass. We also study the influence of the Large Magellanic Cloud (LMC) on the orbits of Draco and Sculptor. Overall, the inclusion of the LMC increases the scatter in the orbital results. Based on our calculations, Draco shows a rather wide range of orbital parameters depending on the MW mass and inclusion/exclusion of the LMC, but Sculptor's orbit is very well constrained with its most recent pericentric approach to the MW being 0.3-0.4 Gyr ago. Our new PMs imply that the orbital trajectories of both Draco and Sculptor are confined within the Disk of Satellites (DoS), better so than implied by earlier PM measurements, and likely rule out the possibility that these two galaxies were accreted together as part of a tightly bound group.

We consider dark matter annihilation into a general set of final states of Standard Model particles, including two-body and four-body final states that result from the decay of intermediate states. For dark matter masses ~10-10^5 GeV, we use updated data from Planck and from high gamma-ray experiments such as Fermi-LAT, MAGIC, and VERITAS to constrain the annihilation cross section for each final state. The Planck constraints are the most stringent over the entire mass range for annihilation into light leptons, and the Fermi-LAT constraints are the most stringent for four-body final states up to masses ~10^4 GeV. We consider these constraints in light of the recent AMS antiproton results, and show that for light mediators it is possible to explain the AMS data with dark matter, and remain consistent with Fermi-LAT Inner Galaxy measurements, for m_\chi ~ 60-100 GeV mass dark matter and mediator masses m_\phi / m_\chi ~< 1.

Non-standard neutrino interactions (NSI) affect both their propagation through matter and their detection, with bounds on NSI parameters coming from various astrophysical and terrestrial neutrino experiments. In this paper, we show that NSI can be probed in future direct dark matter detection experiments through both elastic neutrino-electron scattering and coherent neutrino-nucleus scattering, and that these channels provide complementary probes of NSI. We show NSI can increase the event rate due to solar neutrinos, with a sharp increase for lower nuclear recoil energy thresholds that are within reach for upcoming detectors. We also identify an interference range of NSI parameters for which the rate is reduced by approximately 40\%. Finally, we show that the "dark side" solution for the solar neutrino mixing angle may be discovered at forthcoming direct detection experiments.

We present chemical abundance measurements of three stars in the ultra-faint dwarf galaxy Horologium I, a Milky Way satellite discovered by the Dark Energy Survey. Using high resolution spectroscopic observations we measure the metallicity of the three stars as well as abundance ratios of several $\alpha$-elements, iron-peak elements, and neutron-capture elements. The abundance pattern is relatively consistent among all three stars, which have a low average metallicity of [Fe/H] $\sim -2.6$ and are not $\alpha$-enhanced ([$\alpha$/Fe] $\sim 0.0$). This result is unexpected when compared to other low-metallicity stars in the Galactic halo and other ultra-faint dwarfs and hints at an entirely different mechanism for the enrichment of Hor I compared to other satellites. We discuss possible scenarios that could lead to this observed nucleosynthetic signature including extended star formation, a Population III supernova, and a possible association with the Large Magellanic Cloud.

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.

For models in which dark matter annihilation is Sommerfeld-enhanced, the annihilation cross section increases at low relative velocities. Dwarf spheroidal galaxies (dSphs) have low characteristic dark matter particle velocities and are thus ideal candidates to study such models. In this paper we model the dark matter phase space of dSphs as isotropic and spherically-symmetric, and determine the $J$-factors for several of the most important targets for indirect dark matter searches. For Navarro-Frenk-White density profiles, we quantify the scatter in the $J$-factor arising from the astrophysical uncertainty in the dark matter potential. We show that, in Sommerfeld-enhanced models, the ordering of the most promising dSphs may be different relative to the standard case of velocity-independent cross sections. This result can have important implications for derived upper limits on the annihilation cross section, or on possible signals, from dSphs.

We present Magellan/IMACS spectroscopy of the recently-discovered Milky Way satellite Eridanus II (Eri II). We identify 28 member stars in Eri II, from which we measure a systemic radial velocity of $v_{\rm hel} = 75.6 \pm 1.3~\mbox{(stat.)} \pm 2.0~\mbox{(sys.)}~\mathrm{km\,s^{-1}}$ and a velocity dispersion of $6.9^{+1.2}_{-0.9}~\mathrm{km\,s^{-1}}$. Assuming that Eri~II is a dispersion-supported system in dynamical equilibrium, we derive a mass within the half-light radius of Eri II is $1.2^{+0.4}_{-0.3} \times 10^{7}~\mathrm{M_\odot}$, indicating a mass-to-light ratio of $420^{+210}_{-140}~\mathrm{M_\odot}/\mathrm{L_\odot}$ and confirming that it is a dark matter-dominated dwarf galaxy. From the equivalent width measurements of the CaT lines of 16 red giant member stars, we derive a mean metallicity of ${\rm [Fe/H]} = -2.38 \pm 0.13$ and a metallicity dispersion of $\sigma_{\rm [Fe/H]} = 0.47 ^{+0.12}_{-0.09}$. The velocity of Eri II in the Galactic Standard of Rest frame is $v_{\rm GSR} = -66.6~\mathrm{km\,s^{-1}}$, indicating that either Eri II is falling into the Milky Way potential for the first time or it has passed the apocenter of its orbit on a subsequent passage. At a Galactocentric distance of $\sim$370 kpc, Eri II is one of the Milky Way's most distant satellites known. Additionally, we show that the bright blue stars previously suggested to be a young stellar population are not associated with Eri II. The lack of gas and recent star formation in Eri II is surprising given its mass and distance from the Milky Way, and may place constraints on models of quenching in dwarf galaxies and on the distribution of hot gas in the Milky Way halo. Furthermore, the large velocity dispersion of Eri II can be combined with the existence of a central star cluster to constrain MACHO dark matter with mass $\gtrsim10~\mathrm{M_\odot}$.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

Line-of-sight integrals of the squared density, commonly called the J-factor, are essential for inferring dark matter annihilation signals. The J-factors of dark matter-dominated dwarf spheroidal satellite galaxies (dSphs) have typically been derived using Bayesian techniques, which for small data samples implies that a choice of priors constitutes a non-negligible systematic uncertainty. Here we report the development of a new fully frequentist approach to construct the profile likelihood of the J-factor. Using stellar kinematic data from several classical and ultra-faint dSphs, we derive the maximum likelihood value for the J-factor and its confidence intervals. We validate this method, in particular its bias and coverage, using simulated data from the Gaia Challenge. We find that the method possesses good statistical properties. The J-factors and their uncertainties are generally in good agreement with the Bayesian-derived values, with the largest deviations restricted to the systems with the smallest kinematic datasets. We discuss improvements, extensions, and future applications of this technique.

We investigate constraints on the abundance of primordial black holes (PBHs) in the mass range 10^{15}-10^{17} g using data from the Cosmic Microwave Background (CMB) and MeV extragalactic gamma-ray background (EGB). Hawking radiation from PBHs with lifetime greater than the age of the universe leaves an imprint on the CMB through modification of the ionization history and the damping of CMB anisotropies. Using a model for redshift dependent energy injection efficiencies, we show that a combination of temperature and polarization data from Planck provides the strongest constraint on the abundance of PBHs for masses \sim 10^{15}-10^{16} g, while the EGB dominates for masses \gtrsim 10^{16} g. Both the CMB and EGB now rule out PBHs as the dominant component of dark matter for masses \sim 10^{16}-10^{17} g. Planned MeV gamma-ray observatories are ideal for further improving constraints on PBHs in this mass range.

The Sculptor dwarf spheroidal galaxy contains two distinct stellar populations: one metal-rich and the other metal-poor. Several authors have argued that in order for these two populations to reside in the same gravitational potential, the dark matter halo must have a core similar to that observed in the stellar count profile. This would rule out the cuspy Navarro-Frenk-White (NFW) density profiles predicted for halos and subhalos by dark matter only simulations of the Lambda Cold Dark Matter (Lambda-CDM) cosmological model. We present a new theoretical framework to analyse stellar count and velocity observations in a self-consistent manner based on separable models, f(E,J)=g(J)h(E), for the distribution function of an equilibrium spherical system. We use this machinery to analyse available photometric and kinematic data for the two stellar populations in Sculptor. We find, contrary to some previous claims, that the data are consistent with populations in equilibrium within an NFW dark matter potential with structural parameters in the range expected in Lambda-CDM. Our solution gives a maximum circular velocity for Sculptor between 20 and 35 km/s. We discuss why some previous authors came to a different conclusion.

We use the APOSTLE cosmological hydrodynamic simulations to examine the effects of tidal stripping on cold dark matter (CDM) sub haloes that host three of the most luminous Milky Way (MW) dwarf satellite galaxies: Fornax, Sculptor, and Leo I. We identify simulated satellites that match the observed spatial and kinematic distributions of stars in these galaxies, and track their evolution after infall. We find $\sim$ 30$\%$ of subhaloes hosting satellites with present-day stellar mass $10^6$-$10^8$ $M_{\odot}$ experience $>20\%$ stellar mass loss after infall. Fornax analogues have earlier infall times compared to Sculptor and Leo I analogues. Star formation in Fornax analogues continues for $\sim3$-$6$ Gyr after infall, whereas Sculptor and Leo I analogues stop forming stars $< 2$-$3$ Gyr after infall. Fornax analogues typically show more significant stellar mass loss and exhibit stellar tidal tails, whereas Sculptor and Leo I analogues, which are more deeply embedded in their host DM haloes at infall, do not show substantial mass loss due to tides. When additionally comparing the orbital motion of the host subaloes to the measured proper motion of Fornax we find the matching more difficult; host subhaloes tend to have pericentres smaller than that measured for Fornax itself. From the kinematic and orbital data, we estimate that Fornax has lost $10-20\%$ of its infall stellar mass. Our best estimate for the surface brightness of a stellar tidal stream associated with Fornax is $\Sigma \sim$ 32.6 mag $ {\rm arcsec^{-2}}$, which may be detectable with deep imaging surveys such as DES and LSST.

We search for excess gamma-ray emission coincident with the positions of confirmed and candidate Milky Way satellite galaxies using 6 years of data from the Fermi Large Area Telescope (LAT). Our sample of 45 stellar systems includes 28 kinematically confirmed dark-matter-dominated dwarf spheroidal galaxies (dSphs) and 17 recently discovered systems that have photometric characteristics consistent with the population of known dSphs. For each of these targets, the relative predicted gamma-ray flux due to dark matter annihilation is taken from kinematic analysis if available, and estimated from a distance-based scaling relation otherwise, assuming that the stellar systems are dark-matter-dominated dSphs. LAT data coincident with four of the newly discovered targets show a slight preference (each ~$2 \sigma$ local) for gamma-ray emission in excess of the background. However, the ensemble of derived gamma-ray flux upper limits for individual targets is consistent with the expectation from analyzing random blank-sky regions, and a combined analysis of the population of stellar systems yields no globally significant excess (global significance $<1 \sigma$). Our analysis has increased sensitivity compared to the analysis of 15 confirmed dSphs by Ackermann et al. 2015. The observed constraints on the dark matter annihilation cross section are statistically consistent with the background expectation, improving by a factor of ~2 for large dark matter masses ($m_{{\rm DM},b \bar b} \gtrsim 1$ TeV and $m_{{\rm DM},\tau^{+}\tau^{-}} \gtrsim 70$ GeV) and weakening by a factor of ~1.5 at lower masses relative to previously observed limits.

We present Magellan/IMACS spectroscopy of the recently discovered Milky Way satellite Tucana III (Tuc III). We identify 26 member stars in Tuc III, from which we measure a mean radial velocity of v_hel = -102.3 +/- 0.4 (stat.) +/- 2.0 (sys.) km/s, a velocity dispersion of 0.1^+0.7_-0.1 km/s, and a mean metallicity of [Fe/H] = -2.42^+0.07_-0.08. The upper limit on the velocity dispersion is sigma < 1.5 km/s at 95.5% confidence, and the corresponding upper limit on the mass within the half-light radius of Tuc III is 9.0 x 10^4 Msun. We cannot rule out mass-to-light ratios as large as 240 Msun/Lsun for Tuc III, but much lower mass-to-light ratios that would leave the system baryon-dominated are also allowed. We measure an upper limit on the metallicity spread of the stars in Tuc III of 0.19 dex at 95.5% confidence. Tuc III has a smaller metallicity dispersion and likely a smaller velocity dispersion than any known dwarf galaxy, but a larger size and lower surface brightness than any known globular cluster. Its metallicity is also much lower than those of the clusters with similar luminosity. We therefore tentatively suggest that Tuc III is the tidally-stripped remnant of a dark matter-dominated dwarf galaxy, but additional precise velocity and metallicity measurements will be necessary for a definitive classification. If Tuc III is indeed a dwarf galaxy, it is one of the closest external galaxies to the Sun. Because of its proximity, the most luminous stars in Tuc III are quite bright, including one star at V=15.7 that is the brightest known member star of an ultra-faint satellite.

We analyze future direct data matter detection experiments using Effective Field Theory (EFT) operators with light, $\lesssim 100$ MeV mass mediators. We compare the nuclear recoil energy spectra from these operators to the predicted high energy solar neutrino spectrum. A set of operators that generate spectra similar to the neutrino background is identified, however this set is distinct from those that mimic the neutrino background for heavy, $\gtrsim 100$ MeV mass mediators. We outline a general classification scheme based on momentum dependence of the dark matter-nucleus interaction to determine how strong the discovery limit for a given operator saturates in the presence of the neutrino background. Our results highlight the benefit of considering a general theoretical framework regarding dark matter, and motivate continued experimental progress towards lower nuclear recoil energy thresholds.

By lowering their energy threshold direct dark matter searches can reach the neutrino floor with experimental technology now in development. The 7Be flux can be detected with $\sim 10$ eV nuclear recoil energy threshold and 50 kg-yr exposure. The pep flux can be detected with $\sim 3$ ton-yr exposure, and the first detection of the CNO flux is possible with similar exposure. The pp flux can be detected with threshold of $\sim$ eV and only $\sim$ kg-yr exposure. These can be the first pure neutral current measurements of the low-energy solar neutrino flux. Measuring this flux is important for low mass dark matter searches and for understanding the solar interior.

We discuss prospects for probing Z-prime and non-standard neutrino interactions using neutrino-nucleus coherent scattering with ultra-low energy (~ 10 eV) threshold Si and Ge detectors. The analysis is performed in the context of a specific and contemporary reactor-based experimental proposal, developed in cooperation with the Nuclear Science Center at Texas A&M University, and referencing available technology based upon economical and scalable detector arrays. For expected exposures, we show that sensitivity to the Z-prime mass is on the order of several TeV, and is complementary to the LHC search with low mass detectors in the near term. This technology is also shown to provide sensitivity to the neutrino magnetic moment, at a level that surpasses terrestrial limits, and is competitive with more stringent astrophysical bounds. We demonstrate the benefits of combining silicon and germanium detectors for distinguishing between classes of models of new physics, and for suppressing correlated systematic uncertainties.

Distinguishing a dark matter interaction from an astrophysical neutrino-induced interaction will be major challenge for future direct dark matter searches. In this paper, we consider this issue within non-relativistic Effective Field Theory (EFT), which provides a well-motivated theoretical framework for determining nuclear responses to dark matter scattering events. We analyze the nuclear energy recoil spectra from the different dark matter-nucleon EFT operators, and compare to the nuclear recoil energy spectra that is predicted to be induced by astrophysical neutrino sources. We determine that for 11 of the 14 possible operators, the dark matter-induced recoil spectra can be cleanly distinguished from the corresponding neutrino-induced recoil spectra with moderate size detector technologies that are now being pursued, e.g., these operators would require 0.5 tonne years to be distinguished from the neutrino background for low mass dark matter. Our results imply that in most models detectors with good energy resolution will be able to distinguish a dark matter signal from a neutrino signal, without the need for much larger detectors that must rely on additional information from timing or direction.

We present the first joint analysis of gamma-ray data from the MAGIC Cherenkov telescopes and the Fermi Large Area Telescope (LAT) to search for gamma-ray signals from dark matter annihilation in dwarf satellite galaxies. We combine 158 hours of Segue 1 observations with MAGIC with 6-year observations of 15 dwarf satellite galaxies by the Fermi-LAT. We obtain limits on the annihilation cross-section for dark matter particle masses between 10 GeV and 100 TeV - the widest mass range ever explored by a single gamma-ray analysis. These limits improve on previously published Fermi-LAT and MAGIC results by up to a factor of two at certain masses. Our new inclusive analysis approach is completely generic and can be used to perform a global, sensitivity-optimized dark matter search by combining data from present and future gamma-ray and neutrino detectors.

The launch of the gamma-ray telescope Fermi Large Area Telescope (Fermi-LAT) started a pivotal period in indirect detection of dark matter. By outperforming expectations, for the first time a robust and stringent test of the paradigm of weakly interacting massive particles (WIMPs) is within reach. In this paper, we discuss astrophysical targets for WIMP detection and the challenges they present, review the analysis tools which have been employed to tackle these challenges, and summarize the status of constraints on and the claimed detections in the WIMP parameter space. Methods and results will be discussed in comparison to Imaging Air Cherenkov Telescopes. We also provide an outlook on short term and longer term developments.

We use cosmological simulations to identify dark matter subhalo host candidates of the Fornax dwarf spheroidal galaxy using the stellar kinematic properties of Fornax. We consider cold dark matter (CDM), warm dark matter (WDM), and decaying dark matter (DDM) simulations for our models of structure formation. The subhalo candidates in CDM typically have smaller mass and higher concentrations at z = 0 than the corresponding candidates in WDM and DDM. We examine the formation histories of the ~ 100 Fornax candidate subhalos identified in CDM simulations and, using approximate luminosity-mass relationships for subhalos, we find two of these subhalos that are consistent with both the Fornax luminosity and kinematics. These two subhalos have a peak mass over ten times larger than their z = 0 mass. We suggest that in CDM the dark matter halo hosting Fornax must have been severely stripped of mass and that it had an infall time into the Milky Way of ~ 9 Gyr ago. In WDM, we find that candidate subhalos consistent with the properties of Fornax have a similar infall time and a similar degree of mass loss, while in DDM we find a later infall time of ~ 3 - 4 Gyr ago and significantly less mass loss. We discuss these results in the context of the Fornax star formation history, and show that these predicted subhalo infall times can be linked to different star formation quenching mechanisms. This emphasizes the links between the properties of the dark matter and the mechanisms that drive galaxy evolution.

The search for weakly interacting massive particles (WIMPs) by direct detection faces an encroaching background due to coherent neutrino-nucleus scattering. As the sensitivity of these experiments improves, the question of how to best distinguish a dark matter signal from neutrinos will become increasingly important. A proposed method of overcoming this so-called 'neutrino floor' is to utilize the directional signature that both neutrino and dark matter induced recoils possess. We show that directional experiments can indeed probe WIMP-nucleon cross-sections below the neutrino floor with little loss in sensitivity due to the neutrino background. In particular we find at low WIMP masses (around 6 GeV) the discovery limits for directional detectors penetrate below the non-directional limit by several orders of magnitude. For high WIMP masses (around 100 GeV), the non-directional limit is overcome by a factor of a few. Furthermore we show that even for directional detectors which can only measure 1- or 2-dimensional projections of the 3-dimensional recoil track, the discovery potential is only reduced by a factor of 3 at most. We also demonstrate that while the experimental limitations of directional detectors, such as sense recognition and finite angular resolution, have a detrimental effect on the discovery limits, it is still possible to overcome the ultimate neutrino background faced by non-directional detectors.

Many particle dark matter models predict that the dark matter undergoes cascade annihilations, i.e. the annihilation products are 4-body final states. In the context of model-independent cascade annihilation models, we study the compatibility of the dark matter interpretation of the Fermi- LAT Galactic center gamma-ray emission with null detections from dwarf spheroidal galaxies. For canonical values of the Milky Way density profile and the local dark matter density, we find that the dark matter interpretation to the Galactic center emission is strongly constrained. However, uncertainties in the dark matter distribution weaken the constraints and leave open dark matter interpretations over a wide range of mass scales.

Due to their proximity, high dark-matter content, and apparent absence of non-thermal processes, Milky Way dwarf spheroidal satellite galaxies (dSphs) are excellent targets for the indirect detection of dark matter. Recently, eight new dSph candidates were discovered using the first year of data from the Dark Energy Survey (DES). We searched for gamma-ray emission coincident with the positions of these new objects in six years of Fermi Large Area Telescope data. We found no significant excesses of gamma-ray emission. Under the assumption that the DES candidates are dSphs with dark matter halo properties similar to the known dSphs, we computed individual and combined limits on the velocity-averaged dark matter annihilation cross section for these new targets. If the estimated dark-matter content of these dSph candidates is confirmed, they will constrain the annihilation cross section to lie below the thermal relic cross section for dark matter particles with masses < 20 GeV annihilating via the b-bbar or tau+tau- channels.

We present Magellan/M2FS, VLT/GIRAFFE, and Gemini South/GMOS spectroscopy of the newly discovered Milky Way satellite Reticulum II. Based on the spectra of 25 Ret II member stars selected from Dark Energy Survey imaging, we measure a mean heliocentric velocity of 62.8 +/- 0.5 km/s and a velocity dispersion of 3.3 +/- 0.7 km/s. The mass-to-light ratio of Ret II within its half-light radius is 470 +/- 210 Msun/Lsun, demonstrating that it is a strongly dark matter-dominated system. Despite its spatial proximity to the Magellanic Clouds, the radial velocity of Ret II differs from that of the LMC and SMC by 199 and 83 km/s, respectively, suggesting that it is not gravitationally bound to the Magellanic system. The likely member stars of Ret II span 1.3 dex in metallicity, with a dispersion of 0.28 +/- 0.09 dex, and we identify several extremely metal-poor stars with [Fe/H] < -3. In combination with its luminosity, size, and ellipticity, these results confirm that Ret II is an ultra-faint dwarf galaxy. With a mean metallicity of [Fe/H] = -2.65 +/- 0.07, Ret II matches Segue~1 as the most metal-poor galaxy known. Although Ret II is the third-closest dwarf galaxy to the Milky Way, the line-of-sight integral of the dark matter density squared is log J = 18.8 +/- 0.6 Gev^2/cm^5 within 0.2 degrees, indicating that the predicted gamma-ray flux from dark matter annihilation in Ret II is lower than that of several other dwarf galaxies.

The recent excess in the CMS measurements of $eejj$ and $e\nu jj$ channels and the emergence of PeV comsic neutrino events at the IceCube experiment share an intriguing implication for a leptoquark with a 600-650 GeV mass. We investigate the CMS constraints on the flavor structure of a scenario with the minimal leptoquark Yukawa couplings and correlate such a scenario to the resonant enhancement in the very high energy shower event rates at the IceCube. We find for a single leptoquark, the CMS signals require large couplings to the third generation leptons. This leads to an enhancement in the $\nu_\tau$-nucleon scattering cross-section and subsequently more $\nu_\tau$ events at PeV energies. However, a visible enhancement above the Standard Model scattering would require a leptoquark Yukawa coupling larger than one that can be easily tested at the upcoming LHC runs.

In this work, we sift a simple supersymmetric framework of late invisible decays to/of the gravitino. We investigate two cases where the gravitino is the lightest supersymmetric particle or the next-to-lightest supersymmetric particle. The next-to-lightest supersymmetric particle decays into two dark matter candidates and has a long lifetime due to gravitationally suppressed interactions. However, because of the absence of any hadronic or electromagnetic products, it satisfies the tight bounds set by big bang nucleosynthesis and cosmic microwaved background. One or both of the dark matter candidates produced in invisible decays can contribute to the amount of dark radiation and suppress perturbations at scales that are being probed by the galaxy power spectrum and the Lyman-alpha forest data. We show that these constraints are satisfied in large regions of the parameter space and, as a result, the late invisible decays to/of the gravitino can be responsible for the entire dark matter relic abundance.

We present the results of a search for unidentified emission lines in deep Suzaku X-ray spectra for the central regions of the four X-ray brightest galaxy clusters: Perseus, Coma, Virgo and Ophiuchus. We employ an optimized energy range for our analysis (3.2-5.3 keV) that is relatively free of instrumental features, and a baseline plasma emission model that incorporates the abundances of elements with the strongest expected emission lines at these energies (S, Ar, Ca) as free parameters. For the Perseus Cluster core, employing this baseline model, we find evidence for an additional emission feature at an energy $3.51^{+0.02}_{-0.01}$ keV with a flux of ~$2.87\times10^{-7}$ ph/s/cm^2/arcmin^2. At slightly larger radii, we detect an emission line at 3.59+/-0.02 keV with a flux of ~$4.8\times10^{-8}$ ph/s/cm^2/arcmin^2. The energies and fluxes of these features are broadly consistent with previous claims, although the radial variation of the line strength appears in tension with standard dark matter (DM) model predictions. Assuming a decaying DM origin for the Perseus emission features allows us to predict the energies and line fluxes for the other clusters in our sample. Critically, we do not detect an emission feature at the predicted energy and line flux in the Coma, Virgo and Ophiuchus clusters. The formal 99.5 per cent upper limits on the strengths of an emission line in each cluster are below the decaying DM model predictions, scaling from the Perseus Cluster center, apparently ruling the model out. In the light of these results, we search for other explanations for the ~3.55 keV emission feature in Perseus. Our results suggest that systematic effects associated with modeling the complex spectra for the Perseus Cluster core, details of the assumed ionization balance, and errors in the predicted emissivities of the spectral lines may in part be responsible for the ~3.5 keV feature.

Dark matter detectors will soon be sensitive to Solar neutrinos via two distinct channels: coherent neutrino-nucleus scattering and neutrino electron elastic scattering. We establish an analysis method for extracting Solar model properties and neutrino properties from these measurements, including the possible effects of sterile neutrinos which have been hinted at by some reactor experiments and cosmological measurements. Even including sterile neutrinos, through the coherent scattering channel a 1 ton-year exposure with a low-threshold Germanium detector could improve on the current measurement of the normalization of the $^8$B Solar neutrino flux down to 3% or less. Combining with the elastic scattering data will provide constraints on both the high and low energy survival probability, and will improve on the uncertainty on the active-to-sterile mixing angle by a factor of two. This sensitivity to active-to-sterile transitions is competitive and complementary to forthcoming dedicated short baseline sterile neutrino searches with nuclear decays.

Direct detection dark matter experiments looking for WIMP-nucleus elastic scattering will soon be sensitive to an irreducible background from neutrinos which will drastically affect their discovery potential. Here we explore how the neutrino background will affect future ton-scale experiments considering both spin-dependent and spin-independent interactions. We show that combining data from experiments using different targets can improve the dark matter discovery potential due to target complementarity. We find that in the context of spin-dependent interactions, combining results from several targets can greatly enhance the subtraction of the neutrino background for WIMP masses below 10 GeV/c$^2$ and therefore probe dark matter models to lower cross-sections. In the context of target complementarity, we also explore how one can tune the relative exposures of different target materials to optimize the WIMP discovery potential.

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 dwarf spheroidal satellite galaxies of the Milky Way are some of the most dark-matter-dominated objects known. Due to their proximity, high dark matter content, and lack of astrophysical backgrounds, dwarf spheroidal galaxies are widely considered to be among the most promising targets for the indirect detection of dark matter via gamma rays. Here we report on gamma-ray observations of 25 Milky Way dwarf spheroidal satellite galaxies based on 4 years of Fermi Large Area Telescope (LAT) data. None of the dwarf galaxies are significantly detected in gamma rays, and we present gamma-ray flux upper limits between 500 MeV and 500 GeV. We determine the dark matter content of 18 dwarf spheroidal galaxies from stellar kinematic data and combine LAT observations of 15 dwarf galaxies to constrain the dark matter annihilation cross section. We set some of the tightest constraints to date on the the annihilation of dark matter particles with masses between 2 GeV and 10 TeV into prototypical Standard Model channels. We find these results to be robust against systematic uncertainties in the LAT instrument performance, diffuse gamma-ray background modeling, and assumed dark matter density profile.

Several direct detection experiments, including recently CDMS-II, have reported signals consistent with 5 to 10 GeV dark matter (DM) that appear to be in tension with null results from XENON and LUX experiments; these indicate a careful review of the theoretical basis, including the galactic DM velocity distribution function (VDF). We establish a VDF parameter space from DM-only cosmological simulations and illustrate that seemingly contradictory experimental results can be made consistent within this parameter space. Future experimental limits should be reported after they are marginalized over a range of VDF parameters.

As direct dark matter experiments continue to increase in size, they will become sensitive to neutrinos from astrophysical sources. For experiments that do not have directional sensitivity, coherent neutrino scattering (CNS) from several sources represents an important background to understand, as it can almost perfectly mimic an authentic WIMP signal. Here we explore in detail the effect of neutrino backgrounds on the discovery potential of WIMPs over the entire mass range of 500 MeV to 10 TeV. We show that, given the theoretical and measured uncertainties on the neutrino backgrounds, direct detection experiments lose sensitivity to light (~10 GeV) and heavy (~100 GeV) WIMPs with a spin-independent cross section below 10^{-45} cm^2 and 10^{-49} cm^2, respectively.

We explore systematic biases in the identification of dark matter in future direct detection experiments and compare the reconstructed dark matter properties when assuming a self-consistent dark matter distribution function and the standard Maxwellian velocity distribution. We find that the systematic bias on the dark matter mass and cross-section determination arising from wrong assumptions for its distribution function is of order ~1\sigma. A much larger systematic bias can arise if wrong assumptions are made on the underlying Milky Way mass model. However, in both cases the bias is substantially mitigated by marginalizing over galactic model parameters. We additionally show that the velocity distribution can be reconstructed in an unbiased manner for typical dark matter parameters. Our results highlight both the robustness of the dark matter mass and cross-section determination using the standard Maxwellian velocity distribution and the importance of accounting for astrophysical uncertainties in a statistically consistent fashion.

We examine the Velocity Distribution Function (VDF) in dark matter halos from Milky Way to cluster mass scales. We identify an empirical model for the VDF with a wider peak and a steeper tail than a Maxwell--Boltzmann distribution, and discuss physical explanations. We quantify sources of scatter in the VDF of cosmological halos and their implication for direct detection of dark matter. Given modern simulations and observations, we find that the most significant uncertainty in the VDF of the Milky Way arises from the unknown radial position of the solar system relative to the dark matter halo scale radius.

For nearly a century, more mass has been measured in galaxies than is contained in the luminous stars and gas. Through continual advances in observations and theory, it has become clear that the dark matter in galaxies is not comprised of known astronomical objects or baryonic matter, and that identification of it is certain to reveal a profound connection between astrophysics, cosmology, and fundamental physics. The best explanation for dark matter is that it is in the form of a yet undiscovered particle of nature, with experiments now gaining sensitivity to the most well-motivated particle dark matter candidates. In this article, I review measurements of dark matter in the Milky Way and its satellite galaxies and the status of Galactic searches for particle dark matter using a combination of terrestrial and space-based astroparticle detectors, and large scale astronomical surveys. I review the limits on the dark matter annihilation and scattering cross sections that can be extracted from both astroparticle experiments and astronomical observations, and explore the theoretical implications of these limits. I discuss methods to measure the properties of particle dark matter using future experiments, and conclude by highlighting the exciting potential for dark matter searches during the next decade, and beyond.

I discuss recent progress in dark matter searches, focusing in particular on how rigorous modeling the dark matter distribution in the Galaxy and in its satellite galaxies improves our interpretation of the limits on the annihilation and elastic scattering cross sections. Looking forward to indirect and direct searches that will operate during the next decade, I review methods for extracting the properties of the dark matter in these experiments in the presence of unknown Galactic model parameters.

We study the abundance of satellites akin to the brightest, classical dwarf spheroidals around galaxies similar in magnitude and isolation to the Milky Way and M31 in the Sloan Digital Sky Survey. From a combination of photometric and spectroscopic redshifts, we bound the mean and the intrinsic scatter in the number of satellites down to ten magnitudes fainter than the Milky Way. Restricting to magnitudes brighter than Sagittarius, we show that the Milky Way is not a significant statistical outlier in its population of classical dwarf spheroidals. At fainter magnitudes, we find an upper limit of 13 on the mean number of satellites brighter than the Fornax dwarf spheroidal. Methods to improve these limits that utilize full photometric redshift distributions hold promise, but are currently limited by incompleteness at the very lowest redshifts. Theoretical models are left to explain why the majority of dark matter subhalos that orbit Milky Way-like galaxies are inefficient at making galaxies at the luminosity scale of the brightest dwarf spheroidals, or why these subhalos predicted by Lambda-CDM do not exist.

We derive bounds on the dark matter annihilation cross-section for low-mass (5-20 GeV) dark matter annihilating primarily to up or down quarks, using the Fermi-LAT bound on gamma-rays from Milky Way satellites. For models in which dark matter-Standard Model interactions are mediated by particular contact operators, we show that these bounds can be directly translated into bounds on the dark matter-proton scattering cross-section. For isospin-violating dark matter, these constraints are tight enough to begin to constrain the parameter-space consistent with experimental signals of low-mass dark matter. We discuss possible models that can evade these bounds.

Numerical simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the gamma-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on gamma-ray spectra consistent with WIMP annihilation through the $b \bar b$ channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the $b \bar b$ channel.

We estimate that there may be up to ~10^5 compact objects in the mass range 10^{-8} -10^{-2} solar mass per main sequence star that are unbound to a host star in the Galaxy. We refer to these objects as nomads; in the literature a subset of these are sometimes called free-floating or rogue planets. Our estimate for the number of Galactic nomads is consistent with a smooth extrapolation of the mass function of unbound objects above the Jupiter-mass scale, the stellar mass density limit, and the metallicity of the interstellar medium. We analyze the prospects for detecting nomads via Galactic microlensing. The Wide-Field Infrared Survey Telescope (WFIRST) will measure the number of nomads per main sequence star greater than the mass of Jupiter to ~ 13%, and the corresponding number greater than the mass of Mars to ~25%. All-sky surveys such as GAIA and LSST can identify nomads greater than about the mass of Jupiter. We suggest a dedicated drift scanning telescope that covers approximately 100 square degrees in the Southern hemisphere could identify nomads as small as 10^{-8} solar mass via microlensing of bright stars with characteristic lightcurve timescales of a few seconds.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 1e-26 cm^3 s^-1 at 5 GeV to about 5e-23 cm^3 s^-1 at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (~3e-26 cm^3 s^-1 for a purely s-wave cross section), without assuming additional boost factors.

We present a search for Galactic dark matter (DM) satellites using the Large Area Telescope (LAT). N-body simulations based on the Lambda-CDM model of cosmology predict a large number of as yet unobserved Galactic DM satellites. These satellites could potentially produce gamma rays through the self-annihilation of DM particles. Some DM satellites are expected to have hard gamma-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We searched for LAT sources with these characteristics. We found no candidate DM satellites matching these criteria in one year of LAT data and interpreted this result in the context of N-body simulations.

We investigate the kinematic properties and stellar population of the Galactic satellite Willman 1 (Wil 1) by combining Keck/DEIMOS spectroscopy with KPNO mosaic camera imaging. Wil 1 is an ultra-low luminosity Milky Way companion. This object lies in a region of size-luminosity space (M_V ~ -2 mag, d ~ 38 kpc, r_half ~ 20 pc) also occupied by the Galactic satellites Bo\"otes II and Segue 1 and 2, but no other known old stellar system. We use kinematic and color-magnitude criteria to identify 45 stars as possible members of Wil 1. With a systemic velocity of v_helio = -12.8 +/- 1.0 km/s, Wil 1 stars have velocities similar to those of foreground Milky Way stars. Informed by Monte-Carlo simulations, we identify 5 of the 45 candidate member stars as likely foreground contaminants. We confirm a significant spread in the abundances of the likely Wil 1 red giant branch members ([Fe/H] = -1.73 +/- 0.12 and -2.65 +/- 0.12, [Ca/Fe] = -0.4 +/- 0.18 and +0.13 +/- 0.28). This spread supports the scenario that Wil 1 is an ultra-low luminosity dwarf galaxy rather than a star cluster. Wil 1's innermost stars move with radial velocities offset by 8 km/s from its outer stars and have a velocity dispersion consistent with 0 km/s, suggesting that Wil 1 may not be in dynamical equilibrium. The combination of the foreground contamination and unusual kinematic distribution make it difficult to robustly determine the dark matter mass of Wil 1. As a result, X-ray or gamma-ray observations of Wil 1 that attempt to constrain models of particle dark matter using an equilibrium mass model are strongly affected by the systematics in the observations presented here. We conclude that, despite the unusual features in the Wil 1 kinematic distribution, evidence indicates that this object is, or at least once was, a dwarf galaxy.

We report the results of the observation of the nearby satellite galaxy Segue 1 performed by the MAGIC-I ground-based gamma-ray telescope between November 2008 and March 2009 for a total of 43.2 hours. No significant gamma-ray emission was found above the background. Differential upper limits on the gamma-ray flux are derived assuming various power-law slopes for the possible emission spectrum. Integral upper limits are also calculated for several power-law spectra and for different energy thresholds. The values are of the order of 10^{-11} ph cm^{-2}$ s^{-1} above 100 GeV and 10^{-12} ph cm^{-2} s^{-1} above 200 GeV. Segue 1 is currently considered one of the most interesting targets for indirect dark matter searches. In these terms, the upper limits have been also interpreted in the context of annihilating dark matter particles. For such purpose, we performed a grid scan over a reasonable portion of the parameter space for the minimal SuperGravity model and computed the flux upper limit for each point separately, taking fully into account the peculiar spectral features of each model. We found that in order to match the experimental upper limits with the model predictions, a minimum flux boost of 10^{3} is required, and that the upper limits are quite dependent on the shape of the gamma-ray energy spectrum predicted by each specific model. Finally we compared the upper limits with the predictions of some dark matter models able to explain the PAMELA rise in the positron ratio, finding that Segue 1 data are in tension with the dark matter explanation of the PAMELA spectrum in the case of a dark matter candidate annihilating into tau+tau-. A complete exclusion however is not possible due to the uncertainties in the Segue 1 astrophysical factor.

We investigate the reconstruction capabilities of Dark Matter mass and spin-independent cross-section from future ton-scale direct detection experiments using germanium, xenon or argon as targets. Adopting realistic values for the exposure, energy threshold and resolution of Dark Matter experiments which will come online within 5 to 10 years, the degree of complementarity between different targets is quantified. We investigate how the uncertainty in the astrophysical parameters controlling the local Dark Matter density and velocity distribution affects the reconstruction. For a 50 GeV WIMP, astrophysical uncertainties degrade the accuracy in the mass reconstruction by up to a factor of $\sim 4$ for xenon and germanium, compared to the case when astrophysical quantities are fixed. However, combination of argon, germanium and xenon data increases the constraining power by a factor of $\sim 2$ compared to germanium or xenon alone. We show that future direct detection experiments can achieve self-calibration of some astrophysical parameters, and they will be able to constrain the WIMP mass with only very weak external astrophysical constraints.

WIMPless dark matter provides a framework in which dark matter particles with a wide range of masses naturally have the correct thermal relic density. We show that WIMPless dark matter with mass around 2-10 GeV can explain the annual modulation observed by the DAMA experiment without violating the constraints of other dark matter searches. This explanation implies distinctive and promising signals for other direct detection experiments, GLAST, and the LHC.

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars within 67 pc (10 arcmin, or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes consistent with membership, down to a magnitude limit of r=21.7. Based on photometric, kinematic, and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, we determine a velocity dispersion of 3.7^{+1.4}_{-1.1} km/s, with a corresponding mass within the half-light radius of 5.8^{+8.2}_{-3.1} x 10^5 Msun. The stellar kinematics of Segue 1 require very high mass-to-light ratios unless the system is far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed toward implausibly short periods. With a total luminosity less than that of a single bright red giant and a V-band mass-to-light ratio of 3400 Msun/Lsun, Segue 1 is the darkest galaxy currently known. We critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < -3) of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest measured dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter physics and galaxy formation on small scales.

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

Detecting the dark matter annihilation signal from Galactic substructure, or subhalos, is an important challenge for high-energy gamma-ray experiments. In this paper we discuss detection prospects by combining two different aspects of the gamma-ray signal: the angular distribution and the photon counts probability distribution function (PDF). The true PDF from subhalos has been shown recently (by Lee et al.) to deviate from Poisson; we extend this analysis and derive the signal PDF from a detailed Lambda-CDM-based model for the properties of subhalos. We combine our PDF with a model for Galactic and extra-Galactic diffuse gamma-ray emission to obtain an estimator and projected error on dark matter particle properties (mass and annihilation cross section) using the Fermi Gamma-Ray Space Telescope. We compare the estimator obtained from the true PDF to that obtained from the simpler Poisson analysis. We find that, although both estimators are unbaised in the presence of backgrounds, the error on dark matter properties derived from the true PDF is ~50% smaller than when utilizing the Poisson-based analysis.

We use new kinematic data from the ultra-faint Milky Way satellite Segue 1 to model its dark matter distribution and derive upper limits on the dark matter annihilation cross-section. Using gamma-ray flux upper limits from the Fermi satellite and MAGIC, we determine cross-section exclusion regions for dark matter annihilation into a variety of different particles including charged leptons. We show that these exclusion regions are beginning to probe the regions of interest for a dark matter interpretation of the electron and positron fluxes from PAMELA, Fermi, and HESS, and that future observations of Segue 1 have strong prospects for testing such an interpretation. We additionally discuss prospects for detecting annihilation with neutrinos using the IceCube detector, finding that in an optimistic scenario a few neutrino events may be detected. Finally we use the kinematic data to model the Segue 1 dark matter velocity dispersion and constrain Sommerfeld enhanced models.

We investigate whether the subhalos of Lambda-CDM galaxy halos have potentials consistent with the observed properties of Milky Way satellites, particularly those with high-quality photometric and kinematic data: Fornax, Leo I, Sculptor, Sextans, and Carina. We compare spherical models with isotropic velocity dispersion tensors to the observed, circularly averaged star counts, line-of-sight velocity dispersion profiles and line-of-sight velocity distributions. We identify subhalos within the six high-resolution dark matter halos of the Aquarius Project for which the spherically averaged potentials result in excellent fits to each of the five galaxies. In particular, our simple one-integral models reproduce the observations in the inner regions, proving that these data are fully consistent with Lambda-CDM expectations and do not require cored dark matter distributions. For four of the five satellites the fits require moderately cusped {\it stellar} density profiles. The star count data for Leo I, however, do require a cored distribution of star counts. Current data suggest that these five satellites may be hosted by Lambda-CDM subhalos with maximum circular velocities in the range 10 to 30 km/s.

We present a new method for determining Weakly Interacting Massive Particle (WIMP) properties in future tonne scale direct detection experiments which accounts for uncertainties in the Milky Way (MW) smooth dark matter distribution. Using synthetic data on the kinematics of MW halo stars matching present samples from the Sloan Digital Sky Survey, complemented by local escape velocity constraints, we demonstrate that the local dark matter density can be constrained to approximately 20% accuracy. For low mass WIMPs, we find that a factor of two error in the assumed local dark matter density leads to a severely biased reconstruction of the WIMP spin-independent cross section that is incorrect at the 15-sigma level. We show that this bias may be overcome by marginalizing over parameters that describe the MW potential, and use this formalism to project the accuracy attainable on WIMP properties in future 1 tonne Xenon detectors. Our method can be readily applied to different detector technologies and extended to more detailed MW halo models.

We use the publicly available subhalo catalogs from the Via Lactea simulation along with a Gpc-scale N-body simulation to understand the impact of inhomogeneous reionization on the satellite galaxy population of the Milky Way. The large-volume simulation is combined with a model for reionization that allows us to predict the distribution of reionization times for Milky Way mass halos. Motivated by this distribution, we identify candidate satellite galaxies in the simulation by requiring that any subhalo must grow above a specified mass threshold before it is reionized; after this time the photoionizing background will suppress both the formation of stars and the accretion of gas. We show that varying the reionization time over the range expected for Milky Way mass halos can change the number of satellite galaxies by roughly two orders of magnitude. This conclusion is in contradiction with a number of studies in the literature, and we conclude that this is a result of inconsistent application of the results of Gnedin (2000). We compare our satellite galaxies to observations using both abundance matching and stellar population synthesis methods to assign luminosities to our subhalos and account for observational completeness effects. Additionally, if we assume that the mass threshold is set by the virial temperature Tvir = 8e3K we find that our model accurately matches the vmax distribution, radial distribution, and luminosity function of observed Milky Way satellites for a reionization time zreion = 9.6^{1.0}_{-2.1}, assuming that the Via Lacteasubhalo distribution is representative of the Milky Way. This results in the presence of 119^{+202}_{-50} satellite galaxies.

In the past several years high resolution kinematic data sets from Milky Way satellite galaxies have confirmed earlier indications that these systems are dark matter dominated objects. Further understanding of what these galaxies reveal about cosmology and the small scale structure of dark matter relies in large part on a more detailed interpretation of their internal kinematics. This article discusses a likelihood formalism that extracts important quantities from the kinematic data, including the amplitude of rotation, proper motion, and the mass distribution. In the simplest model the projected error on the rotational amplitude is shown to be $\sim 0.5 $ km s$^{-1}$ with $\sim 10^3$ stars from either classical or ultra-faint satellites. The galaxy Sculptor is analyzed for the presence of a rotational signal; no significant detection of rotation is found, and given this result limits are derived on the Sculptor proper motion. A criteria for model selection is discussed that determines the parameters required to describe the dark matter halo density profiles and the stellar velocity anisotropy. Applied to four data sets with a wide range of velocities, the likelihood is found to be more sensitive to variations in the slope of the dark matter density profile than variations in the velocity anisotropy. Models with variable radial velocity anisotropy are shown to be preferred relative to those in which this quantity is constant at all radii in the galaxy.

Neutrino-induced recoil events may constitute a background to direct dark matter searches, particularly for those detectors that strive to reach the ton-scale and beyond. This paper discusses the expected neutrino-induced background spectrum due to several of the most important sources, including solar, atmospheric, and diffuse supernova neutrinos. The largest rate arises from $^8$B produced solar neutrinos, providing upwards of $\sim 10^3$ events per ton-year over all recoil energies for the heaviest nuclear targets. However the majority of these $^8$B events are expected to be below the recoil threshold of modern detectors. The remaining neutrino sources are found to constitute a background to the WIMP-induced recoil rate only if the WIMP-nucleon cross section is less than $10^{-12}$ pb. Finally the sensitivity to diffuse supernova neutrino flux for non-electron neutrino flavors is discussed, and projected flux limits are compared with existing flux limits.

We provide conservative bounds on the dark matter cross-section and lifetime from final state radiation produced by annihilation or decay into charged leptons, either directly or via an intermediate particle $\phi$. Our analysis utilizes the experimental gamma-ray flux upper limits from four Milky Way dwarf satellites: HESS observations of Sagittarius and VERITAS observations of Draco, Ursa Minor, and Willman 1. Using 90% confidence level lower limits on the integrals over the dark matter distributions, we find that these constraints are largely unable to rule out dark matter annihilations or decays as an explanation of the PAMELA and ATIC/PPB-BETS excesses. However, if there is an additional Sommerfeld enhancement in dwarfs, which have a velocity dispersion ~10 to 20 times lower than that of the local Galactic halo, then the cross-sections for dark matter annihilating through $\phi$'s required to explain the excesses are very close to the cross-section upper bounds from Willman 1. Dark matter annihilation directly into $\tau$'s is also marginally ruled out by Willman 1 as an explanation of the excesses, and the required cross-section is only a factor of a few below the upper bound from Draco. Finally, we make predictions for the gamma-ray flux expected from the dwarf galaxy Segue 1 for the Fermi Gamma-ray Space Telescope. We find that for a sizeable fraction of the parameter space in which dark matter annihilation into charged leptons explains the PAMELA excess, Fermi has good prospects for detecting a gamma-ray signal from Segue 1 after one year of observation.

We present a general methodology for determining the gamma-ray flux from annihilation of dark matter particles in Milky Way satellite galaxies, focusing on two promising satellites as examples: Segue 1 and Draco. We use the SuperBayeS code to explore the best-fitting regions of the Constrained Minimal Supersymmetric Standard Model (CMSSM) parameter space, and an independent MCMC analysis of the dark matter halo properties of the satellites using published radial velocities. We present a formalism for determining the boost from halo substructure in these galaxies and show that its value depends strongly on the extrapolation of the concentration-mass (c(M)) relation for CDM subhalos down to the minimum possible mass. We show that the preferred region for this minimum halo mass within the CMSSM with neutralino dark matter is ~10^-9-10^-6 solar masses. For the boost model where the observed power-law c(M) relation is extrapolated down to the minimum halo mass we find average boosts of about 20, while the Bullock et al (2001) c(M) model results in boosts of order unity. We estimate that for the power-law c(M) boost model and photon energies greater than a GeV, the Fermi space-telescope has about 20% chance of detecting a dark matter annihilation signal from Draco with signal-to-noise greater than 3 after about 5 years of observation.

We consider the prospects for testing the dark matter interpretation of the DAMA/LIBRA signal with the Super-Kamiokande experiment. The DAMA/LIBRA signal favors dark matter with low mass and high scattering cross section. We show that these characteristics imply that the scattering cross section that enters the DAMA/LIBRA event rate determines the annihilation rate probed by Super-Kamiokande. Current limits from Super-Kamiokande through-going events do not test the DAMA/LIBRA favored region. We show, however, that upcoming analyses including fully-contained events with sensitivity to dark matter masses from 5 to 10 GeV may corroborate the DAMA/LIBRA signal. We conclude by considering three specific dark matter candidates, neutralinos, WIMPless dark matter, and mirror dark matter, which illustrate the various model-dependent assumptions entering our analysis.

Over the past five years, searches in Sloan Digital Sky Survey data have more than doubled the number of known dwarf satellite galaxies of the Milky Way, and have revealed a population of ultra-faint galaxies with luminosities smaller than typical globular clusters, L ~ 1000 Lsun. These systems are the faintest, most dark matter dominated, and most metal poor galaxies in the universe. Completeness corrections suggest that we are poised on the edge of a vast discovery space in galaxy phenomenology, with hundreds more of these extreme galaxies to be discovered as future instruments hunt for the low-luminosity threshold of galaxy formation. Dark matter dominated dwarfs of this kind probe the small-scale power-spectrum, provide the most stringent limits on the phase-space packing of dark matter, and offer a particularly useful target for dark matter indirect detection experiments. Full use of dwarfs as dark matter laboratories will require synergy between deep, large-area photometric searches; spectroscopic and astrometric follow-up with next-generation optical telescopes; and subsequent observations with gamma-ray telescopes for dark matter indirect detection.

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data.

The concordance Lambda Cold Dark Matter (Lambda-CDM) model for the formation of structure in the Universe, while remarkably successful at describing observations of structure on large scales, continues to be challenged by observations on galactic scales. Fortunately, CDM models and their various proposed alternatives make a rich variety of testable predictions that make the Local Group and our own Milky Way Galaxy key laboratories for exploring dark matter (DM) in this regime. However, some of the most definitive tests of local DM require microarcsecond astrometry of faint sources, an astrometric regime that is a unique niche of SIM Lite. This chapter explores the important and distinct contributions that can be made by SIM Lite in the exploration of galaxy dynamics and DM on galaxy scales and that have cosmological consequences. Key areas of potential SIM Lite exploration include (1) measuring the shape, orientation, density law, and lumpiness of the dark halo of the Milky Way and other nearby galaxies, (2) determining the orbits of Galactic satellites, which may be representatives of late infall from the hierarchical formation of the Milky Way, (3) ascertaining the distribution of angular momentum and orbital anisotropy of stars and globular clusters to the outer reaches of the Galactic halo, dynamical properties that hold clues to the early hierarchical formation of the Galaxy, (4) measuring the physical nature of DM by placing strong constraints on the phase space density in the cores of nearby dSph galaxies, and (5) reconstructing the dynamical history of the Local Group through the determination of orbits and masses of its constituent galaxies.

We correct the observed Milky Way satellite luminosity function for luminosity bias using published completeness limits for the Sloan Digital Sky Survey DR5. Assuming that the spatial distribution of Milky Way satellites tracks the subhalos found in the Via Lactea LCDM N-body simulation, we show that there should be between ~300 and ~600 satellites within 400 kpc of the Sun that are brighter than the faintest known dwarf galaxies, and that there may be as many as ~ 1000, depending on assumptions. By taking into account completeness limits, we show that the radial distribution of known Milky Way dwarfs is consistent with our assumption that the full satellite population tracks that of subhalos. These results alleviate the primary worries associated with the so-called missing satellites problem in CDM. We show that future, deep wide-field surveys such as SkyMapper, the Dark Energy Survey (DES), PanSTARRS, and the Large Synoptic Survey Telescope (LSST) will deliver a complete census of ultra-faint dwarf satellites out to the Milky Way virial radius, offer new limits on the free-streaming scale of dark matter, and provide unprecedented constraints on the low-luminosity threshold of galaxy formation.

We present Keck/DEIMOS spectroscopy of Segue 1, an ultra-low luminosity (M_V = -1.5) Milky Way satellite companion. While the combined size and luminosity of Segue 1 are consistent with either a globular cluster or a dwarf galaxy, we present spectroscopic evidence that this object is a dark matter-dominated dwarf galaxy. We identify 24 stars as members of Segue 1 with a mean heliocentric recession velocity of 206 +/- 1.3 kms. We measure an internal velocity dispersion of 4.3+/-1.2 kms. Under the assumption that these stars are in dynamical equilibrium, we infer a total mass of 4.5^{+4.7}_{-2.5} x 10^5 Msun in the case where mass-follow-light; using a two-component maximum likelihood model, we determine a similar mass within the stellar radius of 50 pc. This implies a mass-to-light ratio of ln(M/L_V) = 7.2^{+1.1}_{-1.2} or M/L_V = 1320^{+2680}_{-940}. The error distribution of the mass-to-light ratio is nearly log-normal, thus Segue 1 is dark matter-dominated at a high significance. Using spectral synthesis modeling, we derive a metallicity for the single red giant branch star in our sample of [Fe/H]=-3.3 +/- 0.2 dex. Finally, we discuss the prospects for detecting gamma-rays from annihilation of dark matter particles and show that Segue 1 is the most promising satellite for indirect dark matter detection. We conclude that Segue 1 is the least luminous of the ultra-faint galaxies recently discovered around the Milky Way, and is thus the least luminous known galaxy.

NGC 3621 is a late-type (Sd) spiral galaxy with an active nucleus, previously detected through mid-infrared [Ne V] line emission. Archival Hubble Space Telescope (HST) images reveal that the galaxy contains a bright and compact nuclear star cluster. We present a new high-resolution optical spectrum of this nuclear cluster, obtained with the ESI Spectrograph at the Keck Observatory. The nucleus has a Seyfert 2 emission-line spectrum at optical wavelengths, supporting the hypothesis that a black hole is present. The line-of-sight stellar velocity dispersion of the cluster is sigma=43+/-3 km/s, one of the largest dispersions measured for any nuclear cluster in a late-type spiral galaxy. Combining this measurement with structural parameters measured from archival HST images, we carry out dynamical modeling based on the Jeans equation for a spherical star cluster containing a central point mass. The maximum black hole mass consistent with the measured stellar velocity dispersion is 3*10^6 solar masses. If the black hole mass is small compared with the cluster's stellar mass, then the dynamical models imply a total stellar mass of ~1*10^7 solar masses, which is consistent with rough estimates of the stellar mass based on photometric measurements from HST images. From structural decomposition of 2MASS images, we find no clear evidence for a bulge in NGC 3621; the galaxy contains at most a very faint and inconspicuous pseudobulge component (M_K>-17.6 mag). NGC 3621 provides one of the best demonstrations that very late-type spirals can host both active nuclei and nuclear star clusters, and that low-mass black holes can occur in disk galaxies even in the absence of a substantial bulge.

The Milky Way has at least twenty-three known satellite galaxies that shine with luminosities ranging from about a thousand to a billion times that of the Sun. Half of these galaxies were discovered in the past few years in the Sloan Digital Sky Survey, and they are among the least luminous galaxies in the known Universe. A determination of the mass of these galaxies provides a test of galaxy formation at the smallest scales and probes the nature of the dark matter that dominates the mass density of the Universe. Here we use new measurements of the velocities of the stars in these galaxies to show that they are consistent with them having a common mass of about 10^7 M_\odot within their central 300 parsecs. This result demonstrates that the faintest of the Milky Way satellites are the most dark matter-dominated galaxies known, and could be a hint of a new scale in galaxy formation or a characteristic scale for the clustering of dark matter.

Proper motions for several Milky Way dwarf spheroidal (dSph) galaxies have been determined using both ground and space-based imaging. These measurements require long baselines and repeat observations and typical errors are of order ten milli-arcseconds per century. In this paper, we utilize the effect of "perspective rotation" to show that systematic proper motion of some dSphs can be determined to a similar precision using only stellar line-of-sight velocities. We show that including the effects of small intrinsic rotation in dSphs increases the proper motion errors by about a factor of two. We provide error projections for future data sets, and show that proposed thirty meter class telescopes will measure the proper motion of a few dSphs with milli-arcsecond per century precision.

The origin of both the diffuse high-latitude MeV gamma-ray emission and the 511 keV line flux from the Galactic bulge are uncertain. Previous studies have invoked dark matter physics to independently explain these observations, though as yet none has been able to explain both of these emissions within the well-motivated framework of Weakly-Interacting Massive Particles (WIMPs). Here we use an unstable WIMP dark matter model to show that it is in fact possible to simultaneously reconcile both of these observations, and in the process show a remarkable coincidence: decaying dark matter with MeV mass splittings can explain both observations if positrons and photons are produced with similar branching fractions. We illustrate this idea with an unstable branon, which is a standard WIMP dark matter candidate appearing in brane world models with large extra dimensions. We show that because branons decay via three-body final states, they are additionally unconstrained by searches for Galactic MeV gamma-ray lines. As a result, such unstable long-lifetime dark matter particles provide novel and distinct signatures that can be tested by future observations of MeV gamma-rays.

Dark matter particles need not be completely stable, and in fact they may be decaying now. We consider this possibility in the frameworks of universal extra dimensions and supersymmetry with very late decays of WIMPs to Kaluza-Klein gravitons and gravitinos. The diffuse photon background is a sensitive probe, even for lifetimes far greater than the age of the Universe. Remarkably, both the energy spectrum and flux of the observed MeV gamma ray excess may be naturally explained by decaying dark matter with MeV mass splittings. Future observations of continuum and line photon fluxes will test this explanation and may provide novel constraints on cosmological parameters.

Precision astrometry at microarcsecond accuracy has application to a wide range of astrophysical problems. This paper is a study of the science questions that can be addressed using an instrument that delivers parallaxes at about 4 microarcsec on targets as faint as V = 20, differential accuracy of 0.6 microarcsec on bright targets, and with flexible scheduling. The science topics are drawn primarily from the Team Key Projects, selected in 2000, for the Space Interferometry Mission PlanetQuest (SIM PlanetQuest). We use the capabilities of this mission to illustrate the importance of the next level of astrometric precision in modern astrophysics. SIM PlanetQuest is currently in the detailed design phase, having completed all of the enabling technologies needed for the flight instrument in 2005. It will be the first space-based long baseline Michelson interferometer designed for precision astrometry. SIM will contribute strongly to many astronomical fields including stellar and galactic astrophysics, planetary systems around nearby stars, and the study of quasar and AGN nuclei. SIM will search for planets with masses as small as an Earth orbiting in the `habitable zone' around the nearest stars using differential astrometry, and could discover many dozen if Earth-like planets are common. It will be the most capable instrument for detecting planets around young stars, thereby providing insights into how planetary systems are born and how they evolve with time. SIM will observe significant numbers of very high- and low-mass stars, providing stellar masses to 1%, the accuracy needed to challenge physical models. Using precision proper motion measurements, SIM will probe the galactic mass distribution and the formation and evolution of the Galactic halo. (abridged)

Numerical simulations of Milky-Way size Cold Dark Matter (CDM) halos predict a steeply rising mass function of small dark matter subhalos and a substructure count that greatly outnumbers the observed satellites of the Milky Way. Several proposed explanations exist, but detailed comparison between theory and observation in terms of the maximum circular velocity (Vmax) of the subhalos is hampered by the fact that Vmax for satellite halos is poorly constrained. We present comprehensive mass models for the well-known Milky Way dwarf satellites, and derive likelihood functions to show that their masses within 0.6 kpc (M_0.6) are strongly constrained by the present data. We show that the M_0.6 mass function of luminous satellite halos is flat between ~ 10^7 and 10^8 M_\odot. We use the ``Via Lactea'' N-body simulation to show that the M_0.6 mass function of CDM subhalos is steeply rising over this range. We rule out the hypothesis that the 11 well-known satellites of the Milky Way are hosted by the 11 most massive subhalos. We show that models where the brightest satellites correspond to the earliest forming subhalos or the most massive accreted objects both reproduce the observed mass function. A similar analysis with the newly-discovered dwarf satellites will further test these scenarios and provide powerful constraints on the CDM small-scale power spectrum and warm dark matter models.

We use kinematic data from three new, nearby, extremely low-luminosity Milky Way dwarf galaxies (Ursa Major II, Willman 1, and Coma Berenices) to constrain the properties of their dark matter halos, and from these make predictions for the gamma-ray flux from annihilation of dark matter particles in these halos. We show that these 10^3 solar luminosity dwarfs are the most dark matter dominated galaxies in the Universe, with total masses within 100 pc in excess of 10^6 solar masses. Coupled with their relative proximity, their large masses imply that they should have mean gamma-ray fluxes comparable to or greater than any other known satellite galaxy of the Milky Way. Our results are robust to both variations of the inner slope of the density profile and the effect of tidal interactions. The fluxes could be boosted by up to two orders of magnitude if we include the density enhancements caused by surviving dark matter substructure.

The instability of dark matter may produce visible signals in the spectrum of cosmic gamma-rays. We consider this possibility in frameworks with additional spatial dimensions and supersymmetry. Examples of particles include superweakly-interacting massive particles such as gravitinos in supersymmetry models, the lightest Kaluza-Klein (KK) state in models with universal extra dimensions, and weakly-interacting massive particles such as branons in flexible brane-worlds.

We examine the prospects for detecting gamma-rays from dark matter annihilation in the six most promising dwarf spheroidal (dSph) satellite galaxies of the Milky Way. We use recently-measured velocity dispersion profiles to provide a systematic investigation of the dark matter mass distribution of each galaxy, and show that the uncertainty in the gamma-ray flux from mass modeling is less than a factor of ~ 5 for each dSph if we assume a smooth NFW profile. We show that Ursa Minor and Draco are the most promising dSphs for gamma-ray detection with GLAST and other planned observatories. For each dSph, we investigate the flux enhancement resulting from halo substructure, and show that the enhancement factor relative to a smooth halo flux cannot be greater than about 100. This enhancement depends very weakly on the lower mass cut-off scale of the substructure mass function. While the amplitude of the expected flux from each dSph depends sensitively on the dark matter model, we show that the flux ratios between the six Sphs are known to within a factor of about 10. The flux ratios are also relatively insensitive to the current theoretical range of cold dark matter halo central slopes and substructure fractions.

A study of the current status of WIMP dark matter searches has been made in the context of scientific and technical planning for a Deep Underground Science and Engineering Laboratory (DUSEL) in the U.S. The table of contents follows: 1. Overview 2. WIMP Dark Matter: Cosmology, Astrophysics, and Particle Physics 3. Direct Detection of WIMPs 4. Indirect Detection of WIMPs 5. Dark Matter Candidates and New Physics in the Laboratory 6. Synergies with Other Sub-Fields 7. Direct Detection Experiments: Status and Future Prospects 8. Infrastructure 9. International Context 10. Summary and Outlook 11. Acknowledgments

WIMPs with electroweak scale masses (neutralinos, etc.) remain in kinetic equilibrium with other particle species until temperatures approximately in the range of 10 MeV to 1 GeV, leading to the formation of dark matter substructure with masses as small as $10^{-4} M_{\odot}$ to $10^{-12} M_{\odot}$. However, if dark matter consists of particles with MeV scale masses, as motivated by the observation of 511 keV emission from the Galactic Bulge, such particles are naturally expected to remain in kinetic equilibrium with the cosmic neutrino background until considerably later times. This would lead to a strong suppression of small scale structure with masses below about $10^7 M_{\odot}$ to $10^4 M_{\odot}$. This cutoff scale has important implications for present and future searches for faint Local Group satellite galaxies and for the missing satellites problem.

We show that dark matter emerging from late decays (z < 1000) produces a linear power spectrum identical to that of Cold Dark Matter (CDM) on all observationally relavant scales (> 0.1 Mpc), and simultaneously generates observable constant-density cores in small dark matter halos. We refer to this class of models as meta-Cold Dark Matter (mCDM), because it is born with non-relativistic velocities from the decays of cold thermal relics. The constant-density cores are a result of the low phase-space density of mCDM at birth. Warm dark matter cannot produce similar size phase-space limited cores without saturating the Ly-alpha power spectrum bounds. Dark matter dominated galaxy rotation curves and stellar velocity dispersion profiles may provide the best means to discriminate between mCDM and CDM. mCDM candidates are motivated by the particle spectrum of supersymmetric and extra dimensional extensions to the standard model of particle physics.

We show that measurements of stellar proper motions in dwarf spheroidal galaxies provide a powerful probe of the nature of dark matter. Allowing for general dark matter density profiles and stellar velocity anisotropy profiles, we show that the log-slope of the dark matter profile at about twice the stellar core (King) radius can be measured to within \pm 0.2 when the proper motions of 200 stars are added to standard line-of-sight velocity dispersion data. This measurement of the log-slope provides a test of Cold and Warm Dark Matter theories at a sensitivity not possible with line-of-sight velocity dispersion measurements alone. The upcoming SIM PlanetQuest will have the sensitivity to obtain the required number of proper motions in Milky Way dwarf spheroidal galaxies.

Minimal universal extra dimensions (mUED) is often thought to predict that the lightest Kaluza-Klein particle (LKP) is the Kaluza-Klein gauge boson B^1, leading to conventional missing energy signals at colliders and WIMP dark matter. In fact, the implications of mUED are far richer: the B^1, charged Higgs boson H^{\pm 1}, and graviton G^1 are all possible LKPs, leading to many different ``phases'' with distinct signatures. Considering the complete phase diagram, we find predictions for charged or neutral particles with decay lengths of microns to tens of meters; WIMP, superWIMP, or charged relic particles; metastable particles with lifetimes of the order of or in excess of the age of the universe; and scenarios combining two or more of these phenomena. In the cosmologically preferred region, the Higgs boson mass is between 180 and 245 GeV, the LKP mass is between 810 and 1400 GeV, and the maximal splitting between first Kaluza-Klein modes is less than 320 GeV. This region predicts a variety of exotic collider signals, such as slow charged particles, displaced vertices, tracks with non-vanishing impact parameters, track kinks, and even vanishing charged tracks, all of which provide early discovery possibilities at the Large Hadron Collider.

We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its core radius to be no larger than 700 pc. We derive an upper limit core radius of 300 pc by demanding that the central phase space density of Fornax not exceed that directly inferred from the rotation curves of low-mass spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the core to be quite small in order to avoid conservative limits from the Ly alpha forest power spectrum, implying a core radius < 85 pc. We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can be explained by the presence of a large 1.5 kpc core. A self-consistent core of this size would be drastically inconsistent with the expectations of standard warm or cold dark matter models, and would also require an unreasonably massive dark matter halo, with a maximum circular velocity of 200 km/s.

Supernovae are rare nearby, but they are not rare in the Universe, and all past core-collapse supernovae contributed to the Diffuse Supernova Neutrino Background (DSNB), for which the near-term detection prospects are very good. The Super-Kamiokande limit on the DSNB electron {\it antineutrino} flux, $\phi(E_\nu > 19.3 {\rm MeV}) < 1.2$ cm$^{-2}$ s$^{-1}$, is just above the range of recent theoretical predictions based on the measured star formation rate history. We show that the Sudbury Neutrino Observatory should be able to test the corresponding DSNB electron {\it neutrino} flux with a sensitivity as low as $\phi(22.5 < E_\nu < 32.5 {\rm MeV}) \simeq 6 $ cm$^{-2}$ s$^{-1}$, improving the existing Mont Blanc limit by about three orders of magnitude. While conventional supernova models predict comparable electron neutrino and antineutrino fluxes, it is often considered that the first (and forward-directed) SN 1987A event in the Kamiokande-II detector should be attributed to electron-neutrino scattering with an electron, which would require a substantially enhanced electron neutrino flux. We show that with the required enhancements in either the burst or thermal phase $\nu_e$ fluxes, the DSNB electron neutrino flux would generally be detectable in the Sudbury Neutrino Observatory. A direct experimental test could then resolve one of the enduring mysteries of SN 1987A: whether the first Kamiokande-II event reveals a serious misunderstanding of supernova physics, or was simply an unlikely statistical fluctuation. Thus the electron neutrino sensitivity of the Sudbury Neutrino Observatory is an important complement to the electron antineutrino sensitivity of Super-Kamiokande in the quest to understand the DSNB.

We constrain the Cosmic Star Formation Rate (CSFR) by requiring that massive stars produce the observed UV, optical, and IR light while at the same time not overproduce the Diffuse Supernova Neutrino Background as bounded by Super-Kamiokande. With the massive star component so constrained we then show that a reasonable choice of stellar Initial Mass Function and other parameters results in SNIa rates and iron yields in good agreement with data. In this way we define a `concordance' CSFR that predicts the optical SNII rate and the SNIa contribution to the MeV Cosmic Gamma-Ray Background. The CSFR constrained to reproduce these and other proxies of intermediate and massive star formation is more clearly delineated than if it were measured by any one technique and has the following testable consequences: (1) SNIa contribute only a small fraction of the MeV Cosmic Gamma-Ray Background, (2) massive star core-collapse is nearly always accompanied by a successful optical SNII, and (3) the Diffuse Supernova Neutrino Background is tantalizingly close to detectability.

We investigate the effect of dark energy on the density profiles of dark matter haloes with a suite of cosmological N-body simulations and use our results to test analytic models. We consider constant equation of state models, and allow both w>-1 and w<-1. Using five simulations with w ranging from -1.5 to -0.5, and with more than ~1600 well-resolved haloes each, we show that the halo concentration model of Bullock et al. (2001) accurately predicts the median concentrations of haloes over the range of w, halo masses, and redshifts that we are capable of probing. We find that the Bullock et al. (2001) model works best when halo masses and concentrations are defined relative to an outer radius set by a cosmology-dependent virial overdensity. For a fixed power spectrum normalization and fixed-mass haloes, larger values of w lead to higher concentrations and higher halo central densities, both because collapse occurs earlier and because haloes have higher virial densities. While precise predictions of halo densities are quite sensitive to various uncertainties, we make broad comparisons to galaxy rotation curve data. At fixed power spectrum normalization (fixed sigma_8), w>-1 quintessence models seem to exacerbate the central density problem relative to the standard w=-1 model. Meanwhile w<-1 models help to reduce the apparent discrepancy. We confirm that the Jenkins et al. (2001) halo mass function provides an excellent approximation to the abundance of haloes in our simulations and extend its region of validity to include models with w<-1.

We calculate the Supernova Relic Neutrino (SRN) background flux for the KamLAND and Super-Kamiokande (SK) detectors, motivated by the reduction in background at SK and new results for the star formation history (e.g., from the Sloan Digital Sky Survey (SDSS)). Our best estimate for the flux at SK is slightly below, but very close to the current SK upper limit. The SK upper limit is already inconsistent with a range of star formation histories allowed by the SDSS data. We estimate that the SRN background should be detected (at 1-sigma) at SK with a total of about 9 years (including the existing 4 years) of data. While KamLAND is a much smaller detector compared to SK, it profits from being practically background-free and from its sensitivity to the lower energy supernova neutrinos. KamLAND could make a 1-sigma detection of the SRN with a total of about 5 years of data. Given the small expected SRN event rate, we also consider the detection of the SRN in a modified SK detector with a lower threshold and reduced background where the time to detection can be reduced by a factor of 10 relative to the existing SK estimate.

A common property of popular models of quintessence dark energy is the convergence to a common solution from a large range of the initial conditions. We re-examine the popular inverse power-law model of quintessence (where the common solution is dubbed as the 'tracker') with particular attention to the initial conditions for the field and their influence on the evolution. We find that previously derived limits on the parameters of the potential in this model are valid only in a range of initial conditions. A reasonably sharp boundary lies where the initial energy density of the scalar field is equal to that of the background radiation component. An initial quintessence energy density above this equipartition value lead to a solution that will not have joined the tracker solution by the present epoch. These non-tracker solutions possess the property that their present equation of state is very compatible with the observed bounds and independent of the exponent of the potential.