Free Access
Issue
A&A
Volume 632, December 2019
Article Number L8
Number of page(s) 4
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/201937002
Published online 05 December 2019

© ESO 2019

1. Introduction

Pair-instability supernova (Heger & Woosley 2002) signifies the presence of an upper mass gap at 50−135 M of population I/II stellar origin black holes (BHs; Belczynski et al. 2016; Farmer et al. 2019). Up to the time this work is published, no BH mass in binary black hole (BBH) gravitational wave (GW) signals detected by The LIGO Scientific Collaboration & the Virgo Collaboration (2019) exceeds 50 M. Nevertheless, 75% of these BHs have masses greater than 20 M, which is far more than expected from stellar evolution and more than is observed in X-ray binaries. Fishbach & Holz (2017) suggested that because of LIGO’s higher sensitivity to BBH members ≳20 M, an excess of higher mass BHs should be expected, with a maximum at 40 M.

We suggest here that besides LIGO’s bias toward higher BH masses, the BH mass function observed will favor masses that are higher than predicted by stellar evolution because of gas accretion by BHs in dense stellar clusters. This will occur rapidly (Roupas & Kazanas 2019) before the gas is depleted through the first stellar formation event. More importantly, we estimate that this shift is sufficient to fill the upper mass gap. This mechanism is different, but not mutually exclusive, from the repeated mergers scenario (Gerosa & Berti 2017, 2019; Rodriguez et al. 2019; Doctor et al. 2019).

It is plausible that globular clusters (GCs) started life as dense gas clouds, which are referred to as protoglobular or primordial clusters. They underwent prolonged star formation early in their lifetimes (Gratton et al. 2012). A huge gas reservoir, comparable to or higher in mass than that of its stellar component, is available for accretion immediately after the first stellar formation event. Feedback processes from stellar evolution in star-forming regions are believed to clear away the surrounding gas (e.g., Voss et al. 2010; Galván-Madrid et al. 2013; Krumholz et al. 2014). The effectiveness of the process depends on the compactness (Krause et al. 2012, 2016; Silich & Tenorio-Tagle 2017, 2018), that is, the mass over half-mass radius of the cluster,

(1)

For sufficiently compact clouds, feedback processes are expected to become ineffective because they are proportional to the total mass. The gravitational binding energy is proportional to the square of it.

Proposals regarding the precise mechanism for gas depletion in primordial GCs abound in the literature, but a general consensus has not been achieved (Spergel 1991; Thoul et al. 2000; Fender et al. 2005; Moore & Bildsten 2011; Herwig et al. 2012; Krumholz et al. 2014; Krause et al. 2016, 2013, 2012; Silich & Tenorio-Tagle 2017, 2018; Marks et al. 2008; Kruijssen 2012; D’Ercole et al. 2008; Conroy 2012; Renzini et al. 2015; Bagetakos et al. 2011; Jaskot et al. 2011; Fierlinger et al. 2016; Yadav et al. 2017). Leigh et al. (2013) proposed that in any cluster that is able to form massive stars, the primordial gas is depleted exactly due to the accretion onto BHs. They find that accreting BHs can deplete the whole gas reservoir within as short a time as 10 Myr.

Here, we do not investigate this specific scenario and do not focus on the gas reservoir, but on the effect of accretion onto the BHs of the cluster. Calura et al. (2015) estimated that within ∼14 Myr the gas is depleted by 99% by star formation feedback processes in a primordial cluster with an initial total mass of about 107M. Typical globular clusters (GCs) should form about 102–103 BHs within about 3 Myr (Morscher et al. 2013, 2015) and could retain most of them initially if their natal kicks are sufficiently low (see Wong et al. 2012 and references therein). Additional observational evidence suggests that a BH subcluster may be retained even to present-day GCs (Sedda et al. 2019; Abbate et al. 2019).

We assume in our analysis an initial population of 500 BHs, segregated in the core of the cluster (Spitzer 1987; Leigh et al. 2014), and that the primordial gas is depleted by 99% within time tf = 10 Myr, following an exponential law,

(2)

Therefore our analysis applies to any type of depletion mechanism as long as it does not proceed faster than an exponential law and the loss is approximately uniform in the core. To the exponential gas loss, we add the loss by accretion onto the BHs.

In the next section we briefly describe our model, and in Sect. 3 we present the results of our analysis. In the final section we discuss our conclusions.

2. Model

We considered N = 500 BHs moving inside a fixed external potential with gaseous and stellar components. This analysis focuses solely on the effect of accretion onto the BH mass function, and in addition, of accretion within a short timescale tf = 10 Myr. We do not study the general effect on the cluster or BH kinematics. For simplicity and with the aim to provide only statistical estimates, we therefore assumed a spherical distribution of BHs at any t, that is, we did not follow the angular changes of their orbital planes (the BH subsystem dynamics is dominated by the gaseous and stellar components on the timescales considered). We further neglected close encounters. We discuss this further in our conclusions. We intend to include and study both of these effects in a separate more detailed work that will include the estimation of merger rates.

We assumed that the fixed external potential is generated by Plummer density profiles for stellar and gaseous components,

(3)

(4)

where ϵ is the stellar formation efficiency, and τ, given in Eq. (2), corresponds to 99% gas depletion by the final time tf. Following Silich & Tenorio-Tagle (2018), we assumed that the gas is dominated by turbulence (Johnson et al. 2015; Elmegreen 2017), in which case the equation of state is .

The BHs were chosen from a Salpeter initial mass function (Perna et al. 2019) with a cutoff

(5)

with A = 1.35/(5−1.35 − 50−1.35). They were all assumed to be initially bound inside the core radius rc = 0.64a of the initial Plummer sphere. The initial radial distribution and velocities of the BHs were chosen randomly from a Plummer distribution with the same softening radius a of Eqs. (3) and (4). Their initial angular distribution was spherical.

Hot-type accretion, whose typical representative is Bondi-Hoyle accretion, is appropriate in the dense gaseous environment we examined, where the speed of sound is greater than or on the order of the accretor velocity (Merritt 2013). The spherically symmetric accretion rate for a given cross section may be written as (Frank et al. 1985)

(6)

where m is the mass of the BH that accretes the gas with relative velocity

(7)

and v is the BH velocity with respect to the center of mass of the cluster, and cs is the speed of sound of the gas.

As Bonnell et al. (2001) pointed out, choosing the proper accretion radius inside a gaseous stellar cluster depends on the relative amount of gas and on the radial position of the accretor. When the gas dominates the cluster potential, the accretion rates are given by a tidal-lobe accretion radius (Paczyński 1971),

(8)

where M(r <  ri) is the total mass of the cluster within the radial position ri of the ith BH. When the gas has been sufficiently depleted so that stars dominate the potential, the appropriate accretion radius is the Bondi – Hoyle radius,

(9)

Following Bonnell et al. (2001), and with the intention to provide a minimum estimate of the BH mass growth, we chose the accretion radius to be the smaller of the two,

(10)

at any instant of time.

We accounted in the analysis for the dynamical friction generated by the stellar component according to Chandrasekhar’s formula

(11)

where v is the velocity of the BH, and the Coulomb logarithm is lnΛ = ln(bmax/bmin) with

(12)

Lee & Stahler (2011, 2014) reported that the gaseous dynamical friction on an accretor is

(13)

for both subsonic and supersonic accretors. This result seems to agree with the earlier calculation of Hadjidemetriou (1963). It may be understood simply as manifesting angular momentum preservation (see also Roupas & Kazanas 2019). Lee & Stahler (2011, 2014) argued that Eq. (13) encompasses the entire gaseous dynamical friction. In contrast, Tanaka & Haiman (2009) added to this term the gaseous dynamical friction formula proposed by Ostriker (1999; see also Ruffert 1996; Sánchez-Salcedo & Brandenburg 2001; Kim & Kim 2009; Indulekha 2013; Antoni et al. 2019). This issue seems unresolved. We adopted the former view (only the term in Eq. (13), neglecting Ostriker’s term) because it generates the least gaseous dynamical friction. We wish to provide minimum estimates of BH growth, and higher dynamical friction causes BHs to sink deeper into the center where the gas density is higher, leading to more intense BH growth. Including Ostriker’s formula would only lead to results that support our conclusions more strongly. We can also report that when Ostriker’s formula is used instead of Eq. (13), we obtained similar results.

Finally, we describe the equations of motion of the BH-subcluster. As we noted, the BHs were assumed to be initially spherically distributed, and we are interested only in the amount of gas that they can accrete within 10 Myr. Therefore we did not follow the angular redistribution of their orbital planes and also assumed that the background of stars and gas is spatially fixed, as in Eqs. (3)–(4). We further assumed that the gas is depleted by 99% uniformly and exponentially within tf = 10 Myr due to any process. The BHs accrete gas according to Eq. (6) and are subject to dynamical friction generated by stars, Eq. (11), as well as gas, Eq. (13). The system of equations that we solved numerically for i = 1, …, NBH was therefore

(14)

(15)

(16)

(17)

(18)

The total mass within ri, M(r <  ri(t)), includes the enclosed mass of field stars, of the remaining enclosed gas at t, and the enclosed mass within ri of the BH population.

3. Results

We performed the numerical analysis for an initial total cluster mass Mtot = 106M, stellar formation efficiency ϵ = 0.3, and several initial compactness values Cini(t = 0) = {5−10}, corresponding to half-mass radii rhm = {2−1} pc. For each half-mass radius we performed ten simulations with different initial conditions. We also assumed a depletion time tf = 10 Myr.

In Fig. 1 we show the final BH mass function in the case of a half-mass radius rhm = 1 pc, that is, an initial compactness Cini = 10, which results in a final compactness Cfin = 3. It is evident that about 2% of the BHs fill in the theoretical BH mass gap. In addition, the BH mass function increases for mass values higher than ≳35 M. This is also true for any rhm ≤ 1.25 pc. In Fig. 2 we draw the maximum BH mass that is achieved after accretion with respect to the half-mass radius of the cluster.

thumbnail Fig. 1.

Total BH mass distribution of ten samples of N = 500 BHs, with an initial mass cutoff at 50 M, bound in the core of a primordial gaseous stellar cluster with rhm = 1 pc. We assumed the gas to be exponentially depleted by 99% within 10 Myr.

Open with DEXTER

thumbnail Fig. 2.

Maximum BH mass achieved after gas accretion of an initial population of N = 500 BHs, with initial mass cutoff at 50 M, with respect to the half-mass radius rhm of the stellar cluster for ten samples of the BH population at each rhm. The vertical lines represent the range of maximum mass values achieved in all samples for a certain rhm. The gas is assumed to be exponentially depleted by 99% within tf = 10 Myr.

Open with DEXTER

We report that the probability for a BH to exceed the initial cutoff of 50 M is

(19)

where the lower half-mass radius corresponds to the upper probability value. The total accreted gas mass lies in the range ≃{100−900} M (with an initial mass 5400 M of the BH population).

We remark that the more massive BHs sink deeper into the center r ≲ 10−3 pc, where it is most probable that they form hard massive BBHs. About half of the samples at each rhm form a BBH by this process. In particular, for rhm = 1 pc, six hard BBHs were formed in five out of the ten samples, as shown in Table 1. The BBHs will continue to form for t >  10 Myr by dynamical processes.

Table 1.

In five out of the ten samples with rhm = 1 pc, six massive BBHs formed within 10 Myr because they sank to the center of the cluster (≲10−3 pc).

4. Conclusions

We propose that in the early life of dense stellar clusters a rapid accretion process operates. The BHs that have been generated by the more massive stars, most of which are therefore segregated in the core, rapidly accrete primordial gas that is also accumulated in the core within ∼10 Myr before it becomes depleted by stellar formation feedback processes. The amount of accreted gas depends sensitively on the mass and compactness of the cluster and on the accretion timescale.

The BH (and BBH) initial mass function are shifted toward higher values. For sufficiently compact clusters the BH mass limit of 50 M predicted by supernova theory is exceeded. We estimate that 1−2% of the initial BH population exceeds the theoretical limit for an initial cluster mass Mtot = 106M, half-mass radii rhm = {1, 2} pc, depletion time tf = 10 Myr, and stellar formation efficiency ϵ = 0.3.

In addition, the BH mass growth together with dynamical friction causes the more massive BHs to sink deeper into the center of the cluster. There, it is more probable that the BH attains mutual binding energies that are higher than each individual binding energy with the cluster and thus forms hard massive BBHs, as listed in Table 1. This process operates rapidly in the very early life of the cluster, therefore it adds to the dynamical channel of BBH formation and mergers (Abbott et al. 2016), which will continue to operate during the whole life of the cluster (e.g., see Park et al. 2017; Hong et al. 2018). Thus, an observation by LIGO of BBH mergers with member masses above 50 M is to be anticipated.

We assumed an initial BH population bound inside the core of the cluster neglecting close encounters. These are not expected to have a strong effect on our results because their timescale are relatively long. The recoil mechanism of BBH-BH three-body encounters, which ejects BBHs out of the cluster (Sigurdsson & Hernquist 1993), operates on gigayear timescales, but our proposed mechanism works on megayear timescales. In addition, it has become evident during the past decade (Maccarone et al. 2007, 2011; Barnard et al. 2011; Shih et al. 2010; Strader et al. 2012) that recoil is not effective in sufficiently dense GCs, which contain a population (≲1000) of BHs and BBHs in their core (e.g., Morscher et al. 2015 and references therein).

More sophisticated simulations may be required in order to confirm with higher confidence and detail the effects of the rapid accretion process we proposed here and to provide a better estimation of the resulting BH mass function for a wider range of parameter values {Mtot, rhm, tf, ϵ}. Nevertheless, our current analysis strongly supports the idea that the BH upper mass gap can be populated by rapid gas accretion onto the BHs of dense primordial stellar clusters.

References

All Tables

Table 1.

In five out of the ten samples with rhm = 1 pc, six massive BBHs formed within 10 Myr because they sank to the center of the cluster (≲10−3 pc).

All Figures

thumbnail Fig. 1.

Total BH mass distribution of ten samples of N = 500 BHs, with an initial mass cutoff at 50 M, bound in the core of a primordial gaseous stellar cluster with rhm = 1 pc. We assumed the gas to be exponentially depleted by 99% within 10 Myr.

Open with DEXTER
In the text
thumbnail Fig. 2.

Maximum BH mass achieved after gas accretion of an initial population of N = 500 BHs, with initial mass cutoff at 50 M, with respect to the half-mass radius rhm of the stellar cluster for ten samples of the BH population at each rhm. The vertical lines represent the range of maximum mass values achieved in all samples for a certain rhm. The gas is assumed to be exponentially depleted by 99% within tf = 10 Myr.

Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.