Pith. sign in

REVIEW 1 major objections 3 minor 2 cited by

Primordial Black Holes (PBHs) and The Signatures of Cosmic Non-Gaussianity

T0 review · 1 major / 3 minor · reviewed 2026-05-18 · grok-4.3

Challenge this review Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The exact sudden-decay relation in the curvaton scenario yields the full non-Gaussian probability distribution of curvature perturbations from which primordial black hole abundances are computed directly.

desk verdict The paper works out the full non-Gaussian PDF for the curvaton curvature perturbation via explicit branchwise change of variables from the sudden-decay map, then runs PBH and induced-GW calculations from one self-consistent spectrum instead of per-scale fitting. read the letter →

arxiv 2509.10851 v2 submitted 2025-09-13 astro-ph.CO

Owais Farooq , Romana Zahoor , Balungi Francis This is my paper · ORCID
classification astro-ph.CO
keywords primordialblackholescurvatonscenarionon-Gaussianitycurvatureperturbationsuddendecayprobabilitydensityfunctioninducedgravitationalwaves
checked against Cost.FunctionalEquation
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 establishes that primordial black hole formation depends exponentially on the far tail of the curvature perturbation distribution, making it a sensitive test of non-Gaussianity. The authors derive the curvature perturbation exactly from the sudden-decay relation in the curvaton model and obtain the complete probability density function by performing an explicit branchwise change of variables on the underlying Gaussian curvaton field fluctuations. This non-perturbative method is then used to evaluate the black hole formation fraction without Edgeworth truncations or other approximations. The same consistent fluctuation spectrum also determines the induced gravitational wave background, allowing direct comparison with current constraints and future detector sensitivities.

What carries the argument

The exact sudden-decay relation between the curvaton field value and the curvature perturbation, which enables an explicit branchwise change of variables that maps the Gaussian field distribution onto the complete non-Gaussian curvature perturbation probability density function.

What would settle it

A numerical simulation or observation showing that the actual mapping from curvaton field to curvature perturbation deviates substantially from the sudden-decay formula on the relevant scales would invalidate the exact tail used for the black-hole fraction.

Watch

Extended reading notes

Core claim

In the curvaton scenario the curvature perturbation is derived from the exact sudden-decay relation, the full probability density function is obtained through an explicit branchwise change of variables from the Gaussian curvaton-field fluctuation, and the primordial black-hole formation fraction is evaluated from the exact non-perturbative tail; the same model supplies a self-consistent fluctuation spectrum that replaces scale-by-scale variance matching and generates the induced gravitational-wave background.

Load-bearing premise

The curvaton decays suddenly so that the curvature perturbation is given directly by the field value before decay without any gradual transition or smoothing corrections.

Editorial extensions

If this is right

  • The primordial black hole mass function follows directly from the exact non-perturbative tail without perturbative truncation or additional fitting.
  • Induced gravitational wave spectra are generated from the linear curvaton two-point function within the same self-consistent model.
  • Predictions differ measurably from both the pure Gaussian case and the exact local quadratic benchmark.
  • Current conservative constraints on primordial black holes can be confronted with mass functions computed once for the entire spectrum.

Reading between the lines

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

  • The branchwise method could be extended to other sources of non-Gaussianity such as multi-field inflation to obtain similarly exact tails.
  • Precision measurements of the stochastic gravitational wave background by future detectors could distinguish this curvaton signature from other early-universe scenarios.
  • The exact tail calculation may tighten or relax existing limits on curvaton parameters once applied to updated observational envelopes.
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, simulated authors' rebuttal, and a circularity audit.

Referee Report

1 major / 3 minor

Summary. The manuscript studies the curvaton scenario for primordial black hole formation by deriving the curvature perturbation from the exact sudden-decay relation. It obtains the full probability density function of the curvature perturbation through a branchwise change of variables from Gaussian field fluctuations and evaluates the PBH formation fraction from the exact non-perturbative tail. The results are compared to Gaussian and local quadratic non-Gaussian cases. A self-consistent model is used where the dimensionless field fluctuation spectrum is specified once, allowing direct computation of smoothed variances and collapse fractions without additional per-scale fitting. The induced gravitational wave background is generated from the linear spectrum, and the mass functions are compared to constraints with GW spectra shown against observational windows.

Significance. If the central derivations hold, the work provides a non-perturbative analytic treatment of non-Gaussianity in the curvaton model for PBH formation, valuable due to the exponential sensitivity of PBH abundance to the tail. The explicit step-by-step analytic forms for the Jacobian, normalization, support, and small-fluctuation expansion, together with the self-consistent fluctuation model that avoids per-scale variance matching, are strengths. Confrontation with current constraints and display of induced GW spectra against PTA/LISA/DECIGO windows add observational relevance.

major comments (1)
  1. [Derivation of the curvature perturbation from the exact sudden-decay relation] The central construction obtains the PDF of ζ by an explicit change of variables using the sudden-decay mapping ζ = (2/3) ln(1 + (3/2)σ/σ_d) (or equivalent) inserted directly before the Jacobian and support calculation. This assumes instantaneous decay. Realistic models with finite decay rate Γ introduce smoothing from integrated decay dynamics that suppresses the extreme positive tail controlling the PBH collapse fraction. This premise is load-bearing for the claim of an 'exact non-perturbative tail' and the resulting PBH mass functions.
minor comments (3)
  1. [Self-consistent curvaton fluctuation model] The self-consistent model reduces per-scale fitting but still treats the dimensionless field fluctuation spectrum as a free parameter whose value is chosen once; clarify how this choice is motivated or varied to test robustness of the final mass functions and GW spectra.
  2. [Probability density function derivation] Ensure that the support of the distribution and the branchwise handling in the change-of-variables step are illustrated with an explicit equation or figure for the reader to follow the normalization.
  3. [Results and figures] The figures displaying mass functions against the conservative constraint envelope and GW spectra against sensitivity windows are useful; verify that all curves are labeled with the corresponding model (exact curvaton, Gaussian benchmark, quadratic) for clarity.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for the careful reading and constructive comments. We address the major comment below and indicate the revisions we will implement.

read point-by-point responses

  1. Referee: The central construction obtains the PDF of ζ by an explicit change of variables using the sudden-decay mapping ζ = (2/3) ln(1 + (3/2)σ/σ_d) (or equivalent) inserted directly before the Jacobian and support calculation. This assumes instantaneous decay. Realistic models with finite decay rate Γ introduce smoothing from integrated decay dynamics that suppresses the extreme positive tail controlling the PBH collapse fraction. This premise is load-bearing for the claim of an 'exact non-perturbative tail' and the resulting PBH mass functions.

    Authors: We agree that the sudden-decay approximation is central to our analytic derivation and that a finite decay width would introduce smoothing that could suppress the far positive tail. Our manuscript explicitly adopts the sudden-decay relation as a standard analytic limit in the curvaton literature, allowing the exact branchwise change-of-variables PDF to be derived without truncation. We will revise the text to add a dedicated paragraph in Section 2 (or the discussion) that (i) states the approximation explicitly, (ii) cites literature on finite-Γ effects, and (iii) notes that the reported tail represents the limiting case of maximal non-Gaussian enhancement. This clarifies the scope without altering the central derivations. revision: yes

standing simulated objections not resolved
  • Quantitative evaluation of the tail suppression under a finite decay rate Γ, which would require a separate integration over the decay dynamics outside the present sudden-decay framework.

Circularity Check

0 steps flagged · score 0.0 of 10

Derivation proceeds via explicit change-of-variables from stated sudden-decay relation; self-consistent spectrum replaces per-scale fitting without reducing to input by construction

full rationale

The central steps derive the full PDF of ζ from Gaussian σ fluctuations by applying the exact sudden-decay mapping ζ = (2/3) ln(1 + (3/2)σ/σ_d) followed by Jacobian and support calculation. This is a direct mathematical transformation, not a redefinition or fit. The self-consistent curvaton model specifies the dimensionless field fluctuation spectrum once, computes smoothed variances directly, and obtains the collapse fraction with no further per-scale fitting; the resulting mass functions and induced GW spectra are generated from this single pipeline and confronted with external constraints. No self-citations, uniqueness theorems, or ansatzes smuggled via prior work are invoked as load-bearing. The sudden-decay premise is an explicit model assumption whose consequences are computed rather than presupposed, so the derivation chain remains self-contained against the stated benchmarks.

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

The central results rest on the standard assumption that the curvaton field fluctuation is Gaussian at horizon exit and on the sudden-decay approximation; no new free parameters beyond the single spectrum choice are introduced, and no new entities are postulated.

free parameters (1)
  • dimensionless field fluctuation spectrum
    Specified once in the self-consistent model to determine smoothed variances on all scales without further per-scale adjustment.
assumptions (2)
  • domain assumption Curvaton field fluctuations are Gaussian
    Starting distribution for the branchwise change of variables to obtain the curvature perturbation PDF.
  • domain assumption Sudden decay approximation holds
    Used to write the exact relation between curvaton field value and final curvature perturbation before the variable change.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Primordial Black Holes (PBHs) and The Signatures of Cosmic Non-Gaussianity." pith.science (2026). https://pith.science/paper/2509.10851

@misc{pith2026250910851,
  author       = {Pith},
  title        = {Pith review of: Primordial Black Holes (PBHs) and The Signatures of Cosmic Non-Gaussianity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2509.10851}},
  note         = {Machine review of arXiv:2509.10851}
}
read the original abstract

Primordial black-hole formation depends exponentially on the far tail of the primordial curvature-perturbation distribution. That sensitivity makes the small-scale collapse problem a sharp probe of primordial non-Gaussianity. We study the curvaton scenario by deriving the curvature perturbation from the exact sudden-decay relation, obtaining the full probability density function through an explicit branchwise change of variables from the Gaussian curvaton-field fluctuation, and evaluating the primordial black-hole formation fraction from the exact non-perturbative tail. The derivation is written step by step, with the support of the distribution, the Jacobian, the normalization, and the small-fluctuation expansion displayed in analytic form. We place the exact curvaton prediction beside the Gaussian benchmark and beside an exact local quadratic benchmark in which the non-Gaussian probability density is also computed without an Edgeworth truncation. We then replace the scale-by-scale variance-matching ansatz by a self-consistent curvaton fluctuation model in which the dimensionless field fluctuation spectrum is specified once, the smoothed curvaton variance is computed directly, the exact collapse fraction follows with no further fitting on each scale, and the induced gravitational-wave background is generated from the linear curvaton two-point spectrum implied by the same model. The resulting mass functions are confronted with a conservative current constraint envelope motivated by recent primordial-black-hole reviews, and the induced gravitational-wave spectra are displayed against the PTA, LISA, and DECIGO sensitivity windows. The final figures are generated from a single mathematically consistent numerical pipeline.

Figures

Figures reproduced from arXiv: 2509.10851 by the authors.

Figure 1
Figure 1. Primordial black hole (PBH) abundance β(M) as a function of the variance σ. The Gaussian case (black) shows the standard exponential suppression. Positive local non￾Gaussianity (fNL = +5, blue) significantly enhances PBH production, while negative fNL (fNL = −5, orange) suppresses it. This demonstrates the exponential sensitivity of PBH formation to the non-Gaussian tail of the probability distribution [PITH_FULL_I… view at source ↗
Figure 2
Figure 2. Comparison of the probability distribution functions (PDFs) of the curvature pertur [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Exact non-Gaussian PDFs of the curvature perturbation [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: PBH abundance β(M) as a function of the curvaton decay fraction Ωχ,dec. Small values of Ωχ,dec generate large non-Gaussianities and strongly enhance PBH production, whereas in the limit Ωχ,dec → 1 the abundance approaches the Gaussian prediction. 8 [PITH_FULL_IMAGE:fi…
Figure 5
Figure 5. Figure 5: compares β computed in three ways: (i) the Gaussian approximation via Eq. (3) using a target smoothed variance σ 2 , (ii) the perturbative local-fNL correction (Edge￾worth/leading cumulant), and (iii) the exact curvaton PDF Eq. (26) using representative curvaton parame…
Figure 6
Figure 6. Figure 6: Comparison of PBH formation in Gaussian and curvaton scenarios. The [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]
Figure 7
Figure 7. Figure 7: Representative PBH mass spectra β(M), showing illustrative mass fractions computed for different curvature perturbation amplitudes A = 0.01, 0.05, 0.2. The calculation assumes a lognormal power spectrum. The plot demonstrates that increasing the amplitude enhances the …
Figure 8
Figure 8. Figure 8: induced gravitational-wave spectra corresponding to different primordial black hole [PITH_FULL_IMAGE:figures/full_fig_p012_8.png]
Figure 9
Figure 9. Figure 9: induced gravitational-wave spectra ΩGW(f) for PBH mass scales corresponding to PTA, LISA and DECIGO bands. Curve amplitudes were scaled as an illustrative proxy of β(M). 6 Results and Conclusion In this work we have carried out a detailed study of primordial black hole…

Discussion (0). Continue with ORCID to comment.

Lean theorems connected to this paper

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

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    Solving Eq. (15) gives e^{3ζ_χ} = Y(ζ) ≡ 1/Ω (e^{3ζ} + (Ω-1)e^{-ζ}); combining with e^{3ζ_χ} = (1 + δχ/χ̄)^2 yields the mapping δχ(ζ) = χ̄ (±√(Y(ζ)−1)) and the change-of-variables PDF P_ζ(ζ) = Σ_s P_χ(δχ_s(ζ)) |dδχ_s/dζ|.

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 2 Pith papers

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

  1. Long-term 3+1 simulations of primordial black hole formation during radiation domination

    gr-qc 2026-08 conditional novelty 7.0 of 10

    A new time-step gauge makes long 3D primordial black hole formation simulations about 94 times cheaper and reproduces the known spherical collapse threshold and critical exponent.

  2. Numerical simulations of primordial black hole formation via delayed first-order phase transitions

    gr-qc 2026-01 conditional novelty 6.0 of 10

    Spherically symmetric numerical relativity shows false-vacuum domains from delayed first-order phase transitions form type B (baby-universe) or type A (direct-collapse) primordial black holes, separated by a robust t_...

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages · cited by 2 Pith papers

  1. [1]

    K. Ando, M. Kawasaki, and H. Nakatsuka. Formation of primordial black holes in an axion-like curvaton model.Phys. Rev. D, 98(8):083508, 2018

    work page 2018

  2. [2]

    Bartolo, S

    N. Bartolo, S. Matarrese, and A. Riotto. Non-gaussianity from inflation.Phys. Rev. D, 65:103505, 2002

    work page 2002

  3. [3]

    Bean and J

    R. Bean and J. Magueijo. Could supermassive black holes be quintessential primor- dial black holes?Phys. Rev. D, 66:063505, 2002

    work page 2002

  4. [4]

    J. S. Bullock and J. R. Primack. Non-gaussian fluctuations and primordial black holes from inflation.Phys. Rev. D, 55:7423–7439, 1997

    work page 1997

  5. [5]

    C. T. Byrnes, E. J. Copeland, and A. M. Green. Primordial black holes as a tool for constraining non-gaussianity.Phys. Rev. D, 86:043512, 2012

    work page 2012

  6. [6]

    B. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama. Constraints on primordial black holes.Rept. Prog. Phys., 84(11):116902, 2021

    work page 2021

  7. [7]

    B. Carr, F. Kuhnel, and M. Sandstad. Primordial black holes as dark matter.Phys. Rev. D, 94(8):083504, 2016

    work page 2016

  8. [8]

    B. J. Carr and S. W. Hawking. Black holes in the early universe.Mon. Not. Roy. Astron. Soc., 168:399–415, 1974

    work page 1974

Show all 37 references
  1. [9]

    X. Chen. Primordial non-gaussianities from inflation models.Adv. Astron., 2010:638979, 2010

    2010

  2. [10]

    L. S. Collaboration and V. Collaboration. Observation of gravitational waves from a binary black hole merger.Phys. Rev. Lett., 116(6):061102, 2016

    2016

  3. [11]

    Collaboration

    N. Collaboration. The nanograv 15-year data set: Evidence for a gravitational-wave background.Astrophys. J. Lett., 951(1):L8, 2023

    2023

  4. [12]

    Collaboration

    P. Collaboration. Planck 2018 results. x. constraints on inflation.Astron. Astrophys., 641:A10, 2020

    2018

  5. [13]

    Enqvist, A

    K. Enqvist, A. Jokinen, S. Kasuya, and A. Mazumdar. Non-gaussianity from pre- heating.Phys. Rev. Lett., 94:161301, 2005

    2005

  6. [14]

    Enqvist, S

    K. Enqvist, S. Nurmi, O. Taanila, and T. Takahashi. Non-gaussian fingerprints of self-interacting curvaton.JCAP, 11:034, 2013

    2013

  7. [15]

    Garcia-Bellido and E

    J. Garcia-Bellido and E. Ruiz Morales. Primordial black holes from single field models of inflation.Phys. Dark Univ., 18:47–54, 2017

    2017

  8. [16]

    A. M. Green and B. J. Kavanagh. Primordial black holes as a dark matter candidate. J. Phys. G, 48(4):043001, 2021

    2021

  9. [17]

    A. M. Green and A. R. Liddle. Constraints on the density perturbation spectrum from primordial black holes.Phys. Rev. D, 56:6166–6174, 1997

    1997

  10. [18]

    Harada, C

    T. Harada, C. M. Yoo, and K. Kohri. Threshold of primordial black hole formation. Phys. Rev. D, 88(8):084051, 2013. 14

    2013

  11. [19]

    Harada, C

    T. Harada, C. M. Yoo, T. Nakama, and Y. Koga. Cosmological long-wavelength solutions and primordial black hole formation.Phys. Rev. D, 91(8):084057, 2015

    2015

  12. [20]

    P. Ivanov. Nonlinear metric perturbations and production of primordial black holes. Phys. Rev. D, 57:7145–7154, 1998

    1998

  13. [21]

    A. S. Josan, A. M. Green, and K. A. Malik. Generalised constraints on the curvature perturbation from primordial black holes.Phys. Rev. D, 79:103520, 2009

    2009

  14. [22]

    Kawasaki and Y

    M. Kawasaki and Y. Tada. Can primordial black holes as dark matter and second order gravitational waves be probed by pulsar timing array?JCAP, 08:041, 2013

    2013

  15. [23]

    Kohri, C.-M

    K. Kohri, C.-M. Lin, and T. Matsuda. Primordial black holes from the inflating curvaton. 2012

    2012

  16. [24]

    Komatsu and D

    E. Komatsu and D. N. Spergel. Acoustic signatures in the primary microwave back- ground bispectrum.Phys. Rev. D, 63:063002, 2001

    2001

  17. [25]

    LoVerde, A

    M. LoVerde, A. Miller, S. Shandera, and L. Verde. Effects of scale-dependent non- gaussianity on cosmological structures.JCAP, 04:014, 2008

    2008

  18. [26]

    D. H. Lyth and Y. Rodriguez. The inflationary prediction for primordial non- gaussianity.Phys. Rev. Lett., 95:121302, 2005

    2005

  19. [27]

    D. H. Lyth, C. Ungarelli, and D. Wands. The primordial density perturbation in the curvaton scenario.Phys. Rev. D, 67:023503, 2003

    2003

  20. [28]

    Moroi and T

    T. Moroi and T. Takahashi. Effects of the curvaton on the spectrum of cosmological perturbations.Phys. Lett. B, 522:215–221, 2001

    2001

  21. [29]

    Musco, J

    I. Musco, J. C. Miller, and A. G. Polnarev. Primordial black hole formation in the radiative era: Investigation of the critical nature of the collapse.Class. Quant. Grav., 26:235001, 2009

    2009

  22. [30]

    Musco, J

    I. Musco, J. C. Miller, and L. Rezzolla. Computations of primordial black hole formation.Class. Quant. Grav., 22:1405–1424, 2005

    2005

  23. [31]

    J. C. Niemeyer and K. Jedamzik. Dynamics of primordial black hole formation. Phys. Rev. D, 59:124013, 1999

    1999

  24. [32]

    Sasaki, J

    M. Sasaki, J. Valiviita, and D. Wands. Non-gaussianity of the primordial perturba- tion in the curvaton model.Phys. Rev. D, 74:103003, 2006

    2006

  25. [33]

    Seery and J

    D. Seery and J. E. Lidsey. Primordial non-gaussianities in single field inflation. JCAP, 06:003, 2005

    2005

  26. [34]

    Shandera, A

    S. Shandera, A. L. Erickcek, P. Scott, and J. Y. Galarza. Number counts and non- gaussianity. 2012

    2012

  27. [35]

    Shibata and M

    M. Shibata and M. Sasaki. Black hole formation in the friedmann universe: Formu- lation and computation in numerical relativity.Phys. Rev. D, 60:084002, 1999

    1999

  28. [36]

    Young and C

    S. Young and C. T. Byrnes. Primordial black holes in non-gaussian regimes.JCAP, 08:052, 2013

    2013

  29. [37]

    Young and C

    S. Young and C. T. Byrnes. Signatures of non-gaussianity in the isocurvature modes of primordial black holes.JCAP, 04:034, 2015. 15

    2015

Pith tools

Reviewed May 18, 2026 · model on record in the stance chip above.