Search

We study the detectability of GeV-band gamma-ray polarization with the AMS-02 experiment and a proposed successor, from Galactic and extragalactic sources. Characterizing gamma-ray polarization in this energy range could shed light on gamma-ray emission mechanisms in the sources; physics beyond the Standard Model, such as the presence of axion-like particles (ALPs), could also induce a distinctive energy-dependent polarization signal due to propagation effects in magnetic fields. We present estimates for the minimum detectable polarization from bright sources and the forecast reach for axion-like particles (ALPs). We show that for an example bright extragalactic source (NGC1275), AMS-02 will only have sensitivity to ALP-induced polarization from currently open parameter space if the B-field configuration is unusually favorable to a signal. However, a successor experiment such as AMS-100 would be expected to probe new parameter space even for pessimistic B-field models, with prospects to measure the energy-dependence of such a signal. For Galactic sources, polarization measurements could provide a unique test of scenarios where ALPs induce energy-dependent features in the photon intensity. However, in the absence of a bright transient source (such as a Galactic supernova), the parameter space that would be probed by this approach with ten years of AMS-100 data is already nominally excluded by other experiments, although this conflict may be avoided in specific ALP models., Comment: 18+2 pages, 11+2 figures. Comments welcome

A number of new balloon or space-based gamma-ray observatories have been proposed to close a "MeV gap" in sensitivity to gamma rays in the MeV-GeV energy band. One aspect of the science case for these instruments is their ability to constrain or discover decaying or annihilating dark matter. In this work, we forecast the sensitivity of these instruments for dark matter annihilation or decay for a range of possible Standard Model final states and compare to existing bounds., Comment: 12 pages plus appendix, 8 figures, comments welcome

Descriptions of the Galactic Center using Fermi gamma-ray data have so far modeled the Galactic Center Excess (GCE) as a template with fixed spatial morphology or as a linear combination of such templates. Although these templates are informed by various physical expectations, the morphology of the excess is a priori unknown. For the first time, we describe the GCE using a flexible, non-parametric machine learning model -- the Gaussian process (GP). We assess our model's performance on synthetic data, demonstrating that the model can recover the templates used to generate the data. We then fit the \Fermi data with our model in a single energy bin from 2-20 GeV (leaving a spectral GP analysis of the GCE for future work) using a variety of template models of diffuse gamma-ray emission to quantify our fits' systematic uncertainties associated with diffuse emission modeling. We interpret our best-fit GP in terms of GCE templates consisting of an NFW squared template and a bulge component to determine which bulge models can best describe the fitted GP and to what extent the best-fit GP is described better by an NFW squared template versus a bulge template. The best-fit GP contains morphological features that are typically not associated with traditional GCE studies. These include a localized bright source at around $(\ell,b) = (20^{\circ}, 0^{\circ})$ and a diagonal arm extending Northwest from the Galactic Center. In spite of these novel features, the fitted GP is explained best by a template-based model consisting of the bulge presented in Coleman et al. (2020) and a squared NFW component. Our results suggest that the physical interpretation of the GCE in terms of stellar bulge and NFW-like components is highly sensitive to the assumed morphologies, background models, and the region of the sky used for inference., Comment: 60 pages, 39 figures

Long-range attractive interactions between dark matter particles can significantly enhance their annihilation, particularly at low velocities. This ``Sommerfeld enhancement'' is typically computed by evaluating the deformation of the two-particle wavefunction due to the long-range potential, while ignoring the physics associated with the annihilation, and then scaling the appropriate annihilation matrix elements by factors that depend on the wavefunction in the limit where the particles approach zero relative separation. It has long been recognized that this approach is a valid approximation only in the limit where the annihilation rate is small, and breaks down in the regime where the enhanced annihilation rate approaches the unitarity bound, in which case ignoring the impact of the annihilation physics on the two-particle wavefunction cannot be justified and leads to apparent violations of unitarity. In the case where the physics relevant to annihilation occurs at a parametrically shorter distance scale (higher energy scale) compared with the long-range potential, we provide a simple prescription for correcting the Sommerfeld enhancement for the effects of the short-range physics, valid for all partial waves and for systems where multiple states are coupled by the long-range potential., Comment: 70 pages, 17 figures

The temperature and polarization anisotropies of the cosmic microwave background (CMB) have been used to set constraints on decaying dark matter models down to keV masses. In this work, we extend these limits down to the sub-keV mass range. Using principal component analysis, we estimate the lower bound on the decay lifetime for a basis of different dark matter masses and Standard Model final states, from which the bound on an arbitrary model can be calculated. We validate our principal component analysis using Markov chain Monte Carlo methods and Planck 2018 data. We perform a separate analysis for models decaying into photons below the hydrogen ionization threshold. We demonstrate that for these models, the effect of energy deposition can be captured approximately by a single parameter, but the redshift dependence of the effect is very different from higher-energy injections; in particular, the perturbations to CMB anisotropies are more sensitive to energy deposited around the time of recombination., Comment: 15 pages, 10 figures. Ancillary files contain tables of principal components

We investigate the hypothesis that sexaquarks, hypothetical stable six-quark states, could be a significant component of the dark matter. We expand on previous studies of sexaquark cosmology, accounting for the possibility that some relevant interaction cross sections might be strongly suppressed below expectations based on dimensional analysis. We update direct-detection constraints on stable sexaquarks comprising a subdominant fraction of the dark matter, as well as limits on the annihilation of an antisexaquark component from Super-Kamiokande. We argue that the scenario where sexaquarks comprise a $O(1)$ fraction of the dark matter would require either a suppression of $O(10^{-19})$ in sexaquark interactions with baryons, combined with a very high yield of net sexaquark number from the quark-hadron transition, or else a very strong suppression of the cross section for antisexaquark annihilation on nucleons (24+ orders of magnitude below the QCD scale). Independently, we find that a sexaquark component comprising more than $O(10^{-3})$ of the dark matter can be excluded from direct-detection bounds, unless its scattering cross section is severely suppressed compared to the expected scale of strong and even electromagnetic interactions., Comment: 27 pages + appendices, 14 figures, v2: version accepted for publication

The 21-cm signal provides a novel avenue to measure the thermal state of the universe during cosmic dawn and reionization (redshifts $z\sim 5-30$), and thus to probe energy injection from decaying or annihilating dark matter (DM). These DM processes are inherently inhomogeneous: both decay and annihilation are density dependent, and furthermore the fraction of injected energy that is deposited at each point depends on the gas ionization and density, leading to further anisotropies in absorption and propagation. In this work, we develop a new framework for modeling the impact of spatially inhomogeneous energy injection and deposition during cosmic dawn, accounting for ionization and baryon density dependence, as well as the attenuation of propagating photons. We showcase how this first completely inhomogeneous treatment affects the predicted 21-cm power spectrum in the presence of exotic sources of energy injection, and forecast the constraints that upcoming HERA measurements of the 21-cm power spectrum will set on DM decays to photons and to electron/positron pairs. These projected constraints considerably surpass those derived from CMB and Lyman-$\alpha$ measurements, and for decays to electron/positron pairs they exceed all existing constraints in the sub-GeV mass range, reaching lifetimes of $\sim 10^{28}\,\mathrm{s}$. Our analysis demonstrates the unprecedented sensitivity of 21-cm cosmology to exotic sources of energy injection during the cosmic dark ages. Our code, $\mathtt{DM21cm}$, includes all these effects and is publicly available in an accompanying release., Comment: 33 pages, 22 figures, public code at https://github.com/yitiansun/DM21cm

We extend the Effective Field Theory of Heavy Dark Matter to arbitrary odd representations of SU(2) and incorporate the effects of bound states. This formalism is then deployed to compute the gamma-ray spectrum for a 5 of SU(2): quintuplet dark matter. Except at isolated values of the quintuplet mass, the bound state contribution to hard photons with energy near the dark-matter mass is at the level of a few percent compared to that from direct annihilation. Further, compared to smaller representations, such as the triplet wino, the quintuplet can exhibit a strong variation in the shape of the spectrum as a function of mass. Using our results, we forecast the fate of the thermal quintuplet, which has a mass of $\sim$13.6 TeV. We find that existing H.E.S.S. data should be able to significantly test the scenario, however, the final word on this canonical model of minimal dark matter will likely be left to the Cherenkov Telescope Array (CTA)., Comment: 85 pages, 14 figures; v2: appendix and figures added on rapid variation of annihilation rate with mass

The first stars are expected to form through molecular-hydrogen (H$_2$) cooling, a channel that is especially sensitive to the thermal and ionization state of gas, and can thus act as a probe of exotic energy injection from decaying or annihilating dark matter (DM). Here, we use a toy halo model to study the impact of DM-sourced energy injection on the H$_2$ content of the first galaxies, and thus estimate the threshold mass required for a halo to form stars at high redshifts. We find that currently allowed DM models can significantly change this threshold, producing both positive and negative feedback. In some scenarios, the extra heating of the gas raises the halo mass required for collapse, whereas in others, energy injection lowers the threshold by increasing the free-electron fraction and catalyzing H$_2$ formation. The direction of the effect can be redshift-dependent. We also bracket the uncertainties from self-shielding of halos from Lyman-Werner radiation. Hence, exotic energy injection can both delay and accelerate the onset of star formation; we show how this can impact the timing of 21cm signals at cosmic dawn. We encourage detailed simulation follow-ups in the most promising regions of parameter space identified in this work., Comment: 32 pages, 13 figures, 3 tables

We calculate the post-recombination contribution to the Cosmic Microwave Background (CMB) spectral distortion due to general exotic energy injections, including dark matter (DM) decaying or annihilating to Standard Model particles. Upon subtracting the background distortion that would be present even without such energy injections, we find residual distortions that are still potentially large enough to be detectable by future experiments such as PIXIE. The distortions also have a high-energy spectral feature that is a unique signature of the injection of high-energy particles. We present a calculation of the global ionization history in the presence of decaying dark matter with sub-keV masses, and also show that previous calculations of the global ionization history in the presence of energy injection are not significantly modified by these additional spectral distortions. Our improved treatment of low-energy electrons allows us to extend calculations of the CMB anisotropy constraints for decaying DM down to arbitrarily low masses. We also recast these bounds as constraints on the coupling of axion-like particles to photons., Comment: 23 pages, 16 figures

Decaying or annihilating dark matter and other exotic energy injections can modify the spectrum of the universe's photon bath, resulting in e.g. new contributions to spectral distortions of the cosmic microwave background blackbody spectrum and modifications to the temperature and ionization history of the universe. Here, we present an improved version of the $\texttt{DarkHistory}$ code, which is now capable of consistently calculating the spectrum of low-energy photons by properly treating the interactions of these photons with the levels of hydrogen atoms. Other changes to the code include a more detailed treatment of energy deposition by low-energy electrons, and spectral distortions from heating of the intergalactic medium. All of the improvements we have made to $\texttt{DarkHistory}$ are publicly available., Comment: 29 pages, 12 figures

The fundamental nature of Dark Matter is a central theme of the Snowmass 2021 process, extending across all Frontiers. In the last decade, advances in detector technology, analysis techniques and theoretical modeling have enabled a new generation of experiments and searches while broadening the types of candidates we can pursue. Over the next decade, there is great potential for discoveries that would transform our understanding of dark matter. In the following, we outline a road map for discovery developed in collaboration among the Frontiers. A strong portfolio of experiments that delves deep, searches wide, and harnesses the complementarity between techniques is key to tackling this complicated problem, requiring expertise, results, and planning from all Frontiers of the Snowmass 2021 process., Comment: 10 pages, 5 figures. Version prepared for inclusion in the Snowmass Book. Extended version at arXiv:2210.01770. v2: fixed authors and affiliations

Self-interacting dark matter (SIDM) offers the potential to mitigate some of the discrepancies between simulated cold dark matter (CDM) and observed galactic properties. We introduce a physically motivated SIDM model to understand the effects of self interactions on the properties of Milky Way and dwarf galaxy sized haloes. This model consists of dark matter with a nearly degenerate excited state, which allows for both elastic and inelastic scattering. In particular, the model includes a significant probability for particles to up-scatter from the ground state to the excited state. We simulate a suite of zoom-in Milky Way-sized N-body haloes with six models with different scattering cross sections to study the effects of up-scattering in SIDM models. We find that the up-scattering reaction greatly increases the central densities of the main halo through the loss of kinetic energy. However, the physical model still results in significant coring due to the presence of elastic scattering and down-scattering. These effects are not as apparent in the subhalo population compared to the main halo, but the number of subhaloes is reduced compared to CDM., Comment: 19 pages, 13 figures, accepted to MNRAS

The fundamental nature of Dark Matter is a central theme of the Snowmass 2021 process, extending across all frontiers. In the last decade, advances in detector technology, analysis techniques and theoretical modeling have enabled a new generation of experiments and searches while broadening the types of candidates we can pursue. Over the next decade, there is great potential for discoveries that would transform our understanding of dark matter. In the following, we outline a road map for discovery developed in collaboration among the frontiers. A strong portfolio of experiments that delves deep, searches wide, and harnesses the complementarity between techniques is key to tackling this complicated problem, requiring expertise, results, and planning from all Frontiers of the Snowmass 2021 process., Comment: v3: fixed authorlist (on arXiv only)

This report summarizes the findings of the CF1 Topical Subgroup to Snowmass 2021, which was focused on particle dark matter. One of the most important scientific goals of the next decade is to reveal the nature of dark matter (DM). To accomplish this goal, we must delve deep, to cover high priority targets including weakly-interacting massive particles (WIMPs), and search wide, to explore as much motivated DM parameter space as possible. A diverse, continuous portfolio of experiments at large, medium, and small scales that includes both direct and indirect detection techniques maximizes the probability of discovering particle DM. Detailed calibrations and modeling of signal and background processes are required to make a convincing discovery. In the event that a candidate particle is found through different means, for example at a particle collider, the program described in this report is also essential to show that it is consistent with the actual cosmological DM. The US has a leading role in both direct and indirect detection dark matter experiments -- to maintain this leading role, it is imperative to continue funding major experiments and support a robust R\&D program., Comment: Submitted 30 pages, 11 figures, many references, Report of the CF1 Topical Group for Snowmass 2021

We have performed a systematic study of the statistical behavior of non-Poissonian template fitting (NPTF), a method designed to analyze and characterize unresolved point sources in general counts datasets. In this paper, we focus on the properties and characteristics of the Fermi-LAT gamma-ray data set. In particular, we have simulated and analyzed gamma-ray sky maps under varying conditions of exposure, angular resolution, pixel size, energy window, event selection, and source brightness. We describe how these conditions affect the sensitivity of NPTF to the presence of point sources, for inner-galaxy studies of point sources within the Galactic Center excess, and for the simplified case of isotropic emission. We do not find opportunities for major gains in sensitivity from varying these choices, within the range available with current Fermi-LAT data. We provide an analytic estimate of the NPTF sensitivity to point sources for the case of isotropic emission and perfect angular resolution, and find good agreement with our numerical results for that case., Comment: Published in PRD, 24 pages (14 figures) + appendices (14 figures), comments welcome

We show that Dense Neural Networks can be used to accurately model the cooling of high-energy particles in the early universe, in the context of the public code package DarkHistory. DarkHistory self-consistently computes the temperature and ionization history of the early universe in the presence of exotic energy injections, such as might arise from the annihilation or decay of dark matter. The original version of DarkHistory uses large pre-computed transfer function tables to evolve photon and electron spectra in redshift steps, which require a significant amount of memory and storage space. We present a light version of DarkHistory that makes use of simple Dense Neural Networks to store and interpolate the transfer functions, which performs well on small computers without heavy memory or storage usage. This method anticipates future expansion with additional parametric dependence in the transfer functions without requiring exponentially larger data tables., Comment: 12 pages, 10 figures. Code at https://github.com/hongwanliu/DarkHistory . Data files at https://doi.org/10.5281/zenodo.6819281

A near-future detection of the 21cm signal from the epoch of reionization will provide unique opportunities to probe the underlying cosmology, provided that such cosmological information can be extracted with precision. To this end, we further develop effective field theory (EFT) inspired techniques for the 21cm brightness temperature field during the epoch of reionization, incorporating renormalized bias and a treatment of redshift space distortions. Notably, we confirm that in redshift space, measures of the 21cm brightness, e.g the power spectrum, should have irreducible contributions that lack a bias coefficient and therefore contain direct, astrophysics-free information about the cosmological density field; in this work, we study this effect beyond linear order. To validate our theoretical treatment, we fit the predicted EFT Fourier-space shapes to the THESAN suite of hydrodynamical simulations of reionization at the field level, where the considerable number of modes prevents overfitting. We find agreement at the level of a few percent between the 21cm power spectrum from the EFT fits and simulations over the wavenumber range $k \lesssim 0.8$ h/Mpc and neutral fraction $x_\mathrm{HI} \gtrsim 0.4$, which is imminently measurable by the Hydrogen Epoch of Reionization Array (HERA) and future experiments. The ability of the EFT to describe the 21cm signal extends to simulations that have different astrophysical prescriptions for reionization as well as simulations with interacting dark matter., Comment: 25 pages, 10 figures

It has been shown that a first order confinement phase transition can drastically change the relic dark matter abundance in confining dark sectors with only heavy dark quarks. We study the phenomenology of one such model with a Z' portal to Standard Model. We find that dark glueballs are long-lived in this setup and dilute the dark matter abundance after their decay to Standard Model. With this effect, the correct relic abundance is obtained with dark matter masses up to $\mathcal{O}(10^6)$~TeV. We find that while a part of the parameter space is already ruled out by direct detection and collider searches, there is still a broad space of viable scenarios that can be probed by future experiments., Comment: 10 pages, 8 figures, and three Appendices on thermodynamics and Phase Transitions and heavier Glueball states. Version published in JHEP

Intriguing signals with excesses over expected backgrounds have been observed in many astrophysical and terrestrial settings, which could potentially have a dark matter origin. Astrophysical excesses include the Galactic Center GeV gamma-ray excess detected by the Fermi Gamma-Ray Space Telescope, the AMS antiproton and positron excesses, and the 511 and 3.5 keV X-ray lines. Direct detection excesses include the DAMA/LIBRA annual modulation signal, the XENON1T excess, and low-threshold excesses in solid state detectors. We discuss avenues to resolve these excesses, with actions the field can take over the next several years., Comment: 57 pages, solicited white paper submitted to the Proceedings of the US Community Study on the Future of Particle Physics (Snowmass 2021)

A new population of millisecond pulsars is a long-standing proposed explanation for the excess of GeV-scale gamma rays emanating from the region surrounding the center of the Milky Way (the "Galactic Center excess"). We examine several simple parameterizations of possible luminosity functions for this population, as well as several benchmark luminosity functions proposed in the literature, and compare the predicted populations of resolved point sources to the Fermi 4FGL-DR2 point source catalog and a sub-population recently identified using wavelet-based methods. We provide general results that can be used to translate upper limits on the number of resolved point sources associated with the excess, and the fraction of the flux in the excess that can be attributed to resolved sources, into limits on the luminosity function parameter space. We discuss a number of important systematic uncertainties, including in the detection threshold model and the total flux attributed to the excess. We delineate regions of parameter space (containing existing benchmark models) where there is no apparent tension with current data, and the number of total pulsars needed to explain the excess is in the range of O(10,000-100,000). Forecasting the effects of lowered point source detection thresholds, we show that novel analysis methods that probe sub-threshold point source populations can hope to resolve more than 30% of the flux of the excess., Comment: Published in JCAP; updated to fix minor typos. 32 pages plus appendices

We study general freeze-out scenarios where an arbitrary number of initial and final dark matter particles participate in the number-changing freeze-out interaction. We consider a simple sector with two particle species undergoing such a thermal freeze-out; one of the relics is stable and gives rise to the dark matter today, while the other one decays to the Standard Model, injecting significant entropy into the thermal bath that dilutes the dark matter abundance. We show that this setup can lead to a stable relic population that reproduces the observed dark matter abundance without requiring weak scale masses or couplings. The final dark matter abundance is estimated analytically. We carry out this calculation for arbitrary temperature dependence in the freeze-out process and identify the viable dark matter mass and cross section ranges that explain the observed dark matter abundance. This setup can be used to open parameter space for both heavy (above the unitarity bound) or light (sub-GeV) dark matter candidates. We point out that the best strategy for probing most parts of our parameter space is to look for signatures of an early matter-dominant epoch., Comment: 9 pages and technical appendices

These lectures, presented at the 2021 Les Houches Summer School on Dark Matter, provide an introduction to key methods and tools of indirect dark matter searches, as well as a status report on the field circa summer 2021. Topics covered include the possible effects of energy injection from dark matter on the early universe, methods to calculate both the expected energy distribution and spatial distribution of particles produced by dark matter interactions, an outline of theoretical models that predict diverse signals in indirect detection, and a discussion of current constraints and some claimed anomalies. These notes are intended as an introduction to indirect dark matter searches for graduate students, focusing primarily on intuition-building estimates and useful concepts and tools., Comment: Submitted to SciPost Physics Lecture Notes, Les Houches Summer School Series. v2 fixes typos, clarifies a few points, and adds one figure. v3 adds a reference accidentally omitted from v2. arXiv admin note: substantial text overlap with arXiv:1710.05137

We examine the vector-portal inelastic dark matter (DM) model with DM mass $m_\chi$ and dark photon mass $m_{A'}$, in the `forbidden dark matter' regime where $1 \lesssim m_{A'}/m_\chi \lesssim 2$, carefully tracking the dark sector temperature throughout freezeout. The inelastic nature of the dark sector relaxes the stringent cosmic microwave background (CMB) and self-interaction constraints compared to symmetric DM models. We determine the CMB limits on both annihilations involving excited states and annihilation into $e^+e^-$ through initial-state-radiation of an $A'$, as well as limits on the DM self-scattering, which proceeds at the one-loop level. The unconstrained parameter space serves as an ideal target for accelerator $A'$ searches, and provides a DM self-interaction cross section that is large enough to observably impact small-scale structure., Comment: 5 pages, 2 figures + supplemental material 8 pages, 3 figures

We carry out a detailed study of the confinement phase transition in a dark sector with a $SU(N)$ gauge group and a single generation of dark heavy quark. We focus on heavy enough quarks such that their abundance freezes out before the phase transition and the phase transition is of first-order. We find that during this phase transition the quarks are trapped inside contracting pockets of the deconfined phase and are compressed enough to interact at a significant rate, giving rise to a second stage of annihilation that can dramatically change the resulting dark matter abundance. As a result, the dark matter can be heavier than the often-quoted unitarity bound of $\sim100~$TeV. Our findings are almost completely independent of the details of the portal between the dark sector and the Standard Model. We comment briefly on possible signals of such a sector. Our main findings are summarized in a companion letter, while here we provide further details on different parts of the calculation., Comment: 38+30 pages, 14 figures

We study the effect of a first-order phase transition in a confining $SU(N)$ dark sector with heavy dark quarks. The baryons of this sector are the dark matter candidate. During the confinement phase transition the heavy quarks are trapped inside isolated, contracting pockets of the deconfined phase, giving rise to a second stage of annihilation that dramatically suppresses the dark quark abundance. The surviving abundance is determined by the local accidental asymmetry in each pocket. The correct dark matter abundance is obtained for $\mathcal{O}(1-100)$ PeV dark quarks, above the usual unitarity bound., Comment: 5 pages, 3 figures, comments are welcomed

We fully explore the thermal freezeout histories of a vector-portal dark matter model, in the region of parameter space in which the ratio of masses of the dark photon $A^{\prime}$ and dark matter $\chi$ is in the range $1 \lesssim m_{A^{\prime}}/m_{\chi} \lesssim 2$. In this region $2 \rightarrow 2$ and $3 \rightarrow 2$ annihilation processes within the dark sector, as well as processes that transfer energy between the dark sector and the Standard Model, play important roles in controlling the thermal freezeout of the dark matter. We carefully track the temperatures of all species, relaxing the assumption of previous studies that the dark and Standard Model sectors remain in thermal equilibrium throughout dark matter freezeout. Our calculations reveal a rich set of novel pathways which lead to the observed relic density of dark matter, and we develop a simple analytic understanding of these different regimes. The viable parameter space in our model provides a target for future experiments searching for light (MeV-GeV) dark matter, and includes regions where the dark matter self-interaction cross section is large enough to affect the small-scale structure of galaxies., Comment: 36 pages, 18 figures, comments welcome

TeV-scale particles that couple to the standard model through the weak force represent a compelling class of dark matter candidates. The search for such Weakly Interacting Massive Particles (WIMPs) has already spanned multiple decades, and whilst it has yet to provide any definitive evidence for their existence, viable parameter space remains. In this paper, we show that the upcoming Cherenkov Telescope Array (CTA) has significant sensitivity to uncharted parameter space at the TeV mass scale. To do so, we focus on two prototypical dark matter candidates, the Wino and Higgsino. Sensitivity forecasts for both models are performed including the irreducible background from misidentified cosmic rays, as well as a range of estimates for the Galactic emissions at TeV energies. For each candidate, we find substantial expected improvements over existing bounds from current imaging atmospheric Cherenkov telescopes. In detail, for the Wino we find a sensitivity improvement of roughly an order of magnitude in $\langle \sigma v \rangle$, whereas for the Higgsino we demonstrate that CTA has the potential to become the first experiment that has sensitivity to the thermal candidate. Taken together, these enhanced sensitivities demonstrates the discovery potential for dark matter at CTA in the 1-100 TeV mass range., Comment: 20 pages, 8 figures, Accepted for publication in Phys. Rev. D

We derive new constraints on models of decaying and annihilating dark matter (DM) by requiring that the energy injected into the intergalactic medium (IGM) not overheat it at late times, when measurements of the Lyman-$\alpha$ forest constrain the IGM temperature. We improve upon previous analyses by using the recently developed $\texttt{DarkHistory}$ code package, which self-consistently takes into account additional photoionization and photoheating processes due to reionization and DM sources. Our constraints are robust to the uncertainties of reionization and competitive with leading limits on sub-GeV DM that decays preferentially to electrons., Comment: 5+7 pages, 7 figures. v2 updates p-wave constraints, fixes typo in Eq. (4), fixes minor Fig. 1 error, and includes pion and muon constraints in new section of Appendix

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., Comment: 9 pages, 5 figures. Updated to match accepted version

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., Comment: 47 pages, 14 figures, published in PRD. See also companion PRL paper arXiv:2002.12370

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., Comment: 6+6 pages, 3+3 figures, published in PRL. See also companion PRD paper arXiv:2002.12371

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., Comment: 7 pages, 4 figures, 36th International Cosmic Ray Conference

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., Comment: 9+4 pages, 4+4 figures. Comments are welcome

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 ., Comment: 40 pages, 17 figures

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., Comment: 7+32 pages, 2+20 figures

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)., Comment: Astro2020 Decadal Survey science white paper; list of endorsers available at https://github.com/veragluscevic/Astro2020-DM-Cosmology-Endorsers

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., Comment: 11 pages, 2 figures, Science Whitepaper for Astro 2020, more information at https://lsstdarkmatter.github.io

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/ ., Comment: 96 pages, 22 figures, 1 table

We study the bound-state spectrum in a simple model of pseudo-Dirac dark matter, and examine how the rate of bound-state formation through radiative capture compares to Sommerfeld-enhanced annihilation. We use this model as an example to delineate the new features induced by the presence of a mass splitting between the dark matter and a nearly-degenerate partner, compared to the case where only a single dark-matter-like state is present. We provide a simple analytic prescription for estimating the spectrum of bound states in systems containing a mass splitting, which in turn allows characterization of the resonances due to near-zero-energy bound states, and validate this estimate both for pseudo-Dirac dark matter and for the more complex case of wino dark matter. We demonstrate that for pseudo-Dirac dark matter the capture rate into deeply bound states is, to a good approximation, simply related to the Sommerfeld enhancement factor., Comment: 26 pages, 8 figures, Comments welcome

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., Comment: 61 pages, 16 figures

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., Comment: 56 pages, 12 figures; v2, updated to JHEP version

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., Comment: 18 pages, 14 figures

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., Comment: 32 pages, 8 figures, 331 references

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., Comment: 16 pages, 10 figures, references added, accepted for publication in PRD

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., Comment: 16 pages, 12 figures, MNRAS accepted

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., Comment: 22 pages with appendices, 12 figures, comments welcome; v2: references added with comments, typos corrected, minor change to millicharged DM limits