REVIEW 3 major objections 5 minor 15 references
Cosmic Whispers of the Early Universe: Gravitational Waves and Dark Matter from Primordial Black Holes
T0 review · 3 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read Broadest profiles set the critical compaction threshold at 2/5
desk verdict A systematic compilation of PBH abundance methods; the 'always 2/5' threshold for broad spectra is an extrapolation that needs testing, not a settled result. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the compaction function $C(r)$, twice the local mass excess over the areal radius, which on superhorizon scales reads $C(r)=-2\Phi r\,\zeta'(r)[1+\tfrac{r}{2}\zeta'(r)]$. Collapse occurs when its maximum exceeds a threshold $C_{\rm th}$ that depends on the peak curvature $q=-\tfrac{r_m^2 C''(r_m)}{4C(r_m)}$, with $C_{\rm th}(q)$ running from $2/3$ for spiky profiles to $2/5$ for broad ones. The abundance calculation is carried by the joint Gaussian probability of the compaction and the curvature perturbation, with a correlation parameter $\gamma$ that is near unity for peaked spectra and small for broad spectra. Minimizing the effective squared threshold in $q$ pushes the dominant contribution to $q=0$ once $\gamma$ falls below a critical value, which is the mechanism that makes $2/5$ universal. The volume-averaged compaction, $\bar{C}(R_m)=\frac{3}{R_m^3}\int_0^{R_m} dx\,x^2 C(x)$, has threshold $2/5$ for both spiky and broad profiles, providing the profile-independent extension.
What would settle it
Numerical relativity simulations of collapse starting from a set of broad, non-peaked initial profiles with the same compaction maximum but different peak curvatures should find the critical compaction threshold equal to $2/5$ in every case; a spread in thresholds, or values clearly different from $2/5$, would falsify the universal-threshold claim.
Extended reading notes
Core claim
The paper's central claim is that the relevant critical threshold for primordial black hole formation from critical collapse is the one corresponding to the broadest compaction profiles whenever the curvature power spectrum is not sharply peaked. Concretely, the threshold function $C_{\rm th}(q)$, which grows from $2/5$ to $2/3$ as the curvature at the compaction peak increases, enters the abundance integral through an effective threshold; for broad spectra, parametrized by a small correlation parameter $\gamma$, the effective threshold is minimized at $q=0$, forcing $C_{\rm th}=2/5$. The same conclusion holds in the non-Gaussian case for the logarithmic USR-type relation. For very peaked spectra, by contrast, the average profile remains the right choice. The thesis also shows that the volume average of the compaction function has a profile-independent threshold of $2/5$, a step toward an observable with no shape dependence.
Load-bearing premise
The universal $2/5$ threshold holds if the numerically fitted threshold function $C_{\rm th}(q)$ is accurate for every profile shape and if the minimum of the effective threshold really sits at $q=0$ for broad, non-Gaussian spectra; if either fails, the universal value does not follow.
Editorial extensions
If this is right
- For any non-peaked power spectrum, PBH abundance should be computed with the $2/5$ threshold rather than the average-profile threshold; this typically lowers the predicted abundance at a fixed amplitude.
- Peaked spectra remain governed by the average profile, so the universal threshold applies only once the spectrum is sufficiently broad.
- The result extends to the non-Gaussian case for positive, log-type non-Gaussianity, so the $2/5$ rule survives beyond Gaussian statistics.
- Because thresholds enter abundances exponentially, small threshold changes translate into large abundance changes, and retuning the power-spectrum amplitude to compensate changes other predictions such as the induced gravitational-wave signal.
- A profile-independent observable, the volume average of the compaction, can be used with threshold $2/5$; its statistics reduce to computing connected cumulants of the volume-averaged compaction.
Reading between the lines
- If the $2/5$ threshold is universal for broad spectra, published PBH constraints that were derived with average-profile thresholds may need to be re-derived; the direction and size of the shift will depend on each model's spectrum width.
- The volume-averaged compaction route suggests a practical test: compute the one-point statistics of $\bar C(R_m)$ directly from the curvature power spectrum and compare PBH abundances against threshold-statistics results; numerical simulations can then settle whether $2/5$ is truly independent of profile.
- Because the broadest profiles are rare, the claim implies that PBH formation selects atypical fluctuations; if confirmed, this would also affect predictions for the clustering and spin distribution of PBHs, since the collapsing profiles are not the average ones.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This PhD thesis develops a comprehensive framework for computing the primordial black hole (PBH) abundance from inflationary curvature perturbations, with emphasis on primordial non-Gaussianities, threshold statistics on the compaction function, and the dependence of the collapse threshold on the compaction profile. It derives a master formula for the abundance with arbitrary local non-Gaussianity (Eq. 3.58), generalizes it to the PBH mass function (Eq. 3.68), and proposes in Sec. 3.3 that for non-peaked power spectra the abundance is dominated by the broadest compaction profiles, with a claimed universal threshold C_th = 2/5. The thesis then applies these tools to update LVK O3 and PTA constraints, discusses scalar-induced gravitational waves, and analyzes single-field ultra slow-roll and curvaton production models. The abstract and conclusions present PBHs as potential dark matter candidates and as a possible source of the PTA signal.
Significance. If the 'always 2/5' claim survives scrutiny, it is an important result: abundance predictions for broad power spectra would shift by orders of magnitude relative to average-profile prescriptions, with direct consequences for PBH-dark-matter and PTA interpretations. The manuscript is strong in presenting exact resummation-based formulas, in carefully assessing the convergence of perturbative expansions, and in spelling out applied constraints from LVK and PTA data. It also gives credit to the underlying numerical relativity inputs. However, the universal-threshold claim rests on an empirical fit to C_th(q) that has not been validated for the non-Gaussian profiles used in Sec. 3.3.3, and several phenomenological 'predictions' are obtained only after tuning the amplitude A to a target abundance. The core idea is promising, but at present the evidence is not yet at the level required for the strongest version of the claim.
major comments (3)
- [Sec. 3.2.2, Eq. (3.58); Sec. 3.2.3, Eq. (3.68)] The central claim that the threshold is 'always 2/5' for non-peaked power spectra is an extrapolation beyond the calibrated domain of Eq. (3.43). The numerical fit C_th(q) was calibrated on a restricted family of compaction profiles in Refs. [263,265]; the non-Gaussian map ζ = F(ζ_G) changes the relation between C_G, q, and C (Eqs. 3.49 and 3.104), so a profile with q → 0 in the non-Gaussian field is not guaranteed to have the same threshold as the q → 0 profile used to extract 2/5. The manuscript itself notes in Sec. 3.2.1 (footnote 6) that primordial non-Gaussianity can shift δ_c by a few percent for |f_NL| < O(5), yet Sec. 3.3.3 keeps C_th(q) unchanged while using strongly non-Gaussian examples (µ* = 5/2, r_dec = 0.1). Because β depends exponentially on the threshold, even a 10% shift in C_th changes the amplitude A required for a fixed abundance by O(1) and propagates into all derived abundance and constraint results. A dedicated numerical relativity check of the dominant non-Gaussian profiles, or a clearly stated restriction of the universal claim to Gaussian curvature perturbations, is needed before the 'always 2/5' statement can be accepted.
- [Sec. 3.2.2, Eq. (3.58); Sec. 3.2.3, Eq. (3.68)] Several phenomenological conclusions are obtained by fixing the amplitude A to match a reference abundance (e.g., β ≃ 10^-16 in Fig. 3.8, or f_PBH = 1 in Fig. 3.11) and then interpreting the resulting mass function as a prediction. This is a consistency constraint rather than a parameter-free prediction: because β depends exponentially on A and C_th, a different threshold choice is reabsorbed by an O(1) change in A, as the thesis itself acknowledges around Figs. 3.5 and 3.8. The text should state more explicitly in each application which results are robust to this tuning and which are only illustrative; otherwise the distinction between 'prediction' and 'fit' is blurred in the abstract and in Chapter 5.
- [Sec. 3.2.1, Eq. (3.45); Sec. 3.3.4, Fig. 3.17] The manuscript uses two different thresholds for the same broad-power-spectrum case: C_th = 0.56 from the average-profile prescription (Eq. 3.45) and C_th = 2/5 from the broad-profile prescription (Section 3.3). The comparison in Fig. 3.17 is performed for a log-normal spectrum with Δ = 1 and fixed q = 0.5 and shows a ratio β/β_0 of order one, but this is a single benchmark. Because the abundance is exponentially sensitive to the threshold, a systematic scan in Δ, µ* (or r_dec), and q is needed to establish that the two prescriptions agree within the stated accuracy, or to justify replacing the 0.56 threshold by 2/5.
minor comments (5)
- [Sec. 3.2.2] The figure references 'fig.2.10' should read 'fig.3.10' in the three places where they appear in the text.
- [Sec. 4.1] The heading 'Classification of the costraints' contains a typo; it should be 'constraints'.
- [Eq. (3.43)] The formula for C_th(q) is typeset in a way that is ambiguous: the incomplete gamma function should be written with explicit parentheses, e.g., Γ(5/(2q)) − Γ(5/(2q), 1/q), to avoid confusion.
- [Throughout, especially Eqs. (3.27), (3.30), (3.78)] The symbol γ is overloaded: γ_m in Eq. (3.27), the critical-collapse exponent γ in Eq. (3.30), and the correlation coefficient γ in Eq. (3.78). Distinct symbols or a short nomenclature table would improve readability.
- [General structure] The thesis interleaves textbook review material with original computations, and it would benefit from an explicit statement in each chapter separating the author's own results from previously published work with which the thesis is in dialogue.
Circularity Check
No circularity found: the central 2/5 claim is a derived profile-selection rule applied to an externally calibrated threshold function, and the tuned amplitudes are benchmark normalizations rather than hidden predictions.
full rationale
The derivation chain in Chapter 3 is not circular. Section 3.2 introduces threshold statistics on the compaction function with the master formula (Eq. 3.58), and Section 3.3 derives, for broad power spectra, that the PBH abundance is dominated by the smallest-curvature (broadest) compaction profiles; this is done analytically for Gaussian fields (Eqs. 3.90-3.94) and numerically for a non-Gaussian ultra-slow-roll example (Fig. 3.16). The value 2/5 is then taken from the independent numerical-relativity fit C_th(q) in Eq. (3.43), whose small-q limit is given as C_th(q << 1) ≈ 2/5. The statement 'always 2/5' is thus a combination of a derived profile-selection rule and an externally calibrated threshold value, not an identity between output and input. The tuning of the amplitude A to beta ≈ 10^-16 or f_PBH = 1 in Figs. 3.8 and 3.11 is a normalization convention for benchmark models, not a claim that the normalization itself was predicted. The thesis explicitly acknowledges in footnote 6 that primordial non-Gaussianities may shift delta_c by a few percent while keeping C_th(q) unchanged; that is an extrapolation beyond the calibrated domain and a correctness risk, but not a circular reduction because the threshold function is not defined in terms of the abundance being computed. Self-citations are numerous (refs. [1], [6], [7], [8]), but the load-bearing numerical thresholds and the LVK/PTA data are external, and the thesis reproduces the derivations rather than merely citing them. No step of the claimed derivation chain reduces by construction to its own inputs.
Assumptions & free parameters
free parameters (5)
- Amplitude A of curvature power spectrum =
e.g., 9.2e-3 for broad spectrum (Sec 3.2.2)
- Threshold C_th =
0.56 for broad spectra (Sec 3.2); 2/5 in Sec 3.3
- Critical collapse constants K and gamma =
K ~ 3.3 (log-normal) or 4.36 (broad), gamma ~ 0.36
- r_dec (curvaton decay fraction) =
varied in [0.1, 1] in benchmarks
- mu_star (USR parameter) =
taken as free (e.g., 5/2 in Sec 3.3.3)
assumptions (5)
- domain assumption PBH formation occurs when the compaction function exceeds a threshold C_th (Sec 3.2.1).
- domain assumption The curvature perturbation is a local function of a single Gaussian field, ζ = F(ζ_G) (Eq 3.2).
- domain assumption The threshold statistics (Press-Schechter-like) approach is valid for computing PBH abundance (Eq 3.30).
- domain assumption Benchmark power spectra (log-normal, broad) are representative (Eqs 3.25, 3.26).
- standard math Gaussianity of the underlying field ζ_G and conservation of probability (Eqs 3.15, 3.50).
Cite this review
Pith. "Pith review of Cosmic Whispers of the Early Universe: Gravitational Waves and Dark Matter from Primordial Black Holes." pith.science (2026). https://pith.science/paper/RAXSKJSE
@misc{pith2026250103065,
author = {Pith},
title = {Pith review of: Cosmic Whispers of the Early Universe: Gravitational Waves and Dark Matter from Primordial Black Holes},
year = {2026},
howpublished = {\url{https://pith.science/paper/RAXSKJSE}},
note = {Machine review of arXiv:2501.03065}
}
read the original abstract
Various mechanisms have been suggested for the formation of Primordial Black Holes (PBHs), and this thesis focuses on the standard mechanism based on the critical collapse of cosmological fluctuations. The underlying idea is that during inflation, a period of rapid expansion in the early Universe, large cosmological fluctuations could have been generated. After inflation, when the cosmic horizon reached a size comparable to these fluctuations, if the latter were high enough, they could collapse and form a PBH. Beyond the fascinating possibility that these compact objects might make up all or part of the Dark Matter (DM) we observe today, their formation and existence is also associated with the generation of gravitational waves (GWs). These waves could contribute to the merger events observed by the LIGO/Virgo/KAGRA Collaboration (LVK) or account for signals detected by pulsar timing array experiments (PTA). In the first part of this thesis, we investigate the PBH scenario, examining the computation of the abundance beyond the Gaussian paradigm. In the second part, we discuss how PBH formation can produce a stochastic GW background and how observations, related to recent experiments such as LVK and PTA collaborations, can help in distinguishing between different PBH production models. Finally, in the last part, we investigate some aspects of the interplay between black holes and fundamental physics in the early Universe, describing some characteristics and challenges of various PBH production models.
Figures
Figures from the paper (65 more)
Reference graph
Works this paper leans on
-
[1]
IfthemodeiswayoutsidethehorizonatthebeginningoftheUSRphase, itstaysconstant even though its derivative exponentially grows because of the term∼ e−(3−2ηII)N
-
[2]
The curvature perturbation (and its derivative) grows because of the factor e−(3/2−ηII)N
Consider a mode that crosses the Hubble horizon during the USR phase. The curvature perturbation (and its derivative) grows because of the factor e−(3/2−ηII)N. However, it is not immediate to find the exact scaling in time because in this case none of the approximations in eq.(C.4) can be applied. All the above features, even though obtained in the contex...
-
[3]
+ 3W−1(f∆N⋆)] ¯MPl , (C.17) where W−1(z) is the branch withk = −1 of the Lambert W functionWk(z). The scalar spectral index at the CMB pivot scale reads ns = 1 − 16 3[1 + W−1(f∆N⋆)]2 + 8 3[1 + W−1(f∆N⋆)] , (C.18) while for the tensor-to-scalar ratio we find r = 64 3[1 + W−1(f∆N⋆)]2 . (C.19) The amplitude of the scalar power spectrum at the CMB pivot scale...
work page 2024
-
[744]
Could MACHOS be primordial black holes formed during the QCD epoch? Phys
Karsten Jedamzik. Could MACHOS be primordial black holes formed during the QCD epoch? Phys. Rept., 307:155–162, 1998. arXiv:astro-ph/9805147, doi:10.1016/S0370-1573(98) 00067-2
arXiv 1998
-
[745]
Albert Escrivà, Eleni Bagui, and Sebastien Clesse. Simulations of PBH formation at the QCD epoch and comparison with the GWTC-3 catalog.JCAP, 05:004, 2023. arXiv:2209.06196, doi:10.1088/1475-7516/2023/05/004
arXiv 2023
-
[746]
Andrew D. Gow, Christian T. Byrnes, Philippa S. Cole, and Sam Young. The power spectrum on small scales: Robust constraints and comparing PBH methodologies.JCAP, 02:002, 2021. arXiv:2008.03289, doi:10.1088/1475-7516/2021/02/002
arXiv 2021
-
[747]
Gravitational Waves Induced by Scalar Perturbations with a Lognor- mal Peak.JCAP, 09:037, 2020
Shi Pi and Misao Sasaki. Gravitational Waves Induced by Scalar Perturbations with a Lognor- mal Peak.JCAP, 09:037, 2020. arXiv:2005.12306, doi:10.1088/1475-7516/2020/09/037
arXiv 2020
-
[748]
Log-dependent slope of scalar induced gravitational waves in the infrared regions.Phys
Chen Yuan, Zu-Cheng Chen, and Qing-Guo Huang. Log-dependent slope of scalar induced gravitational waves in the infrared regions.Phys. Rev. D, 101(4):043019, 2020.arXiv:1910. 09099, doi:10.1103/PhysRevD.101.043019
Show all 15 references
-
[749]
Constraints on scalar-induced gravitational waves up to third order from a joint analysis of BBN, CMB, and PTA data.Phys
Sai Wang, Zhi-Chao Zhao, and Qing-Hua Zhu. Constraints on scalar-induced gravitational waves up to third order from a joint analysis of BBN, CMB, and PTA data.Phys. Rev. Res., 6(1):013207, 2024. arXiv:2307.03095, doi:10.1103/PhysRevResearch.6.013207
2024 arXiv
-
[750]
Andrei D. Linde. Eternally existing selfreproducing inflationary universe.Phys. Scripta T, 15:169, 1987. doi:10.1088/0031-8949/1987/T15/024
1987 doi
-
[751]
A. S. Goncharov, Andrei D. Linde, and Viatcheslav F. Mukhanov. The Global Struc- ture of the Inflationary Universe. Int. J. Mod. Phys. A, 2:561–591, 1987. doi:10.1142/ S0217751X87000211
1987
-
[752]
Largefieldpolynomialinflation: parameterspace, predictionsand (double) eternal nature.JCAP, 12:005, 2022.arXiv:2209.07545, doi:10.1088/1475-7516/ 2022/12/005
ManuelDreesandYongXu. Largefieldpolynomialinflation: parameterspace, predictionsand (double) eternal nature.JCAP, 12:005, 2022.arXiv:2209.07545, doi:10.1088/1475-7516/ 2022/12/005
2022 arXiv
-
[753]
Gauging Fine-Tuning.Phys
Feraz Azhar and Abraham Loeb. Gauging Fine-Tuning.Phys. Rev. D, 98(10):103018, 2018. arXiv:1809.06220, doi:10.1103/PhysRevD.98.103018
2018 arXiv
-
[754]
Axion isocurvature fluctuations with extremely blue spectrum
Shinta Kasuya and Masahiro Kawasaki. Axion isocurvature fluctuations with extremely blue spectrum. Phys. Rev. D, 80:023516, 2009. arXiv:0904.3800, doi:10.1103/PhysRevD.80. 023516
2009 arXiv
-
[755]
Dark matter and the anthropic principle.Phys
Simeon Hellerman and Johannes Walcher. Dark matter and the anthropic principle.Phys. Rev. D, 72:123520, 2005. arXiv:hep-th/0508161, doi:10.1103/PhysRevD.72.123520. 243 Appendix A Appendix of Chapter 3 A.1 On the convergence of the power-series expansion Consider the power seri...
2005 arXiv
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.