Solar System Evolution

Download as pdf or txt
Download as pdf or txt
You are on page 1of 57

Earth, Moon, and Planets (2006)  Springer 2006

DOI 10.1007/s11038-006-9087-5

3. Solar System Formation and Early Evolution:


the First 100 Million Years

THIERRY MONTMERLE and JEAN-CHARLES AUGEREAU


Laboratoire dAstrophysique de Grenoble, Universite Joseph Fourier, Grenoble, France
(E-mails: [email protected]; [email protected])
MARC CHAUSSIDON
Centre de Recherches Petrographiques et Geochimiques (CRPG), Nancy, France
(E-mail: [email protected])
MATTHIEU GOUNELLE1,2
1
Museum National dHistoire Naturelle, Paris, France; 2Natural History Museum, London, UK
(E-mail: [email protected])
BERNARD MARTY
Ecole Nationale Superieure de Geologie, Nancy, France
(E-mail: [email protected])
ALESSANDRO MORBIDELLI
Observatoire de la Cote dAzur, Nice, France
(E-mail: [email protected])

(Received 1 February 2006; Accepted 4 April 2006)

Abstract. The solar system, as we know it today, is about 4.5 billion years old. It is widely believed
that it was essentially completed 100 million years after the formation of the Sun, which itself took less
than 1 million years, although the exact chronology remains highly uncertain. For instance: which, of
the giant planets or the terrestrial planets, formed rst, and how? How did they acquire their mass?
What was the early evolution of the primitive solar nebula (solar nebula for short)? What is its
relation with the circumstellar disks that are ubiquitous around young low-mass stars today? Is it
possible to dene a time zero (t0), the epoch of the formation of the solar system? Is the solar system
exceptional or common? This astronomical chapter focuses on the early stages, which determine in large
part the subsequent evolution of the proto-solar system. This evolution is logarithmic, being very fast
initially, then gradually slowing down. The chapter is thus divided in three parts: (1) The rst million
years: the stellar era. The dominant phase is the formation of the Sun in a stellar cluster, via accretion of
material from a circumstellar disk, itself fed by a progressively vanishing circumstellar envelope. (2) The
rst 10 million years: the disk era. The dominant phase is the evolution and progressive disappearance
of circumstellar disks around evolved young stars; planets will start to form at this stage. Important
constraints on the solar nebula and on planet formation are drawn from the most primitive objects in
the solar system, i.e., meteorites. (3) The rst 100 million years: the telluric era. This phase is
dominated by terrestrial (rocky) planet formation and dierentiation, and the appearance of oceans and
atmospheres.

Keywords: Star formation: stellar clusters, circumstellar disks, circumstellar dust, jets and outows; solar
nebula: high-energy irradiation, meteorites, short-lived radionuclides, extinct radioactivities, supernovae;
planet formation: planetary embryos, runaway growth, giant planets, migration, asteroid belt, formation
of the Moon; early Earth: atmosphere, core dierentiation, magnetic eld
THIERRY MONTMERLE ET AL.

3.1. The First Million Years: The Stellar Era1

THIERRY MONTMERLE

3.1.1. THE SUNS BIRTHPLACE

Stars are not born in isolation, but in clusters. This is what astronomical
observations of our galaxy (the Milky Way) and other galaxies tell us. The
birthplace of stars are the so-called molecular clouds, i.e., vast, cold volumes
of gas (mostly molecular hydrogen and helium, and also complex organic
molecules, with so far up to 11 C atoms: see Ehrenfreund and Charnley, 2000
for a review). These clouds also contain dust grains (which include heavy
elements in the form of silicates, hydrocarbons, and various ices). The
masses of molecular clouds typically range from 106 to 108 Mx: in principle,
molecular clouds are suciently massive to form millions of stars. However,
somewhat paradoxically, molecular clouds do not naturally tend to form stars:
gravitation, which would tend to generate the free-fall collapse of molecular
clouds in less than 1 Myr,2 appears to be balanced by an internal source of
pressure which keeps them in gravitational equilibrium. The basic answer lies in
the study, in the radio range, of the velocity distribution in the gas. It can be
shown that this distribution corresponds to a state of turbulence, i.e., gaseous
eddies that exchange energy from the large scale (the size of the cloud) to the
small scales (cloudlets of size ~0.1 pc); smaller scales may be present but are
currently beyond the spatial capabilities of existing radiotelescopes. However, on
a large scale, it can be seen from mid- to far-IR observations (which have a better
spatial resolution than in the mm range), that molecular clouds are in fact
lamentary, but these laments are constantly moving, as attested by their
velocity distribution. Figure 3.1 shows a 100-micron image by the IRAS satellite
of the Orion molecular cloud complex, where dense and cold laments are
conspicuous. According to turbulence theories (and dedicated laboratory
experiments), energy is transferred from the large scales to the small scales. So the
question becomes: what drives the turbulence? Or, in other words, where does
the supporting energy come from? Current explanations are still debated. They
focus either on an external energy source like a neighboring supernova, or on an
internal feedback mechanism: as we shall see below (Section 3.1.2), young stars
drive powerful outows of matter (Reipurth and Bally, 2001), in such a way

1
How astronomers determine stellar ages is described in Chapter 2 on Chronometers
(Section 2.1).
2
Here we adopt the usual astronomical convention: 1 Myr = 106 years. This is exactly
synonymous to 1 Ma, as used for instance by geologists elsewhere in the article (where a
stands for annum).
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.1. The Orion complex. Left: image of the Orion nebula M42 in the visible domain (
Anglo-Australian Telescope). Background: far-IR image (100 microns) of the Orion complex,
by the IRAS satellite (1986), covering a very wide area (the angular scale is given). Note the
widespread lamentary structure of the giant molecular cloud. The bright spots are several
star-forming regions belonging to the same complex, the most active one being M42 (box).

that they inate the turbulence cells to keep the cloud from collapsing.
Nevertheless, energy is dissipated at the smallest scales, so that some form of
collapse is inevitable: the idea is that the smallest cloud structures, cloudlets
or prestellar cores, eventually collapse to form stars (e.g., Bate and Bonnell,
2004; Goodwin et al., 2004; Padoan et al., 2004).
Other arguments point to an important role of the interstellar magnetic eld.
In principle, a molecular cloud is by denition cold and entirely neutral, thus
cannot be inuenced by the presence of magnetic elds. But in practice, a
minute fraction of the gas (roughly 10)7) is ionized (electrically charged) by
ambient cosmic rays and also by hard radiation from young stars (UV and
X-rays, see Section 3.2.1.3). The charged particles are tied to the magnetic eld,
and thus, through collisions with them, neutral particles are in turn inuenced
by it: this is called ambipolar diusion. In a way, ambipolar diffusion acts as a
THIERRY MONTMERLE ET AL.

dragnet through which neutral particles ow across magnetic eld lines. This
effect is quantitatively important: measurements of magnetic eld intensity
inside molecular clouds (via the Zeeman effect on molecular lines) show that the
gas pressure and the magnetic pressure are just about equal, with a difference of
at most a factor of 2 in either direction, depending on the clouds (e.g., Padoan
et al., 2004; Crutcher, 2005). This means that in reality we are probably dealing
with magnetically regulated turbulence: the ow of gas in turbulent cells is not
free, but is slowed down by magnetic elds and preferentially proceeds along
laments (e.g., Pety and Falgarone, 2003; Falgarone et al., 2005). In particular,
this means that, at the small scales, gravitational collapse proceeds either on a
short, free-fall time scale (typically ~104 years) if the magnetic eld is on the
weak side (magnetic pressure < gas pressure), or on a long, ambipolar diusion
time scale (which can reach several 105 years or more) if the magnetic elds is
suciently strong (magnetic pressure > gas pressure).
This picture, at least qualitatively, leads to the idea that molecular
clouds are stable, self-supported structures, but on the verge of gravita-
tional collapse. Depending on the intensity of the magnetic eld, star
formation may occur at many places in the cloud, perhaps in sequence,
within a relatively short timescale, of the order of a few 104 years locally,
a few 105 years to 106 years globally. One such global theoretical mech-
anism is called competitive accretion: stars form out of shocks within a
pool of colliding gaseous laments, where they compete to acquire their
mass as they move through the cloud (Bate and Bonnell, 2004, Clark et al.
2005; Figure 3.2). Pure gravitational collapse has also been advocated
(Krumholz et al., 2005).
Whatever the details of the various star formation mechanisms, the net
result is star formation in clusters. Molecular observations (Motte et al.,
1998), in the mass range between ~0.1 and 1 Mx, show that above 0.5 Mx
the core mass distribution and the observed stellar initial mass function
(IMF) are the same, which strongly suggests (but does not prove) that the
stellar mass distribution directly derives from the core mass distribution,
itself linked with the turbulent structure of molecular clouds. (The IMF is the
distribution of stellar masses at formation: it is observed to be a universal
law, expressed as dN*/dlogM*  M*)1.5 for M* 0.5 Mx, where N* is the
number of stars in the mass range M*, M* + dM*; explaining it is one of the
hardest challenges for star formation theories: see Kroupa, 2002 for a
review.) Depending on a number of external conditions, such as the total
molecular cloud mass, the passage of a shock wave of a nearby supernova
explosion from the most massive stars (see below), etc., the high-mass end of
the IMF is observed to be cut-o at some value Mmax. Some clusters have
massive to very massive stars (Mmax up to several tens of Mx), others have
only intermediate-mass stars (Mmax = a few Mx at most), all the way to very
small masses (brown dwarfs, that are not massive enough to eventually
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.2. Numerical three-dimensional simulation of star formation in a 10,000 Mx cloud,


~600,000 yrs after the initial collapse (P.C. Clark, private communication; simulations done at
the UK Astrophysical Fluids Facility in Exeter). The gure is 5 pc on a side. Note the
similarity of the cloud structure with that of the Orion complex shown in the previous gure.
The simulation eventually leads to the formation of ~500 stars. This is less than observed in
Orion (~2000 stars in M42). Complicated eects such as feedback on turbulence from stellar
outows and ionizing radiation have not been included. A better agreement is expected in the
future when these eects are taken into account.

trigger nuclear reactions, M* < 0.08 Mx). Because of the observed univer-
sality of the IMF, a large Mmax implies a large number of stars, a small Mmax
implies a small number of stars. For example, star-forming regions like Orion
display stars up to 20 Mx or more, and contain altogether several thousand
stars, while others like Ophiuchus, Taurus, etc., do not go beyond a few Mx
and harbour only a few tens to a few hundred stars.
Once the stars are formed, what remains of the parental cloud, not yet
condensed into stars, is eventually dispersed, and the stars become optically
visible. At this point, the stellar cluster becomes free from its parent cloud,
and its evolution is regulated by dynamical eects in its own gravitational
potential, leading after a few tens of Myr to open clusters, then to a broad
dispersal of the stars in the galaxy (at typical velocities on the order of a few
km/s), much like a beehive, and thus to a loss of memory of how and
where they were formed individually.
The Sun probably has been one of such stars. The statistics of eld stars
(like the Sun today) vs. the number of stars in star-forming regions leads to a
probability argument drawn from observations: nearly 90% of solar-like stars
must have been born in clusters, of a few tens to a few thousand stars (Adams
and Myers, 2001).
THIERRY MONTMERLE ET AL.

It is therefore impossible at present to know a priori in which type of


stellar cluster the Sun was born: our star is about 50 times older than the
oldest open clusters. But the issue is fundamental for planet formation in
general and for the origin of the solar system in particular, because of the
short time evolution of the circumstellar disks around young stars (see Sec-
tion 3.2.1). To simplify, there are two extreme possibilities for the birthplace
of the Sun:
(i) The Sun was born in a rich, Orion-like environment. (see Hillen-
brand 1997; Figure 3.3) The most massive stars (the Trapezium-like
stars) are very hot, and thus emit UV photons able to strongly ionize
their immediate environment. The disks of the less massive stars then
tend also to be ionized and evaporate, as shown in Figure 3.4
According to most calculations, the disks, which have masses in the
range 10)2 10)4 MSun (i.e., 0.110 Jupiter masses), will disappear in
a few million years only, likely before any terrestrial planet has had
the time to form. One possibility for disks to survive is not to stay
too long in the vicinity of the hottest stars, so as to escape, via
dynamical eects, the original beehive. We observe that most disks
around young stars are typically 10 times larger that the size of the
present-day solar sytem; some disks may be cut-o to the size of the
Kuiper Belt (i.e., the radius of the solar sytem, ~50 AU),3 mainly
because of evaporation processes (Adams et al. 2004, 2006). On the
other hand, the discovery of the existence of Sedna, a relatively
massive, high-eccentricity solar system planet, implies that the
protosolar disk did not suer stellar encounters closer than 1,000 AU
(Morbidelli and Levison, 2004), so retained its initial large size for
some time, or else managed to be not entirely vaporized. This would
be possible in the outskirts of the nebula, or in a less rich cluster like
NGC1333 (Adams et al., 2006). In any case, disk survival in a Orion-
like environment is certainly not easy.4
(ii) The Sun was born in a poor, Ophiuchus-like environment (Fig-
ure 3.5). With no massive star around, stars form deeply embedded
inside the molecular cloud, and once formed they stay protected from
external disturbances. All disks survive, and can remain very large. We
know of some examples of stars with large disks in the immediate
vicinity of molecular clouds (Figure 3.6). It is not clear how in this case

3
1 UA = 1 Astronomical Unit = Sun-Earth distance = 150 million km.
4
To be complete, one should mention the recent work by Throop and Bally (2005), who argue
that dust grains actually grow into planetesimals under the coagulating eect of UV radiation,
hence that planet formation (see below, Section 3.2.1), is favored by evaporation. But even in
this case, the exposure to UV radiation must be ne-tuned for the whole system to survive.
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.3. Near-IR (2 microns) image of the center of M42, revealing the stars of the rich
Orion nebula cluster, and in particular the four central hot stars called the Trapezium,
which excite the nebula. The nebula is 450 pc away. The image is about 10 on a side, which is
250 000 AU. ( ESO, VLT-ISAAC, by M. McCaughrean.).

Figure 3.4. A tear in Orion. This is an evaporating circumstellar disk, 500 AU in diameter.
The central star is clearly visible. The bright spot is oriented towards h1 Ori C, the hottest star
of the Trapezium. The evaporating gas is shaped by the wind from this star ( NASA Hubble
Space Telescope: J. Bally, H. Throop, & C.R. ODell).
THIERRY MONTMERLE ET AL.

Figure 3.5. Left: The poor q Oph cluster in the optical range, hidden in a dark cloud of gas
and dust. It is located 150 pc from the Sun. (The bright stars at the bottom left is the
foreground star Antares, and the fuzzy spot to its right is the distant globular cluster M15. (
Anglo-Australian Telescope.) Right: Mid-IR image of the cloud core, revealing embedded
young stars and protostars invisible in the optical range (ISOCAM, Abergel et al. 1996).

planet-forming disks would shrink by a factor of 10 to be conned


within the size of the solar system. In the end, perhaps the issue of the
existence of a planetary system such as ours is that of the survival of
the ttest: to suer disk truncation in a dense, Orion-like cluster, or at
least to avoid complete evaporation by way of dynamical eects, i.e. to
be ejected from the original cluster quickly enough. The quiet
environment of the Ophiuchus-like, looser clusters, would provide a
safe evolution, but, likely, at the cost of carrying along a large,
massive disk which would lead to planetary systems very dierent from
our own. In this scenario, the solar system would be born in rare,
though not uncommon (10% of the stars), environment, with a key
role played by infrequent dynamical interactions, cutting o its ori-
ginal disk, very early (105 years?) after its birth.

3.1.2. THE SUN AS A FORMING STAR

Let us now zoom on the Sun as a forming star, which we shall assume
isolated for simplicity (knowing from above that it has to be isolated, at some
very early stage, from its cousins in a cluster of forming stars).
The knowledge of what we believe must have been the rst stages of
formation and evolution of the Sun, is drawn from a wealth of observations
of a multitude of star-forming regions that astronomers have been able to
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.6. A lonely, massive egde-on disk in the ourskirts of the q Oph cloud (circle). The
other disk-like object (dotted circle) is a distant galaxy ( ESO, VLT-ISAAC, Grosso et al.,
2003).

obtain. Nowadays, telescopes are being used both on the ground and in
space, covering almost all the electromagnetic spectrum, from mm
wavelengths, to X-rays, across the IR and optical domains. Depending on the
wavelength, it is possible to pierce the darkness of molecular clouds, and
see inside them to watch the hidden birth of solar-like stars, most
importantly in the IR to mm domains, and in the X-ray and gamma-ray
ranges (e.g., Ryter, 1996).
In almost every case, one is able to distinguish three main components,
which simultaneously evolve as star formation proceeds5 (see, e.g., Shu et al.,
1987; Andre and Montmerle, 1994; Andre et al., 2000, a summary and recent
references are given in Feigelson and Montmerle, 1999, and Montmerle,
2005: Figure 3.7). At the so-called protostellar stage, a vast, dense enve-
lope (1,00010,000 AU in radius) is detectable, and from the center emerges a
bipolar outow. The envelope is so dense that its interior is invisible even
at mm wavelengths; only its outer structure can be seen. It is now understood
that the seed of a new star is formed from matter accreted from the
envelope which rains on it under the pull of gravitation. The youngest
observed protostars have an estimated age of ~104 years: this estimate is

5
For a pioneering work, establishing (analytically!) the basic principles of early stellar evo-
lution, see Hayashi (1966).
THIERRY MONTMERLE ET AL.

Figure 3.7. A summary table of the various protostellar and stellar phases, with characteristic
timescales and basic observational properties. (From Feigelson and Montmerle, 1999).

rather uncertain, but is consistent with the number deduced from the
dynamical age of outows (= size/velocity), and from their small number
relative to their more evolved counterparts, like T Tauri stars (see below):
indeed, one nds roughly 1 protostar for every 100 T Tauri stars, aged
110 Myr.
At an age ~105 years, the envelope is much less dense, since most of it
has collapsed onto the disk. It becomes transparent at mm wavelengths,
revealing a dense disk (5001,000 AU in radius), from which the seed star
continues to grow. The source of outows become visible, in the form of
highly collimated jets originating close to the central star, conrming
earlier models in which molecular ourows consist of cold cloud material
entrained by the jet. (In fact, this jet-cloud interaction is believed by some
authors to be the main agent to sustain the turbulent state of the cloud:
this is the feedback mechanism mentioned above, see Matzner and
McKee, 2000.) At this stage is revealed the key three-component structure
that governs the physics of star formation: an outer envelope, an inner
accretion disk, and matter ejected perpendicular to the disk (Figure 3.8).
(The accretion disk probably exists from the start of collapse, but it
cannot be detected because of the opacity of the envelope at the earliest
stage.)
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.8. Sketch of the structure of a protostar, zooming on the star-dik interaction region,
which is dominated by magnetic elds. This region is the seat of the accretionejection
mechanism, by which the majority of the disk mass becomes accreted to form a star at the
center, while the remainder is ejected. (From Feigelson and Montmerle, 1999).

As time passes, this three-component structure considerably evolves. The


envelope is eventually exhausted after ~106 years, leaving a fully developed
real, luminous star, a massive circumstellar disk (Mdisk ~110 MJup), and a
weak, optically visible bipolar jet. This is the start of the so-called classical
T Tauri stage (after the name of the rst-discovered star of this type, which,
as it turns out, is quite atypical for its class; e.g., Bertout, 1989). Then, after a
period which can be as long as 107 years, the disks of T Tauri stars (and jets)
disappear, or at least becomes undetectable, presumably via planet for-
mation (see Section 3.2.1): this is the weak T Tauri stage. Why weak?
Because the classical T Tauri stars have a very unusual optical spectrum,
with very strong emission lines, in particular the Ha line of ionized hydrogen;
in contrast, the weak T Tauri stars have weak emission lines, compa-
rable to solar lines. A weak Ha line can be explained, in analogy with the
Sun, by the presence of active regions, that is, starspots scattered over the
stellar surface, indicative, again as on the Sun, of locally strong magnetic
elds. The strong Ha line of classical T Tauri stars cannot be explained by
magnetic activity alone. Various arguments assign the strength of the Ha line
in this case to matter falling onto the star, fed by the disk: this is the phe-
nomenon known as accretion of matter (e.g., Bertout, 1989) the very
process by which the star grows to reach its nal mass when the disk is
exhausted, perhaps leaving young planetary bodies behind.
Let us now concentrate on young classical T Tauri stars, like HH30 in
Taurus (Figure 3.9). All the observational evidence points to a causal rela-
tionship between the existence of jets and the presence of disks: it is clear, at
least qualitatively, that the jet material is somehow coming from the disk, at
the same time that accretion feeds the central star. This phenomenon,
known as the accretionejection phenomenon, is absolutely central to our
THIERRY MONTMERLE ET AL.

Figure 3.9. Left: HST image of the T Tauri star HH30, showing its edge-on disk and jet. The
young Sun could have been such an object. Right: sketch of the theoretical magnetic structure
used to model the accretionejection mechanism (the background drawing is taken from
Ferreira et al., 2000). The star is an image of the magnetically active Sun, seen in X-rays by
the Yohkoh satellite (see Section 3.2.1.2).

understanding of star formation and early evolution, and to a modern view


of the solar nebula. Such a phenomenon is somewhat paradoxical: it implies
that, for a star to form, it must lose mass! At least, a signicant fraction of
the mass (>10% from observations, e.g., Muzerolle et al., 2001) must
eventually be ejected. The point is that mass accumulation is not the only
necessity to form a star: since (at we see at all scales in the universe, from
clusters of galaxies to planets) rotation is always present, and seen in the form
of circumstellar disks around forming stars, there must exist centrifugal
forces that oppose gravity. For gravity to win and lead to star formation,
angular momentum must be lost. Although the reasons for angular momen-
tum loss in a forming star, and in particular the role of the circumstellar disk,
are not entirely clear because of complex transport processes within them (see
next subsection), it is well known that mass loss in the form of stellar winds
(the solar wind is one example) is very efcient to spin down a star provided
the ejected matter remains coupled to the star. The only way to do it is to link
the star and the wind by a magnetic eld. Along this line of thought, in
current models, accretion (mass gain), and ejection (mass and angular
momentum loss), must be mediated by magnetic elds.

3.1.3. A STELLAR VIEW OF THE PRIMITIVE SOLAR NEBULA

To be more specic, let us now zoom again, this time on the region very close to
the forming star, for instance a young T Tauri star, to which the primitive solar
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

nebula must have been comparable in its rst million years (Figure 3.9, left). A
rich interplay between observations and theory, over the last decade, results in
the following picture (e.g., Matsumoto et al., 2000; Shang et al., 2002; Ferreira
and Casse, 2004, and refs therein; see Figure 3.9, right). Both the star and the
disk are magnetized: (i) the star is surrounded by a dipolar magnetosphere
that surrounds it like a tire, with a closed, loop-like topology of the magnetic
eld; (ii) in contrast, the magnetic eld lines connected to the disk are open,
above and below the disk. Then a special distance, called the corotation ra-
dius Rc, is naturally dened: this is the distance at which the Keplerian (i.e.,
orbital) velocity of a disk particle rotates at exactly the same speed as the star.
(One could say, in analogy with the Earths articial satellites, that this is the
astrostationary orbit.) At distances r < Rc from the star, the intensity of the
magnetic eld is stronger than at Rc, and the magnetosphere rotates in a rigid
fashion, the eld lines being anchored on the stellar surface. At distances
r > Rc, the magnetic eld decreases rapidly and takes an open, spiral form as it
becomes tied to the disk. The point at Rc thus has a very particular magnetic
property: it is the border at which the magnetic eld topology switches from
closed (stellar component) to open (disk component). As such, it is also known
as the X-point (or more exactly the X-ring in three dimensions in view of
the assumed axial symmetry) because of its X-shaped magnetic conguration
(Shu et al. 1997).
The existence of the X-point (in a 2-dimensional cut) holds the key to the
majority of accretionejection theories. There are many discussions among
theorists about its exact status. For instance, it is not clear why the mag-
netically dened X-point (at a distance Rx which depends only on the mag-
netic eld intensity) should be exactly at the same location as the
gravitationally dened corotation radius Rc (which depends only on the
stellar mass and rotation velocity), in other words why should Rx = Rc. It is
not clear either that the X-point should be that: a point (or more precisely an
X-ring, in three dimensions), since this would mean that at Rx there must
be an innite concentration of magnetic elds lines, which is physically
impossible if matter is coupled to it (via some ionization, for instance due to
X-rays), etc. But most theorists (and observers) agree, at least qualitatively,
on the following general accretionejection picture, which will be sucient
for the purpose of this paper.
Because of its special gravitational and magnetic properties, seen in two
dimensions, the X-point de facto behaves like a Lagrange point (Figure 3.9,
right): (i) If a particle, initially located at this point, is pushed towards the
interior (r < Rx), it will start falling freely on the star under the pull of
gravity, along the corresponding rigid magnetic eld line: this is called
magnetospheric accretion. (ii) Conversely, if a particle is pushed outwards,
it will start following an open eld line, and the centrifugal force will push it
even further: this is how centrifugal jets are formed. Thus, the X-point is
THIERRY MONTMERLE ET AL.

intrinsically unstable, analogous to the gravitational Lagrange point L1 be-


tween two celestial bodies (which here would be a three-dimensional mag-
netic Lagrange ring in a star-disk sytem), and owing through it matter
can either fall onto the star or be ejected. In practice, of course, matter does
both: depending on the exact model, calculations predict that 1/10 to 1/3 of
the disk matter should be ejected, the remainder falling on the star and thus
feeding its growth (e.g., Johns-Krull and Gaord, 2002). What is important
for our purpose is the location of the X-point: for a solar-mass star rotating
relatively slowly as T Tauri stars are observed to rotate (period of a few
days), Rx is on the order of several stellar radii. Since a typical T Tauri star
has a radius R* ~ 3 Rx, this means Rx ~ 0.1 AU. In the well-studied T Tauri
star AA Tau, for which the accretion disk is seen edge-on, it has been possible
to reconstruct the magnetic structure of the inner hole, from a study of its
active regions eclipsed by the disk (Bouvier et al., 1999, 2003). Although in
this case the disk is found to be warped, which implies that the magnetic
axis is tilted to the rotation axis, the size of the resulting magnetosphere is
in good agreement with the theoretical picture (Figure 3.10), in particular

Figure 3.10. Sketch of the tilted magnetic structure inside the disk of the T Tauri star AA Tau.
This sketch results from a study of the magnetic activity of the star, eclipsed by its warped
accretion disk seen nearly edge-on (From Bouvier et al., 1999).
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

with the presence of an inner hole ~0.1 AU in radius. As a matter of fact,


the magnetic eld in the inner regions of several young stars has recently
been directly detected, via the Zeeman eect (Donati et al., 2005).
Compared to the present-day solar system, such a hole is well within the
orbit of Mercury (0.4 AU). This is also the size of a pinhole compared to
the disk sizes, which, as mentioned above, can be as large as 1,000 AU at
an early stage.
Yet this pinhole region, which harbors the magnetic central engine
for the accretionejection mechanism, may have played an important role for
the early solar system. Indeed, because it is a region of tangled, unstable
magnetic elds, and located close to the star, itself the seat of its own intense
magnetic activity (as testied by the observed aring X-ray emission, see
Section 3.2.1.2), it is a place where any circumstellar material (gas, grain,
possible small planetary bodies, etc.) will suer a high dose of radiation,
either in the form of hard photons (XUV activity), or energetic particles
(accelerated in stellar ares like on the Sun, and which, while not directly
observable, can be induced from X-rays). We shall return below (Section
3.2.1.2) to this important question, which connects the early evolution of
circumstellar disks with the history of the young solar system.
To summarize, at an age of ~106 years, the young solar system or the
solar nebula can be depicted as follows. At its center, a magnetized solar-
like star has reached most of its nal mass, ~1 Mx This star is surrounded
by a rigid magnetic cavity (its magnetosphere), about 0.1 AU in radius.
At the X-point, corotating with the star, lies the inner edge of the cir-
cumstellar disk, from which most (0.92/3) of the matter continues to fall
on the star. The remainder is ejected on both sides perpendicular to the
disk, in the form of a powerful, highly collimated jet, which evacuates
most, if not all, its angular momentum. Beyond the X-point lies a dense
circumstellar disk, either small (~50 AU, say) if it has been truncated in
the vicinity of hot stars (as in a bright, UV-rich Orion-like environment),
or large, a few hundred AU in radius if it was born in a dark, quiet
Ophiuchus-like environment. Within this disk, matter continues to ow
from the outside to the inside, braked by viscosity, until the disk is
exhausted or because matter starts to assemble to form large grains and,
ultimately, planetesimals and giant planets, as discussed in the next section
below.
In this dynamic picture of the solar nebula, any heavy particle like a
grain can have two fates when it arrives in the vicinity of the X-point: either it
falls onto the star by accretion, or it is entrained outwards by light par-
ticles (the gas), but eventually falls far back onto the disk in a ballistic fashion
because it is too heavy to be carried away along the open magnetic eld lines
by the centrifugal force. In this last case, it will have spent some time very
close to the X-ring, and will have suered there heavy irradiation by hard
THIERRY MONTMERLE ET AL.

photons and energetic particles. These particles will then mix to the disk,
holding specic scars from their passage near the X-ring. As explained in
Section 3.2.2, this is how some models explain the mysterious presence of
extinct radioactivities in meteorites.

3.2. The First 10 Million Years: the Disk Era

3.2.1. THE EVOLUTION OF CIRCUMSTELLAR DISKS AROUND YOUNG STARS AND


IMPLICATIONS FOR THE EARLY SOLAR SYSTEM

3.2.1.1. The path to planets: Astronomical timescale for the growth of dust
grains
JEAN-CHARLES AUGEREAU
Now that a dense circumstellar disk is installed around the central star, it
must evolve: on the one hand, it continues (albeit at a lesser rate) to lose mass
at its inner edge (by way of magnetic accretion), on the other hand, grains
assemble via low relative velocity collisions to form larger, preplanetary
bodies. But how long does it last? Infrared observations of T Tauri stars,
which are sensitive to the presence of circumstellar material, show that disks
disappear on widely dierent timescales. Figure 3.11 (Hillenbrand, 2006),
which collects data from ~30 star-forming regions, shows that so-called
inner disks (i.e., regions warm enough to radiate are near-IR wavelengths)
are ubiquitous at young ages, and tend statistically to disappear after a few
million years only. Quantitatively, the fraction is consistent with 100% at an
age 1 Myr, and drops to less than 10% after 10 Myr, with some clusters
containing no disk at all. Actually, this low fraction of old disks is a lower
limit: mid-IR observations, which are sensitive to cooler disks, hence more
distant from the central star, show that in some cases, like the nearby g Cha
cluster, aged 9 Myr, the fraction of disks is closer to 4060% (resp. Megeath
et al., 2005; Lyo et al. 2003), suggesting that disks may live longer than
previously thought. However, at this stage the disk mass is found to be too
low to form even a Jupiter (10)3 Mx)or perhaps they have already done
soso that the general conclusion is that giant planets, if any, must have
formed on timescales signicantly shorter than 10 Myr.
Therefore, the disk era is a critical period for planet formation, tightly
constrained by astronomical observations. Submicron-sized dust particles
composing young and massive disks constitute the raw material from
which planets form. Tiny dust grains must coagulate to form large dust
aggregates, pebbles, and then larger rocky bodies (planetesimals) before the
dust disk becomes too tenuous. The formation of giant planets through the
core-accretion scenario also requires the formation of planetary cores
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.11. Disk fraction in young stellar clusters, as a function of their age. This fraction is
consistent with 100% at young ages (less than 1 Myr), then declines over a timescale of a few
Myr. After 10 Myr, with a few exceptions, all the disks around young stars have disappeared,
presumably because of planet formation (Hillenbrand, 2006). This puts strong constraints on
the formation of giant planets (which cannot be seen themselves).

before the disk has been mostly depleted of gas. A detailed investigation of
how dust grains grow into large planetary embryos is presently one of the
most important open questions, especially for the formation of the solar
system: this is discussed in detail in Section 3.2.4. Here we give a broad
outline of the results drawn from the study of circumstellar disks around
young stars.
A key conclusion is that the planet formation process is observationally
required to be both fast and common. The disappearance of circumstellar
disks in less than 510 Myr is actually interpreted as a direct consequence of
the formation of larger solid bodies decreasing the opacity and the dust
emission (Haisch et al., 2001; Carpenter et al., 2005; Hillenbrand, 2006).
Unless the circumsolar disk survived dissipation processes longer than usu-
ally observed for circumstellar disks, large solid bodies in the solar system
should have then formed within less than a few million years, which is a
major challenge for terrestrial planet formation theories.
THIERRY MONTMERLE ET AL.

In addition, the growth of solid particles in disks must be suciently


generic for at least two reasons. First, a signicant fraction of solar-like stars
have been found during the last decade to host giant gaseous exoplanets6.
The formation of our solar system may then not be peculiar, and the general
conclusions derived from statistical analysis of star-forming regions may also
apply to the early solar system. Nevertheless, it is not known whether these
giant exoplanets formed through the accretion of gas onto a solid core or not.
But there exists an independent evidence for planetesimal formation to occur
routinely in disks around young stars. A large fraction of main sequence stars
in the solar neighborhood are observed to be surrounded by tenuous disks
composed of short-lived dust grains. The survival of these so-called debris
disks, which contain much less mass that young disks (typically 10)6 Mx),
over hundreds of million years points indirectly towards the presence of
reservoirs of meter-sized (or larger) bodies, similar to the Kuiper Belt in our
solar system, which release grains by mutual collisions. (We see this phe-
nomenon continuing in the present-day solar system in the form of zodiacal
light.) The observation of debris disks is up to now one of the most con-
vincing observational clues indicating that planetesimal formation in young
disks is common.
Various processes contribute to the coagulation and growth of dust grains
in disks of suciently high density. Current theoretical models of circum-
stellar disks (Section 3.2.4 and references therein) include Brownian motion,
vertical settling, radial drift and turbulence. Interstellar like grains (smaller
than about 0.1 micron) stay well mixed with the gas where the gas density is
high enough. In that case, the low-velocity collisions between grains due to
Brownian motion result in aggregates with uy structures. This process is
particularly ecient in the disk mid-plane while collisions at the disk surface
could be destructive. As aggregates grow and reach millimeter sizes, they
settle to the disk mid-plane on short time-scales depending somewhat on the
strength of the turbulence.
Can we test this model? The available tools to witness grain growth in
disks remain limited. Observations can only probe the continuum, quasi-
blackbody emission of grains with sizes of the order of the observing
wavelength. Therefore, astronomers are basically limited to the direct
detection of solid particles smaller than, at most, a few centimeters. Fortu-
nately, disks become optically thin at millimeter wavelengths, which allows
astronomers to probe the denser regions of disks where sand-like grains are
expected to reside. At these wavelengths, the total disk ux directly relies on
the dust opacity that is a function of the grain size, and it can be shown that
grains several orders of magnitude larger than those found in the interstellar

6
Visit The Extrasolar Planets Encyclopaedia: http://www.exoplanets.eu
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

medium are required to explain the (sub-)millimeter observations (Draine,


2006). At shorter wavelengths, the disk mid-plane is opaque and only the
upper layers, mostly consisting of small grains according to the models, are
accessible. But as the opacity of disks inversely depends on the wavelength,
observations should probe deeper and deeper regions as the wavelengths
increases (Duchene et al., 2004). In practice, the interpretation of scattered
light observations of disks at visible and near-infrared wavelengths is not
straightforward as it strongly depends on the light scattering properties of the
individual grains that are, unfortunately, hard to model properly. This type
of study is also limited to a handful of objects, such as the circumbinary disk
of GG Tau (Figure 3.12)
More statistically meaningful conclusions can be derived from the study of
dust mineralogy through infrared spectroscopy. Especially, silicates that are
composed of silicon, oxygen and, in most cases, a cation (often iron or
magnesium), have characteristic solid-state vibrational bands that occur
uniquely at infrared wavelengths. The emission spectra of silicates from
optically thin disk surfaces display characteristic emission features around 10
microns due to the SiO stretching mode, and around 20 microns due to the
OSiO bending modes. Silicates are observed in the interstellar medium and
in many solar system objects, including the Earth mantle where it appears

Figure 3.12. Left panel: observations of the dust ring about the GG Tau young binary system.
In order to properly interpret the almost wavelength-independent appearance of the ring, dust
settling toward the disk midplane must be taking into account in the models (middle and right
panels). Such images may be unveiling the vertical stratied structure of the disk (Duchene
et al., 2004). The cold gas component of the GG Tau ring world has also been observed in
the mm domain (Guilloteau et al., 1999).
THIERRY MONTMERLE ET AL.

mostly in the form of olivine (after its olive-green color). Their presence in
disks around young solar-like stars is thus expected, by analogy with the
solar system. But it is only recently that large fractions of T Tauri stars could
be spectrally studied thanks to highly sensitive infrared space telescopes such
as Spitzer (Figure 3.13). Silicates turn out to be ubiquitous in T Tauri disks in
Myr-old forming regions. But, interestingly, the strength and shape of their
emission features dier in many cases from those observed in the interstellar
medium, indicating signicant processing of silicates in young disks. As the
silicate features strongly depend on grain size, the presence in disk atmo-
spheres of dust particles several orders of magnitude larger than interstellar
grains provides a natural explanation to the observations. Silicate emission
features can thus be used as extremely valuable diagnostics of micron-sized
solid particles in disks, as demonstrated for instance by Kessler-Silacci et al.
(2006). Moreover, as the stars are fairly young (less than a few million years),
this indicates fast grain growth in disks.

Figure 3.13. Ten micron silicate features from dust disks around stars of various masses
(Natta et al. 2006). From left to right: Herbig Ae stars (HAe), a few times more massive than
the Sun), T Tauri stars (TTS), and Brown Dwarfs (BDs). Examples of theoretical silicate
emission features, based on their experimental optical properties, are dispayed in the extreme
right panel (Lab.Sil.). By comparing the shape and strength of these proles to the observa-
tions, one can infer the presence of large grains as well as the degree of crystallinity of the
silicates.
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

The detection in circumstellar disks of crystalline silicates similar to those


observed in solar-system cometary dust is an additional diagnostic of similar
signicant grain processing within the rst 10 Myr. The silicates that are
injected into the primordial stellar nebulae are in the form of amorphous
silicates, as proven by the recurrent non-detections of crystalline silicates in
the interstellar medium. The annealing (crystallization) of amorphous sili-
cates can only happen in the inner regions of the disks where the temperature
is suciently high. Also their incorporation in long-orbital period solar
system objects points towards radial (and vertical) mixing of dust during the
infancy of the solar system. Various mechanisms throughout the disk such as
stellar/disk winds or radial and vertical mixing in (magneto-)turbulent disks
can, in principle, transport dust grains with important implications on the
disk chemistry and mineralogy. The observation of disks around young so-
lar-like stars supports the radial mixing scenario since crystalline silicates are
found at distances where the temperature is too low for in situ annealing. In
contrast, it has been suggested that supernova shock waves in the outer solar
nebula could have rapidly annealed the amorphous silicates in situ prior to
their incorporation into comets, thereby eliminating the need for large-scale
nebular transport processes (Harker and Desch, 2002; see also below, Section
3.2.2.3). But it is currently impossible, even for the solar system, to distin-
guish between the two proposed scenarios.
It is particularly interesting to note that the statistical analysis of large
samples of T Tauri stars shows a poor correlation between the grain size in
the upper layers of disks and the age of the star. The same conclusion can be
drawn when the degree of crystallinity of the silicate grains is considered.
Although star ages are dicult to estimate accurately,7 these results tend to
indicate that age is not the only parameter controlling disk evolution. The
environment (ambient radiation eld, bound companions, etc.) is likely to
impact the disk evolution and hence the planet formation which could
explain star-to-star variations at similar ages. Although the formation of
meter-sized (or larger) rocks is a prerequisite to form terrestrial planets and
planetary cores, one should keep in mind that this step is one of the less
observationally nor theoretically constrained in the planet formation process.
Exoplanet embryos are nevertheless expected to leave characteristic imprints
in disks such as gaps or even to clean up the inner disk regions. The detection
of gaps in young systems due to planet-disk coupling is still pending. In a few
cases, the presence of depleted regions inside disks can be indirectly inferred
from the lack of signicant circumstellar emission at short IR wavelengths
indicating a low amount of warm material close to the star. But such objects
remain exceptional. Thus, at present, astronomers are only able to draft

7
See Chapter 2
THIERRY MONTMERLE ET AL.

general trends without being able yet to further constrain the main steps that
allow to go from sub-micron size grain to km-sized bodies, and this is a major
problem in particular for theories of the formation of the solar system
(Section 3.2.4).

3.2.1.2. X-ray induced irradiation phenomena in circumstellar disks


THIERRY MONTMERLE
A new factor probably plays an important role, for star formation as well as
for planet formation. This is the ubiquitous X-ray emission from young stars,
from the protostar stage (see, e.g., Feigelson and Montmerle, 1999; Micela
and Favata, 2005, for reviews). X-rays are mainly detected in the form of
powerful ares, very similar to those seen in the images of the Sun which the
Yohkoh satellite sent us every day from 1991 to 2001,8 but much more intense
(103105 times stronger than for the quiet Sun).
In the course of these ares, the stellar X-ray luminosity LX amounts to
10 to 10)3 times the total luminosity of the young star (called the bolo-
)4

metric luminosity Lbol), reaching a few percent or more in exceptional cases.


It can be shown that the plasma (a very hot, ionized gas, with a temper-
ature of 107 to 108 K) that emits the X-rays is conned in very large magnetic
loops (up to 23 R, i.e., about 10 Rx for a T Tauri star). By analogy with
the Sun, it is thought that the fast heating results from so-called recon-
nection events, in which magnetic eld lines of opposite polarities get in
contact in a short-circuit, and suddenly release the magnetic energy pre-
viously stored in the course of various motions, for instance when magnetic
footpoints are dragged by convective cells (e.g., Hayashi et al. 2000). The
plasma then cools radiatively typically in a few hours by the emission of
X-ray photons.
Thanks to X-ray satellites (from the 90s), X-ray emission has been detected
from hundreds of T Tauri stars, either concentrated in young stellar clusters
like q Ophiuchi (e.g., Ozawa et al., 2005) or Orion (Getman et al., 2005), or
else more widely dispersed as in the Taurus-Auriga clouds (e.g., Stelzer and
Neuhauser, 2001), and many other star-forming regions. For a given popu-
lation of T Tauri stars (with and without disks), the X-ray detection rate is
consistent with 100%. (Only the fainter stars are not all detected, for lack of
sucient sensitivity.) In other words, observations suggest that all young,
solar-like stars, emit X-rays, at a level LX/Lbol ~ 10)4 to 10)3 with a fairly
large dispersion (see, e.g., Ozawa et al., 2005), to be compared with LX/
Lbol ~ 10)7 for the present-day Sun. This is for instance most vividly illus-
trated by the ~1,500 sources detected in the 10-day Chandra exposure of the
Orion Nebula (Getman et al., 2005; Feigelson et al., 2005) (Figure 3.14).

8
Visit http://www.solar.physics.montana.edu/YPOP/
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.14. The center of the Orion nebula cluster. Left: near-IR image (this is the same as
Figure 3.3). Right: corresponding X-ray image by Chandra (Getman et al., 2005). Note the
excellent identication between the IR and X-ray sources, demonstrating that all young stars
emit X-rays at levels much higher than the Sun.

In addition to providing us information on the magnetic state of young


stars (for instance relevant to ejection phenomena, or to the reconnection
mechanisms, as discussed above), are X-rays are important on another
ground: they induce irradiation eects on the dense surrounding circum-
stellar and interstellar material: ionization and energetic particle interactions
(this last point will be discussed in the next subsection; see Glassgold et al.,
2000, 2005; Feigelson, 2005).
The rst eect of X-ray irradiation is the ionization of the accretion disk.
Direct evidence for irradiation has been found recently in the form of a
uorescence line of neutral iron at 6.4 keV in the spectra of nearly a dozen of
T Tauri stars with disks, mainly in the q Oph and Orion clusters (Favata
et al., 2005; Tsujimoto et al. 2005), and also by the presence of specic
molecules in some disks (e.g., Greaves, 2005; Ilgner and Nelson, 2006, and
references therein). Such a characteristic uorescence line can be excited only
in response to X-ray irradiation of the cold gas. On theoretical grounds,
several authors have studied the case of an accretion disk irradiated by X-
rays (e.g., Glassgold et al., 2000, 2005; Fromang et al., 2002; Matsumura and
Pudritz, 2006). The results depends on details of the adopted disk model, but
the main conclusion, sketched in Figure 3.15, is that the ionization of the
accretion disk is dominated by X-rays, except in the densest equatorial
regions where they cannot penetrate (dead zone), i.e., within a few AU of
the central source. Everywhere else, the disk is partially ionized, and the
ionization fraction xe = ne/nH is roughly comparable to that of the inter-
stellar medium (xe ~ 10)7 or less).
The ionization state of the disk has another important consequence: the
coupling of matter with magnetic elds. Indeed, even very weakly ionized
THIERRY MONTMERLE ET AL.

Figure 3.15. Sketch of the X-ray irradiation of circumstellar disks. Note in particular the
emission of a neutral X-ray uorescence line at 6.4 keV, detected in several systems. The other
line, at 6.7 keV, is characteristic of a hot, X-ray emitting plasma, at temperatures of several
106 K. Note also the dense, neutral dead zone, which, according to some authors, would be
the most favorable to planet formation.

matter sticks to the magnetic eld via collisions between neutral atoms and
charged atoms (ions): this process, called ambipolar diusion, has already
been mentioned above in the context of star formation in molecular clouds
(Section 3.1.1). When this process is coupled with the Keplerian motion of all
atoms in a disk around the central object, it gives rise to a so-called mag-
netorotational instability, discovered by S. Chandrasekhar and extensively
studied by Balbus and Hawley (1991). This instability is invoked to explain the
strong viscosity of accretion disks, and thus regulates accretion itself (e.g.,
Fromang et al., 2004). One interesting consequence is that such strong viscous
coupling would not exist in the (neutral) dead zone, which would then un-
dergo accretion only if some other, non-magnetic mechanism for an ecient
viscosity were at work (which cannot be excluded). In fact, it has even been
suggested that such a protected, neutral region, of size ~20 AU, would be
favorable for planet formation (Glassgold et al., 2000; Matsumura and
Pudritz, 2006). Note, however, that the very existence of an extended dead
zone has been challenged by recent numerical computations, which take into
account turbulent motions within the disk, and which show that the neutral,
dead zone volume tends to mix with the surrounding, ionized material: in the
end, (weak) ionization probably dominates everywhere in the disk (Fleming
and Stone, 2003), although this study does not take grains into account, which
turn out to play an important role in the disk ionization (Ilgner and Nelson,
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

2006). Up to now, such a weak, but widespread ionization, which must also
aect the growth of dust grains since they would be electrically charged, has
not been taken into account in planet formation theories.

3.2.2. THE FIRST FEW MILLION YEARS AS RECORDED BY METEORITE DATA

MARC CHAUSSIDON, MATTHIEU GOUNELLE, THIERRY MONTMERLE

We now turn to the earliest stages of the solar system: the solar nebula, and
the so-called meteoritic record. While most of the extraterrestrial ux to the
Earth is in the form of micrometer size dust (20,000 tons of the so-called
micrometeorites par year), there are about 10 tons of meteorites of centi-
meter to meter size that fall on Earth every year. Most of these meteorites are
likely to come from the asteroidal belt situated between Mars and Jupiter.
Among meteorites, chondrites are primitive objects that have endured no
planetary differenciation. This is attested by:
(i) their mineralogical composition which reects accretion of compo-
nents formed at high temperature (the chondrules) and low temper-
ature (the matrix), components that were never homogenized
chemically and/or isotopically through melting and metamorphism,
(ii) their bulk composition which is similar to the photospheric abundances
of the elements in the Sun (Anders and Grevesse, 1989) and is thus
considered as primitive, i.e. reecting that of the forming solar system,
(iii) their age : they are the oldest rocks of the solar system, with a Rb/Sr age
of 4.55 Gyr (Wasserburg, 1987). (This early pioneering result is now
updated by more precise measurements: see below, Section 3.2.2.2.)

Other meteorites such as irons or some achondrites are the by-products of


planetary dierentiation : they have younger ages (by a few Myr) and have
compositions reecting the interplay of silicatemetal dierentiation (core-
mantle formation) and of silicatesilicate dierentiation (crust-mantle for-
mation). Thus chondrites, which are undierentiated, can help us unravel the
physico-chemical conditions in the solar nebula and the astrophysical context
of the Suns birth, while dierentiated meteorites trace the formation and
evolution of large planetary bodies.

3.2.2.1. Chondrite components and physical conditions in the solar nebula


MARC CHAUSSIDON, MATTHIEU GOUNELLE

Among chondrites, carbonaceous chondrites have a chemical composition


most similar to that of the Sun, and are therefore believed to be our best
proxy for the protosolar nebula (Brearley and Jones, 1998). The mineralogy,
THIERRY MONTMERLE ET AL.

chemistry and isotopic composition of carbonaceous chondrites can help


decipher the formation and history of the solar nebula. Carbonaceous
chondrites are made of CalciumAluminium-rich Inclusions (CAIs),
chondrules and matrix, chondrules representing by far the major component,
7080% in volume (see recent reviews by Zanda, 2004 for chondrules and by
MacPherson and Huss, 2003, for CAIs). (Figure 3.16)
Matrix is rich in volatile elements (e.g., H2O, C,...) : it is a ne-grained
component made of chondrule fragments and of minerals stable at low
temperature. Matrix has endured extensive secondary processing in the
parent-body evidenced by metamorphic transformations (solid state diu-
sion) and hydrothermalism due to uid-rock interactions. Matrix is therefore
not very useful for pinpointing the physical conditions in the solar nebula. At
variance, CAIs and chondrules are high temperature components which were
formed in the solar nebula and thus predate parent-body processes. CAIs are
refractory components made of Al-, Ca-, Ti-rich silicates and oxides.

Figure 3.16. Fragment of the Allende meteorite, revealing (large white area in the bottom
right-hand corner) the so-called CalciumAluminium-rich Inclusions (CAIs), in which
evidence for short-lived, extinct radioactivities have been found.
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Accordingly, they are generally considered from their mineralogy to be the


rst solids formed by condensation from the solar nebula gas at temperatures
higher than 1800 K. Chondrules are spherical objects, made of ironmagne-
sium silicates (mostly olivine and pyroxene), metal and suldes. Chondrules
generally contain a large fraction of glassy mesostasis, implying that they were
once molten and were subsequently quenched at a few to perhaps 1000 K/h
(Hewins, 1997). Experimental studies show that CAIs cooled at 0.1 to
10 K/h (Stolper and Paque, 1986). Both CAIs and chondrules are thought
to have formed in the solar nebula via complex high-temperature processes
including condensation, evaporation, melting, etc.... Despite extensive studies,
the exact mechanisms that led to the formation of CAIs and chondrules are
still elusive, but it is clear that CAIs and chondrules formed at dierent times
or locations. Chondrules probably formed at a higher pressure than CAIs and
it is generally considered that chondrules formed later than CAIs, although
this assumption cannot yet be demonstrated by direct dating of each com-
ponent. In fact, the formation of chondrules and CAIs may overlap and there
exists a variety of chondrules and CAIs which may have formed in dierent
environments and may have complicated histories, involving precursors of
variable composition and more or less extensive exchanges with the nebular
gas (see Hewins et al., 2005, and references therein).

3.2.2.2. Duration of the solar nebula


MARC CHAUSSIDON, MATTHIEU GOUNELLE

Radioactive nuclide abundances in rocks are classically used to date rocks


and to infer timescales for geological processes.9 Among these, the U/Pb
system which combines two parent/daughter couples (238U decays to 206Pb
with a half life of 4.47 Gyr, and 235U decays to 207Pb with a half life of
0.7 Gyr), is the one which has provided the most accurate dating of mete-
oritic components. Accordingly, the CAIs appear to be the oldest component
within chondrites, with a U/Pb age of 4567.2 0.6 Myr (Manhes et al.,
1988; Allegre et al., 1995; Amelin et al., 2002). Slightly dierent values, based
on various methods, have appeared in the recent literature, giving an age of
4568.5 0.4 Myr (Baker, 2005; Bouvier et al., 2005), but for the purpose of
this article these values (which dier by ~1.3 0.5 Myr) are consistent within
the uncertainties. This age can be considered as giving the time zero, t0, for
the start of the formation of the solar system; however this current experi-
mental time resolution of ~1 Myr does not allow to build a precise chro-
nology of the earliest processes which occurred in the solar accretion disk and
which gave birth to the rst solar system solids from the nebula. In fact, the

9
See Chapter 2 on Chronometers (Section 2.2).
THIERRY MONTMERLE ET AL.

typical duration of these processes is most likely less than 1 Myr, so revealing
them requires further improvements in the experimental methods.
Of particular interest for cosmochemistry are the so-called extinct radio-
active nuclides, or short-lived radioactive nuclides, which have half-lives be-
low a few Myr. The ones which have been identied in meteorites up to now
are 7Be (T1/2 = 53 days), 41Ca (T1/2 = 0.1 Myr), 36Cl (T1/2 = 0.3 Myr),
26
Al (T1/2 = 0.74 Myr), 10Be (T1/2 = 1.5 Myr), 60Fe (T1/2 = 1.5 Myr) and
53
Mn (T1/2 = 3.7 Myr). The presence of excesses of their radioactive
daughter elements (e.g., positive correlation between 26Mg/24Mg and
27
Al/24Mg ratios for 26Al, see Figure 3.17) shows that these radioactive
nuclides were present in various amounts both in CAIs and in chondrules
when they formed. Because of their short half-lives, short-lived radioactivities
can be used to constrain very tightly timescales (McKeegan and Davis, 2003
and refs therein): for instance the 26Al/27Al ratio decreases by a factor of two
within 0.74 Myr. In this respect, 10Be, 26Al and 60Fe are of special interest
since they can (i) help constrain the chronology of the rst million years of
evolution of the solar system, and (ii) give nuclear clues to the astrophysical
context of the Suns birth.
Many measurements of the initial content of 26Al in CAIs have led to the
idea of a canonical ratio, 26Al/27Al = 4.5 10)5, that would have dened
the starting time t0 of the protoplanetary solar system (MacPherson et al.,
1995)which must not be confused with the starting time of the formation of
the Sun as a star, which in this context precedes t0 (see Section 2.5). CAIs are
particularly well suited for the determination of the 26Al/27Al ratios, because
since they are refractory objects they are enriched in Al relative to Mg. Less
numerous measurements in chondrules, which are less enriched in Al relative
to Mg and thus more dicult to analyze, suggest that chondrules formed
with 26Al/27Al < 1 10)5 (Mostefaoui et al., 2002). Assuming an homoge-
neous distribution of 26Al in the solar nebula, this yields a ~23 Myr age
dierence between CAIs and chondrules. This timescale has long been con-
sidered as the characteristic timescale of the solar nebula as determined by
chondrite data. If true, such a long duration for the high temperature pro-
cesses in the accretion disk implies that CAIs and chondrules must have been
in some way stored in the nebula for 23 Myr before being accreted to-
gether to form the chondrites.
At this point it is important to stress that this so-called chronological
interpretation of the 26Al/27Al variations relies entirely on the assumption
mentioned above, which is very strong: the existence once in the early solar
system of an homogeneous distribution of 26Al (Gounelle and Russell, 2005).
This assumption is far from being proven and very dicult to demonstrate in
fact, simply because one would need very precise independent absolute ages
for dierent objects to be able to test the homogeneity of their 26Al/27Al
ratios. It could in fact well be that the dierence in the initial content of 26Al
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.17. Variations in CAIs and chondrules of the Mg isotopic compositions (given as
26
Mg/24Mg in the upper panel, or as the permil 26Mg excesses noted d26Mg* in the bottom
panel) due to the radioactive decay of short-lived 26Al (which decays to 26Mg with a half-life of
0.7 Myr). These data show that 26Al was present in the early solar system and that its
abundance (given as 26Al/27Al ratios) can be used to constrain the chronology of the for-
mation of CAIs and chondrules. A canonical 26Al/27Al ratio of 4.5 10)5 has been found by
in situ ion microprobe analysis in most CAIs (data from Podosek et al., 1991, in the upper
panel). Supracanonical 26Al/27Al ratios recently found in CAIs are shown by a grey eld (data
from Young et al., 2005 and Bizzarro et al., 2004) in the bottom panel. In this panel, the
details of the distribution of 26Mg excesses in chondrules (data from Galy et al., 2002; Bizzarro
et al., 2004; Chaussidon et al., 2006) is not yet well understood : it can be interpreted as
reecting either the formation of some chondrules very early (i.e. at the same time than CAIs)
or a late formation of chondrules in the nebula from a mixture of precursors including CAI
material.
between chondrules and CAIs is due in part to spatial heterogeneities in the
solar nebula. Predictions can be made on the existence of such an hetero-
geneity depending on the nucleosynthetic origin of 26Al, either last minute
injection from a nearby massive star or supernova, or production by irra-
diation processes around the early Sun (see next subsection). Note also that
THIERRY MONTMERLE ET AL.

the picture has recently been complicated by results from the Oxford-UCLA
laboratory (Galy et al., 2002; Young et al., 2005) who showed that the
canonical ratio of CAIs might be due to a resetting event of an initially higher
ratio of 7 10)5, and that chondrules might have formed with an initial
26
Al/27Al ratio signicantly higher than 1 10)5. The high 26Al/27Al ratios
inferred for some chondrules can indicate that they formed very early (Biz-
zarro et al., 2004) or be simply due to the fact that they contain an inherited
CAI component (Galy et al., 2002).
Recently, high precision data on Pb isotopes (U/Pb dating system) have
been obtained on CAIs and chondrules (Amelin et al., 2002), providing
absolute ages of the chondrites components. CAIs have ages ranging from
4567.4 1.1 Myr to 4567.17 0.70 Myr, while chondrules have ages rang-
ing from 4564.66 0.3 Myr to 4566.7 1 Myr. This conrms that the solar
nebula could have lasted for several million years, and question the idea that
all chondrules formed 23 Myr after CAIs. Most recent Pb data seem to
indicate that a range of ages do exist for chondrules, from nearly as early as
CAIs to a few Myr later.

3.2.2.3. Extinct radioactivities and the astrophysical context of the birth


of the Sun
MARC CHAUSSIDON, THIERRY MONTMERLE

The presence of short-lived radioactivities in the solar nebula is intriguing.


Because their half-life is shorter than the timescales needed to isolate the future
solar system material from the interstellar medium, they need to have a last
minute origin, as opposed to the steady-state abundance of galactic cosmic-
ray induced isotopes. They could have been injected in the nascent solar
system by a nearby supernova, or made by irradiation of the nebular dust and/
or gas by energetic solar protons, or both. The challenge of any mechanism is
to explain most, if not all, of the extinct radioactivities at the same time.
The ubiquitous, intense and are-like X-ray activity of young stars (Sec-
tion 3.2.1.2) support the idea that irradiation could have been an important
process in the early solar system (see previous sections; also Chaussidon and
Gounelle, 2006, for a review of the traces of early solar system irradiation in
meteoritic components, and Feigelson, 2005 for the relation with X-ray ares
from young stars). Such an origin is supported by the presence of 10Be
(T1/2 = 1.5 Myr) and of 7Be (T1/2 = 53 days) in Allende CAIs (McKeegan
et al., 2000; Chaussidon et al., 2006), beryllium isotopes being made by are-
induced energetic particle irradiation (see Figure 3.18).
On the other hand, the presence of 60Fe in suldes and silicates from
chondrites (Huss and Tachibana, 2004; Mostefaoui et al., 2004) implies the
presence of a nearby supernova, since 60Fe is too neutron-rich to have been
made by in-ight spallation reactions (Lee et al., 1998), and is a signature of
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.18. Variations in one Allende CAI of the 10B/11B ratios vs. the 9Be/11B ratio (data
from Chaussidon et al., 2006). The isochrone-type correlation demonstrates that short-lived
10
Be (which decays to 10B with T1/2=1.5 Myr) was present in the early solar system. The same
CAI contains traces of the in situ decay of short-lived 7Be (T1/2=53 days). Radioactive 10Be
and 7Be can be produced in the early solar system by irradiation processes around the Young
Sun (Chaussidon and Gounelle, 2006).

explosive nucleosynthesis in supernovae. The presence of a nearby supernova


would imply in turn that the Sun was born in a crowded stellar environment,
possibly similar to HII regions as observed in Orion (Hester and Desh, 2005,
and above, Section 3.1.1). 26Al also can be made in massive stars and
supernova explosions (though other stellar processes, such as in novae and
red giant stars, which are old, are also important sources, see Diehl et al.
2006), and it is thus tempting to explain the abundance of both isotopes in
this way. As a matter of fact, one can nd a set of parameters (time interval
between the explosion and the birth of the solar system Dt = 1 Myr, dilution
factor due to transport of supernova material to the presolar core f=10)4)
that explains both abundances, and also other isotopes (Busso et al., 2003;
Gallino et al., 2004). But this may not be completely realistic: as a rule, most
star-forming regions are in some vicinity of OB associations in which massive
stars evolve and explode sequentially. If they explode at a frequency shorter
than ~1 Myr, the resulting 60Fe and 26Al nucleosynthesis will tend to be in a
steady state locally (as long as the massive star formation episode lasts),
rather than in separate bursts, subsequently decaying. For instance, in spite
THIERRY MONTMERLE ET AL.

of the fact that the Orion Nebula Cluster currently contains no supernova,
the 1.809 MeV gamma-rays from the decay of 26Al have been detected by the
GRO satellite, as a signature of past supernovae from massive stars exca-
vating the nearby so-called Eridanus superbubble (Diehl et al., 2004). Present
supernovae, on the other hand, have been found by the RHESSI and
INTEGRAL satellites, by way of the 1.173 and 1.333 MeV emission of 60Fe
in the Cygnus region, which hosts the most massive stars in the Galaxy
(Harris et al. 2005). On the other hand, for the presolar core the values of Dt
and f are loosely constrained since some arbitrary mass cut in the supernova
ejecta has to be invoked for the calculated abundances to be in agreement
with observations in meteorites (Meyer, 2005), and neglects the preceding
contribution of the precursor massive stars. It is thus troubling that the
irradiation model, which takes into account the enhanced stellar energetic
particle ux deduced from X-ray observations, should be able to account for
the 26Al/27Al ratio in CAIs independently from the possible existence of a
supernova.
At this point, it should be noted that, as early as 25 years ago, Montmerle
(1979) identied about 30 massive star-forming regions which were, based on
various observational criteria, tentatively associated with supernova (SN)
remnants. These special regions, dubbed SNOBs (for OB associations
or molecular clouds observed to be associated with supernovae) were at the
time searched in connection with the identication of high-energy gamma-ray
sources. For our purpose, this sample can be taken as examples of the reality
of supernovae exploding in the close vicinity of young stars. It also shows
that only a small fraction of all OB associations are in this situation at any
given time. An illustration of the complexity of the problem is given, again,
by the Orion star-forming region.
In 1895, E.E. Barnard discovered a faint, almost exactly (half-) circular
ring in the outskirts of Orion, spanning several degrees in the sky. This
spectacular structure, now known as Barnards Loop, is shown in Fig-
ure 3.19. It is readily visible in the Ha line of ionized hydrogen. Its exact
nature is still uncertain: proposed to be a supernova remnant by Opik back in
1935, recent radio observations (Heiles et al., 2000) have shown that the
emission is thermal, i.e., radiated only by an ionized gas, just like HII regions
around young stars. The idea then is that this ring has been created by winds
from the Belt stars (not the Trapezium; the Belt stars are the conspicuous
stars visible with the naked eye on clear winter nights, see Figure 3.19). But it
is also possible that the traditional radio signature of supernova remnants,
i.e., the nonthermal emission from high-energy electrons accelerated at the
shock front, is somehow buried in a more intense thermal emission a
classical dilemma in supernova remnant identications. In either case,
however, kinematics indicate that Barnards Loop must have originated
about 3 106 years ago close to the Orion nebula.
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.19. Barnards Loop surrounding the Orion nebula, seen in Ha. It is unclear whether
this extended structure is a the remnant of a supernova explosion or an ionized shell created by
stellar winds from the Belt stars, but kinematical studies give it an age of 3 106 years.
(Photograph by E. Mallart).

When considering the consequences on the possible irradiation of the


proplyds of the Orion nebula (see Figures 3.3 and 3.4), one is therefore
faced with two possibilities: (i) either Barnards Loop is not the remnant of a
supernova explosion having taken place 3 106 years ago, in which case
none of the current young stars in Orion have been peppered with 60Fe, (ii) or
it is, then only a fraction of them have been. This can be seen by looking at
the HR diagram presented in Chapter 2 on chronometers (Figure 2.3). In
this gure, we can easily see that the majority of young stars present in the
Orion Nebula Cluster are in fact younger that 3 106 years, so that any 60Fe
spread by the explosion has disappeared, and these young stars cannot be
contaminated by 60Fe. Quantitatively, one can nd after Figure 3.20 (taken
from Palla and Stahler, 1999), that less than 40% of the older generation of
young suns may have been eectively contaminated by the SN explosion, if
THIERRY MONTMERLE ET AL.

Figure 3.20. Age distribution in four mass ranges for stars in the Orion Nebula Cluster (Palla
and Stahler, 1999). The dotted line indicates the estimated age of Barnards Loop, 3 Myr. For
stars of mass close to 1 Mx, less than 40% are older than 3 Myr and may have been
contaminated by short-lived nucleosynthetic products of a supernova explosion, provided
Barnards Loop is a supernova remnant, which is still unclear.

there was such an explosion. The presence of 60Fe in the early solar system
would then not be the rule, even in an Orion-like birthplace.
In the end, and as discussed in the introductory sections, the birthplace of
the Sun is still an unresolved question, although the birth of the Sun in a rich
cluster seems to be favored by stellar statistics. To understand whether the Sun
was born in a high-mass environment like Orion, or in a low-mass environ-
ment like q Ophiuchi or Taurus has however important implications not only
for the astrophysical conditions for the Suns birth itself, but also for the
chronology of the early solar system. Indeed, depending on their origin, short-
lived radionuclides were or were not homogeneously distributed in the solar
nebula. Usually short-lived radionuclides are expected to be homogeneously
distributed if coming from a supernova and heterogeneously distributed if
originating from in situ irradiation by energetic particles. This comes from the
fact that supernova material must be largely volatilized in the HII region and
homogeneously mixed in the accreting disk where it is injected, though there is
at present no denitive observation demonstrating whether a fraction of the
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

supernova ejecta cannot be in the form of solid grains (Hester and Desch,
2005). At variance, irradiation models based on X-ray are observations and
the X-ring picture (Gounelle et al., 2001; 2004; see above, Section 1.3)
predict possible variations in the production rate of short-lived radionuclides
depending on parameters such as the uence of the accelerated particles, their
composition and the composition of the irradiated target. Irradiation, how-
ever, might also produce rather constant radionuclide abundances if charac-
teristic time scales and compositions are considered for the irradiation. It is
obvious that identifying the source of 26Al and other short lived nuclides
would provide a long-needed basis for the chronological use of these elements.
One could think of solving this problem by looking at young stars with
circumstellar disks, in which particle irradiation is actually taking place as
seen in X-rays, making use of the fact that, as mentioned above, 26Al decays
by the emission of 1.8 MeV gamma-rays, and is thus observable elsewhere
than in the solar system by gamma-ray telescopes such as GRO. The cal-
culation has been done by Montmerle (2002): taking into account that
gamma-ray telescopes have a very wide eld-of-view (several degrees in
diameter), when pointed at a star-forming region they will integrate the ux
of a whole star-forming region, i.e., several hundred young stars at the same
time. As it turns out, even under the most optimistic assumptions the
1.8 MeV ux is undetectable, and dominated by the general 26Al emission in
the Galaxy, which is most conspicuous in massive star-forming regions be-
cause of successive supernova explosions (see Diehl et al. 2006).
In summary, while the presence of 60Fe in CAIs shows that the forming
solar system has been at least polluted by a nearby supernova, it is not clear
whether this supernova has been responsible for the other extinct radioac-
tivities, including 26Al.

3.2.3. INTERMEZZO

THIERRY MONTMERLE

The study of circumstellar disks around solar-like stars, as described in


Section 3.1.2, obviously has strong implications on our current views on the
origin and formation of the solar system. We observe their global evolution
(timescale for the disappearance of disks in stellar clusters), and we begin to
understand their physicochemical evolution (growth and mineralogy of dust
grains via IR spectroscopy) over a few million years, while large quantities of
gas (detectable via mm observations) are still present. However, as already
mentioned, once the grains reach sizes of a few mm, they are not observable
any more. They start to be observable again much later (several tens of
million years), but only as debris of collisions between macroscopic bodies,
presumably planetesimals or asteroids, giving indirect evidence for ongoing
THIERRY MONTMERLE ET AL.

planetary formation. Thus, there is a crucial observational time gap between


the primordial disks (which would correspond to the solar nebula), and
second-generation debris disks (in which planetary formation is well under
way). Translated into sizes, this means that we have essentially no observa-
tional constraints on the transition between cm-sized grains (direct) and
km-sized bodies (indirect).
In principle, the ~200 known exoplanetary systems should also give us
clues about planetary formation in general, and the formation of the solar
system in particular. However, at least as far as the formation of the solar
system is concerned (which is our main concern here in the context of the
origin of life as we know it), there are still many open problems. Indeed, while
current observing techniques are sensitive mostly to super-giant exoplanets
close to the central stars (the so-called hot Jupiters, with masses between 5
and 1015 Jupiter masses or more), and increasingly more to lower-mass
exoplanets, exoplanetary systems with several planets, when they exist, do
not resemble at all the solar system. The recently announced discovery of a
distant 5-Earth-mass exoplanet orbiting a very low-mass, non-solar, star, by
way of gravitational lensing towards the galactic center (Beaulieu et al. 2006),
is not helpful either in this context, since nothing is known about the possible
exoplanetary system it belongs to.
In the end, for most of the phase of planet formation we are essentially left
with theory, backed by rare and dicult laboratory experiments. The goal is
to put together mechanisms for grain agglomeration and destruction, and
mutual dynamical interactions between bodies of various sizes and masses in
the general gravitational eld of the central star, and also in the presence of
gas which aects the motion of small particles. In spite of the relative sim-
plicity of the basic equations, this involves sophisticated numerical N-body
simulations. The following section describes the basic theoretical concepts of
planetary formation restricted to the solar system, with the goal of ultimately
understanding the various steps that led to the formation of our planet, the
Earth.

3.2.4. THE FIRST STAGES OF PLANETARY FORMATION IN THE SOLAR SYSTEM

ALESSANDRO MORBIDELLI

3.2.4.1. From micron-size to kilometer-size bodies (1 Myr)


In the proto-solar nebula, the most refractory materials condense rst,
gradually followed while the local temperature drops by more and more
volatile elements. The dust grains thus formed oat in the gas. Collisions
stick the grains together, forming fractal aggregates, possibly helped
by electrostatic and magnetic forces. Other collisions then rearrange the
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

aggregates, and compact them. When the grains reach a size of about a
centimeter, they begin to rapidly sediment onto the median plane of the disk
in a time
Tsed  R=qp Xa  qv=qp X2 a

where R is the nebula surface density, q is its volume density, v is the r.m.s.
thermal excitation velocity of gas molecules, qp the volume density of the
particles, a is the radius, and X is the local orbital speed of the gas (Goldreich
and Ward, 1973; Weidenschilling, 1980). Assuming Hayashis (1981) minimal
nebula (minimal mass solar composition nebula, with surface density R
proportional to r)3/2, containing materials to create the planets as we know
them), one gets
Tsed  103 =a years
where a is given in cm. This timescale, however, is computed assuming a
quiet, laminar nebula. If the nebula is strongly turbulent, or strongly per-
turbed from the outside or by the ejection of jets in the proximity of the star,
the sedimentation can become much longer.
Once the dust has sedimented on the mid-plane of the nebula, the clock
measuring the timescale of planet accretion starts eectively to tick. Thus, the
time t0 for planetary accretion is not the time t0 usually used in stellar for-
mation theory (start of the collapse of molecular cores) or the time t0 of
cosmochemists (the formation of the rst CAIs, see above, Section 3.2.2.2).
Linking the various times t0 together is one of the major problems in
establishing an absolute chronology for the formation of the solar system. In
addition, notice that Tsed above depends on the heliocentric distance r. This
means that time t0 is dierent from place to place in the disk!
The growth from dust grains to kilometer-size planetesimals is still
unexplained. There are two serious issues that remain unsolved. The rst
issue concerns the physics of collision between such bodies at speed of order
10 m/s (typical collision velocity for Keplerian orbits with eccentricity
e ~ 10)3), which is still poorly understood. For dust grains, current theories
predict that collisions are disruptive at such speeds (Chokshi et al., 1993).
However, recent laboratory experiments on collisions between micrometer
size grains (Poppe and Blum, 1997; Poppe et al., 2000) give a critical velocity
for accretion (velocity below which grains accrete, and above which they
fragment) 10 times larger than predicted by the theory. The reason is prob-
ably that the fractal structure of the dust allows to absorb energy much better
than envisioned a priori. This may help solving the fragmentation paradox
for dustdust collisions. When the agglomeration of dust builds larger bod-
ies, though, the problem of collisional disruption becomes much more severe.
Laboratory experiments cannot be done at this size range and one has to
THIERRY MONTMERLE ET AL.

trust computer models. Specic simulations of this process with SPH


(Smooth Particle Hydrodynamics) techniques have been done by Benz (2000)
for basalt bodies (monolithic or rubble piles) of sizes in the range m to km.
He found that low velocity collisions (540 m/s) are, for equal incoming
kinetic energy per gram of target material, considerably more ecient in
destroying and dispersing bodies than their high velocity counterparts.
Furthermore, planetesimals modeled as rubble piles are found to be char-
acterized by a disruption threshold about ve times smaller than solid bodies.
Thus, unless accretion can proceed avoiding collisions between bodies of
similar masses, the relative weakness of bodies in meter to km size range
creates a serious bottleneck for planetesimal growth. These apparently neg-
ative results, however, may once again depend on our poor understanding on
the internal structure of these primitive small planetesimals. The simulations
assume rocky objects, but it is still unclear how a puy pile of dust becomes a
solid rock. If the planetesimals were not yet solid rocks maybe the impact
energy could be dissipated more eciently.
The second open problem on planetesimal growth concerns radial
migration. Gas drag makes them fall onto the central star: gas being sus-
tained by its own pressure, it behaves as if it felt a central star of lower mass,
and therefore rotates more slowly than a purely Keplerian orbit at the same
heliocentric distance. Solid particles tend to be on Keplerian orbits and
therefore have a larger speed than the gas. So the gas exerts a force (drag) on
the particles, the importance of which is given by the characteristic stopping
time:
Ts mDV=FD ;
where m is the mass of the grain, DV the difference between the particle
velocity and the gas velocity, and FD the gas drag. For small particles, the
Epstein gas drag gives
Ts  qp =R;

while for big particles, the Stokes drag gives


Ts  qp a2 :

(Weideschilling, 1977). We then compare this time with the characteristic


Keplerian time
TK 1=X:
If Ts>>TK, the particle is almost decoupled from the gas. If Ts<<TK, it is
strongly coupled to the gas, and it tends to move with the gas ow. The
maximum eect occurs when Ts ~ TK, which occurs for meter size
particles. For the Hayashi minimal nebula, at about r = 1 AU from the Sun,
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

DV ~ 5055 m/s, and the particles radial velocity induced by the gas drag is
then of order 10100 m/s (Weideschilling, 1977). So meter-sized particles
should fall on the Sun in ~1001,000 years, i.e., before they can grow massive
enough to decouple from the gas.
One way out of this paradox is to have a density of solids larger than that
of the gas. In this case, the growth timescale would be faster than the radial
drift timescale. However, such a composition is not supported by observa-
tions nor by current theories.
Another possibility is the existence of vortices in the protoplanetary disk,
due to the turbulent viscosity of the nebula (Tanga et al., 1996). In this
model, 7090% of the particles are trapped in anti-cyclonic vortices. Once
trapped, the particles do not fall any more towards the Sun, but rather fall
toward the center of the vortex. Such falling timescale varies from a few
tens to a few thousand years, depending mainly on the size of the particle
and on the heliocentric distance (which increases all dynamical timescales).
Once at the center of the vortex, particles dynamics are stable over the
lifetime of the vortex. Their relative velocities are reduced (particles tend to
follow the gas stream lines, so they all tend to have the same velocity) and
the local density is increased, enhancing the accretion process. Therefore,
vortices would help the accretion of kilometer-size planetesimals in two
ways: by stopping the drift of meter-sized bodies towards the Sun, and by
speeding up the accretion process due to the accumulation of the bodies at
the centers of the vortices.
Which of these two possible situations is the real one, profoundly aects
the formation timescale. If there is no way to slow down the fall of growing
planetesimals towards the Sun, then the formation of a multi-km object has
to occur in about a few 1,000 years, probably by gravitational instability. If,
on the contrary, the turbulence of the disks is an eective obstacle to the
inwards drift, then the formation of planetesimals can take much longer. To
add confusion (reality is never easy), it is likely that while the rst plane-
tesimals are building up, new dust is settling on the midplane or drift in from
further heliocentric distances. So, dierent planetesimals can see dierent
times t0. In other words, even if the formation of a planetesimal is locally
very fast, the formation of a population of planetesimals can be ongoing for
much longer.

3.2.4.2. Formation of planetary embryos (110 Myr)


One way or another, we now have a gas disk containing kilometer-size
planetesimals. The dynamics of accretion starts to be dominated by the eect
of the gravitational attraction among the planetesimals, which increases the
collisional cross-sections. A runaway growth phase starts, during which the
big bodies grow faster than the small ones, hence increasing their relative
THIERRY MONTMERLE ET AL.

difference in mass (Greenberg et al., 1978). This process can be summarized


by the equation:
   
d M1 M1 1 dM1 1 dM2
 >0;
dt M2 M2 M1 dt M2 dt
where M1 and M2 are, respectively, the characteristic masses of the big and
of the small bodies, and can be explained as follows.
Generally speaking, accretion is favored by a high collision rate, which
occurs when the relative velocities are large, but also by large collisional
cross-sections and gentle impacts, which occur when the relative velocities are
low. Therefore the relative velocities between the dierent planetesimal
populations govern the growth regime.
At the beginning of the runaway growth phase the large planetesimals
represent only a small fraction of the total mass. Hence the dynamics is
governed by the small bodies, in the sense that the relative velocities among
the bodies is of order the escape velocity of the small bodies Vesc(2). This
velocity is independent of the mass M1 of the big bodies and is smaller than
the escape velocity of the large bodies Vesc(1). For a given body, the collisional
cross-section is enhanced with respect to the geometrical cross-section by the
so-called gravitational focusing factor:
Fg 1 V2esc =V2rel ;
where Vesc is the bodys escape velocity and Vrel is the relative velocity of the
other particles in its environment. Because Vrel ~ Vesc(2), the gravitational
focusing factor of the small bodies [Vsec = Vvesc(2)] is of order unity, while
that of the large bodies [Vesc = Vesc(1)>>Vesc(2)] is much larger. In this
situation one can show that mass growth of a big body is described by the
equation
 
1 dM1 1=3
 M1 V2 rel
M1 dt
(Ida and Makino, 1993) Therefore, the relative growth rate is an increasing
function of the bodys mass, which is the condition for the runaway growth
(Figure 3.21).
The runaway growth stops when the mass of the large bodies becomes
important (Ida and Makino, 1993) and the latter start to govern the
dynamics. The condition for this to occur is:
n1 M21 > n2 M22 ;
where n1 (resp. n2) is the number of big bodies (resp. small bodies). In this
case, Vrel ~ Vesc(1) ~ M11/3, and hence (1/M1)(dM1dt) ~ M11/3. The growing
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.21. A simulation of the runaway growth process for planetary embryos. In a disk of
equal mass planetesimals, two seeds (planetesimals of slightly larger size) are embedded. As
time passes, the two seeds grow in mass much faster than the other planetesimals,, becoming
planetary embryos (the size of each dot is proportional to its mass). While the growing
planetary embryos keep quasi-circular orbits, the remaining planetesimals have their eccen-
tricities (and inclinations) excited by the close encounters with the embryos. Notice also that
the separation between the embryos slowly grows in time (i.e. passing from one panel to the
subsequent one). From Kokubo and Ida (1998).

rate of the embryos gets slower and slower as the bodies grow, and the
relative dierences in mass among the embryos also slowly become smaller.
In principle, one could expect the small bodies themselves to grow, narrowing
their mass dierence with the embryos. But in reality, the now large relative
velocities prevent the small bodies to accrete with each other. The small
THIERRY MONTMERLE ET AL.

bodies can only participate to the growth of the embryos: this phase is called
oligarchic growth.
The runaway growth phase happens throughout the disk, with timescales
that depend on the local dynamical time (Keplerian time) and on the local
density of available solid material. This density will also determine the
maximum size of the embryos and/or planets when the runaway growth ends
(Lissauer, 1987). Assuming a reasonable surface density of solid materials,
the runaway growth process forms planetary embryos of Lunar to Martian
mass at 1 AU in 105106 yr, separated by a few 10)2 AU. Beyond the so-
called snow line at about 4 AU, where condensation of water ice occurred
because of the low temperature, enhancing the surface density of solid
material, runaway growth could produce embryos as large as several Earth
mass in a few million years (Thommes et al., 2003).

3.2.4.3. Formation of the giant planets (10 Myr)


Observations and models of the interior of the giant planets give some
important constraints on the composition and mass of the giant planets (see
Guillot, 1999, and for a review, Guillot, 2005):
(i) Jupiter has a mass of 314 ME (Earth mass), and contains ~1030 ME
of heavy elements;
(ii) Saturn has a mass of 94 ME and contains ~1020 ME of heavy ele-
ments;
(iii) Uranus and Neptune have a mass of 14 and 17 ME respectively, of
which only ~12 ME of hydrogen and helium.

To account for these constraints, the best current models for the formation
of Jupiter and Saturn assume a three-stage formation (Pollack et al., 1996):
(1) the solid core accretes as explained in the previous section; beyond the
snow line, the surface density of solid material is enhanced by a factor of
several, due to the presence of ice grains. This allows embryos to grow to
about 10 ME on a timescale of a million years (Thommes et al., 2003).
(2) The accretion of the solid core slows down (see above), while a slow
accretion of nebular gas begins, due to the gravity of the core. The gas
accretion continues at a roughly constant rate over many million years,
until a total mass of 2030 ME is reached;
(3) When the mass of the protoplanet reaches ~2030 ME the gas gravita-
tionally collapses onto the planet. The mass of the planet grows expo-
nentially and reaches hundreds of Earth masses in ~10,000 years: this is
the runaway phase. (Figure 3.22)

The model explains well the properties of Jupiter and Saturn, in particular
the existence of solid cores of about 515 Earth masses. There are however
four main problems in the above scenario, which have not yet been solved.
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Figure 3.22. The growth of a Jupiter-mass planet. The solid curve gives the mass of metals as a
function of time. The dotted curve gives the mass of the gas, and the dash-dotted curve the
sum of the two, as a function of time. Notice that the growth of the solid core of the
planet almost stalls after the rst 0.5 Myr. During a ~7 Myr timespan, the planet slowly
acquires an atmosphere, and, only when the total mass overcomes a critical threshold, a nal
exponential accretion of gas is possible. The timescale characterizing the slow accretion of the
gas depends on the opacity of the atmosphere. From Pollack et al. (1996).

(I) When the planetary core reaches a mass of several Earth masses, its
tidal interaction with the gas disk forces it to migrate very rapidly towards
the Sun. (This is called Type-I migration). The estimated falling time is
much shorter than the time required for the onset of the exponential accre-
tion of the massive atmosphere. Thus, giant planets should not exist! Two
ways out of this paradox have been proposed.
The rst is that the gas disk was violently turbulent. In this case, the
planetary cores would have suered a random walk, rather than a monotonic
infall towards the star (Nelson, 2005). Some would have collided with the star
even faster than in the absence of turbulence, but others could be lucky
enough to avoid collisions for a time long enough to start phase 3 above.
A second possibility (Masset et al., 2006) is that the gas disk surface
density had a radial discontinuity. For instance, the inner part of the disk
could be depleted by a factor of a few by the ejection of material in the polar
jets, typical of young magnetically active stars (Section 3.1.3). If such a
discontinuity exists, the planetary core would migrate towards the disconti-
nuity, and stop there until the atmosphere is accreted.
(II) The second phase of giant planet accretion (the slow accretion of the
atmosphere, prior to the onset of the runaway growth of the giant planets
mass) is also a problem, as it takes about 1015 Myr, longer then the typical
THIERRY MONTMERLE ET AL.

nebula dissipation time (Haisch et al., 2001; Hillenbrand 2006) (see above).
To shorten the timescale of the second phase, two solutions have been pro-
posed: to have an enveloppe of reduced opacity (Podolack 2003), or the
migration of the planet, which continuously feeds the growth of its core
(Alibert et al., 2005).
(III) The end of the exponential gas accretion is not yet fully understood.
Most likely the growth of the giant planets is slowed down when a gap is
opened in the gaseous disk, and is nally stopped when the nebula is dissi-
pated by photo-evaporation from the central star. If this is true, then the full
formation timescale of Jupiter and Saturn is of the order of the lifetime of the
gas disk, namely of a few Myr. For Uranus and Neptune, it is generally
assumed that the nebula disappeared before that the third phase of accretion
could start. This would explain why these two planets accreted only a few
Earth masses of gas.
(IV) The last problem is that, at the end of stage (III) above, the giant planets
open a gap in the gas disk. Consequently they become locked in the radial
evolution of the disk. As the disks material tends to be accreted by the star, the
giant planets have to migrate inwards. This migration, although slower than
that discussed above for the cores, is nevertheless quite fast. (This is the so-
called Type-II migration.) It is usually invoked to explain the existence of
hot Jupiters, massive extra-solar planets that orbit their star at distances
smaller than the orbital radius of Mercury. But in our solar system this kind of
migration evidently did not happen, or at least did not have a comparably large
radial extent. Again, two solutions to this problem have been advanced.
The rst possibility is that Jupiter and Saturn formed suciently late that
the disk was already in the dissipation phase. Thus, the disk disappeared
before it could signicantly move the planets. In addition, this solution has
the advantage of explaining why Jupiter and Saturn did not grow further in
mass and why their massive atmospheres are enriched in heavy elements
relative to the solar nebula composition (e.g., Guillot and Hueso, 2006).
A second possibility is that Jupiter and Saturn formed almost contem-
poraneously and on orbits that were close to each other. In this case, the gaps
opened by the two planets in the disk would have overlapped (Masset and
Snellgrove, 2001). This would have changed dramatically the migration
evolution of the planets pair, possibly stopping or even reversing it. Of course
the two solutions imply dierent formation timescales. If the planets stopped
because the gas disappeared, the fact that Jupiter and Saturn did not migrate
signicantly implies that their formation timescale is of order of the nebula
dissipation time (310 Myr; Hillenbrand 2006). In the opposite case, they
might have formed even faster.
Finally, to be complete, a model of giant planet formation alternative
to that of Pollack et al. has been proposed by Boss (see Boss, 2003 and ref-
erences herein). In this model, the proto-planetary gas disk was massive
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

enough to be unstable under its own gravity. In this situation, gaseous planets
can form very rapidly, by a process that reproduces in miniature the one that
led to the formation of the central star. There is no need of the presence of a
massive solid core to trigger the capture of the giant planets atmosphere. In
this sense, Bosss model may bring a solution to the timescale problem related
with the Pollack et al. model, discussed above. There is a still open debate in
the literature on whether the clumps of gas observed in numerical simulations
of gravitationally unstable disks are temporary features or would persist
until the disks dissipation. For the case of our solar system, the presence of
massive cores inside all giant planets, and the limited amount of gas in Uranus
and Neptune, make us think that the Pollack et al. model is more appropriate.
It is possible, however, that some or several of the extra-solar planets observed
so far formed through a gravitational instability mechanism.

3.3. The First 100 Million Years: the Telluric Era

3.3.1. FORMATION OF THE TERRESTRIAL PLANETS AND PRIMORDIAL


SCULPTING OF THE ASTEROID BELT

ALESSANDRO MORBIDELLI

After a few 105 years of runaway growth, the embryos in the terrestrial planet
region and in the asteroid belt region have Lunar to Martian masses. They
govern the local dynamics and start perturbing each other. The system be-
comes unstable, and the embryos orbits begin to intersect (Chambers and
Wetherill, 1998). Because of mutual close encounters, the embryos dynam-
ical excitation (increase of eccentricity and inclination) moderately increases,
and accretional collisions among embryos start to occur. The situation
drastically changes when Jupiter and Saturn acquire their current masses.
These two planets strongly perturb the dynamical evolution of the embryos
in the asteroid belt region between ~2 and 5 AU. The embryos acquire a
strong dynamical excitation, begin to cross each other, and cross rather
frequently the orbits of the embryos in the terrestrial planets region. The
collision rate increases. Despite the high relative velocity, these collisions lead
to accretion because of the large mass of the embryos.
The typical result of this highly chaotic phase simulated with several
numerical N-body integrations is the elimination of all the embryos orig-
inally situated in the asteroid belt and the formation a small number of
terrestrial planets on stable orbits in the 0.52 AU region in a timescale
~100 Myr (Figure 3.23).
This scenario has several strong points:
THIERRY MONTMERLE ET AL.

Figure 3.23. The growth of terrestrial planets from a disk of planetary embryos. Each panel
shows the semi-major axis and eccentricity of the bodies in the system, the size of each dot
being proportional to the mass. The color initially reects the starting position of each em-
bryo. When two (or more) embryos collide, the formed object assumes the color corre-
sponding to the embryo population that has mostly contributed to its total mass. A system of
four terrestrial planets, closely resembling our solar system, is formed in 200 Myr. From
Chambers (2001).

(i) Planets are formed on well separated and stable orbits only inside
~2 AU. Their number typically ranges from 2 to 4, depending on the
simulations, and their masses are in the range Mars mass Earth mass
(Chambers and Wetherill, 1998; Agnor et al., 1999).
(ii) Quasi-tangent collisions of Mars-mass embryos onto the proto-
planets are quite frequent (Agnor et al., 1999). These collisions are
expected to generate a disk of ejecta around the proto-planets (Canup
and Asphaug, 2001), from which a satellite is likely to accrete (Canup
and Esposito, 1996). This is the standard, generally accepted, scenario
for the formation of the Moon (see below, Section 3.2)
(iii) The accretion timescale of the terrestrial planets is ~100 Myr. This is
compatible with several constraints on the chronology of accretion
coming from geochemistry (Allegre et al. 1995). On the other hand,
HfW chronology seems to indicate that the formation of the Earths
core occurred within the rst 40 Myr (Yin et al., 2002; Kleine et al.,
2002). This might suggest that the Earth accretion was faster than it
appears in the simulations. However, if the cores of the embryos are
not mixed with the mantles during the collisions as indicated by
SPH simulations (Canup and Asphaug, 2001) this timescale would
measure the mean dierentiation age of the embryos that participated
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

to the formation of the Earth, and not the time required for our planet
to accrete most of its mass.
(iv) All the embryos located beyond 2 AU are eliminated in 2/3 of the
simulations (Chambers and Wetherill, 2001). They either are
dynamically ejected from the solar system, or collide with the Sun, or
are accreted by the forming terrestrial planets.
(v) In the same time, the small planetesimals are subject to the combined
perturbations of the giant planets and of the embryos (Petit et al.,
2001). The dynamical excitation increasing very rapidly (timescale 1
2 Myr), most of the small planetesimals are eliminated in a few million
years by either the ejection from the solar system, or the collision with
the Sun or with a growing planet. In the asteroid belt (24 AU range),
this leads to a remaining population of small bodies (the asteroids) on
stable orbits with quite large eccentricities and inclinations, which
contains only a very small fraction of the total mass initially in the
region. This scenario explains well the current mass decit of the
asteroid belt, the eccentricity and semi-major axis distribution of the
largest asteroids and other more subtle properties of the asteroid belt
population, such as the partial mixing of taxonomic types.
However, this scenario of terrestrial planet formation suers from some
weaknesses:
(i) The nal orbits of the planets formed in the simulations are typically
too eccentric and/or inclined with respect to the real ones. This could
be due to the fact that the current simulations neglect the so-called
phenomenon of dynamical friction, namely the effect of a large pop-
ulation of small bodies, carrying cumulatively a mass comparable to
that of the proto-planets. Dynamical friction should damp the
eccentricities and inclinations of the most massive bodies.
(ii) Obliquities of the terrestrial planets should have random values.
However in reality, only one planet has a retrograde spin (Venus).
Moreover all planetary obliquities are compatible with an initial 0-
degrees obliquity, modied by the subsequent evolution in the
framework of the current architecture of the planetary system (Laskar
and Robutel, 1993).
(iii) The planet formed in the simulations approximately at the location of
Mars is typically too massive.

3.3.2. THE FORMATION OF THE MOON

ALESSANDRO MORBIDELLI

As explained above, the currently accepted model for the Lunar formation is
that of a giant impact occurring during the formation of the Earth. The
THIERRY MONTMERLE ET AL.

current view of terrestrial planet formation implies several giant impacts, and
impacts with an angular momentum similar of that of the EarthMoon
system are not rare, particularly during the end of the accretion process
(Agnor et al., 1999). Simulations of a Moon-forming impact have been done
since 1986, using SPH simulations. The most advanced, recent high-resolu-
tion simulations have been done by Canup and Asphaug (2001) and Canup
(2004a). A very detailed review on the Moon formation can be found in
Canup (2004b).
In the SPH simulations, three critaria has been used to judge the degree of
success: the formation of a circumplanetary disk with about a Lunar mass
outside of the Roche radius of the Earth, a mass of iron in the protolunar
disk that is about 10% of the total mass (the fraction present inside the
Moon), and a total angular momentum of the Earth-disk system of order of
that of the current EarthMoon system. In essence, in case of an early
formation of the Moon, when the proto-Earth was only 60% of the current
Earth mass, the impactor needs to be of about 30% of the total mass (proto-
Earth + impactor). If the impact is late, the impactor can be of about 10%
of the total mass, namely of order of the mass of Mars. The authors privilege
the late impact scenario, because otherwise the Moon would have accreted
too many siderophile elements after its formation, assuming that it accreted
about 10% of the mass that is required to collide with the Earth to complete
the Earths formation and a chondritic compostion of such material.
The dierentiation of the Moon can be dated using the HfW chro-
nometer, and turns out to have occurred at about 40 Myr after CAI for-
mation (Yin et al., 2002; Kleine et al., 2002). Thus, if the Moon-forming
impact was the last one (or close to the last one), the Earth formed (i.e.,
received its last giant impact) in a similar timescale. If on the contrary, the
Moon-forming impact was an early one, the formation of the Earth might
have taken longer. The HfW chronometer seems to indicate an age of 40 My
also for the dierentiation of the Earth, but in this case the interpretation is
less straightforward because, in case there is little equilibration between the
core and the Mantle during the giant impacts, the overall mechanical
accretion process of the Earth could have taken signicantly longer (Sasaki
and Abe, 2004).
In the SPH simulation of the Moon-forming impact, about 80% of the
material that ends in the proto-lunar disk comes from the impactor, rather
than from the proto-Earth. The Moon and the Earth have a very similar
composition under many aspects (for instance the oxygen isotopic compo-
sition is identical). The logical interpretation is then that the proto-Lunar
impactor and the proto-Earth had very similar compositions (and identical
oxygen isotope composition). This however seems in conict with the results
of N-body simulations of the accretion of terrestrial planets. These simula-
tions show that a planet forms by accreting material in a stochastic way from
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

a wide variety of heliocentric distances. Given that the oxygen isotope


compostion of Mars and of all meteorite classes dier from each other, it
then seems unlikely that the proto-Earth and the lunar impactor could have
exactly the same resulting composition. A solution of this apparent paradox
has been proposed by Pahlevan and Stevenson (2005). During the formation
of the proto-lunar disk there might have been enough isotopic exchange with
the proto-Earth to equilibrate the two distributions. Another possibility is
that the fraction of the proto-Lunar disk of terrestrial origin was in reality
much larger than estimated in the SPH simulations. In fact these simulations
assume a non-spinning Earth, so that the totality of the angular momentum
of the Earth-Moon system is carried by the impactor. This has the conse-
quence of having the impactor on a quasi-tangent trajectory relaitve to the
proto-Earth, which maximizes the amount of impactor material that goes in
geocentric orbit. If the Earth was spinning fast at the time of the Moon
formation, then the impactor might have had a more head-on trajectory and
less of its material would have contaminated the proto-Lunar disk. SPH
simulations with a fast rotating Earth have never been done up to now.
The accretion of the Moon from a protolunar disk has been simulated by
Ida et al. (1997) and Kokubo et al. (2000). The Moon forms very quickly, in
about 1 year, at a distance from the Earth of about 1.2 to 1.3 Roche radii. In
many cases, the Moon is accompanied by many moonlets, left-over of the
accretion process in the disk. More rarely, two Moons of similar mass are
formed. The subsequent evolution of these systems, simulated in Canup
et al. (1999), typically leads to an end-state with a single Moon on a stable
orbit.

3.3.3. TOWARDS 1 GYR: THE EARLY EVOLUTION OF THE EARTH

BERNARD MARTY

3.3.3.1. Formation and closure of the Earths atmosphere


The formation of the terrestrial atmosphere is a suite of complex processes
which are not yet fully understood. Indeed its abundance and isotopic
compositions dier from those of the solar nebula and of other known
cosmochemical components like comets or primitive meteorites. The general
consensus is that our atmosphere was formed by contributions of several
volatile-rich components and that it was subject to several episodes of loss to
space throughout its early history. Some of these loss events were
non-fractionating, that is, elements and isotopes were lost in the same pro-
portion. This might have been the case of large impacts. Because the isotopic
composition of some of the atmospheric elements such as noble gases dier
from known end-member compositions, one has to consider also that
THIERRY MONTMERLE ET AL.

isotopic fractionation took place, and that there are several escape processes
for the atmosphere that have been demonstrated to isotopically fractionate
noble gases. Among these, there are the thermal loss, in which the velocity of
a given isotope is higher than that of another given isotope, allowing the
former to escape at a larger rate, and the pick-up ion loss in which atoms at
the top of the atmosphere are ionized during charge transfer from solar wind
ions, light isotopes being statistically more prone for such ionization. Hence
it is convenient to dene an epoch at which the atmosphere became closed to
further loss into space.
During at least the rst tens of Myr of its existence, while the Sun was still
on its way to the main sequence, the Earth was subject to large-scale igneous
(volcanic) events, during which the proto-mantle exchanged volatile elements
with surface reservoirs. It is likely that a primitive atmosphere existed at this
time, but any chemical record of it has been erased owing to the high thermal
state of the Earth evidenced by core formation and magma ocean episodes.
Records of atmospheric processes at this time cannot be found directly at the
Earths surface at Present. The record of extinct radioactivity systems in
which parents dier from daughters by their respective volatilities give strong
clues on the timing of terrestrial dierentiation and the early cycle of volatile
elements. Noble gases are chemically inert and their isotopic composition can
only be modied by kinetic fractionation, or mixing with nucleosynthetic
components, or through nuclear reactions including extinct radioactivity
decays.
Xenon, the heaviest stable noble gas, is of particular interest because some
of its isotopes are the radioactivity products of three dierent decay systems
covering contrasted time intervals. Iodine is a volatile element for which
one isotope, 129I decays with T1/2 = 15.7 Myr to 129Xe. During terrestrial
magma ocean episodes, xenon, which has presumably a higher volatility than
iodine, was degassed preferentially to it. The amount of 129Xe in the atmo-
sphere, in excess of the non-radiogenic xenon composition, corresponds to a
closure interval of about 100 Myr (Allegre et al., 1995). Put in other
words, only a tiny fraction of radiogenic 129Xe has been retained in the
atmosphere, showing evidence that the atmosphere was open to loss in space
for at least several tens of Myr. The terrestrial mantle has kept even less
129
Xe. Heavy xenon isotopes e.g., 136Xe, are produced by the spontaneous
ssion of 244Pu (T1/2 = 82 Myr) so that the combination of both chro-
nometers allows one to date the closure of the terrestrial mantle at
6070 Myr (Kunz et al., 1998). This age could be interpreted as averaging
the period during which the magma ocean episodes declined enough to
quantitatively retain is volatile elements. Notably, heavy Xe isotopes are also
produced by the spontaneous ssion of 238U which decays with a half-life of
4.45 Gyr, so that it is possible to compute closure ages based on
244
Pu238U136Xe. Results indicate a closure age of about 400600 Myr
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

(Yokochi and Marty, 2005), which could correspond to a signicant decline


in mantle geodynamics. Xe isotope variations are therefore consistent with
large-scale dierentiation within the rst 100 Myr and also with prolonged
mantle activity at rates much higher than at Present during the Hadean. A
comparable scenario was derived from other radioactive tracers, notably the
146
Sm142Nd and 147Sm143Nd systems, and this view is fully consistent with
models linking the thermal state of the Earth with the rate of mantle
convection.

3.3.3.2. Geological evidence: Core dierentiation, magnetic eld


The time it took to build the Earth and the other terrestrial planets has been
investigated from two dierent approaches, theoretical modeling on one
hand (see previous sections), and absolute chronology on the other hand. In
short, numerical simulations indicate that the growth of terrestrial planets
was a geologically fast process. In the turbulent nascent solar system, dust
accreted into small bodies, for which gravity became the main coalescent
agent, within 104105 years. Models predict that bodies with sizes up to that
of Mars accreted within a few Myr. Independent evidence for rapid growth
stems from coupled variations of 142Nd, 182W, 129Xe and 131136Xe produced
by extinct radioactivities of 144Sm (T1/2 = 106 Myr), 182Hf (9 Myr), 129I
(15.7 Myr) and 244Pu (82 Myr) observed in Martian meteorites (Halliday
et al., 2001; Marty and Marti, 2002). These geochemical tracers attest that
Martian dierentiation including core formation, crustal development and
mantle degassing took place within 15 Myr, which is remarkably short, and
contemporaneous to the dierentiation of parent bodies of meteorites.
Furthermore, they also indicate that these heterogeneities were not
re-homogenized later on as is the case of the Earth, and therefore attest for a
low mantle convection rate on Mars, if any.
As dicussed previously, it took about 100 Myr to complete the Earth,
considerably longer than the time of Mars-sized object formation (Wether-
hill, 1980). This longer time interval is the result of the decreasing probability
for collisions to occur in an increasingly depleted population of growing
bodies. This modeling approach does not allow one to get more detailed
chronological denition within this time frame, which needs to be con-
strained independently. The tungsten isotopic composition of the Earth
compared to that of the nascent solar system as recorded in primitive
meteorites supports large-scale dierentiation within a few tens of million
years. 182Hf decays into 182W with a half-life T1/2 = 9 Myr. During
metalsilicate dierentiation, tungsten is siderophile, that is, concentrates
preferentially in metal phase whereas hafnium prefers the silicates (litho-
phile). If, in a dierentiated body, this event took place when 182Hf was still
THIERRY MONTMERLE ET AL.

present, the isotopic composition of tungsten is dierent from the solar one
as recorded for example in undierentiated meteorites.
In pioneering attempts to use this chronometer, no dierence in the
tungsten isotopic composition was found between carbonaceous chondrites
(the most primitive meteorites found so far) and terrestrial silicates, leading
the authors to conclude that the last global HfW dierentiation of the Earth
happened after all 182W decayed, in practice 60 Myr after the start of solar
system condensation t0 (Lee and Halliday, 1995). More recently, dierences
between chondrites and the Earth have been found (Kleine et al., 2002; Yin
et al., 2002), implying a mean metalsilicate dierentiation of 30 Myr for the
Earth if terrestrial dierentiation was a single and global event, with a pos-
sible range of 1150 Myr for more realistic accreting conditions. A collision
between a Mars-sized object and a growing proto-Earth is consistent with the
unique EarthMoon angular momentum. Numerical models for such a giant
impact indicate that the proto-Earth was severely disrupted while material
from the proto-mantle and the impactor spiraled at high temperature and
formed the Moon (see above, Section 3.3.2). It is therefore logical, even if it
has not yet been demonstrated, to ascribe to this collision the last global
dierentiation between metal and silicate recorded in the HfW system.
Recent tungsten data for lunar basalts indicate a 182W anomaly for the lunar
mantle, interpreted as record for dierentiation at 45 4 Myr (Kleine et al.,
2005). These authors proposed that this age represented the end of magma
ocean episodes on the Moon.
The earliest record of a geomagnetic eld dates back to the Archean (Hale
and Dunlop, 1984; Yoshihara and Hamano, 2004), some 1 Gyr after the
formation period. The terrestrial magnetic eld has a major role in pre-
venting solar wind ions to interact extensively with the top of the atmosphere,
and create isotope fractionation of atmospheric elements by pick-up charge
exchange. It also tends to preserve the surface of the Earth from cosmic-ray
bombardment that are lethal for the development of organic chemistry.
However, there is no evidence for the existence of magnetic eld induced by
the geodynamo during the rst Gyr.

References

Abergel, A. et al.: 1996, Astron. Astrophys. 315, L329L332.


Adams, F. C. and Myers, P. C.: 2001, Astrophys. J. 553, 744753.
Adams, F. C., Hollenbach, D., Laughlin, G. and Gorti, U.: 2004, Astrophys. J. 611, 360379.
Adams, F. C., Proszkow, E. M., Fatuzzo, M., and Myers, P. C.: 2006, Astrophys. J. 641,
504525.
Agnor, C. B., Canup, R. M. and Levison, H. F.: 1999, Icarus 142, 219237.
Alibert, Y., Mordasini, C., Benz, W. and Winisdoerer, C.: 2005, Astorn. Astrophys. 434,
343353.
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Allegre, C. J., Manhes, G. and Gopel, C.: 1995, Geochim. Cosmochim. Acta. 59, 14451456.
Amelin, Y., Krot, A. N., Hutcheon, I. D. and Ulyanov, A. A.: 2002, Science 297, 16781683.
Anders, E. and Grevesse, N.: 1989, Geochim. Cosmochim. Acta 53, 197214.
Andre, P. and Montmerle, T.: 1994, Astrophys. J. 420, 837862.
Andre, P., Ward-Thompson, D. and Barsony, M.: 2000, in V. Mannings, A. P. Boss and S. S.
Russell (eds.), Protostars and Planets IV, University of Arizona Press, Tucson, pp. 5997.
Baker, J., Bizzarro, M., Wittig, N., Connelly, J. and Haack, H.: 2005, Nature 436, 11271131.
Balbus, S. A. and Hawley, J. F.: 1991, Astrophys. J. 376, 214233.
Bate, M. R. and Bonnell, I. A.: 2004, in H. J. G. L. M. Lamers, L. J. Smith, and A. Nota (eds.),
The Formation and Evolution of Massive Young Star Clusters. ASP Conf. Ser., 322, 289.
Beaulieu, J.-Ph. et al.: 2006, Nature 439, 437440.
Benz, W.: 2000, Sp. Sci. Rev. 92, 279294.
Bertout, C.: 1989, Ann. Rev. Astron. Astrophys. 27, 351395.
Bizzarro, M., Baker, J. A. and Haack, H.: 2004, Nature 431, 275278.
Boss, A. P.: 2003, Astrophys. J. 599, 577581.
Bouvier, A., Blichert-Toft, J., Vervoort, J. D., McClelland, W. and Albarede, F.: 2005,
Geochim. Cosmochim. Acta Suppl. 69(Suppl. 1), A384.
Bouvier, J. et al.: 1999, Astron. Astrophys. 349, 619635.
Bouvier, J. et al.: 2003, Astron. Astrophys. 409, 169192.
Brearley, A. and Jones, R. H.: 1998, in J. J. Papike (ed.), Planetary materials. Rev. Mineral.
36, 3/13/398.
Busso, M., Gallino, R. and Wasserburg, G. J.: 2003, Pub. Astr. Soc. Austr. 20, 356370.
Canup, R. M.: 2004a, Icarus 168, 433456.
Canup, R. M.: 2004b, Ann. Rev. Astron. Astrophys. 42, 441475.
Canup, R. M. and Asphaug, E.: 2001, Nature 412, 708712.
Canup, R. M., Levison, H. F. and Stewart, G. R.: 1999, Astr. J. 117, 603620.
Canup, R. M. and Esposito, L. W.: 1996, Icarus 119, 427446.
Carpenter, J. M., Wolf, S., Schreyer, K., Launhardt, R. and Henning, T.: 2005, Astr. J. 129,
10491062.
Chambers, J. E. and Wetherill, G. W.: 1998, Icarus 136, 304327.
Chambers, J. E. and Wetherill, G. W.: 2001, Met. Plan. Sci. 36, 381399.
Chambers, J. E.: 2001, Icarus 152, 205224.
Chaussidon, M. and Gounelle, M.: 2006, in D. Lauretta and L. Leshin (eds.), Meteorites &
Early Solar System II. Tucson, University of Arizona Press, pp. 323340.
Chaussidon, M., Libourel, G. and Krot, A. N.: 2006, Lunar Planet. Sci. 36, 1335.
Chaussidon, M., Robert, F. and McKeegan, K. D.: 2006, Geochim. Cosmochim. Acta 70, 224
245.
Chokshi, A., Tielens, A. G. G. M. and Hollenbach, D.: 1993, Astrophys. J. 407, 806819.
Clark, P. C., Bonnell, I. A., Zinnecker, H. and Bate, M. R.: 2005, Mon. Not. Roy. Astr. Soc.
359, 809818.
Crutcher R. M.: 2005, in Magnetic elds in the Universe: From Laboratory and Stars to
Primordial Structures, AIP Conf. Proc., vol. 784, pp. 129139.
Diehl, R., Halloin, H., Kretschmer, K., Lichti, G. G., Schonfelder, V., Strong, A. W., von
Kienlin, A., Wang, W., Jean, P., Knodlseder, J., et al.: 2006, Nature 439, 4547.
Diehl, R., Cervino, M., Hartmann, D. H. and Kretschmar, K.: 2004, New Astr. Rev. 48, 8186.
Donati, J.-F., Paletou, F., Bouvier, J. and Ferreira, J.: 2005, Nature 438, 466469.
Draine, B. T.: 2006, Astrophys. J. 636, 11141120.
Duchene, G., McCabe, C., Ghez, A. M. and Macintosh, B. A.: 2004, Astrophys. J. 606, 969
982.
Ehrenfreund, P. and Charnley, S. B.: 2000, Ann. Rev. Astr. Ap. 38, 427483.
THIERRY MONTMERLE ET AL.

Falgarone, E., Hily-Blant, P., Pety, J. and Pineau Des Forets, G.: 2005 in Magnetic elds in
the Universe: From Laboratory and Stars to Primordial Structures, AIP Conf. Proc., vol.
784, pp. 299307.
Favata, F., Micela, G., Silva, B., Sciortino, S. and Tsujimoto, M.: 2005, Astron. Astrophys.
433, 10471054.
Feigelson, E. D.: 2005, Met. Planet. Sci. 40(Suppl.), 5339.
Feigelson, E. D. and Montmerle, T.: 1999, Ann. Rev. Astron. Astrophys. 37, 363408.
Feigelson, E. D., Getman, K., Townsley, L., Garmire, G., Preibisch, T., Grosso, N., Mont-
merle, T., Muench, A. and McCaughrean, M.: 2005, Astrophys. J. Suppl. Ser. 160, 379
389.
Ferreira, J., Pelletier, G. and Appl, S.: 2000, Mon. Not. R. Astron. Soc. 312, 387397.
Ferreira, J. and Casse, F.: 2004, Astrophys. J. 601, L139L142.
Fleming, T. and Stone, J. M.: 2003, Astrophys. J. 585, 908920.
Fromang, S., Balbus, S. A., Terquem, C. and De Villiers, J.-P.: 2004, Astrophys. J. 616, 364
375.
Fromang, S., Terquem, C. and Balbus, S. A.: 2002, Mon. Not. R. Astron. Soc. 329, 1828.
Gallino, R., Busso, M., Wasserburg, G. J. and Straniero, O.: 2004, New Astr. Rev. 48, 133
138.
Galy, A., Young, E. D., Ash, R. D. and ONions, R. K.: 2002, Science 290, 17511753.
Getman, K. V., Feigelson, E. D., Grosso, N., McCaughrean, M. J., Micela, G., Broos, P.,
Garmire, G. and Townsley, L.: 2005, Astrophys. J. Suppl. Ser. 160, 353378.
Glassgold, A. E., Feigelson, E. D. and Montmerle, T.: 2000, in V. Mannings, A. P. Boss and
S. S. Russell (eds.), Protostars and Planets IV, University of Arizona Press, Tucson, pp.
429455.
Glassgold, A. E., Feigelson, E. D., Montmerle, T. and Wolk, S.: 2005, in A. N. Krot, E. R. D.
Scott and B. Repurth (eds.), Chondrites and the Protoplanetary Disk. ASP Conf. Ser., vol.
341, pp. 165182.
Goldreich, P. and Ward, W. R.: 1973, Astrophys. J. 183, 10511061.
Goodwin, S. P., Whitworth, A. P. and Ward-Thompson, D.: 2004, Astron. Astrophys. 423,
169182.
Gounelle, M. and Russell, S. S.: 2005, Geochim. Cosmochim. Acta 69, 31293144.
Gounelle, M., Shu, F. H., Shang, H., Glassgold, A. E., Rehm, K. E. and Lee, T.: 2004, Lunar
Planet. Sci. Conf. 35, 1829.
Gounelle, M., Shu, F. H., Shang, H., Glassgold, A. E., Rehm, K. E. and Lee, T.: 2001,
Astrophys. J. 548, 10511070.
Greaves, J. S.: 2005, Mon. Not. R. Astron. Soc. Letters 364, L47L50.
Greenberg, R., Wacker, J. F., Hartmann, W. K. and Chapman, C. R.: 1978, Icarus 35, 126.
Grosso, N., Alves, J., Wood, K., Neuhauser, R., Montmerle, T. and Bjorkman, J. E.: 2003,
Astrophys. J. 586, 296305.
Guillot, T.: 1999, Science 286, 7277.
Guillot, T.: 2005, Ann. Rev. Earth Plan. Sci. 33, 493530.
Guillot, T. and Hueso, R.: 2006, Mon. Not. R. Astron. Soc. 367, L47L51.
Guilloteau, S., Dutrey, A. and Simon, M.: 1999, Astr. Astrophys. 348, 570578.
Haisch, K. E., Lada, E. A. and Lada, C. J.: 2001, Astrophys. J. Lett. 553, L153L156.
Hale, C. J. and Dunlop, D.: 1984, Geophys. Res. Lett. 11, 97100.
Halliday, A. N., Wanke, H., Birck, J. L. and Clayton, R. N.: 2001, Space Sci. Rev. 96, 134.
Harker, D. E. and Desch, S. J.: 2002, Astrophys. J. Lett. 565, L109L112.
Harris, M. J., Knodlseder, J., Jean, P., Cisana, E., Diehl, R., Lichti, G. G., Roques, J.-P.,
Schanne, S. and Weidenspointner, G.: 2005, Astron. Astrophys. 433, L49L52.
Hayashi, C.: 1966, Ann. Rev. Astron. Astrophys. 4, 171192.
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Hayashi, C.: 1981, Prog. Theor. Phys. 70(Suppl.), 3553.


Hayashi, M., Shibata, K. and Matsumoto, R.: 2000, Adv. Sp. Res. 26, 567570.
Heiles, C., Haner, L. M., Reynolds, R. J. and Tufte, S. L.: 2000, Astrophys. J. 536, 335346.
Hester, J. J. and Desch, S.: 2005, in A. N. Krot and E. R. D. Scott and B. Reipurth (eds.),
Chondrites and the Protoplanetary Disk. ASP Conf series, vol. 341, pp. 107130.
Hewins, R. H.: 1997, Ann. Rev. Earth Planet. Sci. 25, 6183.
Hewins, R. H., Connolly, H. C. Jr., Lofgren, G. E. and Libourel, G.: 2005, in A. N. Krot, E.
R. D. Scot and B. Reipurth (eds.), Chondrites and the Protoplanetary Disk. ASP Conf.
Series, vol. 341, pp. 286316.
Hillenbrand, L. A.: 1997, Astron. J. 113, 17331768.
Hillenbrand, L. A.: 2006, in Livio M. (ed.), A Decade of Discovery: Planets Around Other
Stars. STScI Symposium Series 19, in press.
Huss, G. R. and Tachibana, S.: 2004, Lunar Planet. Sci. Conf. 35, 1811.
Ida, S., Canup, R. M. and Stewart, G. R.: 1997, Nature 389, 353357.
Ida, S. and Makino, J.: 1993, Icarus 106, 210227.
Ilgner, M. and Nelson, R. P.: 2006, Astron. Astrophys. 445, 205222.
Johns-Krull, C. M. and Gaord, A. D.: 2002, Astrophys. J. 573, 685698.
Kessler-Silacci, J. E., Augereau, J.-C. and Dullemond, C. P. et al.: 2006, Astrophys. J. 639,
275291.
Kleine, T., Metzger, K. and Palme, H.: 2005, Lunar Planet. Sci. 36, 1940, CD-ROM.
Kleine, T., Munker, C., Mezfer, K. and Palme, H.: 2002, Nature 418, 952955.
Kokubo, E. and Ida, S.: 1998, Icarus 131, 171178.
Kokubo, E., Ida, S. and Makino, J.: 2000, Icarus 148, 419436.
Kroupa, P.: 2002, Science 295, 8291.
Krumholz, M. R., McKee, C. F. and Klein, R. I.: 2005, Nature 438, 332334.
Kunz, J., Staudacher, T. and Allegre, C. J.: 1998, Science 280(5365), 877880.
Laskar, J. and Robutel, P.: 1993, Nature 361, 608612.
Lee, D. C. and Halliday, A.: 1995, Nature 378, 771774.
Lee, T., Shu, F. H., Shang, H., Glassgold, A. E. and Rhem, K. E.: 1998, Astrophys. J. 506,
898912.
Lissauer, J.: 1987, Icarus 69, 249265.
Lyo, A.-Ran, Lawson, W. A., Mamajek, E. E., Feigelson, E. D., Sung, Eon-Chang and
Crause, L. A.: 2003, Mon. Not. R. Astron. Soc. 338, 616622.
MacPherson, G. J., Davis, A. M. and Zinner, E. K.: 1995, Meteoritics 30, 365386.
MacPherson, G. J. and Huss, G. R.: 2003, Geochim. Cosmochim. Acta 67, 31653179.
Matsumura, S. and Pudritz, R. E.: 2006, Mon. Not. R. Astron. Soc. 365, 572584.
Matsumoto, R., Machida, M., Hayashi, M. and Shibata, K.: 2000, Progr. Th. Phys.
138(Suppl.), 632637.
Matzner, C. D. and McKee, C. F.: 2000, Astrophys. J. 545, 364378.
Megeath, S. T., Hartmann, L., Luhman, K. L. and Fazio, G. G.: 2005, Astrophys. J. Lett. 634,
L113L116.
Manhes G., Gopel C. and Allegre C. J.: 1988, Comptes Rendus de lATP Planetologie, 323
327.
Marty, B. and Marti, K.: 2002, Earth Plan. Sci. Letts. 196(34), 251263.
Masset, F., Morbidelli, A., Crida, A. and Ferreira, J.: 2006, Astrophys. J. 642, 478487.
Masset, F. and Snellgrove, M.: 2001, Mon. Not. R. Astron. Soc. 320, L55L59.
McKeegan, K. D. and Davis, A. M.: 2003, in H. Holl and K. Turekian (eds.), Treatise on
Geochemistry, Elsevier-Pergamon, Oxford, pp. 431460.
McKeegan, K. D., Chaussidon, M. and Robert, F.: 2000, Science 289, 13341337.
THIERRY MONTMERLE ET AL.

Meyer, B. S.: 2005, in A. N. Krot, E. R. D. Scott and B. Reipurth (eds.), Chondrites and the
Protoplanetary Disk, ASP Conf series, vol. 341, pp. 515526.
Micela, G. and Favata, F.: 2005, Sp. Sci. Rev. 108, 577708.
Montmerle, T.: 1979, Astrophys. J. 231, 95110.
Montmerle, T.: 2002, New Astr. Rev. 46, 573583.
Montmerle, T.: 2005, in M. Gargaud, B. Barbier, H. Martin and J. Reisse (eds.), Lectures in
Astrobiology, Springer, Heidelberg, pp. 2759.
Morbidelli, A. and Levison, H. F.: 2004, Astron. J. 128, 25642576.
Mostefaoui, S., Kita, N. T., Togashi, S., Tachibana, S., Nagahara, H. and Morishita, Y.:
2002, Meteoritics Planet. Sci. 37, 421438.
Mostefaoui, S., Lugmair, G. W., Hoppe, P. and El Goresy, A.: 2004, New Astron. Rev. 48,
155159.
Motte, F., Andre, P. and Neri, R.: 1998, Astron. Astrophys. 336, 150172.
Muzerolle, J., Calvet, N. and Hartmann, L.: 2001, Astrophys. J. 550, 944961.
Natta, A., Testi, L., Calvet, N., Henning, Th., Waters, R. and Wilner, D.: 2006, in: B.
Reipurth (ed.), Protostars and Planets V. Tucson, University of Arizona Press, in press
(eprint arXiv:astro-ph/0602041).
Nelson, R. P.: 2005, Astron. Astrophys. 443, 10671085.
Ozawa, H., Grosso, N. and Montmerle, T.: 2005, Astron. Astrophys. 429, 963975.
Padoan, P., Jimenez, R., Juvela, M. and Nordlund, A.: 2004, Astrophys. J. 604, L49L52.
Pahlevan, K. and Stevenson, D. J.: 2005, Lunar Planet. Sci. Conf. 36, 1505.
Palla, F. and Stahler, S. W.: 1999, Astrophys. J. 525, 772783.
Petit, J.-M., Morbidelli, A. and Chambers, J.: 2001, Icarus 153, 338347.
Pety, J. and Falgarone, E.: 2003, Astron. Astrophys. 412, 417430.
Podolak, M.: 2003, Icarus 165, 428437.
Podosek, F. A., Zinner, E. K., MacPherson, G. J., Lundberg, L. L., Brannon, J. C. and Fahey,
A. J.: 1991, Geochim. Cosmochim. Acta 55, 10831110.
Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J, Podolak, M. and Greenzweig, Y.:
1996, Icarus 124, 6285.
Poppe, T. and Blum, J.: 1997, Adv. Space. Res. 20, 15951604.
Poppe, T., Blum, J. and Henning, T.: 2000, Astophys. J. 533, 454471.
Reipurth, B. and Bally, J.: 2001, Ann. Rev. Astron. Astrophys. 39, 403455.
Ryter, C.: 1996, Astr. Sp. Sci. 236, 285291.
Sasaki, T. and Abe, Y.: 2004, Lunar Planet. Sci. Conf. 35, 1505.
Shang, H., Glassgold, A. E., Shu, F. H. and Lizano, S.: 2002, Astrophys. J. 564, 853876.
Shu, F. H., Adams, F. C. and Lizano, S.: 1987, Ann. Rev. Astron. Astrophys. 25, 2381.
Shu, F. H., Shang, H., Glassgold, A. E. and Lee, T.: 1997, Science 277, 14751479.
Stelzer, B. and Neuhauser, R.: 2001, Astron. Astrophys. 377, 538556.
Stolper, E. and Paque, J. M.: 1986, Geochim. Cosmochim. Acta 50, 17851806.
Tanga, P., Babiano, A., Dubrulle, B. and Provenzale, A.: 1996, Icarus 121, 158177.
Thommes, E. W., Duncan, M. J. and Levison, H. F.: 2003, Icarus 161, 431455.
Throop, H. B. and Bally, J.: 2005, Astrophys. J. 623, L149L152.
Tsujimoto, M., Feigelson, E. D., Grosso, N., Micela, G., Tsuboi, Y., Favata, F., Shang, H.
and Kastner, J. H.: 2005, Astrophys. J. Suppl. Ser. 160, 503510.
Wasserburg, G. J.: 1987, Earth Plan. Sci. Lett. 86, 129173.
Weidenschilling, S. J.: 1977, Mon. Not. R. Astron. Soc. 180, 5770.
Weidenschilling, S. J.: 1980, Icarus 44, 172189.
Wetherhill, G. W.: 1980, Ann. Rev. Astron. Astrophys. 18, 77.
Wetherhill, G. W.: 1990, Ann. Rev. of Earth Plan. Sciences 18, 205256.
SOLAR SYSTEM FORMATION AND EARLY EVOLUTION

Yin, Q., Jacobsen, S. B., Yamashita, K., Telouk, P., Blichert-Toft, J. and Albarede, F.: 2002,
Nature 418, 949952.
Yokochi, R. and Marty, B.: 2005, Earth Plan. Sci. Letts 238, 1730.
Young, E. D., Simon, J. I., Galy, A., Russell, S. S., Tonui, E. and Lovera, O.: 2005, Science
308, 223227.
Yoshihara, A. and Hamano, Y.: 2004, Precambr. Res. 131, 111142.
Zanda, B.: 2004, Earth Planet. Sci. Lett. 224, 117.

You might also like