Scalar Deformations of Schwarzschild Holes and Their Stability

We construct two solutions of the minimally coupled Einstein-scalar field equations, representing regular deformations of Schwarzschild black holes by a self-interacting, static, scalar field. One solution features an exponentially decaying scalar field and a triple-well interaction potential; the other one is completely analytic and sprouts Coulomb-like scalar hair. Both evade the no-hair theorem by having partially negative potential, in conflict with the dominant energy condition. The linear perturbation theory around such backgrounds is developed in general, and yields stability criteria in terms of effective potentials for an analog Schr\"odinger problem. We can test for more than half of the perturbation modes, and our solutions prove to be stable against those.Comment: 24 pp, 16 figs, Latex; version published in Int. J. Mod. Phys.

New Critical Behavior in Einstein-Yang-Mills Collapse

We extend the investigation of the gravitational collapse of a spherically symmetric Yang-Mills field in Einstein gravity and show that, within the black hole regime, a new kind of critical behavior arises which separates black holes formed via Type I collapse from black holes formed through Type II collapse. Further, we provide evidence that these new attracting critical solutions are in fact the previously discovered colored black holes with a single unstable mode.Comment: 13 pages, 4 figure

Black Hole--Scalar Field Interactions in Spherical Symmetry

We examine the interactions of a black hole with a massless scalar field using a coordinate system which extends ingoing Eddington-Finkelstein coordinates to dynamic spherically symmetric-spacetimes. We avoid problems with the singularity by excising the region of the black hole interior to the apparent horizon. We use a second-order finite difference scheme to solve the equations. The resulting program is stable and convergent and will run forever without problems. We are able to observe quasi-normal ringing and power-law tails as well an interesting nonlinear feature.Comment: 16 pages, 26 figures, RevTex, to appear in Phys. Rev.

Late-time evolution of nonlinear gravitational collapse

We study numerically the fully nonlinear gravitational collapse of a self-gravitating, minimally-coupled, massless scalar field in spherical symmetry. Our numerical code is based on double-null coordinates and on free evolution of the metric functions: The evolution equations are integrated numerically, whereas the constraint equations are only monitored. The numerical code is stable (unlike recent claims) and second-order accurate. We use this code to study the late-time asymptotic behavior at fixed $r$ (outside the black hole), along the event horizon, and along future null infinity. In all three asymptotic regions we find that, after the decay of the quasi-normal modes, the perturbations are dominated by inverse power-law tails. The corresponding power indices agree with the integer values predicted by linearized theory. We also study the case of a charged black hole nonlinearly perturbed by a (neutral) self-gravitating scalar field, and find the same type of behavior---i.e., quasi-normal modes followed by inverse power-law tails, with the same indices as in the uncharged case.Comment: 14 pages, standard LaTeX, 18 Encapsulated PostScript figures. A new convergence test and a determination of QN ringing were added, in addition to correction of typos and update of reference

On free evolution of self gravitating, spherically symmetric waves

We perform a numerical free evolution of a selfgravitating, spherically symmetric scalar field satisfying the wave equation. The evolution equations can be written in a very simple form and are symmetric hyperbolic in Eddington-Finkelstein coordinates. The simplicity of the system allow to display and deal with the typical gauge instability present in these coordinates. The numerical evolution is performed with a standard method of lines fourth order in space and time. The time algorithm is Runge-Kutta while the space discrete derivative is symmetric (non-dissipative). The constraints are preserved under evolution (within numerical errors) and we are able to reproduce several known results.Comment: 15 pages, 15 figure

Boosted three-dimensional black-hole evolutions with singularity excision

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

Black Hole Data via a Kerr-Schild Approach

We present a new approach for setting initial Cauchy data for multiple black hole spacetimes. The method is based upon adopting an initially Kerr-Schild form of the metric. In the case of non-spinning holes, the constraint equations take a simple hierarchical form which is amenable to direct numerical integration. The feasibility of this approach is demonstrated by solving analytically the problem of initial data in a perturbed Schwarzschild geometry.Comment: 13 pages, RevTeX forma

Bondian frames to couple matter with radiation

A study is presented for the non linear evolution of a self gravitating distribution of matter coupled to a massless scalar field. The characteristic formulation for numerical relativity is used to follow the evolution by a sequence of light cones open to the future. Bondian frames are used to endow physical meaning to the matter variables and to the massless scalar field. Asymptotic approaches to the origin and to infinity are achieved; at the boundary surface interior and exterior solutions are matched guaranteeing the Darmois--Lichnerowicz conditions. To show how the scheme works some numerical models are discussed. We exemplify evolving scalar waves on the following fixed backgrounds: A) an atmosphere between the boundary surface of an incompressible mixtured fluid and infinity; B) a polytropic distribution matched to a Schwarzschild exterior; C) a Schwarzschild- Schwarzschild spacetime. The conservation of energy, the Newman--Penrose constant preservation and other expected features are observed.Comment: 20 pages, 6 figures; to appear in General Relativity and Gravitatio

Numerical Evolution of Black Holes with a Hyperbolic Formulation of General Relativity

We describe a numerical code that solves Einstein's equations for a Schwarzschild black hole in spherical symmetry, using a hyperbolic formulation introduced by Choquet-Bruhat and York. This is the first time this formulation has been used to evolve a numerical spacetime containing a black hole. We excise the hole from the computational grid in order to avoid the central singularity. We describe in detail a causal differencing method that should allow one to stably evolve a hyperbolic system of equations in three spatial dimensions with an arbitrary shift vector, to second-order accuracy in both space and time. We demonstrate the success of this method in the spherically symmetric case.Comment: 23 pages RevTeX plus 7 PostScript figures. Submitted to Phys. Rev.