• Abstract and Figures
• Public Full-text

Icarus, accepted 5/20/2009 A Preprint typeset using LTEX style emulateapj v. 12/14/05 BUILDING THE TERRESTRIAL PLANETS: CONSTRAINED ACCRETION IN THE INNER SOLAR SYSTEM Sean N. Raymond1, David P. O’Brien2, Alessandro Morbidelli3, & Nathan A. Kaib4 (Received; Accepted) Icarus, accepted 5/20/2009 ABSTRACT To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include 1) the orbits, in particular the small eccentricities, and 2) the masses of the terrestrial planets – Mars’ relatively small mass in particular has not been adequately reproduced in previous simulations; 3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; 4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and 5) Earth’s relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N=1000-2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the “Nice model.” We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars’ mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e =0.07 0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely− allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. Our simulations leave us with a problem: we can reasonably satisfy the observed constraints (except for Earth’s water) with a configuration of Jupiter and Saturn that is at best marginally consistent with models of the outer Solar System, as it does not allow for any outer planet migration after a few Myr. Alternately, giant planet configurations which are consistent with the Nice model fail to reproduce Mars’ small size. Subject headings: Terrestrial Planets — Planetary Formation — Accretion — Origin, Solar System 1. INTRODUCTION small because of strong gravitational focusing such that dM/dt M 4/3 (Safronov 1969; Greenberg et al.1978). It is commonly accepted that rocky planets form by ∼ the process of collisional agglomeration of smaller bod- However, viscous stirring by the large bodies increases ies (for recent reviews, see Chambers 2004, Nagasawa the velocity dispersion of planetesimals, thereby reducing et al.2007 or Raymond 2008). This process starts from the growth rate to a roughly geometrical regime, where micron-sized dust grains in young circumstellar disks, dM/dt M 2/3 (Ida & Makino 1993). Dynamical friction ∼ and the current paradigm proceeds as follows. Grains acts on the oligarchs, maintaining small eccentricities settle to a thin disk midplane on a 104 year timescale (Ida & Makino 1992; Kokubo & Ida 1998). The building (Weidenschilling 1980), and grow quickly∼ via sticky col- blocks of the terrestrial planets, Moon-sized planetary 5 6 ∼ lisions until they reach cm- or m- sizes (Dullemond & embryos, form in 10 10 years with a characteristic − Dominik 2004). The time for m-sized bodies to spiral spacing of 5-10 mutual Hill radii (Wetherill & Stewart in to the star is very short ( 100 years) such that 1993; Weidenschilling et al. 1997; Kokubo & Ida 2000, arXiv:0905.3750v1 [astro-ph.EP] 22 May 2009 this size range constitutes a barrier∼ to further growth 2002). Giant collisions between planetary embryos begin (Weidenschilling 1977a). This barrier may be crossed to occur when the local density of planetesimals and em- by rapid accretion (Weidenschilling & Cuzzi 1993; Benz bryos is comparable (Wetherill 1985; Kenyon & Bromley 2000) or by local gravitational instability (Goldreich & 2006). During late-stage accretion, embryo-planetesimal Ward 1973; Youdin & Shu 2002), which can be trig- and embryo-embryo impacts are common and the feeding gered by turbulent concentration (Johansen et al.2007; zones of terrestrial planets can span several AU in width Cuzzi et al.2008). Larger bodies (100 m to 100 km in (Raymond et al.2006a). Late-stage accretion lasts for size), which are more weakly coupled to the gaseous disk, 108 years and sets the final bulk architecture of the ∼ are called planetesimals. Runaway growth of the largest system as well as the composition of the terrestrial plan- planetesimals may occur while the velocity dispersion is ets (e.g., Wetherill 1996). Past simulations of late-stage accretion have succeeded 1 Center for Astrophysics and Space Astronomy, Uni- in reproducing several aspects of the Solar System’s ter- versity of Colorado, UCB 389, Boulder CO 80309-0389; restrial planets. Using only 20-165 particles, Agnor et [email protected] 2 Planetary Science Institute, Tucson, AZ al.(1999) and Chambers (2001) roughly reproduced the 3 Observatoire de la Cˆote d’Azur, Boulevard de l’Observatoire, approximate masses and semimajor axes of Mercury, BP 4229, 06304 Nice Cedex 4, France. Venus, Earth and Mars. Thommes et al.(2008) also re- 4 Department of Astronomy, University of Washington, Seattle, produced the rough mass distribution of the inner Solar WA 98195 System by invoking sweeping secular resonances during 2 Raymond, O’Brien, Morbidelli, & Kaib the depletion of the Solar Nebula. By taking dynam- certain constraints. We present results and analysis of ical friction from remnant planetesimals into account, all simulations in section 5. We discuss these results and O’Brien et al. (2006) and Morishima et al. (2008) repro- present our conclusions in section 6. duced the very low eccentricities of the terrestrial plan- ets. Several groups have succeeded in delivering water to 2. INNER SOLAR SYSTEM CONSTRAINTS Earth from hydrated asteroidal material, following the We consider five broad attributes which we attempt to model of Morbidelli et al. (2000; see also Raymond et reproduce statistically with accretion simulations. Other al. 2004, 2006a, 2007; O’Brien et al. 2006). observations and measurements exist for inner Solar Sys- Despite these achievements, no previous study has ade- tem bodies, but we are limiting ourselves to relatively quately reproduced all aspects of the inner Solar System. broad and well-understood characteristics. These con- Indeed, as pointed out by Wetherill (1991), Mars’ small straints are described below in order from strongest to size remains the most difficult constraint to reproduce weakest. Weaker constraints rely on models or data that (also discussed in Chambers 2001). Agnor et al. (1999), are subject to interpretation, while strong constraints are Chambers (2001) and Morishima et al.(2008) succeeded directly observed. We use several quantities to compare in reproducing Mars’ small size only because their simu- our simulations with the Solar System’s terrestrial plan- lations started from an annulus of material with a fixed ets. These include statistical measures that were intro- width (see also Kominami & Ida 2002). In most cases duced by Chambers (2001). this annulus extended from 0.5-1.5 AU, such that a small planet could form at the outer edge of the initial disk be- 1. The masses and the mass distribution of the terres- cause of spreading. However, no such edge is thought to trial planets. As mentioned above, the mass distri- have existed in the Solar Nebula, so that the assumption bution of the terrestrial planets, and in particular that the planetesimal and embryo population extended the small masses of Mercury and Mars, have not only to 1.5 AU is not justified. Chambers (2001) man- been adequately reproduced in the context of the aged to place Mercury within a planetary mass distri- entire Solar System and its history. In this paper bution but only by adopting an ad-hoc inner disk pro- we do not attempt to reproduce Mercury because file. Thommes et al.(2008) formed a small Mars but its small size and large iron content may be the re- the orbits they assumed for Jupiter and Saturn are in- sult of a mantle-stripping impact (Benz et al. 1988) consistent with any significant late, planetesimal-driven or interesting composition-sorting gaseous effects migration of the giant planets (discussed at length in 6.2 § (Weidenschilling 1978). However, for the case of below). In fact, the scenario of Thommes et al.(2008) is Mars, with its more distant orbit, these effects are incompatible with the two currently viable theories for less likely to be a factor, and it should be repro- the late heavy bombardments because these require ei- ducible in the context of our simulations. In ad- ther a more compact configuration of Jupiter and Saturn dition, the distribution of mass in the inner Solar (Gomes et al.2005) or the formation of a small, sub- System is interesting because the majority is con- Mars-sized planet at 2 AU (Chambers 2007). centrated between the orbits of Venus and Earth. Terrestrial accretion∼ lasts for 108 years, far longer ∼ We therefore use two statistical measures for this than the few Myr lifetimes of the gaseous component constraint: of protoplanetary disks (Haisch et al. 2001; Brice˜no et al.2001; Pascucci et al.2006). Thus, gas giant planets The number of planets formed N . We take must be fully-formed during late-stage accretion and can • p therefore strongly affect terrestrial bodies, especially if Np to represent objects that contain at least the giant planets’ orbits are eccentric (Wetherill 1996; one planetary embryo, that have semimajor Chambers & Cassen 2002; Levison & Agnor 2003; Ray- axes a < 2 AU, and that are on long-term mond et al.

Load more

  21

Copy Link