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.

We revisit a question asked by Dyson: "Is a graviton detectable"We demonstrate that in both Dyson's original sense and in a more modern measurement-theoretic sense, it is possible to construct a detector sensitive to single gravitons, and in fact a variety of existing and near-term gravitational wave detectors can achieve this. However, while such a signal would be consistent with the quantization of the gravitational field, we draw on results from quantum optics to show how the same signal could just as well be explained via classical gravitational waves. We outline the kind of measurements that would be needed to demonstrate quantization of gravitational radiation and explain why these are substantially more difficult than simply counting graviton clicks or observing gravitational noise in an interferometer, and likely impossible to perform in practice.

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 that 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 factorized 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 parameterization 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 parameterization is presented for these priors that 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.

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 g_{aγγ}≲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 g_{aee}≲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.

The dark matter explanation of the 3.5 keV line is strongly disfavored by our work in Dessert et al. (2020). Boyarsky et al. (2004) 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. (2020), though we also showed that the inclusion of such lines is not necessary. (2) Despite suggestions in Boyarsky et al. (2004), even the more conservative limits strongly disfavor a decaying dark matter origin of the 3.5 keV line.

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 instance, we include all relevant electroweak interactions.

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 ^{129}Xe 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.

The Hitomi x-ray satellite mission carried unique high-resolution spectrometers that were set to revolutionize the search for sterile neutrino dark matter (DM) by looking for narrow x-ray lines arising from DM decays. Unfortunately, the satellite was lost shortly after launch, and to date the only analysis using Hitomi for DM decay used data taken towards the Perseus cluster. In this work we present a significantly more sensitive search from an analysis of archival Hitomi data towards blank sky locations, searching for DM decaying in our own Milky Way. The recently launched XRISM satellite has nearly identical soft-x-ray spectral capabilities to Hitomi; we project the full-mission sensitivity of XRISM for analyses of their future blank-sky data, and we find that XRISM will have the leading sensitivity to decaying DM for masses between roughly 1 to 18 keV, with important implications for sterile neutrino and heavy axionlike particle DM scenarios.

Thermal Higgsino dark matter (DM), with a mass near 1.1 TeV, is one of the most well-motivated and untested DM candidates. Leveraging recent hydrodynamic cosmological simulations that give DM density profiles in Milky Way analog galaxies we show that the linelike gamma-ray signal predicted from Higgsino annihilation in the Galactic Center could be detected at high significance with the upcoming Cherenkov Telescope Array (CTA) and Southern Wide-field Gamma-ray Observatory (SWGO) for all but the most pessimistic DM profiles. We perform the most sensitive search to-date for the linelike signal using 15 years of data from the Fermi Large Area Telescope, coming within an order one factor of the necessary sensitivity to detect the Higgsino for some Milky Way analog DM density profiles. We show that H.E.S.S. has subleading sensitivity relative to Fermi for the Higgsino at present. In contrast, we analyze H.E.S.S. inner Galaxy data for the thermal wino model with a mass near 2.8 TeV; we find no evidence for a DM signal and exclude the wino by over a factor of two in cross section for all DM profiles considered. In the process, we identify and attempt to correct what appears to be an inconsistency in previous H.E.S.S. inner Galaxy analyses for DM annihilation related to the analysis effective area, which may weaken the DM cross-section sensitivity claimed in those works by around an order of magnitude.

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.