
In a regime where data are only mildly informative, prior choices can play a
significant role in Bayesian statistical inference, potentially affecting the
inferred physics. We show this is indeed the case for some of the parameters
inferred from current gravitationalwave measurements of binary black hole
coalescences. We reanalyze the first detections performed by the twin LIGO
interferometers using alternative (and astrophysically motivated) prior
assumptions. We find different prior distributions can introduce deviations in
the resulting posteriors that impact the physical interpretation of these
systems. For instance, (i) limits on the $90\%$ credible interval on the
effective black hole spin $\chi_{\rm eff}$ are subject to variations of $\sim
10\%$ if a prior with black hole spins mostly aligned to the binary's angular
momentum is considered instead of the standard choice of isotropic spin
directions, and (ii) under priors motivated by the initial stellar mass
function, we infer tighter constraints on the black hole masses, and in
particular, we find no support for any of the inferred masses within the
putative mass gap $M \lesssim 5 M_\odot$.

Gravitational waves emitted by coalescing compact objects carry information
about the spin of the individual bodies. However, with present detectors only
the massweighted combination of the components of the spin along the orbital
angular momentum can be measured accurately. This quantity, the effective spin
$\chi_{\mathrm{eff}}$, is conserved up to at least the second postNewtonian
order. The measured distribution of $\chi_{\mathrm{eff}}$ values from a
population of detected binaries, and in particular whether this distribution is
symmetric about zero, encodes valuable information about the underlying
compactbinary formation channels. In this paper we focus on two important
complications of using the effective spin to study astrophysical population
properties: (i) an astrophysical distribution for $\chi_{\mathrm{eff}}$ values
which is symmetric does not necessarily lead to a symmetric distribution for
the detected effective spin values, leading to a \emph{selection bias}; and
(ii) the posterior distribution of $\chi_{\mathrm{eff}}$ for individual events
is \emph{asymmetric} and it cannot usually be treated as a Gaussian. We find
that the posterior distributions for $\chi_{\mathrm{eff}}$ systematically show
fatter tails toward larger positive values, unless the total mass is large or
the mass ratio $m_2/m_1$ is smaller than $\sim 1/2$. Finally we show that
uncertainties in the measurement of $\chi_{\mathrm{eff}}$ are systematically
larger when the true value is negative than when it is positive. All these
factors can bias astrophysical inference about the population when we have more
than $\sim 100$ events and should be taken into account when using
gravitationalwave measurements to characterize astrophysical populations.

Binary black holes radiate linear momentum in gravitational waves as they
merge. Recoils imparted to the blackhole remnant can reach thousands of km/s,
thus ejecting black holes from their host galaxies. We exploit recent advances
in gravitational waveform modeling to quickly and reliably extract recoils
imparted to generic, precessing, black hole binaries. Our procedure uses a
numericalrelativity surrogate model to obtain the gravitational waveform given
a set of binary parameters, then from this waveform we directly integrate the
gravitationalwave linear momentum flux. This entirely bypasses the need of
fitting formulae which are typically used to model blackhole recoils in
astrophysical contexts. We provide a thorough exploration of the blackhole
kick phenomenology in the parameter space, summarizing and extending previous
numerical results on the topic. Our extraction procedure is made publicly
available as a module for the Python programming language named SURRKICK. Kick
evaluations take ~0.1s on a standard offtheshelf machine, thus making our
code ideal to be ported to largescale astrophysical studies.

Measurements of blackhole spins are of crucial importance to fulfill the
promise of gravitationalwave astronomy. On the astrophysics side, spins are
perhaps the cleanest indicator of blackhole evolutionary processes, thus
providing a preferred way to discriminate how LIGO's black holes form. On the
relativity side, spins are responsible for peculiar dynamical phenomena (from
precessional modulations in the long inspiral to gravitationalwave recoils at
merger) which encode precious information on the underlying astrophysical
processes. I present some examples to explore this deep and fascinating
interplay between spin dynamics (relativity) and environmental effects
(astrophysics). Blackhole spins indeed hide remarkable surprises on both
fronts: morphologies, resonances, constraints on supernova kicks, multiple
merger generations and more...

In the postNewtonian (PN) regime, the timescale on which the spins of binary
black holes precess is much shorter than the radiationreaction timescale on
which the black holes inspiral to smaller separations. On the precession
timescale, the angle between the total and orbital angular momenta oscillates
with nutation period $\tau$, during which the orbital angular momentum
precesses about the total angular momentum by an angle $\alpha$. This defines
two distinct frequencies that vary on the radiationreaction timescale: the
nutation frequency $\omega \equiv 2\pi/\tau$ and the precession frequency
$\Omega \equiv \alpha/\tau$. We use analytic solutions for generic spin
precession at 2PN order to derive Fourier series for the total and orbital
angular momenta in which each term is a sinusoid with frequency $\Omega 
n\omega$ for integer $n$. As black holes inspiral, they can pass through
nutational resonances ($\Omega = n\omega$) at which the total angular momentum
tilts. We derive an approximate expression for this tilt angle and show that it
is usually less than $10^{3}$ radians for nutational resonances at binary
separations $r > 10M$. The large tilts occurring during transitional precession
(near zero total angular momentum) are a consequence of such states being
approximate $n=0$ nutational resonances. Our new Fourier series for the total
and orbital angular momenta converge rapidly with $n$ providing an intuitive
and computationally efficient approach to understanding generic precession that
may facilitate future calculations of gravitational waveforms in the PN regime.

This letter considers stellar core collapse in massive scalartensor
theories of gravity. The presence of a mass term for the scalar field allows
for dramatic increases in the radiated gravitational wave signal. There are
several potential smoking gun signatures of a departure from general relativity
associated with this process. These signatures could show up within existing
LIGOVirgo searches.

The inferred parameters of the binary black hole GW151226 are consistent with
nonzero spin for the most massive black hole, misaligned from the binary's
orbital angular momentum. If the black holes formed through isolated binary
evolution from an initially aligned binary star, this misalignment would then
arise from a natal kick imparted to the firstborn black hole at its birth
during stellar collapse. We use simple kinematic arguments to constrain the
characteristic magnitude of this kick, and find that a natal kick $v_k \gtrsim
50$ km/s must be imparted to the black hole at birth to produce misalignments
consistent with GW151226. Such large natal kicks exceed those adopted by
default in most of the current supernova and binary evolution models.

Advanced LIGO detectors at Hanford and Livingston made two confirmed and one
marginal detection of binary black holes during their first observing run. The
first event, GW150914, was from the merger of two black holes much heavier that
those whose masses have been estimated so far, indicating a formation scenario
that might differ from "ordinary" stellar evolution. One possibility is that
these heavy black holes resulted from a previous merger. When the progenitors
of a black hole binary merger result from previous mergers, they should (on
average) merge later, be more massive, and have spin magnitudes clustered
around a dimensionless spin ~0.7. Here we ask the following question: can
gravitationalwave observations determine whether merging black holes were born
from the collapse of massive stars ("first generation"), rather than being the
end product of earlier mergers ("second generation")? We construct simple,
observationally motivated populations of black hole binaries, and we use
Bayesian model selection to show that measurements of the masses, luminosity
distance (or redshift), and "effective spin" of black hole binaries can indeed
distinguish between these different formation scenarios.

We analyze the inspiral dynamics of equalmass precessing blackhole binaries
using multitimescale techniques. The orbitaveraged postNewtonian
evolutionary equations admit two constants of motion in the equalmass limit,
namely the magnitude of the total spin $S$ and the effective spin $\xi$. This
feature makes the entire dynamics qualitatively different compared to the
generic unequalmass case, where only $\xi$ is constant while the variable $S$
parametrizes the precession dynamics. For fixed individual masses and spin
magnitudes, an equalmass blackhole inspiral is uniquely characterized by the
two parameters $(S,\xi)$: these two numbers completely determine the entire
evolution under the effect of radiation reaction. In particular, for equalmass
binaries we find that (i) the blackhole binary spin morphology is constant
throughout the inspiral, and that (ii) the precessional motion of the two
blackhole spins about the total spin takes place on a longer timescale than
the precession of the total spin and the orbital plane about the total angular
momentum.

Generic black hole binaries radiate gravitational waves anisotropically,
imparting a recoil, or kick, velocity to the merger remnant. If a component of
the kick along the line of sight is present, gravitational waves emitted during
the final orbits and merger will be gradually Doppler shifted as the kick
builds up. We develop a simple prescription to capture this effect in existing
waveform models, showing that future gravitational wave experiments will be
able to perform direct measurements, not only of the black hole kick velocity,
but also of its accumulation profile. In particular, the eLISA space mission
will measure supermassive black hole kick velocities as low as ~500 km/s, which
are expected to be a common outcome of black hole binary coalescence following
galaxy mergers. Black hole kicks thus constitute a promising new observable in
the growing field of gravitational wave astronomy.

We present the numerical code PRECESSION: a new opensource python module to
study the dynamics of precessing blackhole binaries in the postNewtonian
regime. The code provides a comprehensive toolbox to (i) study the evolution of
the blackhole spins along their precession cycles, (ii) perform
gravitationalwave driven binary inspirals using both orbitaveraged and
precessionaveraged integrations, and (iii) predict the properties of the
merger remnant through fitting formulae obtained from numericalrelativity
simulations. PRECESSION is a readytouse tool to add the blackhole spin
dynamics to largerscale numerical studies such as gravitationalwave parameter
estimation codes, population synthesis models to predict gravitationalwave
event rates, galaxy merger trees and cosmological simulations of structure
formation. PRECESSION provides fast and reliable integration methods to
propagate statistical samples of blackhole binaries from/to large separations
where they form to/from small separations where they become detectable, thus
linking gravitationalwave observations of spinning blackhole binaries to
their astrophysical formation history. The code is also a useful tool to
compute initial parameters for numericalrelativity simulations targeting
specific precessing systems. PRECESSION can be installed from the Python
Package Index and it is freely distributed under version control on Github,
where further documentation is provided.

We present numericalrelativity simulations of spherically symmetric core
collapse and compactobject formation in scalartensor theories of gravity. The
additional scalar degree of freedom introduces a propagating monopole
gravitationalwave mode. Detection of monopole scalar waves with current and
future gravitationalwave experiments may constitute smoking gun evidence for
strongfield modifications of General Relativity. We collapse both polytropic
and more realistic presupernova profiles using a highresolution
shockcapturing scheme and an approximate prescription for the nuclear equation
of state. The most promising sources of scalar radiation are protoneutron stars
collapsing to black holes. In case of a Galactic core collapse event forming a
black hole, Advanced LIGO may be able to place independent constraints on the
parameters of the theory at a level comparable to current SolarSystem and
binarypulsar measurements. In the region of the parameter space admitting
spontaneously scalarised stars, transition to configurations with prominent
scalar hair before blackhole formation further enhances the emitted signal.
Although a more realistic treatment of the microphysics is necessary to fully
investigate the occurrence of spontaneous scalarisation of neutron star
remnants, we speculate that formation of such objects could constrain the
parameters of the theory beyond the current bounds obtained with SolarSystem
and binarypulsar experiments.

Gravitational waves from coalescing binary black holes encode the evolution
of their spins prior to merger. In the postNewtonian regime and on the
precession timescale, this evolution has one of three morphologies, with the
spins either librating around one of two fixed points ("resonances") or
circulating freely. In this work we perform full parameter estimation on
resonant binaries with fixed masses and spin magnitudes, changing three
parameters: a conserved "projected effective spin" $\xi$ and resonant family
$\Delta\Phi=0,\pi$ (which uniquely label the source), the inclination
$\theta_{JN}$ of the binary's total angular momentum with respect to the line
of sight (which determines the strength of precessional effects in the
waveform), and the signal amplitude. We demonstrate that resonances can be
distinguished for a wide range of binaries, except for highly symmetric
configurations where precessional effects are suppressed. Motivated by new
insight into doublespin evolution, we introduce new variables to characterize
precessing black hole binaries which naturally reflects the timescale
separation of the system and therefore better encode the dynamical information
carried by gravitational waves.

One century after its formulation, Einstein's general relativity has made
remarkable predictions and turned out to be compatible with all experimental
tests. Most of these tests probe the theory in the weakfield regime, and there
are theoretical and experimental reasons to believe that general relativity
should be modified when gravitational fields are strong and spacetime curvature
is large. The best astrophysical laboratories to probe strongfield gravity are
black holes and neutron stars, whether isolated or in binary systems. We review
the motivations to consider extensions of general relativity. We present a
(necessarily incomplete) catalog of modified theories of gravity for which
strongfield predictions have been computed and contrasted to Einstein's
theory, and we summarize our current understanding of the structure and
dynamics of compact objects in these theories. We discuss current bounds on
modified gravity from binary pulsar and cosmological observations, and we
highlight the potential of future gravitational wave measurements to inform us
on the behavior of gravity in the strongfield regime.

Binary black holes on quasicircular orbits with spins aligned with their
orbital angular momentum have been testbeds for analytic and numerical
relativity for decades, not least because symmetry ensures that such
configurations are equilibrium solutions to the spinprecession equations. In
this work, we show that these solutions can be unstable when the spin of the
highermass black hole is aligned with the orbital angular momentum and the
spin of the lowermass black hole is antialigned. Spins in these
configurations are unstable to precession to large misalignment when the binary
separation $r$ is between the values $r_{\rm ud\pm}= (\sqrt{\chi_1} \pm \sqrt{q
\chi_2})^4 (1q)^{2} M$, where $M$ is the total mass, $q \equiv m_2/m_1$ is
the mass ratio, and $\chi_1$ ($\chi_2$) is the dimensionless spin of the more
(less) massive black hole. This instability exists for a wide range of spin
magnitudes and mass ratios and can occur in the strongfield regime near
merger. We describe the origin and nature of the instability using recently
developed analytical techniques to characterize fully generic spin precession.
This instability provides a channel to circumvent astrophysical spin alignment
at large binary separations, allowing significant spin precession prior to
merger affecting both gravitationalwave and electromagnetic signatures of
stellarmass and supermassive binary black holes.

The dynamics of precessing binary black holes (BBHs) in the postNewtonian
regime has a strong timescale hierarchy: the orbital timescale is very short
compared to the spinprecession timescale which, in turn, is much shorter than
the radiationreaction timescale on which the orbit is shrinking due to
gravitationalwave emission. We exploit this timescale hierarchy to develop a
multiscale analysis of BBH dynamics elaborating on the analysis of Kesden et
al. (2015). We solve the spinprecession equations analytically on the
precession time and then implement a quasiadiabatic approach to evolve these
solutions on the longer radiationreaction time. This procedure leads to an
innovative "precessionaveraged" postNewtonian approach to studying precessing
BBHs. We use our new solutions to classify BBH spin precession into three
distinct morphologies, then investigate phase transitions between these
morphologies as BBHs inspiral. These precessionaveraged postNewtonian
inspirals can be efficiently calculated from arbitrarily large separations,
thus making progress towards bridging the gap between astrophysics and
numerical relativity.

Gravitational theories with multiple scalar fields coupled to the metric and
each other  a natural extension of the well studied singlescalartensor
theories  are interesting phenomenological frameworks to describe deviations
from general relativity in the strongfield regime. In these theories, the
$N$tuple of scalar fields takes values in a coordinate patch of an
$N$dimensional Riemannian targetspace manifold whose properties are poorly
constrained by weakfield observations. Here we introduce for simplicity a
nontrivial model with two scalar fields and a maximally symmetric targetspace
manifold. Within this model we present a preliminary investigation of
spontaneous scalarization for relativistic, perfect fluid stellar models in
spherical symmetry. We find that the scalarization threshold is determined by
the eigenvalues of a symmetric scalarmatter coupling matrix, and that the
properties of strongly scalarized stellar configurations additionally depend on
the targetspace curvature radius. In preparation for numerical relativity
simulations, we also write down the $3+1$ decomposition of the field equations
for generic tensormultiscalar theories.

Interactions between a supermassive black hole binary and the surrounding
accretion disc can both assist the binary inspiral and align the black hole
spins to the disc angular momentum. While binary migration is due to
angularmomentum transfer within the circumbinary disc, the spinalignment
process is driven by the mass accreting on to each black hole. Mass transfer
between different disc components thus couples the inspiral and the alignment
process together. Mass is expected to leak through the cavity cleared by the
binary, and preferentially accretes on to the lighter (secondary) black hole
which orbits closer to the disc edge. Low accretion rate on to the heavier
(primary) black hole slows the alignment process down. We revisit the problem
and develop a semianalytical model to describe the coupling between gasdriven
inspiral and spin alignment, finding that binaries with mass ratio q<~0.2
approach the gravitationalwave driven inspiral in differential misalignment:
light secondaries prevent primaries from aligning. Binary black holes with
misaligned primaries are ideal candidates for precession effects in the
stronggravity regime and may suffer from moderately large (~1500 km/s) recoil
velocities.

We derive an effective potential for binary blackhole (BBH) spin precession
at second postNewtonian order. This effective potential allows us to solve the
orbitaveraged spinprecession equations analytically for arbitrary mass ratios
and spins. These solutions are quasiperiodic functions of time: after a fixed
period the BBH spins return to their initial relative orientations and jointly
precess about the total angular momentum by a fixed angle. Using these
solutions, we classify BBH spin precession into three distinct morphologies
between which BBHs can transition during their inspiral. We also derive a
precessionaveraged evolution equation for the total angular momentum that can
be integrated on the radiationreaction time and identify a new class of
spinorbit resonances that can tilt the direction of the total angular momentum
during the inspiral. Our new results will help efforts to model and interpret
gravitational waves from generic BBH mergers and predict the distributions of
final spins and gravitational recoils.

We investigate the consequences of superkicks on the population of
supermassive black holes (SMBHs) in the Universe residing in brightest cluster
galaxies (BCGs). There is strong observational evidence that BCGs grew
prominently at late times (up to a factor 24 in mass from z=1), mainly through
mergers with satellite galaxies from the cluster, and they are known to host
the most massive SMBHs ever observed. Those SMBHs are also expected to grow
hierarchically, experiencing a series of mergers with other SMBHs brought in by
merging satellites. Because of the net linear momentum taken away from the
asymmetric gravitational wave emission, the remnant SMBH experiences a kick in
the opposite direction. Kicks may be as large as ~5000 Km/s ("superkicks"),
pushing the SMBHs out in the cluster outskirts for a time comparable to
galaxyevolution timescales. We predict, under a number of plausible
assumptions, that superkicks can efficiently eject SMBHs from BCGs, bringing
their occupation fraction down to a likely range 0.9<f<0.99 in the local
Universe. Future thirtymeterclass telescopes like ELT and TMT will be capable
of measuring SMBHs in hundreds of BCGs up to z=0.2, testing the occurrence of
superkicks in nature and the stronggravity regime of SMBH mergers.

If binary black holes form following the successive core collapses of
sufficiently massive binary stars, precessional dynamics may align their spins
$\mathbf S_1$ and $\mathbf S_2$ and the orbital angular momentum $\mathbf L$
into a plane in which they jointly precess about the total angular momentum
$\mathbf J$. These spin orientations are known as spinorbit resonances since
$\mathbf S_1$, $\mathbf S_2$, and $\mathbf L$ all precess at the same frequency
to maintain their planar configuration. Two families of such spinorbit
resonances exist, alike in dignity but differentiated by whether the components
of the two spins in the orbital plane are either aligned or antialigned. The
fraction of binary black holes in each family is determined by the stellar
evolution of their progenitors, so if gravitationalwave detectors could
measure this fraction they could provide important insights into astrophysical
formation scenarios for binary black holes. In this paper, we show that even
under the conservative assumption that binary black holes are observed along
the direction of $\mathbf J$ (where precessioninduced modulations to the
gravitational waveforms are minimized), the waveforms of many members of each
resonant family can be distinguished from all members of the opposing family in
events with signaltonoise ratios $\rho \simeq 10$, typical of those expected
for the first detections with Advanced LIGO/Virgo. We hope that our preliminary
findings inspire a greater appreciation of the capability of gravitationalwave
detectors to constrain stellar astrophysics and stimulate further studies of
the distinguishability of spinorbit resonant families in more expanded regions
of binary blackhole parameter space.

We study the influence of astrophysical formation scenarios on the
precessional dynamics of spinning blackhole binaries by the time they enter
the observational window of second and thirdgeneration gravitationalwave
detectors, such as Advanced LIGO/Virgo, LIGOIndia, KAGRA and the Einstein
Telescope. Under the plausible assumption that tidal interactions are efficient
at aligning the spins of fewsolar mass blackhole progenitors with the orbital
angular momentum, we find that blackhole spins should be expected to
preferentially lie in a plane when they become detectable by gravitationalwave
interferometers. This "resonant plane" is identified by the conditions
\Delta\Phi=0{\deg} or \Delta\Phi=+/180{\deg}, where \Delta\Phi is the angle
between the components of the blackhole spins in the plane orthogonal to the
orbital angular momentum. If the angles \Delta \Phi can be accurately measured
for a large sample of gravitationalwave detections, their distribution will
constrain models of compact binary formation. In particular, it will tell us
whether tidal interactions are efficient and whether a mechanism such as mass
transfer, stellar winds, or supernovae can induce a massratio reversal (so
that the heavier black hole is produced by the initially lighter stellar
progenitor). Therefore our model offers a concrete observational link between
gravitationalwave measurements and astrophysics. We also hope that it will
stimulate further studies of precessional dynamics, gravitationalwave template
placement and parameter estimation for binaries locked in the resonant plane.