{"id":"fe325877-111a-4189-b06f-855df0a6bf87","arxiv_id":"2412.06889","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"Supercooled radiative symmetry breaking phase transitions generically produce primordial black holes, and the false-vacuum decay rate grows exponentially with time to high accuracy.","lead":"This paper calculates how many primordial black holes can form during supercooled phase transitions with radiative symmetry breaking. It finds that black hole production is generic across a broad range of models and can even explain recent microlensing anomalies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Sec. 3.1 proof of exponential Γ(t) is quantitatively incomplete: the second-order correction to Eq. (2.4) is amplified by S3/T at nucleation, producing an O(e^-1.3) suppression at tmax for β/Hn≈9, so the f_PBH maps inherit an unquantified error.","rationale":"I read the paper as a serious, largely self-contained derivation of PBH abundance, mass, and spin in supercooled RSB phase transitions. The late-blooming mechanism, the model-independent parametrization, and the constraint maps are valuable, and the derivation of Eq. (3.15) for β/H_n is clear. The reader correctly identified the homogeneous-patch collapse mapping as a weakness; that limitation is acknowledged in the text before Eq. (2.14). I focus instead on the exponential-growth claim because it is the paper's advertised new result and because it enters everywhere in Sec. 2 through Eq. (2.4). The back-of-the-envelope estimate in my attack shows that the second-order term in the expansion, while formally suppressed by (H_IΔt/X), is multiplied by the large factor a'/X=S3/T at nucleation. For β/H_n≈9 and the time differences actually used in Figs. 1-3, this yields a ~0.27 multiplicative correction to Γ at t_max. That is not a negligible correction to an exponentiated quantity like Pcoll. The paper's statement that the approximation becomes \"increasingly accurate\" with supercooling is plausible asymptotically, but the quantitative statement \"to a high degree of accuracy\" is not supported by the displayed material. This does not invalidate the framework; it means the ready-to-use formulas should be accompanied by a demonstrated numerical comparison, as the reader's CONDITIONAL verdict already requires. My recommendation therefore leaves the verdict unchanged: CONDITIONAL, pending the concrete test above.","tokens_in":22794,"tokens_out":20497,"duration_ms":217128,"concrete_test":"Recompute the late-blooming dynamics using the exact Γ(t) from Eq. (3.16) with S3/T from Eq. (3.11) (supercool LO) and Eq. (3.12) (improved) instead of the linearized approximation in Eq. (A.8). For a grid of the model-independent parameters used in Figs. 4-10, including β/H_n≈5 and 9 and the HSC point, solve the Appendix A ODE system, recompute Pcoll from Eqs. (2.18)-(2.20), and compare the resulting f_PBH with the published maps. If Pcoll changes by more than a factor of ~2 in the β/H_n≈5-10 region, the Sec. 3.2 formulas need either a corrected β(t) or a quantitative error budget.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central novel step is the justification of Eq. (2.4) in Sec. 3.1. The proof expands 1/(X+H_I(t-t_n)) in Eq. (3.22) and keeps only the linear term. The omitted second-order term in the exponent of Γ is -(a'/X)(H_I(t-t_n)/X)^2, with a'/X = S3/T at t_n. This prefactor is large. For the paper's own example with β/H_n=9 and τ_max-τ_n≈2.08, using H_I(t_max-t_n)=γ(τ_max-τ_n)/√3≈0.916, the RSB relation a'/X=4 ln(T_n/H_n), and β/H_n=a'/X^2-4, one finds X≈8.6 and a'/X≈112 for representative parameters. The quadratic correction is then (a'/X)(0.916/8.6)^2≈1.3, suppressing Γ by e^-1.3≈0.27 relative to Eq. (2.4) at t_max. This is not \"very small\"; the correction is about 15% of the linear exponent β(t-t_n) but is exponentiated in Pcoll via Eqs. (2.18)-(2.21), so f_PBH can shift by an order of magnitude or more. The paper states that a numerical check was performed for the improved supercool expansion but displays no residuals, so this quantitative gap is unresolved and is distinct from the collapse-threshold issue the reader emphasized.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies primordial black hole (PBH) formation via the late-blooming mechanism in supercooled first-order phase transitions with radiative symmetry breaking (RSB). It reformulates the mechanism in Sec. 2, assumes an exponentially growing false-vacuum decay rate Gamma(t), and derives fitting formulas for the collapse probability, PBH abundance, mass, and initial spin. In Sec. 3 it attempts to justify the exponential time dependence of Gamma using the large-supercooling expansion of RSB, and it scans the model-independent parameters chi0, beta-bar, g, and g-tilde to map f_PBH, M_PBH, and the initial spin, including observational constraints. As an application, it identifies a region of a gauged B-L extension of the Standard Model that can fit recently reported HSC microlensing anomalies.","tokens_in":23024,"tokens_out":11452,"duration_ms":120873,"significance":"If the central claims hold, the paper provides a useful model-independent package: for any RSB model with large supercooling one could read off PBH abundance, mass, and spin, and compare with constraints, with a concrete B-L model as a test case. The paper is transparent about some limitations, notably the homogeneous-patch approximation in Sec. 2, and it cites a broad set of PBH constraints. Its main novelty is the attempted justification of Eq. (2.4) and the systematic parameter scan, which extends earlier work in Refs. [20,22,33]. However, the exponential-decay proof is an order-of-magnitude estimate whose controlling second-order term is not small for the parameter range used, and the collapse probability inherits unquantified uncertainties from a fitted function and the collapse threshold. The qualitative conclusion that PBHs are generically produced may survive, but the quantitative predictions do not yet match the strength of the claims made in the abstract and conclusions.","major_comments":[{"comment":"The proof that the non-linear terms in Eq. (3.22) are negligible is quantitatively incomplete. The second-order term in the exponent of Gamma is -(a'/X)(H_I(t-t_n)/X)^2. For the paper's own example with beta/H_n = 9 and tau_max - tau_n = 2.08, one has H_I t_eq = gamma/sqrt(3) = 0.4406, so H_I(t_max-t_n) = 0.916. For representative RSB parameters consistent with Eq. (3.15) and the nucleation condition a'/X = 4 ln(T_n/H_n), one finds X = 8.6 and a'/X = 112, giving a correction at t_max of approximately -1.3 in the exponent. This suppresses Gamma by e^{-1.3} = 0.27 relative to Eq. (2.4), which is about 15% of the linear exponent beta(t_max-t_n). Because this correction enters P_coll through Eqs. (2.18)-(2.21), f_PBH can shift by an order of magnitude or more. The claim that the corrections are 'always very small' is therefore not established by the displayed order-of-magnitude estimate, and the numerical check in the improved supercool expansion is not shown. Please provide residuals or a controlled upper bound on the second- and higher-order terms over the parameter range used in the figures.","section":"Sec. 3.1, Eq. (3.22)"},{"comment":"The collapse criterion is set by requiring the space-averaged homogeneous mass excess delta to reach delta_c = 0.45. The paper candidly states after Eq. (2.13) that 'after t_max one must start including inhomogeneities' and that curvature perturbations cannot be captured. That limitation is load-bearing: Eq. (2.14) fixes the late-blooming time t_PBH_ni through this homogeneous-space criterion, and Eq. (2.18) then converts t_PBH_ni into a collapse probability. Neither the threshold delta_c nor the averaging scale is derived within the mechanism. Consequently, the absolute normalization of f_PBH and the positions of the f_PBH = 1 and constraint contours in Figs. 4-10 carry a systematic uncertainty that is not displayed. A sensitivity analysis over delta_c and over the patch-size choice (for example, Hubble radius versus sound horizon) should be presented, or the central claims should be restricted to relative comparisons between parameter regions.","section":"Sec. 2, Eqs. (2.13)-(2.14)"},{"comment":"Equation (2.21) is a fitting function whose coefficients a_P, b_P, c_P are quoted to four significant figures, but no residuals, fit range, or uncertainties are reported. Since P_coll enters f_PBH exponentially through Eq. (2.22), even a few percent error in b_P or c_P can change f_PBH by orders of magnitude in parts of the parameter space shown in Figs. 5-10. The text says the fit is obtained 'varying delta_c around 0.45,' but it is not clear whether the quoted coefficients are refitted for each delta_c or whether the delta_c dependence is captured by the final factor. Please show the fit residuals, state the validity range of Eq. (2.21), and provide a sensitivity estimate. Without this, the advertised 'ready-to-use' formulas do not have a stated quantitative reliability.","section":"Sec. 2, Eq. (2.21)"}],"minor_comments":[{"comment":"The line defining tau_ni contains a typo: 'tau_ni = t_ni/t_ni' should read 'tau_ni = t_ni/t_eq', and similarly for tau_PBH_ni.","section":"Appendix A, after Eq. (A.5)"},{"comment":"The captions use the notation 'q <a_*^(1/2)>'; according to Eq. (2.23) the plotted quantity is sqrt(<a_*^2>), the RMS dimensionless Kerr parameter. Please make the notation consistent.","section":"Captions of Figs. 5, 6, 8, 9"},{"comment":"The phrase 'estimate of thedelta, which includes non-linear effects' contains a typo ('thedelta') and should read 'estimate of the delta'.","section":"Sec. 2, after Eq. (2.13)"},{"comment":"The abstract and conclusions say the paper 'demonstrates' and provides a 'full justification' of the exponential time dependence of Gamma, while Sec. 3.1 itself describes the estimate as an order-of-magnitude argument. The wording should be softened to match the strength of the evidence actually presented.","section":"Abstract and Sec. 4"},{"comment":"The coefficients beta_2 and beta_3 in Eq. (2.2) are not explicitly defined with the factorial factors of a Taylor expansion. A brief defining sentence would avoid ambiguity.","section":"Eq. (2.2)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about its own limitations, which is to its credit, but the quantitative gaps in Sec. 3.1 and in Eq. (2.21) are central to the advertised ready-to-use predictions. The HSC/B-L application is suggestive and should be retained if the numerical issues are fixed. I see no concerns about attribution; earlier work by Salvio and by Gouttenoire and Volansky is properly cited."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this paper is worth reading and worth sending to a referee, but the headline claim — that the false-vacuum decay rate grows exponentially with time to high accuracy in RSB — is only backed by an order-of-magnitude estimate, and for the paper's own benchmark the next-order correction is not tiny. The abundance maps inherit an unquantified error.\n\nWhat's genuinely new: they attempt to justify the exponential Γ(t) that everyone assumes, and they package the RSB model-independent scan into ready-to-use formulas for f_PBH, M_PBH, and initial spin, including a concrete B-L model fit to the HSC anomalies. The parameter-space plots with constraints are useful, and the paper is clearly written.\n\nThe soft spots: the Sec. 3.1 argument expands 1/(X+H_I Δt) and keeps only the linear term. The paper says the corrections are \"very small,\" but a quick estimate for their β/Hn=9 example gives a quadratic term (a'/X)(H_I Δt/X)^2 ≈ 1.3 at t_max, suppressing Γ by a factor ~0.3. That's not negligible once exponentiated into P_coll. The numerical check is mentioned but no residuals are shown. This doesn't kill the mechanism, but the claimed \"high accuracy\" is not demonstrated. The homogeneous-patch collapse model is a separate limitation, and the authors acknowledge it. Also, the HSC fit is an application with chosen parameters, not a sharp prediction.\n\nThe citation pattern looks fine; prior work by Gouttenoire, Salvio, and others is properly credited. No sign of cooked numbers — the issue is incomplete error quantification, not dishonesty.\n\nWho this is for: people working on PBH production mechanisms and phase-transition cosmology. If I were refereeing it, I'd ask for a quantitative bound on the exponential approximation and a few residual plots. The framework is solid enough to deserve referee time.\n\nRecommendation: send to peer review with major revisions. The core results are useful, but the central accuracy claim needs to be backed by numbers.","headline":"Useful model-independent PBH formulas from RSB phase transitions, but the proof of exponential Γ(t) is backed by an order-of-magnitude argument with an unquantified O(1) correction in the paper's own benchmark.","tokens_in":23714,"tokens_out":5801,"would_cite":true,"duration_ms":54348,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper claims that radiative symmetry breaking with strong supercooling generically produces primordial black holes, with model-independent predictions for their abundance, mass, and initial spin.","keywords":["primordial black holes","radiative symmetry breaking","supercooled phase transition","late-blooming mechanism","false vacuum decay","dark matter abundance","microlensing anomalies","B-L extension"],"falsifier":"A full numerical-relativity simulation of the collapse of a late-blooming Hubble patch in a supercooled RSB transition, computing the compaction-function threshold directly, would settle the central claim: if that threshold departs significantly from $\\delta_c=0.45$, the predicted $f_{\\rm PBH}$ moves by orders of magnitude and the claimed generic-production region in the parameter plane would not hold.","tokens_in":22462,"feed_emoji":"🕳️","tokens_out":12029,"duration_ms":113828,"temperature":0.7,"pith_summary":"The paper works to establish that primordial black hole (PBH) production is a generic, model-independent consequence of radiative symmetry breaking (RSB) whenever the associated first-order phase transition supercools strongly. It shows that in such theories the false-vacuum decay rate grows exponentially with time, $\\Gamma(t)\\approx H_n^4 e^{\\beta(t-t_n)}$, to high accuracy, and that this exponential law makes the late-blooming mechanism the dominant channel for PBH formation. It then gives ready-to-use formulas for the PBH dark matter fraction $f_{\\rm PBH}$, the typical mass $M_{\\rm PBH}$, and the initial spin in terms of the few supercool-expansion parameters, and it maps the region of parameter space where PBHs are produced with observable abundance. A concrete Standard Model extension with gauged $B-L$ and right-handed neutrinos is shown to reproduce reported microlensing anomalies, so the mechanism is testable. A sympathetic reader would care because the mechanism demands no fine-tuning and because $f_{\\rm PBH}\\le1$ becomes a model-independent bound on the parameter space of RSB theories.","feed_headline":"Radiative symmetry breaking generically creates primordial black holes","feed_subtitle":"Ready-to-use formulas predict abundance, mass and spin; one model matches microlensing anomalies","key_machinery":"The load-bearing object is the late-blooming mechanism combined with the supercool expansion of RSB theories. The late-blooming mechanism treats a region that nucleates late as a separate homogeneous patch with its own scale factor $a_l(t)$ and Hubble rate $H_l(t)$, while the background evolves with $a_b(t)$ and $H_b(t)$; the false-vacuum fraction $F(t,t_{ni})=e^{-I(t,t_{ni})}$, with $I$ given by Eq. (2.10), controls how vacuum energy converts into radiation. The supercool expansion rewrites the effective potential using the parameters $\\chi_0$, $\\bar\\beta$, $g$, and $\\tilde g$, so that the nucleation temperature and $\\beta/H_n$ follow from closed formulas such as Eqs. (3.13)-(3.15). The paper's new step is to prove, from this expansion, that $\\Gamma(t)\\approx H_n^4 e^{\\beta(t-t_n)}$ is accurate on the time intervals that matter; that exponential law is what reduces $P_{\\rm coll}$, $f_{\\rm PBH}$, $M_{\\rm PBH}$, and $\\sqrt{\\langle a_*^2\\rangle}$ to simple functions of $\\beta/H_n$, $\\delta_c$, and $T_{\\rm eq}$, and it also explains why smaller $\\beta/H_n$ makes PBH production more efficient.","core_discovery":"On the paper's own terms, the central claim is that PBH production is generic across the model-independent parameter space of perturbative RSB theories with large supercooling, operating through the late-blooming mechanism: some Hubble patches remain in the false vacuum longer than the average background, their radiation density is diluted less, and the resulting mass excess makes them collapse into black holes. The demonstration has two linked parts. First, in the supercool and improved supercool expansions the decay rate $\\Gamma$ is dominated by the time-independent bounce, and the nearly de Sitter expansion before nucleation gives $\\Gamma(t)\\approx H_n^4 e^{\\beta(t-t_n)}$ with corrections suppressed by the supercooling parameter $X=\\log(\\chi_0/T_n)$. Second, solving the two-patch Friedmann equations with this $\\Gamma$ yields the mass excess $\\delta(t,t_{ni})$; when its maximum reaches the threshold $\\delta_c=0.45$ the patch collapses, and the resulting collapse probability gives $f_{\\rm PBH}$, $M_{\\rm PBH}$, and the RMS initial spin $\\sqrt{\\langle a_*^2\\rangle}$ as explicit functions of the model-independent parameters. The paper finds that the PBH abundance scans a broad range up to $f_{\\rm PBH}>1$ over the parameter plane, that the requirement $f_{\\rm PBH}\\le1$ therefore constrains RSB theories whenever the PBHs survive until today, and that a $U(1)_{B-L}$ extension with right-handed neutrinos can fit the HSC microlensing anomaly at $\\chi_0\\sim2\\times10^4$ GeV.","pith_inferences":["Not stated in the paper: if the exponential decay law is accurate, the same parameters that set $f_{\\rm PBH}$ also fix the stochastic gravitational-wave background from the transition, so a future detection of a strongly supercooled RSB transition carries a definite PBH-abundance prediction to cross-check.","Not stated in the paper: the paper's spin result is only the RMS value; applying peak statistics to the shear field would yield the full spin distribution, which binary-merger data could test.","Not stated in the paper: the threshold $\\delta_c=0.45$ is imported from spherical-collapse studies; non-linear collapse simulations of late-blooming patches would sharpen it and could shift the boundaries of the claimed production region.","Not stated in the paper: the same formulas can be applied to other supercooling sectors (axion, dark scalar, or otherwise) to decide whether they produce PBHs without repeating the two-patch integration."],"forward_implications":["Any RSB model with large supercooling can be tested for PBH dark matter without a dedicated simulation: one plugs the model's $(\\chi_0,\\bar\\beta,g)$ into the provided formulas for $f_{\\rm PBH}$, $M_{\\rm PBH}$, and spin.","Requiring $f_{\\rm PBH}\\le1$ yields a model-independent upper bound on the RSB parameter space for $\\chi_0\\lesssim10^9$ GeV, where the produced PBHs survive Hawking evaporation until today.","The predicted mass is concentrated around a value set by $\\chi_0$ through $M_{\\rm PBH}\\approx M_J (2\\pi^2/3\\bar\\beta)^{1/2}(280\\,{\\rm MeV}/\\chi_0)^2$, so PBH searches translate directly into constraints on the symmetry-breaking scale.","The initial PBH spin is predicted to be very small, and the paper notes that accretion, mergers, close hyperbolic encounters, and scalar-driven Hawking evaporation can later spin PBHs up, so a large observed spin would not contradict the mechanism's initial conditions.","A concrete Standard Model extension with gauged $B-L$ and right-handed neutrinos is shown to fit the HSC microlensing anomalies for $\\chi_0\\sim2\\times10^4$ GeV and $g\\simeq0.95$, making the mechanism testable with optical surveys."],"supporting_citations":[{"why":"Provides the model-independent supercooling framework and the concrete $U(1)_{B-L}$ model whose parameter region is tested here.","marker":"[22]"},{"why":"Establishes the model-independent description of RSB phase transitions and the gravitational-wave constraints that the paper relies on.","marker":"[33]"},{"why":"Gives the late-blooming mechanism, the fit for $P_{\\rm coll}$, and the $f_{\\rm PBH}$ and $M_{\\rm PBH}$ formulas that this work re-derives and makes model-independent.","marker":"[20]"},{"why":"Lays the foundation of radiative symmetry breaking as the origin of spontaneous symmetry breaking, the mechanism under study.","marker":"[38]"},{"why":"Extends radiative symmetry breaking to general flat directions of the scalar potential, justifying the parameterization used in the supercool expansion.","marker":"[40]"},{"why":"Supplies the survival probability of a Hubble patch against bubble nucleation, used in Eq. (2.18).","marker":"[7]"},{"why":"Determines the false-vacuum fraction $F(t,t_{ni})$ that converts vacuum energy to radiation in the two-patch equations.","marker":"[54,55]"},{"why":"Provides the mass-excess threshold for gravitational collapse, the value $\\delta_c=0.45$ adopted for the late-blooming patch.","marker":"[52]"},{"why":"Supplies the RMS initial-spin formula for PBHs from an early phase transition, combined here with the RSB mass and abundance.","marker":"[58]"}],"fun_headline_variants":["Supercooled RSB generically produces black holes","Ready-to-use PBH predictions from radiative symmetry breaking","RSB supercooling explains microlensing anomaly via PBHs","Generic PBH mass, spin, and abundance from supercooled RSB","Late-blooming vacuum patches become black holes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation treats the late-blooming patch as a homogeneous region and assumes it collapses into a black hole when its space-averaged mass excess $\\delta$ reaches $\\delta_c=0.45$, a threshold taken from spherical-collapse studies; the true collapse threshold and the inhomogeneities that develop after $t_{\\max}$ are not computed, so the predicted $f_{\\rm PBH}$ shifts if that threshold is inaccurate.","fun_headline_variants_meta":{"raw":{"variants":["Supercooled RSB generically produces black holes","Ready-to-use PBH predictions from radiative symmetry breaking","RSB supercooling explains microlensing anomaly via PBHs","Generic PBH mass, spin, and abundance from supercooled RSB","Late-blooming vacuum patches become black holes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00116,"raw_usage":{"total_tokens":4862,"prompt_tokens":1065,"completion_tokens":3797,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":681,"completion_tokens_details":{"reasoning_tokens":3714}},"tokens_in":681,"tokens_out":3797,"duration_ms":25277,"temperature":1.0,"reasoning_tokens":3714,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T19:22:08.059784+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A full numerical-relativity simulation of the collapse of a late-blooming Hubble patch in a supercooled RSB transition, computing the compaction-function threshold directly, would settle the central claim: if that threshold departs significantly from $\\delta_c=0.45$, the predicted $f_{\\rm PBH}$ moves by orders of magnitude and the claimed generic-production region in the parameter plane would not hold.","supporting_citations":[],"review_version":1}