Wellposed constrained evolution of 3+1 formulations of general relativity
Abstract
We present an analysis of wellposedness of constrained evolution of 3+1 formulations of general relativity. In this analysis we explicitly take into account the energy and momentum constraints as well as possible algebraic constraints on the evolution of highfrequency perturbations of solutions of Einstein's equations. In this respect, our approach is principally different from standard analyses of wellposedness of free evolution in general relativity. Our study reveals the existence of subsets of the linearized Einstein's equations that control the wellposedness of constrained evolution. It is demonstrated that the wellposedness of ArnowittDeserMisner (ADM), BaumgarteShapiroShibataNakamura and other 3+1 formulations derived from the ADM formulation by adding combinations of constraints to the righthand side of the ADM formulation and/or by linear transformation of the dynamical ADM variables depends entirely on the properties of the gauge. For certain classes of gauges we formulate conditions for wellposedness of constrained evolution. This provides a new basis for constructing stable numerical integration schemes for a classical ADM and many other 3+1 formulations of general relativity.
 Authors:
 Department of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 (United States)
 (United States)
 Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 (United States)
 (Denmark)
 (Russian Federation)
 Publication Date:
 OSTI Identifier:
 21010896
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Physical Review. D, Particles Fields; Journal Volume: 75; Journal Issue: 2; Other Information: DOI: 10.1103/PhysRevD.75.024026; (c) 2007 The American Physical Society; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; COSMOLOGY; DISTURBANCES; EINSTEIN FIELD EQUATIONS; FOURDIMENSIONAL CALCULATIONS; GAUGE INVARIANCE; GENERAL RELATIVITY THEORY; MATHEMATICAL SOLUTIONS; TRANSFORMATIONS
Citation Formats
Paschalidis, Vasileios, Khokhlov, Alexei, Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, Novikov, Igor, Niels Bohr Institute, Blegdamsvej 17, DK2100 Copenhagen, and Astro Space Center of P.N. Lebedev Physical Institute, Profsoyuznaya 84/32, Moscow, 117810. Wellposed constrained evolution of 3+1 formulations of general relativity. United States: N. p., 2007.
Web. doi:10.1103/PHYSREVD.75.024026.
Paschalidis, Vasileios, Khokhlov, Alexei, Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, Novikov, Igor, Niels Bohr Institute, Blegdamsvej 17, DK2100 Copenhagen, & Astro Space Center of P.N. Lebedev Physical Institute, Profsoyuznaya 84/32, Moscow, 117810. Wellposed constrained evolution of 3+1 formulations of general relativity. United States. doi:10.1103/PHYSREVD.75.024026.
Paschalidis, Vasileios, Khokhlov, Alexei, Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637, Novikov, Igor, Niels Bohr Institute, Blegdamsvej 17, DK2100 Copenhagen, and Astro Space Center of P.N. Lebedev Physical Institute, Profsoyuznaya 84/32, Moscow, 117810. Mon .
"Wellposed constrained evolution of 3+1 formulations of general relativity". United States.
doi:10.1103/PHYSREVD.75.024026.
@article{osti_21010896,
title = {Wellposed constrained evolution of 3+1 formulations of general relativity},
author = {Paschalidis, Vasileios and Khokhlov, Alexei and Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, Illinois 60637 and Novikov, Igor and Niels Bohr Institute, Blegdamsvej 17, DK2100 Copenhagen and Astro Space Center of P.N. Lebedev Physical Institute, Profsoyuznaya 84/32, Moscow, 117810},
abstractNote = {We present an analysis of wellposedness of constrained evolution of 3+1 formulations of general relativity. In this analysis we explicitly take into account the energy and momentum constraints as well as possible algebraic constraints on the evolution of highfrequency perturbations of solutions of Einstein's equations. In this respect, our approach is principally different from standard analyses of wellposedness of free evolution in general relativity. Our study reveals the existence of subsets of the linearized Einstein's equations that control the wellposedness of constrained evolution. It is demonstrated that the wellposedness of ArnowittDeserMisner (ADM), BaumgarteShapiroShibataNakamura and other 3+1 formulations derived from the ADM formulation by adding combinations of constraints to the righthand side of the ADM formulation and/or by linear transformation of the dynamical ADM variables depends entirely on the properties of the gauge. For certain classes of gauges we formulate conditions for wellposedness of constrained evolution. This provides a new basis for constructing stable numerical integration schemes for a classical ADM and many other 3+1 formulations of general relativity.},
doi = {10.1103/PHYSREVD.75.024026},
journal = {Physical Review. D, Particles Fields},
number = 2,
volume = 75,
place = {United States},
year = {Mon Jan 15 00:00:00 EST 2007},
month = {Mon Jan 15 00:00:00 EST 2007}
}

Two new formulations of general relativity are introduced. The first one is a parabolization of the ArnowittDeserMisner formulation and is derived by the addition of combinations of the constraints and their derivatives to the righthand side of the ArnowittDeserMisner evolution equations. The desirable property of this modification is that it turns the surface of constraints into a local attractor because the constraint propagation equations become secondorder parabolic independently of the gauge conditions employed. This system may be classified as mixed hyperbolicsecondorder parabolic. The second formulation is a parabolization of the KidderScheelTeukolsky formulation and is a manifestly mixed strongly hyperbolicsecondorderparabolic setmore »

Boundary Conditions for Constrained Systems of Evolution Equations in Numerical Relativity
In numerical relativity, to study astrophysical events such as blackhole collisions, the Einstein equations for geometry of spacetime are solved as system of partial differential equations. Current simulations are usually based on socalled BSSN system of 3+1 constrained hyperbolic evolution equations for tensorial fields of various ranks. The boundary conditions for evolved fields are given by the fact that the simulated events happen in an empty space and that far from the center the waves should propagate outwards. We analyze the usual outgoing radiation boundary conditions for the BSSN variables using potentials. Introduction of these potentials simplifies the linearized BSSNmore »