Measuring scattering variations in pulsar timing observations: a test of the fidelity of current methods

ABSTRACT The turbulent nature of the ionized interstellar medium (IISM) causes dispersion measure (DM) and scattering variations in pulsar timing measurements. To improve precision of gravitational wave measurements, pulsar timing array (PTA) collaborations have begun the use of sophisticated and intricate noise modelling techniques such as modelling stochastic variations induced by the turbulent IISM and quasi-deterministic processes attributed to discrete structures. However, the reliability of these techniques has not been studied in detail, and it is unclear whether the recovered processes are physical or if they are impacted by misspecification. In this work, we present an analysis to test the efficacy of IISM noise models based on the data from the MeerKAT Pulsar Timing Array (MPTA) 4.5-yr data release. We first performed multifrequency, long-length (500 refractive length scale) simulations of multipath propagation in the IISM to study the properties of scattering variations under a variety of scattering conditions. The results of our simulations show the possibility of significant radio-frequency decorrelation in the scattering variations, particularly for the anisotropic scattering medium. Our analysis of the observed DM and scattering variations using the MPTA 4.5-yr data set shows that there can be apparent anticorrelations between DM and scattering variations, which we attribute to the model fitting methods. We also report a possibility that plasma underdensities might exist along the sight lines of PSR J1431–5740 and PSR J1802–2124. Finally, using simulations, we show that the IISM noise models can result in the apparent measurement of strong frequency dependence of scattering variations observed in the MPTA data set. Our analysis shows that improvements in the IISM noise modelling techniques are necessary to accurately measure the IISM properties.