pith. machine review for the scientific record. sign in

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

Recognition: 2 theorem links

· Lean Theorem

The NANOGrav 15-year Data Set: Search for Signals from New Physics

Adeela Afzal , Gabriella Agazie , Akash Anumarlapudi , Anne M. Archibald , Zaven Arzoumanian , Paul T. Baker , Bence B\'ecsy , Jose Juan Blanco-Pillado
show 115 more authors
Laura Blecha Kimberly K. Boddy Adam Brazier Paul R. Brook Sarah Burke-Spolaor Rand Burnette Robin Case Maria Charisi Shami Chatterjee Katerina Chatziioannou Belinda D. Cheeseboro Siyuan Chen Tyler Cohen James M. Cordes Neil J. Cornish Fronefield Crawford H. Thankful Cromartie Kathryn Crowter Curt J. Cutler Megan E. DeCesar Dallas DeGan Paul B. Demorest Heling Deng Timothy Dolch Brendan Drachler Richard von Eckardstein Elizabeth C. Ferrara William Fiore Emmanuel Fonseca Gabriel E. Freedman Nate Garver-Daniels Peter A. Gentile Kyle A. Gersbach Joseph Glaser Deborah C. Good Lydia Guertin Kayhan G\"ultekin Jeffrey S. Hazboun Sophie Hourihane Kristina Islo Ross J. Jennings Aaron D. Johnson Megan L. Jones Andrew R. Kaiser David L. Kaplan Luke Zoltan Kelley Matthew Kerr Joey S. Key Nima Laal Michael T. Lam William G. Lamb T. Joseph W. Lazio Vincent S. H. Lee Natalia Lewandowska Rafael R. Lino dos Santos Tyson B. Littenberg Tingting Liu Duncan R. Lorimer Jing Luo Ryan S. Lynch Chung-Pei Ma Dustin R. Madison Alexander McEwen James W. McKee Maura A. McLaughlin Natasha McMann Bradley W. Meyers Patrick M. Meyers Chiara M. F. Mingarelli Andrea Mitridate Jonathan Nay Priyamvada Natarajan Cherry Ng David J. Nice Stella Koch Ocker Ken D. Olum Timothy T. Pennucci Benetge B. P. Perera Polina Petrov Nihan S. Pol Henri A. Radovan Scott M. Ransom Paul S. Ray Joseph D. Romano Shashwat C. Sardesai Ann Schmiedekamp Carl Schmiedekamp Kai Schmitz Tobias Schr\"oder Levi Schult Brent J. Shapiro-Albert Xavier Siemens Joseph Simon Magdalena S. Siwek Ingrid H. Stairs Daniel R. Stinebring Kevin Stovall Peter Stratmann Jerry P. Sun Abhimanyu Susobhanan Joseph K. Swiggum Jacob Taylor Stephen R. Taylor Tanner Trickle Jacob E. Turner Caner Unal Michele Vallisneri Sonali Verma Sarah J. Vigeland Haley M. Wahl Qiaohong Wang Caitlin A. Witt David Wright Olivia Young Kathryn M. Zurek (for the NANOGrav Collaboration)
Authors on Pith no claims yet

Pith reviewed 2026-05-16 05:33 UTC · model grok-4.3

classification 🌌 astro-ph.HE astro-ph.COgr-qchep-ph
keywords gravitational wave backgroundpulsar timing arraysNANOGravcosmological sourcesphase transitionscosmic stringsultralight dark matter
0
0 comments X

The pith

The NANOGrav 15-year pulsar timing data can be reproduced by several cosmological gravitational wave sources, which often fit the observations better than supermassive black hole binaries.

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

The paper analyzes whether the low-frequency gravitational wave background detected in 15 years of NANOGrav pulsar timing observations could originate from cosmological processes. It tests models of cosmic inflation, scalar-induced gravitational waves, first-order phase transitions, cosmic strings, and domain walls against the data. All these models except stable cosmic strings of field theory origin are able to match the observed signal amplitude and spectrum. Several cosmological models yield better statistical fits than the standard interpretation in terms of inspiraling supermassive black hole binaries, with Bayes factors ranging from 10 to 100. The work also excludes regions of parameter space for these models and reports no evidence for signals from ultralight dark matter or Milky Way dark matter substructures.

Core claim

The 15-year NANOGrav data set shows positive evidence for a stochastic low-frequency gravitational-wave background that can be reproduced by cosmological models including cosmic inflation, scalar-induced gravitational waves, first-order phase transitions, and domain walls. Stable cosmic strings of field theory origin cannot reproduce the signal. When compared to the standard supermassive black hole binary interpretation, many cosmological models provide a better fit with Bayes factors from 10 to 100. Parameter regions are excluded where the predicted cosmological signals would exceed the observed strength, independent of the signal origin. No evidence is found for deterministic signals from

What carries the argument

Bayesian model comparison of cosmological gravitational wave production mechanisms against astrophysical supermassive black hole binary populations, applied to the NANOGrav 15-year pulsar timing array dataset.

If this is right

  • Cosmological sources such as inflation and phase transitions can account for the full observed gravitational wave background without requiring black hole binaries.
  • Certain ranges of parameters in these cosmological models are ruled out because they would generate a stronger background than detected.
  • No stochastic or deterministic signals from ultralight dark matter or galactic dark matter substructures are present in the data.
  • Constraints on ultralight dark matter couplings to electrons, muons, and gluons are tightened beyond previous limits from torsion balances and atomic clocks.

Where Pith is reading between the lines

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

  • Independent data from other pulsar timing arrays could test the cosmological interpretations without relying on black hole binary modeling choices.
  • Pulsar timing arrays now function as a direct probe of early-universe physics at energy scales far above those accessible to colliders.
  • If the signal is confirmed to be cosmological, it would require new early-universe mechanisms operating at temperatures around 10^9 to 10^12 GeV.

Load-bearing premise

The statistical preference for cosmological models over supermassive black hole binaries depends on specific assumptions about the black hole population and its evolution.

What would settle it

A future pulsar timing measurement that reveals a spectral shape or amplitude inconsistent with both the tested cosmological models and the black hole binary population would falsify the current interpretations.

read the original abstract

The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons.

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.

Referee Report

2 major / 3 minor

Summary. The manuscript analyzes the NANOGrav 15-year pulsar timing array data for cosmological interpretations of the detected stochastic gravitational-wave background. It examines models including cosmic inflation, scalar-induced gravitational waves, first-order phase transitions, cosmic strings, and domain walls, finding that all but stable field-theory cosmic strings can reproduce the signal. Several cosmological models yield Bayes factors of 10-100 relative to the standard supermassive black hole binary (SMBHB) interpretation. The paper also derives excluded parameter regions independent of signal origin and reports null results with updated constraints on ultralight dark matter and Milky Way dark matter substructures.

Significance. If the results hold, this work shows that pulsar timing data can place meaningful constraints on a wide range of new-physics models, including parameter exclusions that do not assume the origin of the NANOGrav signal. The explicit caveat on SMBHB modeling dependence and the competitive ultralight dark matter limits are strengths. The analysis follows standard Bayesian methods for timing residuals and provides falsifiable predictions for future data sets.

major comments (2)
  1. [Section 5 (Bayesian model comparison)] The Bayes-factor comparisons (reported as 10-100 in favor of several cosmological models) are computed against one fixed SMBHB population model. While the dependence on SMBHB assumptions is noted in the abstract and discussion, the manuscript does not quantify how the Bayes factors vary under alternative priors on merger rates, mass distributions, or eccentricities; this sensitivity analysis is load-bearing for the comparative claim.
  2. [Section 6 (parameter constraints)] In the parameter-exclusion analysis, the claim that excluded regions are independent of signal origin is central. The text should explicitly demonstrate that the upper limits remain unchanged when the NANOGrav spectrum is fixed to the best-fit SMBHB model versus a joint cosmological-plus-SMBHB fit.
minor comments (3)
  1. [Introduction] The introduction would benefit from a one-paragraph summary of the key 15-year data characteristics (e.g., number of pulsars, timing precision, and frequency range) to make the manuscript self-contained.
  2. [Figures 4-7] Figure captions for the posterior plots should explicitly state the credible-interval levels (68 % and 95 %) and the priors used for each cosmological parameter.
  3. [Section 3 (methods)] A few instances of inconsistent notation appear between the text and equations for the GW energy-density spectrum; a short notation table would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major point below and will revise the manuscript accordingly to strengthen the presentation of the results.

read point-by-point responses

  1. Referee: [Section 5 (Bayesian model comparison)] The Bayes-factor comparisons (reported as 10-100 in favor of several cosmological models) are computed against one fixed SMBHB population model. While the dependence on SMBHB assumptions is noted in the abstract and discussion, the manuscript does not quantify how the Bayes factors vary under alternative priors on merger rates, mass distributions, or eccentricities; this sensitivity analysis is load-bearing for the comparative claim.

    Authors: We agree that a quantitative sensitivity analysis would improve the robustness of the Bayes-factor comparisons. In the revised manuscript we will add a dedicated subsection (or appendix) that recomputes the Bayes factors under several alternative SMBHB priors, including variations in merger-rate evolution, black-hole mass functions, and eccentricity distributions drawn from the literature. We will report the resulting range of Bayes factors and discuss how the comparative preference for cosmological models changes (or does not change) under these alternatives. This will make the dependence on SMBHB modeling explicit and quantitative while preserving the existing caveat that the results should not yet be interpreted as evidence for new physics. revision: yes

  2. Referee: [Section 6 (parameter constraints)] In the parameter-exclusion analysis, the claim that excluded regions are independent of signal origin is central. The text should explicitly demonstrate that the upper limits remain unchanged when the NANOGrav spectrum is fixed to the best-fit SMBHB model versus a joint cosmological-plus-SMBHB fit.

    Authors: We thank the referee for this clarifying suggestion. In the revised version we will add an explicit comparison (new figure or table in Section 6) showing the excluded parameter regions obtained when the spectrum is fixed to the best-fit SMBHB model versus those obtained from a joint cosmological-plus-SMBHB fit. Because the exclusion criterion is that the predicted cosmological signal exceeds the observed NANOGrav power at any frequency, the limits are expected to be essentially unchanged; the new comparison will demonstrate this directly and thereby strengthen the claim of origin independence. revision: yes

Circularity Check

1 steps flagged

Bayes factors favoring cosmological models over SMBHBs depend on specific SMBHB population modeling assumptions

specific steps
  1. fitted input called prediction [Abstract]
    "When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics."

    The Bayes factor is obtained by comparing each cosmological model likelihood against a single SMBHB population model whose amplitude, spectral index, and other parameters are chosen to reproduce the observed NANOGrav signal. Because the reference model is fitted to the same data, the numerical preference (BF 10-100) is forced by the particular SMBHB parameterization adopted; alternative priors on merger rates or eccentricity distributions would alter the reference spectrum and thus the reported Bayes factors.

full rationale

The paper's central comparison computes Bayes factors between cosmological sources and a fixed SMBHB population model. The abstract explicitly states that cosmological models yield Bayes factors of 10-100 but immediately caveats that results 'strongly depend on modeling assumptions about the cosmic SMBHB population' and 'should not be regarded as evidence for new physics.' This matches the fitted-input-called-prediction pattern: the reference SMBHB spectrum is itself parameterized and fitted to the identical NANOGrav dataset, so the reported preference is conditional on that choice rather than an independent benchmark. No self-citation chain or definitional loop is present; the derivation of individual cosmological spectra is independent. The score is therefore moderate (4) rather than high, reflecting acknowledged model dependence without full reduction of the claim to tautology.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The analysis rests on standard stochastic GW background spectra for each cosmological source and on the assumption that the observed common-spectrum process is the only signal present; no new entities are postulated.

free parameters (2)
  • SMBHB population parameters
    Parameters describing the supermassive black hole binary population are chosen to match the observed amplitude and are varied to test robustness of cosmological model preference.
  • cosmological model parameters
    Parameters such as phase transition temperature, bubble wall velocity, string tension, and inflation scale are fitted or scanned to reproduce the NANOGrav signal spectrum.
axioms (2)
  • domain assumption The observed common-spectrum process is a stochastic gravitational wave background
    The paper assumes the timing residuals arise from a GW background rather than other noise processes when performing model comparisons.
  • standard math Standard forms for GW spectra from inflation, phase transitions, cosmic strings, and domain walls
    The paper adopts the usual analytic or semi-analytic expressions for the energy density spectra of each source class.

pith-pipeline@v0.9.0 · 6239 in / 1498 out tokens · 53425 ms · 2026-05-16T05:33:46.487952+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 19 Pith papers

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

  1. Probing the small-scale primordial power spectrum via relic neutrinos and acoustic reheating

    hep-ph 2026-05 unverdicted novelty 7.0

    Dissipation of small-scale primordial perturbations after neutrino decoupling cools relic neutrinos and reduces their abundance, enabling PTOLEMY to constrain the primordial curvature power spectrum to O(0.1) on scale...

  2. Tracing Primordial Gravitational Waves via non-Gaussian Signatures of Halo Bias

    astro-ph.CO 2026-05 unverdicted novelty 7.0

    Tensor-induced non-Gaussianity from primordial gravitational waves generates a unique scale-dependent halo bias correction that can reach order-one amplitude for rare high-redshift halos at z=7.

  3. 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.

  4. 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.

  5. Stone Skipping Black Holes in Ultralight Dark Matter Solitons

    astro-ph.CO 2026-02 unverdicted novelty 7.0

    Black holes in ultralight dark matter solitons undergo quasi-periodic stone-skipping orbits driven by soliton dipole excitations, modifying inspiral dynamics when the black hole is much lighter than the soliton.

  6. 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...

  7. Massive Scalar Quasinormal Modes, Greybody Factors, and Absorption Cross Section of a Parity-Symmetric Beyond-Horndeski Black Hole

    gr-qc 2026-05 unverdicted novelty 6.0

    Increasing the mass of a scalar field around a parity-symmetric beyond-Horndeski black hole strongly reduces the damping rate of quasinormal modes while suppressing low-frequency absorption and shifting efficient abso...

  8. 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.

  9. Imprint of domain wall annihilation on induced gravitational waves

    hep-ph 2026-04 unverdicted novelty 6.0

    Domain wall annihilation imprints a two-peaked spectrum on induced gravitational waves via an early matter-dominated phase and entropy dilution.

  10. Irreducible Gravitational Wave Background as a Particle Detector

    hep-ph 2026-04 unverdicted novelty 6.0

    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.

  11. 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...

  12. Two types of quasinormal modes of Casadio-Fabbri-Mazzacurati brane-world black holes

    gr-qc 2026-02 unverdicted novelty 6.0

    Quasinormal modes of massive scalars in CFM brane-world black holes split into two types, with modes disappearing at critical masses where real or imaginary frequency parts reach zero.

  13. Long-lived massive scalar modes, grey-body factors, and absorption cross sections of the Reissner--Nordstr\"om-like brane-world black hole

    gr-qc 2026-05 unverdicted novelty 5.0

    Positive tidal charge in this brane-world black hole lowers the effective potential barrier, pushes massive scalar quasinormal modes toward arbitrarily long lifetimes, and increases transmission and absorption.

  14. 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...

  15. Long-lived quasinormal modes, shadows and particle motion in four-dimensional quasi-topological gravity

    gr-qc 2026-03 unverdicted novelty 5.0

    Massive scalar quasinormal modes in quasi-topological black holes become long-lived as scalar mass grows, while photon-sphere radius, shadow size, and ISCO exhibit moderate deviations from Schwarzschild.

  16. F-Term Hybrid Inflation with T-Model K\"ahler Geometry and Beyond

    hep-ph 2026-05 unverdicted novelty 4.0

    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...

  17. Metastable strings at PTAs: classical stability analysis

    hep-ph 2026-05 unverdicted novelty 4.0

    Classical instabilities in metastable strings from two-step symmetry breaking can restrict the viable parameter space for explaining the PTA gravitational wave signal.

  18. 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.

  19. Massive scalar quasinormal modes of an asymptotically flat regular black hole supported by a phantom Dirac--Born--Infeld field

    gr-qc 2026-04 unverdicted novelty 4.0

    Massive scalar quasinormal modes in this DBI-supported regular black hole show higher oscillation frequencies and lower damping as field mass increases, with larger regularity scales producing softer and longer-lived ringing.