Pith. sign in

REVIEW 3 major objections 5 minor 1 cited by

Super-exponential Primordial Black Hole Production via Delayed Vacuum Decay

T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Delayed vacuum decay makes primordial black hole abundance super-exponential in the Euclidean action-to-temperature ratio.

desk verdict Clean derivation of the double-exponential f_pbh formula, but the 'any framework' generality claim is not supported by the evidence. read the letter →

arxiv 2412.10666 v2 pith:MJL445VH submitted 2024-12-14 hep-ph astro-ph.COastro-ph.HEgr-qc

Yanda Wu , Stefano Profumo This is my paper · ORCID
classification hep-phastro-ph.COastro-ph.HEgr-qc
keywords primordialblackholesfirst-orderphasetransitiondelayedvacuumdecaybubblenucleationEuclideanactionsuper-exponentialsensitivitymodifiedexpansionratedarkmattercandidate
topics Dark Matter
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

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.

What carries the argument

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).

What would settle it

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$.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

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

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

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.

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 (3)
  1. [Section II.B, Eq. (22), Table II] 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.
  2. [Section II.C, Eqs. (23)–(26), Fig. 4] 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.
  3. [Section II, Eq. (6), footnote 1] 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.
minor comments (5)
  1. [Footnotes and text] 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.
  2. [Table II] 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.
  3. [Fig. 3] 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.
  4. [Section II.C, Eq. (24)] 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.
  5. [Section II.B, Eq. (16)] 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.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: Eq. (22) is derived from the nucleation-rate definition and saddle-point evaluation, with M and Q computed from model inputs and validated numerically; the only self-citation (Ref. [18]) is a minor, non-load-bearing ansatz.

full rationale

The central result f_pbh ≈ M exp(-Q exp(-S3(Tp)/Tp)) is not obtained by fitting the target f_pbh curves. M is defined in Eq. (21) from m_pbh, V_H,nor, s0/s, ρ0, and Ωdm; Q is defined in Eq. (19) from N(tpbh), A(Tp), and S3''(Tp); both are computed from the model parameters and numerical integrals, with representative values listed in Table II. The double-exponential form follows mathematically from the survival-probability definition P(td) = exp(-Pint) in Eq. (9) combined with the saddle-point evaluation of Pint in Eqs. (14)–(20), where the nucleation rate Γ itself contains exp(-S3/T). This is a derivation, not a hidden restatement of the conclusion, and Fig. 3 shows the analytic formula against independently computed numerical f_pbh values. The paper explicitly flags its main robustness assumption: 'If other models also feature a large, nearly parameter-independent N(tpbh), then variation in remaining components of Q or in M will likewise be subdominant to variations in S3(Tp)/Tp.' This is a generality caveat about N(tpbh) near-constancy, not circularity, because the paper demonstrates O(10^85) values for its benchmarks rather than assuming them. The only self-citation is Ref. [18] (D'Eramo, Fernandez, Profumo) for the ρφ ∝ a^{-(4+n)} parameterization and the BBN constraint on Tr; this is a standard phenomenological ansatz used to set up the modified-expansion scenario, it is not a uniqueness theorem, and it does not carry the main derivation. Accordingly, the paper is self-contained and the circularity score is at the low end.

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

The central claim rests on standard tunneling physics plus benchmark-specific assumptions. The free parameters are the collapse threshold and the modified-expansion reference temperature; the main unproven premise is the largeness and near-constancy of N(tpbh). The φ component is an invented phenomenological ingredient for the expansion-rate study.

free parameters (2)
  • δc = 0.45
    Critical density contrast for PBH collapse, adopted from the literature. The normalization of f_pbh depends on this choice; footnote 1 acknowledges the criterion is disputed.
  • Tr = 50 MeV
    Reference temperature where the extra component φ has the same energy density as radiation. Chosen by hand, constrained only by BBN to be above O(10) MeV; affects the modified-expansion results.
assumptions (5)
  • standard math Nucleation rate Γ(T) = T^4 (S3/(2πT))^(3/2) exp(-S3/T) (Eq. 2).
    Standard semiclassical bubble nucleation rate (Coleman, Callan, De Luccia); used as an input to the no-nucleation probability.
  • ad hoc to paper A delayed-decay patch collapses into a PBH when δ(t)=ρ_del/ρ_nor - 1 exceeds δc=0.45 (Eq. 6).
    Adopted from previous literature but flagged as disputed in footnote 1; the quantitative f_pbh values depend on this criterion.
  • domain assumption N(tpbh) is O(10^85) and nearly independent of model parameters (Table II).
    Supported numerically for five benchmarks but not proven; the generality argument for super-exponential dominance rests on it.
  • domain assumption The extra component φ has ρφ ~ a^{-(4+n)} and λ(T)=4+n T^n/(T_r^n+T^n) (Eqs. 23 to 26).
    Phenomenological parametrization for superfast expansion, not derived from a specific particle physics model.
  • standard math The S3(T)/T integral is dominated by the saddle point at Tp, so Gaussian integration with extended limits applies (Eqs. 15 to 18).
    Standard saddle-point approximation; validated numerically by δ(Tp) < O(0.01) in Table II.
invented entities (1)
  • Energy density component φ with ρφ ~ a^{-(4+n)}
    purpose: Model an early-universe superfast expansion rate to test the robustness of the super-exponential PBH abundance scaling.
    Introduced as a phenomenological component with no specified particle physics realization; only constrained by Tr > O(10) MeV from BBN.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Super-exponential Primordial Black Hole Production via Delayed Vacuum Decay." pith.science (2026). https://pith.science/paper/MJL445VH

@misc{pith2026241210666,
  author       = {Pith},
  title        = {Pith review of: Super-exponential Primordial Black Hole Production via Delayed Vacuum Decay},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MJL445VH}},
  note         = {Machine review of arXiv:2412.10666}
}
read the original abstract

If a cosmological first-order phase transition occurs sufficiently slowly, delayed vacuum decay may lead to the formation of primordial black holes. Here we consider a simple model as a case study of how the abundance of the produced black holes depends on the model's input parameters. We demonstrate, using both numerical and analytical arguments and methods, that the black hole abundance is controlled by a double, ``super''-exponential dependence on the three-dimensional Euclidean action over temperature at its minimal value. We show that a modified expansion rate during the phase transition, such as one driven by an additional energy density component, leads to a weaker dependence on the underlying model parameters, but maintains the same super-exponential structure. We argue that our findings generalize to any framework of black hole production via delayed vacuum decay.

Figures

Figures reproduced from arXiv: 2412.10666 by the authors.

Figure 1
Figure 1. fpbh with µ3 under BMs given in Table.I. Upper figure: fpbh with µ3 in BM1; Bottom figure: fpbh with µ3/µ∗ 3 in five BMs, where µ ∗ 3 is chosen when fpbh = 1 within each BM. that because of the extreme sensitivity of the abundance of PBH on the parameters in the effective potential, for convenience we define BM2-5 off of the values of the ref￾erence BM1. We show our numerical results for the abundance fpbh as a func… view at source ↗
Figure 2
Figure 2. Approximation of the integral of Peint, given in Eq.(14). Tp denotes the temperature when Se3(T) reaches its minimum. The grey shadow region represents the full integral from Td to Tc, while the red shadow region denotes the ap￾proximated integral from [Ta, Tb], which is centered around Tp with a typical full width at half maximum of the integrand, denoted by D. we approximate the original integral range (from Td to… view at source ↗
Figure 3
Figure 3. fpbh with S3(Tp)/Tp, where Tp is the tempera￾ture when S3(Tp)/Tp reaches its minimum value. The solid blue line refers to our analytic approximation of fpbh, given in Eq.(22), with M = 4 × 107 and Q = 1 × 1077. The five BMs corresponding to the our numerical values of fpbh of model parameters given in Table.I, where µ ∗ 3 denotes the value of µ3 with fpbh = 1. For each BM, we first shift µ3 away from µ ∗ 3, thereby … view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: fpbh with µ3 under mϕ = 300 MeV, ω = 860 MeV, and c = 0.14 for different cosmological scenarios. Eq.(19) that receives large modification is Hdel(Tp), which increases after introducing the ϕ component; larger values of n naturally lead to larger values of Hdel(Tp). We …

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

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

  1. 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

19 extracted references · 6 canonical work pages · cited by 1 Pith paper

  1. [1]

    Escrivà, F

    A. Escrivà, F. Kuhnel, and Y. Tada, (2022), 10.1016/B978-0-32-395636-9.00012-8, arXiv:2211.05767 [astro-ph.CO]

    arXiv 2022

  2. [2]

    B. Carr, K. Kohri, Y. Sendouda, and J. Yokoyama, Rept. Prog. Phys.84, 116902 (2021), arXiv:2002.12778 [astro- ph.CO]

    arXiv 2021

  3. [3]

    J. Liu, L. Bian, R.-G. Cai, Z.-K. Guo, and S.-J. Wang, Phys. Rev. D 105, L021303 (2022), arXiv:2106.05637 [astro-ph.CO]

    arXiv 2022

  4. [4]

    Pbh formation from overdensities in delayed vacuum transitions,

    K. Kawana, T. Kim, and P. Lu, “Pbh formation from overdensities in delayed vacuum transitions,” (2022)

    work page 2022

  5. [5]

    Revisiting formation of primordial black holes in a supercooled first- order phase transition,

    M. M. Flores, A. Kusenko, and M. Sasaki, “Revisiting formation of primordial black holes in a supercooled first- order phase transition,” (2024)

    work page 2024

  6. [6]

    Primordial black hole for- mation during first-order phase transitions,

    K. Jedamzik and J. Niemeyer, “Primordial black hole for- mation during first-order phase transitions,” (1999)

    work page 1999

  7. [7]

    M. Y. Khlopov, R. V. Konoplich, S. G. Rubin, and A. S. Sakharov, (1998), arXiv:hep-ph/9807343

    arXiv 1998

  8. [8]

    Gonçalves, A

    D. Gonçalves, A. Kaladharan, and Y. Wu, Phys. Rev. D 111, 035009 (2025), arXiv:2406.07622 [hep-ph]

    arXiv 2025

Show all 19 references
  1. [9]

    Connecting the higgs potential and primordial black holes,

    D. Dai, R. Gregory, and D. Stojkovic, “Connecting the higgs potential and primordial black holes,” (2019). 9

    2019

  2. [10]

    X. Wang, F. P. Huang, and X. Zhang, JCAP05, 045 (2020), arXiv:2003.08892 [hep-ph]

    2020 arXiv

  3. [11]

    M. Y. Khlopov, R. V. Konoplich, S. G. Rubin, and A. S. Sakharov, Grav. Cosmol.2, S1 (1999), arXiv:hep- ph/9912422

    1999

  4. [12]

    M. J. Baker, M. Breitbach, J. Kopp, and L. Mittnacht, Phys. Rev. D 111, 063544 (2025), arXiv:2110.00005 [astro-ph.CO]

    2025 arXiv

  5. [13]

    Kanemura, M

    S. Kanemura, M. Tanaka, and K.-P. Xie, JHEP06, 036 (2024), arXiv:2404.00646 [hep-ph]

    2024 arXiv

  6. [14]

    C. G. Callan, Jr. and S. R. Coleman, Phys. Rev. D16, 1762 (1977)

    1977

  7. [15]

    S. R. Coleman and F. De Luccia, Phys. Rev. D21, 3305 (1980)

    1980

  8. [16]

    strength

    clarifies that the original critical collapse criterion still ap- plies, as in inflationary PBH production scenarios. We postpone a full discussion of this issue to future work. Table I. Details of the five BMs. In all BMs, we takemϕ = 300 MeV, and the µ∗ 3 is chosen for fpbh ...

  9. [17]

    Cai, Y.-S

    R.-G. Cai, Y.-S. Hao, and S.-J. Wang, Sci. China Phys. Mech. Astron.67, 290411 (2024), arXiv:2404.06506 [astro-ph.CO]

    2024 arXiv

  10. [18]

    Navaset al

    S. Navaset al. (Particle Data Group), Phys. Rev. D110, 030001 (2024)

    2024

  11. [19]

    D’Eramo, N

    F. D’Eramo, N. Fernandez, and S. Profumo, JCAP05, 012 (2017), arXiv:1703.04793 [hep-ph]. Appendix: Details of the simplified model The toy model we consider consists of a single real scalar field, with an effective potential at finite tempera- ture V (ϕ, T) = 1 2 µ3ω − m2 ϕ 2 ...

    2017 arXiv

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.