Bayesian Forecasts for Dark Matter Substructure Searches with Mock Pulsar Timing Data

Understand

Dark matter substructure, such as primordial black holes (PBHs) and axion miniclusters, can induce phase shifts in pulsar timing arrays (PTAs) measurements due to gravitational effects.

• In order to gain a more realistic forecast for the detectability of such models of dark matter with PTAs, we propose a Bayesian inference framework to search for phase shifts generated by PBHs and perform the analysis on mock PTA data.
• For most PBH masses the constraints on the dark matter abundance agree with previous (frequentist) analyses (without mock data) to $\mathcal{O}(1)$ factors.
• This further motivates a dedicated search for PBHs (and dense small scale structures) in the mass range from $10^{-8}\,M_{\odot}$ to well above $10^2\,M_{\odot}$ with the Square Kilometer Array.