Pith. sign in

REVIEW 3 major objections 6 minor 4 cited by

Primordial Black Hole Formation in a Scalar Field Dominated Universe: Investigation of the Critical nature of the Collapse

T0 review · 3 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash

Pith's one-line read A quartic scalar field and a radiation fluid collapse into primordial black holes with nearly identical critical exponents.

desk verdict Useful first measurement of the quartic-scalar critical exponent, but the abstract overstates the difference from radiation and the gauge-defining fluid is never tested. read the letter →

arxiv 2509.10431 v1 pith:Y3J7Z3DH submitted 2025-09-12 astro-ph.CO

classification astro-ph.CO
keywords primordialblackholescriticalgravitationalcollapseexponentquarticscalarfieldtypeIIMisner-Sharpformalismradiationfluidnearuniversality
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper asks whether a universe dominated by a self-interacting scalar field with a quartic potential forms primordial black holes in the same critical way as a radiation-filled universe. Using fully nonlinear, spherically symmetric simulations in the Misner-Sharp formalism, it finds that the black hole mass near threshold follows the same type II power law in both cases, with critical exponents agreeing to within about $2\sigma$. If correct, this extends the common shortcut of modelling oscillating scalar fields as perfect fluids into the strongly nonlinear collapse regime and supports the near universality of the critical exponent in primordial black hole formation. The paper reports a small residual difference between the exponents rather than claiming exact equality.

What carries the argument

The central object is the power-law scaling of the black hole mass near threshold, $M_{\rm BH}\propto(p-p_{\rm th})^\gamma$, in type II critical collapse. The machinery that extracts it is the Misner-Sharp formulation of spherical relativistic collapse, evolved with a fourth-order Runge-Kutta method and finite differences, with a scalar field plus an extremely diluted homogeneous perfect fluid that defines the comoving frame. Black hole formation is detected through the apparent-horizon condition $\Theta_+=0$; the final mass is read from the Misner-Sharp mass once the apparent-horizon velocity satisfies $v_{\rm AH}\approx 1$, and $\gamma$ is obtained by fitting the mass against $p-p_{\rm th}$.

What would settle it

Run the same scalar-field collapse with the homogeneous background fluid density reduced by factors of 10 while holding the scalar perturbation fixed; if the fitted exponent $\gamma$ moves, or if the fluid's share of the Misner-Sharp mass near horizon formation grows, the quoted scalar-field exponent is contaminated by the frame-defining fluid.

Watch

Extended reading notes

Core claim

The paper's central claim is that a quartic self-interacting scalar field, $V(\psi)\propto\psi^4$, exhibits type II critical collapse (black hole mass becoming arbitrarily small near threshold and scaling as a power law) with critical exponent $\gamma = 0.3401 \pm 0.0071$, statistically close to the radiation-fluid value $\gamma = 0.3474 \pm 0.004$ measured with the same code. The threshold amplitudes are also close: $p_{\rm th}^{\rm sf}=0.0271205$ versus $p_{\rm th}^{\rm rad}=0.0270895$. The authors present this as the first explicit verification that the scalar-field/perfect-fluid correspondence extends into the critical collapse regime, at least for the quartic potential, while emphasising that the two exponents differ by about $2\sigma$ and that the quadratic case departs from dust-like behaviour.

Load-bearing premise

The scalar-field simulations carry a tiny homogeneous perfect fluid purely to define the comoving frame, and the paper does not show that this fluid stays dynamically negligible at the critical threshold.

Editorial extensions

If this is right

  • For quartic-potential scenarios, primordial black hole abundance calculations that approximate the field as a radiation fluid should capture the near-threshold scaling, with the $\sim2\%$ exponent difference entering as a small systematic uncertainty.
  • The same code reproduces the known radiation critical exponent, so the comparison baseline is not an artefact of the numerical setup.
  • The quadratic scalar-field case departs from dust-like behaviour, so the fluid analogy is potential-dependent and cannot be assumed for all scalar-field models.
  • The $2\sigma$ separation between the two exponents defines a concrete target: higher-resolution runs will either close the gap or establish a genuine matter-model dependence of $\gamma$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Varying the potential exponent between the quadratic and quartic cases, or changing the quartic coupling, could map how $\gamma$ drifts with the effective equation of state and reveal whether the $2\sigma$ gap is a physical trend rather than a numerical artefact.
  • The frame-defining homogeneous fluid could be eliminated by reformulating the evolution in a different gauge; if the same exponent survives, the near-universality claim would no longer rest on that technical crutch.
  • The compaction function oscillates near threshold in the scalar-field runs, so peak-theory estimates of primordial black hole abundance that rely on a single maximum of the compaction function may need oscillation-aware prescriptions.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 6 minor

Summary. This paper uses a Misner–Sharp, spherically symmetric, fully nonlinear relativistic code to study critical collapse of a quartic self-interacting scalar field and a radiation fluid in an expanding universe. It reports type II critical behavior in both cases, with critical exponents γ = 0.3474 ± 0.004 (radiation) and γ = 0.3401 ± 0.0071 (scalar field), and interprets the difference as about 2σ, supporting the near universality of the critical exponent in primordial black hole formation. The numerical methodology includes horizon detection, excision, Hamiltonian-constraint monitoring, and resolution tests.

Significance. If the result is robust, this is a useful first direct numerical test of the scalar-field/perfect-fluid correspondence in the critical collapse regime, with implications for PBH mass calculations in scalar-field-dominated early-universe scenarios. The paper's strengths are that it uses a fully nonlinear evolution with constraint monitoring, reports resolution dependence, and compares two matter models with the same code. However, the central quantitative claim depends on an unquantified artificial fluid in the scalar-field runs, and the stated 2σ separation is not supported by the quoted errors.

major comments (3)
  1. [II.C.2 and Eqs. (12e), (12h)] The diluted perfect fluid is not demonstrated to be dynamically negligible. The fluid is not a passive marker: it contributes to the stress-energy tensor and, through Eq. (12h), it directly fixes the lapse gradient. Eq. (12e) sources δρ_pf/ρ_pf through terms independent of the fluid's initial amplitude because the source is proportional to (ρ_pf+P_pf), so fractional perturbations of order unity can develop even for an 'extremely diluted' initial fluid. The manuscript gives no value for ρ_b0pf, does not include fluid perturbations in the initial-data list, and shows no diagnostic for the fluid density in Fig. 4. Please quantify ρ_b0pf and demonstrate by explicit convergence tests, for example by varying ρ_b0pf over several orders of magnitude, that γ_sf is unaffected; otherwise the measured exponent may be contaminated by the gauge fluid.
  2. [Abstract and Sections V.B, VI] The claimed 'about 2σ' difference between the exponents is not consistent with the quoted errors. The difference is 0.3474 − 0.3401 = 0.0073; the quadrature-summed error is sqrt(0.004^2 + 0.0071^2) ≈ 0.0082, giving about 0.9σ. Please correct this statement and avoid interpreting the result as evidence for a statistically significant difference.
  3. [IV.B and V.A] The quoted uncertainties in γ are fit errors only. The conclusion of 'near universality' rests on a small difference between exponents, so systematic uncertainties from the choice of fitting range, resolution, threshold bracket in Eq. (36), horizon-mass extrapolation, and the excision criterion need to be estimated. Without such estimates, the comparison between γ_sf and γ_rad cannot be assigned a reliable significance.
minor comments (6)
  1. [Throughout] There are typos: 'dominateed' in the Introduction, 'collaps' in the Introduction, and 'threashold' in Section V.A.
  2. [II.C.2, Eq. (16)] The initial condition for the scalar field is normalized as ~Π_b0^2/2 = 1 − ρ_b0pf, but ρ_b0pf is never specified; please state its value and justify that it is too small to alter the background dynamics.
  3. [V.A, Eq. (36)] The threshold is defined as the midpoint between p_min_bh and p_max_no-bh, but the bracketing procedure and how close to threshold the simulations are run are not described; please specify the bisection or bracketing algorithm.
  4. [Figs. 2 and 5] The set of resolutions used for the convergence tests is not stated in the captions; please list the dA values and the corresponding line styles.
  5. [III.C and IV.B] The notation for the horizon velocity criterion is inconsistent: Section IV.B uses v_c − 1 < 3 × 10^-3, while Section III.C and Appendix A define v_AH; please unify the notation.
  6. [IV.A, footnote 2] The footnote acknowledges that w = 1/3 is only an effective equation of state for the quartic scalar field; near criticality the field may not oscillate rapidly, so please discuss how this could affect the interpretation of the comparison with radiation.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the critical exponents are dynamical outputs of independent numerical evolutions, and self-citations are methodological rather than load-bearing.

full rationale

The central claim—gamma_sf = 0.3401 +/- 0.0071 versus gamma_rad = 0.3474 +/- 0.004—is obtained by evolving the full Misner-Sharp equations (12) for different initial amplitudes p, locating apparent horizons, measuring PBH masses, and fitting the scaling law M proportional to (p - p_th)^gamma. The critical exponent is not an input to the code; it is a fitting output. The initial conditions (16)-(19) are physical setup choices (purely kinetic scalar field, Gaussian curvature profile, background equation of state w = 1/3, quartic coupling lambda4 = 10) and do not encode the target scaling exponent. The paper's reliance on Ref. [21] for the derivation of initial conditions and the numerical code is a self-citation, but it is methodological rather than load-bearing: the scalar-field result is compared with a radiation-fluid run performed with the same code and with the external value of Musco et al. The 'inherently circular' remark about determining lambda by iteration is a numerical fixed-point procedure, not a logical circularity. The homogeneous, extremely diluted fluid used to define the Misner-Sharp comoving frame is an admitted modeling choice; whether it contaminates the dynamics is a robustness/correctness concern, not a case of the prediction being equivalent to its inputs by construction. No fitted parameter is renamed as a prediction, and no uniqueness theorem from the authors' prior work is invoked to forbid alternatives. Therefore no circular step is identified.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The critical exponent measurement does not rely on fitting the target result, but the physical setup contains several hand-fixed parameters (coupling, profile width, diluted fluid density, initial phase) that could in principle affect the exponent. The paper does not vary these to demonstrate universality.

free parameters (4)
  • p (initial amplitude of Gaussian overdensity) = p_rad_th = 0.0270895, p_sf_th = 0.0271205
    Control parameter tuned to bracket the collapse threshold; the critical exponent is fitted as a function of p - p_th.
  • Sigma (Gaussian width) = 5 R_H
    Fixed for all runs; the critical exponent is assumed independent of initial profile shape, but this is not tested.
  • lambda_hat_4 (quartic coupling) = 10
    Sets the strength of the scalar field self-interaction relative to the Hubble scale; only one value is simulated.
  • rho_b0pf (diluted fluid density) = not specified
    A diluted perfect fluid is added in scalar field runs to define the comoving frame; its value is not reported and its effect is not tested.
assumptions (5)
  • domain assumption The Misner-Sharp formalism with spherical symmetry is sufficient to capture critical collapse in an expanding universe.
    The evolution equations are solved in spherical symmetry using a comoving fluid frame; no non-spherical or back-reaction effects are considered.
  • domain assumption Initial conditions are valid to first order in gradient and amplitude perturbations, selecting the growing mode.
    Initial data are taken from gradient expansion results (Refs. [21,34,35]) without testing higher-order corrections for the scalar field case.
  • domain assumption The scalar field starts purely kinetic at the initial time (psi0 = 0, chi0 = 0).
    Fixes the initial phase of the scalar field oscillations; no exploration of phase dependence is performed.
  • domain assumption The background universe obeys an effective barotropic equation of state w = 1/3 for the quartic scalar field.
    Uses the averaged equation of state for a rapidly oscillating quartic scalar field, while the full dynamics include oscillations.
  • ad hoc to paper The diluted perfect fluid in scalar field runs is dynamically negligible.
    The fluid is required for the comoving Misner-Sharp frame but its dynamical effect near criticality is not validated by tests.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Primordial Black Hole Formation in a Scalar Field Dominated Universe: Investigation of the Critical nature of the Collapse." pith.science (2026). https://pith.science/paper/Y3J7Z3DH

@misc{pith2026250910431,
  author       = {Pith},
  title        = {Pith review of: Primordial Black Hole Formation in a Scalar Field Dominated Universe: Investigation of the Critical nature of the Collapse},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Y3J7Z3DH}},
  note         = {Machine review of arXiv:2509.10431}
}
abstract

In this paper, we investigate the critical collapse leading to primordial black hole (PBH) formation in a universe dominated by a self-interacting scalar field with a quartic potential, comparing it to the well-known radiation-dominated case. Using fully relativistic nonlinear numerical simulations in spherical symmetry, based on the Misner--Sharp formalism, we analyze the dynamics near the collapse threshold and track the scaling of the black hole mass. Our results confirm that both the scalar field and radiation cases exhibit type II critical behavior with similar -- though not identical -- critical exponents, differing by about $2\sigma$. This suggests that, while a quartic scalar field effectively mimics a radiation fluid even in the nonlinear collapse regime, small differences in the critical exponent persist. Our findings provide direct numerical evidence for the near universality of the critical exponent in PBH formation, with only mild dependence on whether the collapse is driven by a scalar field or a perfect fluid.

Figures

Figures reproduced from arXiv: 2509.10431 by the authors.

Figure 1
Figure 1. FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 3
Figure 3. FIG. 3 [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. shows the evolution of the Hamiltonian con￾straint and the maximum of the compaction function for different initial conditions. As before, the compaction function clearly distinguishes between collapsing and dis￾persing configurations, while the Hamiltonian constraint remains well controlled, showing significant growth only near criticality. We note, however, that in the scalar field case the maximum of the compacti… view at source ↗
Figures from the paper (2 more)
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6 [PITH_FULL_IMAGE:figures/full_fig_p011_6.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 4 Pith papers

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

  1. Long-term 3+1 simulations of primordial black hole formation during radiation domination

    gr-qc 2026-08 conditional novelty 7.0 of 10

    A new time-step gauge makes long 3D primordial black hole formation simulations about 94 times cheaper and reproduces the known spherical collapse threshold and critical exponent.

  2. The statistics of curvature-profile dispersion in primordial black hole formation

    astro-ph.CO 2026-07 conditional novelty 7.0 of 10

    Shape dispersion around the average peak profile is a genuine statistical ingredient: rare deformed curvature profiles can dominate primordial black hole formation when the power spectrum is broad or non-Gaussianity i...

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

    astro-ph.CO 2026-07 conditional novelty 6.0 of 10

    In one non-attractor inflation model, initially similar field fluctuations produce three distinct black-hole formation channels—trapped, overdense, void—with collapse thresholds determined by the full density profile,...

  4. Numerical simulations of primordial black hole formation via delayed first-order phase transitions

    gr-qc 2026-01 conditional novelty 6.0 of 10

    Spherically symmetric numerical relativity shows false-vacuum domains from delayed first-order phase transitions form type B (baby-universe) or type A (direct-collapse) primordial black holes, separated by a robust t_...

Reference graph

Works this paper leans on

44 extracted references · 40 canonical work pages · cited by 4 Pith papers

  1. [1]

    Initial conditions for a perfect fluid-dominated scenario The idea behind the derivation of our initial conditions is to linearize the equations of motion by expressing the 5 variables as ˜X= ˜Xb +ϵδ X,(13) where ˜Xstands for the relevant dynamical quantities, such as ˜mpf, ˜U, ˜R, ˜ρpf, etc.; ˜Xb denotes the correspond- ing background value1, andϵis a fo...

  2. [2]

    Laboratorio de Mod- elos y Datos (LAMOD) para proyectos de Investigaci´ on Cient´ ıfica: Censos Astrof´ ısicos

    Initial conditions for a scalar field-dominated scenario Since our system is evolved in coordinates comoving with a perfect fluid, even in this scalar field-dominated 1 Note that most variables are unity at the background level. scenario, we include a homogeneous, isotropic, and ex- tremely diluted fluid component to regularize the origin and maintain the...

  3. [3]

    Universality and scaling in gravita- tional collapse of a massless scalar field,

    M. W. Choptuik, “Universality and scaling in gravita- tional collapse of a massless scalar field,” Phys. Rev. Lett., vol. 70, pp. 9–12, Jan 1993

  4. [4]

    Dynamics of pri- mordial black hole formation,

    J. C. Niemeyer and K. Jedamzik, “Dynamics of pri- mordial black hole formation,” Phys. Rev. D, vol. 59, p. 124013, May 1999

  5. [5]

    Near-critical gravita- tional collapse and the initial mass function of primordial black holes,

    J. C. Niemeyer and K. Jedamzik, “Near-critical gravita- tional collapse and the initial mass function of primordial black holes,” Phys. Rev. Lett., vol. 80, pp. 5481–5484, Jun 1998

  6. [6]

    The dynamics of primordial black-hole formation,

    I. Hawke and J. M. Stewart, “The dynamics of primordial black-hole formation,” Classical and Quantum Gravity, vol. 19, p. 3687, jun 2002

  7. [7]

    Computa- tions of primordial black-hole formation,

    I. Musco, J. C. Miller, and L. Rezzolla, “Computa- tions of primordial black-hole formation,” Classical and Quantum Gravity, vol. 22, p. 1405, mar 2005

  8. [8]

    Curvature profiles as initial conditions for primordial black hole formation,

    A. G. Polnarev and I. Musco, “Curvature profiles as initial conditions for primordial black hole formation,” Classical and Quantum Gravity, vol. 24, p. 1405, mar 2007

Show all 44 references
  1. [9]

    Primordial black hole formation in the radiative era: Investigation of the critical nature of the collapse,

    I. Musco, J. C. Miller, and A. G. Polnarev, “Primordial black hole formation in the radiative era: Investigation of the critical nature of the collapse,” Class. Quant. Grav., vol. 26, p. 235001, 2009

  2. [10]

    Critical be- havior in gravitational collapse of a yang-mills field,

    M. W. Choptuik, T. Chmaj, and P. Bizo´ n, “Critical be- havior in gravitational collapse of a yang-mills field,” 13 Phys. Rev. Lett., vol. 77, pp. 424–427, Jul 1996

  3. [11]

    Phases of massive scalar field collapse,

    P. R. Brady, C. M. Chambers, and S. M. C. V. Gon¸ calves, “Phases of massive scalar field collapse,” Phys. Rev. D, vol. 56, pp. R6057–R6061, Nov 1997

  4. [12]

    Boson stars driven to the brink of black hole formation,

    S. H. Hawley and M. W. Choptuik, “Boson stars driven to the brink of black hole formation,” Phys. Rev. D, vol. 62, p. 104024, Oct 2000

  5. [13]

    Critical gravi- tational collapse of a massive complex scalar field,

    E. Jimenez-Vazquez and M. Alcubierre, “Critical gravi- tational collapse of a massive complex scalar field,”Phys. Rev. D, vol. 106, no. 4, p. 044071, 2022

  6. [14]

    Critical phenom- ena in gravitational collapse,

    C. Gundlach and J. M. Martin-Garcia, “Critical phenom- ena in gravitational collapse,” Living Rev. Rel., vol. 10, p. 5, 2007

  7. [15]

    Collapse threshold for a cosmological Klein Gordon field,

    J. C. Hidalgo, J. De Santiago, G. German, N. Barbosa- Cendejas, and W. Ruiz-Luna, “Collapse threshold for a cosmological Klein Gordon field,” Phys. Rev. D, vol. 96, no. 6, p. 063504, 2017

  8. [16]

    Gravother- malizing into primordial black holes, boson stars, and cannibal stars,

    P. Ralegankar, D. Perri, and T. Kobayashi, “Gravother- malizing into primordial black holes, boson stars, and cannibal stars,” 10 2024

  9. [17]

    Complex Scalar Field Reheating and Primordial Black Hole production,

    K. Carrion, J. C. Hidalgo, A. Montiel, and L. E. Padilla, “Complex Scalar Field Reheating and Primordial Black Hole production,” JCAP, vol. 07, p. 001, 2021

  10. [18]

    Primordial black holes from the preheating instability in single-field inflation,

    J. Martin, T. Papanikolaou, and V. Vennin, “Primordial black holes from the preheating instability in single-field inflation,” JCAP, vol. 01, p. 024, 2020

  11. [19]

    Black hole threshold solutions in stiff fluid col- lapse,

    P. R. Brady, M. W. Choptuik, C. Gundlach, and D. W. Neilsen, “Black hole threshold solutions in stiff fluid col- lapse,” Class. Quant. Grav., vol. 19, pp. 6359–6376, 2002

  12. [20]

    Growth of primordial black holes in a universe containing a massless scalar field,

    T. Harada and B. J. Carr, “Growth of primordial black holes in a universe containing a massless scalar field,” Phys. Rev. D, vol. 71, p. 104010, 2005

  13. [21]

    Long- wavelength nonlinear perturbations of a complex scalar field,

    L. E. Padilla, J. C. Hidalgo, and D. N´ u˜ nez, “Long- wavelength nonlinear perturbations of a complex scalar field,” Phys. Rev. D, vol. 104, no. 8, p. 083513, 2021

  14. [22]

    Coherent scalar-field oscillations in an expanding universe,

    M. S. Turner, “Coherent scalar-field oscillations in an expanding universe,” Phys. Rev. D, vol. 28, pp. 1243– 1247, Sep 1983

  15. [23]

    Primordial Black Hole Formation in a Scalar Field Dominated Universe,

    E. Milligan, L. E. Padilla, D. J. Mulryne, and J. C. Hidalgo, “Primordial Black Hole Formation in a Scalar Field Dominated Universe,” 4 2025

  16. [24]

    New mechanism for primordial black hole formation during reheating,

    L. E. Padilla, J. C. Hidalgo, and K. A. Malik, “New mechanism for primordial black hole formation during reheating,” Phys. Rev. D, vol. 106, no. 2, p. 023519, 2022

  17. [25]

    Produc- tion of PBHs from inflaton structures,

    J. C. Hidalgo, L. E. Padilla, and G. German, “Produc- tion of PBHs from inflaton structures,” Phys. Rev. D, vol. 107, no. 6, p. 063519, 2023

  18. [26]

    Constrain- ing inflationary potentials with inflaton PBHs,

    L. E. Padilla, J. C. Hidalgo, and G. German, “Constrain- ing inflationary potentials with inflaton PBHs,” Phys. Rev. D, vol. 108, no. 6, p. 063511, 2023

  19. [27]

    Detecting the stochastic gravitational wave background from primordial black holes in slow-reheating scenarios,

    L. E. Padilla, J. C. Hidalgo, K. A. Malik, and D. Mulryne, “Detecting the stochastic gravitational wave background from primordial black holes in slow-reheating scenarios,” JCAP, vol. 12, p. 011, 2024

  20. [28]

    Inflaton clusters and inflaton stars,

    J. C. Niemeyer and R. Easther, “Inflaton clusters and inflaton stars,” JCAP, vol. 07, p. 030, 2020

  21. [29]

    Primordial black hole formation during slow-reheating: a review,

    L. E. Padilla, J. C. Hidalgo, T. D. Gomez-Aguilar, K. A. Malik, and G. German, “Primordial black hole formation during slow-reheating: a review,” Front. Astron. Space Sci., vol. 11, p. 1361399, 2024

  22. [30]

    Relativistic equations for adiabatic, spherically symmetric gravitational collapse,

    C. W. Misner and D. H. Sharp, “Relativistic equations for adiabatic, spherically symmetric gravitational collapse,” Phys. Rev., vol. 136, pp. B571–B576, Oct 1964

  23. [31]

    Dynamics of pri- mordial black hole formation,

    J. C. Niemeyer and K. Jedamzik, “Dynamics of pri- mordial black hole formation,” Phys. Rev. D, vol. 59, p. 124013, 1999

  24. [32]

    Simulation of primordial black hole forma- tion using pseudo-spectral methods,

    A. Escriv` a, “Simulation of primordial black hole forma- tion using pseudo-spectral methods,” Phys. Dark Univ., vol. 27, p. 100466, 2020

  25. [33]

    Computations of primordial black hole formation,

    I. Musco, J. C. Miller, and L. Rezzolla, “Computations of primordial black hole formation,” Class. Quant. Grav., vol. 22, pp. 1405–1424, 2005

  26. [34]

    Identifying the most crucial parameters of the initial curvature profile for primordial black hole formation,

    T. Nakama, T. Harada, A. G. Polnarev, and J. Yokoyama, “Identifying the most crucial parameters of the initial curvature profile for primordial black hole formation,” JCAP, vol. 01, p. 037, 2014

  27. [36]

    Black hole formation in the friedmann universe: Formulation and computation in nu- merical relativity,

    M. Shibata and M. Sasaki, “Black hole formation in the friedmann universe: Formulation and computation in nu- merical relativity,” Physical Review D, vol. 60, Sept. 1999

  28. [37]

    Formalism for Primordial Black Hole Formation in Spherical Symme- try,

    J. Bloomfield, D. Bulhosa, and S. Face, “Formalism for Primordial Black Hole Formation in Spherical Symme- try,” 4 2015

  29. [38]

    Self- consistent initial conditions for primordial black hole formation,

    A. Polnarev, T. Nakama, and J. Yokoyama, “Self- consistent initial conditions for primordial black hole formation,” Journal of Cosmology and Astroparticle Physics, vol. 2012, p. 027–027, Sept. 2012

  30. [39]

    Identifying the most crucial parameters of the ini- tial curvature profile for primordial black hole forma- tion,

    T. Nakama, T. Harada, A. Polnarev, and J. Yokoyama, “Identifying the most crucial parameters of the ini- tial curvature profile for primordial black hole forma- tion,” Journal of Cosmology and Astroparticle Physics, vol. 2014, p. 037–037, Jan. 2014

  31. [40]

    Causal nature and dynamics of trapping horizons in black hole collapse,

    A. Helou, I. Musco, and J. C. Miller, “Causal nature and dynamics of trapping horizons in black hole collapse,” Classical and Quantum Gravity, vol. 34, p. 135012, June 2017

  32. [41]

    Causal Nature and Dynamics of Trapping Horizons in Black Hole Collapse,

    A. Helou, I. Musco, and J. C. Miller, “Causal Nature and Dynamics of Trapping Horizons in Black Hole Collapse,” Class. Quant. Grav., vol. 34, no. 13, p. 135012, 2017

  33. [42]

    Towards a singularity-proof scheme in numerical relativity,

    E. Seidel and W.-M. Suen, “Towards a singularity-proof scheme in numerical relativity,”Phys. Rev. Lett., vol. 69, pp. 1845–1848, Sep 1992

  34. [43]

    Horizon boundary condition for black hole space- times,

    P. Anninos, G. Daues, J. Mass´ o, E. Seidel, and W.-M. Suen, “Horizon boundary condition for black hole space- times,” Phys. Rev. D, vol. 51, pp. 5562–5578, May 1995

  35. [44]

    A new approach for simulating PBH forma- tion from generic curvature fluctuations with the Misner- Sharp formalism,

    A. Escriv` a, “A new approach for simulating PBH forma- tion from generic curvature fluctuations with the Misner- Sharp formalism,” 4 2025

  36. [45]

    Comparison of accurate methods for the integration of hyperbolic equations,

    H.-O. Kreiss and J.Oliger, “Comparison of accurate methods for the integration of hyperbolic equations,” Tellus, vol. 24, pp. 199–125, 1972

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.