No Arabic abstract
We describe an explicit in time, finite-difference code designed to simulate black holes by using the excision method. The code is based upon the harmonic formulation of the Einstein equations and incorporates several features regarding the well-posedness and numerical stability of the initial-boundary problem for the quasilinear wave equation. After a discussion of the equations solved and of the techniques employed, we present a series of testbeds carried out to validate the code. Such tests range from the evolution of isolated black holes to the head-on collision of two black holes and then to a binary black hole inspiral and merger. Besides assessing the accuracy of the code, the inspiral and merger test has revealed that individual apparent horizons can touch and even intersect. This novel feature in the dynamics of the marginally trapped surfaces is unexpected but consistent with theorems on the properties of apparent horizons.
We analyze the excision strategy for simulating black holes. The problem is modeled by the propagation of quasi-linear waves in a 1-dimensional spatial region with timelike outer boundary, spacelike inner boundary and a horizon in between. Proofs of well-posed evolution and boundary algorithms for a second differential order treatment of the system are given for the separate pieces underlying the finite difference problem. These are implemented in a numerical code which gives accurate long term simulations of the quasi-linear excision problem. Excitation of long wavelength exponential modes, which are latent in the problem, are suppressed using conservation laws for the discretized system. The techniques are designed to apply directly to recent codes for the Einstein equations based upon the harmonic formulation.
Binary black hole interactions provide potentially the strongest source of gravitational radiation for detectors currently under development. We present some results from the Binary Black Hole Grand Challenge Alliance three- dimensional Cauchy evolution module. These constitute essential steps towards modeling such interactions and predicting gravitational radiation waveforms. We report on single black hole evolutions and the first successful demonstration of a black hole moving freely through a three-dimensional computational grid via a Cauchy evolution: a hole moving ~6M at 0.1c during a total evolution of duration ~60M.
We present, in an explicit form, the metric for all spherically symmetric Schwarzschild-Bach black holes in Einstein-Weyl theory. In addition to the black hole mass, this complete family of spacetimes involves a parameter that encodes the value of the Bach tensor on the horizon. When this additional non-Schwarzschild parameter is set to zero the Bach tensor vanishes everywhere and the Schwa-Bach solution reduces to the standard Schwarzschild metric of general relativity. Compared with previous studies, which were mainly based on numerical integration of a complicated form of field equations, the new form of the metric enables us to easily investigate geometrical and physical properties of these black holes, such as specific tidal effects on test particles, caused by the presence of the Bach tensor, as well as fundamental thermodynamical quantities.
We consider the black hole information problem in an explicitly defined spacetime modelling black hole evaporation. Using this context we review basic aspects of the problem, with a particular effort to be unambiguous about subtle topics, for instance precisely what is meant by entropy in various circumstances. We then focus on questions of unitarity, and argue that commonly invoked semiclassical statements of long term, evaporation time, and Page time unitarity may all be violated even if physics is fundamentally unitary. This suggests that there is no horizon firewall. We discuss how the picture is modified for regular (nonsingular) evaporation models. We also compare our conclusions to recent holographic studies, and argue that they are mutually compatible.
We analyze the effect of gravitational radiation reaction on generic orbits around a body with an axisymmetric mass quadrupole moment Q to linear order in Q, to the leading post-Newtonian order, and to linear order in the mass ratio. This system admits three constants of the motion in absence of radiation reaction: energy, angular momentum, and a third constant analogous to the Carter constant. We compute instantaneous and time-averaged rates of change of these three constants. For a point particle orbiting a black hole, Ryan has computed the leading order evolution of the orbits Carter constant, which is linear in the spin. Our result, when combined with an interaction quadratic in the spin (the coupling of the black holes spin to its own radiation reaction field), gives the next to leading order evolution. The effect of the quadrupole, like that of the linear spin term, is to circularize eccentric orbits and to drive the orbital plane towards antialignment with the symmetry axis. In addition we consider a system of two point masses where one body has a single mass multipole or current multipole. To linear order in the mass ratio, to linear order in the multipole, and to the leading post-Newtonian order, we show that there does not exist an analog of the Carter constant for such a system (except for the cases of spin and mass quadrupole). With mild additional assumptions, this result falsifies the conjecture that all vacuum, axisymmetric spacetimes posess a third constant of geodesic motion.