Stable long-term evolution in numerical relativity

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