Pith. Sign in

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

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 7 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDgr-qc/9905064 (version 1)
Pith IDRVESMHCO, Pith's permanent record of this paper and its review
CitationMasaru Shibata, Misao Sasaki. "Black hole formation in the Friedmann universe: Formulation and computation in numerical relativity." arXiv:gr-qc/9905064 (1999). https://pith.science/paper/RVESMHCO
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

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

The paper asks how large a primordial density fluctuation must be, in the very early universe, before it collapses into a black hole. Its answer is that there is no single number: the threshold depends on the shape of the fluctuation's neighbourhood. A peak sitting in an average-density flat universe needs a metric distortion of about $1.79$ at the initial time, while a peak surrounded by a low-density region can collapse from as little as about $1.43$. To reach that conclusion the authors build a single computational scheme that covers both ends of the problem, an analytic treatment while the fluctuation is larger than the Hubble horizon and a full numerical relativity simulation once it is well inside. The practical stake is that primordial black hole abundance is inferred from the high-amplitude tail of a Gaussian distribution, so a shift in the threshold of a few tenths changes the predicted number of black holes by orders of magnitude.

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

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

  • 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.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 5,137 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.04 · states limitation 0.60 · human participants 0.01 · controlled comparison 0.20 · randomized assignment 0.01 · states quantitative result 0.98 · compares with prior results 0.10 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

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

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
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 kind of paper the classifier reports this is: simulation study 0.98 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. FIG. 1. The spectrum shapes of the density fluctuation
Figure 4. FIG. 4. The same as Fig. 3, but for
Figure 7. It is found that the maximum value C(r)max is about 0.4 irrespective of σ. We note that for each σ, if the maximum value is large than the value shown in

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

Evidence: substantial. 184 of the citing papers share no author with it, and 36 of those carry a confident label from Pith's citation classifier.

Pith read 17 of the 374 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 3 of 23 papers with this label
How the literature uses it

164 arXiv papers cite this in 324 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 238 sentences without a confident label

Contests it 3 papers
Uses its method or tool 18 papers
Uses its data or values 13 papers
Compares with it 1 paper
Cites it as support 1 paper
Cites it as background 42 papers

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.