We report the first plausible optical electromagnetic (EM) counterpart to a (candidate) binary black hole (BBH) merger. Detected by the Zwicky Transient Facility (ZTF), the EM flare is consistent with expectations for a kicked BBH merger in the accretion disk of an active galactic nucleus (AGN), and is unlikely ($<O(0.01\%$)) due to intrinsic variability of this source. The lack of color evolution implies that it is not a supernovae and instead is strongly suggestive of a constant temperature shock. Other false-positive events, such as microlensing or a tidal disruption event, are ruled out or constrained to be $<O(0.1\%$). If the flare is associated with S190521g, we find plausible values of: total mass $ M_{\rm BBH} \sim 100 M_{\odot}$, kick velocity $v_k \sim 200\, {\rm km}\, {\rm s}^{-1}$ at $\theta \sim 60^{\circ}$ in a disk with aspect ratio $H/a \sim 0.01$ (i.e., disk height $H$ at radius $a$) and gas density $\rho \sim 10^{-10}\, {\rm g}\, {\rm cm}^{-3}$. The merger could have occurred at a disk migration trap ($a \sim 700\, r_{g}$; $r_g \equiv G M_{\rm SMBH} / c^2$, where $M_{\rm SMBH}$ is the mass of the AGN supermassive black hole). The combination of parameters implies a significant spin for at least one of the black holes in S190521g. The timing of our spectroscopy prevents useful constraints on broad-line asymmetry due to an off-center flare. We predict a repeat flare in this source due to a re-encountering with the disk in $\sim 1.6\, {\rm yr}\, (M_{\rm SMBH}/10^{8}M_{\odot})\, (a/10^{3}r_{g})^{3/2}$.

Active galactic nuclei (AGN) are powered by the accretion of disks of gas onto supermassive black holes (SMBHs). Stars and stellar remnants orbiting the SMBH in the nuclear star cluster (NSC) will interact with the AGN disk. Orbiters plunging through the disk experience a drag force and, through repeated passage, can have their orbits captured by the disk. A population of embedded objects in AGN disks may be a significant source of binary black hole mergers, supernovae, tidal disruption events and embedded gamma-ray bursts. For two representative AGN disk models we use geometric drag and Bondi-Hoyle-Littleton drag to determine the time to capture for stars and stellar remnants. We assume a range of initial inclination angles and semi-major axes for circular Keplerian prograde orbiters. Capture time strongly depends on the density and aspect ratio of the chosen disk model, the relative velocity of the stellar object with respect to the disk, and the AGN lifetime. We expect that for an AGN disk density $\rho \gtrsim 10^{-11}\rm g/cm^3$ and disk lifetime $\geq 1$Myr, there is a significant population of embedded stellar objects, which can fuel mergers detectable in gravitational waves with LIGO-Virgo and LISA.

The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

The mass, rate, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and Advanced Virgo provide challenges to traditional BBH formation and merger scenarios. An active galactic nucleus (AGN) disk provides a promising additional merger channel, because of the powerful influence of the gas that drives orbital evolution, makes encounters dissipative, and leads to migration. Previous work showed that stellar mass black holes (sBHs) in an AGN disk migrate to regions of the disk, known as migration traps, where positive and negative gas torques cancel out, leading to frequent BBH formation. Here we build on that work by simulating the evolution of additional sBHs that enter the inner disk by either migration or inclination reduction. We also examine whether the BBHs formed in our models have retrograde or prograde orbits around their centers of mass with respect to the disk, determining the orientation, relative to the disk, of the spin of the merged BBHs. We find that orbiters entering the inner disk form BBHs, that are commonly asymmetric in mass, with sBHs on resonant orbits near the migration trap. When these black holes reach 80 solar masses, they form BBHs with sBHs in the migration trap, which over 10 Myr reach about 1000 solar masses. We also find that 68% of the BBHs in our simulation orbit in the retrograde direction, which implies our merger channel could easily produce the low values of the dimensionless aligned spin, chi_eff, observed in a majority of BBH mergers detected thus far.

Accretion disks around supermassive black holes are promising sites for stellar mass black hole mergers detectable with LIGO. Here we present the results of Monte-Carlo simulations of black hole mergers within 1-d AGN disk models. For the spin distribution in the disk bulk, key findings are: (1) The distribution of $\chi_{\rm eff}$ is naturally centered around $\tilde{\chi}_{\rm eff} \approx 0.0$, (2) the width of the $\chi_{\rm eff}$ distribution is narrow for low natal spins. For the mass distribution in the disk bulk, key findings are: (3) mass ratios $\tilde{q} \sim 0.5-0.7$, (4) the maximum merger mass in the bulk is $\sim 100-200M_{\odot}$, (5) $\sim 1\%$ of bulk mergers involve BH $>50M_{\odot}$ with (6) $\simeq 80\%$ of bulk mergers are pairs of 1st generation BH. Additionally, mergers at a migration trap grow an IMBH with typical merger mass ratios $\tilde{q}\sim 0.1$. Ongoing LIGO non-detections of black holes $>10^{2}M_{\odot}$ puts strong limits on the presence of migration traps in AGN disks (and therefore AGN disk density and structure) as well as median AGN disk lifetime. The highest merger rate occurs for this channel if AGN disks are relatively short-lived ($\leq 1$Myr) so multiple AGN episodes can happen per Galactic nucleus in a Hubble time.

Advanced LIGO \& Advanced Virgo are detecting a large number of binary stellar origin black hole (BH) mergers. A promising channel for accelerated BH merger lies in active galactic nucleus (AGN) disks of gas around super-masssive black holes. Here we investigate the relative number of compact object mergers in AGN disk models, including BH, neutron stars (NS) and white dwarfs, via Monte Carlo simulations. We find the number of all merger types in the bulk disk grows $\propto t^{1/3}$ which is driven by the Hill sphere of the more massive merger component. Median mass ratios of NS-BH mergers in AGN disks are $\tilde{q}=0.07\pm 0.06(0.14\pm 0.07)$ for mass functions (MF) $M^{-1(-2)}$. If a fraction $f_{AGN}$ of the observed rate of BH-BH mergers (${\cal}{R}_{BH-BH}$) come from AGN, we expect a rate of NS-BH (NS-NS) mergers from the AGN channel of ${\cal}{R}_{BH-NS} \sim f_{AGN}[10,300]\rm{Gpc}^{-3} \rm{yr}^{-1}$,(${\cal}{R}_{NS-NS} \leq f_{AGN}400\rm{Gpc}^{-3} \rm{yr}^{-1}$). Allowing for the ratio of NS-NS/BH-BH LIGO search volumes, from preliminary O3 results we can rule out the AGN channel as a dominant contribution to observed NS-NS mergers. The number of lower mass gap events in this channel is a strong function of the nuclear MF and mass segregation efficiency. Compact object merger ratios derived from LIGO can restrict models of MF, mass segregation and populations embedded in AGN disks. The expected number of EM counterparts to NS-BH mergers in AGN disks at $z<1$ is $\sim [30,900]{\rm{yr}}^{-1}(f_{AGN}/0.1)$. EM searches for flaring due to embedded events in large AGN surveys will complement LIGO constraints on AGN models and the embedded populations that must live in them.

We present analysis of Chandra X-ray observations of seven quasars that were identified as candidate sub-parsec binary supermassive black hole (SMBH) systems in the Catalina Real-Time Transient Survey (CRTS) based on apparent periodicity in their optical light curves. Simulations predict close-separation accreting SMBH binaries will have different X-ray spectra than single accreting SMBHs, including harder or softer X-ray spectra, ripple-like profiles in the Fe K-$\alpha$ line, and distinct peaks in the spectrum due to the separation of the accretion disk into a circumbinary disk and mini-disks around each SMBH. We obtained Chandra observations to test these models and assess whether these quasars could contain binary SMBHs. We instead find that the quasar spectra are all well fit by simple absorbed power law models, with the rest frame 2-10 keV photon indices, $\Gamma$, and the X-ray-to-optical power slopes, $\alpha_{\rm OX}$, indistinguishable from the larger quasar population. This may indicate that these seven quasars are not truly sub-parsec binary SMBH systems, or it may simply reflect that our sample size was too small to robustly detect any differences. Alternatively, the X-ray spectral changes might only be evident at higher energies than probed by Chandra. Given the available models and current data, no firm conclusions are drawn. These observations will help motivate and direct further work on theoretical models of binary SMBH systems, such as modeling systems with thinner accretion disks and larger binary separations.

We present the results of a multi-wavelength follow up campaign for the luminous nuclear transient Gaia16aax, which was first identified in January 2016. The transient is spatially consistent with the nucleus of an active galaxy at z=0.25, hosting a black hole of mass $\rm \sim6\times10^8M_\odot$. The nucleus brightened by more than 1 magnitude in the Gaia G-band over a timescale of less than one year, before fading back to its pre-outburst state over the following three years. The optical spectra of the source show broad Balmer lines similar to the ones present in a pre-outburst spectrum. During the outburst, the $\rm H\alpha$ and $\rm H\beta$ emission lines develop a secondary peak. We also report on the discovery of two transients with similar light curve evolution and spectra: Gaia16aka and Gaia16ajq. We consider possible scenarios to explain the observed outbursts. We exclude that the transient event could be caused by a microlensing event, variable dust absorption or a tidal encounter between a neutron star and a stellar mass black hole in the accretion disk. We consider variability in the accretion flow in the inner part of the disk, or a tidal disruption event of a star $\geq 1 M_{\odot}$ by a rapidly spinning supermassive black hole as the most plausible scenarios. We note that the similarity between the light curves of the three Gaia transients may be a function of the Gaia alerts selection criteria.

We report on three redshift $z>2$ quasars with dramatic changes in their C IV emission lines, the first sample of changing-look quasars (CLQs) at high redshift. This is also the first time the changing-look behaviour has been seen in a high-ionisation emission line. SDSS J1205+3422, J1638+2827, and J2228+2201 show interesting behaviour in their observed optical light curves, and subsequent spectroscopy shows significant changes in the C IV broad emission line, with both line collapse and emergence being displayed on rest-frame timescales of $\sim$240-1640 days. These are rapid changes, especially when considering virial black hole mass estimates of $M_{\rm BH} > 10^{9} M_{\odot}$ for all three quasars. Continuum and emission line measurements from the three quasars show changes in the continuum-equivalent width plane with the CLQs seen to be on the edge of the full population distribution, and showing indications of an intrinsic Baldwin effect. We put these observations in context with recent state-change models, and note that even in their observed low-state, the C IV CLQs are generally above $\sim$5\% in Eddington luminosity.

The merger rate of stellar-mass black hole binaries (sBHBs) inferred by the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO) suggests the need for an efficient source of sBHB formation. Active galactic nucleus (AGN) disks are a promising location for the formation of these sBHBs, as well as binaries of other compact objects, because of powerful torques exerted by the gas disk. These gas torques cause orbiting compact objects to migrate towards regions in the disk where inward and outward torques cancel, known as migration traps. We simulate the migration of stellar mass black holes in an example of a model AGN disk, using an augmented N-body code that includes analytic approximations to migration torques, stochastic gravitational forces exerted by turbulent density fluctuations in the disk, and inclination and eccentricity dampening produced by passages through the gas disk, in addition to the standard gravitational forces between objects. We find that sBHBs form rapidly in our model disk as stellar-mass black holes migrate towards the migration trap. These sBHBs are likely to subsequently merge on short time-scales. The process continues, leading to the build-up of a population of over-massive stellar-mass black holes. The formation of sBHBs in AGN disks could contribute significantly to the sBHB merger rate inferred by LIGO.

Accretion disks around supermassive black holes (SMBH) are promising sites for stellar mass black hole (BH) mergers due to mass segregation and merger acceleration by disk gas torques. Here we show that a GW-kick at BH merger causes ram-pressure stripping of gas within the BH Hill sphere. If $R_{H}\geq H$, the disk height, an off-center UV flare at $a_{\rm BH} \sim 10^{3}r_{g}$ emerges within $t_{\rm UV} \sim \rm{O}(2 {\rm days})(a_{\rm BH}/10^{3}r_{g})(M_{\rm SMBH}/10^{8}M_{\odot})(v_{\rm kick}/10^{2}\rm{km/s})$ post-merger and lasts O$(R_{H}/v_{\rm kick}) \sim \rm{O}(5 t_{\rm UV}$). The flare emerges with luminosity O($10^{42}{\rm erg/s})(t_{\rm UV}/2{\rm days})^{-1}(M_{\rm Hill}/1M_{\odot})(v_{\rm kick}/10^{2}{\rm km/s})^{2}$. AGN optical/UV photometry alters and asymmetric broad emission line profiles can develop after weeks. If $R_{H}<H$, detectability depends on disk optical depth. Follow-up by large optical sky surveys is optimized for small GW error volumes and for LIGO/Virgo triggers $>50M_{\odot}$.

Low ionization emission line regions (LINERs) are a heterogeneous collection of up to $1/3$ of galactic nuclei in the local Universe. It is unclear whether LINERs are simply the result of low accretion rates onto supermassive black holes or whether they include a large number of optically thick radiatively inefficient but super-Eddington accretion flows (RIAFs). Optically thick RIAFs are typically disks of large scale-height or quasi-spherical gas flows. These should be dense enough to trap and merge a large number of the stellar mass black holes, which we expect to exist in galactic nuclei. Electromagnetic observations of photospheres of accretion flows do not allow us to break model degeneracies. However, gravitational wave observations probe the interior of accretion flows where the merger of stellar mass black holes can be greatly accelerated over the field rate. Here we show that the upper limits on the rate of black hole mergers observed with LIGO demonstrate that most LINERs cannot be optically thick RIAFs.

In this paper, we continue our study on the evolution of black holes (BHs) that receive velocity kicks at the origin of their host star cluster potential. We now focus on BHs in rotating clusters that receive a range of kick velocities in different directions with respect to the rotation axis. We perform N-body simulations to calculate the trajectories of the kicked BHs and develop an analytic framework to study their motion as a function of the host cluster and the kick itself. Our simulations indicate that for a BH that is kicked outside of the cluster's core, as its orbit decays in a rotating cluster the BH will quickly gain angular momentum as it interacts with stars with high rotational frequencies. Once the BH decays to the point where its orbital frequency equals that of local stars, its orbit will be circular and dynamical friction becomes ineffective since local stars will have low relative velocities. After circularization, the BH's orbit decays on a longer timescale than if the host cluster was not rotating. Hence BHs in rotating clusters will have longer orbital decay times. The timescale for orbit circularization depends strongly on the cluster's rotation rate and the initial kick velocity, with kicked BHs in slowly rotating clusters being able to decay into the core before circularization occurs. The implication of the circularization phase is that the probability of a BH undergoing a tidal capture event increases, possibly aiding in the formation of binaries and high-mass BHs.

We present the results of a systematic search for quasars in the Catalina Real-time Transient Survey exhibiting both strong photometric and spectroscopic variability over a decadal baseline. We identify 73 sources with specific patterns of optical and mid-IR photometric behavior and a defined spectroscopic change. These "Changing-State" quasars (CSQs) form a higher luminosity sample to complement existing sets of "Changing-Look" AGN and quasars in the literature. The CSQs (by selection) exhibit larger photometric variability than the CLQs. The spectroscopic variability is marginally stronger in the CSQs than CLQs as defined by the change in H$\beta$/[OIII] ratio. We find 36 sources with declining H$\beta$ flux, 37 sources with increasing H$\beta$ flux and discover seven sources with $z > 0.8$, further extending the redshift arm. Our CSQ sample compares to the literature CLQ objects in similar distributions of H$\beta$ flux ratios and differential Eddington ratios between high (bright) and low (dim) states. Taken as a whole, we find that this population of extreme varying quasars is associated with changes in the Eddington ratio and the timescales imply cooling/heating fronts propagating through the disk.

The conventional view of stellar-mass binary black hole (sBBH) mergers is that there should not be enough matter present to produce a detectable electromagnetic transient. However, there ARE a number of mechanisms for producing such a counterpart through accretion of matter from various reserves, charged black holes, or interactions with magnetic or exotic fields. After reviewing these mechanisms, we describe what we can learn from multi-messenger observations of sBBH mergers in the areas of stellar evolution and compact binary formation, cosmological measurements, searches for charged black holes, and tests of general relativity and fundamental physics. We end with some remarks about the need to support both observing capabilities and modeling in the next decade and beyond.

The stellar mass binary black hole (sBBH) mergers presently detected by LIGO may originate wholly or in part from binary black hole mergers embedded in disks of gas around supermassive black holes. Determining the contribution of these active galactic nucleus (AGN) disks to the sBBH merger rate enables us to uniquely measure important parameters of AGN disks, including their typical density, aspect ratio, and lifetime, thereby putting unique limits on an important element of galaxy formation. For the first time, gravitational waves will allow us to reveal the properties of the hidden interior of AGN disks, while electromagnetic radiation (EM) probes the disk photosphere. The localization of sBBH merger events from LIGO is generally insufficient for association with a single EM counterpart. However, the contribution to the LIGO event rate from rare source types (such as AGNs) can be determined on a statistical basis. To determine the contribution to the sBBH rate from AGNs in the next decade requires: {\em 1) a complete galaxy catalog for the LIGO search volume, 2) strategic multi-wavelength EM follow-up of LIGO events and 3) significant advances in theoretical understanding of AGN disks and the behavior of objects embedded within them.}

Electromagnetic observations of the sky have been the basis for our study of the Universe for millennia, cosmic ray studies are now entering their second century, the first neutrinos from an astrophysical source were identified three decades ago, and gravitational waves were directly detected only four years ago. Detections of these messengers are now common. Astrophysics will undergo a revolution in the 2020s as multimessenger detections become routine. The 8th Astro2020 Thematic Area is Multimessenger Astronomy and Astrophysics, which includes the identification of the sources of gravitational waves, astrophysical and cosmogenic neutrinos, cosmic rays, and gamma-rays, and the coordinated multimessenger and multiwavelength follow-ups. Identifying and characterizing multimessenger sources enables science throughout and beyond astrophysics. Success in the multimessenger era requires: (i) sensitive coverage of the non-electromagnetic messengers, (ii) full coverage of the electromagnetic spectrum, with either fast-response observations or broad and deep high-cadence surveys, and (iii) improved collaboration, communication, and notification platforms.

Changing-look quasars are a recently identified class of active galaxies in which the strong UV continuum and/or broad optical hydrogen emission lines associated with unobscured quasars either appear or disappear on timescales of months to years. The physical processes responsible for this behaviour are still debated, but changes in the black hole accretion rate or accretion disk structure appear more likely than changes in obscuration. Here we report on four epochs of spectroscopy of SDSS J110057.70-005304.5, a quasar at a redshift of $z=0.378$ whose UV continuum and broad hydrogen emission lines have faded, and then returned over the past $\approx$20 years. The change in this quasar was initially identified in the infrared, and an archival spectrum from 2010 shows an intermediate phase of the transition during which the flux below rest-frame $\approx$3400\AA\ has decreased by close to an order of magnitude. This combination is unique compared to previously published examples of changing-look quasars, and is best explained by dramatic changes in the innermost regions of the accretion disk. The optical continuum has been rising since mid-2016, leading to a prediction of a rise in hydrogen emission line flux in the next year. Increases in the infrared flux are beginning to follow, delayed by a $\sim$3 year observed timescale. If our model is confirmed, the physics of changing-look quasars are governed by processes at the innermost stable circular orbit (ISCO) around the black hole, and the structure of the innermost disk. The easily identifiable and monitored changing-look quasars would then provide a new probe and laboratory of the nuclear central engine.

We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 \mu$m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of fifteen decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 which, when compared with the 2006 archival SDSS spectrum, reveals that the broad H$\beta$ emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes in the quasar luminosity over this timescale, and favor changes in the innermost regions of the accretion disk that occur on the thermal and heating/cooling front timescales. We discuss possible physical triggers that could cause these changes, and predict the multiwavelength signatures that could distinguish these physical scenarios.

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.

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 consider how dynamical friction acts on black holes that receive a velocity kick while located at the center of a gravitational potential, analogous to a star cluster, due to either a natal kick or the anisotropic emission of gravitational waves during a black hole-black hole merger. Our investigation specifically focuses on how well various Chandrasekhar-based dynamical friction models can predict the orbital decay of kicked black holes with $m_{bh} \lesssim 100 M_\odot$ due to an inhomogeneous background stellar field. In general, the orbital evolution of a kicked black hole follows that of a damped oscillator where two-body encounters and dynamical friction serve as sources of damping. However, we find models for approximating the effects of dynamical friction do not accurately predict the amount of energy lost by the black hole if the initial kick velocity $v_{k}$ is greater than the stellar velocity dispersion $\sigma$. For all kick velocities, we also find that two-body encounters with nearby stars can cause the energy evolution of a kicked BH to stray significantly from standard dynamical friction theory as encounters can sometimes lead to an energy gain. For larger kick velocities, we find the orbital decay of a black hole departs from classical theory completely as the black hole's orbital amplitude decays linearly with time as opposed to exponentially. Therefore, we have developed a linear decay formalism which scales linearly with black hole mass and $\frac{v_{k}}{\sigma}$ in order to account for the variations in the local gravitational potential.

White dwarfs (WDs) are believed to detonate via explosive Carbon-fusion in a Type Ia Supernova when their temperature and/or density reach the point where Carbon is ignited in a runaway reaction. Observations of the Type Ia supernova (SN) rate imply all WD binaries that merge through the emission of gravitational radiation within a Hubble time should result in SNe, regardless of total mass. Here we investigate the conditions under which a single WD in a binary system might extract energy from its orbit, depositing enough energy into a resonant mode such that it detonates before merger. We show that, ignoring non-linear effects, in a WD binary in tidal lock at small binary separations, the sustained tidal forcing of a low-order quadrupolar g-mode or a harmonic of a low-order quadrupolar p-mode could in principle drive the average temperature of Carbon nuclei in the mode over the runaway fusion threshold. If growing mode energy is thermalized at a core/atmosphere boundary, rapid Helium burning and inward-travelling p-waves may result in core detonation. Thermalization at a boundary in the core can also result in detonation. If energy can be efficiently transferred from the orbit to modes as the WD binary passes through resonances, the WD merger timescale will be shortened by Myr-Gyr compared to expected timescales from GW-emission alone and GW detectors will observe deviations from predicted chirp profiles in resolved WD binaries. Future work in this area should focus on whether tidal locking in WD binaries is naturally driven towards low-order mode frequencies.

Accretion disks around supermassive black holes (SMBHs) in active galactic nuclei contain stars, stellar mass black holes, and other stellar remnants, which perturb the disk gas gravitationally. The resulting density perturbations in turn exert torques on the embedded masses causing them to migrate through the disk in a manner analogous to the behavior of planets in protoplanetary disks. We determine the strength and direction of these torques using an empirical analytic description dependent on local disk gradients, applied to two different analytic, steady-state disk models of SMBH accretion disks. We find that there are radii in such disks where the gas torque changes sign, trapping migrating objects. Our analysis shows that major migration traps generally occur where the disk surface density gradient changes sign from positive to negative, around 20--300$R_{\rm g}$, where $R_{\rm g}=2GM/c^{2}$ is the Schwarzschild radius. At these traps, massive objects in the AGN disk can accumulate, collide, scatter, and accrete. Intermediate mass black hole formation is likely in these disk locations, which may lead to preferential gap and cavity creation at these radii. Our model thus has significant implications for SMBH growth as well as gravitational wave source populations.

The barycenter of a massive black hole binary will lie outside the event horizon of the primary black hole for modest values of mass ratio and binary separation. Analagous to radial velocity shifts in stellar emission lines caused by the tug of planets, the radial velocity of the primary black hole around the barycenter can leave a tell-tale oscillation in the broad component of Fe K$\alpha$ emission from accreting gas. Near-future X-ray telescopes such as Astro-H and Athena will have the energy resolution ($\delta E/E \lesssim 10^{-3}$) to search nearby active galactic nuclei (AGN) for the presence of binaries with mass ratios $q \gtrsim 0.01$, separated by several hundred gravitational radii. The general-relativistic and Lense-Thirring precession of the periapse of the secondary orbit imprints a detectable modulation on the oscillations. The lowest mass binaries in AGN will oscillate many times within typical X-ray exposures, leading to a broadening of the line wings and an over-estimate of black hole spin in these sources. Detection of periodic oscillations in the AGN line centroid energy will reveal a massive black hole binary close to merger and will provide an early warning of gravitational radiation emission.

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.

JWST's Near-Infrared Imager and Slitless Spectrograph (NIRISS) includes an Aperture Masking Interferometry (AMI) mode designed to be used between 2.7{\mu}m and 4.8{\mu}m. At these wavelengths, it will have the highest angular resolution of any mode on JWST, and, for faint targets, of any existing or planned infrastructure. NIRISS AMI is uniquely suited to detect thermal emission of young massive planets and will permit the characterization of the mid-IR flux of exoplanets discovered by the GPI and SPHERE adaptive optics surveys. It will also directly detect massive planets found by GAIA through astrometric accelerations, providing the first opportunity ever to get both a mass and a flux measurement for non-transiting giant planets. NIRISS AMI will also enable the study of the nuclear environment of AGNs.

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.

Due to feedback from accretion onto supermassive black holes (SMBHs), Active Galactic Nuclei (AGNs) are believed to play a key role in LambdaCDM cosmology and galaxy formation. However, AGNs' extreme luminosities and the small angular size of their accretion flows create a challenging imaging problem. We show James Webb Space Telescope's Near Infrared Imager and Slitless Spectrograph (JWST-NIRISS) Aperture Masking Interferometry (AMI) mode will enable true imaging (i.e. without any requirement of prior assumptions on source geometry) at ~65 mas angular resolution at the centers of AGNs. This is advantageous for studying complex extended accretion flows around SMBHs, and in other areas of angular-resolution-limited astrophysics. By simulating data sequences incorporating expected sources of noise, we demonstrate that JWST-NIRISS AMI mode can map extended structure at a pixel-to-pixel contrast of ~10^{-2} around an L=7.5 point source, using short exposure times (minutes). Such images will test models of AGN feedback, fuelling and structure (complementary with ALMA observations), and are not currently supported by any ground-based IR interferometer or telescope. Binary point source contrast with NIRISS is ~10^{-4} (for observing binary nuclei in merging galaxies), significantly better than current ground-based optical or IR interferometry. JWST-NIRISS' seven-hole non-redundant mask has a throughput of 15%, and utilizes NIRISS' F277W (2.77\micron), F380M (3.8\micron), F430M (4.3\micron), and F480M (4.8\micron) filters. NIRISS' square pixels are 65mas per side, with a field of view ~2\arcmin x 2\arcmin. We also extrapolate our results to AGN science enabled by non-redundant masking on future 2.4m and 16m space telescopes working at long-UV to near-IR wavelengths.

(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.

Here we propose a mechanism for efficiently growing intermediate mass black holes (IMBH) in disks around supermassive black holes. Stellar mass objects can efficiently agglomerate when facilitated by the gas disk. Stars, compact objects and binaries can migrate, accrete and merge within disks around supermassive black holes. While dynamical heating by cusp stars excites the velocity dispersion of nuclear cluster objects (NCOs) in the disk, gas in the disk damps NCO orbits. If gas damping dominates, NCOs remain in the disk with circularized orbits and large collision cross-sections. IMBH seeds can grow extremely rapidly by collisions with disk NCOs at low relative velocities, allowing for super-Eddington growth rates. Once an IMBH seed has cleared out its feeding zone of disk NCOs, growth of IMBH seeds can become dominated by gas accretion from the AGN disk. However, the IMBH can migrate in the disk and expand its feeding zone, permitting a super-Eddington accretion rate to continue. Growth of IMBH seeds via NCO collisions is enhanced by a pile-up of migrators. We highlight the remarkable parallel between the growth of IMBH in AGN disks with models of giant planet growth in protoplanetary disks. If an IMBH becomes massive enough it can open a gap in the AGN disk. IMBH migration in AGN disks may stall, allowing them to survive the end of the AGN phase and remain in galactic nuclei. Our proposed mechanisms should be more efficient at growing IMBH in AGN disks than the standard model of IMBH growth in stellar clusters. Dynamical heating of disk NCOs by cusp stars is transferred to the gas in a AGN disk helping to maintain the outer disk against gravitational instability. Model predictions, observational constraints and implications are discussed in a companion paper (Paper II).

Galactic nuclei should contain a cluster of stars and compact objects in the vicinity of the central supermassive black hole due to stellar evolution, minor mergers and gravitational dynamical friction. By analogy with protoplanetary migration, nuclear cluster objects (NCOs) can migrate in the accretion disks that power active galactic nuclei by exchanging angular momentum with disk gas. Here we show that an individual NCO undergoing runaway outward migration comparable to Type III protoplanetary migration can generate an accretion rate corresponding to Seyfert AGN or quasar luminosities. Multiple migrating NCOs in an AGN disk can dominate traditional viscous disk accretion and at large disk radii, ensemble NCO migration and accretion could provide sufficient heating to prevent the gravitational instability from consuming disk gas in star formation. The magnitude and energy of the X-ray soft excess observed at ~0.1-1keV in Seyfert AGN could be explained by a small population of ~10^{2}-10^{3} accreting stellar mass black holes or a few ULXs. NCO migration and accretion in AGN disks are therefore extremely important mechanisms to add to realistic models of AGN disks.

Small and intermediate mass black holes should be expected in galactic nuclei as a result of stellar evolution, minor mergers and gravitational dynamical friction. If these minor black holes accrete as X-ray binaries or ultra-luminous X-ray sources, and are associated with star formation, they could account for observations of many low luminosity AGN or LINERs. Accreting and inspiralling intermediate mass black holes could provide a crucial electromagnetic counterpart to strong gravitational wave signatures, allowing tests of strong gravity. Here we discuss observational signatures of minor black holes in galactic nuclei and we demonstrate that optical line ratios observed in LINERs or transition-type objects can be produced by an ionizing radiation field from ULXs. We conclude by discussing constraints from existing observations as well as candidates for future study.

We propose a new mechanism for the delivery of gas to the heart of galactic nuclei. We show that warm halo clouds must periodically impact galactic centers and potentially deliver a large (~10^{4-6} M_{solar}) mass of gas to the galactic nucleus in a singular event. The impact of an accreting warm halo cloud originating far in the galactic halo can, depending on mixing, produce a nuclear starburst of low metallicity stars as well as low luminosity accretion onto the central black hole. Based on multiphase cooling around a LambdaCDM distribution of halos we calculate the nuclear impact rate, the mass captured by the central black hole and the fraction of active nuclei for impacting cloud masses in the range ~10^{4}-10^{6}M_{solar}. If there is moderate braking during cloud infall, our model predicts an average fraction of low luminosity active nuclei consistent with observations.

We investigate the role of host galaxy classification and black hole mass in a heterogeneous sample of 276 mostly nearby (z<0.1) X-ray and IR selected AGN. Around 90% of Seyfert 1 AGN in bulge-dominated host galaxies (without disk contamination) span a very narrow range in the observed 12um to 2-10keV luminosity ratio (1<R_{IR/X}<7). This narrow dispersion incorporates all possible variations among AGN central engines, including accretion mechanism and efficiency, disk opening angle, orientation to sightline, covering fraction of absorbing material, patchiness of X-ray corona and measured variability. As a result, all models of X-ray and IR production in AGN are very strongly constrained. Among Seyfert 1 AGN, median X-ray and IR luminosities increase with black hole mass at >99% confidence. Using ring morphology of the host galaxy as a proxy for lack of tidal interaction, we find that AGN luminosity in host galaxies within 70Mpc is independent of host galaxy interaction for $\sim$ Gyrs, suggesting that the timescale of AGN activity due to secular evolution is much shorter than that due to tidal interactions. We find that LINER hosts have lower 12um luminosity than the median 12um luminosity of normal disk- and bulge-dominated galaxies which may represent observational evidence for past epochs of feedback that supressed star formation in LINER host galaxies. We propose that nuclear ULXs may account for the X-ray emission from LINER 2s without flat-spectrum, compact radio cores. We confirmed the robustness of our results in X-rays by comparing them with the 14-195keV 22-month BAT survey of AGN, which is all-sky and unbiased by photoelectric absorption.

Active galactic nuclei (AGN) produce vast amounts of high energy radiation deep in their central engines. X-rays either escape the AGN or are absorbed and re-emitted mostly as IR. By studying the dispersion in the ratio of observed mid-IR luminosity to observed 2-10keV X-ray luminosity (R_{ir/x}) in AGN we can investigate the reprocessing material (possibly a torus or donut of dust) in the AGN central engine, independent of model assumptions. We studied the ratio of observed mid-IR and 2-10keV X-ray luminosities in a heterogeneous sample of 245 AGN from the literature. We found that when we removed AGN with prominent jets, ~90% of Type I AGN lay within a very tight dispersion in luminosity ratio (1<R_{ir/x}<30). This implies that the AGN central engine is extremely uniform and models of the physical AGN environment (e.g. cloud cover, turbulent disk, opening angle of absorbing structures such as dusty tori) must span a very narrow range of parameters. We also found that the far-IR(100um) to mid-IR (12um) observed luminosity ratio is an effective descriminator between heavily obscured AGN and relatively unobscured AGN.

We report the detection of H2CO (formaldehyde) around the carbon-rich AGB star, IRC+10216. We find a fractional abundance with respect to molecular hydrogen of x(H2CO)= (1.3 {+1.5}{-0.8}) x 10^{-8}. This corresponds to a formaldehyde abundance with respect to water vapor of x(H2CO)/x(H2O)=(1.1 +/- 0.2) x 10^{-2}, in line with the formaldehyde abundances found in Solar System comets, and indicates that the putative extrasolar cometary system around IRC+10216 may have a similar chemical composition to Solar System comets. However, we also failed to detect CH3OH (methanol) around IRC+10216 and our upper limit of x(CH3OH)/x(H2O) < 7.7 x 10^{-4}, (3 sigma), indicates that methanol is substantially underabundant in IRC+10216, compared to Solar System comets. We also conclude, based on offset observations, that formaldehyde has an extended source in the envelope of IRC+10216 and may be produced by the photodissociation of a parent molecule, similar to the production mechanism for formaldehyde in Solar System comet comae. Preliminary mapping observations also indicate a possible asymmetry in the spatial distribution of formaldehyde around IRC+10216, but higher signal-to-noise observations are required to confirm this finding. This study is based on observations carried out with the IRAM 30m telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain). (abridged)

We report the detection of the 1665 and 1667 MHz main lines of OH (hydroxyl) and upper limits on the 1612 MHz satellite line of OH towards the carbon-rich AGB star, IRC+10216. We find a beam averaged fractional abundance x(OH)~4x10^{-8}. This detection supports the identification by Melnick et al. (2001) of the 1_10-1_01 transition of water vapor with a 556.936 GHz emission feature detected towards IRC+10216, since OH is the expected photodissociation product of water vapor. The shape of the OH lines, however, differs significantly from the shape expected based on the observations of Melnick et al. (2001). Possible explanations for the anomalous shapes of the 1665 and 1667 MHz lines are discussed. The most likely explanations for the unexpected OH line shapes are either masing or an asymmetric distribution of OH molecules around IRC+10216.

We argue that the presence of water vapor in the circumstellar outflow of a carbon-rich AGB star is potentially a distinctive signature of extra-solar cometary systems. Detailed models show that at suitable distances from the star, water ice can survive well into the carbon-rich AGB phase; water vapor abundances as large as 10^-6 could result from the vaporization of a collection of orbiting icy bodies with a total mass comparable to what might have been originally present in the solar system's Kuiper Belt. In particular, the recently-reported detection by the Submillimeter Wave Astronomy Satellite of water vapor in the circumstellar outflow of IRC+10216 can be explained if ~10 Earth masses of ice is present at a distance ~300 AU from that carbon-rich star. Future observations with the Herschel Space Observatory (HSO, formerly known as FIRST) will facilitate sensitive multi-transition observations of water, yielding line ratios that can establish the radial distribution of water vapor in IRC+10216. The greater sensitivity of HSO will also allow searches for water vapor to be carried out in a much larger sample of carbon-rich AGB stars.

Since 1995, astronomers have discovered planets with masses comparable to that of Jupiter (318 times Earth's mass) in orbit around approximately 60 stars. Although unseen directly, the presence of these planets is inferred by the small reflex motions that they gravitationally induce on the star they orbit; these result in small periodic wavelength shifts in the stellar spectrum. Since this method favors the detection of massive objects orbiting in close proximity to the star, the question of whether these systems also contain analogs of the smaller constituents of our Solar System has remained unanswered. Using an alternative approach, we report here observations of an aging carbon-star, IRC+10216, that reveal the presence of circumstellar water vapor, a molecule not expected in measurable abundances around such a star and thus a distinctive signature of an orbiting cometary system. The only plausible explanation for this water vapor is that the recent evolution of IRC+10216 - which is accompanied by a prodigious increase in its luminosity - is now causing the vaporization of a collection of orbiting icy bodies, a process first considered in a previous theoretical study.