The astrophysical origin of gravitational wave (GW) events is one of the most timely problems in the wake of the LIGO/Virgo discoveries. In active galactic nuclei (AGN), binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Previous studies have suggested that binary black hole (BBH) mergers in AGN disks can contribute significantly to BBH mergers observed by GW interferometers. Here we examine the distribution of the effective spin parameter $\chi_\mathrm{eff}$ of this GW source population. We extend our semi-analytical model of binary formation and evolution in AGN disks by following the evolution of the binary orbital angular momenta and black hole (BH) spins. BH spins change due to gas accretion and BH mergers, while the binary orbital angular momenta evolve due to gas accretion and binary-single interactions. We find that the distribution of $\chi_\mathrm{eff}$ predicted by our AGN model is similar to the distribution observed during LIGO/Virgo O1 and O2. On the other hand, if radial migration of BHs is inefficient, $\chi_\mathrm{eff}$ is skewed toward higher values compared with the observed distribution, because of the paucity of scattering events that would randomize spin directions relative to the orbital plane. We suggest that high binary masses and the positive correlation between binary mass and the standard deviation of $\chi_\mathrm{eff}$ for chirp masses up to $\approx 20$ $\mathrm{M}_\odot$, can be possible signatures for mergers originating in AGN disks. Finally, hierarchical mergers in AGN disks naturally produce properties of the recent GW event GW190412, including a low mass ratio, a high primary BH spin, and a significant spin component in the orbital plane.

The astrophysical origin of gravitational wave (GW) events discovered by LIGO/VIRGO remains an outstanding puzzle. In active galactic nuclei (AGN), compact-object binaries form, evolve, and interact with a dense star cluster and a gas disk. An important question is whether and how binaries merge in these environments. To address this question, we have performed one-dimensional $N$-body simulations combined with a semi-analytical model which includes the formation, disruption, and evolution of binaries self-consistently. We point out that binaries can form in single-single interactions by the dissipation of kinetic energy in a gaseous medium. This ``gas capture'' binary formation channel contributes up to $97\,\%$ of gas-driven mergers and leads to a high merger rate in AGN disks even without pre-existing binaries. We find the merger rate to be in the range $\sim 0.02-60\,\mathrm{Gpc^{-3}yr^{-1}}$. The results are insensitive to the assumptions on gaseous hardening processes: we find that once they are formed, binaries merge efficiently via binary-single interactions even if these gaseous processes are neglected. We find that the average number of mergers per BH is $0.4$, and the probability for repeated mergers in 30 Myr is $\sim 0.21-0.45$. High BH masses due to repeated mergers, high eccentricities, and a significant Doppler drift of GWs are promising signatures which distinguish this merger channel from others. Furthermore, we find that gas-capture binaries reproduce the distribution of LMXBs in the Galactic center, including an outer cutoff at $\sim1$ pc due to the competition between migration and hardening by gas torques.

A third object in the vicinity of a binary system causes variations in the eccentricity and the inclination of the binary through the Kozai-Lidov effect. We examine if such variations leave a detectable imprint on the gravitational waves of a binary consisting of intermediate mass black holes and stellar mass objects. As a proof of concept, we present an example where LISA may detect the Kozai-Lidov modulated gravitational wave signals of such sources from at least a distance of 1Mpc if the perturbation is caused by a supermassive black hole tertiary. Although the quick pericenter precession induced by general relativity significantly reduces the appropriate parameter space for this effect by quenching the Kozai-Lidov oscillations, we still find reasonable parameters where the Kozai-Lidov effect may be detected with high signal-to-noise ratios.

Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies.

Approximately two hundred supermassive black holes (SMBHs) have been discovered within the first $\sim$Gyr after the Big Bang. One pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs). A possible obstacle to this scenario is that the collapsing gas fragments and forms a cluster of main-sequence stars. Here we raise the possibility that stellar collisions may be sufficiently frequent and energetic to inhibit the contraction of the massive protostar, avoiding strong UV radiation driven outflows, and allowing it to continue growing into an SMS. We investigate this scenario with semianalytic models incorporating star formation, gas accretion, dynamical friction from stars and gas, stellar collisions, and gas ejection. We find that when the collapsing gas fragments at a density of $\lesssim 3\times 10^{10}\,\mathrm{cm^{-3}}$, the central protostar contracts due to infrequent stellar mergers, and in turn photoevaporates the remaining collapsing gas, resulting in the formation of a $\lesssim 10^4~{\rm M_\odot}$ object. On the other hand, when the collapsing gas fragments at higher densities (expected for a metal-poor cloud with $Z\lesssim10^{-5}\,{\rm Z_\odot}$ with suppressed ${\rm H_2}$ abundance) the central protostar avoids contraction and keeps growing via frequent stellar mergers, reaching masses as high as $\sim 10^5-10^6\,{\rm M_\odot}$. We conclude that frequent stellar mergers represent a possible pathway to form massive BHs in the early universe.

Binary black hole mergers encode information about their environment and the astrophysical processes that led to their formation. Measuring the redshift dependence of their merger rate will help probe the formation and evolution of galaxies and the evolution of the star formation rate. Here we compute the cosmic evolution of the merger rate for stellar-mass binaries in the disks of Active Galactic Nuclei (AGNs). We focus on recent evolution out to redshift $z=2$, covering the accessible range of current Earth-based gravitational-wave observatories. On this scale, the AGN population density is the main contributor to redshift-dependence. We find that the AGN-assisted merger rate does not meaningfully evolve with redshift, differentiating this channel from field binaries and some other dynamical formation scenarios.

Although their existence is not yet confirmed observationally, intermediate mass black holes (IMBHs) may play a key role in the dynamics of galactic nuclei. In this paper, we neglect the effect of the nuclear star cluster itself and investigate only how a small reservoir of IMBHs influences the secular dynamics of stellar-mass black hole binaries, using N-body simulations. We show that our simplifications are valid and that the IMBHs significantly enhance binary evaporation by pushing the binaries into the Hill-unstable region of parameter space, where they are separated by the SMBH's tidal field. For binaries in the S-cluster region of the Milky Way, IMBHs drive the binaries to merge in up to 1-6% of cases, assuming five IMBHs within 5 pc of mass 10,000 solar masses each. Observations of binaries in the Galactic center may strongly constrain the population of IMBHs therein.

Many astrophysical scenarios have been proposed to explain the several black hole (BH) and neutron star binary mergers observed via gravitational waves (GWs) by the LIGO-Virgo collaboration. Contributions from various channels can be statistically disentangled by mass, spin, eccentricity and redshift distributions of merging binaries. In this paper, we investigate the signatures of BH-BH binary mergers induced by a third companion through the Lidov-Kozai mechanism in triple systems. We adopt different prescriptions for the supernovae natal kicks and consider different progenitor metallicities and initial orbital parameters. We show that the typical eccentricity in the LIGO band is $0.01$-$0.1$ and that the merger rate is in the range $0.008-9 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, depending on the natal kick prescriptions and progenitor metallicity. Furthermore, we find that the typical distribution of effective projected spin is peaked at $\chi_{\rm eff}\sim 0$ with significant tails. We show that the triple scenario could reproduce the distribution of $\chi_{\rm eff}$. We find that the triple channel may be strongly constrained by the misalignment angle between the binary component spins in future detections with spin precession.

Despite the rapidly growing number of stellar-mass binary black hole mergers discovered through gravitational waves, the origin of these binaries is still not known. In galactic centers, black holes can be brought to each others' proximity by dynamical processes, resulting in mergers. It is also possible that black holes formed in previous mergers encounter new black holes, resulting in so-called hierarchical mergers. Hierarchical events carry signatures such as higher-than usual black hole mass and spin. Here we show that the recently reported gravitational-wave candidate, GW170817A, could be the result of such a hierarchical merger. In particular, its chirp mass $\sim40$ M$_\odot$ and effective spin of $\chi_{\rm eff}\sim0.5$ are the typically expected values from hierarchical mergers within the disks of active galactic nuclei. We find that the reconstructed parameters of GW170817A strongly favor a hierarchical merger origin over having been produced by an isolated binary origin (with an Odds ratio of $>10^3$, after accounting for differences between the expected rates of hierarchical versus isolated mergers)

We investigate the internal dynamics of anisotropic, rotating globular clusters with a multimass stellar population by performing new direct N-body simulations. In addition to the well-known radial mass segregation effect, where heavy stars and stellar remnants sink toward the center of the cluster, we find a mass segregation in the distribution of orbital inclinations as well. This newly discovered anisotropic mass segregation leads to the formation of a disk-like structure of massive objects near the equatorial plane of a rotating cluster. This result has important implications on the expected spatial distribution of black holes in globular clusters.

We consider the formation of binary intermediate black holes (BIMBH) in globular clusters (GC), which could happen either in situ or due to the mergers between clusters. We simulate the evolution of the BIMBH orbit (and its subsequent merger) due to stellar ejections. We also take into account the evaporation of GCs due to the tidal field of the host galaxy and two-body relaxation. Our results show that if at least $10^{-3}$ of all GCs become BIMBH hosts and the BIMBH masses are $\sim1\%$ of the GC mass, at least one of the inspiralling (or merging) BIMBHs will be detected by LISA during its 4-year mission lifetime. Most of the detected BIMBHs come 1) from heavy GCs ($\gtrsim 3\times10^5M_\odot$), as lower-mass GCs end up being disrupted before their BIMBHs have time to merge, and 2) from redshifts $1<z<3$, assuming that most of GCs form around $z\sim4$ and given that the merger timescale for most BIMBHs is $\sim1$ Gyr. If the BIMBH to GC mass ratio is lower ($\sim10^{-3}$) but the fraction of BIMBH hosts among GCs is higher ($\gtrsim10^{-2}$), some of their mergers will also be detected by LIGO, VIRGO, and KAGRA and the proposed Einstein Telescope.

The origins of the stellar-mass black hole mergers discovered by LIGO/Virgo are still unknown. Here we show that, if migration traps develop in the \add{accretion} disks of Active Galactic Nuclei (AGNs) and promote the mergers of their captive black holes, the majority of black holes within disks will undergo hierarchical mergers---with one of the black holes being the remnant of a previous merger. 40% of AGN-assisted mergers detected by LIGO/Virgo will include a black hole with mass $\gtrsim 50$ M$_\odot$, the mass limit from stellar core collapse. Hierarchical mergers at traps in AGNs will exhibit black hole spins (anti-)aligned with the binary's orbital axis, a distinct property from other hierarchical channels. Our results are suggestive, although not definitive (with Odds ratio of $\sim 1$), that LIGO's heaviest merger so far, GW170729, could have originated from this channel.

We consider a black hole (BH) density cusp in a nuclear star cluster (NSC) hosting a supermassive back hole (SMBH) at its center. Assuming the stars and BHs inside the SMBH sphere of influence are mass-segregated, we calculate the number of BHs that sink into this region under the influence of dynamical friction. We find that the total number of BHs increases significantly in this region due to this process for lower mass SMBHs by up to a factor of 5, but there is no increase in the vicinity of the highest mass SMBHs. Due to the high BH number density in the NSC, BH-BH binaries form during close approaches due to GW emission. We update the previous estimate of O'Leary et al. for the rate of such GW capture events by estimating the $\langle n^2\rangle/\langle n\rangle^2$ parameter where $n$ is the number density. We find a BH merger rate for this channel to be in the range $\sim0.01-0.1 \, \mathrm{Gpc^{-3}yr^{-1}}$. The total merger rate is dominated by the smallest galaxies hosting SMBHs, and the number of heaviest BHs in the NSC. It is also exponentially sensitive to the radial number density profile exponent, reaching $>100 \, \mathrm{Gpc^{-3}yr^{-1}}$ when the BH mass function is $m^{-2.3}$ or shallower and the heaviest BH radial number density is close to $r^{-3}$. Even if the rate is much lower than the range constrained by the current LIGO detections, the GW captures around SMBHs can be distinguished by their high eccentricity in the LIGO band.

Stars passing too close to a black hole can produce tidal disruption events (TDEs), when the tidal force across the star exceeds the gravitational force that binds it. TDEs have usually been discussed in relation to massive black holes that reside in the centres of galaxies or lurk in star clusters. In this paper, we investigate the possibility that triple stars hosting a stellar black hole (SBH) may be sources of TDEs. We start from a triple system made up of three main sequence (MS) stars and model the supernova (SN) kick event that led to the production of an inner binary comprised of a SBH. We evolve these triples in isolation with a high precision $N$-body code and study their TDEs as a result of Kozai-Lidov oscillations. We explore a variety of distributions of natal kicks imparted during the SN event, various maximum initial separations for the triples, and different distributions of eccentricities. We show that the main parameter that governs the properties of the SBH-MS binaries which produce a TDE in triples is the mean velocity of the natal kick distribution. Smaller $\sigma$'s lead to larger inner and outer semi-major axes of the systems that undergo a TDE, smaller SBH masses, and longer timescales. We find that the fraction of systems that produce a TDE is roughly independent of the initial conditions, while estimate a TDE rate of $7.3\times 10^{-5}-4.1 \ \mathrm{yr}^{-1}$, depending on the prescriptions adopted for the SBH natal kicks. This rate is almost comparable to the expected TDE rate for massive black holes.

The localization of stellar-mass binary black hole mergers using gravitational waves is critical in understanding the properties of the binaries' host galaxies, observing possible electromagnetic emission from the mergers, or using them as a cosmological distance ladder. The precision of this localization can be substantially increased with prior astrophysical information about the binary system. In particular, constraining the inclination of the binary can reduce the distance uncertainty of the source. Here we present the first realistic set of localizations for binary black hole mergers, including different prior constraints on the binaries' inclinations. We find that prior information on the inclination can reduce the localization volume by a factor of 3. We discuss two astrophysical scenarios of interest: (i) follow-up searches for beamed electromagnetic/neutrino counterparts and (ii) mergers in the accretion disks of active galactic nuclei.

With the hundreds of merging binary black hole (BH) signals expected to be detected by LIGO/Virgo, LISA and other instruments in the next few years, the modeling of astrophysical channels that lead to the formation of compact-object binaries has become of fundamental importance. In this paper, we carry out a systematic statistical study of quadruple BHs consisting of two binaries in orbit around their center of mass, by means of high-precision direct $N$-body simulations including Post-Newtonian (PN) terms up to 2.5PN order. We found that most merging systems have high initial inclinations and the distributions peak at $\sim 90^\circ$ as for triples, but with a more prominent broad distribution tail. We show that BHs merging through this channel have a significant eccentricity in the LIGO band, typically much larger than BHs merging in isolated binaries and in binaries ejected from star clusters, but comparable to that of merging binaries formed via the GW capture scenario in clusters, mergers in hierarchical triples, or BH binaries orbiting intermediate-mass black holes in star clusters. We show that the merger fraction can be up to $\sim 3$--$4\times$ higher for quadruples than for triples. Thus even if the number of quadruples is $20\%$--$25\%$ of the number of triples, the quadruple scenario can represent an important contribution to the events observed by LIGO/VIRGO.

The growing number of stellar-mass binary black hole mergers discovered by Advanced LIGO and Advanced Virgo are starting to constrain the binaries' origin and environment. However, we still lack sufficiently accurate modeling of binary formation channels to obtain strong constraints, or to identify sub-populations. One promising formation mechanism that could result in different black hole properties is binaries merging within the accretion disks of Active Galactic Nuclei (AGN). Here we show that the black holes' orbital alignment with the AGN disks preferentially selects heavier black holes. We carry out Monte Carlo simulations of orbital alignment with AGN disks, and find that AGNs harden the initial black hole mass function. Assuming an initial power law mass distribution $M_{\rm bh}^{-\beta}$, we find that the power law index changes by $\Delta \beta\sim1.3$, resulting in a more top-heavy population of merging black holes. This change is independent of the mass of, and accretion rate onto, the supermassive black hole in the center of the AGN. Our simulations predict an AGN-assisted merger rate of $\sim4$Gpc$^{-3}$yr$^{-1}$. With its hardened mass spectra, the AGN channel could be responsible for $10-50$% of gravitational-wave detections.

Stellar-mass black hole binaries (BHBs) near supermassive black holes (SMBH) in galactic nuclei undergo eccentricity oscillations due to gravitational perturbations from the SMBH. Previous works have shown that this channel can contribute to the overall BHB merger rate detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo Interferometer. Significantly, the SMBH gravitational perturbations on the binary's orbit may produce eccentric BHBs which are expected to be visible using the upcoming Laser Interferometer Space Antenna (LISA) for a large fraction of their lifetime before they merge in the LIGO/Virgo band. For a proof-of-concept, we show that the eccentricity oscillations of these binaries can be detected with LISA for BHBs in the local universe up to a few Mpcs, with observation periods shorter than the mission lifetime, thereby disentangling this merger channel from others. The approach presented here is straightforward to apply to a wide variety of compact object binaries with a tertiary companion.

LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

In a recent paper, we determined the measurement accuracy of physical parameters for eccentric, precessing, non-spinning, inspiraling, stellar-mass black hole - black hole (BH-BH) binaries for the upcoming second-generation LIGO/VIRGO/KAGRA detector network at design sensitivity using the Fisher matrix method. Here we extend that study to a wide range of binary masses including neutron star - neutron star (NS-NS), NS-BH, and BH-BH binaries with BH masses up to $110 \, M_{\odot}$. The measurement error of eccentricity $e_{10 \,\rm Hz}$ at a gravitational-wave (GW) frequency of $10 \, {\rm Hz}$ is in the range $(10^{-4}-10^{-3}) \times (D_L/ 100\,\rm Mpc)$ for NS-NS, NS-BH, and BH-BH binaries at a luminosity distance of $D_L$ if $e_{10 \,\rm Hz} \gtrsim 0.1 $. For events with masses and distances similar to the detected 10 GW transients, we show that nonzero orbital eccentricities may be detected if $0.081 \lesssim e_{10 \,\rm Hz}$. Consequently, the LIGO/VIRGO/KAGRA detector network at design sensitivity will have the capability to distinguish between eccentric waveforms and circular waveforms. In comparison to circular inspirals, we find that the chirp mass measurement precision can improve by up to a factor of $\sim 20$ and $\sim 50-100$ for NS-NS and NS-BH binaries with BH companion masses $\lesssim 40 \, M_{\odot}$, respectively. The identification of eccentric sources may give information on their astrophysical origin; it would indicate merging binaries in triple or higher multiplicity systems or dynamically formed binaries in dense stellar systems such as globular clusters or galactic nuclei.

We investigate the evolution of triple, non-hierarchical, black hole (BH) systems making use of $2.9\times10^4$ 3-body simulations. Varying the mutual orbital inclination, the three BH masses and the inner and outer eccentricities, we show that retrograde, nearly planar configurations lead to a significant shrinkage of the inner binary. We find an universal trend of triple systems, that they tend to evolve toward prograde configurations, Moreover, we demonstrate that the orbital flip, driven by the torque exerted on the inner BH binary (BHB) by the outer BH, leads in general to tighter inner orbits. In some cases, the resulting BHB undergoes coalescence within a Hubble time, releasing gravitational waves (GWs). Frequently, the inner BHB merger occurs after a component swap between one of its components and the outer BH. The mass spectrum of the BHBs that underwent the component exchange differs significantly from the case in which the BHB merge without any swap. A large fraction of merging BHBs with initial separation $1$ AU enter the $10^{-3}-10^{-1}$ Hz frequency band with large eccentricities, thus representing potential LISA sources. Mergers originating from initially tighter BHB ($a\sim 0.01$ AU), instead, have a large probability to have eccentricities above 0.7 in the $1$ Hz band. We find that the mergers' mass distribution in this astrophysical channel maps the original BH binary spectrum. This might have interesting consequences in light of the growing population of BH mergers detected by LIGO.

Mergers of binary black holes on eccentric orbits are among the targets for second-generation ground-based gravitational-wave detectors. These sources may commonly form in galactic nuclei due to gravitational-wave emission during close flyby events of single objects. We determine the distributions of initial orbital parameters for a population of these gravitational-wave sources. Our results show that the initial dimensionless pericenter distance systematically decreases with the binary component masses and the mass of the central supermassive black hole, and its distribution depends sensitively on the highest possible black hole mass in the nuclear star cluster. For a multi-mass black hole population with masses between 5 Msun and 80 Msun, we find that between 43-69% (68-94%) of 30 Msun - 30 Msun (10 Msun - 10 Msun) sources have an eccentricity greater than 0.1 when the gravitational-wave signal reaches 10 Hz, but less than 10% of the sources with binary component masses less than 30 Msun remain eccentric at this level near the last stable orbit (LSO). The eccentricity at LSO is typically between 0.005-0.05 for the lower-mass BHs, and 0.1 - 0.2 for the highest-mass BHs. Thus, due to the limited low-frequency sensitivity, the six currently known quasi-circular LIGO/Virgo sources could still be compatible with this originally highly eccentric source population. However, at the design sensitivity of these instruments, the measurement of the eccentricity and mass distribution of merger events may be a useful diagnostic to identify the fraction of GW sources formed in this channel.

We present a semi-analytic model for self-consistently evolving a population of globular clusters (GCs) in a given host galaxy across cosmic time. We compute the fraction of GCs still hosting intermediate-mass black holes (IMBHs) at a given redshift in early and late type galaxies of different masses and sizes, and the corresponding rate of tidal disruption events (TDEs), both main-sequence (MS) and white dwarf (WD) stars. We find that the integrated TDE rate for the entire GC population can exceed the corresponding rate in a given galactic nucleus and that $\sim 90$% of the TDEs reside in GCs within a maximum radius of $\sim 2-15$ kpc from the host galaxy's center. This suggests that observational efforts designed to identify TDEs should not confine themselves to galactic nuclei alone, but should also consider the outer galactic halo where massive old GCs hosting IMBHs would reside. Indeed, such off-centre TDEs as predicted here may already have been observed. MS TDE rates are more common than WD TDE rates by a factor 30 (100) at $z\leq 0.5$ ($z=2$). We also calculate the rate of IMBH-SBH mergers across cosmic time, finding that the typical IMRI rate at low redshift is of the order of $\sim 0.5-3$ Gpc$^{-3}$ yr$^{-1}$, which becomes as high as $\sim 100$ Gpc$^{-3}$ yr$^{-1}$ near the peak of GC formation. Advanced LIGO combined with VIRGO, KAGRA, ET and LISA will be able to observe the bottom-end and top-end of the IMBH population, respectively.

The high rate of black hole (BH) mergers detected by LIGO/Virgo opened questions on their astrophysical origin. One possibility is the dynamical channel, in which binary formation and hardening is catalyzed by dynamical encounters in globular clusters (GCs). Previous studies have shown that the BH merger rate from the present day GC density in the Universe is lower than the observed rate. In this \textit{Letter}, we study the BH merger rate by accounting for the first time for the evolution of GCs within their host galaxies. The mass in GCs was initially $\sim 8\times$ higher, which decreased to its present value due to evaporation and tidal disruption. Many BH binaries that were ejected long before their merger, originated in GCs that no longer exist. We find that the comoving merger rate in the dynamical channel from GCs varies between $18$ to $35\,{\rm Gpc}^{-3}\,{\rm yr}^{-1}$ between redshift $z=0.5$ to $2$, and the total rate is $1$, $5$, $24$ events per day within $z=0.5$, $1$, and $2$, respectively. The cosmic evolution and disruption of GCs systematically increases the present-day merger rate by a factor $\sim 2$ relative to isolated clusters. Gravitational wave detector networks offer an unique observational probe of the initial number of GC populations and their subsequent evolution across cosmic time.

Gravitational torques among objects orbiting a supermassive black hole drive the rapid reorientation of orbital planes in nuclear star clusters (NSCs), a process known as vector resonant relaxation. In this Letter, we determine the statistical equilibrium of systems with a distribution of masses, semimajor axes, and eccentricities. We average the interaction over the apsidal precession time and construct a Monte Carlo Markov chain method to sample the microcanonical ensemble of the NSC. We examine the case of NSCs formed by 16 episodes of star formation or globular cluster infall. We find that the massive stars and stellar mass black holes form a warped disk, while low mass stars resemble a spherical distribution with a possible net rotation. This explains the origin of the clockwise disk in the Galactic center and predicts a population of black holes (BHs) embedded within this structure. The rate of mergers among massive stars, tidal disruption events of massive stars by BHs, and BH-BH mergers are highly increased in such disks. The first two may explain the origin of the observed G1 and G2 clouds, the latter may be important for gravitational wave detections with LIGO and VIRGO. More generally, black holes are expected to settle in disks in all dense spherical stellar systems assembled by mergers of smaller systems including globular clusters.

Recently several gravitational wave detections have shown evidence for compact object mergers. However, the astrophysical origin of merging binaries is not well understood. Stellar binaries are typically at much larger separations than what is needed for the binaries to merge due to gravitational wave emission, which leads to the so-called final AU problem. In this letter we propose a new channel for mergers of compact object binaries which solves the final AU problem. We examine the binary evolution following gas expansion due to a weak failed supernova explosion, neutrino mass loss, core disturbance, or envelope instability. In such situations the binary is possibly hardened by ambient gas. We investigate the evolution of the binary system after a shock has propagated by performing smoothed particle hydrodynamics simulations. We find that significant binary hardening occurs when the gas mass bound to the binary exceeds that of the compact objects. This mechanism represents a new possibility for the pathway to mergers for gravitational wave events.

Resonant relaxation has been discussed as an efficient process that changes the angular momenta of stars orbiting around a central supermassive black hole due to the fluctuating gravitational field of the stellar cluster. Other spherical stellar systems, such as globular clusters, exhibit a restricted form of this effect where enhanced relaxation rate only occurs in the directions of the angular momentum vectors, but not in their magnitudes; this is called vector resonant relaxation (VRR). To explore this effect, we performed a large set of direct N-body simulations, with up to 512k particles and ~500 dynamical times. Contrasting our simulations with Spitzer-style Monte Carlo simulations, that by design only exhibit 2-body relaxation, we show that the temporal behavior of the angular momentum vectors in $N$-body simulations cannot be explained by 2-body relaxation alone. VRR operates efficiently in globular clusters with $N>10^4$. The fact that VRR operates in globular clusters may open way to use powerful tools in statistical physics for their description. In particular, since the distribution of orbital planes relaxes much more rapidly than the distribution of the magnitude of angular momentum and the radial action, the relaxation process reaches an internal statistical equilibrium in the corresponding part of phase space while the whole cluster is generally out of equilibrium, in a state of quenched disorder. We point out the need to include effects of VRR in Monte Carlo simulations of globular clusters.

The Milky Way centre exhibits an intense flux in the gamma and X-ray bands, whose origin is partly ascribed to the possible presence of a large population of millisecond pulsars (MSPs) and cataclysmic variables (CVs), respectively. However, the number of sources required to generate such an excess is much larger than what is expected from in situ star formation and evolution, opening a series of questions about the formation history of the Galactic nucleus. In this paper we make use of direct $N$-body simulations to investigate whether these sources could have been brought to the Galactic centre by a population of star clusters that underwent orbital decay and formed the Galactic nuclear star cluster (NSC). Our results suggest that the gamma ray emission is compatible with a population of MSPs that were mass segregated in their parent clusters, while the X-ray emission is consistent with a population of CVs born via dynamical interactions in dense star clusters. Combining observations with our modelling, we explore how the observed $\gamma$ ray flux can be related to different NSC formation scenarios. Finally, we show that the high-energy emission coming from the galactic central regions can be used to detect black holes heavier than $10^5\Ms$ in nearby dwarf galaxies.

The recent discovery of gravitational waves has opened new horizons for physics. Current and upcoming missions, such as LIGO, VIRGO, KAGRA, and LISA, promise to shed light on black holes of every size from stellar mass (SBH) sizes up to supermassive black holes which reside in galactic nuclei. The intermediate mass black hole (IMBH) family has not been detected beyond any reasonable doubt neither directly nor indirectly. Recent analyses suggest observational evidence for the presence of IMBHs in the centers of two Galactic globular clusters. In this paper, we investigate the possibility that globular clusters were born with a central IMBH, which undergo repeated merger events with SBHs in the cluster core. By means of a semi-analytical method, we follow the evolution of the primordial cluster population in the galactic potential and the Gravitational Wave (GW) mergers of the binary IMBH-SBH systems. Our models predict $\approx 1000$ IMBHs within $1$ kpc from the Galactic Center. Our results show that the IMBH-SBH merger rate density changes from $\mathcal{R}\approx 1000$ Gpc$^{-3}$ yr$^{-1}$ beyond $z\approx 2$ to $\mathcal{R}\approx 1-10$ Gpc$^{-3}$ yr$^{-1}$ at $z\approx 0$. The rates at low redshifts may be significantly higher if young massive star clusters host IMBHs. The merger rates are dominated by IMBHs with masses between $10^3$ and $10^4\,\mathrm{M}_{\odot}$. Currently there are no LIGO/VIRGO upper limits for GW sources in this mass range, but our results show that at design sensitivity these instruments may detect these IMBH-SBH mergers in the coming years. \textit{LISA} and the Einstein Telescope will be best suited to detect these GW events. The inspirals of IMBH-SBH systems may also generate an unresolved GW background.

The gravitational interaction among bodies orbiting in a spherical potential leads to the rapid relaxation of the orbital planes' distribution, a process called vector resonant relaxation. We examine the statistical equilibrium of this process for a system of bodies with similar semimajor axes and eccentricities. We extend the previous model of Roupas et al. (2017) by accounting for the multipole moments beyond the quadrupole, which dominate the interaction for radially overlapping orbits. Nevertheless, we find no qualitative differences between the behavior of the system with respect to the model restricted to the quadrupole interaction. The equilibrium distribution resembles a counterrotating disk at low temperature and a spherical structure at high temperature. The system exhibits a first order phase transition between the disk and the spherical phase in the canonical ensemble if the total angular momentum is below a critical value. We find that the phase transition erases the high order multipoles, i.e. small-scale structure in angular momentum space, most efficiently. The system admits a maximum entropy and a maximum energy, which lead to the existence of negative temperature equilibria.

Nuclear star clusters around massive black holes are expected to be abundant in stellar mass black holes and black hole binaries. These binaries form a hierarchical triple system with the massive black hole at the center. Gravitational perturbations from the massive black hole can cause high eccentricity excitation. During this process, the eccentricity may approach unity, and the pericenter distance may become sufficiently small that gravitational wave emission drives the binary to merge. In this paper, we consider a simple proof of concept and explore the effect of the eccentric Kozai-Lidov mechanism for unequal mass binaries. We perform a set of Monte Carlo simulations on BH-BH binaries in galactic nuclei with quadrupole and octupole-level secular perturbations, general relativistic precession, and gravitational wave emission. For a nominal number of steady-state BH-BH binaries, our model gives a total merger rate $\sim 1 - 3$$Gpc^{-3} yr^{-1}$, depending on the assumed density profile. Thus, our model potentially competes with other dynamical mechanisms, such as the dynamical formations and mergers of BH binaries in globular clusters or dense nuclear clusters without a massive black hole. We provide predictions for the distributions of these LIGO sources in galactic nuclei.

Collisional relaxation describes the stochastic process with which a self-gravitating system near equilibrium evolves in phase space due to the fluctuating gravitational field of the system. The characteristic timescale of this process is called the relaxation time. In this paper, we highlight the difference between two measures of the relaxation time in globular clusters: (i) the diffusion time with which the isolating integrals of motion (i.e. energy E and angular momentum magnitude L) of individual stars change stochastically and (ii) the asymptotic timescale required for a family of orbits to mix in the cluster. More specifically, the former corresponds to the instantaneous rate of change of a star's E or L, while the latter corresponds to the timescale for the stars to statistically forget their initial conditions. We show that the diffusion timescales of E and L vary systematically around the commonly used half-mass relaxation time in different regions of the cluster by a factor of ~10 and ~100, respectively, for more than 20% of the stars. We define the mixedness of an orbital family at any given time as the correlation coefficient between its E or L probability distribution functions and those of the whole cluster. Using Monte Carlo simulations, we find that mixedness converges asymptotically exponentially with a decay timescale that is ~10 times the half-mass relaxation time.

Mergers of stellar-mass black holes on highly eccentric orbits are among the targets for ground-based gravitational-wave detectors, including LIGO, VIRGO, and KAGRA. These sources may commonly form through gravitational-wave emission in high velocity dispersion systems or through the secular Kozai-Lidov mechanism in triple systems. Gravitational waves carry information about the binaries' orbital parameters and source location. Using the Fisher matrix technique, we determine the measurement accuracy with which the LIGO-VIRGO-KAGRA network could measure the source parameters of eccentric binaries using a matched filtering search of the repeated burst and eccentric inspiral phases of the waveform. We account for general relativistic precession and the evolution of the orbital eccentricity and frequency during the inspiral. We find that the signal-to-noise ratio and the parameter measurement accuracy may be significantly higher for eccentric sources than for circular sources. This increase is sensitive to the initial pericenter distance, the initial eccentricity, and component masses. For instance, compared to a 30 Msun-30 Msun non-spinning circular binary, the chirp mass and sky localization accuracy can improve for an initially highly eccentric binary by a factor of ~129 (38) and ~2 (11) assuming an initial pericenter distance of 20 Mtot (10 Mtot).

Black hole mergers detectable with LIGO can occur in active galactic nucleus (AGN) disks. Here we parameterize the merger rates, the mass spectrum and the spin spectrum of black holes (BH) in AGN disks. The predicted merger rate spans $\sim 10^{-4}-10^{4} \rm{Gpc}^{-1} \rm{yr}^{-1}$, so upper limits from LIGO ($<212\rm{Gpc}^{-1}\rm{yr}^{-1}$) already constrain it. The predicted mass spectrum has the form of a broken power-law consisting of a pre-existing BH powerlaw mass spectrum and a harder powerlaw mass spectrum resulting from mergers. The predicted spin spectrum is multi-peaked with the evolution of retrograde spin BH in the gas disk playing a key role. We outline the large uncertainties in each of these LIGO observables for this channel and we discuss ways in which they can be constrained in the future.

It has been proposed that primordial black holes (PBHs) form binaries in the radiation dominated era. Once formed, some fraction of them may merge within the age of the Universe by gravitational radiation reaction. We investigate the merger rate of the PBH binaries when the PBHs have a distribution of masses around ${\cal O}(10) M_\odot$, which is a generalization of the previous studies where the PBHs are assumed to have the same mass. After deriving a formula for the merger time probability distribution in the PBH mass plane, we evaluate it under two different approximations. We identify a quantity constructed from the mass-distribution of the merger rate density per unit cosmic time and comoving volume $\mathcal{R}(m_1,m_2)$, $\alpha = -{(m_1+m_2)}^2\partial^2 \ln\mathcal{R}/\partial m_1\partial m_2 $, which universally satisfies $0.97 \lesssim \alpha \lesssim 1.05$ for all binary masses independently of the PBH mass function. This result suggests that the measurement of this quantity is useful for testing the PBH scenario.

The advent of time domain astronomy is revolutionizing our understanding of the Universe. Programs such as the Catalina Real-time Transient Survey (CRTS) or the Palomar Transient Factory (PTF) surveyed millions of objects for several years, allowing variability studies on large statistical samples. The inspection of $\approx$250k quasars in CRTS resulted in a catalogue of 111 potentially periodic sources, put forward as supermassive black hole binary (SMBHB) candidates. A similar investigation on PTF data yielded 33 candidates from a sample of $\approx$35k quasars. Working under the SMBHB hypothesis, we compute the implied SMBHB merger rate and we use it to construct the expected gravitational wave background (GWB) at nano-Hz frequencies, probed by pulsar timing arrays (PTAs). After correcting for incompleteness and assuming virial mass estimates, we find that the GWB implied by the CRTS sample exceeds the current most stringent PTA upper limits by almost an order of magnitude. After further correcting for the implicit bias in virial mass measurements, the implied GWB drops significantly but is still in tension with the most stringent PTA upper limits. Similar results hold for the PTF sample. Bayesian model selection shows that the null hypothesis (whereby the candidates are false positives) is preferred over the binary hypothesis at about $2.3\sigma$ and $3.6\sigma$ for the CRTS and PTF samples respectively. Although not decisive, our analysis highlights the potential of PTAs as astrophysical probes of individual SMBHB candidates and indicates that the CRTS and PTF samples are likely contaminated by several false positives.

We assess the contribution of dynamical hardening by direct three-body scattering interactions to the rate of stellar-mass black hole binary (BHB) mergers in galactic nuclei. We derive an analytic model for the single-binary encounter rate in a nucleus with spherical and disk components hosting a super-massive black hole (SMBH). We determine the total number of encounters $N_{\rm GW}$ needed to harden a BHB to the point that inspiral due to gravitational wave emission occurs before the next three-body scattering event. This is done independently for both the spherical and disk components. Using a Monte Carlo approach, we refine our calculations for $N_{\rm GW}$ to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically $N_{\rm GW} \lesssim$ 10. We find two separate regimes for the efficient dynamical hardening of BHBs: (1) spherical star clusters with high central densities, low velocity dispersions and no significant Keplerian component; and (2) migration traps in disks around SMBHs lacking any significant spherical stellar component in the vicinity of the migration trap, which is expected due to effective orbital inclination reduction of any spherical population by the disk. We also find a weak correlation between the ratio of the second-order velocity moment to velocity dispersion in galactic nuclei and the rate of BHB mergers, where this ratio is a proxy for the ratio between the rotation- and dispersion-supported components. Because disks enforce planar interactions that are efficient in hardening BHBs, particularly in migration traps, they have high merger rates that can contribute significantly to the rate of BHB mergers detected by the advanced Laser Interferometer Gravitational-Wave Observatory.

We examine dense self-gravitating stellar systems dominated by a central potential, such as nuclear star clusters hosting a central supermassive black hole. Different dynamical properties of these systems evolve on vastly different timescales. In particular, the orbital-plane orientations are typically driven into internal thermodynamic equilibrium by vector resonant relaxation before the orbital eccentricities or semimajor axes relax. We show that the statistical mechanics of such systems exhibit a striking resemblance to liquid crystals, with analogous ordered-nematic and disordered-isotropic phases. The ordered phase consists of bodies orbiting in a disk in both directions, with the disk thickness depending on temperature, while the disordered phase corresponds to a nearly isotropic distribution of the orbit normals. We show that below a critical value of the total angular momentum, the system undergoes a first-order phase transition between the ordered and disordered phases. At the critical point the phase transition becomes second-order while for higher angular momenta there is a smooth crossover. We also find metastable equilibria containing two identical disks with mutual inclinations between $90^{\circ}$ and $180^\circ$.

The Laser Interferometer Gravitational-Wave Observatory, LIGO, found direct evidence for double black hole binaries emitting gravitational waves. Galactic nuclei are expected to harbor the densest population of stellar-mass black holes. A significant fraction ($\sim30\%$) of these black holes can reside in binaries. We examine the fate of the black hole binaries in active galactic nuclei, which get trapped in the inner region of the accretion disk around the central supermassive black hole. We show that binary black holes can migrate into and then rapidly merge within the disk well within a Salpeter time. The binaries may also accrete a significant amount of gas from the disk, well above the Eddington rate. This could lead to detectable X-ray or gamma-ray emission, but would require hyper-Eddington accretion with a few percent radiative efficiency, comparable to thin disks. We discuss implications for gravitational wave observations and black hole population studies. We estimate that Advanced LIGO may detect $\sim20$ such, gas-induced binary mergers per year.

The dynamical formation of stellar-mass black hole-black hole binaries has long been a promising source of gravitational waves for the Laser Interferometer Gravitational-Wave Observatory (LIGO). Mass segregation, gravitational focusing, and multibody dynamical interactions naturally increase the interaction rate between the most massive black holes in dense stellar systems, eventually leading them to merge. We find that dynamical interactions, particularly three-body binary formation, enhance the merger rate of black hole binaries with total mass M_tot roughly as ~M_tot^beta, with beta >~ 4. We find that this relation holds mostly independently of the initial mass function, but the exact value depends on the degree of mass segregation. The detection rate of such massive black hole binaries is only further enhanced by LIGO's greater sensitivity to massive black hole binaries with M_tot <~ 80 solar masses. We find that for power-law BH mass functions dN/dM ~ M^-alpha with alpha <~ 2, LIGO is most likely to detect black hole binaries with a mass twice that of the maximum initial black hole mass and a mass ratio near one. Repeated mergers of black holes inside the cluster result in about ~5% of mergers being observed between two and three times the maximum initial black hole mass. Using these relations, one may be able to invert the observed distribution to the initial mass function with multiple detections of merging black hole binaries.

Most, if not all, stars in the field are born in binary configurations or higher multiplicity systems. In dense stellar environment such as the Galactic Center (GC), many stars are expected to be in binary configurations as well. These binaries form hierarchical triple body systems, with the massive black hole (MBH) as the third, distant object. The stellar binaries are expected to undergo large amplitude eccentricity and inclination oscillations via the so-called "eccentric Kozai-Lidov" (EKL) mechanism. These eccentricity excitations, combined with post main sequence stellar evolution, can drive the inner stellar binaries to merge. We study the mergers of stellar binaries in the inner 0.1 pc of the GC caused by gravitational perturbations due to the MBH. We run a large set of Monte Carlo simulations that include the secular evolution of the orbits, general relativistic precession, tides, and post-main-sequence stellar evolution. We find that about 13 % of the initial binary population will have merged after a few million years and about 29 % after a few billion years. These expected merged systems represent a new class of objects at the GC and we speculate that they are connected to G2-like objects and the young stellar population.

The Laser Interferometer Gravitational Wave Observatory (LIGO) has recently discovered gravitational waves (GWs) emitted by merging black hole binaries. We examine whether future GW detections may identify triple companions of merging binaries. Such a triple companion causes variations in the GW signal due to (1) the varying path length along the line of sight during the orbit around the center of mass, (2) relativistic beaming, Doppler, and gravitational redshift, (3) the variation of the "light"-travel time in the gravitational field of the triple companion, and (4) secular variations of the orbital elements. We find that the prospects for detecting the triple companion are the highest for low-mass compact object binaries which spend the longest time in the LIGO frequency band. In particular, for merging neutron star binaries, LIGO may detect a white dwarf or M-dwarf perturber at signal to noise ratio of 8, if it is within 0.4 solar radius distance from the binary and the system is within a distance of 100 Mpc. Stellar mass (supermassive) black hole perturbers may be detected at a factor 5x (1000x) larger separations. Such pertubers in orbit around the merging binary emit GWs at frequencies above 1 mHz detectable by the Laser Interferometer Space Antenna (LISA) in coincidence.

The Fermi satellite has recently detected gamma ray emission from the central regions of our Galaxy. This may be evidence for dark matter particles, a major component of the standard cosmological model, annihilating to produce high-energy photons. We show that the observed signal may instead be generated by millisecond pulsars that formed in dense star clusters in the Galactic halo. Most of these clusters were ultimately disrupted by evaporation and gravitational tides, contributing to a spherical bulge of stars and stellar remnants. The gamma ray amplitude, angular distribution, and spectral signatures of this source may be predicted without free parameters, and are in remarkable agreement with the observations. These gamma rays are from fossil remains of dispersed clusters, telling the history of the Galactic bulge.

An enhanced rate of stellar tidal disruption events (TDEs) may be an important characteristic of supermassive black hole (SMBH) binaries at close separations. Here we study the evolution of the distribution of stars around a SMBH binary due to the eccentric Kozai-Lidov (EKL) mechanism, including octupole effects and apsidal precession caused by the stellar mass distribution and general relativity. We identify a region around one of the SMBHs in the binary where the EKL mechanism drives stars to high eccentricities, which ultimately causes the stars to either scatter off the second SMBH or get disrupted. For SMBH masses 10^7 Msun and 10^8 Msun, the TDE rate can reach 10^{-2} yr and deplete a region of the stellar cusp around the secondary SMBH in ~0.5 Myr. As a result, the final geometry of the stellar distribution between 0.01 and 0.1 pc around the secondary SMBH is a torus. These effects may be even more prominent in nuclear stellar clusters hosting a supermassive and an intermediate mass black hole.

Stars bound to a supermassive black hole interact gravitationally. Persistent torques acting between stellar orbits lead to the rapid resonant relaxation of the orbital orientation vectors ("vector" resonant relaxation) and slower relaxation of the eccentricities ("scalar" resonant relaxation), both at rates much faster than two-body or non-resonant relaxation. We describe a new parallel symplectic integrator, N-ring, which follows the dynamical evolution of a cluster of N stars through vector resonant relaxation, by averaging the pairwise interactions over the orbital period and periapsis-precession timescale. We use N-ring to follow the evolution of clusters containing over 10^4 stars for tens of relaxation times. Among other results, we find that the evolution is dominated by torques among stars with radially overlapping orbits, and that resonant relaxation can be modelled as a random walk of the orbit normals on the sphere, with angular step size ranging from 0.5-1 radian. The relaxation rate in a cluster with a fixed number of stars is proportional to the RMS mass of the stars. The RMS torque generated by the cluster stars is reduced below the torque between Kepler orbits due to apsidal precession and declines weakly with the eccentricity of the perturbed orbit. However since the angular momentum of an orbit also decreases with eccentricity, the relaxation rate is approximately eccentricity-independent for e<0.7 and grows rapidly with eccentricity for e>0.8. We quantify the relaxation using the autocorrelation function of the spherical multipole moments; this decays exponentially and the e-folding time may be identified with the vector resonant relaxation timescale.

Quadrupole oscillation modes in stars can resonate with incident gravitational waves (GWs), and grow non-linear at the expense of GW energy. Stars near massive black hole binaries (MBHB) can act as GW-charged batteries, cooling radiatively. Mass-loss from these stars can prompt MBHB accretion at near-Eddington rates. GW opacity is independent of amplitude, so distant resonating stars can eclipse GW sources. Absorption by the Sun of GWs from Galactic white dwarf binaries may be detectable with second-generation space-based GW detectors as a shadow within a complex diffraction pattern.

If intermediate mass black holes (IMBHs) grow efficiently in gas disks around supermassive black holes, their host active galactic nucleus (AGN) disks should exhibit myriad observational signatures. Gap-opening IMBHs in AGN disks can exhibit spectral features and variability analagous to gapped protoplanetary disks. A gap-opening IMBH in the innermost disk imprints ripples and oscillations on the broad Fe K$\alpha$ line which may be detectable with future X-ray missions. A non-gap-opening IMBH will accrete and produce a soft X-ray excess relative to continuum emission. An IMBH on a retrograde orbit in an AGN disk will not open a gap and will generate soft X-rays from a bow-shock 'headwind'. Accreting IMBH in a large cavity can generate ULX-like X-ray luminosities and LINER-like optical line ratios from local ionized gas. We propose that many LINERs house a weakly accreting MBH binary in a large central disk cavity and will be luminous sources of gravitational waves (GW). IMBHs in galactic nuclei may also be detected via intermittent observational signatures including: UV/X-ray flares due to tidal disruption events, asymmetric X-ray intensity distributions as revealed by AGN transits, quasi-periodic oscillations and underluminous Type Ia supernovae. GW emitted during IMBH inspiral and collisions may be detected with eLISA and LIGO, particularly from LINERs. We summarize observational signatures and compare to current data where possible or suggest future observations.

The dynamical evolution of a hierarchical three body system is well characterized by the eccentric Kozai-Lidov mechanism, where the inner orbit can undergo large eccentricity and inclination oscillations. It was shown before that starting with a circular inner orbit, large mutual inclination (40-140 degree) can produce long timescale modulations that drives the eccentricity to extremely large value and can flip the orbit. Here, we demonstrate that starting with an almost coplanar configuration, for eccentric inner and outer orbits, the eccentricity of the inner orbit can still be excited to high values, and the orbit can flip by ~180 degree, rolling over its major axis. The ~180 degree flip criterion and the flip timescale are described by simple analytic expressions that depend on the initial orbital parameters. With tidal dissipation, this mechanism can produce counter-orbiting exo-planetary systems. In addition, we also show that this mechanism has the potential to enhance the tidal disruption or collision rates for different systems.

Black holes orbiting the Super Massive Black Hole (SMBH) Sgr A* in the Milky-way galaxy center (GC) generate gravitational waves. The spectrum, due to stars and black holes, is continuous below 40 nHz while individual BHs within about 200 AU of the central SMBH stick out in the spectrum at higher frequencies. The GWs can be detected by timing radio pulsars within a few parsecs of this region. Future observations with the Square Kilometer Array of such pulsars with sufficient timing accuracy may be sensitive to signals from intermediate mass BHs (IMBH) in a 3 year observation baseline. The recent detection of radio pulsations from the magnetar SGR J1745-29 very near the GC opens up the possibilities of detecting millisecond pulsars (which can be used as probes of the GWs) through lines of sight with only moderate pulse and angular broadening due to scattering.

We study the secular, hierarchical three-body problem to first-order in a post-Newtonian expansion of General Relativity. We expand the first-order post-Newtonian Hamiltonian to leading-order in the ratio of the semi-major axis of the two orbits. In addition to the well-known terms that correspond to the GR precession of the inner and outer orbits, we find a new secular post-Newtonian interaction term that can affect the long-term evolution of the triple. We explore the parameter space for highly inclined and eccentric systems, where the Kozai-Lidov mechanism can produce large-amplitude oscillations in the eccentricities. The standard lore, i.e., that General Relativity effects suppress eccentricity, is only consistent with the parts of phase space where the General Relativity timescales are several orders of magnitude shorter than the secular Newtonian one. In other parts of phase space, however, post-Newtonian corrections combined with the three body ones, can excite eccentricities. In particular, for systems where the General Relativity timescale is comparable to the secular Newtonian timescales, the three-body interactions give rise to a resonant-like eccentricity excitation. Furthermore, for triples with a comparable-mass inner binary, where the eccentric Kozai-Lidov mechanism is suppressed, post-Newtonian corrections can further increase the eccentricity and lead to orbital flips even when the timescale of the former is much longer than the timescale of the secular Kozai-Lidov quadrupole perturbations.

Black holes are the ultimate prisons of the Universe, regions of spacetime where the enormous gravity prohibits matter or even light to escape to infinity. Yet, matter falling toward the black holes may shine spectacularly, generating the strongest source of radiation. These sources provide us with astrophysical laboratories of extreme physical conditions that cannot be realized on Earth. This chapter offers a review of the basic menus for feeding matter onto black holes and discusses their observational implications.

Galactic nuclei are expected to be densely populated with stellar and intermediate mass black holes. Exploring this population will have important consequences for the observation prospects of gravitational waves as well as understanding galactic evolution. The gas cloud G2 currently approaching Sgr A* provides an unprecedented opportunity to probe the black hole and neutron star population of the Galactic nucleus. We examine the possibility of a G2-black hole encounter and its detectability with current X-ray satellites, such as Chandra and NuSTAR. We find that multiple encounters are likely to occur close to the pericenter, which may be detectable upon favorable circumstances. This opportunity provides an additional, important science case for leading X-ray observatories to closely follow G2 on its way to the nucleus.

(abridged) When a sufficiently massive satellite (or secondary) black hole is embedded in a gas disk around a (primary) supermassive black hole, it can open an empty gap in the disk. A gap-opening secondary close to the primary will leave an imprint in the broad component of the FeKa emission line, which varies in a unique and predictable manner. If the gap persists into the innermost disk, the effect consists of a pair of dips in the broad line which ripple blue-ward and red-ward from the line centroid energy respectively, as the gap moves closer to the primary. This ripple effect could be unambiguously detectable and allow an electromagnetic monitoring of massive black hole mergers as they occur. As the mass ratio of the secondary to primary black hole increases to q>0.01, we expect the gap to widen, possibly clearing a central cavity in the inner disk, which shows up in the broad FeKa line component. If the secondary stalls at >100r_{g} in its in-migration, due to low co-rotating gas mass, a detectable ripple effect occurs in the broad line component on the disk viscous timescale as the inner disk drains and the outer disk is dammed. If the secondary maintains an accretion disk within a central cavity, due to dam bursting or leakage, a periodic 'see-saw' oscillation effect is exhibited in the observed line profile. Here we demonstrate the range of ripple effect signatures potentially detectable with Astro-H and IXO/Athena, and oscillation effects potentially detectable with XMM or LOFT for a wide variety of merger and disk conditions. Observations of the ripple effect and periodic oscillations can be used to provide an early warning of gravitational radiation emission from the AGN.

It is commonly assumed that ground-based gravitational wave (GW) instruments will not be sensitive to supermassive black holes (SMBHs) because the characteristic GW frequencies are far below the ~ 10 - 1000 Hz sensitivity bands of terrestrial detectors. Here, however, we explore the possibility of SMBH gravitational waves to leak to higher frequencies. In particular, if the high frequency spectral tail asymptotes to h(f) ~ f^{-alpha}, where alpha<=2, then the spectral amplitude is a constant or increasing function of the mass M at a fixed frequency f>>c^3/GM. This will happen if the time domain waveform or its derivative exhibits a discontinuity. Ground based instruments could search for these universal spectral tails to detect or rule out such features irrespective of their origin. We identify the following processes which may generate high frequency signals: (i) gravitational bremsstrahlung of ultrarelativistic objects in the vicinity of a SMBH, (ii) ringdown modes excited by an external process that has a high frequency component or terminates abruptly, (iii) gravitational lensing echos and diffraction. In particular for (iii), SMBHs produce GW echos of inspiraling stellar mass binaries in galactic nuclei with a delay of a few minutes to hours. The lensed primary signal and GW echo are amplified if the binary is within a ~10 deg (r/100M)^{-1/2} cone behind the SMBH relative to the line of sight at distance r from the SMBH. For the rest of the binaries near SMBHs, the amplitude of the GW echo is ~0.1 (r/100M)^{-1} of the primary signal on average.

Inspiraling supermassive black hole binary systems with high orbital eccentricity are important sources for space-based gravitational wave (GW) observatories like the Laser Interferometer Space Antenna (LISA). Eccentricity adds orbital harmonics to the Fourier transform of the GW signal and relativistic pericenter precession leads to a three-way splitting of each harmonic peak. We study the parameter estimation accuracy for such waveforms with different initial eccentricity using the Fisher matrix method and a Monte Carlo sampling of the initial binary orientation. The eccentricity improves the parameter estimation by breaking degeneracies between different parameters. In particular, we find that the source localization precision improves significantly for higher-mass binaries due to eccentricity. The typical sky position errors are $\sim1 $deg for a nonspinning, $10^7\,M_{\odot}$ equal-mass binary at redshift $z=1$, if the initial eccentricity 1 yr before merger is $e_0\sim 0.6$. Pericenter precession does not affect the source localization accuracy significantly, but it does further improve the mass and eccentricity estimation accuracy systematically by a factor of 3--10 for masses between $10^6$ and $10^7\,M_{\odot}$ for $e_0 \sim 0.3$.

Supermassive black holes (SMBH) are typically surrounded by a dense stellar population in galactic nuclei. Stars crossing the line of site in active galactic nuclei (AGN) produce a characteristic transit lightcurve, just like extrasolar planets do when they transit their host star. We examine the possibility of finding such AGN transits in deep optical, UV, and X-ray surveys. We calculate transit lightcurves using the Novikov--Thorne thin accretion disk model, including general relatistic effects. Based on the expected properties of stellar cusps, we find that around 10^6 solar mass SMBHs, transits of red giants are most common for stars on close orbits with transit durations of a few weeks and orbital periods of a few years. We find that detecting AGN transits requires repeated observations of thousands of low mass AGNs to 1% photometric accuracy in optical, or ~ 10% in UV bands or soft X-ray. It may be possible to identify stellar transits in the Pan-STARRS and LSST optical and the eROSITA X-ray surveys. Such observations could be used to constrain black hole mass, spin, inclination and accretion rate. Transit rates and durations could give valuable information on the circumnuclear stellar clusters as well. Transit lightcurves could be used to image accretion disks with unprecedented resolution, allowing to resolve the SMBH silhouette in distant AGNs.

Many astrophysical binaries, from planets to black holes, exert strong torques on their circumbinary accretion disks, and are expected to significantly modify the disk structure. Despite the several decade long history of the subject, the joint evolution of the binary + disk system has not been modeled with self-consistent assumptions for arbitrary mass ratios and accretion rates. Here we solve the coupled binary-disk evolution equations analytically in the strongly perturbed limit, treating the azimuthally-averaged angular momentum exchange between the disk and the binary and the modifications to the density, scale-height, and viscosity self-consistently, including viscous and tidal heating, diffusion limited cooling, radiation pressure, and the orbital decay of the binary. We find a solution with a central cavity and a migration rate similar to those previously obtained for Type-II migration, applicable for large masses and binary separations, and near-equal mass ratios. However, we identify a distinct new regime, applicable at smaller separations and masses, and mass ratio in the range 0.001< q < 0.1. For these systems, gas piles up outside the binary's orbit, but rather than creating a cavity, it continuously overflows as in a porous dam. The disk profile is intermediate between a weakly perturbed disk (producing Type-I migration) and a disk with a gap (with Type-II migration). However, the migration rate of the secondary is typically slower than both Type-I and Type-II rates. We term this new regime "Type-1.5" migration.

We study the interaction of a supermassive black hole (SMBH) binary and a standard radiatively efficient thin accretion disk. We examine steady-state configurations of the disk and migrating SMBH system, self-consistently accounting for tidal and viscous torques and heating, radiative diffusion limited cooling, gas and radiation pressure, and the decay of the binary's orbit. We obtain a "phase diagram" of the system as a function of binary parameters, showing regimes in which both the disk structure and migration have a different character. Although massive binaries can create a central gap in the disk at large radii, the tidal barrier of the secondary causes a significant pile-up of gas outside of its orbit, which can lead to the closing of the gap. We find that this spillover occurs at an orbital separation as large as ~200 M_7^{-1/2} gravitational radii, where M = 10^7 M_7 Msun is the total binary mass. If the secondary is less massive than ~10^6 Msun, then the gap is closed before gravitational waves (GWs) start dominating the orbital decay. In this regime, the disk is still strongly perturbed, but the piled-up gas continuously overflows as in a porous dam, and crosses inside the secondary's orbit. The corresponding migration rate, which we label Type 1.5, is slower than the usual limiting cases known as Type I and II migration. Compared to an unperturbed disk, the steady-state disk in the overflowing regime is up to several hundred times brighter in the optical bands. Surveys such as PanSTARRS or LSST may discover the periodic variability of this population of binaries. Our results imply that the circumbinary disks around SMBHs can extend to small radii during the last stages of their merger, when they are detectable by LISA, and may produce coincident electromagnetic (EM) emission similar to active galactic nuclei (AGN).

Galactic nuclei are densely populated by stellar mass compact objects such as black holes and neutron stars. Bound, highly eccentric binaries form as a result of gravitational wave (GW) losses during close flybys between these objects. We study the evolution of these systems using 2.5 and 3.5 order post-Newtonian equations of motion. The GW signal consists of many thousand repeated bursts (RB) for minutes to days (depending on the impact parameter and masses), followed by a powerful GW chirp and an eccentric merger. We show that a significant signal to noise ratio (SNR) accumulates already in the RB phase, corresponding to a detection limit around 200-300 Mpc and 300--600 Mpc for Advanced LIGO for an average orientation BH/NS or BH/BH binary, respectively. The theoretical errors introduced by the inaccuracy of the PN templates are typically much less severe for the RB phase than in the following eccentric merger. The GW signal in the RB phase is broadband; we show that encounters involving intermediate mass black holes are detectable in multiple frequency bands coincidentally using LIGO and LISA.

We investigate the electromagnetic (EM) counterpart of gravitational waves (GWs) emitted by a supermassive black hole binary (SMBHB) through the viscous dissipation of the GW energy in an accretion disk and stars surrounding the SMBHB. We account for the suppression of the heating rate if the forcing period is shorter than the turnover time of the largest turbulent eddies. We find that the viscous heating luminosity in 0.1 solar mass stars can be significantly higher than their intrinsic luminosity, but still too low to be detected for extragalactic sources. The relative brightening is small for accretion disks.

We examine the nHz gravitational wave (GW) foreground of stars and black holes (BHs) orbiting SgrA* in the Galactic Center. A cusp of stars and BHs generates a continuous GW spectrum below 40 nHz; individual BHs within 1 mpc to SgrA* stick out in the spectrum at higher GW frequencies. The GWs and gravitational near-field effects can be resolved by timing pulsars within a few pc of this region. Observations with the Square Kilometer Array (SKA) may be especially sensitive to intermediate mass black holes (IMBHs) in this region, if present. A 100ns-10mus timing accuracy is sufficient to detect BHs of mass 1000 Msun with pulsars at distance 0.1-1 pc in a 3 yr observation baseline. Unlike electromagnetic imaging techniques, the prospects for resolving individual objects through GW measurements improve closer to SgrA*, even if the number density of objects increases inwards steeply. Scattering by the interstellar medium will pose the biggest challenge for such observations.

We examine the electromagnetic (EM) and gravitational wave (GW) signatures of stellar-mass compact objects (COs) spiraling into a supermassive black hole (extreme mass-ratio inspirals or EMRIs), embedded in a thin, radiation-pressure dominated, accretion disk. At large separations, the tidal effect of the secondary CO clears a gap. We show that the gap refills during the late GW-driven phase of the inspiral, leading to a sudden EM brightening of the source. The accretion disk leaves an imprint on the GW through its angular momentum exchange with the binary, the mass increase of the binary members due to accretion, and its gravity. We compute the disk-modified GWs both in an analytical Newtonian approximation and in a numerical effective-one-body approach. We find that disk-induced migration provides the dominant perturbation to the inspiral, with weaker effects from the mass accretion onto the CO and hydrodynamic drag. Depending on whether a gap is present, the perturbation of the GW phase is between 10 and 1000 radians per year, detectable with the future Laser Interferometer Space Antenna (LISA) at high significance. The Fourier transform of the disk-modified GW in the stationary phase approximation is sensitive to disk parameters with a frequency trend different from post-Newtonian vacuum corrections. Our results suggest that observations of EMRIs may place new sensitive constraints on the physics of accretion disks.

We study the effects of a thin gaseous accretion disk on the inspiral of a stellar--mass black hole into a supermassive black hole. We construct a phenomenological angular momentum transport equation that reproduces known disk effects. Disk torques modify the gravitational wave phase evolution to detectable levels with LISA for reasonable disk parameters. The Fourier transform of disk-modified waveforms acquires a correction with a different frequency trend than post-Newtonian vacuum terms. Such inspirals could be used to detect accretion disks with LISA and to probe their physical parameters.

Observations of the spatial distribution and kinematics of young stars in the Galactic centre can be interpreted as showing that the stars occupy one, or possibly two, discs of radii ~0.05-0.5 pc. The most prominent (`clockwise') disc exhibits a strong warp: the normals to the mean orbital planes in the inner and outer third of the disc differ by ~60 deg. Using an analytical model based on Laplace-Lagrange theory, we show that such warps arise naturally and inevitably through vector resonant relaxation between the disc and the surrounding old stellar cluster.

Pulsar timing arrays (PTAs) measure nHz frequency gravitational waves (GWs) generated by orbiting massive black hole binaries (MBHBs) with periods between 0.1-10 yr. Previous studies on the nHz GW background assumed that the inspiral is purely driven by GWs. However, torques generated by a gaseous disk can shrink the binary much more efficiently than GW emission, reducing the number of binaries at these separations. We use simple disk models for the circumbinary gas and for the binary-disk interaction to follow the orbital decay of MBHBs through physically distinct regions of the disk, until GWs take over their evolution. We extract MBHB cosmological merger rates from the Millennium simulation, generate Monte Carlo realizations of a population of gas driven binaries, and calculate the corresponding GW amplitudes of the most luminous individual binaries and the stochastic GW background. For steady state alpha-disks with alpha>0.1 we find that the nHz GW background can be significantly modified. The number of resolvable binaries is however not changed by the presence of gas; we predict 1-10 individually resolvable sources to stand above the noise for a 1-50 ns timing precision. Gas driven migration reduces predominantly the number of small total mass or unequal mass ratio binaries, which leads to the attenuation of the mean stochastic GW--background, but increases the detection significance of individually resolvable binaries. The results are sensitive to the model of binary--disk interaction. The GW background is not attenuated significantly for time-dependent models of Ivanov, Papaloizou, & Polnarev (1999).

In order to attract and retain excellent researchers and diverse individuals in astrophysics, we recommend action be taken in several key areas impacting young scientists: (1) Maintain balance between large collaborations and individual projects through distribution of funding; encourage public releases of observational and simulation data for use by a broader community. (2) Improve the involvement of women, particularly at leading institutions. (3) Address the critical shortage of child care options and design reasonable profession-wide parental leave policies. (4) Streamline the job application and hiring process. We summarize our reasons for bringing these areas to the attention of the committee, and we suggest several practical steps that can be taken to address them.

Supermassive black hole binaries (SMBHBs) in galactic nuclei are thought to be a common by-product of major galaxy mergers. We use simple disk models for the circumbinary gas and for the binary-disk interaction to follow the orbital decay of SMBHBs with a range of total masses (M) and mass ratios (q), through physically distinct regions of the disk, until gravitational waves (GWs) take over their evolution. Prior to the GW-driven phase, the viscous decay is in the stalled "secondary-dominated" regime. SMBHBs spend a non-negligible fraction of 10^7 years at orbital periods t_var between a day and a year. A dedicated optical or X-ray survey could identify coalescing SMBHBs statistically, as a population of periodically variable quasars, whose abundance N_var is proportional to t_var^alpha, in a range of periods t_var around tens of weeks. SMBHBs with M < 10^7 M_sun, with 0.5 < alpha < 1.5, would probe the physics of viscous orbital decay, whereas the detection of a population of higher-mass binaries, with alpha=8/3, would confirm that their decay is driven by GWs. The lowest mass SMBHBs (M < 10^{5-6} M_sun) enter the GW-driven regime at short orbital periods, in the frequency band of the Laser Interferometric Space Antenna (LISA). While viscous processes are strongly sub-dominant in the last few years of coalescence, they could reduce the amplitude of any unresolved background of near-stationary LISA sources. We discuss constraints on the SMBHB population available from existing data, and the sensitivity and sky coverage requirements for a detection in future surveys. SMBHBs may also be identified from velocity shifts in their spectra; we discuss the expected abundance of SMBHBs as a function of their orbital velocity.

Supermassive black hole binaries (SMBHBs) produced in galaxy mergers are thought to complete their coalescence, below separations of r_GW=10^{-3} (M_BH/10^8 M_sun)^{3/4} pc, as their orbit decays due to the emission of gravitational waves (GWs). It may be possible to identify such GW-driven inspirals statistically in an electromagnetic (EM) survey for variable sources. A GW-driven binary spends a characteristic time T_GW at each orbital separation r_orb < r_GW that scales with the corresponding orbital time t_orb as T_GW = (const) t_orb^{8/3}. If the coalescing binary produces variations in the EM emission on this timescale, then it could be identified as a variable source with a characteristic period t_var = t_orb. The incidence rate of sources with similar inferred BH masses, showing near-periodic variability on the time-scale t_var, would then be proportional to t_var^{8/3}. Luminosity variations corresponding to a fraction f_Edd<0.01 of the Eddington luminosity would have been missed in current surveys. However, if the binary inspirals are associated with quasars, we show that a dedicated survey could detect the population of SMBHBs with a range of periods around tens of weeks. The discovery of a population of periodic sources whose abundance obeys N_var = (const) t_var^{8/3} would confirm (i) that the orbital decay is indeed driven by GWs, and (ii) that circumbinary gas is present at small orbital radii and is being perturbed by the BHs. Deviations from the t_var^{8/3} power-law could constrain the structure of the circumbinary gas disk and viscosity-driven orbital decay. We discuss constraints from existing data, and quantify the sensitivity and sky coverage that could yield a detection in future surveys.

Stellar mass black holes (BHs) are expected to segregate and form a steep density cusp around supermassive black holes (SMBHs) in galactic nuclei. We follow the evolution of a multi-mass system of BHs and stars by numerically integrating the Fokker-Planck energy diffusion equations for a variety of BH mass distributions. We find that the BHs "self-segregate'', and that the rarest, most massive BHs dominate the scattering rate closest to the SMBH (< 0.1 pc). BH--BH binaries form out of gravitational wave emission during BH encounters. We find that the expected rate of BH coalescence events detectable by Advanced LIGO is ~1 - 100/yr, depending on the initial mass function of stars in galactic nuclei and the mass of the most massive BHs. We find that the actual merger rate is likely ~10 times larger than this due to the intrinsic scatter of stellar densities in many different galaxies. The BH binaries that form this way in galactic nuclei have significant eccentricities as they enter the LIGO band (90% with e > 0.9), and are therefore distinguishable from other binaries, which circularise before becoming detectable. We also show that eccentric mergers can be detected to larger distances and greater BH masses than circular mergers, up to ~ 700 solar masses. Future ground-based gravitational wave observatories will be able to constrain both the mass function of BHs and stars in galactic nuclei.

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

Supermassive black hole binaries (SMBHBs) with masses in the range 10^4-10^7 M_sun/(1+z), produced in galaxy mergers, are thought to complete their coalescence due to the emission of gravitational waves (GWs). The anticipated detection of the GWs by the LISA will constitute a milestone for fundamental physics and astrophysics. While the GW signatures themselves will provide a treasure trove of information, if the source can be securely identified in electromagnetic (EM) bands, this would open up entirely new scientific opportunities, to probe fundamental physics, astrophysics, and cosmology. We discuss several ideas, involving wide-field telescopes, that may be useful in locating electromagnetic counterparts to SMBHBs detected by LISA. In particular, the binary may produce a variable electromagnetic flux, such as a roughly periodic signal due to the orbital motion prior to coalescence, or a prompt transient signal caused by shocks in the circumbinary disk when the SMBHB recoils and "shakes" the disk. We discuss whether these time-variable EM signatures may be detectable, and how they can help in identifying a unique counterpart within the localization errors provided by LISA. We also discuss a possibility of identifying a population of coalescing SMBHBs statistically, in a deep optical survey for periodically variable sources, before LISA detects the GWs directly. The discovery of such sources would confirm that gas is present in the vicinity and is being perturbed by the SMBHB - serving as a proof of concept for eventually finding actual LISA counterparts.

Mergers of supermassive black hole binaries release peak power of up to ~10^57 erg/s in gravitational waves (GWs). As the GWs propagate through ambient gas, they induce shear and a small fraction of their power is dissipated through viscosity. The dissipated heat appears as electromagnetic (EM) radiation, providing a prompt EM counterpart to the GW signal. For thin accretion disks, the GW heating rate exceeds the accretion power at distances farther than 10^3 Schwarzschild radii, independently of the accretion rate and viscosity coefficient.

Transiting exoplanetary systems are surpassingly important among the planetary systems since they provide the widest spectrum of information for both the planet and the host star. If a transiting planet is on an eccentric orbit, the duration of transits T_D is sensitive to the orientation of the orbital ellipse relative to the line of sight. The precession of the orbit results in a systematic variation in both the duration of individual transit events and the observed period between successive transits, P_obs. The periastron of the ellipse slowly precesses due to general relativity and possibly the presence of other planets in the system. This secular precession can be detected through the long-term change in P_obs (transit timing variations, TTV) or in T_D (transit duration variations, TDV). We estimate the corresponding precession measurement precision for repeated future observations of the known eccentric transiting exoplanetary systems (XO-3b, HD 147506b, GJ 436b and HD 17156b) using existing or planned space-borne instruments. The TDV measurement improves the precession detection sensitivity by orders of magnitude over the TTV measurement. We find that TDV measurements over a ~4 year period can typically detect the precession rate to a precision well exceeding the level predicted by general relativity.

The notion that microparsec-scale black holes can be used to probe gigaparsec-scale physics may seem counterintuitive, at first. Yet, the gravitational observatory LISA will detect cosmologically-distant coalescing pairs of massive black holes, accurately measure their luminosity distance and help identify an electromagnetic counterpart or a host galaxy. A wide variety of new black hole studies and a gravitational version of Hubble's diagram become possible if host galaxies are successfully identified. Furthermore, if dark energy is a manifestation of large-scale modified gravity, deviations from general relativistic expectations could become apparent in a gravitational signal propagated over cosmological scales, especially when compared to the electromagnetic signal from a same source. Finally, since inspirals of white dwarfs into massive black holes at cosmological distances may permit pre-merger localizations, we suggest that careful monitoring of these events and any associated electromagnetic counterpart could lead to high-precision cosmological measurements with LISA.

Electromagnetic (EM) counterparts to supermassive black hole binary mergers observed by LISA can be localized to within the field of view of astronomical instruments ~10 deg^2 hours to weeks prior to coalescence. The temporal coincidence of any prompt EM counterpart with a gravitationally-timed merger may offer the best chance of identifying a unique host galaxy. We discuss the challenges posed by searches for prompt EM counterparts and propose novel observational strategies to address them. In particular, we discuss the size and shape evolution of the LISA localization error ellipses on the sky, and quantify the requirements for dedicated EM surveys of the area prior to coalescence. A triggered EM counterpart search campaign will require monitoring a several-square degree area. It could aim for variability at the 24-27 mag level in optical bands, for example, which corresponds to 1-10% of the Eddington luminosity of the prime LISA sources of 10^6-10^7 Msun BHs at z=1-2, on time-scales of minutes to hours, the orbital time-scale of the binary in the last 2-4 weeks. A cross-correlation of the period of any variable EM signal with the quasi-periodic gravitational waveform over 10-1000 cycles may aid the detection. Alternatively, EM searches can detect a transient signal accompanying the coalescence. We highlight the measurement of differences in the arrival times of photons and gravitons from the same cosmological source as a valuable independent test of the massive character of gravity, and of possible violations of Lorentz invariance in the gravity sector.

When a source emits a gravity-wave (GW) pulse over a short period of time, the leading edge of the GW signal is redshifted more than the inner boundary of the pulse. The GW pulse is distorted by the gravitational effect of the self-energy residing in between these shells. We illustrate this distortion for GW pulses from the final plunge of black hole (BH) binaries, leading to the evolution of the GW profile as a function of the radial distance from the source. The distortion depends on the total GW energy released and the duration of the emission, scaled by the total binary mass, M. The effect should be relevant in finite box simulations where the waveforms are extracted within a radius of <~ 100M. For characteristic emission parameters at the final plunge between binary BHs of arbitrary spins, this effect could distort the simulated GW templates for LIGO and LISA by a fraction of 0.001. Accounting for the wave distortion would significantly decrease the waveform extraction errors in numerical simulations.

The continuous improvement in localization errors (sky position and distance) in real time as LISA observes the gradual inspiral of a supermassive black hole (SMBH) binary can be of great help in identifying any prompt electromagnetic counterpart associated with the merger. We develop a new method, based on a Fourier decomposition of the time-dependent, LISA-modulated gravitational-wave signal, to study this intricate problem. The method is faster than standard Monte Carlo simulations by orders of magnitude. By surveying the parameter space of potential LISA sources, we find that counterparts to SMBH binary mergers with total mass M~10^5-10^7 M_Sun and redshifts z<~3 can be localized to within the field of view of astronomical instruments (~deg^2) typically hours to weeks prior to coalescence. This will allow targeted searches for variable electromagnetic counterparts as the merger proceeds, as well as monitoring of the most energetic coalescence phase. A rich set of astrophysical and cosmological applications would emerge from the identification of electromagnetic counterparts to these gravitational-wave standard sirens.

The rapid advance of gravitational-wave (GW) detector facilities makes it very important to estimate the event rates of possible detection candidates. We consider an additional possibility of GW bursts produced during parabolic encounters (PEs) of stellar mass compact objects. We estimate the rate of successful detections for specific detectors: the initial Laser Interferometric Gravitational-Wave Observatory (InLIGO), the French-Italian gravitational-wave antenna VIRGO, the near-future Advanced-LIGO (AdLIGO), the space-based Laser Interferometric Space Antenna (LISA), and the Next Generation LISA (NGLISA). Simple GC models are constructed to account for the compact object mass function, mass segregation, number density distribution, and velocity distribution. We calculate encounters both classically and account for general relativistic corrections by extrapolating the results for infinite mass ratios. We also include the cosmological redshift of waveforms and event rates. We find that typical PEs with masses m_1=m_2=40 Msun are detectable with matched filtering over a signal to noise ratio of 5 within a distance d_L~200Mpc for InLIGO and VIRGO, z=1 for AdLIGO, 0.4Mpc for LISA, and 1Gpc for NGLISA. We estimate single datastream total detection rates of 5.5 x 10^{-5} for InLIGO, 7.2 x 10^{-5} for VIRGO, 0.063 for AdLIGO, 2.9 x 10^{-6} for LISA, and 1.0 for NGLISA per year, for reasonably conservative assumptions. These estimates are subject to uncertainties in the GC parameters, most importantly the total number and mass-distribution of black holes (BHs) in the cluster core. In reasonably optimistic cases, we get >~1 detections for AdLIGO per year. The regular detection of GWs during PEs would provide a unique observational probe for constraining the stellar BH mass function of dense clusters. (abridged)

The gravitational waves (GW) emitted during the coalescence of supermassive black holes (SMBHs) in the mass range 10^4-10^7 M_sun will be detectable out to high redshifts with LISA. We calculate the size and orientation of the three-dimensional error ellipse in solid angle and redshift within which the LISA event could be localized using the GW signatures alone. We take into account uncertainties in LISA's measurements of the luminosity distance and direction to the source, in the background cosmology, in weak gravitational lensing magnification due to inhomogeneities along the line of sight, and potential source peculiar velocities. We find that weak lensing errors exceed other sources of uncertainties by nearly an order of magnitude. Under the plausible assumption that BH mergers are accompanied by gas accretion leading to Eddington-limited quasar activity, we then compute the number of quasars that would be found in a typical LISA error volume, as a function of BH mass and redshift. We find that low redshifts offer the best opportunities to identify quasar counterparts to cosmological standard sirens, and that the LISA error volume will typically contain a single near-Eddington quasar at z=1. This will allow a straightforward test of the hypothesis that BH mergers are accompanied by bright quasar activity and, if the hypothesis proves correct, will guarantee the identification of a unique quasar counterpart. This would yield unprecedented tests of the physics of SMBH accretion, and offer an alternative method to precisely constrain cosmological parameters [abridged].

In cosmological structure formation models, massive non-linear objects in the process of formation, such as galaxy clusters, are surrounded by large-scale shocks at or around the expected virial radius. Direct observational evidence for such virial shocks is currently lacking, but we show here that their presence can be inferred from future, high resolution, high-sensitivity observations of the Sunyaev-Zel'dovich (SZ) effect in galaxy clusters. We study the detectability of virial shocks in mock SZ maps, using simple models of cluster structure (gas density and temperature distributions) and noise (background and foreground galaxy clusters projected along the line of sight, as well as the cosmic microwave background anisotropies). We find that at an angular resolution of 2'' and sensitivity of 10 micro K, expected to be reached at ~ 100 GHz frequencies in a ~ 20 hr integration with the forthcoming ALMA instrument, virial shocks associated with massive M ~ 10^15 M_Sun clusters will stand out from the noise, and can be detected at high significance. More generally, our results imply that the projected SZ surface brightness profile in future, high-resolution experiments will provide sensitive constraints on the density profile of cluster gas.