Pith. sign in

REVIEW 3 major objections 5 minor 9 cited by

A gauge-invariant calculation shows that curvature perturbations from slow first-order phase transitions are far smaller than earlier separate-universe estimates, making primordial black hole formation from this mechanism very unlikely.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-03 09:12 UTC pith:P5OK6JWT

load-bearing objection Solid confirmation that slow-FOPT super-horizon inhomogeneities don't form PBHs once gauge issues are handled; the new fitting template and PTA analysis are useful, but the linearization assumption needs a check. the 3 major comments →

arxiv 2601.14412 v2 pith:P5OK6JWT submitted 2026-01-20 hep-ph

How large are curvature perturbations from slow first-order phase transitions? A gauge-invariant analysis

classification hep-ph
keywords first-order phase transitionscurvature perturbationsprimordial black holesgauge invarianceseparate universesupercooled transitionsscalar-induced gravitational wavesmulti-fluid perturbation theory
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper argues that the super-horizon density contrasts produced when a slow, strongly supercooled cosmological first-order phase transition nucleates bubbles are much smaller than previously thought. Working with a two-fluid system of radiation and vacuum, the authors compute the gauge-invariant comoving curvature perturbation and find its amplitude is suppressed relative to earlier separate-universe results. Consequently, the density contrast distribution lies below the typical black-hole formation threshold, so primordial black holes are unlikely to form through this mechanism. The paper also provides a fitting formula for the curvature power spectrum and shows that scalar-induced gravitational waves from these perturbations are subdominant.

Core claim

The central claim is that the comoving curvature perturbation generated by super-horizon inhomogeneities from slow first-order phase transitions is much smaller than implied by previous separate-universe studies, and the resulting density contrast does not exceed the conventional primordial black hole formation threshold. Using gauge-invariant multi-fluid perturbation theory, the authors evolve the curvature perturbation on super-Hubble scales, accounting for the relative entropy perturbation between radiation and vacuum. They derive a fitting formula for the curvature power spectrum that scales as k^3 at small k and beta^-5 at large beta, but with an amplitude suppressed compared to earlier

What carries the argument

The central object is the gauge-invariant comoving curvature perturbation R (equivalently zeta), evolved on super-Hubble scales via d(zeta)/dt = -H/(rho+P) deltaP_nad, where deltaP_nad is the total non-adiabatic pressure perturbation. For a two-fluid system of radiation and vacuum, the intrinsic non-adiabatic pressure vanishes because each fluid has a definite barotropic equation of state, leaving only the relative entropy perturbation between the fluids. The authors avoid gauge ambiguities by treating the density contrasts from separate-universe-style nucleation simulations as linear-order perturbations on a specific but unspecified time slicing, then use the gauge-invariant equations to co

Load-bearing premise

The calculation assumes that the density contrasts from separate-universe simulations can be approximated, at linear order, as perturbations on a specific but unspecified time slicing, so that the gauge-invariant multi-fluid formalism applies.

What would settle it

A fully nonlinear simulation of bubble nucleation in an expanding universe that computes the comoving curvature perturbation without invoking the linear-slicing assumption would settle the issue; if it produced a density contrast distribution with a tail above about 0.45 for some parameters, the central suppression claim would be refuted.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Primordial black hole formation from super-horizon inhomogeneities at slow first-order phase transitions is strongly disfavored; the density contrast distribution stays below the critical threshold.
  • Scalar-induced gravitational waves from these curvature perturbations are subdominant compared to primary bubble-collision or sound-wave signals.
  • The derived fitting formula for the curvature power spectrum can be used to translate current and future constraints on primordial perturbations into constraints on phase transition parameters.
  • Future radio and dark-matter-halo observations may constrain both slow and fast phase transitions, with the suppressed spectrum shifting the reachable parameter space.
  • The gauge-invariant treatment provides a template for evaluating other cosmological observables arising from large-scale phase-transition inhomogeneities.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the suppression holds, the window for primordial black hole dark matter from this specific mechanism closes, though other formation channels such as false vacuum islands remain unaffected.
  • The linear-slicing assumption could be tested by comparing against a fully nonlinear simulation of bubble nucleation; if nonlinear corrections are large, the suppression may change.
  • The fitting formula offers a ready way to forecast LISA-era constraints on supercooled transitions, which the paper only sketches via existing datasets.
  • The result underscores that any claim of black hole production from phase transitions must specify the gauge in which the density contrast is defined.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper claims that super-horizon curvature perturbations generated by slow, strongly supercooled first-order phase transitions are much smaller than earlier separate-universe estimates suggested, so primordial black hole production from these perturbations is unlikely. The authors use the multi-fluid gauge-invariant perturbation formalism of Refs. [32,33], feeding in density contrasts δ_NG obtained from the DeltaPT separate-universe simulation of Ref. [24]. They compute the comoving curvature perturbation R via the relative entropy perturbation S_Vr, present its variance for several β/H_n and k/k_max, derive a fitting formula for P_R(k) (Eq. 20), and use this template to discuss constraints from CMB, Lyman-α, FIRAS, PTA, SKA, LISA, and UCMHs, as well as scalar-induced gravitational waves and a NANOGrav 15-year fit. The central conclusion is that the δ_C distribution lies far below the adopted PBH threshold δ_crit ≈ 0.45.

Significance. If the main claim is correct, the paper resolves a gauge ambiguity that has affected recent estimates of PBHs and SIGWs from slow FOPTs, and it provides a useful fitting formula for phenomenological studies. The use of the standard multi-fluid gauge-invariant formalism is appropriate, and the paper is explicit about its key assumption. The numerical analysis with 50,000 realizations per parameter point is a strength, as is the inclusion of current and projected observational constraints and a PTArcade-based posterior analysis. However, the central result is load-bearing on a linearization assumption that is not quantitatively validated, and the simulation data/code are not shipped, so the core suppression claim is not yet fully demonstrable to the standards of a journal publication.

major comments (3)
  1. [Sec. IV.A; Eq. (18); Fig. 1] The central result depends on the assumption, stated in Sec. IV.A, that the separate-universe density contrasts δ_NG can be approximated at linear order as density contrasts on a specific but unknown time-slicing. While a linear time shift cancels in S_Vr, the linearization itself is not justified: the manuscript nowhere shows δ_NG ≪ 1. Figure 1 plots |δ_NG| and |9/4 δ_NG| on an axis extending to order unity, and no convergence test in δ is reported. Since the final |R| values are orders of magnitude smaller (~10^-4 or less), second-order corrections to the gauge transformation and to the mapping from a finite patch to a global perturbation enter S_Vr and δP_rel at O(δ^2) and could change the sign or magnitude of R. The conclusion that PBHs are unlikely is therefore not robust against this most obvious failure mode. A quantitative bound on δ_NG in the relevant (β/H_n, k/k_max) regime, or
  2. [Sec. IV.B; Fig. 2; Eq. (20)] No statistical uncertainties are given for σ_R^2. The claim that deviations from k^3 scaling at large β are 'likely attributable to the limited number of bubbles' is not supported by any convergence test or error estimate. Since Eq. (20) is fitted to these same σ_R^2 values, the fitting coefficients A=0.0038, a=0.043, b=1.77 inherit these unquantified uncertainties. The observational constraints in Sec. V and Table I are computed using this template, so the presented exclusion regions are not accompanied by a systematic error budget. Please provide error bars (e.g., jackknife or bootstrap over histories) and a convergence test in the number of bubbles/realizations, or state clearly that the template is a central-value fit.
  3. [Sec. V; Eq. (20); Table I] The fitting parameters in Eq. (20) are fitted to the same DeltaPT simulations that produce the central quantities, and the same template is then used to derive all observational constraints and the SIGW spectra in Sec. V. This is not circular for the direct PBH conclusion, which is computed from the simulated R distribution, but it limits the independent validation of the template. The manuscript would be stronger if the simulation input (or at least the σ_R^2 data points) were provided as supplementary material, and if the template fit were validated against out-of-sample parameter points. Without this, the Table I constraints should be seen as dependent on an unvalidated functional form.
minor comments (5)
  1. [Sec. IV.B; Fig. 2] The horizontal axis label in the right panel appears garbled: '1016 × 100 2 × 101' should read '10^0, 2×10^1' (or similar). Please fix this and ensure all axis labels are clear.
  2. [Sec. IV.A] The text states 'For each super-horizon wave number k > kmax', but the figures and earlier discussion concern modes k < k_max that re-enter after percolation. Please correct this apparent typo.
  3. [Eq. (18) and surrounding text] The cancellation of the divergence in S_Vr when ρ̇_r = 0 is described only in words. Since the prefactor in Eq. (17) also contains 1/(6H ρ̇), it would be helpful to show the explicit cancellation or state the limiting behavior when both ρ̇_r and ρ̇_V approach zero.
  4. [Eq. (20)] Typo: 'heaviside' should be 'Heaviside'. Also, the phrase 'Θ(k_max − k)' should be written with a space or consistent notation.
  5. [Table I] The ranges for β/H_n (e.g., '>~ 35−80') are not explained. Are these ranges due to different constraints inside a column, or to the adopted reference values? Please clarify in the text or caption.

Circularity Check

0 steps flagged

No significant circularity: the central PBH-suppression claim is computed directly from the simulated density contrasts via standard gauge-invariant multi-fluid perturbation theory, not from the fitted template.

full rationale

The derivation chain is self-contained. The paper takes DeltaPT separate-universe density contrasts δ_NG (Ref. [24]) and, using the standard gauge-invariant multi-fluid formalism (Eqs. 9–18, citing Refs. [31–33]), computes the relative entropy perturbation S_Vr and integrates Eq. (13) to obtain ζ and R ≈ −ζ. The key assumption—that δ_NG can be treated, at linear order, as density contrasts on some unknown slicing—is an explicit physical premise (Sec. IV.A), not a reduction of the output to the input; the linear-order gauge invariance of S_Vr is what makes the computation well-defined. The PBH-unlikely conclusion is read directly off the simulated δ_C distribution in Fig. 3, not off the fitted template. Equation (20)'s coefficients A, a, b are fitted to the same simulations, but the paper labels it a fitting formula and uses it only for observational projections (Sec. V); it is not used to derive the central suppression claim. No load-bearing self-citation or imported uniqueness theorem appears; the authors' own previous work (Refs. [9–12]) is cited only for background on SGWB and phase-transition dynamics. The main risk is the unvalidated linearization premise, which is a correctness/physics concern, not circularity.

Axiom & Free-Parameter Ledger

3 free parameters · 7 axioms · 0 invented entities

The central calculation relies on standard multi-fluid perturbation theory plus a stated but unvalidated mapping from separate-universe simulation outputs to a time slicing. The only free parameters introduced for the paper's own purposes are the three fitting coefficients in Eq. (20).

free parameters (3)
  • A = 0.0038
    Amplitude of the fitting formula for P_R in Eq. (20), fitted to the DeltaPT simulation variances in Sec. IV.B.
  • a = 0.043
    Shape parameter in Eq. (20), fitted to the simulated σ_R^2(k, β/H_n) data.
  • b = 1.77
    Shape parameter in Eq. (20), fitted to the simulated σ_R^2(k, β/H_n) data.
axioms (7)
  • ad hoc to paper Separate-universe density contrasts δ_NG are approximated, at linear order, as density contrasts on some specific but unknown time-slicing.
    Sec. IV.A explicitly states this assumption. It is necessary to map simulation outputs into the gauge-invariant multi-fluid formalism, and it is the main unresolved premise.
  • domain assumption The bubble nucleation rate follows Γ = H_n^4 exp(β t).
    Eq. (1), standard exponential nucleation model used throughout FOPT phenomenology.
  • domain assumption The false-vacuum fraction F(t) obeys the standard expression involving the integral of the nucleation rate.
    Eq. (3), taken from Refs. [28,29]; standard for radiation-filled backgrounds with vacuum domination.
  • domain assumption Vacuum energy behaves as a barotropic fluid with P_V = -ρ_V and c_V^2 = -1.
    Used to set the vacuum sound speed and to cancel the divergence in S_Vr via background quantities; standard in multi-fluid treatments of vacuum.
  • standard math Linear perturbation theory and the super-horizon relation R ≈ -ζ are valid for the computed density contrasts.
    Invoked in Sec. III and IV; the paper notes the assumption is 'at linear order', but does not quantify the size of neglected nonlinear terms.
  • domain assumption The PBH formation threshold is δ_crit ≈ 0.45.
    Adopted in Sec. IV.B from Refs. [44-46]; the conclusion that PBHs are unlikely depends on this threshold value.
  • domain assumption UCMH-based constraints assume WIMP dark matter annihilating into b bbar with m_DM = 1 TeV and s-wave thermal relic cross section.
    Stated in Sec. V.A and Fig. 4 caption; the authors note other assumptions would relax these limits.

pith-pipeline@v1.3.0-alltime-deepseek · 11979 in / 12623 out tokens · 132624 ms · 2026-08-03T09:12:07.163330+00:00 · methodology

0 comments
read the original abstract

When strongly supercooled cosmological first-order phase transitions (FOPTs) are sufficiently slow, super-horizon inhomogeneities can be generated. We compute these super-horizon curvature perturbations by employing a gauge-invariant, multi-fluid formalism. By resolving the gauge ambiguities inherent in conventional separate-universe simulations, we demonstrate that Primordial Black Holes are unlikely to be produced by these super-horizon inhomogeneities. We also derive a fitting formula for the resulting curvature perturbations and discuss potential observational constraints on FOPTs imposed by limits on primordial curvature perturbations and associated scalar-induced gravitational waves.

Figures

Figures reproduced from arXiv: 2601.14412 by Chi Tian, Csaba Bal\'azs, Ran Ding, Xiao Wang.

Figure 1
Figure 1. Figure 1: FIG. 1. Time evolution of the absolute value of the false vac [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Variance of the comoving curvature perturbations [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Distributions of density contrast of independently [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Various benchmark spectra [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. GWs spectra from [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Posterior distributions of phase transition parameters from NANOGrav-15 data, based on [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 9 Pith papers

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

  1. Can the universe be matter-dominated after a supercooled first-order phase transition?

    hep-ph 2026-07 conditional novelty 7.0

    After a supercooled first-order phase transition, the scalar field's equation of state is set by the bubble-wall Lorentz factor γ*, and matter domination is delayed until a/a* ≃ γ* in the free-streaming limit.

  2. Are Primordial Black Holes a Natural Dark Matter Candidate?

    hep-ph 2026-06 unverdicted novelty 7.0

    PBH dark matter spans all naturalness tiers, with some mechanisms as natural as WIMPs or freeze-in particles, determined by abundance map structure rather than candidate type.

  3. Gravitational Waves from Black Hole Reheating: The Scalar-Induced Component

    hep-ph 2026-05 unverdicted novelty 7.0

    Accounting for the minimal mass spread of primordial black holes from gravitational collapse suppresses the Poltergeist GW background to the level of generic scalar-induced signals and reopens ultra-light PBH parameter space.

  4. Reviving primordial black hole formation in slow first-order phase transitions

    hep-ph 2026-05 unverdicted novelty 7.0

    Primordial black hole formation from slow first-order phase transitions remains viable if slow reheating after supercooling creates an early matter-dominated era allowing small overdensities to grow, producing black h...

  5. Primordial Black Hole from Tensor-induced Density Fluctuation: First-order Phase Transitions and Domain Walls

    astro-ph.CO 2026-05 unverdicted novelty 6.0

    Tensor perturbations from first-order phase transitions and domain wall annihilation induce curvature fluctuations at second order that form primordial black holes, allowing asteroid-mass PBHs to comprise all dark mat...

  6. Primordial Black Holes from Slow Phase Transitions with Delayed Reheating: A Peak-Theory Approach

    hep-ph 2026-06 unverdicted novelty 5.0

    PBH production from slow phase transitions with delayed reheating is modeled via peak theory and Monte Carlo simulations, showing extreme sensitivity to reheating efficiency and potential to explain all dark matter.

  7. Primordial Black Hole from Tensor-induced Density Fluctuation: First-order Phase Transitions and Domain Walls

    astro-ph.CO 2026-05 reject novelty 5.0

    Tensor perturbations from FOPT and domain-wall sources are claimed to induce second-order scalar perturbations large enough to form primordial black holes, potentially all of the dark matter.

  8. Reviving primordial black hole formation in slow first-order phase transitions

    hep-ph 2026-05 unverdicted novelty 5.0

    Slow reheating after a supercooled first-order phase transition allows an early matter-dominated era in which small curvature perturbations grow sufficiently to form primordial black holes.

  9. Phenomenology of Vector Dark Matter produced by a First Order Phase Transition

    hep-ph 2026-04 unverdicted novelty 5.0

    Dark sector first-order phase transitions near 10 MeV can substantially modify vector dark matter relic densities away from standard thermal freeze-out predictions, with distinct mass windows and calculable gravitatio...

Reference graph

Works this paper leans on

56 extracted references · 1 canonical work pages · cited by 7 Pith papers

  1. [1]

    V. A. Kuzmin, V. A. Rubakov, and M. E. Shaposhnikov, On the Anomalous Electroweak Baryon Number Non- conservation in the Early Universe, Phys. Lett. B155, 36 (1985)

  2. [2]

    A. G. Cohen, D. B. Kaplan, and A. E. Nelson, Progress in electroweak baryogenesis, Ann. Rev. Nucl. Part. Sci. 43, 27 (1993), arXiv:hep-ph/9302210

  3. [3]

    V. A. Rubakov and M. E. Shaposhnikov, Electroweak baryon number nonconservation in the early universe and in high-energy collisions, Usp. Fiz. Nauk166, 493 (1996), arXiv:hep-ph/9603208

  4. [4]

    Riotto and M

    A. Riotto and M. Trodden, Recent progress in baryogen- esis, Ann. Rev. Nucl. Part. Sci.49, 35 (1999), arXiv:hep- ph/9901362

  5. [5]

    D. E. Morrissey and M. J. Ramsey-Musolf, Elec- troweak baryogenesis, New J. Phys.14, 125003 (2012), arXiv:1206.2942 [hep-ph]

  6. [6]

    Hindmarsh, S

    M. Hindmarsh, S. J. Huber, K. Rummukainen, and D. J. Weir, Gravitational waves from the sound of a first order phase transition, Phys. Rev. Lett.112, 041301 (2014), arXiv:1304.2433 [hep-ph]

  7. [7]

    Capriniet al., Detecting gravitational waves from cos- mological phase transitions with LISA: an update, JCAP 03, 024, arXiv:1910.13125 [astro-ph.CO]

    C. Capriniet al., Detecting gravitational waves from cos- mological phase transitions with LISA: an update, JCAP 03, 024, arXiv:1910.13125 [astro-ph.CO]

  8. [8]

    Athron, C

    P. Athron, C. Bal´ azs, A. Fowlie, L. Morris, and L. Wu, Cosmological phase transitions: From perturbative parti- cle physics to gravitational waves, Prog. Part. Nucl. Phys. 135, 104094 (2024), arXiv:2305.02357 [hep-ph]

  9. [9]

    C. Tian, X. Wang, and C. Bal´ azs, Gravitational waves from cosmological first-order phase transitions with pre- cise hydrodynamics, Eur. Phys. J. C85, 1091 (2025), arXiv:2409.14505 [hep-ph]

  10. [10]

    X. Wang, C. Tian, and C. Bal´ azs, Self-consistent pre- diction of gravitational waves from cosmological phase transitions, (2024), arXiv:2409.06599 [hep-ph]

  11. [11]

    C. Tian, X. Wang, and C. Bal´ azs, DeepSSM: an emulator of gravitational wave spectra from sound waves during cosmological first-order phase transitions, JCAP08, 060, arXiv:2501.10244 [astro-ph.CO]

  12. [12]

    Balazs, X

    C. Balazs, X. Wang, and C. Tian, Gravitational waves from cosmological first-order phase transitions with self-consistent hydrodynamics, PoSCORFU2024, 261 (2025)

  13. [13]

    Agazieet al.(NANOGrav), The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background, Astrophys

    G. Agazieet al.(NANOGrav), The NANOGrav 15 yr Data Set: Evidence for a Gravitational-wave Background, Astrophys. J. Lett.951, L8 (2023), arXiv:2306.16213 [astro-ph.HE]

  14. [14]

    Xuet al., Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I, Res

    H. Xuet al., Searching for the Nano-Hertz Stochastic Gravitational Wave Background with the Chinese Pulsar Timing Array Data Release I, Res. Astron. Astrophys. 23, 075024 (2023), arXiv:2306.16216 [astro-ph.HE]

  15. [15]

    Antoniadiset al.(EPTA, InPTA:), The second data release from the European Pulsar Timing Array - III

    J. Antoniadiset al.(EPTA, InPTA:), The second data release from the European Pulsar Timing Array - III. Search for gravitational wave signals, Astron. Astrophys. 678, A50 (2023), arXiv:2306.16214 [astro-ph.HE]

  16. [16]

    D. J. Reardonet al., Search for an Isotropic Gravitational-wave Background with the Parkes Pul- sar Timing Array, Astrophys. J. Lett.951, L6 (2023), arXiv:2306.16215 [astro-ph.HE]

  17. [17]

    M. T. Mileset al., The MeerKAT Pulsar Timing Array: the first search for gravitational waves with the MeerKAT radio telescope, Mon. Not. Roy. Astron. Soc.536, 1489 9 (2024), arXiv:2412.01153 [astro-ph.HE]

  18. [18]

    Afzalet al.(NANOGrav), The NANOGrav 15 yr Data Set: Search for Signals from New Physics, Astrophys

    A. Afzalet al.(NANOGrav), The NANOGrav 15 yr Data Set: Search for Signals from New Physics, Astrophys. J. Lett.951, L11 (2023), [Erratum: Astrophys.J.Lett. 971, L27 (2024), Erratum: Astrophys.J. 971, L27 (2024)], arXiv:2306.16219 [astro-ph.HE]

  19. [19]

    J. Liu, L. Bian, R.-G. Cai, Z.-K. Guo, and S.-J. Wang, Constraining First-Order Phase Transitions with Curva- ture Perturbations, Phys. Rev. Lett.130, 051001 (2023), arXiv:2208.14086 [astro-ph.CO]

  20. [20]

    Franciolini, Y

    G. Franciolini, Y. Gouttenoire, and R. Jinno, Curva- ture Perturbations from First-Order Phase Transitions: Implications to Black Holes and Gravitational Waves, (2025), arXiv:2503.01962 [hep-ph]

  21. [21]

    Cai, Y.-S

    R.-G. Cai, Y.-S. Hao, and S.-J. Wang, Primordial black holes and curvature perturbations from false vacuum is- lands, Sci. China Phys. Mech. Astron.67, 290411 (2024), arXiv:2404.06506 [astro-ph.CO]

  22. [22]

    Fu and J

    C. Fu and J. Liu, The power of SKA to Constrain cos- mological gravitational-wave backgrounds below the as- trophysical foreground noise, (2025), arXiv:2506.14366 [astro-ph.CO]

  23. [23]

    Gouttenoire, WIMPs and new physics interpreta- tions of the PTA signal are incompatible, (2025), arXiv:2503.03857 [hep-ph]

    Y. Gouttenoire, WIMPs and new physics interpreta- tions of the PTA signal are incompatible, (2025), arXiv:2503.03857 [hep-ph]

  24. [24]

    Lewicki, P

    M. Lewicki, P. Toczek, and V. Vaskonen, Black Holes and Gravitational Waves from Slow First-Order Phase Transitions, Phys. Rev. Lett.133, 221003 (2024), arXiv:2402.04158 [astro-ph.CO]

  25. [25]

    J. Liu, L. Bian, R.-G. Cai, Z.-K. Guo, and S.-J. Wang, Primordial black hole production during first- order phase transitions, Phys. Rev. D105, L021303 (2022), arXiv:2106.05637 [astro-ph.CO]

  26. [26]

    Kawana, T

    K. Kawana, T. Kim, and P. Lu, PBH formation from overdensities in delayed vacuum transitions, Phys. Rev. D108, 103531 (2023), arXiv:2212.14037 [astro-ph.CO]

  27. [27]

    Gouttenoire and T

    Y. Gouttenoire and T. Volansky, Primordial black holes from supercooled phase transitions, Phys. Rev. D110, 043514 (2024), arXiv:2305.04942 [hep-ph]

  28. [28]

    M. S. Turner, E. J. Weinberg, and L. M. Widrow, Bubble nucleation in first order inflation and other cosmological phase transitions, Phys. Rev. D46, 2384 (1992)

  29. [29]

    X. Wang, F. P. Huang, and X. Zhang, Phase transition dynamics and gravitational wave spectra of strong first- order phase transition in supercooled universe, JCAP05, 045, arXiv:2003.08892 [hep-ph]

  30. [30]

    G. Elor, R. Jinno, S. Kumar, R. McGehee, and Y. Tsai, Finite Bubble Statistics Constrain Late Cosmological Phase Transitions, Phys. Rev. Lett.133, 211003 (2024), arXiv:2311.16222 [hep-ph]

  31. [31]

    Wands, K

    D. Wands, K. A. Malik, D. H. Lyth, and A. R. Liddle, A New approach to the evolution of cosmological pertur- bations on large scales, Phys. Rev. D62, 043527 (2000), arXiv:astro-ph/0003278

  32. [32]

    K. A. Malik, D. Wands, and C. Ungarelli, Large scale cur- vature and entropy perturbations for multiple interact- ing fluids, Phys. Rev. D67, 063516 (2003), arXiv:astro- ph/0211602

  33. [33]

    K. A. Malik and D. Wands, Adiabatic and entropy per- turbations with interacting fluids and fields, JCAP02, 007, arXiv:astro-ph/0411703

  34. [34]

    J. M. Bardeen, Gauge Invariant Cosmological Perturba- tions, Phys. Rev. D22, 1882 (1980)

  35. [35]

    V. F. Mukhanov, H. A. Feldman, and R. H. Branden- berger, Theory of cosmological perturbations. Part 1. Classical perturbations. Part 2. Quantum theory of per- turbations. Part 3. Extensions, Phys. Rept.215, 203 (1992)

  36. [36]

    A. M. Green, A. R. Liddle, K. A. Malik, and M. Sasaki, A New calculation of the mass fraction of primordial black holes, Phys. Rev. D70, 041502 (2004), arXiv:astro- ph/0403181

  37. [37]

    Aghanimet al.(Planck), Planck 2018 results

    N. Aghanimet al.(Planck), Planck 2018 results. VI. Cosmological parameters, Astron. Astrophys.641, A6 (2020), [Erratum: Astron.Astrophys. 652, C4 (2021)], arXiv:1807.06209 [astro-ph.CO]

  38. [38]

    S. Bird, H. V. Peiris, M. Viel, and L. Verde, Mini- mally Parametric Power Spectrum Reconstruction from the Lyman-alpha Forest, Mon. Not. Roy. Astron. Soc. 413, 1717 (2011), arXiv:1010.1519 [astro-ph.CO]

  39. [39]

    Chluba, A

    J. Chluba, A. L. Erickcek, and I. Ben-Dayan, Probing the inflaton: Small-scale power spectrum constraints from measurements of the CMB energy spectrum, Astrophys. J.758, 76 (2012), arXiv:1203.2681 [astro-ph.CO]

  40. [40]

    C. T. Byrnes, P. S. Cole, and S. P. Patil, Steepest growth of the power spectrum and primordial black holes, JCAP 06, 028, arXiv:1811.11158 [astro-ph.CO]

  41. [41]

    Inomata and T

    K. Inomata and T. Nakama, Gravitational waves in- duced by scalar perturbations as probes of the small-scale primordial spectrum, Phys. Rev. D99, 043511 (2019), arXiv:1812.00674 [astro-ph.CO]

  42. [42]

    M. S. Delos, A. L. Erickcek, A. P. Bailey, and M. A. Al- varez, Density profiles of ultracompact minihalos: Impli- cations for constraining the primordial power spectrum, Phys. Rev. D98, 063527 (2018), arXiv:1806.07389 [astro- ph.CO]

  43. [43]

    Franco Abell´ an and G

    G. Franco Abell´ an and G. Facchinetti, Minihalos as probes of the inflationary spectrum: accurate boost fac- tor calculation and new CMB constraints, JCAP06, 032, arXiv:2304.02996 [astro-ph.CO]

  44. [44]

    Escriv` a, PBH Formation from Spherically Symmetric Hydrodynamical Perturbations: A Review, Universe8, 66 (2022), arXiv:2111.12693 [gr-qc]

    A. Escriv` a, PBH Formation from Spherically Symmetric Hydrodynamical Perturbations: A Review, Universe8, 66 (2022), arXiv:2111.12693 [gr-qc]

  45. [45]

    Musco, Threshold for primordial black holes: Depen- dence on the shape of the cosmological perturbations, Phys

    I. Musco, Threshold for primordial black holes: Depen- dence on the shape of the cosmological perturbations, Phys. Rev. D100, 123524 (2019), arXiv:1809.02127 [gr- qc]

  46. [46]

    Escriv` a, C

    A. Escriv` a, C. Germani, and R. K. Sheth, Universal threshold for primordial black hole formation, Phys. Rev. D101, 044022 (2020), arXiv:1907.13311 [gr-qc]

  47. [47]

    Saikawa and S

    K. Saikawa and S. Shirai, Precise WIMP Dark Mat- ter Abundance and Standard Model Thermodynamics, JCAP08, 011, arXiv:2005.03544 [hep-ph]

  48. [48]

    A. S. Josan and A. M. Green, Gamma-rays from ultra- compact minihalos: potential constraints on the primor- dial curvature perturbation, Phys. Rev. D82, 083527 (2010), arXiv:1006.4970 [astro-ph.CO]

  49. [49]

    Bringmann, P

    T. Bringmann, P. Scott, and Y. Akrami, Improved con- straints on the primordial power spectrum at small scales from ultracompact minihalos, Phys. Rev. D85, 125027 (2012), arXiv:1110.2484 [astro-ph.CO]

  50. [50]

    Nakama, T

    T. Nakama, T. Suyama, K. Kohri, and N. Hiroshima, Constraints on small-scale primordial power by annihi- lation signals from extragalactic dark matter minihalos, Phys. Rev. D97, 023539 (2018), arXiv:1712.08820 [astro- ph.CO]

  51. [51]

    S. J. Huber and T. Konstandin, Gravitational Wave Pro- duction by Collisions: More Bubbles, JCAP09, 022, 10 arXiv:0806.1828 [hep-ph]

  52. [52]

    Caprini, R

    C. Caprini, R. Durrer, and G. Servant, The stochastic gravitational wave background from turbulence and mag- netic fields generated by a first-order phase transition, JCAP12, 024, arXiv:0909.0622 [astro-ph.CO]

  53. [53]

    Kohri and T

    K. Kohri and T. Terada, Semianalytic calculation of gravitational wave spectrum nonlinearly induced from primordial curvature perturbations, Phys. Rev. D97, 123532 (2018), arXiv:1804.08577 [gr-qc]

  54. [54]

    J. R. Espinosa, D. Racco, and A. Riotto, A Cosmologi- cal Signature of the SM Higgs Instability: Gravitational Waves, JCAP09, 012, arXiv:1804.07732 [hep-ph]

  55. [55]

    Mitridate, D

    A. Mitridate, D. Wright, R. von Eckardstein, T. Schr¨ oder, J. Nay, K. Olum, K. Schmitz, and T. Trickle, PTArcade, (2023), arXiv:2306.16377 [hep-ph]

  56. [56]

    Mitridate, Ptarcade 10.5281/zenodo.7876430 (2023)

    A. Mitridate, Ptarcade 10.5281/zenodo.7876430 (2023)