Remnant mass, spin, and recoil from spin aligned black-hole binaries
  We perform a set of 36 nonprecessing black-hole binary simulations with spins
either aligned or counteraligned with the orbital angular momentum in order to
model the final mass, spin, and recoil of the merged black hole as a function
of the individual black hole spin magnitudes and the mass ratio of the
progenitors. We find that the maximum recoil for these configurations is
$V_{max}=526\pm23\,km/s$, which occurs when the progenitor spins are maximal,
the mass ratio is $q_{max}=m_1/m_2=0.623\pm0.038$, the smaller black-hole spin
is aligned with the orbital angular momentum, and the larger black-hole spin is
counteraligned ($\alpha_1=-\alpha_2=1$). This maximum recoil is about
$80\,km/s$ larger than previous estimates, but most importantly, because the
maximum occurs for smaller mass ratios, the probability for a merging binary to
recoil faster than $400\,km/s$ can be as large as $17\%$, while the probability
for recoils faster than $250\, km/s$ can be as large as $45\%$. We provide
explicit phenomenological formulas for the final mass, spin, and recoil as a
function of the individual BH spins and the mass difference between the two
black holes. Here we include terms up through fourth-order in the initial spins
and mass difference, and find excellent agreement (within a few percent) with
independent results available in the literature. The maximum radiated energy is
$E_{\rm rad}/m\approx11.3\%$ and final spin $\alpha_{\rm rem}^{\rm
max}\approx0.952$ for equal mass, aligned maximally spinning binaries.
