Old pre-main-sequence Stars: Disc reformation by Bondi-Hoyle ...

c ESO 2014

Astronomy & Astrophysics manuscript no. BHaccretion May 26, 2014

Letter to the Editor

Old pre-main-sequence stars Disc reformation by Bondi-Hoyle accretion P. Scicluna1,2 , G. Rosotti3,4,5 , J.E. Dale4 , and L. Testi2,4,6 1 2 3 4

arXiv:1405.6051v1 [astro-ph.SR] 23 May 2014

5 6

ITAP, Universit¨at zu Kiel, Leibnizstr. 15, 24118 Kiel, Germany European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching b. M¨unchen, Germany Max-Planck-Institut f¨ur extraterrestrische Physik, Giessenbachstraße, D-85748 Garching, Germany Excellence Cluster Universe, Boltzmannstr. 2, D-85748 Garching, Germany Universitats-Sternwarte M¨unchen, Scheinerstraße 1, D-81679 M¨unchen, Germany INAF-Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, I-50125 Firenze, Italy

Received 17 February 2014; Accepted 21 May 2014 ABSTRACT

Young stars show evidence of accretion discs which evolve quickly and disperse with an e-folding time of ∼ 3Myr. This is in striking contrast with recent observations that suggest evidence for numerous > 30 Myr old stars with an accretion disc in large star-forming complexes. We consider whether these observations of apparently old accretors could be explained by invoking Bondi-Hoyle accretion to rebuild a new disc around these stars during passage through a clumpy molecular cloud. We combine a simple Monte Carlo model to explore the capture of mass by such systems with a viscous evolution model to infer the levels of accretion that would be observed. We find that a significant fraction of stars may capture enough material via the Bondi-Hoyle mechanism to rebuild a disc of mass & 1 minimum-mass solar nebula, and . 10% accrete at observable levels at any given time. A significant fraction of the observed old accretors may be explained with our proposed mechanism. Such accretion may provide a chance for a second epoch of planet formation, and have unpredictable consequences for planetary evolution. Key words. Accretion, accretion discs – circumstellar matter – Stars:formation – Stars: pre-main sequence – Protoplanetary discs

1. Introduction Circumstellar discs form around protostars as a result of angular momentum conservation during gravitational collapse (e.g. Shu et al. 1987). In the early phases of star formation, disc material loses angular momentum and is accreted onto the central star. The most direct observational signature of the presence of a protoplanetary disc is the excess emission, on top of the expected naked stellar photosphere, at infrared and millimetre wavelengths, in the ultraviolet and in optical/infrared emission lines. The long wavelength emission is produced by a dusty disc, heated by internal dissipation processes or reprocessing of stellar radiation (e.g. Dullemond et al. 2007). The short wavelength excess and the optical/infrared emission lines are thought to be produced by the disc-star interaction as matter accretes onto the star or is ejected in a wind/jet (Hartmann 2009). Strong observational evidence shows that both the inner dusty disc and accretion onto the central star quickly disappear during the early stages of pre-main-sequence evolution; the fractions of stars with near infrared excess and with accretion signatures decay with an efolding time of 2-3 Myr (Fedele et al. 2010; Hern´andez et al. 2007). This disc dissipation timescale, even considering the pos-

Send offprint requests to: [email protected]

sible revision by Bell et al. (2013), sets a stringent constraint on the timescales for planet formation. Recent work has challenged this paradigm. Sensitive, wide field Hα surveys of large star-forming complexes in the Magellanic Clouds and our own Galaxy have revealed a population of pre-main-sequence stars that appear to be older than 10 Myr but still show prominent Hα emission and/or infrared excess (Beccari et al. 2010; De Marchi et al. 2013a, 2011a,c, 2013b). Although some of these “old” accretor candidates in nearby star-forming regions have been shown to be misclassified young stellar objects (Manara et al. 2013), it is difficult to believe that this is the case for all the candidates; these populations of old accretors are not as centrally condensed as the young stellar clusters in the same fields (e.g. De Marchi et al. 2011b). If the line emission is interpreted as due to accretion as in young pre-main-sequence stars, the implied accretion rates are similar to those derived at early ages, and typically higher than nearby transitional discs1 . These findings are hard to understand in a framework in which the primordial disc is still the reservoir of accreting material at such old ages; even one disc of age > 30Myr implies an initial population > 105 (assuming exponential decay with an e-folding timescale of 3Myr). 1

although these are systematically lower mass objects 1

P. Scicluna et al.: Old pre-main-sequence stars

between 0.7M and 3.2M 2 from a Salpeter IMF (M ∝ M −α Salpeter 1955) and velocities generated assuming a velocity dispersion of σv = 1 km s−1 , we sample the parameters required in Eq. 2 from the values given in Table 1. The model simulates 10Myr treated as a series of quasi-static time steps of length tst = 2Rcl /v∗ , assuming that each star is independent. For each star, we calculate RBH , the volume swept out per time-step Vst = v∗ tst × π (Rcl + RBH )2 , and hence the probability of encountering a dense clump

Table 1. Parameters for Monte Carlo models Parameter

Values

Parameter

Values

fV Nstars Rcl ncl

10−2 , 10−3 , 10−4 , 10−5 105 , 105 , 106 , 107 0.1 pc 104 cm−3

cs σv α

0.3 km s−1 1 km s−1 2.35

In this paper we explore the possibility that the old accretors do not have a primordial disc, but a disc that they re-accreted after the primordial disc had dissipated. Previous studies (Moeckel & Throop 2009; Padoan et al. 2005; Throop & Bally 2008) have investigated the influence of Bondi-Hoyle accretion on premain-sequence mass-accretion rates and the protoplanetary disc at earlier phases, during the initial evolution of the disc-star system within the progenitor cloud. Here we investigate the possibility that a star older than 5-10Myr happens to travel through a clumpy molecular cloud, typically unrelated to that in which the star formed, and is able to accrete enough material to form a new accretion disc.

2. Modelling 2.1. Bondi-Hoyle accretion

Hoyle & Lyttleton (1939), Bondi & Hoyle (1944) and Bondi (1952) proposed a mechanism by which objects can capture matter from the interstellar medium. A massive object moving through the ISM causes a perturbation, pulling material toward the object. As the capture of material is roughly symmetrical with respect to the direction of motion of the star, much of the angular momentum of the material cancels out, and hence it is captured by the star to eventually be accreted (Davies & Pringle 1980). The rate at which material is captured is given by ˙ BH = µnvπR2BH , M

(1)

where v is the relative velocity between the star and the ISM, n is the number density of the ISM, and µ is the mean molecular weight (usually taken as 2.3mH ). The gravitational cross-section is given by πR2BH , where RBH is the Bondi-Hoyle radius RBH =

2GM∗ ; v2 + c2s

(2)

cs is the sound-speed of the ISM, typically 0.3 km s−1 . For a 1M star moving at 1km s−1 , RBH ∼ 1500 au. To explore the effect of this process in reconstituting discs around young stars, we build a simple Monte Carlo model to treat interactions between stars and clumps with densities typical for molecular clouds. We assume a stationary clumpy molecular cloud, which we model as a collection of identical spherical clumps with radius Rcl and density ncl . We parametrise the density of clumps through a volume filling factor of dense gas fV . We assume a population of “old” young stars that has lost their primordial disc enters the cloud and moves through the clumpy medium. By randomly generating stars with masses 2

p=

Vst × fV . (4/3) πR3cl

(3)

In each time-step a uniform random number ζ is drawn, and the star encounters a clump when ζ ≤ p; the impact parameter b of the encounter is given by drawing a second random number ζ2 from the same generator such that b = (Rcl + RBH ) ζ21/2 . We then determine the accretion rate (Eq. 1) and resolve the stellar accretion and the clump-mass depletion on a finer time-grid of 1000 sub-steps to accurately determine the accreted mass. Interactions where RBH > Rcl and grazing encounters are treated correctly by taking the projected area of intersection. By repeating this process for > 105 stars we build up meaningful statistics about the range of possible BH accretion histories and their probabilities. Note that each star is modelled independently, and mass accreted by a star does not influence the mass-budget available to later stars. The accretion histories determined by this model are then passed to a viscous evolution model (Sect. 2.2) to estimate the rate at which material is accreted by the star. Our choice of fV is based on a reanalysis of SPH simulations of star-forming regions including feedback mechanisms presented in Dale et al. (2012, 2013) to determine the filling factor of gas at densities higher than 104 cm−3 . We find that for bound clouds of similar stellar mass to the regions observed by Beccari et al. (2010); De Marchi et al. (2013b), 10−6 < fV . 10−3 irrespective of whether feedback from massive stars is included. While this provides a useful estimate of the amount of mass captured in this way, it somewhat overestimates the total as we neglect a number of physical processes. First, we neglect the motion of the clumps and assume that v = v∗ in Eq. 2. Correct treatment of the relative motions would in general reduce RBH and hence the accretion rates. Second, stars above 2M have significant wind and radiation pressure that will depress the accretion rate (Edgar & Clarke 2004). Similarly, we do not include the possible influence of the X-ray photoevaporation on the accretion, which may have an analagous effect for lower mass stars. We also ignore the possible influence of magnetic fields, which recent studies (e.g. Lee et al. 2014) have shown may reduce accretion rates by a factor of a few. Likewise, we neglect structure on scales smaller than a single clump; such structure is required for a disc to form, and would reduce accretion rates relative to the homogeneous clump case treated here. Finally, we do not include binaries. However, the only influence of binarity in the context of Bondi-Hoyle accretion is to increase RBH , since binaries behave as a single object of mass M = M1 + M2 . 2 Stars above ∼ 3M have strong winds which make a simple model inappropriate, while observations of old accretors are incomplete for stars below 0.7-1M depending on the distance to the observed region.

P. Scicluna et al.: Old pre-main-sequence stars

2.2. Viscous evolution modelling

1

fV = 10-2

Fraction

Cumulative fraction

Due to the angular momentum of the material accreted from the fV = 10-3 fV = 10-4 clump, which may be due to a density gradient within the clump -1 fV = 10-5 10 or the rotation of the clump itself, accretion cannot proceed directly onto the star (Ruffert 1997). Therefore, the formation of a thin accretion disc is expected as the result of the viscous spread10-2 ing of a thin ring. Throop & Bally (2008) described the “buffer” effect of an accretion disc, but did not directly model it. We assume that the material accreted from the medium circularises at a radius r0 = 0.1RBH . After a single impulse of accretion onto the 10-3 disc, the surface density is described by Σ(r) = M0 /(2π)δ(r − r0 ), where M0 is the deposited mass. Under the influence of an effective viscosity ν that redistributes the angular momentum in the -3 -2 -1 0 1 disc, the spreading ring solution (Lynden-Bell & Pringle 1974) Accreted mass (Log (M/MO•)) describes the evolution in time of this initial surface density,   1/2  Fig. 1. Cumulative fraction of the stellar population that has ac  (r0 − r1/2 )2 r  GM∗ (rr0 )1/4 mass. The solid blue Σ(r, t) = exp −   exp(−λ)I1/2 (λ), (4) creted mass as a function of total accreted line indicates a filling-factor of 10−2 , the dotted magenta line 3tν 3πr2 νΩ 10−3 , the dashed red line 10−4 , and the dot-dashed green line where ν is the kinematic viscosity of the gas, Ω the keplerian 10−5 . angular speed, I1/2 the modified Bessel function of order 1/2, λ = 2r3/2 /(3(GM∗ )3/2 νtr0 ), and we have specialized the expres1 sion for the ν ∝ r case. From this analytical solution, it is possi˙ ble to compute the mass accretion rate onto the star Mkernel . To derive the mass accretion rate history onto the star, we convolve 10-1 this function with the mass accretion rate history onto the disc: Z 10-2 ˙ ˙ BH (t0 ) M ˙ kernel (t − t0 )dt0 . M∗ (t) = M (5) Given a stellar mass, the loading radius, and a law for viscosity, the evolution in time is now completely determined. We fix the viscosity by using the well-known Shakura & Sunyaev (1973) prescription, ν = α(h/r)2 r2 Ω, where α is the Shakura-Sunyaev parameter and h/r the aspect ratio of the disc. We choose typical values of α = 0.01 and h/r = 0.05(r/1AU)1/4 (Armitage 2011). Operationally, we sample Eq. 4 numerically on a space and time grid. We integrate over space to get the mass of the disc and we numerically differentiate the result to get the mass accretion rate kernel, which can be convolved with the Bondi-Hoyle history (Sect. 2.1).

3. Results Our model indicates that a fraction of the population ∼ 40 − 50 × fV encounter dense regions and accrete more than 0.001M material by the end of the simulation (Fig. 1). The median accreted mass is typically ∼ 0.01M , similar to the mass of discs around young pre-main-sequence stars, with strong dependence on the stellar mass. In extreme cases, however, more massive stars (>2M ) with low v∗ that encounter several clumps can capture ≥M . Our treatment of the disc formation and evolution is probably inadequate for these extreme cases. Converting the Bondi-Hoyle accretion into stellar accretion ˙ ∗ . 10−6 M yr−1 after the formation of the rates, we find M disc. Owing to the assumptions inherent in our model, this rate declines from the peak as a power law as in primordial discs. By calculating the time each star spends accreting above a certain threshold accretion rate, one can derive a mean time per

10-3 10-4 10-9

10-8 10-7 Accretion threshold (MO• yr-1)

10-6

Fig. 2. Fraction of the population that would be detected as an old accretor at a given time, plotted as a function of the instantaneous accretion rate. The models are indicated using the same colours and line-styles as Fig. 1. star as a function of the threshold and hence an estimate of the fraction of the population which one expects to observe accreting at a given time. As shown in Fig. 2, for a threshold rate of 10−8 M yr −1 we typically find that the cumulative probability is ∼ 20 fV , i.e. the fraction of a stellar population that one expects to observe as old accretors at a given time is an order of magnitude larger than the volume filling-factor of dense clumps.

4. Discussion Our primary goal is to assess whether the Bondi-Hoyle mechanism can contribute significantly to observations of old accretors in regions with ongoing star formation, under a number of simple assumptions. This involves stars from a previous star-formation episode, after their primordial discs have dispersed, interacting with a clumpy molecular cloud. Our model 3

P. Scicluna et al.: Old pre-main-sequence stars

indicates that up to several percent of the population passing through a region containing dense clumps may accrete more than 0.001M of material. Because of the factors indicated above (Sect. 2.1), the model is likely to overestimate the total accreted mass. However, since the Bondi-Hoyle accretion is a well-understood process, the largest sources of uncertainty derive from the parameters assumed as input to the model, and in particular the clump geometry and filling factor, as well as the assumption that the accreted material will form a thin disc. Our initial choice of filling factor was based on a reanalysis of the simulations of Dale et al. (2012, 2013) for clouds similar to those observed to host old accretors. A further estimate can be obtained from the high-resolution sub-mm maps of the 30 Dor region from Indebetouw et al. (2013). These reveal a wealth of clumpy structures, similar in scale and density to the clumps in the Monte Carlo model used here. Assuming that the clumps are uniform spheres with an average radius Rcl = 0.15 pc and distributed in a cube whose depth is equal to the projected size of the observed region (10 × 10 × 10 pc3 ) yields a filling factor of fV = 1.5 × 10−3 , at the upper end of our parameter range. The behaviour of the accretion disc depends strongly on the viscous timescale τν , as parametrised in terms of r0 and α. An order of magnitude change in τν has little effect on the observable old-accretor fractions at low thresholds, but the fractions at high thresholds decline approximately in proportion to 1/τν . For larger changes in viscosity, this also affects the lowest thresholds explored in Sect. 3 Since we do not include stars down to the peak of IMF (∼0.3M ) and Bondi-Hoyle accretion rates are ∝ M 2 , we may overestimate the total fraction of old accretors by a factor ∼ 3 for the Salpeter IMF assumed here. However, Eq. 3 is dominated by Rcl for low-mass stars, so one would expect a similar fraction of ˙ is a factor of 4 lower. old accretors when M Comparisons between our model and the observations of old accretors are difficult, as there are no firm constraints on the size of the old population (including non-accretors). Nevertheless, from Fig. 2 one can see that without an unrealistically large filling factor ( 10−3 ) of dense clumps, the small, nearby starforming regions are unlikely to produce more than one old accretor, as their typical mass is a few hundred M . As no old accretors have been identified in these regions, this is consistent with our model. From the recent identification of a large (∼ 3 × 103 M ) diffuse population with ages & 10 Myr toward Orion (Bouy et al. 2014) one expects a few tens of reformed discs, although it is unclear whether there is any overlap between this population and the Orion molecular clouds. Observations of old accretors in large star-forming complexes typically detect up to several hundred such sources in each observed region. Given the formation efficiency we have computed and our assumed filling factors, this requires a total population at least of the order of 104 stars in the mass range of the observed old accretors, or ∼ 3 × 104 stars correcting for the IMF, which must have passed through the regions in which the clumps are distributed. In the case of NGC3603, which is inferred to have a population ∼ 104.2 M (Rahman et al. 2013) and ∼ 100 old accretors, this implies either that the old population was significantly richer, or that fV is or was very high. The 30 Doradus region, on the other hand, shows a similar total of old accretors, although the total population is likely ∼ 100 times 4

larger than NGC3603. Only a small fraction (1%) of the stars in 30Dor need to pass through regions containing dense clumps to produce the observed numbers. In reality, fV will evolve with time, and it is possible that the difference we observe between these regions may be due to 30Dor being more evolved, or having evolved more rapidly, than NGC3603. In our model, a significant fraction (up to several tens of percent) of stars capture enough material to form a circumstellar disc of mass similar to primordial protoplanetary discs. This raises a number of interesting questions, such as whether a second epoch of planet formation is possible, and how the interaction between inflowing material and an existing planetary system might alter the accretion or the planetary evolution. The answers to these queries depend strongly on how the inflowing material interacts with the existing system, which we have not treated. Nevertheless, Bondi-Hoyle accretion presents a mechanism by which a new reservoir of potentially planetforming material may be built by up to a few percent of stars. This gives them a second chance to form planets, from material that is potentially of different composition from the material that formed the star. Another possibility is that these stars are already surrounded by a planetary system formed out of the primordial disc. If they accrete new material, typically with an angular momentum different from that of the original planetary system, the interaction of the new material and the existing planets may have a range of outcomes. Understanding the range of possible outcomes will require detailed simulations of the accretion process and of the dynamical interactions with the planetary systems which are beyond the scope of the present paper.

5. Conclusions We have presented a model in which Bondi-Hoyle accretion by stars passing through dense clumps in the outer regions of their natal molecular cloud leads to the re-formation of a circumstellar disc. As a result, these stars may masquerade as pre-mainsequence objects due to ongoing accretion and the presence of infrared excess emission. A significant part of the observed populations of old accretors in large star-forming regions may be explained by this mechanism. As it may have wide-ranging consequences for the early evolution of planetary systems in rich stellar environments, we believe that further investigation of this mechanism is warranted. Acknowledgements. We wish to thank the anonymous referee for their careful reading of the text. The idea explored in this paper came up during discussions at the ESO science days and star formation coffee as well as the Munich Star Formation workshops. We thank the ESO Office for Science and all the institutes in the Munich area for providing a stimulating environment. We thank P. Armitage, G. Beccari, G. Costigan, B. Ercolano, G. De Marchi, C. Manara, N. Moeckel, A. Natta, P. Padoan, R. Siebenmorgen and S. Wolf for discussions and insights on the various aspects discussed in this paper. PS is supported under DFG programme no. WO 857/10-1. GR acknowledges the support of the International Max Planck Research School (IMPRS). This research was supported by the DFG cluster of excellence ‘Origin and Structure of the Universe’ (JED).

References Armitage, P. J. 2011, ARA&A, 49, 195 Beccari, G., Spezzi, L., De Marchi, G., et al. 2010, ApJ, 720, 1108

P. Scicluna et al.: Old pre-main-sequence stars Bell, C. P. M., Naylor, T., Mayne, N. J., Jeffries, R. D., & Littlefair, S. P. 2013, MNRAS, 434, 806 Bondi, H. 1952, MNRAS, 112, 195 Bondi, H. & Hoyle, F. 1944, MNRAS, 104, 273 Bouy, H., Alves, J., Bertin, E., Sarro, L. M., & Barrado, D. 2014, A&A, 564, A29 Dale, J. E., Ercolano, B., & Bonnell, I. A. 2012, MNRAS, 424, 377 Dale, J. E., Ercolano, B., & Bonnell, I. A. 2013, MNRAS, 430, 234 Davies, R. E. & Pringle, J. E. 1980, MNRAS, 191, 599 De Marchi, G., Beccari, G., & Panagia, N. 2013a, ApJ, 775, 68 De Marchi, G., Panagia, N., Guarcello, M. G., & Bonito, R. 2013b, MNRAS, 435, 3058 De Marchi, G., Panagia, N., Romaniello, M., et al. 2011a, ApJ, 740, 11 De Marchi, G., Panagia, N., & Sabbi, E. 2011b, ApJ, 740, 10 De Marchi, G., Paresce, F., Panagia, N., et al. 2011c, ApJ, 739, 27 Dullemond, C. P., Hollenbach, D., Kamp, I., & D’Alessio, P. 2007, Protostars and Planets V, 555 Edgar, R. & Clarke, C. 2004, MNRAS, 349, 678 Fedele, D., van den Ancker, M. E., Henning, T., Jayawardhana, R., & Oliveira, J. M. 2010, A&A, 510, A72 Hartmann, L. 2009, Accretion Processes in Star Formation: Second Edition (Cambridge University Press) Hern´andez, J., Hartmann, L., Megeath, T., et al. 2007, ApJ, 662, 1067 Hoyle, F. & Lyttleton, R. A. 1939, Proceedings of the Cambridge Philosophical Society, 35, 405 Indebetouw, R., Brogan, C., Chen, C.-H. R., et al. 2013, ApJ, 774, 73 Lee, A. T., Cunningham, A. J., McKee, C. F., & Klein, R. I. 2014, ApJ, 783, 50 Lynden-Bell, D. & Pringle, J. E. 1974, MNRAS, 168, 603 Manara, C. F., Beccari, G., Da Rio, N., et al. 2013, A&A, 558, A114 Moeckel, N. & Throop, H. B. 2009, ApJ, 707, 268 Padoan, P., Kritsuk, A., Norman, M. L., & Nordlund, Å. 2005, ApJ, 622, L61 Rahman, M., Matzner, C. D., & Moon, D.-S. 2013, ApJ, 766, 135 Ruffert, M. 1997, A&A, 317, 793 Salpeter, E. E. 1955, ApJ, 121, 161 Shakura, N. I. & Sunyaev, R. A. 1973, A&A, 24, 337 Shu, F. H., Adams, F. C., & Lizano, S. 1987, ARA&A, 25, 23 Throop, H. B. & Bally, J. 2008, AJ, 135, 2380

5