Axions and other very weakly interacting slim particles (WISPs), with masses below 1 GeV, arise naturally in many extensions of the Standard Model of particle physics. In particular, they could offer a new framework to explain the nature of dark matter and may help address a range of puzzling observations in astrophysics and particle physics. This review provides an overview of ongoing WISP searches and outlines the prospects for the next decade, spanning their theoretical motivation, indirect signatures in astrophysical observations, and dedicated laboratory experiments. It is based on the work carried on by the EU-funded COST Action “Cosmic WISPers in the Dark Universe: Theory, astrophysics, and experiments” (CA21106, https://www.cost.eu/actions/CA21106). This network plays a key role in coordinating and supporting WISP searches across Europe, while also contributing to the development of a roadmap aimed at securing European leadership in this research area. It is emphasized that Europe is currently pursuing a rich, diverse, and cost-effective experimental program, with the potential to deliver one or more transformative discoveries.

The quantum chromodynamics (QCD) axion arises as the pseudo-Goldstone mode of a spontaneously broken Abelian Peccei-Quinn (PQ) symmetry. If the scale of PQ symmetry breaking occurs below the inflationary reheat temperature and the domain wall number is unity, then there is a unique axion mass that gives the observed dark matter (DM) abundance. Computing this mass has been the subject of intensive numerical simulations for decades since the mass prediction informs laboratory experiments. Axion strings develop below the PQ symmetry-breaking temperature, and as the string network evolves, it emits axions that go on to become the DM. A key ingredient in the axion mass prediction is the spectral index of axion radiation emitted by the axion strings. We compute this index in this Letter using the most precise and accurate large-scale simulations to date of the axion-string network leveraging adaptive mesh refinement to achieve the precision that would, otherwise, require a static lattice with 262,1443 lattice sites. We find a scale-invariant axion radiation spectrum to within 1% precision and find no evidence that the spectral index of radiation evolves with time. Accounting for axion production from strings prior to the QCD phase transition leads us to predict that the axion mass should be approximately ma∈(45,65) μeV. However, we provide preliminary evidence that axions are produced in greater quantities from the string-domain-wall network collapse during the QCD phase transition, potentially increasing the mass prediction to as much as 300 μeV.

Understanding the nonlinear dynamics of coupled scalar fields often necessitates simulations on a 3D mesh. These simulations can be computationally expensive if a large scale separation is involved. A common solution is adaptive mesh refinement, which, however, greatly increases a simulation's complexity. In this work, we present sledgehamr , an AMReX -based code package to make the simulation of coupled scalar fields on an adaptive mesh more accessible. Compatible with both GPU and CPU clusters, sledgehamr offers a flexible and customizable framework. While the code had been primarily developed to evolve axion string networks, this framework enables various other applications, such as the study of gravitational waves sourced by the dynamics of scalar fields.

Axion strings are horizon-size topological defects that may be produced in the early Universe. Ultra-light axion-like particles may form strings that persist to temperatures below that of big bang nucleosynthesis. Such strings have been considered previously as sources of gravitational waves and cosmic microwave background (CMB) polarization rotation. In this work we show, through analytic arguments and dedicated adaptive mesh refinement cosmological simulations, that axion strings deposit a sub-dominant fraction of their energy into high-energy Standard Model (SM) final states, for example, by the direct production of heavy radial modes that subsequently decay to SM particles. This high-energy SM radiation is absorbed by the primordial plasma, leading to novel signatures in precision big bang nucleosynthesis, the CMB power spectrum, and gamma-ray surveys. In particular, we show that CMB power spectrum data constrains axion strings with decay constants $f_a \lesssim 10^{12}$ GeV, up to model dependence on the ultraviolet completion, for axion masses $m_a \lesssim 10^{-29}$ eV; future CMB surveys could find striking evidence of axion strings with lower decay constants.

The quantum chromodynamics (QCD) axion may solve the strong CP problem and explain the dark matter (DM) abundance of our Universe. The axion was originally proposed to arise as the pseudo-Nambu Goldstone boson of global $\mathrm{U}(1)_{\rm PQ}$ Peccei-Quinn (PQ) symmetry breaking, but axions also arise generically in string theory as zero modes of higher-dimensional gauge fields. In this work we show that string theory axions behave fundamentally differently from field theory axions in the early Universe. Field theory axions may form axion strings if the PQ phase transition takes place after inflation. In contrast, we show that string theory axions do not generically form axion strings. In special inflationary paradigms, such as D-brane inflation, string theory axion strings may form; however, their tension is parametrically larger than that of field theory axion strings. We then show that such QCD axion strings overproduce the DM abundance for all allowed QCD axion masses and are thus ruled out, except in scenarios with large warping. A loop-hole to this conclusion arises in the axiverse, where an axion string could be composed of multiple different axion mass eigenstates; a heavier eigenstate could collapse the network earlier, allowing for the QCD axion to produce the correct DM abundance and also generating observable gravitational wave signals.

This is the final technical report of the DOE Early Career award for PI Safdi. The award was active from 09/01/2018 - 05/01/2021. The Early Career award funded theoretical work to investigate the particle nature of dark matter across a range of astrophysical and particle physics scales.

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.

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.

We summarize progress made in theoretical astrophysics and cosmology over the past decade and areas of interest for the coming decade. This Report is prepared as the TF09 "Astrophysics and Cosmology" topical group summary for the Theory Frontier as part of the Snowmass 2021 process.

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.

The non-linear process of cosmic structure formation produces gravitationally bound overdensities of dark matter known as halos. The abundances, density profiles, ellipticities, and spins of these halos can be tied to the underlying fundamental particle physics that governs dark matter at microscopic scales. Thus, macroscopic measurements of dark matter halos offer a unique opportunity to determine the underlying properties of dark matter across the vast landscape of dark matter theories. This white paper summarizes the ongoing rapid development of theoretical and experimental methods, as well as new opportunities, to use dark matter halo measurements as a pillar of dark matter physics.

Exploring dark matter via observations of extreme astrophysical environments — defined here as heavy compact objects such as white dwarfs, neutron stars, and black holes, as well as supernovae and compact object merger events — has been a major field of growth since the last Snowmass process. Theoretical work has highlighted the util-ity of current and near-future observatories to constrain novel dark matter parameter space across the full mass range. This includes gravitational wave instruments and observatories spanning the electromagnetic spectrum, from radio to gamma-rays. While recent searches already provide leading sensitivity to various dark matter models, this work also highlights the need for theoretical astrophysics research to better constrain the properties of these extreme astrophysical systems. The unique potential of these search signatures to probe dark matter adds motivation to proposed next-generation astronomical and gravitational wave instruments.

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.

Axions may be produced thermally inside the cores of neutron stars (NSs), escape the stars due to their feeble interactions with matter, and subsequently convert into x rays in the magnetic fields surrounding the stars. We show that a recently discovered excess of hard x-ray emission in the 2-8 keV energy range from the nearby magnificent seven isolated NSs could be explained by this emission mechanism. These NSs are unique in that they had previously been expected to only produce observable flux in the UV and soft x-ray bands from thermal surface emission at temperatures ∼100  eV. No conventional astrophysical explanation of the magnificent seven hard x-ray excess exists at present. We show that the hard x-ray excess may be consistently explained by an axionlike particle with mass m_{a}≲2×10^{-5}  eV and g_{aγγ}×g_{ann}∈(2×10^{-21},10^{-18})  GeV^{-1} at 95% confidence, accounting for both statistical and theoretical uncertainties, where g_{aγγ} (g_{ann}) is the axion-photon (axion-neutron) coupling constant.

Malte Buschmann, ∗ Benjamin R. Safdi, 3, † and Katelin Schutz ‡ Department of Physics, Princeton University, Princeton, NJ 08544, USA Berkeley Center for Theoretical Physics, University of California, Berkeley, CA 94720, USA Theoretical Physics Group, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA (Dated: March 10, 2021)

We present an approach to measure the Milky Way (MW) potential using the angular accelerations of stars in aggregate as measured by astrometric surveys like Gaia. Accelerations directly probe the gradient of the MW potential, as opposed to indirect methods using, e.g., stellar velocities. We show that end-of-mission Gaia stellar acceleration data may be used to measure the potential of the MW disk at approximately 3σ significance and, if recent measurements of the solar acceleration are included, the local dark matter density at ∼2σ significance. Since the significance of detection scales steeply as t^{5/2} for observing time t, future surveys that include angular accelerations in the astrometric solutions may be combined with Gaia to precisely measure the local dark matter density and shape of the density profile.

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.

Ultracompact dark matter (DM) minihalos at masses at and below 10^{-12}  M_{⊙} arise in axion DM models where the Peccei-Quinn (PQ) symmetry is broken after inflation. The minihalos arise from density perturbations that are generated from the nontrivial axion self-interactions during and shortly after the collapse of the axion-string and domain-wall network. We perform high-resolution simulations of this scenario starting at the epoch before the PQ phase transition and ending at matter-radiation equality. We characterize the spectrum of primordial perturbations that are generated and comment on implications for efforts to detect axion DM. We also measure the DM density at different simulated masses and argue that the correct DM density is obtained for m_{a}=25.2±11.0  μeV.

This report summarises the physics opportunities in the search and study of physics beyond the Standard Model at a 100 TeV pp collider.