We conduct a descriptive analysis of the multipolar structure of gravitational-radiation waveforms from equal-mass aligned-spin mergers, following an approach first presented in the complementary context of nonspinning black holes of varying mass ratio [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)]. We find that, as with the nonspinning mergers, the dominant waveform mode phases evolve together in lock-step through inspiral and merger, supporting the previous waveform description in terms of an adiabatically rigid rotator driving gravitational-wave emission-an implicit rotating source. We further apply the late-time merger-ringdown model for the rotational frequency introduced in [J.G. Baker et al. Phys. Rev. D 78 044046 (2008)], along with an improved amplitude model appropriate for the dominant (2, \pm2) modes. This provides a quantitative description of the merger-ringdown waveforms, and suggests that the major features of these waveforms can be described with reference only to the intrinsic parameters associated with the state of the final black hole formed in the merger. We provide an explicit model for the merger-ringdown radiation, and demonstrate that this model agrees to fitting factors better than 95% with the original numerical waveforms for system masses above \sim150M\odot. This model may be directly applicable to gravitational-wave detection of intermediate-mass black-hole mergers.

Although the gravitational wave kick velocity in the orbital plane of coalescing black holes has been understood for some time, apparently conflicting formulae have been proposed for the dominant out-of-plane kick, each a good fit to different data sets. This is important to resolve because it is only the out-of-plane kicks that can reach more than 500 km/s and can thus eject merged remnants from galaxies. Using a different ansatz for the out-of-plane kick, we show that we can fit almost all existing data to better than 5 %. This is good enough for any astrophysical calculation, and shows that the previous apparent conflict was only because the two data sets explored different aspects of the kick parameter space.

Coalescing massive black hole binaries are produced by the mergers of galaxies. The final stages of the black hole coalescence produce strong gravitational radiation that can be detected by the space-borne LISA. In cases where the black hole merger takes place in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. Modeling such electromagnetic counterparts of the final merger requires evolving the behavior of both gas and fields in the strong-field regions around the black holes. We have taken a step towards solving this problem by mapping the flow of pressureless matter in the dynamic, 3-D general relativistic spacetime around the merging black holes. We find qualitative differences in collision and outflow speeds, including a signature of the merger when the net angular momentum of the matter is low, between the results from single and binary black holes, and between nonrotating and rotating holes in binaries. If future magnetohydrodynamic results confirm these differences, it may allow assessment of the properties of the binaries as well as yielding an identifiable electromagnetic counterpart to the attendant gravitational wave signal.

Recent numerical relativistic results demonstrate that the merger of comparable-mass spinning black holes has a maximum ``recoil kick'' of up to $\sim 4000 \kms$. However the scaling of these recoil velocities with mass ratio is poorly understood. We present new runs showing that the maximum possible kick perpendicular to the orbital plane does not scale as $\sim\eta^2$ (where $\eta$ is the symmetric mass ratio), as previously proposed, but is more consistent with $\sim\eta^3$, at least for systems with low orbital precession. We discuss the effect of this dependence on galactic ejection scenarios and retention of intermediate-mass black holes in globular clusters.

We present a multipolar analysis of the gravitational recoil computed in recent numerical simulations of binary black hole (BH) coalescence, for both unequal masses and non-zero, non-precessing spins. We show that multipole moments up to and including l=4 are sufficient to accurately reproduce the final recoil velocity (within ~2%) and that only a few dominant modes contribute significantly to it (within ~5%). We describe how the relative amplitudes, and more importantly, the relative phases, of these few modes control the way in which the recoil builds up throughout the inspiral, merger, and ringdown phases. We also find that the numerical results can be reproduced by an ``effective Newtonian'' formula for the multipole moments obtained by replacing the radial separation in the Newtonian formulae with an effective radius computed from the numerical data. Beyond the merger, the numerical results are reproduced by a superposition of three Kerr quasi-normal modes (QNMs). Analytic formulae, obtained by expressing the multipole moments in terms of the fundamental QNMs of a Kerr BH, are able to explain the onset and amount of ``anti-kick'' for each of the simulations. Lastly, we apply this multipolar analysis to help explain the remarkable difference between the amplitudes of planar and non-planar kicks for equal-mass spinning black holes.

Several groups have recently computed the gravitational radiation recoil produced by the merger of two spinning black holes. The results suggest that spin can be the dominant contributor to the kick, with reported recoil speeds of hundreds to even thousands of kilometers per second. The parameter space of spin kicks is large, however, and it is ultimately desirable to have a simple formula that gives the approximate magnitude of the kick given a mass ratio, spin magnitudes, and spin orientations. As a step toward this goal, we perform a systematic study of the recoil speeds from mergers of black holes with mass ratio $q\equiv m_1/m_2=2/3$ and dimensionless spin parameters of $a_1/m_1$ and $a_2/m_2$ equal to 0 or 0.2, with directions aligned or anti-aligned with the orbital angular momentum. We also run an equal-mass $a_1/m_1=-a_2/m_2=0.2$ case, and find good agreement with previous results. We find that, for currently reported kicks from aligned or anti-aligned spins, a simple kick formula inspired by post-Newtonian analyses can reproduce the numerical results to better than $\sim$10%.

Coalescing black-hole binaries are expected to be the strongest sources of gravitational waves for ground-based interferometers as well as the space-based interferometer LISA. Recent progress in numerical relativity now makes it possible to calculate the waveforms from the strong-field dynamical merger and is revolutionizing our understanding of these systems. We review these dramatic developments, emphasizing applications to issues in gravitational wave observations. These new capabilities also make possible accurate calculations of the recoil or kick imparted to the final remnant black hole when the merging components have unequal masses, or unequal or unaligned spins. We highlight recent work in this area, focusing on results of interest to astrophysics.