We study the association between active galactic nuclei (AGN) and environment at scales of $0.01-1\ h^{-1}$Mpc in the IllustrisTNG (TNG100) simulated universe. We identify supermassive black hole (BH) pairs and multiples within scales of 0.01, 0.1, & 1 $h^{-1}$Mpc and examine their AGN activity in relation to randomly-selected pairs and multiples. The number density of BHs in TNG100 is $n=0.06\,h^3$Mpc$^{-3}$ at $z\lesssim1.5$ ($n=0.02\,h^3$ Mpc$^{-3}$ at $z=3$). About $\sim10$% and $\sim1$% of them live in pairs and multiples, respectively, within 0.1 $h^{-1}$Mpc scales. We find that BH systems have enhanced likelihood (up to factors of 3-6) of containing high Eddington ratio ($\eta\gtrsim0.7$) AGN compared to random pairs and multiples. Conversely, the likelihood of an AGN to live in 0.1$h^{-1}$Mpc scale systems is also higher (by factors $\sim4$ for $\eta\gtrsim0.7$) compared to random pairs and multiples. We also estimate that $\sim10$% of ultra-hard X-ray selected AGN in TNG100 have detectable 2-10 keV AGN companions on $0.1\ h^{-1}$Mpc scales, in agreement with observations. On larger spatial scales ($\sim 1$ $h^{-1}$Mpc), however, no significant enhancement is associated with BH pairs and multiples, even at high Eddington ratios. The enhancement of AGN activity in rich, small-scale ($\lesssim0.1$ $h^{-1}$Mpc) environments is therefore likely to be driven by galaxy interactions and mergers. Nonetheless, the overall percentage of AGN that live in $\lesssim0.1$ $h^{-1}$Mpc scale multiples is still subdominant (at most $\sim40$% for the highest Eddington ratio AGN). Furthermore, the enhancement in Eddington ratios of BH systems(as well as merging BHs) is only up to factors of $\sim2-3$. Thus, our results support the existence of a merger-AGN connection, but they also suggest that mergers and interactions play a relatively minor role in fueling the AGN population as a whole.

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

Galaxy mergers are predicted to trigger accretion onto the central supermassive black holes, with the highest rates occurring during final coalescence. Previously, we have shown elevated rates of both optical and mid-IR selected active galactic nuclei (AGN) in post-mergers, but to date the prevalence of X-ray AGN has not been examined in the same systematic way. We present XMM-Newton data of 43 post-merger galaxies selected from the Sloan Digital Sky Survey along with 430 non-interacting control galaxies matched in stellar mass, redshift, and environment in order to test for an excess of hard X-ray (2-10 keV) emission in post-mergers attributable to triggered AGN. We find 2 X-ray detections in the post-mergers (4.7^{+9.3}_{-3.8}%) and 9 in the controls (2.1^{+1.5}_{-1.0}%), an excess of 2.22^{+4.44}_{-2.22}, where the confidence intervals are 90%. While we therefore do not find statistically significant evidence for an X-ray AGN excess in post-mergers (p = 0.26), we find a factor of ~17 excess of mid-IR AGN in our sample, consistent with past work and inconsistent with the observed X-ray excess (p = 2.7 x 10^{-4}). Dominant, luminous AGN are therefore more frequent in post-mergers, and the lack of a comparable excess of 2-10 keV X-ray AGN suggests that AGN in post-mergers are more likely to be heavily obscured. Our results are consistent with the post-merger stage being characterised by enhanced AGN fueling, heavy AGN obscuration, and more intrinsically luminous AGN, in line with theoretical predictions.

Theory predicts that a supermassive black hole binary (SMBHB) could be observed as a luminous active galactic nucleus (AGN) that periodically varies on the order of its orbital timescale. In X-rays, periodic variations could be caused by mechanisms including relativistic Doppler boosting and shocks. Here we present the first systematic search for periodic AGNs using $941$ hard X-ray light curves (14-195 keV) from the first 105 months of the Swift Burst Alert Telescope (BAT) survey (2004-2013). We do not find evidence for periodic AGNs in Swift-BAT, including the previously reported SMBHB candidate MCG+11$-$11$-$032. We find that the null detection is consistent with the combination of the upper-limit binary population in AGNs in our adopted model, their expected periodic variability amplitudes, and the BAT survey characteristics. We have also investigated the detectability of SMBHBs against normal AGN X-ray variability in the context of the eROSITA survey. Under our assumptions of a binary population and the periodic signals they produce which have long periods of hundreds of days, up to $13$% true periodic binaries can be robustly distinguished from normal variable AGNs with the ideal uniform sampling. However, we demonstrate that realistic eROSITA sampling is likely to be insensitive to long-period binaries because longer observing gaps reduce their detectability. In contrast, large observing gaps do not diminish the prospect of detecting binaries of short, few-day periods, as 19% can be successfully recovered, the vast majority of which can be identified by the first half of the survey.

Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing.

In contrast to massive galaxies with Solar or super-Solar gas phase metallicities, very few Active Galactic Nuclei (AGN) are found in low-metallicity dwarf galaxies. Such a population could provide insight into the origins of supermassive black holes. Here we report near-infrared spectroscopic and X-ray observations of SDSS J105621.45+313822.1, a low-mass, low-metallicity galaxy with optical narrow line ratios consistent with star forming galaxies but a broad H$\alpha$ line and mid-infrared colors consistent with an AGN. We detect the [Si VI] 1.96$\mu$m coronal line and a broad Pa$\alpha$ line with a FWHM of $850 \pm 25$~km~s$^{-1}$. Together with the optical broad lines and coronal lines seen in the SDSS spectrum, we confirm the presence of a highly accreting black hole with mass $(2.2 \pm 1.3) \times 10^{6}$~M$_{\odot}$, with a bolometric luminosity of $\approx10^{44}$~erg~s$^{-1}$ based on the coronal line luminosity, implying a highly accreting AGN. Chandra observations reveal a weak nuclear point source with $L_{\textrm{X,2-10 keV}} = (2.3 \pm 1.2) \times 10^{41}$~erg~s$^{-1}$, $\sim 2$ orders of magnitude lower than that predicted by the mid-infrared luminosity, suggesting that the AGN is highly obscured despite showing broad lines in the optical spectrum. The low X-ray luminosity and optical narrow line ratios of J1056+3138 highlight the limitations of commonly employed diagnostics in the hunt for AGNs in the low metallicity low mass regime.

ESA and NASA are moving forward with plans to launch LISA around 2034. With data from the Illustris cosmological simulation, we provide analysis of LISA detection rates accompanied by characterization of the merging massive black hole population. Massive black holes of total mass $\sim10^5-10^{10} M_\odot$ are the focus of this study. We evolve Illustris massive black hole mergers, which form at separations on the order of the simulation resolution ($\sim$kpc scales), through coalescence with two different treatments for the binary massive black hole evolutionary process. The coalescence times of the population, as well as physical properties of the black holes, form a statistical basis for each evolutionary treatment. From these bases, we Monte Carlo synthesize many realizations of the merging massive black hole population to build mock LISA detection catalogs. We analyze how our massive black hole binary evolutionary models affect detection rates and the associated parameter distributions measured by LISA. With our models, we find massive black hole binary detection rates with LISA of $\sim0.5-1$ yr$^{-1}$ for massive black holes with masses greater than $10^5M_\odot$. This should be treated as a lower limit primarily because our massive black hole sample does not include masses below $10^5M_\odot$, which may significantly add to the observed rate. We suggest reasons why we predict lower detection rates compared to much of the literature.

The co-evolution of galaxies and the supermassive black holes (SMBHs) at their centers via hierarchical galaxy mergers is a key prediction of $\Lambda$CDM cosmology. As gas and dust are funneled to the SMBHs during the merger, the SMBHs light up as active galactic nuclei (AGNs). In some cases, a merger of two galaxies can encounter a third galaxy, leading to a triple merger, which would manifest as a triple AGN if all three SMBHs are simultaneously accreting. Using high-spatial resolution X-ray, near-IR, and optical spectroscopic diagnostics, we report here a compelling case of an AGN triplet with mutual separations < 10 kpc in the advanced merger SDSS J084905.51+111447.2 at z = 0.077. The system exhibits three nuclear X-ray sources, optical spectroscopic line ratios consistent with AGN in each nucleus, a high excitation near-IR coronal line in one nucleus, and broad Pa$\alpha$ detections in two nuclei. Hard X-ray spectral fitting reveals a high column density along the line of sight, consistent with the picture of late-stage mergers hosting heavily absorbed AGNs. Our multiwavelength diagnostics support a triple AGN scenario, and we rule out alternative explanations such as star formation activity, shock-driven emission, and emission from fewer than three AGN. The dynamics of gravitationally bound triple SMBH systems can dramatically reduce binary SMBH inspiral timescales, providing a possible means to surmount the "Final Parsec Problem." AGN triplets in advanced mergers are the only observational forerunner to bound triple SMBH systems and thus offer a glimpse of the accretion activity and environments of the AGNs prior to the gravitationally-bound triple phase.

Observations and theoretical simulations suggest that a significant fraction of merger-triggered accretion onto supermassive black holes is highly obscured, particularly in late-stage galaxy mergers, when the black hole is expected to grow most rapidly. Starting with the Wide-Field Infrared Survey Explorer all-sky survey, we identified a population of galaxies whose morphologies suggest ongoing interaction and which exhibit red mid-infrared colors often associated with powerful active galactic nuclei (AGNs). In a follow-up to our pilot study, we now present Chandra/ACIS and XMM-Newton X-ray observations for the full sample of the brightest 15 IR-preselected mergers. All mergers reveal at least one nuclear X-ray source, with 8 out of 15 systems exhibiting dual nuclear X-ray sources, highly suggestive of single and dual AGNs. Combining these X-ray results with optical line ratios and with near-IR coronal emission line diagnostics, obtained with the near-IR spectrographs on the Large Binocular Telescope, we confirm that 13 out of the 15 mergers host AGNs, two of which host dual AGNs. Several of these AGNs are not detected in the optical. All X-ray sources appear X-ray weak relative to their mid-infrared continuum, and of the nine X-ray sources with sufficient counts for spectral analysis, eight reveal strong evidence of high absorption with column densities of $N_\mathrm{H} \gtrsim 10^{23}$~cm$^{-2}$. These observations demonstrate that a significant population of single and dual AGNs are missed by optical studies, due to high absorption, adding to the growing body of evidence that the epoch of peak black hole growth in mergers occurs in a highly obscured phase.

The discoveries made over the past 20 years by Chandra and XMM-Newton surveys in conjunction with multiwavelength imaging and spectroscopic data available in the same fields have significantly changed the view of the supermassive black hole (SMBH) and galaxy connection. These discoveries have opened up several exciting questions that are beyond the capabilities of current X-ray telescopes and will need to be addressed by observatories in the next two decades. As new observatories peer into the early Universe, we will begin to understand the physics and demographics of SMBH infancy (at $z>6$) and investigate the influence of their accretion on the formation of the first galaxies ($\S$ 2.1). We will also be able to understand the accretion and evolution over the cosmic history (at $z\sim$1-6) of the full population of black holes in galaxies, including low accretion rate, heavily obscured AGNs at luminosities beyond the reach of current X-ray surveys ($\S$2.2 and $\S$2.3), enabling us to resolve the connection between SMBH growth and their environment.

Gravitational wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. Multi-wavelength observations with next-generation telescopes, including systematic large-area surveys, will provide unprecedented opportunities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales. We highlight ways in which these observations will open a new avenue for multi-messenger studies in the dawn of low-frequency (~ nHz - mHz) GW astronomy.

Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts.

Hierarchical models of galaxy formation predict that galaxy mergers represent a significant transitional stage of rapid supermassive black hole (SMBH) growth. Yet, the connection between the merging process and enhanced active galactic nuclei (AGN) activity as well as the timescale of SMBH mergers remains highly uncertain. The breakthrough in reconciling the importance of galaxy mergers with black hole growth lies in a thoroughly-studied census of dual AGN across cosmic history, which will be enabled by next-generation observational capabilities, theoretical advances, and simulations. This white paper outlines the key questions in galaxy mergers, dual and offset AGN, and proposes multiwavelength solutions using future high-resolution observatories in the X-rays (AXIS, Lynx), near and mid-infrared (30 meter class telescopes, JWST), and submillimeter (ALMA).

Intermediate-mass black holes (IMBHs) with masses between $100 - 10^5M_{\odot}$ are crucial to our understanding of black hole seed formation and are the prime targets for LISA, yet black holes in this mass range have eluded detection by traditional optical spectroscopic surveys aimed at finding active galactic nuclei (AGNs). In this paper, we have modeled for the first time the dependence of the optical narrow emission line strengths on the black hole mass of accreting AGN over the range of $100-10^8M_{\odot}$. We show that as the black hole mass decreases, the hardening of the spectral energy distribution from the accretion disk changes the ionization structure of the nebula. The enhanced high energy emission from IMBHs results in a more extended partially ionized zone compared with models for higher mass black holes. This effect produces a net decrease in the predicted [OIII]/H$\beta$ and [NII]/H$\alpha$ emission line ratios. Based on this model, we demonstrate that the standard optical narrow emission line diagnostics used to identify massive black holes fail when black hole mass falls below $\approx10^4M_{\odot}$ for highly accreting IMBHs and for radiatively inefficient IMBHs with active star formation. Our models call into question the ability of common optical spectroscopic diagnostics to confirm AGN candidates in dwarf galaxies, and indicate that the low-mass black hole occupation fraction inferred from such diagnostics will be severely biased.

Major galaxy mergers are thought to play an important part in fuelling the growth of supermassive black holes. However, observational support for this hypothesis is mixed, with some studies showing a correlation between merging galaxies and luminous quasars and others showing no such association. Recent observations have shown that a black hole is likely to become heavily obscured behind merger-driven gas and dust, even in the early stages of the merger, when the galaxies are well separated (5 to 40 kiloparsecs). Merger simulations further suggest that such obscuration and black-hole accretion peaks in the final merger stage, when the two galactic nuclei are closely separated (less than 3 kiloparsecs). Resolving this final stage requires a combination of high-spatial-resolution infrared imaging and high-sensitivity hard-X-ray observations to detect highly obscured sources. However, large numbers of obscured luminous accreting supermassive black holes have been recently detected nearby (distances below 250 megaparsecs) in X-ray observations. Here we report high-resolution infrared observations of hard-X-ray-selected black holes and the discovery of obscured nuclear mergers, the parent populations of supermassive-black-hole mergers. We find that obscured luminous black holes (bolometric luminosity higher than 2x10^44 ergs per second) show a significant (P<0.001) excess of late-stage nuclear mergers (17.6 per cent) compared to a sample of inactive galaxies with matching stellar masses and star formation rates (1.1 per cent), in agreement with theoretical predictions. Using hydrodynamic simulations, we confirm that the excess of nuclear mergers is indeed strongest for gas-rich major-merger hosts of obscured luminous black holes in this final stage.

Gravitational-wave (GW) and gravitational slingshot recoil kicks, which are natural products of SMBH evolution in merging galaxies, can produce active galactic "nuclei" that are offset from the centers of their host galaxies. Detections of offset AGN would provide key constraints on SMBH binary mass and spin evolution and on GW event rates. Although numerous offset AGN candidates have been identified, none have been definitively confirmed. The ngVLA offers unparalleled capabilities to identify and confirm candidate offset AGN from sub-parsec to kiloparsec scales, opening a new avenue for multi-messenger studies in the dawn of low-frequency GW astronomy.

The Next-Generation Very Large Array (ngVLA) has the potential to be a workhorse for the discovery and study of paired supermassive black holes either at large separations (dual) or in tightly bound systems (binary). In this chapter, we outline the science case for the study of these supermassive pairs, and summarize discovery methods that can be used at radio wavelengths to discover them: including morphological, spectral, and time-domain identifications. One critical aspect of this work is that multi-messenger binary black hole studies may be possible with the ngVLA when combined with gravitational-wave searches using pulsar timing array techniques. However, long-baseline interferometery (>>1000 km) will make this possibility more likely by expanding the redshift range at which radio emission arising from two separate black holes may be resolved and studied.

The role of galaxy mergers in fueling active galactic nuclei (AGN) is still debated, owing partly to selection effects inherent to studies of the merger/AGN connection. In particular, luminous AGN are often heavily obscured in late-stage mergers. Mid-infrared (IR) color selection of dust-enshrouded AGN with, e.g., the Wide-field Infrared Survey Explorer (WISE) has uncovered large new populations of obscured AGN. However, this method is sensitive mainly to AGN that dominate emission from the host. To understand how these selection biases affect mid-IR studies of the merger/AGN connection, we simulate the evolution of obscured AGN throughout galaxy mergers. Although mid-IR colors closely trace luminous, obscured AGN, we show that nearly half of merger-triggered AGN are missed with common mid-IR selection criteria, even in late-stage, gas-rich major mergers. At z < 0.5, where merger signatures and dual nuclei can most easily be detected, we find that a more lenient W1-W2 > 0.5 cut greatly improves completeness without significantly decreasing reliability. Extreme nuclear starbursts are briefly able to mimic this AGN signature, but this is largely irrelevant in mergers, where such starbursts are accompanied by AGN. We propose a two-color cut that yields high completeness and reliability even in starbursting systems. Further, we show that mid-IR color selection very effectively identifies dual AGN hosts, with the highest fraction at the smallest separations (< 3 kpc). Thus, many merger hosts of mid-IR AGN should contain unresolved dual AGN; these are ideal targets for high-resolution follow-up, particularly with the James Webb Space Telescope.

Intermediate mass black holes (IMBHs), with masses between 100 to 10^5 M_\odot, represent the link between stellar mass black holes and the supermassive black holes that reside in galaxy centers. While IMBHs are crucial to our understanding of black hole seed formation, black holes of less than \approx 10^4 M_\odot eluded detection by traditional searches. Observations of the infrared coronal lines (CLs) offer us one of the most promising tools to discover IMBHs in galaxies. We have modeled the infrared emission line spectrum that is produced by gas photoionized by an AGN radiation field and explored for the first time the dependence of the infrared CL spectrum on black hole mass over the range of 10^2 M_\odot to 10^8 M_\odot. We show that infrared CLs are expected to be prominent in the spectra of accreting IMBHs and can potentially be a powerful probe of the black hole mass in AGNs. We identify key emission line ratios in the 1-30 \mu m range that are most sensitive to black hole mass. While variations in accretion rate and the physical parameters of the gas can also affect the CL spectrum, we demonstrate that the effect of black hole mass is likely to be the most dramatic over the mass range explored in our models. With the unprecedented sensitivity of JWST, a large number of CLs will be detectable for the first time, providing important insight into the existence and properties of IMBHs in the local universe, potentially revolutionizing our understanding of this class of object.

We calculate the properties, occurrence rates and detection prospects of individually resolvable 'single sources' in the low frequency gravitational wave (GW) spectrum. Our simulations use the population of galaxies and massive black hole binaries from the Illustris cosmological hydrodynamic simulations, coupled to comprehensive semi-analytic models of the binary merger process. Using mock pulsar timing arrays (PTA) with, for the first time, varying red-noise models, we calculate plausible detection prospects for GW single sources and the stochastic GW background (GWB). Contrary to previous results, we find that single sources are at least as detectable as the GW background. Using mock PTA, we find that these 'foreground' sources (also 'deterministic'/'continuous') are likely to be detected with $\sim 20 \, \textrm{yr}$ total observing baselines. Detection prospects, and indeed the overall properties of single sources, are only moderately sensitive to binary evolution parameters---namely eccentricity & environmental coupling, which can lead to differences of $\sim 5 \,\textrm{yr}$ in times to detection. Red noise has a stronger effect, roughly doubling the time to detection of the foreground between a white-noise only model ($\sim 10$ - $15 \, \textrm{yr}$) and severe red noise ($\sim 20$ - $30 \, \textrm{yr}$). The effect of red noise on the GWB is even stronger, suggesting that single source detections may be more robust. We find that typical signal-to-noise ratios for the foreground peak near $f = 0.1 \, \textrm{yr}^{-1}$, and are much less sensitive to the continued addition of new pulsars to PTA.

A direct consequence of hierarchical galaxy formation is the existence of dual supermassive black holes (SMBHs), which may be preferentially triggered as active galactic nuclei (AGN) during galaxy mergers. Despite decades of searching, however, dual AGNs are extremely rare, and most have been discovered serendipitously. Using the all-sky WISE survey, we identified a population of over 100 morphologically identified interacting galaxies or mergers that display red mid-infrared colors often associated in extragalactic sources with powerful AGNs. The vast majority of these advanced mergers are optically classified as star-forming galaxies suggesting that they may represent an obscured population of AGNs that cannot be found through optical studies. In this work, we present Chandra/ACIS observations and near-infrared spectra with the Large Binocular Telescope of six advanced mergers with projected pair separations less than ~ 10 kpc. The combined X-ray, near-infrared, and mid-infrared properties of these mergers provide confirmation that four out of the six mergers host at least one AGN, with four of the mergers possibly hosting dual AGNs with projected separations less than ~10 kpc, despite showing no firm evidence for AGNs based on optical spectroscopic studies. Our results demonstrate that 1) optical studies miss a significant fraction of single and dual AGNs in advanced mergers, and 2) mid-infrared pre-selection is extremely effective in identifying dual AGN candidates in late-stage mergers. Our multi-wavelength observations suggest that the buried AGNs in these mergers are highly absorbed, with intrinsic column densities in excess of N_H >10^24cm^-2, consistent with hydrodynamic simulations.

Pulsar Timing Arrays (PTA) around the world are using the incredible consistency of millisecond pulsars to measure low frequency gravitational waves from (super)Massive Black Hole (MBH) binaries. We use comprehensive MBH merger models based on cosmological hydrodynamic simulations to predict the spectrum of the stochastic Gravitational-Wave Background (GWB). We use real Time-of-Arrival (TOA) specifications from the European, NANOGrav, Parkes, and International PTA (IPTA) to calculate realistic times to detection of the GWB across a wide range of model parameters. In addition to exploring the parameter space of environmental hardening processes (in particular: stellar scattering efficiencies), we have expanded our models to include eccentric binary evolution which can have a strong effect on the GWB spectrum. Our models show that strong stellar scattering and high characteristic eccentricities enhance the GWB strain amplitude near the PTA sensitive "sweet-spot" (near the frequency $f = 1 \, \mathrm{yr}^{-1}$), slightly improving detection prospects in these cases. While the GWB $amplitude$ is degenerate between cosmological and environmental parameters, the location of a spectral turnover at low frequencies ($f \lesssim 0.1 \, \mathrm{yr}^{-1}$) is strongly indicative of environmental coupling. At high frequencies ($f\gtrsim 1 \, \mathrm{yr}^{-1}$), the GWB spectral index can be used to infer the number density of sources and possibly their eccentricity distribution. Even with merger models that use pessimistic environmental and eccentricity parameters, if the current rate of PTA expansion continues, we find that the International PTA is highly likely to make a detection within about 10 years.

Mergers of galaxies are thought to cause significant gas inflows to the inner parsecs, which can activate rapid accretion onto supermassive black holes (SMBHs), giving rise to Active Galactic Nuclei (AGN). During a significant fraction of this process, SMBHs are predicted to be enshrouded by gas and dust. Studying 52 galactic nuclei in infrared-selected local Luminous and Ultra-luminous infrared galaxies in different merger stages in the hard X-ray band, where radiation is less affected by absorption, we find that the amount of material around SMBHs increases during the last phases of the merger. We find that the fraction of Compton-thick (CT, $N_{\rm\,H}\geq 10^{24}\rm\,cm^{-2}$) AGN in late merger galaxies is higher ($f_{\rm\,CT}=65^{+12}_{-13}\%$) than in local hard X-ray selected AGN ($f_{\rm\,CT}=27\pm 4\%$), and that obscuration reaches its maximum when the nuclei of the two merging galaxies are at a projected distance of $D_{12}\simeq0.4-10.8$ kiloparsecs ($f_{\rm\,CT}=77_{-17}^{+13}\%$). We also find that all AGN of our sample in late merger galaxies have $N_{\rm\,H}> 10^{23}\rm\,cm^{-2}$, which implies that the obscuring material covers $95^{+4}_{-8}\%$ of the X-ray source. These observations show that the material is most effectively funnelled from the galactic scale to the inner tens of parsecs during the late stages of galaxy mergers, and that the close environment of SMBHs in advanced mergers is richer in gas and dust with respect to that of SMBHs in isolated galaxies, and cannot be explained by the classical AGN unification model in which the torus is responsible for the obscuration.

We present a combined morphological and X-ray analysis of Was 49, an isolated, dual AGN system notable for the presence of a dominant AGN Was 49b in the disk of the primary galaxy Was 49a, at a projected radial distance of 8 kpc from the nucleus. Using X-ray data from Chandra, NuSTAR, and Swift, we find that this AGN has a bolometric luminosity of L_bol ~ 2 x 10^45 erg/s, with a black hole mass of M_BH=1.3^{+2.9}_{-0.9} x 10^8 M_Sol. Despite its large mass, our analysis of optical data from the Discovery Channel Telescope shows that the supermassive black hole is hosted by a stellar counterpart with a mass of only 5.6^{+4.9}_{-2.6} x 10^9 M_Sol, making the SMBH potentially larger than expected from SMBH-galaxy scaling relations, and the stellar counterpart exhibits a morphology that is consistent with dwarf elliptical galaxies. Our analysis of the system in the r and K bands indicates that Was 49 is a minor merger, with a mass ratio of Was 49a to Was 49b between 1:7 and 1:15. This is in contrast with findings that the most luminous merger-triggered AGNs are found in major mergers, and that minor mergers predominantly enhance AGN activity in the primary galaxy.

Gravitational Waves (GW) have now been detected from stellar-mass black hole binaries, and the first observations of GW from Massive Black Hole (MBH) Binaries are expected within the next decade. Pulsar Timing Arrays (PTA), which can measure the years long periods of GW from MBHB, have excluded many standard predictions for the amplitude of a stochastic GW Background (GWB). We use coevolved populations of MBH and galaxies from hydrodynamic, cosmological simulations ('Illustris') to calculate a predicted GWB. The most advanced predictions so far have included binary hardening mechanisms from individual environmental processes. We present the first calculation including all of the environmental mechanisms expected to be involved: dynamical friction, stellar 'loss-cone' scattering, and viscous drag from a circumbinary disk. We find that MBH binary lifetimes are generally multiple gigayears, and only a fraction coalesce by redshift zero. For a variety of parameters, we find all GWB amplitudes to be below the most stringent PTA upper limit of $A_{\textrm{yr}^{-1}} \approx 10^{-15}$. Our fairly conservative fiducial model predicts an amplitude of $A_{\textrm{yr}^{-1}} \approx 0.4\times 10^{-15}$---less than a factor of three below the current limit. At lower frequencies, we find $A_{0.1\,\textrm{yr}^{-1}} \approx 1.5\times 10^{-15}$ with spectral indices between $-0.4$ and $-0.6$---significantly flatter than the canonical value of $-2/3$ due to purely GW-driven evolution. Typical MBHB driving the GWB signal come from redshifts around $0.3$, with total masses of a few times $10^9\,M_\odot$, and in host galaxies with very large stellar masses. Even without GWB detections, our results can be connected to observations of dual AGN to constrain binary evolution.

Supermassive black hole (BH) mergers produce powerful gravitational wave (GW) emission. Asymmetry in this emission imparts a recoil kick to the merged BH, which can eject the BH from its host galaxy altogether. Recoiling BHs could be observed as offset active galactic nuclei (AGN). Several candidates have been identified, but systematic searches have been hampered by large uncertainties regarding their observability. By extracting merging BHs and host galaxy properties from the Illustris cosmological simulations, we have developed a comprehensive model for recoiling AGN. Here, for the first time, we model the effects of BH spin alignment and recoil dynamics based on the gas-richness of host galaxies. We predict that if BH spins are not highly aligned, seeing-limited observations could resolve offset AGN, making them promising targets for all-sky surveys. For randomly-oriented spins, less than about 10 spatially-offset AGN may be detectable in HST-COSMOS, and > 10^3 could be found with Pan-STARRS, LSST, Euclid, and WFIRST. Nearly a thousand velocity-offset AGN are predicted within the SDSS footprint; the rarity of large broad-line offsets among SDSS quasars is likely due in part to selection effects but suggests that spin alignment plays a role in suppressing recoils. Nonetheless, in our most physically motivated model where alignment occurs only in gas-rich mergers, hundreds of offset AGN should be found in all-sky surveys. Our findings strongly motivate a dedicated search for recoiling AGN.

We present the full public release of all data from the Illustris simulation project. Illustris is a suite of large volume, cosmological hydrodynamical simulations run with the moving-mesh code Arepo and including a comprehensive set of physical models critical for following the formation and evolution of galaxies across cosmic time. Each simulates a volume of (106.5 Mpc)^3 and self-consistently evolves five different types of resolution elements from a starting redshift of z=127 to the present day, z=0. These components are: dark matter particles, gas cells, passive gas tracers, stars and stellar wind particles, and supermassive black holes. This data release includes the snapshots at all 136 available redshifts, halo and subhalo catalogs at each snapshot, and two distinct merger trees. Six primary realizations of the Illustris volume are released, including the flagship Illustris-1 run. These include three resolution levels with the fiducial "full" baryonic physics model, and a dark matter only analog for each. In addition, we provide four distinct, high time resolution, smaller volume "subboxes". The total data volume is ~265 TB, including ~800 full volume snapshots and ~30,000 subbox snapshots. We describe the released data products as well as tools we have developed for their analysis. All data may be directly downloaded in its native HDF5 format. Additionally, we release a comprehensive, web-based API which allows programmatic access to search and data processing tasks. In both cases we provide example scripts and a getting-started guide in several languages: currently, IDL, Python, and Matlab. This paper addresses scientific issues relevant for the interpretation of the simulations, serves as a pointer to published and on-line documentation of the project, describes planned future additional data releases, and discusses technical aspects of the release.

While performing a survey to detect recoiling supermassive black holes, we have identified an unusual source having a projected offset of 800 pc from a nearby dwarf galaxy. The object, SDSS J113323.97+550415.8, exhibits broad emission lines and strong variability. While originally classified as a supernova (SN) because of its nondetection in 2005, we detect it in recent and past observations over 63 yr and find over a magnitude of rebrightening in the last 2 years. Using high-resolution adaptive optics observations, we constrain the source emission region to be <12 pc and find a disturbed host-galaxy morphology indicative of recent merger activity. Observations taken over more than a decade show narrow [O III] lines, constant ultraviolet emission, broad Balmer lines, a constant putative black hole mass over a decade of observations despite changes in the continuum, and optical emission-line diagnostics consistent with an active galactic nucleus (AGN). However, the optical spectra exhibit blueshifted absorption, and eventually narrow Fe II and [Ca II] emission, each of which is rarely found in AGN spectra. While this peculiar source displays many of the observational properties expected of a potential black hole recoil candidate, some of the properties could also be explained by a luminous blue variable star (LBV) erupting for decades since 1950, followed by a Type IIn SN in 2001. Interpreted as an LBV followed by a SN analogous to SN 2009ip, the multi-decade LBV eruptions would be the longest ever observed, and the broad Halpha emission would be the most luminous ever observed at late times (>10 yr), larger than that of unusually luminous supernovae such as SN 1988Z, suggesting one of the most extreme episodes of pre-SN mass loss ever discovered.

The unusual morphologies of the Andromeda spiral galaxy (M31) and its dwarf companion M32 have been characterized observationally in great detail. The two galaxies' apparent proximity suggests that Andromeda's prominent star-forming ring as well as M32's compact elliptical structure may result from a recent collision. Here we present the first self-consistent model of the M31-M32 interaction that simultaneously reproduces observed positions, velocities, and morphologies for both galaxies. Andromeda's spiral structure is resolved in unprecedented detail, showing that a rare head-on orbit is not necessary to match Andromeda's ring-like morphology. The passage of M32 through Andromeda's disk perturbs the disk velocity structure. We find tidal stripping of M32's stars to be inefficient during the interaction, suggesting that some cEs are intrinsically compact. Additionally, the orbital solution implies that M32 is currently closer to the Milky Way than models have typically assumed, a prediction that may be testable with upcoming observations.

We present a first attempt to model the narrow-line (NL) region of active galactic nuclei (AGN) in hydrodynamic simulations of galaxy mergers, using a novel physical prescription. This model is used to determine the origin of double-peaked NL (dNL) AGN in merging galaxies and their connection to supermassive black hole (SMBH) pairs, motivated by recent observations of such objects. We find that dNL AGN induced by the relative motion of SMBH pairs are a generic but short-lived feature of gaseous major mergers. dNL AGN should often be observed in late-stage mergers, during the kpc-scale phase of SMBH inspiral or soon after the SMBH merger. However, even within the kpc-scale phase, only a minority of dNL AGN are directly induced by SMBH motion; their lifetimes are typically a few Myr. Most double peaks arise from gas kinematics near the SMBH, although prior to the SMBH merger up to 80% of all dNL profiles may be influenced by SMBH motion via altered peak ratios or velocity offsets. The total lifetimes of dNL AGN depend strongly on viewing angle and on properties of the merging galaxies. Also, in a typical merger, at least 10-40% of the double peaks induced by SMBH motion have small projected separations, 0.1-1 kpc, such that dual peaks of stellar surface brightness are not easily resolved. Diffuse tidal features can indicate late-stage galaxy mergers, although they do not distinguish SMBH pairs from merged SMBHs. We show that dNL profiles with peak velocity splittings > 500 km s^-1 or with measurable overall velocity shifts are often associated with SMBH pairs. Our results support the notion that selection of dNL AGN is a promising method for identifying dual SMBH candidates, but demonstrate the critical importance of high-resolution, multi-wavelength follow-up observations, and the use of multiple lines of evidence, for confirming the dual nature of candidate SMBH pairs. (Abridged)

The galaxy CXOC J100043.1+020637, also known as CID-42, is a highly unusual object. An apparent galaxy merger remnant, it displays signatures of both an inspiraling, kiloparsec-scale active galactic nucleus (AGN) pair and of a recoiling AGN with a kick velocity > 1300 km s^-1. Among recoiling AGN candidates, CID-42 alone has both spatial offsets (in optical and X-ray bands) and spectroscopic offsets. In order to constrain the relative likelihood of both scenarios, we develop models using hydrodynamic galaxy merger simulations coupled with radiative transfer calculations. Our gas-rich, major merger models are generally well matched to the galactic morphology and to the inferred stellar mass and star formation rate. We show that a recoiling supermassive black hole (SMBH) in CID-42 should be observable as an AGN at the time of observation. However, in order for the recoiling AGN to produce narrow-line emission, it must be observed shortly after the kick while it still inhabits a dense gaseous region, implying a large total kick velocity (v_k > 2000 km s^-1). For the dual AGN scenario, an unusually large broad-line offset is required, and the best match to the observed morphology requires a galaxy that is less luminous than CID-42. Further, the lack of X-ray emission from one of the two optical nuclei is not easily attributed to an intrinsically quiescent SMBH or to a Compton-thick galactic environment. While the current data do not allow either the recoiling or the dual AGN scenario for CID-42 to be excluded, our models highlight the most relevant parameters for distinguishing these possibilities with future observations. In particular, high-quality, spatially-resolved spectra that can pinpoint the origin of the broad and narrow line features will be critical for determining the nature of this unique source.

We present Chandra High Resolution Camera observations of CID-42, a candidate recoiling supermassive black hole (SMBH) at z=0.359 in the COSMOS survey. CID-42 shows two optical compact sources resolved in the HST/ACS image embedded in the same galaxy structure and a velocity offset of ~1300 km/s between the H\beta\ broad and narrow emission line, as presented by Civano et al. (2010). Two scenarios have been proposed to explain the properties of CID-42: a GW recoiling SMBH and a double Type 1/ Type 2 AGN system, where one of the two is recoiling because of slingshot effect. In both scenario, one of the optical nuclei hosts an unobscured AGN, while the other one, either an obscured AGN or a star forming compact region. The X-ray Chandra data allow to unambiguously resolve the X-ray emission, and unveil the nature, of the two optical sources in CID-42. We find that only one of the optical nuclei is responsible for the whole X-ray unobscured emission observed and a 3sigma upper limit on the flux of the second optical nucleus is measured. The upper limit on the X-ray luminosity plus the analysis of the multiwavelength spectral energy distribution indicate the presence of a starforming region in the second source rather than an obscured SMBH, thus favoring the GW recoil scenario. However the presence of a very obscured SMBH cannot be fully ruled-out. A new X-ray feature, in a SW direction with respect to the main source, is discovered and discussed.

Central supermassive black holes (SMBHs) are a ubiquitous feature of locally-observed galaxies, and ample evidence suggests that the growth of SMBHs and their host galaxies is closely linked. However, in the event of a merger, gravitational-wave (GW) recoil may displace a SMBH from its galactic center, or eject it entirely. To explore the consequences of this phenomenon, we use hydrodynamic simulations of gaseous galaxy mergers that include a range of BH recoil velocities. We have generated a suite of over 200 simulations with more than 60 merger models, enabling us to identify systematic trends in the behavior of recoiling BHs -- specifically (i) their dynamics, (ii) their observable signatures, and (iii) their effects on BH/galaxy co-evolution. (i) Recoiling BH trajectories depend heavily on the gas content of the host galaxy; maximal BH displacements from the center may vary by up to an order of magnitude between gas-rich and gas-poor mergers. In some cases, recoil trajectories also depend on the timing of the BH merger relative to the formation of the galaxy merger remnant. (ii) Recoiling BHs may be observable as offset active galactic nuclei (AGN) via either kinematic offsets (v > 800 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of about 1 - 100 Myr. In addition, recoil events affect the total AGN lifetime. GW recoil generally reduces the lifetimes of bright AGN, but may extend lower-luminosity AGN lifetimes. (iii) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. Furthermore, recoil events displace AGN feedback from the galactic center, which enhances central star formation rates. This results in longer starburst phases and higher central stellar densities in merger remnants.

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. With our suite of over 200 merger simulations, we identify systematic trends in the behavior of recoiling BHs. Our main results are as follows. (1) While BHs kicked at nearly the central escape speed (vesc) are essentially "lost" to the galaxy, in gas rich mergers, BHs kicked with up to about 0.7 vesc may be confined to the central few kpc of the galaxy. (2) The inflow of cold gas during a gas-rich major merger may cause a rapid increase in central escape speed; in such cases recoil trajectories will depend on the timing of the BH merger relative to the change in vesc. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN) but may actually extend AGN lifetimes at lower luminosities. (4) Recoiling AGN may be observable via kinematic offsets (v > 500 km s^-1) or spatial offsets (R > 1 kpc) for lifetimes of up to about 10 - 100 Myr. (5) Rapidly-recoiling BHs may be up to about 5 times less massive than their stationary counterparts. These mass deficits lower the normalization of the M - sigma relation and contribute to both intrinsic and overall scatter. (6) Finally, the displacement of AGN feedback by a recoil event causes higher central star formation rates in the merger remnant, thereby extending the starburst phase of the merger and creating a denser, more massive stellar cusp.

Gravitational-wave (GW) recoil of merging supermassive black holes (SMBHs) may influence the co-evolution of SMBHs and their host galaxies. We examine this possibility using SPH/N-body simulations of gaseous galaxy mergers in which the merged BH receives a recoil kick. This enables us to follow recoiling BHs in self-consistent, evolving merger remnants. In contrast to recent studies on similar topics, we conduct a large parameter study, generating a suite of over 200 simulations with more than 60 merger models and a range of recoil velocities (vk). Our main results are as follows. (1) BHs kicked at nearly the central escape speed (vesc) may oscillate on large orbits for up to a Hubble time, but in gas-rich mergers, BHs kicked with up to ~ 0.7 vesc may be confined to the central few kpc of the galaxy, owing to gas drag and steep central potentials. (2) vesc in gas-rich mergers may increase rapidly during final coalescence, in which case trajectories may depend on the timing of the BH merger relative to the formation of the potential well. (3) Recoil events generally reduce the lifetimes of bright active galactic nuclei (AGN), but may actually extend AGN lifetimes at lower luminosities. (4) Kinematically-offset AGN (v > 800 km s^-1) may be observable for up to ~ 10 Myr either immediately after the recoil or during pericentric passages through a gas-rich remnant. (5) Spatially-offset AGN (R > 1 kpc) generally have low luminosities and lifetimes of ~ 1 - 100 Myr. (6) Rapidly-recoiling BHs may be up to ~ 5 times less massive than their stationary counterparts. This lowers the normalization of the M-sigma relation and contributes to both intrinsic and overall scatter. (7) Finally, the displacement of AGN feedback after a recoil event enhances central star formation rates, thereby extending the starburst phase of the merger and creating a denser stellar cusp. [Abridged.]

We present a detailed study of a peculiar source in the COSMOS survey at z=0.359. Source CXOCJ100043.1+020637 (CID-42) presents two compact optical sources embedded in the same galaxy. The distance between the 2, measured in the HST/ACS image, is 0.495" that, at the redshift of the source, corresponds to a projected separation of 2.46 kpc. A large (~1200 km/s) velocity offset between the narrow and broad components of Hbeta has been measured in three different optical spectra from the VLT/VIMOS and Magellan/IMACS instruments. CID-42 is also the only X-ray source having in its X-ray spectra a strong redshifted broad absorption iron line, and an iron emission line, drawing an inverted P-Cygni profile. The Chandra and XMM data show that the absorption line is variable in energy by 500 eV over 4 years and that the absorber has to be highly ionized, in order not to leave a signature in the soft X-ray spectrum. That these features occur in the same source is unlikely to be a coincidence. We envisage two possible explanations: (1) a gravitational wave recoiling black hole (BH), caught 1-10 Myr after merging, (2) a Type 1/ Type 2 system in the same galaxy where the Type 1 is recoiling due to slingshot effect produced by a triple BH system. The first possibility gives us a candidate gravitational waves recoiling BH with both spectroscopic and imaging signatures. In the second case, the X-ray absorption line can be explained as a BAL-like outflow from the foreground nucleus (a Type 2 AGN) at the rearer one (a Type 1 AGN), which illuminates the otherwise undetectable wind, giving us the first opportunity to show that fast winds are present in obscured AGN.

Simulations of binary black hole mergers indicate that asymmetrical gravitational wave (GW) emission can cause black holes to recoil at speeds up to thousands of km/s. These GW recoil events can dramatically affect the coevolution of recoiling supermassive black holes (SMBHs) and their host galaxies. However, theoretical studies of SMBH-galaxy evolution almost always assume a stationary central black hole. In light of the numerical results on GW recoil velocities, we relax that assumption here and consider the consequences of recoil for SMBH evolution. We follow the trajectories of SMBHs ejected in a smooth background potential that includes both a stellar bulge and a multi-component gaseous disk. In addition, we calculate the accretion rate onto the SMBH as a function of time using a hybrid prescription of viscous (alpha-disk) and Bondi accretion. We find that recoil kicks between 100 km/s and the escape speed cause SMBHs to wander through the galaxy and halo for about 1 Myr - 1 Gyr before settling back to the galactic center. However, the mass accreted during this time is roughly constant at about 10% of the initial mass, independent of the recoil velocity. This indicates that while large recoils may disrupt active galactic nuclei feedback processes, recoil itself is an effective means of regulating SMBH growth. Recoiling SMBHs may be observable as spatially or kinematically offset quasars, but finding such systems could be challenging, because the largest offsets correspond to the shortest quasar lifetimes.

While many observed ultra-luminous X-ray sources (ULXs, Lx > 10^39 erg s^-1) could be extragalactic X-ray binaries (XRBs) emitting close to the Eddington limit, the highest-luminosity ULXs (Lx > 3x10^39 erg s^-1) exceed the isotropic Eddington luminosity for even high-stellar-mass accreting black hole XRBs. It has been suggested that these highest-luminosity ULXs may contain accreting intermediate-mass black hole (IMBH) binaries. We consider this hypothesis for dense, young (about 100 Myr) stellar clusters where we assume that a 50-500 solar mass central IMBH has formed through runaway growth of a massive star. Using numerical simulations of the dynamics and evolution of the central black hole's captured companions, we obtain estimates of the incidence of mass transfer phases and possible ULX activity throughout the IMBH's evolutionary history. We find that, although it is common for the central black hole to acquire binary companions, there is a very low probability that these interacting binaries will become observable ULX sources.