Massive black hole (MBH) binary inspiral time scales are uncertain, and their spins are even more poorly constrained. Spin misalignment, along with unequal mass ratios and spin magnitudes, introduces asymmetry in the gravitational radiation, which imparts a recoil kick to the merged MBH. Understanding how MBH binary spins evolve is crucial for determining their recoil velocities, their gravitational wave (GW) waveforms detectable with LISA, as well as their post-merger retention rate in galaxies and thus their subsequent merger rate. Here we present a novel study that introduces a sub-resolution model for gas- and GW-driven MBH binary spin evolution using a population of accreting MBHs from the Illustris cosmological hydrodynamics simulations. We also model sub-resolution binary inspiral via dynamical friction, stellar scattering, viscous gas drag, and GW emission. Our model assumes differential accretion, which causes greater alignment of the secondary MBH spin in unequal-mass mergers. We find that 47% of the MBHs in our population merge by $z=0$. Of these, 19% have misaligned primaries and 10% have misaligned secondaries at the time of merger in our (conservative) fiducial model. The MBH misalignment fraction depends strongly on the accretion disc parameters, however. Reducing accretion rates by a factor of 100, in a thicker disc, yields 79% and 42% misalignment for primaries and secondaries, respectively. Even in the fiducial model, more than 12% of binaries experience recoils of $>500$ km s$^{-1}$, which could displace them at least temporarily from galactic nuclei. We additionally find that a significant number of systems experience strong precession.

Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astrophysical questions by estimating tidal disruption event rates from observational samples of flares? These are the questions we aim to address in this Chapter, which summarizes current theoretical knowledge about rates of stellar tidal disruption, and compares theoretical predictions to the current state of observations.

Within the framework of 2PN black-hole binary spin precession, we explore configurations where one of the two spins oscillates from being completely aligned with the orbital angular momentum to being completely anti-aligned with it during a single precession cycle. This "wide nutation" is the extreme limit of the generic phenomenon of spin nutation in black-hole binaries. Crucially, wide nutation happens on the short precession time scale and it is not a secular effect due to gravitational-wave radiation reaction. The spins of these binaries, therefore, flip repeatedly as one of these special configurations is entered. Binaries with total mass $M$, mass ratio $q$, and dimensionless spin $\chi_1$ ($\chi_2$) of the more (less) massive black hole are allowed to undergo wide nutation at binary separations $r \leq r_{\rm wide} \equiv [(q \chi_2 - \chi_1)/(1-q)]^2 M$. Sources that are more likely to nutate widely have similar masses and effective spins close to zero.

We study the expected spin misalignments of merging binary black holes (BHs) formed in isolation by combining state-of-the-art population-synthesis models with efficient post-Newtonian evolutions, thus tracking sources from stellar formation to gravitational-wave detection. We present extensive predictions of the properties of sources detectable by both current and future interferometers. We account for the fact that detectors are more sensitive to spinning BH binaries with suitable spin orientations and find that this significantly impacts the population of sources detectable by LIGO, while this is not the case for 3rd-generation detectors. We find that three formation pathways, differentiated by the order of core collapse and common-envelope phases, dominate the observed population, and that their relative importance critically depends on the recoils imparted to BHs at birth. Our models suggest that measurements of the "effective-spin" parameter $\chi_{\rm eff}$ will allow for powerful constraints. For instance, we find that the role of spin magnitudes and spin directions in $\chi_{\rm eff}$ can be largely disentangled, and that the symmetry of the effective-spin distribution is a robust indicator of the binary's formation history. Our predictions for individual spin directions and their precessional morphologies confirm and extend early toy models, while exploring substantially more realistic and broader sets of initial conditions. Our main conclusion is that specific subpopulations of BH binaries will exhibit distinctive precessional dynamics: these classes include (but are not limited to) sources where stellar tidal interactions act on sufficiently short timescales, and massive binaries produced in pulsational pair-instability supernovae. Measurements of BH spin orientations have enormous potential to constrain specific evolutionary processes in the lives of massive binary stars.

We compare binary evolution models with different assumptions about black-hole natal kicks to the first gravitational-wave observations performed by the LIGO detectors. Our comparisons attempt to reconcile merger rate, masses, spins, and spin-orbit misalignments of all current observations with state-of-the-art formation scenarios of binary black holes formed in isolation. We estimate that black holes (BHs) should receive natal kicks at birth of the order of $\sigma\simeq 200$ (50) km/s if tidal processes do (not) realign stellar spins. Our estimate is driven by two simple factors. The natal kick dispersion $\sigma$ is bounded from above because large kicks disrupt too many binaries (reducing the merger rate below the observed value). Conversely, the natal kick distribution is bounded from below because modest kicks are needed to produce a range of spin-orbit misalignments. A distribution of misalignments increases our models' compatibility with LIGO's observations, if all BHs are likely to have natal spins. Unlike related work which adopts a concrete BH natal spin prescription, we explore a range of possible BH natal spin distributions. Within the context of our models, for all of the choices of $\sigma$ used here and within the context of one simple fiducial parameterized spin distribution, observations favor low BH natal spin.

A tidal disruption event (TDE) ensues when a star passes too close to the supermassive black hole (SMBH) in a galactic center and is ripped apart by the tidal field of the SMBH. The gaseous debris produced in a TDE can power a bright electromagnetic flare as it is accreted by the SMBH; so far, several dozen TDE candidates have been observed. For SMBHs with masses above $\sim 10^7 M_\odot$, the tidal disruption of solar-type stars occurs within ten gravitational radii of the SMBH, implying that general relativity (GR) is needed to describe gravity. Three promising signatures of GR in TDEs are: (1) a super-exponential cutoff in the volumetric TDE rate for SMBH masses above $\sim 10^8 M_\odot$ due to direct capture of tidal debris by the event horizon, (2) delays in accretion disk formation (and a consequent alteration of the early-time light curve) caused by the effects of relativistic precession on stream circularization, and (3) quasi-periodic modulation of X-ray emission due to global precession of misaligned accretion disks and the jets they launch. We review theoretical models and simulations of TDEs in Newtonian gravity, then describe how relativistic modifications give rise to these proposed observational signatures, as well as more speculative effects of GR. We conclude with a brief summary of TDE observations and the extent to which they show indications of these predicted relativistic signatures.

We present the numerical code PRECESSION: a new open-source python module to study the dynamics of precessing black-hole binaries in the post-Newtonian regime. The code provides a comprehensive toolbox to (i) study the evolution of the black-hole spins along their precession cycles, (ii) perform gravitational-wave driven binary inspirals using both orbit-averaged and precession-averaged integrations, and (iii) predict the properties of the merger remnant through fitting formulae obtained from numerical-relativity simulations. PRECESSION is a ready-to-use tool to add the black-hole spin dynamics to larger-scale numerical studies such as gravitational-wave parameter estimation codes, population synthesis models to predict gravitational-wave event rates, galaxy merger trees and cosmological simulations of structure formation. PRECESSION provides fast and reliable integration methods to propagate statistical samples of black-hole binaries from/to large separations where they form to/from small separations where they become detectable, thus linking gravitational-wave observations of spinning black-hole binaries to their astrophysical formation history. The code is also a useful tool to compute initial parameters for numerical-relativity 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.

In the post-Newtonian (PN) regime, the timescale on which the spins of binary black holes precess is much shorter than the radiation-reaction 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 radiation-reaction 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.

The recent discovery of gravitational waves (GWs) from stellar-mass binary black holes (BBHs) provided direct evidence of the existence of these systems. BBH lenses would have gravitational microlensing signatures that are distinct from single-lens signals. We apply Bayesian statistics to examine the distinguishability of BBH microlensing events from single-lens events under ideal observing conditions, using the photometric capabilities of the Korean Microlensing Telescope Network. Given one year of observations, a source star at the Galactic Centre, a GW150914-like BBH lens (total mass 65M$_\odot$, mass ratio 0.8) at half that distance, and an impact parameter of 0.4 Einstein radii, we find that binarity is detectable for BBHs with separations down to 0.0250 Einstein radii, which is nearly 3.5 times greater than the maximum separation for which such BBHs would merge within the age of the Universe. Microlensing searches are thus sensitive to more widely separated BBHs than GW searches, perhaps allowing the discovery of BBH populations produced in different channels of binary formation.

A supermassive black hole moving through a field of stars will gravitationally scatter the stars, inducing a backreaction force on the black hole known as dynamical friction. In Newtonian gravity, the axisymmetry of the system about the black hole's velocity $\mathbf{v}$ implies that the dynamical friction must be anti-parallel to $\mathbf{v}$. However, in general relativity the black hole's spin $\mathbf{S}$ need not be parallel to $\mathbf{v}$, breaking the axisymmetry of the system and generating a new component of dynamical friction similar to the Lorentz force $\mathbf{F} = q\mathbf{v} \times \mathbf{B}$ experienced by a particle with charge $q$ moving in a magnetic field $\mathbf{B}$. We call this new force gravitomagnetic dynamical friction and calculate its magnitude for a spinning black hole moving through a field of stars with Maxwellian velocity dispersion $\sigma$, assuming that both $v$ and $\sigma$ are much less than the speed of light $c$. We use post-Newtonian equations of motion accurate to $\mathcal{O}(v^3/c^3)$ needed to capture the effect of spin-orbit coupling and also include direct stellar capture by the black hole's event horizon. Gravitomagnetic dynamical friction will cause a black hole with uniform speed to spiral about the direction of its spin, similar to a charged particle spiraling about a magnetic field line, and will exert a torque on a supermassive black hole orbiting a galactic center, causing the angular momentum of this orbit to slowly precess about the black-hole spin. As this effect is suppressed by a factor $(\sigma/c)^2$ in nonrelativistic systems, we expect it to be negligible in most astrophysical contexts but provide this calculation for its theoretical interest and potential application to relativistic systems.

Stars on orbits with pericenters sufficiently close to the supermassive black hole at the center of their host galaxy can be ripped apart by tidal stresses. Some of the resulting stellar debris becomes more tightly bound to the hole and can potentially produce an observable flare called a tidal-disruption event (TDE). We provide a self-consistent, unified treatment of TDEs by non-spinning (Schwarzschild) black holes, investigating several effects of general relativity including changes to the boundary in phase space that defines the loss-cone orbits on which stars are tidally disrupted or captured. TDE rates decrease rapidly at large black-hole masses due to direct stellar capture, but this effect is slightly countered by the widening of the loss cone due to the stronger tidal fields in general relativity. We provide a new mapping procedure that translates between Newtonian gravity and general relativity, allowing us to better compare predictions in both gravitational theories. Partial tidal disruptions in relativity will strip more material from the star and produce more tightly bound debris than in Newtonian gravity for a stellar orbit with the same angular momentum. However, for deep encounters leading to full disruption in both theories, the stronger tidal forces in relativity imply that the star is disrupted further from the black hole and that the debris is therefore less tightly bound, leading to a smaller peak fallback accretion rate. We also examine the capture of tidal debris by the horizon and the relativistic pericenter precession of tidal debris, finding that black holes of $10^6$ solar masses and above generate tidal debris precessing by $10^\circ$ or more per orbit.

Gravitational waves from coalescing binary black holes encode the evolution of their spins prior to merger. In the post-Newtonian 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 double-spin 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.

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 spin-precession equations. In this work, we show that these solutions can be unstable when the spin of the higher-mass black hole is aligned with the orbital angular momentum and the spin of the lower-mass black hole is anti-aligned. 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 (1-q)^{-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 strong-field 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 gravitational-wave and electromagnetic signatures of stellar-mass and supermassive binary black holes.

The dynamics of precessing binary black holes (BBHs) in the post-Newtonian regime has a strong timescale hierarchy: the orbital timescale is very short compared to the spin-precession timescale which, in turn, is much shorter than the radiation-reaction timescale on which the orbit is shrinking due to gravitational-wave emission. We exploit this timescale hierarchy to develop a multi-scale analysis of BBH dynamics elaborating on the analysis of Kesden et al. (2015). We solve the spin-precession equations analytically on the precession time and then implement a quasi-adiabatic approach to evolve these solutions on the longer radiation-reaction time. This procedure leads to an innovative "precession-averaged" post-Newtonian 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 precession-averaged post-Newtonian inspirals can be efficiently calculated from arbitrarily large separations, thus making progress towards bridging the gap between astrophysics and numerical relativity.

We derive an effective potential for binary black-hole (BBH) spin precession at second post-Newtonian order. This effective potential allows us to solve the orbit-averaged spin-precession 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 precession-averaged evolution equation for the total angular momentum that can be integrated on the radiation-reaction time and identify a new class of spin-orbit 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.

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 spin-orbit 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 spin-orbit 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 gravitational-wave 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 precession-induced 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 signal-to-noise 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 gravitational-wave detectors to constrain stellar astrophysics and stimulate further studies of the distinguishability of spin-orbit resonant families in more expanded regions of binary black-hole parameter space.

We study the influence of astrophysical formation scenarios on the precessional dynamics of spinning black-hole binaries by the time they enter the observational window of second- and third-generation gravitational-wave detectors, such as Advanced LIGO/Virgo, LIGO-India, KAGRA and the Einstein Telescope. Under the plausible assumption that tidal interactions are efficient at aligning the spins of few-solar mass black-hole progenitors with the orbital angular momentum, we find that black-hole spins should be expected to preferentially lie in a plane when they become detectable by gravitational-wave 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 black-hole spins in the plane orthogonal to the orbital angular momentum. If the angles \Delta \Phi can be accurately measured for a large sample of gravitational-wave 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 mass-ratio reversal (so that the heavier black hole is produced by the initially lighter stellar progenitor). Therefore our model offers a concrete observational link between gravitational-wave measurements and astrophysics. We also hope that it will stimulate further studies of precessional dynamics, gravitational-wave template placement and parameter estimation for binaries locked in the resonant plane.

A star orbiting a supermassive black hole can be tidally disrupted if the black hole's gravitational tidal field exceeds the star's self gravity at pericenter. Some of this stellar tidal debris can become gravitationally bound to the black hole, leading to a bright electromagnetic flare with bolometric luminosity proportional to the rate at which material falls back to pericenter. In the Newtonian limit, this flare will have a light curve that scales as t^-5/3 if the tidal debris has a flat distribution in binding energy. We investigate the time dependence of the black-hole mass accretion rate when tidal disruption occurs close enough the black hole that relativistic effects are significant. We find that for orbits with pericenters comparable to the radius of the marginally bound circular orbit, relativistic effects can double the peak accretion rate and halve the time it takes to reach this peak accretion rate. The accretion rate depends on both the magnitude of the black-hole spin and its orientation with respect to the stellar orbit; for orbits with a given pericenter radius in Boyer-Lindquist coordinates, a maximal black-hole spin anti-aligned with the orbital angular momentum leads to the largest peak accretion rate.

Recent numerical relativity simulations have shown that the final black hole produced in a binary merger can recoil with a velocity as large as 5,000 km/s. Because of enhanced gravitational-wave emission in the so-called "hang-up" configurations, this maximum recoil occurs when the black-hole spins are partially aligned with the orbital angular momentum. We revisit our previous statistical analysis of post-Newtonian evolutions of black-hole binaries in the light of these new findings. We demonstrate that despite these new configurations with enhanced recoil velocities, spin alignment during the post-Newtonian stage of the inspiral will still significantly suppress (or enhance) kick magnitudes when the initial spin of the more massive black hole is more (or less) closely aligned with the orbital angular momentum than that of the smaller hole. We present a preliminary study of how this post-Newtonian spin alignment affects the ejection probabilities of supermassive black holes from their host galaxies with astrophysically motivated mass ratio and initial spin distributions. We find that spin alignment suppresses (enhances) ejection probabilities by ~ 40% (20%) for an observationally motivated mass-dependent galactic escape velocity, and by an even greater amount for a constant escape velocity of 1,000 km/s. Kick suppression is thus at least a factor two more efficient than enhancement.

A supermassive black hole can disrupt a star when its tidal field exceeds the star's self-gravity, and can directly capture stars that cross its event horizon. For black holes with mass M > 10^7 solar masses, tidal disruption of main-sequence stars occurs close enough to the event horizon that a Newtonian treatment of the tidal field is no longer valid. The fraction of stars that are directly captured is also no longer negligible. We calculate generically oriented stellar orbits in the Kerr metric, and evaluate the relativistic tidal tensor at pericenter for those stars not directly captured by the black hole. We combine this relativistic analysis with previous calculations of how these orbits are populated to determine tidal-disruption rates for spinning black holes. We find, consistent with previous results, that black-hole spin increases the upper limit on the mass of a black hole capable of tidally disrupting solar-like stars to ~7 x 10^8 solar masses. More quantitatively, we find that direct stellar capture reduces tidal-disruption rates by a factor 2/3 (1/10) at M = 10^7 (10^8) solar masses. The strong dependence of tidal-disruption rates on black-hole spin for M > 10^8 solar masses implies that future surveys like LSST that discover thousands of tidal disruption events can constrain supermassive black-hole spin demographics.

Stars are transparent to the passage of primordial black holes (PBHs) and serve as seismic detectors for such objects. The gravitational field of a PBH squeezes a star and causes it to ring acoustically. We calculate the seismic signature of a PBH passing through the Sun. The background for this signal is the observed spectrum of solar oscillations excited by supersonic turbulence. We predict that PBHs more massive than 10^21 g (comparable in mass to an asteroid) are detectable by existing solar observatories. The oscillations excited by PBHs peak at large scales and high frequencies, making them potentially detectable in other stars. The discovery of PBHs would have profound implications for cosmology and high-energy physics.

A test particle of mass mu on a bound geodesic of a Kerr black hole of mass M >> mu will slowly inspiral as gravitational radiation extracts energy and angular momentum from its orbit. This inspiral can be considered adiabatic when the orbital period is much shorter than the timescale on which energy is radiated, and quasi-circular when the radial velocity is much less than the azimuthal velocity. Although the inspiral always remains adiabatic provided mu << M, the quasi-circular approximation breaks down as the particle approaches the innermost stable circular orbit (ISCO). In this paper, we relax the quasi-circular approximation and solve the radial equation of motion explicitly near the ISCO. We use the requirement that the test particle's 4-velocity remain properly normalized to calculate a new contribution to the difference between its energy and angular momentum. This difference determines how a black hole's spin changes following a test-particle merger, and can be extrapolated to help predict the mass and spin of the final black hole produced in finite-mass-ratio black-hole mergers. Our new contribution is particularly important for nearly maximally spinning black holes, as it can affect whether a merger produces a naked singularity.

Black holes of mass M must have a spin angular momentum S below the Kerr limit chi = S/M^2 < 1, but whether astrophysical black holes can attain this limiting spin depends on their accretion history. Gas accretion from a thin disk limits the black-hole spin to chi_gas < 0.9980 +- 0.0002, as electromagnetic radiation from this disk with retrograde angular momentum is preferentially absorbed by the black hole. Extrapolation of numerical-relativity simulations of equal-mass binary black-hole mergers to maximum initial spins suggests these mergers yield a maximum spin chi_eq < 0.95. Here we show that for smaller mass ratios q = m/M << 1, the superradiant extraction of angular momentum from the larger black hole imposes a fundamental limit chi_lim < 0.9979 +- 0.0001 on the final black-hole spin even in the test-particle limit q -> 0 of binary black-hole mergers. The nearly equal values of chi_gas and chi_lim imply that measurement of supermassive black-hole spins cannot distinguish a black hole built by gas accretion from one assembled by the gravitational inspiral of a disk of compact stellar remnants. We also show how superradiant scattering alters the mass and spin predicted by models derived from extrapolating test-particle mergers to finite mass ratios.

Numerical-relativity simulations indicate that the black hole produced in a binary merger can recoil with a velocity up to v_max ~ 4,000 km/s with respect to the center of mass of the initial binary. This challenges the paradigm that most galaxies form through hierarchical mergers, yet retain supermassive black holes at their centers despite having escape velocities much less than v_max. Interaction with a circumbinary disk can align the binary black hole spins with their orbital angular momentum, reducing the recoil velocity of the final black hole produced in the subsequent merger. However, the effectiveness of this alignment depends on highly uncertain accretion flows near the binary black holes. In this Letter, we show that if the spin S_1 of the more massive binary black hole is even partially aligned with the orbital angular momentum L, relativistic spin precession on sub-parsec scales can align the binary black hole spins with each other. This alignment significantly reduces the recoil velocity even in the absence of gas. For example, if the angle between S_1 and L at large separations is 10 degrees while the second spin S_2 is isotropically distributed, the spin alignment discussed in this paper reduces the median recoil from 864 km/s to 273 km/s for maximally spinning black holes with a mass ratio of 9/11. This reduction will greatly increase the fraction of galaxies retaining their supermassive black holes.

The inspiral of binary black holes is governed by gravitational radiation reaction at binary separations r < 1000 M, yet it is too computationally expensive to begin numerical-relativity simulations with initial separations r > 10 M. Fortunately, binary evolution between these separations is well described by post-Newtonian equations of motion. We examine how this post-Newtonian evolution affects the distribution of spin orientations at separations r ~ 10 M where numerical-relativity simulations typically begin. Although isotropic spin distributions at r ~ 1000 M remain isotropic at r ~ 10 M, distributions that are initially partially aligned with the orbital angular momentum can be significantly distorted during the post-Newtonian inspiral. Spin precession tends to align (anti-align) the binary black hole spins with each other if the spin of the more massive black hole is initially partially aligned (anti-aligned) with the orbital angular momentum, thus increasing (decreasing) the average final spin. Spin precession is stronger for comparable-mass binaries, and could produce significant spin alignment before merger for both supermassive and stellar-mass black hole binaries. We also point out that precession induces an intrinsic accuracy limitation (< 0.03 in the dimensionless spin magnitude, < 20 degrees in the direction) in predicting the final spin resulting from the merger of widely separated binaries.

The unique dynamics of the tidal disruption of satellite galaxies is an extremely sensitive probe of long-range interactions between dark-matter particles. Dark-matter forces that are several percent the strength of gravity will lead to order unity changes in the ratio of the number of stars in the leading and trailing tidal streams of a satellite galaxy. The approximate symmetry of the stellar tidal streams of the Sagittarius dwarf galaxy would thus exclude attractive dark-matter forces greater than 10% the strength of gravity which would entirely eliminate the leading stream. However, recent simulations suggest that dark-matter forces 100% the strength of gravity could completely strip the stellar component of Sagittarius of its dark matter, allowing for the subsequent development of symmetric tidal streams. Here we argue that these simulations use inconsistent initial conditions corresponding to separate pure stellar and pure dark-matter satellites moving independently in the host galaxy's halo, rather than a single disrupting composite satellite as had been intended. A new simulation with different initial conditions, in particular a much more massive satellite galaxy, might demonstrate a scenario in which symmetric tidal streams develop in the presence of large dark-matter forces. This scenario must satisfy several highly restrictive criteria described in this paper.

Gravitational waves carry away both energy and angular momentum as binary black holes inspiral and merge. The relative efficiency with which they are radiated determines whether the final black hole of mass $M_f$ and spin $S_f$ saturates the Kerr limit ($\chi_f \equiv S_f/M_f^2 \leq 1$). Extrapolating from the test-particle limit, we propose expressions for $S_f$ and $M_f$ for mergers with initial spins aligned or anti-aligned with the orbital angular momentum. We predict the the final spin at plunge for equal-mass non-spinning binaries to better than 1%, and that equal-mass maximally spinning aligned mergers lead to nearly maximally spinning final black holes ($\chi_f \simeq 0.9988$). We also find black holes can always be spun up by aligned mergers provided the mass ratio is small enough.

In a recent paper arXiv:0709.0299, we introduced a spin expansion that provides a simple yet powerful way to understand aspects of binary black hole (BBH) merger. This approach relies on the symmetry properties of initial and final quantities like the black hole mass m, kick velocity {\bf k}, and spin vector {\bf s}, rather than a detailed understanding of the merger dynamics. In this paper, we expand on this proposal, examine how well its predictions agree with current simulations, and discuss several future directions that would make it an even more valuable tool. The spin expansion yields many new predictions, including several exact results that may be useful for testing numerical codes. Some of these predictions have already been confirmed, while others await future simulations. We explain how a relatively small number of simulations -- 10 equal-mass simulations, and 16 unequal-mass simulations -- may be used to calibrate all of the coefficients in the spin expansion up to second order at the minimum computational cost. For a more general set of simulations of given covariance, we derive the minimum-variance unbiased estimators for the spin expansion coefficients. We discuss how this calibration would be interesting and fruitful for general relativity and astrophysics. Finally, we sketch the extension to eccentric orbits.

We regard binary black hole (BBH) merger as a map from a simple initial state (two Kerr black holes, with dimensionless spins {\bf a} and {\bf b}) to a simple final state (a Kerr black hole with mass m, dimensionless spin {\bf s}, and kick velocity {\bf k}). By expanding this map around {\bf a} = {\bf b} = 0 and applying symmetry constraints, we obtain a simple formalism that is remarkably successful at explaining existing BBH simulations. It also makes detailed predictions and suggests a more efficient way of mapping the parameter space of binary black hole merger. Since we rely on symmetry rather than dynamics, our expansion complements previous analytical techniques.

Satellite galaxies currently undergoing tidal disruption offer a unique opportunity to constrain an effective violation of the equivalence principle in the dark sector. Theories in which cold dark matter (CDM) couples to a light scalar field naturally lead to a long-range force between dark matter particles. An inverse-square-law force of this kind would manifest itself as a violation of the equivalence principle in the dynamics of CDM compared to baryons in the form of gas or stars. In a previous paper, we showed that an attractive force would displace stars outwards from the bottom of the satellite's gravitational potential well, leading to a higher fraction of stars being disrupted from the tidal bulge further from the Galactic center. Since stars disrupted from the far (near) side of the satellite go on to form the trailing (leading) tidal stream, an attractive dark-matter force will produce a relative enhancement of the trailing stream compared to the leading stream. This distinctive signature of a dark-matter force might be detected through detailed observations of the tidal tails of a disrupting satellite, such as those recently performed by the Two-Micron All-Sky Survey (2MASS) and Sloan Digital Sky Survey (SDSS) on the Sagittarius (Sgr) dwarf galaxy. Here we show that this signature is robust to changes in our models for both the satellite and Milky Way, suggesting that we might hope to search for a dark-matter force in the tidal features of other recently discovered satellite galaxies in addition to the Sgr dwarf.

Satellite galaxies are tidally disrupted as they orbit the Milky Way. If dark matter (DM) experiences a stronger self-attraction than baryons, stars will preferentially gain rather than lose energy during tidal disruption leading to an enhancement in the trailing compared to the leading tidal stream. The Sgr dwarf galaxy is seen to have roughly equal streams, challenging models in which DM and baryons accelerate differently by more than 10%. Future observations and a better understanding of DM distribution should allow detection of equivalence violation at the percent level.

Event horizons are among the most intriguing of general relativity's predictions. Although on firm theoretical footing, direct indications of their existence have yet to be observed. With this motivation in mind, we explore here the possibility of finding a signature for event horizons in the gravitational waves (GWs) produced during the inspiral of stellar-mass compact objects (COs) into the supermassive ($\sim 10^6 M_\odot$) objects that lie at the center of most galaxies. Such inspirals will be a major source for LISA, the future space-based GW observatory. We contrast supermassive black holes with models in which the central object is a supermassive boson star (SMBS). Provided the COs interact only gravitationally with the SMBS, stable orbits exist not just outside the Schwarzschild radius but also inside the surface of the SMBS as well. The absence of an event horizon allows GWs from these orbits to be observed. Here we solve for the metric in the interior of a fairly generic class of SMBS and evolve the trajectory of an inspiraling CO from the Schwarzschild exterior through the plunge into the exotic SMBS interior. We calculate the approximate waveforms for GWs emitted during this inspiral. Geodesics within the SMBS surface will exhibit extreme pericenter precession and other features making the emitted GWs readily distinguishable from those emitted during an inspiral into a black hole.

The lowest multipole moments of the cosmic microwave background (CMB) are smaller than expected for a scale-invariant power spectrum. One possible explanation is a cutoff in the primordial power spectrum below a comoving scale of $k_c \simeq 5.0 \times 10^{-4}$ Mpc$^{-1}$. This would affect not only the CMB but also the cosmic-shear (CS) distortion of the CMB. Such a cutoff increases significantly the cross-correlation between the large-angle CMB and cosmic-shear patterns. The cross-correlation may be detectable at $> 2\sigma$ which, when combined with the low CMB moments, may tilt the balance between a $2\sigma$ result and a firm detection of a large-scale power-spectrum cutoff. As an aside, we also note that the cutoff increases the large-angle cross-correlation between the CMB and low-redshift tracers of the mass distribution.

Weak gravitational lensing by intervening large-scale structure induces a distinct signature in the cosmic microwave background (CMB) that can be used to reconstruct the weak-lensing displacement map. Estimators for individual Fourier modes of this map can be combined to produce an estimator for the lensing-potenial power spectrum. The naive estimator for this quantity will be biased upwards by the uncertainty associated with reconstructing individual modes; we present an iterative scheme for removing this bias. The variance and covariance of the lensing-potenial power spectrum estimator are calculated and evaluated numerically in a $\Lambda$CDM universe for Planck and future polarization-sensitive CMB experiments.

We discuss the real-space moments of temperature anisotropies in the cosmic microwave background (CMB) due to weak gravitational lensing by intervening large-scale structure. We show that if the probability distribution function of primordial temperature anisotropies is Gaussian, then it remains unchanged after gravitational lensing. With finite resolution, however, non-zero higher-order cumulants are generated both by lensing autocorrelations and by cross-correlations between the lensing potential and secondary anisotropies in the CMB such as the Sunayev-Zel'dovich (SZ) effect. Skewness is produced by these lensing-SZ correlations, while kurtosis receives contributions from both lensing alone and lensing-SZ correlations. We show that if the projected lensing potential is Gaussian, all cumulants of higher-order than the kurtosis vanish. While recent results raise the possibility of detection of the skewness in upcoming data, the kurtosis will likely remain undetected.

Inflationary gravitational waves (GW) contribute to the curl component in the polarization of the CMB. Cosmic shear--gravitational lensing of the CMB-- converts a fraction of the dominant gradient polarization to the curl component. Higher-order correlations can be used to map the cosmic shear and subtract this contribution to the curl. Arcminute resolution will be required to pursue GW amplitudes smaller than those accessible by Planck. The finite cutoff in CMB power at small scales leads to a minimum detectable GW amplitude corresponding to an inflation energy near 10^15 GeV.

We discuss the four-point correlation function, or the trispectrum in Fourier space, of CMB temperature and polarization anisotropies due to the weak gravitational lensing effect by intervening large scale structure. We discuss the squared temperature power spectrum as a probe of this trispectrum and, more importantly, as an observational approach to extracting the power spectrum of the deflection angle associated with the weak gravitational lensing effect on the CMB. We extend previous discussions on the trispectrum and associated weak lensing reconstruction from CMB data by calculating non-Gaussian noise contributions, beyond the previously discussed dominant Gaussian noise. Non-Gaussian noise contributions are generated by lensing itself and by the correlation between the lensing effect and other foreground secondary anisotropies in the CMB such as the Sunyaev-Zel'dovich (SZ) effect. When the SZ effect is removed from temperature maps using its spectral dependence, we find these additional non-Gaussian noise contributions to be an order of magnitude lower than the dominant Gaussian noise. If the noise-bias due to the dominant Gaussian part of the temperature squared power spectrum is removed, then these additional non-Gaussian contributions provide the limiting noise level for the lensing reconstruction. The temperature squared power spectrum allows a high signal-to-noise extraction of the lensing deflections and a confusion-free separation of the curl (or B-mode) polarization due to inflationary gravitational waves from that due to lensed gradient (or E-mode) polarization. The small angular scale temperature and polarization anisotropy measurements provide a novel approach to weak lensing studies, complementing the approach based on galaxy ellipticities.