Numerical stability of a new conformal-traceless 3 + 1 formulation of the Einstein equation

There is strong evidence indicating that the particular form used to recast the Einstein equation as a 3 + 1 set of evolution equations has a fundamental impact on the stability properties of numerical evolutions involving black holes and/or neutron stars. Presently, the longest lived evolutions have been obtained using a parametrized hyperbolic system developed by Kidder, Scheel and Teukolsky or a conformal-traceless system introduced by Baumgarte, Shapiro, Shibata and Nakamura. We present a new conformal-traceless system. While this new system has some elements in common with the Baumgarte Shapiro Shibata Nakamura system, it differs in both the type of conformal transformations and how the nonlinear terms involving the extrinsic curvature are handled. We show results from 3D numerical evolutions of a single, non-rotating black hole in which we demonstrate that this new system yields a significant improvement in the lifetime of the simulations.