Axions with couplings g_{aγγ}∼few×10^{-11}  GeV^{-1} to electromagnetism may resolve a number of astrophysical anomalies, such as unexpected ∼TeV transparency, anomalous stellar cooling, and x-ray excesses from nearby neutron stars. We show, however, that such axions are severely constrained by the nonobservation of x rays from the magnetic white dwarf (MWD) RE J0317-853 using ∼40  ks of data acquired from a dedicated observation with the Chandra X-ray Observatory. Axions may be produced in the core of the MWD through electron bremsstrahlung and then convert to x rays in the magnetosphere. The nonobservation of x rays constrains the axion-photon coupling to g_{aγγ}≲5.5×10^{-13}sqrt[C_{aγγ}/C_{aee}]  GeV^{-1} at 95% confidence for axion masses m_{a}≲5×10^{-6}  eV, with C_{aee} and C_{aγγ} the dimensionless coupling constants to electrons and photons. Considering that C_{aee} is generated from the renormalization group, our results robustly disfavor g_{aγγ}≳4.4×10^{-11}  GeV^{-1} even for models with no ultraviolet contribution to C_{aee}.

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.

We outline the unique opportunities and challenges in the search for “ultraheavy” dark matter candidates with masses between roughly 10 TeV and the Planck scale mpl ≈ 1016 TeV. This mass range presents a wide and relatively unexplored dark matter parameter space, with a rich space of possible models and cosmic histories. We emphasize that both current detectors and new, targeted search techniques, via both direct and indirect detection, are poised to contribute to searches for ultraheavy particle dark matter in the coming decade. We highlight the need for new developments in this space, including new analyses of current and imminent direct and indirect experiments targeting ultraheavy dark matter and development of new, ultra-sensitive detector technologies like next-generation liquid noble detectors, neutrino experiments, and specialized quantum sensing techniques.