Primordial black holes spin from cosmological first-order phase transitions
Pith reviewed 2026-06-28 13:33 UTC · model grok-4.3
The pith
Cosmological first-order phase transitions generate spin in primordial black holes, with the Kerr parameter scaling with latent heat strength and transition rate.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
By integrating the nucleation history without assuming a Gaussian distribution, the expectation values and variances of the semi-axis lengths of overdense ellipsoidal regions are calculated; combined with the statistical properties of the velocity shear tensor, this yields a quantitative relationship in which the Kerr parameter a* increases with latent heat strength α and decreases with phase transition rate β, reaching a typical magnitude of 10^{-3}.
What carries the argument
The nucleation history integration method within the false vacuum island model and accumulation mechanism, which tracks nonspherical collapse of overdense regions seeded by delayed bubble nucleation.
If this is right
- The Kerr parameter of these primordial black holes increases with the latent heat strength α of the phase transition.
- The Kerr parameter decreases with the phase transition rate β.
- Typical values reach 10^{-3}, exceeding those for peak-theory primordial black holes formed in the radiation-dominated era.
- Spin serves as a distinguishing feature compared to matter-era estimates.
Where Pith is reading between the lines
- If correct, spin measurements of primordial black holes could indirectly constrain the strength and duration of early-universe phase transitions.
- The mechanism might be extended to predict the full probability distribution of spins rather than only average values and variances.
- This formation channel suggests a route to separate phase-transition-seeded black holes from those produced by density fluctuations alone.
Load-bearing premise
The false vacuum island model and accumulation mechanism correctly describe the non-Gaussian density perturbations that seed primordial black hole formation from delayed bubble nucleation.
What would settle it
A measurement of primordial black hole spins in the relevant mass range that shows no dependence on phase transition parameters α and β or values far from 10^{-3}.
read the original abstract
The stochastic bubble nucleation during cosmological first-order phase transitions leads to variations in the phase transition initiation times across different Hubble volumes, thereby generating non-Gaussian density perturbations in regions with delayed transitions. Based on the accumulation mechanism and the false vacuum island model . This paper investigates the spin angular momentum of primordial black holes formed from nonspherical collapse. By introducing the nucleation history integration method, without assuming a Gaussian distribution, we calculate the expectation values and variances of the semi-axis lengths of overdense ellipsoidal regions, combined with the statistical properties of the velocity shear tensor, we derive the quantitative relationship between the Kerr parameter $a_*$ describing black hole spin and the phase transition parameters , latent heat strength $\alpha$ and phase transition rate $\beta$. The study finds that the Kerr parameter increases with $\alpha$ and decreases with $\beta$; estimate the typical the magnitude of $a_*$ can reach $10^{-3}$, which is significantly higher than that of primordial black holes formed in the radiation-dominated era under peak theory, but still lower than that in a matter-dominated era.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that primordial black holes formed via nonspherical collapse from delayed bubble nucleation in first-order phase transitions acquire spin parameterized by the Kerr parameter a*. Using the false vacuum island model together with an accumulation mechanism, the authors introduce a nucleation-history integration method to compute the expectation values and variances of the semi-axes of overdense ellipsoids (without assuming Gaussian statistics) and combine these with the velocity shear tensor to obtain an explicit dependence a*(α, β). They report that a* increases with latent-heat strength α, decreases with transition rate β, and reaches typical values ~10^{-3}, larger than radiation-era peak-theory estimates but smaller than matter-era ones.
Significance. If the false-vacuum-island statistics are reliable, the work supplies a concrete, non-Gaussian route from FOPT parameters to PBH spin distributions. This could furnish an independent observable channel (distinct from stochastic GW backgrounds) for constraining α and β, and the avoidance of Gaussian assumptions is a methodological strength.
major comments (2)
- [§2–3] The quantitative relation between a* and (α, β) is derived entirely inside the false vacuum island + accumulation model (§2–3). No comparison is presented to lattice bubble-nucleation simulations or to the peak-theory limits cited in the abstract; because the claimed non-Gaussian overdensity statistics and the resulting 10^{-3} magnitude rest on this unvalidated model, the central scaling result cannot be assessed for robustness.
- [§3–4] The nucleation-history integration yields expectation values and variances of the ellipsoid semi-axes, which are then combined with the shear tensor to produce a*. The manuscript does not show how uncertainties in the integration (or in the assumed ranges of α and β) propagate into the final a* distribution, so the statement that the “typical magnitude can reach 10^{-3}” lacks a quantified error estimate.
minor comments (2)
- [Abstract] The abstract contains a repeated definite article (“estimate the typical the magnitude”); this should be corrected for clarity.
- [§1] Notation for the phase-transition parameters α and β is introduced without an explicit reference to the standard definitions used in the FOPT literature; a brief reminder equation would aid readers.
Simulated Author's Rebuttal
We thank the referee for the constructive feedback on our manuscript. We address each major comment below and outline the revisions we will make to strengthen the presentation of our results while remaining within the analytic scope of the work.
read point-by-point responses
-
Referee: [§2–3] The quantitative relation between a* and (α, β) is derived entirely inside the false vacuum island + accumulation model (§2–3). No comparison is presented to lattice bubble-nucleation simulations or to the peak-theory limits cited in the abstract; because the claimed non-Gaussian overdensity statistics and the resulting 10^{-3} magnitude rest on this unvalidated model, the central scaling result cannot be assessed for robustness.
Authors: The false-vacuum-island plus accumulation framework is an established analytic approximation in the literature on FOPT-induced perturbations; we will add a dedicated paragraph in §2 that explicitly contrasts our computed semi-axis expectation values and variances against the Gaussian peak-theory predictions referenced in the abstract. This will include a brief derivation showing how the non-Gaussian tail from delayed nucleation increases the variance by a factor that yields a* ~ 10^{-3} rather than the smaller radiation-era peak-theory values. Direct comparison to lattice bubble-nucleation simulations is not feasible within the present analytic study, as no existing lattice data sets match the required parameter space and new simulations would constitute a separate numerical project. revision: partial
-
Referee: [§3–4] The nucleation-history integration yields expectation values and variances of the ellipsoid semi-axes, which are then combined with the shear tensor to produce a*. The manuscript does not show how uncertainties in the integration (or in the assumed ranges of α and β) propagate into the final a* distribution, so the statement that the “typical magnitude can reach 10^{-3}” lacks a quantified error estimate.
Authors: We agree that an explicit propagation of uncertainties will improve the robustness statement. In the revised manuscript we will insert a new subsection at the end of §4 that performs a sensitivity analysis: we vary the integration cutoffs and scan α ∈ [0.05, 1] and β/H ∈ [1, 100] over physically allowed ranges, recompute the ellipsoid moments, and propagate the resulting spread through the shear-tensor combination to obtain a 1σ interval on a* (approximately 3×10^{-4}–4×10^{-3}). This will replace the order-of-magnitude statement with a quantified range while preserving the central scaling with α and β. revision: yes
- Direct quantitative comparison to lattice bubble-nucleation simulations, which would require new large-scale numerical work outside the analytic scope of the present study.
Circularity Check
No circularity; forward derivation of a* from α, β via explicit integration
full rationale
The paper states it computes semi-axis lengths of overdense ellipsoids by nucleation history integration (without Gaussian assumption), then combines with velocity shear statistics to obtain a* as a function of the input parameters α and β. This is a model-based calculation whose output is not equivalent to its inputs by construction, nor does it rename a fit as a prediction. No self-citation is shown as load-bearing for the central result, and the derivation chain remains self-contained once the false-vacuum-island and accumulation assumptions are granted.
Axiom & Free-Parameter Ledger
free parameters (2)
- α (latent heat strength)
- β (phase transition rate)
axioms (2)
- domain assumption False vacuum island model correctly captures the statistics of delayed bubble nucleation regions
- domain assumption Accumulation mechanism maps nucleation time variations to overdense ellipsoidal regions
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.