2010/01/31 by Miguel Zilhão, Miguel Zilhao, Helvi Witek +6 · 65 citations
Mathematics · Physics and Astronomy · #Algorithm #Black Holes and Theoretical Physics #Cosmology and Gravitation Theories #Einstein field equations #General relativity #Mathematical physics #Mathematics #Numerical relativity #Particle physics theoretical and experimental studies #Physics #Theoretical physics #astro-ph.HE #gr-qc #hep-ph #hep-th
paper · pdf · doi:10.1103/physrevd.81.084052
published in Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields 81(8) (American Physical Society) · 31 pages, 6 figures; v2 Minor changes and added two references. Matches the published version in PRD.
openalex publication_date 2010/04/29 · arxiv created 2010/05/03 · arxiv updated 2014/11/20 · openalex created_date 2016/06/24 · openalex updated_date 2026/08/05
The numerical evolution of Einstein's field equations in a generic background has the potential to answer a variety of important questions in physics: from applications to the gauge-gravity duality, to modeling black hole production in TeV gravity scenarios, to analysis of the stability of exact solutions, and to tests of cosmic censorship. In order to investigate these questions, we extend numerical relativity to more general space-times than those investigated hitherto, by developing a framework to study the numerical evolution of D dimensional vacuum space-times with an SO(D\ensuremath-2) isometry group for D\ensuremath≥5, or SO(D\ensuremath-3) for D\ensuremath≥6. Performing a dimensional reduction on a (D\ensuremath-4) sphere, the D dimensional vacuum Einstein equations are rewritten as a 3+1 dimensional system with source terms, and presented in the Baumgarte, Shapiro, Shibata, and Nakamura formulation. This allows the use of existing 3+1 dimensional numerical codes with small adaptations. Brill-Lindquist initial data are constructed in D dimensions and a procedure to match them to our 3+1 dimensional evolution equations is given. We have implemented our framework by adapting the Lean code and perform a variety of simulations of nonspinning black hole space-times. Specifically, we present a modified moving puncture gauge, which facilitates long-term stable simulations in D=5. We further demonstrate the internal consistency of the code by studying convergence and comparing numerical versus analytic results in the case of geodesic slicing for D=5, 6.