Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://www.arcetri.astro.it/~palla/webpage/fp_SMBH.pdf Äàòà èçìåíåíèÿ: Sun Mar 2 19:03:56 2014 Äàòà èíäåêñèðîâàíèÿ: Sun Apr 10 02:05:57 2016 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: massive stars

A& A 5 5 8 , A5 9 ( 2 0 1 3 ) DOI: 10.1051/0004-6361/201321949
c ESO 2013

Astronomy & Astrophysics

Massive black hole factories: Supermassive and quasi-star formation in primordial halos
Dominik R. G. Schleicher1 , Francesco Palla2 , Andrea Ferrara3,4 , Daniele Galli2 , and Muhammad Latif1
1 2 3 4

Institut fÝr Astrophysik, Georg-August-UniversitÄt GÆttingen, Friedrich-Hund-Platz 1, 37077 GÆttingen, Germany e-mail: [dschleic;mlatif]@astro.physik.uni-goettingen.de INAF ­ Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, 50125 Firenze, Italy e-mail: [palla;galli]@arcetri.astro.it Scuola Normale Superiore, Piazza dei Cavalieri 7, 56126 Pisa, Italy e-mail: andrea.ferrara@sns.it Kavli IPMU (WPI), the University of Tokyo, 5-1-5 Kashiwanoha, 277-8583 Kashiwa, Japan

Received 23 May 2013 / Accepted 14 August 2013
ABSTRACT

Context. Supermassive stars and quasi-stars (massive stars with a central black hole) are both considered as potential progenitors for the formation of supermassive black holes. They are expected to form from rapidly accreting protostars in massive primordial halos. Aims. We explore how long rapidly accreting protostars remain on the Hayashi track, implying large protostellar radii and weak accretion luminosity feedback. We assess the potential role of energy production in the nuclear core, and determine what regulates the evolution of such protostars into quasi-stars or supermassive stars. Methods. We followed the contraction of characteristic mass shells in rapidly accreting protostars, and inferred the timescales for them to reach nuclear densities. We compared the characteristic timescales for nuclear burning with those for which the extended protostellar envelope can be maintained. Results. We find that the extended envelope can be maintained up to protostellar masses of 3.6 â 108 m3 M , where m denotes the accretion rate in solar masses per year. We expect the nuclear core to exhaust its hydrogen content in 7 â 106 yr. If accretion rates m 0.14 can still be maintained at this point, a black hole may form within the accreting envelope, leading to a quasi-star. Alternatively, the accreting object will gravitationally contract to become a main-sequence supermassive star. Conclusions. Due to the limited gas reservoir in typical 107 M dark matter halos, the accretion rate onto the central object may drop at late times, implying the formation of supermassive stars as the typical outcome of direct collapse. However, if high accretion rates are maintained, a quasi-star with an interior black hole may form.
Key words. cosmology: theory ­ dark ages, reionization, first stars ­ stars: protostars ­ stars: formation

1. Introduction
Supermassive black holes (SMBHs) with more than 109 M have been observed at z > 6(Fan et al. 2001, 2003, 2004, 2006), and recently even at z = 7.085 (Mortlock et al. 2011). While theoretical considerations indicate a variety of potential pathways, including direct collapse to a black hole, collapse to a supermassive star, collapse of stellar clusters or clusters of black holes (Rees 1984; Haiman 2006; Volonteri & Bellovary 2012), substantial efforts have been made to assess the validity of such scenarios.
1.1. Black hole formation scenarios

As the most straightforward possibility, one could imagine these SMBHs to originate from remnants of the first stars, which were proposed to be very massive with 100-1000 M (Abel et al. 2002; Bromm & Larson 2004). Stellar evolution calculations by Heger & Woosley (2002) indeed show that black holes may form for stellar masses between 10-30 M , or 100-300 M , thus yielding the seeds for further growth to SMBHs.

New calculations following the evolution of primordial star formation beyond the first core, however, show that fragmentation can occur efficiently (Stacy et al. 2010, 2012; Turk et al. 2009; Clark et al. 2011; Greif et al. 2011; Smith et al. 2011). While the effect of the accretion luminosity on the protostellar core seems to be minor (Smith et al. 2011, 2012), UV feedback during the main sequence phase appears to set the upper limit to the stellar mass to 50-100 M (Hosokawa et al. 2011, 2012b; Susa 2013). Even if very massive stars could form, the resulting HII region is likely to inhibit accretion for at ´ least 108 yrs (Milosavljevic et al. 2009a,b). An additional problem has been identified by Whalen & Fryer (2012), who showed that black holes resulting from such stars would be born with strong dynamical kicks, implying that they should readily be expelled from the halo. On the other hand, dense stellar clusters may form at high redshift in the presence of trace amounts of dust (Schneider et al. 2003, 2008, 2012a,b; Omukai et al. 2005, 2012; Clark et al. 2008; Dopcke et al. 2011, 2013; Klessen et al. 2012). The relativistic instability of such clusters has been explored in a semi-analytic framework by Devecchi & Volonteri (2009) and A5 9 , p a g e 1 o f 9

Article published by EDP Sciences


A&A 558, A59 (2013)

Devecchi et al. (2010, 2012), finding characteristic black hole masses of up to 3000 M . The latter is already substantially higher than the masses of the first stars, though it may still be difficult to obtain masses as high as 109 M at z = 7. In the presence of a significant black hole spin, which will be rapidly obtained in the presence of efficient accretion, one typically expects an additional growth by about four orders of magnitude, thus requiring seed masses on the order of 105 M (Shapiro 2005). The formation of such massive seeds has been considered in the context of the direct collapse model, where gas in a 107 ­108 M halo is expected to collapse without fragmentation into a single central object (Koushiappas et al. 2004; Begelman et al. 2006; Lodato & Natarajan 2006; Spaans & Silk 2006; Volonteri et al. 2008; Begelman & Shlosman 2009). In order to avoid fragmentation, molecular hydrogen needs to be efficiently dissociated, requiring a strong ambient UV field (Omukai 2001; Bromm & Loeb 2003; Shang et al. 2010; Schleicher et al. 2010; Latif et al. 2011; Van Borm & Spaans 2013). Such radiation backgrounds may have been provided by nearby starburst galaxies, which may indeed occur frequently enough to explain the observed abundance of SMBHs at z 6 (Dijkstra et al. 2008; Agarwal et al. 2012). As an alternative, Inayoshi & Omukai (2012) suggested that black hole formation may occur in the presence of strong shocks, leading to collisional dissociation of molecular hydrogen. Using numerical simulations, Wise et al. (2008) have modeled the gravitational collapse of massive primordial halos cooling via atomic hydrogen, reporting an isothermal density profile and angular momentum transport via bar-like instabilities during the formation of the first core. Regan & Haehnelt (2009) presented the first study following the evolution beyond the first core, and modeling the formation of self-gravitating disks on parsec-scales. While these studies employed a typical resolution of 16 cells per Jeans length, it was recently demonstrated that turbulent structures can only be resolved with a resolution of at least 32 cells per Jeans length, preferably more (Sur et al. 2010; Federrath et al. 2011; Turk et al. 2012). In a set of highresolution simulations following the formation of the first core, Latif et al. (2013a,b) found extended turbulent structures in the center of massive primordial halos, but no signs of simple disks or bar-like instabilities. Disks do, however, form at later stages of the evolution, with masses of 1000 M on scales of 30 AU, and characteristic accretion rates of 1 M yr-1 (Latif et al. 2013c). As a result, very massive central objects may indeed form in rather short cosmic times. The dynamics at early times have also been explored by Choi et al. (2013), finding a somewhat different result including the formation of toroidal structures on scales of several parsec. It is, however, not fully clear to which extent the latter is a result of the simplified initial conditions employed in their calculation. To provide more quantitative predictions regarding the conditions where black holes may form, Prieto et al. (2013) explored the correlations of the baryon spin with the spin of the dark matter, showing that the resulting correlation is, however, weak, and that the baryon properties cannot be naively extrapolated from the dark matter. In a similar study, de Souza et al. (2013) report that the dynamical properties of the dark matter halo may leave an imprint on the resulting stellar mass.
1.2. Quasi-stars vs. supermassive stars

While the accumulation of high masses seems feasible from a hydrodynamical point of view, the resulting object is still considerably more uncertain. As a first step, we expect the formation of a massive protostar, as the gas becomes optically thick at high
A59, page 2 of 9

densities (Omukai & Palla 2001, 2003; Hosokawa et al. 2012a). It is, however, unclear how these objects are going to evolve, and whether they will form a supermassive star or a quasi-star. In the present context, we define a supermassive main sequence star as a conventional star with very high mass ( M = 103 -106 M ) (see also Shapiro & Teukolsky 1986), which predominantly generates its energy via nuclear burning. A quasistar refers to an object with mass similar to that of a supermassive star, but whose central core has collapsed into a black hole (Begelman et al. 2006; Begelman 2010). The progenitors of both types of objects are rapidly accreting high-mass protostars, and we aim to explore in this paper under which conditions they evolve into a quasi-star or a supermassive main-sequence star. To explore the evolution of such protostars, Hosokawa et al. (2012a) recently followed the evolution of rapidly accreting protostars, showing that they expand as cool supergiants, thus inhibiting the feedback from accretion luminosity and potentially allowing accretion to proceed for very long times. Hosokawa et al. also found that nuclear burning starts at stellar masses of about 50 M , and is maintained until the end of their calculation at 103 M . Therefore, their model indicates no transition toward a quasi-star at least during the evolutionary phase considered. It is, however, important to assess how long this phase o f e ffi c i e n t a c c r e t i o n c a n b e m a i n t a i n e d , a n d u n d e r wh i c h c o n ditions a transition to a supermassive main-sequence star can be expected, which can potentially influence the accretion flow via radiative feedback (Omukai & Inutsuka 2002; Johnson et al. 2011, 2012). While the majority of such supermassive stars directly collapses into a black hole (Fryer & Heger 2011), a small mass window exists around 55 000 M where violent supernova explosions of up to 1055 erg may occur (Johnson et al. 2013). Following a different approach, Begelman et al. (2006) and Begelman (2010) have proposed the formation of a black hole in the interior of rapidly accreting objects, leading to quasi-stars as the progenitors of SMBHs. Begelman et al. (2006) argued that for quasi-stars with 104 ­105 M , the typical accretion timescale is considerably shorter than the timescale for nuclear burning, and as a result, the latter may not be able to stop the collapse of the central core. A potential advantage of such a configuration is that the accretion onto the central object is not limited by the Eddington accretion rate of the black hole, but by the Eddington accretion rate of the more massive quasistar (Begelman et al. 2008). Begelman (2010) considers objects with more than 106 M , which first evolve toward the main sequence as a supermassive star, but then form a black hole after an extended phase of nuclear burning. Ball et al. (2011) followed the evolution of black holes in a quasi-star employing the Cambridge STARS stellar evolution package (Eggleton 1971; Pols et al. 1995), finding that the black hole may efficiently accrete 10% of the stellar mass before hydrostatic equilibrium breaks down. They also report that the results are sensitive to the boundary condition in the interior. The contraction of an embedded isothermal core in a stellar envelope was further explored by Ball et al. (2012) in the framework of the SchÆnberg-Chandrasekhar limit. While the models for the quasi-stars typically employed a Thomson-scattering opacity, Hosokawa et al. (2012a) reported that H- rather than Thomson scattering dominates the opacity in the protostellar atmosphere. They find luminosities close to the Eddington luminosity, LE = 4G M mp c = 3.8 â 104 L T M , M (1)

dd


D. R. G. Schleicher et al.: Massive black hole factories

where G is the gravitational constant, mp the proton mass, c the speed of light, T the Thomson scattering cross section and M the mass of the star. The temperatures in the atmosphere are rather cool, 5000 K, as the protostars are on the Hayashi track. As a result, the characteristic radii evolve as Rin = 2.6 â 103 R M M
1/2

2. The interplay of mass accretion and Kelvin-Helmholtz contraction
In the following, we will sketch why the protostar may maintain large envelopes even if the Kelvin-Helmholtz timescale becomes shorter than the timescale for accretion. For this purpose, we consider mass shells of enclosed mass M , located at radii R( M, t). A mass shell M forms at the time M/ M , with an initial radius as given by Eq. (2). The radii of these mass shells evolve on their Kelvin-Helmholtz timescale, given as tKH ( M, R) = GM 2 · R( M, t)LEdd (7)

·

(2)

The protostars are thus considerably more extended and show an increasing stellar radius as a function of mass, while the models of Begelman et al. (2006) indicate a constant radius as a function of stellar mass. Begelman (2010) consider a mass-dependent stellar radius, which is, however, still more compact due to the assumption of Thomson scattering opacities. The difference between these scenarios may, however, turn out to be highly relevant, as the stellar radius regulates the temperature on the surface, and thus the strength of protostellar feedback. In fact, the behavior reported by Hosokawa et al. (2012a) is well-known (Omukai & Palla 2001, 2003) in cases where the accretion timescale, tacc = M , M (3)

As a result, we have dR R( M, t) =- · dt tKH ( M, R) The equation is the mass shell time tin ( M ) = integration yield (8)

integrated between the initial radius Rin when forms, given in (2), corresponding to the M/ M , to the current radius R at time t. The s (9)

is much smaller than the Kelvin-Helmholtz timescale, tKH GM 2 · = RL (4) behavwh e r e via acto conximum

4 mp c 1 1 = (t - tin ( M ))· + R Rin T M

Inserting the initial radius as given in Eq. (2) and expressing the result in astrophysical units, we obtain 1000 R R = 1000 R 2.6 â 103 ( M/M )1/2 1.16 â 10-1 yr-1 (t - tin ( M )) + · M/M

However, Hosokawa et al. (2012a) report that this ior extends into the regime where tKH < tacc , i.e. Kelvin-Helmholtz contraction is faster than mass growth cretion. The physical mechanism which allows this phase tinue will be discussed in detail in Sect. 2. Taking the ma luminosity of the star Lmax 0.6 L M M
11/2

(10)

R R

-1/2

,

(5)

which assumes Kramers' opacity T -3.5 (Hayashi et al. 1962) and an initial radius as given in Eq. (2), they show that both timescales are equal at protostellar masses of Meq = 14.9 M m 0.01
0.26

,

(6)

where we parametrized the accretion rate as M m M yr-1 . On the other hand, the accreting star remains on the Hayashi track up to stellar masses of 1000 M , the highest mass reached in the calculation of Hosokawa et al. (2012a). The extended envelope locked at T eff 5000 K allows efficient accretion with only moderate feedback for a longer period. In this paper, we aim to assess how long this efficient accretion phase can be maintained beyond M 1000 M without strong feedback from the protostar. In addition, we will discuss the potential impact of nuclear burning for the evolution of the accreting objects, and the conditions under which a quasi-star as opposed to a supermassive star may form. For this purpose, we explore the interplay of mass accretion and Kelvin-Helmholtz contraction in Sect. 2. We calculate the impact of nuclear burning in Sect. 3. A final discussion of our results is provided in Sect. 4. The impact of additional processes such as deuterium shell burning and hydrogen burning via the pp-chain is assessed in the Appendix.

The first term on the right-hand side thus dominates right after the formation of a mass shell and determines its initial radius, while the second term describes its evolution due to Kelvin-Helmholtz contraction. It is remarkable that at late times, the second term dominates and the evolution of the radii appears to be independent of their initial position. This can be understood, as our expression for the Kelvin-Helmholtz timescale, Eq. (7), scales with the inverse radius of the mass shell. As a result, the evolution slows down during the contraction, implying that the Kelvin-Helmholtz timescale may increase significantly in the interior. An example for the resulting structure in a 105 M star is given in Fig. 1, including both the regime where the first and the second term of the equation dominate. We now aim to estimate when a given mass shell reaches the densities of nuclear burning, which we take as nuc 1 g cm-3 (Hosokawa et al. 2012a). We obtain a radius of nuclear burning, which is given as R
nuc



3M 4nuc

1/3

1 .2 R

M M

1/3

·

(11)

We note that the above provides only a rough estimate of the radius of the nuclear core. If the latter were more compact, the timescale for a given mass scale to reach nuclear burning would become even longer, implying that more mass could be accreted during that time. A comparison with Eq. (2) shows that this radius always remains smaller than the radius of the star by at least a factor of 1000, and in fact increases more gradually with stellar mass. A given mass shell M will only reach nuclear densities on timescales much longer than the initial Kelvin-Helmholtz timescale tKH , which increases as R-1 . We can thus determine the
A59, page 3 of 9


A&A 558, A59 (2013)

during the contraction. This behavior will only change when the timescale of accretion becomes comparable to the timescale on which the outer mass shell is able to reach nuclear densities. As we show here, the latter is significantly larger than the Kelvin-Helmholtz timescale of the star, and scales as R-1 . The nuc transition occurs when a critical mass scale of 3.6 â 108 m3 M is reached, implying a transition toward a supermassive main sequence star. We note that the transition may even occur earlier if m decreases with time.

3. The nuclear evolution
In this section, we will assess the potential role of nuclear burning. For this purpose, we show that the transition to the CNO-cycle rapidly occurs and regulates the production of helium in the nuclear core. We will then calculate under which conditions the fuel in the nuclear core will be exhausted while the protostar maintains its extended envelope. The latter implies the potential formation of a black hole in the interior, thus a transition into a quasi-star.
3.1. Impor tance of the CNO-cycle

Fig. 1. Enclosed mass as a function of radius in the protostar calculated from Eq. (10), assuming a total stellar mass of 105 M and mass accretion rates m = 0.1, 1., 10. The current time t is calculated as t = ( M/M )/m. For a given mass shell, higher accretion rates imply larger radii, as the star had less time to contract.

time when nuclear densities are reached by equating (11) with the second term in (10). As a result, we obtain t = t - tin ( M ) = 710 yrs M M
2/3

·

(12)

During the time interval t, the star will accrete an additional mass M = M t. Considering a given mass shell, we can calculate the ratio M M M t = 710 m M M M
-1/3

While the pp-cycle will dominate in the very beginning (see appendix), we show here that the transition to the CNO-cycle rapidly occurs as a result of helium burning via the triple- process. Following Padmanabhan (2000), the energy production rate via triple- is given as
3

=

5.1 â 108 2 Y 3 e 3 T9

-4.4027/T

9

erg g-1 s-1 ,

(16)

,

(13)

which describes how much additional mass is accreted before the shell reaches nuclear densities. The ratio drops below 1 only for protostellar masses of M 3.6 â 108 m3 M . (14)

While for typical accretion rates of m 10-3 , this happens al ready at mass scales of order unity, this transition occurs only at substantially higher masses above 1000 M for m > 10-2 . While a given shell evolves toward nuclear densities, the amount of additional mass that is accreted is thus substantially higher than the mass in that shell. From Eq. (12), we can further derive the mass in the nuclear core in the limit that t tin ( M ). We obtain Mnuc = t 710 yr
3/2 3/2 M = t710 M ,

(15)

where Y 0.25 denotes the mass fraction of helium and T 9 denotes the gas temperature in units of 109 K. It is straightforward to check that the energy production rate changes by several orders of magnitude for temperature changes on the order of 10%. As shown by Hosokawa et al. (2012a), the highest temperatures with T 9 0.15 are expected in the center of the nuclear core, which has an enhanced nuclear density of 10 g cm-3 . With a volume filling factor fill 0.01-0.1, the contribution of that region to helium burning is significant due to the steep temperature dependence. For the following estimates, we will assume that efficient mixing occurs throughout the nuclear core, implying that the heavy elements produced in the central region will be available throughout the core. On the other hand, if mixing is inefficient, nuclear burning via the CNO cycle will dominate even more in the central region due to its higher metallicity. Noting that the energy released by a single triple- reaction corresponds to 1.166 â 10-5 erg, the number of heavy elements produced per unit time is given as N3 = Mnuc · 1.166 â 10-5 erg
3 fill

where we defined t710 = t/(710 yr). As we neglected the term M/ M in Eq. (12), we note that the derivation here implicitly assumes a sufficient mass supply to the protostar in order to feed the core. Clearly, this approximation will break down at the mass scale derived in (14). At that point, the protostellar mass may become approximately constant, and we expect a transition toward supermassive main sequence stars. From this model, we already obtain a set of relevant conclusions concerning the evolution of the protostar. In particular, if the accretion rate is high, a substantial amount of additional matter is accreted before a given shell reaches nuclear densities. As a result, the evolution of the protostar is dominated not by the interior shells, but by the additional matter which is accreted
A59, page 4 of 9

(17)

The CNO mass in the core thus evolves as and an integration yields MCNO = 7.1 â 10
-17 5/2 fill t710

MCNO = 12mp â N3

(18)

M.

(19)

Assuming efficient mixing, the resulting metallicity in the nuclear core is then MCNO = 6.3 â 10-9 t710 fill . (20) Z= Mnuc


D. R. G. Schleicher et al.: Massive black hole factories

We recall that the energy production rate in the CNO cycle is given as (Padmanabhan 2000)
CNO

=

4.4 â 1025 XZ T
2/3 9

e

-15.228/T

1/3 9

erg g

-1 -1

s.

(21)

Since the energy production rate depends sensitively on the temperature, we expect its contribution in the innermost core to be dominant. For the temperature T 9 = 0.15, it is straightforward to show that energy production via the CNO cycle becomes comparable to the pp cycle for a metallicity of Zc = 2 â 10-12 . A comparison with Eq. (20) shows that the CNO cycle thus dominates after a short time of tCNO = 2.25 where we introduced fill = 0.1 fill,-1 . We note here that this timescale is probably not accurate, as our assumptions in Sect. 2 (in particular concerning the Eddington luminosity) only become valid at later times. Nevertheless, this result implies that the CNO cycle can be expected to be relevant early on, and thus needs to be considered. As the production of elements becomes increasingly efficient at late times, we expect that the expression (20) will yield a reasonable estimate in the period of interest. With the above assumptions, the CNO luminosity due to nuclear burning is given as LCNO =
CNO,Zc -1 fill,-1

yr,

(22)
Fig. 2. Ratio of nuclear luminosity provided via the CNO-cycle over the Eddington luminosity of the star as a function of protostellar mass for different mass accretion rates. The CNO luminosity is calculated here as the minimum of the estimate in Eq. (24) and the Eddington luminosity of the nuclear core in Eq. (25). At high protostellar masses, one can clearly recognize the transition to the regime where nuclear burning is limited by the Eddington luminosity of the core. Higher accretion rates imply that less time was available for contraction, thus reducing the nuclear luminosities.

Z Mnuc Zc

An integration yields MHe
,CNO

fill

,

(23)

= 1.0 â 10

-3 2 7/2 fill t710

M.

(29)

where CNO,Zc denotes the energy production rate via the CNO cycle evaluated at the critical metallicity Zc . Inserting (20), we obtain The Eddington luminosity of the nuclear core is given as LE
dd,core

A comparison with the mass of the core indicates that they become comparable at tHe
,CNO

LCNO = 1.3 â 10

5 2 5/2 fill t710

L.

(24)

= 2.25 â 10

5 -1 fill,-1

yr,

(30)

The ratio of these luminosities to the Eddington luminosity of the star is given in Fig. 2. Even here, the luminosity produced by nuclear burning approaches the Eddington luminosity of the star only around stellar masses of 108 M , implying no relevant impact on the stellar evolution during the earlier stages. A comparison of these expressions yields the timescale For fill,-1 1, as indicated by Hosokawa et al. (2012a) for m = 0.1, the CNO luminosity would exceed the Eddington luminosity of the core after 2 â 104 yr, shortly after the end of their simulations. However, when the Eddington luminosity is reached, one would expect an expansion of the nuclear core, implying lower densities and an adiabatically decreasing temperature. Due to this thermostat, it is likely that the core will adjust to a state maintaining its Eddington luminosity. We will now demonstrate that the nuclear core is not converted to helium before the critical timescale tc after which the core radiates at its Eddington luminosity. We consider the production of helium by the CNO process, given as NHe
,CNO

= 3.8 â 104 mt710 L . 3/2

(25)

implying that a pure helium core will only be formed once the nuclear luminosity is equal to the Eddington luminosity of the core.
3.2. Formation of a helium core

We will now calculate the helium production rate under the assumption that the formation of a helium core occurs once the nuclear luminosity equals the Eddington luminosity of the core. Thus, NHe = LEdd,core · 4.3 â 10-5 erg (31)

tc = 2.1 â 10

4 -2 fill,-1

yr.

(26)

The helium mass production rate is thus given as MHe = 4mp â NHe and an integration yields MHe = 1.05 â 10
-4 5/2 t710

(32)

M.

(33)

=

and the helium mass production rate MHe
,CNO

LCNO , 4.3 â 10-5 erg

(27)

= 12mp â NHe

,CNO

.

(28)

Equating Eq. (33) with Eq. (15), we find that the nuclear core becomes a pure helium core after t f 6.8 â 106 yr. We note that at this point, the mass of the nuclear core is given as M f 9.4 â 105 M . In order to illustrate how the latter translates into total stellar masses, the ratio of helium mass to total mass in the core is displayed in Fig. 3 for different accretion rates, showing that the hydrogen fuel is indeed exhausted for stellar masses of 3 â 105 -3 â 107 M depending on the accretion rate. After that point, the burning of heavier elements will occur, but only last for a short period. As the mass of the core is
A59, page 5 of 9


A&A 558, A59 (2013) Table 1. Pathways to the formation of main-sequence supermassive stars and quasi-stars depending on the critical mass scale Mcr ( ) = m 3.6 â 108 m3 M and the accretion rate m in solar masses per year. M < Mcr ( ) m protostar with extended envelope, cool atmospheres M > Mcr ( ), m m 0.14 transition to main-sequence supermassive star M > Mcr ( ), m m 0.14 collapse of the nuclear core, quasi-star

Fig. 3. Fraction of helium in the nuclear core, corresponding amount of exhausted fuel for nuclear burning, as a function of stellar mass for different accretion rates. Higher accretion rates less time for protostellar contraction. As a result, nuclear densit reached later, implying less time for nuclear burning.

to the protoimply ies are

considerably higher than the Chandrasekhar mass scale and the Tolman-Oppenheimer-Volkoff mass scale, the core may collapse and form a black hole within the star, a so-called quasi-star, as proposed by Begelman et al. (2006) and Begelman (2010). The interior black hole may then continue accreting from the stellar envelope, while the accretion rate of the quasi-star will only be limited by the Eddington rate corresponding to the total mass of the configuration. The derivation of t f and M f , however, implicitly assumes that enough mass has been accreted in order to build up the core, requiring that M
f

t f M.

(34)

To form a quasi-star, we thus require an accretion rate of M M f /t f , i.e. m 0.14. (35)

For lower accretion rates, we expect that the extended envelope of the star cannot be maintained at late times, implying that the star will contract and form a supermassive main sequence star. During this period, both the stellar mass and the mass in the nuclear core will be approximately constant, implying a characteristic hydrogen burning timescale of 2 â 106 yr. For stellar masses above 300 M , these are expected to collapse to massive black holes (Fryer et al. 2001; Fryer & Heger 2011).

4. Discussion and outlook
To investigate the interplay of accretion and Kelvin-Helmholtzcontraction of the protostar, we have presented an analytical model following the evolution of the mass shells from the extended envelope down to nuclear densities. Based on this model, we are able to explain the fundamental result by Hosokawa et al. (2012a) that the accreting protostar keeps its extended envelope beyond the adiabatic accretion phase. The reason is that the timescale on which a given mass shell may reach nuclear densities is much longer than the Kelvin-Helmholtz timescale of the star, as the protostellar radius is considerably larger than the radius of the nuclear core. Instead of the Kelvin-Helmholtz timescale evaluated at the protostellar radius, the relevant timescale is thus the Kelvin-Helmholtz timescale
A59, page 6 of 9

at the expected core radius after contraction, providing an appropriate estimate of the evolutionary timescale. As a consequence, rapidly accreting protostars are expected to keep their extended envelopes until stellar masses of 3.6 â 108 m3 M . At this point, the accretion timescale becomes longer than the Kelvin-Helmholtz timescale in the stellar interior, implying that the star will subsequently contract toward the main sequence. Then, we expect it to follow the typical evolution of a supermassive star as sketched by Shapiro & Teukolsky (1986). We further assessed the impact of nuclear burning for the evolution of these stars. While deuterium burning will never be relevant during the protostellar evolution, hydrogen burning via the pp-cycle could potentially affect the stellar evolution after 1015 years, which is, however, practically irrelevant. Due to the efficient metal production via the triple- process, a transition to the CNO cycle is expected early on. The energy production by the CNO cycle is considerably more efficient, implying that the core can radiate at its Eddington luminosity. As a result, a helium core forms after 7 â 106 yr. The expected mass of the core is then 106 M , which provides an upper limit on the initial mass of the interior black hole. Therefore, if the available gas reservoir is sufficient to maintain accretion rates m 0.14 until this point, the central core may collapse during the protostellar evolution phase, and the object becomes a quasi-star. Alternatively, the protostar will evolve toward the main sequence to become a supermassive star. Such an object is then expected to collapse into a SMBH at the end of its lifetime (Fryer & Heger 2011). A summary of our results regarding the evolution of rapidly accreting protostars is given in Table 1. Our results have been derived based on a comparison of the mass accumulated by accretion with the characteristic mass scale where we expect a pure helium core. This mass scale is based on the reasonable assumption that the luminosity of the core is given by its Eddington luminosity, in agreement with the numerical results by Hosokawa et al. (2012a). Uncertainties are, however, present, as our calculation for the evolution of the core mass assumes that the luminosity of each mass shell is given by the Eddington luminosity. Following Hosokawa et al. (2012a), the latter is a good but not precise approximation. Further corrections can be expected in particular when approaching nuclear densities, when feedback from nuclear burning becomes significant. We further note that in the regime where deviations from the hydrostatic equilibrium are significant, as during the collapse of the nuclear core or during the accretion onto the protostellar envelope, it is conceivable that the luminosity released from the gravitational potential may locally exceed the Eddington luminosity if the accretion rate is high enough. Due to the momentum of the flow, we, however, expect that such accretion rates could still overcome the radiative feedback and effectively trap the photons in the optically thick material. We further note that Dotan & Shaviv (2012) propose a scenario where porosity induced by instabilities may lead to super-Eddington luminosities.


D. R. G. Schleicher et al.: Massive black hole factories

The results obtained differ from the scenario proposed by Begelman et al. (2006) and Begelman (2010), as we make use of new results on the protostellar evolution that were obtained in the recent study by Hosokawa et al. (2012a). In particular, we note that Begelman et al. (2006) still assumed a timeindependent radius of the quasi-stars, and both Begelman et al. (2006) and Begelman (2010) adopted Thomson scattering as the only opacity source. However, Hosokawa et al. (2012a) have shown that H- may significantly contribute to the opacity in the atmosphere. By considering the timescale required for a given mass shell to reach nuclear densities, we further have shown here that the star may remain on the Hayashi track for a much longer period than expected before. The latter gives rise to considerable differences in the stellar evolution, and allows for significant a m o u n t s o f a c c r e t i o n wh i l e f e e d b a c k i s w e a k . Additional effects could be introduced as a result of rotation, which was not included in this model. We thus expect our results to provide an order of magnitude estimate for the critical accretion r ate, which m ay be de termined more accurately employing stellar evolution calculations, but also a more realistic model for the time-dependent accretion rate. We also note here that Hosokawa et al. (2012a) reported a slight dependence of their results on the employed boundary conditions for the protostar. While their main calculations adopted shock boundary conditions, they also explored the effect of lower-entropy accretion provided by photospheric boundary conditions. In that case, a slightly higher accretion rate of 0.3 M yr-1 appears to be required to maintain the extended envelopes, but the overall evolution remains very similar. For comparison, we note that recent simulations by Latif et al. (2013c) reported accretion rates of 1 M yr-1 in halos with 107 M . Adopting a constant accretion rate, the transition to a supermassive main sequence star would occur at a stellar mass of 3.6 â 108 M . As the nuclear fuel is, however, exhausted earlier, a central black hole may form before this transition, giving rise to an evolution as sketched by Ball et al. (2011). Due to the limited gas reservoir, and as here the halo mass is in fact comparable to the mass scale of the corresponding quasi-star, it seems likely that the accretion rate will substantially decrease at late times, implying that the transition toward the supermassive star should occur at an earlier stage which is then determined by the time evolution of the accretion rate, and that supermassive stars with 104 -105 M might be the most generic outcome of the collapse. More massive objects appear to be possible at least in principle, if a larger gas reservoir is available in more massive dark matter halos. In order to assess the potential influence of rotation, we estimated the amount of rotational support in the central 103 M clumps reported by Latif et al. (2013c), which varied between 5-20% in different simulations. With densities of 10-10 g cm-3 , these are still orders of magnitudes below the characteristic densities within the protostar, such that we cannot yet draw strong conclusions regarding the final amount of rotational energy. However, we note that Stacy et al. (2012) reported a significant amount of rotation in primordial protostars based on numerical simulations of Greif et al. (2012), and a quite similar case can be expected here. As recently shown by Reisswig et al. (2013), the latter may have a substantial impact on the collapse of supermassive stars, implying the potential formation of a black hole binary and a subsequent merger, accompanied with efficient emission of gravitational waves. The latter provides a potential pathway of probing black hole formation scenarios

with LISA1 . In a narrow mass range around 55 000 M , one may further expect the occurence of highly energetic supernovae with energies up to 1044 erg, which can be potentially detected with JWST2 (Johnson et al. 2013). Returning to the fate of supermassive stars, their evolution in the presence of UV feedback was assessed by Omukai & Inutsuka (2002) and Johnson et al. (2012) in the case of spherical symmetry. These authors find that UV feedback is unable to stop accretion for rates above 0.1 M yr-1 . As a result, the expected stellar mass is MUV 103 M m 10-
8/7 3

·

(36)

For accretion rates of 0.1 M yr-1 , our results exceed this value, as the protostars remain on the Hayashi track up to a mass of 3 â 105 M , making feedback inefficient. Therefore, taking protostellar evolution into account favors the formation of more massive objects. In all cases, a critical question concerns the time evolution of the accretion rate, since the transition to supermassive stars is regulated by the accretion rate at late times. While this paper provides a first assessment for the case of constant accretion rates and spherical symmetry, the potential implications of time-dependent accretion rates need to be addressed in more detail in the future, along with the implications of rotation during protostellar evolution.
Acknowledgements. D.R.G.S. and M.L. thank for funding from the Deutsche Forschungsgemeinschaft (DFG) via the SFB 963/1 "Astrophysical flow instabilities and turbulence" (project A12). DRGS further acknowledges financial support via the Schwerpunktprogramm SPP 1573 "Physics of the Interstellar Medium" under grant SCHL 1964/1-1. We thank the anonymous referee for valuable comments that helped to improve our manuscript.

Appendix A: Additional nuclear processes
In Sect. 3.1, we have shown that the transition to the CNO cycle rapidly occurs, implying that the latter will regulate the formation of a helium core as discussed in Sect. 3.2. However, additional nuclear processes are expected to simultaneously occur, which we discuss here for completeness. In the two sub-sections below, we assess the role of deuterium shell burning as well as the impact of hydrogen burning via the pp-chain during the stellar evolution, showing that these processes will only have a minor impact on the stellar evolution.
A.1. Deuterium burning

As reported by Hosokawa & Omukai (2009), deuterium shell burning may occur even before nuclear burning starts in the core. The luminosity from deuterium shell burning is given as LD = 1.5 â 105 L m 0 .1 [D/H] 2.5 â 10-
5

,

( A .1 )

where [D/H] denotes the deuterium abundance relative to hydrogen. For an accretion rate of m = 0.1, Hosokawa et al. (2012a) demonstrated that deuterium shell burning starts only at protostellar masses of 80 M , and presumably even higher masses in case of higher accretion rates. A comparison with the Eddington accretion rate in Eq. (1) thus shows that the luminosity from deuterium shell burning will never become dominant in this mass
1 2

LISA webpage: http://lisa.nasa.gov/ JWST webpage: http://www.jwst.nasa.gov/ A59, page 7 of 9


A&A 558, A59 (2013)

range, but in fact becomes increasingly less relevant for higher masses. It is also clear that deuterium shell burning is not going to occur close to the outer surface. While the characteristic temperature for a given mass shell scales as T GM/R, we note that, from Eq. (10), GM/R t = M/ M at times t M/ M . In this regime, the central temperature is almost spatially constant, consistent with the results of Hosokawa et al. (2012a), but already considerably larger than the deuterium burning temperature of 106 K. Deuterium burning will thus occur when the first term in Eq. (10) is still relevant. From Fig. 1, it is evident that the radius of the shells changes considerably in that regime, while the mass in the shells is almost unchanged. Adopting a scaling relation of T GM/R with almost constant M implies that the radius has to change by 2­3 orders of magnitude for the atmospheric temperature of 5000 K to increase to a value above 106 K. We therefore expect that the radius of deuterium burning will correspond to a fixed fraction of the protostellar radius as long as the protostar maintains its bloated envelope. From a comparison with Hosokawa et al. (2012a) at a stellar mass of 1000 M , we obtain the normalization of this relation to be RD = 10 R M M
1/3

Fig. A.1. Ratio of nuclear luminosity provided via the pp-cycle (Eq. (A.4)) over the Eddington luminosity of the star (Eq. (A.5)) as a function of protostellar mass for different mass accretion rates. Higher accretion rates imply that less time was available for contraction, thus reducing the nuclear luminosities.

·

( A .2 )

As long as the protostar remains on the Hayashi track, the deuterium burning shell will not be able to move toward the atmosphere, and, therefore, it does not have an impact on the evolution of the protostar.
A.2. The pp-cycle

For practical accretion rates of m = 10-3 - 10 and protostellar masses between 103 -108 M , the contribution from the pp-chain will never be relevant, and can thus be neglected for the overall evolution of the stars. We now estimate the amount of helium produced in the core. During one fusion event, an average energy of 4.3 â 10-5 erg is released, implying a helium production rate of NHe
,pp

=

We consider first the energy production in the nuclear core via the pp-cycle, as the composition of the star is initially primordial. As shown by Hosokawa et al. (2012a), the typical temperature in the core is 108 K, and the nuclear density 1 g cm-3 . We adopt here the expression of Padmanabhan (2000) for the energy production rate, noting that
pp

The helium mass in the core thus evolves as MHe
,pp

Lpp · 4.3 â 10-5 erg

( A .6 )

= 4mp â NHe,pp . = 6.0 â 10
-8 5/2 t710

( A .7 )

An integration yields MHe
,pp

=

2.4 â 10 X T
2/3 9

4

2

M.

( A .8 )

e

-3.38/T

1/3 9

erg s

-1

g,

-1

( A .3 )

where T 9 denotes the temperature in units of 109 K and X 0.75 the mass fraction of hydrogen. While the pp-cycle considered here has only a moderate temperature dependence, we note that the CNO cycle scales as T 20 around temperatures of 106 K, and the triple- process as T 40 around temperatures of 108 K. For the latter cases, we will thus need to take into account the higher temperatures within the center of the core, as they substantially contribute to the energy production via nuclear burning. Now, with nuc 1 g cm-3 and T 9 = 0.1, we have pp = 43 erg g-1 s-1 . The luminosity provided by the pp-cycle is then given as Lpp =
pp 3/2 Mnuc = 21.5 t710 L .

Equating the resulting mass with the total mass in the core, a helium core can be expected after a time of 1.2 â 1010 yr. In summary, when only the pp-cycle is considered, nuclear burning is rather inefficient, and never going to have a significant impact on the evolution of the star.

References
Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 295, 93 Agarwal, B., Khochfar, S., Johnson, J. L., et al. 2012, MNRAS, 425, 2854 Ball, W. H., Tout, C. A., Zytkow, A. N., & Eldridge, J. J. 2011, MNRAS, 414, 2751 Ball, W. H., Tout, C. A., & Zytkow, A. N. 2012, MNRAS, 421, 2713 Begelman, M. C. 2010, MNRAS, 402, 673 Begelman, M. C., & Shlosman, I. 2009, ApJ, 702, L5 Begelman, M. C., Volonteri, M., & Rees, M. J. 2006, MNRAS, 370, 289 Begelman, M. C., Rossi, E. M., & Armitage, P. J. 2008, MNRAS, 387, 1649 Bromm, V., & Loeb, A. 2003, ApJ, 596, 34 Bromm, V., & Larson, R. B. 2004, ARA&A, 42, 79 Choi, J.-H., Shlosman, I., & Begelman, M. C. 2013, ApJ, 774, 149 Clark, P. C., Glover, S. C. O., & Klessen, R. S. 2008, ApJ, 672, 757 Clark, P. C., Glover, S. C. O., Smith, R. J., et al. 2011, Science, 331, 1040 de Souza, R. S., Ciardi, B., Maio, U., & Ferrara, A. 2013, MNRAS, 428, 2109 Devecchi, B., & Volonteri, M. 2009, ApJ, 694, 302 Devecchi, B., Volonteri, M., Colpi, M., & Haardt, F. 2010, MNRAS, 409, 1057 Devecchi, B., Volonteri, M., Rossi, E. M., Colpi, M., & Portegies Zwart, S. 2012, MNRAS, 421, 1465

( A .4 )

For comparison, the Eddington luminosity of the protostar is given as LE
dd

The ratio of these luminosities is given in Fig. A.1, showing that it remains considerably smaller than unity for stellar masses up to 108 M . Equating the two expressions, we find that the luminosity resulting from the pp-burning is relevant only at very late times tpp = 1.1 â 1015 yr.
A59, page 8 of 9

= 3.8 â 104 M t L = 2.7 â 107 m t710 L .

( A .5 )


D. R. G. Schleicher et al.: Massive black hole factories
Dijkstra, M., Haiman, Z., Mesinger, A., & Wyithe, J. S. B. 2008, MNRAS, 391, 1961 Dopcke, G., Glover, S. C. O., Clark, P. C., & Klessen, R. S. 2011, ApJ, 729, L3 Dopcke, G., Glover, S. C. O., Clark, P. C., & Klessen, R. S. 2013, ApJ, 766, 103 Dotan, C., & Shaviv, N. J. 2012, MNRAS, 427, 3071 Eggleton, P. P. 1971, MNRAS, 151, 351 Fan, X., Naraynan, V. K., Lupton, R. H., et al. 2001, AJ, 122, 2833 Fan, X., Strauss, M. A., Schneider, D. P., et al. 2003, AJ, 125, 1649 Fan, X., Hennawi, J. F., Richards, G. T., et al. 2004, AJ, 128, 515 Fan, X., Strauss, M. A., Richards, G. T., et al. 2006, AJ, 131, 1203 Federrath, C., Sur, S., Schleicher, D. R. G., Banerjee, R., & Klessen, R. S. 2011, ApJ, 731, 62 Fryer, C. L., & Heger, A. 2011, Astron. Nachr., 332, 408 Fryer, C. L., Woosley, S. E., & Heger, A. 2001, ApJ, 550, 372 Greif, T. H., Springel, V., White, S. D. M., et al. 2011, ApJ, 737, 75 Greif, T. H., Bromm, V., Clark, P. C., et al. 2012, MNRAS, 424, 399 Haiman, Z. 2006, New Astron. Rev., 50, 672 Hayashi, C., Hoshi, R., & Sugimoto, D. 1962, Progr. Theoret. Phys. Suppl., 22, ¯ 1 Heger, A., & Woosley, S. E. 2002, ApJ, 567, 532 Hosokawa, T., & Omukai, K. 2009, ApJ, 703, 1810 Hosokawa, T., Omukai, K., Yoshida, N., & Yorke, H. W. 2011, Science, 334, 1250 Hosokawa, T., Omukai, K., & Yorke, H. W. 2012a, ApJ, 756, 93 Hosokawa, T., Yoshida, N., Omukai, K., & Yorke, H. W. 2012b, ApJ, 760, L37 Inayoshi, K., & Omukai, K. 2012, MNRAS, 422, 2539 Johnson, J. L., Khochfar, S., Greif, T. H., & Durier, F. 2011, MNRAS, 410, 919 Johnson, J. L., Whalen, D. J., Fryer, C. L., & Li, H. 2012, ApJ, 750, 66 Johnson, J. L., Whalen, D. J., Even, W., et al. 2013, ApJ, submitted [arXiv:1304.4601] Klessen, R. S., Glover, S. C. O., & Clark, P. C. 2012, MNRAS, 421, 3217 Koushiappas, S. M., Bullock, J. S., & Dekel, A. 2004, MNRAS, 354, 292 Latif, M. A., Schleicher, D. R. G., Spaans, M., & Zaroubi, S. 2011, A&A, 532, A66 Latif, M. A., Schleicher, D. R. G., Schmidt, W., & Niemeyer, J. 2013a, MNRAS, 430, 588 Latif, M. A., Schleicher, D. R. G., Schmidt, W., & Niemeyer, J. 2013b, MNRAS, 432, 668 Latif, M. A., Schleicher, D. R. G., Schmidt, W., & Niemeyer, J. 2013c, MNRAS, 433, 1607 Lodato, G., & Natarajan, P. 2006, MNRAS, 371, 1813 ´ Milosavljevic, M., Bromm, V., Couch, S. M., & Oh, S. P. 2009a, ApJ, 698, 766 ´ Milosavljevic, M., Couch, S. M., & Bromm, V. 2009b, ApJ, 696, L146 Mortlock, D. J., Warren, S. J., Venemans, B. P., et al. 2011, Nature, 474, 616 Omukai, K. 2001, ApJ, 546, 635 Omukai, K. 2012, PASJ, 64, 114 Omukai, K., & Inutsuka, S.-I. 2002, MNRAS, 332, 59 Omukai, K., & Palla, F. 2001, ApJ, 561, L55 Omukai, K., & Palla, F. 2003, ApJ, 589, 677 Omukai, K., Tsuribe, T., Schneider, R., & Ferrara, A. 2005, ApJ, 626, 627 Omukai, K., Schneider, R., & Haiman, Z. 2008, ApJ, 686, 801 Padmanabhan, T. 2000, Theoretical Astrophysics ­ Vol. 1, Astrophys. Proc. (Cambridge University Press) Pols, O. R., Tout, C. A., Eggleton, P. P., & Han, Z. 1995, MNRAS, 274, 964 Prieto, J., Jimenez, R., & Haiman, Z. 2013 [arXiv:1301.5567] Rees, M. J. 1984, ARA&A, 22, 471 Regan, J. A., & Haehnelt, M. G. 2009, MNRAS, 396, 343 Reisswig, C., Ott, C. D., Abdikamalov, E., et al. 2013 [arXiv:1304.7787] Schleicher, D. R. G., Spaans, M., & Glover, S. C. O. 2010, ApJ, 712, L69 Schneider, R., Ferrara, A., Salvaterra, R., Omukai, K., & Bromm, V. 2003, Nature, 422, 869 Schneider, R., Omukai, K., Bianchi, S., & Valiante, R. 2012a, MNRAS, 419, 1566 Schneider, R., Omukai, K., Limongi, M., et al. 2012b, MNRAS, 423, L60 Shang, C., Bryan, G. L., & Haiman, Z. 2010, MNRAS, 402, 1249 Shapiro, S. L. 2005, ApJ, 620, 59 Shapiro, S. L., & Teukolsky, S. A. 1986, Black Holes, White Dwarfs and Neutron Stars: The Physics of Compact Objects Smith, R. J., Glover, S. C. O., Clark, P. C., Greif, T., & Klessen, R. S. 2011, MNRAS, 414, 3633 Smith, R. J., Hosokawa, T., Omukai, K., Glover, S. C. O., & Klessen, R. S. 2012, MNRAS, 424, 457 Spaans, M., & Silk, J. 2006, ApJ, 652, 902 Stacy, A., Greif, T. H., & Bromm, V. 2010, MNRAS, 403, 45 Stacy, A., Greif, T. H., & Bromm, V. 2012, MNRAS, 422, 290 Sur, S., Schleicher, D. R. G., Banerjee, R., Federrath, C., & Klessen, R. S. 2010, ApJ, 721, L134 Susa, H. 2013, ApJ, 773, 185 Turk, M. J., Abel, T., & O'Shea, B. 2009, Science, 325, 601 Turk, M. J., Oishi, J. S., Abel, T., & Bryan, G. L. 2012, ApJ, 745, 154 Van Borm, C., & Spaans, M. 2013, A&A, 553, L9 Volonteri, M., & Bellovary, J. 2012, Rep. Progr. Phys., 75, 124901 Volonteri, M., Lodato, G., & Natarajan, P. 2008, MNRAS, 383, 1079 Whalen, D. J., & Fryer, C. L. 2012, ApJ, 756, L19 Wise, J. H., Turk, M. J., & Abel, T. 2008, ApJ, 682, 745

A59, page 9 of 9