If dark matter consists of ultralight bosons, on galactic scales, it can be effectively described as a coherent classical field experiencing oscillations. Such a field could perturb the dynamics of celestial bodies via a direct coupling to ordinary matter, introducing signatures detectable through high-precision pulsar timing analysis. In this work, we extend a two-step Bayesian inference framework, originally developed for linearly coupled scalar ultralight dark matter (ULDM), to probe a quadratic scalar coupling and spin-1 vector dark matter. By explicitly marginalizing over nuisance orbital parameters, our approach provides robust sensitivity limits that avoid the artificial overestimation often associated with direct fitting techniques. For quadratic scalar ULDM, we establish new constraints on the coupling $\beta$ for masses between $2\times {10}^{-22}$ and $2\times {10}^{-21}\mathrm{eV}$ inaccessible to other experiments, while identifying mass regimes where the sensitivity is dominated by the orbital phase ${\mathrm{\Psi }}^{\prime }$ or the projected semimajor axis $x$. For vector ULDM, we characterize resonant signatures present even in circular orbits and obtain bounds on the coupling $g$ within the ${10}^{-23}$ to ${10}^{-18}\mathrm{eV}$ range, yielding results within the same orders of magnitude as current laboratory and space-based experiments.

Close

PrevNext

View figure in article

• Cosmology
• Dark matter
• Gravitation
• Particle astrophysics
• Particle dark matter
• Axions
• Neutron stars & pulsars