We present a novel implementation of supermassive black hole (SMBH) formation, dynamics, and accretion in the massively parallel tree+SPH code, ChaNGa. This approach improves the modeling of SMBHs in fully cosmological simulations, allowing for a more detailed analysis of SMBH-galaxy co-evolution throughout cosmic time. Our scheme includes novel, physically motivated models for SMBH formation, dynamics and sinking timescales within galaxies, and SMBH accretion of rotationally supported gas. The sub-grid parameters that regulate star formation (SF) and feedback from SMBHs and SNe are optimized against a comprehensive set of z = 0 galaxy scaling relations using a novel, multi-dimensional parameter search. We have incorporated our new SMBH implementation and parameter optimization into a new set of high resolution, large-scale cosmological simulations called Romulus. We present initial results from our flagship simulation, Romulus25, showing that our SMBH model results in SF efficiency, SMBH masses, and global SF and SMBH accretion histories at high redshift that are consistent with observations. We discuss the importance of SMBH physics in shaping the evolution of massive galaxies and show how SMBH feedback is much more effective at regulating star formation compared to SNe feedback in this regime. Further, we show how each aspect of our SMBH model impacts this evolution compared to more common approaches. Finally, we present a science application of this scheme studying the properties and time evolution of an example dual AGN system, highlighting how our approach allows simulations to better study galaxy interactions and SMBH mergers in the context of galaxy-BH co-evolution.

We study supermassive black holes (BHs) in merging galaxies, using a suite of hydrodynamical simulations with very high spatial (∼10 pc) and temporal (∼1 Myr) resolution, where we vary the initial mass ratio, the orbital configuration, and the gas fraction. (i) We address the question of when and why, during a merger, increased BH accretion occurs, quantifying gas inflows and BH accretion rates. (ii) We also quantify the relative effectiveness in inducing AGN activity of merger-related versus secularrelated causes, by studying different stages of the encounter: the stochastic (or early) stage, the (proper) merger stage, and the remnant (or late) stage. (iii) We assess which galaxy mergers preferentially enhance BH accretion, finding that the initial mass ratio is the most important factor. (iv) We study the evolution of the BH masses, finding that the BH mass contrast tends to decrease in minor mergers and to increase in major mergers. This effect hints at the existence of a preferential range of mass ratios for BHs in the final pairing stages. (v) In both merging and dynamically quiescent galaxies, the gas accreted by the BH is not necessarily the gas with low angular momentum, but the gas that loses angular momentum.

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 exert torques on the embedded masses causing them to migrate through the disk in a manner analogous to 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-300R g , where R 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.

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 ∼ 10 −4 − 10 4 Gpc −1 yr −1 , so upper limits from LIGO (< 212Gpc −1 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.

Supermassive black holes (SMBHs) have been detected in the centers of most nearby massive galaxies. Galaxies today are not only the products of billions of years of galaxy mergers, but also billions of years of SMBH activity as active galactic nuclei (AGNs) that is connected to galaxy mergers. In this context, detection of AGN pairs should be relatively common. Observationally, however, dual AGNs are scant, being just a few percent of all AGNs. In this Letter, we investigate the triggering of AGN activity in merging galaxies via a suite of highresolution hydrodynamical simulations. We follow the dynamics and accretion onto the SMBHs as they move from separations of tens of kiloparsecs to tens of parsecs. Our resolution, cooling, and star formation implementation produce an inhomogeneous, multi-phase interstellar medium, allowing us to accurately trace star formation and accretion onto the SMBHs. We study the impact of gas content, morphology, and mass ratio, focusing on AGN activity and dynamics across a wide range of relevant conditions. We test when the two AGNs are simultaneously detectable, for how long and at which separations. We find that strong dual AGN activity occurs during the late phases of the mergers, at small separations (<1-10 kpc) below the resolution limit of most surveys. Much of the SMBH accretion is not simultaneous, limiting the dual AGN fraction detectable through imaging and spectroscopy to a few percent, in agreement with observational samples.

Recent discoveries of massive black holes (MBHs) in dwarf galaxies suggest that they may have a more common presence than once thought. Systematic searches are revealing more candidates, but this process could be accelerated by predictions from simulations. We perform a study of several high-resolution, cosmological, zoom-in simulations focusing on dwarf galaxies that host massive black holes at z = 0, with the aim of determining when the black holes are most observable. Larger dwarf galaxies are more likely to host MBHs than those of lower mass. About 50% of the MBHs in dwarfs are not centrally located, but rather are wandering within a few kpc of the galaxy center. The accretion luminosities of MBHs in dwarfs are low throughout cosmic time, rendering them extremely difficult to detect. However, the merger history of these MBHs is optimal for gravitational wave detection by LISA.

We perform a suite of high-resolution smoothed particle hydrodynamics simulations to investigate the orbital decay and mass evolution of massive black hole (MBH) pairs down to scales of ∼ 30 pc during minor mergers of disk galaxies. Our simulation set includes star formation and accretion onto the MBHs, as well as feedback from both processes. We consider 1:10 merger events starting at z ∼ 3, with MBH masses in the sensitivity window of the Laser Interferometer Space Antenna, and we follow the coupling between the merger dynamics and the evolution of the MBH mass ratio until the satellite galaxy is tidally disrupted. While the more massive MBH accretes in most cases as if the galaxy were in isolation, the satellite MBH may undergo distinct episodes of enhanced accretion, owing to strong tidal torques acting on its host galaxy and to orbital circularization inside the disk of the primary galaxy. As a consequence, the initial 1:10 mass ratio of the MBHs changes by the time the satellite is disrupted. Depending on the initial fraction of cold gas in the galactic disks and the geometry of the encounter, the mass ratios of the MBH pairs at the time of satellite disruption can stay unchanged or become as large as 1:2. Remarkably, the efficiency of MBH orbital decay correlates with the final mass ratio of the pair itself: MBH pairs that increase significantly their mass ratio are also expected to inspiral more promptly down to nuclear-scale separations. These findings indicate that the mass ratios of MBH pairs in galactic nuclei do not necessarily trace the mass ratios of their merging host galaxies, but are determined by the complex interplay between gas accretion and merger dynamics.

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 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 GW to include gravitational wave emission between scattering events. For astrophysically plausible models we find that typically N GW 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.