No Arabic abstract
This paper gives a detailed pedagogic presentation of the central concepts underlying a new algorithm for the numerical solution of Einsteins equations for gravitation. This approach incorporates the best features of the two leading approaches to computational gravitation, carving up spacetime via Cauchy hypersurfaces within a central worldtube, and using characteristic hypersurfaces in its exterior to connect this region with null infinity and study gravitational radiation. It has worked well in simplified test problems, and is currently being used to build computer codes to simulate black hole collisions in 3-D.
From Einsteins theory we know that besides the electromagnetic spectrum, objects like quasars, active galactic nuclei, pulsars and black holes also generate a physical signal of purely gravitational nature. The actual form of the signal is impossible to determine analytically, which lead to use of numerical methods. Two major approaches emerged. The first one formulates the gravitational radiation problem as a standard Cauchy initial value problem, while the other approach uses a Characteristic Initial value formulation. In the strong field region, where caustics in the wavefronts are likely to form, the Cauchy formulation is more advantageous. On the other side, the Characteristic formulation is uniquely suited to study radiation problems because it describes space-time in terms of radiation wavefronts. The fact that the advantages and disadvantages of these two systems are complementary suggests that one may want to use the two of them together. In a full nonlinear problem it would be advantageous to evolve the inner (strong field) region using Cauchy evolution and the outer (radiation) region with the Characteristic approach. Cauchy Characteristic Matching enables one to evolve the whole space-time matching the boundaries of Cauchy and Characteristic evolution. The methodology of Cauchy Characteristic Matching has been successful in numerical evolution of the spherically symmetric Klein-Gordon-Einstein field equations as well as for 3-D non-linear wave equations. In this thesis the same methodology is studied in the context of the Einstein equations.
We present several improvements to the Cauchy-characteristic evolution procedure that generates high-fidelity gravitational waveforms at $mathcal{I}^+$ from numerical relativity simulations. Cauchy-characteristic evolution combines an interior solution of the Einstein field equations based on Cauchy slices with an exterior solution based on null slices that extend to $mathcal{I}^+$. The foundation of our improved algorithm is a comprehensive method of handling the gauge transformations between the arbitrarily specified coordinates of the interior Cauchy evolution and the unique (up to BMS transformations) Bondi-Sachs coordinate system of the exterior characteristic evolution. We present a reformulated set of characteristic evolution equations better adapted to numerical implementation. In addition, we develop a method to ensure that the angular coordinates used in the volume during the characteristic evolution are asymptotically inertial. This provides a direct route to an expanded set of waveform outputs and is guaranteed to avoid pure-gauge logarithmic dependence that has caused trouble for previous spectral implementations of the characteristic evolution equations. We construct a set of Weyl scalars compatible with the Bondi-like coordinate systems used in characteristic evolution, and determine simple, easily implemented forms for the asymptotic Weyl scalars in our suggested set of coordinates.
We present a new method of extracting gravitational radiation from three-dimensional numerical relativity codes and providing outer boundary conditions. Our approach matches the solution of a Cauchy evolution of Einsteins equations to a set of one-dimensional linear wave equations on a curved background. We illustrate the mathematical properties of our approach and discuss a numerical module we have constructed for this purpose. This module implements the perturbative matching approach in connection with a generic three-dimensional numerical relativity simulation. Tests of its accuracy and second-order convergence are presented with analytic linear wave data.
We extract gravitational waveforms from numerical simulations of black hole binaries computed using the Spectral Einstein Code. We compare two extraction methods: direct construction of the Newman-Penrose (NP) scalar $Psi_4$ at a finite distance from the source and Cauchy-characteristic extraction (CCE). The direct NP approach is simpler than CCE, but NP waveforms can be contaminated by near-zone effects---unless the waves are extracted at several distances from the source and extrapolated to infinity. Even then, the resulting waveforms can in principle be contaminated by gauge effects. In contrast, CCE directly provides, by construction, gauge-invariant waveforms at future null infinity. We verify the gauge invariance of CCE by running the same physical simulation using two different gauge conditions. We find that these two gauge conditions produce the same CCE waveforms but show differences in extrapolated-$Psi_4$ waveforms. We examine data from several different binary configurations and measure the dominant sources of error in the extrapolated-$Psi_4$ and CCE waveforms. In some cases, we find that NP waveforms extrapolated to infinity agree with the corresponding CCE waveforms to within the estimated error bars. However, we find that in other cases extrapolated and CCE waveforms disagree, most notably for $m=0$ memory modes.
Recently, folk questions on the smoothability of Cauchy hypersurfaces and time functions of a globally hyperbolic spacetime M, have been solved. Here we give further results, applicable to several problems: (1) Any compact spacelike acausal submanifold H with boundary can be extended to a spacelike Cauchy hypersurface S. If H were only achronal, counterexamples to the smooth extension exist, but a continuous extension (in fact, valid for any compact achronal subset K) is still possible. (2) Given any spacelike Cauchy hypersurface S, a Cauchy temporal function T (i.e., a smooth function with past-directed timelike gradient everywhere, and Cauchy hypersurfaces as levels) with S equal to one of the levels, is constructed -thus, the spacetime splits orthogonally as $R times S$ in a canonical way. Even more, accura