William Krivan

Dynamics of perturbations of rotating black holes

We present a numerical study of the time evolution of perturbations of rotating black holes. The solutions are obtained by integrating the Teukolsky equation written as a first order in time, coupled system of equations. We address the numerical difficulties of solving the equation in its original form. We follow the propagation of generic initial data through the burst, quasinormal ringing and power-law tail phases. In particular, we calculate the effects due to the rotation of the black hole on the scattering of incident gravitational wave pulses. These effects include the amplitude enhancement due to so-called super-radiance. The results may help explain how the angular momentum of the black hole affects the gravitational waves that are generated during the final stages of black hole coalescence.

Initial data for superposed rotating black holes

The standard approach to initial data for both analytic and numerical computations of black hole collisions has been to use conformally-flat initial geometry. Among other advantages, this choice allows the simple superposition of holes with arbitrary mass, location and spin. The conformally flat restriction, however, is inappropriate to the study of Kerr holes, for which the standard constant-time slice is not conformally flat. Here we point out that for axisymmetric arrangements of rotating holes, a nonconformally flat form of the 3-geometry can be chosen which allows fairly simple superposition of Kerr holes with arbitrary mass and spin. We present initial data solutions representing locally Kerr holes at large separation, and representing rotating holes close enough so that outside a common horizon the spacetime geometry is a perturbation of a single Kerr hole.

Dynamics of scalar fields in the background of rotating black holes

A numerical study of the evolution of a massless scalar field in the background of rotating black holes is presented. First, solutions to the wave equation are obtained for slowly rotating black holes. In this approximation, the background geometry is treated as a perturbed Schwarzschild spacetime with the angular momentum per unit mass playing the role of a perturbative parameter. To first order in the angular momentum of the black hole, the scalar wave equation yields two coupled one-dimensional evolution equations for a function representing the scalar field in the Schwarzschild background and a second field that accounts for the rotation. Solutions to the wave equation are also obtained for rapidly rotating black holes. In this case, the wave equation does not admit complete separation of variables and yields a two-dimensional evolution equation. The study shows that, for rotating black holes, the late time dynamics of a massless scalar field exhibit the same power-law behavior as in the case of a Schwarzschild background independently of the angular momentum of the black hole.

Quasi-stationary binary inspiral: II. Radiation-balanced boundary conditions

The quasi-stationary method for black hole binary inspiral is an approximation for studying strong-field effects while suppressing radiation reaction. In this paper we use a nonlinear scalar field toy model to: (a) explain the underlying method of approximating binary motion by periodic orbits with radiation; (b) show how the fields in such a model are found by the solution of a boundary value problem; (c) demonstrate how a good approximation to the outgoing radiation can be found by finding fields with a balance of ingoing and outgoing radiation (a generalization of standing waves).

Late-time dynamics of scalar fields on rotating black hole backgrounds

Motivated by results of recent analytic studies, we present a numerical investigation of the late-time dynamics of scalar test fields on Kerr backgrounds. We pay particular attention to the issue of mixing of different multipoles and their fall-off behavior at late times. Confining ourselves to the special case of axisymmetric modes with equatorial symmetry, we show that, in agreement with the results of previous work, the late-time behavior is dominated by the lowest allowed l-multipole. However the numerical results imply that, in general, the late-time falloff of the dominating multipole is different from that in the Schwarzschild case, and seems to be incompatible with a result of a recently published analytic study.

The Imposition of Cauchy data to the Teukolsky equation. 2. Numerical comparison with the Zerilli-Moncrief approach to black hole perturbations

We revisit the question of the imposition of initial data representing astrophysical gravitational perturbations of black holes. We study their dynamics for the case of nonrotating black holes by numerically evolving the Teukolsky equation in the time domain. In order to express the Teukolsky function Psi explicitly in terms of hypersurface quantities, we relate it to the Moncrief waveform phi_M through a Chandrasekhar transformation in the case of a nonrotating black hole. This relation between Psi and phi_M holds for any constant time hypersurface and allows us to compare the computation of the evolution of Schwarzschild perturbations by the Teukolsky and by the Zerilli and Regge-Wheeler equations. We explicitly perform this comparison for the Misner initial data in the close limit approach. We evolve numerically both, the Teukolsky (with the recent code of Ref. [1]) and the Zerilli equations, finding complete agreement in resulting waveforms within numerical error. The consistency of these results further supports the correctness of the numerical code for evolving the Teukolsky equation as well as the analytic expressions for Psi in terms only of the three-metric and the extrinsic curvature.