Shear-free Null Quasi-Spherical Spacetimes

The residual gauge freedom within the null quasi-spherical coordinate condition is studied, for spacetimes admitting an expanding, shear-free null foliation. The freedom consists of a boost and rotation at each coordinate sphere, corresponding to a specification of inertial frame at each sphere. Explicit formulae involving arbitrary functions of two variables are obtained for the accelerated Minkowski, Schwarzschild, and Robinson-Trautman spacetimes. These examples will be useful as test metrics in numerical relativity.Comment: 20 pages, revte

Einstein equations in the null quasi-spherical gauge III: numerical algorithms

We describe numerical techniques used in the construction of our 4th order evolution for the full Einstein equations, and assess the accuracy of representative solutions. The code is based on a null gauge with a quasi-spherical radial coordinate, and simulates the interaction of a single black hole with gravitational radiation. Techniques used include spherical harmonic representations, convolution spline interpolation and filtering, and an RK4 "method of lines" evolution. For sample initial data of "intermediate" size (gravitational field with 19% of the black hole mass), the code is accurate to 1 part in 10^5, until null time z=55 when the coordinate condition breaks down.Comment: Latex, 38 pages, 29 figures (360Kb compressed

On a Localized Riemannian Penrose Inequality

Consider a compact, orientable, three dimensional Riemannian manifold with boundary with nonnegative scalar curvature. Suppose its boundary is the disjoint union of two pieces: the horizon boundary and the outer boundary, where the horizon boundary consists of the unique closed minimal surfaces in the manifold and the outer boundary is metrically a round sphere. We obtain an inequality relating the area of the horizon boundary to the area and the total mean curvature of the outer boundary. Such a manifold may be thought as a region, surrounding the outermost apparent horizons of black holes, in a time-symmetric slice of a space-time in the context of general relativity. The inequality we establish has close ties with the Riemannian Penrose Inequality, proved by Huisken and Ilmanen, and by Bray.Comment: 16 page

Einstein equations in the null quasi-spherical gauge

The structure of the full Einstein equations in a coordinate gauge based on expanding null hypersurfaces foliated by metric 2-spheres is explored. The simple form of the resulting equations has many applications -- in the present paper we describe the structure of timelike boundary conditions; the matching problem across null hypersurfaces; and the propagation of gravitational shocks.Comment: 12 pages, LaTeX (revtex, amssymb), revision 18 pages, contains expanded discussion and explanations, updated references, to appear in CQ

Trapped Surfaces in Vacuum Spacetimes

An earlier construction by the authors of sequences of globally regular, asymptotically flat initial data for the Einstein vacuum equations containing trapped surfaces for large values of the parameter is extended, from the time symmetric case considered previously, to the case of maximal slices. The resulting theorem shows rigorously that there exists a large class of initial configurations for non-time symmetric pure gravitational waves satisfying the assumptions of the Penrose singularity theorem and so must have a singularity to the future.Comment: 14 page

Far-from-constant mean curvature solutions of Einstein's constraint equations with positive Yamabe metrics

In this article we develop some new existence results for the Einstein constraint equations using the Lichnerowicz-York conformal rescaling method. The mean extrinsic curvature is taken to be an arbitrary smooth function without restrictions on the size of its spatial derivatives, so that it can be arbitrarily far from constant. The rescaled background metric belongs to the positive Yamabe class, and the freely specifiable part of the data given by the traceless-transverse part of the rescaled extrinsic curvature and the matter fields are taken to be sufficiently small, with the matter energy density not identically zero. Using topological fixed-point arguments and global barrier constructions, we then establish existence of solutions to the constraints. Two recent advances in the analysis of the Einstein constraint equations make this result possible: A new type of topological fixed-point argument without smallness conditions on spatial derivatives of the mean extrinsic curvature, and a new construction of global super-solutions for the Hamiltonian constraint that is similarly free of such conditions on the mean extrinsic curvature. For clarity, we present our results only for strong solutions on closed manifolds. However, our results also hold for weak solutions and for other cases such as compact manifolds with boundary; these generalizations will appear elsewhere. The existence results presented here for the Einstein constraints are apparently the first such results that do not require smallness conditions on spatial derivatives of the mean extrinsic curvature.Comment: 4 pages, no figures, accepted for publication in Physical Review Letters. (Abstract shortenned and other minor changes reflecting v4 version of arXiv:0712.0798

Numerical Bifurcation Analysis of Conformal Formulations of the Einstein Constraints

The Einstein constraint equations have been the subject of study for more than fifty years. The introduction of the conformal method in the 1970's as a parameterization of initial data for the Einstein equations led to increased interest in the development of a complete solution theory for the constraints, with the theory for constant mean curvature (CMC) spatial slices and closed manifolds completely developed by 1995. The first general non-CMC existence result was establish by Holst et al. in 2008, with extensions to rough data by Holst et al. in 2009, and to vacuum spacetimes by Maxwell in 2009. The non-CMC theory remains mostly open; moreover, recent work of Maxwell on specific symmetry models sheds light on fundamental non-uniqueness problems with the conformal method as a parameterization in non-CMC settings. In parallel with these mathematical developments, computational physicists have uncovered surprising behavior in numerical solutions to the extended conformal thin sandwich formulation of the Einstein constraints. In particular, numerical evidence suggests the existence of multiple solutions with a quadratic fold, and a recent analysis of a simplified model supports this conclusion. In this article, we examine this apparent bifurcation phenomena in a methodical way, using modern techniques in bifurcation theory and in numerical homotopy methods. We first review the evidence for the presence of bifurcation in the Hamiltonian constraint in the time-symmetric case. We give a brief introduction to the mathematical framework for analyzing bifurcation phenomena, and then develop the main ideas behind the construction of numerical homotopy, or path-following, methods in the analysis of bifurcation phenomena. We then apply the continuation software package AUTO to this problem, and verify the presence of the fold with homotopy-based numerical methods.Comment: 13 pages, 4 figures. Final revision for publication, added material on physical implication