Showing 1–50 of 64 results for author: Bellovary, J

An Enhanced Massive Black Hole Occupation Fraction Predicted in Cluster Dwarf Galaxies

Authors: Michael Tremmel, Angelo Ricarte, Priyamvada Natarajan, Jillian Bellovary, Ramon Sharma, Thomas R. Quinn

Abstract: The occupation fraction of massive black holes (MBHs) in dwarf galaxies offers interesting insights into initial black hole seeding mechanisms and their mass assembly history, though disentangling these two effects remains challenging. Using the {\sc Romulus} cosmological simulations we examine the impact of environment on the occupation fraction of MBHs in low mass galaxies. Unlike most modern co… ▽ More The occupation fraction of massive black holes (MBHs) in dwarf galaxies offers interesting insights into initial black hole seeding mechanisms and their mass assembly history, though disentangling these two effects remains challenging. Using the {\sc Romulus} cosmological simulations we examine the impact of environment on the occupation fraction of MBHs in low mass galaxies. Unlike most modern cosmological simulations, {\sc Romulus} seeds MBHs based on local gas properties, selecting dense ($n>3$ cm$^{-3}$), pristine ($Z<3e-4Z_{\odot}$), and rapidly collapsing regions in the early Universe as sites to host MBHs without assuming anything about MBH occupation as a function of galaxy stellar mass, or halo mass, {\it a priori}. The simulations predict that dwarf galaxies with M$_{\star}<10^9$ M$_{\odot}$ in cluster environments are $\sim2$ times more likely to host a MBH compared to those in the field. The predicted occupation fractions are remarkably consistent with those of nuclear star clusters. Across cluster and field environments, dwarf galaxies with earlier formation times are more likely to host a MBH. While the MBH occupation function is similar between cluster and field environments at high redshift ($z>3$), a difference arises as late-forming dwarfs -- which do not exist in the cluster environment -- begin to dominate in the field and pull the MBH occupation fraction down for low mass galaxies. Additionally, prior to in-fall some cluster dwarfs are similar to progenitors of massive, isolated galaxies, indicating that they might have grown to higher masses had they not been impeded by the cluster environment. While the population of MBHs in dwarf galaxies is already widely understood to be important for understanding MBH formation, this work demonstrates that environmental dependence is important to consider as future observations search for low mass black holes in dwarf galaxies. △ Less

Submitted 15 April, 2024; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: 16 pages, 7 figures, Accepted to the Open Journal of Astrophysics

AGN quenching in simulated dwarf galaxies

Authors: Ray S. Sharma, Alyson M. Brooks, Michael Tremmel, Jillian Bellovary, Thomas R. Quinn

Abstract: We examine the quenching characteristics of $328$ isolated dwarf galaxies $\left(10^{8} < M_{\rm star}/M_\odot < 10^{10} \right)$ within the \Rom{} cosmological hydrodynamic simulation. Using mock observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological sim… ▽ More We examine the quenching characteristics of $328$ isolated dwarf galaxies $\left(10^{8} < M_{\rm star}/M_\odot < 10^{10} \right)$ within the \Rom{} cosmological hydrodynamic simulation. Using mock observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological simulations, we find a population of quenched, isolated dwarf galaxies below $M_{\rm star} < 10^{9} M_\odot$ not detected within the NSA. We find that the presence of massive black holes (MBHs) in \Rom{} is largely responsible for the quenched, isolated dwarfs, while isolated dwarfs without an MBH are consistent with quiescent fractions observed in the field. Quenching occurs between $z=0.5-1$, during which the available supply of star-forming gas is heated or evacuated by MBH feedback. Mergers or interactions seem to play little to no role in triggering the MBH feedback. At low stellar masses, $M_{\rm star} \lesssim 10^{9.3} M_\odot$, quenching proceeds across several Gyr as the MBH slowly heats up gas in the central regions. At higher stellar masses, $M_{\rm star} \gtrsim 10^{9.3} M_\odot$, quenching occurs rapidly within $1$ Gyr, with the MBH evacuating gas from the central few kpc of the galaxy and driving it to the outskirts of the halo. Our results indicate the possibility of substantial star formation suppression via MBH feedback within dwarf galaxies in the field. On the other hand, the apparent over-quenching of dwarf galaxies due to MBH suggests higher resolution and/or better modeling is required for MBHs in dwarfs, and quenched fractions offer the opportunity to constrain current models. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 23 pages, 11 figures; submitted to ApJ

Aligning Retrograde Nuclear Cluster Orbits with an Active Galactic Nucleus Accretion Disc

Authors: Syeda S. Nasim, Gaia Fabj, Freddy Caban, Amy Secunda, K. E. Saavik Ford, Barry McKernan, Jillian M. Bellovary, Nathan W. C. Leigh, Wladimir Lyra

Abstract: Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $i<90^{\circ}$) decrease inclination, as well as semi-major axis $(a)$ and eccentricity $(e)$ until orbital alignment with the gas disc ('disc capture'). Captured stellar-origin black holes (sBH) add to the embedded AGN population w… ▽ More Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination $i<90^{\circ}$) decrease inclination, as well as semi-major axis $(a)$ and eccentricity $(e)$ until orbital alignment with the gas disc ('disc capture'). Captured stellar-origin black holes (sBH) add to the embedded AGN population which drives sBH-sBH mergers detectable in gravitational waves using LIGO-Virgo-KAGRA (LVK) or sBH-SMBH mergers detectable with LISA (Laser Interferometer Space Antenna). Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive 'immortal' stars. Here, we investigate the behaviour of polar and retrograde orbiters $(i \geq 90^{\circ})$ interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation $(i<90^{\circ})$ during capture, and decrease $a$ dramatically towards the SMBH. For sBH, we find a critical angle $i_{\rm ret} \sim 113^{\circ}$, below which retrograde sBH decay towards embedded prograde orbits $(i \to 0^{\circ})$, while for $i_{\rm o}>i_{\rm ret}$ sBH decay towards embedded retrograde orbits $(i \to 180^{\circ})$. sBH near polar orbits $(i \sim 90^{\circ})$ and stars on nearly embedded retrograde orbits $(i \sim 180^{\circ})$ show the greatest decreases in $a$. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial $a$. For sBH, disc capture-time is longest for polar orbits, low mass sBH and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the embedded sBH spin distribution. △ Less

Submitted 17 May, 2023; v1 submitted 19 July, 2022; originally announced July 2022.

Comments: 9 pages, 7 figures, 1 table

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 4, July 2023, Pages 5393-5401

Dwarf AGNs from Variability for the Origins of Seeds (DAVOS): Intermediate-mass black hole demographics from optical synoptic surveys

Authors: Colin J. Burke, Yue Shen, Xin Liu, Priyamvada Natarajan, Neven Caplar, Jillian M. Bellovary, Z. Franklin Wang

Abstract: We present a phenomenological forward Monte Carlo model for forecasting the population of active galactic nuclei (AGNs) in dwarf galaxies observable via their optical variability. Our model accounts for expected changes in the spectral energy distribution of AGNs in the intermediate-mass black hole (IMBH) mass range and uses observational constraints on optical variability as a function of black h… ▽ More We present a phenomenological forward Monte Carlo model for forecasting the population of active galactic nuclei (AGNs) in dwarf galaxies observable via their optical variability. Our model accounts for expected changes in the spectral energy distribution of AGNs in the intermediate-mass black hole (IMBH) mass range and uses observational constraints on optical variability as a function of black hole (BH) mass to generate mock light curves. Adopting several different models for the BH occupation function, including one for off-nuclear IMBHs, we quantify differences in the predicted local AGN mass and luminosity functions in dwarf galaxies. As a result, we are able to model the variable fraction of AGNs as a function of physical host properties, such as host galaxy stellar mass, in the presence of complex selection effects. We find that our adopted occupation fractions for the "heavy" and "light" initial BH seeding scenarios can be distinguished with variability data at the $2-3 σ$ level for galaxy host stellar masses below $\sim 10^8 M_\odot$ with the Vera C. Rubin Observatory. We demonstrate the prevalence of a selection bias whereby recovered IMBH masses fall, on average, above the predicted value from the local host galaxy - BH mass scaling relation with the strength of the bias dependent on the survey sensitivity. The methodology developed in this work can be used more broadly to forecast and correct for selection effects for AGN demographic studies in synoptic surveys. Finally, we show that a targeted $\sim$ hourly cadence program over a few nights with the Rubin Observatory can provide strong constraints on IMBH masses given their expected rapid variability timescales. △ Less

Submitted 10 August, 2023; v1 submitted 8 July, 2022; originally announced July 2022.

Comments: 26 pages, 16 figures incl. 5 appendices; minor revisions and typo corrections to match published version in MNRAS

Journal ref: MNRAS, 518, 1880 (2023)

RESOLVE and ECO: Finding Low-Metallicity $z\sim0$ Dwarf AGN Candidates Using Optimized Emission-Line Diagnostics

Authors: Mugdha S. Polimera, Sheila J. Kannappan, Chris T. Richardson, Ashley S. Bittner, Carlynn Ferguson, Amanda J. Moffett, Kathleen D. Eckert, Jillian M. Bellovary, Mark A. Norris

Abstract: Existing star-forming vs. active galactic nucleus (AGN) classification schemes using optical emission-line diagnostics mostly fail for low-metallicity and/or highly star-forming galaxies, missing AGN in typical $z\sim0$ dwarfs. To recover AGN in dwarfs with strong emission lines (SELs), we present a classification scheme optimizing the use of existing optical diagnostics. We use SDSS emission-line… ▽ More Existing star-forming vs. active galactic nucleus (AGN) classification schemes using optical emission-line diagnostics mostly fail for low-metallicity and/or highly star-forming galaxies, missing AGN in typical $z\sim0$ dwarfs. To recover AGN in dwarfs with strong emission lines (SELs), we present a classification scheme optimizing the use of existing optical diagnostics. We use SDSS emission-line catalogs overlapping the volume- and mass-limited RESOLVE and ECO surveys to determine the AGN percentage in SEL dwarfs. Our photoionization grids show that the [O III]/H$β$ versus [S II]/H$α$ diagram (SII plot) and [O III]/H$β$ versus [O I]/H$α$ diagram (OI plot) are less metallicity sensitive and more successful in identifying dwarf AGN than the popular [O III]/H$β$ versus [N II]/H$α$ diagnostic (NII plot or "BPT diagram"). We identify a new category of "star-forming AGN" (SF-AGN) classified as star-forming by the NII plot but as AGN by the SII and/or OI plots. Including SF-AGN, we find the $z\sim0$ AGN percentage in dwarfs with SELs to be $\sim$3-16\%, far exceeding most previous optical estimates ($\sim$1\%). The large range in our dwarf AGN percentage reflects differences in spectral fitting methodologies between catalogs. The highly complete nature of RESOLVE and ECO allows us to normalize strong emission-line galaxy statistics to the full galaxy population, reducing the dwarf AGN percentage to $\sim$0.6-3.0\%. The newly identified SF-AGN are mostly gas-rich dwarfs with halo mass $ < 10^{11.5} M_\odot$, where highly efficient cosmic gas accretion is expected. Almost all SF-AGN also have low metallicities (Z $\lesssim 0.4$ Z$_\odot$), demonstrating the advantage of our method. △ Less

Submitted 20 May, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 31 pages, 16 figures, 3 tables. Accepted by ApJ on April 3 2022

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

A hidden population of massive black holes in simulated dwarf galaxies

Authors: Ray S. Sharma, Alyson M. Brooks, Michael Tremmel, Jillian Bellovary, Angelo Ricarte, Thomas R. Quinn

Abstract: We explore the characteristics of actively accreting MBHs within dwarf galaxies in the \textsc{Romulus25} cosmological hydrodynamic simulation. We examine the MBH occupation fraction, x-ray active fractions, and AGN scaling relations within dwarf galaxies of stellar mass $10^{8} < M_{\rm star} < 10^{10} M_\odot$ out to redshift $z=2$. In the local universe, the MBH occupation fraction is consisten… ▽ More We explore the characteristics of actively accreting MBHs within dwarf galaxies in the \textsc{Romulus25} cosmological hydrodynamic simulation. We examine the MBH occupation fraction, x-ray active fractions, and AGN scaling relations within dwarf galaxies of stellar mass $10^{8} < M_{\rm star} < 10^{10} M_\odot$ out to redshift $z=2$. In the local universe, the MBH occupation fraction is consistent with observed constraints, dropping below unity at $M_{\rm star} < 3\times10^{10} M_{\odot}$, $M_{\rm 200} < 3\times10^{11} M_\odot$. Local dwarf AGN in \textsc{Romulus25} follow observed scaling relations between AGN x-ray luminosity, stellar mass, and star formation rate, though they exhibit slightly higher active fractions and number densities than comparable x-ray observations. Since $z=2$, the MBH occupation fraction has decreased, the population of dwarf AGN has become overall less luminous, and as a result, the overall number density of dwarf AGN has diminished. We predict the existence of a large population of MBHs in the local universe with low x-ray luminosities and high contamination from x-ray binaries and the hot interstellar medium that are undetectable by current x-ray surveys. These hidden MBHs make up $76\%$ of all MBHs in local dwarf galaxies, and include many MBHs that are undermassive relative to their host galaxy's stellar mass. Their detection relies not only on greater instrument sensitivity but on better modeling of x-ray contaminants or multi-wavelength surveys. Our results indicate dwarf AGN were substantially more active in the past despite being low-luminosity today, and indicate future deep x-ray surveys may uncover many hidden MBHs in dwarf galaxies out to at least $z=2$. △ Less

Submitted 1 August, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Optical and JWST Mid-IR Emission Line Diagnostics for Simultaneous IMBH and Stellar Excitation in z~0 Dwarf Galaxies

Authors: C. T. Richardson, C. Simpson, M. S. Polimera, S. J. Kannappan, J. M. Bellovary, C. Greene, S. Jenkins

Abstract: Current observational facilities have yet to conclusively detect $10^3 - 10^4 M_{\odot}$ intermediate mass black holes (IMBHs) that fill in the evolutionary gap between early universe seed black holes and $z \sim 0$ supermassive black holes. Dwarf galaxies present an opportunity to reveal active IMBHs amidst persistent star formation. We introduce photoionization simulations tailored to address ke… ▽ More Current observational facilities have yet to conclusively detect $10^3 - 10^4 M_{\odot}$ intermediate mass black holes (IMBHs) that fill in the evolutionary gap between early universe seed black holes and $z \sim 0$ supermassive black holes. Dwarf galaxies present an opportunity to reveal active IMBHs amidst persistent star formation. We introduce photoionization simulations tailored to address key physical uncertainties: coincident vs. non-coincident mixing of IMBH and starlight excitation, open vs. closed surrounding gas cloud geometries, and different AGN SED shapes. We examine possible AGN emission line diagnostics in the optical and mid-IR, and find that the diagnostics are often degenerate with respect to the investigated physical uncertainties. In spite of these setbacks, and in contrast to recent work, we are able to show that [O III]/H$β$ typically remains bright for dwarf AGN powered by IMBHs down to $10^3 M_{\odot}$. Dwarf AGN are predicted to have inconsistent star-forming and Seyfert/LINER classifications using the most common optical diagnostics. In the mid-IR, [O IV] 25.9$μ$m and [Ar II] 6.98$μ$m are less sensitive to physical uncertainties than are optical diagnostics. Based on these emission lines, we provide several mid-IR emission line diagnostic diagrams with demarcations for separating starbursts and AGN with varying levels of activity. The diagrams are valid over a wide range of ionization parameters and metallicities out to $z\sim0.1$, so will prove useful for future JWST observations of local dwarf AGN in the search for IMBHs. We make our photoionization simulation suite freely available. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: 24 pages, 13 figures, accepted to ApJ

Merging Black Holes in Dwarf Galaxies: Calculating Binary Black Hole Coalescence Timescales from Simulations for LISA Detection

Authors: Victoria I. De Cun, Jillian M. Bellovary, Michael L. Katz

Abstract: Supermassive black holes (SMBHs) merging in dwarf galaxies will be detectable by the Laser Interferometer Space Antenna (LISA) in the mid-2030s. Previous cosmological hydrodynamic simulations have shown the prediction of massive black holes merging in dwarf galaxies, but these simulations are limited by their resolution and cannot follow black hole pairs all the way to coalescence. We calculate th… ▽ More Supermassive black holes (SMBHs) merging in dwarf galaxies will be detectable by the Laser Interferometer Space Antenna (LISA) in the mid-2030s. Previous cosmological hydrodynamic simulations have shown the prediction of massive black holes merging in dwarf galaxies, but these simulations are limited by their resolution and cannot follow black hole pairs all the way to coalescence. We calculate the delay time between black hole pairing and merger based on the properties of the black holes and their host galaxies, and use these properties to calculate gravitational wave strains for eleven different binary black holes that merge inside dwarf galaxies from eight cosmological simulations. This delay time calculation accounts for dynamical friction due to gas and stars, loss-cone scattering, and hardening of the binary due to gravitational radiation. Out of the eleven black hole mergers in the simulations, five black hole pairs will merge within 0.8 - 8 Gyr of forming a close pair and could be observed by LISA, and the remaining six are unresolved due to resolution limitations of the simulation. As all five of the resolved close pairs merge within a Hubble time, we make the broad estimate that close SMBH pairs in dwarf galaxies will merge and be detectable by LISA, but this estimate depends on either the presence of gas during orbital decay or a solution to the dynamical buoyancy problem in cored potentials. △ Less

Submitted 15 February, 2023; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: 8 pages, 3 figures

The Origins of Off-Centre Massive Black Holes in Dwarf Galaxies

Authors: Jillian M. Bellovary, Sarra Hayoune, Katheryn Chafla, Donovan Vincent, Alyson Brooks, Charlotte Christensen, Ferah Munshi, Michael Tremmel, Thomas R. Quinn, Jordan Van Nest, Serena K. Sligh, Michelle Luzuriaga

Abstract: Massive black holes often exist within dwarf galaxies, and both simulations and observations have shown that a substantial fraction of these may be off-center with respect to their hosts. We trace the evolution of off-center massive black holes (MBHs) in dwarf galaxies using cosmological hydrodynamical simulations, and show that the reason for off-center locations is mainly due to galaxy-galaxy me… ▽ More Massive black holes often exist within dwarf galaxies, and both simulations and observations have shown that a substantial fraction of these may be off-center with respect to their hosts. We trace the evolution of off-center massive black holes (MBHs) in dwarf galaxies using cosmological hydrodynamical simulations, and show that the reason for off-center locations is mainly due to galaxy-galaxy mergers. We calculate dynamical timescales and show that off-center MBHs are unlikely to sink to their galaxys' centers within a Hubble time, due to the shape of the hosts' potential wells and low stellar densities. These wandering MBHs are unlikely to be detected electromagnetically, nor is there a measurable dynamical effect on the galaxy's stellar population. We conclude that off-center MBHs may be common in dwarfs, especially if the mass of the MBH is small or the stellar mass of the host galaxy is large. However detecting them is extremely challenging, because their accretion luminosities are very low and they do not measurably alter the dynamics of their host galaxies. △ Less

Submitted 22 October, 2021; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: updated to reflect MNRAS accepted version. Includes carbon footprint estimate due to supercomputer use

Getting Ready for LISA: The Data, Support and Preparation Needed to Maximize US Participation in Space-Based Gravitational Wave Science

Authors: Kelly Holley-Bockelmann, :, Jillian Bellovary, Peter Bender, Emanuele Berti, Warren Brown, Robert Caldwell, Neil Cornish, Jeremy Darling, Matthew Digman, Mike Eracleous, Kayhan Gultekin, Zoltan Haiman, Kelly Holley-Bockelmann, Joey Key, Shane Larson, Xin Liu, Sean McWilliams, Priyamvada Natarajan, David Shoemaker, Deirdre Shoemaker, Krista Lynne Smith, Marcelle Soares-Santos, Robin Stebbins

Abstract: The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US scienc… ▽ More The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US science center and an open access data model. As other major missions have demonstrated, a science center acts as both a locus and an amplifier of research innovation, data analysis, user support, user training and user interaction. In its most basic function, a US Science Center could facilitate entry into LISA science by hosting a Data Processing Center and a portal for the US community to access LISA data products. However, an enhanced LISA Science Center could: support one of the parallel independent processing pipelines required for data product validation; stimulate the high level of research on data analysis that LISA demands; support users unfamiliar with a novel observatory; facilitate astrophysics and fundamental research; provide an interface into the subtleties of the instrument to validate extraordinary discoveries; train new users; and expand the research community through guest investigator, postdoc and student programs. Establishing a US LISA Science Center well before launch can have a beneficial impact on the participation of the broader astronomical community by providing training, hosting topical workshops, disseminating mock catalogs, software pipelines, and documentation. Past experience indicates that successful science centers are established several years before launch; this early adoption model may be especially relevant for a pioneering mission like LISA. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: 93 pages with a lovely cover page thanks to Bernard Kelly and Elizabeth Ferrara

Aligning Nuclear Cluster Orbits with an Active Galactic Nucleus Accretion Disk

Authors: Gaia Fabj, Syeda S. Nasim, Freddy Caban, K. E. Saavik Ford, Barry McKernan, Jillian M. Bellovary

Abstract: 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 d… ▽ More 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 $ρ\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. △ Less

Submitted 30 October, 2020; v1 submitted 19 June, 2020; originally announced June 2020.

Comments: Submitted to MNRAS; 10 pages, 6 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 499, Issue 2, December 2020, Pages 2608-2616

Orbital Migration of Interacting Stellar Mass Black Holes in Disks around Supermassive Black Holes II. Spins and Incoming Objects

Authors: Amy Secunda, Jillian Bellovary, Mordecai-Mark Mac Low, K. E. Saavik Ford, Barry McKernan, Nathan W. C. Leigh, Wladimir Lyra, Zsolt Sandor, Jose I. Adorno

Abstract: The masses, rates, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and 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. Previou… ▽ More The masses, rates, and spins of merging stellar-mass binary black holes (BBHs) detected by aLIGO and 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. Orbiters entering the inner disk form BBHs with sBHs on resonant orbits near the migration trap. When these sBHs reach ~80 Msun, they form BBHs with sBHs in the migration trap, which over 10 Myr reach ~1000 Msun. We find 68% of the BBHs in our simulation orbit in the retrograde direction, which implies BBHs in our merger channel will have small dimensionless aligned spins, χ_eff. Overall, our models produce BBHs that resemble both the majority of BBH mergers detected thus far (0.66 to 120 Gpc^-3 yr^-1 ) and two recent unusual detections, GW190412 (~0.3 Gpc^-3 yr^-1 ) and GW190521 (~0.1 Gpc^-3 yr^-1 ). △ Less

Submitted 25 September, 2020; v1 submitted 24 April, 2020; originally announced April 2020.

Comments: 16 pages, 9 figures, accepted to ApJ

The role of gravitational recoil in the assembly of massive black hole seeds

Authors: Glenna Dunn, Kelly Holley-Bockelmann, Jillian Bellovary

Abstract: When two black holes merge, the asymmetric emission of gravitational waves provides an impulse to the merged system; this gravitational wave recoil velocity can be up to 4000 km s$^{-1}$, easily fast enough for the black hole to escape its host galaxy. We combine semi-analytic modeling with cosmological zoom-in simulations of a Milky Way-type galaxy to investigate the role of black hole spin and g… ▽ More When two black holes merge, the asymmetric emission of gravitational waves provides an impulse to the merged system; this gravitational wave recoil velocity can be up to 4000 km s$^{-1}$, easily fast enough for the black hole to escape its host galaxy. We combine semi-analytic modeling with cosmological zoom-in simulations of a Milky Way-type galaxy to investigate the role of black hole spin and gravitational recoil in the epoch of massive black hole seeding. We sample four different spin distributions (random, aligned, anti-aligned, and zero spin), and compare the resulting merger rates, occupation fractions, and MBH-host relations with what is expected by excluding the effect of recoil. The inclusion of gravitational recoil and MBH spin in the assembly of MBH seeds can reduce the final $z=5$ MBH mass by up to an order of magnitude. The MBH occupation fraction, however, remains effectively unaltered due to episodes of black hole formation following a recoil event. While electromagnetic detections of these events are unlikely, LISA is ideally suited to detect gravitational wave signals from such events. △ Less

Submitted 11 February, 2020; originally announced February 2020.

Comments: 8 pages, 6 figures. Submitted to ApJ

Building A Field: The Future of Astronomy with Gravitational Waves, A State of The Profession Consideration for Astro2020

Authors: Kelly Holley-Bockelmann, Joey Shapiro Key, Brittany Kamai, Robert Caldwell, Warren Brown, Bill Gabella, Karan Jani, Quentin Baghi, John Baker, Jillian Bellovary, Pete Bender, Emanuele Berti, T. J. Brandt, Curt Cutler, John W. Conklin, Michael Eracleous, Elizabeth C. Ferrara, Bernard J. Kelly, Shane L. Larson, Jeff Livas, Maura McLaughlin, Sean T. McWilliams, Guido Mueller, Priyamvada Natarajan, Norman Rioux , et al. (6 additional authors not shown)

Abstract: Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scienti… ▽ More Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field. △ Less

Submitted 16 December, 2019; originally announced December 2019.

Black Hole Growth and Feedback in Isolated Romulus25 Dwarf Galaxies

Authors: Ray Sharma, Alyson Brooks, Rachel S. Somerville, Michael Tremmel, Jillian Bellovary, Anna Wright, Thomas Quinn

Abstract: We investigate the effects of massive black hole growth on the structural evolution of dwarf galaxies within the Romulus25 cosmological hydrodynamical simulation. We study a sample of 228 central, isolated dwarf galaxies with stellar masses $M_{star} < 10^{10} M_\odot$ and a central BH. We find that the local $M_{BH} - M_{star}$ relation exhibits a high degree of scatter below… ▽ More We investigate the effects of massive black hole growth on the structural evolution of dwarf galaxies within the Romulus25 cosmological hydrodynamical simulation. We study a sample of 228 central, isolated dwarf galaxies with stellar masses $M_{star} < 10^{10} M_\odot$ and a central BH. We find that the local $M_{BH} - M_{star}$ relation exhibits a high degree of scatter below $M_{star} < 10^{10} M_\odot$, which we use to classify BHs as overmassive or undermassive relative to their host $M_{star}$. Overmassive BHs grow through a mixture of BH mergers and relatively high average accretion rates, while undermassive BHs grow slowly through accretion. We find that isolated dwarf galaxies that host overmassive BHs also follow different evolutionary tracks relative to their undermassive BH counterparts, building up their stars and dark matter earlier and experiencing star formation suppression starting around $z=2$. By $z=0.05$, overmassive BH hosts above $M_{star} > 10^{9} M_\odot$ are more likely to exhibit lower central stellar mass density, lower HI gas content, and lower star formation rates than their undermassive BH counterparts. Our results suggest that overmassive BHs in isolated galaxies above $M_{star} > 10^{9} M_\odot$ are capable of driving feedback, in many cases suppressing and even quenching star formation by late times. △ Less

Submitted 27 September, 2021; v1 submitted 13 December, 2019; originally announced December 2019.

Comments: 19 pages, 12 figures

Journal ref: ApJ 897 103 (2020)

The Social Sciences Interdisciplinarity for Astronomy and Astrophysics -- Lessons from the History of NASA and Related Fields

Authors: Anamaria Berea, Kathryn Denning, Monica Vidaurri, Kimberly Arcand, Michael P. Oman-Reagan, Jillian Bellovary, Arsev Umur Aydinoglu, Mark Lupisella

Abstract: In this paper we showcase the importance of understanding and measuring interdisciplinarity and other -disciplinarity concepts for all scientists, the role social sciences have historically played in NASA research and missions, the sparsity of social science interdisciplinarity in space and planetary sciences, including astronomy and astrophysics, while there is an imperative necessity for it, and… ▽ More In this paper we showcase the importance of understanding and measuring interdisciplinarity and other -disciplinarity concepts for all scientists, the role social sciences have historically played in NASA research and missions, the sparsity of social science interdisciplinarity in space and planetary sciences, including astronomy and astrophysics, while there is an imperative necessity for it, and the example of interdisciplinarity between social sciences and astrobiology. Ultimately we give voice to the scientists across all fields with respect to their needs, aspirations and experiences in their interdisciplinary work with social sciences through an ad-hoc survey we conducted within the Astro2020 Decadal Survey scientific community. △ Less

Submitted 17 July, 2019; originally announced July 2019.

The Laser Interferometer Space Antenna: Unveiling the Millihertz Gravitational Wave Sky

Authors: John Baker, Jillian Bellovary, Peter L. Bender, Emanuele Berti, Robert Caldwell, Jordan Camp, John W. Conklin, Neil Cornish, Curt Cutler, Ryan DeRosa, Michael Eracleous, Elizabeth C. Ferrara, Samuel Francis, Martin Hewitson, Kelly Holley-Bockelmann, Ann Hornschemeier, Craig Hogan, Brittany Kamai, Bernard J. Kelly, Joey Shapiro Key, Shane L. Larson, Jeff Livas, Sridhar Manthripragada, Kirk McKenzie, Sean T. McWilliams , et al. (17 additional authors not shown)

Abstract: The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the… ▽ More The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps. △ Less

Submitted 26 July, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: White Paper submitted to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics). v2: fixed a reference

The evolution of kicked stellar-mass black holes in star cluster environments II. Rotating star clusters

Authors: Jeremy J. Webb, Nathan W. C. Leigh, Roberto Serrano, Jillian Bellovary, K. E. Saavik Ford, Barry McKernan, Mario Spera, Alessandro A. Trani

Abstract: 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 fra… ▽ More 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. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 13 pages, 11 figures, Accepted for publication in MNRAS

Disentangling nature from nurture: tracing the origin of seed black holes

Authors: Priyamvada Natarajan, Angelo Ricarte, Vivienne Baldassare, Jillian Bellovary, Peter Bender, Emanuele Berti, Nico Cappelluti, Andrea Ferrara, Jenny Greene, Zoltan Haiman, Kelly Holley-Bockelmann, Guido Mueller, Fabio Pacucci, David Shoemaker, Deirdre Shoemaker, Michael Tremmel, C. Meg Urry, Alexey Vikhlinin, Marta Volonteri

Abstract: The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains chal… ▽ More The origin and properties of black hole seeds that grow to produce the detected population of supermassive black holes are unconstrained at present. Despite the existence of several potentially feasible channels for the production of initial seeds in the high redshift universe, since even actively growing seeds are not directly observable at these epochs, discriminating between models remains challenging. Several new observables that encapsulate information about seeding have been proposed in recent years, and these offer exciting prospects for truly unraveling the nature of black hole seeds in the coming years. One of the key challenges for this task lies in the complexity of the problem, the required disentangling of the confounding effects of accretion physics and mergers, as mergers and accretion events over cosmic time stand to erase these initial conditions. Nevertheless, some unique signatures of seeding do survive and still exist in: local scaling relations between black holes and their galaxy hosts at low-masses; in high-redshift luminosity functions of accreting black holes; and in the total number and mass functions of gravitational wave coalescence events from merging binary black holes. One of the clearest discriminants for seed models are these high redshift gravitational wave detections of mergers from space detectable in the milliHertz range. These predicted event rates offer the most direct constraints on the properties of initial black hole seeds. Improving our theoretical understanding of black hole dynamics and accretion will also be pivotal in constraining seeding models in combination with the wide range of multi-messenger data. △ Less

Submitted 19 April, 2019; originally announced April 2019.

Comments: White paper submitted to the Astro2020 decadal survey

AGN (and other) astrophysics with Gravitational Wave Events

Authors: K. E. Saavik Ford, Imre Bartos, Barry McKernan, Zoltan Haiman, Alessandra Corsi, Azadeh Keivani, Szabolcs Marka, Rosalba Perna, Matthew Graham, Nicholas P. Ross, Daniel Stern, Jillian Bellovary, Emanuele Berti, Matthew O'Dowd, Wladimir Lyra, Mordecai-Mark Mac Low, Zsuzsanna Marka

Abstract: 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, aspec… ▽ More 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.} △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: 5 pages, 2 figures, Astro2020 Decadal Survey white paper

Astro2020 Science White Paper: The Local Relics of of Supermassive Black Hole Seeds

Authors: Jenny E Greene, Aaron Barth, Andrea Bellini, Jillian Bellovary, Kelly Holley-Bockelmann, Tuan Do, Elena Gallo, Karl Gebhardt, Kayhan Gultekin, Zoltan Haiman, Matthew Hosek Jr., Dongwon Kim, Mattia Libralato, Jessica Lu, Kristina Nyland, Matthew Malkan, Amy Reines, Anil Seth, Tommaso Treu, Jonelle Walsh, Joan Wrobel

Abstract: We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observe… ▽ More We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Science White Paper Submitted for the Astro2020 Decadal Survey on Astronomy and Astrophysics

Electromagnetic Window into the Dawn of Black Holes

Authors: Zoltán Haiman, William N. Brandt, Alexey Vikhlinin, Jillian Bellovary, Elena Gallo, Jenny Greene, Kohei Inayoshi, Joseph Lazio, Bret Lehmer, Bin Luo, Piero Madau, Priya Natarajan, Feryal Özel, Fabio Pacucci, Alberto Sesana, Daniel Stern, Christian Vignali, Eli Visbal, Fabio Vito, Marta Volonteri, Joan Wrobel, Emanuele Berti, Volker Bromm, Greg Bryan, Nico Cappelluti , et al. (6 additional authors not shown)

Abstract: Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanyi… ▽ More Massive 10^6-10^10 Msun black holes (BHs) are ubiquitous in local galactic nuclei. They were common by the time the Universe is several Gyr old, and many of them were in place within the first 1~Gyr after the Big Bang. Their quick assembly has been attributed to mechanisms such as the rapid collapse of gas into the nuclei of early protogalaxies, accretion and mergers of stellar-mass BHs accompanying structure formation at early times, and the runaway collapse of early, ultra-dense stellar clusters. The origin of the early massive BHs remains an intriguing and long-standing unsolved puzzle in astrophysics. Here we discuss strategies for discerning between BH seeding models using electromagnetic observations. We argue that the most direct answers will be obtained through detection of BHs with masses M<10^5 Msun at redshifts z>10, where we expect them to first form. Reaching out to these redshifts and down to these masses is crucial, because BHs are expected to lose the memory of their initial assembly by the time they grow well above 10^5 Msun and are incorporated into higher-mass galaxies. The best way to detect 10^4-10^5 Msun BHs at high redshifts is by a sensitive X-ray survey. Critical constraining power is augmented by establishing the properties and the environments of their host galaxies in deep optical/IR imaging surveys. Required OIR data can be obtained with the JWST and WFIRST missions. The required X-ray flux limits (down to 10^{-19} erg/s/cm^2) are accessible only with a next-generation X-ray observatory which has both high (sub-1") angular resolution and high throughput. A combination of deep X-ray and OIR surveys will be capable of probing several generic markers of the BH seed scenarios, and resolving the long-stanging puzzle of their origin. These electromagnetic observations are also highly synergistic with the information from LISA on high-z BH mergers. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Science white paper submitted to the Astro2020 Decadal Survey

Where are the Intermediate Mass Black Holes?

Authors: Jillian Bellovary, Alyson Brooks, Monica Colpi, Michael Eracleous, Kelly Holley-Bockelmann, Ann Hornschemeier, Lucio Mayer, Priya Natarajan, Jacob Slutsky, Michael Tremmel

Abstract: Observational evidence has been mounting for the existence of intermediate mass black holes (IMBHs, 10^2-10^5 Msun), but observing them at all, much less constraining their masses, is very challenging. In one theorized formation channel, IMBHs are the seeds for supermassive black holes in the early universe. As a result, IMBHs are predicted to exist in the local universe in dwarf galaxies, as well… ▽ More Observational evidence has been mounting for the existence of intermediate mass black holes (IMBHs, 10^2-10^5 Msun), but observing them at all, much less constraining their masses, is very challenging. In one theorized formation channel, IMBHs are the seeds for supermassive black holes in the early universe. As a result, IMBHs are predicted to exist in the local universe in dwarf galaxies, as well as wandering in more massive galaxy halos. However, these environments are not conducive to the accretion events or dynamical signatures that allow us to detect IMBHs. The Laser Interferometer Space Antenna (LISA) will demystify IMBHs by detecting the mergers of these objects out to extremely high redshifts, while measuring their masses with extremely high precision. These observations of merging IMBHs will allow us to constrain the formation mechanism and subsequent evolution of massive black holes, from the 'dark ages' to the present day, and reveal the role that IMBHs play in hierarchical galaxy evolution. △ Less

Submitted 22 March, 2019; v1 submitted 19 March, 2019; originally announced March 2019.

Comments: Science white paper #196 submitted to the Astro2020 Decadal Survey

Astro2020 science white paper: The gravitational wave view of massive black holes

Authors: Monica Colpi, Kelly Holley-Bockelmann, Tamara Bogdanovic, Priya Natarajan, Jillian Bellovary, Alberto Sesana, Michael Tremmel, Jeremy Schnittman, Julia Comerford, Enrico Barausse, Emanuele Berti, Marta Volonteri, Fazeel Khan, Sean McWilliams, Sarah Burke-Spolaor, Jeff Hazboun, John Conklin, Guido Mueller, Shane Larson

Abstract: Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic… ▽ More Coalescing, massive black-hole (MBH) binaries are the most powerful sources of gravitational waves (GWs) in the Universe, which makes MBH science a prime focus for ongoing and upcoming GW observatories. The Laser Interferometer Space Antenna (LISA) -- a gigameter scale space-based GW observatory -- will grant us access to an immense cosmological volume, revealing MBHs merging when the first cosmic structures assembled in the Dark Ages. LISA will unveil the yet unknown origin of the first quasars, and detect the teeming population of MBHs of $10^4 - 10^7$ solar masses. forming within protogalactic halos. The Pulsar Timing Array, a galactic-scale GW survey, can access the largest MBHs the Universe, detecting the cosmic GW foreground from inspiraling MBH binaries of about 10^9 solar masses. LISA can measure MBH spins and masses with precision far exceeding that from electromagnetic (EM) probes, and together, both GW observatories will provide the first full census of binary MBHs, and their orbital dynamics, across cosmic time. Detecting the loud gravitational signal of these MBH binaries will also trigger alerts for EM counterpart searches, from decades (PTAs) to hours (LISA) prior to the final merger. By witnessing both the GW and EM signals of MBH mergers, precious information will be gathered about the rich and complex environment in the aftermath of a galaxy collision. The unique GW characterization of MBHs will shed light on the deep link between MBHs of $10^4-10^{10}$ solar masses and the grand design of galaxy assembly, as well as on the complex dynamics that drive MBHs to coalescence. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: 5 pages, 2 figures submitted to the Astro2020 Decadal

Opportunities for Multimessenger Astronomy in the 2020s

Authors: E. Burns, A. Tohuvavohu, J. M. Bellovary, E. Blaufuss, T. J. Brandt, S. Buson, R. Caputo, S. B. Cenko, N. Christensen, J. W. Conklin, F. D'Ammando, K. E. S. Ford, A. Franckowiak, C. Fryer, C. M. Hui, K. Holley-Bockelmann, T. Jaffe, T. Kupfer, M. Karovska, B. D. Metzger, J. Racusin, B. Rani, M. Santander, J. Tomsick, C. Wilson-Hodge

Abstract: 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… ▽ More 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. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Astro2020 White Paper for the 8th Thematic Area of Multimessenger Astronomy and Astrophysics

Orbital Migration of Interacting Stellar Mass Black Holes in Disks around Supermassive Black Holes

Authors: Amy Secunda, Jillian Bellovary, Mordecai-Mark Mac Low, K. E. Saavik Ford, Barry McKernan, Nathan Leigh, Wladimir Lyra, Zsolt Sandor

Abstract: 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. T… ▽ More 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. △ Less

Submitted 12 November, 2019; v1 submitted 8 July, 2018; originally announced July 2018.

Comments: 18 pages, 13 figures, Accepted to ApJ

Multimessenger Signatures of Massive Black Holes in Dwarf Galaxies

Authors: Jillian Bellovary, Colleen Cleary, Ferah Munshi, Michael Tremmel, Charlotte Christensen, Alyson Brooks, Thomas Quinn

Abstract: 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… ▽ More 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. △ Less

Submitted 19 October, 2018; v1 submitted 1 June, 2018; originally announced June 2018.

Comments: replaced with accepted version

Sowing black hole seeds: Direct collapse black hole formation with realistic Lyman-Werner radiation in cosmological simulations

Authors: Glenna Dunn, Jillian Bellovary, Kelly Holley-Bockelmann, Charlotte Christensen, Thomas Quinn

Abstract: We study the birth of supermassive black holes from the direct collapse process and characterize the sites where these black hole seeds form. In the pre-reionization epoch, molecular hydrogen (H$_2$) is an efficient coolant, causing gas to fragment and form Population III stars, but Lyman-Werner radiation can suppress H$_2$ formation and allow gas to collapse directly into a massive black hole. Th… ▽ More We study the birth of supermassive black holes from the direct collapse process and characterize the sites where these black hole seeds form. In the pre-reionization epoch, molecular hydrogen (H$_2$) is an efficient coolant, causing gas to fragment and form Population III stars, but Lyman-Werner radiation can suppress H$_2$ formation and allow gas to collapse directly into a massive black hole. The critical flux required to inhibit H$_2$ formation, $J_{\rm crit}$, is hotly debated, largely due to the uncertainties in the source radiation spectrum, H$_2$ self-shielding, and collisional dissociation rates. Here, we test the power of the direct collapse model in a self-consistent, time-dependant, non-uniform Lyman-Werner radiation field -- the first time such has been done in a cosmological volume -- using an updated version of the SPH+N-body tree code Gasoline with H$_2$ non-equilibrium abundance tracking, H$_2$ cooling, and a modern SPH implementation. We vary $J_{\rm crit} $ from $30$ to $10^3$ in units of $J_{21}$ to study how this parameter impacts the number of seed black holes and the type of galaxies which host them. We focus on black hole formation as a function of environment, halo mass, metallicity, and proximity of the Lyman-Werner source. Massive black hole seeds form more abundantly with lower $J_{\rm crit}$ thresholds, but regardless of $J_{\rm crit}$, these seeds typically form in halos that have recently begun star formation. Our results do not confirm the proposed atomic cooling halo pair scenario; rather black hole seeds predominantly form in low-metallicity pockets of halos which already host star formation. △ Less

Submitted 11 June, 2018; v1 submitted 2 March, 2018; originally announced March 2018.

Comments: 10 pages, 8 figures, accepted by ApJ

On the rate of black hole binary mergers in galactic nuclei due to dynamical hardening

Authors: Nathan W. C. Leigh, Aaron M. Geller, B. McKernan, K. E. S. Ford, M. -M. Mac Low, J. Bellovary, Z. Haiman, W. Lyra, J. Samsing, M. O'Dowd, B. Kocsis, S. Endlich

Abstract: 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 ha… ▽ More 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. △ Less

Submitted 28 November, 2017; originally announced November 2017.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

The Evolution of Kicked Stellar-Mass Black Holes in Star Cluster Environments

Authors: Jeremy J. Webb, Nathan W. C. Leigh, Abhishek Singh, K. E. Saavik Ford, Barry McKernan, Jillian Bellovary

Abstract: 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 pre… ▽ More 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 $σ$. 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}}σ$ in order to account for the variations in the local gravitational potential. △ Less

Submitted 24 November, 2017; originally announced November 2017.

Comments: 13 pages, 9 figures, Accepted for publication in MNRAS

The Contribution Of Outer HI Disks To The Merging Binary Black Hole Population

Authors: Sukanya Chakrabarti, Philip Chang, Richard O'Shaughnessy, Alyson Brooks, Sijing Shen, Jillian Bellovary, Wojciech Gladysz, Chris Belczynski

Abstract: We investigate the contribution of outer HI disks to the observable population of merging black hole binaries. Like dwarf galaxies, the outer HI disks of spirals have low star formation rates and lower metallicities than the inner disks of spirals. Since low-metallicity star formation can produce more detectable compact binaries than typical star formation, the environments in the outskirts of spi… ▽ More We investigate the contribution of outer HI disks to the observable population of merging black hole binaries. Like dwarf galaxies, the outer HI disks of spirals have low star formation rates and lower metallicities than the inner disks of spirals. Since low-metallicity star formation can produce more detectable compact binaries than typical star formation, the environments in the outskirts of spiral galaxies may be conducive to producing a rich population of massive binary black holes. We consider here both detailed controlled simulations of spirals and cosmological simulations, as well as the current range of observed values for metallicity and star formation in outer disks. We find that outer HI disks contribute at least as much as dwarf galaxies do to the observed LIGO/Virgo detection rates. Identifying the host galaxies of merging massive black holes should provide constraints on cosmological parameters and insights into the formation channels of binary mergers. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Comments: accepted to ApJL, 5 pages, 2 figures

The birth of a supermassive black hole binary

Authors: Hugo Pfister, Alessandro Lupi, Pedro R. Capelo, Marta Volonteri, Jillian M. Bellovary, Massimo Dotti

Abstract: We study the dynamical evolution of supermassive black holes, in the late stage of galaxy mergers, from kpc to pc scales. In particular, we capture the formation of the binary, a necessary step before the final coalescence, and trace back the main processes causing the decay of the orbit. We use hydrodynamical simulations of galaxy mergers with different resolutions, from $20\,\rm pc$ down to… ▽ More We study the dynamical evolution of supermassive black holes, in the late stage of galaxy mergers, from kpc to pc scales. In particular, we capture the formation of the binary, a necessary step before the final coalescence, and trace back the main processes causing the decay of the orbit. We use hydrodynamical simulations of galaxy mergers with different resolutions, from $20\,\rm pc$ down to $1\,\rm pc$, in order to study the effects of the resolution on our results, remove numerical effects, and assess that resolving the influence radius of the orbiting black hole is a minimum condition to fully capture the formation of the binary. Our simulations include the relevant physical processes, namely star formation, supernova feedback, accretion onto the black holes and the ensuing feedback. We find that, in these mergers, dynamical friction from the smooth stellar component of the nucleus is the main process that drives black holes from kpc to pc scales. Gas does not play a crucial role and even clumps do not induce scattering or perturb the orbits. We compare the time needed for the formation of the binary to analytical predictions and suggest how to apply such analytical formalism to obtain estimates of binary formation times in lower resolution simulations. △ Less

Submitted 13 June, 2017; originally announced June 2017.

Comments: 12 pages, 12 Figures, submitted to MNRAS

On stellar-mass black hole mergers in AGN disks detectable with LIGO

Authors: B. McKernan, K. E. S. Ford, J. Bellovary, N. W. C. Leigh, Z. Haiman, B. Kocsis, W. Lyra, M. -M. MacLow, B. Metzger, M. O'Dowd, S. Endlich, D. J. Rosen

Abstract: 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 t… ▽ More 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. △ Less

Submitted 30 January, 2018; v1 submitted 24 February, 2017; originally announced February 2017.

Comments: 9pages, 2 figures, revised version submitted to ApJ. Comments welcome

Journal ref: 2018ApJ...866...66M

A survey of dual active galactic nuclei in simulations of galaxy mergers: frequency and properties

Authors: Pedro R. Capelo, Massimo Dotti, Marta Volonteri, Lucio Mayer, Jillian M. Bellovary, Sijing Shen

Abstract: We investigate the simultaneous triggering of active galactic nuclei (AGN) in merging galaxies, using a large suite of high-resolution hydrodynamical simulations. We compute dual-AGN observability time-scales using bolometric, X-ray, and Eddington-ratio thresholds, confirming that dual activity from supermassive black holes (BHs) is generally higher at late pericentric passages, before a merger re… ▽ More We investigate the simultaneous triggering of active galactic nuclei (AGN) in merging galaxies, using a large suite of high-resolution hydrodynamical simulations. We compute dual-AGN observability time-scales using bolometric, X-ray, and Eddington-ratio thresholds, confirming that dual activity from supermassive black holes (BHs) is generally higher at late pericentric passages, before a merger remnant has formed, especially at high luminosities. For typical minor and major mergers, dual activity lasts ~20-70 and ~100-160 Myr, respectively. We also explore the effects of X-ray obscuration from gas, finding that the dual-AGN time decreases at most by a factor of ~2, and of contamination from star formation. Using projected separations and velocity differences rather than three-dimensional quantities can decrease the dual-AGN time-scales by up to ~4, and we apply filters which mimic current observational-resolution limitations. In agreement with observations, we find that, for a sample of major and minor mergers hosting at least one AGN, the fraction harbouring dual AGN is ~20-30 and ~1-10 per cent, respectively. We quantify the effects of merger mass ratio (0.1 to 1), geometry (coplanar, prograde, retrograde, and inclined), disc gas fraction, and BH properties, finding that the mass ratio is the most important factor, with the difference between minor and major mergers varying between factors of a few to orders of magnitude, depending on the luminosity and filter used. We also find that a shallow imaging survey will require very high angular resolution, whereas a deep imaging survey will be less resolution-dependent. △ Less

Submitted 29 April, 2017; v1 submitted 28 November, 2016; originally announced November 2016.

Comments: Accepted for publication in MNRAS; added online-only material

The Pairing of Accreting Massive Black Holes in Multiphase Circumnuclear Disks: the Interplay between Radiative Cooling, Star Formation, and Feedback Processes

Authors: Rafael Souza Lima, Lucio Mayer, Pedro R. Capelo, Jillian M. Bellovary

Abstract: We study the orbital decay of a pair of massive black holes (BHs) with masses 5 * 10^5 and 10^7 M_sun, using hydrodynamical simulations of circumnuclear disks (CNDs) with the alternating presence of sub-grid physics such as radiative cooling, star formation, supernova feedback, BH accretion and feedback. In the absence of such processes, the orbit of the secondary BH decays over timescales of ~10… ▽ More We study the orbital decay of a pair of massive black holes (BHs) with masses 5 * 10^5 and 10^7 M_sun, using hydrodynamical simulations of circumnuclear disks (CNDs) with the alternating presence of sub-grid physics such as radiative cooling, star formation, supernova feedback, BH accretion and feedback. In the absence of such processes, the orbit of the secondary BH decays over timescales of ~10 Myr to the center of the CND, where the primary BH resides. When strong dissipation operates in CNDs, fragmentation into massive objects the size of giant molecular clouds and with densities in the range 10^4 - 10^7 amu / cm^3 occurs, causing stochastic torques and hits that can eject the secondary BH from the midplane. Outside the plane, the low-density medium provides only weak drag, and the BH return is governed by inefficient dynamical friction. In rare cases, clump-BH interactions can lead to a faster decay. Feedback processes lead to outflows, but do not change significantly the overall density of the CND midplane. However, with a spherically distributed BH feedback a hot bubble is generated behind the secondary, which almost shuts off dynamical friction, a phenomenon we dub "wake evacuation", leading to delays in the decay of possibly ~0.3 Gyr. We discuss the non-trivial implications on the discovery space of the eLISA telescope. Our results suggest the largest uncertainty in predicting BH merger rates lies in the potentially wide variety of galaxy host systems, with different degrees of gas dissipation and heating, yielding decay timescales from ~10 to ~300 Myr. △ Less

Submitted 16 February, 2017; v1 submitted 5 October, 2016; originally announced October 2016.

Comments: 20. pages, 14 figures. Accepted for publication in ApJ. Comments are welcome

Cosmological Hydrodynamic Simulations of Preferential Accretion in the SMBH of Milky Way Size Galaxies

Authors: N. Nicole Sanchez, Jillian M. Bellovary, Kelly Holley-Bockelmann, Michael Tremmel, Alyson Brooks, Fabio Governato, Tom Quinn, Marta Volonteri, James Wadsley

Abstract: Using a new, high-resolution cosmological hydrodynamic simulation of a Milky Way-type (MW-type) galaxy, we explore how a merger-rich assembly history affects the mass budget of the central supermassive black hole (SMBH). We examine a MW-mass halo at the present epoch whose evolution is characterized by several major mergers to isolate the importance of merger history on black hole accretion. This… ▽ More Using a new, high-resolution cosmological hydrodynamic simulation of a Milky Way-type (MW-type) galaxy, we explore how a merger-rich assembly history affects the mass budget of the central supermassive black hole (SMBH). We examine a MW-mass halo at the present epoch whose evolution is characterized by several major mergers to isolate the importance of merger history on black hole accretion. This study is an extension of Bellovary et. al. 2013, which analyzed the accretion of high mass, high redshift galaxies and their central black holes, and found that the gas content of the central black hole reflects what is accreted by the host galaxy halo. In this study, we find that a merger-rich galaxy will have a central SMBH preferentially fed by merger gas. Moreover, we find that nearly 30$\%$ of the accreted mass budget of the SMBH enters the galaxy through the two major mergers in its history, which may account for the increase of merger-gas fueling the SMBH. Through an investigation of the angular momentum of the gas entering the host and its SMBH, we determine that merger gas enters the galaxy with lower angular momentum compared to smooth accretion, partially accounting for the preferential fueling witnessed in the SMBH. In addition, the presence of mergers, particularly major mergers, also helps funnel low angular momentum gas more readily to the center of the galaxy. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories. Our results imply that galaxy mergers play an important role in feeding the SMBH in MW-type galaxies with merger-rich histories. △ Less

Submitted 5 June, 2018; v1 submitted 4 October, 2016; originally announced October 2016.

Comments: 7 pages, 5 figures, Accepted by ApJ

The effects of host galaxy properties on merging compact binaries detectable by LIGO

Authors: Richard O'Shaughnessy, Jillian Bellovary, Alyson Brooks, Sijing Shen, Fabio Governato, Charlotte Christensen

Abstract: Cosmological simulations of galaxy formation can produce present-day galaxies with a large range of assembly and star formation histories. A detailed study of the metallicity evolution and star formation history of such simulations can assist in predicting LIGO-detectable compact object binary mergers. Recent simulations of compact binary evolution suggest the compact object merger rate depends se… ▽ More Cosmological simulations of galaxy formation can produce present-day galaxies with a large range of assembly and star formation histories. A detailed study of the metallicity evolution and star formation history of such simulations can assist in predicting LIGO-detectable compact object binary mergers. Recent simulations of compact binary evolution suggest the compact object merger rate depends sensitively on the progenitor's metallicity. Rare low-metallicity star formation during galaxy assembly can produce more detected compact binaries than typical star formation. Using detailed simulations of galaxy and chemical evolution, we determine how sensitively the compact binary populations of galaxies with similar present-day appearance depend on the details of their assembly. We also demonstrate by concrete example the extent to which dwarf galaxies overabundantly produce compact binary mergers, particularly binary black holes, relative to more massive galaxies. We discuss the implications for transient multimessenger astronomy with compact binary sources. △ Less

Submitted 21 September, 2016; originally announced September 2016.

Comments: 9 pages, re-submitted to MNRAS after responding to referee comments

The Romulus Cosmological Simulations: A Physical Approach to the Formation, Dynamics and Accretion Models of SMBHs

Authors: Michael Tremmel, Michael Karcher, Fabio Governato, Marta Volonteri, Tom Quinn, Andrew Pontzen, Lauren Anderson, Jillian Bellovary

Abstract: 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 de- tailed analysis of SMBH-galaxy co-evolution throughout cosmic time. Our scheme includes novel, physically motivated models for SMBH formation, d… ▽ More 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 de- tailed 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. △ Less

Submitted 27 June, 2017; v1 submitted 7 July, 2016; originally announced July 2016.

Comments: 21 pages, 17 figures, Accepted to MNRAS, in press. Updated references

Bar-driven evolution and quenching of spiral galaxies in cosmological simulations

Authors: Daniele Spinoso, Silvia Bonoli, Massimo Dotti, Lucio Mayer, Piero Madau, Jillian Bellovary

Abstract: We analyse the output of the hi-res cosmological zoom-in simulation ErisBH to study self-consistently the formation of a strong stellar bar in a Milky Way-type galaxy and its effect on the galactic structure, on the central gas distribution and on star formation. The simulation includes radiative cooling, star formation, SN feedback and a central massive black hole which is undergoing gas accretio… ▽ More We analyse the output of the hi-res cosmological zoom-in simulation ErisBH to study self-consistently the formation of a strong stellar bar in a Milky Way-type galaxy and its effect on the galactic structure, on the central gas distribution and on star formation. The simulation includes radiative cooling, star formation, SN feedback and a central massive black hole which is undergoing gas accretion and is heating the surroundings via thermal AGN feedback. A large central region in the ErisBH disk becomes bar-unstable after z~1.4, but a clear bar-like structure starts to grow significantly only after z~0.4, possibly triggered by the interaction with a massive satellite. At z~0.1 the bar reaches its maximum radial extent of l~2.2 kpc. As the bar grows, it becomes prone to buckling instability, which we quantify based on the anisotropy of the stellar velocity dispersion. The actual buckling event is observable at z~0.1, resulting in the formation of a boxy-peanut bulge clearly discernible in the edge-on view of the galaxy at z=0. The bar in ErisBH does not dissolve during the formation of the bulge but remains strongly non-axisymmetric down to the resolution limit of ~100 pc at z=0. During its early growth, the bar exerts a strong torque on the gas within its extent and drives gas inflows that enhance the nuclear star formation on sub-kpc scales. Later on the infalling gas is nearly all consumed into stars and, to a lesser extent, accreted onto the central black hole, leaving behind a gas-depleted region within the central ~2 kpc. Observations would more likely identify a prominent, large-scale bar at the stage when the galactic central region has already been quenched. Bar-driven quenching may play an important role in disk-dominated galaxies at all redshift. [Abridged] △ Less

Submitted 7 July, 2016; originally announced July 2016.

Comments: 13 pages, 12 figures, MNRAS submitted

Supermassive black hole pairs in clumpy galaxies at high redshift: delayed binary formation and concurrent mass growth

Authors: Valentina Tamburello, Pedro R. Capelo, Lucio Mayer, Jillian M. Bellovary, James Wadsley

Abstract: Massive gas-rich galaxy discs at $z \sim 1-3$ host massive star-forming clumps with typical baryonic masses in the range $10^7-10^8$~M$_{\odot}$ which can affect the orbital decay and concurrent growth of supermassive black hole (BH) pairs. Using a set of high-resolution simulations of isolated clumpy galaxies hosting a pair of unequal-mass BHs, we study the interaction between massive clumps and… ▽ More Massive gas-rich galaxy discs at $z \sim 1-3$ host massive star-forming clumps with typical baryonic masses in the range $10^7-10^8$~M$_{\odot}$ which can affect the orbital decay and concurrent growth of supermassive black hole (BH) pairs. Using a set of high-resolution simulations of isolated clumpy galaxies hosting a pair of unequal-mass BHs, we study the interaction between massive clumps and a BH pair at kpc scales, during the early phase of the orbital decay. We find that both the interaction with massive clumps and the heating of the cold gas layer of the disc by BH feedback tend to delay significantly the orbital decay of the secondary, which in many cases is ejected and then hovers for a whole Gyr around a separation of 1--2 kpc. In the envelope, dynamical friction is weak and there is no contribution of disc torques: these lead to the fastest decay once the orbit of the secondary BH has circularised in the disc midplane. In runs with larger eccentricities the delay is stronger, although there are some exceptions. We also show that, even in discs with very sporadic transient clump formation, a strong spiral pattern affects the decay time-scale for BHs on eccentric orbits. We conclude that, contrary to previous belief, a gas-rich background is not necessarily conducive to a fast BH decay and binary formation, which prompts more extensive investigations aimed at calibrating event-rate forecasts for ongoing and future gravitational-wave searches, such as with Pulsar Timing Arrays and the future evolved Laser Interferometer Space Antenna. △ Less

Submitted 5 October, 2016; v1 submitted 29 February, 2016; originally announced March 2016.

Comments: Accepted by MNRAS

Migration Traps in Disks Around Supermassive Black Holes

Authors: Jillian Bellovary, Mordecai-Mark Mac Low, Barry McKernan, K. E. Saavik Ford

Abstract: 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 determ… ▽ More 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. △ Less

Submitted 26 February, 2016; v1 submitted 30 October, 2015; originally announced November 2015.

Comments: replaced with ApJL accepted version

Black Hole Starvation and Bulge Evolution in a Milky Way-like Galaxy

Authors: Silvia Bonoli, Lucio Mayer, Stelios Kazantzidis, Piero Madau, Jillian Bellovary, Fabio Governato

Abstract: We present a new zoom-in hydrodynamical simulation, "Erisbh", which follows the cosmological evolution and feedback effects of a supermassive black hole at the center of a Milky Way-type galaxy. ErisBH shares the same initial conditions, resolution, recipes of gas cooling, star formation and feedback, as the close Milky Way-analog "Eris", but it also includes prescriptions for the formation, growt… ▽ More We present a new zoom-in hydrodynamical simulation, "Erisbh", which follows the cosmological evolution and feedback effects of a supermassive black hole at the center of a Milky Way-type galaxy. ErisBH shares the same initial conditions, resolution, recipes of gas cooling, star formation and feedback, as the close Milky Way-analog "Eris", but it also includes prescriptions for the formation, growth and feedback of supermassive black holes. We find that the galaxy's central black hole grows mainly through mergers with other black holes coming from infalling satellite galaxies. The growth by gas accretion is minimal because very little gas reaches the sub-kiloparsec scales. The final black hole is, at z=0, about 2.6 million solar masses and it sits closely to the position of SgrA* on the MBH-MBulge and MBH-sigma planes, in a location consistent with what observed for pseudobulges. Given the limited growth due to gas accretion, we argue that the mass of the central black hole should be above 10^5 solar masses already at z~8. The effect of AGN feedback on the host galaxy is limited to the very central few hundreds of parsecs. Despite being weak, AGN feedback seems to be responsible for the limited growth of the central bulge with respect to the original Eris, which results in a significantly flatter rotation curve in the inner few kiloparsecs. Moreover, the disk of ErisBH is more prone to instabilities, as its bulge is smaller and its disk larger then Eris. As a result, the disk of ErisBH undergoes a stronger dynamical evolution relative to Eris and around z=0.3 a weak bar grows into a strong bar of a few disk scale lengths in size. The bar triggers a burst of star formation in the inner few hundred parsecs, provides a modest amount of new fuel to the central black hole, and causes the bulge of ErisBH to have, by z=0, a box/peanut morphology.(Abridged) △ Less

Submitted 28 August, 2015; originally announced August 2015.

Comments: 16 pages, 16 figures. Submitted to MNRAS

A comparison of black hole growth in galaxy mergers with Gasoline and Ramses

Authors: J. M. Gabor, Pedro R. Capelo, Marta Volonteri, Frédéric Bournaud, Jillian Bellovary, Fabio Governato, Thomas Quinn

Abstract: Supermassive black hole dynamics during galaxy mergers is crucial in determining the rate of black hole mergers and cosmic black hole growth. As simulations achieve higher resolution, it becomes important to assess whether the black hole dynamics is influenced by the treatment of the interstellar medium in different simulation codes. We here compare simulations of black hole growth in galaxy merge… ▽ More Supermassive black hole dynamics during galaxy mergers is crucial in determining the rate of black hole mergers and cosmic black hole growth. As simulations achieve higher resolution, it becomes important to assess whether the black hole dynamics is influenced by the treatment of the interstellar medium in different simulation codes. We here compare simulations of black hole growth in galaxy mergers with two codes: the Smoothed Particle Hydrodynamics code Gasoline, and the Adaptive Mesh Refinement code Ramses. We seek to identify predictions of these models that are robust despite differences in hydrodynamic methods and implementations of sub-grid physics. We find that the general behavior is consistent between codes. Black hole accretion is minimal while the galaxies are well-separated (and even as they "fly-by" within 10 kpc at first pericenter). At late stages, when the galaxies pass within a few kpc, tidal torques drive nuclear gas inflow that triggers bursts of black hole accretion accompanied by star formation. We also note quantitative discrepancies that are model-dependent: our Ramses simulations show less star formation and black hole growth, and a smoother gas distribution with larger clumps and filaments, than our Gasoline simulations. We attribute these differences primarily to the sub-grid models for black hole fueling and feedback and gas thermodynamics. The main conclusion is that differences exist quantitatively between codes, and this should be kept in mind when making comparisons with observations. However, reassuringly, both codes capture the same dynamical behaviors in terms of triggering of black hole accretion, star formation, and black hole dynamics. △ Less

Submitted 10 August, 2015; originally announced August 2015.

Comments: 11 pages, 7 figures. Submitted to A&A. Comments welcome

Black hole accretion versus star formation rate: theory confronts observations

Authors: Marta Volonteri, Pedro R. Capelo, Hagai Netzer, Jillian Bellovary, Massimo Dotti, Fabio Governato

Abstract: We use a suite of hydrodynamical simulations of galaxy mergers to compare star formation rate (SFR) and black hole accretion rate (BHAR) for galaxies before the interaction ('stochastic' phase), during the `merger' proper, lasting ~0.2-0.3 Gyr, and in the `remnant' phase. We calculate the bi-variate distribution of SFR and BHAR and define the regions in the SFR-BHAR plane that the three phases occ… ▽ More We use a suite of hydrodynamical simulations of galaxy mergers to compare star formation rate (SFR) and black hole accretion rate (BHAR) for galaxies before the interaction ('stochastic' phase), during the `merger' proper, lasting ~0.2-0.3 Gyr, and in the `remnant' phase. We calculate the bi-variate distribution of SFR and BHAR and define the regions in the SFR-BHAR plane that the three phases occupy. No strong correlation between BHAR and galaxy-wide SFR is found. A possible exception are galaxies with the highest SFR and the highest BHAR. We also bin the data in the same way used in several observational studies, by either measuring the mean SFR for AGN in different luminosity bins, or the mean BHAR for galaxies in bins of SFR. We find that the apparent contradiction or SFR versus BHAR for observed samples of AGN and star forming galaxies is actually caused by binning effects. The two types of samples use different projections of the full bi-variate distribution, and the full information would lead to unambiguous interpretation. We also find that a galaxy can be classified as AGN-dominated up to 1.5 Gyr after the merger-driven starburst took place. Our study is consistent with the suggestion that most low-luminosity AGN hosts do not show morphological disturbances. △ Less

Submitted 27 May, 2015; v1 submitted 25 February, 2015; originally announced February 2015.

Comments: MNRAS Letters, in press

Growing black holes and galaxies: black hole accretion versus star formation rate

Authors: Marta Volonteri, Pedro R. Capelo, Hagai Netzer, Jillian Bellovary, Massimo Dotti, Fabio Governato

Abstract: We present a new suite of hydrodynamical simulations and use it to study, in detail, black hole and galaxy properties. The high time, spatial and mass resolution, and realistic orbits and mass ratios, down to 1:6 and 1:10, enable us to meaningfully compare star formation rate (SFR) and BH accretion rate (BHAR) timescales, temporal behaviour and relative magnitude. We find that (i) BHAR and galaxy-… ▽ More We present a new suite of hydrodynamical simulations and use it to study, in detail, black hole and galaxy properties. The high time, spatial and mass resolution, and realistic orbits and mass ratios, down to 1:6 and 1:10, enable us to meaningfully compare star formation rate (SFR) and BH accretion rate (BHAR) timescales, temporal behaviour and relative magnitude. We find that (i) BHAR and galaxy-wide SFR are typically temporally uncorrelated, and have different variability timescales, except during the merger proper, lasting ~0.2-0.3 Gyr. BHAR and nuclear (<100 pc) SFR are better correlated, and their variability are similar. Averaging over time, the merger phase leads typically to an increase by a factor of a few in the BHAR/SFR ratio. (ii) BHAR and nuclear SFR are intrinsically proportional, but the correlation lessens if the long-term SFR is measured. (iii) Galaxies in the remnant phase are the ones most likely to be selected as systems dominated by an active galactic nucleus (AGN), because of the long time spent in this phase. (iv) The timescale over which a given diagnostic probes the SFR has a profound impact on the recovered correlations with BHAR, and on the interpretation of observational data. △ Less

Submitted 23 February, 2015; originally announced February 2015.

Comments: Accepted for publication in MNRAS

Massive Black Hole Science with eLISA

Authors: Enrico Barausse, Jillian Bellovary, Emanuele Berti, Kelly Holley-Bockelmann, Brian Farris, Bangalore Sathyaprakash, Alberto Sesana

Abstract: The evolving Laser Interferometer Space Antenna (eLISA) will revolutionize our understanding of the formation and evolution of massive black holes along cosmic history by probing massive black hole binaries in the $10^3-10^7$ solar mass range out to redshift $z\gtrsim 10$. High signal-to-noise ratio detections of $\sim 10-100$ binary coalescences per year will allow accurate measurements of the pa… ▽ More The evolving Laser Interferometer Space Antenna (eLISA) will revolutionize our understanding of the formation and evolution of massive black holes along cosmic history by probing massive black hole binaries in the $10^3-10^7$ solar mass range out to redshift $z\gtrsim 10$. High signal-to-noise ratio detections of $\sim 10-100$ binary coalescences per year will allow accurate measurements of the parameters of individual binaries (such as their masses, spins and luminosity distance), and a deep understanding of the underlying cosmic massive black hole parent population. This wealth of unprecedented information can lead to breakthroughs in many areas of physics, including astrophysics, cosmology and fundamental physics. We review the current status of the field, recent progress and future challenges. △ Less

Submitted 8 January, 2015; v1 submitted 10 October, 2014; originally announced October 2014.

Comments: 22 pages, 3 figures. Minor change (added two references) to match version accepted in the Proceedings of LISA Symposium X (Journal of Physics: Conference Series)

Growth and activity of black holes in galaxy mergers with varying mass ratios

Authors: Pedro R. Capelo, Marta Volonteri, Massimo Dotti, Jillian M. Bellovary, Lucio Mayer, Fabio Governato

Abstract: 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. (… ▽ More 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 secular-related 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. △ Less

Submitted 27 November, 2014; v1 submitted 29 August, 2014; originally announced September 2014.

Comments: Accepted for publication in MNRAS, 23 pages, 22 figures, 3 tables

Effects of Inclination on Measuring Velocity Dispersion and Implications for Black Holes

Authors: Jillian Bellovary, Kelly Holley-Bockelmann, Kayhan Gültekin, Charlotte Christensen, Fabio Governato, Alyson Brooks, Sarah Loebman, Ferah Munshi

Abstract: The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$σ$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are… ▽ More The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$σ$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are essentially straightforward. We examine five disk galaxies from cosmological SPH simulations and find that line-of-sight effects due to galaxy orientation can affect the measured $σ$ by 30%, and consequently black hole mass predictions by up to 1.0 dex. Face-on orientations correspond to systematically lower velocity dispersion measurements, while more edge-on orientations give higher velocity dispersions, due to contamination by disk stars when measuring line of sight quantities. We caution observers that the uncertainty of velocity dispersion measurements is at least 20 km/s, and can be much larger for moderate inclinations. This effect may account for some of the scatter in the locally measured M-$σ$ relation, particularly at the low-mass end. We provide a method for correcting observed $σ_{\rm los}$ values for inclination effects based on observable quantities. △ Less

Submitted 22 September, 2014; v1 submitted 1 May, 2014; originally announced May 2014.

Comments: 11 pages, 7 figures, replaced with accepted version

Nuclear coups: dynamics of black holes in galaxy mergers

Authors: Sandor Van Wassenhove, Pedro R. Capelo, Marta Volonteri, Massimo Dotti, Jillian M. Bellovary, Lucio Mayer, Fabio Governato

Abstract: We study the dynamical evolution of supermassive black holes (BHs) in merging galaxies on scales of hundreds of kpc to 10 pc, to identify the physical processes that aid or hinder the orbital decay of BHs. We present hydrodynamical simulations of galaxy mergers with a resolution of $\leq$20 pc, chosen to accurately track the motion of the nuclei and provide a realistic environment for the evolutio… ▽ More We study the dynamical evolution of supermassive black holes (BHs) in merging galaxies on scales of hundreds of kpc to 10 pc, to identify the physical processes that aid or hinder the orbital decay of BHs. We present hydrodynamical simulations of galaxy mergers with a resolution of $\leq$20 pc, chosen to accurately track the motion of the nuclei and provide a realistic environment for the evolution of the BHs. We find that, during the late stages of the merger, tidal shocks inject energy in the nuclei, causing one or both nuclei to be disrupted and leaving their BH `naked', without any bound gas or stars. In many cases, the nucleus that is ultimately disrupted is that of the larger galaxy (`nuclear coup'), as star formation grows a denser nuclear cusp in the smaller galaxy. We supplement our simulations with an analytical estimate of the orbital-decay time required for the BHs to form a binary at unresolved scales, due to dynamical friction. We find that, when a nuclear coup occurs, the time-scale is much shorter than when the secondary's nucleus is disrupted, as the infalling BH is more massive, and it also finds itself in a denser stellar environment. △ Less

Submitted 27 January, 2014; v1 submitted 28 October, 2013; originally announced October 2013.

Comments: Accepted for publication in MNRAS, 16 pages, 13 figures, 2 tables