REVIEW
Black hole formation in the Friedmann universe: Formulation and computation in numerical relativity
General Relativity and Quantum Cosmology · Submitted 18 May 1999 · 18 references
Pith's verdict: The constant-mean-curvature formulation is the real contribution and it holds up; the profile-dependence physics is right, but the abstract leans on the variable that maximizes the spread and on a mass asymptote the authors themselves say they didn't reach.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
Whether a dense peak becomes a black hole depends on what surrounds it, not just how tall it is.
The reading
What carries the argument
The load-bearing device is the constant mean curvature slicing condition $K=-3ȧ/a$, which makes the equations for the lapse $α$ and the conformal factor $ψ$ resemble their maximal-slicing counterparts in asymptotically flat space and lets the authors factor the cosmological expansion out of the dynamical variables. This is what allows one set of variables, one gauge, and one code to run continuously from the epoch when a fluctuation is far outside the horizon, where the long-wavelength expansion gives an exact analytic answer, to the epoch when it is far inside and a black hole forms. The analytic piece supplies the initial data; the numerical piece supplies the threshold and the final mass, with the Hamiltonian constraint and the Kodama mass used as accuracy checks on the apparent-horizon mass.
Extended reading notes
Core claim
The authors set up the Einstein and hydrodynamic equations in a Friedmann background using constant mean curvature time-slicing, chosen so that the geometric equations take a form close to maximal slicing and so that the same variables remain valid from the superhorizon epoch through collapse. In the long-wavelength limit they obtain an analytic solution in which the density contrast and the lapse perturbation are proportional to $t^{2-4/3Γ}$ while the conformal factor $ψ$ is frozen to leading order, and they use that solution to build initial data for a spherically symmetric, radiation-dominated ($Γ=4/3$) simulation. Scanning a two-parameter family of profiles, a Gaussian peak of amplitude $C_δ$ minus a broader negative Gaussian of width $σ$, they find that a black hole forms only when the initial metric perturbation is already non-linear, with $ψ_0$ at $t=0$ above roughly $1.4$ in the most favourable case. The threshold is not universal: it is about $1.79$ when the peak is surrounded by a flat, average-density universe ($σ=∞$), falls for smaller $σ$ where the peak sits inside a low-density region, and reaches about $1.43$ in the flat-spacetime time-symmetric limit, which they treat as approximately the smallest possible threshold. The final apparent-horizon mass varies just as strongly, exceeding $8r_0$ for $σ=∞$ but reaching only about $1.1r_0$ for $σ=2$. As a profile-independent alternative, they define a compaction function $C(r)=δM_K/[rψ(r)^2a]$ from the Kodama mass and report that its maximum is about $0.4$ at threshold for every $σ$ they tried.
Load-bearing premise
The claim that the threshold depends on the surrounding density rests on one two-parameter family of profiles, a Gaussian peak minus a broader negative Gaussian of width $σ$ with the scale $r_0$ fixed to unity, and the paper does not derive that family from an inflationary spectrum or show that other compensated shapes give the same spread in thresholds.
Editorial extensions
If this is right
- Threshold estimates for primordial black hole formation that use a single amplitude criterion, independent of the surrounding density field, will misstate the abundance by large factors.
- The number of primordial black holes predicted from an inflationary power spectrum must be computed from a joint distribution of amplitude and local environment, not from the one-point amplitude distribution alone.
- Final black hole masses span a wide range for the same peak amplitude, so a claimed mass function cannot be inferred from the threshold alone.
- The compaction function maximum, about $0.4$ at threshold across the profiles tested, offers a profile-independent collapse criterion in spherical symmetry.
- The same slicing and variable choice can be carried over to non-spherical simulations, which the authors name as the next step.
Reading between the lines
- If the environment matters this much in spherical symmetry, then in three dimensions the tidal field of neighbouring peaks is a second environmental variable of comparable weight, and the threshold should be treated as a functional of the surrounding density field rather than a number.
- The reported near-constancy of the compaction maximum suggests the right invariant is a mass-excess-to-radius ratio at horizon crossing; a general-relativistic version of a Jeans-type criterion written in that variable would unify the $σ$ dependence the authors see.
- Because the threshold shift is modest while the abundance is exponentially sensitive to it, the dominant uncertainty in primordial black hole dark matter bounds may be the shape of the fluctuation profile rather than the amplitude of the power spectrum.
- The authors' remark that a peak surrounded by a void collapses most easily implies that primordial black hole formation should correlate spatially with underdense regions, which is in principle a signature distinguishable from a purely local threshold model.
Original abstract
We study formation of black holes in the Friedmann universe. We present a formulation of the Einstein equations under the constant mean curvature time-slicing condition. Our formalism not only gives us the analytic solution of the perturbation equations for non-linear density and metric fluctuations on superhorizon scales, but also allows us to carry out a numerical relativity simulation for black hole formation after the scale of the density fluctuations is well within the Hubble horizon scale. We perform a numerical simulation of spherically symmetric black hole formation in the radiation-dominated, spatially flat background universe for a realistic initial condition supplied from the analytic solution. It is found that the initial metric perturbation has to be non-linear (the maximum value of 3D conformal factor $\psi_0$ at $t=0$ should be larger than $\sim 1.4$) for a black hole to be formed, but the threshold amplitude for black hole formation and the final black hole mass considerably depend on the initial density (or metric) profile of the perturbation: The threshold value of $\psi_0$ at $t=0$ for formation of a black hole is smaller for a high density peak surrounded by a low density region than for that surrounded by the average density region of the flat universe. This suggests that it is necessary to take into account the spatial correlation of density fluctuations in the study of primordial black hole formation.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
We study formation of black holes in the Friedmann universe. We present a formulation of the Einstein equations under the constant mean curvature time-slicing condition. Our formalism not only gives us the analytic solution of the perturbation equations for non-linear density and metric fluctuations on superhorizon scales, but also allows us to carry out a numerical relativity simulation for black hole formation after the scale of the density fluctuations is well within the Hubble horizon scale. We perform a numerical simulation of spherically symmetric black hole formation in the radiation-dominated, spatially flat background universe for a realistic initial condition supplied from the analytic solution. It is found that the initial metric perturbation has to be non-linear (the maximum value of 3D conformal factor $\psi_0$ at $t=0$ should be larger than $\sim 1.4$) for a black hole to be formed, but the threshold amplitude for black hole formation and the final black hole mass considerably depend on the initial density (or metric) profile of the perturbation: The threshold value of $\psi_0$ at $t=0$ for formation of a black hole is smaller for a high density peak surrounded by a low density region than for that surrounded by the average density region of the flat universe. This suggests that it is necessary to take into account the spatial correlation of density fluctuations in the study of primordial black hole formation.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | We study formation of black holes in the Friedmann universe. | contribution 0.69 | 0.90 | 0.03 | asserted 0.70 |
| s1 | We present a formulation of the Einstein equations under the constant mean curvature time-slicing condition. | method 0.71 | 0.93 | 0.03 | asserted 1.00 |
| s2 | Our formalism not only gives us the analytic solution of the perturbation equations for non-linear density and metric fluctuations on superhorizon scales, but also allows us to carry out a numerical relativity simulation for black hole formation after the scale of the density fluctuations is well within the Hubble horizon scale. | contribution 0.80 | 0.93 | 0.10 | asserted 1.00 |
| s3 | We perform a numerical simulation of spherically symmetric black hole formation in the radiation-dominated, spatially flat background universe for a realistic initial condition supplied from the analytic solution. | method 0.79 | 0.96 | 0.06 | asserted 0.99 |
| s4 | It is found that the initial metric perturbation has to be non-linear (the maximum value of 3D conformal factor $\psi_0$ at $t=0$ should be larger than $\sim 1.4$) for a black hole to be formed, but the threshold amplitude for black hole formation and the final black hole mass considerably depend on the initial density (or metric) profile of the perturbation: The threshold value of $\psi_0$ at $t=0$ for formation of a black hole is smaller for a high density peak surrounded by a low density region than for that surrounded by the average density region of the flat universe. | result 1.00 | 0.96 | 0.63 | asserted 1.00 |
| s5 | This suggests that it is necessary to take into account the spatial correlation of density fluctuations in the study of primordial black hole formation. | interpretation 1.00 | 0.66 | 0.63 | hedged possibility 0.99 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Showing 3 of the paper's 7 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:gr-qc/9905064 by the authors.
References
How later papers used it
192 later papers cite it
3 later papers by other authors are listed below: 3 used its method or tool.
Uses its method or tool 3 of 23 papers with this label
-
More precisely, we consider an optimized criterion for PBH formation which is given by a threshold for the compaction function proposed in Ref. [ref] .
PBH abundance from random Gaussian curvature perturbations and a local density threshold (2018) · Summary and Discussion · 0.89
-
Used the compaction function introduced in this paper as the criterion for primordial black hole formation.
Effect of nonlinearity between density and curvature perturbations on the primordial black hole formation (2019) · 0.80 · paraphrased: the licence does not allow quoting
-
Used the criterion proposed in this paper, based on the compaction function.
From Primordial Black Holes Abundance to Primordial Curvature Power Spectrum (and back) (2019) · 0.99 · paraphrased: the licence does not allow quoting
How the literature uses it
164 arXiv papers cite this in 324 sentences
Papers citing it · 238 sentences without a confident label
Contests it 3 papers
-
In this paper, we show that the Shibata-Sasaki compaction function introduced in Ref. [ref] is not what they intended.
Revisiting compaction functions for primordial black hole formation · Introduction · contest 0.87
-
However, it should be noted that, for the CMC slice, the contribution from the second term in the integrand of Eq. [ref] cannot be neglected[In Refs. [ref] , the second term in the integrand of Eq. [ref] has been erroneously dropped, and the authors got the wrong relation $𝒞_CMC=𝒞_SS$.
The basics of primordial black hole formation and abundance estimation · PBH formation criterion · contest 0.93
-
Disregard for this issue affected early work on this subject [ref] , as pointed out in Ref. [ref] .
Formalism for Primordial Black Hole Formation in Spherical Symmetry · Initial Conditions · contest 0.77
Uses its method or tool 18 papers
-
Therefore, in this section, we introduce the criterion for PBH formation assuming spherical symmetry originally proposed in Ref. [ref] .
PBH abundance from random Gaussian curvature perturbations and a local density threshold · perturbation variables and threshold for PBH formation · 0.84
-
This construction combines the long-wavelength treatment of nonlinear cosmological perturbations with Gaussian peak statistics and subsequent refinements accounting for the nonlinear curvature–density relation, the profile dependence of the collapse threshold, and critical scaling near the threshold [ref] .
Projected Constraints on Primordial Black Holes from Scalar-Induced Gravitational Waves with Taiji · PBH formation from the curvature spectrum · 0.94
-
Here as the criterion for the PBH formation, we use the Compaction function $𝒞(r)$ introduced in Ref. [ref] .
Effect of nonlinearity between density and curvature perturbations on the primordial black hole formation · Introduction · 0.80
Uses its data or values 13 papers
-
In the constant mean curvature slicing, the numerical simulation gives the threshold value of the PBH formation as $𝒞_m;CMC^th≃0.4$ during the radiation dominated eta [ref] .
Effect of nonlinearity between density and curvature perturbations on the primordial black hole formation · PBH formation criterion and Press-Schechter formalism · 1.00
-
We use $δ_c ≃0.42$ [ref] (see also Refs. [ref] ), $k=3.3$ [ref] and $γ=0.36$ [ref] .
Primordial Black Holes from Thermal Inflation · Primordial black holes · 1.00
-
Shibata and Sasaki [ref] obtained the threshold value $≃0.4$ for radiation $Γ=4/3$.
Revisiting compaction functions for primordial black hole formation · Shibata-Sasaki compaction function and its significance · 0.92
Compares with it 1 paper
-
Further, this analytical estimate of $δ_c$ is in rough agreement with the range specified from the numerical analysis in [ref] .
Primordial black holes dark matter from inflection point models of inflation and the effects of reheating · Critical value of the density contrast · 0.81
Cites it as support 1 paper
-
We do not include the effect of the critical collapse [ref] for simplicity because Ref. [ref] implies that it does not modify the PBH abundance so much.
NANOGrav results and LIGO-Virgo primordial black holes in axion-like curvaton model · Introduction · 0.97
Cites it as background 42 papers
-
Another quantity, the maximum value of the so-called compaction function has been proposed as the one giving a robust threshold of primordial black hole formation by Shibata and Sasaki [ref] through general relativistic numerical simulations.
Revisiting compaction functions for primordial black hole formation · Introduction · 0.84
-
After these pioneering works, the collapse of initial perturbations in the form of primordial curvature fluctuations, was numerically studied in Refs. [ref] , followed some years afterwards by an extensive numerical analysis by one of us using an explicit Lagrangian hydrodynamics code developed and used in Refs. [ref] and, recently, by other authors in Refs. [ref] .
From Primordial Black Holes Abundance to Primordial Curvature Power Spectrum (and back) · Primordial black holes formation · 0.91
-
Our present considerations do not apply to such situation.] A rough criterion for the formation of PBHs was first proposed by Carr [ref] , and much numerical work has been done in search of a more accurate criterion [ref] .
PBH abundance from random Gaussian curvature perturbations and a local density threshold · Introduction · 0.96
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.