Quick Search:

uzh logo
Browse by:

Zurich Open Repository and Archive

Maintenance: Tuesday, 5.7.2016, 07:00-08:00

Maintenance work on ZORA and JDB on Tuesday, 5th July, 07h00-08h00. During this time there will be a brief unavailability for about 1 hour. Please be patient.

Permanent URL to this publication: http://dx.doi.org/10.5167/uzh-34360

Morishima, R; Stadel, J; Moore, B (2010). From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets. Icarus, 207(2):517-535.

Accepted Version
View at publisher


We present results from a suite of N-body simulations that follow the formation and accretion history of the terrestrial planets using a new parallel treecode that we have developed. We initially place 2000 equal size planetesimals between 0.5 and 4.0 AU and the collisional growth is followed until the completion of planetary accretion (>100 Myr). A total of 64 simulations were carried out to explore sensitivity to the key parameters and initial conditions. All the important effect of gas in laminar disks are taken into account: the aerodynamic gas drag, the disk–planet interaction including Type I migration, and the global disk potential which causes inward migration of secular resonances as the gas dissipates. We vary the initial total mass and spatial distribution of the planetesimals, the time scale of dissipation of nebular gas (which dissipates uniformly in space and exponentially in time), and orbits of Jupiter and Saturn. We end up with 1–5 planets in the terrestrial region. In order to maintain sufficient mass in this region in the presence of Type I migration, the time scale of gas dissipation needs to be 1–2 Myr. The final configurations and collisional histories strongly depend on the orbital eccentricity of Jupiter. If today’s eccentricity of Jupiter is used, then most of bodies in the asteroidal region are swept up within the terrestrial region owing to the inward migration of the secular resonance, and giant impacts between protoplanets occur most commonly around 10 Myr. If the orbital eccentricity of Jupiter is close to zero, as suggested in the Nice model, the effect of the secular resonance is negligible and a large amount of mass stays for a long period of time in the asteroidal region. With a circular orbit for Jupiter, giant impacts usually occur around 100 Myr, consistent with the accretion time scale indicated from isotope records. However, we inevitably have an Earth size planet at around 2 AU in this case. It is very difficult to obtain spatially concentrated terrestrial planets together with very late giant impacts, as long as we include all the above effects of gas and assume initial disks similar to the minimum mass solar nebular.
From planetesimals to terrestrial planets: N-body simulations including the effects of nebular gas and giant planets
Morishima, Ryuji; Stadel, Joachim; Moore, Ben
AA(Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA; Institute of Geophysics and Planetary Physics, UCLA, Los Angels, CA 90065, USA), AB(Institute for Theoretical Physics, University of Zurich, Switzerland), AC(Institute for Theoretical Physics, University of Zurich, Switzerland)
Icarus, Volume 207, Issue 2, p. 517-535. (Icarus Homepage)
Publication Date:
Abstract Copyright:
Elsevier Inc.
Bibliographic Code:


37 citations in Web of Science®
39 citations in Scopus®
Google Scholar™



48 downloads since deposited on 03 Mar 2011
14 downloads since 12 months

Detailed statistics

Additional indexing

Item Type:Journal Article, refereed, original work
Communities & Collections:07 Faculty of Science > Institute for Computational Science
Dewey Decimal Classification:530 Physics
Date:June 2010
Deposited On:03 Mar 2011 16:15
Last Modified:05 Apr 2016 14:09
Publisher DOI:10.1016/j.icarus.2009.11.038
Related URLs:http://arxiv.org/abs/1007.0579

Users (please log in): suggest update or correction for this item

Repository Staff Only: item control page