Generalized ray tracing for axions in astrophysical plasmas

Ray tracing plays a vital role in black hole imaging, modeling the emission mechanisms of pulsars, and deriving signatures from physics beyond the Standard Model. In this work we focus on one specific application of ray tracing, namely, predicting radio signals generated from the resonant conversion of axion dark matter in the strongly magnetized plasma surrounding neutron stars. The production and propagation of low-energy photons in these environments are sensitive to both the anisotropic response of the background plasma and curved spacetime; here, we employ a fully covariant framework capable of treating both effects. We implement this both via forward and backward ray tracing. In forward ray tracing, photons are sampled at the point of emission and propagated to infinity, whilst in the backward-tracing approach, photons are traced backwards from an image plane to the point of production. We explore various approximations adopted in prior work, quantifying the importance of gravity, plasma anisotropy, the neutron star mass and radius, and imposing the proper kinematic matching of the resonance. Finally, using a more realistic model for the charge distribution of magnetar magnetospheres, we revisit the sensitivity of current and future radio and sub-mm telescopes to spectral lines emanating from the Galactic Center magnetar, showing such observations may extend sensitivity to axion masses m_a ∼ O(few) × 10⁻³  eV, potentially even probing parameter space of the QCD axion.