Perturbations and Stability of Black Ellipsoids

We study the perturbations of two classes of static black ellipsoid solutions of four dimensional vacuum Einstein equations. Such solutions are described by generic off--diagonal metrics which are generated by anholonomic transforms of diagonal metrics. The analysis is performed in the approximation of small eccentricity deformations of the Schwarzschild solution. We conclude that such anisotropic black hole objects may be stable with respect to the perturbations parametrized by the Schrodinger equations in the framework of the one--dimensional inverse scattering theory.Comment: Published variant in IJMD with small modifications in formulas and new reference

Understanding initial data for black hole collisions

Numerical relativity, applied to collisions of black holes, starts with initial data for black holes already in each other's strong field. The initial hypersurface data typically used for computation is based on mathematical simplifying prescriptions, such as conformal flatness of the 3-geometry and longitudinality of the extrinsic curvature. In the case of head on collisions of equal mass holes, there is evidence that such prescriptions work reasonably well, but it is not clear why, or whether this success is more generally valid. Here we study these questions by considering the ``particle limit'' for head on collisions of nonspinning holes. Einstein's equations are linearized in the mass of the small hole, and described by a single gauge invariant spacetime function psi, for each multipole. The resulting equations have been solved by numerical evolution for collisions starting from various initial separations, and the evolution is studied on a sequence of hypersurfaces. In particular, we extract hypersurface data, that is psi and its time derivative, on surfaces of constant background Schwarzschild time. These evolved data can then be compared with ``prescribed'' data, evolved data can be replaced by prescribed data on any hypersurface, and evolved further forward in time, a gauge invariant measure of deviation from conformal flatness can be evaluated, etc. The main findings of this study are: (i) For holes of unequal mass the use of prescribed data on late hypersurfaces is not successful. (ii) The failure is likely due to the inability of the prescribed data to represent the near field of the smaller hole. (iii) The discrepancy in the extrinsic curvature is more important than in the 3-geometry. (iv) The use of the more general conformally flat longitudinal data does not notably improve this picture.Comment: 20 pages, REVTEX, 26 PS figures include

Evolving the Bowen-York initial data for spinning black holes

The Bowen-York initial value data typically used in numerical relativity to represent spinning black hole are not those of a constant-time slice of the Kerr spacetime. If Bowen-York initial data are used for each black hole in a collision, the emitted radiation will be partially due to the ``relaxation'' of the individual holes to Kerr form. We compute this radiation by treating the geometry for a single hole as a perturbation of a Schwarzschild black hole, and by using second order perturbation theory. We discuss the extent to which Bowen-York data can be expected accurately to represent Kerr holes.Comment: 10 pages, RevTeX, 4 figures included with psfi

Cold Plasma Dispersion Relations in the Vicinity of a Schwarzschild Black Hole Horizon

We apply the ADM 3+1 formalism to derive the general relativistic magnetohydrodynamic equations for cold plasma in spatially flat Schwarzschild metric. Respective perturbed equations are linearized for non-magnetized and magnetized plasmas both in non-rotating and rotating backgrounds. These are then Fourier analyzed and the corresponding dispersion relations are obtained. These relations are discussed for the existence of waves with positive angular frequency in the region near the horizon. Our results support the fact that no information can be extracted from the Schwarzschild black hole. It is concluded that negative phase velocity propagates in the rotating background whether the black hole is rotating or non-rotating.Comment: 27 pages, 11 figures accepted for publication in Gen. Relat. & Gravi

Gravitational waveforms from a point particle orbiting a Schwarzschild black hole

We numerically solve the inhomogeneous Zerilli-Moncrief and Regge-Wheeler equations in the time domain. We obtain the gravitational waveforms produced by a point-particle of mass $\mu$ traveling around a Schwarzschild black hole of mass M on arbitrary bound and unbound orbits. Fluxes of energy and angular momentum at infinity and the event horizon are also calculated. Results for circular orbits, selected cases of eccentric orbits, and parabolic orbits are presented. The numerical results from the time-domain code indicate that, for all three types of orbital motion, black hole absorption contributes less than 1% of the total flux, so long as the orbital radius r_p(t) satisfies r_p(t)> 5M at all times.Comment: revtex4, 24 pages, 23 figures, 3 tables, submitted to PR

Fourth order indirect integration method for black hole perturbations: even modes

On the basis of a recently proposed strategy of finite element integration in time domain for partial differential equations with a singular source term, we present a fourth order algorithm for non-rotating black hole perturbations in the Regge-Wheeler gauge. Herein, we address even perturbations induced by a particle plunging in. The forward time value at the upper node of the $(r^*,t)$ grid cell is obtained by an algebraic sum of i) the preceding node values of the same cell, ii) analytic expressions, related to the jump conditions on the wave function and its derivatives, iii) the values of the wave function at adjacent cells. In this approach, the numerical integration does not deal with the source and potential terms directly, for cells crossed by the particle world line. This scheme has also been applied to circular and eccentric orbits and it will be object of a forthcoming publication.Comment: This series of papers deals with EMRI for LISA. With the respect to the v1 version, the algorithm has been improved; convergence tests and references have been added; v2 is composed by 23 pages, and 6 figures. Paper accepted by Class. Quantum Gravity for the special issue on Theory Meets Data Analysis at Comparable and Extreme Mass Ratios (Capra and NRDA) at Perimeier Institute in June 201

Axial perturbations of general spherically symmetric spacetimes

The aim of this paper is to present a governing equation for first order axial metric perturbations of general, not necessarily static, spherically symmetric spacetimes. Under the non-restrictive assumption of axisymmetric perturbations, the governing equation is shown to be a two-dimensional wave equation where the wave function serves as a twist potential for the axisymmetry generating Killing vector. This wave equation can be written in a form which is formally a very simple generalization of the Regge-Wheeler equation governing the axial perturbations of a Schwarzschild black hole, but in general the equation is accompanied by a source term related to matter perturbations. The case of a viscous fluid is studied in particular detail.Comment: 16 pages, no figures, minor correction

Head-on collisions of black holes: the particle limit

We compute gravitational radiation waveforms, spectra and energies for a point particle of mass $m_0$ falling from rest at radius $r_0$ into a Schwarzschild hole of mass $M$. This radiation is found to lowest order in $(m_0/M)$ with the use of a Laplace transform. In contrast with numerical relativity results for head-on collisions of equal-mass holes, the radiated energy is found not to be a monotonically increasing function of initial separation; there is a local radiated-energy maximum at $r_0\approx4.5M$. The present results, along with results for infall from infinity, provide a complete catalog of waveforms and spectra for particle infall. We give a representative sample from that catalog and an interesting observation: Unlike the simple spectra for other head-on collisions (either of particle and hole, or of equal mass holes) the spectra for $\infty>r_0>\sim5M$ show a series of evenly spaced bumps. A simple explanation is given for this. Lastly, our energy vs. $r_0$ results are compared with approximation methods used elsewhere, for small and for large initial separation.Comment: 15 pages, REVTeX, 25 figure

Shell sources as a probe of relativistic effects in neutron star models

A perturbing shell is introduced as a device for studying the excitation of fluid motions in relativistic stellar models. We show that this approach allows a reasonably clean separation of radiation from the shell and from fluid motions in the star, and provides broad flexibility in the location and timescale of perturbations driving the fluid motions. With this model we compare the relativistic and Newtonian results for the generation of even parity gravitational waves from constant density models. Our results suggest that relativistic effects will not be important in computations of the gravitational emission except possibly in the case of excitation of the neutron star on very short time scales.Comment: 16 pages LaTeX with 6 eps figures; submitted to Phys. Rev.

Energy and angular momentum of the weak gravitational waves on the Schwarzschild background -- quasilocal gauge-invariant formulation

It is shown that the axial and polar perturbations of the spherically symmetric black hole can be described in a gauge-invariant way. The reduced phase space describing gravitational waves outside of the horizon is described by the gauge-invariant quantities. Both degrees of freedom fulfill generalized scalar wave equation. For the axial degree of freedom the radial part of the equation corresponds to the Regge-Wheeler result (Phys. Rev. 108, 1063-1069 (1957)) and for the polar one we get Zerilli result (Phys. Rev. D2, 2141-2160 (1970)), see also Chandrasekhar (The Mathematical Theory of Black Holes,(Clarendon Press Oxford, 1983)), Moncrief (Annals of Physics 88, 323-342 (1974)) for both. An important ingredient of the analysis is the concept of quasilocality which does duty for the separation of the angular variables in the usual approach. Moreover, there is no need to represent perturbations by normal modes (with time dependence $\exp(-ikt)$), we have fields in spacetime and the Cauchy problem for them is well defined outside of the horizon. The reduced symplectic structure explains the origin of the axial and polar invariants. It allows to introduce an energy and angular momentum for the gravitational waves which is invariant with respect to the gauge transformations. Both generators represent quadratic approximation of the ADM nonlinear formulae in terms of the perturbations of the Schwarzschild metric. We also discuss the boundary-initial value problem for the linearized Einstein equations on a Schwarzschild background outside of the horizon.Comment: 23 page