(gr-qc/0204002) to appear in Phys. Rev. D.)Several numerical relativity groups are using a modified ADM formulation for their simulations, which was developed by Nakamura et al (and widely cited as Baumgarte-Shapiro-Shibata-Nakamura system). This so-called BSSN formulation is shown to be more stable than the standard ADM formulation in many cases, and there have been many attempts to explain why this re-formulation has such an advantage. We try to explain the background mechanism of the BSSN equations by using eigenvalue analysis of constraint propagation equations. This analysis has been applied and has succeeded in explaining other systems in our series of works. We derive the full set of the constraint propagation equations, and study it in the flat background space-time. We carefully examine how the replacements and adjustments in the equations change the propagation structure of the constraints, i.e. whether violation of constraints (if it exists) will decay or propagate away. We conclude that the better stability of the BSSN system is obtained by their adjustments in the equations, and that the combination of the adjustments is in a good balance, i.e. a lack of their adjustments might fail to obtain the present stability. We further propose other adjustments to the equations, which may offer more stable features than the current BSSN equations.

The current important issue in numerical relativity is to determine which formulation of the Einstein equations provides us with stable and accurate simulations. Based on our previous work on "asymptotically constrained" systems, we here present constraint propagation equations and their eigenvalues for the Arnowitt-Deser-Misner (ADM) evolution equations with additional constraint terms (adjusted terms) on the right hand side. We conjecture that the system is robust against violation of constraints if the amplification factors (eigenvalues of Fourier-component of the constraint propagation equations) are negative or pure-imaginary. We show such a system can be obtained by choosing multipliers of adjusted terms. Our discussion covers Detweiler's proposal (1987) and Frittelli's analysis (1997), and we also mention the so-called conformal-traceless ADM systems.

revised version) gr-qc/0007034 to appear in Class. Quant. Grav.)We study asymptotically constrained systems for numerical integration of the Einstein equations, which are intended to be robust against perturbative errors for the free evolution of the initial data. First, we examine the previously proposed "λ-system", which introduces artificial flows to constraint surfaces based on the symmetric hyperbolic formulation. We show that this system works as expected for the wave propagation problem in the Maxwell system and in general relativity using Ashtekar's connection formulation. Second, we propose a new mechanism to control the stability, which we call the "adjusted system". This is simply obtained by adding constraint terms in the dynamical equations and adjusting its multipliers. We explain why a particular choice of multiplier reduces the numerical errors from non-positive or pure-imaginary eigenvalues of the adjusted constraint propagation equations. This "adjusted system" is also tested in the Maxwell system and in the Ashtekar's system. This mechanism affects more than the system's symmetric hyperbolicity.

In order to perform accurate and stable long-time numerical integration of the Einstein equation, several hyperbolic systems have been proposed. We here present numerical comparisons between weakly hyperbolic, strongly hyperbolic, and symmetric hyperbolic systems based on Ashtekar's connection variables. The primary advantage for using this connection formulation in this experiment is that we can keep using the same dynamical variables for all levels of hyperbolicity. Our numerical code demonstrates gravitational wave propagation in plane symmetric spacetimes, and we compare the accuracy of the simulation by monitoring the violation of the constraints. By comparing with results obtained from the weakly hyperbolic system, we observe the strongly and symmetric hyperbolic system show better numerical performance (yield less constraint violation), but not so much difference between the latter two. Rather, we find that the symmetric hyperbolic system is not always the best in numerical performances.This study is the premier to present full numerical simulations using Ashtekar's variables. We also describe our procedures in detail.gr-qc/0005003, revised version

In order to find a way to have a better formulation for numerical evolution of the Einstein equations, we study the propagation equations of the constraints based on the Arnowitt-Deser-Misner formulation. By adjusting constraint terms in the evolution equations, we try to construct an "asymptotically constrained system" which is expected to be robust against violation of the constraints, and to enable a long-term stable and accurate numerical simulation. We first provide useful expressions for analyzing constraint propagation in a general spacetime, then apply it to Schwarzschild spacetime. We search when and where the negative real or non-zero imaginary eigenvalues of the homogenized constraint propagation matrix appear, and how they depend on the choice of coordinate system and adjustments. Our analysis includes the proposal of Detweiler (1987), which is still the best one according to our conjecture but has a growing mode of error near the horizon. Some examples are snapshots of a maximally sliced Schwarzschild black hole. The predictions here may help the community to make further improvements. * Electronic address: shinkai@atlas.riken.go.jp † Electronic address: yoneda@mn.waseda.ac.jp