Arxiv:2009.11326V2 [Astro-Ph.HE] 25 Dec 2020
Total Page:16
File Type:pdf, Size:1020Kb
Draft version December 29, 2020 Typeset using LATEX twocolumn style in AASTeX63 Growth of stellar mass black holes in dense molecular clouds and GW190521 Jared R. Rice 1 and Bing Zhang 2 1Department of Physics, Texas State University, San Marcos, TX 78666, USA 2Department of Physics and Astronomy, University of Nevada Las Vegas, Las Vegas, NV 89154, USA ABSTRACT A stellar mass black hole can grow its mass noticeably through Bondi accretion, if it is embedded in an extremely dense and massive molecular cloud with slow motion with respect to the ambient medium for an extended period of time. This provides a novel, yet challenging channel for the formation of massive stellar-mass black holes. We discuss how this channel may account for the massive binary black hole merger system GW190521 as observed by LIGO/Virgo gravitational wave detectors as well as the claimed massive black hole candidate LB-1. Keywords: Black holes (162), Bondi accretion (174), Gravitational waves (678) 1. INTRODUCTION cretes gas inside of a nuclear star cluster. Previous Recently, the LIGO/Virgo Scientific Collaboration re- works on single and binary black hole accretion, such ported the discovery Abbott et al.(2020) of a bi- as (Roupas & Kazanas 2019a,b), have suggested that nary black hole (BH) merger gravitational wave event, massive black holes with masses in the pair-instability GW190521. The masses of the two pre-merger BHs are gap or lower may form via accretion inside star clusters. +21 +17 More exotic explanations for the massive black holes ob- 85−14M and 66−18M , respectively, and the mass of +28 served in GW190521 may be found in Sakstein et al. the final BH is 142−16M . The probability that at least one of the merging black holes is in the mass range of (2020), where several possible explanations beyond the the so-called pulsational pair-instability supernova gap Standard Model are discussed. An in-depth account of black hole accretion inside nuclear clusters is discussed (65 120M ) is 99.0%, suggesting that such systems are not− easy to form based on standard stellar evolution in Kroupa et al.(2020). models. Abbott et al.(2020) discussed several possi- Another massive stellar-mass BH candidate is in a bi- ble astrophysical channels for forming such massive bi- nary system called LB-1, which was initially reported +11 nary BH systems, such as hierarchical mergers (in glob- to have a mass M = 68−13 M (Liu et al. 2019). If 1 ular clusters), stellar mergers (direct collapse of merged this mass is true , then its existence also challenges the stars with an oversized hydrogen envelope), mergers in standard stellar evolution theories. AGN disks, as well as primordial BHs. They also disfa- Here we propose an alternative astrophysical channel vored the possibility of strong gravitational lensing for to form massive stellar-mass BHs. We suggest that if GW190521. De Luca et al.(2020) proposed that the molecular clouds (MCs) are formed in the regions where BH masses observed in GW190521 may be explained stellar-mass BHs (from the deaths of the previous gener- ation of massive stars) are located, the BHs may undergo arXiv:2009.11326v2 [astro-ph.HE] 25 Dec 2020 by single or binary primordial black hole systems that have undergone a period of early (pre-reionization) ac- significant growth through accretion from the ambient cretion. Safarzadeh & Haiman(2020) suggested that gas under certain conditions. We derive these conditions the BH is a remnant of a Population III star, accret- and discuss the possible formation of the GW190521 ing ambient high-density gas inside of a high-redshift progenitor system and LB-1 within such a framework. minihalo. Natarajan(2020) suggested that the BH ac- 1 This suggested mass was questioned by some groups (e.g. Abdul- Masih et al. 2020; Shenar et al. 2020; Eldridge et al. 2020). Re- [email protected] cent updated analyses by the original group (Liu et al. 2020a,b) derived a lower BH mass for LB-1. [email protected] 2 Rice and Zhang 2. BLACK HOLES IN MOLECULAR CLOUDS model to explain ultraluminous X-ray sources (ULXs) as Isolated BHs in MCs have been discussed in the accreting intermediate mass BHs. In his model, the BHs past, mostly within the context of radiation signatures accrete at the standard BHL rate, but since radiative that reveal their existence (McDowell 1985; Campana feedback significantly reduces the local number density & Pardi 1993). Campana & Pardi(1993) estimated and increases the ambient ISM temperature, the actual −4 −1 accretion rate is reduced by a few orders of magnitude. that each generation of stars produce 3:5 10 M BHs. MCs have a typical diameter ∼40 pc× and mass Park & Ricotti(2013) performed 2-D axisymmetric hy- 5 6 ∼ drodynamical simulations to arrive at an analytical for- 10 10 M . Considering a filling factor of MCs in the Milky− Way, they estimated there are 2 detectable mulation of BHL accretion, regulated by radiative feed- accreting BHs per MC. In these previous∼ studies, the back, in the case of a moving BH. They found that for growth of BHs was not investigated. a BH supersonically traversing a high density medium, Following our previous study on the growth of pri- the accretion rate onto the BH is steady. Interestingly, if mordial BHs Rice & Zhang(2017), we study the possi- the BH's relative velocity is about a few times the sound ble growth of BHs embedded in MCs within the lifetime speed, the accretion rate maintains or even rises above of MCs. Jeffreson & Kruijssen(2018); Jeffreson et al. the Bondi rate (Bondi 1952) even if radiative feedback (2020) proposed an analytical model of MC lifetimes is considered. This is because the dense ionization front that incorporates galactic dynamics. In their model (see and bow shock decrease the downstream relative gas ve- Jeffreson & Kruijssen(2018) Eq. 24), the lifetime of a locity. Additionally, Park & Ricotti(2013) discovered molecular cloud is the sum that the accretion rate in this parameter regime is in- dependent of the ambient ISM temperature so that the −1 −1 −1 −1 −1 −1 τ = τ + τ + τ + τ τ ; (1) growth rate of moving BHs is actually higher than that κ ΩP ff,g cc − β for stationary BHs. In our problem, the MCs have low where τ is the epicyclic perturbation timescale, τ is κ ΩP sound speeds (see Eq.(6) below), so that the supersonic the spiral arm crossing timescale, τff,g is the ISM freefall condition in the Park & Ricotti(2013) simulations is timescale, τcc is the cloud-cloud collision timescale, and readily satisfied. τβ is the galactic shear timescale. In this analytical In the following, we adopt the following MC fiducial model, the predicted MC lifetimes vary with the galac- 5 −3 parameters: electron number density n = (10 cm ) n5, tocentric radius. For the Milky Way, their prediction is 6 mass Mmc = (10 M ) Mmc;6, so that the radius that the MC lifetimes should vary between 20 Myr 1=3 of the cloud is Rmc = [(3=4π)Mmc=nmp] = and 60 Myr. Beyond the Solar neighborhood∼ at dis- (4:6 pc) M 1=3 n−1=3, where m is the proton mass. In tances∼ greater than 8:0 kpc, the predicted MC lifetime mc;6 5 p general, the accretion rate of an isolated BH can be ex- increases to its maximum∼ value. In our model the BHs pressed as exist prior to the formation of the MCs, so one may use dM 20 60 Myr as the duration for BH mass growth. = minM_ ; M_ ; (2) − dt BHL E 3. BLACK HOLE ACCRETION which is bound by the Eddington accretion rate We consider quasi-spherical accretion of BHs in MCs _ 4πGmp for two reasons. This allows us to estimate the ac- ME = M εσT c cretion rate using the simplest Bondi-Hoyle-Lyttleton −7 −1 −1 M (hereafter BHL; see Bondi(1952)) rate. Also a low an- (2:2 10 M yr ) " ; (3) ' × −1 10 M gular momentum cloud can avoid an extended accretion disk so that the inflow can maintain a low sound speed. where G is Newton's gravitation constant, σT is the Numerical simulations (Waters et al. 2020) show that Thomson scattering cross section for electrons, M is the inhomogeneous distributions of rotational velocity and mass of the BH, c is the speed of light, and " denotes density along the outer radius marking the sphere of in- the efficiency of converting the accretion power to lumi- nosity, which has a typical value of 0:1. The BHL fluence of a BH would modify Bondi-like accretion to ∼ form outflows along with inflows, but with the accretion accretion rate is 2 rate remaining close to the Bondi rate value. Simula- 4πλcG 2 M_ BHL = mpnM tions of choked black hole accretion with the accretion V 3 −3 rate remaining near the Bondi rate value were shown in −8 −1 V (1:5 10 M yr ) Aguayo-Ortiz et al.(2019, 2020); Tejeda et al.(2020). ' × 10 km s−1 One factor that will slow down BH growth is the feed- 2 M back effect. For example, Krolik(2004) proposed a n5 ; (4) × 10 M Growth of stellar mass black holes via Bondi accretion 3 where massive BH can undergo significant growth if v is small 2 2 1=2 5 −3 V (cs + v ) ; (5) enough. For our nominal parameters (n = 10 cm , ≡ 2 −1 and T = 10 K, and v = 10:0 km s ), a black hole with ini- tial mass mi = 10 M ; 40 M would gain a fractional 1=2 γk^ BT −1 1=2 mass ∆M=M 10%; 57%, respectively.