⚛️ general relativity

Probing Picohertz Gravitational Waves with Pulsars

This paper presents the first Bayesian search for continuous picohertz gravitational waves using the drifts of pulsar binary orbital and spin periods, achieving a tenfold sensitivity improvement over previous efforts and demonstrating that future Square Kilometre Array observations will be capable of detecting signals from supermassive black hole mergers and probing new physics.

Original authors: Qinyuan Zheng, Chiara M. F. Mingarelli, William DeRocco, Jonathan Nay, Kimberly K. Boddy, Jeff A. Dror

Published 2026-02-09

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Qinyuan Zheng, Chiara M. F. Mingarelli, William DeRocco, Jonathan Nay, Kimberly K. Boddy, Jeff A. Dror

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe is filled with a constant, low hum of ripples in space-time, known as gravitational waves. Scientists have already learned to "hear" these ripples using different tools: giant laser detectors on Earth catch the high-pitched "cracks" from colliding stars, and future space missions will listen for the mid-range "thumps" from massive black holes.

But there is a huge gap in the music: a very deep, slow "bass note" called the picohertz (pHz) regime. These waves are so slow that a single ripple takes hundreds or even thousands of years to pass by. Our current detectors, which have only been listening for a few decades, are like trying to hear a song by listening to just a single second of it. You can't hear the melody; you only hear a tiny, static slice.

The New Approach: Listening to the "Drift"
This paper proposes a clever new way to hear these ultra-slow waves using pulsars. Pulsars are cosmic lighthouses—dead stars that spin incredibly fast and beam radio waves at us with perfect regularity. They are nature's most precise clocks.

Usually, scientists look for "timing residuals," which are tiny glitches in the clock's rhythm caused by a passing wave. But for these super-slow waves, the rhythm doesn't glitch; it just drifts very slowly over time.

The authors suggest we stop looking for the glitch and start measuring the drift itself. They focus on two specific ways the pulsar's "clock" changes:

The Orbital Drift ( $\dot{P}_b$ ): If a pulsar has a companion star, the time it takes to orbit that star changes slightly.
The Spin Drift ( $\ddot{P}$ ): The speed at which the pulsar itself spins changes slightly (specifically, how the change in speed is changing).

Think of it like this: If you are driving a car and a giant, slow-moving wave of air pushes you, you won't feel a sudden bump. Instead, your speedometer will slowly creep up or down. By measuring that slow creep over many years, you can prove the wave is there, even if you never saw the wave itself.

The Detective Work
The researchers used a sophisticated statistical tool (a "Bayesian search") to sift through data from 30 pulsars (more than double the number used in previous attempts). They looked for a specific pattern in the drift of these pulsars that would match the signature of a continuous gravitational wave.

The Results

No Signal Found: Just like a detective searching a crime scene and finding no fingerprints, they did not detect a specific gravitational wave signal in the current data.
Better Limits: However, they didn't come up empty-handed. They established the strictest rules yet for what these waves cannot be. They improved the sensitivity of the search by ten times compared to previous studies. This means if these waves exist, they must be even quieter than we previously thought.

The Future: The Square Kilometre Array (SKA)
The paper is optimistic about the future. They predict that the upcoming Square Kilometre Array (SKA), a massive new radio telescope, will be a game-changer.

More Clocks: The SKA will find hundreds of new pulsars, giving us a much larger "choir" to listen to.
Sharper Ears: It will measure the pulsars' timing with much greater precision.

Their simulations show that with the SKA, we might finally be able to detect the "drift" caused by the very early stages of supermassive black holes merging. This would allow us to study how these giant black holes evolve long before they actually crash into each other.

Why It Matters
Even though they didn't find a wave yet, this paper builds a new "microphone" for the deepest bass notes in the universe. It opens a door to studying:

Supermassive Black Holes: How they pair up and dance before merging.
New Physics: It could reveal clues about the very early universe, cosmic strings, or other exotic phenomena that happened right after the Big Bang.

In short, this work teaches us how to listen to the universe's slowest, deepest songs, and promises that with the next generation of telescopes, we might finally hear the melody.

Technical Summary: Probing Picohertz Gravitational Waves with Pulsars

Problem Statement
Gravitational wave (GW) astronomy currently operates in distinct frequency bands: ground-based interferometers cover Hz to kHz, space-based missions target mHz, and Pulsar Timing Arrays (PTAs) operate in the nHz regime. A significant observational gap exists between $10^{-16}$ Hz and 1 nHz, known as the picohertz (pHz) regime. In this band, continuous GWs (CWs) from supermassive black hole binaries (SMBHBs) are expected to exist as progenitors to the nHz population. However, traditional PTA analysis, which relies on timing residuals, cannot detect these signals because their periods (hundreds to thousands of years) far exceed current observation baselines (decades). Consequently, these low-frequency signals are typically absorbed into the timing model noise. The authors posit that while the waveform itself is not resolved, pHz GWs induce secular drifts in pulsar timing parameters that persist over observation times, offering a potential detection channel.

Methodology
The authors develop a novel Bayesian framework to search for continuous pHz GWs by analyzing the drift of specific pulsar timing parameters rather than timing residuals.

Analytical Model: The method focuses on two parameters: the time derivative of the binary orbital period ( $\dot{P}_b$ ) and the second-order time derivative of the pulsar spin period ( $\ddot{P}$ ). In the pHz regime, a passing GW induces an effective velocity between the pulsar and the Solar System Barycenter (SSB). This manifests as a constant offset in the observed $\dot{P}_b$ (interpreted as an acceleration, $a_{GW}$ ) and a constant offset in $\ddot{P}$ (interpreted as a jerk, $j_{GW}$ ). The authors derive analytic expressions linking these observed derivatives to the GW strain amplitude ( $h_0$ ) and frequency ( $f_{GW}$ ), accounting for antenna beam patterns and the separation between the Earth and pulsar terms.
Bayesian Framework: A likelihood function is constructed assuming uncorrelated uncertainties across $N$ pulsars. The analysis employs Markov Chain Monte Carlo (MCMC) to determine the posterior distribution of signal parameters ( $h_0, f_{GW}$ , and source geometry $\Theta$ ). Uniform priors are used for sky location and frequency, while a logarithmic uniform prior is applied to the strain amplitude to ensure sufficient sampling of small values. Detection significance is quantified using the Bayes Factor (BF) via the Savage-Dickey density ratio method.
Population Modeling: To assess detection prospects, the authors extend an existing SMBHB population synthesis model (calibrated to NANOGrav 15-year data) into the pHz regime. This model incorporates environmental effects such as dynamical friction and stellar hardening, which dominate over GW emission at large separations, predicting a turnover in the GWB spectrum around 60 pHz.
Datasets: The study utilizes an updated dataset of 30 pulsars for the $\dot{P}_b$ analysis (doubling the sample size of previous work) and 13 pulsars for the $\ddot{P}$ analysis. Distance measurements are critical; for $\dot{P}_b$ , only pulsars with distance precision better than the GW wavelength ( $f_{GW} \Delta d_p \leq 1$ ) are used, while all pulsars are retained for $\ddot{P}$ to avoid information loss, despite lower distance precision.

Key Results

Search Outcome: No continuous wave signal was detected in the current data. The Bayes Factors for both $\dot{P}_b$ ( $BF \approx 0.69$ ) and $\ddot{P}$ ( $BF \approx 0.47$ ) searches were consistent with the null hypothesis.
Upper Limits: The authors set new upper limits on the GW strain amplitude ( $h_0$ ). The analysis improves upon previous constraints (DeRocco and Dror [28]) by approximately an order of magnitude. At a reference frequency of 100 pHz, the upper limit is improved by a factor of 9.53.
Scaling Laws: The sensitivity scaling follows theoretical predictions: $h_0 \propto f_{GW}^{-1}$ for $\dot{P}_b$ and $h_0 \propto f_{GW}^{-2}$ for $\ddot{P}$ . At very low frequencies ( $\lesssim 10$ pHz), sensitivity degrades as the Earth and pulsar terms become "in phase."
Population Constraints: The extended SMBHB model predicts a spectral turnover near 60 pHz due to environmental hardening, resulting in a suppressed GWB power at the lowest frequencies compared to a simple power-law extrapolation.
Future Prospects (SKA): Simulations for the Square Kilometre Array (SKA) indicate that increasing the number of pulsars to $\sim 800$ will improve sensitivity by another order of magnitude solely due to sample size. However, detecting individual SMBHB sources requires a significant reduction in timing parameter uncertainty. The authors find that detecting a realistic CW source would require either a uniform $\sim 2$ -order-of-magnitude improvement in $\ddot{P}$ measurement precision across all SKA pulsars, or a $\sim 3$ -order-of-magnitude improvement in just 10 of the 795 pulsars.

Significance and Claims
The paper claims to present the first Bayesian search for continuous pHz GWs using pulsar timing parameter drifts. Its primary contribution is establishing a robust method that achieves an order-of-magnitude improvement in strain limits over previous studies in this frequency band. The authors argue that this technique opens a new observational frontier for:

Black Hole Astrophysics: Probing the early stages of SMBHB evolution and the environmental mechanisms (dynamical friction, stellar hardening) that drive them before they enter the nHz PTA band.
Fundamental Physics: Testing beyond-the-Standard-Model scenarios, including cosmic strings, first-order phase transitions, and primordial GWs from inflation.

The authors emphasize that while current data yields only upper limits, the methodology is poised to transform into a genuine discovery tool with the advent of SKA, provided that timing precision for higher-order derivatives ( $\ddot{P}$ ) can be sufficiently enhanced. They note that the framework is adaptable for targeted searches using multi-messenger priors (e.g., from dual AGN observations) and for searching for other exotic physics, such as ultra-light dark matter.

More like this