Ultimate Black Hole Recoil: What is the Maximum High-Energy Collision Kick?

We performed a series of 1381 full numerical simulations of high energy collision of black holes to search for the maximum recoil velocity after their merger. We consider equal mass binaries with opposite spins pointing along their orbital plane and perform a search of spin orientations, impact parameters, and initial linear momenta to find the maximum recoil for a given spin magnitude s. This spin sequence for s=0.4, 0.7, 0.8, 0.85, 0.9 is then extrapolated to the extreme case, s=1, to obtain an estimated maximum recoil velocity of 28,562±342 km/s, thus approximately bounded by 10% of the speed of light.

Exploring the Small Mass Ratio Binary Black Hole Merger via Zeno's Dichotomy Approach.

We perform a sequence of binary black hole simulations with increasingly small mass ratios, reaching to a 128:1 binary that displays 13 orbits before merger. Based on a detailed convergence study of the q=m_{1}/m_{2}=1/15 nonspinning case, we apply additional mesh refinement levels around the smaller hole horizon [30] to reach successively the q=1/32, q=1/64, and q=1/128 cases. Roughly a linear computational resources scaling with 1/q is observed on eight-nodes simulations. We compute the remnant properties of the merger: final mass, spin, and recoil velocity, finding precise consistency between horizon and radiation measures. We also compute the gravitational waveforms: their peak frequency, amplitude, and luminosity. We compare those values with predictions of the corresponding phenomenological formulas, reproducing the particle limit within 2%, and we then use the new results to improve their fitting coefficients.

Flip-flopping binary black holes.

We study binary spinning black holes to display the long term individual spin dynamics. We perform a full numerical simulation starting at an initial proper separation of d≈25M between equal mass holes and evolve them down to merger for nearly 48 orbits, 3 precession cycles, and half of a flip-flop cycle. The simulation lasts for t=20 000M and displays a total change in the orientation of the spin of one of the black holes from an initial alignment with the orbital angular momentum to a complete antialignment after half of a flip-flop cycle. We compare this evolution with an integration of the 3.5 post-Newtonian equations of motion and spin evolution to show that this process continuously flip flops the spin during the lifetime of the binary until merger. We also provide lower order analytic expressions for the maximum flip-flop angle and frequency. We discuss the effects this dynamics may have on spin growth in accreting binaries and on the observational consequences for galactic and supermassive binary black holes.

Hangup kicks: still larger recoils by partial spin-orbit alignment of black-hole binaries.

We revisit the scenario of the gravitational radiation recoil acquired by the final remnant of a black-hole-binary merger by studying a set of configurations that have components of the spin both aligned with the orbital angular momentum and in the orbital plane. We perform a series of 42 new full numerical simulations for equal-mass and equal-spin-magnitude binaries. We extend previous recoil fitting formulas to include nonlinear terms in the spins and successfully include both the new and known results. The new predicted maximum velocity approaches 5000 km/s for spins partially aligned with the orbital angular momentum, which leads to an important increase of the probabilities of large recoils in generic astrophysical mergers. We find non-negligible probabilities for recoils of several thousand km/s from accretion-aligned binaries.

Orbital evolution of extreme-mass-ratio black-hole binaries with numerical relativity.

We perform the first fully nonlinear numerical simulations of black-hole binaries with mass ratios 100∶1. Our technique is based on the moving puncture formalism with a new gauge condition and an optimal choice of the mesh refinement. The evolutions start with a small nonspinning black hole just outside the ISCO that orbits twice before plunging. We compute the gravitational radiation, as well as the final remnant parameters, and find close agreement with perturbative estimates. We briefly discuss the relevance of these simulations for Advanced LIGO, third-generation ground-based detectors, LISA observations, and self-force computations.

Intermediate-mass-ratio black-hole binaries: numerical relativity meets perturbation theory.

We study black-hole binaries in the intermediate-mass-ratio regime 0.01≲q≲0.1 with a new technique that makes use of nonlinear numerical trajectories and efficient perturbative evolutions to compute waveforms at large radii for the leading and nonleading (ℓ, m) modes. As a proof-of-concept, we compute waveforms for q=1/10. We discuss applications of these techniques for LIGO and VIRGO data analysis and the possibility that our technique can be extended to produce accurate waveform templates from a modest number of fully nonlinear numerical simulations.

Maximum gravitational recoil.

Recent calculations of gravitational radiation recoil generated during black-hole binary mergers have reopened the possibility that a merged binary can be ejected even from the nucleus of a massive host galaxy. Here we report the first systematic study of gravitational recoil of equal-mass binaries with equal, but counteraligned, spins parallel to the orbital plane. Such an orientation of the spins is expected to maximize the recoil. We find that recoil velocity (which is perpendicular to the orbital plane) varies sinusoidally with the angle that the initial spin directions make with the initial linear momenta of each hole and scales up to a maximum of approximately 4000 km s-1 for maximally rotating holes. Our results show that the amplitude of the recoil velocity can depend sensitively on spin orientations of the black holes prior to merger.