The International Gamma-Ray Astrophysics Laboratory (INTEGRAL) satellite has yielded unprecedented measurements of the soft gamma-ray spectrum of our Galaxy. Here we use those measurements to set constraints on dark matter (DM) that decays or annihilates into photons with energies $E\approx 0.02-2$ MeV. First, we revisit the constraints on particle DM that decays or annihilates to photon pairs. In particular, for decaying DM, we find that previous limits were overstated by roughly an order of magnitude. Our new, conservative analysis finds that the DM lifetime must satisfy $\tau\gtrsim 5\times 10^{26}\,{\rm s}\times (m_{\chi}/\rm MeV)^{-1}$ for DM masses $m_{\chi}=0.054-3.6$ MeV. For MeV-scale DM that annihilates into photons INTEGRAL sets the strongest constraints to date. Second, we target ultralight primordial black holes (PBHs) through their Hawking radiation. This makes them appear as decaying DM with a photon spectrum peaking at $E\approx 5.77/(8\pi G M_{\rm PBH})$, for a PBH of mass $M_{\rm PBH}$. We use the INTEGRAL data to demonstrate that, at 95\% C.L., PBHs with masses less than $1.2\times 10^{17}$ g cannot comprise all of the DM, setting the tightest bound to date on ultralight PBHs.

The Galactic Center GeV excess (GCE) has garnered great interest as a possible signal of either dark matter annihilation or some novel astrophysical phenomenon, such as a new population of gamma-ray emitting pulsars. In a companion paper, we showed that in a $10^\circ$ radius region of interest (ROI) surrounding the Galactic Center, apparent evidence for GCE point sources (PSs) from non-Poissonian template fitting (NPTF) is actually an artifact of unmodeled north-south asymmetry of the GCE. In this work, we develop a simplified analytic description of how signal mismodeling can drive an apparent preference for a PS population, and demonstrate how the behavior pointed out in the companion paper also appears in simpler simulated datasets that contain no PS signals at all. We explore the generality of this behavior in the real gamma-ray data, and discuss the implications for past and future studies using NPTF techniques. While the drop in PS preference once north-south asymmetry is included is not ubiquitous in larger ROIs, we show that any overly-rigid signal model is expected to yield a spurious PS signal that can appear very convincing: as well as apparent significance comparable to what one would expect from a true PS population, the signal can exhibit stability against a range of variations in the analysis, and a source count function that is very consistent with previous apparent NPTF-based detections of a GCE PS population. This contrasts with previously-studied forms of systematic mismodeling which are unlikely to mimic a PS population in the same way. In the light of this observation, and its explicit realization in the region where the GCE is brightest, we argue that a dominantly smooth origin for the GCE is not in tension with existing NPTF analyses.

We re-examine evidence that the Galactic Center Excess (GCE) originates primarily from point sources (PSs). We show that in our region of interest, non-Poissonian template fitting (NPTF) evidence for GCE PSs is an artifact of unmodeled north-south asymmetry of the GCE. This asymmetry is strongly favored by the fit (although it is unclear if this is physical), and when it is allowed, the preference for PSs becomes insignificant. We reproduce this behavior in simulations, including detailed properties of the spurious PS population. We conclude that NTPF evidence for GCE PSs is highly susceptible to certain systematic errors, and should not at present be taken to robustly disfavor a dominantly smooth GCE.

We calculate the one-loop corrections to TeV scale dark matter annihilation in a model where the dark matter is described by an SU(2)$_L$ triplet of Majorana fermions, such as the wino. We use this framework to determine the high and low-scale MS-bar matching coefficients at both the dark matter and weak boson mass scales at one loop. Part of this calculation has previously been performed in the literature numerically; we find our analytic result differs from the earlier work and discuss potential origins of this disagreement. Our result is used to extend the dark matter annihilation rate to NLL' (NLL+O($\alpha_2$) corrections) which enables a precise determination of indirect detection signatures in present and upcoming experiments.

A TeV scale electroweak particle is a well motivated candidate for the dark matter (DM) of our Universe. Yet such a particle may only be detectable using indirect detection instruments sensitive to TeV-scale gamma rays that can result from dark matter annihilations. We present a mock analysis of the sensitivity for the present ground-based Cherenkov telescope array H.E.S.S. (High Energy Spectroscopic System) to detect TeV scale DM in the Galactic Center region. The work combines next-to-leading-logarithmic order calculations for the annihilation photon spectrum, as well as a comprehensive treatment of detector effects and expected backgrounds. Forecast limits on the sensitivity of H.E.S.S. have been derived across the important TeV mass range, assuming different DM density profiles and focusing on the canonical WIMP dark matter candidate Wino.These limits test our present and future ability to probe the predicted thermal cross section for some of the most promising DM candidates that could be discovered in the coming decade.

Development of the hardware, data analysis, and simulation techniques for large compact radio arrays dedicated to mapping the 21 cm line of neutral hydrogen gas has proven to be more difficult than imagined twenty years ago when such telescopes were first proposed. Despite tremendous technical and methodological advances, there are several outstanding questions on how to optimally calibrate and analyze such data. On the positive side, it has become clear that the outstanding issues are purely technical in nature and can be solved with sufficient development activity. Such activity will enable science across redshifts, from early galaxy evolution in the pre-reionization era to dark energy evolution at low redshift.

We consider dark matter models in which the mass splitting between the dark matter particles and their annihilation products is tiny. Compared to the previously proposed Forbidden Dark Matter scenario, the mass splittings we consider are much smaller, and are allowed to be either positive or negative. To emphasize this modification, we dub our scenario "Impeded Dark Matter". We demonstrate that Impeded Dark Matter can be easily realized without requiring tuning of model parameters. For negative mass splitting, we demonstrate that the annihilation cross-section for Impeded Dark Matter depends linearly on the dark matter velocity or may even be kinematically forbidden, making this scenario almost insensitive to constraints from the cosmic microwave background and from observations of dwarf galaxies. Accordingly, it may be possible for Impeded Dark Matter to yield observable signals in clusters or the Galactic center, with no corresponding signal in dwarfs. For positive mass splitting, we show that the annihilation cross-section is suppressed by the small mass splitting, which helps light dark matter to survive increasingly stringent constraints from indirect searches. As specific realizations for Impeded Dark Matter, we introduce a model of vector dark matter from a hidden $SU(2)$ sector, and a composite dark matter scenario based on a QCD-like dark sector.

We construct an effective field theory (EFT) description of the hard photon spectrum for heavy WIMP annihilation. This facilitates precision predictions relevant for line searches, and allows the incorporation of non-trivial energy resolution effects. Our framework combines techniques from non-relativistic EFTs and soft-collinear effective theory (SCET), as well as its multi-scale extensions that have been recently introduced for studying jet substructure. We find a number of interesting features, including the simultaneous presence of SCET$_{\text{I}}$ and SCET$_{\text{II}}$ modes, as well as collinear-soft modes at the electroweak scale. We derive a factorization formula that enables both the resummation of the leading large Sudakov double logarithms that appear in the perturbative spectrum, and the inclusion of Sommerfeld enhancement effects. Consistency of this factorization is demonstrated to leading logarithmic order through explicit calculation. Our final result contains both the exclusive and the inclusive limits, thereby providing a unifying description of these two previously-considered approximations. We estimate the impact on experimental sensitivity, focusing for concreteness on an SU(2)$_{W}$ triplet fermion dark matter - the pure wino - where the strongest constraints are due to a search for gamma-ray lines from the Galactic Center. We find numerically significant corrections compared to previous results, thereby highlighting the importance of accounting for the photon spectrum when interpreting data from current and future indirect detection experiments.

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 provide precise predictions for the hard photon spectrum resulting from neutral SU$(2)_W$ triplet (wino) dark matter annihilation. Our calculation is performed utilizing an effective field theory expansion around the endpoint region where the photon energy is near the wino mass. This has direct relevance to line searches at indirect detection experiments. We compute the spectrum at next-to-leading logarithmic (NLL) accuracy within the framework established by a factorization formula derived previously by our collaboration. This allows simultaneous resummation of large Sudakov logarithms (arising from a restricted final state) and Sommerfeld effects. Resummation at NLL accuracy shows good convergence of the perturbative series due to the smallness of the electroweak coupling constant - scale variation yields uncertainties on our NLL prediction at the level of $5\%$. We highlight a number of interesting field theory effects that appear at NLL associated with the presence of electroweak symmetry breaking, which should have more general applicability. We also study the importance of using the full spectrum as compared with a single endpoint bin approximation when computing experimental limits. Our calculation provides a state of the art prediction for the hard photon spectrum that can be easily generalized to other DM candidates, allowing for the robust interpretation of data collected by current and future indirect detection experiments.

We present results from N-body simulations of self-interacting dark matter (SIDM) subhalos, which could host ultra-faint dwarf spheroidal galaxies, inside a Milky-Way-like main halo. We find that high-concentration subhalos are driven to gravothermal core collapse, while low-concentration subhalos develop large (kpc-sized) low-density cores, with both effects depending sensitively on the satellite's orbit and the self-interaction cross section over mass $\sigma/m$. The overall effect for $\sigma/m \gtrsim 3 \ \rm cm^2/g$ is to increase the range of inner densities, potentially explaining the observed diversity of Milky Way satellites, which include compact systems like Draco and Segue 1 that are dense in dark matter, and less dense, diffuse systems like Sextans and Crater II. We discuss possible ways of distinguishing SIDM models from collisionless dark matter models using the inferred dark matter densities and stellar sizes of the dwarf spheroidal galaxies.

We present a new public Python package, DarkHistory, for computing the effects of dark matter annihilation and decay on the temperature and ionization history of the early universe. DarkHistory simultaneously solves for the evolution of the free electron fraction and gas temperature, and for the cooling of annihilation/decay products and the secondary particles produced in the process. Consequently, we can self-consistently include the effects of both astrophysical and exotic sources of heating and ionization, and automatically take into account backreaction, where modifications to the ionization/temperature history in turn modify the energy-loss processes for injected particles. We present a number of worked examples, demonstrating how to use the code in a range of different configurations, in particular for arbitrary dark matter masses and annihilation/decay final states. Possible applications of DarkHistory include mapping out the effects of dark matter annihilation/decay on the global 21cm signal and the epoch of reionization, as well as the effects of exotic energy injections other than dark matter annihilation/decay. The code is available at https://github.com/hongwanliu/DarkHistory with documentation at https://darkhistory.readthedocs.io . Data files required to run the code can be downloaded at https://doi.org/10.7910/DVN/DUOUWA .

Statistical evidence has previously suggested that the Galactic Center GeV Excess (GCE) originates largely from point sources, and not from annihilating dark matter. We examine the impact of unmodeled source populations on identifying the true origin of the GCE using non-Poissonian template fitting (NPTF) methods. In a proof-of-principle example with simulated data, we discover that unmodeled sources in the Fermi Bubbles can lead to a dark matter signal being misattributed to point sources by the NPTF. We discover striking behavior consistent with a mismodeling effect in the real Fermi data, finding that large artificial injected dark matter signals are completely misattributed to point sources. Consequently, we conclude that dark matter may provide a dominant contribution to the GCE after all.

Self-interacting dark matter provides a promising alternative for the cold dark matter paradigm to solve potential small-scale galaxy formation problems. Nearly all self-interacting dark matter simulations so far have considered only elastic collisions. Here we present simulations of a galactic halo within a generic inelastic model using a novel numerical implementation in the Arepo code to study arbitrary multi-state inelastic dark matter scenarios. For this model we find that inelastic self-interactions can: (i) create larger subhalo density cores compared to elastic models for the same cross section normalisation; (ii) lower the abundance of satellites without the need for a power spectrum cutoff; (iii) reduce the total halo mass by about 10%; (iv) inject the energy equivalent of O(100) million Type II supernovae in galactic haloes through level de-excitation; (v) avoid the gravothermal catastrophe due to removal of particles from halo centers. We conclude that a ~5 times larger elastic cross section is required to achieve the same central density reduction as the inelastic model. This implies that well-established constraints on self-interacting cross sections have to be revised if inelastic collisions are the dominant mode. In this case significantly smaller cross sections can achieve the same core density reduction thereby increasing the parameter space of allowed models considerably.

Cosmological observations offer unique and robust avenues for probing the fundamental nature of dark matter particles-they broadly test a range of compelling theoretical scenarios, often surpassing or complementing the reach of terrestrial and other experiments. We discuss observational and theoretical advancements that will play a pivotal role in realizing a strong program of cosmological searches for the identity of dark matter in the coming decade. Specifically, we focus on measurements of the cosmic-microwave-background anisotropy and spectral distortions, and tracers of structure (such as the Lyman-$\alpha$ forest, galaxies, and the cosmological 21-cm signal).

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 examine the possible indirect signatures of dark matter annihilation processes with a non-standard scaling with the dark matter density, and in particular the case where more than two dark matter particles participate in the annihilation process. We point out that such processes can be strongly enhanced at low velocities without violating unitarity, similar to Sommerfeld enhancement in the standard case of two-body annihilation, potentially leading to visible signals in indirect searches. We study in detail the impact of such multi-body annihilations on the ionization history of the universe and consequently the cosmic microwave background, and find that unlike in the two-body case, the dominant signal can naturally arise from the end of the cosmic dark ages, after the onset of structure formation. We examine the complementary constraints from the Galactic Center, Galactic halo, and galaxy clusters, and outline the circumstances under which each search would give rise to the strongest constraints. We also show that if there is a population of ultra-compact dense dark matter clumps present in the Milky Way with sufficiently steep density profile, then it might be possible to detect point sources illuminated by multi-body annihilation, even if there is no large low-velocity enhancement. Finally, we provide a case study of a model where 3-body annihilation dominates the freezeout process, and in particular the resonant regime where a large low-velocity enhancement is naturally generated.

Observing gamma rays using ground-based atmospheric Cherenkov telescopes provides one of the only probes of heavy weakly interacting dark matter. A canonical target is the thermal wino, for which the strongest limits come from searches for photon lines from annihilations in the Galactic Center. Irreducible finite energy resolution effects motivate refining the prediction for a wino signal beyond the photon line approximation; recently, modern effective field theory techniques have been utilized to obtain a precise calculation of the full photon energy spectrum from wino annihilation. In this paper, we investigate the implications for a realistic mock H.E.S.S.-like line search. We emphasize the impact of including the non-trivial spectral shape, and we carefully treat the region of interest, presenting results for choices between $1^{\circ}$ and $4^{\circ}$ from the Galactic Center. Projected limits for wino masses from $1$-$70$ TeV are interpreted as a constraint on the wino annihilation rate, or alternatively as the minimum core size required such that the wino is not excluded. If there is a thermal wino, H.E.S.S. will be able to probe cores of several kpc, which would begin to cause tension between this dark matter candidate and astrophysical observations/simulations.

The search for the particle nature of dark matter has given rise to a number of experimental, theoretical and statistical challenges. Here, we report on a number of these statistical challenges and new techniques to address them, as discussed in the DMStat workshop held Feb 26 - Mar 3 2018 at the Banff International Research Station for Mathematical Innovation and Discovery (BIRS) in Banff, Alberta.

Weakly Interacting Massive Particles (WIMPs) have long reigned as one of the leading classes of dark matter candidates. The observed dark matter abundance can be naturally obtained by freezeout of weak-scale dark matter annihilations in the early universe. This "thermal WIMP" scenario makes direct predictions for the total annihilation cross section that can be tested in present-day experiments. While the dark matter mass constraint can be as high as $m_\chi\gtrsim100$ GeV for particular annihilation channels, the constraint on the total cross section has not been determined. We construct the first model-independent limit on the WIMP total annihilation cross section, showing that allowed combinations of the annihilation-channel branching ratios considerably weaken the sensitivity. For thermal WIMPs with s-wave $2\rightarrow2$ annihilation to visible final states, we find the dark matter mass is only known to be $m_\chi\gtrsim20$ GeV. This is the strongest largely model-independent lower limit on the mass of thermal-relic WIMPs, together with the upper limit on the mass from the unitarity bound ($m_\chi\lesssim 100$ TeV), it defines what we call the "WIMP window". To probe the remaining mass range, we outline ways forward.

Measurements of the temperature of the baryons at the end of the cosmic dark ages can potentially set very precise constraints on energy injection from exotic sources, such as annihilation or decay of the dark matter. However, additional effects that lower the gas temperature can substantially weaken the expected constraints on exotic energy injection, whereas additional radiation backgrounds can conceal the effect of an increased gas temperature in measurements of the 21-cm hyperfine transition of neutral hydrogen. Motivated in part by recent claims of a detection of 21-cm absorption from a redshift of 17 by the EDGES experiment, we derive the constraints on dark matter annihilation and decay that can be placed in the presence of extra radiation backgrounds or effects that modify the gas temperature, such as dark matter-baryon scattering and early baryon-photon decoupling. We find that if the EDGES observation is confirmed, then constraints on light dark matter decaying or annihilating to electrons will in most scenarios be stronger than existing state-of-the-art limits from the cosmic microwave background, potentially by several orders of magnitude. More generally, our results allow mapping any future measurement of the global 21-cm signal into constraints on dark matter annihilation and decay, within the broad range of scenarios we consider.

We present and compare several cosmological constraints on the cross section for elastic scattering between dark matter (DM) and baryons, for cross sections with a range of power-law dependences on the DM-baryon relative velocity $v$, especially focusing on the case of $\sigma \propto v^{-4}$. We study constraints spanning a wide range of epochs in cosmological history, from pre-recombination distortions to the blackbody spectrum and anisotropies of the cosmic microwave background (CMB), to modifications to the intergalactic medium temperature and the resulting 21cm signal, and discuss the allowed signals in the latter channels given the constraints from the former. We improve previous constraints on DM-baryon scattering from the CMB anisotropies, demonstrate via principal component analysis that the effect on the CMB can be written as a simple function of DM mass, and map out the redshifts dominating this signal. We show that given high-redshift constraints on DM-baryon scattering, a $v^{-4}$ scaling of the cross section for light DM would be sufficient to explain the deep 21cm absorption trough recently claimed by the EDGES experiment, if 100% of the DM scatters with baryons. For millicharged DM models proposed to explain the observation, where only a small fraction of the DM interacts, we estimate that a PIXIE-like future experiment measuring CMB spectral distortion could test the relevant parameter space.

We explore the possibility that bound states involving dark matter particles could be detected by resonance searches at the LHC, and the generic implications of such scenarios for indirect and direct detection. We demonstrate that resonance searches are complementary to mono-jet searches and can probe dark matter masses above 1 TeV with current LHC data. We argue that this parameter regime, where the bound-state resonance channel is the most sensitive probe of the dark sector, arises most naturally in the context of non-trivial dark sectors with large couplings, nearly-degenerate dark-matter-like states, and multiple force carriers. The presence of bound states detectable by the LHC implies a minimal Sommerfeld enhancement that is appreciable, and potentially also radiative bound state formation in the Galactic halo, leading to large signals in indirect searches. We calculate these complementary constraints, which favor either models where the bound-state-forming dark matter constitutes a small fraction of the total density, or models where the late-time annihilation is suppressed at low velocities or late times. We present concrete examples of models that satisfy all these constraints and where the LHC resonance search is the most sensitive probe of the dark sector.

The $\textit{Fermi}$-LAT Collaboration recently presented a new catalog of gamma-ray sources located within the $40^{\circ} \times 40^{\circ}$ region around the Galactic Center~(Ajello et al. 2017) -- the Second Fermi Inner Galaxy (2FIG) catalog. Utilizing this catalog, they analyzed models for the spatial distribution and luminosity function of sources with a pulsar-like gamma-ray spectrum. Ajello et al. 2017 v1 also claimed to detect, in addition to a disk-like population of pulsar-like sources, an approximately 7$\sigma$ preference for an additional centrally concentrated population of pulsar-like sources, which they referred to as a "Galactic Bulge" population. Such a population would be of great interest, as it would support a pulsar interpretation of the gamma-ray excess that has long been observed in this region. In an effort to further explore the implications of this new source catalog, we attempted to reproduce the results presented by the $\textit{Fermi}$-LAT Collaboration, but failed to do so. Mimicking as closely as possible the analysis techniques undertaken in Ajello et al. 2017, we instead find that our likelihood analysis favors a very different spatial distribution and luminosity function for these sources. Most notably, our results do not exhibit a strong preference for a "Galactic Bulge" population of pulsars. Furthermore, we find that masking the regions immediately surrounding each of the 2FIG pulsar candidates does $\textit{not}$ significantly impact the spectrum or intensity of the Galactic Center gamma-ray excess. Although these results refute the claim of strong evidence for a centrally concentrated pulsar population presented in Ajello et al. 2017, they neither rule out nor provide support for the possibility that the Galactic Center excess is generated by a population of low-luminosity and currently largely unobserved pulsars.

These lectures, presented at TASI 2016: Anticipating the Next Discoveries in Particle Physics, provide an introduction to some key methods and tools of indirect dark matter searches. Topics covered include estimation of dark matter signals, thermal freezeout and related scenarios, potential effects of dark matter annihilation on the early universe, modeling photon signals from annihilation or decay, and a brief and qualitative introduction to diffusive propagation of cosmic rays. The second half of the notes gives a status report (circa summer 2016) on selected experimental searches, the resulting constraints and some potential signal candidates. These notes are intended as an introduction to indirect dark matter searches for graduate students, focusing on back-of-the-envelope estimates and useful concepts rather than detailed quantitative computations.

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.

Observations by the Fermi-LAT have uncovered a bright, spherically symmetric excess surrounding the center of the Milky Way galaxy. The spectrum of the gamma-ray excess peaks sharply at an energy ~2 GeV, exhibiting a hard spectrum at lower energies, and falls off quickly above an energy ~5 GeV. The spectrum of the excess above ~10 GeV is potentially an important discriminator between different physical models for its origin. We focus our study on observations of the gamma-ray excess at energies exceeding 10 GeV, finding: (1) a statistically significant excess remains in the energy range 9.5-47.5 GeV, which is not degenerate with known diffuse emission templates such as the Fermi Bubbles, (2) the radial profile of the excess at high energies remains relatively consistent with data near the spectral peak, (3) the data above ~5 GeV prefer a slightly greater ellipticity with a major axis oriented perpendicular to the Galactic plane. Using the recently developed non-Poissonian template fit, we find mild evidence for a point-source origin for the high-energy excess, although given the statistical and systematic uncertainties we show that a smooth origin of the high-energy emission cannot be ruled out. We discuss the implication of these findings for pulsar and dark matter models of the gamma-ray excess. Finally we provide a number of updated measurements of the gamma-ray excess, utilizing novel diffuse templates and the Pass 8 dataset.

The spectrum of Weakly-Interacting-Massive-Particle (WIMP) dark matter generically possesses bound states when the WIMP mass becomes sufficiently large relative to the mass of the electroweak gauge bosons. The presence of these bound states enhances the annihilation rate via resonances in the Sommerfeld enhancement, but they can also be produced directly with the emission of a low-energy photon. In this work we compute the rate for SU(2) triplet dark matter (the wino) to bind into WIMPonium -- which is possible via single-photon emission for wino masses above 5 TeV for relative velocity v < O(10^{-2}) -- and study the subsequent decays of these bound states. We present results with applications beyond the wino case, e.g. for dark matter inhabiting a nonabelian dark sector; these include analytic capture and transition rates for general dark sectors in the limit of vanishing force carrier mass, efficient numerical routines for calculating positive and negative-energy eigenstates of a Hamiltonian containing interactions with both massive and massless force carriers, and a study of the scaling of bound state formation in the short-range Hulthen potential. In the specific case of the wino, we find that the rate for bound state formation is suppressed relative to direct annihilation, and so provides only a small correction to the overall annihilation rate. The soft photons radiated by the capture process and by bound state transitions could permit measurement of the dark matter's quantum numbers; for wino-like dark matter, such photons are rare, but might be observable by a future ground-based gamma-ray telescope combining large effective area and a low energy threshold.

Dark matter (DM) halos have long been known to be triaxial, but in studies of possible annihilation and decay signals they are often treated as approximately spherical. In this work, we examine the asymmetry of potential indirect detection signals of DM annihilation and decay, exploiting the large statistics of the hydrodynamic simulation Illustris. We carefully investigate the effects of the baryons on the sphericity of annihilation and decay signals for both the case where the observer is at 8.5 kpc from the center of the halo (exemplified in the case of Milky Way-like halos), and for an observer situated well outside the halo. In the case of Galactic signals, we find that both annihilation and decay signals are expected to be quite symmetric, with axis ratios very different from 1 occurring rarely. In the case of extragalactic signals, while decay signals are still preferentially spherical, the axis ratio for annihilation signals has a much flatter distribution, with elongated profiles appearing frequently. Many of these elongated profiles are due to large subhalos and/or recent mergers. Comparing to gamma-ray emission from the Milky Way and X-ray maps of clusters, we find that the gamma-ray background appears less spherical/more elongated than the expected DM signal from the large majority of halos, and the Galactic gamma ray excess appears very spherical, while the X- ray data would be difficult to distinguish from a DM signal by elongation/sphericity measurements alone.

Precise measurements of the temperature and polarization anisotropies of the cosmic microwave background can be used to constrain the annihilation and decay of dark matter. In this work, we demonstrate via principal component analysis that the imprint of dark matter decay on the cosmic microwave background can be approximately parameterized by a single number for any given dark matter model. We develop a simple prescription for computing this model-dependent detectability factor, and demonstrate how this approach can be used to set model-independent bounds on a large class of decaying dark matter scenarios. We repeat our analysis for decay lifetimes shorter than the age of the universe, allowing us to set constraints on metastable species other than the dark matter decaying at early times, and decays that only liberate a tiny fraction of the dark matter mass energy. We set precise bounds and validate our principal component analysis using a Markov Chain Monte Carlo approach and Planck 2015 data.

Dark matter annihilation or decay could have a significant impact on the ionisation and thermal history of the universe. In this paper, we study the potential contribution of dark matter annihilation (s-wave- or p-wave-dominated) or decay to cosmic reionisation, via the production of electrons, positrons and photons. We map out the possible perturbations to the ionisation and thermal histories of the universe due to dark matter processes, over a broad range of velocity-averaged annihilation cross-sections/decay lifetimes and dark matter masses. We have employed recent numerical studies of the efficiency with which annihilation/decay products induce heating and ionization in the intergalactic medium, and in this work extended them down to a redshift of $1+z = 4$ for two different reionisation scenarios. We also improve on earlier studies by using the results of detailed structure formation models of dark matter haloes and subhaloes that are consistent with up-to-date $N$-body simulations, with estimates on the uncertainties that originate from the smallest scales. We find that for dark matter models that are consistent with experimental constraints, a contribution of more than 10% to the ionisation fraction at reionisation is disallowed for all annihilation scenarios. Such a contribution is possible only for decays into electron/positron pairs, for light dark matter with mass $m_\chi \lesssim $ 100 MeV, and a decay lifetime $\tau_\chi \sim 10^{24} - 10^{25}$ s.

If dark matter inhabits an expanded "hidden sector", annihilations may proceed through sequential decays or multi-body final states. We map out the potential signals and current constraints on such a framework in indirect searches, using a model-independent setup based on multi-step hierarchical cascade decays. While remaining agnostic to the details of the hidden sector model, our framework captures the generic broadening of the spectrum of secondary particles (photons, neutrinos, e+e- and antiprotons) relative to the case of direct annihilation to Standard Model particles. We explore how indirect constraints on dark matter annihilation limit the parameter space for such cascade/multi-particle decays. We investigate limits from the cosmic microwave background by Planck, the Fermi measurement of photons from the dwarf galaxies, and positron data from AMS-02. The presence of a hidden sector can change the constraints on the dark matter annihilation cross section by up to an order of magnitude in either direction (although the effect can be much smaller). We find that generally the bound from the Fermi dwarfs is most constraining for annihilations to photon-rich final states, while AMS-02 is most constraining for electron and muon final states; however in certain instances the CMB bounds overtake both, due to their approximate independence of the details of the hidden sector cascade. We provide the full set of cascade spectra considered here as publicly available code with examples at http://web.mit.edu/lns/research/CascadeSpectra.html.

The Fermi Bubbles are giant Galactic structures observed in both gamma-rays and microwaves. Recent studies have found support for the hypothesis that the gamma-ray and microwave emission can both be understood as arising from a hard cosmic-ray electron population within the volume of the Bubbles, via inverse Compton scattering and synchrotron radiation respectively. The relative rates of these processes are set by the relative energy density of the interstellar radiation field and the magnetic field within the Bubbles; consequently, under the hypothesis of a common origin, the combination of the gamma-ray and microwave measurements can be used to estimate the magnetic field within the Bubbles. We revisit the consistency of this hypothesis on a latitude-by-latitude basis, using data from Fermi, WMAP and Planck; estimate the variation of the electron spectrum within the Bubbles; and infer bounds on the magnetic field within the Bubbles as a function of distance from the Galactic plane. We find that while the microwave and gamma-ray spectra are generally consistent with the leptonic hypothesis for few-microGauss magnetic fields, there appears to be a preference for spectral hardening in the microwaves at mid-latitudes (especially in the |b|~ 25-35 degree range) that is not mirrored in the gamma rays. This result may hint at a non-leptonic contribution to the gamma-ray spectra; however, the discrepancy can be reconciled in purely leptonic models if the cutoff energy for the electrons is lower in this latitude range and the spectrum below the cutoff is harder.

We present a new method to characterize unresolved point sources (PSs), generalizing traditional template fits to account for non-Poissonian photon statistics. We apply this method to Fermi Large Area Telescope gamma-ray data to characterize PS populations at high latitudes and in the Inner Galaxy. We find that PSs (resolved and unresolved) account for ~50% of the total extragalactic gamma-ray background in the energy range ~1.9 to 11.9 GeV. Within 10$^\circ$ of the Galactic Center with $|b| \geq 2^\circ$, we find that ~5-10% of the flux can be accounted for by a population of unresolved PSs, distributed consistently with the observed ~GeV gamma-ray excess in this region. The excess is fully absorbed by such a population, in preference to dark-matter annihilation. The inferred source population is dominated by near-threshold sources, which may be detectable in future searches.

Any injection of electromagnetically interacting particles during the cosmic dark ages will lead to increased ionization, heating, production of Lyman-alpha photons and distortions to the energy spectrum of the cosmic microwave background, with potentially observable consequences. In this note we describe numerical results for the low-energy electrons and photons produced by the cooling of particles injected at energies from keV to multi-TeV scales, at arbitrary injection redshifts (but focusing on the post-recombination epoch). We use these data, combined with existing calculations modeling the cooling of these low-energy particles, to estimate the resulting contributions to ionization, excitation and heating of the gas, and production of low-energy photons below the threshold for excitation and ionization. We compute corrected deposition-efficiency curves for annihilating dark matter, and demonstrate how to compute equivalent curves for arbitrary energy-injection histories. These calculations provide the necessary inputs for the limits on dark matter annihilation presented in the accompanying Paper I, but also have potential applications in the context of dark matter decay or de-excitation, decay of other metastable species, or similar energy injections from new physics. We make our full results publicly available at http://nebel.rc.fas.harvard.edu/epsilon, to facilitate further independent studies. In particular, we provide the full low-energy electron and photon spectra, to allow matching onto more detailed codes that describe the cooling of such particles at low energies.

Recent measurements of the cosmic microwave background (CMB) anisotropies by Planck provide a sensitive probe of dark matter annihilation during the cosmic dark ages, and specifically constrain the annihilation parameter $f_\mathrm{eff} \langle \sigma v \rangle/m_\chi$. Using new results (Paper II) for the ionization produced by particles injected at arbitrary energies, we calculate and provide $f_\mathrm{eff}$ values for photons and $e^+e^-$ pairs injected at keV-TeV energies; the $f_\mathrm{eff}$ value for any dark matter model can be obtained straightforwardly by weighting these results by the spectrum of annihilation products. This result allows the sensitive and robust constraints on dark matter annihilation presented by the Planck Collaboration to be applied to arbitrary dark matter models with $s$-wave annihilation. We demonstrate the validity of this approach using principal component analysis. As an example, we integrate over the spectrum of annihilation products for a range of Standard Model final states to determine the CMB bounds on these models as a function of dark matter mass, and demonstrate that the new limits generically exclude models proposed to explain the observed high-energy rise in the cosmic ray positron fraction. We make our results publicly available at http://nebel.rc.fas.harvard.edu/epsilon.

If dark matter is embedded in a non-trivial dark sector, it may annihilate and decay to lighter dark-sector states which subsequently decay to the Standard Model. Such scenarios - with annihilation followed by cascading dark-sector decays - can explain the apparent excess GeV gamma-rays identified in the central Milky Way, while evading bounds from dark matter direct detection experiments. Each 'step' in the cascade will modify the observable signatures of dark matter annihilation and decay, shifting the resulting photons and other final state particles to lower energies and broadening their spectra. We explore, in a model-independent way, the effect of multi-step dark-sector cascades on the preferred regions of parameter space to explain the GeV excess. We find that the broadening effects of multi-step cascades can admit final states dominated by particles that would usually produce too sharply peaked photon spectra; in general, if the cascades are hierarchical (each particle decays to substantially lighter particles), the preferred mass range for the dark matter is in all cases 20-150 GeV. Decay chains that have nearly-degenerate steps, where the products are close to half the mass of the progenitor, can admit much higher DM masses. We map out the region of mass/cross-section parameter space where cascades (degenerate, hierarchical or a combination) can fit the signal, for a range of final states. In the current work, we study multi-step cascades in the context of explaining the GeV excess, but many aspects of our results are general and can be extended to other applications.

Past studies have identified a spatially extended excess of $\sim$1-3 GeV gamma rays from the region surrounding the Galactic Center, consistent with the emission expected from annihilating dark matter. We revisit and scrutinize this signal with the intention of further constraining its characteristics and origin. By applying cuts to the \textit{Fermi} event parameter CTBCORE, we suppress the tails of the point spread function and generate high resolution gamma-ray maps, enabling us to more easily separate the various gamma-ray components. Within these maps, we find the GeV excess to be robust and highly statistically significant, with a spectrum, angular distribution, and overall normalization that is in good agreement with that predicted by simple annihilating dark matter models. For example, the signal is very well fit by a 36-51 GeV dark matter particle annihilating to $b\bar{b}$ with an annihilation cross section of $\sigma v = (1-3)\times 10^{-26}$ cm$^3$/s (normalized to a local dark matter density of 0.4 GeV/cm$^3$). Furthermore, we confirm that the angular distribution of the excess is approximately spherically symmetric and centered around the dynamical center of the Milky Way (within $\sim$$0.05^{\circ}$ of Sgr A$^*$), showing no sign of elongation along the Galactic Plane. The signal is observed to extend to at least $\simeq10^{\circ}$ from the Galactic Center, disfavoring the possibility that this emission originates from millisecond pulsars.

Updated constraints on dark matter cross section and mass are presented combining CMB power spectrum measurements from Planck, WMAP9, ACT, and SPT as well as several low-redshift datasets (BAO, HST, supernovae). For the CMB datasets, we combine WMAP9 temperature and polarization data for l <= 431 with Planck temperature data for 432 < l < 2500, ACT and SPT data for l > 2500, and Planck CMB four-point lensing measurements. We allow for redshift-dependent energy deposition from dark matter annihilation by using a `universal' energy absorption curve. We also include an updated treatment of the excitation, heating, and ionization energy fractions, and provide updated deposition efficiency factors (f_eff) for 41 different dark matter models. Assuming perfect energy deposition (f_eff = 1) and a thermal cross section, dark matter masses below 26 GeV are excluded at the 2-sigma level. Assuming a more generic efficiency of f_eff = 0.2, thermal dark matter masses below 5 GeV are disfavored at the 2-sigma level. These limits are a factor of ~2 improvement over those from WMAP9 data alone. These current constraints probe, but do not exclude, dark matter as an explanation for reported anomalous indirect detection observations from AMS-02/PAMELA and the Fermi Gamma-ray Inner Galaxy data. They also probe relevant models that would explain anomalous direct detection events from CDMS, CRESST, CoGeNT, and DAMA, as originating from a generic thermal WIMP. Projected constraints from the full Planck release should improve the current limits by another factor of ~2, but will not definitely probe these signals. The proposed CMB Stage IV experiment will more decisively explore the relevant regions and improve upon the Planck constraints by another factor of ~2.

Self-interacting dark matter scenarios have recently attracted much attention, as a possible means to alleviate the tension between N-body simulations and observations of the dark matter distribution on galactic and sub-galactic scales. The presence of internal structure for the dark matter --- for example, a nearly-degenerate state in the spectrum that could decay, or be collisionally excited or de-excited --- has also been proposed as a possible means to address these discrepancies. Such internal structure can be a source of interesting signatures in direct and indirect dark matter searches, for example providing a novel explanation for the 3.5 keV line recently observed in galaxies and galaxy clusters. We analyze a simple model of dark matter self-scattering including a nearly-degenerate excited state, and develop an accurate analytic approximation for the elastic and inelastic s-wave cross sections, which is valid outside the perturbative regime provided the particle velocity is sufficiently low (this condition is also required for the s-wave to dominate over higher partial waves). We anticipate our results will be useful in incorporating inelastic self-scattering into N-body simulations, in order to study the quantitative impact of nearly-degenerate states in the dark matter spectrum on galactic structure and dynamics, and in computing the indirect signatures of multi-state dark matter.

We formulate an effective field theory description for SU(2)$_L$ triplet fermionic dark matter by combining nonrelativistic dark matter with gauge bosons in the soft-collinear effective theory. For a given dark matter mass, the annihilation cross section to line photons is obtained with 5% precision by simultaneously including Sommerfeld enhancement and the resummation of electroweak Sudakov logarithms at next-to-next-to-leading logarithmic order. Using these results, we present more accurate and precise predictions for the gamma-ray line signal from annihilation, updating both existing constraints and the reach of future experiments.

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

Dark sectors, consisting of new, light, weakly-coupled particles that do not interact with the known strong, weak, or electromagnetic forces, are a particularly compelling possibility for new physics. Nature may contain numerous dark sectors, each with their own beautiful structure, distinct particles, and forces. This review summarizes the physics motivation for dark sectors and the exciting opportunities for experimental exploration. It is the summary of the Intensity Frontier subgroup "New, Light, Weakly-coupled Particles" of the Community Summer Study 2013 (Snowmass). We discuss axions, which solve the strong CP problem and are an excellent dark matter candidate, and their generalization to axion-like particles. We also review dark photons and other dark-sector particles, including sub-GeV dark matter, which are theoretically natural, provide for dark matter candidates or new dark matter interactions, and could resolve outstanding puzzles in particle and astro-particle physics. In many cases, the exploration of dark sectors can proceed with existing facilities and comparatively modest experiments. A rich, diverse, and low-cost experimental program has been identified that has the potential for one or more game-changing discoveries. These physics opportunities should be vigorously pursued in the US and elsewhere.

Anisotropies of the cosmic microwave background (CMB) have proven to be a very powerful tool to constrain dark matter annihilation at the epoch of recombination. However, CMB constraints are currently derived using a number of reasonable but yet un-tested assumptions that could potentially lead to a misestimation of the true bounds. In this paper we examine the potential impact of these systematic effects. In particular, we separately study the propagation of the secondary particles produced by annihilation in two energy regimes; first following the shower from the initial particle energy to the keV scale, and then tracking the resulting secondary particles from this scale to the absorption of their energy as heat, ionization, or excitation of the medium. We improve both the high and low energy parts of the calculation, in particular finding that our more accurate treatment of losses to sub-10.2 eV photons produced by scattering of high-energy electrons weakens the constraints on particular DM annihilation models by up to a factor of two. On the other hand, we find that the uncertainties we examine for the low energy propagation do not significantly affect the results for current and upcoming CMB data. We include the evaluation of the precise amount of excitation energy, in the form of Lyman-alpha photons, produced by the propagation of the shower, and examine the effects of varying the Helium fraction and Helium ionization fraction. In the recent literature, simple approximations for the fraction of energy absorbed in different channels have often been used to derive CMB constraints: we assess the impact of using accurate versus approximate energy fractions. Finally we check that the choice of recombination code (between RECFAST v1.5 and COSMOREC), to calculate the evolution of the free electron fraction in the presence of dark matter annihilation, introduces negligible differences.

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.

We study the variation of the spectrum of the Fermi Bubbles with Galactic latitude. Far from the Galactic plane (|b| > 30 degrees), the observed gamma-ray emission is nearly invariant with latitude, and is consistent with arising from inverse Compton scattering of the interstellar radiation field by cosmic-ray electrons with an approximately power-law spectrum. The same electrons in the presence of microgauss-scale magnetic fields can also generate the the observed microwave "haze". At lower latitudes (b < 20 degrees), in contrast, the spectrum of the emission correlated with the Bubbles possesses a pronounced spectral feature peaking at 1-4 GeV (in E^2 dN/dE) which cannot be generated by any realistic spectrum of electrons. Instead, we conclude that a second (non-inverse-Compton) emission mechanism must be responsible for the bulk of the low-energy, low-latitude emission. This second component is spectrally similar to the excess GeV emission previously reported from the Galactic Center (GC), and also appears spatially consistent with a luminosity per volume falling approximately as r^-2.4, where r is the distance from the GC. We argue that the spectral feature visible in the low-latitude Bubbles is the extended counterpart of the GC excess, now detected out to at least 2-3 kpc from the GC. The spectrum and angular distribution of the signal is consistent with that predicted from ~10 GeV dark matter particles annihilating to leptons, or from ~50 GeV dark matter particles annihilating to quarks, following a distribution similar to the canonical Navarro-Frenk-White (NFW) profile. We also consider millisecond pulsars as a possible astrophysical explanation for the signal, as observed millisecond pulsars possess a spectral cutoff at approximately the required energy. Any such scenario would require a large population of unresolved millisecond pulsars extending at least 2-3 kpc from the GC.

Dark matter annihilation or de-excitation, decay of metastable species, or other new physics may inject energetic electrons and photons into the photon-baryon fluid during and after recombination. As such particles cool, they partition their energy into a large number of efficiently ionizing electrons and photons, which in turn modify the ionization history. Recent work has provided a simple method for constraining arbitrary energy deposition histories using the cosmic microwave background (CMB); in this note, we present results describing the energy deposition histories for photons and electrons as a function of initial energy and injection redshift. With these results, the CMB bounds on any process injecting some arbitrary spectrum of electrons, positrons and/or photons with arbitrary redshift dependence can be immediately computed.

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

In models of dark matter (DM) with Sommerfeld-enhanced annihilation, where the annihilation rate scales as the inverse velocity, N-body simulations of DM structure formation suggest that the local annihilation signal may be dominated by small, dense, cold subhalos. This contrasts with the usual assumption of a signal originating from the smooth DM halo, with much higher velocity dispersion. Accounting for local substructure modifies the parameter space for which Sommerfeld-enhanced annihilating DM can explain the PAMELA and Fermi excesses. Limits from the inner galaxy and the cosmic microwave background are weakened, without introducing new tension with substructure-dependent limits, such as from dwarf galaxies or isotropic gamma-ray studies. With substructure, previously excluded parameter regions with mediators of mass ~1-200 MeV are now easily allowed. For O(MeV) mediators, subhalos in a specific range of host halo masses may be evaporated, further suppressing diffuse signals without affecting substructure in the Milky Way.

Observations by ARCADE-2 and other telescopes sensitive to low frequency radiation have revealed the presence of an isotropic radio background with a hard spectral index. The intensity of this observed background is found to exceed the flux predicted from astrophysical sources by a factor of approximately 5-6. In this article, we consider the possibility that annihilating dark matter particles provide the primary contribution to the observed isotropic radio background through the emission of synchrotron radiation from electron and positron annihilation products. For reasonable estimates of the magnetic fields present in clusters and galaxies, we find that dark matter could potentially account for the observed radio excess, but only if it annihilates mostly to electrons and/or muons, and only if it possesses a mass in the range of approximately 5-50 GeV. For such models, the annihilation cross section required to normalize the synchrotron signal to the observed excess is sigma v ~ (0.4-30) x 10^-26 cm^3/s, similar to the value predicted for a simple thermal relic (sigma v ~ 3 x 10^-26 cm^3/s). We find that in any scenario in which dark matter annihilations are responsible for the observed excess radio emission, a significant fraction of the isotropic gamma ray background observed by Fermi must result from dark matter as well.

High-precision measurements of the temperature and polarization anisotropies of the cosmic microwave background radiation have been previously employed to set robust constraints on dark matter annihilation during recombination. In this work we improve and generalize these constraints to apply to energy deposition during the recombination era with arbitrary redshift dependence. Our approach also provides more rigorous and model-independent bounds on dark matter annihilation and decay scenarios. We employ principal component analysis to identify a basis of weighting functions for the energy deposition. The coefficients of these weighting functions parameterize any energy deposition model and can be constrained directly by experiment. For generic energy deposition histories that are currently allowed by WMAP7 data, up to 3 principal component coefficients are measurable by Planck and up to 5 coefficients are measurable by an ideal cosmic variance limited experiment. For WIMP dark matter, our analysis demonstrates that the effect on the CMB is described well by a single (normalization) parameter and a "universal" redshift dependence for the energy deposition history. We give WMAP 7 constraints on both generic energy deposition histories and the universal WIMP case.

(Abridged) The extragalactic background light (EBL) observed at multiple wavelengths is a promising tool to probe the nature of dark matter since it might contain a significant contribution from gamma-rays produced promptly by dark matter annihilation. Additionally, the electrons and positrons produced in the annihilation give energy to the CMB photons to populate the EBL with X-rays and gamma-rays. We here create full-sky maps of the radiation from both of these contributions using the high-resolution Millennium-II simulation. We use upper limits on the contributions of unknown sources to the EBL to constrain the intrinsic properties of dark matter using a model-independent approach that can be employed as a template to test different particle physics models (including those with a Sommerfeld enhancement). These upper limits are based on observations spanning eight orders of magnitude in energy (from soft X-rays measured by CHANDRA to gamma-rays measured by Fermi), and on expectations for the contributions from blazars and star forming galaxies. To exemplify this approach, we analyze a set of benchmark Sommerfeld-enhanced models that give the correct dark matter abundance, satisfy CMB constraints, and fit the cosmic ray spectra measured by PAMELA and Fermi without any contribution from local subhalos. We find that these models are in conflict with the EBL constraints unless the contribution of unresolved subhalos is small and the annihilation signal dominates the EBL. We conclude that provided the collisionless cold dark matter paradigm is accurate, even for conservative estimates of the contribution from unresolved subhalos and astrophysical backgrounds, the EBL is at least as sensitive a probe of these types of scenarios as the CMB. Our results disfavor an explanation of the positron excess measured by PAMELA based only on dark matter annihilation in the smooth Galactic halo.

Anomalies in direct and indirect detection have motivated models of dark matter consisting of a multiplet of nearly-degenerate states, coupled by a new GeV-scale interaction. We perform a careful analysis of the thermal freezeout of dark matter annihilation in such a scenario. We compute the range of "boost factors" arising from Sommerfeld enhancement in the local halo for models which produce the correct relic density, and show the effect of including constraints on the saturated enhancement from the cosmic microwave background (CMB). We find that boost factors from Sommerfeld enhancement of up to ~800 are possible in the local halo. When the CMB bounds on the saturated enhancement are applied, the maximal boost factor is reduced to ~400 for 1-2 TeV dark matter and sub-GeV force carriers, but remains large enough to explain the observed Fermi and PAMELA electronic signals. We describe regions in the DM mass-boost factor plane where the cosmic ray data is well fit for a range of final states, and show that Sommerfeld enhancement alone is enough to provide the large annihilation cross sections required to fit the data, although for light mediator masses (less than ~200 MeV) there is tension with the CMB constraints in the absence of astrophysical boost factors from substructure. Additionally, we consider the circumstances under which WIMPonium formation is relevant and find for heavy WIMPs (greater than ~2 TeV) and soft-spectrum annihilation channels it can be an important consideration; we find regions with dark matter mass greater than 2.8 TeV that are consistent with the CMB bounds and have ~600-700 present-day boost factors.

Data from the Fermi-LAT reveal two large gamma-ray bubbles, extending 50 degrees above and below the Galactic center, with a width of about 40 degrees in longitude. The gamma-ray emission associated with these bubbles has a significantly harder spectrum (dN/dE ~ E^-2) than the IC emission from electrons in the Galactic disk, or the gamma-rays produced by decay of pions from proton-ISM collisions. There is no significant spatial variation in the spectrum or gamma-ray intensity within the bubbles, or between the north and south bubbles. The bubbles are spatially correlated with the hard-spectrum microwave excess known as the WMAP haze; the edges of the bubbles also line up with features in the ROSAT X-ray maps at 1.5-2 keV. We argue that these Galactic gamma-ray bubbles were most likely created by some large episode of energy injection in the Galactic center, such as past accretion events onto the central massive black hole, or a nuclear starburst in the last ~10 Myr. Dark matter annihilation/decay seems unlikely to generate all the features of the bubbles and the associated signals in WMAP and ROSAT; the bubbles must be understood in order to use measurements of the diffuse gamma-ray emission in the inner Galaxy as a probe of dark matter physics. Study of the origin and evolution of the bubbles also has the potential to improve our understanding of recent energetic events in the inner Galaxy and the high-latitude cosmic ray population.

We describe a simple test of the spatial uniformity of an ensemble of discrete events. Given an estimate for the point source luminosity function and an instrumental point spread function (PSF), a robust upper bound on the fractional point source contribution to a diffuse signal can be found. We verify with Monte Carlo tests that the statistic has advantages over the two-point correlation function for this purpose, and derive analytic estimates of the statistic's mean and variance as a function of the point source contribution. As a case study, we apply this statistic to recent gamma-ray data from the Fermi Large Area Telescope (LAT), and demonstrate that at energies above 10 GeV, the contribution of unresolved point sources to the diffuse emission is small in the region relevant for study of the WMAP Haze.

The Fermi Gamma-Ray Space Telescope reveals a diffuse inverse Compton signal in the inner Galaxy with a similar spatial morphology to the microwave haze observed by WMAP, supporting the synchrotron interpretation of the microwave signal. Using spatial templates, we regress out pi0 gammas, as well as IC and bremsstrahlung components associated with known soft-synchrotron counterparts. We find a significant gamma-ray excess towards the Galactic center with a spectrum that is significantly harder than other sky components and is most consistent with IC from a hard population of electrons. The morphology and spectrum are consistent with it being the IC counterpart to the electrons which generate the microwave haze seen at WMAP frequencies. In addition, the implied electron spectrum is hard; electrons accelerated in supernova shocks in the disk which then diffuse a few kpc to the haze region would have a softer spectrum. We describe the full sky Fermi maps used in this analysis and make them available for download.

We present an analysis of the Sommerfeld enhancement to dark matter annihilation in the presence of an excited state, where the interaction inducing the enhancement is purely off-diagonal, such as in models of exciting or inelastic dark matter. We derive a simple and accurate semi-analytic approximation for the s-wave enhancement, which is valid provided the mass splitting between the ground and excited states is not too large, and discuss the cutoff of the enhancement for large mass splittings. We reproduce previously derived results in the appropriate limits, and demonstrate excellent agreement with numerical calculations of the enhancement. We show that the presence of an excited state leads to generically larger values of the Sommerfeld enhancement, larger resonances, and shifting of the resonances to lower mediator masses. Furthermore, in the presence of a mass splitting the enhancement is no longer a monotonic function of velocity: the enhancement where the kinetic energy is close to that required to excite the higher state can be up to twice as large as the enhancement at zero velocity.

We present the case for a dark matter detector with directional sensitivity. This document was developed at the 2009 CYGNUS workshop on directional dark matter detection, and contains contributions from theorists and experimental groups in the field. We describe the need for a dark matter detector with directional sensitivity; each directional dark matter experiment presents their project's status; and we close with a feasibility study for scaling up to a one ton directional detector, which would cost around $150M.

Recently, a preliminary spectrum from the Fermi Gamma-ray Space Telescope has been presented for the inner galaxy (-30 < l < 30, -5 < b < 5), as well as the galactic center (-1 < l < 1, -1< b < 1). We consider the implications of these data for dark matter annihilation models, especially models capable of producing the cosmic-ray excesses previously observed by PAMELA and Fermi. These local cosmic-ray excesses, when extrapolated to the inner galaxy, imply inverse Compton scattering (ICS) gamma-ray signals largely consistent with the preliminary Fermi gamma-ray spectrum. For specific halos and models, particularly those with prompt photons, the data have begun to constrain the allowed parameter space. Although significant modeling and background uncertainties remain, we explore how large a signal is permitted by the current data. Based upon this, we make predictions for what signal could be present in other regions of the sky where dark matter signals may be easier to isolate from the astrophysical backgrounds.

Models of dark matter with ~ GeV scale force mediators provide attractive explanations of many high energy anomalies, including PAMELA, ATIC, and the WMAP haze. At the same time, by exploiting the ~ MeV scale excited states that are automatically present in such theories, these models naturally explain the DAMA/LIBRA and INTEGRAL signals through the inelastic dark matter (iDM) and exciting dark matter (XDM) scenarios, respectively. Interestingly, with only weak kinetic mixing to hypercharge to mediate decays, the lifetime of excited states with delta < 2 m_e is longer than the age of the universe. The fractional relic abundance of these excited states depends on the temperature of kinetic decoupling, but can be appreciable. There could easily be other mechanisms for rapid decay, but the consequences of such long-lived states are intriguing. We find that CDMS constrains the fractional relic population of ~100 keV states to be <~ 10^-2, for a 1 TeV WIMP with sigma_n = 10^-40 cm^2. Upcoming searches at CDMS, as well as xenon, silicon, and argon targets, can push this limit significantly lower. We also consider the possibility that the DAMA excitation occurs from a metastable state into the XDM state, which decays via e+e- emission, which allows lighter states to explain the INTEGRAL signal due to the small kinetic energies required. Such models yield dramatic signals from down-scattering, with spectra peaking at high energies, sometimes as high as ~1 MeV, well outside the usual search windows. Such signals would be visible at future Ar and Si experiments, and may be visible at Ge and Xe experiments. We also consider other XDM models involving ~ 500 keV metastable states, and find they can allow lighter WIMPs to explain INTEGRAL as well.

We compute in detail the rate at which energy injected by dark matter annihilation heats and ionizes the photon-baryon plasma at z ~ 1000, and provide accurate fitting functions over the relevant redshift range for a broad array of annihilation channels and DM masses. The resulting perturbations to the ionization history can be constrained by measurements of the CMB temperature and polarization angular power spectra. We show that models which fit recently measured excesses in 10-1000 GeV electron and positron cosmic rays are already close to the 95% confidence limits from WMAP. The recently launched Planck satellite will be capable of ruling out a wide range of DM explanations for these excesses. In models of dark matter with Sommerfeld-enhanced annihilation, where sigma v rises with decreasing WIMP velocity until some saturation point, the WMAP5 constraints imply that the enhancement must be close to saturation in the neighborhood of the Earth.

We propose a comprehensive theory of dark matter that explains the recent proliferation of unexpected observations in high-energy astrophysics. Cosmic ray spectra from ATIC and PAMELA require a WIMP with mass M_chi ~ 500 - 800 GeV that annihilates into leptons at a level well above that expected from a thermal relic. Signals from WMAP and EGRET reinforce this interpretation. Taken together, we argue these facts imply the presence of a GeV-scale new force in the dark sector. The long range allows a Sommerfeld enhancement to boost the annihilation cross section as required, without altering the weak scale annihilation cross section during dark matter freezeout in the early universe. If the dark matter annihilates into the new force carrier, phi, its low mass can force it to decay dominantly into leptons. If the force carrier is a non-Abelian gauge boson, the dark matter is part of a multiplet of states, and splittings between these states are naturally generated with size alpha m_phi ~ MeV, leading to the eXciting dark matter (XDM) scenario previously proposed to explain the positron annihilation in the galactic center observed by the INTEGRAL satellite. Somewhat smaller splittings would also be expected, providing a natural source for the parameters of the inelastic dark matter (iDM) explanation for the DAMA annual modulation signal. Since the Sommerfeld enhancement is most significant at low velocities, early dark matter halos at redshift ~10 potentially produce observable effects on the ionization history of the universe, and substructure is more detectable than with a conventional WIMP. Moreover, the low velocity dispersion of dwarf galaxies and Milky Way subhalos can greatly increase the substructure annihilation signal.

We consider the nuclear scattering cross section for the eXciting Dark Matter (XDM) model. In XDM, the Weakly Interacting Massive Particles (WIMPs) couple to the Standard Model only via an intermediate light scalar which mixes with the Higgs: this leads to a suppression in the nuclear scattering cross section relative to models in which the WIMPs couple to the Higgs directly. We estimate this suppression factor to be of order 10^(-5). The elastic nuclear scattering cross section for XDM can also be computed directly: we perform this computation for XDM coupled to the Higgs sector of the Standard Model and find a spin-independent cross section in the order of 4 x 10^(-13) pb in the decoupling limit, which is not within the range of any near-term direct detection experiments. However, if the XDM dark sector is instead coupled to a two-Higgs-doublet model, the spin-independent nuclear scattering cross section can be enhanced by up to four orders of magnitude for large tan(beta), which should be observable in the upcoming SuperCDMS and ton-scale xenon experiments.