Scalar-Induced Gravitational Waves and PBHs
Scrollable formula. Use Left and Right arrow keys to reveal hidden content.
Scrollable table. Use Left and Right arrow keys to reveal additional columns.
Scrollable figure. Use Left and Right arrow keys to inspect hidden content.
Enhanced primordial scalar perturbations source gravitational waves at second order and can also collapse into primordial black holes (PBHs). The shared origin makes the pair attractive: a stochastic spectrum probes relatively low moments of the scalar statistics, while PBH abundance probes a rare nonlinear tail. It also makes careless claims easy. Through 2026-08-20, the production mechanism and band-specific constraints are mature, and scalar-induced templates can fit the nanohertz background; neither the background’s origin nor a linked PBH population is identified. Non-Gaussianity, collapse physics, thermal history, and astrophysical source competition keep the evidence contested.
Helpful background. Scalar-Induced Gravitational Waves derives the second-order tensor source and transfer kernels. This assessment asks whether one specified primordial probability distribution predicts both the observed-background likelihood and an allowed PBH population.
One primordial distribution, two observables
Section titled “One primordial distribution, two observables”The finite proposition is whether a declared scalar probability distribution and expansion history jointly predict a scalar-induced gravitational-wave spectrum and a PBH mass function that fit a measured or bounded stochastic-background likelihood and every independent abundance constraint. Fitting a free peaked spectrum is not enough; the same parameters, non-Gaussian statistics, collapse threshold, and thermal history must generate the PBHs.
The decisive axes are power-spectrum shape and width; Gaussian versus fully specified non-Gaussian statistics; collapse threshold, profile, smoothing window, and critical scaling; radiation domination versus a nonstandard expansion or reheating phase; PBH mass and abundance evolution; PTA or interferometer response; pulsar correlations and red-noise models; astrophysical foregrounds; and compatibility versus model evidence.
| Assessment at 2026-08-20 | Meaning |
|---|---|
| Question status | Partially resolved: production theory and band-specific constraints are mature, while observational attribution and joint PBH abundance remain open |
| Evidence status | Contested: PTA data admit scalar-induced spectra but also credible astrophysical and other primordial sources, and PBH consistency is assumption sensitive |
| Confidence | Moderate; high for second-order production, moderate for band-limited constraints, low for PBH-linked source attribution |
| Domain | Primordial scalar-induced backgrounds in PTA and ground-interferometer bands, with explicitly stated scalar statistics, expansion history, and PBH formation model |
| What could change the conclusion | A new PTA or interferometer release, a joint non-Gaussian formation calculation tied to public chains, or a correction to the abundance–background mapping |
A starting project: reproduce the abundance bridge
Section titled “A starting project: reproduce the abundance bridge”This direction suits students comfortable with cosmological perturbations, numerical integrals, and Gaussian tail probabilities. Begin with the linked induced-wave derivation, then read the scalar-induced model in Afzal et al. 2023, §V.2 and Appendix C.2, Eqs. (C7)–(C9). Extract the scalar spectrum, smoothing window, radiation-era transfer function, threshold, and formation-to-present conversion as separate inputs. Keep the paper’s PBH estimate distinct from a collapse simulation.
For one of its specified scalar-spectrum shapes, implement the variance integral and compare the exact Gaussian tail integral with its asymptotic approximation. Reproduce the paper’s abundance prescription for a declared parameter point before scanning amplitude or width. Report quadrature and integration-range convergence independently of uncertainty in the collapse threshold. Then display how changing the threshold within a stated diagnostic range affects abundance while holding the input scalar spectrum fixed. This is a sensitivity calculation, not a consistent change of collapse physics unless the profile and window are recalibrated together.
The output should be an executable abundance calculation and a plot explaining when the tail approximation ceases to be adequate. A laptop suffices for this baseline; a joint timing analysis needs pinned PTA products and inference expertise, and calibrated collapse thresholds need a specialist numerical collaborator. A candidate extension is a fully specified non-Gaussian distribution propagated consistently into both observables. Check the literature first: inserting one skewness correction into an exponentially rare tail is not, by itself, a new controlled prediction.
Pivot to a narrower uncertainty study if the tail expansion is uncontrolled or the required collapse calibration is missing. Do not proceed to interpreting PTA posterior overlap as PBH evidence while these inputs remain unresolved. The successful first project identifies the weakest step in the abundance mapping and a concrete calculation needed to strengthen it.
What theory and present data establish
Section titled “What theory and present data establish”At second perturbative order, products of scalar perturbations source transverse tensor modes. The spectrum is a convolution of scalar modes weighted by time-dependent transfer kernels, so its shape depends on the primordial spectrum and expansion history. Domènech reviews the gauge issues, radiation and nonstandard eras, observables, and PBH connection in §§2–8 2021, §§2–8. This establishes a mechanism, not an observed source identity.
The PBH bridge is more nonlinear. Abundance depends exponentially on a rare tail and can change drastically under non-Gaussianity even when the induced-wave spectrum changes modestly. Pi treats non-Gaussian statistics, abundance, and the joint induced-wave interpretation in §§2–7 2025, §§2–7. Specifying only a power spectrum and one perturbative non-Gaussian coefficient is insufficient when the tail expansion itself fails.
NANOGrav’s 15-year correlation analysis provides direct evidence for a Hellings–Downs-correlated nanohertz background Agazie et al. 2023, §§III–VI and Fig. 3. A separate new-physics analysis presents the source comparisons, scalar-induced sector, and PBH caveats in §3, §5.2, §7, and Appendix C.2 Afzal et al. 2023, §3, §5.2, §7, and Appendix C.2. Selected scalar-induced templates can fit the data, but that paper does not claim evidence for new physics; supermassive-black-hole binaries and other templates remain credible.
In the ground band, the current LVK early-Universe analysis uses O1–O4a cross-correlation data, finds no stochastic detection, and constrains prescribed second-order scalar spectra LVK 2025, §II, scalar-perturbation section and Fig. 11. It probes scalar scales far from the PTA range. The null cannot be transported across frequency without a physical spectrum and thermal transfer.
Results and open questions at 2026-08-20
Section titled “Results and open questions at 2026-08-20”| Claim or evidence | What it supports | What it does not establish |
|---|---|---|
| Second-order tensor derivation | A robust scalar-induced production mechanism in a specified background | That an observed stochastic signal has this origin |
| Non-Gaussian PBH analyses | Strong dependence of abundance–wave consistency on the full scalar statistics | A universal conversion from wave amplitude to PBH fraction |
| NANOGrav 15-year Hellings–Downs correlations | Direct evidence for a nanohertz gravitational-wave background | Its source, primordial origin, or a linked PBH population |
| NANOGrav 15-year source fits | Compatibility of selected scalar-induced templates with that background likelihood | Independent detection evidence, preference over all astrophysical sources, or proof of PBHs |
| LVK O1–O4a ground-band null | Bounds on prescribed scalar spectra mapped into the ground band | Exclusion of nanohertz models with different scale support |
| Current PBH synthesis | A conditional joint landscape and forecast targets | Primordial attribution of the PTA background |
| PBH-reheating adverse recalculation | Large suppression of one forecast signal after including irreducible mass spread | Refutation of ordinary radiation-era scalar-induced waves |
The rows are complementary, not competing measurements of one amplitude. PTA and LVK sample vastly separated scalar scales. A source comparison on one PTA dataset supplies alternative explanations but not independent observations; a distinct frequency band provides independent data but only after a model maps between bands.
Source degeneracy and the rare-tail problem
Section titled “Source degeneracy and the rare-tail problem”The optimistic interpretation says that the observed nanohertz shape can be produced by an enhanced primordial scalar spectrum and that the same enhancement naturally predicts PBHs. The current PBH synthesis describes this promising connection while emphasizing non-Gaussianity, formation, and competing sources Carr et al. 2026, §VII.A and Figs. 7–8. A fit demonstrates compatibility in its parameterization.
The skeptical interpretation has two independent parts. First, source degeneracy: supermassive-black-hole binaries and several primordial templates can fit the same timing data, with results sensitive to red-noise, ephemeris, spectral, and population priors. Second, rare-tail degeneracy: two scalar distributions with similar two-point power can give radically different PBH abundances. The missing discriminator is a joint model that predicts additional spectral structure, non-Gaussian correlations, or an independently observed PBH population.
Counterexamples, null results, and failed routes
Section titled “Counterexamples, null results, and failed routes”The LVK null is a counterexample to arbitrarily extending a large scalar spectrum into the ground band, but not to a narrow nanohertz feature. Conversely, a PTA fit is not a counterexample to an astrophysical background because source models reuse the same timing data and remain viable.
A concrete failed forecast appears in the PBH-dominated reheating or “Poltergeist” scenario. Gouttenoire, Leister, and Schwaller find that an irreducible PBH mass spread can suppress a previously forecast induced component by orders of magnitude 2026, §§2–7. This is an unrefereed calculation in a distinct mechanism; it warns that idealized monochromatic assumptions can fail, but it does not refute standard radiation-era scalar-induced production.
What would change the assessment
Section titled “What would change the assessment”The attribution would strengthen if a new PTA release resolved spectral or correlation structure predicted before unblinding and disfavored flexible supermassive-binary populations. A joint calculation using a full non-Gaussian distribution, validated collapse simulations, and public timing chains could turn compatibility into a sharper PBH statement. Detection or exclusion of the same physical spectrum in a second frequency band would be decisive only if the thermal transfer is fixed. A material revision to PTA noise or the collapse/abundance mapping could instead remove the current overlap.
Research checklist
Section titled “Research checklist”- Specify the complete scalar probability distribution, not only the two-point spectrum or one tail coefficient.
- Pin the expansion and reheating history, transfer functions, gauge convention, and perturbative validity range.
- State collapse threshold, profile, smoothing window, critical scaling, equation of state, and QCD treatment.
- Map formation abundance to present mass and fraction, including accretion, evaporation, and mergers.
- Pin PTA timing files, clocks, ephemeris, pulsar-noise and correlation model, astrophysical prior, and chains.
- For LVK, pin cross spectra, calibration, spectral template, foreground model, and likelihood revision.
- Compare source models with evidence or predictive checks and expose prior sensitivity.
- Separate replaying chains from reproducing timing residuals and the stochastic likelihood.
Source and data map
Section titled “Source and data map”Domènech supplies the production framework; Pi supplies the current non-Gaussian abundance bridge. Agazie et al. provide direct evidence for the nanohertz background, while Afzal et al. provide the NANOGrav source-model likelihood and explicit caution. Carr et al. provide current PBH synthesis; LIGO-P2500150-v13 supplies the independent ground-band translation of the O1–O4a null; Gouttenoire et al. provide adverse evidence for the distinct PBH-reheating forecast.
NANOGrav 15-year v1.0.1, DOI 10.5281/zenodo.8104459, includes timing files, templates, clocks, and documented timing software. The official NANOGrav data page links PTArcade, chains, and model files. The LIGO-G2502630-v1 release provides the O1–O4a cross spectra and isotropic-search products used by model translations; no separate durable public model-likelihood environment was verified for LIGO-P2500150-v13. No end-to-end rerun was performed. Reproduction must pin timing and noise products, clocks and ephemeris, correlations, astrophysical priors, scalar statistics, thermal history, collapse physics, sampler, and external PBH constraints. IPTA DR3 was not verified as a final public product, and no durable code archive was found for arXiv:2605.21474v1.
References
Section titled “References”- Afzal, A., et al. 2023, “The NANOGrav 15-year Data Set: Search for Signals from New Physics,” Astrophysical Journal Letters 951, L11, DOI 10.3847/2041-8213/acdc91.
- Agazie, G., et al. 2023, “The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background,” Astrophysical Journal Letters 951, L8, DOI 10.3847/2041-8213/acdac6.
- Carr, B., Iovino, A. J., Perna, G., Vaskonen, V., and Veermäe, H. 2026, “Primordial Black Holes: Constraints, Potential Evidence and Prospects,” Rivista del Nuovo Cimento 49, 225–274, arXiv:2601.06024v2, DOI 10.1007/s40766-026-00080-z.
- Domènech, G. 2021, “Scalar Induced Gravitational Waves Review,” Universe 7, 398, DOI 10.3390/universe7110398.
- Gouttenoire, Y., Leister, M., and Schwaller, P. 2026, “Gravitational Waves from Black Hole Reheating: The Scalar-Induced Component,” arXiv:2605.21474v1, preprint, arXiv.
- LIGO Scientific Collaboration, Virgo Collaboration, and KAGRA Collaboration. 2025, “Cosmological and High Energy Physics Implications from Gravitational-Wave Background Searches in LIGO–Virgo–KAGRA’s O1–O4a Runs,” LIGO-P2500150-v13, arXiv:2510.26848v2, official record.
- LIGO Scientific Collaboration, Virgo Collaboration, and KAGRA Collaboration. 2025, “Data for Upper Limits on the Isotropic Gravitational-Wave Background from the First Part of LIGO, Virgo, and KAGRA’s Fourth Observing Run,” LIGO-G2502630-v1, official data release.
- Pi, S. 2025, “Non-Gaussianities and Primordial Black Holes,” in Primordial Black Holes, arXiv:2404.06151v3, DOI 10.1007/978-981-97-8887-3_7.