Vacuum Critical Collapse in Axisymmetric Spacetimes

• Vacuum critical collapse is the study of black-hole threshold phenomena in matter-free Einstein equations, focusing on axisymmetric gravitational-wave initial data.
• It exhibits key characteristics such as approximate discrete self-similarity, power-law scaling of black-hole mass, and recurring curvature echoes near threshold.
• Recent numerical studies show robust threshold scaling while highlighting the absence of a single universal critical solution across different initial data families.

Searching arXiv for papers on vacuum critical collapse and closely related threshold phenomena. arXiv search query: "vacuum critical collapse axisymmetric gravitational waves critical phenomena" Vacuum critical collapse is the study of the black-hole threshold in the matter-free Einstein equations, most prominently for axisymmetric gravitational-wave initial data, because in $3+1$ dimensions there is no nontrivial spherically symmetric vacuum collapse. In the standard critical-collapse framework, one considers a one-parameter family of smooth asymptotically flat initial data labeled by $p$, with dispersal for one range of $p$ and black-hole formation for another, and a threshold value $p_*$ separating the two. The classical paradigm predicts universality, power-law scaling such as $M_{\rm BH}\sim C(p-p_*)^\gamma$, and self-similar threshold solutions, often discretely self-similar (DSS). In vacuum axisymmetry, however, the modern picture is more qualified: critical behavior is supported by several numerical studies, but a single firmly established universal critical solution is not (Gundlach et al., 10 Jul 2025, Baumgarte et al., 2023).

1. Conceptual framework and historical benchmark

In the dynamical-systems description of critical collapse, the black-hole threshold is a codimension-one boundary in phase space, and a critical solution is an attractor within that boundary. If it has exactly one unstable mode, then near-threshold evolutions first approach that solution and only later depart along the unique growing direction. For Type II collapse, the threshold solution is self-similar, there is no intrinsic mass scale in the near-critical dynamics, and the black-hole mass obeys a scaling law. For DSS, the defining relation is

$\Phi_* g_{ab}=e^{-2\Delta} g_{ab},$

with $\Delta$ the echoing period (Gundlach et al., 10 Jul 2025).

The foundational numerical benchmark for the vacuum problem is the work of Abrahams and Evans, who studied axisymmetric vacuum gravitational waves using a partially constrained formulation, maximal slicing, and quasi-isotropic spatial gauge. They reported power-law scaling of black-hole mass, $\gamma \simeq 0.36$, tentative evidence for DSS echoing with $\Delta \simeq 0.6$, and some evidence for universality from a second Teukolsky-wave family. Later work, including recent reviews and reanalyses, treats those results as historically central but not yet decisively reproduced in the original setting (Gundlach et al., 10 Jul 2025, Rostworowski, 21 Mar 2025).

The modern significance of the subject lies precisely in this tension. Vacuum critical collapse is the cleanest setting in which to ask what in critical phenomena is intrinsic to gravity rather than to a matter model, yet it is also the setting in which the standard universal-attractor picture is least secure. Current evidence supports threshold scaling and repeated near-threshold structure for individual families, but it also suggests that the critical exponent and echoing period are not universal across families (Gundlach et al., 10 Jul 2025).

2. Classical axisymmetric vacuum setting

Most contemporary work concerns axisymmetric, asymptotically flat vacuum spacetimes, often in the twist-free sector. In the formulation revisited by Rinne, the physical system is vacuum axisymmetric gravitational collapse in the vanishing-twist sector, meaning axisymmetry generated by a hypersurface-orthogonal Killing field. In perturbative language, this is a consistent truncation to the polar or even-parity sector, so there is only one genuine gravitational-wave degree of freedom. The spacetime is written in quasi-isotropic spherical-polar coordinates $(t,r,\theta,\varphi)$, with $p$0, and supplemented by maximal slicing (Rostworowski, 21 Mar 2025).

That formulation is fully constrained. The evolved pair is $p$1, while the remaining variables are obtained on each slice from the momentum constraints, the nonlinear Hamiltonian constraint, a linear elliptic slicing equation, and linear equations for the shift. A central technical point is the use of conformally scaled extrinsic-curvature variables rather than bare $p$2; with the scaled variables, the Hamiltonian constraint decouples properly and is expected to have a unique solution, whereas unscaled choices can lead to nonuniqueness and related pathologies (Rostworowski, 21 Mar 2025).

Several initial-data constructions are used across the literature. Brill-wave data are time-symmetric and encode the wave in the spatial metric, for example through

$p$3

with amplitude $p$4 as the tuning parameter (Baumgarte et al., 2023). Teukolsky-wave and Nakamura-wave families instead begin from linearized gravitational-wave solutions and then solve the nonlinear constraints. In the 2026 Nakamura study, the wave content is encoded in the extrinsic curvature rather than the spatial curvature, which simplifies the nonlinear dressing and allowed slightly better fine-tuning than in the authors’ previous Teukolsky-wave study (Baumgarte et al., 25 Jun 2026).

The threshold itself is defined in the usual way. One chooses a one-parameter family of quadrupolar or otherwise axisymmetric vacuum-wave data, with subcritical data dispersing and supercritical data forming a black hole. In the axisymmetric vacuum problem, however, the strong-field region need not remain centered at the origin, and collapse may occur on the symmetry axis away from the center, or in multiple centers. This feature already distinguishes the vacuum problem from the canonical spherical matter models (Gundlach et al., 10 Jul 2025).

3. Numerical evidence and the question of universality

Recent numerical work substantially improved the evidentiary basis for vacuum critical collapse. A major cross-code comparison used three independent implementations—a BSSN finite-difference code based on the Einstein Toolkit, a first-order generalized harmonic pseudo-spectral code, and a spherical-polar BSSN code—and found excellent agreement on the existence of critical behavior in axisymmetric vacuum gravitational-wave collapse. For negative-amplitude Brill waves, the three codes agreed on

$p$5

and on a threshold solution that is approximately, but not exactly, DSS with an accumulation point at the center; the estimated echoing period was $p$6. For positive-amplitude Brill waves, by contrast, they found no evidence for a central DSS solution, and the near-threshold structure was qualitatively different. The principal conclusion was explicit: there is no single, universal critical solution for vacuum gravitational-wave collapse (Baumgarte et al., 2023).

A complementary study based on curvature invariants argued for a weaker form of universality. Near-threshold vacuum spacetimes generated repeated, approximate, scaled, irregular echoes of the same local curvature structure across several asymptotically flat initial-data families. The principal invariants were the Kretschmann scalar,

$p$7

and the axis-adapted scalar

$p$8

with on-axis relations

$p$9

That work did not establish a unique global DSS vacuum critical solution, nor a family-independent critical exponent, but it did support a recurring local spacetime template in the strong-curvature region (Ledvinka et al., 2021).

Further evidence against universality came from comparisons of different vacuum multipoles. Simulations of quadrupolar and hexadecapolar Teukolsky-wave initial data found that both threshold solutions appeared to display at least approximate DSS with an accumulation event at the center, but the hexadecapolar threshold solution appeared distinct from the quadrupolar one. The quadrupolar family showed characteristics consistent with earlier quadrupolar studies, while the hexadecapolar family showed a much shorter approximate echoing period, around $p$0, compared with $p$1 for the quadrupolar case (Baumgarte et al., 2023).

The 2026 Nakamura-wave study strengthened that picture rather than reversing it. Using new families of vacuum initial data and slightly better fine-tuning, it found one additional echo in the approximately self-similar threshold solution. The conclusions were again consistent with earlier studies: threshold solutions are approximately DSS, the self-similarity is not exact, and there is no evidence for a unique critical solution. The paper also highlighted recurring geometric features shared across families, including alternating maxima in the direction of the poles and the equator (Baumgarte et al., 25 Jun 2026).

4. Local geometry, diagnostics, and nearby wave-driven systems

A persistent theme in the vacuum literature is that the most useful diagnostics are invariant or nearly invariant curvature measures and adapted coordinates rather than coordinate radii or raw gauge variables. For DSS tests in the three-code Brill-wave comparison, the key adapted time was

$p$2

where $p$3 is proper time at the center and $p$4 is the accumulation time, together with outgoing null geodesics labeled by an affine parameter $p$5. The rescaled curvature $p$6 was then examined in $p$7 coordinates. For the $p$8 Brill family this showed approximate periodicity; for $p$9 it did not (Baumgarte et al., 2023).

The spatial structure of the strong-field region is also diagnostically important. The review literature emphasizes that subcritical maxima of the Kretschmann scalar can lie on the symmetry axis but away from the origin, and that supercritical evolutions can form disjoint apparent horizons around separated curvature peaks. This is one reason the standard picture of a single centered self-similar attractor is less constraining in vacuum than in spherical matter collapse (Gundlach et al., 10 Jul 2025).

Several nonvacuum but wave-driven systems have clarified why this may happen. In axisymmetric Einstein-Maxwell collapse, fine-tuning to threshold gave approximate power-law scaling

$p_*$0

and approximate DSS with $p_*$1, but neither DSS nor universality was exact. Dipole and quadrupole electromagnetic waves then showed different approximate echoing periods and different critical exponents, with quadrupole data producing two separate collapse centers on the symmetry axis rather than one central center. The authors explicitly discussed the implications for vacuum gravitational-wave collapse, arguing that the absence of a spherically symmetric radiative critical solution may naturally weaken or destroy global uniqueness (Baumgarte et al., 2019, Mendoza et al., 2021).

A related bridge is provided by twist-free axisymmetric collapse of a massless complex scalar field. In spherical symmetry, that system appears to share the same critical spacetime metric as the Choptuik solution, with

$p_*$2

But as the degree of asphericity increases, the extracted $p_*$3 and $p_*$4 drift, periodicity becomes less clean, and at sufficiently high asphericity the center of collapse bifurcates into two centers on the symmetry axis away from the origin. The authors interpret this as a competition between matter-driven and gravitational-wave-driven collapse, with increasing Weyl or gravitational-wave contribution making the dynamics more vacuum-like. This suggests that at least some of the anomalies seen in vacuum collapse are part of a broader class of highly aspherical threshold phenomena (Marouda et al., 2024).

5. Reproducibility problem and formulation obstacles

A major recent development is the recognition that the vacuum critical-collapse problem is as much a formulation problem as a brute-force resolution problem. Rinne’s re-entry into the original Abrahams–Evans problem was explicitly not a new demonstration of vacuum critical collapse. Its goal was to reconstruct the original constrained axisymmetric formulation faithfully, derive the corresponding linearized theory in a form suited to code validation, and diagnose the specific numerical obstacles that prevent stable long-time evolutions near the black-hole threshold (Rostworowski, 21 Mar 2025).

One of the strongest results of that study is the linearized benchmark around Minkowski spacetime,

$p_*$5

with master scalar obeying

$p_*$6

For $p_*$7, explicit formulas were derived for the metric variables in the Abrahams–Evans gauge. These exact solutions provide a stringent weak-field benchmark, and the corresponding pseudospectral code reproduced them with spectral accuracy in space and fourth-order convergence in time (Rostworowski, 21 Mar 2025).

The same work isolated the main failure mode near threshold. In the nonlinear equations, mixed derivative terms such as

$p_*$8

are analytically finite because the relevant scalars behave like $p_*$9 near the origin, but numerically the sequence “differentiate, then divide by $M_{\rm BH}\sim C(p-p_*)^\gamma$0” is not spectrally benign. In spectral space this seeds a growing high-frequency plateau of numerical contamination. In the linearized regime the low modes remain accurate and the contamination is largely harmless; in long-lived nonlinear near-threshold evolutions, however, the high-mode noise couples back into the low modes and destroys the evolution on exactly the timescales needed for critical-collapse studies. The resulting conclusion is direct: simply rerunning Abrahams and Evans with modern hardware is not straightforward (Rostworowski, 21 Mar 2025).

A related formulation paper, motivated explicitly by difficulties encountered in the vacuum setting, addressed two generic obstacles: coordinate singularities and large constraint violations. Working in first-order generalized harmonic gauge, it derived a necessary condition for a gauge to respect DSS and proposed a DSS-compatible gauge source

$M_{\rm BH}\sim C(p-p_*)^\gamma$1

together with a modified reduction-constraint damping prescription

$M_{\rm BH}\sim C(p-p_*)^\gamma$2

In spherical scalar-field tests these changes improved threshold tuning from 6 up to 11 digits and clearly observed up to 3 echoes (Cors et al., 2023). This does not by itself solve the vacuum problem, but it suggests that gauge design and reduction-constraint control are likely central to any successful high-precision vacuum code.

6. Broader uses of the term and semiclassical extensions

The phrase “vacuum critical collapse” is sometimes used more broadly than the axisymmetric vacuum gravitational-wave problem. One neighboring usage concerns thin-wall vacuum bubbles. In that setting, a spherical thin-wall bubble separating two constant-vacuum-energy regions obeys

$M_{\rm BH}\sim C(p-p_*)^\gamma$3

or, after rescaling,

$M_{\rm BH}\sim C(p-p_*)^\gamma$4

The critical structure is then the separatrix at the single maximum of $M_{\rm BH}\sim C(p-p_*)^\gamma$5, together with distinguished masses such as $M_{\rm BH}\sim C(p-p_*)^\gamma$6, $M_{\rm BH}\sim C(p-p_*)^\gamma$7, $M_{\rm BH}\sim C(p-p_*)^\gamma$8, and $M_{\rm BH}\sim C(p-p_*)^\gamma$9. This is a threshold problem in vacuum bubble dynamics, but it is not Choptuik-type critical collapse with universal self-similar PDE attractors (Ng et al., 2010).

A different extension concerns semiclassical backreaction near the classical black-hole threshold. In analytically tractable Einstein-scalar critical spacetimes, one-loop vacuum polarization in a regular Boulware-like state produces a universal quantum growing mode with exponent

$\Phi_* g_{ab}=e^{-2\Delta} g_{ab},$0

In the examples studied, this quantum mode shifts the threshold, generates a finite mass gap, and changes the effective behavior from classical Type II to quantum-modified Type I, thereby cloaking the classical naked singularity behind a horizon (Tomašević et al., 3 Sep 2025). This is not classical vacuum gravitational-wave collapse, but it shows that quantum vacuum effects can qualitatively change near-threshold dynamics.

In the standard $\Phi_* g_{ab}=e^{-2\Delta} g_{ab},$1 axisymmetric vacuum problem, the present status remains unsettled but sharply defined. Critical behavior in pure gravity is now numerically real enough to show threshold scaling and repeated near-threshold structure, yet the field does not presently support a single universal axisymmetric vacuum critical solution. What is established is approximate DSS for some families, approximate scaling for individual families, and recurring local curvature structure; what remains open is whether deeper fine-tuning, different gauges, or different formulations will reveal a stricter universality class, or whether family dependence is itself a genuine property of vacuum critical collapse (Gundlach et al., 10 Jul 2025)