Olivier Sarbach

Instability of wormholes supported by a ghost scalar field: I. Linear stability analysis

We examine the linear stability of static, spherically symmetric wormhole solutions of Einsteins field equations coupled to a massless ghost scalar field. These solutions are parametrized by the areal radius of their throat and the product of the masses at their asymptotically flat ends. We prove that all these solutions are unstable with respect to linear fluctuations and possess precisely one unstable, exponentially in time growing mode. The associated time scale is shown to be of the order of the wormhole throat divided by the speed of light. The nonlinear evolution is analyzed in a subsequent article.

Turduckening black holes: An Analytical and computational study

We provide a detailed analysis of several aspects of the turduckening technique for evolving black holes. At the analytical level we study the constraint propagation for a family of formulations of Einsteins field equations and identify under what conditions the turducken procedure is rigorously justified and under what conditions constraint violations will propagate to the outside of the black holes. We present high resolution spherically symmetric studies which verify our analytical predictions. Then we present three-dimensional simulations of single distorted black holes using different variations of the turduckening method and also the puncture method. We study the effect that these different methods have on the coordinate conditions, constraint violations, and extracted gravitational waves. We find that the waves agree up to small but nonvanishing differences, caused by escaping superluminal gauge modes. These differences become smaller with increasing detector location.

Well Posed Constraint-Preserving Boundary Conditions for the Linearized Einstein Equations

In this paper we address the problem of specifying boundary conditions for Einsteins equations when linearized around Minkowski space using the generalized Einstein-Christoffel symmetric hyperbolic system of evolution equations. The boundary conditions we work out guarantee that the constraints are satisfied provided they are satisfied on the initial slice and ensures a well posed initial-boundary value formulation. We consider the case of a manifold with a non-smooth boundary, as is the usual case of the cubic boxes commonly used in numerical relativity. The techniques discussed should be applicable to more general cases, as linearizations around more complicated backgrounds, and may be used to establish well posedness in the full non-linear case.

Hyperbolicity of the Baumgarte-Shapiro-Shibata-Nakamura system of Einstein evolution equations

We discuss an equivalence between the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formulation of the Einstein evolution equations, a subfamiliy of the Kidder--Scheel--Teukolsky formulation, and other strongly or symmetric hyperbolic first order systems with fixed shift and densitized lapse. This allows us to show under which conditions the BSSN system is, in a sense to be discussed, hyperbolic. This desirable property may account in part for the empirically observed better behavior of the BSSN formulation in numerical evolutions involving black holes.

Instability of wormholes supported by a ghost scalar field. II. Nonlinear evolution

We analyze the nonlinear evolution of spherically symmetric wormhole solutions coupled to a massless ghost scalar field using numerical methods. In a previous article, we have shown that static wormholes with these properties are unstable with respect to linear perturbations. Here, we show that depending on the initial perturbation the wormholes either expand or decay to a Schwarzschild black hole. We estimate the time scale of the expanding solutions and those collapsing to a black hole, and show that they are consistent in the regime of small perturbations with those predicted from perturbation theory. In the collapsing case, we also present a systematic study of the final black hole horizon and discuss the possibility for a luminous signal to travel from one universe to the other and back before the black hole forms. In the expanding case, the wormholes seem to undergo an exponential expansion, at least during the run time of our simulations.

On the well posedness of the Baumgarte-Shapiro-Shibata-Nakamura formulation of Einstein's field equations

We give a well posed initial value formulation of the Baumgarte-Shapiro-Shibata-Nakamura form of Einsteins equations with gauge conditions given by a Bona-Masso-like slicing condition for the lapse and a frozen shift. This is achieved by introducing extra variables and recasting the evolution equations into a first order symmetric hyperbolic system. We also consider the presence of artificial boundaries and derive a set of boundary conditions that guarantee that the resulting initial-boundary value problem is well posed, though not necessarily compatible with the constraints. In the case of dynamical gauge conditions for the lapse and shift we obtain a class of evolution equations which are strongly hyperbolic and so yield well posed initial value formulations.

Continuum and Discrete Initial-Boundary-Value Problems and Einstein's Field Equations

Many evolution problems in physics are described by partial differential equations on an infinite domain; therefore, one is interested in the solutions to such problems for a given initial dataset. A prominent example is the binary black-hole problem within Einstein’s theory of gravitation, in which one computes the gravitational radiation emitted from the inspiral of the two black holes, merger and ringdown. Powerful mathematical tools can be used to establish qualitative statements about the solutions, such as their existence, uniqueness, continuous dependence on the initial data, or their asymptotic behavior over large time scales. However, one is often interested in computing the solution itself, and unless the partial differential equation is very simple, or the initial data possesses a high degree of symmetry, this computation requires approximation by numerical discretization. When solving such discrete problems on a machine, one is faced with a finite limit to computational resources, which leads to the replacement of the infinite continuum domain with a finite computer grid. This, in turn, leads to a discrete initial-boundary value problem. The hope is to recover, with high accuracy, the exact solution in the limit where the grid spacing converges to zero with the boundary being pushed to infinity.The goal of this article is to review some of the theory necessary to understand the continuum and discrete initial boundary-value problems arising from hyperbolic partial differential equations and to discuss its applications to numerical relativity; in particular, we present well-posed initial and initial-boundary value formulations of Einstein’s equations, and we discuss multi-domain high-order finite difference and spectral methods to solve them.

Excision without excision

to turducken (turduckens, turduckening, turduckened, turduckened) [math.]: To stuff a black hole. We analyze and apply an alternative to black hole excision based on smoothing the interior of black holes with arbitrary initial data, and solving the vacuum Einstein evolution equations everywhere. By deriving the constraint propagation system for our hyperbolic formulation of the BSSN evolution system we rigorously prove that the constraints propagate causally and so any constraint violations introduced inside the black holes cannot affect the exterior spacetime. We present evolutions of Cook-Pfeiffer binary black hole initial configurations showing that these techniques appear to work robustly for generic data. We also present evidence from spherically symmetric evolutions that for the gauge conditions used the same stationary end-state is approached irrespective of the choice of initial data and smoothing procedure.

Convergence and stability in numerical relativity

It is often the case in numerical relativity that schemes that are known to be convergent for well posed systems are used in evolutions of weakly hyperbolic (WH) formulations of Einsteins equations. Here we explicitly show that with several of the discretizations that have been used through out the years, this procedure leads to non-convergent schemes. That is, arbitrarily small initial errors are amplified without bound when resolution is increased, independently of the amount of numerical dissipation introduced. The lack of convergence introduced by this instability can be particularly subtle, in the sense that it can be missed by several convergence tests, especially in 3+1 dimensional codes. We propose tests and methods to analyze convergence that may help detect these situations.

No naked singularities in homogeneous, spherically symmetric bubble spacetimes?

We study the evolution of bubble spacetimes in vacuum and electrovac scenarios by numerical means. We find strong evidence against the formation of naked singularities in (i) scenarios with negative masses displaying initially collapsing conditions and (ii) scenarios with negative masses displaying initially expanding conditions, previously reported to give rise to such singularities. Additionally, we show that the presence of strong gauge fields implies that an initially collapsing bubble bounces back and expands. By fine-tuning the strength of the gauge field we find that the solution approaches a static bubble solution.