Critical bubbles and implications for critical black strings

We demonstrate the existence of gravitational critical phenomena in higher dimensional electrovac bubble spacetimes. To this end, we study linear fluctuations about families of static, homogeneous spherically symmetric bubble spacetimes in Kaluza-Klein theories coupled to a Maxwell field. We prove that these solutions are linearly unstable and posses a unique unstable mode with a growth rate that is universal in the sense that it is independent of the family considered. Furthermore, by a double analytical continuation this mode can be seen to correspond to marginally stable stationary modes of perturbed black strings whose periods are integer multiples of the Gregory-Laflamme critical length. This allow us to rederive recent results about the behavior of the critical mass for large dimensions and to generalize them to the charged black string case.Comment: A reference to unpublished work for the case q=2, by J. Hovdebo adde

The discrete energy method in numerical relativity: Towards long-term stability

The energy method can be used to identify well-posed initial boundary value problems for quasi-linear, symmetric hyperbolic partial differential equations with maximally dissipative boundary conditions. A similar analysis of the discrete system can be used to construct stable finite difference equations for these problems at the linear level. In this paper we apply these techniques to some test problems commonly used in numerical relativity and observe that while we obtain convergent schemes, fast growing modes, or ``artificial instabilities,'' contaminate the solution. We find that these growing modes can partially arise from the lack of a Leibnitz rule for discrete derivatives and discuss ways to limit this spurious growth.Comment: 18 pages, 22 figure

3D simulations of Einstein's equations: symmetric hyperbolicity, live gauges and dynamic control of the constraints

We present three-dimensional simulations of Einstein equations implementing a symmetric hyperbolic system of equations with dynamical lapse. The numerical implementation makes use of techniques that guarantee linear numerical stability for the associated initial-boundary value problem. The code is first tested with a gauge wave solution, where rather larger amplitudes and for significantly longer times are obtained with respect to other state of the art implementations. Additionally, by minimizing a suitably defined energy for the constraints in terms of free constraint-functions in the formulation one can dynamically single out preferred values of these functions for the problem at hand. We apply the technique to fully three-dimensional simulations of a stationary black hole spacetime with excision of the singularity, considerably extending the lifetime of the simulations.Comment: 21 pages. To appear in PR

Towards the Final Fate of an Unstable Black String

Black strings, one class of higher dimensional analogues of black holes, were shown to be unstable to long wavelength perturbations by Gregory and Laflamme in 1992, via a linear analysis. We revisit the problem through numerical solution of the full equations of motion, and focus on trying to determine the end-state of a perturbed, unstable black string. Our preliminary results show that such a spacetime tends towards a solution resembling a sequence of spherical black holes connected by thin black strings, at least at intermediate times. However, our code fails then, primarily due to large gradients that develop in metric functions, as the coordinate system we use is not well adapted to the nature of the unfolding solution. We are thus unable to determine how close the solution we see is to the final end-state, though we do observe rich dynamical behavior of the system in the intermediate stages.Comment: 17 pages, 7 figure

Relativistic MHD and black hole excision: Formulation and initial tests

A new algorithm for solving the general relativistic MHD equations is described in this paper. We design our scheme to incorporate black hole excision with smooth boundaries, and to simplify solving the combined Einstein and MHD equations with AMR. The fluid equations are solved using a finite difference Convex ENO method. Excision is implemented using overlapping grids. Elliptic and hyperbolic divergence cleaning techniques allow for maximum flexibility in choosing coordinate systems, and we compare both methods for a standard problem. Numerical results of standard test problems are presented in two-dimensional flat space using excision, overlapping grids, and elliptic and hyperbolic divergence cleaning.Comment: 22 pages, 8 figure

Geometrically motivated hyperbolic coordinate conditions for numerical relativity: Analysis, issues and implementations

We study the implications of adopting hyperbolic driver coordinate conditions motivated by geometrical considerations. In particular, conditions that minimize the rate of change of the metric variables. We analyze the properties of the resulting system of equations and their effect when implementing excision techniques. We find that commonly used coordinate conditions lead to a characteristic structure at the excision surface where some modes are not of outflow-type with respect to any excision boundary chosen inside the horizon. Thus, boundary conditions are required for these modes. Unfortunately, the specification of these conditions is a delicate issue as the outflow modes involve both gauge and main variables. As an alternative to these driver equations, we examine conditions derived from extremizing a scalar constructed from Killing's equation and present specific numerical examples.Comment: 9 figure

On computations of angular momentum and its flux in numerical relativity

The purpose of this note is to point out ambiguities that appear in the calculation of angular momentum and its radiated counterpart when some simple formulae are used to compute them. We illustrate, in two simple different examples, how incorrect results can be obtained with them. Additionally, we discuss the magnitude of possible errors in well known situations.Comment: 8 pages. Minor improvements . To appear in Class. Quantum Gra

Numerical stability of a new conformal-traceless 3+1 formulation of the Einstein equation

There is strong evidence indicating that the particular form used to recast the Einstein equation as a 3+1 set of evolution equations has a fundamental impact on the stability properties of numerical evolutions involving black holes and/or neutron stars. Presently, the longest lived evolutions have been obtained using a parametrized hyperbolic system developed by Kidder, Scheel and Teukolsky or a conformal-traceless system introduced by Baumgarte, Shapiro, Shibata and Nakamura. We present a new conformal-traceless system. While this new system has some elements in common with the Baumgarte-Shapiro-Shibata-Nakamura system, it differs in both the type of conformal transformations and how the non-linear terms involving the extrinsic curvature are handled. We show results from 3D numerical evolutions of a single, non-rotating black hole in which we demonstrate that this new system yields a significant improvement in the life-time of the simulations.Comment: 7 pages, 2 figure

Numerical Simulations of Hyperfine Transitions of Antihydrogen

One of the ASACUSA (Atomic Spectroscopy And Collisions Using Slow Antiprotons) collaboration's goals is the measurement of the ground state hyperfine transition frequency in antihydrogen, the antimatter counterpart of one of the best known systems in physics. This high precision experiment yields a sensitive test of the fundamental symmetry of CPT. Numerical simulations of hyperfine transitions of antihydrogen atoms have been performed providing information on the required antihydrogen events and the achievable precision