Fluctuations in ultralight dark matter produce significant metric perturbations, which may be detected by monitoring the arrival times of light from millisecond pulsars. While searches using this technique are already underway, they do not consistently account for the statistical properties of the dark matter field. The statistics of this field depend on the velocity dispersion of dark matter and, consequently, its coherence length. In the mass range relevant for pulsar timing arrays, the coherence length is comparable to separations between pulsars, making it crucial to incorporate its effects into the analysis. This Letter presents a consistent statistical method for gravitational direct detection of ultralight dark matter. Our key result is the derivation of the two-point function of the metric fluctuations, which we apply to pulsar timing and discuss its implementation in future searches.
Boddy et al. (Tue,) studied this question.