{"id":"a085316b-f576-4293-87d8-b64bc246f9c8","arxiv_id":"2501.11040","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Delayed first-order electroweak phase transitions can form super-critical primordial black holes, and a timescale ratio t_H/t_V captures the threshold better than the standard density contrast.","lead":"Using the equations for a thin bubble wall, this paper shows numerically that a delayed electroweak phase transition can leave vacuum regions that expand forever and look like black holes to outside observers. The authors propose a timescale-based formation criterion and update the model regions that black-hole searches can test.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The claim that the timescale criterion is more appropriate than the density-contrast criterion is tested against δ(t_H), not δ_max, so the paper's headline comparison is not yet established.","rationale":"The central claim has two parts: (1) super-critical PBHs can form under the stated assumptions, and (2) the timescale criterion t_V <~ t_H is superior to δ > δ_C. The δ_max issue attacks (2) directly and is an internal inconsistency: the authors themselves identify δ_max > δ_C as the conventional criterion in Sec. III B, then test δ(t_H) in Fig. 2. Since δ is time-dependent and grows during radiation domination until vacuum energy dominates, δ(t_H) is not the quantity the conventional criterion uses. The paper does not compute δ_max anywhere, so its conclusion about appropriateness rests on comparing the timescale to a suboptimal version of the density criterion. This is more load-bearing than the no-friction assumption (A18), which, while physically simplifying, is stated explicitly and affects both criteria equally through the EoM. Spherical symmetry is also acknowledged as a limitation in Sec. V. The missing δ_max baseline is a testable, self-contained flaw that can be resolved from the paper's own numerical setup. I therefore maintain the reader's CONDITIONAL verdict: the formation mechanism is plausible and numerically supported under the stated assumptions, but the headline criterion claim should be verified against δ_max before being accepted.","tokens_in":18001,"tokens_out":18577,"duration_ms":200145,"concrete_test":"For the two benchmark points in Fig. 1, locate the sub/super-critical threshold χ_in/χ_h by bisection in the EoM (16), then compute δ(t) from Eq. (25) throughout the evolution and record δ_max = max_t δ(t) up to the classification time t_end. Repeat for the Λ range and v_w values scanned in Fig. 2. If δ_max at threshold is approximately flat (within, say, 10%) and close to δ_C ≈ 0.45, the conclusion that Eq. (24) is more appropriate than the conventional criterion is not supported; if δ_max varies as much as δ(t_H), the timescale criterion is genuinely more robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central new claim—that the timescale criterion t_V <~ t_H is 'more appropriate' than the conventional density-contrast criterion—is supported in Fig. 2 by comparing t_H/t_V with δ(t_H), not with δ_max. But Sec. III B defines the conventional criterion as δ_max > δ_C (ref. [19]) and notes that δ depends on the time at which it is evaluated. Since δ grows as the background radiation density decays while ρ_V stays constant, δ(t_H) can be systematically smaller than δ_max; applying a fixed δ_C to δ(t_H) will under-predict super-criticality, making the timescale criterion look better than it is. The paper's own ground truth is the EoM threshold from Eq. (16), so the correct test is to evaluate δ_max for the same threshold solutions. If δ_max turns out to be approximately constant (~0.45) across the scanned (Λ, v_w) values, the conventional criterion would work just as well and the headline claim would be undermined. This is an internal mismatch between the stated criterion and the tested quantity, not a disagreement with external consensus. The no-friction junction condition (A18) is a further simplification that changes the EoM, but the δ_max comparison is the most direct threat to the paper's distinctive conclusion.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies super-critical primordial black hole (PBH) formation from delayed first-order electroweak phase transitions. It derives an ordinary differential equation for the comoving area radius of a spherically symmetric false vacuum domain boundary from Israel's junction conditions under the thin-wall and no-friction approximations (Eq. (16), Appendix A), solves it numerically for benchmarks in the nearly-aligned Higgs effective field theory (naHEFT), and identifies supercritical versus subcritical solutions. The paper then proposes the timescale condition t_V ≲ t_H (Eq. (24)) as a criterion that is more appropriate than the conventional density-contrast criterion δ > δ_C, and uses the numerical solutions to estimate PBH masses and abundances and to map naHEFT parameter regions probed by microlensing, gravitational-wave, and collider experiments.","tokens_in":18327,"tokens_out":5290,"duration_ms":61099,"significance":"If the central comparison is established, this work would provide a quantitative validation of the timescale criterion for super-critical PBH formation and would strengthen the case for using PBH observations to probe supercooled electroweak phase transitions. The paper has clear strengths: the derivation of Eq. (16) is explicit and self-contained; the benchmark scans are internally consistent; the final parameter regions agree with earlier work (Ref. [23]), which is a useful cross-check; and the authors explicitly acknowledge the spherical-symmetry and thin-wall limitations in Sec. V. However, as detailed below, the headline claim that t_H/t_V is more appropriate than the conventional density-contrast criterion is not yet supported because the comparison uses δ(t_H) rather than δ_max.","major_comments":[{"comment":"The paper defines the conventional criterion as δ_max > δ_C and explicitly notes that δ depends on the time at which it is evaluated, but the numerical test in Fig. 2 compares t_H/t_V with δ(t_H), not with δ_max. During supercooling, ρ_R decays while ρ_V remains constant, so δ(t_H) is systematically smaller than δ_max; applying a fixed δ_C to δ(t_H) biases the comparison in favor of the timescale criterion. The authors should compute δ_max for the same threshold solutions obtained from Eq. (16) and test whether δ_max is approximately constant (~0.45) across the scanned (Λ, v_w) values. If it is approximately constant, the conventional criterion would perform comparably and the paper's central claim would be undermined; if not, the reasons should be stated explicitly.","section":"Sec. III B/C, Eq. (25), Fig. 2"},{"comment":"The no-friction junction condition [u^μ ξ_μ]_0 = 0 forces [ρ_R]_0 = 0 and makes the density contrast come entirely from vacuum energy. For a supercooled electroweak phase transition with fast-moving walls, wall-plasma friction is not obviously negligible, and v_w is treated as a free parameter in Sec. II rather than being derived from the same wall dynamics. Since Eq. (16) and the extracted t_H/t_V threshold depend on ρ_V and σ, the authors should either justify this assumption for the parameter range considered or quantify the sensitivity of the threshold to wall-fluid friction. Without this, the quantitative validity of Eq. (24) for realistic electroweak phase transitions remains open.","section":"Appendix A, Eq. (A18), Eq. (16)"}],"minor_comments":[{"comment":"Equation (A21) contains \"H ˙r\" where Eq. (16) has \"H ˙χ\"; this appears to be a typographical error and should be corrected.","section":"Appendix A, Eq. (A21)"},{"comment":"The text states that the gray region corresponds to the requirement for the validity of naHEFT Λ > v, while the Fig. 4 caption says the gray region is where the requirement (Λ < v) is not satisfied; these statements are inconsistent and should be reconciled.","section":"Sec. IV and Fig. 4"},{"comment":"The caption says the red solid and dashed lines \"indicate the requirement of the time ratio t_H/t_V\" but does not state whether these are the threshold values at the supercritical boundary for each wall velocity; the threshold values and their extraction from the numerical solutions should be stated explicitly.","section":"Fig. 2"},{"comment":"The quantity P(χ_in) is described as the probability for a given value of χ_in, but it is not clear whether this is a probability density or a probability for a fixed comoving radius; the authors should clarify the normalization and explain explicitly how χ_th is chosen from the numerical threshold when evaluating f_PBH.","section":"Sec. IV, Eqs. (28)-(29)"},{"comment":"The text notes that the wall velocity v_w and χ̇ are independent variables, but a reader may wonder whether the FVD boundary speed should be related to the bubble wall speed; a brief clarification of the physical distinction would be helpful.","section":"Sec. III A"}],"recommendation":"major_revision","confidential_remarks":"The main unresolved issue is the δ_max comparison in Fig. 2; this is a fixable but load-bearing point because it directly concerns the paper's distinctive conclusion. The no-friction assumption in Appendix A also deserves a dedicated sensitivity test. The use of the EoM as ground truth is acceptable and does not amount to circularity in my view. If the authors can provide the δ_max analysis and a discussion of the no-friction assumption, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this paper does something genuinely useful—it solves the thin-wall false vacuum domain boundary equation numerically for two naHEFT benchmark points and shows that the t_H/t_V timescale tracks the super-critical threshold. As far as I know, that explicit check is new. The derivation of the EoM in Appendix A is clean, the scans are internally consistent, and the authors honestly state that their updated parameter regions are almost unchanged from Ref. [23]. Credit where it is due.\n\nThe soft spot is the headline comparison. The paper claims the timescale criterion is \"more appropriate\" than the conventional density-contrast criterion, but Fig. 2 compares t_H/t_V against δ evaluated at t_H, not against δ_max. The paper itself cites the conventional criterion as δ_max > δ_C (Ref. [19]). Since δ grows as radiation dilutes while vacuum energy stays constant, δ(t_H) is systematically smaller than δ_max, which stacks the deck in favor of the timescale criterion. The right test is to evaluate δ_max for the same super-critical EoM solutions and check whether it is roughly constant around 0.45. If it is, the conventional criterion works just as well and the central methodological claim weakens substantially. As it stands, the claim is not established.\n\nOther concerns are minor relative to that. Only two benchmark points are scanned, so the claimed insensitivity of t_H/t_V to model parameters rests on thin evidence. The no-friction junction condition (A18) forces the density contrast to come entirely from vacuum energy; wall-plasma friction could change the EoM. Spherical asymmetry is acknowledged as potentially making formation easier. None of these are fatal for the numerical result itself; they just limit how much can be concluded.\n\nWho should read this: people working on PBH formation from phase transitions and model-builders in naHEFT. It deserves a serious referee—the numerical EoM solution is a step beyond earlier work—but the \"more appropriate criterion\" sentence should be rewritten or supported with a δ_max comparison before I would trust it. I would send it to review with a request for that comparison.","headline":"Useful numerical check of super-critical PBH formation from delayed EWPTs, but the headline claim that the timescale criterion beats δ > δ_C is tested against δ(t_H), not δ_max, so that claim is not yet established.","tokens_in":18852,"tokens_out":2197,"would_cite":true,"duration_ms":23506,"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":"A delayed electroweak phase transition can leave false-vacuum regions that collapse into super-critical primordial black holes, and the threshold is set by a timescale comparison rather than by the conventional density contrast.","keywords":["primordial black holes","electroweak phase transition","super-critical PBH","false vacuum domains","thin-wall approximation","nearly aligned Higgs effective field theory","bubble wall dynamics","PBH abundance constraints"],"falsifier":"Recompute Eq. (16) with a nonzero fluid-wall coupling, dropping the condition $[u^\\mu\\xi_\\mu]_0=0$, or run a non-spherical numerical-relativity simulation of an ellipsoidal false-vacuum domain of comparable initial radius; if the threshold shifts substantially or super-critical solutions disappear even for $t_V\\lesssim t_H$, the timescale criterion and the formation claim would be falsified.","tokens_in":17810,"feed_emoji":"🕳️","tokens_out":8410,"duration_ms":93521,"temperature":0.7,"pith_summary":"This paper argues that a sufficiently delayed first-order electroweak phase transition can leave isolated regions of the early universe stuck in the unbroken (false) vacuum, and that such regions can end up as primordial black holes of a special 'super-critical' kind. The authors solve the equation of motion for the boundary of such a false-vacuum domain, under a thin-wall, spherically symmetric, frictionless approximation, and find numerically that super-critical solutions exist when the initial domain radius is large enough. They conclude that the appropriate formation criterion is a timescale comparison, $t_V \\lesssim t_H$, rather than the conventional density-contrast bound $\\delta>\\delta_C$. Using the nearly aligned Higgs effective field theory, they update the new-physics parameter regions that PBH observations can probe, obtaining typical super-critical PBH masses near $4\\times10^{-5}M_\\odot$.","feed_headline":"A delayed electroweak phase transition can form black holes","feed_subtitle":"A vacuum-energy timescale, not density contrast, best predicts collapse; new-model search regions are updated.","key_machinery":"The central object is Eq. (16), the equation of motion for the comoving area radius $\\chi$ of the false-vacuum domain boundary. It comes from Israel's junction conditions for a thin shell: the wall has surface energy $\\sigma$, and the matching is simplified by the condition $[u^\\mu\\xi_\\mu]_0=0$, which means the radiation fluid crosses the wall with no friction and therefore has no density jump across the boundary, so the density contrast is sourced entirely by vacuum energy. The same machinery supplies the Misner-Sharp mass used to set the initial comoving horizon radius and to estimate the PBH mass at formation. The dynamics are classified by whether the boundary radius grows (super-critical) or collapses (sub-critical), and the threshold is distilled into the timescale condition $t_V\\lesssim t_H$, where $t_V$ is the vacuum-domination time inside the domain and $t_H$ is the horizon-crossing time.","core_discovery":"On the paper's own terms, the central discovery is that super-critical PBH formation through a delayed electroweak phase transition happens whenever the false-vacuum domain's vacuum energy takes over the local expansion before the domain crosses the horizon. Solving the boundary equation (16) for the comoving area radius, the authors find an initial-radius threshold separating sub-critical solutions, whose boundary radius collapses to zero, from super-critical solutions, whose boundary radius grows without bound and which outside observers see as a baby universe inside a wormhole, i.e. a super-critical PBH. The threshold is numerically far better described by $t_V \\lesssim t_H$ than by the widely used $\\delta>\\delta_C$ condition, because the timescale ratio is almost insensitive to the model parameters and wall velocity while the density fluctuation is not. For their naHEFT benchmarks the resulting PBHs have mass about $4\\times10^{-5}M_\\odot$, and the parameter regions accessible to PBH observations nearly coincide with those found in Ref. [23].","pith_inferences":["If wall-plasma friction is not actually negligible, the no-friction junction condition overstates the driving pressure, so the super-critical threshold would shift toward larger initial domain radii and lower PBH abundances than reported.","The same timescale comparison should apply to other deeply supercooled phase transitions beyond the electroweak one, so PBH bounds could be mapped onto any model with a sufficiently suppressed nucleation rate.","A numerical-relativity study of non-spherical false-vacuum domains would show how much of the favourable threshold comes from the spherical-symmetry assumption; until then, the formation rates should be read as upper estimates.","The near-invariance of the naHEFT parameter regions across formation criteria suggests that PBH observability is controlled mainly by the nucleation rate rather than by the boundary dynamics, so future probes of PBHs would be measuring the phase-transition action more directly than the geometry of collapse."],"forward_implications":["Super-critical PBHs from a delayed EWPT would have masses around $4\\times10^{-5}M_\\odot$, lying in the window probed by HSC, OGLE, EROS and future Roman/PRIME microlensing searches.","The $t_V\\lesssim t_H$ criterion gives model-builders a fast, parameter-insensitive check for whether a proposed phase transition can produce PBHs, without solving the boundary dynamics in full.","Because the naHEFT parameter regions are nearly unchanged when the super-critical dynamics are solved explicitly, earlier EFT-based predictions for PBH observability are corroborated rather than displaced.","PBH abundance is exponentially sensitive to the bubble-nucleation action, so PBH observations would constrain the nucleation rate strongly even though the PBH mass is fixed near the electroweak Hubble scale.","Combining PBH constraints with gravitational-wave spectra (LISA, DECIGO) and triple-Higgs-coupling measurements covers complementary regions of the naHEFT parameter space."],"supporting_citations":[{"why":"introduced super-critical PBHs and the $t_V\\lesssim t_H$ heuristic that the paper verifies numerically","marker":"[8]"},{"why":"classified super-critical versus sub-critical vacuum-domain solutions in the wormhole/PBH picture","marker":"[9]"},{"why":"derived the vacuum-bubble junction-condition EoM that the paper extends to the EWPT setting","marker":"[10]"},{"why":"supplied the boundary EoM and Misner-Sharp mass expressions used in the numerical solutions","marker":"[11]"},{"why":"gave the spherical domain-wall collapse formalism and the sub-critical mass estimate the paper contrasts with super-critical formation","marker":"[7]"},{"why":"Israel's junction conditions are the basis of the thin-wall matching and of Eq. (16)","marker":"[87]"},{"why":"set the original PBH-as-probe-of-EWPT framework and the mass estimate the paper updates","marker":"[20]"},{"why":"provided the earlier naHEFT parameter-region analysis whose results the paper confirms with explicit boundary dynamics","marker":"[23]"},{"why":"supports the neglect of surface-energy contributions and the super-critical formation picture in supercooled transitions","marker":"[29]"},{"why":"CosmoTransitions is used to compute bounce solutions and the surface energy $\\sigma$","marker":"[97]"}],"fun_headline_variants":["Delayed electroweak phase transition births black holes","Timing, not density, sets black hole formation threshold","False-vacuum pockets become black holes in delayed transition","Electroweak delay redefines how black holes form","Vacuum energy timescale predicts black hole collapse"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the unbroken region is perfectly spherical at its boundary and that the boundary wall moves through the surrounding plasma with no friction; if either assumption fails, black hole formation could turn out to be harder than predicted.","fun_headline_variants_meta":{"raw":{"variants":["Delayed electroweak phase transition births black holes","Timing, not density, sets black hole formation threshold","False-vacuum pockets become black holes in delayed transition","Electroweak delay redefines how black holes form","Vacuum energy timescale predicts black hole collapse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000384,"raw_usage":{"total_tokens":2010,"prompt_tokens":901,"completion_tokens":1109,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":517,"completion_tokens_details":{"reasoning_tokens":1032}},"tokens_in":517,"tokens_out":1109,"duration_ms":11753,"temperature":1.0,"reasoning_tokens":1032,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T18:41:57.306236+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute Eq. (16) with a nonzero fluid-wall coupling, dropping the condition $[u^\\mu\\xi_\\mu]_0=0$, or run a non-spherical numerical-relativity simulation of an ellipsoidal false-vacuum domain of comparable initial radius; if the threshold shifts substantially or super-critical solutions disappear even for $t_V\\lesssim t_H$, the timescale criterion and the formation claim would be falsified.","supporting_citations":[],"review_version":1}