Gauge fixing for the simulation of black hole spacetimes

I consider the initial-boundary-value-problem of nonlinear general relativistic vacuum spacetimes, which today cannot yet be evolved numerically in a satisfactory manner. Specifically, I look at gauge conditions, classifying them into gauge evolution conditions and gauge fixing conditions. In this terminology, a gauge fixing condition is a condition that removes all gauge degrees of freedom from a system, whereas a gauge evolution condition determines only the time evolution of the gauge condition, while the gauge condition itself remains unspecified. I find that most of today's gauge conditions are only gauge evolution conditions. I present a system of evolution equations containing a gauge fixing condition, and describe an efficient numerical implementation using constrained evolution. I examine the numerical behaviour of this system for several test problems, such as linear gravitational waves or nonlinear gauge waves. I find that the system is robustly stable and second-order convergent. I then apply it to more realistic configurations, such as Brill waves or single black holes, where the system is also stable and accurate.

[1]  Ken-ichi Oohara,et al.  General Relativistic Collapse to Black Holes and Gravitational Waves from Black Holes , 1987 .

[2]  S. Shapiro,et al.  On the numerical integration of Einstein's field equations , 1998, gr-qc/9810065.

[3]  Marsa,et al.  Black-hole-scalar-field interactions in spherical symmetry. , 1996, Physical review. D, Particles and fields.

[4]  Olaf Dreyer,et al.  Introduction to isolated horizons in numerical relativity , 2003 .

[5]  Saul A. Teukolsky.  Stability of the iterated Crank-Nicholson method in numerical relativity , 2000 .

[6]  J. Thornburg.  Numerical relativity in black hole spacetimes , 1993 .

[7]  Seidel,et al.  Three-dimensional numerical relativity: The evolution of black holes. , 1995, Physical review. D, Particles and fields.

[8]  Jacques E. Romain,et al.  Gravitation: An Introduction to Current Research , 1963 .

[9]  Cauchy boundaries in linearized gravitational theory , 1999, gr-qc/9912030.

[10]  Lawrence E. Kidder,et al.  Black hole evolution by spectral methods , 2000, gr-qc/0005056.

[11]  A cure for unstable numerical evolutions of single black holes: adjusting the standard ADM equations , 2001, gr-qc/0103099.

[12]  Choptuik.  Consistency of finite-difference solutions of Einstein's equations. , 1991, Physical review. D, Particles and fields.

[13]  Gabrielle Allen,et al.  Cactus Grid Computing: Review of Current Development , 2001, Euro-Par.

[14]  Kerr?Schild-type initial data for black holes with angular momenta , 2002, gr-qc/0205060.

[15]  Bernd Bruegmann,et al.  A Simple Construction of Initial Data for Multiple Black Holes , 1997 .

[16]  D. Shoemaker,et al.  Grazing collisions of black holes via the excision of singularities. , 2000, Physical review letters.

[17]  Boundary conditions in linearized harmonic gravity , 2001, gr-qc/0106026.

[18]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[19]  Long Chen.  INTRODUCTION TO MULTIGRID METHODS , 2005 .

[20]  C. Lousto,et al.  The Lazarus project : A pragmatic approach to binary black hole , 2001, gr-qc/0104063.

[21]  William Gropp,et al.  Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries , 1997, SciTools.

[22]  Dieter R. Brill,et al.  On the positive definite mass of the Bondi-Weber-Wheeler time-symmetric gravitational waves , 1959 .

[23]  Finding apparent horizons in numerical relativity. , 1995, Physical review. D, Particles and fields.

[24]  John Shalf,et al.  Cactus Tools for Grid Applications , 2001, Cluster Computing.

[25]  Generic tracking of multiple apparent horizons with level flow , 2000, gr-qc/0004062.

[26]  J. York.  Conformally invariant orthogonal decomposition of symmetric tensors on Riemannian manifolds and the initial‐value problem of general relativity , 1973 .

[27]  Richard A. Matzner,et al.  Initial data and coordinates for multiple black hole systems , 1998 .

[28]  J. Frauendiener.  Numerical treatment of the hyperboloidal initial value problem for the vacuum Einstein equations. III. On the determination of radiation , 1998, gr-qc/9808072.

[29]  M. F.,et al.  Bibliography , 1985, Experimental Gerontology.

[30]  Brian D. Fisher,et al.  University of British Columbia , 2002, INTR.

[31]  J. York.  Covariant decompositions of symmetric tensors in the theory of gravitation , 1974 .

[32]  Larry Smarr,et al.  Kinematical conditions in the construction of spacetime , 1978 .

[33]  Locating boosted Kerr and Schwarzschild apparent horizons , 2000, gr-qc/0002076.

[34]  Apparent horizon finder for general three-dimensional spaces , 2000 .

[35]  Nakamura,et al.  Evolution of three-dimensional gravitational waves: Harmonic slicing case. , 1995, Physical review. D, Particles and fields.

[36]  Re-formulating the Einstein equations for stable numerical simulations -F ormulation Problem in Numerical Relativity - , 2002, gr-qc/0209111.

[37]  D. Raine.  General relativity , 1980, Nature.