Search

The stochastic gravitational-wave background is imprinted on the times of arrival of radio pulses from millisecond pulsars. Traditional pulsar timing analyses fit a timing model to each pulsar and search the residuals of the fit for a stationary time correlation. This method breaks down at gravitational-wave frequencies below the inverse observation time of the array; therefore, existing analyses restrict their searches to frequencies above 1 nHz. An effective method to overcome this challenge is to study the correlation of secular drifts of parameters in the pulsar timing model itself. In this paper, we show that timing model correlations are sensitive to sub-nanohertz stochastic gravitational waves and perform a search using existing measurements of binary spin-down rates and pulsar spin-decelerations. We do not observe a signal at our present sensitivity, constraining the stochastic gravitational-wave relic energy density to $\Omega_\text{GW} ( f ) < 3.9 \times 10 ^{ - 9} $ at 450 pHz with sensitivity which scales as the frequency squared until approximately 10 pHz. We place additional limits on the amplitude of a power-law spectrum of $A_\star \lesssim 8\times10^{-15}$ for the spectral index expected from supermassive black hole binaries, $\gamma = 13/3$. If a detection of a supermassive black hole binary signal above 1 nHz is confirmed, this search method will serve as a critical complementary probe of the dynamics of galaxy evolution., Comment: 13 pages, 2 figures

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 supermassive black hole mergers. We do not see an excess in the data, setting a limit on the strain of $ 7.1 \times 10 ^{ - 14} $ at 1 nHz 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., Comment: 12 pages, 2 figures, 3 appendices