Upper Limit on the QCD Axion Mass from Isolated Neutron Star Cooling.

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.

Dark matter from axion strings with adaptive mesh refinement.

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.

Galactic Potential and Dark Matter Density from Angular Stellar Accelerations.

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.

Axion Emission Can Explain a New Hard X-Ray Excess from Nearby Isolated Neutron Stars.

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.

Early-Universe Simulations of the Cosmological Axion.

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.

Probing Dark Matter Using Precision Measurements of Stellar Accelerations.

Dark matter comprises the bulk of the matter in the Universe, but its particle nature and cosmological origin remain mysterious. Knowledge of the dark matter density distribution in the Milky Way Galaxy is crucial both to our understanding of the standard cosmological model and for grounding direct and indirect searches for the particles comprising dark matter. Current measurements of Galactic dark matter content rely on model assumptions to infer the forces acting upon stars from the distribution of observed velocities. Here, we propose to apply the precision radial velocity method, optimized in recent years for exoplanet astronomy, to measure the change in the velocity of stars over time, thereby providing a direct probe of the local gravitational potential in the Galaxy. Using numerical simulations, we develop a realistic strategy to observe the differential accelerations of stars in our Galactic neighborhood with next-generation telescopes, at the level of 10^{-8} cm/s^{2}. Our simulations show that detecting accelerations at this level with an ensemble of 10^{3} stars requires the effect of stellar noise on radial velocity measurements to be reduced to <10 cm/s. The measured stellar accelerations may then be used to extract the local dark matter density and morphological parameters of the density profile.

Stellar Wakes from Dark Matter Subhalos.

We propose a novel method utilizing stellar kinematic data to detect low-mass substructure in the Milky Way's dark matter halo. By probing characteristic wakes that a passing dark matter subhalo leaves in the phase-space distribution of ambient halo stars, we estimate sensitivities down to subhalo masses of ∼10^{7} M_{⊙} or below. The detection of such subhalos would have implications for dark matter and cosmological models that predict modifications to the halo-mass function at low halo masses. We develop an analytic formalism for describing the perturbed stellar phase-space distributions, and we demonstrate through idealized simulations the ability to detect subhalos using the phase-space model and a likelihood framework. Our method complements existing methods for low-mass subhalo searches, such as searches for gaps in stellar streams, in that we can localize the positions and velocities of the subhalos today.