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 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 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 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 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. A larger Mi cs = (1.2 km s ) T2 (6) ∼ mp ' or a smaller v would enhance the mass growth of the BH. This provides a novel channel for producing mas- is the sound speed in the MC far from the hole, v is sive stellar-mass BHs. the relative velocity between the BH and the ambient In our calculation, we have assumed that the MC is gas, λ = 1/4 (ifγ ˆ = 5/3) is the non-dimensional ac- c not destroyed due to radiative feedback throughout the cretion eigenvalue given in Bondi(1952), ˆγ = 5/3 is epoch of BH growth. Even though radiative feedback the adiabatic index, k is the Boltzmann constant, and B was found not to be important locally in the simula- T = 100 K T is the medium temperature far from the 2 tions of Park & Ricotti(2013), one still needs to justify hole. this over a much longer timescale. With our fiducial The temperature range of MCs is T (10 102) K, ∼ − parameters, the total binding energy of the MC is which corresponds to a sound speed range of cs (0.4 −1 ∼ − 1.2) km s . The proper motion velocity v of the BH GM 2  R −1 E = mc (1.9 1052 erg) M 2 mc . is likely greater than cs, so the BHL accretion rate is b mc,6 Rmc ' × 4.6 pc mainly defined by V v cs. Hereafter our nominal (9) parameter is v 10 km∼ s−1. From the scaling in Eq.(4), ∼ The dissociation energy per molecule of H2 is 4.5 eV and we immediately infer that a relatively massive BH (say, the ionization energy per atom of H is 13.6 eV, so one can M > 10M ) embedded in a MC could grow significantly estimate the total energy needed to break up molecules within the lifetime of the MC if V is small enough and n is large enough. Comparing Eqs.(3) and (4), we also M 51 Emol = 4.5 eV 4.3 10 erg Mmc,6, (10) conclude that the Eddington limit does not kick in as 2mp ' × long as the BH mass satisfies and the energy needed to break up atoms  V 3 −1 −1 M M < 150 M ε−1n5 −1 . (7) 52 10 km s Eatm = 13.6 eV 2.6 10 erg Mmc,6. (11) mp ' × A supersonic BH will remain inside a MC as long as 6 the escape velocity of the MC is higher than the BH As a result, the total energy to destroy a 10 M MC is velocity. The escape velocity of our nominal MC is 52 Etot = Eb + Emol + Eatm 4.9 10 erg. (12) ∼ × r2GM v = mc esc With the Eddington luminosity LEdd = 1.3 Rmc 1039 erg s−1(M/10 M ) of the accreting black hole, the×  R −1/2 = 43 km s−1M 1/2 mc , (8) time it takes to destroy the MC is2 mc,6 4.6 pc 6 tdestroy = Etot/LEdd = 1.2 10 yr. (13) Our nominal BH is traveling at v = 10 km s−1, so v < × vesc and the BH remains within the MC for the entire This is shorter than the BH mass growth time (20-60 lifetime of the MC as it grows its mass. Notice that Myr) in our calculation, but is much longer than the we envisage a scenario where the MC is formed at the timescale tesc for the radiation to escape the MC. The location of BHs formed from the previous generation latter can be calculated as stars. If the BH is captured by the MC, its initial speed  2 should exceed the escape velocity. A significant mass Rmc 4 Rmc tesc τ = 2.1 10 yr n5σ−21 , (14) growth is possible only if the BH can slow down due to ' c × 4.6 pc dynamical friction and remain within the MC to gain mass. 2 This estimate is likely an upper limit, since depositing heat to the Figure1 left panel shows the growth of BHs with ini- MC would alter the structure of the MC, which may decrease the ˙ tial masses m = 10 M , 40 M embedded in a MC with accretion rate. In any case, MBHL in our calculation (Eq.(4)) is i dominantly defined by the proper motion velocity V ∼ 10 km s−1, n5 = 1. The two black dashed vertical lines bracket which is one order of magnitude greater than the sound speed for the range of the lifetime of a typical MC from Jeffre- fiducial parameters. The results would not change significantly 4 son & Kruijssen(2018). One can see that a relatively unless the temperature of cloud reaches the order of 10 K. 4 Rice and Zhang

2 100

] 80

1 M [ M/M ∆

M 60

0 40 106 107 108 106 107 108 t [yr] t [yr]

Figure 1. Left: Fractional mass increase of Bondi-accreting BHs embedded in a molecular cloud, as a function of time. The two dot-dashed vertical lines represent upper and lower bound estimates from Jeffreson & Kruijssen(2018) on the lifetime of MCs in the Milky Way, from t = 19.9 Myr to t = 63.1 Myr. In each case, the BH is embedded in a molecular cloud with the following physical parameters: n = 105 cm−3, T = 102 K, and v = 10.0 km s−1 (solid) or v = 7.5 km s−1 (dashed). Blue: a BH accreting at the Eddington rate given in Eq.3. Green: a BH of mass M = 40 M accreting at the BHL rate given in Eq.4. Yellow: a BH of mass M = 10 M accreting at the BHL rate given in Eq.4. The two green dotted horizontal lines represent the fractional mass gain of ∆M/M = 45/40 = 1.125 (top) and ∆M/M = 26/40 = 0.65 (bottom) of a 40 M BH accreting to the estimated BH masses in the GW190521 system (see Abbott et al.(2020)) of 85 M and 66 M , respectively. Right: Mass (in solar masses) of Bondi-accreting BHs embedded in a molecular cloud as a function of time. The two dot-dashed vertical lines are the same as in the left panel. The two dotted horizontal lines represent the estimated GW190521 masses of 85 M and 66 M from Abbott et al.(2020). The solid lines represent the evolution track of two BHs with initial critical masses that allow the masses to grow to the GW190521 masses within the upper-bound timescale t = 63.1 Myr. The initial masses of such BHs are 41.3 M (purple) and 48.0 M (blue), respectively. The dashed lines represent BHs with those initial masses that are accreting at the Eddington rate. where 4.1. GW190521   Both pre-merger BHs in GW190521 are massive 3 Rmc τ = nσRmc 1.4 10 n5σ−21 . (15) +21 +17 ' × 4.6 pc (85−14 M and 66−18 M ). The hierarchical BH merger scenario (e.g. Fragione et al. 2020) would require that is the optical depth of X-ray photons from the disk, with each of the two components has undergone one merger in −21 2 the absorption cross section σ 10 cm σ−21 (Wilms stellar clusters prior to the detected merger. Within that ∼ et al. 2000) normalized to the cross section at 0.5 scenario, one would expect some other merging systems ∼ keV, roughly the typical photon energy of a Shakura- should have one massive component (after one merger) Sunyaev accretion disk around a 10 M BH (Shakura and one less massive component (without merger yet). ∼ & Sunyaev 1973). If the photon energy is lower (e.g. by In our scenario, the two BHs likely grow in the same MC getting reprocessed by absorption and re-emission), the and gain mass simultaneously. It is natural to obtain cross section would be higher, which would increase the two heavy BHs in one merging system. Figure1 right photon trapping time. The MC may be then destroyed panel shows that for our nominal parameters, two BHs to halt BH growth. As a result, radiative feedback would with initial masses 41.3M and 48M (both consistent not be important only under extreme conditions. It with formation from massive, especially first generation, can be indeed important for most MCs with less ex- stars; Heger et al. 2003) could grow to masses consistent treme parameters (e.g. lower n and lower Mmc), so BHs with GW190521 within the lifetime of a dense MC. can grow significantly in MCs only in a small parameter The requirement for significant mass growth is rel- space. This is consistent with the rarity of massive BHs atively slow motion between the accreting BH and the in the MW. Further detailed numerical simulations are MC. For a binary system with a total mass m m1+m2, ≡ needed to verify our suggestion for the fiducial parame- reduced mass µ m1m2/m, and symmetric mass ratio ters adopted in this paper. η µ/m, to lowest≡ order the dimensionless orbital ve- ≡ 4. CASE STUDIES Growth of stellar mass black holes via Bondi accretion 5 locity parameter is given by model would eliminate systems that might over-produce massive BH binaries. v Gm 1/3 4.2. LB-1 = πf , (16) c c3 LB-1 is a controversial system. Liu et al.(2019) ini- tially reported this black hole candidate via radial veloc- where f is the gravitational wave frequency defined ity measurements of the B-type star LB-1 and estimated through Kepler’s third law (πf)2r3 = Gm. For a sys- +11 its mass to be M1 = 68 M , with a companion star −1 −13 tem to merge within the Hubble time tH H = +0.9 +13 ≡ 0 mass and age of M2 = 8.2−1.2 M and t = 35−7 Myr, 14.4 Gyr (using the Planck parameters, Planck Collab- respectively. Abdul-Masih et al.(2020) disfavored a oration et al. 2016), the initial velocity of the two BHs black hole with mass M1 > 50 M while favoring a lower should be (assuming m1 = m2 so that η = 1/4) +0.8 mass companion star with M2 = 4.2−0.7 M . In a reply, Liu et al.(2020a) updated their black hole and com-  −1/8 v0 5 1 Gm 1 panion star mass estimates to M1 (23 65) M and = 3 ∼ − c 256 η c t M2 (5 8) M , respectively. In a follow-up study Liu ∼ −  1/8 −1/8 et al.(2020b) favor a mass ratio of M1/M2 = 5.1 0.8, −3 1/8 m t ± = 1.74 10 η0.25 ,(17) thus estimating a black hole of mass M1 (4 36) M × 100 M tH ∼ − and a companion star of mass M2 (1 6) M . The initial reported BH mass challenged∼ the stellar− evolution or v0 522 km/s. theories, and the scenario discussed in this paper would This' value is much greater than the nominal value v = help to account for its existence. Though the latest up- 10 km s−1 in our calculations. However, if the orbital per bound on M no longer challenges the stellar evolu- separation, a, of the binary is much smaller than the 1 tion theories, we suggest that the existence of 70 M Bondi radius: BHs in the Milky Way is not impossible according∼ to the 2Gm scenario discussed in this paper; but such holes should R B 2 be rare given the demanding conditions discussed above. ≡ cs   16 m −2 5. CONCLUSIONS AND DISCUSSION = 2.65 10 cm cs,6, (18) × 100 M We outlined a physical scenario for growing stellar- mass BHs if they are embedded in a dense MC and have then the BH system will accrete at a rate similar to a slow motion with respect to the MC medium. For nom- system consisting of a single BH with a mass equal to the inal parameters of n = 105 cm−3 and v = 10 km s−1 (in- total mass of the binary; see Comerford et al.(2019) and sensitive to the MC temperature), we find that BHs with Antoni et al.(2019). For a binary to merge within t H mass 10 M and 40 M can grow by ∆M/M 10% and in our nominal MC, it will always be true that a R . B 57%, respectively. We apply this scenario to GW190521,∼ For GW190521, there are two possibilities to allow and found that two relatively massive BHs with masses our envisaged scenario to happen. The first possibil- 41 M and 48 M can grow to the observed masses ity is that two relatively massive BHs (probably formed ∼from the same∼ MC, with the condition that the rela- from first generation stars) were not initially in a bi- tive velocities with respect to the MC for both BHs are nary system. After significant mass growth in an MC v 10 km s−1. Our scenario provides a novel channel for a period of < 108 yr, they form a binary system for∼ producing heavy BHs in the universe. According through dynamical interactions with other massive bod- to this scenario, the existence of a 70 M BH in the ies in the MC. The binary system subsequently loses Milky Way is not impossible, even∼ though the claimed orbital energy and angular momentum through gravita- LB-1 BH no longer needs such a large mass. Searches tional waves and merges after billions of years. The sec- for such massive BHs in the Galaxy and beyond are en- ond possibility is that the two BHs were in a binary sys- couraged. tem initially within the common Bondi radius and grow together with the Bondi accretion rate. The relative speeds between the BHs and the MC are much smaller ACKNOWLEDGMENTS than the orbital speeds so that the nominal parameters We thank an anonymous referee for constructive com- used in our calculations can be realized. Both scenarios ments, Daniel Proga and Tim Waters for discussions on require relatively contrived conditions. However, in light numerical simulations of BH accretion, and Jifeng Liu of the evidence that GW190521-like massive BH merg- for discussion on LB-1. ers are rare, the demanding conditions imposed in our 6 Rice and Zhang

REFERENCES

Abbott, R., Abbott, T. D., Abraham, S., et al. 2020, Phys. Krolik, J. H. 2004, ApJ, 615, 383, doi: 10.1086/424034 Rev. Lett., 125, 101102, Kroupa, P., Subr, L., Jerabkova, T., & Wang, L. 2020, doi: 10.1103/PhysRevLett.125.101102 MNRAS, 498, 5652, doi: 10.1093/mnras/staa2276 Abbott, R., Abbott, T. D., Abraham, S., et al. 2020, ApJL, Liu, J., Soria, R., Zheng, Z., et al. 2020a, Nature, 580, E16, 900, L13, doi: 10.3847/2041-8213/aba493 doi: 10.1038/s41586-020-2217-9 Abdul-Masih, M., Banyard, G., Bodensteiner, J., et al. Liu, J., Zhang, H., Howard, A. W., et al. 2019, Nature, 575, 2020, Nature, 580, E11, doi: 10.1038/s41586-020-2216-x 618, doi: 10.1038/s41586-019-1766-2 Aguayo-Ortiz, A., Sarbach, O., & Tejeda, E. 2020, arXiv Liu, J., Zheng, Z., Soria, R., et al. 2020b, ApJ, 900, 42, doi: 10.3847/1538-4357/aba49e e-prints, arXiv:2009.06653. McDowell, J. 1985, MNRAS, 217, 77, https://arxiv.org/abs/2009.06653 doi: 10.1093/mnras/217.1.77 Aguayo-Ortiz, A., Tejeda, E., & Hernandez, X. 2019, Natarajan, P. 2020, MNRAS, doi: 10.1093/mnras/staa3724 MNRAS, 490, 5078, doi: 10.1093/mnras/stz2989 Park, K., & Ricotti, M. 2013, ApJ, 767, 163, Antoni, A., MacLeod, M., & Ramirez-Ruiz, E. 2019, ApJ, doi: 10.1088/0004-637X/767/2/163 884, 22, doi: 10.3847/1538-4357/ab3466 Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. Bondi, H. 1952, MNRAS, 112, 195, 2016, A&A, 594, A13, doi: 10.1051/0004-6361/201525830 doi: 10.1093/mnras/112.2.195 Rice, J. R., & Zhang, B. 2017, Journal of High Energy Campana, S., & Pardi, M. C. 1993, A&A, 277, 477 , 13, 22, doi: 10.1016/j.jheap.2017.02.002 Comerford, T. A. F., Izzard, R. G., Booth, R. A., & Roupas, Z., & Kazanas, D. 2019a, A&A, 621, L1, Rosotti, G. 2019, MNRAS, 490, 5196, doi: 10.1051/0004-6361/201834609 doi: 10.1093/mnras/stz2977 —. 2019b, A&A, 632, L8, De Luca, V., Desjacques, V., Franciolini, G., Pani, P., & doi: 10.1051/0004-6361/201937002 Riotto, A. 2020, arXiv e-prints, arXiv:2009.01728. Safarzadeh, M., & Haiman, Z. 2020, ApJL, 903, L21, https://arxiv.org/abs/2009.01728 doi: 10.3847/2041-8213/abc253 Eldridge, J. J., Stanway, E. R., Breivik, K., et al. 2020, Sakstein, J., Croon, D., McDermott, S. D., Straight, M. C., MNRAS, 495, 2786, doi: 10.1093/mnras/staa1324 & Baxter, E. J. 2020, arXiv e-prints, arXiv:2009.01213. Fragione, G., Loeb, A., & Rasio, F. A. 2020, ApJL, 902, https://arxiv.org/abs/2009.01213 L26, doi: 10.3847/2041-8213/abbc0a Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 500, 33 Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Shenar, T., Bodensteiner, J., Abdul-Masih, M., et al. 2020, Hartmann, D. H. 2003, ApJ, 591, 288, A&A, 639, L6, doi: 10.1051/0004-6361/202038275 doi: 10.1086/375341 Tejeda, E., Aguayo-Ortiz, A., & Hernandez, X. 2020, ApJ, Jeffreson, S. M. R., & Kruijssen, J. M. D. 2018, MNRAS, 893, 81, doi: 10.3847/1538-4357/ab7ffe 476, 3688, doi: 10.1093/mnras/sty594 Waters, T., Aykutalp, A., Proga, D., et al. 2020, MNRAS, 491, L76, doi: 10.1093/mnrasl/slz168 Jeffreson, S. M. R., Kruijssen, J. M. D., Keller, B. W., Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914, Chevance, M., & Glover, S. C. O. 2020, MNRAS, doi: 10.1086/317016 doi: 10.1093/mnras/staa2127