{"id":"01d2e7f5-c946-4b36-8e29-6e02b59665e1","arxiv_id":"2412.07709","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A broad power spectrum of primordial curvature perturbations produces a bimodal primordial black hole mass function whose dominant peak is near the infrared scale, not the ultraviolet scale.","lead":"Primordial black holes could form when the early universe's density ripples span many sizes at once. This paper shows that such broad ripples make heavy black holes, tied to the largest scale, much more abundant than the standard picture assumed.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The dominant IR peak rests on the unvalidated large-w extrapolation of the collapse threshold gc(w); Fig. 2 inherits Eq. (7) and cannot independently confirm it.","rationale":"I read the paper as claiming that, for a nearly scale-invariant plateau, peak-theory statistics produce a dominant heavy-PBH population from large-w profiles, overturning the usual UV-dominated expectation. The calculation is internally consistent, and the authors provide both a direct numerical integration of Eq. (23) and analytic estimates (Eqs. (20), (24)); the agreement in Figs. 1 and 2, and the smooth-cutoff robustness check, are real supporting evidence. The load-bearing point is that the entire large-w contribution sits at the edge of the regime where the threshold gc(w) has been tested. The reader's weakest-assumption analysis identifies the same point, and the authors' own Outlook flags it. This is a quantitative fragility rather than a demonstrated error: if the threshold saturates at a value below 4/3, the IR peak could be even stronger; if it approaches 4/3 much faster than 32/(9w), the heavy-peak contribution would shrink and could lose dominance. Therefore the honest verdict remains CONDITIONAL, pending dedicated simulations.","tokens_in":11107,"tokens_out":25832,"duration_ms":250755,"concrete_test":"Perform spherically symmetric numerical-relativity collapse simulations for initial profiles generated from Eq. (3) with Pζ of Eq. (15), α=100 and As=10^{-2}, selecting realisations whose g(r) peaks at r≈1.69/kIR with w≈wmax from Eq. (19) (w∼10^3). Measure the critical value of g separating collapse from dispersal and compare with Eq. (7); also test K and γcr in Eq. (14) for such thin-shell profiles. If the measured threshold deviates from Eq. (7) by more than a few percent, or does not saturate near 4/3, recompute Eq. (23) and check whether the IR peak remains dominant.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central prediction—a dominant heavy peak associated with r≈1.69/kIR—is obtained by integrating Eq. (12) over g in the window [gc(w), 4/3]. For the values of w that dominate the IR peak (w≳10^3 for α=100), the lower limit uses Eq. (7), gc(w)≈4/3−32/(9w), an analytic fit whose large-w behavior is extrapolated from simulations that do not cover the thin-shell regime. This is not a peripheral detail: the exact numerical mass functions in Fig. 2 use the same gc(w), so they cannot independently validate the extrapolation. The authors explicitly state in the Outlook that 'dedicated simulations ... are needed to confirm the extrapolation ... of a threshold saturation in the Type I case.' If the true threshold approaches 4/3 more slowly, or with a different coefficient, the amplitude and mass scale of the IR peak—Eqs. (20) and (24), Mheavy/MUV∝α^{2−4γcr}—shift. If it approaches 4/3 much faster, the 1/w window over g shrinks and the large-w contribution may no longer outrank the UV peak. Since Type II fluctuations (g>4/3) are excluded from Eq. (12), the claim of Type II relevance is an additional extrapolation beyond the computed abundance.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter computes the PBH mass function for a flat, enhanced curvature power spectrum bounded by infrared and ultraviolet scales. Using compaction-function peak theory together with the critical-collapse mass relation, the authors integrate the abundance over smoothing radius, peak curvature w, and amplitude g, and find that the mass function is bimodal: a UV feature at r≈4/k_UV and a dominant heavy feature at r≈1.69/k_IR, with an analytic saddle-point estimate reproducing the numerical result. They interpret the heavy peak as arising from thin spherical shells whose profiles have large curvature w, close to the Type I/Type II boundary, and they discuss consequences for PBH overproduction bounds and for Type II collapse.","tokens_in":11368,"tokens_out":7479,"duration_ms":70641,"significance":"If correct, the main result overturns the standard expectation that a broad enhancement is dominated by light PBHs associated with the UV cutoff, and it would strengthen overproduction constraints in the asteroid-mass window. The paper's strengths are its transparent analytic estimates (Eqs. 17-22 and 24), which agree with the direct numerical integration, and the explicit condition in Eq. (16) that delimits the bimodal regime. The bimodality is not put in by hand; it emerges from the statistics. The principal caveat is that the dominant IR peak depends on an extrapolated form of the collapse threshold at large w, Eq. (7), and this dependence is acknowledged but not quantified. On balance the result is significant and plausible, but it needs a robustness or sensitivity analysis before publication.","major_comments":[{"comment":"The dominant IR feature rests on the large-w behavior of the collapse threshold. For α=100 and r=1.69, the relevant values of w are of order γνcσw≃(2/3)(rα)^2/ln α, i.e. thousands, far beyond the regime in which Eq. (7) has been calibrated. The exact numerical curves in Fig. 2 use the same analytic fit gc(w) in the lower integration limit of Eq. (12), so they cannot independently confirm the extrapolation. As the authors note in the Outlook, dedicated simulations are required. If the true threshold saturates more slowly, or with a different coefficient, the amplitude and scale of the IR peak (Eqs. 20 and 24, with Mheavy/MUV∝α^{2−4γcr}) shift; if it saturates more quickly, the [gc(w),4/3] window narrows and the IR contribution may no longer dominate. I therefore regard the central prediction as conditional on Eq. (7) and ask for an explicit sensitivity analysis, for example using alternative asymptotic forms for gc(w), in the revised manuscript.","section":"Sec. III, Eq. (7) and Fig. 2"},{"comment":"The claim that Type II fluctuations may contribute significantly is not backed by a computation. Equation (12) restricts the g integral to g≤4/3, so all Type II fluctuations are excluded by construction. The fact that the dominant Type I configurations are close to the boundary is suggestive, but the statement in the Abstract that the results imply a \"higher-than-expected abundance of PBH originating from Type II initial fluctuations\" goes beyond the quantitative content of the paper. Please either include an explicit estimate of a Type II contribution with suitable modeling, or clearly soften the claim to a qualitative conjecture.","section":"Sec. V and Eq. (12)"},{"comment":"The robustness of the result for power spectra without sharp cutoffs is asserted in footnote 2 but not demonstrated. Since Eq. (22) contains sharp-boundary terms such as CosIntegral(2r), and since the IR peak is located at r≈1.69, it is important to check that this feature is not a cutoff artifact. A figure or a short quantitative statement showing that smooth cutoffs produce the same bimodal structure would remove this uncertainty.","section":"Footnote 2 and Eq. (22)"}],"minor_comments":[{"comment":"There is a duplicated word in \"As shown in in Figure 2\"; please correct it.","section":"Sec. IV"},{"comment":"There is a duplicated word in \"a wide class of of profiles\"; please correct it.","section":"Sec. II"},{"comment":"Equation (16) is quoted as the bimodality condition, but the derivation of the numerical coefficient 0.05 is not shown; a brief derivation or reference would help the reader assess its range of validity.","section":"Sec. III, Eq. (16)"},{"comment":"The phrase \"withv(r)=0\" should read \"with v(r)=0\".","section":"Sec. II, below Eq. (3)"},{"comment":"Reference [3] is printed with a DOI but without a visible publication venue; completing the bibliographic entry would be helpful.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the paper is within scope and the authors are transparent about the main limitation. My concern is not internal consistency but the external validity of the asymptotic threshold fit. A sensitivity analysis under different gc(w) forms would significantly strengthen the paper; without it, the main quantitative claim is conditional. I saw no issue of attribution or novelty disclosure, and no sign that the result is fitted to produce the bimodality."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a real result, with a real caveat. The paper shows that for a broad, flat power spectrum, the PBH mass function is bimodal, and the heavier IR peak can dominate. That contradicts the standard expectation (Byrnes et al., De Luca et al.) of a single UV peak with an M^{-1/2} tail. The analytic estimate M_heavy/M_UV proportional to alpha^{2-4gamma_cr} is new, and the numerics in Eqs. (12)-(13) back it. The claim is not circular: bimodality emerges from the calculation rather than being put in by hand.\n\nWhat the paper does well: the peak-theory machinery is applied carefully, the saddle-point estimate matches the numerical curves, and the authors flag their own weak point explicitly. The Appendix's narrow-spectrum limit is a useful consistency check. Credit also goes to the authors for engaging with refs. [18] and [23], and for spelling out where their result departs from recent claims about broad profiles.\n\nThe soft spot is the one the authors themselves admit in the Outlook: the IR peak rests on extrapolating the Type I threshold gc(w) to w >> 1 via Eq. (7). The exact numerical mass functions in Fig. 2 use that same fit, so they cannot independently validate the extrapolation. If the true threshold approaches 4/3 more slowly, or with a different coefficient, the amplitude and location of the IR peak shift (Eqs. 20 and 24). If it saturates much faster, the IR peak may no longer outrank the UV peak. That is a load-bearing assumption, and only dedicated simulations can settle it. The Type II relevance claim is even more speculative, since the calculation explicitly excludes g > 4/3.\n\nThat said, I do not see a fatal flaw. The core structural result — that thin shells at the IR scale dominate the statistics — may well survive threshold corrections. The paper is honest, the math checks out internally, and the significance for PBH overproduction bounds is genuine if the prediction holds.\n\nWho this is for: PBH phenomenologists, especially people using peak theory to derive mass functions. It deserves a serious referee, and the referee should push hard on the large-w threshold behavior. I would send it to review rather than desk reject, and I would bring it to reading group.","headline":"Real new result on bimodal PBH mass functions, with a load-bearing but clearly flagged extrapolation of the collapse threshold; worth a serious referee.","tokens_in":11917,"tokens_out":1563,"would_cite":true,"duration_ms":14239,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["98.80.-k","04.70.-s"],"model":"deepseek-v4-flash","headline":"Broad enhancements of the primordial power spectrum make heavy black holes, not light ones, dominate.","keywords":["primordial black holes","mass function","compaction function","critical collapse","peak statistics","broad power spectrum","Type II fluctuations","dark matter"],"falsifier":"A set of numerical relativity simulations that initialize thin spherical-shell compaction profiles with $w\\gg1$ and scale separation $\\alpha\\gtrsim100$, measuring the critical threshold $g_c(w)$ and the resulting PBH mass, would settle the claim. If $g_c(w)$ does not keep rising toward $4/3$, or if such shells do not collapse, the predicted infrared peak at $r\\simeq1.69/k_{\\rm IR}$ and the scaling $M_{\\rm heavy}/M_{\\rm UV}\\propto\\alpha^{2-4\\gamma_{\\rm cr}}$ fail. A cheaper partial test is to evaluate the full mass-function integral for a smooth broad spectrum and check whether the bimodal shape and infrared dominance survive.","tokens_in":10915,"feed_emoji":"🕳️","tokens_out":11300,"duration_ms":90238,"temperature":0.7,"pith_summary":"This paper establishes a counterintuitive result for primordial black holes (PBHs) formed from nearly Gaussian curvature perturbations: when the curvature power spectrum is enhanced over a broad range of scales, between an infrared scale $k_{\\rm IR}$ and an ultraviolet scale $k_{\\rm UV}$, the PBH mass function is bimodal, and the dominant, heavier peak sits near the infrared scale rather than the ultraviolet one. Standard estimates for a broad, nearly scale-invariant spectrum predict that light black holes formed at the ultraviolet scale dominate, because they form earlier during radiation domination. The authors show analytically and numerically that the statistics of peaks select thin spherical shells with infrared radius and ultraviolet thickness, pushing the typical collapsing configurations toward large values of the curvature parameter $w$ and toward the Type I/Type II boundary. If this is right, predicted PBH abundances are higher at large masses, so overproduction bounds on the power spectrum amplitude tighten and the mass range in which PBHs could be all of the dark matter narrows.","feed_headline":"Broad early-universe ripples make heavy black holes dominate","feed_subtitle":"If right, overproduction bounds tighten and the asteroid-mass PBH dark-matter window narrows.","key_machinery":"The engine is the joint peak statistics of the compaction function $C=g(1-3g/8)$, with $g=-(4/3)r\\,\\partial_r\\zeta$, conditioned on a radial extremum $v=0$; the curvature parameter $w\\equiv-r^2\\partial_r^2g$ controls the collapse threshold $g_c(w)$. For a broad spectrum the correlators simplify to $\\tilde\\gamma\\simeq\\gamma\\simeq1/\\sqrt{\\ln\\alpha}\\ll1$ and $\\sigma_w^2\\simeq(2/9)A_s(r\\alpha)^4$, and the saddle point of the abundance integrand, Eq.~(17), sits at $w_{\\max}/\\sigma_w=(\\gamma\\nu_c+\\sqrt{\\gamma^2\\nu_c^2+4\\lambda})/2$ with $\\lambda\\simeq2.28$ set by the large-$w$ peak-theory shape factor $f\\simeq(w/\\sigma_w)^3$. Combined with $\\sigma_g^2\\simeq(8/9)A_s[\\ln\\alpha+{\\rm CosIntegral}(2r)+(\\sin^2r-r\\sin2r)/r^2]$, whose maximum is at $r\\simeq1.69/k_{\\rm IR}$, this yields the infrared peak and the analytic approximation Eq.~(24). The condition $A_s(\\ln\\alpha)^2\\gtrsim0.05$ separates the bimodal regime from a single-peak regime whose abundance is negligible.","core_discovery":"On the paper's own terms, the central claim is that the PBH mass function from a flat plateau $P_\\zeta(k)=A_s\\,\\theta(k-k_{\\rm IR})\\theta(k_{\\rm UV}-k)$ with $\\alpha=k_{\\rm UV}/k_{\\rm IR}\\gtrsim25$ develops two peaks: a subdominant ultraviolet peak near $r\\simeq4/k_{\\rm UV}$ and a dominant infrared peak near $r\\simeq1.69/k_{\\rm IR}$, with the heavy peak scaling as $M_{\\rm heavy}/M_{\\rm UV}\\propto\\alpha^{2-4\\gamma_{\\rm cr}}$ for $\\gamma_{\\rm cr}\\simeq0.36$. The infrared peak arises because, conditioned on meeting the collapse threshold, the most probable compaction function is a thin spherical shell of radius $\\sim1/k_{\\rm IR}$ and thickness $\\sim1/k_{\\rm UV}$, corresponding to large $w$ and to $g$ close to $4/3$, i.e. near the Type I/Type II boundary. The paper derives this result with a saddle-point estimate of the peak-theory integral and verifies it by numerically evaluating the full mass-function integral, for power-spectrum amplitudes satisfying $A_s(\\ln\\alpha)^2\\gtrsim0.05$.","pith_inferences":["A direct test would be numerical relativity simulations that evolve thin-shell compaction profiles with $w\\gg1$ and scale separation $\\alpha\\gtrsim100$; the paper explicitly states that such simulations are needed, and their outcome would confirm or refute the infrared peak.","The same peak-statistics mechanism should operate for smooth broad spectra, not just the sharp plateau used for analytic control; the paper reports numerical robustness, so a lognormal or running-power-law enhancement is a clean place to look for the same bimodal signature.","If the infrared peak is real, mapping a fixed gravitational-wave stochastic background onto $A_s$ should be revised: a given background would be compatible with a lower amplitude than single-peak estimates suggest, because the heavy population is more abundant.","The result implies that replacing a broad spectrum by a Dirac-delta spectrum of the same logarithmic area, a common shortcut, systematically underestimates the heavy-mass tail; only the ultraviolet peak is captured by that shortcut."],"forward_implications":["If the central claim holds, standard estimates for broad spectra miss the dominant heavy-PBH population, so overproduction bounds on $A_s$ tighten because the infrared peak adds a large-mass contribution that was previously neglected.","The viable window for PBHs as all of the dark matter narrows, especially in the asteroid-mass range: the subdominant ultraviolet population must avoid Hawking-evaporation bounds, which pushes the dominant infrared population to larger masses.","Gravitational-wave and lensing constraints that assume a single-peaked mass function need to be re-evaluated for enhancements spanning $N\\gtrsim3$ e-folds, where $\\alpha\\gtrsim25$ already produces bimodality.","Typical collapsing configurations are thin shells with large $w$, so Type II fluctuations, or the Type I/Type II boundary, become statistically relevant even though their per-profile probability is small.","The narrow-spectrum approximation is not valid for the overall mass function once $\\alpha\\gtrsim25$, and even the ultraviolet peak shifts to $r\\simeq4/k_{\\rm UV}$ once $\\alpha\\gtrsim1.5$, according to the paper's numerical checks."],"supporting_citations":[{"why":"Defines the compaction function and the trapped-surface criterion that sets the collapse threshold.","marker":"[6]"},{"why":"Supplies the peak-theory phase-space factors $f(\\chi/\\sigma_\\chi)$ used in the abundance integral.","marker":"[7]"},{"why":"Establishes the threshold for PBH formation and its dependence on the radial profile.","marker":"[11]"},{"why":"Provides the peak-statistics formalism and correlators for the compaction function that the paper generalizes to broad spectra.","marker":"[12]"},{"why":"Gives the analytic approximation to $g_c(w)$ and the mass-function inversion used for Eq. (23).","marker":"[13]"},{"why":"Matches the analytic threshold $g_c(w)$ to numerical simulations to within a few percent.","marker":"[14]"},{"why":"Provides the numerical simulations used to validate the threshold behavior that the paper extrapolates to large $w$.","marker":"[15]"},{"why":"Supplies the critical-scaling mass relation $M\\simeq K M_H(C-C_c)^{\\gamma_{\\rm cr}}$ that converts profile parameters into PBH mass.","marker":"[16]"},{"why":"Represents the standard expectation of a UV-dominated mass function with a $M^{-1/2}$ dilution tail that the paper overturns.","marker":"[20]"},{"why":"Also represents the standard broad-spectrum expectation against which the bimodal result is compared.","marker":"[21]"}],"fun_headline_variants":["Broad primordial ripples yield twin black hole mass peaks","Heavy black holes dominate when spectrum broadens","Bimodal PBH masses from a plateau in curvature power","Broad curvature enhancements favor heavy black holes","Plateau ripples generate two black hole mass peaks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The heavy, infrared peak rests on extrapolating the analytic collapse threshold to very narrow, strongly curved profiles with $w\\gg1$, where the paper notes dedicated simulations are still missing; if the threshold saturates at a different value there, the peak's height and location shift.","fun_headline_variants_meta":{"raw":{"variants":["Broad primordial ripples yield twin black hole mass peaks","Heavy black holes dominate when spectrum broadens","Bimodal PBH masses from a plateau in curvature power","Broad curvature enhancements favor heavy black holes","Plateau ripples generate two black hole mass peaks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000399,"raw_usage":{"total_tokens":2101,"prompt_tokens":973,"completion_tokens":1128,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":589,"completion_tokens_details":{"reasoning_tokens":1054}},"tokens_in":589,"tokens_out":1128,"duration_ms":8555,"temperature":1.0,"reasoning_tokens":1054,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:34:26.980899+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A set of numerical relativity simulations that initialize thin spherical-shell compaction profiles with $w\\gg1$ and scale separation $\\alpha\\gtrsim100$, measuring the critical threshold $g_c(w)$ and the resulting PBH mass, would settle the claim. If $g_c(w)$ does not keep rising toward $4/3$, or if such shells do not collapse, the predicted infrared peak at $r\\simeq1.69/k_{\\rm IR}$ and the scaling $M_{\\rm heavy}/M_{\\rm UV}\\propto\\alpha^{2-4\\gamma_{\\rm cr}}$ fail. A cheaper partial test is to evaluate the full mass-function integral for a smooth broad spectrum and check whether the bimodal shape and infrared dominance survive.","supporting_citations":[],"review_version":1}