On June 20, 2020, the LIGO-Virgo collaboration announced the discovery of GW190814, a gravitational wave event originating from a binary system merger between a black hole of mass $M_1 = 23.2^{+1.1} _ {-1.0}M_\odot$ and an unidentified object with a mass of $M_2 = 2.59^{+0.08} _ {-0.09}M_\odot$. This second object would be either the heaviest neutron star or lightest black hole observed to date. Here we investigate the possibility of the $\sim 2.6M_\odot$ object being a primordial black hole (PBH). We find that a primordial black hole explanation to GW190814 is unlikely as it is limited by the formation rate of the primary stellar progenitor and the observed merger rates of $\mathcal{O}(20)M_\odot$ massive black hole pairs.

It has been suggested that late-universe dark matter decays can alleviate the tension between measurements of $H_0$ in the local universe and its value inferred from cosmic microwave background fluctuations. Decaying dark matter can potentially account for this discrepancy as it reshuffles the energy density between matter and radiation and as a result allows dark energy to become dominant at earlier times. We show that the low multipoles amplitude of the cosmic microwave background anisotropy power spectrum severely constrains the feasibility of late-time decays as a solution to the $H_0$ tension.

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

We study the cosmological effects of two-body dark matter decays where the products of the decay include a massless and a massive particle. We show that if the massive daughter particle is slightly warm it is possible to relieve the tension between distance ladder measurements of the present day Hubble parameter with measurements from the cosmic microwave background.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

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 present methods to assess whether gamma-ray excesses towards Milky Way dwarf galaxies can be attributed to astrophysical sources rather than to dark matter annihilation. As a case study we focus on Reticulum II, the dwarf which shows the strongest evidence for a gamma-ray signal in Fermi data. Dark matter models and those with curved energy spectra provide good fits to the data, while a simple power law is ruled out at 97.5% confidence. We compare RetII's spectrum to known classes of gamma-ray sources and find a useful representation in terms of spectral curvature and the energy at which the spectral energy distribution peaks. In this space the blazar classes appear segregated from the confidence region occupied by RetII. Pulsars have similar gamma-ray spectra to RetII but we show that RetII is unlikely to host a pulsar population detectable in gamma rays. Tensions with astrophysical explanations are stronger when analyzing 6.5 years of Pass 7 than with the same amount of Pass 8 data, where the excess is less significant. These methods are applicable to any dwarf galaxy which is a promising dark matter target and shows signs of gamma-ray emission along its line of sight.

A thick dark matter disk is predicted in cold dark matter simulations as the outcome of the interaction between accreted satellites and the stellar disk in Milky Way sized halos. We study the effects of a self-interacting thick dark disk on the energetic neutrino flux from the Sun. We find that for particle masses between 100 GeV and 1 TeV and dark matter annihilation to heavy leptons either the self-interaction may not be strong enough to solve the small scale structure motivation or a dark disk cannot be present in the Milky Way.

Future generation of gravitational wave detectors will have the sensitivity to detect gravitational wave events at redshifts far beyond any detectable electromagnetic sources. We show that if the observed event rate is greater than one event per year at redshifts z > 40, then the probability distribution of primordial density fluctuations must be significantly non-Gaussian or the events originate from primordial black holes. The nature of the excess events can be determined from the redshift distribution of the merger rate.

We study the effects of black hole dark matter on the dynamical evolution of stars in dwarf galaxies. We find that mass segregation leads to a depletion of stars in the center of dwarf galaxies and the appearance of a ring in the projected stellar surface density profile. Using Segue 1 as an example we show that current observations of the projected surface stellar density rule out at the 99.9% confidence level the possibility that more than 6% of the dark matter is composed of black holes with a mass of few tens of solar masses.

We present constraints on the annihilation cross section of WIMP dark matter based on the joint statistical analysis of four dwarf galaxies with VERITAS. These results are derived from an optimized photon weighting statistical technique that improves on standard imaging atmospheric Cherenkov telescope (IACT) analyses by utilizing the spectral and spatial properties of individual photon events. We report on the results of $\sim$230 hours of observations of five dwarf galaxies and the joint statistical analysis of four of the dwarf galaxies. We find no evidence of gamma-ray emission from any individual dwarf nor in the joint analysis. The derived upper limit on the dark matter annihilation cross section from the joint analysis is $1.35\times 10^{-23} {\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the bottom quark ($b\bar{b}$) final state, $2.85\times 10^{-24}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the tau lepton ($\tau^{+}\tau^{-}$) final state and $1.32\times 10^{-25}{\mathrm{ cm^3s^{-1}}}$ at 1 TeV for the gauge boson ($\gamma\gamma$) final state.

We present a study of cosmological implications of generic dark matter decays. We consider two-body and many-body decaying scenarios. In the two-body case the massive particle has a possibly relativistic kick velocity and thus possesses a dynamical equation of state. This has implications to the expansion history of the universe. We use recent observational data from the cosmic microwave background, baryon acoustic oscillations and supernovae Type Ia to obtain constraints on the lifetime of the dark matter particle. We find that for an energy splitting where more than 40% of the dark matter particle energy is transferred to massless, relativistic particles in the two-body case, or more than 50% in the many-body case, lifetimes less than the age of the universe are excluded at more than 95% confidence. When the energy splitting falls to 10% the lifetime is constrained to be more than roughly half the age.

The dwarf spheroidal galaxies (dSph) of the Milky Way are among the most attractive targets for indirect searches of dark matter. In this work, we reconstruct the dark matter annihilation (J-factor) and decay profiles for the newly discovered dSph Reticulum II. Using an optimized spherical Jeans analysis of kinematic data obtained from the Michigan/Magellan Fiber System (M2FS), we find Reticulum II's J-factor to be among the largest of any Milky Way dSph. We have checked the robustness of this result against several ingredients of the analysis. Unless it suffers from tidal disruption or significant inflation of its velocity dispersion from binary stars, Reticulum II may provide a unique window on dark matter particle properties.

We present a search for gamma-ray emission from the direction of the newly discovered dwarf galaxy Reticulum II. Using Fermi-LAT data, we detect a signal that exceeds expected backgrounds between ~2-10 GeV and is consistent with annihilation of dark matter for particle masses less than a few x 10^2 GeV. Modeling the background as a Poisson process based on Fermi-LAT diffuse models, and taking into account trials factors, we detect emission with p-value less than 9.8 x 10^-5 (>3.7 sigma). An alternative, model-independent treatment of background reduces the significance, raising the p-value to 9.7 x 10^-3 (2.3 sigma). Even in this case, however, Reticulum II has the most significant gamma-ray signal of any known dwarf galaxy. If Reticulum II has a dark matter halo that is similar to those inferred for other nearby dwarfs, the signal is consistent with the s-wave relic abundance cross section for annihilation.

We present a new formalism designed to discover dark matter annihilation occurring in the Milky Way's dwarf galaxies. The statistical framework extracts all available information in the data by simultaneously combining observations of all the dwarf galaxies and incorporating the impact of particle physics properties, the distribution of dark matter in the dwarfs, and the detector response. The method performs maximally powerful frequentist searches and produces confidence limits on particle physics parameters. Probability distributions of test statistics under various hypotheses are constructed exactly, without relying on large sample approximations. The derived limits have proper coverage by construction and claims of detection are not biased by imperfect background modeling. We implement this formalism using data from the Fermi Gamma-ray Space Telescope to search for an annihilation signal in the complete sample of Milky Way dwarfs whose dark matter distributions can be reliably determined. We find that the observed data is consistent with background for each of the dwarf galaxies individually as well as in a joint analysis. The strongest constraints are at small dark matter particle masses. Taking the median of the systematic uncertainty in dwarf density profiles, the cross section upper limits are below the pure s-wave weak scale relic abundance value (2.2 x 10^-26 cm^3/s) for dark matter masses below 26 GeV (for annihilation into b quarks), 29 GeV (tau leptons), 35 GeV (up, down, strange, and charm quarks and gluons), 6 GeV (electrons/positrons), and 114 GeV (two-photon final state). For dark matter particle masses less than 1 TeV, these represent the strongest limits obtained to date using dwarf galaxies.

Gamma-ray searches for dark matter annihilation and decay in dwarf galaxies rely on an understanding of the dark matter density profiles of these systems. Conversely, uncertainties in these density profiles propagate into the derived particle physics limits as systematic errors. In this paper we quantify the expected dark matter signal from 20 Milky Way dwarfs using a uniform analysis of the most recent stellar-kinematic data available. Assuming that the observed stellar populations are equilibrium tracers of spherically-symmetric gravitational potentials that are dominated by dark matter, we find that current stellar-kinematic data can predict the amplitudes of annihilation signals to within a factor of a few for the ultra-faint dwarfs of greatest interest. On the other hand, the expected signal from several classical dwarfs (with high-quality observations of large numbers of member stars) can be localized to the ~20% level. These results are important for designing maximally sensitive searches in current and future experiments using space and ground-based instruments.

We present a decaying dark matter scenario where the daughter products are a single massless relativistic particle and a single, massive but possibly relativistic particle. We calculate the velocity distribution of the massive daughter particle and its associated equation of state and derive its dynamical evolution in an expanding universe. In addition, we present a model of decaying dark matter where there are many massless relativistic daughter particles together with a massive particle at rest. We place constraints on these two models using supernovae type Ia observations. We find that for a daughter relativistic fraction of 1% and higher, lifetimes of at least less than 10 Gyrs are excluded, with larger relativistic fractions constraining longer lifetimes.

In this white paper, submitted as a part of Snowmass 2013 (subgroup CF2), we examine the current status and future prospects of the VERITAS indirect dark matter detection program. The VERITAS array of imaging atmospheric Cherenkov telescopes (IACTs), sensitive in the 0.1-50 TeV regime, is in the process of completing a multi-year program aimed at detecting signatures of neutralino dark matter. This program is spread out over a range of astrophysical targets which can potentially yield definitive signatures of neutralino self-annihilation such as dwarf spheroidal galaxies (dSphs) and the center of the Milky Way galaxy. While the program is still in progress, initial results on dSphs have produced very competitive upper limits on the thermally averaged cross-section of neutralino self-annihilation as well as strongly constraining leptophillic dark matter models (such as those inferred from the PAMELA and AMS positron fraction results). As the program continues over the next 5 years, VERITAS observations will yield some of the strongest constraints available from IACTs, pushing down the limit on the thermally averaged cross-section of neutralino self-annihilation to the expected natural thermal relic scale.

We explore the cosmological halo-to-halo scatter of the distribution of mass within dark matter halos utilizing a well-resolved statistical sample of clusters from the cosmological Millennium simulation. We find that at any radius, the spherically-averaged dark matter density of a halo (corresponding to the "smooth-component") and its logarithmic slope are well-described by a Gaussian probability distribution. At small radii (within the scale radius), the density distribution is fully determined by the measured Gaussian distribution in halo concentrations. The variance in the radial distribution of mass in dark matter halos is important for the interpretation of direct and indirect dark matter detection efforts. The scatter in mass profiles imparts approximately a 25 percent cosmological uncertainty in the dark matter density at the Solar neighborhood and a factor of ~3 uncertainty in the expected Galactic dark matter annihilation flux. The aggregate effect of halo-to-halo profile scatter leads to a small (few percent) enhancement in dark matter annihilation background if the Gaussian concentration distribution holds for all halo masses versus a 10 percent enhancement under the assumption of a log-normal concentration distribution. The Gaussian nature of the cluster profile scatter implies that the technique of "stacking" halos to improve signal to noise should not suffer from bias.

We perform a joint analysis of dwarf galaxy data from the Fermi Gamma-ray Space Telescope in search of dark matter annihilation into a gamma-ray line. We employ a novel statistical method that takes into account the spatial and spectral information of individual photon events from a sample of seven dwarf galaxies. Dwarf galaxies show no evidence of a gamma-ray line between 10 GeV and 1 TeV. The subsequent upper limit on the annihilation cross section to a two-photon final state is 3.9(+7.1)(-3.7) x 10^-26 cm^3/s at 130 GeV, where the errors reflect the systematic uncertainty in the distribution of dark matter within the dwarf galaxies.

We present a statistical technique which can be used to detect the presence and properties of moving sources contributing to a diffuse background. The method is a generalization of the 2-point correlation function to include temporal as well as spatial information. We develop a formalism which allows for a derivation of the spacetime 2-point function in terms of the properties of the contributing sources. We test this technique in simulated sky maps, and demonstrate its robustness in identifying the presence of moving and stationary sources. Applications of this formalism to the diffuse gamma-ray background include searches for solar system bodies, fast moving primordial black holes, and dense cores of dark matter proto-halos in the solar neighborhood. Other applications include detecting the contribution of energetic neutrinos originating in the solar system, as well as probing compact objects in long timeline lensing experiments.

Dwarf spheroidal galaxies are known to be excellent targets for the detection of annihilating dark matter. We present new limits on the annihilation cross section of Weakly Interacting Massive Particles (WIMPs) based on the joint analysis of seven Milky Way dwarfs using a frequentist Neyman construction and Pass 7 data from the Fermi Gamma-ray Space Telescope. We exclude generic WIMP candidates annihilating into b-bbar with mass less than 40 GeV that reproduce the observed relic abundance. To within 95% systematic errors on the dark matter distribution within the dwarfs, the mass lower limit can be as low as 19 GeV or as high as 240 GeV. For annihilation into tau+tau- these limits become 19 GeV, 13 GeV, and 80 GeV respectively.

We explore the halo-to-halo variation of dark matter substructure in galaxy-sized dark matter halos, focusing on its implications for strongly gravitational lensed systems. We find that the median value for projected substructure mass fractions within projected radii of 3% of the host halo virial radius is approximately fsub ~ 0.25%, but that the variance is large with a 95-percentile range of 0 <= fsub <= 1%. We quantify possible effects of substructure on quadruply-imaged lens systems using the cusp relation and the simple statistic, Rcusp. We estimate that the probability of obtaining the large values of the Rcusp which have been observed from substructure effects is roughly ~10^-3 to ~10^-2. We consider a variety of possible correlations between host halo properties and substructure properties in order to probe possible sample biases. In particular, low-concentration host dark matter halos have more large substructures and give rise to large values of Rcusp more often. However, there is no known observational bias that would drive observed quadruply-imaged quasars to be produced by low-concentration lens halos. Finally, we show that the substructure mass fraction is a relatively reliable predictor of the value of Rcusp.

We present a statistical framework which can be used to determine the contribution of an unresolved population of pulsars to the gamma-ray background. This formalism is based on the joint analysis of photon time series over extended regions of the sky. We demonstrate the robustness of this technique in controlled simulations of pulsar populations, and show that the Fermi Gamma-ray Space Telescope can be used to detect a pulsar contribution as small as 0.1% of the gamma-ray background. This technique is sensitive to pulsar populations with photon fluxes greater than ~1E-10 cm^-2 s^-1. The framework is extensible to arbitrarily complex searches for periodicity and can therefore be tailored to specific applications such as all-sky surveys and studies of the Galactic center and globular clusters.

We calculate the probability distribution function (PDF) of the expected annihilation luminosities of dark matter subhalos as a function of subhalo mass and distance from the Galactic center using a semi-analytical model of halo evolution. We find that the PDF of luminosities is relatively broad, exhibiting a spread of as much as an order of magnitude at fixed subhalo mass and halo-centric distance. The luminosity PDF allows for simple construction of mock samples of gamma-ray luminous subhalos and assessment of the variance in among predicted gamma-ray signals from dark matter annihilation. Other applications include quantifying the variance among the expected luminosities of dwarf spheroidal galaxies, assessing the level at which dark matter annihilation can be a contaminant in the expected gamma-ray signal from other astrophysical sources, as well as estimating the level at which nearby subhalos can contribute to the antimatter flux.

We investigate the feasibility of extracting the gravitational nanolensing signal due to the presence of subsolar mass halos within galaxy-sized dark matter halos. We show that subsolar mass halos in a lensing galaxy can cause strong nanolensing events with shorter durations and smaller amplitudes than microlensing events caused by stars. We develop techniques that can be used in future surveys such as Pan-STARRS, LSST and OMEGA to search for the nanolensing signal from subsolar mass halos.

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 study the effects of substructure on the rate of dark-matter annihilation in the Galactic halo. We use an analytic model for substructure that can extend numerical simulation results to scales too small to be resolved by the simulations. We first calibrate the analytic model to numerical simulations, and then determine the annihilation boost factor, for standard WIMP models as well as those with Sommerfeld (or other) enhancements, as a function of Galactocentric radius in the Milky Way. We provide an estimate of the dependence of the gamma-ray intensity of WIMP annihilation as a function of angular distance from the Galactic center. This methodology, coupled with future numerical simulation results can be a powerful tool that can be used to constrain WIMP properties using Fermi all-sky data.

Dark matter halos contain a wealth of substructure in the form of subhalos and tidal streams. Enhancements in the dark matter density of these regions leads to enhanced rates in direct detection experiments, as well as enhanced dark matter capture rates in the Sun and the Earth. Direct detection experiments probe the present-day dark matter density, while energetic neutrinos probe the past history of the dark matter density along the solar system's orbit about the Galactic center. We discuss how an elevated energetic neutrino flux can be used to probe the level of substructure present at the Galactic radius of the solar system.

We consider the effects of Galactic substructure on energetic neutrinos from annihilation of weakly-interacting massive particles (WIMPs) that have been captured by the Sun and Earth. Substructure gives rise to a time-varying capture rate and thus to time variation in the annihilation rate and resulting energetic-neutrino flux. However, there may be a time lag between the capture and annihilation rates. The energetic-neutrino flux may then be determined by the density of dark matter in the Solar System's past trajectory, rather than the local density. The signature of such an effect may be sought in the ratio of the direct- to indirect-detection rates.

Dark matter halos of sub-solar mass are the first bound objects to form in cold dark matter theories. In this article, I discuss the present understanding of "microhalos'', their role in structure formation, and the implications of their potential presence, in the interpretation of dark matter experiments.

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.

This is a report on the findings of the dark matter science working group for the white paper on the status and future of TeV gamma-ray astronomy. The white paper was commissioned by the American Physical Society, and the full white paper can be found on astro-ph (arXiv:0810.0444). This detailed section discusses the prospects for dark matter detection with future gamma-ray experiments, and the complementarity of gamma-ray measurements with other indirect, direct or accelerator-based searches. We conclude that any comprehensive search for dark matter should include gamma-ray observations, both to identify the dark matter particle (through the charac- teristics of the gamma-ray spectrum) and to measure the distribution of dark matter in galactic halos.

One of the goals of NASA's Fermi Gamma-ray Space Telescope (formerly GLAST) will be detection of gamma rays from dark-matter annihilation in the Galactic halo. Theoretical arguments suggest that dark matter may be bound into subhalos with masses as small as 10^{-4}--10^2 Earth mass. If so, it may be possible to detect individual subhalos as point sources in the Fermi Telescope. It has further been argued that some of these point sources may exhibit proper motions. Here we show that upper limits to the diffuse gamma-ray background constrain the number of subhalos close enough to exhibit proper motions to be less than one.

We study the effects of substructure in the Galactic halo on direct detection of dark matter, on searches for energetic neutrinos from WIMP annihilation in the Sun and Earth, and on the enhancement in the WIMP annihilation rate in the halo. Our central result is a probability distribution function (PDF) P(\rho) for the local dark-matter density. This distribution must be taken into account when using null dark-matter searches to constrain the properties of dark-matter candidates. We take two approaches to calculating the PDF. The first is an analytic model that capitalizes on the scale-invariant nature of the structure--formation hierarchy in order to address early stages in the hierarchy (very small scales; high densities). Our second approach uses simulation-inspired results to describe the PDF that arises from lower-density larger-scale substructures which formed in more recent stages in the merger hierarchy. The distributions are skew positive, and they peak at densities lower than the mean density. The local dark-matter density may be as small as 1/10th the canonical value of ~ 0.4 GeV/cm^3, but it is probably no less than 0.2 GeV/cm^3.

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.

Notwithstanding the advent of the Gamma-ray Large Area Telescope, theoretical models predict that a significant fraction of the cosmic gamma-ray background (CGB), at the level of 20% of the currently measured value, will remain unresolved. The angular power spectrum of intensity fluctuations of the CGB contains information on its origin. We show that probing the latter from a few tens of arcmin to several degree scales, together with complementary GLAST observations of gamma-ray emission from galaxy clusters and the blazars luminosity function, can discriminate between a background that originates from unresolved blazars or cosmic rays accelerated at structure formation shocks.

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.

I discuss the prospects for detecting the dark matter via the proper motion of sub-solar mass dark matter halos in the vicinity of the solar neighbourhood. Microhalos that survive tidal disruption could exhibit proper motion of order few arcminutes per year. For dark matter particles that couple to photons, such as the lightest supersymmetric or Kaluza-Klein particles, microhalos could be detected via their Gamma-ray photon emission from annihilations. A detection of proper motion of a microhalo in the Gamma-ray part of the spectrum contains not only information about the particle physics properties of the dark matter particle, but also provides an insight into hierarchical structure formation at very early times.

We present upper limits on line emission in the Cosmic X-ray background (CXB) that would be produced by decay of sterile neutrino dark matter. We employ the spectra of the unresolved component of the CXB in the Chandra Deep Fields North and South obtained with the Chandra CCD detector in the E=0.8-9 keV band. The expected decay flux comes from the dark matter on the lines of sight through the Milky Way galactic halo. Our constraints on the sterile neutrino decay rate are sensitive to the modeling of the Milky Way halo. The highest halo mass estimates provide a limit on the sterile neutrino mass of m_s<2.9 keV in the Dodelson-Widrow production model, while the lowest halo mass estimates provide the conservative limit of m_s<5.7 keV (2-sigma). We also discuss constraints from a short observation of the softer (E<1 keV) X-ray background with a rocket-borne calorimeter by McCammon and collaborators.

I discuss the prospects of detecting the smallest dark matter bound structures present in the Milky Way by searching for the proper motion of $\gamma$-ray sources in the upcoming GLAST all sky map. I show that for dark matter particle candidates that couple to photons the detection of at least one $\gamma$-ray microhalo source with proper motion places a constraint on the couplings and mass of the dark matter particle. For SUSY dark matter, proper motion detection implies that the mass of the particle is less than 500 GeV and the kinetic decoupling temperature is in the range of [4-100] MeV.

We present a comprehensive analysis of constraints on the sterile neutrino as a dark matter candidate. The minimal production scenario with a standard thermal history and negligible cosmological lepton number is in conflict with conservative radiative decay constraints from the cosmic X-ray background in combination with stringent small-scale structure limits from the Lyman-alpha forest. We show that entropy release through massive particle decay after production does not alleviate these constraints. We further show that radiative decay constraints from local group dwarf galaxies are subject to large uncertainties in the dark matter density profile of these systems. Within the strongest set of constraints, resonant production of cold sterile neutrino dark matter in non-zero lepton number cosmologies remains allowed.

We study the gravitational wave background produced from the formation and assembly of supermassive black holes within the cosmological paradigm of hierarchical structure formation. In particular, we focus on a supermassive black hole formation scenario in which the present-day population of supermassive black holes is built from high-mass seed black holes and we compute the concomitant spectrum of gravitational radiation produced by mergers of the seed black holes. We find that this scenario predicts a large, gravitational wave background that should be resolved into individual sources with space interferometers such as the proposed Laser Interferometric Space Antenna (LISA). The number of inspiral, merger and ringdown events above a signal to noise ratio of 5 that result from massive black hole seeds is of order 10^3. This prediction is robust and insensitive to several of the details of the model. We conclude that an interferometer such as LISA will be able to effectively rule out or confirm a class of models where supermassive black holes grow from high-mass seed black holes and may be able to place strong limits on the role of mergers as a channel for supermassive black hole growth. Supermassive black hole seeds likely form in the earliest proto-galactic structures at high redshift and the masses of supermassive black holes are known to be strongly correlated with the potentials of the spheroids in which they reside, therefore these results imply that space interferometers can be used as a powerful probe of the physics of galaxy formation and proto-galaxy formation at very high redshift.

In this proceeding, we briefly review three recent results. First, we show that halos formed in simulations with gas cooling are significantly rounder than halos formed in dissipationless $N$-body simulations. The increase in principle axis ratios is $\delta (c/a) ~ 0.2 - 0.4$ in the inner halo and remains significant at large radii. Second, we discuss the CDM substructure crisis and demonstrate the sensitivity of the crisis to the spectrum of primordial density fluctuations on small scales. Third, we assess the ability of experiments like VERITAS and GLAST to detect $\gamma$-rays from neutralino dark matter annihilation in dark subhalos about the MW.

(ABRIDGED) We present a model in which the seeds of supermassive black holes form from the lowest angular momentum gas in proto-galaxies at high redshift. We show that this leads to a correlation between black hole masses and spheroid properties, as observed today. We assume that gas in early-forming, rare halos has a distribution of specific angular momentum similar to that expected for dark matter halos. This distribution has a significant low angular momentum tail which implies that every proto-galaxy should contain gas that ends up in a high-density disc. In halos more massive than a critical threshold of \~10^8Msun, the discs are gravitationally unstable, and experience an efficient Lin-Pringle viscosity that allows mass inflow. This process continues until the first massive stars disrupt the disc. Black holes are created with a characteristic mass of ~10^5 Msun, independent of the redshift of formation. This serves as a lower bound for black-hole masses in galaxies today. The comoving mass density of black hole seeds grows with time, and saturates when cosmic reionization stops gas cooling in these low-mass systems. By z~15, the comoving mass density becomes comparable to that inferred from observations, with room for appropriate additional luminous growth during a later quasar accretion phase. Hierarchical merging after z~15 naturally leads to a linear correlation between black-hole mass and stellar spheroid mass, with negligible black hole masses in disc-dominated galaxies.

We estimate the probability of detecting gamma-rays from the annihilation of neutralino dark matter in dense, central regions of Milky Way substructure. We characterize Galactic substructure statistically based on Monte Carlo realizations of the formation of a Milky Way-like halo using a semi-analytic method that has been calibrated against the results of high-resolution N-body simulations. We find that it may be possible for the upcoming experiments GLAST and VERITAS, working in concert, to detect gamma-rays from dark matter substructure if the neutralino is relatively light ($\Mchi \lsim 100 \GeV$), while for $\Mchi \gsim 500 \GeV$ such a detection would be unlikely. We perform most of our calculations within the framework of the standard $\Lambda$CDM cosmological model; however, we also investigate the robustness of our results to various assumptions and find that the probability of detection is sensitive to poorly-constrained input parameters, particularly those that characterize the primordial power spectrum. Specifically, the best-fitting power spectrum of the WMAP team, with a running spectral index, predicts roughly a factor of fifty fewer detectable subhalos compared to the standard $\LCDM$ cosmological model with scale-invariant power spectrum. We conclude that the lack of a detected gamma-ray signal gives very little information about the supersymmetric parameter space due to uncertainties associated with both the properties of substructure and cosmological parameters.

We estimate the probability of detecting gamma-rays from the annihilation of neutralino dark matter in the substructure of the Milky Way. We characterize substructure statistically based on Monte Carlo realizations of the formation of a Milky Way-like halo using semi-analytic methods that have been calibrated against N-body simulations. We find that it may be possible for the upcoming GLAST and VERITAS experiments, working in concert, to detect gamma-rays from dark matter substructure if $\Mchi \lsim 100 \GeV$, while for $\Mchi \gsim 500 \GeV$ such a detection seems unlikely. We investigate the effects of the underlying cosmological model and find that the probability of detection is sensitive to the primordial power spectrum of density fluctuations on small (galactic and sub-galactic) scales. We conclude that the lack of such a detection reveals little about the supersymmetric parameter space due to the uncertainties associated with the properties of substructure and cosmological parameters.