pith. sign in

arxiv: 2606.27431 · v1 · pith:JR652WDHnew · submitted 2026-06-25 · 🌀 gr-qc

Critical collapse of vacuum spacetimes: Nakamura wave initial data

Pith reviewed 2026-06-29 01:46 UTC · model grok-4.3

classification 🌀 gr-qc
keywords critical collapsevacuum spacetimesgravitational wavesdiscrete self-similarityblack hole formationaxisymmetric initial dataNakamura waves
0
0 comments X

The pith

Simulations of critical collapse in vacuum gravitational waves find only approximate discrete self-similarity.

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

This paper examines critical phenomena in the collapse of axisymmetric vacuum gravitational waves by adopting new families of initial data constructed via Nakamura's approach. The data encode wave content in the extrinsic curvature, which simplifies satisfying the nonlinear constraint equations. Improved fine-tuning to the black hole threshold reveals one additional echo in the threshold solution. The results confirm that these solutions are approximately discretely self-similar but not exactly so, with no evidence for a unique critical solution, and they share common features such as alternating maxima toward the poles and equator.

Core claim

Adopting Nakamura wave initial data allows fine-tuning to the onset of black hole formation slightly better than before, making one more echo visible in the approximately self-similar threshold solution. The self-similarity remains inexact, no unique critical solution appears, and alternating maxima in the polar and equatorial directions recur across different threshold solutions.

What carries the argument

Nakamura wave initial data, which encode the gravitational wave content in the extrinsic curvature and are adjusted to satisfy the nonlinear Einstein constraints.

If this is right

  • Threshold solutions exhibit approximate discrete self-similarity.
  • The self-similarity is not exact.
  • No evidence supports the existence of a unique critical solution.
  • Different threshold solutions share the feature of alternating maxima toward the poles and the equator.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The degree of self-similarity may depend on the specific family of initial data chosen.
  • Further improvements in resolution could test whether the observed deviations from exact self-similarity remain stable.
  • Recurring features across families point to some aspects of the critical approach that are independent of the initial data construction.

Load-bearing premise

The numerical resolution, gauge choices, and evolution methods suffice to distinguish approximate from exact discrete self-similarity without introducing mimicking artifacts.

What would settle it

A higher-resolution simulation or different gauge choice that produces exact discrete self-similarity together with a unique critical solution would contradict the reported observations.

Figures

Figures reproduced from arXiv: 2606.27431 by Carsten Gundlach, David Hilditch, Thomas W. Baumgarte.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of the curvature invariants [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Scaling for our three families of initial data. For each [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The threshold solution for the Nakamura family [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The maximum curvature [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Contour plots in the [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Spacetime profiles of the curvature invariants [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Same as Figs [PITH_FULL_IMAGE:figures/full_fig_p012_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Spacetime profiles of the curvature invariants [PITH_FULL_IMAGE:figures/full_fig_p012_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Spacetime profiles of the curvature invariant [PITH_FULL_IMAGE:figures/full_fig_p013_12.png] view at source ↗
read the original abstract

We report on numerical simulations of critical phenomena in the collapse of axisymmetric vacuum gravitational waves, adopting families of initial data that, to the best of our knowledge, have not been used in this context before. Like Teukolsky waves, the data are based on linear wave solutions to the Einstein equations. We follow Nakamura's construction and encode the wave content in the extrinsic curvature rather than the spatial curvature, which leads to several simplifications when the data are "dressed up" so that they satisfy the nonlinear constraint equations. We are able to fine-tune these data to the onset of black hole formation slightly better than in our previous simulations, allowing us to observe and examine one more echo in the approximately self-similar threshold solution. Our findings are consistent with earlier studies: while we find threshold solutions that are approximately discretely self-similar, the self-similarity is not exact, and we find no evidence for a unique critical solution. We discuss common features between the different threshold solutions, including the appearance of alternating maxima in the direction of the poles and the equator.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript reports numerical simulations of critical collapse for axisymmetric vacuum gravitational waves using a new family of initial data based on Nakamura's construction, in which the wave content is encoded in the extrinsic curvature. The authors achieve slightly improved fine-tuning to the black-hole formation threshold relative to prior work, permitting observation of one additional echo. They conclude that the resulting threshold solutions are approximately but not exactly discretely self-similar, exhibit no evidence for a unique critical solution, and display common features such as alternating maxima along the poles and equator; these findings are stated to be consistent with earlier studies that used Teukolsky-wave data.

Significance. If the reported numerical observations hold, the work supplies an independent cross-check on the approximate discrete self-similarity of vacuum critical solutions using a distinct initial-data construction that simplifies the constraint-solving step. The modest advance in fine-tuning (one extra echo) is a concrete technical contribution that strengthens the empirical case against exact self-similarity and uniqueness. The manuscript also supplies the expected details on initial-data construction, constraint solving, evolution, and standard convergence checks, so the stress-test concern about numerical resolution and gauge artifacts does not appear to land.

minor comments (3)
  1. The abstract states consistency with prior results and improved tuning but supplies no details on numerical methods, convergence tests, or error estimates; a single sentence summarizing these aspects would improve readability for readers who consult only the abstract.
  2. Figure captions and axis labels should explicitly state the grid resolution and the number of echoes shown so that the visual evidence for approximate (but not exact) discrete self-similarity can be assessed without returning to the text.
  3. A short table or paragraph comparing the achieved fine-tuning parameter value and the number of observed echoes with the corresponding quantities from the authors' previous Teukolsky-wave runs would make the claimed improvement quantitative.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of our manuscript and for recommending minor revision. The report raises no specific major comments, so we have nothing to address point by point. We are pleased that the work is recognized as supplying an independent cross-check on approximate discrete self-similarity using a distinct initial-data family.

Circularity Check

0 steps flagged

No significant circularity

full rationale

This is a numerical relativity study that constructs initial data from linear wave solutions (Nakamura waves encoded in extrinsic curvature), solves constraints, evolves the system with standard methods, and reports observations of approximate discrete self-similarity in threshold solutions. The central claims rest on direct simulation output and comparison to prior literature rather than any derivation that reduces a result to its own fitted inputs or self-citations by construction. No equations or steps exhibit self-definitional, fitted-prediction, or uniqueness-imported circularity; the work is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on numerical evolution of the vacuum Einstein equations using tuned initial data families; no new physical constants or entities are introduced beyond standard GR assumptions.

free parameters (1)
  • Wave amplitude tuning parameter
    The strength of the initial gravitational wave is adjusted to approach the threshold of black hole formation.
axioms (2)
  • standard math Vacuum Einstein equations govern the spacetime evolution
    The simulations solve the vacuum Einstein field equations.
  • domain assumption Axial symmetry of the spacetime
    The setup assumes axisymmetry around one coordinate axis.

pith-pipeline@v0.9.1-grok · 5721 in / 1467 out tokens · 78348 ms · 2026-06-29T01:46:13.536278+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

52 extracted references · 19 linked inside Pith

  1. [1]

    thresh- old solutions

    for scalar fields. Given these observations, it became generally expected that the properties of critical collapse in spherical symmetry – in particular the existence of a unique self-similar critical solution with an associated unique critical exponent – should continue to hold even in the absence of spherical symmetry. As we will discuss in more detail ...

  2. [2]

    dressed up

    but are then “dressed up” to satisfy the non-linear constraint equations (see Section II A below). The lat- ter step typically involves inverting curved and non-linear Laplace operators, which can be complicated (see also the discussion in [29, 34]). Many authors therefore adopted Brill wave initial data [35] instead, which entails solving a flat and line...

  3. [3]

    off-centered

    for details). We chose the seed functionFto be a linear combination ofF + andF − so that the time deriva- tive of the combination vanishes att= 0, meaning that the initial slice represents a moment of time symmetry. 1 We use a hat to distinguish this choice of conformal rescaling, which is commonly adopted in the context of the initial value problem, from...

  4. [4]

    Here we focus on axisymmetric solutions and therefore construct new solutions form= 0

    Family A Explicit solutions for ˆANak ij are provided in [36, 37] for quadrupolar waves (ℓ= 2) withm= 2. Here we focus on axisymmetric solutions and therefore construct new solutions form= 0. Adopting seed functions similar to 5 those of [36, 37], P20 =A 4π 5 1/2 ∂ ∂t exp −(r−t) 2/2 (6a) Q20 =−A 4π 5 1/2 ∂ ∂t exp −(r+t) 2/2 ,(6b) whereAis the amplitude, w...

  5. [5]

    bean-like

    Family B We build a second family of quadrupolar waves from the functions (compare [36, 37]) P20 =A 4π 5 1/2 ∂ ∂t (r−t) 2 + 1 exp −(r−t) 2/2 (8a) Q20 =−A 4π 5 1/2 ∂ ∂t (r+t) 2 −1 exp −(r+t) 2/2 ,(8b) in which case the extrinsic curvature becomes ˆANak rr =Ae −r2/2(5−r 2)(1−3 cos 2 θ) (9a) ˆANak rθ =Ae −r2/2r(15−10r 2 +r 4) sinθcosθ(9b) ˆANak rφ = 0 (9c) ˆ...

  6. [6]

    M. W. Choptuik, Universality and scaling in gravita- tional collapse of a massless scalar field, Phys. Rev. Lett. 70, 9 (1993)

  7. [7]

    Garfinkle and G

    D. Garfinkle and G. C. Duncan, Scaling of curvature in subcritical gravitational collapse, Phys. Rev. D58, 064024 (1998), arXiv:gr-qc/9802061 [gr-qc]

  8. [8]

    Koike, T

    T. Koike, T. Hara, and S. Adachi, Critical Behavior in Gravitational Collapse of Radiation Fluid: A Renormal- ization Group (Linear Perturbation) Analysis, Phys. Rev. Lett.74, 5170 (1995), arXiv:gr-qc/9503007 [gr-qc]

  9. [9]

    Maison, Non-universality of critical behaviour in spherically symmetric gravitational collapse, Physics Let- ters B366, 82 (1996), arXiv:gr-qc/9504008 [gr-qc]

    D. Maison, Non-universality of critical behaviour in spherically symmetric gravitational collapse, Physics Let- ters B366, 82 (1996), arXiv:gr-qc/9504008 [gr-qc]

  10. [10]

    C. R. Evans and J. S. Coleman, Critical phenomena and self-similarity in the gravitational collapse of radi- ation fluid, Phys. Rev. Lett.72, 1782 (1994), arXiv:gr- qc/9402041 [gr-qc]

  11. [11]

    Hod and T

    S. Hod and T. Piran, Fine structure of Choptuik’s mass- scaling relation, Phys. Rev. D55, R440 (1997), arXiv:gr- qc/9606087 [gr-qc]

  12. [12]

    Gundlach, Understanding critical collapse of a scalar field, Phys

    C. Gundlach, Understanding critical collapse of a scalar field, Phys. Rev. D , 695 (1997), arXiv:gr-qc/9604019 [gr- qc]

  13. [13]

    Gundlach, D

    C. Gundlach, D. Hilditch, and J. M. Mart´ ın-Garc´ ıa, Crit- ical Phenomena in Gravitational Collapse, arXiv e-prints , arXiv:2507.07636 (2025), arXiv:2507.07636 [gr-qc]. 4 As in Sect. II A we adopt the notation and convention of [38] here, so that some expressions differ from those of [33] by constant factors

  14. [14]

    A. M. Abrahams and C. R. Evans, Critical behavior and scaling in vacuum axisymmetric gravitational collapse, Phys. Rev. Lett.70, 2980 (1993)

  15. [15]

    A. M. Abrahams and C. R. Evans, Universality in ax- isymmetric vacuum collapse, Phys. Rev. D49, 3998 (1994)

  16. [16]

    M. W. Choptuik, E. W. Hirschmann, S. L. Liebling, and F. Pretorius, Critical collapse of the massless scalar field in axisymmetry, Phys. Rev. D68, 044007 (2003), arXiv:gr-qc/0305003 [gr-qc]

  17. [17]

    T. W. Baumgarte, Aspherical deformations of the Choptuik spacetime, Phys. Rev. D98, 084012 (2018), arXiv:1807.10342 [gr-qc]

  18. [18]

    Marouda, D

    K. Marouda, D. Cors, H. R. R¨ uter, F. Atteneder, and D. Hilditch, Twist-free axisymmetric critical collapse of a complex scalar field, Phys. Rev. D109, 124042 (2024), arXiv:2402.06724 [gr-qc]

  19. [19]

    Gundlach, T

    C. Gundlach, T. W. Baumgarte, and D. Hilditch, Sim- ulations of gravitational collapse in null coordinates. II. Critical collapse of an axisymmetric scalar field, Phys. Rev. D110, 024019 (2024), arXiv:2404.15839 [gr-qc]

  20. [20]

    T. W. Baumgarte, C. Gundlach, and D. Hilditch, Criti- cal Phenomena in the Gravitational Collapse of Electro- magnetic Waves, Phys. Rev. Lett.123, 171103 (2019), arXiv:1909.00850 [gr-qc]

  21. [21]

    M. F. Perez Mendoza and T. W. Baumgarte, Critical phenomena in the gravitational collapse of electromag- netic dipole and quadrupole waves, Phys. Rev. D103, 124048 (2021), arXiv:2104.03980 [gr-qc]

  22. [22]

    G. D. Reid and M. W. Choptuik, Universality in the crit- ical collapse of the Einstein-Maxwell system, Phys. Rev. 15 D108, 104021 (2023), arXiv:2308.03943 [gr-qc]

  23. [23]

    Alcubierre, G

    M. Alcubierre, G. Allen, B. Br¨ ugmann, G. Lanfermann, E. Seidel, W.-M. Suen, and M. Tobias, Gravitational col- lapse of gravitational waves in 3D numerical relativity, Phys. Rev. D61, 041501 (2000), arXiv:gr-qc/9904013 [gr-qc]

  24. [24]

    Garfinkle and G

    D. Garfinkle and G. C. Duncan, Numerical evolution of Brill waves, Phys. Rev. D63, 044011 (2001), arXiv:gr- qc/0006073 [gr-qc]

  25. [25]

    Rinne, Constrained evolution in axisymmetry and the gravitational collapse of prolate Brill waves, Classical and Quantum Gravity25, 135009 (2008), arXiv:0802.3791 [gr-qc]

    O. Rinne, Constrained evolution in axisymmetry and the gravitational collapse of prolate Brill waves, Classical and Quantum Gravity25, 135009 (2008), arXiv:0802.3791 [gr-qc]

  26. [26]

    Sorkin, On critical collapse of gravitational waves, Classical and Quantum Gravity28, 025011 (2011), arXiv:1008.3319 [gr-qc]

    E. Sorkin, On critical collapse of gravitational waves, Classical and Quantum Gravity28, 025011 (2011), arXiv:1008.3319 [gr-qc]

  27. [27]

    Hilditch, T

    D. Hilditch, T. W. Baumgarte, A. Weyhausen, T. Di- etrich, B. Br¨ ugmann, P. J. Montero, and E. M¨ uller, Collapse of nonlinear gravitational waves in moving- puncture coordinates, Phys. Rev. D88, 103009 (2013), arXiv:1309.5008 [gr-qc]

  28. [28]

    C. Bona, J. Mass´ o, E. Seidel, and J. Stela, New Formal- ism for Numerical Relativity, Phys. Rev. Lett.75, 600 (1995), arXiv:gr-qc/9412071 [gr-qc]

  29. [29]

    Alcubierre, Appearance of coordinate shocks in hy- perbolic formalisms of general relativity, Phys

    M. Alcubierre, Appearance of coordinate shocks in hy- perbolic formalisms of general relativity, Phys. Rev. D 55, 5981 (1997), arXiv:gr-qc/9609015 [gr-qc]

  30. [30]

    Alcubierre and J

    M. Alcubierre and J. Mass´ o, Pathologies of hyperbolic gauges in general relativity and other field theories, Phys. Rev. D57, R4511 (1998), arXiv:gr-qc/9709024 [gr-qc]

  31. [31]

    Alcubierre, Are gauge shocks really shocks?, Clas- sical and Quantum Gravity22, 4071 (2005), arXiv:gr- qc/0503030 [gr-qc]

    M. Alcubierre, Are gauge shocks really shocks?, Clas- sical and Quantum Gravity22, 4071 (2005), arXiv:gr- qc/0503030 [gr-qc]

  32. [32]

    T. W. Baumgarte, C. Gundlach, and D. Hilditch, Critical phenomena in the collapse of quadrupolar and hexade- capolar gravitational waves, Phys. Rev. D107, 084012 (2023), arXiv:2303.05530 [gr-qc]

  33. [33]

    Khirnov and T

    A. Khirnov and T. Ledvinka, Slicing conditions for axisymmetric gravitational collapse of Brill waves, Classical and Quantum Gravity35, 215003 (2018), arXiv:1908.06034 [gr-qc]

  34. [34]

    Ledvinka and A

    T. Ledvinka and A. Khirnov, Universality of Curvature Invariants in Critical Vacuum Gravitational Collapse, Phys. Rev. Lett.127, 011104 (2021), arXiv:2102.09579 [gr-qc]

  35. [35]

    Su´ arez Fern´ andez, S

    I. Su´ arez Fern´ andez, S. Renkhoff, D. Cors Agull´ o, B. Br¨ ugmann, and D. Hilditch, Evolution of Brill waves with an adaptive pseudospectral method, Phys. Rev. D 106, 024036 (2022)

  36. [36]

    Alcubierre, Hyperbolic slicings of spacetime: singular- ity avoidance and gauge shocks, Classical and Quantum Gravity20, 607 (2003), arXiv:gr-qc/0210050 [gr-qc]

    M. Alcubierre, Hyperbolic slicings of spacetime: singular- ity avoidance and gauge shocks, Classical and Quantum Gravity20, 607 (2003), arXiv:gr-qc/0210050 [gr-qc]

  37. [37]

    T. W. Baumgarte, B. Br¨ ugmann, D. Cors, C. Gund- lach, D. Hilditch, A. Khirnov, T. Ledvinka, S. Renkhoff, and I. S. Fern´ andez, Critical Phenomena in the Collapse of Gravitational Waves, Phys. Rev. Lett.131, 181401 (2023), arXiv:2305.17171 [gr-qc]

  38. [38]

    S. A. Teukolsky, Linearized quadrupole waves in general relativity and the motion of test particles, Phys. Rev. D 26, 745 (1982)

  39. [39]

    Rostworowski, Vacuum axisymmetric gravitational collapse revisited: Preliminary investigation, Phys

    A. Rostworowski, Vacuum axisymmetric gravitational collapse revisited: Preliminary investigation, Phys. Rev. D111, 104067 (2025), arXiv:2503.17549 [gr-qc]

  40. [40]

    D. R. Brill, On the positive definite mass of the Bondi- Weber-Wheeler time-symmetric gravitational waves, An- nals of Physics7, 466 (1959)

  41. [41]

    T. Nakamura, General Solutions to the Linearized Ein- stein Equations and Initial Data for Three Dimensional Time Evolution of Pure Gravitational Waves, Progress of Theoretical Physics72, 746 (1984)

  42. [42]

    Shibata and T

    M. Shibata and T. Nakamura, Evolution of three- dimensional gravitational waves: Harmonic slicing case, Phys. Rev. D52, 5428 (1995)

  43. [43]

    Rinne, Explicit solution of the linearized Einstein equations in TT gauge for all multipoles, Classical and Quantum Gravity29, 048003 (2009)

    O. Rinne, Explicit solution of the linearized Einstein equations in TT gauge for all multipoles, Classical and Quantum Gravity29, 048003 (2009)

  44. [44]

    T. W. Baumgarte and S. L. Shapiro, Numerical inte- gration of Einstein’s field equations, Phys. Rev. D59, 024007 (1998), arXiv:gr-qc/9810065 [gr-qc]

  45. [45]

    T. W. Baumgarte, P. J. Montero, I. Cordero-Carri´ on, and E. M¨ uller, Numerical relativity in spherical polar co- ordinates: Evolution calculations with the BSSN formu- lation, Phys. Rev. D87, 044026 (2013), arXiv:1211.6632 [gr-qc]

  46. [46]

    Nakamura, K

    T. Nakamura, K. Oohara, and Y. Kojima, General Rela- tivistic Collapse to Black Holes and Gravitational Waves from Black Holes, Progress of Theoretical Physics Sup- plement90, 1 (1987)

  47. [47]

    J. D. Brown, Covariant formulations of Baumgarte, Shapiro, Shibata, and Nakamura and the standard gauge, Phys. Rev. D79, 104029 (2009), arXiv:0902.3652 [gr-qc]

  48. [48]

    Gourgoulhon,3+1 Formalism in General Relativity, Vol

    E. Gourgoulhon,3+1 Formalism in General Relativity, Vol. 846 (Springer, 2012)

  49. [49]

    Ruchlin, Z

    I. Ruchlin, Z. B. Etienne, and T. W. Baumgarte, SENR /NRPy+: Numerical relativity in singular curvilinear coordinate systems, Phys. Rev. D97, 064036 (2018), arXiv:1712.07658 [gr-qc]

  50. [50]

    T. W. Baumgarte and D. Hilditch, Shock-avoiding slicing conditions: Tests and calibrations, Phys. Rev. D106, 044014 (2022), arXiv:2207.06376 [gr-qc]

  51. [51]

    R. Beig, P. T. Chru´ sciel, and W. Cong, Who’s afraid of a negative lapse?, arXiv e-prints , arXiv:2505.03516 (2025), arXiv:2505.03516 [gr-qc]

  52. [52]

    Stephani, D

    H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers, and E. Herlt,Exact solutions of Einstein’s field equations (Cambridge University Press, 2003)