Radio-loud active galactic nuclei (AGN) are expected to contribute substantially to both the intensity and anisotropy of the isotropic gamma-ray background (IGRB). In turn, the measured properties of the IGRB can be used to constrain the characteristics of proposed contributing source classes. We consider individual subclasses of radio-loud AGN, including low-, intermediate-, and high-synchrotron-peaked BL Lacertae objects, flat-spectrum radio quasars, and misaligned AGN. Using updated models of the gamma-ray luminosity functions of these populations, we evaluate the energy-dependent contribution of each source class to the intensity and anisotropy of the IGRB. We find that collectively radio-loud AGN can account for the entirety of the IGRB intensity and anisotropy as measured by the Fermi Large Area Telescope (LAT). Misaligned AGN provide the bulk of the measured intensity but a negligible contribution to the anisotropy, while high-synchrotron-peaked BL Lacertae objects provide the dominant contribution to the anisotropy. In anticipation of upcoming measurements with the Fermi-LAT and the forthcoming Cherenkov Telescope Array, we predict the anisotropy in the broader energy range that will be accessible to future observations.

(abridged) The Galactic Center is one of the most promising targets for indirect detection of dark matter with gamma rays. We investigate the sensitivity of the upcoming Cherenkov Telescope Array (CTA) to dark matter annihilation and decay in the Galactic Center. As the inner density profile of the Milky Way's dark matter halo is uncertain, we study the impact of the slope of the Galactic density profile, inwards of the Sun, on the prospects for detecting a dark matter signal with CTA. We find that the sensitivity achieved by CTA to annihilation signals is strongly dependent on the inner profile slope, whereas the dependence is more mild in the case of dark matter decay. Surprisingly, we find that the optimal choice of signal and background regions is virtually independent of the assumed density profile. For the fiducial case of a Navarro-Frenk-White profile, we find that CTA will be able to probe annihilation cross sections well below the canonical thermal relic value for dark matter masses from a few tens of GeV up to $\sim 5$ TeV for annihilation to $\tau^{+}\tau^{-}$, and will achieve only a slightly weaker sensitivity for annihilation to $b\bar{b}$ or $\mu^{+}\mu^{-}$. CTA will improve significantly on current sensitivity to annihilation signals for dark matter masses above $\sim 100$ GeV, covering parameter space that is complementary to that probed by searches with the Fermi Large Area Telescope. The interpretation of apparent excesses in the measured cosmic-ray electron and positron spectra as signals of dark matter decay will also be testable with CTA observations of the Galactic Center. We demonstrate that both for annihilation and for decay, including spectral information for hard channels (such as $\mu^{+}\mu^{-}$ and $\tau^{+}\tau^{-}$) leads to enhanced sensitivity for dark matter masses above $m_{\rm DM}\sim 200$ GeV.

Indirect searches for products of dark matter annihilation and decay face the challenge of identifying an uncertain and subdominant signal in the presence of uncertain backgrounds. Two valuable approaches to this problem are (1) using analysis methods which take advantage of different features in the energy spectrum and angular distribution of the signal and backgrounds, and (2) more accurate characterization of backgrounds, which allows for more robust identification of possible signals. These two approaches are complementary and can be significantly strengthened when used together. I review the status of indirect searches with gamma rays using two promising targets, the Inner Galaxy and the Isotropic Gamma-Ray Background. For both targets, uncertainties in the properties of backgrounds is a major limitation to the sensitivity of indirect searches. I then highlight approaches which can enhance the sensitivity of indirect searches using these targets.

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

Dark matter particles captured by the Sun through scattering may annihilate and produce neutrinos, which escape. Current searches are for the few high-energy neutrinos produced in the prompt decays of some final states. We show that interactions in the solar medium lead to a large number of pions for nearly all final states. Positive pions and muons decay at rest, producing low-energy neutrinos with known spectra, including nuebar through neutrino mixing. We demonstrate that Super-Kamiokande can thereby provide a new probe of the spin-dependent WIMP-proton cross section. Compared to other methods, the sensitivity is competitive and the uncertainties are complementary.

Using data from the Fermi Gamma-Ray Space Telescope, a spatially extended component of gamma rays has been identified from the direction of the Galactic Center, peaking at energies of ~2-3 GeV. More recently, it has been shown that this signal is not confined to the innermost hundreds of parsecs of the Galaxy, but instead extends to at least ~3 kpc from the Galactic Center. While the spectrum, intensity, and angular distribution of this signal is in good agreement with predictions from annihilating dark matter, it has also been suggested that a population of unresolved millisecond pulsars could be responsible for this excess GeV emission from the Inner Galaxy. In this paper, we consider this later possibility in detail. Comparing the observed spectral shape of the Inner Galaxy's GeV excess to the spectrum measured from 37 millisecond pulsars by Fermi, we find that these sources exhibit a spectral shape that is much too soft at sub-GeV energies to accommodate this signal. We also construct population models to describe the spatial distribution and luminosity function of the Milky Way's millisecond pulsars. After taking into account constraints from the observed distribution of Fermi sources (including both sources known to be millisecond pulsars, and unidentified sources which could be pulsars), we find that millisecond pulsars can account for no more than ~10% of the Inner Galaxy's GeV excess. Each of these arguments strongly disfavor millisecond pulsars as the source of this signal.

The total anisotropy of a diffuse background composed of two or more sources, such as the Fermi-LAT--measured gamma-ray background, is set by the anisotropy of each source population and the contribution of each population to the total intensity. The total anisotropy as a function of energy (the anisotropy energy spectrum) will modulate as the relative contributions of the sources change, implying that the anisotropy energy spectrum also encodes the intensity spectrum of each source class. We develop techniques, applicable to any such diffuse background, for unraveling the intensity spectrum of each component source population given a measurement of the total intensity spectrum and the total anisotropy energy spectrum, without introducing \emph{a priori} assumptions about the spectra of the source classes. We demonstrate the potential of these methods by applying them to example scenarios for the composition of the Fermi-LAT gamma-ray background consistent with current data and feasible within 10 years of observation.

The detailed origin of the diffuse gamma-ray background is still unknown. However, the contribution of unresolved sources is expected to induce small-scale anisotropies in this emission, which may provide a way to identify and constrain the properties of its contributors. Recent studies have predicted the contributions to the angular power spectrum (APS) from extragalactic and galactic dark matter (DM) annihilation or decay. The Fermi-LAT collaboration reported detection of angular power with a significance larger than $3\sigma$ in the energy range from 1 GeV to 10 GeV on 22 months of data [Ackermann et al. 2012]. For these preliminary results the already published Fermi-LAT APS measurements [Ackermann et al. 2012] are compared to the accurate predictions for DM anisotropies from state-of-the-art cosmological simulations as presented in [Fornasa et al. 2013] to derive constraints on different DM candidates.

The Fermi-LAT collaboration has recently reported the detection of angular power above the photon noise level in the diffuse gamma-ray background between 1 and 50 GeV. Such signal can be used to constrain a possible contribution from Dark-Matter-induced photons. We estimate the intensity and features of the angular power spectrum (APS) of this potential Dark Matter (DM) signal, for both decaying and annihilating DM candidates, by constructing template all-sky gamma-ray maps for the emission produced in the galactic halo and its substructures, as well as in extragalactic (sub)halos. The DM distribution is given by state-of-the-art N-body simulations of cosmic structure formation, namely Millennium-II for extragalactic (sub)halos, and Aquarius for the galactic halo and its subhalos. We use a hybrid method of extrapolation to account for (sub)structures that are below the resolution limit of the simulations, allowing us to estimate the total emission all the way down to the minimal self-bound halo mass. We describe in detail the features appearing in the APS of our template maps and we estimate the effect of various uncertainties such as the value of the minimal halo mass, the fraction of substructures hosted in a halo and the shape of the DM density profile. Our results indicate that the fluctuation APS of the DM-induced emission is of the same order as the Fermi-LAT APS, suggesting that one can constrain this hypothetical emission from the comparison with the measured anisotropy. We also quantify the uncertainties affecting our results, finding "theoretical error bands" spanning more than two orders of magnitude and dominated (for a given particle physics model) by our lack of knowledge of the abundance of low-mass (sub)halos.

We place new constraints on the contribution of blazars to the large-scale isotropic gamma-ray background (IGRB) by jointly analyzing the measured source count distribution (logN-logS) of blazars and the measured intensity and anisotropy of the IGRB. We find that these measurements point to a consistent scenario in which unresolved blazars make less than 20% of the IGRB intensity at 1-10 GeV while accounting for the majority of the measured anisotropy in that energy band. These results indicate that the remaining fraction of the IGRB intensity is made by a component with a low level of intrinsic anisotropy. We determine upper limits on the anisotropy from non-blazar sources, adopting the best-fit parameters of the measured source count distribution to calculate the unresolved blazar anisotropy. In addition, we show that the anisotropy measurement excludes some recently proposed models of the unresolved blazar population.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

The small angular scale fluctuations of the (on large scale) isotropic gamma-ray background (IGRB) carry information about the presence of unresolved source classes. A guaranteed contribution to the IGRB is expected from the unresolved gamma-ray AGN while other extragalactic sources, Galactic gamma-ray source populations and dark matter Galactic and extragalactic structures (and sub-structures) are candidate contributors. The IGRB was measured with unprecedented precision by the Large Area Telescope (LAT) on-board of the Fermi gamma-ray observatory, and these data were used for measuring the IGRB angular power spectrum (APS). Detailed Monte Carlo simulations of Fermi-LAT all-sky observations were performed to provide a reference against which to compare the results obtained for the real data set. The Monte Carlo simulations are also a method for performing those detailed studies of the APS contributions of single source populations, which are required in order to identify the actual IGRB contributors. We present preliminary results of an anisotropy search in the IGRB. At angular scales <2deg (e.g. above multipole 155), angular power above the photon noise level is detected, at energies between 1 and 10 GeV in each energy bin, with statistical significance between 7.2 and 4.1 sigmas. The energy not dependence of the fluctuation anisotropy is pointing to the presence of one or more unclustered source populations, while the energy dependence of the intensity anisotropy is consistent with source populations having average photon index 2.40\pm0.07.

Pulsars emerge in the Fermi era as a sizable population of gamma-ray sources. Millisecond pulsars (MSPs) constitute an older subpopulation whose sky distribution extends to high Galactic latitudes, and it has been suggested that unresolved members of this class may contribute a significant fraction of the measured large-scale isotropic gamma-ray background (IGRB). We investigate the possible energy-dependent contribution of unresolved MSPs to the anisotropy of the Fermi-measured IGRB. For observationally-motivated MSP population models, we show that the preliminary Fermi anisotropy measurement places an interesting constraint on the abundance of MSPs in the Galaxy and the typical MSP flux, about an order of magnitude stronger than constraints on this population derived from the intensity of the IGRB alone. We also examine the possibility of a MSP component in the IGRB mimicking a dark matter signal in anisotropy-based searches, and conclude that the energy dependence of an anisotropy signature would distinguish MSPs from all but very light dark matter candidates.

Although most proposed dark matter candidates are stable, in order for dark matter to be present today, the only requirement is that its lifetime is longer than the age of the Universe, t_U ~ 4 10^17 s. Moreover, the dark matter particle could be produced via non-thermal processes and have a larger annihilation cross section from the canonical value for thermal dark matter, <sigma v> ~ 3 10^{-26} cm3/s. We propose a strategy to distinguish between dark matter annihilation and/or decay in the case that a clear signal is detected in future gamma-ray observations of Milky Way dwarf galaxies with gamma-ray experiments. The discrimination between these cases would not be possible in the case of the measurement of only the energy spectrum. We show that by studying the dependence of the intensity and energy spectrum on the angular distribution of the signal, the origin of the signal could be identified, and some information about the presence of substructure might be extracted.

The contribution of unresolved sources to the diffuse gamma-ray background could produce anisotropies in this emission on small angular scales. Recent studies have considered the angular power spectrum and other anisotropy metrics as tools for identifying contributions to diffuse emission from unresolved source classes, such as extragalactic and Galactic dark matter as well as various astrophysical gamma-ray source populations. We present preliminary results of an anisotropy analysis of the diffuse emission measured by the Fermi-LAT.

The energy-dependence of the anisotropy (the anisotropy energy spectrum) of the large-scale diffuse gamma-ray background can reveal the presence of multiple source populations. Annihilating dark matter in the substructure of the Milky Way halo could give rise to a modulation in the anisotropy energy spectrum of the diffuse gamma-ray emission measured by Fermi, enabling the detection of a dark matter signal. We determine the detectability of a dark-matter-induced modulation for scenarios in which unresolved blazars are the primary contributor to the measured emission above ~1 GeV and find that in some scenarios pair-annihilation cross sections of order the value expected for thermal relic dark matter can produce a detectable feature. We anticipate that the sensitivity of this technique to specific dark matter models could be improved by tailored likelihood analysis methods.

Most proposed dark matter candidates are stable and are produced thermally in the early Universe. However, there is also the possibility of unstable (but long-lived) dark matter, produced thermally or otherwise. We propose a strategy to distinguish between dark matter annihilation and/or decay in the case that a clear signal is detected in gamma-ray observations of Milky Way dwarf spheroidal galaxies with gamma-ray experiments. The sole measurement of the energy spectrum of an indirect signal would render the discrimination between these cases impossible. We show that by examining the dependence of the intensity and energy spectrum on the angular distribution of the emission, the origin could be identified as decay, annihilation, or both. In addition, once the type of signal is established, we show how these measurements could help to extract information about the dark matter properties, including mass, annihilation cross section, lifetime, dominant annihilation and decay channels, and the presence of substructure. Although an application of the approach presented here would likely be feasible with current experiments only for very optimistic dark matter scenarios, the improved sensitivity of upcoming experiments could enable this technique to be used to study a wider range of dark matter models.

Due to their higher concentrations and small internal velocities, Milky Way subhalos can be at least as important as the smooth halo in accounting for the GeV positron excess via dark matter annihilation. After showing how this can be achieved in various scenarios, including in Sommerfeld models, we demonstrate that, in this case, the diffuse inverse-Compton emission resulting from electrons and positrons produced in substructure leads to a nearly-isotropic signal close to the level of the isotropic GeV gamma-ray background seen by Fermi. Moreover, we show that HESS cosmic-ray electron measurements can be used to constrain multi-TeV internal bremsstrahlung gamma rays arising from annihilation to charged leptons.

Dark matter annihilation in Galactic substructure produces diffuse gamma-ray emission of remarkably constant intensity across the sky, and in general this signal dominates over the smooth halo signal at angles greater than a few tens of degrees from the Galactic Center. The large-scale isotropy of the emission from substructure suggests that it may be difficult to extract this Galactic dark matter signal from the extragalactic gamma-ray background. I show that dark matter substructure induces characteristic small-scale anisotropies in the diffuse emission which may provide a robust means of distinguishing this component. I discuss predictions for the angular power spectrum of the diffuse emission from various extragalactic source classes as well as from Galactic dark matter, and show that the energy dependence of the angular power spectrum of the total measured emission could be used to confidently identify a gamma-ray signal from Galactic dark matter.

Dark matter annihilation in Galactic substructure will produce diffuse gamma-ray emission of remarkably constant intensity across the sky, making it difficult to disentangle this Galactic dark matter signal from the extragalactic gamma-ray background. We show that if Galactic dark matter contributes a modest fraction of the measured emission in an energy range accessible to the Fermi Gamma-ray Space Telescope, the energy dependence of the angular power spectrum of the total measured emission could be used to confidently identify gamma rays from Galactic dark matter substructure.

The identification of point sources poses a great challenge for the high energy community. We present a new approach to evaluate the likelihood of a set of sources being a Galactic population based on the simple assumption that galaxies similar to the Milky Way host comparable populations of gamma-ray emitters. We propose a luminosity constraint on Galactic source populations which complements existing approaches by constraining the abundance and spatial distribution of any objects of Galactic origin, rather than focusing on the properties of a specific candidate emitter. We use M31 as a proxy for the Milky Way, and demonstrate this technique by applying it to the unidentified EGRET sources. We find that it is highly improbable that the majority of the unidentified EGRET sources are members of a Galactic halo population (e.g., dark matter subhalos), but that current observations do not provide any constraints on all of these sources being Galactic objects if they reside entirely in the disk and bulge. Applying this method to upcoming observations by the Fermi Gamma-ray Space Telescope has the potential to exclude association of an even larger number of unidentified sources with any Galactic source class.

The majority of gamma-ray emission from Galactic dark matter annihilation is likely to be detected as a contribution to the diffuse gamma-ray background. I show that dark matter substructure in the halo of the Galaxy induces characteristic anisotropies in the diffuse background that could be used to determine the small-scale dark matter distribution. I calculate the angular power spectrum of the emission from dark matter substructure for several models of the subhalo population, and show that features in the power spectrum can be used to infer the presence of substructure. The shape of the power spectrum is largely unaffected by the subhalo radial distribution and mass function, and for many scenarios I find that a measurement of the angular power spectrum by Fermi will be able to constrain the abundance of substructure. An anti-biased subhalo radial distribution is shown to produce emission that differs significantly in intensity and large-scale angular dependence from that of a subhalo distribution which traces the smooth dark matter halo, potentially impacting the detectability of the dark matter signal for a variety of targets and methods.

In the past decade, surveys of the stellar component of the Galaxy have revealed a number of streams from tidally disrupted dwarf galaxies and globular clusters. Simulations of hierarchical structure formation in LCDM cosmologies predict that the dark matter halo of a galaxy like the Milky Way contains hundreds of subhalos with masses of ~10^8 solar masses and greater, and it has been suggested that the existence of coherent tidal streams is incompatible with the expected abundance of substructure. We investigate the effects of dark matter substructure on tidal streams by simulating the disruption of a self-gravitating satellite on a wide range of orbits in different host models both with and without substructure. We find that the halo shape and the specific orbital path more strongly determine the overall degree of disruption of the satellite than does the presence or absence of substructure, i.e., the changes in the large-scale properties of the tidal debris due to substructure are small compared to variations in the debris from different orbits in a smooth potential. Substructure typically leads to an increase in the degree of clumpiness of the tidal debris in sky projection, and in some cases a more compact distribution in line-of-sight velocity. Substructure also leads to differences in the location of sections of debris compared to the results of the smooth halo model, which may have important implications for the interpretation of observed tidal streams. A unique signature of the presence of substructure in the halo which may be detectable by upcoming surveys is identified. We conclude, however, that predicted levels of substructure are consistent with a detection of a coherent tidal stream from a dwarf galaxy.

The large majority of EGRET point sources remain without an identified low-energy counterpart, and a large fraction of these sources are most likely extragalactic. Whatever the nature of the extragalactic EGRET unidentified sources, faint unresolved objects of the same class must have a contribution to the diffuse extragalactic gamma-ray background (EGRB). Understanding this component of the EGRB, along with other guaranteed contributions from known sources, is essential if we are to use this emission to constrain exotic high-energy physics. Here, we follow an empirical approach to estimate whether a potential contribution of unidentified sources to the EGRB is likely to be important, and we find that it is. Additionally, we show how upcoming GLAST observations of EGRET unidentified sources, as well as of their fainter counterparts, can be combined with GLAST observations of the Galactic and extragalactic diffuse backgrounds to shed light on the nature of the EGRET unidentified sources even without any positional association of such sources with low-energy counterparts.

The third EGRET catalog contains a large number of unidentified sources. This subset of objects is expected to include known gamma-ray emitters of Galactic origin such as pulsars and supernova remnants, in addition to an extragalactic population of blazars. However, current data allows the intriguing possibility that some of these objects may represent a new class of yet undiscovered gamma-ray sources. Many theoretically motivated candidate emitters (e.g. clumps of annihilating dark matter particles) have been suggested to account for these detections. We take a new approach to determine to what extent this population is Galactic and to investigate the nature of the possible Galactic component. By assuming that galaxies similar to the Milky Way should host comparable populations of objects, we constrain the allowed Galactic abundance and distribution of various classes of gamma-ray sources using the EGRET data set. We find it is highly improbable that a large number of the unidentified sources are members of a Galactic halo population, but that a distribution of the sources entirely in the disk and bulge is plausible. Finally, we discuss the additional constraints and new insights that GLAST will provide.

The large majority of EGRET point sources remain to this day without an identified low-energy counterpart. Whatever the nature of the EGRET unidentified sources, faint unresolved objects of the same class must have a contribution to the diffuse gamma-ray background: if most unidentified objects are extragalactic, faint unresolved sources of the same class contribute to the background, as a distinct extragalactic population; on the other hand, if most unidentified sources are Galactic, their counterparts in external galaxies will contribute to the unresolved emission from these systems. Understanding this component of the gamma-ray background, along with other guaranteed contributions from known sources, is essential in any attempt to use gamma-ray observations to constrain exotic high-energy physics. Here, we follow an empirical approach to estimate whether a potential contribution of unidentified sources to the extragalactic gamma-ray background is likely to be important, and we find that it is. Additionally, we comment on how the anticipated GLAST measurement of the diffuse gamma-ray background will change, depending on the nature of the majority of these sources.