Pith. sign in

REVIEW 3 major objections 6 minor 1 cited by

Scalar Induced Gravitational Waves signaling Primordial Black Hole Dark Matter

T0 review · 3 major / 6 minor · reviewed 2026-08-16 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper's central claim is that the reported pulsar-timing-array gravitational wave background cannot be explained by primordial black hole dark matter if the curvature power spectrum is monochromatic, because matching the signal would…

desk verdict Short proceedings restating the author's JCAP work; the new finite-width scan is preliminary and the monochromatic-exclusion claim needs a posterior scan over A_G. read the letter →

arxiv 2504.18237 v1 pith:L6LAYAME submitted 2025-04-25 astro-ph.CO gr-qc

Cristian Joana This is my paper · ORCID
classification astro-ph.COgr-qc
keywords primordialblackholesscalar-inducedgravitationalwavesnon-Gaussianitypulsartimingarraysdarkmatterwavebackgroundcurvaturepowerspectrumlognormalpeak
topics Dark Matter
open problems Dark Matter
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper argues that the nHz gravitational-wave background reported by pulsar timing arrays, if interpreted as scalar-induced gravitational waves from primordial black hole (PBH) formation, cannot come from a monochromatic curvature power spectrum. With the spectrum amplitude fixed so that PBHs make up all of the dark matter, the signal amplitude needed to match the pulsar data implies more black holes than the dark matter budget allows. The proposed way out is a power spectrum with finite width, such as a lognormal peak, which suppresses the overproduction while keeping the signal visible. The analysis includes local non-Gaussianity through a logarithmic duality that changes both the collapse threshold and the probability tail, and computes the gravitational-wave spectrum to third order in the scalar amplitude. If the claim is right, finite-width peaks become the viable route for PBH dark matter and the resulting spectra give sharp targets for mHz and nHz observatories.

What carries the argument

The load-bearing machinery is the logarithmic duality, a relation that ties the non-Gaussian parameter $\gamma$ to both the curvature perturbation profile $\zeta(r)$ and its probability distribution function. The profile is written $\zeta(r;\mu,\gamma)=-\frac{1}{\gamma}\ln[1-\gamma\zeta_G(r;\mu)]$ with $\zeta_G(r;\mu)=\mu\,\mathrm{sinc}(kr)$, and the PDF takes a shifted exponential form. This duality carries the argument because it lets the same calculation track how non-Gaussianity changes the collapse threshold $\mu_c(\gamma)$ (obtained from full numerical relativity simulations) and reshapes the rare tail that controls PBH abundance. The method assumes a monochromatic power spectrum $P_\zeta(k)=A_G\delta(k_\star-k)$, fixes $A_G$ by demanding that all dark matter is PBHs, and computes the scalar-induced gravitational wave spectrum perturbatively to order $O(A_S^3)$ with non-Gaussian corrections. The finite-width lognormal peak is the ingredient that evades the overproduction bound.

What would settle it

Recompute the PBH abundance and the scalar-induced gravitational-wave spectrum for a strictly monochromatic peak with $A_G$ chosen so the nHz gravitational-wave amplitude matches the reported signal; if the resulting dark-matter fraction $f_{\rm tot}^{\rm PBH}$ is less than one, the overproduction claim is refuted. An independent measurement that PBHs constitute, say, at most 1% of the dark matter would also lower the normalization enough to test whether the incompatibility survives.

Watch

Extended reading notes

Core claim

The core discovery is an incompatibility: under the standard assumption of a monochromatic curvature power spectrum, the pulsar-timing-array gravitational wave background cannot be consistently explained by PBH dark matter. Fixing the spectrum amplitude so that PBHs account for the entire dark matter abundance ($f_{\rm tot}^{\rm PBH}=1$), the amplitude required to reproduce the detected nHz background implies a PBH abundance that overshoots the dark matter density. The author reports that this tension is removed when the power spectrum has finite width, e.g. a lognormal peak of width $\Delta$, where peak amplitudes around $A_s/\sqrt{2\pi\Delta}\gtrsim 3\times10^{-2}$ can explain the PTA signal without overproducing PBHs, even for nearly Gaussian perturbations ($\gamma\approx f_{\rm NL}\approx 0$). Non-Gaussianity is treated through a logarithmic duality relating the non-Gaussian parameter $\gamma$ to both the curvature profile and the probability density, which shifts the collapse threshold and the abundance. The gravitational-wave spectrum is computed perturbatively including non-Gaussian corrections up to third order in the scalar amplitude.

Load-bearing premise

The argument assumes PBHs make up all of the dark matter ($f_{\rm tot}^{\rm PBH}=1$), which fixes the curvature power-spectrum amplitude; if PBHs are only a fraction of the dark matter, the required amplitude falls and the claimed overproduction may disappear.

Editorial extensions

If this is right

  • A monochromatic curvature power spectrum cannot be the source of the reported nHz gravitational-wave background if PBHs are all of the dark matter, because the required amplitude would overproduce PBHs.
  • A finite-width, lognormal peak in the curvature power spectrum can restore compatibility; for widths $\Delta>0.1$, peak amplitudes around $A_s/\sqrt{2\pi\Delta}\gtrsim 3\times10^{-2}$ can explain the PTA signal even for nearly Gaussian perturbations.
  • Local non-Gaussianity shifts the collapse threshold and the abundance: positive $\gamma$ enhances PBH production, negative $\gamma$ down to about $-3.1$ suppresses it, and strongly negative values enhance it again.
  • The scalar-induced gravitational-wave spectrum includes non-Gaussian corrections up to third order, giving PBH dark matter models concrete predictions in both the mHz band and the nHz band.
  • If the incompatibility holds, future pulsar-timing data can discriminate monochromatic from finite-width sources through the shape and amplitude of the gravitational-wave background.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Inference: The overproduction argument depends as much on the dark matter normalization as on the signal; if independent measurements show PBHs are only a small fraction of the dark matter, the required curvature amplitude drops and a monochromatic peak may fit the nHz data after all.
  • Inference: The finite-width resolution shifts model-building preferences: inflationary scenarios producing broad curvature peaks become natural candidates, and the width of the gravitational-wave bump could be used to infer the width $\Delta$ of the peak.
  • Inference: The logarithmic duality turns the gravitational-wave background into a possible probe of primordial non-Gaussianity on small scales, where CMB measurements cannot reach, if the spectrum shape can be separated from the finite-width degeneracy.
  • Inference: The compatibility of finite-width PBH dark matter with pulsar-timing data is preliminary; a full treatment should check whether the finite-width peak also satisfies CMB spectral-distortion bounds and the many PBH abundance constraints across masses.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. This is a two-page proceedings contribution summarizing recent work by the author and collaborators on primordial black hole (PBH) formation and scalar-induced gravitational waves (SIGWs) in the presence of local non-Gaussianities. The paper describes a methodology that combines numerical-relativity collapse thresholds, the logarithmic duality linking the non-Gaussian parameter γ to the curvature profile and PDF, and a monochromatic curvature power spectrum normalized so that PBHs constitute all dark matter (f_tot_PBH = 1). The central claim, stated in Section 3, is that the NANOGrav signal is incompatible with monochromatic PBH formation because matching the detected amplitude would overproduce PBHs, and that preliminary finite-width (lognormal) peaks may alleviate this tension. The abstract highlights consequences for LISA and PTA experiments.

Significance. If correct, the incompatibility claim would be a strong and useful constraint on monochromatic PBH dark matter models, redirecting attention to finite-width spectra and providing a concrete, falsifiable target for LISA and PTA. The framework itself, based on the logarithmic duality and full numerical relativity thresholds, is sophisticated and goes beyond simpler Press-Schechter estimates. The paper also gives credit to machine-checkable or reproducible elements by referencing the companion work [1] and the peak-theory method [4]. However, the present manuscript is a short proceedings summary: the central negative claim is asserted rather than derived in the text, and the finite-width resolution is explicitly labeled preliminary. As a standalone paper, the evidence presented is insufficient to fully support its conclusions; its significance would be substantially higher if the derivations or posterior scans were included.

major comments (3)
  1. [Section 3, Fig. 1 (left)] The statement that 'the PTA signal is incompatible with PBH formation scenarios under the assumption of a monochromatic power spectrum, as the detected amplitude would imply an overproduction of PBHs' is not demonstrated in this manuscript. Section 2 fixes A_G by requiring f_tot_PBH = 1, but the incompatibility claim requires scanning the NANOGrav posterior over the SIGW amplitude, inverting the monochromatic SIGW prediction to obtain A_G for each allowed amplitude, and then computing f_PBH(A_G) with the same thresholds and PDF. No such scan is reported here, and the left panel of Fig. 1 is only described as 'adapted from Ref. [1]' without showing axes, posterior contours, or the quantitative relation between A_G and the PTA amplitude. This missing step is load-bearing because the overproduction argument depends on the actual A_G inferred from the detected signal, not on the A_G chosen by the f=1 normalization.
  2. [Section 3, right panel and Section 2] The finite-width resolution is explicitly preliminary, but it is still used to draw the paper's main positive conclusion. The text states 'For widths larger than Δ > 0.1 we find peak amplitudes above As/√(2πΔ) ≳ 3·10^-2 which can explain the reported PTA signal', but it does not define the normalization convention for the lognormal peak, the value of A_s, the details of the peak-theory computation, the assumed value of γ (given as γ ≈ 0), or any error bars. Without this information, a reader cannot check whether the finite-width mechanism indeed resolves the overproduction tension or whether the quoted amplitude threshold is a fit to the data rather than a prediction. The preliminary nature of this part should either be acknowledged in the conclusions as an outlook or the computation must be specified in enough detail to be reproducible.
  3. [Section 2, bullet list] The normalization 'The amplitude A_G is fixed by that all dark matter in the form of PBHs (i.e. f_tot_PBH = 1)' is an assumption, not a prediction, and it is structurally load-bearing for the central claim in Section 3. If PBHs are subdominant, a smaller A_G would produce fewer PBHs and a weaker SIGW background, so the claimed overproduction and the resulting incompatibility with the NANOGrav signal would not necessarily hold. The manuscript should state explicitly in the conclusions that the incompatibility is conditional on PBHs being all of the dark matter, and it should discuss the regime f_PBH < 1, where the constraint does not follow.
minor comments (6)
  1. [Title] The title contains a typo: 'W aves' should be 'Waves'.
  2. [Eq. (1)] Equation (1) is typeset in a garbled way; the fraction and exponent are not readable in the provided text. Please ensure the PDF is written with clear notation, e.g., P(ζ) = (1/√(2πσ^2)) exp[- (e^{-γζ} - 1)^2/(2γ^2σ^2) - γζ].
  3. [Section 2, bullet list] The phrase 'The amplitude AG is fixed by that all dark matter' is grammatically incomplete; it should read 'The amplitude A_G is fixed by requiring that all dark matter be in the form of PBHs (i.e., f_tot_PBH = 1).'
  4. [Section 3, right panel] The expression 'As/√(2πΔ)' appears to be a typo: for a lognormal power spectrum, the peak value is typically A_s/(√(2π)Δ), not A_s/√(2πΔ). Please define A_s and the normalization convention explicitly.
  5. [Section 2, text near Eq. (2)] The notation 'γ ≈ fNL ≈ 0' identifies γ with the conventional non-Gaussianity parameter f_NL, but γ is introduced as the logarithmic-duality parameter. Please clarify the relationship between γ and f_NL or avoid equating them without definition.
  6. [Fig. 1] Figure 1 lacks axis labels and a legend in the displayed version; the left panel's posterior contours cannot be interpreted without knowing what parameters are plotted (presumably A_G and γ) and the corresponding confidence levels.

Circularity Check

0 steps flagged · score 0.0 of 10

No circular reduction found: the f_PBH=1 normalization is an explicit input assumption, and the PTA statements are conditional constraints rather than re-labeled fitted parameters.

full rationale

The proceedings is a summary of Ref [1] and does not attempt to derive its inputs from its conclusions. A_G is explicitly fixed by assuming PBHs constitute all dark matter (f_tot_PBH=1), and the subsequent SIGW spectrum is a computed output of that assumption. The statement that monochromatic PBH scenarios are incompatible with the PTA signal is a constraint obtained by comparing that computed SIGW amplitude with the NANOGrav data; it is not equivalent by construction to the input f_PBH=1, because the comparison could in principle yield either overproduction or underproduction depending on the numerical values. The finite-width panel is explicitly labeled 'preliminary' and merely scans Δ>0.1 to show that the tension can be alleviated; it is not presented as a parameter-free prediction. The cited methods [2-4] are prior external computations; although several are by overlapping authors, the proceedings does not invoke an unverified uniqueness theorem or define its quantities in terms of its conclusions. The notable gap that the 'detected amplitude would imply an overproduction' sentence is not backed by a displayed full-posterior inversion is a support/completeness concern, not a circularity under the hard rules.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The central claim rests on a chain of assumptions and model choices: the logarithmic duality, third-order SIGW computation, the f_PBH=1 normalization, peak theory, and the interpretation of the NANOGrav signal. No new physical entities are introduced.

free parameters (2)
  • gamma (non-Gaussianity parameter) = not fitted, scanned
    Parameter of the local non-Gaussian model; the paper reports how PBH abundance and SIGW spectrum depend on gamma, with strong suppression for 0 > gamma > -3.1 and enhancement for gamma < -3.1. It is an input to the model, not determined by the analysis.
  • delta (lognormal peak width) = greater than 0.1
    Width of the lognormal power spectrum used in the preliminary finite-width analysis; the text states that for widths larger than delta > 0.1 the peak amplitudes can explain the PTA signal, so delta is scanned to match the data.
assumptions (5)
  • domain assumption The logarithmic duality of Ref [2] correctly links the non-Gaussian parameter gamma to the curvature profile and PDF.
    Equations (1) and (2) adopt this relation from Pi and Sasaki, with authorship overlapping the present work; its validity is assumed.
  • domain assumption Scalar-induced gravitational waves can be computed perturbatively to third order in A_S using the formalism of Ref [3].
    The spectrum computation is taken from Abe et al., which shares authors with the present paper.
  • domain assumption Primordial black holes make up all of the dark matter, f_tot_PBH = 1, fixing the amplitude A_G.
    Section 2 states that A_G is fixed by this condition; the overproduction argument and the comparison with the PTA signal depend on this assumption.
  • domain assumption The peak theory methodology of Ref [4] correctly predicts the PBH mass function for finite-width spectra.
    The right panel of Fig 1 relies on this method from Pi et al., with overlapping authorship; no independent validation is shown in this paper.
  • domain assumption The NANOGrav 15-year signal is a stochastic gravitational wave background suitable for comparison with SIGW predictions.
    Ref [5] is used as the observational data set; the paper does not model possible astrophysical foregrounds or alternative origins.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Scalar Induced Gravitational Waves signaling Primordial Black Hole Dark Matter." pith.science (2026). https://pith.science/paper/L6LAYAME

@misc{pith2026250418237,
  author       = {Pith},
  title        = {Pith review of: Scalar Induced Gravitational Waves signaling Primordial Black Hole Dark Matter},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/L6LAYAME}},
  note         = {Machine review of arXiv:2504.18237}
}
read the original abstract

Primordial black holes are a unique probe of the early Universe and offer a potential link between inflationary dynamics and dark matter. In a recent work \cite{Inui:2024fgk}, we analyse PBH formation and the associated stochastic gravitational wave background from scalar induced gravitational waves, while taking into account the contributions from local non-Gaussianities. We highlight the observable consequences for the LISA and PTA experiments.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. Primordial black holes forming during kination: the trapped, the overdense, and the void

    astro-ph.CO 2026-07 conditional novelty 6.0 of 10

    In one non-attractor inflation model, initially similar field fluctuations produce three distinct black-hole formation channels—trapped, overdense, void—with collapse thresholds determined by the full density profile,...

Reference graph

Works this paper leans on

9 extracted references · 3 canonical work pages · cited by 1 Pith paper

  1. [1]

    R. Inui, C. Joana, H. Motohashi, S. Pi, Y. Tada and S. Yokoyama, JCAP 03 (2025), 021

    work page 2025

  2. [4]

    S. Pi, M. Sasaki, V. Takhistov and J. Wang, [arXiv:2501.00295 [astro-ph.CO]]

    arXiv

  3. [2]

    Pi and M

    S. Pi and M. Sasaki, Phys. Rev. Lett. 131, no.1, 011002 (2023)

    2023

  4. [3]

    K. T. Abe, R. Inui, Y. Tada and S. Yokoyama, JCAP 05, 044 (2023)

    work page 2023

  5. [5]

    Afzal et al

    A. Afzal et al. [NANOGrav], Astrophys. J. Lett. 951, no.1, L11 (2023)

    work page 2023

  6. [6]

    C Jarlskog (World Scientific, Singapore, 1988)

    C Jarlskog in CP Violation , ed. C Jarlskog (World Scientific, Singapore, 1988)

    1988

  7. [7]

    Maiani, 62 183 1976

    L. Maiani, 62 183 1976

    1976

  8. [8]

    Bjorken and I

    J.D. Bjorken and I. Dunietz, 36 2109 1987

    1987

Show all 9 references
  1. [9]

    Buchanan et al , 45 4088 1992

    C.D. Buchanan et al , 45 4088 1992

    1992

Pith tools

Reviewed August 16, 2026 · model on record in the stance chip above.