Using Pulsar Parameter Drifts to Detect Sub-Nanohertz Gravitational Waves

Abstract

Gravitational waves with frequencies below 1~nHz are notoriously difficult to detect. With periods exceeding current experimental lifetimes, they induce slow drifts in observables rather than periodic correlations. Observables with well-known intrinsic contributions provide a means to probe this regime. In this work, we demonstrate the viability of using observed pulsar timing parameters to discover such "ultralow" frequency gravitational waves, presenting two complementary observables for which the systematic shift induced by ultralow-frequency gravitational waves can be extracted. Using existing data for these parameters, we search the ultralow frequency regime for continuous-wave signals, finding a sensitivity near the expected prediction from inspirals of supermassive black holes. We do not see an excess in the data, setting a limit on the strain of $ 1.3 \times 10 ^{ - 12} $ at $ 450~{\rm pHz} $ with a sensitivity dropping approximately quadratically with frequency until $ 10~{\rm pHz}$. Our search method opens a new frequency range for gravitational wave detection and has profound implications for astrophysics, cosmology, and particle physics.