pith. machine review for the scientific record. sign in

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

Recognition: 2 theorem links

· Lean Theorem

The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

Aaron D. Johnson, Abhimanyu Susobhanan, Adam Brazier, Akash Anumarlapudi, Alexander McEwen, Andrea Lommen, Andrea Mitridate, Andrew R. Kaiser, Anne M. Archibald, Ann Schmiedekamp, Belinda D. Cheeseboro, Bence Becsy, Benetge B. P. Perera, Bradley W. Meyers, Brendan Drachler, Brent J. Shapiro-Albert, Caitlin A. Witt, Caner Unal, Carl Schmiedekamp, Cherry Ng, Chiara M. F. Mingarelli, Chung-Pei Ma, Curt J. Cutler, Dallas DeGan, Daniel R. Stinebring, David J. Nice, David L. Kaplan, Deborah C. Good, Duncan R. Lorimer, Dustin R. Madison, Elizabeth C. Ferrara, Emmanuel Fonseca, Fronefield Crawford, Gabriel E. Freedman, Gabriella Agazie, Haley M. Wahl, Heling Deng, Henri A. Radovan, H. Thankful Cromartie, Ingrid H. Stairs, Jacob E. Turner, Jacob Taylor, James M. Cordes, James W. McKee, Jeffrey S. Hazboun, Jerry P. Sun, Jing Luo, Joey S. Key, Joseph D. Romano, Joseph Glaser, Joseph K. Swiggum, Joseph Simon, Justin A. Ellis, Kai Schmitz, Katerina Chatziioannou, Kathryn Crowter, Kayhan Gultekin, Ken D. Olum, Kevin Stovall, Kristina Islo, Kyle A. Gersbach, Laura Blecha, Levi Schult, Luke Zoltan Kelley, Magdalena S. Siwek, Margaret A. Mattson, Maria Charisi, Matthew Kerr, Maura A. McLaughlin, Megan E. DeCesar, Megan L. Jones, Michael T. Lam, Michele Vallisneri, Natalia Lewandowska, Natasha McMann, Nate Garver-Daniels, Neil J. Cornish, Nihan S. Pol, Nima Laal, Olivia Young (The NANOGrav Collaboration), Patrick M. Meyers, Paul B. Demorest, Paul R. Brook, Paul S. Ray, Paul T. Baker, Peter A. Gentile, Polina Petrov, Priyamvada Natarajan, Qiaohong Wang, Rand Burnette, Robin Case, Ross J. Jennings, Rutger van Haasteren, Ryan S. Lynch, Sarah Burke-Spolaor, Sarah J. Vigeland, Scott M. Ransom, Shami Chatterjee, Shashwat C. Sardesai, Siyuan Chen, Sophie Hourihane, Stella Koch Ocker, Stephen R. Taylor, Timothy Dolch, Timothy T. Pennucci, Tingting Liu, T. Joseph W. Lazio, Tonia C. Klein, Tyler Cohen, Tyson B. Littenberg, William Fiore, William G. Lamb, Xavier Siemens, Zaven Arzoumanian

Authors on Pith no claims yet

Pith reviewed 2026-05-12 04:03 UTC · model grok-4.3

classification 🌌 astro-ph.HE gr-qc
keywords gravitational wavespulsar timing arraysstochastic backgroundsupermassive black holesHellings-Downs correlationsBayesian model selectionnanohertz frequencies
0
0 comments X

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.

The paper examines 15 years of arrival-time data from 67 millisecond pulsars to search for a common low-frequency signal. It finds that the measured correlations between different pulsars closely follow the angular dependence predicted for an isotropic gravitational-wave background. Statistical model comparison yields Bayes factors exceeding 10^14 in favor of a correlated power-law spectrum over models containing only independent pulsar noise, and factors of several hundred over an uncorrelated common spectrum. The inferred strain amplitude at a one-year reference frequency is consistent with the expected contribution from a large population of distant supermassive black-hole binaries. A sympathetic reader would care because the result supplies the first direct evidence that nanohertz gravitational waves permeate the universe and can be accessed through pulsar timing.

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

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

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

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

0 major / 2 minor

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)
  1. [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.
  2. 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

0 responses · 0 unresolved

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

0 steps flagged

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

1 free parameters · 1 axioms · 0 invented entities

The claim rests on the standard assumption that the Hellings-Downs pattern is the unique signature of an isotropic gravitational-wave background and on the completeness of the pulsar noise models; the amplitude is a fitted parameter under a fiducial spectral index.

free parameters (1)
  • characteristic strain amplitude = 2.4e-15
    Fitted to the data under the assumed f^{-2/3} spectrum; reported as 2.4^{+0.7}_{-0.6} times 10^{-15} at 1/yr.
axioms (1)
  • standard math An isotropic stochastic gravitational-wave background produces the Hellings-Downs spatial correlation pattern in pulsar timing residuals.
    Invoked as the distinguishing observable that identifies the signal as gravitational-wave in origin.

pith-pipeline@v0.9.0 · 6237 in / 1503 out tokens · 53941 ms · 2026-05-12T04:03:56.095337+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.

Forward citations

Cited by 29 Pith papers

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

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

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

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

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

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

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

  7. Gravitational wave signal and noise response of an optically levitated sensor in a Fabry-P\'erot cavity

    gr-qc 2026-04 unverdicted novelty 7.0

    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.

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

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

  10. High-Power AM-CW Lunar Laser Ranging as a $\mu$Hz SGWB Detector

    gr-qc 2026-05 unverdicted novelty 6.0

    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.

  11. Gravitational Waves from a Black Hole Falling Radially into a Thin-Shell Traversable Wormhole

    gr-qc 2026-05 unverdicted novelty 6.0

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

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

  13. Spectral Butterfly Effect and Resilient Ringdown in Thick Braneworlds

    gr-qc 2026-04 unverdicted novelty 6.0

    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.

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

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

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

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

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

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

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

  21. Reconstructing inflationary features on large scales using genetic algorithm

    astro-ph.CO 2026-04 unverdicted novelty 5.0

    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.

  22. Grand Unified Origin of Enhanced Scalar Couplings: Connecting Radiative Electroweak Symmetry Breaking to SO(10) Dynamics

    hep-ph 2026-04 unverdicted novelty 5.0

    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.

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

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

  25. Constraints on Einstein-aether gravity from the precision timing of PSR J1738+0333

    gr-qc 2026-05 unverdicted novelty 4.0

    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.

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

  27. Colloquium: Radio astronomy with the Arecibo 305-m telescope: In contemporaneous context

    astro-ph.GA 2026-04 unverdicted novelty 2.0

    The Arecibo telescope advanced radio astronomy via serendipitous discoveries and successive upgrades to its reflector, optics, receivers, and data systems over its operational lifetime.

  28. The CosmoVerse White Paper: Addressing observational tensions in cosmology with systematics and fundamental physics

    astro-ph.CO 2025-04 accept novelty 2.0

    The CosmoVerse White Paper compiles observational tensions in cosmology and maps strategies using improved systematics checks and tests of fundamental physics to resolve them.

  29. A short course in general relativity

    gr-qc 2026-04 unverdicted

    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

17 extracted references · 17 canonical work pages · cited by 29 Pith papers

  1. [1]

    P., et al

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

    http://www.jstor.org/stable/24306036 Gelman, A., & Rubin, D. B. 1992, Statistical Science, 7,

  5. [5]

    http://www.jstor.org/stable/2246093 Godsill, S. J. 2001, Journal of Computational and Graphical Statistics, 10,

  6. [6]

    M., Reardon, D

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

    J., et al

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

  9. [9]

    In NG12gwb we use an analytic approximation for the uncertainty of marginalized-posterior statistics (Wilcox 2012)

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

  10. [10]

    freezing

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

  11. [11]

    curn13/3 Bayes factor

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

  12. [12]

    The TPS model is spread more widely and deviates from the perfect power law at bins f1, f6, f7, and f8, as expected

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

  13. [13]

    Left: reconstructed spectra, compared to free-spectral bin-variance posteriors ( curnTPS, violin plots)

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

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

  15. [15]

    We show distri- butions computed using a match threshold of ¯M < 0.17 (blue), ¯M < 0.1 (orange), and ¯M < 0.08 (green)

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

  16. [16]

    astrophysical,

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

  17. [17]

    Each point represents a draw η(k) from hd13/3 posterior, which is used to simulate δtsim,(k) and to compute both S/Ns

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