Constraint propagation in N+1-dimensional space-time

Higher dimensional space-time models provide us an alternative interpretation of nature, and give us different dynamical aspects than the traditional four-dimensional space-time models. Motivated by such recent interests, especially for future numerical research of higher-dimensional space-time, we study the dimensional dependence of constraint propagation behavior. The $N+1$ Arnowitt-Deser-Misner evolution equation has matter terms which depend on $N$, but the constraints and constraint propagation equations remain the same. This indicates that there would be problems with accuracy and stability when we directly apply the $N+1$ ADM formulation to numerical simulations as we have experienced in four-dimensional cases. However, we also conclude that previous efforts in re-formulating the Einstein equations can be applied if they are based on constraint propagation analysis.Comment: 4 pages, to appear in Gen. Rel. Gra

Diagonalizability of Constraint Propagation Matrices

In order to obtain stable and accurate general relativistic simulations, re-formulations of the Einstein equations are necessary. In a series of our works, we have proposed using eigenvalue analysis of constraint propagation equations for evaluating violation behavior of constraints. In this article, we classify asymptotical behaviors of constraint-violation into three types (asymptotically constrained, asymptotically bounded, and diverge), and give their necessary and sufficient conditions. We find that degeneracy of eigenvalues sometimes leads constraint evolution to diverge (even if its real-part is not positive), and conclude that it is quite useful to check the diagonalizability of constraint propagation matrices. The discussion is general and can be applied to any numerical treatments of constrained dynamics.Comment: 4 pages, RevTeX, one figure, added one paragraph in concluding remarks. The version to appear in Class. Quant. Grav. (Lett

Formulations of the Einstein equations for numerical simulations

We review recent efforts to re-formulate the Einstein equations for fully relativistic numerical simulations. In order to complete a long-term and accurate simulations of binary compact objects, people seek a robust set of equations against the violation of constraints. Many trials have revealed that mathematically equivalent sets of evolution equations show different numerical stability in free evolution schemes. In this article, we overview the efforts of the community, categorizing them into three directions: (1) modifications of the standard Arnowitt-Deser-Misner equations initiated by the Kyoto group (the so-called Baumgarte-Shapiro-Shibata-Nakamura equations), (2) rewriting the evolution equations in a hyperbolic form, and (3) construction of an "asymptotically constrained" system. We then introduce our series of works that tries to explain these evolution behaviors in a unified way using eigenvalue analysis of the constraint propagation equations. The modifications of (or adjustments to) the evolution equations change the character of constraint propagation, and several particular adjustments using constraints are expected to damp the constraint-violating modes. We show several set of adjusted ADM/BSSN equations, together with their numerical demonstrations.Comment: 15 pages, 12 figures; Invited Lecture at APCTP Winter School on Black Hole Astrophysics, Daejeon and Pohang, Korea, January 24-29, 200