Critical Phenomena in Gravitational Collapse

Critical phenomena organize the boundary between black-hole formation and a competing endpoint. Near that boundary, nonlinear evolution can forget most details of the initial data, approach a universal intermediate solution, and depart along one unstable mode. This chapter derives the lifetime and mass laws, explains continuous and discrete self-similarity, and marks the boundary between established spherical results and the open nonspherical problem.

Helpful background. The Einstein equation supplies the constrained nonlinear flow, and Ordinary Differential Equations supplies eigenmodes, growth rates, and asymptotic scaling.

Enter this chapter

The common setup is a smooth one-parameter family of regular initial data that crosses a black-hole threshold transversely at $p=p_*$. The matter model, asymptotics, symmetry class, and endpoint criterion must be fixed. The phase-space picture assumes a gauge prescription and physical mode spectrum; gauge and constraint-violating directions do not count as instabilities.

The pages form one argument:

1. Type I Critical Collapse derives logarithmic lifetime scaling near a finite-scale metastable solution.
2. Type II Critical Collapse derives a shrinking departure scale and arbitrarily small classical black holes.
3. Continuous and Discrete Self-Similarity distinguishes homothetic fixed points from echoing limit cycles.
4. Mass Scaling and Critical Exponents connects observable powers to the growing eigenvalue and explains reliable fitting.
5. Universality from Perturbation Modes derives codimension and initial-family independence.
6. Critical Collapse Beyond Spherical Symmetry assesses the dated evidence for angular modes, rotation, and vacuum waves.

One phase-space picture, two types

Let $Z_*$ lie on the boundary between basins and have one growing physical mode,

$$\delta Z=C_0(p)e^{\lambda_0 T}Z_0+\text{decaying modes}, \qquad C_0(p)\simeq\kappa(p-p_*).$$

Here $Z_0$ is the growing physical eigenfunction, $C_0(p)$ its family-dependent amplitude, and $\lambda_0>0$ its growth exponent after gauge and constraint directions are removed. Take $p$ to be dimensionless locally; otherwise replace $p-p_*$ inside logarithms by $(p-p_*)/p_0$ for a fixed parameter scale. For type I, set $T=t/t_0$ using a declared physical time scale $t_0$, and let $Z_*$ be stationary or periodic. Departure at fixed amplitude gives

$$t_p=-\frac{t_0}{\lambda_0}\ln\lvert p-p_*\rvert+\text{constant}.$$

The intermediate solution has finite mass, so the black-hole branch can have a mass gap.

For type II, $T=\tau$ is logarithmic scale time and $Z_*$ is continuously or discretely self-similar. Departure generates

$$L_*\propto\lvert p-p_*\rvert^{1/\lambda_0}, \qquad GM_{\rm BH}\propto(p-p_*)^\gamma, \qquad \gamma=\lambda_0^{-1}.$$

The same matter theory can display both types in different regions if a dimensionful scale is relevant in one regime and negligible in another. Gundlach and Martín-García derive the common phase-space structure and the distinction in Gundlach and Martín-García 2007, §§2.1–2.4, pp. 4–8.

Self-similarity and measurable scaling

Continuous self-similarity is generated by a homothety $\mathcal L_\xi g=2g$. Discrete self-similarity repeats only after $\tau\mapsto\tau+\Delta$ and a corresponding metric rescaling. In scale-quotiented phase space, these appear as a fixed point and a limit cycle. Neither symmetry is usually global: regular asymptotically flat data approach it in a shrinking core for a finite number of echoes.

Observable dimensions then turn $L_*$ into a hierarchy of exponents. Area scales with $2\gamma$, curvature with negative powers, and quantities controlled by additional decaying modes can carry other exponents. DSS adds periodic fine structure to the logarithmic mass law. A consistent analysis compares supercritical mass, subcritical curvature, local echoes, and an independently calculated eigenvalue.

Choptuik’s scalar-field calculation established the prototype in Choptuik 1993, pp. 9–12. Koike, Hara, and Adachi connected the perfect-fluid exponent to a perturbative eigenvalue in Koike, Hara, and Adachi 1995, pp. 5170–5173.

What universality includes

Decaying modes remove dependence on the initial pulse shape, leaving family-specific threshold locations, amplitudes, and DSS phases but a shared intermediate profile and exponent. This is universality across families that reach the same critical solution. It is not a claim that all equations of state, matter fields, or symmetry classes share one value of $\gamma$.

The phase space must be stated. Each additional independent real growing amplitude raises the codimension by one; a nonspherical angular multiplet can contribute several amplitudes despite sharing one eigenvalue. Conversely, the massless scalar spherical DSS solution has strong perturbative evidence that its nonspherical modes decay Martín-García and Gundlach 1999, §VIII and Table I, pp. 15–19, PDF. Axisymmetric vacuum collapse and fully three-dimensional universality remain less settled.

The frontier evidence cutoff is 2026-08-12. Updating it requires a broader stability theorem, a converged counterexample, independent reproduction of a disputed critical solution, or discovery of a new one-mode nonspherical attractor.

Testing critical scaling

A numerical threshold study must converge $p_*$, repeat several initial-data families, resolve the same echoes at at least three resolutions, and separate truncation, finite-domain, and fitting error. It should quote the actual scaling window and show sensitivity to shifting $p_*$. Dimensionless profiles and invariant curvature or horizon quantities are stronger evidence than a single log–log slope.

At the far end of the window, family-specific decaying modes matter. At the near end, grid scale, roundoff, another physical length, or quantum physics intervenes. Classical type II scaling toward zero mass is therefore an asymptotic statement within a declared continuum model, not a trans-Planckian prediction.

Review the chapter

A data family shows $M_{\rm BH}\propto(p-p_*)^{0.36}$ and dimensionless lifetime $t_p/t_0\propto-2.78\ln\lvert p-p_*\rvert$. Can both laws arise from the same single eigenvalue in one regime?

Solution

No. The mass exponent $\gamma=0.36$ implies a logarithmic-scale eigenvalue $\lambda_{\rm II}=1/0.36\simeq2.78$. Using this chapter’s dimensionless type I clock $T=t/t_0$, the lifetime coefficient instead satisfies $t_p/t_0=-\lambda_{\rm I}^{-1}\ln\lvert p-p_*\rvert+\mathrm{constant}$, so a coefficient $2.78$ implies $\lambda_{\rm I}\simeq0.360$. Its dimensional growth rate is $\lambda_{\rm I}/t_0$. These are reciprocal coefficients of different kinds, not equal eigenvalues. Moreover, physical time and logarithmic scale time are different clocks: a finite-scale type I state has a mass gap, whereas type II has a shrinking core. Check the intermediate solution and its clock before comparing any fitted numbers.

References

• Choptuik, Matthew W. “Universality and Scaling in Gravitational Collapse of a Massless Scalar Field.” Physical Review Letters 70 (1993): 9–12. doi:10.1103/PhysRevLett.70.9.
• Gundlach, Carsten, and José M. Martín-García. “Critical Phenomena in Gravitational Collapse.” Living Reviews in Relativity 10 (2007): 5. doi:10.12942/lrr-2007-5.
• Koike, Tatsuhiko, Takashi Hara, and Satoshi Adachi. “Critical Behavior in Gravitational Collapse of Radiation Fluid: A Renormalization Group Analysis.” Physical Review Letters 74 (1995): 5170–5173. doi:10.1103/PhysRevLett.74.5170.
• Martín-García, José M., and Carsten Gundlach. “All Nonspherical Perturbations of the Choptuik Spacetime Decay.” Physical Review D 59 (1999): 064031. doi:10.1103/PhysRevD.59.064031.

Dated Research assessment. For open questions and recent results, follow the link below. Each assessment states the date through which it considers results and the limits of its conclusions.

• Universality in Critical Gravitational Collapse