Abstract : The quantum chromodynamics axion with a decay constant near the Grand Unification (GUT) scale has an ultralight mass near a neV. We show, however, that axion-like particles with masses near the keV–PeV range with GUT-scale decay constants are also well motivated in that they naturally arise from axiverse theories with dark c gauge groups. We demonstrate that the correct dark matter abundance may be achieved by the heavy axions in these models through the misalignment mechanism in combination with a period of early matter domination from the long-lived dark glueballs of the same gauge group. Heavy axion dark matter may decay to two photons, yielding mono-energetic electromagnetic signatures that may be detectable by current or next-generation space-based telescopes. We project the sensitivity of next-generation telescopes including Athena, AMEGO, and e-ASTROGAM to such decaying axion dark matter. If the dark sector contains multiple confining gauge groups, then the observed primordial baryon asymmetry may also be achieved in this scenario through spontaneous baryogenesis. We present explicit orbifold constructions where the dark gauge groups unify with the SM at the GUT scale and axions emerge as the fifth components of dark gauge fields with bulk Chern-Simons terms.

Axions may be produced in abundance inside stellar cores and then convert into observable x rays in the Galactic magnetic fields. We focus on the Quintuplet and Westerlund 1 super star clusters, which host large numbers of hot, young stars including Wolf-Rayet stars; these stars produce axions efficiently through the axion-photon coupling. We use Galactic magnetic field models to calculate the expected x-ray flux locally from axions emitted from these clusters. We then combine the axion model predictions with archival Nuclear Spectroscopic Telescope Array (NuSTAR) data from 10-80 keV to search for evidence of axions. We find no significant evidence for axions and constrain the axion-photon coupling g_{aγγ}≲3.6×10^{-12}  GeV^{-1} for masses m_{a}≲5×10^{-11}  eV at 95% confidence.

Abstract: The 3.5 keV line is a purported emission line observed in galaxies, galaxy clusters, and the Milky Way whose origin is inconsistent with known atomic transitions and has previously been suggested to arise from dark matter decay. We systematically reexamine the bulk of the evidence for the 3.5 keV line, attempting to reproduce six previous analyses that found evidence for the line. We only reproduce one of the analyses; in the other five, we find no significant evidence for a 3.5 keV line when following the described analysis procedures on the original data sets. For example, previous results claimed 4σ evidence for a 3.5 keV line from the Perseus cluster; we dispute this claim, finding no evidence for a 3.5 keV line. We find evidence for background mismodeling in multiple analyses. We show that analyzing these data in narrower energy windows diminishes the effects of mismodeling but returns no evidence for a 3.5 keV line. We conclude that there is little robust evidence for the existence of the 3.5 keV line. Some of the discrepancy of our results from those of the original works may be due to the earlier reliance on local optimizers, which we demonstrate can lead to incorrect results. For ease of reproducibility, all code and data are publicly available.

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 wavelike 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 wavelike 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 subdegree accuracy, or the direction of a putative cold DM stream to the subarcminute 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 three-dimensional DM velocity distribution.

Dark matter (DM) in the Milky Way halo may annihilate or decay to photons, producing monochromatic gamma rays. We search for DM-induced spectral lines using 14 years of data from the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope (Fermi-LAT) between 10 GeV and 2 TeV in the inner Milky Way leveraging both the spatial and spectral morphology of an expected signal. We present new constraints as strong as (σv) ≲6×10-30 cm3/s for the two-to-two annihilations and τ≳1030 s for one-to-two decays, representing leading sensitivity between 10 GeV and ∼500 GeV. We consider the implications of our line-constraints on the Galactic Center excess (GCE), which is a previously-observed excess of continuum ∼GeV gamma-rays that may be explained by DM annihilation. The Higgs portal and neutralinolike DM scenarios, which have been extensively discussed as possible origins of the GCE, are constrained by our work because of the lack of observed one-loop decays to two photons. More generally, we interpret our null results in a variety of annihilating and decaying DM models, such as neutralinos, gravitinos, and glueballs, showing that in many cases the line search is more powerful than the continuum, despite the continuum annihilation being at tree level.

Thermal Higgsino dark matter (DM), with mass around 1 TeV, is a well-motivated, minimal DM scenario that arises in supersymmetric extensions of the standard model. Higgsinos may naturally be the lightest superpartners in split-supersymmetry models that decouple the scalar superpartners while keeping Higgsinos and gauginos close to the TeV scale. Higgsino DM may annihilate today to give continuum γ-ray emission at energies less than a TeV in addition to a linelike signature at energies equal to the mass. Previous searches for Higgsino DM, for example with the H.E.S.S. γ-ray telescope, have not reached the necessary sensitivity to probe the Higgsino annihilation cross section. In this work we make use of 14 years of data from the Fermi Large Area Telescope at energies above ∼10  GeV to search for the continuum emission near the Galactic Center from Higgsino annihilation. We interpret our results using DM profiles from Milky Way analog galaxies in the FIRE-2 hydrodynamic cosmological simulations. We set the strongest constraints to date on Higgsino-like DM. Our results show a mild, ∼2σ preference for Higgsino DM with a mass near the thermal Higgsino mass and, depending on the DM density profile, the expected cross section.

The quantum chromodynamics (QCD) axion may modify the cooling rates of neutron stars (NSs). The axions are produced within the NS cores from nucleon bremsstrahlung and, when the nucleons are in superfluid states, Cooper pair breaking and formation processes. We show that four of the nearby isolated magnificent seven NSs along with PSR J0659 are prime candidates for axion cooling studies because they are coeval, with ages of a few hundred thousand years known from kinematic considerations, and they have well-measured surface luminosities. We compare these data to dedicated NS cooling simulations incorporating axions, profiling over uncertainties related to the equation of state, NS masses, surface compositions, and superfluidity. Our calculations of the axion and neutrino emissivities include high-density suppression factors that also affect SN 1987A and previous NS cooling limits on axions. We find no evidence for axions in the isolated NS data, and within the context of the Kim-Shifman-Vainshtein-Zakharov QCD axion model, we constrain m_{a}≲16  meV at 95% confidence level. An improved understanding of NS cooling and nucleon superfluidity could further improve these limits or lead to the discovery of the axion at weaker couplings.

A number of proposals have been put forward for detecting axion dark matter (DM) with grand unification scale decay constants that rely on the conversion of coherent DM axions to oscillating magnetic fields in the presence of static, laboratory magnetic fields. Crucially, such experiments - including ABRACADABRA - have to date worked in the limit that the axion Compton wavelength is larger than the size of the experiment, which allows one to take a magnetoquasistatic (MQS) approach to characterize the detector apparatus and model the axion signal. We use finite element methods to solve the coupled axion-electromagnetism equations of motion without assuming the MQS approximation. We show that the MQS approximation becomes a poor approximation at frequencies 2 orders of magnitude lower than the naive MQS limit frequency commonly defined by the inverse diameter of a lumped-element detector. Radiation losses diminish the quality factor of an otherwise high-Q resonant readout circuit, though this may be mitigated through shielding and minimizing lossy materials. Additionally, self-resonances associated with the detector geometry change the reactive properties of the pickup system, leading to two generic features beyond MQS: There are frequencies that require an inductive rather than capacitive tuning to maintain resonance, and the detector itself becomes a multipole resonator at high frequencies. Accounting for these features, competitive sensitivity to the axion-photon coupling may be extended well beyond the naive MQS limit.

Axions are hypothetical particles that may explain the observed dark matter density and the non-observation of a neutron electric dipole moment. An increasing number of axion laboratory searches are underway worldwide, but these efforts are made difficult by the fact that the axion mass is largely unconstrained. If the axion is generated after inflation there is a unique mass that gives rise to the observed dark matter abundance; due to nonlinearities and topological defects known as strings, computing this mass accurately has been a challenge for four decades. Recent works, making use of large static lattice simulations, have led to largely disparate predictions for the axion mass, spanning the range from 25 microelectronvolts to over 500 microelectronvolts. In this work we show that adaptive mesh refinement simulations are better suited for axion cosmology than the previously-used static lattice simulations because only the string cores require high spatial resolution. Using dedicated adaptive mesh refinement simulations we obtain an over three order of magnitude leap in dynamic range and provide evidence that axion strings radiate their energy with a scale-invariant spectrum, to within ~5% precision, leading to a mass prediction in the range (40,180) microelectronvolts.

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 Letter, 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.