Using Pulsar Parameter Drifts to Detect Subnanohertz Gravitational Waves

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 Letter, 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×10^{-12} at 450 pHz with a sensitivity dropping approximately quadratically with frequency until 10 pHz. Our search method opens a new frequency range for gravitational wave detection and has profound implications for astrophysics, cosmology, and particle physics.