MNRAS 000, ??{?? (2017) Preprint 22 December 2020 Compiled using MNRAS LATEX style file v3.0 A new channel to form IMBHs throughout cosmic time Priyamvada Natarajan1;2;3 1 Department of Astronomy, Yale University, New Haven, CT 06511 2 Department of Physics, Yale University, New Haven, CT 06511 3 Yale Center for Astronomy & Astrophysics, Yale University, New Haven, CT 06520 22 December 2020 ABSTRACT While the formation of the first black holes at high redshift is reasonably well under- stood though debated, massive black hole (BH) formation at later cosmic epochs has not been adequately explored. We present a gas accretion driven mechanism that can build up black hole BH masses rapidly in dense, gas-rich nuclear star clusters (NSCs). Wind-fed supra-exponential accretion in these environments under the assumption of net zero angular momentum for the gas, can lead to extremely rapid growth, scal- ing stellar mass remnant seed BHs up to the intermediate mass black hole (IMBH) range. This new long-lived channel for IMBH formation permits growth to final masses 5 ranging from 50 − 10 M . Growth is modulated by the gas supply, and premature termination can result in the formation of BHs with masses between 50 - a few 100M filling in the so-called mass gap. Typically, growth is unimpeded and will result in the 5 formation of IMBHs with masses ranging from ∼ 100 − 10 M . New detections from 5 the LIGO-VIRGO source GW190521 to the emerging population of ∼ 10 M BHs harbored in low-mass dwarf galaxies - are revealing this elusive population. Naturally accounting for the presence of off-center BHs in low-mass dwarfs, this new pathway also predicts the existence of a population of wandering non-central BHs in more mas- sive galaxies detectable via tidal disruption events and as GW coalescences. Gas rich NSCs could therefore serve as incubators for the continual formation of BHs over a wide range in mass throughout cosmic time. Key words: early universe|black hole physics|accretion|Stars: dynamics and kinematics|Galaxy: formation|galaxies: high-redshift 1 INTRODUCTION end of the IMBH mass range (The LIGO Scientific Collab- oration et al.(2020)). Intriguingly, the mass of the primary Our universe appears to be littered with black holes, near in GW190521 lies in the so-called pair instability supernova and far, active and dormant, with masses ranging from a mass gap, the mass range between 50 - 150 M in which few M to supermassive ones weighing 106−10 M . Bridg- stellar processes are not expected to leave behind any com- ing these mass ranges, an expected population of interme- pact remnants Woosley(2017). The large mass of the pri- diate mass black holes (IMBHs) with masses in the range mary black hole in the GW190521 merger puts it right in 100 − 105 M are predicted to exist (Rees 1984; Bond et al. the mass gap and therefore necessitates alternate formation 1984; Madau & Rees 2001). This population of IMBHs has mechanisms. been challenging to detect, but they are now slowly start- ing to be uncovered. The more massive IMBHs with masses Observational detection of accreting black holes has pri- 5 arXiv:2009.09156v2 [astro-ph.GA] 19 Dec 2020 of ∼ 10 M hosted in low-mass and low-luminosity galax- marily been in electromagnetic wavelengths, while more re- ies have been detected across a range of wavelengths from cently detection via gravitational wave (GW) emission dur- X-ray to radio (Reines & Deller 2012; Baldassare et al. ing BH-BH coalescences by the LIGO-VIRGO collaboration 2015; Lemons et al. 2015; Pardo et al. 2016; Nucita et al. has opened up a new window into the mass function of stellar 2017, 2018; Chilingarian et al. 2018; Baldassare et al. 2018; mass black holes and potentially IMBHs (The LIGO Scien- Reines et al. 2020; Mezcua & Dom´ınguezS´anchez 2020). In tific Collaboration et al. 2016, 2020). Accreting supermas- addition, the recent detection of the LIGO-VIRGO source sive black holes (SMBHs) have been typically detected as GW190521 - a merger of two black holes with masses ∼ 85 quasars and AGN at high and moderate redshifts, while ac- M and 66 M and a final merged mass of ∼ 142 M - creting stellar mass black holes, are typically detected in has revealed the existence of black holes at the low mass the nearby universe as X-ray binaries. In our own Galaxy, © 2017 The Authors 2 Natarajan a total of 106 − 108 stellar mass black holes are predicted et al.(2008) and many other authors subsequently). This, it to exist, however, only a hundred or so X-ray binaries are is argued, happens when a major merger between galaxies seen, therefore it is anticipated that a vast population of dor- funnels gas to the centre of the galaxy, that both triggers mant stellar mass black holes might be lurking (see review by an active galactic nucleus (AGN) and a burst of star forma- Remillard & McClintock(2006)). Obtaining a full census of tion (Di Matteo et al. 2005). One empirical clue supporting non-accreting stellar mass black holes has been challenging, this picture is the observed \AGN main sequence" a linear though the recent detection of unusual triple systems like correlation between SMBH accretion rates and star forma- HR 6819 are likely revealing a new repository of dormant tion rates inferred from fitting spectral energy distributions, stellar mass remnants (Rivinius et al. 2020). The prospects detected by stacking star-forming galaxies (Mullaney et al. for a more complete multi-messenger census of stellar mass 2012). black holes looks extremely promising. It has been challenging observationally to unravel the While for dormant SMBHs indirect detection via their link between galaxy mergers and AGN activity, see for ex- dynamical influence on stars orbiting in their vicinity in ample (Mechtley et al. 2016) and recent review by Alexan- galactic centers is viable for nearby galaxies, we need to der & Hickox(2012). It has been proposed that perhaps it rely on occasional luminous tidal disruption events (TDEs) is only the most luminous AGN that are triggered via ma- - stripping of stars that stray close - for more distant galax- jor mergers (Treister et al. 2012; Hong et al. 2015; Marian ies (Rees 1988). Once again, with advances in time-domain et al. 2019), although this matter is far from settled (Vill- astronomy - more sophisticated survey strategies and dedi- forth et al. 2017). Some authors have suggested that these cated instrumentation - there has been tremendous progress empirical correlations seen locally could simply be a conse- in the detection of transient tidal disruption flares. At the quence of the central limit theorem: a correlation can result present time, several dozen TDEs have been observed with simply from many mergers of initially uncorrelated SMBHs multi-wavelength surveys and analyzed (c.f. recent reviews and host galaxies that occur over cosmic time (Peng 2007; by Komossa(2015); Stone et al.(2017)). Jahnke & Macci`o 2011). The local demography of SMBHs suggests that most Whether this correlation extends with the same ob- galaxies in the universe harbor them at their centers (Kor- servationally determined slope all the way down to IMBH mendy & Richstone 1995). SMBHs are detected in galactic masses or whether it flattens at some characteristic mass centers using a multiplicity of observational techniques us- scale, is yet to be determined (Greene 2012). Till recently, ing gas, stellar, or maser dynamics as tracers (Saglia et al. the issue of whether low surface brightness, low mass galax- 2016; van den Bosch 2016), including Sgr A* at the cen- ies harbor central black holes at all was unsettled. New ob- tre of our own galaxy (Ghez et al. 2008; Genzel et al. 2010; servational studies have been steadily uncovering accreting 5 Gillessen et al. 2016). Studies have also revealed correlations IMBHs ∼ 10 M black holes - faint AGN - in low mass, between BH mass and various host galaxy properties, like dwarf galaxies (Reines & Deller 2012; Baldassare et al. 2015; the stellar mass and luminosity of the bulge, however, the Lemons et al. 2015; Pardo et al. 2016; Nucita et al. 2017, tightest correlation is the so-called M• − σ relation, that re- 2018; Baldassare et al. 2018; Reines et al. 2020). These AGN lates the mass of the central SMBH to the velocity dispersion detections out to z ∼ 3, are hosted in either disky or dwarf 9 of stars in inner-most regions of the host galaxy (Ferrarese galaxies with stellar masses M∗ 6 10 M and have been 2002; Tremaine et al. 2002; Kormendy & Ho 2013). These detected via their narrow emission line diagnostics as well relations, it was surmised, are likely the result of SMBH as broad lines that have enabled their mass estimates; X- and host galaxy coevolution in the context of the currently ray and radio emission signatures as well as the presence accepted LCDM model. While the idea of co-evolution has of higher ionization emission lines and on the basis of opti- increasingly garnered support, it remains unknown how and cal variability (Mezcua et al.(2018); Mezcua & Dom´ınguez precisely when, in terms of cosmic epoch, these relations are S´anchez(2020)). In several instances, as shown in this most established; how far down in black hole mass they extend recent compilation, Reines et al.(2020) find that these ac- to and if they evolve with redshift. The origin of these cor- creting sources lie off-center, displaced from the galactic nu- relations, and whether they belie underlying physical cau- cleus. sation mechanisms, perhaps an imprint of the initial condi- In cosmological simulations tailored to track the growth tions that reflect the formation of initial black hole seeds or history of black holes in dwarf galaxies, Bellovary et al.