Pith. sign in

REVIEW 3 major objections 4 minor 94 references

A single family of inflationary fluctuations yields three distinct routes to primordial black hole formation during kination, each with a different collapse mechanism and its own threshold condition.

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 →

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

T0 review · deepseek-v4-flash

2026-08-01 09:51 UTC pith:OMT4LY3O

load-bearing objection A careful numerical proof-of-principle showing three distinct PBH formation channels from one family of inflaton fluctuations; the spherical-symmetry caveat is real but explicitly acknowledged, and the paper deserves a serious referee. the 3 major comments →

arxiv 2607.20423 v1 pith:OMT4LY3O submitted 2026-07-22 astro-ph.CO

Primordial black holes forming during kination: the trapped, the overdense, and the void

Cristian Joana This is my paper
classification astro-ph.CO
keywords primordial black holeskinationultra-slow-roll inflationnumerical relativitycollapse thresholdcurvature profilevoid collapsenon-Gaussianity
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 shows that primordial black hole formation in a single-field model transitioning from slow roll through ultra-slow roll into kination cannot be captured by a single amplitude criterion. Using full numerical relativity, it follows localized sub-Hubble field fluctuations continuously through inflation and into collapse, finding three qualitatively different outcomes: a patch trapped in inflation, a positive-curvature overdensity, and a negative-curvature void bounded by an overdense shell. Each can form a black hole by a distinct mechanism. The paper measures collapse thresholds for the generated profiles and for analytic sinc profiles, finding that neither central curvature nor peak linear compaction provides a universal collapse criterion. The implication is that PBH abundance calculations in such models must work with complete profile distributions, treating the trapped, overdense, and void channels separately.

Core claim

On its own terms, the paper demonstrates a full end-to-end general-relativistic evolution of a localized sub-Hubble inflaton fluctuation through slow roll, ultra-slow roll, and kination to apparent-horizon formation, with no intermediate curvature profile prescribed or matched. Within this setup, three channels emerge: a sufficiently delayed patch remains trapped in inflation and becomes hidden behind a horizon in a manner akin to false-vacuum-bubble collapse; a more moderate positive displacement produces a positive-curvature overdensity that collapses in the usual way; and an advanced patch generates a negative-curvature void whose compensating overdense shell implodes, compressing the cor

What carries the argument

The central methodological object is the continuous numerical-relativity evolution of a localized field fluctuation through the SR–USR–kination transition, from which the conserved super-Hubble curvature profile is extracted; collapse thresholds are then scanned by rescaling that fixed profile's central amplitude. The paper also uses the linear compaction function C_lin(x) = -2 f(w) x ∂_x R, whose peak value is a more stable diagnostic than central curvature but still not universal. The compensating shell of the void channel is the load-bearing feature carrying collapse information, in contrast to the central curvature.

Load-bearing premise

All results assume the initial fluctuations and collapse are spherically symmetric; if non-spherical perturbations fragment the trapped patch or destabilize the void shell, the channel classification and measured thresholds may not survive.

What would settle it

A 3D numerical-relativity simulation of the same initial field profiles: if the trapped patch fragments or the void's overdense shell breaks apart before collapse, the three-channel taxonomy fails. Alternatively, a calculation that adds a short-wavelength non-spherical perturbation to the void profile and checks whether the shell still collapses would settle it.

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

If this is right

  • PBH abundance predictions in non-attractor models must be based on the distribution of complete curvature profiles, with separate treatment of trapped, overdense, and void channels.
  • The trapped channel requires a different abundance calculation, set by the probability of trapping and the horizon mass when the locally inflating patch becomes hidden, not by near-threshold scaling.
  • Collapse thresholds measured with analytic sinc profiles do not transfer to profiles generated nonlinearly by the SR–USR–kination dynamics; the shift can be tens of percent.
  • Scalar-field kination thresholds lie below those of an ideal stiff fluid, by about 10% in peak linear compaction, due to anisotropic scalar-gradient stress softening the effective response.
  • If PBHs form during kination, even a small abundance can come to dominate the universe because black holes redshift as matter while the background dilutes as a^{-6}, potentially reheating the universe via Hawking evaporation.

Where Pith is reading between the lines

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

  • If the three channels are stable in 3D, the trapped patch mechanism may resemble false-vacuum-bubble collapse and could produce PBHs in any model with a shallow local minimum, not just this one.
  • The large threshold shift between sinc and generated profiles suggests that peak-theory estimates based on the compaction function may need to be replaced by shape-aware statistics; this is an editorial extension.
  • The three channels likely have distinct gravitational-wave signatures, since each involves different accretion histories and shell dynamics; a search for such signatures could test the model.
  • A natural next calculation is to map the boundaries between the three channels in the space of initial amplitudes and shapes; the paper's three example configurations suggest such a phase diagram exists.

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 / 4 minor

Summary. This Letter reports spherically symmetric numerical-relativity simulations of primordial black hole (PBH) formation in a single-field inflation model whose background passes from slow roll through ultra-slow roll into kination. Localized, initially sub-Hubble scalar-field fluctuations are evolved continuously through the non-attractor transition, into the kination epoch, and through collapse to apparent-horizon formation. The paper identifies three qualitatively distinct outcomes within the explored initial configurations: a classically trapped inflating patch, a positive-curvature overdensity, and a negative-curvature void surrounded by an overdense shell. For one representative profile in each of the overdense and void channels, the author measures the critical central amplitude for collapse and compares it with the same quantity for analytic sinc profiles. The comparison shows that neither the central curvature nor the peak linear compaction is a profile-independent criterion, with threshold shifts of about 22% and 35% in central amplitude while peak compaction changes by only 6% and 2%. The paper concludes that PBH abundance calculations in non-attractor models should be formulated as a distribution over complete profiles, with separate treatment of the trapped, overdense, and void channels.

Significance. If the three-channel taxonomy and the associated thresholds survive further scrutiny, this is a substantial step forward for PBH formation in non-attractor models. The paper is, to my knowledge, among the first to follow an initially sub-Hubble field fluctuation through the full SR-USR-kination transition and into apparent-horizon formation without stitching together separate stages. The explicit extraction of nonlinear curvature profiles and the use of those profiles in dedicated threshold scans is a clear strength, as is the monitoring of the Hamiltonian constraint. The author is also transparent about the main limitations: the thresholds are conditional on a fixed profile shape, the trapped channel is studied only as an existence example, and the stability of the channels beyond spherical symmetry is explicitly deferred to future work. The paper does not ship machine-checked proofs, public code, or data, and it reports no convergence tests, so the quantitative claims rest on the reliability of the numerical implementation. Given those caveats, the work is a valuable but preliminary contribution that will be of interest to the numerical-relativity and PBH communities.

major comments (3)
  1. [Discussion, final paragraph; also Abstract] The concluding recommendation that PBH abundances in non-attractor models must treat the trapped, overdense, and void channels separately is a generalization beyond spherical symmetry, but all simulations are one-dimensional: the initial data are radially symmetric (Eq. 2) and the evolution code is 1D. The stability of the void-shell channel is especially nontrivial: an imploding overdense shell is susceptible to fragmentation, and an aspherical trapped-patch boundary could alter the horizon structure and the relation between patch size and PBH mass. The final paragraph explicitly defers this check to future three-dimensional simulations. This is a load-bearing caveat for the abundance prescription. I recommend either softening the abundance statements to the spherically symmetric setting or adding a stability analysis / preliminary 3D evidence. The existence of the three channels in sph
  2. [Table I and Supplemental Material, 'Code validation'] The quantitative claims—thresholds quoted to ±0.01 and the 22%/35% central-amplitude shifts—are central to the paper. The only numerical-quality indicator reported is the normalized Hamiltonian-constraint residual (Fig. 8), which is small for the representative evolutions. However, no convergence test across resolutions is shown for either the inflationary evolutions or the threshold scans. In near-critical collapse, threshold values can be sensitive to grid resolution and numerical dissipation. I request a convergence study (at least two or three resolutions) for the central evolution and for the threshold scans, and a statement of how the ±0.01 uncertainties in Table I were obtained. Without this, the quoted precision and the resulting percentage shifts are not fully supported.
  3. [Eq. (5), Table I, and 'Collapse thresholds and PBH masses'] The paper's main qualitative conclusion—profile- and channel-dependent collapse thresholds—rests on exactly one extracted profile per channel and one analytic reference profile. The extracted shape depends on the initial amplitude, comoving scale k_star, window function, and phase, as acknowledged. The statement that 'neither the central curvature nor the peak linear compaction alone provides a profile- or channel-independent collapse criterion' would be considerably strengthened by a small family of profiles (e.g., varying k_star or the oscillation phase within a channel) showing the range of R_c and C_max_lin. If such a scan is not feasible for a Letter, I suggest rephrasing the claim explicitly as valid for the specific profiles studied here and moving the general abundance prescription to the discussion as a forward-looking research direction.
minor comments (4)
  1. [Eq. (41)] The symbol R is used both for the comoving curvature perturbation and for the areal radius. This becomes confusing in a paper where R_0 also denotes the central curvature amplitude. Please use a different symbol for the areal radius, e.g., r_areal or r_m.
  2. [Eq. (30)] The gauge condition is written as \partial_t \alpha = -\mu_L \alpha p (K_phys - <K_phys>). Since 1 <= p <= 2 is a parameter, this is presumably \mu_L \alpha^p. Please fix the typographical omission of the exponent.
  3. [Abstract and Fig. 3] The abstract says the trapped configuration 'may form a PBH', but no threshold scan or horizon mass is reported for this channel. A sentence specifying that this is an existence demonstration for one configuration, with the relevant mass being the initial horizon mass rather than a re-entry mass, would avoid over-interpretation. The same clarification would help the caption of Fig. 3, since the reported masses are not asymptotic: the horizons are still accreting at the end of the simulations.
  4. [Supplemental Material, 'Numerical implementation'] No code or data release is mentioned. Given the detail of the supplemental material, publishing the grid parameters, the resolution set, and a minimal run script would significantly aid reproducibility and verification by independent groups.

Circularity Check

0 steps flagged

No circularity: thresholds are measured in independent dedicated scans, the three-channel results come from direct numerical evolutions, and self-citations are used for comparison or method, not as load-bearing derivation.

full rationale

The paper's central claims are supported by full numerical-relativity evolutions that are not fitted to the claimed results. The initial data are specified independently (Eq. 2), and the curvature profiles are extracted from continuous SR–USR–kination evolution: “with only the gauge adapted during collapse and nothing prescribed or matched between stages.” The threshold measurements are explicitly separated from the original evolutions: “The threshold measurements in Table I, however, come from dedicated scans,” and the profiles are normalized and scanned in amplitude only, with the caveat that thresholds are “conditional on a fixed shape P(x)”. The sinc reference is an external analytic benchmark, and the paper emphasizes that the extracted profiles “differ appreciably from the sinc shape,” so the comparison is not a self-fulfilling fit. Citations to the author's prior work, including the void mechanism [61] and the ideal stiff-fluid baseline [61], are used for comparison and interpretation, not as the derivation of the new measured thresholds; the void and overdense channels are directly simulated to apparent-horizon formation in this paper. The final-paragraph caveat that “three-dimensional simulations should determine the stability of the trapped, overdense, and void-shell channels beyond spherical symmetry” is an honest limitation on external validity, not a circular step: it does not make any equation in the paper reduce to its own inputs. No load-bearing self-citation chain, fitted-input-called-prediction, or definitional equivalence was found.

Axiom & Free-Parameter Ledger

3 free parameters · 5 axioms · 0 invented entities

The central claims rest on a handful of chosen model parameters, a spherical-symmetry assumption, and the validity of separate-universe/profile extraction. The ledger shows that the three-channel result is demonstrated for hand-picked configurations rather than derived from a statistical distribution.

free parameters (3)
  • Potential parameters (V0, φ, φ0, σ, λ) = V0=2.9e-10 Mpl^4, φ=0.2 Mpl, φ0=0.45 Mpl, σ=0.02 Mpl, λ=2.9685e-3
    Chosen by hand to match the CMB amplitude (As≈2.1e-9) at ~50 e-folds and to produce a USR peak Aζ≈0.08; the Gaussian feature creates the shallow local minimum needed for the trapped channel. These are inputs, not predictions.
  • Initial fluctuation amplitude φ_amp per channel = trapped −6e-4 Mpl; overdense −2.5e-4 Mpl; void +5e-4 Mpl
    Selected so that each of the three channels is realized; the sign and magnitude determine whether the patch is delayed, overdense, or void-forming. The paper states these are illustrative rather than representative.
  • Comoving scale k⋆ and window W(r) = k⋆ chosen so Hubble exit occurs near the SR–USR transition; W=1 inside 2r_max/3 with smooth falloff
    Sets the initial perturbation size and phase; the channel classification and frozen profiles depend on it. Chosen by hand to maximize nonlinear amplification.
axioms (5)
  • domain assumption Spherical symmetry is preserved and sufficient for the three channels
    All simulations are 1D radial (Eq. 9); the paper defers 3D stability to future work. If non-spherical dynamics fragment the shells or trapped patch, thresholds and channel taxonomy may change.
  • domain assumption Separate-universe expression (3)/(37) is valid on super-Hubble scales at extraction checkpoints
    Frozen curvature profiles are read from Eq. (37) assuming k⋆/(aH)≪1 and the background trajectory approximation for H/Π; the linearised version is checked to agree after kination, but the nonlinear integral is not independently validated against a gauge-invariant definition.
  • standard math The numerical relativity formulation converges to the Einstein–Klein–Gordon system
    BSSN/Z4c equations are standard; the authors monitor the normalized Hamiltonian residual |Hrel|<1e-2 but provide no convergence tests. The central claim assumes the discrete equations approximate GR.
  • domain assumption Sinc profile is a relevant analytic benchmark for USR-generated peaks
    Motivated by the logarithmic map ζ=−γ^{-1} log(1−γ log ζG) [24,37]; used as the reference template for the threshold comparison. The conclusion about profile dependence relies on this choice being representative of template-based abundance calculations.
  • domain assumption Initial conditions are classical localised fluctuations, ignoring quantum diffusion and backreaction
    The perturbations are seeded as classical field profiles (Eq. 2); stochastic effects that generate the actual distribution of fluctuations are not simulated.

pith-pipeline@v1.3.0-alltime-deepseek · 15404 in / 17013 out tokens · 139607 ms · 2026-08-01T09:51:21.052645+00:00 · methodology

0 comments
read the original abstract

We study primordial black-hole (PBH) formation in a single-field model which passes from slow roll, through a transient ultra-slow-roll phase, into kination. Nonlinear numerical-relativity simulations are used to follow localised, initially sub-Hubble field fluctuations across the non-attractor transition, to determine the resulting super-Hubble curvature profiles, and subsequently to evolve their re-entry during kination. Within the set of initial configurations explored here, we find three qualitatively distinct outcomes, namely a patch which remains trapped in inflation, a positive-curvature overdensity, and a negative-curvature void bounded by an overdense shell. Each may form a PBH, although by a different physical mechanism. For one selected profile in each of the overdense and void channels, we measure collapse thresholds and compare them with thresholds for sinc profiles. The comparison shows that neither the central curvature nor the peak linear compaction alone provides a profile- or channel-independent collapse criterion. PBH abundances in non-attractor models must therefore be formulated in terms of the distribution of complete profiles, with the trapped, overdense, and void channels treated separately.

Figures

Figures reproduced from arXiv: 2607.20423 by Cristian Joana.

Figure 1
Figure 1. Figure 1: FIG. 1. The left panel shows the scalar potential (1), with the SR (green), USR (blue), and kination (red) portions of the [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Comoving-curvature profiles at successive e-fold numbers for the trapped (left, black), overdense (centre, red), and void [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Apparent-horizon mass [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Critical curvature profiles (left) and their linear compaction functions (right). Solid curves are the nonlinear (NL) [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Inflationary potential (upper left), linear curvature spectrum (upper right), and the first two Hubble-flow parameters [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Central comoving curvature for the trapped (left), overdense (centre), and void (right) solutions, evaluated with the [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Misner–Sharp compactness during collapse for the three PBH-forming channels. Darker curves denote later times; the [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Relative Hamiltonian-constraint residual (44) for the selected overdense (left) and void (right) simulations at successive [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

94 extracted references · 76 linked inside Pith

  1. [1]

    Carr and F

    B. Carr and F. Kuhnel, Primordial Black Holes as Dark Matter: Recent Developments, Ann. Rev. Nucl. Part. Sci. 70, 355 (2020), arXiv:2006.02838 [astro-ph.CO]

    Pith/arXiv arXiv 2020

  2. [2]

    A. M. Green and B. J. Kavanagh, Primordial Black Holes as a dark matter candidate, J. Phys. G48, 043001 (2021), arXiv:2007.10722 [astro-ph.CO]

    Pith/arXiv arXiv 2021

  3. [3]

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

    Pith/arXiv arXiv 2021

  4. [4]

    Carr and J

    B. Carr and J. Silk, Primordial Black Holes as Generators of Cosmic Structures, Mon. Not. Roy. Astron. Soc.478, 3756 (2018), arXiv:1801.00672 [astro-ph.CO]

    Pith/arXiv arXiv 2018

  5. [5]

    Sasaki, T

    M. Sasaki, T. Suyama, T. Tanaka, and S. Yokoyama, Primordial black holes—perspectives in gravitational wave astronomy, Class. Quant. Grav.35, 063001 (2018), arXiv:1801.05235 [astro-ph.CO]

    Pith/arXiv arXiv 2018

  6. [6]

    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), arXiv:2306.16219 [astro-ph.HE]

    Pith/arXiv arXiv 2023

  7. [7]

    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]

    Pith/arXiv arXiv 2023

  8. [8]

    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]

    Pith/arXiv arXiv 2023

  9. [9]

    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]

    Pith/arXiv arXiv 2023

  10. [10]

    Baguiet al.(LISA Cosmology Working Group), Pri- mordial black holes and their gravitational-wave signa- tures, Living Rev

    E. Baguiet al.(LISA Cosmology Working Group), Pri- mordial black holes and their gravitational-wave signa- tures, Living Rev. Rel.28, 1 (2025), arXiv:2310.19857 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  11. [11]

    Shankaranarayanan, S

    S. Shankaranarayanan, S. Bhattacharya, and A. Vid- yarthi, Primordial Black Holes: A Review of Formation and Evolution (2026), arXiv:2606.23846 [gr-qc]

    Pith/arXiv arXiv 2026

  12. [12]

    Hawking, Gravitationally collapsed objects of very low mass, Mon

    S. Hawking, Gravitationally collapsed objects of very low mass, Mon. Not. Roy. Astron. Soc.152, 75 (1971)

    1971

  13. [13]

    B. J. Carr and S. W. Hawking, Black holes in the early Universe, Mon. Not. Roy. Astron. Soc.168, 399 (1974)

    1974

  14. [14]

    B. J. Carr, The Primordial black hole mass spectrum, Astrophys. J.201, 1 (1975)

    1975

  15. [15]

    M. Y. Khlopov, B. A. Malomed, I. B. Zeldovich, and Y. B. Zeldovich, Gravitational instability of scalar fields and formation of primordial black holes, Mon. Not. Roy. Astron. Soc.215, 575 (1985)

    1985

  16. [16]

    Motohashi, A

    H. Motohashi, A. A. Starobinsky, and J. Yokoyama, Inflation with a constant rate of roll, JCAP09, 018, arXiv:1411.5021 [astro-ph.CO]

    Pith/arXiv arXiv

  17. [17]

    Motohashi and A

    H. Motohashi and A. A. Starobinsky, Constant-roll infla- tion: confrontation with recent observational data, EPL 117, 39001 (2017), arXiv:1702.05847 [astro-ph.CO]

    Pith/arXiv arXiv 2017

  18. [18]

    Motohashi, S

    H. Motohashi, S. Mukohyama, and M. Oliosi, Con- stant Roll and Primordial Black Holes, JCAP03, 002, arXiv:1910.13235 [gr-qc]

    Pith/arXiv arXiv 1910

  19. [19]

    Pi and M

    S. Pi and M. Sasaki, Gravitational Waves Induced by Scalar Perturbations with a Lognormal Peak, JCAP09, 037, arXiv:2005.12306 [gr-qc]

    Pith/arXiv arXiv 2005

  20. [20]

    Tomberg, Stochastic constant-roll inflation and pri- mordial black holes, Phys

    E. Tomberg, Stochastic constant-roll inflation and pri- mordial black holes, Phys. Rev. D108, 043502 (2023), arXiv:2304.10903 [astro-ph.CO]

    Pith/arXiv arXiv 2023

  21. [21]

    Pattison, V

    C. Pattison, V. Vennin, H. Assadullahi, and D. Wands, Quantum diffusion during inflation and primordial black holes, JCAP10, 046, arXiv:1707.00537 [hep-th]

    Pith/arXiv arXiv

  22. [22]

    J. M. Ezquiaga, J. Garc ´ ıa-Bellido, and V. Vennin, The exponential tail of inflationary fluctuations: con- sequences for primordial black holes, JCAP03, 029, arXiv:1912.05399 [astro-ph.CO]

    Pith/arXiv arXiv 1912

  23. [23]

    D. G. Figueroa, S. Raatikainen, S. Rasanen, and E. Tomberg, Non-Gaussian Tail of the Curvature Per- turbation in Stochastic Ultraslow-Roll Inflation: Im- plications for Primordial Black Hole Production, Phys. Rev. Lett.127, 101302 (2021), arXiv:2012.06551 [astro- ph.CO]

    arXiv 2021

  24. [24]

    Pi and M

    S. Pi and M. Sasaki, Logarithmic Duality of the Curva- ture Perturbation, Phys. Rev. Lett.131, 011002 (2023), arXiv:2211.13932 [astro-ph.CO]

    Pith/arXiv arXiv 2023

  25. [25]

    Caravano, G

    A. Caravano, G. Franciolini, and S. Renaux-Petel, Ultraslow-roll inflation on the lattice: Backreaction and nonlinear effects, Phys. Rev. D111, 063518 (2025), arXiv:2410.23942 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  26. [26]

    Caravano, G

    A. Caravano, G. Franciolini, and S. Renaux-Petel, Ultraslow-roll inflation on the lattice. II. Nonperturba- tive curvature perturbation, Phys. Rev. D112, 083508 (2025), arXiv:2506.11795 [astro-ph.CO]

    arXiv 2025

  27. [27]

    Musco, J

    I. Musco, J. C. Miller, and L. Rezzolla, Computations of primordial black hole formation, Class. Quant. Grav.22, 1405 (2005), arXiv:gr-qc/0412063

    Pith/arXiv arXiv 2005

  28. [28]

    Musco and J

    I. Musco and J. C. Miller, Primordial black hole for- mation in the early universe: critical behaviour and self-similarity, Class. Quant. Grav.30, 145009 (2013), arXiv:1201.2379 [gr-qc]

    Pith/arXiv arXiv 2013

  29. [29]

    Harada, C.-M

    T. Harada, C.-M. Yoo, T. Nakama, and Y. Koga, Cosmological long-wavelength solutions and primordial black hole formation, Phys. Rev. D91, 084057 (2015), arXiv:1503.03934 [gr-qc]

    Pith/arXiv arXiv 2015

  30. [30]

    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]

    Pith/arXiv arXiv 2020

  31. [31]

    Escriv` a, C

    A. Escriv` a, C. Germani, and R. K. Sheth, Analytical thresholds for black hole formation in general cosmolog- ical backgrounds, JCAP01, 030, arXiv:2007.05564 [gr- qc]

    Pith/arXiv arXiv 2007

  32. [32]

    Escriv` a, F

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

    Pith/arXiv arXiv 2022

  33. [33]

    Escriv` a, The threshold for PBH formation in the type- II region and its analytical estimation, arXiv e-prints (2025), arXiv:2504.05814 [astro-ph.CO]

    A. Escriv` a, The threshold for PBH formation in the type- II region and its analytical estimation, arXiv e-prints (2025), arXiv:2504.05814 [astro-ph.CO]

    arXiv 2025

  34. [34]

    Yuwen, C

    Z.-Y. Yuwen, C. Joana, S.-J. Wang, and R.-G. Cai, Pri- mordial black hole formation in bulk-viscous cosmology, arXiv e-prints (2026), arXiv:2606.26532 [gr-qc]

    Pith/arXiv arXiv 2026

  35. [35]

    V. Atal, J. Garriga, and A. Marcos-Caballero, Primor- dial black hole formation with non-Gaussian curvature 7 perturbations, JCAP09, 073, arXiv:1905.13202 [astro- ph.CO]

    Pith/arXiv arXiv 1905

  36. [36]

    V. Atal, J. Cid, A. Escriv` a, and J. Garriga, PBH in single field inflation: the effect of shape dispersion and non- Gaussianities, JCAP05, 022, arXiv:1908.11357 [astro- ph.CO]

    Pith/arXiv arXiv 1908

  37. [37]

    R. Inui, C. Joana, H. Motohashi, S. Pi, Y. Tada, and S. Yokoyama, Primordial black holes and induced gravi- tational waves from logarithmic non-Gaussianity, JCAP 03, 021, arXiv:2411.07647 [astro-ph.CO]

    Pith/arXiv arXiv

  38. [38]

    Young, Peaks and primordial black holes: the effect of non-Gaussianity, JCAP05(05), 037, arXiv:2201.13345 [astro-ph.CO]

    S. Young, Peaks and primordial black holes: the effect of non-Gaussianity, JCAP05(05), 037, arXiv:2201.13345 [astro-ph.CO]

    Pith/arXiv arXiv

  39. [39]

    Clough, E

    K. Clough, E. A. Lim, B. S. DiNunno, W. Fischler, R. Flauger, and S. Paban, Robustness of Inflation to Inhomogeneous Initial Conditions, JCAP09, 025, arXiv:1608.04408 [hep-th]

    Pith/arXiv arXiv

  40. [40]

    Clough, R

    K. Clough, R. Flauger, and E. A. Lim, Robustness of Inflation to Large Tensor Perturbations, JCAP05, 065, arXiv:1712.07352 [hep-th]

    Pith/arXiv arXiv

  41. [41]

    Joana and S

    C. Joana and S. Clesse, Inhomogeneous preinflation across Hubble scales in full general relativity, Phys. Rev. D103, 083501 (2021), arXiv:2011.12190 [astro-ph.CO]

    Pith/arXiv arXiv 2021

  42. [42]

    Joana, Gravitational dynamics in Higgs inflation: Pre- inflation and preheating with an auxiliary field, Phys

    C. Joana, Gravitational dynamics in Higgs inflation: Pre- inflation and preheating with an auxiliary field, Phys. Rev. D106, 023504 (2022), arXiv:2202.07604 [astro- ph.CO]

    Pith/arXiv arXiv 2022

  43. [43]

    Joana, Beginning inflation in conformally curved spacetimes, Phys

    C. Joana, Beginning inflation in conformally curved spacetimes, Phys. Rev. D110, 063534 (2024), arXiv:2406.00811 [astro-ph.CO]

    Pith/arXiv arXiv 2024

  44. [44]

    J. C. Aurrekoetxea, K. Clough, R. Flauger, and E. A. Lim, The Effects of Potential Shape on Inhomoge- neous Inflation, JCAP05, 030, arXiv:1910.12547 [astro- ph.CO]

    Pith/arXiv arXiv 1910

  45. [45]

    Elley, J

    M. Elley, J. C. Aurrekoetxea, K. Clough, R. Flauger, P. Giannadakis, and E. A. Lim, Robustness of in- flation to kinetic inhomogeneities, JCAP01, 050, arXiv:2405.03490 [astro-ph.CO]

    Pith/arXiv arXiv

  46. [46]

    S. E. Brady, K. Clough, P. Figueras, and ´A. D. Kov´ acs, Inflaton Dynamics in Higher-Derivative Scalar- Tensor Theories of Gravity, arXiv e-prints (2025), arXiv:2505.17986 [gr-qc]

    Pith/arXiv arXiv 2025

  47. [47]

    S. E. Brady, T. W. Baumgarte, and K. Clough, Starting inflation in asymptotically flat spacetimes, arXiv e-prints (2026), arXiv:2607.06441 [gr-qc]

    Pith/arXiv arXiv 2026

  48. [48]

    J. C. Aurrekoetxea, K. Clough, and E. A. Lim, Cosmol- ogy using numerical relativity, Living Rev. Rel.28, 5 (2025), arXiv:2409.01939 [gr-qc]

    Pith/arXiv arXiv 2025

  49. [49]

    de Jong, J

    E. de Jong, J. C. Aurrekoetxea, and E. A. Lim, Primor- dial black hole formation with full numerical relativity, JCAP03(03), 029, arXiv:2109.04896 [astro-ph.CO]

    Pith/arXiv arXiv

  50. [50]

    X.-X. Kou, C. Tian, and S.-Y. Zhou, Oscillon Preheat- ing in Full General Relativity, Class. Quant. Grav.38, 045005 (2021), arXiv:1912.09658 [gr-qc]

    Pith/arXiv arXiv 2021

  51. [51]

    Cheng, P

    C. Cheng, P. Giannadakis, L. Heurtier, and E. A. Lim, Non-linear Dynamics and Primordial Black Hole Formation During Kination, arXiv e-prints (2025), arXiv:2507.19166 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  52. [52]

    L. E. Padilla, E. Milligan, D. J. Mulryne, and J. C. Hi- dalgo, Primordial Black Hole Formation in a Scalar Field Dominated Universe: Investigation of the Critical nature of the Collapse, arXiv e-prints (2025), arXiv:2509.10431 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  53. [53]

    Milligan, L

    E. Milligan, L. E. Padilla, D. J. Mulryne, and J. C. Hi- dalgo, Primordial black hole formation in a scalar field dominated universe, JCAP10, 025, arXiv:2504.02600 [astro-ph.CO]

    Pith/arXiv arXiv

  54. [54]

    S. K. Blau, E. I. Guendelman, and A. H. Guth, The Dy- namics of False Vacuum Bubbles, Phys. Rev. D35, 1747 (1987)

    1987

  55. [55]

    Garriga, A

    J. Garriga, A. Vilenkin, and J. Zhang, Black holes and the multiverse, JCAP02, 064, arXiv:1512.01819 [hep-th]

    Pith/arXiv arXiv

  56. [56]

    Deng and A

    H. Deng and A. Vilenkin, Primordial black hole formation by vacuum bubbles, JCAP12, 044, arXiv:1710.02865 [gr-qc]

    Pith/arXiv arXiv

  57. [57]

    Deng, Primordial black hole formation by vacuum bubbles

    H. Deng, Primordial black hole formation by vacuum bubbles. Part II, JCAP09, 023, arXiv:2006.11907 [astro- ph.CO]

    Pith/arXiv arXiv 2006

  58. [58]

    Escriv` a, V

    A. Escriv` a, V. Atal, and J. Garriga, Formation of trapped vacuum bubbles during inflation, and consequences for PBH scenarios, JCAP10, 035, arXiv:2306.09990 [astro- ph.CO]

    arXiv

  59. [59]

    Wang, Y.-l

    H. Wang, Y.-l. Zhang, and T. Suyama, Nearly Monochro- matic Primordial Black Holes as total Dark Mat- ter from Bubble Collapse, arXiv e-prints (2025), arXiv:2510.19233 [astro-ph.CO]

    arXiv 2025

  60. [60]

    Franciolini, M

    G. Franciolini, M. Peloso, and A. Riotto, Dark Matter from Eternity, arXiv e-prints (2026), arXiv:2602.08338 [astro-ph.CO]

    arXiv 2026

  61. [61]

    Joana and Z.-Y

    C. Joana and Z.-Y. Yuwen, Primordial black holes from primordial voids, Phys. Rev. D113, 023518 (2026), arXiv:2510.11611 [astro-ph.CO]

    arXiv 2026

  62. [62]

    Martin, C

    J. Martin, C. Ringeval, and V. Vennin, Encyclopædia Inflationaris: Opiparous Edition, Phys. Dark Univ.5-6, 75 (2014), arXiv:1303.3787 [astro-ph.CO]

    Pith/arXiv arXiv 2014

  63. [63]

    Kallosh and A

    R. Kallosh and A. Linde, Universality Class in Conformal Inflation, JCAP07, 002, arXiv:1306.5220 [hep-th]

    Pith/arXiv arXiv

  64. [64]

    T. W. Baumgarte and S. L. Shapiro, Numerical integra- tion of einstein’s field equations, Physical Review D59, 10.1103/physrevd.59.024007 (1998)

    doi:10.1103/physrevd.59.024007 1998

  65. [65]

    Shibata and T

    M. Shibata and T. Nakamura, Evolution of three- dimensional gravitational waves: Harmonic slicing case, Phys. Rev. D52, 5428 (1995)

    1995

  66. [66]

    D. Alic, C. Bona-Casas, C. Bona, L. Rezzolla, and C. Palenzuela, Conformal and covariant formulation of the Z4 system with constraint-violation damping, Phys. Rev. D85, 064040 (2012), arXiv:1106.2254 [gr-qc]

    Pith/arXiv arXiv 2012

  67. [67]

    Bernuzzi and D

    S. Bernuzzi and D. Hilditch, Constraint violation in free evolution schemes: comparing BSSNOK with a con- formal decomposition of Z4, Phys. Rev. D81, 084003 (2010), arXiv:0912.2920 [gr-qc]

    Pith/arXiv arXiv 2010

  68. [68]

    A. G. Polnarev and I. Musco, Curvature profiles as ini- tial conditions for primordial black hole formation, Class. Quant. Grav.24, 1405 (2007), arXiv:gr-qc/0605122

    Pith/arXiv arXiv 2007

  69. [69]

    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]

    Pith/arXiv arXiv 2019

  70. [70]

    Germani and L

    C. Germani and L. Montell` a, The trichotomy of primor- dial black holes initial conditions, arXiv e-prints (2025), arXiv:2510.02006 [gr-qc]

    arXiv 2025

  71. [71]

    Germani and R

    C. Germani and R. K. Sheth, The Statistics of Pri- mordial Black Holes in a Radiation-Dominated Uni- verse: Recent and New Results, Universe9, 421 (2023), arXiv:2308.02971 [astro-ph.CO]

    Pith/arXiv arXiv 2023

  72. [72]

    Ning, X.-X

    Z. Ning, X.-X. Zeng, R.-G. Cai, and S.-J. Wang, Nu- merical simulations of primordial black hole formation 8 via delayed first-order phase transitions, arXiv e-prints (2026), arXiv:2601.21878 [gr-qc]

    Pith/arXiv arXiv 2026

  73. [73]

    R. K. Sheth and R. van de Weygaert, A hierarchy of voids: Much ado about nothing, Mon. Not. Roy. Astron. Soc.350, 517 (2004), arXiv:astro-ph/0311260

    Pith/arXiv arXiv 2004

  74. [74]

    Young, Computation of the Abundance of Primordial Black Holes (2025) arXiv:2405.13259 [astro-ph.CO]

    S. Young, Computation of the Abundance of Primordial Black Holes (2025) arXiv:2405.13259 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  75. [75]

    Escriv` a, The statistics of curvature-profile dispersion in primordial black hole formation, arXiv e-prints (2026), arXiv:2607.08738 [astro-ph.CO]

    A. Escriv` a, The statistics of curvature-profile dispersion in primordial black hole formation, arXiv e-prints (2026), arXiv:2607.08738 [astro-ph.CO]

    Pith/arXiv arXiv 2026

  76. [76]

    Dom` enech and J

    G. Dom` enech and J. Tr¨ ankle, From formation to evapora- tion: Induced gravitational wave probes of the primordial black hole reheating scenario, Phys. Rev. D111, 063528 (2025), arXiv:2409.12125 [astro-ph.CO]

    Pith/arXiv arXiv 2025

  77. [77]

    Martin, T

    J. Martin, T. Papanikolaou, and V. Vennin, Primordial black holes from the preheating instability in single-field inflation, JCAP01, 024, arXiv:1907.04236 [astro-ph.CO]

    Pith/arXiv arXiv 1907

  78. [78]

    Auclair and V

    P. Auclair and V. Vennin, Primordial black holes from metric preheating: mass fraction in the excursion-set ap- proach, JCAP02, 038, arXiv:2011.05633 [astro-ph.CO]

    Pith/arXiv arXiv 2011

  79. [79]

    Dom` enech, C

    G. Dom` enech, C. Lin, and M. Sasaki, Gravitational wave constraints on the primordial black hole dominated early universe, JCAP04, 062, [Erratum: JCAP 11, E01 (2021)], arXiv:2012.08151 [gr-qc]

    Pith/arXiv arXiv 2021

  80. [80]

    Balaji, G

    S. Balaji, G. Dom` enech, G. Franciolini, A. Ganz, and J. Tr¨ ankle, Probing modified Hawking evaporation with gravitational waves from the primordial black hole domi- nated universe, JCAP11, 026, arXiv:2403.14309 [gr-qc]

    Pith/arXiv arXiv

Showing first 80 references.