Search

Physical review letters 130(18), 181801 (2023). doi:10.1103/PhysRevLett.130.181801, We study the signal and background that arise in nuclear magnetic resonance searches for axion dark matter, finding key differences with the existing literature. We find that spin-precession instruments are much more sensitive than what has been previously estimated in a sizable range of axion masses, with sensitivity improvement of up to a factor of 100 using a Xe129 sample. This improves the detection prospects for the QCD axion, and we estimate the experimental requirements to reach this motivated target. Our results apply to both the axion electric and magnetic dipole moment operators., Published by APS, College Park, Md.

Searches for dark matter decaying into photons constrain its lifetime to be many orders of magnitude larger than the age of the Universe. A corollary statement is that the abundance of any particle that can decay into photons over cosmological timescales is constrained to be much smaller than the cold dark-matter density. We show that an irreducible freeze-in contribution to the relic density of axions is in violation of that statement in a large portion of the parameter space. This allows us to set stringent constraints on axions in the mass range 100 eV–100 MeV. At 10 keV our constraint on a photophilic axion is gaγγ≲8.1×10-14 GeV-1, almost 3 orders of magnitude stronger than the bounds established using horizontal branch stars; at 100 keV our constraint on a photophobic axion coupled to electrons is gaee≲8.0×10-15, almost 4 orders of magnitude stronger than the present results. Although we focus on axions, our argument is more general and can be extended to, for instance, sterile neutrinos. Searches for dark matter decaying into photons constrain its lifetime to be many orders of magnitude larger than the age of the Universe. A corollary statement is that the abundance of any particle that can decay into photons over cosmological timescales is constrained to be much smaller than the cold dark-matter density. We show that an $\textit{irreducible}$ freeze-in contribution to the relic density of axions is in violation of that statement in a large portion of the parameter space. This allows us to set stringent constraints on axions in the mass range $100\rm \;eV - 100\; MeV$. At $10\rm \; keV$ our constraint on a photophilic axion is $g_{a\gamma \gamma} \lesssim 8.1 \times 10^{-14}~{\rm GeV}^{-1}$, almost three orders of magnitude stronger than the bounds established using horizontal branch stars; at $100~{\rm keV}$ our constraint on a photophobic axion coupled to electrons is $g_{aee} \lesssim 8.0 \times 10^{-15}$, almost four orders of magnitude stronger than present results. Although we focus on axions, our argument is more general and can be extended to, for instance, sterile neutrinos.

While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligencebased statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade. While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade. While astrophysical and cosmological probes provide a remarkably precise and consistent picture of the quantity and general properties of dark matter, its fundamental nature remains one of the most significant open questions in physics. Obtaining a more comprehensive understanding of dark matter within the next decade will require overcoming a number of theoretical challenges: the groundwork for these strides is being laid now, yet much remains to be done. Chief among the upcoming challenges is establishing the theoretical foundation needed to harness the full potential of new observables in the astrophysical and cosmological domains, spanning the early Universe to the inner portions of galaxies and the stars therein. Identifying the nature of dark matter will also entail repurposing and implementing a wide range of theoretical techniques from outside the typical toolkit of astrophysics, ranging from effective field theory to the dramatically evolving world of machine learning and artificial-intelligence-based statistical inference. Through this work, the theory frontier will be at the heart of dark matter discoveries in the upcoming decade.

Physical review letters 129(4), 041101 (2022). doi:10.1103/PhysRevLett.129.041101, Gravitational waves (GWs) generate oscillating electromagnetic effects in the vicinity of external electric and magnetic fields. We discuss this phenomenon with a particular focus on reinterpreting the results of axion haloscopes based on lumped-element detectors, which probe GWs in the 100 kHz–100 MHz range. Measurements from ABRACADABRA and SHAFT already place bounds on GWs, although the present strain sensitivity is weak. However, we demonstrate that the sensitivity scaling with the volume of such instruments is significant—faster than for axions—and so rapid progress will be made in the future. With no modifications, DMRadio-m$^3$ will have a GW strain sensitivity of h∼10$^{-20}$ at 200 MHz. A simple modification of the pickup loop used to readout the induced magnetic flux can parametrically enhance the GW sensitivity, particularly at lower frequencies., Published by APS, College Park, Md.

We compute the decay spectrum for dark matter (DM) with masses above the scale of electroweak symmetry breaking, all the way to the Planck scale. For an arbitrary hard process involving a decay to the unbroken standard model, we determine the prompt distribution of stable states including photons, neutrinos, positrons, and antiprotons. These spectra are a crucial ingredient in the search for DM via indirect detection at the highest energies as being probed in current and upcoming experiments including IceCube, HAWC, CTA, and LHAASO. Our approach improves considerably on existing methods. For example, we include all relevant electroweak interactions. The importance of these effects grow with DM mass, and by an EeV our spectra can differ by orders of magnitude from existing results., 8+24 pages, 3+12 figures, updated to published version. Spectra available at https://github.com/nickrodd/HDMSpectra

Indirect detection opens a unique window for probing thermal dark matter (DM): the same annihilation process that determined the relic abundance in the early Universe drives the present day astrophysical signal. While TeV-scale particles weakly coupled to the Standard Model face undoubted challenges from decades of null searches, the scenario remains compelling, and simple realizations such as Higgsino DM remain largely unexplored. The fate of such scenarios could be determined by gamma-ray observations of the centre of the Milky Way with Imaging Atmospheric Cherenkov Telescopes (IACTs). We consider the ultimate sensitivity of current IACTs to a broad range of TeV-scale DM candidates - including specific ones such as the Wino, Higgsino, and Quintuplet. To do so, we use realistic mock H.E.S.S.-like observations of the inner Milky Way halo, and provide a careful assessment of the impact of recent Milky Way mass modeling, instrumental and astrophysical background uncertainties in the Galactic Center region, and the theoretical uncertainty on the predicted signal. We find that the dominant systematic for IACT searches in the inner Galaxy is the unknown distribution of DM in that region, however, beyond this the searches are currently statistically dominated indicating a continued benefit from more observations. For two-body final states at $1~{\rm TeV}$, we find a H.E.S.S.-like observatory is sensitive to $\langle \sigma v \rangle \sim 3 \times 10^{-26}-4 \times 10^{-25}~{\rm cm}^3{\rm s}^{-1}$, except for neutrino final states, although we find results competitive with ANTARES. In addition, the thermal masses for the Wino and Quintuplet can be probed; the Higgsino continues to be out of reach by at least a factor of a few. Our conclusions are also directly relevant to the next generation Cherenkov Telescope Array, which remains well positioned to be the discovery instrument for thermal DM., Comment: 23 pages, 13 figures, 3 tables, including appendix. Accepted in Phys. Rev. D

Existing searches for cosmic axions relics have relied heavily on the axion being nonrelativistic and making up dark matter. However, light axions can be copiously produced in the early Universe and remain relativistic today, thereby constituting a Cosmic axion Background (CaB). As prototypical examples of axion sources, we consider thermal production, dark-matter decay, parametric resonance, and topological defect decay. Each of these has a characteristic frequency spectrum that can be searched for in axion direct detection experiments. We focus on the axion-photon coupling and study the sensitivity of current and future versions of ADMX, HAYSTAC, DMRadio, and ABRACADABRA to a CaB, finding that the data collected in search of dark matter can be repurposed to detect axion energy densities well below limits set by measurements of the energy budget of the Universe. In this way, direct detection of relativistic relics offers a powerful new opportunity to learn about the early Universe and, potentially, discover the axion.

The identification and description of point sources is one of the oldest problems in astronomy; yet, even today the correct statistical treatment for point sources remains one of the field's hardest problems. For dim or crowded sources, likelihood based inference methods are required to estimate the uncertainty on the characteristics of the source population. In this work, a new parametric likelihood is constructed for this problem using Compound Poisson Generator (CPG) functionals which incorporate instrumental effects from first principles. We demonstrate that the CPG approach exhibits a number of advantages over Non-Poissonian Template Fitting (NPTF) - an existing method - in a series of test scenarios in the context of X-ray astronomy. These demonstrations show that the effect of the point-spread function, effective area, and choice of point-source spatial distribution cannot, generally, be factorised as they are in NPTF, while the new CPG construction is validated in these scenarios. Separately, an examination of the diffuse-flux emission limit is used to show that most simple choices of priors on the standard parameterisation of the population model can result in unexpected biases: when a model comprising both a point-source population and diffuse component is applied to this limit, nearly all observed flux will be assigned to either the population or to the diffuse component. A new parametrisation is presented for these priors which properly estimates the uncertainties in this limit. In this choice of priors, CPG correctly identifies that the fraction of flux assigned to the population model cannot be constrained by the data., 30+10 pages, 15 figures, comments welcome; v2 updated to published version. Accompanying code available at https://github.com/ghcollin/cpg_likelihood

The next generation of ultralight dark matter (DM) direct detection experiments, which could confirm sub-eV bosons as the dominant source of DM, will feature multiple detectors operating at various terrestrial locations. As a result of the wave-like nature of ultralight DM, spatially separated detectors will each measure a unique DM phase. When the separation between experiments is comparable to the DM coherence length, the spatially-varying phase contains information beyond that which is accessible at a single detector. We introduce a formalism to extract this information, which performs interferometry directly on the DM wave. In particular, we develop a likelihood-based framework that combines data from multiple experiments to constrain directional information about the DM phase space distribution. We show that the signal in multiple detectors is subject to a daily modulation effect unique to wave-like DM. Leveraging daily modulation, we illustrate that within days of an initial discovery multiple detectors acting in unison could localize directional parameters of the DM velocity distribution such as the direction of the solar velocity to sub-degree accuracy, or the direction of a putative cold DM stream to the sub-arcminute level. We outline how to optimize the locations of multiple detectors with either resonant cavity (such as ADMX or HAYSTAC) or quasistatic (such as ABRACADABRA or DM-Radio) readouts to have maximal sensitivity to the full 3-dimensional DM velocity distribution., 25 pages, 9 figures, animation available at https://github.com/joshwfoster/DM_Interferometry, updated to published version

Sterile neutrinos with masses in the keV range are well-motivated extensions to the Standard Model that could explain the observed neutrino masses while also making up the dark matter (DM) of the Universe. If sterile neutrinos are DM then they may slowly decay into active neutrinos and photons, giving rise to the possibility of their detection through narrow spectral features in astrophysical X-ray data sets. In this work, we perform the most sensitive search to date for this and other decaying DM scenarios across the mass range from 5 to 16 keV using archival XMM-Newton data. We reduce 547 Ms of data from both the MOS and PN instruments using observations taken across the full sky and then use this data to search for evidence of DM decay in the ambient halo of the Milky Way. We determine the instrumental and astrophysical baselines with data taken far away from the Galactic Center, and use Gaussian Process modeling to capture additional continuum background contributions. No evidence is found for unassociated X-ray lines, leading us to produce the strongest constraints to date on decaying DM in this mass range., 8+29 pages, 3+31 figures, Supplementary Data at https://github.com/bsafdi/XMM_BSO_DATA ; v2 additional cross-checks, expanded discussion

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.

A fundamental question regarding the Galactic Center excess (GCE) is whether the underlying structure is pointlike or smooth, often framed in terms of a millisecond pulsar or annihilating dark matter (DM) origin for the emission. We show that Bayesian neural networks (NNs) have the potential to resolve this debate. In simulated data, the method is able to predict the flux fractions from inner Galaxy emission components to on average $\ensuremath{\sim}0.5%$. When applied to the Fermi photon-count map, the NN identifies a smooth GCE in the data, suggestive of the presence of DM, with the estimates for the background templates being consistent with existing results.

The Fermi Large Area Telescope has observed an excess of ~GeV energy gamma rays from the center of the Milky Way, which may arise from near-thermal dark matter annihilation. Firmly establishing the dark matter origin for this excess is however complicated by challenges in modeling diffuse cosmic-ray foregrounds as well as unresolved astrophysical sources, such as millisecond pulsars. Non-Poissonian Template Fitting (NPTF) is one statistical technique that has previously been used to show that at least some fraction of the GeV excess is likely due to a population of dim point sources. These results were recently called into question by Leane and Slatyer (2019), who showed that a synthetic dark matter annihilation signal injected on top of the real Fermi data is not recovered by the NPTF procedure. In this work, we perform a dedicated study of the Fermi data and explicitly show that the central result of Leane and Slatyer (2019) is likely driven by the fact that their choice of model for the Galactic foreground emission does not provide a sufficiently good description of the data. We repeat the NPTF analyses using a state-of-the-art model for diffuse gamma-ray emission in the Milky Way and introduce a novel statistical procedure, based on spherical-harmonic marginalization, to provide an improved description of the Galactic diffuse emission in a data-driven fashion. With these improvements, we find that the NPTF results continue to robustly favor the interpretation that the Galactic Center excess is due, in part, to unresolved astrophysical point sources across the analysis variations that we have explored., Comment: 27 pages, 18 figures, comments welcome

We use unitarity and analyticity of scattering amplitudes to constrain fermionic operators in the standard model effective field theory. For four-fermion operators at mass dimension 8, we scatter flavor superpositions in fixed standard model representations and find the Wilson coefficients to be constrained so that their contraction with any pair of pure density matrices is positive. These constraints imply that flavor-violating couplings are upper-bounded by their flavor-conserving cousins. For instance, LEP data already appears to preclude certain operators in upcoming $\mu \to 3e$ measurements., Comment: 7 pages, 1 figure; v3 corrected a sign in the definition of several operators

The Galactic Center Excess (GCE) of GeV gamma rays can be explained as a signal of annihilating dark matter or of emission from unresolved astrophysical sources, such as millisecond pulsars. Evidence for the latter is provided by a statistical procedure---referred to as Non-Poissonian Template Fitting (NPTF)---that distinguishes the smooth distribution of photons expected for dark matter annihilation from a "clumpy" photon distribution expected for point sources. In this paper, we perform an extensive study of the NPTF on simulated data, exploring its ability to recover the flux and luminosity function of unresolved sources at the Galactic Center. When astrophysical background emission is perfectly modeled, we find that the NPTF successfully distinguishes between the dark matter and point source hypotheses when either component makes up the entirety of the GCE. When the GCE is a mixture of dark matter and point sources, the NPTF may fail to reconstruct the correct contribution of each component. We further study the impact of mismodeling the Galactic diffuse backgrounds, finding that while a dark matter signal could be attributed to point sources in some outlying cases for the scenarios we consider, the significance of a true point source signal remains robust. Our work enables us to comment on a recent study by Leane and Slatyer (2019) that questions prior NPTF conclusions because the method does not recover an artificial dark matter signal injected on actual Fermi data. We demonstrate that the failure of the NPTF to extract an artificial dark matter signal can be natural when point sources are present in the data---with the effect further exacerbated by the presence of diffuse mismodeling---and does not on its own invalidate the conclusions of the NPTF analysis in the Inner Galaxy.

Indirect searches for dark matter (DM) have conventionally been applied to the products of DM annihilation or decay. If DM couples to light force carriers, however, it can be captured into bound states via dissipation of energy that may yield detectable signals. We extend the indirect searches to DM bound state formation and transitions between bound levels, and constrain the emission of unstable dark photons. Our results significantly refine the predicted signal flux that could be observed in experiments. As a concrete example, we use Fermi-LAT dwarf spheroidal observations to obtain constraints in terms of the dark photon mass and energy which we use to search for the formation of stable or unstable bound states., SCOPUS: ar.j, info:eu-repo/semantics/published

The dark matter explanation of the 3.5 keV line is strongly disfavored by our work in Dessert et al. 2020. Boyarsky et al. 2020 questions that conclusion: modeling additional background lines is claimed to weaken the limit sufficiently to re-allow a dark matter interpretation. We respond as follows. 1) A more conservative limit is obtained by modeling additional lines; this point appeared in its entirety in our work in Dessert et al., though we also showed that the inclusion of such lines is not necessary. 2) Despite suggestions in Boyarsky et al., even the more conservative limits strongly disfavor a decaying dark matter origin of the 3.5 keV line., Comment: 7 pages, 1 figure, response to technical comment solicited by Science. Version published in PDU. Full analysis available at https://github.com/bsafdi/BlankSkyfor3p5

The Fermi-LAT Collaboration recently presented a new catalog of gamma-ray sources located within the 40 ° × 40 ° 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 σ 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 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 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. In a spirit of maximal openness and transparency, we have made our analysis code available at https://github.com/bsafdi/GCE-2FIG .

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

Observations of nearby galaxies and galaxy clusters have reported an unexpected X-ray emission line around 3.5 kilo-electron volts (keV). Proposals to explain this line include decaying dark matter$-$in particular, that the decay of sterile neutrinos with a mass around 7 keV could match the available data. If this interpretation is correct, the 3.5 keV line should also be emitted by dark matter in the halo of the Milky Way. We used more than 30 megaseconds of XMM-Newton (X-ray Multi-Mirror Mission) blank-sky observations to test this hypothesis, finding no evidence of the 3.5-keV line emission from the Milky Way halo. We set an upper limit on the decay rate of dark matter in this mass range, which is inconsistent with the possibility that the 3.5-keV line originates from dark matter decay., Comment: v1, 6+13 pages, 3 + 11 figures; v3, updated to version published in Science. Supplementary data at https://github.com/nickrodd/XMM-DM and supplementary analysis examples at https://github.com/bsafdi/BlankSkyfor3p5

The axion is a promising dark matter candidate, which was originally proposed to solve the strong-CP problem in particle physics. To date, the available parameter space for axion and axion-like particle dark matter is relatively unexplored, particularly at masses $m_a\lesssim1\,\mu$eV. ABRACADABRA is a new experimental program to search for axion dark matter over a broad range of masses, $10^{-12}\lesssim m_a\lesssim10^{-6}$ eV. ABRACADABRA-10 cm is a small-scale prototype for a future detector that could be sensitive to the QCD axion. In this Letter, we present the first results from a 1 month search for axions with ABRACADABRA-10 cm. We find no evidence for axion-like cosmic dark matter and set 95% C.L. upper limits on the axion-photon coupling between $g_{a\gamma\gamma}, Comment: 7 pages, 3 Figures. Submitted to PRL

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

The next generation of axion direct detection experiments may rule out or confirm axions as the dominant source of dark matter. We develop a general likelihood-based framework for studying the time-series data at such experiments, with a focus on the role of dark-matter astrophysics, to search for signatures of the QCD axion or axion like particles. We illustrate how in the event of a detection the likelihood framework may be used to extract measures of the local dark matter phase-space distribution, accounting for effects such as annual modulation and gravitational focusing, which is the perturbation to the dark matter phase-space distribution by the gravitational field of the Sun. Moreover, we show how potential dark matter substructure, such as cold dark matter streams or a thick dark disk, could impact the signal. For example, we find that when the bulk dark matter halo is detected at 5$\sigma$ global significance, the unique time-dependent features imprinted by the dark matter component of the Sagittarius stream, even if only a few percent of the local dark matter density, may be detectable at $\sim$2$\sigma$ significance. A co-rotating dark disk, with lag speed $\sim$50 km$/$s, that is $\sim$20$\%$ of the local DM density could dominate the signal, while colder but as-of-yet unknown substructure may be even more important. Our likelihood formalism, and the results derived with it, are generally applicable to any time-series based approach to axion direct detection., Comment: 34 pages, 12 figures, code available at: https://github.com/bsafdi/AxiScan v2, updated to PRD version v3, minor improvements

Author(s): Keeley, Ryan E; Abazajian, Kevork N; Kwa, Anna; Rodd, Nicholas L; Safdi, Benjamin R | Abstract: The Milky Way's Galactic Center harbors a gamma-ray excess that is a candidate signal of annihilating dark matter. Dwarf galaxies remain predominantly dark in their expected commensurate emission. In this work we quantify the degree of consistency between these two observations through a joint likelihood analysis. In doing so we incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center Extended gamma-ray excess (GCE) detected by the Fermi Gamma-Ray Space Telescope. The preferred range of annihilation rates and masses expands when including these unknowns. Even so, using two recent determinations of the Milky Way halo's local density leave the GCE preferred region of single-channel dark matter annihilation models to be in strong tension with annihilation searches in combined dwarf galaxy analyses. A third, higher Milky Way density determination, alleviates this tension. Our joint likelihood analysis allows us to quantify this inconsistency. We provide a set of tools for testing dark matter annihilation models' consistency within this combined dataset. As an example, we test a representative inverse Compton sourced self-interacting dark matter model, which is consistent with both the GCE and dwarfs.

The High Altitude Water Cherenkov (HAWC) gamma-ray observatory is a wide field-of-view observatory sensitive to 500 GeV - 100 TeV gamma rays and cosmic rays. With its observations over 2/3 of the sky every day, the HAWC observatory is sensitive to a wide variety of astrophysical sources, including possible gamma rays from dark matter. Dark matter annihilation and decay in the Milky Way Galaxy should produce gamma-ray signals across many degrees on the sky. The HAWC instantaneous field-of-view of 2 sr enables observations of extended regions on the sky, such as those from dark matter in the Galactic halo. Here we show limits on the dark matter annihilation cross-section and decay lifetime from HAWC observations of the Galactic halo with 15 months of data. These are some of the most robust limits on TeV and PeV dark matter, largely insensitive to the dark matter morphology. These limits begin to constrain models in which PeV IceCube neutrinos are explained by dark matter which primarily decays into hadrons., 24 pages, 11 figures

We propose two new analysis strategies for studying charm and beauty quarks at colliders. The first strategy is aimed at testing the kinematics of heavy-flavor quarks within an identified jet. Here, we use the SoftDrop jet-declustering algorithm to identify two subjets within a large-radius jet, using subjet flavor tagging to test the heavy-quark splitting functions of QCD. For subjets containing a $J / \psi$ or $\Upsilon$, this declustering technique can also help probe the mechanism for quarkonium production. The second strategy is aimed at isolating heavy-flavor production from gluon splitting. Here, we introduce a new FlavorCone algorithm, which smoothly interpolates from well-separated heavy-quark jets to the gluon-splitting regime where jets overlap. Because of its excellent ability to identify charm and beauty hadrons, the LHCb detector is ideally suited to pursue these strategies, though similar measurements should also be possible at ATLAS and CMS. Together, these SoftDrop and FlavorCone studies should clarify a number of aspects of heavy-flavor physics at colliders, and provide crucial information needed to improve heavy-flavor modeling in parton-shower generators., Comment: 22 pages, 14 figures; v2: updated figures with new z_tag condition, references and discussion added

Dark matter in the halos surrounding galaxy groups and clusters can annihilate to high-energy photons. Recent advancements in the construction of galaxy group catalogs provide many thousands of potential extragalactic targets for dark matter. In this paper, we outline a procedure to infer the dark matter signal associated with a given galaxy group. Applying this procedure to a catalog of sources, one can create a full-sky map of the brightest extragalactic dark matter targets in the nearby Universe ($z\lesssim 0.03$), supplementing sources of dark matter annihilation from within the Local Group. As with searches for dark matter in dwarf galaxies, these extragalactic targets can be stacked together to enhance the signals associated with dark matter. We validate this procedure on mock $\textit{Fermi}$ gamma-ray data sets using a galaxy catalog constructed from the $\texttt{DarkSky}$ $N$-body cosmological simulation and demonstrate that the limits are robust, at $\mathcal{O}(1)$ levels, to systematic uncertainties on halo mass and concentration. We also quantify other sources of systematic uncertainty arising from the analysis and modeling assumptions. Our results suggest that a stacking analysis using galaxy group catalogs provides a powerful opportunity to discover extragalactic dark matter and complements existing studies of Milky Way dwarf galaxies., 17+7 pages, 9+4 figures; v2, updated to PRD version with several updates, results and conclusions unchanged

We use 413 weeks of publicly-available $\textit{Fermi}$ Pass 8 gamma-ray data, combined with recently-developed galaxy group catalogs, to search for evidence of dark matter annihilation in extragalactic halos. In our study, we use luminosity-based mass estimates and mass-to-concentration relations to infer the $J$-factors and associated uncertainties for hundreds of galaxy groups within a redshift range $z \lesssim 0.03$. We employ a conservative substructure boost-factor model, which only enhances the sensitivity by an $\mathcal{O}(1)$ factor. No significant evidence for dark matter annihilation is found and we exclude thermal relic cross sections for dark matter masses below $\sim$30 GeV to 95% confidence in the $b\bar{b}$ annihilation channel. These bounds are comparable to those from Milky Way dwarf spheroidal satellite galaxies. The results of our analysis increase the tension, but do not rule out, the dark matter interpretation of the Galactic Center excess. We provide a catalog of the galaxy groups used in this study and their inferred properties, which can be broadly applied to searches for extragalactic dark matter., 5+18 pages, 1+14 figures, catalog available at: https://github.com/bsafdi/DMCat; v2 updated to journal version with several updates, results and conclusions unchanged

Utilizing the Fermi measurement of the gamma-ray spectrum toward the Galactic Center, we derive some of the strongest constraints to date on the dark matter (DM) lifetime in the mass range from hundreds of MeV to above an EeV. Our profile-likelihood based analysis relies on 413 weeks of Fermi Pass 8 data from 200 MeV to 2 TeV, along with up-to-date models for diffuse gamma-ray emission within the Milky Way. We model Galactic and extragalactic DM decay and include contributions to the DM-induced gamma-ray flux resulting from both primary emission and inverse-Compton scattering of primary electrons and positrons. For the extragalactic flux, we also calculate the spectrum associated with cascades of high-energy gamma-rays scattering off of the cosmic background radiation. We argue that a decaying DM interpretation for the 10 TeV-1 PeV neutrino flux observed by IceCube is disfavored by our constraints. Our results also challenge a decaying DM explanation of the AMS-02 positron flux. We interpret the results in terms of individual final states and in the context of simplified scenarios such as a hidden-sector glueball model., Comment: Updated to correct an issue in Fig. S12, 9+19 pages, 2+14 figures, Supplementary Data available at https://dspace.mit.edu/handle/1721.1/105550

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., Comment: 27 pages, 5 figures. Corrected counter-term contributions, fixing an issue first noted in 1805.07367. Impact on final results in plots is imperceptible, leaving the conclusions and numerical analysis unchanged

The standard quantization formalism in spacetimes with event horizons implies a non-unitary evolution of quantum states, as initial pure states may evolve into thermal states. This phenomenon is behind the famous black hole information loss paradox which provoked long-standing debates on the compatibility of quantum mechanics and gravity. In this paper we demonstrate that within an alternative time-symmetric quantization formalism thermal radiation is absent and states evolve unitarily in spacetimes with event horizons. We also discuss the theoretical consistency of the proposed formalism. We explicitly demonstrate that the theory preserves the microcausality condition and suggest a "reinterpretation postulate" to resolve other apparent pathologies associated with negative energy states. Accordingly as there is a consistent alternative, we argue that choosing to use time-asymmetric quantization is a necessary condition for the black hole information loss paradox., Comment: 9 pages

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 p̄ p) 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 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 on 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., United States. Department of Energy (grant contract number DE-SC00012567), United States. Department of Energy (grant contract number DE−SC0013999)

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., 40 pages, 36 figures

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., 18 pages, 15 figures, 2 tables; comments welcome. Updated to published version