Recognition: 2 theorem links
· Lean TheoremThe NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background
Pith reviewed 2026-05-12 04:03 UTC · model grok-4.3
The pith
The NANOGrav 15-year pulsar timing data set shows evidence for a stochastic gravitational-wave background whose spatial correlations match the Hellings-Downs pattern.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Analysis of the NANOGrav 15-year data set reveals a stochastic signal that is correlated among the pulsars according to the Hellings-Downs pattern expected for a gravitational-wave background. The presence of a power-law-spectrum background is favored over independent pulsar noises with a Bayes factor in excess of 10^14 and over an uncorrelated common power-law spectrum with Bayes factors of 200-1000. A frequentist test based on a weighted sum of inter-pulsar correlations gives a p-value between 5 times 10^-5 and 1.9 times 10^-4. Under the assumption of a characteristic-strain spectrum proportional to f^-2/3, the median amplitude is 2.4 times 10^-15 at a reference frequency of one cycle per
What carries the argument
The Hellings-Downs pattern of inter-pulsar correlations, which is the unique angular dependence produced by an isotropic stochastic gravitational-wave background; it is used to distinguish the signal from other common processes through both Bayesian model comparison and a direct frequentist correlation statistic.
If this is right
- The amplitude and spectral index are consistent with the superposition of gravitational waves from an ensemble of inspiraling supermassive black-hole binaries.
- More exotic cosmological sources remain possible but are not required by the current data.
- Continued monitoring will tighten constraints on the background spectrum and begin to resolve individual sources.
- The detection establishes pulsar timing arrays as a functioning observatory for nanohertz gravitational waves.
Where Pith is reading between the lines
- Independent confirmation by other pulsar timing arrays would raise the significance and help isolate any array-specific systematics.
- If the background is dominated by black-hole binaries, the measured amplitude supplies a new integral constraint on the high-redshift black-hole mass function and merger rate.
- With more pulsars the same data set could begin to search for anisotropy or for deviations from the pure power-law spectrum that would indicate additional source populations.
Load-bearing premise
That the measured correlations between pulsars arise from gravitational waves rather than from any unidentified instrumental systematic or incomplete description of the individual pulsar noise processes.
What would settle it
A future data release with additional pulsars in which the measured correlation as a function of angular separation deviates from the Hellings-Downs curve by more than the combined statistical and systematic uncertainty, or the identification of a common non-gravitational process that reproduces the same correlation pattern.
read the original abstract
We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3\sigma$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4\sigma$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports multiple lines of evidence for a stochastic gravitational-wave background in the NANOGrav 15-year pulsar-timing data set from 67 pulsars. The observed inter-pulsar correlations follow the Hellings-Downs pattern. A power-law GW background model is favored over independent pulsar noise models with Bayes factor >10^{14} and over an uncorrelated common power-law model with Bayes factors 200-1000. A data-driven null distribution (inter-pulsar correlations removed) yields p≈10^{-3} (∼3σ); a frequentist weighted correlation statistic yields p=5×10^{-5} to 1.9×10^{-4} (∼3.5-4σ). The characteristic strain amplitude is reported as 2.4^{+0.7}_{-0.6}×10^{-15} at 1 yr^{-1}, consistent with expectations from a supermassive black-hole binary population.
Significance. If the central claim holds, this constitutes the first observational evidence for a nanohertz stochastic gravitational-wave background. The result opens a new observational window complementary to LIGO/Virgo and has direct implications for supermassive black-hole binary demographics and possible cosmological sources. Strengths include the use of multiple independent statistical tests (Bayes factors, constructed null distribution, and frequentist statistic) and explicit comparison against both noise-only and uncorrelated common-spectrum models. The manuscript also provides a falsifiable amplitude prediction tied to an astrophysical population.
minor comments (2)
- [Abstract] Abstract and § on model comparison: the reported Bayes-factor range (200-1000) is stated to depend on spectral modeling choices; explicit listing of the exact models (e.g., broken power-law, free spectral, etc.) and their associated Bayes factors would improve reproducibility.
- The construction of the null distribution by removing inter-pulsar correlations is central to the p-value claim; a brief statement confirming that the procedure preserves the marginal posterior widths on the common-signal hyperparameters would further address potential covariance concerns.
Simulated Author's Rebuttal
We thank the referee for their positive review of our manuscript and their recommendation to accept. We appreciate the recognition of the multiple independent statistical tests and the astrophysical implications of the result.
Circularity Check
No significant circularity; derivation uses external theoretical template and data-driven null
full rationale
The paper's evidence chain compares observed pulsar correlations against the Hellings-Downs pattern (an independent GR prediction) and against models with no or uncorrelated common noise. Bayes factors are computed directly from the likelihoods under these fixed correlation structures. The null distribution is generated by explicitly removing inter-pulsar correlations from the same dataset, yielding an empirical p-value. No equation, parameter fit, or model definition reduces the reported Bayes factors or amplitude to a tautology of the inputs. Prior NANOGrav citations supply only analysis methods and are not invoked to establish uniqueness or to substitute for the present data's statistical tests. The central result therefore remains independent of its own fitted values.
Axiom & Free-Parameter Ledger
free parameters (1)
- characteristic strain amplitude =
2.4e-15
axioms (1)
- standard math An isotropic stochastic gravitational-wave background produces the Hellings-Downs spatial correlation pattern in pulsar timing residuals.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/DimensionForcingalexander_duality_circle_linking echoesThe correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background.
Forward citations
Cited by 29 Pith papers
-
Gravity Echoes from Supermassive Black Hole Binaries
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.
-
Phase-resolved field-space distance bounds in ekpyrotic, bouncing and cyclic cosmologies
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.
-
Self-acceleration of Hardening Binaries
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 ...
-
Forecasting graviton-mass constraints from the full covariance of PTA-astrometry ORF estimators
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
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.
-
Detecting Chiral Gravitational Wave Background with a Dipole Pulsar Timing Array
A dipole pulsar timing array detects chiral nanohertz gravitational waves and extends PTA sensitivity into the microhertz regime.
-
Gravitational wave signal and noise response of an optically levitated sensor in a Fabry-P\'erot cavity
A general relativistic derivation of gravitational wave response in an optically levitated cavity sensor reveals position-dependent strain sensitivity and suppressed input-mirror noise coupling.
-
Testing General Relativity with Individual Supermassive Black Hole Binaries
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?
PTA statistical tests lose sensitivity to non-Gaussian GW features after decorrelation and cannot distinguish them model-agnostically.
-
High-Power AM-CW Lunar Laser Ranging as a $\mu$Hz SGWB Detector
AM-CW lunar laser ranging achieves μHz SGWB sensitivity of 5.29×10^{-9} D_cov (80 μm range uncertainty) or 2.07×10^{-9} D_cov (50 μm) over 5 years, with discovery possible if covariance degradation stays below ~3.6-13.7.
-
Gravitational Waves from a Black Hole Falling Radially into a Thin-Shell Traversable Wormhole
Analytic gravitational waveforms from radial test-particle infall into a thin-shell traversable wormhole exhibit a characteristic pulse-gap structure from repeated throat crossings and lie within reach of ground-based...
-
Imprint of domain wall annihilation on induced gravitational waves
Domain wall annihilation imprints a two-peaked spectrum on induced gravitational waves via an early matter-dominated phase and entropy dilution.
-
Spectral Butterfly Effect and Resilient Ringdown in Thick Braneworlds
Thick braneworlds feature fragile quasinormal mode spectra due to a butterfly effect but maintain a resilient early ringdown, keeping the standard gravitational wave fingerprint usable.
-
Irreducible Gravitational Wave Background as a Particle Detector
Spectral features imprinted by long-lived BSM particles on any primordial GWB directly determine the particles' mass and decay rate once the model and initial abundance are specified.
-
The Heavy Tailed Non-Gaussianity of the Supermassive Black Hole Gravitational Wave Background
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
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
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
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.
-
Purely Quadratic Non-Gaussianity from Tachyonic Instability: Primordial Black Holes and Scalar-Induced Gravitational Waves
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
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...
-
Reconstructing inflationary features on large scales using genetic algorithm
Genetic algorithm reconstructs single-field inflationary models with features in the scalar power spectrum that fit Planck 2018 CMB data better by Δχ² ≲ -10 and suggest alternative background parameters.
-
Grand Unified Origin of Enhanced Scalar Couplings: Connecting Radiative Electroweak Symmetry Breaking to SO(10) Dynamics
SO(10) scalar threshold corrections generate the Higgs quartic enhancement factor k≈6 required by radiative electroweak symmetry breaking, placing the Landau pole at 1.5-2×10^16 GeV near the GUT scale.
-
F-Term Hybrid Inflation with T-Model K\"ahler Geometry and Beyond
F-term hybrid inflation with SU(1,1)/U(1) or SU(2)/U(1) Kähler geometry in GUTs can be realized without inflationary extrema for broad parameters, matching ACT/SPT data via curvature and tadpole adjustments while pred...
-
Constraints on Ultralight Scalar and Dark Photon Dark Matter from PPTA-DR3 and EPTA-DR2
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.
-
Constraints on Einstein-aether gravity from the precision timing of PSR J1738+0333
Precision timing of PSR J1738+0333 from EPTA and NANOGrav data yields the tightest strong-field constraints on Einstein-aether parameters from any single binary pulsar.
-
Precision Analysis for $\boldsymbol{H_0}$ Using Upcoming Multi-band Gravitational Wave Observations
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.
-
Colloquium: Radio astronomy with the Arecibo 305-m telescope: In contemporaneous context
The Arecibo telescope advanced radio astronomy via serendipitous discoveries and successive upgrades to its reflector, optics, receivers, and data systems over its operational lifetime.
-
The CosmoVerse White Paper: Addressing observational tensions in cosmology with systematics and fundamental physics
The CosmoVerse White Paper compiles observational tensions in cosmology and maps strategies using improved systematics checks and tests of fundamental physics to resolve them.
-
A short course in general relativity
Lecture notes that derive the nonlinear Einstein field equations from the gravitational action and apply them to black holes, cosmology, and gravitational waves at the advanced undergraduate level.
Reference graph
Works this paper leans on
-
[1]
Abbott, B. P., Abbott, R., Abbott, T. D., et al. 2016, PhRvL, 116, 061102, doi: 10.1103/PhysRevLett.116.061102 Afzal, A., et al. 2023, in preparation, doi: 10.3847/2041-8213/acdc91 Agazie, G., et al. 2023a, in preparation —. 2023b, in preparation, doi: 10.3847/2041-8213/acda9a —. 2023c, in preparation, doi: 10.3847/2041-8213/acda88 Aggarwal, K., Arzoumani...
-
[2]
2022, MNRAS, 509, 5538, doi: 10.1093/mnras/stab3283 Chamberlin, S
http://www.jstor.org/stable/2346151 Chalumeau, A., Babak, S., Petiteau, A., et al. 2022, MNRAS, 509, 5538, doi: 10.1093/mnras/stab3283 Chamberlin, S. J., Creighton, J. D. E., Siemens, X., et al. 2015, PhRvD, 91, 044048, doi: 10.1103/PhysRevD.91.044048 Chen, S., Caballero, R. N., Guo, Y. J., et al. 2021, MNRAS, 508, 4970, doi: 10.1093/mnras/stab2833 Cordes...
-
[3]
2021, Universe, 7, 398, doi: 10.3390/universe7110398 DuPlain, R., Ransom, S., Demorest, P., et al
http://www.jstor.org/stable/2958475 Dom` enech, G. 2021, Universe, 7, 398, doi: 10.3390/universe7110398 DuPlain, R., Ransom, S., Demorest, P., et al. 2008, in Proc. SPIE, Vol. 7019, Advanced Software and Control for Astronomy II, 70191D, doi: 10.1117/12.790003 Einstein, A. 1916, Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys.), 1916, 1 Ellis, J., & v...
- [4]
- [5]
-
[6]
http://www.jstor.org/stable/1391010 Goncharov, B., Shannon, R. M., Reardon, D. J., et al. 2021a, ApJL, 917, L19, doi: 10.3847/2041-8213/ac17f4 Goncharov, B., Reardon, D. J., Shannon, R. M., et al. 2021b, MNRAS, 502, 478, doi: 10.1093/mnras/staa3411 Goncharov, B., Thrane, E., Shannon, R. M., et al. 2022, ApJL, 932, L22, doi: 10.3847/2041-8213/ac76bb Guzzet...
-
[7]
https://arxiv.org/abs/1908.05356 Reardon, D. J., et al. 2023, in preparation Roebber, E. 2019, ApJ, 876, 55, doi: 10.3847/1538-4357/ab100e Roebber, E., & Holder, G. 2017, ApJ, 835, 21, doi: 10.3847/1538-4357/835/1/21 Romani, R. W. 1989, in Timing Neutron Stars, ed. H. ¨Ogelman & E. P. J. Heuvel (Springer), 113–117 Romano, J. D., Hazboun, J. S., Siemens, X...
-
[8]
Markers indicate which telescopes observed the pulsar
Top: Sky locations of the 67 pulsars used in the 15-year GWB analysis. Markers indicate which telescopes observed the pulsar. Bottom: Distribution of angular sepa- rations probed by the pulsars in the full data set (orange), the Arecibo data set (blue), and the GBT data set (red). Because Arecibo and GBT mostly observed pulsars at differ- ent declinations...
work page 1992
-
[9]
and postprocessed sam- ples with chainconsumer (Hinton 2016). In NG12gwb we use an analytic approximation for the uncertainty of marginalized-posterior statistics (Wilcox 2012). Here we instead adopt a boostrap approach: we resample the original MCMC samples (with replace- ment) to generate new sets that act as independent sam- pling realizations. We then...
work page 2016
-
[10]
fixed fbend [Hz] broken–power-law bend frequency log-Uniform [ −8.7,−7] one parameter for PTA ℓ broken–power-law high-freq. transition sharpness delta function ( ℓ = 0.1) fixed all common processes, t-process spectrum A power-law amplitude log-Uniform [ −18, −11] one parameter for PTA γ power-law spectral index Uniform [0 , 7] one parameter per PTA xi mod...
work page 1989
-
[11]
shows that the frequency bins at f1, f6, f7, and f8 appear to be in tension with a pure power law, skew- ing the estimation of γ and reducing the hd13/3 vs. curn13/3 Bayes factor. Assuming that those frequency components reflect unmodeled systematics or stronger- than-expected statistical fluctuations, we can make our inference more robust to such outlier...
work page 2022
-
[12]
and compared to curnfree bin variances. The TPS model is spread more widely and deviates from the perfect power law at bins f1, f6, f7, and f8, as expected. The right panel of Figure 13 shows the joint log10 A, γ posteriors for curnγ and curnTPS. The latter is more consistent with steeper power laws, and it includes γ = 13/3 at 1 σ credibility. E. TURNOVE...
work page 2002
-
[13]
Power-law (curnγ, blue) and t-process power-law ( curnTPS, orange) spectral posteriors. Left: reconstructed spectra, compared to free-spectral bin-variance posteriors ( curnTPS, violin plots). Right: joint (log 10 A, γ) posteriors. The “fuzzy” t-process allows local deviations from a perfect power law, producing wider constraints that are more consistent ...
work page 2016
-
[14]
The linear-model co- efficients are the squared amplitudes of the components
fits the inter-pulsar corre- lation coefficients ρab with a linear model that includes multiple components with different correlation patterns, but with the same spectral shape. The linear-model co- efficients are the squared amplitudes of the components. Within such a model, the significance of each component can be quoted as a S/N given by its best-fit ...
work page 2018
-
[15]
Comparison between empirical background dis- tributions for the noise-marginalized optimal statistic, as computed by the sky-scramble technique. We show distri- butions computed using a match threshold of ¯M < 0.17 (blue), ¯M < 0.1 (orange), and ¯M < 0.08 (green). Dotted lines indicate Gaussian-equivalent 2σ, 3σ, and 4σ thresholds. The dashed vertical lin...
work page 2023
-
[16]
T op: MCOS HD S/N values recovered in the three simulations described in App. H.1, compared to the MCOS HD S/N measured in the real data set (vertical dashed red line), which has p-values of < 10−2 for simula- tions i and iii, and 0 .64 for simulations ii. Bottom: MCOS monopole S/N values recovered in the three simulations, compared to the real-data MCOS ...
work page 2021
-
[17]
Distribution of real-data and replicated MCOS monopole S/Ns. Each point represents a draw η(k) from hd13/3 posterior, which is used to simulate δtsim,(k) and to compute both S/Ns. The replicated monopole S/N is greater for 11% of the simulations. lations: (i) injecting no spatially correlated power-law GWB or excess uncorrelated common-spectrum noise; (ii...
work page 1996
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.