Research / MBH Recoil
Javiera Guedes, Piero Madau, Michael Kuhlen, Jurg Diemand, Marcel Zemp, Lucio Mayer and Simone Callegari.
Numerical relativistic codes are now able follow the orbits of black holes from inspiral to merger. Back in the 60’s, Peres suggested that asymmetries in the configuration black hole binaries (different masses or spins) would have the net effect of beaming gravitational wave radiation in some preferred direction, a problem theoretically analog to electromagnetic radiation recoil. In the gravitational case the recoil arises from the interference of the mass quadrupole and the mass octopole (or alternatively, the flow quadrupole radiation.) The magnitude of the kick velocity depends on: a) the mass ratio of the binary, b) the spin magnitudes and c) spin orientations. The Baker et al. 2008 numerical fit looks like this:
![]() |
![]() |
a) In spherical bulges the black hole will oscillate right through the center, and return to the center after a few dynamical times, or deplete or "heat up" the stars in the nucleus, reducing the efficiency of dynamical and extending the return time.
b) If the black hole is able to leave the bulge and enters a triaxial dark matter halo, it may not return at all within a Hubble time (we recently published a paper that looks at this in detail.)
c) High resolution simulations of gas mergers (e.g. Mayer et al. 2007) show that up to 60% of the initial gas is funneled to the center of the remnant. In that particular simulation, the mass of the nuclear disk is Md = 3e9 Msun, and recoling black holes with velocities of ~500 km/s are not able to escape the central region (scale length L=75 pc).
A movie of a black hole recoiling with initial kick velocity V=1200 km/s in the direction perpendicular to the nuclear disk is shown below. Click for movie.
![]() |