Dynamics of stellar black holes in young star clusters with different metallicities - I. Implications for X-ray binaries
We present body simulations of intermediate-mass ( M) young star clusters (SCs) with three different metallicities (, 0.1 and 1 Z), including metal-dependent stellar evolution recipes and binary evolution. Following recent theoretical models of wind mass loss and core collapse supernovae, we assume that the mass of the stellar remnants depends on the metallicity of the progenitor stars. In particular, massive metal-poor stars ( Z) are enabled to form massive stellar black holes (MSBHs, with mass M) through direct collapse. We find that three-body encounters, and especially dynamical exchanges, dominate the evolution of the MSBHs formed in our simulations. In SCs with and 0.1 Z, about 75 per cent of simulated MSBHs form from single stars and become members of binaries through dynamical exchanges in the first 100 Myr of the SC life. This is a factor of more efficient than in the case of low-mass ( M) stellar black holes. A small but non-negligible fraction of MSBHs power wind-accreting ( per cent) and Roche lobe overflow (RLO, per cent) binary systems. The vast majority of MSBH binaries that undergo wind accretion and/or RLO were born from dynamical exchange. This result indicates that MSBHs can power X-ray binaries in low-metallicity young SCs, and is very promising to explain the association of many ultraluminous X-ray sources with low-metallicity and actively star forming environments.
keywords:black hole physics – stars: binaries: general – galaxies: star clusters: general – X-rays: binaries – methods: numerical – stars: kinematics and dynamics.
The mass spectrum of black holes (BHs) that form from the collapse of massive stars is highly uncertain. An accurate dynamical mass estimate has been derived only for stellar BHs (see table 2 of Özel et al. 2010 for one of the most updated compilations, but see also Lee, Brown & Wijers 2002; Orosz 2003; Narayan & McClintock 2005). Most of the derived BH masses are in the range , with an apparent absence (in X-ray binaries) of BHs with M, difficult to explain with observational biases (Özel et al. 2010). In the Milky Way (MW), the most massive BHs in X-ray binaries do not seem to significantly exceed M (e.g. M for GS2023+338, Charles & Coe 2006; M for GRS 1915+105, Harlaftis & Greiner 2004), whereas a few BHs in nearby galaxies may have higher masses: M33 X-7 ( M, Orosz et al. 2007), IC 10 X-1 ( M, Prestwich et al. 2007; Silverman & Filippenko 2008) and NGC 300 X-1 ( M, Carpano et al. 2007; Crowther et al. 2007; M, under reasonable assumptions for the inclination and for the mass of the companion, Crowther et al. 2010).
Interestingly, these three relatively massive stellar BHs are hosted in regions with relatively low metallicity. A metallicity Z is estimated for the dwarf irregular galaxy IC 10 from an electron-temperature based calibration of spectra of HII regions (, Garnett 1990), assuming Z. The metallicity of M 33 in proximity of X-7 (i.e. at , where is the Holmberg radius) is Z, and that of NGC 300 in the vicinity of X-1 (i.e. at ) is Z (derived from the metallicity gradients of M 33 and of NGC 300, respectively, provided by Pilyugin, Vílchez & Contini 2004).
From a theoretical perspective, our knowledge of the mass spectrum of stellar BHs is hampered by the uncertainties about two issues: (i) mass losses by stellar winds in massive stars (during and especially after the main sequence, MS); (ii) the hydrodynamics of core-collapse supernova (SN) explosions.
The rate of mass loss by stellar winds during the MS likely increases with the metallicity of the star (, where , depending on the model, e.g. Kudritzki, Pauldrach & Puls 1987; Leitherer, Robert & Drissen 1992; Maeder 1992; Kudritzki & Puls 2000; Vink, de Koter & Lamers 2001; Kudritzki 2002; Belkus, Van Bever & Vanbeveren 2007; Pauldrach, Vanbeveren & Hoffmann 2012). The behaviour of post-MS massive stars, such as luminous blue variable stars (LBVs) and Wolf-Rayet stars (WRs), is much more uncertain (e.g. Vink & de Koter 2005).
According to models of stellar evolution and SN explosion (Fryer 1999; Fryer & Kalogera 2001; Woosley, Heger & Weaver 2002; Heger et al. 2003a), a star with a final mass
In the case of a failed SN, it is reasonable to expect that the mass of the remnant is comparable to the final mass of the progenitor star (or at least more than half of it, see e.g. Heger et al. 2003b). Therefore, BHs may form, via this channel, with mass higher than in the case of SN explosion. According to Belczynski et al. (2010; hereafter B10), the mass of a BH formed via direct collapse may be as high as M. The models by B10 consider only the evolution of single, non-rotating stars. Rotation and binarity can affect the final mass of the remnant (e.g. Vanbeveren 2009; Maeder & Meynet 2010, 2012).
In the following, we will refer to massive stellar BHs (MSBHs) to indicate BHs with mass M formed via direct collapse. The existence of MSBHs in the nearby Universe may be crucial for our understanding of X-ray sources. The scenario of X-ray binaries powered by MSBHs was recently proposed to explain a large fraction of the ultra-luminous X-ray sources (ULXs, i.e. X-ray sources with luminosity, assumed isotropic, higher than erg s), without requiring excessive super-Eddington factors or more exotic mechanisms (e.g. Mapelli, Colpi & Zampieri 2009; Zampieri & Roberts 2009; Mapelli et al. 2010, 2011a).
In this paper, we present new body simulations of dense intermediate-mass (a few M) young ( Myr) star clusters (SCs), including an accurate treatment of dynamics and updated recipes for metal-dependent stellar evolution, stellar winds and failed SNe. Our aim is to study the formation and dynamical evolution of stellar BHs and MSBHs in young SCs with different metallicity. In the current paper (which is the first of a series), we will focus on the effects of the dynamics of stellar BHs and MSBHs on the population of X-ray sources. In the next papers of the series we will consider also other effects of BH dynamics (e.g. the consequences for the population of gravitational wave sources).
We simulate SCs, as most stars ( per cent) are expected to form in SCs (e.g. Lada & Lada 2003). We restrict our analysis to intermediate-mass (a few M) young SCs. These form with a higher frequency than larger clusters in the nearby Universe (as the mass function of star clusters is , with , Lada & Lada 2003), and are often sites of an intense X-ray activity: many bright high-mass X-ray binaries (HMXBs) and ULXs are associated with OB associations and with young intermediate-mass SCs (e.g. Goad et al. 2002; Zezas et al. 2002; Liu, Bregman & Seitzer 2004; Soria et al. 2005; Ramsey et al. 2006; Terashima, Inoue & Wilson 2006; Abolmasov et al. 2007; Berghea 2009; Swartz, Tennant & Soria 2009; Tao et al. 2011; Grisé et al. 2011, 2012). Tens of intermediate-mass young SCs have been discovered in the MW in the last few years (e.g. Bica et al. 2003; Mercer et al. 2005; Borissova et al. 2011; Richards et al. 2012)
The paper is organized as follows. In Section 2, we briefly review the previous work on this topic, to summarize the state-of-art and to highlight the differences with the present analysis. In Section 3, we describe the method adopted for the simulations. In Section 4, we present the results, focusing on the mass spectrum of BHs and on the effects of dynamics on accreting BHs in X-ray binaries. In Section 5, we summarize the most relevant results and discuss future challenges for body simulations of SCs.
2 Short review of previous work: the importance of combining dynamics and stellar evolution
Most population synthesis codes study the formation of BHs from single stars and from stars in primordial binaries (i.e. stars that are in the same binary since their formation). This method has been widely used to investigate the population of X-ray binaries (e.g. Portegies Zwart, Verbunt & Ergma 1997; Hurley, Tout & Pols 2002; Podsiadlowski, Rappaport & Pfahl 2002; Podsiadlowski, Rappaport & Han 2003; Belczynski et al. 2004a; Belczynski, Sadowski & Rasio 2004b; Rappaport, Podsiadlowski & Pfahl 2005; Dray 2006; Madhusudhan et al. 2006, 2008; Belczynski et al. 2008; Linden et al. 2010, hereafter L10). Some studies include also recipes for the effects of metallicity on stellar evolution (e.g. Hurley et al. 2002; Belczynski et al. 2004a; Dray 2006; Belczynski et al. 2008; L10).
In particular, L10 adopt recipes for stellar winds and failed SNe that are very similar to B10, and find that MSBHs are unlikely to power bright X-ray binaries. This occurs because primordial binaries merge if they are sufficiently tight to enter a common envelope (CE) phase before the SN of the primary and if they are still in Roche lobe overflow (RLO) at the end of the CE phase. If they survive the CE phase, they may be tight enough to later enter a RLO phase when the secondary evolves, but typically their BHs are small, as a consequence of the mass lost during CE. In both cases, the binary does not evolve into a RLO HMXB with a MSBH.
On the other hand, L10 and the large majority of studies that take into account metallicity do not include the effects of dynamics on the evolution of primordial binaries and on the formation of new binaries. This is a strong limitation, as most stars ( per cent) likely form in SCs (e.g. Lada & Lada 2003), and most SCs are collisional environments, that is sites where close encounters between stars and binaries (three-body encounters) are extremely frequent and have important consequences (see, e.g. Bonnell & Kroupa 1998, for a review dedicated to young SCs). Binaries are important for close encounters, because they have a larger cross-section than single stars and because they have an energy reservoir (their internal energy), which can be exchanged with single stars (e.g. Heggie 1975; Heggie & Hut 1993; Davies 1995; Colpi, Mapelli & Possenti 2003). Close encounters with single stars statistically lead to the increase (decrease) of the binding energy of a hard (soft) binary, defined as a binary with binding energy higher (lower) than the average kinetic energy of a star in the SC (Heggie 1975). Close encounters can even unbind binaries (ionization). Dynamical exchanges are also possible, that is interactions where one of the former members of the binary is replaced by the single star (for a description of the possible outcomes of a binary-single star interaction, see, e.g. Sigurdsson & Phinney 1993). Finally, recoil velocities due to three-body encounters
can cause the ejection of the star and/or of the binary from the parent cluster (see Section 4.1 of Sigurdsson & Phinney 1993 for a general definition of recoil velocity in three-body encounters). Recent body simulations of dense young SCs with MSBH binaries
Blecha et al. (2006, hereafter B06) is one of the few studies of massive BHs where a stellar evolution code (although not accounting for different metallicities) is combined with recipes for dynamics (in part semi-analytical prescriptions and in part a body integrator for three-body encounters). B06 study the evolution of an intermediate-mass BH (IMBH) born via runaway collapse (Portegies Zwart & McMillan 2002) at the centre of a 5 M dense young SC. The IMBH mass in the simulations of B06 is in the range M, i.e. partially overlapping with our definition of MSBHs. B06 find that the dynamical effects are very important for the evolution of the IMBH, as captures of stars and three-body encounters allow it to have a companion star for most of the SC lifetime.
The main differences between this paper and B06 are the following. Firstly, our treatment of dynamics is fully body, rather than based on a hybrid code. Secondly, in our simulations MSBHs form self-consistently, as a consequence of metallicity-dependent stellar evolution and dynamics, whereas B06 generate one IMBH per cluster, following recipes from runaway-merger simulations. Thus, ‘our’ MSBHs are generally smaller than those by B06, and do not form necessarily in the core. Thirdly, we included a metallicity-dependent treatment of stellar evolution and wind-mass losses. At last, the total mass of each of the SCs we simulate is a factor of smaller than in B06 simulations.
Finally, Monte Carlo codes are suitable for the study of the largest collisional systems ( objects), such as globular clusters (e.g. Hénon 1971a, 1971b; Stodólkiewicz 1982, 1986; Giersz 1998; Joshi et al. 2000, 2001; Fregeau et al. 2003; Fregeau & Rasio 2007). The number of objects in these systems makes prohibitive to run wide grids of body models so far. Some Monte Carlo codes also include accurate recipes for metal-dependent stellar and binary evolution (e.g. Chatterjee et al. 2010; Downing et al. 2010, 2011; Pattabiraman et al. 2012). Monte Carlo codes were recently used to study the population of X-ray binaries in globular clusters (e.g. Ivanova et al. 2006, 2008), the evolution of possible sources of gravitational waves (e.g. Downing et al. 2010, 2011), and the formation of IMBHs by runaway collapse (e.g. Gürkan, Freitag & Rasio 2004). These studies do not include recipes for the formation of MSBHs in the M range.
3 Method and simulations
The simulations were done using the Starlab
3.1 Single star evolution
We added the metallicity dependence of stellar radius, temperature and luminosity, using the polynomial fitting formulas by Hurley, Pols & Tout (2000). We changed the recipes for mass loss by winds for MS stars, by using the metal-dependent fitting formulas given by Vink et al. (2001; see also B10).
We added a very approximate treatment for LBV stars and for WR mass losses by stellar winds, according to the recipes by B10. In particular, we assume that a post-MS star becomes a LBV when its luminosity and radius satisfy the requirement that and (Humphreys & Davidson 1994). The mass-loss rate by stellar winds for a LBV is then calculated as M yr, where (chosen by B10 because it allows to reproduce the most massive known stellar BHs, see the Introduction and Özel et al. 2010).
Naked helium giants coming from stars with zero age MS (ZAMS) mass M (e.g., van der Hucht 1991 and references therein) are labelled as WR stars in the new version of the code and undergo a mass-loss rate by stellar winds defined by M yr, where . This formula was first used by B10, and is a combination of the Hamann & Koesterke (1998) wind rate estimate (taking into account WR wind clumping) and Vink & de Koter (2005) wind -dependence for WRs.
We assume that all stars with final mass M collapse quietly to a BH, without SN explosion (Fryer 1999; Fryer & Kalogera 2001), and that all the BHs born in this way do not receive any natal kicks (see Fryer et al. 2012, for a discussion of this assumption). If M, the mass of the BH is equal to , where is the final mass of the Carbon Oxygen (CO) content of the progenitor, while and are the residual mass of Helium (He) and of Hydrogen (H), respectively. is the fraction of He and H mass which collapses to the BH in the failed SN scenario. The value of is uncertain, and can range between 0 and 1 (depending whether a faint SN or completely no SN occurs). We assume to match the maximum values of at low derived by B10.
If M, the SN takes place and we adopt the standard routine of Starlab (Portegies Zwart & Verbunt 1996), including a SN explosion and partial fallback of H, He and CO, depending on the ratio between energy released by the SN and binding energy of each shell.
Fig. 1 shows the mass of the BH as a function of the ZAMS mass for four different metallicities (from 0.01 to 1 Z, where we adopt Z). The values of shown in Fig. 1 have been obtained from our code when binary evolution is switched off. The main feature is the increase of the maximum allowed for decreasing metallicity. It is interesting to note the abrupt step in for and M, in the case of Z and Z, respectively. The step is produced by the minimum value of for which M and the direct collapse takes place, at a given metallicity.
The main differences between our treatment of the collapse and that by B10 are that (i) B10 do not use the total final mass to discriminate between direct collapse and ‘standard’ SN, but the mass of the CO core, assuming that all stars with final CO mass M collapse directly to a BH; (ii) B10 adopt slightly different recipes for fallback (see Portegies Zwart & Verbunt 1996 and Belczynski et al. 2008 for details).
Because of these differences, in our simulations fallback is less efficient for ZAMS masses in the M range. Furthermore, the maximum mass of BHs for is higher (by a factor of 1.4) in our simulations with respect to B10 results. Our choice of using to discriminate between SN and direct collapse produces the abrupt step visible in Fig. 1 for low , which is much less pronounced in B10. Our results are consistent with those of B10 in light of the uncertainty in the BH mass spectrum. In general, our stellar evolution recipes are less sophisticated than those adopted by B10, but are the most accurate that can be presently implemented within the framework of a complete body calculation, maintaining an acceptable computational time.
One of the simplifying approximations in our model (as well as in B10) is that stars do not rotate. Rotation may have important effects both on the mass loss of the progenitor star, inducing rotational mixing (e.g., Maeder & Meynet 2010, 2012, and references therein), and on the formation of the BH, producing asymmetries in the collapse (e.g. Fryer & Heger 2000; Akiyama et al. 2003; Fryer & Warren 2004; Ardeljan, Bisnovatyi-Kogan & Moiseenko 2005; Wheeler, Akiyama& Williams 2005; Thompson, Quataert & Burrows 2005; Woosley & Bloom 2006). On the other hand, how the effects of rotation affect the mass spectrum of BHs is debated and highly uncertain (e.g. Dessart, O’Connor, & Ott 2012; Maeder & Meynet 2012). Including the effects of rotation in our model requires a significant revision of the adopted recipes for stellar evolution and wind mass losses, which goes beyond the aims of the current paper. We will investigate the effects of rotation in a forthcoming paper.
3.2 Binary evolution
For the evolution of binaries we maintain the recipes already present in the original version of SeBa (Portegies Zwart & Verbunt 1996).
In particular, SeBa distinguishes among the evolution of detached binaries (where both stars are smaller than their Roche lobes), that of semi-detached binaries (one of the components fills its Roche lobe) and that of contact binaries (both stars fill their Roche lobe). In the case of detached binaries, mass loss and accretion can occur via stellar winds (according to the formulation by Livio & Warner 1984). In the case of semi-detached binaries, either stable or unstable RLO takes place, according to the formalism described in appendix C of Portegies Zwart & Verbunt (1996). If the accretor is a BH, the maximum mass-accretion rate of the BH is constrained by the requirement that the luminosity does not exceed the Eddington limit.
A further critical assumption about binary evolution is the CE efficiency (we define as CE efficiency the product of the two degenerate parameters and , see equation 3 of Podsiadlowski et al. 2003 for a standard definition). For the runs presented in this paper, we adopt , which is a rather standard value and tends to favour the formation of BH binaries (e.g. Podsiadlowski et al. 2003). Test runs with different values of show that the choice of this parameter does not significantly affect the results for values , in agreement with Podsiadlowski et al. (2003) and with L10.
|0.01, 0.1, 1.0|
: central adimensional potential in the King (1966) model; : number of stars per cluster; : initial core radius; : concentration ( is the initial tidal radius); and : minimum and maximum simulated stellar mass, respectively; : fraction of primordial binaries, defined as the number of primordial binaries in each SC divided by the number of ‘centres of mass’ (CMs) in the SC. In each simulated SC, there are initially 5000 CMs, among which 500 are designated as ‘binaries’ and 4500 are ‘single stars’ (see Downing et al. 2010 for a description of this formalism). Thus, 1000 stars per SC are initially in binaries.
3.3 Initial conditions and simulation grid
In this paper, we focus on moderately dense SCs, adopting a spherical King profile with central adimensional potential (King 1966). Each simulated SC is initially composed of stars, corresponding to a total mass M per SC. The resulting core density, at the beginning of the simulation, is M pc. The main parameters adopted for the initial conditions are reported in Table 1.
For the runs presented in this paper, we fix the primordial binary fraction to . The single stars and the primary stars of each binary are generated according to a Kroupa initial mass function (IMF, Kroupa 2001), with minimum and maximum mass equal to 0.1 and 150 M, respectively
The simulated SCs have half-mass relaxation time , where is the initial half-mass radius of the SC (in our simulations pc). Thus, the core collapse time (Portegies Zwart & McMillan 2002) is . We integrate the evolution of these SCs for the first 100 Myr, therefore studying a phase of the life of the cluster in which dynamical interactions are particularly intense.
We make three sets of runs corresponding to three different metallicities: 0.01 Z, 0.1 Z and 1 Z. For each of these metallicities we simulate 100 different clusters, for a total of 300 SCs.
The properties of the simulated SCs (total mass, number of stars, core density, core and half-mass radius) are consistent with the properties of observed young intermediate-mass SCs (see e.g. the recent review by Portegies, McMillan & Gieles 2010; see also Hillenbrand & Hartmann 1998; Dias et al. 2002; Portegies Zwart 2004; Pfalzner 2009; Kuhn et al. 2012).
Finally, our simulations do not include recipes for the tidal field of the host galaxy. Accounting for the tidal field may increase the fraction of mass lost, and even transform the SCs into unbound associations (e.g. Gieles & Portegies Zwart 2011). The effect of tidal fields will be added and discussed in forthcoming papers.
The simulations were run on the graphics processing unit (GPU) cluster IBM-PLX at CINECA. The processors available on PLX are six-cores Intel Westmere 2.40 GHz (two per node), while the GPUs are NVIDIA Tesla M2070 and M2070Q (two per node). Each single job ran over one processor and two GPUs and required 50 CPU hours on average. Starlab runs on GPUs through the SAPPORO library (Gaburov, Harfst & Portegies Zwart 2009).
4 Results and discussion
4.1 Effects of binary evolution of the progenitor star on BH mass
Fig. 2 shows the effects of binary evolution of the progenitor star on the final mass of the BH. In particular, the points in Fig. 2 show all the simulated BHs (originated from stars in binaries) whose mass differs by more than per cent from the BH mass calculated if the progenitor was a single star evolving in isolation. In the following, we will simply refer to them as BHs with . All the BHs with form from progenitors that were in a primordial binary and underwent a mass transfer (MT) phase before collapsing to a BH. The effects of MT before the formation of the first BH can be dramatic, especially when the two stars undergo a CE phase before the first SN. Depending on the binding energy of the envelope, a CE phase can end either with the ejection of the envelope or with the merger of the two stars (see Portegies Zwart & Verbunt 1996 for more details). If the envelope is ejected, the final BH mass will likely be smaller than expected for the evolution of an unperturbed progenitor; whereas if the two stars merge, the final BH mass will be higher than expected for a single star. For simplicity, the code assumes that, if two stars merge during a CE phase, no mass is lost during the merger and that the merger remnant evolves as a MS star (the merger product is a blue straggler star, if its mass is higher than the turn-off mass).
For example, a CE phase followed by envelope ejection is responsible for the two BHs formed by stars with M and Z, which have M, i.e. a factor of less than their analogues born from single stars or wider binaries.
Interestingly, for five BHs form with mass M (one of them with mass close to 40 M). These five BHs all form from the merger of a primordial binary before the SN explosion of the primary star
Table 2 allows to understand the statistical importance of binary evolution on the final mass of BHs. We define as the fraction of BHs with with respect to the total number of simulated BHs. From Table 2, for all the metallicities. Since the primordial binary fraction is , approximately all BHs that formed in primordial binaries have . This depends on the chosen initial distribution of the semi-major axes and on the adopted recipes for CE, but also on the effects of dynamical encounters, as we will discuss in the next Section. This result indicates that the initial binary fraction is an essential ingredient to shape the mass distribution of BHs.
The resulting mass distribution of BHs for the three considered metallicities is shown in Fig. 3. The low-mass tails of the three distributions do not differ significantly, showing a peak at M. This is in fair agreement with the observational limits for MW BHs discussed by Özel et al. (2010), especially considering the uncertainties on fallback models.
The most significant difference among different metallicities appears at the high-mass tail: the most massive BHs at Z have a cut-off at M, whereas a significant fraction of BHs form with M at Z and with M at Z.
: average number of BHs per cluster; : fraction of BHs that are members of a binary at least once in the simulations. This and all the other fractions reported in this Table (except for ) are calculated with respect to the total number of simulated BHs; : fraction of BHs that form from single stars and become members of a binary at least once in the simulations; : average number of exchanges per binary that hosts at least one BH (hereafter BH binary); : average BH binary lifetime; : fraction of BHs that undergo wind accretion at least once in the simulated time interval; : fraction of BHs that undergo wind accretion with an exchanged companion; : fraction of BHs that undergo RLO at least once in the simulated time interval; : fraction of BHs that undergo RLO with an exchanged companion; : average number of strong interactions per BH binary (considering only binaries that will undergo RLO); : fraction of BHs whose mass is affected by per cent by binary evolution and dynamics; : fraction of MSBHs (i.e. BHs with ); : fraction of MSBHs that undergo RLO (RLO MSBHs) with respect to the total number of BHs undergoing RLO (RLO BHs).
: average number of MSBHs per cluster; : fraction of MSBHs that are members of a binary at least once in the simulations. This and all the fractions reported in this Table are calculated with respect to the total number of simulated MSBHs; : fraction of MSBHs that form from single stars and become members of a binary at least once in the simulations; : average number of exchanges per binary that hosts at least one MSBH (hereafter MSBH binary); : average MSBH binary lifetime; : fraction of MSBHs that undergo wind accretion at least once in the simulated time interval; : fraction of MSBHs that undergo wind accretion with an exchanged companion; : fraction of MSBHs that undergo RLO at least once in the simulated time interval; : fraction of MSBHs that undergo RLO with an exchanged companion; : average number of strong interactions per MSBH binary (considering only binaries that will undergo RLO).
4.2 Statistical properties of single and binary BHs
Table 2 provides a striking evidence of the importance of dynamics on the evolution of BHs in binaries, although the statistics is still quite low. The average number of BHs per cluster () is consistent with the expectations for a Kroupa IMF and does not appreciably depend on the metallicity.
The fraction of BHs that became members of a binary at least once in the simulated time interval is , regardless of the metallicity. This fraction is higher (by a factor of 2) than the fraction we expect if all the BH binaries come from primordial binaries. Furthermore, the fraction of BHs that formed from single stars and then became members of a binary as a consequence of a dynamical exchange is very high (). In particular, is very similar to , indicating that most of the BHs in binaries formed from single stars. In addition, the average number of exchanges per BH binary () is more than one during the simulated time interval. Therefore, dynamical exchanges dominate the life of binary BHs in the simulated young SCs.
The average lifetime of BH binaries is Myr in Table 2. In calculating , we assume that a binary survives even if one of its members exchanges, and dies only when it is completely ionized or when the simulation is stopped.
It is worth noting that BHBH and BHneutron star (NS) binaries (i.e. binaries where both the primary and the secondary member are compact objects) live longer than other BH binaries (on average), because they are not perturbed by the stellar evolution of the secondary and are sufficiently massive to avoid ionization. In particular, the average lifetime of BHBH binaries is 46, 37 and 32 Myr at , 0.1 and 1 Z, respectively. The average lifetime of BHNS binaries is similar: 49, 38 and 65 Myr at , 0.1 and 1 Z, respectively
The fraction of BHs that undergo wind accretion () indicates that there is about one wind-accreting system every two young clusters in 100 Myr. Interestingly, about half of these systems are consequences of dynamical exchanges (see in Table 2). This is true for all the considered metallicities (for Z the fraction of wind-accreting systems formed by exchanges is slightly lower than for lower metallicities, but this might be a fluctuation due to low statistics).
The fraction of BHs in RLO () is lower than . A large fraction of RLO systems are originated from a dynamical exchange (in Table 2, ). Most of the donor stars in RLO systems are post-MS stars ( per cent post-MS stars, versus per cent MS stars). A non-negligible fraction of such post-MS stars in RLO systems are LBV and WR stars ( per cent of the total donor stars in RLO systems). These results are in agreement with B06. We note that a non-negligible fraction of MS companions are blue straggler stars, i.e. stars rejuvenated by stellar mergers or by MT (e.g. Mapelli et al. 2004, 2006): they behave as MS stars, although they have mass higher than the turn-off mass.
In Table 2, is defined as the average number of strong resonant interactions per binary (i.e. interactions that lead to the formation of an unstable triple system and that change significantly the orbital period of the binary), calculated only for those binaries that will undergo RLO. These interactions are mostly three-body encounters and, in a few cases, four-body encounters (i.e. binary-binary encounters). The number of strong interactions is quite high (), confirming the importance of three-body encounters.
Finally, the last two columns of Table 2 show the number of MSBHs normalized to the total number of BHs (), and the number of RLO MSBHs normalized to the total number of RLO BHs (), respectively. and are similar, suggesting that the incidence of RLO systems among MSBHs is comparable to the incidence of RLO systems among low-mass BHs.
We stress that the main results presented in this Section, and especially the statistics of accreting systems, strongly depend on the assumed fraction of primordial binaries. In fact, primordial binaries have at least two important effects on accreting systems. Firstly, a number of primordial binary systems are expected to be sufficiently tight to start RLO at early times, as a consequence of stellar evolution. If there are no primordial binaries, the first RLO systems will appear later in the evolution of the SC, as an effect of dynamical interactions.
Secondly, primordial binaries (and especially hard primordial binaries) represent an initial reservoir of binding energy. This means that, if a BH becomes member of a hard primordial binary after an exchange, this binary might be sufficiently tight to start RLO immediately after the first exchange. On the contrary, non-primordial binaries that form from the encounter of three single stars are initially quite soft, and even the hardest among them become sufficiently tight to start RLO only after a lot of dynamical encounters (e.g. Hut et al. 1992). These arguments are discussed quantitatively (through a supplementary set of runs) in Appendix A.
We also note that the results presented in this Section strongly depend on the stellar evolution recipes and on the simplified model of failed SN adopted for our simulations. In Appendix B, we highlight the differences with respect to a scenario where the formation of MSBHs is strongly suppressed. In the forthcoming investigations, we will consider in more detail different models of stellar evolution.
4.3 Statistical properties of single and binary MSBHs
Table 3 reports the same quantities as shown in the first 11 columns of Table 2, but calculated only for the MSBHs (rather than for all the BHs). Table 3 shows that dynamical interactions are crucial for MSBHs, even more important than for low-mass stellar BHs.
First, the fraction of MSBHs that became members of a binary at least once in the simulated time interval (with respect to the total number of simulated MSBHs) is , i.e. a factor of higher than in the overall BH sample. At low metallicity ( and 0.1 Z), most of the MSBHs in binaries formed from single stars, as the fraction of MSBHs that formed from single stars and then became members of a binary as a consequence of a dynamical exchange is . Instead, for Z, no MSBHs in binaries form from single stars (). This is naturally explained by the fact that the only channel to form MSBHs at Z is the merger between two massive stars in a primordial binary (see Section 4.1).
Furthermore, the average number of exchanges per MSBH binary is , higher than for the overall BH sample (). The lifetime of MSBH binaries is slightly longer than that of the entire BH binary sample ( Myr). This is in fair agreement with B06, finding that a BH with M, born from runaway collapse, spends about 25 to 60 per cent of the simulated time interval with a companion.
Columns of Table 3 provide a key result to understand the formation of X-ray binaries powered by MSBHs. The comparison between columns and indicates that the large majority of wind-accreting MSBH binaries at low ( and 0.1 Z) were originated by a dynamical exchange (). The comparison between columns and indicates that all RLO MSBH binaries at low metallicity ( and 0.1 Z) were originated by a dynamical exchange (). Furthermore, all but one of the RLO MSBH binaries in our simulations host MSBHs that formed from single stars (this is not shown in Table 3, but we checked it from our simulations). The only MSBH (among RLO MSBH binaries) that did not form from a single star was originated from the merger of a primordial binary. Thus, even this MSBH formed as a single object and acquired a stellar companion after dynamical exchange.
The fraction of wind-accreting MSBH binaries (with respect to the total number of simulated MSBHs) is higher than the corresponding fraction for all BH binaries. The fraction of RLO MSBH binaries is similar to the corresponding fraction for all BH binaries, but the low statistics makes fluctuations very important (we find 18 RLO MSBH binaries over 300 simulations). Therefore, to study RLO MSBH binaries in detail, it will be essential to run a much larger grid of simulations (Mapelli et al., in preparation).
We find no wind-accretion MSBH binaries and no RLO MSBH binaries for Z, but this is likely an effect of the low statistics, as only five MSBHs form in our simulations at Z.
The results reported in Table 3 are in fair agreement with the findings by L10, based on population synthesis calculations. In particular, L10 find that no or very few RLO systems powered by MSBHs form from primordial binaries. This is consistent with our finding that no RLO MSBH binaries come from primordial binaries. On the other hand, the code used by L10 does not account for dynamical interactions. All RLO MSBH binaries found in our simulations were originated from dynamical interactions, in the sense that the MSBH became a member of the binary as a consequence of a dynamical exchange.
4.4 Orbital properties of accreting BHs
Figs. 4, 5 and 6 show some important characteristics of MT BHs. In these Figures, the filled circles are wind-accreting systems, while the crosses are RLO systems. Fig. 4 shows the mass of the BH versus the mass of the companion star (). At all considered metallicities, ranges from relatively low values ( M, mostly asymptotic giant branch stars) to very high values ( M or more) for both wind-accreting systems and RLO systems. The masses of accreting BHs span all the possible BH masses for a given metallicity (see Fig. 3, for comparison), even the higher masses. Not only wind-accreting systems but also RLO systems can host MSBHs ( M). This is an important difference with respect to L10 results, and it is mainly the consequence of dynamical interactions (which were not accounted by L10, see the discussion in Section 4.3). We note that for Z there is even one RLO system powered by a M MSBH. This MSBH is more massive than the maximum mass ( M) that can be reached at this metallicity through the adopted models of stellar evolution. In fact, such MSBH comes from the merger between a M BH and a M MS companion.
Fig. 5 shows the mass of the BH versus the orbital period . For all the considered metallicities, wind-accreting systems can form with periods as long as a few days, in agreement with previous studies (e.g. L10).
RLO systems have periods spanning from less than one day ( hours for one system at Z) up to yr. We stress that the periods of RLO systems shown in Fig. 5 are the values of the period at the first Roche lobe approach, as the code does not trace with accuracy the late stages of RLO and the time interval between snapshots is not sufficiently short to follow the evolution of all the systems. For example, if a system evolves into tidal instability and the companion is a MS star, the code assumes that the system is undergoing merger and removes the binary from calculation (Portegies Zwart & Verbunt 1996). Therefore, the plotted periods of RLO systems must be considered upper limits. We note that the companion star is an evolved star, with a very large radius (of the order of R), in most of the wide RLO systems ( yr).
Fig. 6 shows the mass of the BH versus the time elapsed since the beginning of the simulation. This plot gives information about the duty cycle (), defined as the lapse of time for which a binary is wind accreting and/or in RLO, divided by the total elapsed time in the simulation. Fig. 6 indicates that most systems are in RLO for less than one snapshot (corresponding to Myr). This is in fair agreement with the results by B06, which find that RLO systems powered by IMBHs are on average short-lived ( per cent of the simulated time interval). Wind accretion can last for a longer time (a few Myr).
Relatively low-mass BHs ( M) tend to start the RLO phase short after their formation. These systems start RLO as a consequence of stellar evolution of the secondary and/or as an effect of the first SN kick (see L10). Instead, more massive BHs start RLO at later times ( Myr after the beginning of the simulation). This means that the most massive BHs were single or in relatively wide binaries and can start RLO only as a consequence of the hardening of the binary by three-body encounters and/or of a dynamical exchange. Therefore, dynamical interactions are essential to allow MSBHs to power RLO X-ray binaries. This difference is also important to understand which X-ray sources are associated with low-mass BHs and which X-ray sources might be powered by MSBHs. In fact, from Fig. 6 we expect that the donor stars in X-ray binaries powered by MSBHs are on average older than those in X-ray binaries powered by low-mass BHs. In particular, most MSBHs enter the RLO phase with companions that are Myr old. Interestingly, most of the ULXs for which information about the stellar environment is available are associated with Myr old stellar populations (see e.g. Soria et al. 2005; Liu et al. 2007; Grisé et al. 2008; Swartz et al. 2009; Grisé et al. 2011; Voss et al. 2011).
Further hints about the importance of dynamics come from the variation of the orbital period. Fig. 7 shows the final period of a RLO binary (defined as the period at the beginning of RLO) versus the initial period of the same binary (defined as the period at the beginning of the simulation). For consistency, in Fig. 7 we show only primordial binaries that do not undergo dynamical exchanges before starting RLO. For most systems , as it was reasonable to expect. The shrinking is due to the joint effect of stellar evolution (e.g. a CE phase before the formation of the first BH forces the semi-major axis to shrink) and of three-body encounters (especially when the period changes by a factor of or more).
The large majority of systems have yr. Systems with an initial period below this threshold merge before the formation of the first BH. There are only two systems with yr, both for Z. These systems undergo at least one strong three-body encounter before the formation of the first BH. The three-body encounters widen the semi-major axis of these binaries, allowing them to avoid merger and to survive till the formation of the first BH
4.5 BH ejections
: fraction of BHs (including MSBHs) that are ejected from the SC (i.e. that have distance from the centre of the cluster ). This and all the other fractions reported in this Table are normalized to the total number of simulated BHs; : fraction of BHs that are ejected by natal kick (i.e. after the SN explosion); : fraction of MSBHs that are ejected from the SC; : fraction of BHs that are ejected from the SC together with their companion star; : fraction of BHs that are ejected from the SC together with their companion star and start RLO after the ejection.
The possibility that a BH is ejected from the parent SC is relevant for various astrophysical issues. For example, bright HMXBs and ULXs are often close to young SCs and star forming regions, but displaced by pc with respect to their centre (e.g. Zezas et al. 2002; Kaaret et al. 2004; Berghea 2009; Swartz et al. 2009; Swartz 2010; Rangelov, Prestwich & Chandar 2011; Voss et al. 2011; Poutanen et al. 2012). This fact has been generally interpreted as the indication that bright X-ray sources are powered by runaway binaries, that is, by binaries that were ejected from the parent cluster because of a natal kick (e.g. Sepinsky, Kalogera & Belczynski 2005; Zuo & Li 2010) or because of a close encounter (e.g. Kaaret et al. 2004; Berghea 2009; Mapelli et al. 2011b).
In our simulations, BHs can be ejected both through SN explosion (natal kick) and through three-body encounters. Table 4 shows that per cent of simulated BHs are ejected from the SC, almost independently of the metallicity (we classify a BH as ejected when its distance from the centre of mass of the SC is ). About half of the ejections are consequences of the natal kick, while the remaining half is due to three body encounters. SN explosions can also unbind a binary system. We estimate that primordial binaries per SC are ionized by a SN explosion leading to the formation of a BH (we do not include NSs in this estimate).
At and 0.1 Z, the fraction of ejected MSBHs (with respect to the total number of simulated BHs) is . Since MSBHs are per cent of all the BHs at these metallicities, this means that about per cent of all the simulated MSBHs are ejected as a consequence of three-body encounters (we recall that in our simulations we assume that MSBHs receive no natal kick). This percentage is moderately lower, but still in agreement with Mapelli et al. (2011b), who find that per cent of all MSBHs are ejected as a consequence of three-body encounters. The difference between these two estimates can be explained by the fact that Mapelli et al. (2011b) generate the MSBHs already in the initial conditions (rather than letting them form later, through stellar evolution), and assume that all MSBHs are members of primordial binaries. Thus, the available time for MSBHs to undergo thee-body encounters is longer in the simulations by Mapelli et al. (2011b). This assumption increases the number of interactions involving MSBH binaries, especially during the first stage of core collapse (which occurs as early as Myr in our simulated SCs). Furthermore, the simulations in Mapelli et al. (2011b) do not include stellar evolution.
Table 4 also shows that about 10 per cent of all the simulated BHs are ejected together with their companion star, regardless of the metallicity. About one tenth of these ejected binaries enter RLO after leaving the cluster (). Since the total fraction of RLO systems is (see Table 2), this means that about 20 per cent of all RLO systems enter the RLO phase after being ejected from the SC (and in most of the cases the RLO is triggered by a dynamical interaction). Only one of these ejected RLO systems is powered by a MSBH. Thus, most of the simulated RLO systems are inside the parent SC. This result is apparently at odds with observations (e.g. Zezas et al. 2002; Kaaret et al. 2004; Berghea 2009; Swartz et al. 2009; Swartz 2010), which indicate that a significant fraction of bright X-ray sources are offset from the parent cluster by more than 10 pc. On the other hand, we stress that the simulations presented in this paper do not include the tidal field of the host galaxy: as our simulated SCs have a relatively small initial mass, a large fraction of binaries can be stripped from the SC by external tidal forces.
Finally, in the current paper we have assumed that the MSBHs are born without natal kick. This assumption depends crucially on the model of direct collapse and failed SN. Thus, the fraction of ejected MSBHs reported in this paper must be regarded as a lower limit.
4.6 A schematic interpretation of MSBH behaviour
In the previous sections, we highlighted the differences between the behaviour of MSBHs and that of low-mass stellar BHs, discussing the results of the body simulations. In this section, we show how the results of the body simulations can be intuitively understood in light of the interplay between stellar evolution and some of the basic properties of three-body encounters. The three most relevant aspects from the physics of three-body encounters are the following.
(i) The probability of a single star (or stellar remnant) with mass to exchange into a binary () is higher if or , where and are the masses of the two components of the binary (e.g. Hills & Fullerton 1980; Sigurdsson & Phinney 1993). In particular, is very close to one, if the binary is hard (so that it cannot be easily ionized) and if or (Hills & Fullerton 1980). After the first Myr of the cluster lifetime (when the turn-off mass goes below M), single MSBHs are among the most massive objects in the SC and are likely to replace the lowest-mass members of primordial binaries.
(ii) The cross section for three-body interactions is larger for more massive binaries (see e.g. Sigurdsson & Phinney 1993; Davies, Benz & Hills 1994; Davies 2002; Miller & Hamilton 2002). After the first Myr of the cluster lifetime, MSBH binaries become significantly more massive than the other binaries in the SC. Thus, their rate of three-body encounters is higher than that of low-mass binaries.
(iii) Hard binaries tend to harden as a consequence of three-body encounters, i.e. to reduce their semi-major axis (e.g. Heggie 1975; Heggie & Hut 1993; Davies 1995; Quinlan 1996; Merritt 2001). The decrease of can start RLO into a BH binary. A MSBH binary in our simulated SCs is hard if A.U. (corresponding to a period yr). Therefore, most of the MSBH binaries in our simulations are hard.
Combining stellar evolution with the above notions from the theory of three-body encounters, we can summarize the evolution of the simulated BHs as shown in Fig. 8. In the following, we discuss the main differences between MSBHs and low-mass BHs, relatively to the flow chart in Fig. 8.
Case A: the BH forms from a single star. If the BH is a MSBH, then the exchange probability is very high already after Myr since the SC formation, and the MSBH becomes soon a binary member via dynamical exchange. The resulting MSBH binary undergoes efficiently three-body encounters and hardens rapidly.
If the BH has relatively low mass, its is low for most of the SC life: the low-mass BH can remain single for a much longer time, or even for the entire lifetime of the cluster. Furthermore, low-mass BHs are more easily ejected out of the SC.
Case B: the BH forms from a primordial binary. Depending on the initial semi-major axis of the binary and on the stellar evolution, the binary can merge or avoid merger before the formation of the BH. If the binary merges before the formation of the BH (Case B1), then the BH forms as a single object and its evolution is the same as in the Case A.
If the binary does not merge before the formation of the BH (Case B2), then there are at least three possibilities.
Case B2.1: The SN kick of the primary unbinds the binary. In this case, the BH is substantially a single object and behaves as in Case A. In our simulations, this case can occur only for low-mass BHs, as MSBHs are assumed to form without natal kick.
Case B2.2: The binary remains bound and the stellar evolution of the secondary is faster than three-body encounters. This occurs especially when the binary is very close and/or the mass of the secondary is similar to the mass of the primary. For binaries where the radius of the secondary star is close to the Roche lobe at the time of formation of the BH, the binary enters a RLO phase because of the evolution of the secondary. As already shown by L10, MSBHs can hardly form in such close binaries, because the first CE phase leads generally to the formation of small BHs. Therefore, this case is more frequent for low-mass BHs.
Case B2.3: The binary remains bound and three-body encounters are more efficient than the stellar evolution of the secondary. This occurs especially for relatively wide binaries (where stellar evolution is not sufficient to drive RLO), and/or for low-mass secondary stars, which evolve much more slowly than the primary. In this case, the evolution is very different depending on the mass of the BH and depending on whether the binary is hard or soft. In the following, we mention only three of the possible cases (the most relevant for our simulations).
(i) If the binary is soft (unlikely if a MSBH is member of the binary), three-body encounters are expected to ionize it.
(ii) If the BH is a MSBH and the binary is relatively hard, then three-body encounters harden the binary. Exchanges can occur, but are unlikely to remove the MSBH from the binary, as it is more massive than most single stars. The MSBH binary survives and may undergo MT.
(iii) If the BH has a low mass (lower than the companion mass) and the binary is relatively hard, then three-body encounters harden the binary, but exchanges can remove the BH from the binary. If it is expelled by an exchange, the BH becomes single again.
We stress that the flow chart in Fig. 8 is schematic and somehow simplistic, as it neglects some further effects that can take place (e.g. the SN explosion of the companion star, the formation of a binary composed of two compact objects, the definitive ejection of the BH from the SC). Despite this, Fig. 8 allows to understand why exchanges are so important to enhance the formation and evolution of MSBH binaries (see Table 3). In fact, MSBHs evolve predominantly from Case A (i.e. they form from single stars) to Case B2.3, where exchanges and three-body encounters dominate the evolution of a hard and massive MSBH binary.
It clarifies also why most low-mass BHs start RLO immediately after their formation (see Fig. 6), whereas MSBHs enter RLO much later than their formation ( Myr after the beginning of the simulation). In fact, it is unlikely for MSBHs to evolve through Case B2.2 and to start RLO because of the stellar evolution of the secondary star in a primordial binary. MSBHs evolve mainly through Case A and Case B2.3: they can efficiently acquire (new) companions because of dynamical exchanges, and their binaries harden because of three-body encounters. On the other hand, dynamical interactions occur on a longer timescale (a few tens Myr) than massive star evolution.
Instead, low-mass BHs evolve mainly through Case A, Case B2.2 and Case B2.3. In Case A, low-mass BHs remain likely single objects, as the probability of an exchange into a binary is low. Case B2.2 can lead to MT on the timescale of massive star evolution. Case B2.3 can result in the ejection of the low-mass BH from the binary, as a consequence of dynamical exchanges. Therefore, most of the RLO systems powered by low-mass BHs evolve as described in Case B2.2, and only a few of them through Case B2.3.
4.7 RLO systems and Eddington luminosity
To describe the evolution of mass accretion and the emission properties of RLO binaries is beyond the scope of this paper. On the other hand, from the simulations described in the previous sections, we can extract some basic hints about the evolution of RLO binaries.
The code we use limits the accretion rate so that the luminosity of the binary cannot exceed the Eddington limit for BHs. We prefer not to change this assumption and not to introduce arbitrary rules for super-Eddington accretion, as the dynamical code is not sufficiently accurate to distinguish between different accretion models. In a forthcoming paper, we will re-simulate the MT systems obtained from this study, without accounting for three-body encounters but with a more accurate recipe for the MT rate (e.g. Patruno & Zampieri 2008).
Fig. 9 shows the Eddington luminosity () for the simulated RLO systems. There are important differences between the solar-metallicity simulations and the two sub-solar environments, although the statistics is quite low. No RLO systems formed at have erg s (corresponding to a BH mass M). RLO systems formed at and reach erg s and erg s, respectively (corresponding to BH mass and 58 M, respectively).
The mass spectrum of BHs born from the collapse of massive stars is very uncertain. Recent theoretical studies (e.g. Mapelli et al. 2009; Zampieri & Roberts 2009; B10) indicate that the mass of the BH depends on the metallicity of the progenitor star and that metal-poor massive stars can produce MSBHs with a mass as high as M (for ). Observational results suggest that low-metallicity environments host more massive BHs than the MW (e.g. Prestwich et al. 2007) and highlight a possible anti-correlation between the population of bright X-ray sources (especially ULXs) and the metallicity of the host galaxy (e.g. Mapelli et al. 2010, 2011a).
However, it is still unclear whether theoretically predicted MSBHs are efficient in powering X-ray sources. Studies based on population synthesis codes (e.g. L10), including metallicity effects but neglecting dynamical interactions, find that MSBHs can hardly power bright X-ray sources, because of various effects connected with binary evolution.
In this paper, we present preliminary simulations run with a body plus stellar and binary evolution code, based on the public version of Starlab. These simulations include metallicity-dependent stellar evolution, binary evolution and an accurate treatment of three-body encounters. We simulate dense young SCs with three different metallicities: , and .
Our simulations show that the dynamics is essential to map the population of BH binaries and that of X-ray binaries in young SCs. Firstly, the rate of BH binaries formed through dynamical exchanges is very high for all the considered metallicities. Secondly, the hardening of the binary due to three-body encounters is an important mechanism to induce RLO in binaries hosting BHs.
The role of dynamics is particularly important for MSBHs. In SCs with and 0.1 Z, about 75 per cent of simulated MSBHs form from single stars and become members of binaries through dynamical exchanges in the first 100 Myr of the SC life. This is a factor of more efficient than in the case of low-mass ( M) stellar BHs. We show that the vast majority of MSBHs in RLO binaries originated from single stars and went through a dynamical exchange (Table 3 and Section 4.3). This is consistent with the qualitative predictions from the theory of three-body encounters combined with stellar evolution (Section 4.6).
In fact, the higher is the mass of the BH involved, the more important the effects of dynamics on X-ray binaries, as the rate of three-body encounters scales approximately with the total mass of the binary (e.g. Sigurdsson & Phinney 1993). Furthermore, the probability for a single BH to become a binary member through a dynamical exchange is higher if the mass of the single BH is higher than the mass of one of the two components of the binary (e.g. Hills & Fullerton 1980). Therefore, MSBHs are even more influenced by dynamics than their low-mass analogues.
Our results agree with the basic conclusion by L10: it is very hard that MSBHs power RLO X-ray binaries, if they form in primordial binaries that evolve unperturbed. On the other hand, we find that the situation is completely different for MSBHs that form in dynamically active environments (such as dense young SCs). In these environments, MSBHs efficiently power RLO X-ray binaries, as a consequence of dynamical exchanges and three-body encounters, which dramatically alter the orbital properties of primordial binaries. Our results are also in fair agreement with those by B06, although there are some important differences, as B06 study IMBHs born from the runaway collapse in denser SCs, do not account for metallicity-dependent stellar evolution, and have a much simplified dynamical treatment. Both our paper and B06 find that BHs with mass M spend a large fraction of the simulated time with a companion star and that these massive BHs can power RLO binaries as a consequence of three-body encounters.
Our preliminary results show that the study of MSBHs is very promising, but there is a lot of work still to be done. First, the sample of MSBHs obtained from our simulations is statistically small (112 and 118 MSBHs for the simulations with and Z, respectively). Although our main results are physically well motivated and in agreement with previous studies about three-body encounters (e.g., Hills & Fullerton 1980; Sigurdsson & Phinney 1993) and population synthesis models (e.g. L10), a much larger sample is required to study in detail the properties of MSBHs in RLO binaries.
From the theoretical point of view, different models of stellar evolution and especially of mass loss by stellar winds need to be investigated with the same approach. Different environments (not only dense young SCs, but also globular clusters and open clusters) need to be studied with larger statistics. Furthermore, in this paper we assume a primordial binary fraction . Higher binary fractions are not unrealistic (e.g. Delgado-Donate et al. 2004; Ivanova et al. 2005; Sollima et al. 2010) and must be considered in future simulations. The code described in this paper traces accurately the dynamics of close encounters, but adopts simplified recipes for stellar and binary evolution, as well as for MT. The MT systems individuated by our simulations will be re-simulated with more accurate binary evolution codes, to trace how they evolve, which are their expected emission features and whether they can explain a fraction of the ULXs.
A fundamental question is where we can search for MSBHs. MW globular clusters can reach metallicities as low as . Thus, they may host MSBHs, single or with low-mass companions. Metallicities of the order of are not infrequent in nearby irregular galaxies (e.g. IC 10), which are good candidates to further investigate the metallicityBH mass connection.
We thank the anonymous referee for her/his critical reading of the manuscript, which improved it. We made use of the public software package Starlab (version 4.4.4) and of the SAPPORO library (Gaburov, Harfst & Portegies Zwart 2009) to run Starlab on graphics processing units (GPUs). We acknowledge all the developers of Starlab, and especially its primary authors: Piet Hut, Steve McMillan, Jun Makino, and Simon Portegies Zwart. We thank the authors of SAPPORO, and in particular E. Gaburov, S. Harfst and S. Portegies Zwart. We also thank Paola Marigo, Alessia Gualandris, Antonella Vallenari, Rosanna Sordo, Laura Greggio, Mauro Barbieri and Monica Colpi for stimulating discussions. We acknowledge the CINECA Award N. HP10CXB7O8 and HP10C894X7, 2011 for the availability of high performance computing resources and support. MM, LZ and ER acknowledge financial support from INAF through grant PRIN-2011-1. LZ acknowledges financial support from ASI/INAF grant no. I/009/10/0. AB acknowledges financial support from MIUR through grant PRIN-2009-1. MM thanks the Aspen Center for Physics, where part of this work was done.
Appendix A The importance of primordial binaries
An initial binary fraction was adopted for all the simulations presented in the main text. In this Appendix, we show the results of a set of 100 supplementary runs without primordial binaries (). While these simulations are quite unrealistic, as primordial binaries are expected to exist in SCs, they are useful to illustrate the importance of primordial binaries. All the other properties of the SCs are the same as in the main text. In order to better see the effects on MSBHs, we maximize their formation probability by taking Z.
The quantities shown in this Table are the same as defined in columns and of Table 2.
Table A1 shows some of the most relevant properties of the BH population in the runs with and Z. Some interesting considerations can be drawn from the comparison of Table A1 with Table 2. First, the number of BHs per cluster is about 10 per cent lower than in the case with . This is because a number of BHs form from the merger of relatively small stars (13-25 M) in primordial binaries, if . This is evident from Fig. 10, where the initial mass spectrum of BHs in the runs with and with are compared. No BHs with mass M form in the case with , whereas all the BHs with mass M formed in the runs with are remnants of merged progenitors. We stress that the mass of BHs born from merged progenitors depends strongly on the assumptions about mass loss during the merger.
The fraction of BHs that become members of a binary at least once in the simulation () in Table A1 is almost the same as in Table 2. On the other hand, by construction in the case with . Thus, the number of exchanges that produce BH binaries is higher if . Similarly, the average number of exchanges per binary () is almost double in the case of than in the case of . These results indicate that exchanges involving BHs are more numerous if .
Furthermore, the properties of accreting BH binaries are considerably affected by the assumed fraction of primordial binaries. Wind-accreting systems and especially RLO systems are strongly suppressed if . This can be easily understood in light of the formation mechanism of binaries. If primordial binaries are allowed to form, then a number of primordial binary systems are expected to be sufficiently close to start RLO at early times, as a consequence of stellar evolution. In the absence of primordial binaries, binaries can form only through the dynamical interaction of three single stars (hereafter, three-body capture) and through a tidal interaction between two single stars (hereafter, tidal capture). Three-body capture binaries are generally wide and eccentric (Hut et al. 1992), whereas tidal-capture binaries are very hard and generally merge (e.g., Portegies Zwart et al. 1997). These two mechanisms are unlikely for most of the SC life, but are enhanced during core collapse (Spitzer 1987). In our simulations, we do not include any recipes for tidal-capture binaries. Thus we have only three-body capture binaries. Only a small fraction () of BH binaries born from three-body capture survives ionization and becomes hard enough (through three-body encounters) to start RLO within 100 Myr since the beginning of the simulation. The right-hand panel of Fig. 11 confirms this consideration, by showing that all RLO systems switch on at relatively late times.
The quantities shown in this Table are the same as defined in columns of Table 3.
The last two columns of Table A1 suggest another important consideration: if , MSBHs are much more efficient than low-mass BHs in powering RLO systems. In fact, the fraction of RLO systems powered by MSBHs () is about twice as high as the fraction of MSBHs with respect to all BHs (). The bias towards high BH masses in RLO systems appears even more evident if we look at Fig. 11: all the BHs powering RLO systems have masses M. This is a strong confirmation that the more a BH is massive, the more its chances of entering a binary by dynamical exchange are high.
Finally, Table A2 shows some of the most relevant properties of MSBHs in the runs with . From the comparison with Table 3, it is apparent that more MSBHs can form in the absence of primordial binaries ( and 1.18 if and 0.1, respectively). If , the mass loss following the CE phase in a primordial binary tends to produce smaller BH masses. Finally, the fraction of RLO systems powered by MSBHs is a factor of 10 smaller if than if .
In summary, all the aspects discussed in this Appendix indicate that primordial binaries have a crucial importance for the population of X-ray binaries in SCs.
Appendix B Different stellar evolution and SN explosion recipes
The quantities shown in this Table are the same as defined in columns and of Table 2.
The quantities shown in this Table are the same as defined in columns of Table 3.
In this paper, we follow the stellar evolution recipes described in Section 3.1 and we assume that stars with M directly collapse to MSBHs. Different stellar evolution recipes and different assumptions for the end of massive star life have important effects on the presented results. A complete comparison between different stellar evolution recipes will be done in the next papers of the series. In this Appendix, we just show a simplified, rather extreme case, in which the formation of MSBHs is strongly suppressed. In particular, we assume that the maximum BH mass is M, even for Z. Thus, in these runs the MSBHs can have only M.
We performed 50 runs of SC evolution assuming M, Z and leaving all the other properties of the SCs as described in Section 3. The main results for the BH population and for X-ray binaries powered by BHs are summarized in Tables B1 and B2. Table B1 shows that , , and are very similar to the case with Z described in Table 2. The most significant difference between Table B1 and Table 2 is the value of . The fraction of MSBHs that power RLO systems is much smaller in Table B1 than in Table 2 (0.04 versus 0.24). This confirms that MSBHs in the mass range M are much more efficient in powering RLO systems than BHs with M. The same conclusion can be derived from the comparison between Table 3 (for Z) and Table B2, where the statistics for MSBHs is shown. The MSBHs in Table B2 (whose mass is M by construction) are members of binaries and power RLO systems less often than the MSBHs in Table 3.
- We name ‘final mass’, , of a star the mass bound to the star immediately before the collapse.
- We stress that our results cannot be easily generalized to more massive SCs and in particular to globular clusters, as the masses and the relevant timescales are too different (in globular clusters the half-mass relaxation time is orders of magnitude longer than the timescale for the evolution of massive stars, while in intermediate-mass SCs these two timescales are comparable). The study of such larger systems requires dedicated simulations.
- In the following, we call MSBH binary (BH binary) a binary hosting at least one MSBH (BH).
- Recent studies show that the IMF might be top-heavy in dense low-metallicity regions (Marks et al. 2012). Thus, our choice of a Kroupa IMF for all the considered metallicities is quite conservative, as it reduces the differences in the number of BHs among different metallicities.
- In the simulated SCs, binaries are soft if erg, where is the gravitational constant. For M this corresponds to R
- The formation of MSBHs from the merger of two massive stars at depends strongly on the treatment of the star (which is a blue straggler star) formed from the merger, and in particular on the mass lost by this star before the collapse. This issue is delicate and deserves further study.
- The long average lifetime ( Myr) of BHNS binaries at Z is explained by statistical fluctuations: only four BHNS binaries form at Z in our simulations.
- According to Heggie’s law (Heggie 1975), hard binaries tend to harden as a consequence of three-body encounters. This is true in a statistical sense. Single interactions can widen even hard binaries.
- Abolmasov P. K., Swartz D. A., Fabrika S., Ghosh K. K., Sholukhova O., Tennant A. F., 2007, ApJ, 668, 124
- Akiyama Sh., Wheeler J. C., Meier D. L., Lichtenstadt I., 2003, ApJ, 584, 954
- Ardeljan N. V., Bisnovatyi-Kogan G. S., Moiseenko S. G., 2005, MNRAS, 359, 333
- Belczynski K., Kalogera V., Zezas A., Fabbiano G., 2004a, ApJ, 601, L147
- Belczynski K., Sadowski A., Rasio F. A., 2004b, ApJ, 611, 1068
- Belczynski K., Kalogera V., Rasio F. A., Taam R. E., Zezas A., Bulik T., Maccarone Th. J., Ivanova N., 2008, ApJS, 174, 223
- Belczynski K., Bulik T., Fryer C. L., Ruiter A., Valsecchi F., Vink J. S., Hurley J. R., 2010, ApJ, 714, 1217 (B10)
- Belkus H., Van Bever J., Vanbeveren D. 2007, ApJ, 659, 1576
- Berghea C., 2009, Ultraluminous X-ray sources and their Environment, PhD thesis, http://students.cua.edu/79berghea/papers/dissertation.pdf
- Bica E., Dutra C. M., Soares J., Barbuy B., 2003, A&A, 404, 223
- Blecha L., Ivanova N., Kalogera V., Belczynski K., Fregeau J., Rasio F., 2006, ApJ, 642, 427 (B06)
- Bonnell I., Kroupa P., 1998, to appear in The Orion Complex Revisited, eds. M.J.McCaughrean, A.Burkert, ASP Conference Series, arXiv:astro-ph/9802306
- Borissova J., et al., 2011, A&A, 532A, 131
- Carpano S., Pollock A. M. T., Prestwich A., Crowther P., Wilms J., Yungelson L., Ehle M., 2007, A&A, 466, L17
- Charles P. A., Coe M. J., 2006, In: Compact stellar X-ray sources. Edited by Walter Lewin & Michiel van der Klis. Cambridge Astrophysics Series, No. 39. Cambridge, UK: Cambridge University Press, 215
- Chatterjee S., Fregeau J. M., Umbreit S., Rasio F. A., 2010, ApJ, 719, 915
- Colpi M., Mapelli M., Possenti A., 2003, ApJ, 599, 1260
- Crowther P. A., Carpano S., Hadfield L. J., Pollock A. M. T., 2007, A&A, 469, L31
- Crowther P. A., Barnard R., Carpano S., Clark J. S., Dhillon V. S., Pollock A. M. T., 2010, MNRAS, 403, L41
- Davies M. B., Benz W., Hills J. G., 1994, ApJ, 424, 870
- Davies M. B., 1995, MNRAS, 276, 887
- Davies M. B., 2002, in ‘Omega Centauri, A Unique Window into Astrophysics’. ASP Conference Proceedings, Vol. 265. Edited by Floor van Leeuwen, Joanne D. Hughes, and Giampaolo Piotto. San Francisco, Astronomical Society of the Pacific, p.215
- Delgado-Donate E. J., Clarke C. J., Bate M. R., Hodgkin S. T., 2004, MNRAS, 351, 617
- Dessart L., O’Connor E., Ott Ch., 2012, ApJ, 754, 76
- Dias W. S., Alessi B. S., Moitinho A., Lépine J. R. D., 2002, A&A, 389, 871
- Downing J. M. B., Benacquista M. J., Giersz M., Spurzem R., 2010, MNRAS, 407, 1946
- Downing J. M. B., Benacquista M. J., Giersz M., Spurzem R., 2011, MNRAS, 416, 133
- Dray L. M., 2006, MNRAS, 370, 2079
- Duquennoy A., & Mayor M., 1991, A&A, 248, 485
- Fregeau J. M., Gürkan M. A., Joshi K. J., Rasio F. A., 2003, ApJ, 593, 772
- Fregeau J. M., & Rasio F. A., 2007, ApJ, 658, 1047
- Fryer Ch. L., 1999, ApJ, 522, 413
- Fryer Ch. L., & Heger A., 2000, ApJ, 541, 1033
- Fryer Ch. L., Kalogera V., 2001, ApJ, 554, 548
- Fryer Ch. L., & Warren M. S., 2004, ApJ, 601, 391
- Fryer Ch. L., Belczynski K., Wiktorowicz G., Dominik M., Kalogera V., Holz D. E., 2012, ApJ, 749, 91
- Gaburov E., Harfst S., Portegies Zwart S., 2009, New Astronomy, 14, 630
- Garnett D. R., 1990, ApJ, 363, 142
- Gieles M., Portegies Zwart S. F., 2011, MNRAS, 410, L6
- Giersz M., 1998, MNRAS, 298, 1239
- Goad M. R., Roberts T. P., Knigge C., Lira P., 2002, MNRAS, 335, L67
- Grisé F., Pakull M. W., Soria R., Motch C., Smith I. A., Ryder S. D., Böttcher M., 2008, A&A, 486, 151
- Grisé F., Kaaret P., Pakull M. W., Motch C., 2011, ApJ, 734, 23
- Grisé F., Kaaret P., Corbel S., Feng H., Cseh D., Tao L., 2012, ApJ, 745, 123
- Gürkan M. A., Freitag M., Rasio F. A., 2004, ApJ, 604, 632
- Hamann W.-R., Koesterke L., 1998, A&A, 335, 1003
- Harlaftis E. T., Greiner J., 2004, A&A, 414L, 13
- Heger A., Fryer C.L., Woosley S.E., Langer N., Hartmann D.H., 2003a, ApJ, 591, 288
- Heger A., Woosley S.E., Fryer C.L., Langer N., 2003b, in From Twilight to Highlight: The Physics of Supernovae: Proceedings of the ESO/MPA/MPE Workshop Held at Garching, Germany, 29-31 Juli 2002, ESO ASTROPHYSICS SYMPOSIA. ISBN 3-540-00483-1. Edited by W. Hillebrandt and B. Leibundgut. Springer-Verlag, p. 3
- Heggie D. C., 1975, MNRAS, 173, 729
- Heggie D. C., Hut P., 1993, ApJS, 85, 347
- Hénon M. H., 1971a, Ap&SS, 13, 284
- Hénon M. H., 1971b, Ap&SS, 14, 151
- Hillenbrand L. A., Hartmann L. W., 1998, ApJ, 492, 540
- Hills J. G., Fullerton L. W., 1980, AJ, 85, 1281
- Humphreys R. M., Davidson K., 1994, Astronomical Society of the Pacific, Publications, 106, 1025
- Hurley J. R., Pols O. R., Tout C. A., 2000, MNRAS, 315, 543
- Hurley J. R., Tout C. A., Pols O. R., 2002, MNRAS, 329, 897
- Hut P., McMillan S., Goodman J., Mateo M., Phinney E. S., Pryor C., Richer H. B., Verbunt F., Weinberg M., 1992, Astronomical Society of the Pacific, Publications, 104, 981
- Ivanova N., Belczynski K., Fregeau J. M., Rasio F., 2005, MNRAS, 358, 572
- Ivanova N., Heinke C. O., Rasio F. A., Taam R. E., Belczynski K., Fregeau J., 2006, MNRAS, 372, 1043
- Ivanova N., Heinke C. O., Rasio F. A., Belczynski K., Fregeau J., 2008, MNRAS, 386, 553
- Joshi K. J., Rasio F. A., Portegies Zwart S., 2000, ApJ, 540, 969
- Joshi K. J., Nave C. P., Rasio F. A., 2001, ApJ, 550, 691
- Kaaret P., Alonso-Herrero A., Gallagher J. S., Fabbiano G., Zezas A., Rieke M. J., 2004, MNRAS, 348L, 28
- King I. R., 1966, AJ, 71, 64
- Kraicheva Z. T., Popova E. I., Tutukov A. V., Iungelson L. R., 1978, Astronomicheskii Zhurnal, 55, 1176
- Kroupa P., 2001, MNRAS, 322, 231
- Kudritzki R. P., Pauldrach A., Puls J., 1987, A&A, 173, 293
- Kudritzki R.-P., Puls J., 2000, ARA&A, 38, 613
- Kudritzki R. P., 2002, ApJ, 577, 389
- Kuhn M. A., Baddeley A., Feigelson E. D., Getman K. V., Broos P. S., Townsley L. K., Povich M. S., Naylor T., King R. R., Busk H. A., Luhman K. L., 2012, arXiv:1208.3492, to appear in The Labyrinth of Star Formation, (eds.) D. Stamatellos, S. Goodwin, and D. Ward-Thompson, Springer, in press
- Lada Ch. J., Lada E. A., 2003, ARA&A, 41, 57
- Lee C.-H., Brown G. E., Wijers R. A. M. J., 2002, ApJ, 575, 996
- Leitherer C., Robert C., Drissen L., 1992, ApJ, 401, 596
- Linden T., Kalogera V., Sepinsky J. F., Prestwich A., Zezas A., Gallagher J., 2010, ApJ, 725, 1984 (L10)
- Liu J.-F., Bregman J., Miller J., Kaaret Ph., 2007, ApJ, 661, 165
- Liu J.-F., Bregman J., Seitzer P., 2004, ApJ, 602, 249
- Livio M., Warner B., 1984, The Observatory, 104, 152
- Madhusudhan N., Justham S., Nelson L., Paxton B., Pfahl E., Podsiadlowski Ph., Rappaport S., 2006, ApJ, 640, 918
- Madhusudhan N., Rappaport S., Podsiadlowski Ph., Nelson L., 2008, ApJ, 688, 1235
- Maeder A., Meynet G., 2010, New Astronomy, 54, 32
- Maeder A., Meynet G., 2012, Reviews of Modern Physics, 84, 25
- Maeder A., 1992, A&A, 264, 105 (M92)
- Mapelli M., Sigurdsson S., Colpi M., Ferraro F. R., Possenti A., Rood R. T., Sills A., Beccari G., 2004, ApJ, 605L, 29
- Mapelli M., Sigurdsson S., Ferraro F. R., Colpi M., Possenti A., Lanzoni B., 2006, MNRAS, 373, 361
- Mapelli M., Colpi M., Zampieri L., 2009, MNRAS, 395L, 71
- Mapelli M., Ripamonti E., Zampieri L., Colpi M., Bressan A., 2010, MNRAS, 408, 234
- Mapelli M., Ripamonti E., Zampieri L., Colpi M., 2011a, AN, 332, 414
- Mapelli M., Ripamonti E., Zampieri L., Colpi M., 2011b, MNRAS, 416, 1756
- Marks M., Kroupa P., Dabringhausen J., Pawlowski M. S., 2012, MNRAS, 422, 2246
- Mercer E. P., et al., 2005, ApJ, 635, 560
- Merritt D., 2001, ApJ, 556, 245
- Miller M. C., Hamilton D. P., 2002, MNRAS, 330, 232
- Narayan R., McClintock J. E., 2005, ApJ, 623, 1017
- Nelemans G., Yungelson L. R., Portegies Zwart S. F., Verbunt F., 2001, A&A, 365, 491
- Orosz J. A., 2003, in van der Hucht K. A., Herrero A., Esteban, C., eds, Proc. IAU Symp. 212, A Massive Star Odyssey: From Main Sequence to Supernova, p. 365
- Orosz J. A., et al., 2007, Nature, 449, 872
- Özel F., Psaltis D., Narayan R., McClintock J. E., 2010, ApJ, 725, 1918
- Patruno A., Zampieri L., 2008, MNRAS, 386, 543
- Pattabiraman Bh., Umbreit S., Liao W.-K., Choudhary A., Kalogera V., Memik G., Rasio F. A., 2012, submitted to ApJ Supplements, arXiv:1206.5878
- Pauldrach A. W. A., Vanbeveren D., Hoffmann T. L. 2012,A&A, 538, 75
- Pfalzner S., 2009, A&A, 498, L37
- Pilyugin L. S., Vílchez J. M., Contini T., 2004, A&A, 425, 849
- Podsiadlowski Ph., Rappaport S., Pfahl E. D., 2002, ApJ, 565, 1107
- Podsiadlowski Ph., Rappaport S., Han Z., 2003, MNRAS, 341, 385
- Portegies Zwart S. F., Verbunt F., 1996, A&A, 309, 179
- Portegies Zwart S. F., Verbunt F., Ergma E., 1997, A&A, 321, 207
- Portegies Zwart S. F., McMillan S. L. W., Hut P., Makino J., 2001, MNRAS, 321, 199
- Portegies Zwart S. F., McMillan S. L. W., 2002, ApJ, 576, 899
- Portegies Zwart S. F., 2004, arXiv:astro-ph/0406550, Lecture note to appear in “Joint Evolution of Black Holes and Galaxies” of the Series in High Energy Physics, Cosmology and Gravitation. IOP Publishing, Bristol and Philadelphia, 2005, eds M. Colpi, V.Gorini, F.Haardt and U.Moschella
- Portegies Zwart S. F., McMillan S. L. W., Makino J., 2007, MNRAS, 374, 95
- Portegies Zwart S. F., McMillan S. L. W., Gieles M., 2010, ARA&A, 48, 431
- Poutanen J., Fabrika S., Valeev A. F., Sholukhova O., Greiner J., 2012, MNRAS, submitted, arXiv:1210.1210
- Prestwich A. H., Kilgard R., Crowther P. A., Carpano S., Pollock A. M. T., Zezas A., Saar S. H., Roberts T. P., Ward M. J., 2007, ApJ, 669, L21
- Quinlan G. D., 1996, NewAstronomy, 1, 35
- Ramsey C. J., Williams R. M. Gruendl R. A., Chen C.-H. R., Chu Y.-H., Wang Q. D., 2006, ApJ, 641, 241
- Rangelov B., Prestwich A. H., Chandar R., 2011, ApJ, 741, 86
- Rappaport S. A., Podsiadlowski Ph., Pfahl E., 2005, MNRAS, 356, 401
- Richards E. E., Lang C. C., Trombley Ch., Figer D. F., 2012, AJ, 144, 89
- Sepinsky J., Kalogera V., Belczynski K., 2005, ApJ, 621L, 37
- Sigurdsson S., Phinney E. S., 1993, ApJ, 415, 631
- Sigurdsson S., Phinney E. S., 1995, ApJS, 99, 609
- Silverman J. M., Filippenko A. V., 2008, ApJ, 678, L17
- Smartt S. J., 2009, ARA&A, 47, 63
- Sollima A., Carballo-Bello J. A., Beccari G., Ferraro F. R., Fusi Pecci F., Lanzoni B., 2010, MNRAS, 401, 577
- Soria R., Cropper M., Pakull M., Mushotzky R., Wu K., 2005, MNRAS, 356, 12
- Soria R., 2006, in Populations of High Energy Sources in Galaxies, Edited by Meurs E. J. A.; Fabbiano G.. Cambridge: Cambridge University Press, 473
- Soria R., Kuntz K. D., Winkler P. F., Blair W. P., Long K. S., Plucinsky P. P., Whitmore B. C., 2012, ApJ, 750, 152
- Spitzer L., 1987 Dynamical evolution of globular clusters, Princeton University Press
- Stodolkiewicz J. S., 1982, Acta Astronomica, 32, 63
- Stodolkiewicz J. S., 1986, Acta Astronomica, 36, 19
- Swartz D. A., Tennant A. F., Soria R., 2009, ApJ, 703, 159
- Swartz D. A., 2010, Galaxy Wars: Stellar Populations and Star Formation in Interacting Galaxies ASP Conference Series Vol. 423. Edited by Smith B., Bastian N., Higdon S. J. U., and Higdon J. L.. San Francisco: Astronomical Society of the Pacific, p.277
- Tao L., Feng H., Grisé F., Kaaret Ph., 2011, ApJ, 737, 81
- Terashima Y., Inoue H., Wilson A. S., 2006, ApJ, 645, 264
- Thompson T. A., Quataert E., Burrows A., 2005, ApJ, 620, 861
- Vanbeveren D., 2009, New Astronomy Reviews, 53, 27
- van der Hucht K. A., 1991, in ‘Wolf-Rayet Stars and Interrelations with Other Massive Stars in Galaxies’: Proceedings of the 143rd Symposium of the International Astronomical Union. Edited by Karel A. van der Hucht and Bambang Hidayat. International Astronomical Union. Symposium no. 143, Kluwer Academic Publishers, Dordrecht, p.19
- Vink J. S., de Koter A., Lamers H. J. G. L. M. , 2001, A&A, 369, 574
- Vink J. S., de Koter A., 2005, A&A, 442, 587
- Voss R., Nielsen M. T. B., Nelemans G., Fraser M., Smartt S. J., 2011, MNRAS, 418, L124
- Wheeler J. C., Akiyama S., Williams P. T., 2005, Ap&SS, 298, 3
- Woosley S. E., Heger A., Weaver T. A., 2002, Reviews of Modern Physics, 74, 1015
- Woosley S. E., Bloom J. S., 2006, ARA&A, 44, 507
- Zampieri L., Roberts T., 2009, MNRAS, 400, 677
- Zezas A., Fabbiano G., Rots A. H., Murray S. S., 2002, ApJ, 577, 710
- Zuo Z.-Y., Li X.-D., 2010, MNRAS, 405, 2768