Pith. sign in

REVIEW 5 major objections 4 minor 52 references

Stable long-term evolution in numerical relativity

T0 review · 5 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Undamped momentum-constraint evolution suppresses a late-time instability in black hole simulations, reaching $10^5 M$.

desk verdict A useful, clearly reported numerical recipe for long-term black hole runs, but the mechanism is inferred rather than measured and the evidence is entirely 1D. read the letter →

arxiv 2501.01055 v2 pith:6QB3VPK7 submitted 2025-01-02 gr-qc hep-th

Sebastian Garcia-Saenz , Guangzhou Guo , Peng Wang , Xinmiao Wang This is my paper · ORCID
classification gr-qchep-th MSC 83C0583C5765M0683-08
keywords numericalrelativityBSSNformulationCCZ4constraintviolationsmomentumlate-timeinstabilityblackholeevolutionlong-termstability
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper reports that the standard BSSN formulation of numerical relativity develops a late-time numerical instability even for the simplest black hole, Schwarzschild, and traces the culprit to accumulated violations of the momentum constraint. To fix it, the authors modify the conformal covariant Z4 scheme so that momentum constraint violations are evolved as dynamical variables but are not damped. Their two proposed schemes, CCZ4' and CCZ3, keep black hole simulations stable out to times of order $10^5 M$ for Schwarzschild and for charged black holes undergoing spontaneous scalarization. If correct, the result gives a practical recipe for studying weak, slow instabilities like superradiance that current codes cannot follow long enough.

What carries the argument

The central object is the spatial momentum-constraint-violation vector $Z^i$, the projection of the Z4 four-vector onto the spatial hypersurface. The key move is to evolve $Z^i$ without damping: CCZ4' splits the single damping parameter $\kappa_1$ into $\kappa_\Theta$ (damping $\Theta$) and $\kappa_\Gamma$ (damping $Z^i$) and sets $\kappa_\Gamma = 0$; CCZ3 removes the evolution of $\Theta$ entirely (sets $\Theta = 0$) while keeping the $Z^i$-carrying equation for $\tilde{\Lambda}^i$. This separation lets the Hamiltonian constraint violation be damped strongly without inducing the nonlinear instabilities that damping of momentum violations causes.

What would settle it

A three-dimensional evolution of a Kerr black hole using CCZ3 with $\kappa_\Gamma = 0$ that develops the same late-time blow-up in the apparent horizon area would falsify the general claim; alternatively, a spherical BSSN run with a different outer boundary condition that removes the instability would show the effect is a boundary artifact, not a formulation failure.

Watch

Extended reading notes

Core claim

The central discovery is that the late-time instability seen in BSSN evolutions of black holes is driven by violations of the momentum constraint $Z^i$, not by the Hamiltonian constraint violation $\Theta$, and that the cure is to let $Z^i$ propagate freely without damping. The paper introduces two schemes, CCZ4' and CCZ3, in which the momentum constraint violation is evolved through the conformal connection variable $\tilde{\Lambda}^i$ without a damping term. In these schemes the damping of the Hamiltonian constraint violation can be made strong, which keeps the Hamiltonian constraint under control, while the absence of momentum damping avoids the nonlinear instabilities that strong damping otherwise triggers. The schemes are demonstrated in spherical symmetry for a Schwarzschild black hole, a Reissner-Nordström black hole, and black hole spontaneous scalarization in the Einstein-Maxwell-scalar model, with stable evolutions reaching times of order $10^5 M$.

Load-bearing premise

The numerical evidence is one-dimensional: every simulation is spherically symmetric, so the claim that these schemes resolve the instability in general black hole spacetimes rests on the assumption that the same behavior holds in full three-dimensional settings such as rotating or binary black holes.

Editorial extensions

If this is right

  • Long-term black hole simulations with the CCZ3 and CCZ4' schemes remain stable to at least $10^5 M$, a regime where BSSN and unmodified CCZ4 break down.
  • The momentum constraint violation, not the Hamiltonian one, is the driver of the late-time instability; disabling momentum propagation (CCZ0) reproduces the BSSN failure.
  • Strong damping of the Hamiltonian constraint violation is beneficial once momentum damping is removed, so constraint control can be improved without destabilizing the simulation.
  • In matter spacetimes, the CCZ3 scheme with non-propagating electromagnetic constraints gives the most accurate and robust evolution of spontaneous scalarization.

Reading between the lines

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

  • If the mechanism is generic, the 'no damping of $Z^i$' rule may also improve other free-evolution formulations, including generalized harmonic evolutions, wherever momentum constraint violations accumulate.
  • The same principle could be tested directly in full three-dimensional evolutions of Kerr black holes with superradiant scalar clouds; the predicted requirement is that simulations remain stable to times of order $10^5$ to $10^6 M$ without momentum-constraint damping.
  • The result suggests that the choice of which constraint to damp is more consequential than the overall strength of damping, which may inform future constraint-damping designs.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

5 major / 4 minor

Summary. The paper investigates a late-time numerical instability observed in long-term black hole evolutions with the BSSN formulation, even for a Schwarzschild spacetime. The authors attribute this instability to accumulated violations of the momentum constraint and propose two modified versions of the conformal covariant Z4 system, CCZ4' and CCZ3, distinguished by propagating momentum-constraint violations Z^i without damping while allowing strong damping of the Hamiltonian-constraint violation Θ. Using the BlackHoles@Home platform in spherical symmetry, they show that CCZ4' and CCZ3 maintain apparent-horizon areas close to the expected Schwarzschild value for times up to 10^5 M, whereas BSSN and a comparison scheme CCZ0 (with Z^i set to zero) develop late-time instabilities. They extend the tests to Reissner-Nordström black holes and to spontaneous scalarization in the Einstein-Maxwell-scalar model, reporting that CCZ3 is the most robust, while CCZ4' fails to converge accurately in the scalarization case. The paper's central claim is that undamped propagation of momentum-constraint violations is the key to removing the late-time instability.

Significance. If the central claim holds, the paper offers a simple and practically useful modification for numerical relativity: turning off momentum-constraint damping in CCZ4-type systems, which could enable long-term simulations of weakly unstable black hole systems (e.g., superradiance, light-ring instabilities). The strength of the paper lies in the breadth of parameter variation (outer boundary, resolution, dissipation, CFL) and the use of a known benchmark, the Schwarzschild horizon area, so the stability result is not fitted to a target. However, the significance is currently tempered by three gaps: (i) the causal mechanism is inferred from scheme-by-scheme comparisons without direct measurement of constraint violations; (ii) the evidence is entirely one-dimensional, while the abstract claims resolution of the instability in 'black hole spacetimes with matter fields' in general; and (iii) no convergence study is reported, leaving open the possibility that the observed stability is partly a numerical artifact. The paper also contains a printed inconsistency in the evolution equation for Θ that affects the reproducibility of the CCZ4' scheme.

major comments (5)
  1. [Section III.A, Figs. 1-3] The paper's central causal claim—that the late-time instability stems from accumulated violations of the momentum constraint—is not directly evidenced. No figure or table reports the norm of Z^i, Θ, or the residuals of the Hamiltonian and momentum constraints for any scheme. The inference relies on comparing the stability of BSSN/CCZ0 (no Z^i propagation) with CCZ4'/CCZ3 (Z^i propagation without damping). Without direct diagnostic data, alternative explanations such as gauge drift, outer-boundary noise, or resolution/dissipation effects cannot be excluded. Please plot, for representative runs, the L2 norms of Z^i and Θ as functions of time for BSSN, CCZ4', CCZ3, and CCZ0, and show that the growth of Z^i correlates with the onset of the late-time instability.
  2. [Section II, Eq. (17)] The evolution equation for Θ in the conformal CCZ4 system, Eq. (17), is inconsistent with the non-conformal version, Eq. (6). Equation (17) lacks the matter source term (-16παρ), the damping term (-ακ1(2+κ2)Θ), and the -Z^i∂_iα term present in Eq. (6). This is not a minor typo because CCZ4' is defined by replacing κ1 in Eq. (15) with κΘ and in Eq. (18) with κΓ, leaving the damping of Eq. (17) unspecified. As printed, κΘ does not damp Θ in its own evolution equation, yet the paper attributes the stability of CCZ4' to Hamiltonian constraint damping. The authors must clarify the actual damping structure used in their CCZ4' implementation and correct Eq. (17) so that the scheme is reproducible and the interpretation is sound.
  3. [Section III, Fig. 3] No convergence study is reported. The resolution tests in Fig. 3 use NR=200, 300, and 400, but the plots do not show whether the error in the apparent horizon area decreases with resolution, and no convergence order is quoted. Given that the central evidence is numerical stability over 10^5 M, a convergence analysis is necessary to ensure the results are not dominated by the large Kreiss-Oliger dissipation (ϵKO=0.2) or other numerical artifacts. Please provide a convergence test, e.g., Richardson extrapolation of Ah at selected times for at least three resolutions.
  4. [Section IV.C, Fig. 6 and Abstract] The abstract states that CCZ4' and CCZ3 'effectively resolve the late-time numerical instability not only in Schwarzschild spacetimes but also in black hole spacetimes with matter fields.' However, Fig. 6 (upper-left panel) and the text state that the CCZ4' scheme does not converge in the spontaneous scalarization case, at least as far as the apparent horizon is concerned. Thus the abstract overstates the success of CCZ4' in matter spacetimes. Please revise the abstract and conclusions to distinguish the performance of CCZ3 from that of CCZ4', or restrict the general claim to CCZ3, which is the only scheme that converges in the matter-field tests.
  5. [Section III and Abstract] The numerical evidence is restricted to spherical symmetry: the paper states in Section III that 'we restrict our attention to spherically symmetric systems in this paper.' The abstract's claim that the schemes resolve the late-time instability 'not only in Schwarzschild spacetimes but also in black hole spacetimes with matter fields' goes beyond the tested domain. In spherical symmetry, Zi has only a radial component and the vector structure of the Einstein equations is degenerate. The paper should explicitly qualify the abstract and conclusions as applying to spherically symmetric spacetimes, or, if the general claim is intended, provide at least one non-spherical test (e.g., a Kerr or binary black hole run) to support it.
minor comments (4)
  1. [Section III, Eq. (27) vs Fig. 1] The benchmark parameter R0 is given as R0=0.00012 in Eq. (27) but as R0=0.0012 in the caption of Fig. 1 for the rmax=60000M case; please correct this typo.
  2. [Section III] The text says 'we use fourth-order finite differential on the spatial direction'; this should read 'fourth-order finite differences.'
  3. [Section II, Eq. (19)] The paper never explains how Zi is reconstructed from the evolved variable Λ̃i in the CCZ4' and CCZ3 implementations. Since Eq. (19) defines Λ̃i ≡ Λ̄i + 2γ̄ij Zj, please clarify the reconstruction step used in the code.
  4. [Figs. 4 and 5] In the RN and scalarization figures, the BSSN and CCZ0 panels are plotted only up to t=1000M, whereas the CCZ4' and CCZ3 panels extend to 10^5M; using the same time range in all panels would make the comparison more straightforward.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the schemes are validated against externally fixed benchmarks, and no fitted quantity is renamed as a prediction.

full rationale

The paper's central claim is an empirical numerical result: BSSN evolutions of Schwarzschild develop a late-time instability, while CCZ4' and CCZ3, defined by propagating Z^i without damping, remain stable up to about 10^5 M. The target observable, the apparent horizon area A_h, is an external benchmark with known expected value 4πM^2 for Schwarzschild; no parameter in the proposed schemes is fitted to that value, and the damping parameters are varied and reported openly (Figs. 2-3). The scheme definitions in Section II are explicit algebraic choices about which constraint variables are evolved and damped, and the conclusion that undamped Z^i propagation is key is inferred from pairwise comparisons (CCZ4' with κΓ=0 versus CCZ4; CCZ0 versus CCZ3), not from the definitions themselves. Citations to the authors' prior work (Refs. [7] and [44]) are background for the EMS scalarization model and are not load-bearing for the numerical stability claim, which is supported by code runs against known horizon areas. The paper's limitations—spherical symmetry, no direct constraint-violation norm reported, and the discrepancy between Eq. (17) and Eq. (6) in the Θ damping term—affect evidentiary strength but do not make any prediction reduce to its inputs by construction. The derivation and numerical evidence are self-contained against external benchmarks, so the circularity score is 0.

Assumptions & free parameters 6 free parameters · 5 assumptions · 0 invented entities

The central stability results depend on the adopted Z4 and CCZ4 formulation, the moving puncture gauge, the spherical reduction, and the choice to switch off momentum damping. These are pulled from prior literature or chosen by hand; none is derived in the paper. The numerical parameter choices are shown to be non-critical in most cases, but the spherical-to-3D extrapolation and the lack of direct constraint diagnostics remain unverified assumptions.

free parameters (6)
  • kappa_Gamma (CCZ4' momentum damping) = 0
    Set to zero in CCZ4' to achieve stable long-term evolution; the paper argues damping of Zi causes instabilities. This is a hand-chosen value, not fitted to data, but central to the proposed scheme.
  • kappa_Theta (CCZ4' Hamiltonian damping) = 1.0
    Chosen to show that strong Hamiltonian damping is compatible with stability once momentum damping is off; stability is claimed for sufficiently large kappa_Theta.
  • kappa_1 (CCZ3 damping) = 0
    CCZ3 is stable only with kappa_1 = 0 in the tested setup; nonzero kappa_1 reintroduces late-time instabilities (Fig. 2 lower-right).
  • Gauge damping eta = 1
    Standard Gamma-driver damping used in the main text; variations change BSSN stability but not CCZ3, so it is not load-bearing.
  • Benchmark grid parameters (rmax, R0, a, NR, eps_KO, CFL) = 60000M, 0.00012, 0.07, 300, 0.2, 1.0
    Default numerical resolution and dissipation settings; the qualitative instability and cure persist across variations shown in Figs. 1 and 3.
  • Scalarization perturbation amplitude p and coupling alpha_0 = 1e-4 and 1
    Chosen so the scalarization grows slowly and demands long-term evolution; not fitted to the stability claim.
assumptions (5)
  • domain assumption The damped Z4 extension of Einstein's equations (Eq. 1) with Z_mu = 0 reducing to GR is a valid starting point.
    Adopted from Refs. [25,34]; the paper does not derive it.
  • domain assumption The moving puncture gauge, 1+log lapse and Gamma-driver (Eqs. 20-21), is suitable for long-term black hole evolutions.
    Standard gauge from Refs. [35,36]; the stability conclusions depend on it.
  • domain assumption Fourth-order finite differencing with Kreiss-Oliger dissipation of strength 0.2 is sufficiently accurate and does not mask the instability.
    Parameter variations are shown, but no formal convergence test or error estimate is reported.
  • domain assumption Spherical symmetry is representative of the instability mechanism in full 3D numerical relativity.
    The paper restricts to spherically symmetric systems and does not test rotating or binary black holes.
  • domain assumption The apparent horizon area is a reliable diagnostic for numerical instability.
    All stability statements are based on Ah time series; no constraint norms are shown.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Stable long-term evolution in numerical relativity." pith.science (2026). https://pith.science/paper/6QB3VPK7

@misc{pith2026250101055,
  author       = {Pith},
  title        = {Pith review of: Stable long-term evolution in numerical relativity},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6QB3VPK7}},
  note         = {Machine review of arXiv:2501.01055}
}
read the original abstract

We report on the potential occurrence of a numerical instability in the long-time simulation of black holes using the Baumgarte-Shapiro-Shibata-Nakamura formulation of numerical relativity, even in the simple set-up of a Schwarzschild black hole. Through extensive numerical experiments, we identify that this "late-time instability" arises from accumulated violations of the momentum constraint. To address this issue, we propose two modified versions of the so-called conformal covariant Z4 scheme, designed to propagate momentum constraint violations without damping. Our results demonstrate that these alternative formulations, which we refer to as CCZ4' and CCZ3, effectively resolve the late-time numerical instability not only in Schwarzschild spacetimes but also in black hole spacetimes with matter fields. Notably, by preventing damping of the momentum constraint violation, the Hamiltonian constraint damping can be significantly increased, which plays a crucial role in stabilizing long-term evolution in our proposed schemes.

Figures

Figures reproduced from arXiv: 2501.01055 by the authors.

Figure 1
Figure 1. FIG. 1. Schwarzschild black hole evolution, as measured by the area of the apparent horizon, [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Numerical evolution of a Schwarzschild black hole, as measured by the area of the apparent hori [PITH_FULL_IMAGE:figures/full_fig_p011_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Numerical evolution of a Schwarzschild black hole, as measured by the area of the apparent horizon, [PITH_FULL_IMAGE:figures/full_fig_p012_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Evolution of a RN black hole, as measured by the area of the apparent horizon, [PITH_FULL_IMAGE:figures/full_fig_p015_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. The area of the apparent horizon, [PITH_FULL_IMAGE:figures/full_fig_p017_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. The evolution of the subtracted apparent horizon area ∆ [PITH_FULL_IMAGE:figures/full_fig_p018_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Evolution of Schwarzschild ( [PITH_FULL_IMAGE:figures/full_fig_p021_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

52 extracted references · 12 canonical work pages

  1. [1]

    Abbott et al

    B.P. Abbott et al. Observation of Gravitational Waves from a Binary Black Hole Merger. Phys. Rev. Lett., 116(6):061102, 2016. arXiv:1602.03837, doi:10.1103/PhysRevLett.116.061102

    arXiv 2016

  2. [2]

    First M87 Event Horizon Telescope Results

    Kazunori Akiyama et al. First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole. Astrophys. J. Lett. , 875:L1, 2019. arXiv:1906.11238, doi:10.3847/2041-8213/ab0ec7. 21 BSSN (p = 0.75) BSSN (p = 0.5) BSSN (p = 0.2) 0 5000 10 000 15 000 20 000 25 000 3.9980 3.9985 3.9990 3.9995 4.0000 4.0005 4.0010 4.0015 4.0020 t/M Ah/4π M2 ...

    arXiv 2019

  3. [3]

    First Sagittarius A* Event Horizon Telescope Results

    Kazunori Akiyama et al. First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way. Astrophys. J. Lett. , 930(2):L12, 2022. arXiv:2311.08680, doi:10.3847/2041-8213/ac6674

    arXiv 2022

  4. [4]

    late-time instability

    Similarly to the vacuum and electro-vacuum cases, our results show that the BSSN and CCZ0 formulations do not provide robust results in the long term, failing in this case to fully capture the formation of the scalarized black hole, with the evolution diverging in a comparatively short time. In contrast, both the CCZ4’ and CCZ3 formulations successfully e...

    work page 2000

  5. [5]

    Manuela Campanelli, C. O. Lousto, P. Marronetti, and Y. Zlochower. Accurate evolutions of orbiting black-hole binaries without excision. Phys. Rev. Lett. , 96:111101, 2006. arXiv:gr-qc/0511048, doi: 10.1103/PhysRevLett.96.111101

    arXiv 2006

  6. [6]

    Evolution of binary black hole spacetimes

    Frans Pretorius. Evolution of binary black hole spacetimes. Phys. Rev. Lett., 95:121101, 2005. arXiv: gr-qc/0507014, doi:10.1103/PhysRevLett.95.121101

    arXiv 2005

  7. [7]

    Baker, Joan Centrella, Dae-Il Choi, Michael Koppitz, and James van Meter

    John G. Baker, Joan Centrella, Dae-Il Choi, Michael Koppitz, and James van Meter. Gravitational wave extraction from an inspiraling configuration of merging black holes. Phys. Rev. Lett. , 96:111102,

    work page

  8. [8]

    Damour, N

    T. Damour, N. Deruelle, and R. Ruffini. On Quantum Resonances in Stationary Geometries. Lett. Nuovo Cim. , 15:257–262, 1976. doi:10.1007/BF02725534

    work page doi:10.1007/bf02725534 1976

Show all 52 references
  1. [9]

    Nonlinear Stability of Black Holes with a Stable Light Ring

    Guangzhou Guo, Peng Wang, and Yupeng Zhang. Nonlinear Stability of Black Holes with a Stable Light Ring. 3 2024. arXiv:2403.02089. 22

    2024 arXiv

  2. [10]

    Detweiler

    Steven L. Detweiler. KLEIN-GORDON EQUATION AND ROTATING BLACK HOLES. Phys. Rev. D, 22:2323–2326, 1980. doi:10.1103/PhysRevD.22.2323

    1980 doi

  3. [11]

    String Axiverse

    Asimina Arvanitaki, Savas Dimopoulos, Sergei Dubovsky, Nemanja Kaloper, and John March-Russell. String Axiverse. Phys. Rev. D , 81:123530, 2010. arXiv:0905.4720, doi:10.1103/PhysRevD.81. 123530

    2010 arXiv

  4. [12]

    Sam R. Dolan. Instability of the massive Klein-Gordon field on the Kerr spacetime. Phys. Rev. D , 76:084001, 2007. arXiv:0705.2880, doi:10.1103/PhysRevD.76.084001

    2007 arXiv

  5. [13]

    Vitor Cardoso, Oscar J. C. Dias, Jose P. S. Lemos, and Shijun Yoshida. The Black hole bomb and superradiant instabilities. Phys. Rev. D , 70:044039, 2004. [Erratum: Phys.Rev.D 70, 049903 (2004)]. arXiv:hep-th/0404096, doi:10.1103/PhysRevD.70.049903

    2004 arXiv

  6. [14]

    Montero, Jos´ e A

    Nicolas Sanchis-Gual, Juan Carlos Degollado, Pedro J. Montero, Jos´ e A. Font, and Carlos Herdeiro. Explosion and Final State of an Unstable Reissner-Nordstr¨ om Black Hole. Phys. Rev. Lett. , 116(14):141101, 2016. arXiv:1512.05358, doi:10.1103/PhysRevLett.116.141101

    2016 arXiv

  7. [15]

    Sam R. Dolan. Superradiant instabilities of rotating black holes in the time domain. Phys. Rev. D , 87(12):124026, 2013. arXiv:1212.1477, doi:10.1103/PhysRevD.87.124026

    2013 arXiv

  8. [16]

    Horowitz

    Gary T. Horowitz. Introduction to Holographic Superconductors. Lect. Notes Phys., 828:313–347, 2011. arXiv:1002.1722, doi:10.1007/978-3-642-04864-7_10

    2011 arXiv

  9. [17]

    East and Frans Pretorius

    William E. East and Frans Pretorius. Superradiant Instability and Backreaction of Massive Vector Fields around Kerr Black Holes. Phys. Rev. Lett. , 119(4):041101, 2017. arXiv:1704.04791, doi: 10.1103/PhysRevLett.119.041101

    2017 arXiv

  10. [18]

    Pedro V. P. Cunha, Carlos Herdeiro, Eugen Radu, and Nicolas Sanchis-Gual. Exotic Compact Objects and the Fate of the Light-Ring Instability. Phys. Rev. Lett., 130(6):061401, 2023. arXiv:2207.13713, doi:10.1103/PhysRevLett.130.061401

    2023 arXiv

  11. [19]

    Vitor Cardoso, Lu ´ ıs C. B. Crispino, Caio F. B. Macedo, Hirotada Okawa, and Paolo Pani. Light rings as observational evidence for event horizons: long-lived modes, ergoregions and nonlinear instabilities of ultracompact objects. Phys. Rev. D , 90(4):044069, 2014. arXiv:1406....

    2014 arXiv

  12. [20]

    Baumgarte and Stuart L

    Thomas W. Baumgarte and Stuart L. Shapiro. On the numerical integration of Einstein’s field equa- tions. Phys. Rev. D , 59:024007, 1998. arXiv:gr-qc/9810065, doi:10.1103/PhysRevD.59.024007

    1998 arXiv

  13. [21]

    Evolution of three-dimensional gravitational waves: Harmonic slicing case

    Masaru Shibata and Takashi Nakamura. Evolution of three-dimensional gravitational waves: Harmonic slicing case. Phys. Rev. D , 52:5428–5444, 1995. doi:10.1103/PhysRevD.52.5428

    1995 doi

  14. [22]

    Etienne, and Thomas W

    Ian Ruchlin, Zachariah B. Etienne, and Thomas W. Baumgarte. SENR/NRPy+: Numerical Relativity in Singular Curvilinear Coordinate Systems. Phys. Rev. D , 97(6):064036, 2018. arXiv:1712.07658, 23 doi:10.1103/PhysRevD.97.064036

    2018 arXiv

  15. [23]

    David Brown

    J. David Brown. Covariant formulations of BSSN and the standard gauge. Phys. Rev. D , 79:104029,

  16. [24]

    On the stability of covariant BSSN formulation

    Ryosuke Urakawa, Takuya Tsuchiya, and Gen Yoneda. On the stability of covariant BSSN formulation. Class. Quant. Grav. , 39(16):165002, 2022. arXiv:2206.13944, doi:10.1088/1361-6382/ac7e16

    2022 arXiv

  17. [25]

    Conformal and covariant formulation of the Z4 system with constraint-violation damping.Phys

    Daniela Alic, Carles Bona-Casas, Carles Bona, Luciano Rezzolla, and Carlos Palenzuela. Conformal and covariant formulation of the Z4 system with constraint-violation damping.Phys. Rev. D, 85:064040,

  18. [26]

    Baumgarte and Stuart L

    Thomas W. Baumgarte and Stuart L. Shapiro. Relativistic radiation hydrodynamics in a reference- metric formulation. Phys. Rev. D, 102(10):104001, 2020. arXiv:2009.08990, doi:10.1103/PhysRevD. 102.104001

    2020 arXiv

  19. [27]

    Compact binary evolutions with the Z4c formulation

    David Hilditch, Sebastiano Bernuzzi, Marcus Thierfelder, Zhoujian Cao, Wolfgang Tichy, and Bernd Bruegmann. Compact binary evolutions with the Z4c formulation. Phys. Rev. D , 88:084057, 2013. arXiv:1212.2901, doi:10.1103/PhysRevD.88.084057

    2013 arXiv

  20. [28]

    Constraint violation in free evolution schemes: Comparing BSSNOK with a conformal decomposition of Z4

    Sebastiano Bernuzzi and David Hilditch. Constraint violation in free evolution schemes: Comparing BSSNOK with a conformal decomposition of Z4. Phys. Rev. D , 81:084003, 2010. arXiv:0912.2920, doi:10.1103/PhysRevD.81.084003

    2010 arXiv

  21. [29]

    Montero, Jose A

    Nicolas Sanchis-Gual, Pedro J. Montero, Jose A. Font, Ewald M¨ uller, and Thomas W. Baumgarte. Fully covariant and conformal formulation of the Z4 system in a reference-metric approach: comparison with the BSSN formulation in spherical symmetry. Phys. Rev. D , 89(10):104033, 2...

    2014 arXiv

  22. [30]

    Constraint preserving boundary conditions for the Z4c formulation of general relativity

    Milton Ruiz, David Hilditch, and Sebastiano Bernuzzi. Constraint preserving boundary conditions for the Z4c formulation of general relativity. Phys. Rev. D , 83:024025, 2011. arXiv:1010.0523, doi: 10.1103/PhysRevD.83.024025

    2011 arXiv

  23. [31]

    Reid and Matthew W

    Gray D. Reid and Matthew W. Choptuik. Reference metric approach to the Z4 system. Phys. Rev. D , 108(12):124070, 2023. arXiv:2309.05094, doi:10.1103/PhysRevD.108.124070

    2023 arXiv

  24. [32]

    Numerical stability of the Z4c formulation of general relativity

    Zhoujian Cao and David Hilditch. Numerical stability of the Z4c formulation of general relativity. Phys. Rev. D , 85:124032, 2012. arXiv:1111.2177, doi:10.1103/PhysRevD.85.124032

    2012 arXiv

  25. [33]

    Apples with Apples comparison of 3+1 conformal numerical relativity schemes

    David Daverio, Yves Dirian, and Ermis Mitsou. Apples with Apples comparison of 3+1 conformal numerical relativity schemes. 10 2018. arXiv:1810.12346

    2018 arXiv

  26. [34]

    Baumgarte, Zachariah B

    Vassilios Mewes, Yosef Zlochower, Manuela Campanelli, Thomas W. Baumgarte, Zachariah B. Etienne, Federico G. Lopez Armengol, and Federico Cipolletta. Numerical relativity in spherical coordinates: A new dynamical spacetime and general relativistic MHD evolution framework for t...

    2020 arXiv

  27. [35]

    A New formalism for numerical relativity

    Carles Bona, Joan Masso, Edward Seidel, and Joan Stela. A New formalism for numerical relativity. Phys. Rev. Lett. , 75:600–603, 1995. arXiv:gr-qc/9412071, doi:10.1103/PhysRevLett.75.600. 24

    1995 arXiv

  28. [36]

    Gauge conditions for long term numerical black hole evolutions without excision

    Miguel Alcubierre, Bernd Bruegmann, Peter Diener, Michael Koppitz, Denis Pollney, Edward Seidel, and Ryoji Takahashi. Gauge conditions for long term numerical black hole evolutions without excision. Phys. Rev. D , 67:084023, 2003. arXiv:gr-qc/0206072, doi:10.1103/PhysRevD.67.084023

    2003 arXiv

  29. [37]

    Time Step Size Limitation Introduced by the BSSN Gamma Driver

    Erik Schnetter. Time Step Size Limitation Introduced by the BSSN Gamma Driver. Class. Quant. Grav., 27:167001, 2010. arXiv:1003.0859, doi:10.1088/0264-9381/27/16/167001

    2010 arXiv

  30. [38]

    Constraint damping of the conformal and co- variant formulation of the Z4 system in simulations of binary neutron stars.Phys

    Daniela Alic, Wolfgang Kastaun, and Luciano Rezzolla. Constraint damping of the conformal and co- variant formulation of the Z4 system in simulations of binary neutron stars.Phys. Rev. D, 88(6):064049,

  31. [39]

    Press, Saul A

    William H. Press, Saul A. Teukolsky, William T. Vetterling, and B. P. Flannery. Numerical Recipes: The Art of Scientific Computing (Third Edition) . Cambridge University Press, 2007

    2007

  32. [40]

    David Brown

    J. David Brown. Probing the puncture for black hole simulations. Phys. Rev. D , 80:084042, 2009. arXiv:0908.3814, doi:10.1103/PhysRevD.80.084042

    2009 arXiv

  33. [41]

    Zachariah B. Etienne. Improved Moving-Puncture Techniques for Compact Binary Simulations. 4 2024. arXiv:2404.01137

    2024 arXiv

  34. [42]

    Martin-Garcia, Gioel Calabrese, and Ian Hinder

    Carsten Gundlach, Jose M. Martin-Garcia, Gioel Calabrese, and Ian Hinder. Constraint damping in the Z4 formulation and harmonic gauge. Class. Quant. Grav., 22:3767–3774, 2005. arXiv:gr-qc/0504114, doi:10.1088/0264-9381/22/17/025

    2005 arXiv

  35. [43]

    Etienne and Ian Ruchlin et al

    Zachariah B. Etienne and Ian Ruchlin et al. BlackHoles@Home, 2022, To find out more, visit https://blackholesathome.net/

    2022

  36. [44]

    Black hole accretion of scalar clouds with spontaneous symmetry breaking

    Sebastian Garcia-Saenz, Guangzhou Guo, Peng Wang, and Xinmiao Wang. Black hole accretion of scalar clouds with spontaneous symmetry breaking. Phys. Rev. D , 110(12):124045, 2024. arXiv: 2409.13184, doi:10.1103/PhysRevD.110.124045

    2024 arXiv

  37. [45]

    Hirschmann, Luis Lehner, Steven L

    Eric W. Hirschmann, Luis Lehner, Steven L. Liebling, and Carlos Palenzuela. Black Hole Dynamics in Einstein-Maxwell-Dilaton Theory. Phys. Rev. D , 97(6):064032, 2018. arXiv:1706.09875, doi: 10.1103/PhysRevD.97.064032

    2018 arXiv

  38. [46]

    The Einstein-Maxwell system in 3+1 form and initial data for multiple charged black holes.Phys

    Miguel Alcubierre, Juan Carlos Degollado, and Marcelo Salgado. The Einstein-Maxwell system in 3+1 form and initial data for multiple charged black holes.Phys. Rev. D, 80:104022, 2009. arXiv:0907.1151, doi:10.1103/PhysRevD.80.104022

    2009 arXiv

  39. [47]

    Bubbles kick off primordial black holes to form more binaries

    Zi-Yan Yuwen, Cristian Joana, Shao-Jiang Wang, and Rong-Gen Cai. Bubbles kick off primordial black holes to form more binaries. 6 2024. arXiv:2406.05838

    2024 arXiv

  40. [48]

    Herdeiro, Eugen Radu, Nicolas Sanchis-Gual, and Jos´ e A

    Carlos A.R. Herdeiro, Eugen Radu, Nicolas Sanchis-Gual, and Jos´ e A. Font. Spontaneous Scalarization of Charged Black Holes. Phys. Rev. Lett. , 121(10):101102, 2018. arXiv:1806.05190, doi:10.1103/ PhysRevLett.121.101102

    2018 arXiv

  41. [2006]

    arXiv:gr-qc/0511103, doi:10.1103/PhysRevLett.96.111102

    arXiv

  42. [2009]

    arXiv:0902.3652, doi:10.1103/PhysRevD.79.104029

    arXiv

  43. [2012]

    arXiv:1106.2254, doi:10.1103/PhysRevD.85.064040

    arXiv

  44. [2013]

    arXiv:1307.7391, doi:10.1103/PhysRevD.88.064049

    arXiv

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.