REVIEW 3 major objections 7 minor 1 cited by
Rare profile shapes can dominate primordial black hole production
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · glm-5.2
2026-07-10 02:01 UTC pith:IAVIU3H6
load-bearing objection Finite-action framework for PBH profile dispersion is well-built, but orders-of-magnitude abundance claims rest on a single shape direction the 3 major comments →
The statistics of curvature-profile dispersion in primordial black hole formation
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The dominant contribution to PBH abundance is selected by minimizing a combined cost: the Gaussian action of realizing a coherent shape deformation plus the squared collapse threshold along that branch. This competition can favor rare, several-sigma profile deformations over the mean profile, because the exponential sensitivity of PBH abundance to the collapse threshold can overwhelm the Gaussian suppression of the deformation. The effect is strongest for broad power spectra (where radial shape freedom is large) and for negative non-Gaussianity (where the threshold curve develops strongly shape-dependent branches). In the examples studied, the abundance enhancement reaches factors of 2 to 3,
What carries the argument
The finite-action deformation framework: the power spectrum defines a Gaussian-action metric on the space of curvature profiles, and coherent deformations are expanded in action-normalized modes orthogonal to the BBKS peak variables. Each mode carries a Gaussian cost n² (its action) and shifts the collapse threshold μ_c(n). The quantity being implicitly optimized is the critical action W_c(n) = n² + ν_c²(n), where ν_c = μ_c/σ₀. A deformed branch dominates when its lower threshold compensates its higher Gaussian cost.
Load-bearing premise
The dominance of rare shape deformations is demonstrated only within a one-parameter family of spherical radial deformations. The full shape space includes ellipticity, prolateness, and higher angular multipoles, all held at their mean values. If threshold variation along additional directions is smaller, or if the combined Gaussian cost of multi-mode deformations is higher, the enhancement from shape dispersion could be reduced.
What would settle it
If the collapse threshold were found to depend only weakly on residual radial shape modes once the BBKS height and curvature are fixed (i.e., if μ_c(n) were nearly flat across all action-normalized deformation directions), then the Gaussian cost n² would always dominate and the mean profile would always give the dominant contribution, falsifying the central claim.
If this is right
- PBH abundance calculations that fix only the mean or representative profile may underestimate the true abundance by orders of magnitude, particularly for broad power spectra or negative non-Gaussianity.
- The power-spectrum amplitude inferred from a given PBH abundance target can be significantly lower than no-dispersion estimates suggest (by factors of 0.3–0.6 in the examples studied).
- For scenarios using negative non-Gaussianity to suppress PBH production while retaining scalar-induced gravitational waves, profile dispersion may strengthen the tension between PBH constraints and gravitational-wave signals.
- The monochromatic-spectrum approximation eliminates radial shape dispersion entirely and may therefore be particularly misleading for abundance estimates.
- Extending the framework to non-spherical modes (ellipticity, prolateness, higher multipoles) could further modify the dominant branch, especially at low peak heights where the high-peak near-spherical approximation is less restrictive.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper develops a finite-action framework for describing curvature-profile dispersion in primordial black hole (PBH) formation. Using a multipolar Fourier–Bessel decomposition, the authors separate BBKS peak variables (height, curvature, ellipsoidal Hessian) from residual coherent deformations, each normalized by their Gaussian action cost. The framework is applied to two spherical numerical-collapse examples: (A) a sharply peaked finite-width spectrum with logarithmic local non-Gaussianity, and (B) a finite-band scale-invariant spectrum with tunable bandwidth. The central result is that the dominant contribution to the PBH abundance is selected by a competition between the Gaussian cost of realizing a coherent shape deformation and the exponential gain from a reduced collapse threshold—not necessarily by the mean profile or the profile with the lowest threshold. For broad spectra or negative non-Gaussianity, rare several-sigma shape deformations can enhance the integrated abundance by orders of magnitude. The theoretical framework is internally consistent: the Fourier–Bessel decomposition, action normalization, and recovery of the standard BBKS profile are clean. The numerical simulations use a public code (SPriBHoS) and threshold curves are compared against analytical estimates (EGS, HYK) with stated deviations (~4%).
Significance. The paper addresses a genuinely important gap in PBH abundance calculations: the statistical treatment of profile shape dispersion around reference peak profiles. The finite-action framework is a principled contribution that connects the Gaussian probability cost of coherent deformations to their effect on the collapse threshold, and the key insight—that the dominant profile is selected by a competition between statistical cost and threshold reduction—is well-motivated and potentially impactful. The use of a public, reproducible numerical code (SPriBHoS) and the comparison with analytical threshold prescriptions are commendable. The identification of specific regimes (broad spectra, negative non-Gaussianity) where shape dispersion can dominate the abundance is a falsifiable and quantitative prediction. However, the quantitative dominance claims rest on a restricted slice of the full shape space (a single radial deformation direction in spherical symmetry), which limits the scope of the conclusions at present.
major comments (3)
- §5.1.1, Fig. 13, and Table 2: The orders-of-magnitude abundance enhancements (e.g., β_disp/β0 ≈ 2.78×10^5 for β_NG = −2) depend entirely on the collapse threshold μ_c(s) computed along one specific deformation direction: the equal-variance split mode (Eq. 5.13). The paper states this mode is chosen to 'maximize the real-space profile dispersion,' not because it is a generic or typical direction in shape space. Since the threshold depends on the compaction function (involving radial derivatives of ζ), the direction that maximizes pointwise profile dispersion need not maximize threshold variation per unit Gaussian cost. The dominant branches occur at large deviations (n* ≈ −5.87 for β_NG = −2; n* ≈ −4.5 for R_F = 50), where the deformed profiles develop split-mode-specific features (oscillatory tails, branch switching). A different action-normalized shape direction could produce qualitativ
- §5.2.1, Eq. (5.104): The only consistency check provided is the projection of the split mode against the BBKS curvature direction q_x. This removes a component of the same split mode; it does not test an independent shape direction. The full integral (Eq. 2.91) is over all residual modes, but only one is sampled. The paper should more prominently acknowledge this limitation in the quantitative claims, or ideally test at least one additional independent radial deformation direction to assess whether the enhancement is a robust feature of the shape space or an artifact of the split-mode choice.
- §5.1.2, Eq. (5.46) and Table 2: For β_NG = −2 and β_NG = −3, the dominant branches occur at large negative s (n* ≈ −5.87 and −3.59 respectively), where the profiles develop features such as type-I/type-II transitions and central underdensities (Fig. 11). The logarithmic non-Gaussian map (Eq. 5.46) is used as a phenomenological template, and the paper notes (footnote 6) that deviations from the exact δN result can be ~8%. At these large deformations, the validity of the logarithmic template itself may be questionable. The paper should discuss whether the dominance result is robust to the non-Gaussian map prescription, or at minimum state this as a caveat on the quantitative enhancements.
minor comments (7)
- §2.4, Eq. (2.58): The standard BBKS profile is recovered, but the notation switches between γ and γ_BBKS in the surrounding equations. Using one symbol consistently would help the reader.
- §5.1, Eq. (5.34): The definition of the split mode G(x) in Eq. (5.31) involves the error function with complex argument. A brief note on how this is evaluated numerically would improve readability.
- Table 2: The column header 'f_shape_PBH / f^(n=0)_PBH' is somewhat ambiguous. Clarifying whether this is the mass-function-integrated ratio or the peak-abundance ratio would help.
- Fig. 10, upper-left panel: The threshold curves for different β_NG values are color-coded but the vertical dotted lines marking s_div are only shown for β_NG = −3 and −2. Adding these for all relevant β_NG values, or stating which curves have them, would improve clarity.
- §5.2, Eq. (5.60): The top-hat spectrum is introduced with Θ(k−k_min)Θ(k_max−k), but the relationship between R_F = k_max/k_min and the e-fold width L = ln R_F is only stated later (in the Conclusions). Stating this relationship at first introduction would help.
- §6, final paragraph: The statement that 'profile dispersion can become quantitatively important for enhanced power-spectrum features extending over several e-folds' could be made more concrete by specifying the approximate threshold width above which the effect becomes significant.
- The paper uses both 'logarithmic non-Gaussianity' and 'logarithmic local non-Gaussianity' in different places. Standardizing the terminology would improve consistency.
Circularity Check
No significant circularity; framework is derived from first principles and thresholds come from independent numerical simulations
full rationale
The paper's central framework is constructed from standard, independently verifiable ingredients: the Gaussian action (Eq. 2.5), the Fourier-Bessel decomposition (Eq. 2.2), and BBKS peak theory (Eqs. 2.80-2.84). The collapse thresholds μ_c(s) are obtained from independent relativistic numerical simulations using the SPriBHoS code, not from the statistical framework itself. The key quantitative results (Tables 2-3) combine these independently determined thresholds with the Gaussian action cost n², which is a direct consequence of the orthonormality condition (Eq. 2.8) — a standard construction, not a fitted or self-defined relation. The split-mode deformation (Eq. 5.13) is an explicit, action-normalized spectral direction satisfying the orthonormality constraint; it is not defined in terms of the abundance result it is used to compute. Self-citations (Refs. [37], [52], [58], [62]) provide the numerical code and threshold fitting formulas, but these are independently reproducible tools, not circular dependencies. The projected mode in Sec. 5.2.1 provides a partial consistency check using a different (though related) shape direction, yielding qualitatively similar results. The paper's limitation — that only one shape direction is sampled — is a scope restriction honestly acknowledged in Sec. 6, not a circularity. The derivation chain does not reduce to its inputs by construction. The skeptic's concern that the split mode may be atypical is a validity/scope issue, not a circularity issue: the threshold curves are computed from simulations, not assumed or fitted to reproduce the abundance enhancement. Score 2 reflects minor self-citation of the author's own numerical code and threshold formulas, which are not load-bearing for the logical structure of the argument but provide tools used in the computation.
Axiom & Free-Parameter Ledger
free parameters (4)
- K_eff =
6
- k_peak =
1.1e13 Mpc^-1
- nu_ref =
8.45
- beta_NG =
various (-3 to 3)
axioms (4)
- domain assumption The primordial curvature field is Gaussian to leading order, with non-Gaussianity introduced via a local logarithmic map (Eq. 5.46).
- domain assumption PBH formation threshold is determined by the maximum of the compaction function and can be found via spherical numerical collapse for the radial shape modes studied.
- ad hoc to paper The split-mode direction is a representative direction in shape space that captures the dominant radial dispersion effect.
- standard math Critical collapse scaling M ~ (mu - mu_c)^gamma_cr with gamma_cr ~ 0.356 applies to the shape-dispersed profiles.
read the original abstract
In the standard curvature-perturbation scenario, PBHs form from the collapse of superhorizon curvature fluctuations after horizon re-entry. The predicted abundance is exponentially sensitive to the collapse threshold and hence to the shape of the primordial curvature profile. In this work we develop a finite-action framework to describe curvature-profile dispersion around representative peak profiles. Using a multipolar Fourier-Bessel decomposition, we separate the local peak variables of the Gaussian field from residual radial and angular deformations, normalized by their Gaussian action. We apply the formalism to spherical numerical-collapse examples in order to isolate the effect of radial shape dispersion. For finite-width spectra, and in the presence of logarithmic local non-Gaussianity, we compute the collapse threshold as a function of a coherent shape variable and combine the result with peak statistics. We find that the dominant contribution to the PBH abundance is not necessarily the conditional-mean reference profile, nor simply the profile with the lowest threshold. Instead, it is selected by a competition between the Gaussian cost of realizing a coherent deformation and the exponential gain associated with lowering the collapse threshold. Broad spectra and negative non-Gaussianity can make rare shape deformations dominate the abundance. In the examples studied here, the dominant branches can correspond to several-sigma coherent shape fluctuations while enhancing the integrated abundance by orders of magnitude. Equivalently, including shape dispersion can reduce the power-spectrum amplitude required to obtain a fixed PBH abundance. Our results show that residual profile dispersion is a genuine statistical ingredient in PBH formation and can be quantitatively important for accurate abundance estimates.
Forward citations
Cited by 1 Pith paper
-
Primordial black holes forming during kination: the trapped, the overdense, and the void
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,...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.