{"id":"bfb6b060-b307-426b-8f41-8b33a0182c2e","arxiv_id":"2412.10666","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"PBH abundance from delayed vacuum decay follows f_pbh ≈ M exp(-Q exp(-S3(Tp)/Tp)), so S3(Tp)/Tp super-exponentially controls how many black holes form.","lead":"This paper studies how many primordial black holes form from delayed vacuum decay in the early universe, finding that the abundance is extremely sensitive to one ratio: the bubble-nucleation action divided by temperature. The result gives model-builders a precise target quantity to compute when designing particle physics models that make dark matter from black holes.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Generality of Eq. (22) rests on unproven near-constancy of Q and M; Table II already shows ~25% Q variation, and no analytic argument explains why N(tpbh) should be parameter-independent.","rationale":"The reader identified the near-constancy of N(tpbh) as the weakest assumption, and I agree that this is the most load-bearing element of the paper's central claim. The double-exponential form itself (f_pbh = M exp(-Q e^{-x})) is internally consistent: it follows from exponentiating the expected number of nucleations, and the saddle-point approximation is checked numerically for the benchmarks. The collapse criterion δc = 0.45, while uncertain, does not threaten the super-exponential structure because the saddle-point integral is dominated by Tp and is largely insensitive to the exact position of the cutoff Td (as long as Td < Tp), so changes in δc would mostly be absorbed into M and Q. The modified-expansion-rate section is a valid perturbation: it changes H_del and hence Q, but preserves the structure. The remaining vulnerable point is exactly the claim that Q and M are so weakly parameter-dependent that S3(Tp)/Tp is the unique control parameter. The paper provides real evidence for this within the toy model—five benchmark models, three of which are tabulated—but no proof. Since the abstract and conclusions generalize to 'any framework of black hole production via delayed vacuum decay,' this unproven constancy is the gap on which the generality claim rests. My proposed test directly quantifies whether Q variations are subdominant to x variations across the allowed parameter space; if they are, the paper's claim is strengthened, and if not, the claim should be weakened to apply only to models where N(tpbh) is demonstrably large and stable. I therefore keep the reader's CONDITIONAL verdict unchanged: the core scaling is plausible and well-supported in the specific model, but the broader generality requires the additional verification I propose. I do not see grounds for rejection or acceptance without the conditional check.","tokens_in":10249,"tokens_out":12729,"duration_ms":126820,"concrete_test":"Perform a dense Monte Carlo scan over the allowed parameter region of the toy model (vary µ3, ω, and c across the interval in Eq. A.2, e.g., 20×20×20 = 8000 points). For each point, compute numerically x = S3(Tp)/Tp, Q from Eq. (19), M from Eq. (21), and f_pbh. Restrict to the observationally relevant window f_pbh ∈ [10^-10, 1]. Measure the local sensitivity ratio R = |(d ln f_pbh / d ln Q)| / |(d ln f_pbh / d x)| = |Q d ln Q / (Q e^{-x} d x)| using finite differences between neighboring points. If the median R exceeds 0.1, then variations in Q compete with variations in x and the claim that x 'super-exponentially controls' the abundance is not supported. A complementary check is to derive N(tpbh) analytically in radiation domination and show it is independent of potential parameters at leading order, rather than inferring this from a handful of tuned benchmarks.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that f_pbh ≈ M exp(-Q exp(-S3(Tp)/Tp)) with Q and M effectively constant, so that S3(Tp)/Tp 'super-exponentially controls' the PBH abundance (Eq. 22, Section II.B). The load-bearing step is the assertion that N(tpbh), and hence Q, is O(10^85) and nearly parameter-independent (Section II.B, Table II). This is verified for only three benchmark points in Table II and implicitly for five in Fig. 3; no analytic derivation is given for why N(tpbh) should be nearly constant across the parameter space of this model, let alone 'any framework of black hole production via delayed vacuum decay.' The danger is concrete: near the observable window, Pint = Q e^{-x} is O(100) (since x≈172–175 and Q≈10^77), so a fractional change in Q of, say, 25% changes ln f_pbh by roughly 25% of Pint, i.e., tens of e-folds. In the same window, a change Δx of only 0.5 changes ln f_pbh by ~Q e^{-x} Δx ≈ 50 e-folds. Thus x dominates over Q only because the scanned Δx is relatively large (~2.5 between the benchmarks). If in other parameter regions Q varies substantially on the same scale as the x variation required to traverse the observable window, the claimed 'overwhelming' control by S3(Tp)/Tp would fail, and Eq. (22) would lose most of its predictive power. The paper's own conditional sentence ('If other models also feature a large, nearly parameter-independent N(tpbh)...') acknowledges this, but the generality claim in the abstract and conclusions goes beyond what is demonstrated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies primordial black hole (PBH) formation from delayed vacuum decay during a first-order phase transition in a single-field toy model. It derives an analytic approximation for the PBH-to-dark-matter abundance ratio, f_pbh ≈ M exp(−Q exp(−S3(Tp)/Tp)) (Eq. 22), using a saddle-point evaluation of the no-decay probability and treating the prefactors M and Q as effectively constant. The approximation is checked against numerical integration for several benchmark models, and the analysis is extended to a modified, superfast expansion rate driven by an additional energy component, where the same super-exponential structure is claimed to persist with a reduced value of Q. The paper argues that the result generalizes to any delayed-vacuum-decay PBH production framework.","tokens_in":10620,"tokens_out":7868,"duration_ms":75165,"significance":"If Eq. (22) is valid, it provides a practical and striking criterion: the PBH abundance is controlled overwhelmingly by the single number S3(Tp)/Tp, with a double-exponential sensitivity. A strength of the paper is that M and Q are not fitted to the target f_pbh curve; they are computed from model parameters and numerical integrals, and the analytic curve reproduces the numerical results in Fig. 3 with fixed M and Q. The saddle-point expansion is also checked, with the correction δ(Tp) below O(0.01) in Table II. The modified-expansion section is an interesting extension, but its quantitative claims are less fully validated. The main limitation is that the near-constancy of N(tpbh) and hence Q is demonstrated for only a few benchmarks, while the abstract and conclusions make a broader generality claim that goes beyond what is shown. The paper would be a useful contribution to the PBH/FOPT literature after this robustness issue is addressed or the claims are appropriately scoped.","major_comments":[{"comment":"The central claim that f_pbh is overwhelmingly more sensitive to x ≡ S3(Tp)/Tp than to M or Q rests on the assertion that N(tpbh), and hence Q, is large and nearly parameter-independent. Table II reports only three benchmark points, with Q varying from 9.4×10^76 to 1.2×10^77 (about 25%). Since Pint = Q e^{−x} is O(100) in the observable window, a 25% fractional change in Q changes ln f_pbh by tens of e-folds; this is smaller than, but not negligible compared with, the effect of the Δx ≈ 0.5 shifts emphasized in Fig. 3. The paper provides no analytic argument for why N(tpbh) should remain nearly constant over the full parameter space, and the text itself states this only conditionally: “If other models also feature a large, nearly parameter-independent N(tpbh)…” (Section II.B). Yet the abstract and conclusions assert the generality more strongly. I request either (i) a denser numerical scan that decomposes the variation in f_pbh into the contribution from x and the contributions from M and Q, or (ii) a softening of the generality claims to explicitly limit Eq. (22) to parameter regions where near-constancy of N(tpbh) is demonstrated.","section":"Section II.B, Eq. (22), Table II"},{"comment":"The modified-expansion analysis claims that the only quantity receiving a large modification is Hdel(Tp) in Q, and that N(tpbh) is not sensitive to the additional ϕ component. No analogue of Table II is provided for the n = 2 and n = 4 cases, so this assertion is not quantified. Fig. 4 shows f_pbh curves but not the values of Q or N(tpbh) under the modified cosmology. Please provide the modified-expansion analogue of Table II, or otherwise show explicitly that the O(10) reduction of Q is sufficient to explain the enhancement and that the super-exponential structure remains numerically validated in those cases.","section":"Section II.C, Eqs. (23)–(26), Fig. 4"},{"comment":"The collapse criterion δc = 0.45 is an input assumption, and footnote 1 acknowledges that this criterion has been questioned in the literature and defers a full discussion. Because f_pbh changes by many orders of magnitude for small changes in the effective collapse condition, the quantitative statements in the paper—such as the window S3(Tp)/Tp ∈ [171.5,175.5] for successful PBH formation—are conditional on this choice. Please provide a robustness check for δc over a plausible range, or state explicitly in the abstract and conclusions that the quantitative window and the values of f_pbh assume δc = 0.45.","section":"Section II, Eq. (6), footnote 1"}],"minor_comments":[{"comment":"There are several typos: “Plank Mass” should be “Planck mass” in footnote 2; “Big Bag Nucleosynthesis” should be “Big Bang Nucleosynthesis” in Section II.C; “M read” should be “M reads” in Eq. (21); “as been shown in Table.II” should be “as shown in Table II” in Section II.B.","section":"Footnotes and text"},{"comment":"The labels BMa, BMb, and BMb do not directly correspond to BM1–BM5 in Table I, and the listed μ3 values appear to be off-peak values rather than the μ3* values defined in Table I. Please clarify which benchmark configurations are used and why these particular μ3 values were chosen.","section":"Table II"},{"comment":"The solid blue line in Fig. 3 uses fixed values M = 4×10^7 and Q = 1×10^77. The text should state explicitly that this is not a fit but the analytic prediction with these computed values; this would strengthen the validation claim.","section":"Fig. 3"},{"comment":"The choice Tr = 50 MeV is not varied, although the BBN bound only requires Tr ≳ O(10) MeV. A sentence on the sensitivity of the modified-expansion results to Tr would improve the robustness discussion.","section":"Section II.C, Eq. (24)"},{"comment":"In the saddle-point expansion, the A'(Tp)Dy term is dropped because it integrates to zero. The text could state this explicitly to avoid the impression that the term was omitted without justification.","section":"Section II.B, Eq. (16)"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for a hep-ph journal and makes a worthwhile contribution. The main unresolved issue is the robustness of the near-constancy of Q and M; the current evidence is limited to three benchmark points in Table II, while the abstract and conclusions claim broader generality. This is fixable by additional numerical scans or by carefully limiting the claims, so I recommend major revision rather than rejection. The paper's own conditional sentence in Section II.B is an important caveat that should be reflected in the abstract and conclusions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe one thing to know: this paper derives a compact double-exponential formula for PBH production from delayed vacuum decay, f_pbh ≈ M exp(-Q exp(-S3(Tp)/Tp)), and then checks it numerically rather than just asserting it. The check is genuine. The saddle-point approximation is validated with δ(Tp) below O(0.01) for five benchmarks, and the analytic curve with fixed M and Q reproduces the numerical f_pbh across a wide dynamic range. M and Q are computed from the model parameters, not fitted to the target curve, so there is no circularity. That part is solid and genuinely useful for model builders.\n\nWhat is actually new: the explicit analytic expression itself, and the observation that a faster-than-radiation expansion phase weakens but preserves the super-exponential structure. Prior work by Kawana et al. and Kanemura et al. had the mechanism; this paper turns the numerical story into a one-line scaling law.\n\nThe soft spot, and it is the main one, is the generality claim. The abstract says the findings generalize to any framework of black hole production via delayed vacuum decay. What is demonstrated is that Q and N(tpbh) are nearly constant across three benchmark points in one toy model. That is not nothing, but it is not a general proof. The paper itself hedges later: if other models also feature a large, nearly parameter-independent N(tpbh), then variations in Q and M are subdominant. That conditional sentence is doing real work, because a 25% change in Q shifts ln f_pbh by tens of e-folds in the observable window. The reason S3/Tp appears to dominate in the scans is that the scanned range of S3/Tp (about 2.5 units) is large compared to the Q variation. In a different model where Q varies more steeply with the same parameters that move S3/Tp, the double-exponential form might still hold but the overwhelming control by S3/Tp alone would not. The abstract and conclusions should be moderated to match what is actually shown.\n\nTwo lesser points. The collapse criterion δc = 0.45 is taken as given; the paper footnotes the debate and defers to Cai et al., which is honest but means the quantitative predictions inherit that uncertainty. And the modified-expansion section is thinner than the rest—one benchmark, n=2 and 4, no table of the same quantities that were checked in the main analysis. Plausible, but not at the same verification level.\n\nWho should read this: anyone computing S3/T for a first-order phase transition and wanting a quick estimate of whether delayed vacuum decay can produce observable PBHs, plus people interpreting PBH constraints. It deserves a serious referee. My recommendation: send it to review, but ask the authors to either prove or soften the generalization, and to scan a wider parameter range showing Q and M stability.\n\nBest,","headline":"Clean derivation of the double-exponential f_pbh formula, but the 'any framework' generality claim is not supported by the evidence.","tokens_in":11142,"tokens_out":4452,"would_cite":true,"duration_ms":36305,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Delayed vacuum decay makes primordial black hole abundance super-exponential in the Euclidean action-to-temperature ratio.","keywords":["primordial black holes","first-order phase transition","delayed vacuum decay","bubble nucleation","Euclidean action","super-exponential sensitivity","modified expansion rate","dark matter candidate"],"falsifier":"Recompute $f_{\\mathrm{pbh}}$ without truncating the integral in Eq. (14) for a benchmark where the integrand's peak is broad or where $T_d$ and $T_c$ are close; if the full numerical integral differs from $Q\\exp(-S_3(T_p)/T_p)$ by more than the prefactor variation, the super-exponential approximation fails. A second, parameter-space test: choose a model with much smaller $N(t_{\\mathrm{pbh}})$ and see whether $f_{\\mathrm{pbh}}$ still collapses onto the universal curve versus $S_3(T_p)/T_p$.","tokens_in":10005,"feed_emoji":"🕳️","tokens_out":8675,"duration_ms":74272,"temperature":0.7,"pith_summary":"In a cosmological first-order phase transition, patches of false vacuum that decay late can become overdense and collapse into primordial black holes. This paper studies a minimal single-scalar model of that delayed vacuum decay and shows that the resulting relic abundance is controlled by a double exponential, $f_{\\mathrm{pbh}}\\simeq M\\exp\\!\\left(-Q\\,e^{-S_3(T_p)/T_p}\\right)$, where $S_3$ is the three-dimensional Euclidean action of the nucleating bubble, $T_p$ is the temperature at which $S_3(T)/T$ attains its minimum, and $M$ and $Q$ are slowly varying prefactors. The consequence is that percent-level shifts in underlying potential parameters swing $f_{\\mathrm{pbh}}$ through dozens of orders of magnitude, because the exponent itself depends exponentially on $S_3(T_p)/T_p$. The paper also shows that a faster-than-radiation expansion phase, as from an extra energy-density component, enhances the PBH abundance and softens the parameter sensitivity while preserving the same super-exponential structure.","feed_headline":"One phase-transition ratio super-exponentially sets PBH abundance","feed_subtitle":"Small shifts in the tunneling action swing the predicted black hole abundance by dozens of orders of magnitude.","key_machinery":"The load-bearing object is the probability that a Hubble patch has not nucleated by a delayed time, $P(t_d)=\\exp(-P_{\\mathrm{int}})$, whose exponent is a temperature integral of $A(T)\\exp[-S_3(T)/T]$ with $A(T)$ collecting the prefactors from the nucleation rate and the Hubble volume. The machinery is the saddle-point approximation of that integral around the temperature $T_p$ where $S_3(T)/T$ is minimal; it turns the integral into $Q\\exp(-S_3(T_p)/T_p)$, with $Q\\sim 1\\times10^{77}$ dominated by the large, nearly constant factor $N(t_{\\mathrm{pbh}})\\sim10^{85}$. A prefactor $M\\sim4\\times10^7$ gathers the mass, volume, entropy and abundance factors outside the integral. The combination is what converts a tunnelling-rate suppression into the double-exponential form of Eq. (22).","core_discovery":"The central claim is that delayed vacuum decay converts the tunneling suppression into a super-exponential control of PBH production. In the paper's own terms, Eq. (22), $f_{\\mathrm{pbh}}\\simeq M\\exp(-Q\\exp(-S_3(T_p)/T_p))$, is the outcome of a saddle-point evaluation of the no-nucleation probability; for the five benchmarks the prefactors take values $M\\sim 4\\times 10^7$ and $Q\\sim 1\\times 10^{77}$, while $S_3(T_p)/T_p$ lies in the narrow interval $[171.5,175.5]$. Small changes in the cubic coupling therefore move $f_{\\mathrm{pbh}}$ over more than a hundred orders of magnitude. The authors claim the same structure persists when the expansion rate is modified by an extra energy component: the modified cosmology mainly reduces $Q$ by about an order of magnitude, which enhances $f_{\\mathrm{pbh}}$ and weakens the dependence on model parameters.","pith_inferences":["If Eq. (22) survives in other models, then observational bounds on $f_{\\mathrm{pbh}}$ can be inverted into constraints on the tunnelling action ratio $S_3(T_p)/T_p$, making PBH searches a direct probe of the nucleation barrier.","A natural numerical check is to apply the same saddle-point reduction to multi-field or multi-stage phase transitions; the super-exponential form should survive whenever the integrand is sharply peaked and $N(t_{\\mathrm{pbh}})$ stays large.","The paper itself flags that the collapse criterion $\\delta_c=0.45$ has been questioned in the literature; if a different collapse prescription were adopted, the prefactor $M$ and the required $S_3(T_p)/T_p$ window would shift, even though the super-exponential dependence might persist.","Because $Q$ is set partly by the Hubble rate at PBH formation, a measured abundance that is far above the standard-cosmology prediction would be a hint of kination-like early-universe expansion."],"forward_implications":["In the delayed-vacuum-decay scenario, the observable PBH abundance pins $S_3(T_p)/T_p$ to a narrow window; for the benchmarks here, the window is roughly $171.5$ to $175.5$.","Percent-level changes in a potential parameter such as the cubic coupling $\\mu_3$ sweep $f_{\\mathrm{pbh}}$ over more than one hundred orders of magnitude in the numerical examples.","A superfast expansion phase with an extra component scaling as $a^{-(4+n)}$ for $n=2$ or $4$ raises $f_{\\mathrm{pbh}}$ and relaxes fine-tuning while keeping the super-exponential relation intact.","The result gives a criterion for scanning first-order phase transition models: compute $S_3(T_p)/T_p$ and compare with the required window to decide whether delayed vacuum decay can produce a significant PBH population.","The authors argue that because the saddle-point structure relies on general features of the integrand, Eq. (22) extends beyond the toy model to any delayed-vacuum-decay framework."],"supporting_citations":[{"why":"Supplies the single-scalar finite-temperature effective potential used as the case study throughout the paper.","marker":"[13]"},{"why":"Initiates the delayed-vacuum-decay PBH mechanism that the paper extends and generalizes.","marker":"[3]"},{"why":"Provides the earlier analytical treatment of collapse probability and phase-transition duration that this work augments.","marker":"[4]"},{"why":"Gives the thermal tunneling rate formula whose exponent contains $S_3(T)/T$.","marker":"[14, 15]"},{"why":"Parametrizes the superfast expansion with an extra energy density scaling as $a^{-(4+n)}$ and sets the reference temperature used in the modified-cosmology analysis.","marker":"[18]"},{"why":"Defines PBH abundance conventions, the critical density contrast, and the Hubble-mass conversion used in the abundance formulas.","marker":"[1, 2]"}],"fun_headline_variants":["Super-exponential PBH yield from delayed vacuum decay","One tunneling ratio swings PBH abundance by 100+ orders","Slow phase transition: PBH count super-exponentially sensitive","Delayed vacuum decay makes PBH abundance super-exponential"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation requires that the number of relevant Hubble patches, and hence the prefactor $Q$, is enormous around $10^{85}$ and nearly independent of the model parameters; if $Q$ were small or varied strongly, changes in $M$ and $Q$ could compete with the exponential of $S_3(T_p)/T_p$ and the claimed universal super-exponential dominance would fail.","fun_headline_variants_meta":{"raw":{"variants":["Super-exponential PBH yield from delayed vacuum decay","One tunneling ratio swings PBH abundance by 100+ orders","Slow phase transition: PBH count super-exponentially sensitive","Delayed vacuum decay makes PBH abundance super-exponential"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000188,"raw_usage":{"total_tokens":1298,"prompt_tokens":878,"completion_tokens":420,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":494,"completion_tokens_details":{"reasoning_tokens":352}},"tokens_in":494,"tokens_out":420,"duration_ms":4189,"temperature":1.0,"reasoning_tokens":352,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T15:44:46.118446+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute $f_{\\mathrm{pbh}}$ without truncating the integral in Eq. (14) for a benchmark where the integrand's peak is broad or where $T_d$ and $T_c$ are close; if the full numerical integral differs from $Q\\exp(-S_3(T_p)/T_p)$ by more than the prefactor variation, the super-exponential approximation fails. A second, parameter-space test: choose a model with much smaller $N(t_{\\mathrm{pbh}})$ and see whether $f_{\\mathrm{pbh}}$ still collapses onto the universal curve versus $S_3(T_p)/T_p$.","supporting_citations":[{"cited_title":"Pbh formation from overdensities in delayed vacuum transitions,","cited_arxiv_id":null,"evidence_quote":"Provides the earlier analytical treatment of collapse probability and phase-transition duration that this work augments."}],"review_version":1}