arxiv: 2306.16215 · v1 · submitted 2023-06-28 · 🌌 astro-ph.HE · astro-ph.CO· astro-ph.GA· gr-qc

Recognition: 2 theorem links

· Lean Theorem

Search for an isotropic gravitational-wave background with the Parkes Pulsar Timing Array

Daniel J. Reardon , Andrew Zic , Ryan M. Shannon , George B. Hobbs , Matthew Bailes , Valentina Di Marco , Agastya Kapur , Axl F. Rogers

show 21 more authors

Eric Thrane Jacob Askew N. D. Ramesh Bhat Andrew Cameron Ma{\l}gorzata Cury{\l}o William A. Coles Shi Dai Boris Goncharov Matthew Kerr Atharva Kulkarni Yuri Levin Marcus E. Lower Richard N. Manchester Rami Mandow Matthew T. Miles Rowina S. Nathan Stefan Os{\l}owski Christopher J. Russell Ren\'ee Spiewak Songbo Zhang Xing-Jiang Zhu

Authors on Pith no claims yet

Pith reviewed 2026-05-12 17:01 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.COastro-ph.GAgr-qc

keywords gravitational wave backgroundpulsar timing arraysupermassive black hole binariesnanohertz gravitational wavesParkes Pulsar Timing Arraystochastic gravitational wavesspatial correlations

0 comments

The pith

Parkes Pulsar Timing Array data recovers a common-spectrum process and spatial correlations consistent with an isotropic gravitational-wave background at about 2 sigma.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper reports results from an 18-year search for nanohertz gravitational waves using timing observations of 30 millisecond pulsars. A Bayesian analysis identifies a common red noise process across the pulsars with an amplitude and spectral index that can be interpreted as arising from inspiraling supermassive black hole binaries. When the spectral index is fixed to the value expected for such binaries, the amplitude is measured as 2.04^{+0.25}_{-0.22} × 10^{-15}. Hierarchical analysis of pulsar pairs reveals spatial correlations matching the predicted pattern for an isotropic background, with a false-alarm probability of at most 0.02. However, the strength of the signal depends on the data segment used, creating tension with upper limits from the first half of the observations.

Core claim

We measure a common-spectrum noise process represented as a strain spectrum with amplitude A = 3.1^{+1.3}_{-0.9} × 10^{-15} and spectral index α = -0.45 ± 0.20. For the expected index of α = -2/3 we recover A = 2.04^{+0.25}_{-0.22} × 10^{-15}. Through analysis of individual pulsar pairs we measure spatial correlations consistent with a gravitational-wave background, with an estimated false-alarm probability of p ≲ 0.02.

What carries the argument

The common-spectrum process modeled as a power-law gravitational-wave strain spectrum combined with hierarchical Bayesian analysis of pairwise spatial correlations in pulsar timing residuals, validated by randomizing pulsar sky positions.

If this is right

The measured amplitude implies a population of supermassive black hole binaries radiating gravitational waves at nanohertz frequencies.
The long 18-year baselines and access to southern pulsars improve constraints on the isotropic background compared to northern-only arrays.
Continued observations with the Parkes array will help clarify the time dependence of the signal and contribute to the International Pulsar Timing Array.
The 2-sigma spatial correlation supports the interpretation that the common noise is due to gravitational waves rather than other processes.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the signal is real, it provides the first direct evidence for the existence of a stochastic gravitational wave background at nanohertz frequencies.
Cross-validation with data from other pulsar timing arrays would be needed to confirm the detection.
The observed tension between data halves suggests that refined noise modeling may be required before claiming a definitive detection.
Future increases in timing precision could allow measurement of the background's spectrum in more detail.

Load-bearing premise

The common-spectrum process and the observed spatial correlations are due to an isotropic gravitational-wave background and not to unmodeled time-varying noise, systematics, or selection effects in the pulsar timing residuals.

What would settle it

If a reanalysis of the full dataset using more comprehensive noise models eliminates the common-spectrum process or the spatial correlations while fitting the data equally well, the gravitational-wave background interpretation would be falsified.

read the original abstract

Pulsar timing arrays aim to detect nanohertz-frequency gravitational waves (GWs). A background of GWs modulates pulsar arrival times and manifests as a stochastic process, common to all pulsars, with a signature spatial correlation. Here we describe a search for an isotropic stochastic gravitational-wave background (GWB) using observations of 30 millisecond pulsars from the third data release of the Parkes Pulsar Timing Array (PPTA), which spans 18 years. Using current Bayesian inference techniques we recover and characterize a common-spectrum noise process. Represented as a strain spectrum $h_c = A(f/1 {\rm yr}^{-1})^{\alpha}$, we measure $A=3.1^{+1.3}_{-0.9} \times 10^{-15}$ and $\alpha=-0.45 \pm 0.20$ respectively (median and 68% credible interval). For a spectral index of $\alpha=-2/3$, corresponding to an isotropic background of GWs radiated by inspiraling supermassive black hole binaries, we recover an amplitude of $A=2.04^{+0.25}_{-0.22} \times 10^{-15}$. However, we demonstrate that the apparent signal strength is time-dependent, as the first half of our data set can be used to place an upper limit on $A$ that is in tension with the inferred common-spectrum amplitude using the complete data set. We search for spatial correlations in the observations by hierarchically analyzing individual pulsar pairs, which also allows for significance validation through randomizing pulsar positions on the sky. For a process with $\alpha=-2/3$, we measure spatial correlations consistent with a GWB, with an estimated false-alarm probability of $p \lesssim 0.02$ (approx. $2\sigma$). The long timing baselines of the PPTA and the access to southern pulsars will continue to play an important role in the International Pulsar Timing Array.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

PPTA DR3 recovers a common-spectrum process at A ~ 2e-15 for alpha=-2/3 but shows clear time dependence that undercuts a stationary GWB claim, with spatial correlations only at ~2 sigma via pair randomization.

read the letter

The main things to know are that the third Parkes data release finds a common red-noise process whose amplitude for a -2/3 spectrum sits around 2.04 x 10^-15, and that a hierarchical pulsar-pair analysis with position randomization gives spatial correlations at roughly 2 sigma. The paper also flags that the first half of the 18-year span yields an upper limit in tension with the full-dataset amplitude. That tension is the central caveat they report openly in the abstract. They apply standard Bayesian inference to the timing residuals of 30 millisecond pulsars and add the pair-wise approach to test the Hellings-Downs pattern without fitting the entire array at once. The long baseline and southern-sky coverage are real assets for the broader IPTA effort, and the randomization test supplies a straightforward way to estimate false-alarm probability. The time dependence is the clearest soft spot. A nanohertz background from inspiraling supermassive black holes should be stationary across these timescales. When excess power appears concentrated in the later half, it opens the door to unmodeled red noise, calibration effects, or selection biases that could produce both the common spectrum and the apparent correlations. The abstract notes the tension but leaves open how the spatial statistic behaves once that non-stationarity is modeled. The work is aimed at the PTA and low-frequency gravitational-wave community. Anyone following the accumulating IPTA results will want the full methods on noise modeling and the individual pair contributions. It is incremental rather than decisive, but the dataset and the explicit caveat make it worth referee time. I would send it to peer review. The data volume is substantial and the analysis choices, especially around time variation, need external scrutiny before the spatial-correlation claim can be taken as firm evidence.

Referee Report

2 major / 2 minor

Summary. The manuscript reports a Bayesian search for an isotropic stochastic gravitational-wave background using 18 years of PPTA DR3 timing data on 30 millisecond pulsars. It recovers a common-spectrum process with amplitude A=3.1^{+1.3}_{-0.9}×10^{-15} and spectral index α=-0.45±0.20; fixing α=-2/3 yields A=2.04^{+0.25}_{-0.22}×10^{-15}. A hierarchical pair-wise analysis finds spatial correlations consistent with the Hellings-Downs pattern at FAP p≲0.02 (≈2σ) via pulsar-position randomization. The paper explicitly notes that the common-spectrum amplitude is time-dependent, with the first half of the dataset yielding an upper limit in tension with the full-dataset value.

Significance. If the spatial correlations are shown to arise from a stationary isotropic GWB rather than time-varying noise or systematics, the result would constitute an important incremental step toward nanohertz GW detection and would support the expected SMBHB background. The long baseline and southern-sky coverage add value to IPTA efforts. The explicit demonstration of time-dependence, however, limits the strength of the current claim.

major comments (2)

[Abstract and results on time-dependence] Abstract and the results section on time-dependence: the reported tension between the full-dataset amplitude A=2.04^{+0.25}_{-0.22}×10^{-15} (α=-2/3) and the upper limit obtained from only the first half of the 18-year span directly challenges the stationarity assumption required to interpret the pair-wise spatial correlations as evidence for a GWB. A true SMBHB background must be stationary; if the excess common power is concentrated in the later half, the position-randomization FAP (p≲0.02) may be contaminated by non-stationary red noise, calibration effects, or selection biases.
[Spatial-correlation analysis section] Section describing the hierarchical pair-wise spatial-correlation analysis: the randomization test for false-alarm probability must be shown to remain valid when the common-spectrum process exhibits the reported time-dependence. If the excess power is localized in time, randomizing only sky positions (without also resampling or splitting the time series) does not fully test whether the observed correlations could arise from time-localized systematics rather than the isotropic Hellings-Downs pattern.

minor comments (2)

[Methods] Clarify in the methods how the common-spectrum process is parameterized relative to the timing-residual power spectral density; the conversion from the reported strain amplitude A to the residual spectrum should be stated explicitly.
[Figures] Figure showing the posterior for A and α (or the time-split comparison) should include the first-half upper limit as an explicit overlay or separate panel for direct visual comparison.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their insightful comments on our manuscript. We have addressed each of the major comments point by point below, making revisions to the manuscript where necessary to improve the discussion of limitations associated with the time-dependent common-spectrum process.

read point-by-point responses

Referee: Abstract and results on time-dependence: the reported tension between the full-dataset amplitude A=2.04^{+0.25}_{-0.22}×10^{-15} (α=-2/3) and the upper limit obtained from only the first half of the 18-year span directly challenges the stationarity assumption required to interpret the pair-wise spatial correlations as evidence for a GWB. A true SMBHB background must be stationary; if the excess common power is concentrated in the later half, the position-randomization FAP (p≲0.02) may be contaminated by non-stationary red noise, calibration effects, or selection biases.

Authors: We agree that the time-dependence of the common-spectrum process challenges the stationarity assumption and thus requires careful interpretation of the spatial correlations as evidence for a GWB. In the revised manuscript, we have modified the abstract and the results section to more clearly highlight this tension and its potential implications for non-stationary effects. We now stress that the FAP from the position randomization should be viewed with this caveat in mind, and that additional data will be needed to confirm the stationarity of the process. revision: yes
Referee: Section describing the hierarchical pair-wise spatial-correlation analysis: the randomization test for false-alarm probability must be shown to remain valid when the common-spectrum process exhibits the reported time-dependence. If the excess power is localized in time, randomizing only sky positions (without also resampling or splitting the time series) does not fully test whether the observed correlations could arise from time-localized systematics rather than the isotropic Hellings-Downs pattern.

Authors: The position randomization test is designed to assess the significance of the spatial correlations by destroying the Hellings-Downs pattern while preserving the individual pulsar timing properties, including any time-dependent features in the data. We maintain that this provides a valid assessment of the false-alarm probability under the observed data properties. Nevertheless, we recognize the referee's point that time-localized systematics could potentially influence the result. We have added text to the manuscript discussing this assumption and the limitations of the current test, recommending that future analyses incorporate time-domain splits or resampling to further validate the findings. revision: partial

Circularity Check

0 steps flagged

No significant circularity; results are direct data-driven measurements

full rationale

The paper performs Bayesian inference on pulsar timing residuals to recover parameters of a common-spectrum process (A and alpha) and then conducts a separate hierarchical pairwise analysis to test for spatial correlations with position randomization for false-alarm probability. These steps are explicit fits and statistical tests applied to the observations, with no reduction of outputs to inputs by construction, no self-definitional equations, and no load-bearing self-citations or ansatzes invoked in the provided text. The noted time-dependence of the signal is an empirical observation about the dataset splits rather than a circular step in the analysis chain. The work is self-contained as an observational search against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The analysis rests on standard PTA assumptions about noise processes and gravitational-wave signal morphology drawn from prior literature; no new entities are introduced.

free parameters (2)

A = 3.1^{+1.3}_{-0.9} × 10^{-15}
Amplitude of the strain spectrum h_c(f) fitted via Bayesian inference to the common-spectrum process in the timing residuals.
alpha = -0.45 ± 0.20
Spectral index of the power-law strain spectrum also fitted from the data.

axioms (1)

domain assumption An isotropic stochastic gravitational-wave background produces a common-spectrum process with Hellings-Downs spatial correlations across pulsar pairs.
Invoked when interpreting the recovered common noise as a GWB and when performing the hierarchical pair analysis.

pith-pipeline@v0.9.0 · 5815 in / 1708 out tokens · 70485 ms · 2026-05-12T17:01:08.473729+00:00 · methodology

discussion (0)

Forward citations

Cited by 29 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Gravity Echoes from Supermassive Black Hole Binaries
astro-ph.HE 2026-04 unverdicted novelty 8.0

Future microhertz detections combined with nanohertz pulsar terms can serve as gravity echoes to measure supermassive black hole binary inspiral rates from hundreds to thousands of years in the past.
Gravitational Memory from Hairy Binary Black Hole Mergers
gr-qc 2026-04 unverdicted novelty 8.0

Gravitational memory from hairy binary black hole mergers in scalar-Gauss-Bonnet gravity differs from GR by a few percent due to altered nonlinear dynamics, with direct scalar contributions suppressed, and including m...
Phase-resolved field-space distance bounds in ekpyrotic, bouncing and cyclic cosmologies
gr-qc 2026-05 unverdicted novelty 7.0

Phase-resolved scalar distance bounds are derived for ekpyrotic, bouncing, and cyclic cosmologies, yielding a master condition that lower-bounds ε_ek from remaining distance after conversion and bounce.
Phase-resolved field-space distance bounds in ekpyrotic, bouncing and cyclic cosmologies
gr-qc 2026-05 unverdicted novelty 7.0

Phase-resolved field-space distance bounds for non-inflationary smoothing yield a master lower bound on ε_ek and imply ultra-fast-roll ekpyrosis or modified bounces to match observed red-tilted perturbations.
Gravitational wave emission from nonspherical collapse in an early matter-dominated era using N-body simulations
astro-ph.CO 2026-05 unverdicted novelty 7.0

Full numerical N-body treatment is required for reliable gravitational wave predictions from nonspherical collapse in early matter-dominated eras, with resulting spectra mappable to detector sensitivities via horizon ...
Self-acceleration of Hardening Binaries
astro-ph.GA 2026-05 unverdicted novelty 7.0

Hardening binaries experience deterministic self-acceleration of their center of mass, induced precession, and plane rotation in uniform isotropic media, driving outward spiraling and eccentricity growth in all cases ...
New gravitational-wave templates for metastable cosmic strings: Loop breaking versus network collapse
hep-ph 2026-04 unverdicted novelty 7.0

Metastable cosmic strings produce a gravitational wave background that is best modeled with three parameters (string tension Gμ plus independent time scales t_LB and t_NC), yielding a compact analytical spectrum when ...
Forecasting graviton-mass constraints from the full covariance of PTA-astrometry ORF estimators
gr-qc 2026-04 unverdicted novelty 7.0

A full-covariance formalism for PTA-astrometry ORF estimators forecasts graviton-mass upper limits of 4.41e-24 eV/c2 for current-like setups and 0.48e-24 eV/c2 for SKA/Theia-like future setups, with astrometry adding ...
Probing Supermassive Black Hole Mergers with Pulsar Timing Arrays
astro-ph.HE 2026-04 unverdicted novelty 7.0

Pulsar timing arrays can probe supermassive black hole binaries that merged prior to observations via the pulsar term, with SKA potentially detecting a few such zombie binaries at SNR > 3.
A practical theorem on gravitational-wave background statistics
astro-ph.CO 2026-04 unverdicted novelty 7.0

For large but finite source counts, the PDF of rescaled GWB characteristic strain squared follows the universal form N^{1/3} times the reflected map-Airy distribution evaluated at N^{1/3}(y-1), fully determined by the...
Gravitational-wave lensing beyond rays: a disordered-system approach
astro-ph.CO 2026-04 unverdicted novelty 7.0

A quenched-disorder approach with Schwinger-Keldysh path integrals produces an averaged density matrix for gravitational waves that separates phase-suppressing exponential terms from oscillatory corrections to coheren...
Detecting Chiral Gravitational Wave Background with a Dipole Pulsar Timing Array
gr-qc 2026-04 unverdicted novelty 7.0

A dipole pulsar timing array detects chiral nanohertz gravitational waves and extends PTA sensitivity into the microhertz regime.
Primordial Black Hole from Tensor-induced Density Fluctuation: First-order Phase Transitions and Domain Walls
astro-ph.CO 2026-05 unverdicted novelty 6.0

Tensor perturbations from first-order phase transitions and domain wall annihilation induce curvature fluctuations at second order that form primordial black holes, allowing asteroid-mass PBHs to comprise all dark mat...
Testing General Relativity with Individual Supermassive Black Hole Binaries
gr-qc 2026-05 unverdicted novelty 6.0

A framework is developed to test beyond-GR effects in nanohertz continuous waves from individual SMBHBs, deriving modified inter-pulsar correlations, antenna responses, and phase delays for three deviation classes, va...
Are PTA measurements sensitive to gravitational wave non-Gaussianities?
astro-ph.CO 2026-05 unverdicted novelty 6.0

PTA statistical tests lose sensitivity to non-Gaussian GW features after decorrelation and cannot distinguish them model-agnostically.
Exploring the statistical anisotropy of primordial curvature perturbations with pulsar timing arrays
gr-qc 2026-04 unverdicted novelty 6.0

A phenomenological dipole anisotropy in primordial perturbations induces dipolar and quadrupolar anisotropies in SIGW energy density spectra, producing frequency-dependent PTA overlap reduction functions that depend o...
Neutron Portal and Dark Matter-Baryon Coincidence: from UV Completion to Phenomenology
hep-ph 2026-04 unverdicted novelty 6.0

A UV-complete neutron portal model dynamically solves the dark matter-baryon coincidence via a supercooled dark confinement transition that generates GeV-scale asymmetric DM and links to observed gravitational waves.
The Heavy Tailed Non-Gaussianity of the Supermassive Black Hole Gravitational Wave Background
astro-ph.CO 2026-04 unverdicted novelty 6.0

The gravitational wave background from supermassive black hole binaries has a universal heavy-tailed amplitude distribution with power-law index -4, causing divergent higher moments and dominance of the strongest sign...
Fixing the Renormalization of Inflationary Loops via Ward Identities
gr-qc 2026-05 unverdicted novelty 5.0

Ward identities from large gauge symmetry impose model-independent constraints on renormalizing inflationary loops and non-perturbatively govern the infrared power spectrum evolution.
One Merge to Rule Them All: From Galaxy Interactions to Black Hole Mergers Using Horizon-AGN
astro-ph.GA 2026-05 unverdicted novelty 5.0

Horizon-AGN shows galaxy and black hole merger rates both rise with stellar mass and fall with redshift, peaking near z=2-3, establishing a direct evolutionary link from galaxy interactions to black hole coalescences.
Reviving Motivated Inflationary Potentials with $K$-inflation in the light of ACT
gr-qc 2026-05 unverdicted novelty 5.0

K-inflation with non-canonical kinetic term G(φ) shifts α-attractor T-models and natural inflation into the Planck-ACT-LB-BK18 allowed region while satisfying Swampland conjectures and producing testable GW spectra.
Transient Parity Violation during Inflation: Implications for PTA Gravitational Waves
astro-ph.CO 2026-04 unverdicted novelty 5.0

A transient parity-violating phase during inflation generates a robust blue-tilted (n_T ≃ 2) primordial gravitational wave spectrum at small scales with nearly maximal helicity coherence and linear polarization, offer...
Primordial Magnetogenesis and Gravitational Waves from ALP-assisted Phase Transition
hep-ph 2026-04 unverdicted novelty 5.0

ALP-assisted first-order phase transitions can explain observed intergalactic magnetic fields and produce detectable gravitational waves, linking cosmology with particle physics searches.
Purely Quadratic Non-Gaussianity from Tachyonic Instability: Primordial Black Holes and Scalar-Induced Gravitational Waves
astro-ph.CO 2026-04 unverdicted novelty 5.0

Purely quadratic non-Gaussianity from tachyonic instability allows narrow curvature spectra to exponentially suppress primordial black hole overproduction via correlation coefficient ρ approaching -1 while retaining s...
Hawking area law in quantum gravity
gr-qc 2026-04 unverdicted novelty 5.0

Exact Hawking area law from black hole mergers restricts quantum gravity to singular Ricci-flat or specific regular black holes in Stelle and nonlocal theories, derives the standard entropy-area law, and realizes Barr...
Gravitational Waves from Matter Perturbations of Spectator Scalar Fields
hep-ph 2026-04 unverdicted novelty 5.0

A spectator scalar field with strong portal coupling to the inflaton sources a stochastic gravitational wave background reaching Ω_GW h² ∼ 10^{-11} at frequencies 10^7-10^8 Hz for benchmark parameters σ/λ ≃ 10^4 and T...
Constraints on Ultralight Scalar and Dark Photon Dark Matter from PPTA-DR3 and EPTA-DR2
astro-ph.CO 2026-05 unverdicted novelty 4.0

Bayesian analysis of PPTA-DR3 and EPTA-DR2 finds no statistically significant ULDM signals and sets 95% CL upper limits on scalar and dark photon dark matter, improving prior bounds in most mass ranges.
Precision Analysis for $\boldsymbol{H_0}$ Using Upcoming Multi-band Gravitational Wave Observations
astro-ph.CO 2026-04 unverdicted novelty 4.0

Multi-band GW observations of PBHs can reduce H0 uncertainty to ≲2 km/s/Mpc (conservative) or O(0.1) km/s/Mpc (optimistic) via Fisher forecasts on M_PBH and f_PBH.
Probing non-Gaussianity during reheating with SIGW in the LISA band
astro-ph.CO 2026-04 unverdicted novelty 4.0

Non-standard reheating imprints detectable features on SIGW spectra via non-Gaussianity, with dynamics that can suppress or boost the signal amplitude for LISA.