Making the planetary material diversity during the early assembling of the Solar System
VERSION: ArXiv DATE: July 19, 2019
Chondritic meteorites, the building blocks of terrestrial planets, are made of an out-of-equilibrium assemblage of solids formed at high and low temperatures, either in our Solar system or previous generations of stars. This was considered for decades to result from large scale transport processes in the Sun’s isolated accretion disk. However, mounting evidences suggest that refractory inclusions in chondrites formed contemporaneously with the disk building.
Here we numerically investigate, using a 1D model and several physical and chemical processes, the formation and transport of rocky materials during the collapse of the Sun’s parent cloud and the consequent Solar Nebula assembling.
The interplay between the cloud collapse, the dynamics of gas and dust, vaporization, recondensation and thermal processing of different species in the disk, results in a local mixing of solids with different thermal histories. Moreover, our results also explains the overabundance of refractory materials far from the Sun and their short formation timescales, during the first tens of kyr of the Sun, corresponding to class 0-I, opening new windows into the origin of the compositional diversity of chondrites.
Francesco C. Pignatale
0000-0003-0902-7421]Francesco C. Pignatale
Chondritic meteorites are fragments of planetesimals formed in the Solar Nebula within the first 4 Myr of Solar System history (the time 0 being the formation of the highly refractory inclusions, Ca-, Al- rich, CAIs) (Chaussidon & Liu, 2015). They are made of a mixture of materials that formed at high and low temperatures (Scott & Krot, 2003). This peculiar assemblage has always been puzzling and is one key ingredient in the composition of terrestrial planets. A striking paradox is that carbonaceous chondrites, which presumably accreted in the outer solar system at temperatures lower than 150 K (Warren, 2011), are the richest in CAIs formed at temperatures higher than 1600 K, presumably close to the Sun (Scott & Krot, 2003). The chemical composition of CAIs links them to precursors condensed from the solar gas (Richter et al., 2006).
Absolute Pb-Pb ages show that CAIs in carbonaceous chondrites were formed over a brief interval (time, kyr) at the very beginning of the Solar System (Connelly et al., 2012) from refractory precursors that, according to short-lived chronologies, could have condensed over an even shorter period of 10-20 kyr (Jacobsen et al., 2008; Larsen et al., 2011; Mishra & Chaussidon, 2014). Such a timescale is shorter than the assembling time of a protoplanetary disk from the infalling molecular cloud (Williams & Cieza, 2011). It seems thus unavoidable that the formation of the earliest solar system solids (CAIs at least) occurred concurrently with the building of the Solar Nebula, in agreement with observations of young stellar objects showing that dust grows and becomes largely crystalline as soon as the star begins to form (Williams & Cieza, 2011; Chiang et al., 2012). Since increasing observational and theoretical evidence shows that disks start off compact (Dullemond et al., 2006; Jacquet et al., 2011; Hennebelle et al., 2016; Cox et al., 2017) most of the infalling material should be located initially close to the protostar and undergoes thermal processing before incorporation into planetary material, like chondrites or planetesimals.
To understand the origin of this material, we have developed a 1D simulation to track the formation of the disk and the dust thermal evolution during the collapse of the parent cloud.
2 Methods and models
The code aims at tracking the 1D evolution of a protoplanetary disks, with different gaseous and solid species. We use a simplified treatment of the gas and dust dynamics and chemistry (nevertheless including a variety of relevant coupled processes) in order to unveil key physical processes associated to the coupling of transport and thermal processing in a connected cloud-disk system. This is the reason the infalling envelope is described as a spherical isothermal shell collapsing while conserving angular momentum, in order to make the calculations computationally tractable (see e.g. Hueso & Guillot (2005)). Although more detailed models are available for the cloud collapse, e.g. Hennebelle et al. (2016), and numerous are the studies investigating the disk chemistry (Henning & Semenov, 2013), our approach represents the first attempt to implement in one single simulation several dynamical and chemical aspects of cloud collapse and disk formation and evolution.
The code is a classical alpha-disk model built on the model of Hueso & Guillot (2005) that considers cloud collapse and radiative and viscous heating. Our 1D disk is described by a logarithmic grid (100 cells, between 0.01 and 1000 au). We include gravitational instability, dead zone in the form of a parameterized description (Terquem, 2008). We use a layered accretion model of Zhu et al. (2010) with a low- midplane and high- surface layers, dust transport through advection and diffusion, and growth/fragmentation using the method desribed by Dra̧żkowska et al. (2016) and Birnstiel et al. (2012). The opacity is calculated using the opacity table from Baillié et al. (2016) and procedures described in Helling et al. (2000) and in Semenov et al. (2003). Starting from the work of Yang & Ciesla (2012) on refractory dust, we now include multiple chemical species allowing material condensation and sublimation using a simple chemical model.
The infall of the parent cloud is described as a source term and here below we describe the main aspects. The cloud collapses onto the protostar-disk system with a constant accretion rate (Shu, 1977):
where is the gravitational constant and is the isothermal sound speed111 where is the Boltzmann constant, is the cloud temperature, is the molecular weight of the gas in terms of the proton mass, . in the cloud. Assuming that the angular momentum of the cloud material is conserved during the infall, the surface density accreted below a threshold radius, called “centrifugal radius”, is (Hueso & Guillot, 2005; Yang & Ciesla, 2012):
where the centrifugal radius is the radius under which the cloud collapses onto the disk. is initially so small that all infalling material falls directly onto the protostar. As the star mass increases, increases and the cloud material feeds progressively the disk. is determined as the place where the specific angular momentum of the infalling material is equal to the specific angular momentum of the keplerian disk. Thus, deeper cloud envelopes fall close to the protostar while exterior cloud envelopes fall far from the protostar. The expression for the centrifugal radius is
where is the constant angular velocity of the cloud, is the temperature of the cloud and is the total mass of the protostar+disk system at time .
Concerning the chemistry, a simplified approach is used where simple phase transitions are considered. Using a thermochemical code to compute in real time the exact chemical evolution of the different species in the simulation is very challenging. Indeed, at each time step and location the local gas and dust composition changes. This implies the implementation and tracking of at least several dozens of compounds if only a âshortâ list of elements were considered (see for example the shortened list of 60 species that can be stable in protoplanetary disk conditions derived in Pignatale et al. (2011) using the fifteen most abundant Sun’s elements (Asplund et al., 2009)).
We, thus, introduce a simplified chemical model (but nevertheless physically motivated, as the temperatures of condensation/processing of numerous species are well established (Barshay & Lewis, 1976; Ebel, 2006)), in order to preserve the main trends of the chemistry. In our simulations, the bulk chemical composition of the parent cloud is assumed to be solar (Asplund et al., 2009; Pignatale et al., 2011) with all the major rock-forming elements (Si, Mg, Al, Ca, Fe) hosted in interstellar unprocessed solid phases. We allow the infalling interstellar materials to be either vaporised, âprocessedâ (i.e. heated at temperatures K enough to thermally alter the dust and destroy presolar grains but below vaporisation) or unaltered upon its arrival in the disk.
The model tracks the formation and transport in the disk of gaseous species and different solid species with their own vaporisation temperatures (See Table 1: (i) condensed refractory phases, which are the potential precursors of CAIs present in meteorites, named hereafter, for sake of simplicity, CAIs, (ii) condensed metal named hereafter metallic iron, (iii) condensed silicates, (iv) thermally processed and unprocessed silicate dust, (v) thermally processed and unprocessed iron dust, (vi) water ice, and (vii) CO ice. Molecular hydrogen is then the main gas species.
Chondrules formation is neglected, as we essentially model their precursors. However, chondrule formation is not expected to affect large-scale compositional fractionations in the disk (Goldberg et al., 2015). Heating from accretion shock is also neglected because of the low temperatures involved (Dominik & Dullemond, 2008).
|species||T (K)||become||T (K)||become|
|unprocessed refractories||considered as silicates|
|unprocessed silicates||processed silicates||silicate(g)|
|unprocessed iron||processed iron||iron(g)|
|metallic iron||processed iron|
|water ice||water vapour||water ice|
|CO ice||CO(g)||CO ice|
3 Results and Discussion
Figure 1 shows the time evolution of the centrifugal radius and the locations of disk outer edge, different condensation fronts, and dead zone (dark grey). With our initial conditions ( K, rad/s, , K, , , , , , ), the cloud infall and diskâs assemblage lasts kyr. While being fed, the disks also expands due to viscous spreading (Hueso & Guillot, 2005; Yang & Ciesla, 2012). After kyr the material starts flowing mostly inwards as the disk accretes onto the proto-Sun. During the infall, a dead zone (DZ), a region of low diffusivity (Turner & Sano, 2008; Terquem, 2008), develops between and AU.
As the dust is incorporated into the disk, its composition changes because of heating of minerals with different sublimation temperatures. Our simulations show that the final distribution of minerals is controlled by the interplay between two key parameters: (i) the centrifugal radius, (t), i.e. the distance (increasing with time) within which the parent cloud material is injected and (ii) the location of the different condensation fronts and properties of the different mineral species, parameters never fully considered previously. Because of angular momentum constraints, (t) is initially close to the forming Sun fed from the parent cloud’s innermost regions (Kimura et al., 2016). Then (t) increases with time due to the arrival of originally more distant cloud materials with higher specific angular momentum. Under our initial conditions, (t) increases from 0 to AU when the infall ends (Fig. 1). The location of the condensation fronts of the different solid species considered depends on the disk temperature profile which is, itself, controlled by viscous and radiative heating (Charnoz et al., 2018). For a given species, if is inside the condensation front of that species, then the dust will be totally vaporised upon its injection in the disk. However, the disk spreads outward concurrently, so that gas injected at high temperature cools down while moving outward, allowing its condensation at larger heliocentric distance. Moreover, dust injected at later stage, in a lower temperature environment, can be just processed without being vaporised, or survive thermally unaltered. These processes are sketched in Fig. 1.
Figures 2 and 3 show four time-snapshots of the surface density profiles of gas, dust and ices and the radial distribution of the mass fraction of the considered rocky material. Our results show that at kyr (Fig. 3(a)) the disk is only populated by condensed dust (CAIs, iron and silicates) and only one kind of processed dust: condensed metallic iron that could experience low temperature chemical alteration ( K, see Table 1). Later, as moves outward, parent cloud dust can be preserved if a given mineral is injected outside its condensation front. At kyr, ) penetrates the DZ (Fig. 3(b)) where dust and gas can be stored because of low turbulence. This lowers dramatically the concentration of refractory condensates within the DZ. At the end of the infall, most of the DZ has been fed by parent cloud material, leading to a DZ poor in refractory condensates and a CAIs-enriched outer disk beyond the DZ (Fig. 3(c)). Inside the inner half of the DZ, material can be further thermally processed as the temperature could reach K while the outer half of the DZ contains mostly thermally unaltered dust. Our simulations also show that this heterogeneous dust distribution will be qualitatively preserved over the lifetime of the disk (Fig. 3(d)). However, inward the DZ inner edge ( AU), the disk becomes CAIs-poor because of strong inward accretion of gas that removes solids.
This distribution predicted for the various dust components in the disk (Fig. 3(d)) shows striking similarities with the composition of the three main families of chondritic meteorites (enstatite, ordinary, carbonaceous chondrites) as defined by their fractions of CAIs, metal, chondrules and matrix (Scott & Krot, 2003). To do this exercise of comparison, we consider to first order the unprocessed dust as a proxy for the fine-grained matrix containing different kinds of circumstellar grains (Nittler, 2003) and GEMS (glass with embedded metal and sulfides) which, for a fraction of them, are likely surviving interstellar grains (Keller & Messenger, 2011). Within 1-10 AU, dust is dominated by disk-born silicates (processed or condensed), iron (processed or condensed, inside 2 AU essentially) and poor in CAIs and pristine interstellar dust. This is reminiscent of (matrix-poor, chondrule-rich) non-carbonaceous chondrites, i.e. ordinary chondrites and enstatite chondrites (the latter likely formed well inside the water snow line). Further away in the disk, the dust is dominated by unprocessed ISM dust that was injected late when AU, processed dust and condensed silicates. This region, though beyond the snow-line and thus having water ice-bearing dust, has the highest fraction of CAIs of the whole disk. This bears similarities to carbonaceous chondrites. Note the presence, in addition to CAIs, of condensed silicates formed from the same parent gas reservoir. These silicates could then correspond to the Amoeboid Olivine Aggregates (AOAs) of carbonaceous chondrites that are known to have strong links with CAIs (from their oxygen isotopic composition, see Krot et al. (2004)).
Interestingly, formation of CAIs is found to be a transient event. In Fig 4(a) we report the time evolution of the CAIs mass in the disk and the time evolution of the accretion rates. In Fig. 4(b) we show the time evolution of the temperature in a given location of the disk and the time evolving location of the CAIs condensation front.
Massive production of CAIs happens, in the simulation, from kyr to kyr, in very good agreement with Pb–Pb and – chronologies of CAIs (Connelly et al., 2012; Mishra & Chaussidon, 2014) (Fig. 4(a)). This corresponds to a short time window, when interstellar solids are injected in the disk at temperatures high enough for them to be fully evaporated (i.e. when is inward the condensation radius of CAIs, see Fig.1), so that upon spreading outward, the gas can condense into CAIs. This time is shorter by a factor of a few than that found by Yang & Ciesla (2012) because of their lower condensation temperature for refractories.
Still, the production of CAIs appears to occur by steps, linked to accretion bursts onto the star following mass loading and instabilities of the DZ (Zhu et al., 2009; Owen & Armitage, 2014) (Fig. 4(a)). This may help to explain part of the diversity of CAI populations across chondrite groups (MacPherson, 2014) despite overall derivation from the same reservoir (the inner disk in our model). Possibly, some of the reheating episodes experienced by many CAIs (Richter et al., 2006; Simon et al., 2017) may be linked to such bursts of accretion, re-heating, partial vaporisation condensations (as suggested in Fig. 4(a,b)) although many of those episodes were likely of much shorter duration. Our results set the formation of the main bulk of CAIs precursors found in chondrites as early as the formation of the disk itself (class 0 to I).
In this work we investigated the effects of the cloud collapse and disk building in forming and distributing different solids species in the disk. We find that coupling dust transport and thermal processing during the cloud collapse and disk formation naturally results in local mixing of materials with different thermal histories, in qualitative agreement with the three main chondritic families. Our model also explains the overabundance of refractory materials far from the Sun and their short formation timescales, during the first tens of kyr of the Sun. Our finding can reconcile all these long-existing cosmochemical questions within one single picture and place the time 0 of the Solar System during the class 0 to I of the Sun’s formation history.
Some additional processes (e.g. disk winds acting at the disk surface) may redistribute material at large distance, and a more refined chemistry can give more precise predictions. However, this simple model unveils new ways, never described before, to interpret the distribution of chondritic material in our Solar System, that may have kept the imprint of the could infall. This works implies that it is now critical to study the protoplanetary disk, not in isolation, but in the context of the early environment of the Sun which has left fingerprints in the diversity of our Solar System material.
- Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
- Baillié et al. (2016) Baillié, K., Charnoz, S., & Pantin, E. 2016, A&A, 590, A60
- Barshay & Lewis (1976) Barshay, S. S., & Lewis, J. S. 1976, ARA&A, 14, 81
- Birnstiel et al. (2012) Birnstiel, T., Klahr, H., & Ercolano, B. 2012, A&A, 539, A148
- Charnoz et al. (2018) Charnoz, S., Pignatale, F.C., Hyodo, R., et al. 2018, A&A, submitted
- Chaussidon & Liu (2015) Chaussidon, M., & Liu, M.-C. 2015, Washington DC American Geophysical Union Geophysical Monograph Series, 212, 1
- Chiang et al. (2012) Chiang, H.-F., Looney, L. W., & Tobin, J. J. 2012, ApJ, 756, 168
- Connelly et al. (2012) Connelly, J. N., Bizzarro, M., Krot, A. N., et al. 2012, Science, 338, 651
- Cox et al. (2017) Cox, E. G., Harris, R. J., Looney, L. W., et al. 2017, ApJ, 851, 83
- Dominik & Dullemond (2008) Dominik, C., & Dullemond, C. P. 2008, A&A, 491, 663
- Dra̧żkowska et al. (2016) Dra̧żkowska, J., Alibert, Y., & Moore, B. 2016, A&A, 594, A105
- Dullemond et al. (2006) Dullemond, C. P., Apai, D., & Walch, S. 2006, ApJ, 640, L67
- Ebel (2006) Ebel, D. S. 2006, Meteorites and the Early Solar System II, 253
- Goldberg et al. (2015) Goldberg, A. Z., Owen, J. E., & Jacquet, E. 2015, MNRAS, 452, 4054
- Helling et al. (2000) Helling, C., Winters, J. M., & Sedlmayr, E. 2000, A&A, 358, 651
- Hennebelle et al. (2016) Hennebelle, P., Commerçon, B., Chabrier, G., & Marchand, P. 2016, ApJ, 830, L8
- Henning & Semenov (2013) Henning, T., & Semenov, D. 2013, Chemical Reviews, 113, 9016
- Hueso & Guillot (2005) Hueso, R., & Guillot, T. 2005, A&A, 442, 703
- Jacobsen et al. (2008) Jacobsen, B., Yin, Q.-z., Moynier, F., et al. 2008, Earth and Planetary Science Letters, 272, 353
- Jacquet et al. (2011) Jacquet, E., Fromang, S., & Gounelle, M. 2011, A&A, 526, L8
- Keller & Messenger (2011) Keller, L. P., & Messenger, S. 2011, Geochim. Cosmochim. Acta, 75, 5336
- Kimura et al. (2016) Kimura, S. S., Kunitomo, M., & Takahashi, S. Z. 2016, MNRAS, 461, 2257
- Krot et al. (2004) Krot, A. N., Petaev, M. I., Russell, S. S., et al. 2004, Chemie der Erde / Geochemistry, 64, 185
- Larsen et al. (2011) Larsen, K. K., Trinquier, A., Paton, C., et al. 2011, ApJ, 735, L37
- MacPherson (2014) MacPherson, G. J. 2014, Meteorites and Cosmochemical Processes, 139
- Mendybaev et al. (2013) Mendybaev, R. A., Richter, F. M., Georg, R. B., et al. 2013, Geochim. Cosmochim. Acta, 123, 368
- Mishra & Chaussidon (2014) Mishra, R. K., & Chaussidon, M. 2014, Earth and Planetary Science Letters, 390, 318
- Nittler (2003) Nittler, L. R. 2003, Earth and Planetary Science Letters, 209, 259
- Owen & Armitage (2014) Owen, J. E., & Armitage, P. J. 2014, MNRAS, 445, 2800
- Pignatale et al. (2011) Pignatale, F. C., Maddison, S. T., Taquet, V., Brooks, G., & Liffman, K. 2011, MNRAS, 414, 2386
- Richter et al. (2006) Richter, F. M., Mendybaev, R. A., & Davis, A. M. 2006, Meteoritics and Planetary Science, 41, 83
- Simon et al. (2017) Simon, J. I., Jordan, M. K., Tappa, M. J., et al. 2017, Earth and Planetary Science Letters, 472, 277
- Scott & Krot (2003) Scott, E. R. D., & Krot, A. N. 2003, Treatise on Geochemistry, 1, 711
- Semenov et al. (2003) Semenov, D., Henning, T., Helling, C., Ilgner, M., & Sedlmayr, E. 2003, A&A, 410, 611
- Shu (1977) Shu, F. H. 1977, ApJ, 214, 488
- Terquem (2008) Terquem, C. E. J. M. L. J. 2008, ApJ, 689, 532
- Turner & Sano (2008) Turner, N. J., & Sano, T. 2008, ApJ, 679, L131
- Warren (2011) Warren, P. H. 2011, Geochim. Cosmochim. Acta, 75, 6912
- Williams & Cieza (2011) Williams, J. P., & Cieza, L. A. 2011, ARA&A, 49, 67
- Yang & Ciesla (2012) Yang, L., & Ciesla, F. J. 2012, Meteoritics and Planetary Science, 47, 99
- Zhu et al. (2009) Zhu, Z., Hartmann, L., Gammie, C., & McKinney, J. C. 2009, ApJ, 701, 620
- Zhu et al. (2010) Zhu, Z., Hartmann, L., & Gammie, C. 2010, ApJ, 713, 1143