Source author record

Alessandro Morbidelli

Alessandro Morbidelli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

84works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

84 published item(s)

preprint2026arXiv

Two-source terrestrial planet formation with a sweeping secular resonance

The models that most successfully reproduce the orbital architecture of the Solar System terrestrial planets start from a narrow annulus of material that grows into embryos and then planets. However, it is not clear how this ring model can be made consistent with the chemical structure of the inner Solar System, which shows a reduced-to-oxidized gradient from Mercury to Mars and a parallel gradient in the asteroid belt. We propose that there were two primary reservoirs in the early inner Solar System: a narrow, refractory enriched ring inside of 1 au, and a less massive, extended planetesimal disk outside of 1 au with oxidation states ranging from enstatite chondrites to ordinary chondrites. We show through a suite of N-body simulations that an inwardly sweeping secular resonance, caused by aerodynamic drag and perturbations from a mean-motion resonant Jupiter and Saturn, gathers the outer planetesimal disk into a narrow ring that migrates radially, forms Mars, and contributes oxidized material to proto-Earth. Remaining unaccreted planetesimals can be implanted into the asteroid belt as the parent bodies of aubrites and non-carbonaceous iron meteorites, while the most reduced material is not implanted and thus unsampled in the meteorite collection. This model explains the oxidation and isotopic gradients within the inner Solar System within the context of a low-viscosity, magnetic wind-driven disk.

preprint2023arXiv

Formation of Rocky Super-Earths From A Narrow Ring of Planetesimals

The formation of super-Earths, the most abundant planets in the Galaxy, remains elusive. These planets have masses that typically exceed that of the Earth by a factor of a few; appear to be predominantly rocky, although often surrounded by H/He atmospheres; and frequently occur in multiples. Moreover, planets that encircle the same star tend to have similar masses and radii, whereas those belonging to different systems exhibit remarkable overall diversity. Here, we advance a theoretical picture for rocky planet formation that satisfies the aforementioned constraints: building upon recent work - which demonstrates that planetesimals can form rapidly at discrete locations in the disk - we propose that super-Earths originate inside rings of silicate-rich planetesimals at approximately ~1 AU. Within the context of this picture, we show that planets grow primarily through pairwise collisions among rocky planetesimals, until they achieve terminal masses that are regulated by isolation and orbital migration. We quantify our model with numerical simulations and demonstrate that our synthetic planetary systems bear a close resemblance to compact, multi-resonant progenitors of the observed population of short-period extrasolar planets. Our results thus indicate that the absence of short-period super-Earths within the solar system can simply be attributed to the comparatively low mass of the primordial planetesimal ring within the protosolar nebula.

preprint2022arXiv

A Criterion for the Stability of Planets in Chains of Resonances

Uncovering the formation process that reproduces the distinct properties of compact super-Earth exoplanet systems is a major goal of planet formation theory. The most successful model argues that non-resonant systems begin as resonant chains of planets that later experience a dynamical instability. However, both the boundary of stability in resonant chains and the mechanism of the instability itself are poorly understood. Previous work postulated that a secondary resonance between the fastest libration frequency and a difference in synodic frequencies destabilizes the system. Here, we use that hypothesis to produce a simple and general criterion for resonant chain stability that depends only on planet orbital periods and masses. We show that the criterion accurately predicts the maximum mass of planets in synthetic resonant chains up to six planets. More complicated resonant chains produced in population synthesis simulations are found to be less stable than expected, although our criterion remains useful and superior to machine learning models.

preprint2022arXiv

Dynamical origin of the Dwarf Planet Ceres

The Dwarf Planet Ceres revealed the presence of ammonia and other unique properties compared to other asteroids in the main belt which suggests that it was not formed in situ. We model the early dynamical evolution of the outer Solar System to study possible dynamical mechanisms to implant a Ceres-sized planetesimal in the asteroid belt from the trans-Saturnian region. We calculate that the fraction of the population of Ceres-sized planetesimals that are captured in the asteroid belt is in the range of 2.8e-5 to 1.2e-3 depending on the initial location in the outer planetesimal disk. The captured bodies have a 70% probability to have a semimajor axis between 2.5 and 3 au, a 33% probability to have an eccentricity smaller than 0.2 and a 45% probability to have an orbital inclination smaller than 10 degrees. Assuming the existence of 3,600 Ceres-size planetesimals in the inner part of the trans-Saturnian disk, consistent with the estimate of Nesvorny & Vokrouhlicky (2016) for the trans-Neptunian disk, our estimated capture probability and a final 80% depletion of the asteroid belt during the subsequent giant planet instability, lead to capture 1 Ceres in the asteroid belt, with a probability of 15%, 34%, and 51% to be located in the inner, middle and outer belt respectively.

preprint2022arXiv

Stochastic accretion of the Earth

Earth is depleted in volatile elements relative to chondritic meteorites, its possible building blocks. The extent of this depletion increases with decreasing condensation temperature, and is approximated by a cumulative normal distribution, unlike that in any chondrite. However, moderately volatile elements, occupying the mid-range of the distribution, have chondritic isotope ratios, contrary to that expected from loss by partial vaporisation/condensation. Here we reconcile these observations by showing, using N-body simulations, that Earth accreted stochastically from many precursor bodies whose variable compositions reflect the temperatures at which they formed. Impact-induced atmospheric loss was efficient only when the proto-Earth was small, and elements that accreted thereafter retain near-chondritic isotope ratios. Earth's composition is reproduced when initial temperatures of planetesimal- to embryo-sized bodies are set by disk accretion rates of (1.08 $\pm$ 0.17) $\times$ 10$^{-7}$ solar masses/yr, although they may be perturbed by $^{26}$Al heating on bodies formed at different times. The model implies a heliocentric gradient in composition and rapid planetesimal formation within $\sim$ 1 Myr, in accord with radiometric volatile depletion ages of Earth.

preprint2022arXiv

Terrestrial planet formation from lost inner solar system material

Two fundamentally different processes of rocky planet formation exist, but it is unclear which one built the terrestrial planets of the solar system. Either they formed by collisions among planetary embryos from the inner solar system, or by accreting sunward-drifting millimeter-sized 'pebbles' from the outer solar system. We show that the isotopic compositions of Earth and Mars are governed by two-component mixing among inner solar system materials, including material from the innermost disk unsampled by meteorites, whereas the contribution of outer solar system material is limited to a few percent by mass. This refutes a pebble accretion origin of the terrestrial planets, but is consistent with collisional growth from inner solar system embryos. The low fraction of outer solar system material in Earth and Mars indicates the presence of a persistent dust-drift barrier in the disk, highlighting the specific pathway of rocky planet formation in the solar system.

preprint2021arXiv

Contemporary formation of early solar system planetesimals at two distinct radial locations

The formation of planetesimals is expected to occur via particle-gas instabilities that concentrate dust into self-gravitating clumps. Triggering these instabilities requires the prior pileup of dust in the protoplanetary disk. Until now, this has been successfully modeled exclusively at the disk's snowline, whereas rocky planetesimals in the inner disk were obtained only by assuming either unrealistically large particle sizes or an enhanced global disk metallicity. However, planetesimal formation solely at the snowline is difficult to reconcile with the early and contemporaneous formation of iron meteorite parent bodies with distinct oxidation states and isotopic compositions, indicating formation at different radial locations in the disk. Here, by modeling the evolution of a disk with ongoing accretion of material from the collapsing molecular cloud, we show that planetesimal formation may have been triggered within the first 0.5 million years by dust pileup at both the snowline (at approximately 5 au) and the silicate sublimation line (at approximately 1 au), provided turbulent diffusion was low. Particle concentration at approximately 1 au is due to the early outward radial motion of gas and is assisted by the sublimation and recondensation of silicates. Our results indicate that, although the planetesimals at the two locations formed about contemporaneously, those at the snowline accreted a large fraction of their mass (approximately 60 percent) from materials delivered to the disk in the first few 10,000 yr, whereas this fraction is only 30 percent for the planetesimals formed at the silicate line. Thus, provided that the isotopic composition of the delivered material changed with time, these two planetesimal populations should have distinct isotopic compositions, consistent with observations.

preprint2020arXiv

Characterization of Temporarily-Captured Minimoon 2020 CD$_3$ by Keck Time-resolved Spectrophotometry

We present time-resolved visible spectrophotometry of minimoon 2020 CD$_3$, the second asteroid known to become temporarily captured by the Earth-Moon system's gravitational field. The spectrophotometry was taken with Keck I/LRIS between wavelengths 434 nm and 912 nm in $B$, $g$, $V$, $R$, $I$ and RG850 filters as it was leaving the Earth-Moon system on 2020 March 23 UTC. The spectrophotometry of 2020 CD$_3$ most closely resembles the spectra of V-type asteroids and some Lunar rock samples with a reddish slope of ~18$\%$/100 nm between 434 nm and 761 nm corresponding to colors of $g$-$r$ = 0.62$\pm$0.08, $r$-$i$ = 0.21 $\pm$ 0.06 and an absorption band at ~900 nm corresponding to $i$-$z$ = -0.54$\pm$0.10. Combining our measured 31.9$\pm$0.1 absolute magnitude with an albedo of 0.35 typical for V-type asteroids, we determine 2020 CD$_3$'s diameter to be ~0.9$\pm$0.1 m making it the first minimoon and one of the smallest asteroids to be spectrally studied. We use our time-series photometry to detect periodic lightcurve variations with a $<$10$^{-4}$ false alarm probability corresponding to a lightcurve period of ~573 s and a lightcurve amplitude of ~1 mag implying 2020 CD$_3$ possesses a $b/a$ axial ratio of ~2.5. In addition, we extend the observational arc of 2020 CD$_3$ to 37 days between 2020 February 15 UTC and 2020 March 23 UTC. From the improved orbital solution for 2020 CD$_3$, we estimate its likely duration of its capture to be ~2 y, and we measure the non-gravitation perturbation on its orbit due to radiation pressure with an area-to-mass ratio of 6.9$\pm$2.4$\times$10$^{-4}$ m$^2$/kg implying a density of 2.3$\pm$0.8 g/cm$^3$, broadly compatible with the densities of other meter-scale asteroids and Lunar rock. We searched for pre-discovery detections of 2020 CD$_3$ in the ZTF archive as far back as 2018 October, but were unable to locate any positive detections.

preprint2020arXiv

Characterization of the Nucleus, Morphology and Activity of Interstellar Comet 2I/Borisov by Optical and Near-Infrared GROWTH, Apache Point, IRTF, ZTF and Keck Observations

We present visible and near-infrared photometric and spectroscopic observations of interstellar object 2I/Borisov taken from 2019 September 10 to 2019 November 29 using the GROWTH, the APO ARC 3.5 m and the NASA/IRTF 3.0 m combined with post and pre-discovery observations of 2I obtained by ZTF from 2019 March 17 to 2019 May 5. Comparison with imaging of distant Solar System comets shows an object very similar to mildly active Solar System comets with an out-gassing rate of $\sim$10$^{27}$ mol/sec. The photometry, taken in filters spanning the visible and NIR range shows a gradual brightening trend of $\sim0.03$ mags/day since 2019 September 10 UTC for a reddish object becoming neutral in the NIR. The lightcurve from recent and pre-discovery data reveals a brightness trend suggesting the recent onset of significant H$_2$O sublimation with the comet being active with super volatiles such as CO at heliocentric distances $>$6 au consistent with its extended morphology. Using the advanced capability to significantly reduce the scattered light from the coma enabled by high-resolution NIR images from Keck adaptive optics taken on 2019 October 04, we estimate a diameter of 2I's nucleus of $\lesssim$1.4 km. We use the size estimates of 1I/'Oumuamua and 2I/Borisov to roughly estimate the slope of the ISO size-distribution resulting in a slope of $\sim$3.4$\pm$1.2, similar to Solar System comets and bodies produced from collisional equilibrium.

preprint2020arXiv

Chemical Diversity of Super-Earths As a Consequence of Formation

Recent observations of rocky super-Earths have revealed an apparent wider distribution of Fe/Mg ratios, or core to mantle ratios, than the planets in our Solar System. This study aims to understand how much of the chemical diversity in the super-Earth population can arise from giant impacts during planetary formation. Planet formation simulations have only recently begun to treat collisions more realistically in an attempt to replicate the planets in our Solar System. We investigate planet formation more generally by simulating the formation of rocky super-Earths with varying initial conditions using a version of SyMBA, a gravitational N-body code, that incorporates realistic collisions. We track the maximum plausible change in composition after each impact. The final planets span a range of Fe/Mg ratios similar to the Solar System planets, but do not completely match the distribution in super-Earth data. We only form a few planets with minor iron-depletion, suggesting other mechanisms are at work. The most iron-rich planets have a lower Fe/Mg ratio than Mercury, and are less enriched than planets such as Kepler-100b. This indicates that further work on our understanding of planet formation and further improvement of precision of mass and radius measurements are required to explain planets at the extremes of this Fe/Mg distribution.

preprint2020arXiv

Debiased albedo distribution for Near Earth Objects

We extend the most recent orbital and absolute magnitude Near Earth Object (NEO) model (Granvik et al., 2018) to provide a statistical description of NEO geometric albedos. Our model is calibrated on NEOWISE albedo data for the NEO population and reproduces these data very well once a simple model for the NEOWISE observational biases is applied. The results are consistent with previous estimates. There are about 1,000 NEOs with diameter D>1km and the mean albedo to convert absolute magnitude into diameter is 0.147. We do not find any statistically significant evidence that the albedo distribution of NEOs depends on NEO size. Instead, we find evidence that the disruption of NEOs at small perihelion distances found in Granvik et al. (2016) occurs preferentially for dark NEOs. The interval between km-sized bodies striking the Earth should occur on average once every 750,000 years. Low and high albedo NEOs are well mixed in orbital space, but a trend remains with higher albedo objects being at smaller semimajor axes and lower albedo objects more likely found at larger semimajor axes.

preprint2020arXiv

Effects of early intense bombardment on megaregolith evolution and on lunar (and planetary) surface samples

Impact rates in the first 500 Myr of the solar system are critical to an understanding of lunar geological history, but they have been controversial. The widely accepted, post-Apollo paradigm of early lunar impact cratering (ca. 1975-2014) proposed very low or negligible impact cratering in the period from accretion (>4.4 Ga) to about 4.0 Ga ago, followed by a 170-million-year-long spike of cataclysmic cratering, during which most prominent multi-ring impact basins formed at age of about 3.9 Ga. More recent dynamical models suggest very early intense impact rates, declining throughout the period from accretion until an age of about 3.0 Ga. These models remove the basin-forming spike. This shift has important consequences on megaregolith evolution and properties of rock samples that can be collected on the lunar surface today. We adopt the Morbidelli et al. (2018) "accretion tail" model of early intense bombardment, declining as a function of time. We find effects differing from the previous models: early crater saturation and supersaturation; disturbance of magma ocean solidification; deep early megaregolith; and erosive destruction of the earliest multi-ring basins, their impact melts, and their ejecta blankets. Our results explain observations such as differences in numbers of early lunar impact melts vs. numbers of early igneous crustal rocks, highland breccias containing impact melts as old as 4.35 Ga, absence of a 170 Myr-long spike in impact melt ages at 3.9 Ga among lunar and asteroidal meteorites, and GRAIL observations of lunar crustal structure.

preprint2020arXiv

Formation of Giant Planet Satellites

Recent analyses have shown that the concluding stages of giant planet formation are accompanied by the development of large-scale meridional flow of gas inside the planetary Hill sphere. This circulation feeds a circumplanetary disk that viscously expels gaseous material back into the parent nebula, maintaining the system in a quasi-steady state. Here we investigate the formation of natural satellites of Jupiter and Saturn within the framework of this newly outlined picture. We begin by considering the long-term evolution of solid material, and demonstrate that the circumplanetary disk can act as a global dust trap, where $s_{\bullet}\sim0.1-10\,$mm grains achieve a hydrodynamical equilibrium, facilitated by a balance between radial updraft and aerodynamic drag. This process leads to a gradual increase in the system's metallicity, and eventually culminates in the gravitational fragmentation of the outer regions of the solid sub-disk into $\mathcal{R}\sim100\,$km satellitesimals. Subsequently, satellite conglomeration ensues via pairwise collisions, but is terminated when disk-driven orbital migration removes the growing objects from the satellitesimal feeding zone. The resulting satellite formation cycle can repeat multiple times, until it is brought to an end by photo-evaporation of the parent nebula. Numerical simulations of the envisioned formation scenario yield satisfactory agreement between our model and the known properties of the Jovian and Saturnian moons.

preprint2020arXiv

Isotopic evolution of the inner Solar System inferred from molybdenum isotopes in meteorites

The fundamentally different isotopic compositions of non-carbonaceous (NC) and carbonaceous (CC) meteorites reveal the presence of two distinct reservoirs in the solar protoplanetary disk that were likely separated by Jupiter. However, the extent of material exchange between these reservoirs, and how this affected the composition of the inner disk are not known. Here we show that NC meteorites display broadly correlated isotopic variations for Mo, Ti, Cr, and Ni, indicating the addition of isotopically distinct material to the inner disk. The added material resembles bulk CC meteorites and Ca-Al-rich inclusions in terms of its enrichment in neutron-rich isotopes, but unlike the latter materials is also enriched in s-process nuclides. The comparison of the isotopic composition of NC meteorites with the accretion ages of their parent bodies reveals that the isotopic variations within the inner disk do not reflect a continuous compositional change through the addition of CC dust, indicating an efficient separation of the NC and CC reservoirs and limited exchange of material between the inner and outer disk. Instead, the isotopic variations among NC meteorites more likely record a rapidly changing composition of the disk during infall from the Sun's parental molecular cloud, where each planetesimal locks the instant composition of the disk when it forms. A corollary of this model is that late-formed planetesimals in the inner disk predominantly accreted from secondary dust that was produced by collisions among pre-existing NC planetesimals.

preprint2020arXiv

Planet formation by pebble accretion in ringed disks

Context: Pebble accretion is expected to be the dominant process for the formation of massive solid planets, such as the cores of giant planets and super-Earths. So, far, this process has been studied under the assumption that dust coagulates and drifts throughout the full protoplanetary disk. However, observations show that many disks are structured in rings that may be due to pressure maxima, preventing the global radial drift of the dust. Aims: We study how the pebble-accretion paradigm changes if the dust is confined in a ring. Results: Planet Type-I migration is stopped in a ring, but not necessarily at its center. If the entropy-driven corotation torque is desaturated, the planet is located in a region with a low density of dust, which severely limits its accretion rate. If instead the planet is near the ring's center, its accretion rate can be similar to the one it would have in a classic (ring-less) disk of equivalent dust density. However, the growth rate of the planet is limited by the diffusion of dust in the ring and the final planet's mass is bounded by the total ring's mass. The DSHARP rings are too far from the star to allow the formation of massive planets within the disk's lifetime. However, a similar ring rescaled to 5 AU could lead to the formation of a planet incorporating the full ring's mass in less than 1/2 My. Conclusions: The existence of rings may not be an obstacle to planet formation by pebble-accretion. However, for accretion to be effective the resting position of the planet has to be relatively near the ring's center and the ring needs to be not too far from the central star. The formation of planets in rings can explain the existence of giant planets with core masses smaller than the so-called pebble isolation mass.

preprint2020arXiv

Planet formation: key mechanisms and global models

Models of planet formation are built on underlying physical processes. In order to make sense of the origin of the planets we must first understand the origin of their building blocks. This review comes in two parts. The first part presents a detailed description of six key mechanisms of planet formation: 1) The structure and evolution of protoplanetary disks 2) The formation of planetesimals 3) Accretion of protoplanets 4) Orbital migration of growing planets 5) Gas accretion and giant planet migration 6) Resonance trapping during planet migration. While this is not a comprehensive list, it includes processes for which our understanding has changed in recent years or for which key uncertainties remain. The second part of this review shows how global models are built out of planet formation processes. We present global models to explain different populations of known planetary systems, including close-in small/low-mass planets (i.e., super-Earths), giant exoplanets, and the Solar System's planets. We discuss the different sources of water on rocky exoplanets, and use cosmochemical measurements to constrain the origin of Earth's water. We point out the successes and failings of different models and how they may be falsified. Finally, we lay out a path for the future trajectory of planet formation studies.

preprint2020arXiv

The onset of instability in resonant chains

There is evidence that most chains of mean motion resonances of type $k$:$k-1$ among exoplanets become unstable once the dissipative action from the gas is removed from the system, particularly for large $N$ (the number of planets) and $k$ (indicating how compact the chain is). We present a novel dynamical mechanism that can explain the origin of these instabilities and thus the dearth of resonant systems in the exoplanet sample. It relies on the emergence of secondary resonances between a fraction of the synodic frequency $2 π(1/P_1-1/P_2)$ and the libration frequencies in the mean motion resonance. These secondary resonances excite the amplitudes of libration of the mean motion resonances thus leading to an instability. We detail the emergence of these secondary resonances by carrying out an explicit perturbative scheme to second order in the planetary masses and isolating the harmonic terms that are associated with them. Focusing on the case of three planets in the 3:2 -- 3:2 mean motion resonance as an example, a simple but general analytical model of one of these resonances is obtained which describes the initial phase of the activation of one such secondary resonance. The dynamics of the excited system is also briefly described. This scheme shows how one can obtain analytical insight into the emergence of these resonances, and into the dynamics that they trigger. Finally, a generalisation of this dynamical mechanism is obtained for arbitrary $N$ and $k$. This leads to an explanation of previous numerical experiments on the stability of resonant chains, showing why the critical planetary mass allowed for stability decreases with increasing $N$ and $k$.

preprint2019arXiv

A record of the final phase of giant planet migration fossilized in the asteroid belt's orbital structure

The asteroid belt is characterized by an extreme low total mass of material on dynamically excited orbits. The Nice Model explains many peculiar qualities of the solar system, including the belt's excited state, by invoking an orbital instability between the outer planets. However, previous studies of the Nice Model's effect on the belt's structure struggle to reproduce the innermost asteroids' orbital inclination distribution. Here, we show how the final phase of giant planet migration sculpts the asteroid belt, in particular its inclination distribution. As interactions with leftover planetesimals cause Saturn to move away from Jupiter, its rate of orbital precession slows as the two planets' mutual interactions weaken. When the planets approach their modern separation, where Jupiter completes just short of five orbits for every two of Saturn's, Jupiter's eccentric forcing on Saturn strengthens. We use numerical simulations to show that the absence of asteroids with orbits that precess between 24-28 arcsec/yr is related to the inclination problem. As Saturn's precession speeds back up, high inclination asteroids are excited on to planet crossing orbits and removed from the inner main belt. Through this process, the asteroid belt's orbital structure is reshaped, leading to markedly improved simulation outcomes.

preprint2019arXiv

Formation of compact systems of super-Earths via dynamical instabilities and giant impacts

NASA's Kepler mission discovered $\sim700$ planets in multi-planet systems containing 3 or more transiting bodies, many of which are super-Earths and mini-Neptunes in compact configurations. Using $N$-body simulations, we examine the in situ, final stage assembly of multi-planet systems via the collisional accretion of protoplanets. Our initial conditions are constructed using a subset of the Kepler 5-planet systems as templates. Two different prescriptions for treating planetary collisions are adopted. The simulations address numerous questions: do the results depend on the accretion prescription?; do the resulting systems resemble the Kepler systems, and do they reproduce the observed distribution of planetary multiplicities when synthetically observed?; do collisions lead to significant modification of protoplanet compositions, or to stripping of gaseous envelopes?; do the eccentricity distributions agree with those inferred for the Kepler planets? We find the accretion prescription is unimportant in determining the outcomes. The final planetary systems look broadly similar to the Kepler templates adopted, but the observed distributions of planetary multiplicities or eccentricities are not reproduced, because scattering does not excite the systems sufficiently. In addition, we find that $\sim 1$% of our final systems contain a co-orbital planet pair in horseshoe or tadpole orbits. Post-processing the collision outcomes suggests they would not significantly change the ice fractions of initially ice-rich protoplanets, but significant stripping of gaseous envelopes appears likely. Hence, it may be difficult to reconcile the observation that many low mass Kepler planets have H/He envelopes with an in situ formation scenario that involves giant impacts after dispersal of the gas disc.

preprint2019arXiv

Kuiper belt: formation and evolution

This chapter reviews accretion models for Kuiper belt objects (KBOs), discussing in particular the compatibility of the observed properties of the KBO population with the streaming instability paradigm. Then it discusses how the dynamical structure of the KBO population, including the formation of its 5 sub-components (cold, hot, resonant, scattered and fossilized), can be quantitatively understood in the framewok of the giant planet instability. We also establish the connections between the KBO population and the Trojans of Jupiter and Neptune, the irregular satellites of all giant planets, the Oort cloud and the D-type main belt asteroids. Finally, we discuss the collisional evolution of the KBO population, arguing that the current size-frequency distribution below 100~km in size has been achieved as a collisional equilibrium in a few tens of My inside the original massive trans-Neptunain disk, possibly with the exception of the cold population sub-component.

preprint2018arXiv

Solar System Formation in the Context of Extra-Solar Planets

Exoplanet surveys have confirmed one of humanity's (and all teenagers') worst fears: we are weird. If our Solar System were observed with present-day Earth technology -- to put our system and exoplanets on the same footing -- Jupiter is the only planet that would be detectable. The statistics of exo-Jupiters indicate that the Solar System is unusual at the ~1% level among Sun-like stars (or ~0.1% among all stars). But why are we different? Successful formation models for both the Solar System and exoplanet systems rely on two key processes: orbital migration and dynamical instability. Systems of close-in super-Earths or sub-Neptunes require substantial radial inward motion of solids either as drifting mm- to cm-sized pebbles or migrating Earth-mass or larger planetary embryos. We argue that, regardless of their formation mode, the late evolution of super-Earth systems involves migration into chains of mean motion resonances, generally followed by instability when the disk dissipates. This pattern is likely also ubiquitous in giant planet systems. We present three models for inner Solar System formation -- the low-mass asteroid belt, Grand Tack, and Early Instability models -- each invoking a combination of migration and instability. We identify bifurcation points in planetary system formation. We present a series of events to explain why our Solar System is so weird. Jupiter's core must have formed fast enough to quench the growth of Earth's building blocks by blocking the flux of inward-drifting pebbles. The large Jupiter/Saturn mass ratio is rare among giant exoplanets but may be required to maintain Jupiter's wide orbit. The giant planets' instability must have been gentle, with no close encounters between Jupiter and Saturn, also unusual in the larger (exoplanet) context. Our Solar System system is thus the outcome of multiple unusual, but not unheard of, events.

preprint2016arXiv

Collisionless Encounters and the Origin of the Lunar Inclination

The Moon is generally thought to have formed from the debris ejected by the impact of a planet-sized object with the proto-Earth towards the end of planetary accretion. Modeling of the impact process predicts that the lunar material was disaggregated into a circumplanetary disk and that lunar accretion subsequently placed the Moon in a near equatorial orbit. Forward integration of the lunar orbit from this initial state predicts a modern inclination at least an order of magnitude smaller than the lunar value, a long-standing discrepancy known as the lunar inclination problem. Here we show that the modern lunar orbit provides a sensitive record of gravitational interactions with Earth-crossing planetesimals not yet accreted at the time of the Moon-forming event. The excited lunar orbit can naturally be reproduced via interaction with a small quantity of mass (corresponding to 0.0075-0.015 ME eventually accreted to the Earth) carried by a few bodies, consistent with constraints and models of late accretion. While the process has a stochastic element, the observed value of the lunar inclination is among the most likely outcomes for a wide range of parameters. The excitation of the lunar orbit is most readily reproduced via collisionless encounters with strong tidal dissipation on the early Earth. This mechanism obviates the necessity of previously proposed (but idealized) excitation mechanisms and can constrain the planet formation context of the Moon-forming event and the pristineness of the dynamical state of the Earth-Moon system.

preprint2016arXiv

Evolution of Protoplanetary Discs with Magnetically Driven Disc Winds

Aims: We investigate the evolution of protoplanetary discs (PPDs hereafter) with magnetically driven disc winds and viscous heating. Methods: We consider an initially massive disc with ~0.1 Msun to track the evolution from the early stage of PPDs. We solve the time evolution of surface density and temperature by taking into account viscous heating and the loss of the mass and the angular momentum by the disc winds within the framework of a standard alpha model for accretion discs. Our model parameters, turbulent viscosity, disc wind mass loss, and disc wind torque, which are adopted from local magnetohydrodynamical simulations and constrained by the global energetics of the gravitational accretion, largely depends on the physical condition of PPDs, particularly on the evolution of the vertical magnetic flux in weakly ionized PPDs. Results: Although there are still uncertainties concerning the evolution of the vertical magnetic flux remaining, surface densities show a large variety, depending on the combination of these three parameters, some of which are very different from the surface density expected from the standard accretion. When a PPD is in a "wind-driven accretion" state with the preserved vertical magnetic field, the radial dependence of the surface density can be positive in the inner region <1-10 au. The mass accretion rates are consistent with observations, even in the very low level of magnetohydrodynamical turbulence. Such a positive radial slope of the surface density gives a great impact on planet formation because (i)it inhibits the inward drift or even results in the outward drift of pebble/boulder-sized solid bodies, and (ii) it also makes the inward type-I migration of proto-planets slower or even reversed. Conclusions: The variety of our calculated PPDs should yield a wide variety of exoplanet systems.

preprint2016arXiv

Highly siderophile elements were stripped from Earth's mantle by iron sulfide segregation

Highly siderophile elements (HSEs) are strongly depleted in the bulk silicate Earth (BSE) but are present in near-chondritic relative abundances. The conventional explanation is that the HSEs were stripped from the mantle by the segregation of metal during core formation but were added back in near-chondritic proportions by late accretion, after core formation had ceased. Here we show that metal-silicate equilibration and segregation during Earth's core formation actually increased HSE mantle concentrations because HSE partition coefficients are relatively low at the high pressures of core formation within Earth. The pervasive exsolution and segregation of iron sulfide liquid from silicate liquid (the "Hadean matte") stripped magma oceans of HSEs during cooling and crystallization, before late accretion, and resulted in slightly suprachondritic palladium/iridium and ruthenium/iridium ratios.

preprint2016arXiv

Impact-induced melting during accretion of the Earth

Because of the high energies involved, giant impacts that occur during planetary accretion cause large degrees of melting. The depth of melting in the target body after each collision determines the pressure and temperature conditions of metal-silicate equilibration and thus geochemical fractionation that results from core-mantle differentiation. The accretional collisions involved in forming the terrestrial planets of the inner Solar System have been calculated by previous studies using N-body accretion simulations. Here we use the output from such simulations to determine the volumes of melt produced and thus the pressure and temperature conditions of metal-silicate equilibration, after each impact, as Earth-like planets accrete. For these calculations a parametrised melting model is used that takes impact velocity, impact angle and the respective masses of the impacting bodies into account. The evolution of metal-silicate equilibration pressures (as defined by evolving magma ocean depths) during Earth's accretion depends strongly on the lifetime of impact-generated magma oceans compared to the time interval between large impacts. In addition, such results depend on starting parameters in the N-body simulations, such as the number and initial mass of embryos. Thus, there is the potential for combining the results, such as those presented here, with multistage core formation models to better constrain the accretional history of the Earth.

preprint2016arXiv

Magnitude and timing of the giant planet instability: A reassessment from the perspective of the asteroid belt

It is generally accepted today that our solar system has undergone a phase during which the orbits of the giant planets became very unstable. In recent years, many studies have identified traces of this event and have provided reasonable justification for this occurrence. The magnitude (in terms of orbital variation) and the timing of the instability though (early or late with respect to the dispersal of the gas disk) still remains an open debate. The terrestrial planets seem to set a strict constraint: either the giant planet instability happened early, while the terrestrial planets were still forming, or the orbits of Jupiter and Saturn had to separate from each other impulsively, with a large enough `jump' in semimajor axis (Brasser et al. 2009; Kaib and Chambers 2016) for the terrestrial planets to remain stable. Because a large orbital jump is a low probability event, the early instability hypothesis seems to be favored. However, the asteroid belt would also evolve in a different way, assuming different instability amplitudes. These two constraints need to match each other in order to favor one scenario over the other. Considering an initially dynamically cold disk of asteroids, Morbidelli et al. (2010) concluded that a comparably large jump is needed to reconstruct the current asteroid belt. Here we confirm the same conclusion, but considering an asteroid population already strongly excited in eccentricity, such as that produced in the Grand Tack scenario (Walsh et al. 2011). Because the asteroids existed since the time of removal of the gas disk, unlike the terrestrial planets, this constraint on the width of the giant planet jump is valid regardless of whether the instability happened early or late. Hence, at this stage, assuming an early instability does not appear to provide any advantage in terms of the probabilistic reconstruction of the solar system structure.

preprint2016arXiv

Oxygen isotopic evidence for vigorous mixing during the Moon-forming Giant Impact

Earth and Moon are shown here to be composed of oxygen isotope reservoirs that are indistinguishable, with a difference in Δ"17O of -1 +/- 5ppm (2se). Based on these data and our new planet formation simulations that include a realistic model for oxygen isotopic reservoirs, our results favor vigorous mixing during the giant impact and therefore a high-energy high- angular-momentum impact. The results indicate that the late veneer impactors had an average Δ"17O within approximately 1 per mil of the terrestrial value, suggesting that these impactors were water rich.

preprint2016arXiv

Planet Formation Imager (PFI): science vision and key requirements

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

preprint2016arXiv

The Asteroid Belt as a Relic From a Chaotic Early Solar System

The orbital structure of the asteroid belt holds a record of the Solar System's dynamical history. The current belt only contains ${\rm \sim 10^{-3}}$ Earth masses yet the asteroids' orbits are dynamically excited, with a large spread in eccentricity and inclination. In the context of models of terrestrial planet formation, the belt may have been excited by Jupiter's orbital migration. The terrestrial planets can also be reproduced without invoking a migrating Jupiter; however, as it requires a severe mass deficit beyond Earth's orbit, this model systematically under-excites the asteroid belt. Here we show that the orbits of the asteroids may have been excited to their current state if Jupiter and Saturn's early orbits were chaotic. Stochastic variations in the gas giants' orbits cause resonances to continually jump across the main belt and excite the asteroids' orbits on a timescale of tens of millions of years. While hydrodynamical simulations show that the gas giants were likely in mean motion resonance at the end of the gaseous disk phase, small perturbations could have driven them into a chaotic but stable state. The gas giants' current orbits were achieved later, during an instability in the outer Solar System. Although it is well known that the present-day Solar System exhibits chaotic behavior, our results suggest that the early Solar System may also have been chaotic.

preprint2016arXiv

The inclination of the planetary system relative to the solar equator may be explained by the presence of Planet 9

We evaluate the effects of a distant planet, commonly known as planet 9, on the dynamics of the giant planets of the Solar System. We find that, given the large distance of planet 9, the dynamics of the inner giant planets can be decomposed into a classic Lagrange-Laplace dynamics relative to their own invariant plane (the plane orthogonal to their total angular momentum vector) and a slow precession of said plane relative to the total angular momentum vector of the Solar System, including planet 9. Under some specific configurations for planet 9, this precession can explain the current tilt of approximately 6 degrees between the invariant plane of the giant planets and the solar equator. An analytical model is developed to map the evolution of the inclination of the inner giant planets' invariant plane as a function of the planet 9's mass, inclination, eccentricity and semimajor axis, and some numerical simulations of the equations of motion of the giant planets and planet 9 are performed to validate our analytical approach. The longitude of the ascending node of planet 9 is found to be linked to the longitude of the ascending node of the giant planets' invariant plane, which also constrain the longitude of the node of planet 9 on the ecliptic. Some of the planet 9 configurations that allow explaining the current solar tilt are compatible with those proposed to explain the orbital confinement of the most distant Kuiper belt objects. Thus, this work on the one hand gives an elegant explanation for the current tilt between the invariant plane of the inner giant planets and the solar equator and, on the other hand, adds new constraints to the orbital elements of planet 9.

preprint2016arXiv

The radial dependence of pebble accretion rates: A source of diversity in planetary systems I. Analytical formulation

Context. The classical "planetesimal" accretion scenario for the formation of planets has recently evolved with the idea that "pebbles", centimeter- to meter-sized icy grains migrating in protoplanetary disks, can control planetesimal and/or planetary growth. Aims. We investigate how pebble accretion depends on disk properties and affects the formation of planetary systems Methods. We construct analytical models of pebble accretion onto planetary embryos that consistently account for the mass and orbital evolution of the pebble flow and reflect disk structure. Results. We derive simple formulas for pebble accretion rates in the so-called "settling" regime for planetary embryos with more than 100 km in size. For relatively smaller embryos or in outer disk regions, the accretion mode is 3D, meaning that the thickness of the pebble flow must be taken into account, and resulting in an accretion rate that is independent of the embryo mass. For larger embryos or in inner regions, the accretion is in a 2D mode, i.e., the pebble disk may be considered to be infinitely thin. We show that the radial dependence of the pebble accretion rata is different (even the sign of the power-law exponent changes) for different disk conditions such as the disk heating source (viscous heating or stellar irradiation), the drag law (Stokes or Epstein, and weak or strong coupling), and in the 2D or 3D accretion modes. We also discuss the effect of the sublimation/destruction of icy pebbles inside the snow line. Conclusions. Pebble accretion easily produces a large diversity of planetary systems. In other words, to correctly infer the results of planet formation through pebble accretion, detailed prescriptions of disk evolution and pebble growth, sublimation/destruction and migration are required.

preprint2016arXiv

Yarkovsky V-shape identification of asteroid families

There are only a few known main belt (MB) asteroid families with ages greater than 2 Gyr (Brož et al. 2013, Spoto et al. 2015). Estimates based on the family producing collision rate suggest that the lack of >2 Gyr-old families may be due to a selection bias in current techniques used to identify families. Family fragments disperse in their orbital elements, semi-major axis, $a$, eccentricity, $e$, and inclination, $i$, due to secular resonances, close encounters with massive asteroids and the non-gravitational Yarkovsky force. This causes the family fragments to be indistinguishable from the background of the main belt making them more difficult to identify with the hierarchical clustering method (HCM) with increasing family age. The discovery of the Eulalia and new Polana families in the inner belt relied on new techniques because Yarkovsky spreading made them too disperse to be identified using the classical HCM. The technique used to discover the new Polana and Eulalia families is modified to identify asteroid families by searching for correlations between $a$ and asteroid diameter, $D$, or absolute magnitude, $H$. A group of asteroids is identified as a collisional family if its boundary in the $a$ vs. 1/$D$ or $a$ vs. $H$ planes has a characteristic V-shape which is due to the size dependent Yarkovsky spreading. The V-shape boundary is identified with two separate techniques. The first identifies a border by measuring the drop between the number of objects inside and outside of the border. The second identifies the V-shape border by measuring the number density of objects in $a$ vs. 1/$D$,$H$ space. The V-shape techniques are demonstrated on known families Erigone, Vesta, Koronis, and families difficult to identify by HCM such as Flora, Baptistina, new Polana, Eulalia and Karin. Application of the technique in a search for >2 Gyr-old families throughout the MB is discussed.

preprint2015arXiv

A reassessment of the in situ formation of close-in super-Earths

A large fraction of stars host one or multiple close-in super-Earth planets. There is an active debate about whether these planets formed in situ or at greater distances from the central star and migrated to their current position. It has been shown that part of their observed properties (e.g., eccentricity distribution) can be reproduced by N-body simulations of in situ formation starting with a population of protoplanets of high masses and neglecting the effects of the disk gas. We plan to reassess the in situ formation of close-in super-Earths through more complete simulations. We performed N-body simulations of a population of small planetary embryos and planetesimals that include the effects of disk-planet interactions (e.g., eccentricity damping, type I migration). In addition, we also consider the accretion of a primitive atmosphere from a protoplanetary disk. We find that planetary embryos grow very quickly well before the gas dispersal, and thus undergo rapid inward migration, which means that one cannot neglect the effects of a gas disk when considering the in-situ formation of close-in super-Earths. Owing to their rapid inward migration, super-Earths reach a compact configuration near the disk's inner edge whose distribution of orbital parameters matches the observed close-in super-Earths population poorly. On the other hand, simulations including eccentricity damping, but no type I migration, reproduce the observed distributions better. Including the accretion of an atmosphere does not help reproduce the bulk architecture of observations. Interestingly, we find that the massive embryos can migrate inside the disk edge while capturing only a moderately massive hydrogen/helium atmosphere. By this process they avoid becoming giant planets. The bulk of close-in super-Earths cannot form in situ, unless type I migration is suppressed in the entire disk inside 1 AU.

preprint2015arXiv

Accretion of Uranus and Neptune from inward-migrating planetary embryos blocked by Jupiter and Saturn

Reproducing Uranus and Neptune remains a challenge for simulations of solar system formation. The ice giants' peculiar obliquities suggest that they both suffered giant collisions during their formation. Thus, there must have been an epoch of accretion dominated by collisions among large planetary embryos in the primordial outer solar system. We test this idea using N-body numerical simulations including the effects of a gaseous protoplanetary disk. One strong constraint is that the masses of the ice giants are very similar -- the Neptune/Uranus mass ratio is $\sim1.18$. We show that similar-size ice giants do indeed form by collisions between planetary embryos beyond Saturn. The fraction of successful simulations varies depending on the initial number of planetary embryos in the system, their individual and total masses. Similar-sized ice giants are consistently reproduced in simulations starting with 5-10 planetary embryos with initial masses of $\sim$3-6 ${\rm M_\oplus}$. We conclude that accretion from a population of planetary embryos is a plausible scenario for the origin of Uranus and Neptune.

preprint2015arXiv

Comets as collisional fragments of a primordial planetesimal disk

The Rosetta mission and its exquisite measurements have revived the debate on whether comets are pristine planetesimals or collisionally evolved objects. We investigate the collisional evolution experienced by the precursors of current comet nuclei during the early stages of the Solar System, in the context of the so-called "Nice Model". We consider two environments for the collisional evolution: (1) the trans-planetary planetesimal disk, from the time of gas removal until the disk was dispersed by the migration of the ice giants, and (2) the dispersing disk during the time that the scattered disk was formed. Simulations have been performed, using different methods in the two cases, to find the number of destructive collisions typically experienced by a comet nucleus of 2km radius. In the widely accepted scenario, where the dispersal of the planetesimal disk occurred at the time of the Late Heavy Bombardment about 4Gy ago, comet-sized planetesimals have a very small chance to survive against destructive collisions in the disk. On the extreme assumption that the disk was dispersed directly upon gas removal, there is a chance for a significant fraction of the planetesimals to remain intact. However, these survivors would still bear the marks of many non-destructive impacts. Thus, the Nice Model of Solar System evolution predicts that typical km-sized comet nuclei are predominantly fragments resulting from collisions experienced by larger parent bodies. An important goal for further research is to investigate, whether the observed properties of comet nuclei are compatible with such a collisional origin.

preprint2015arXiv

Gas giant planets as dynamical barriers to inward-migrating super-Earths

Planets of 1-4 times Earth's size on orbits shorter than 100 days exist around 30-50% of all Sun-like stars. In fact, the Solar System is particularly outstanding in its lack of "hot super-Earths" (or "mini-Neptunes"). These planets -- or their building blocks -- may have formed on wider orbits and migrated inward due to interactions with the gaseous protoplanetary disk. Here, we use a suite of dynamical simulations to show that gas giant planets act as barriers to the inward migration of super-Earths initially placed on more distant orbits. Jupiter's early formation may have prevented Uranus and Neptune (and perhaps Saturn's core) from becoming hot super-Earths. Our model predicts that the populations of hot super-Earth systems and Jupiter-like planets should be anti-correlated: gas giants (especially if they form early) should be rare in systems with many hot super-Earths. Testing this prediction will constitute a crucial assessment of the validity of the migration hypothesis for the origin of close-in super-Earths.

preprint2015arXiv

Highly Siderophile Elements in the Earth's Mantle as a Clock for the Moon-forming Impact

According to the generally accepted scenario, the last giant impact on the Earth formed the Moon and initiated the final phase of core formation by melting the Earth's mantle. A key goal of geochemistry is to date this event, but different ages have been proposed. Some argue for an early Moon-forming event, approximately 30 million years (Myr) after the condensation of the first solids in the Solar System, whereas others claim a date later than 50 Myr (and possibly as late as around 100 My) after condensation. Here we show that a Moon-forming event at 40 Myr after condensation, or earlier, is ruled out at a 99.9 per cent confidence level. We use a large number of N-body simulations to demonstrate a relationship between the time of the last giant impact on an Earth-like planet and the amount of mass subsequently added during the era known as Late Accretion. As the last giant impact is delayed, the late-accreted mass decreases in a predictable fashion. This relationship exists within both the classical scenario and the Grand Tack scenario of terrestrial planet formation, and it holds across a wide range of disk conditions. The concentration of highly siderophile elements (HSEs) in Earth's mantle constrains the mass of chondritic material added to Earth during Late Accretion. Using HSE abundance measurements, we determine a Moon-formation age of 95 +/- 32 Myr since condensation. The possibility exists that some late projectiles were differentiated and left an incomplete HSE record in Earth's mantle. Even in this case, various isotopic constraints strongly suggest that the late-accreted mass did not exceed 1 per cent of Earth's mass, and so the HSE clock still robustly limits the timing of the Moon-forming event to significantly later than 40 My after condensation.

preprint2015arXiv

New Paradigms For Asteroid Formation

Asteroids and meteorites provide key evidence on the formation of planetesimals in the Solar System. Asteroids are traditionally thought to form in a bottom-up process by coagulation within a population of initially km-scale planetesimals. However, new models challenge this idea by demonstrating that asteroids of sizes from 100 to 1000 km can form directly from the gravitational collapse of small particles which have organised themselves in dense filaments and clusters in the turbulent gas. Particles concentrate passively between eddies down to the smallest scales of the turbulent gas flow and inside large-scale pressure bumps and vortices. The streaming instability causes particles to take an active role in the concentration, by piling up in dense filaments whose friction on the gas reduces the radial drift compared to that of isolated particles. In this chapter we review new paradigms for asteroid formation and compare critically against the observed properties of asteroids as well as constraints from meteorites. Chondrules of typical sizes from 0.1 to 1 mm are ubiquitous in primitive meteorites and likely represent the primary building blocks of asteroids. Chondrule-sized particles are nevertheless tightly coupled to the gas via friction and are therefore hard to concentrate in large amounts in the turbulent gas. We review recent progress on understanding the incorporation of chondrules into the asteroids, including layered accretion models where chondrules are accreted onto asteroids over millions of years. We highlight in the end ten unsolved questions in asteroid formation where we expect that progress will be made over the next decade.

preprint2015arXiv

Spin-Spin Coupling in the Solar System

The richness of dynamical behavior exhibited by the rotational states of various solar system objects has driven significant advances in the theoretical understanding of their evolutionary histories. An important factor that determines whether a given object is prone to exhibiting non-trivial rotational evolution is the extent to which such an object can maintain a permanent aspheroidal shape, meaning that exotic behavior is far more common among the small body populations of the solar system. Gravitationally bound binary objects constitute a substantial fraction of asteroidal and TNO populations, comprising systems of triaxial satellites that orbit permanently deformed central bodies. In this work, we explore the rotational evolution of such systems with specific emphasis on quadrupole-quadrupole interactions, and show that for closely orbiting, highly deformed objects, both prograde and retrograde spin-spin resonances naturally arise. Subsequently, we derive capture probabilities for leading order commensurabilities and apply our results to the illustrative examples of (87) Sylvia and (216) Kleopatra asteroid systems. Cumulatively, our results suggest that spin-spin coupling may be consequential for highly elongated, tightly orbiting binary objects.

preprint2015arXiv

Suppression of type I migration by disk winds

Planets less massive than Saturn tend to rapidly migrate inward in protoplanetary disks. This is the so-called type I migration. Simulations attempting to reproduce the observed properties of exoplanets show that type I migration needs to be significantly reduced over a wide region of the disk for a long time. However, the mechanism capable of suppressing type I migration over a wide region has remained elusive. The recently found turbulence-driven disk winds offer new possibilities. We investigate the effects of disk winds on the disk profile and type I migration for a range of parameters that describe the strength of disk winds. We also examine the in situ formation of close-in super-Earths in disks that evolve through disk winds. The disk profile, which is regulated by viscous diffusion and disk winds, was derived by solving the diffusion equation. We carried out a number of simulations and plot here migration maps that indicate the type I migration rate. We also performed N-body simulations of the formation of close-in super-Earths from a population of planetesimals and planetary embryos. We define a key parameter, Kw, which determines the ratio of strengths between the viscous diffusion and disk winds. For a wide range of Kw, the type I migration rate is presented in migration maps. These maps show that type I migration is suppressed over the whole close-in region when the effects of disk winds are relatively strong (Kw < 100). From the results of N-body simulations, we see that type I migration is significantly slowed down assuming Kw = 40. We also show that the results of N-body simulations match statistical orbital distributions of close-in super-Earths.

preprint2015arXiv

Terrestrial Planet Formation Constrained by Mars and the Structure of the Asteroid Belt

Reproducing the large Earth/Mars mass ratio requires a strong mass depletion in solids within the protoplanetary disk between 1 and 3 AU. The Grand Tack model invokes a specific migration history of the giant planets to remove most of the mass initially beyond 1 AU and to dynamically excite the asteroid belt. However, one could also invoke a steep density gradient created by inward drift and pile-up of small particles induced by gas-drag, as has been proposed to explain the formation of close-in super Earths. Here we show that the asteroid belt's orbital excitation provides a crucial constraint against this scenario for the Solar System. We performed a series of simulations of terrestrial planet formation and asteroid belt evolution starting from disks of planetesimals and planetary embryos with various radial density gradients and including Jupiter and Saturn on nearly circular and coplanar orbits. Disks with shallow density gradients reproduce the dynamical excitation of the asteroid belt by gravitational self-stirring but form Mars analogs significantly more massive than the real planet. In contrast, a disk with a surface density gradient proportional to $r^{-5.5}$ reproduces the Earth/Mars mass ratio but leaves the asteroid belt in a dynamical state that is far colder than the real belt. We conclude that no disk profile can simultaneously explain the structure of the terrestrial planets and asteroid belt. The asteroid belt must have been depleted and dynamically excited by a different mechanism such as, for instance, in the Grand Tack scenario.

preprint2015arXiv

The Complex History of Trojan Asteroids

The Trojan asteroids provide a unique perspective on the history of Solar System. As a large population of small bodies, they record important gravitational interactions and dynamical evolution of the Solar System. In the past decade, significant advances have been made in understanding physical properties, and there has been a revolution in thinking about the origin of Trojans. The ice and organics generally presumed to be a significant part of Trojan compositions have yet to be detected directly, though low density of the binary system Patroclus (and possibly low density of the binary/moonlet system Hektor) is consistent with an interior ice component. By contrast, fine-grained silicates that appear to be similar to cometary silicates in composition have been detected, and a color bimodality may indicate distinct compositional groups among the Trojans. Whereas Trojans had traditionally been thought to have formed near 5 AU, a new paradigm has developed in which the Trojans formed in the proto-Kuiper Belt, and they were scattered inward and captured in the Trojan swarms as a result of resonant interactions of the giant planets. Whereas the orbital and population distributions of current Trojans are consistent with this origin scenario, there are significant differences between current physical properties of Trojans and those of Kuiper Belt objects. These differences may be indicative of surface modification due to the inward migration of objects that became the Trojans, but understanding of appropriate modification mechanisms is poor and would benefit from additional laboratory studies. Many open questions remain, and the future promises significant strides in our understanding of Trojans. The time is ripe for a spacecraft mission to the Trojans, to turn these objects into geologic worlds that can be studied in detail to unravel their complex history.

preprint2015arXiv

The Dynamical Evolution of the Asteroid Belt

The asteroid belt is the leftover of the original planetesimal population in the inner solar system. However, currently the asteroids have orbits with all possible values of eccentricities and inclinations compatible with long-term dynamical stability, whereas the initial planetesimal orbits should have been quasi-circular and almost co-planar. The total mass in the asteroid population is a small fraction of that existing primordially. Also, asteroids with different chemical/mineralogical properties are not ranked in an orderly manner with mean heliocentric distance as one could expect from the existence of a radial gradient of the temperature in the proto-planetary disk, but they are partially mixed. These properties show that the asteroid belt has been severely sculpted by one or a series of processes during its lifetime. This paper reviews the processes that have been proposed so far, discussing the properties that they explain and the problems that they are confronted with. Emphasis is paid to the interplay between the dynamical and the collisional evolution of the asteroid population, which allows the use of the size distribution to constrain the dynamical models. We divide the asteroid belt evolution into three phases. The first phase started during the lifetime of the gaseous proto-planetary disk, when the giant planets formed and presumably experienced large-scale migrations, and continued after the removal of the gas, during the build-up of the terrestrial planets. The second phase occurred after the removal of the gaseous proto-planetary disk and it became particularly lively for the asteroid belt when the giant planets suddenly changed their orbits, as a result of a mutual dynamical instability and the interaction with the trans-Neptunian planetesimal disk. The third phase covers the aftermath of the giant planet instability, until today.

preprint2014arXiv

Accretion and differentiation of the terrestrial planets with implications for the compositions of early-formed Solar System bodies and accretion of water

In order to test planetary accretion and differentiation scenarios, we integrated a multistage core-mantle differentiation model with N-body accretion simulations. Impacts between embryos and planetesimals result in magma ocean formation and episodes of core formation. The core formation model combines rigorous chemical mass balance with metal-silicate element partitioning data. The primary constraint on the combined model is the composition of the Earth's primitive mantle, the composition of the Martian mantle, and the mass fractions of the metallic cores of Earth and Mars. The model is refined by least squares minimization with up to five fitting parameters that consist of the metal-silicate equilibrium pressure and 1-4 parameters that define the starting compositions of primitive bodies. This integrated model has been applied to 6 Grand Tack simulations. Investigations of a broad parameter space indicate that: accretion of Earth was heterogeneous, metal-silicate equilibration pressures increase as accretion progresses and are 60-70% of core-mantle boundary pressures at the time of each impact, and a large fraction (70-100%) of the metal of impactor cores equilibrates with a small fraction of the silicate mantles of protoplanets during each core formation event. Acceptable fits to the Earth's mantle composition are obtained only when bodies that originated close to the Sun, at <0.9-1.2 AU, are highly reduced and those beyond this distance are increasingly oxidized. The FeO content of the Martian mantle depends critically on the heliocentric distance at which the Mars-forming embryo originated. Finally, the Earth's core is predicted to contain 8-9 wt% silicon, 2-4 wt% oxygen and 10-60 ppm hydrogen, whereas the Martian core is predicted to contain low concentrations (<1 wt%) of Si and O.

preprint2014arXiv

Chaotic Disintegration of the Inner Solar System

On timescales that greatly exceed an orbital period, typical planetary orbits evolve in a stochastic yet stable fashion. On even longer timescales, however, planetary orbits can spontaneously transition from bounded to unbound chaotic states. Large-scale instabilities associated with such behavior appear to play a dominant role in shaping the architectures of planetary systems, including our own. Here we show how such transitions are possible, focusing on the specific case of the long-term evolution of Mercury. We develop a simple analytical model for Mercury's dynamics and elucidate the origins of its short term stochastic behavior as well as of its sudden progression to unbounded chaos. Our model allows us to estimate the timescale on which this transition is likely to be triggered, i.e. the dynamical lifetime of the Solar System as we know it. The formulated theory is consistent with the results of numerical simulations and is broadly applicable to extrasolar planetary systems dominated by secular interactions. These results constitute a significant advancement in our understanding of the processes responsible for sculpting of the dynamical structures of generic planetary systems.

preprint2014arXiv

Earth-like Habitats in Planetary Systems

Understanding the concept of habitability is related to an evolutionary knowledge of the particular planet-in-question. Additional indications so-called "systemic aspects" of the planetary system as a whole governs a particular planet's claim on habitability. Here we focus on such systemic aspects and discuss their relevance to the formation of an 'Earth-like' habitable planet. We summarize our results obtained by lunar sample work and numerical models within the framework of the Research Alliance "Planetary Evolution and Life". We consider various scenarios which simulate the dynamical evolution of the Solar System and discuss the likelihood of forming an Earth-like world orbiting another star. Our model approach is constrained by observations of the modern Solar System and the knowledge of its history. Results suggest that the long-term presence of terrestrial planets is jeopardized due to gravitational interactions if giant planets are present. But habitability of inner rocky planets may be supported in those planetary systems hosting giant planets. Gravitational interactions within a complex multiple-body structure including giant planets may supply terrestrial planets with materials which formed in the colder region of the proto-planetary disk. During these processes, water, the prime requisite for habitability, is delivered to the inner system. This may occur either during the main accretion phase of terrestrial planets or via impacts during a post-accretion bombardment. Results for both processes are summarized and discussed with reference to the lunar crater record. Starting from a scenario involving migration of the giant planets this contribution discusses the delivery of water to Earth, the modification of atmospheres by impacts in a planetary system context and the likelihood of the existence of extrasolar Earth-like habitable worlds.

preprint2014arXiv

Late Accretion and the Late Veneer

The concept of Late Veneer has been introduced by the geochemical community to explain the abundance of highly siderophile elements in the Earth's mantle and their chondritic proportions relative to each other. However, in the complex scenario of Earth accretion, involving both planetesimal bombardment and giant impacts from chondritic and differentiated projectiles, it is not obvious what the "Late Veneer" actually corresponds to. In fact, the process of differentiation of the Earth was probably intermittent and there was presumably no well-defined transition between an earlier phase where all metal sunk into the core and a later phase in which the core was a closed entity separated from the mantle. In addition, the modellers of Earth accretion have introduced the concept of "Late Accretion", which refers to the material accreted by our planet after the Moon-forming event. Characterising Late Veneer, Late Accretion and the relationship between the two is the major goal of this chapter.

preprint2014arXiv

Lunar and Terrestrial Planet Formation in the Grand Tack Scenario

We present conclusions from a large number of N-body simulations of the giant impact phase of terrestrial planet formation. We focus on new results obtained from the recently proposed Grand Tack model, which couples the gas-driven migration of giant planets to the accretion of the terrestrial planets. The giant impact phase follows the oligarchic growth phase, which builds a bi-modal mass distribution within the disc of embryos and planetesimals. By varying the ratio of the total mass in the embryo population to the total mass in the planetesimal population and the mass of the individual embryos, we explore how different disc conditions control the final planets. The total mass ratio of embryos to planetesimals controls the timing of the last giant (Moon forming) impact and its violence. The initial embryo mass sets the size of the lunar impactor and the growth rate of Mars. After comparing our simulated outcomes with the actual orbits of the terrestrial planets (angular momentum deficit, mass concentration) and taking into account independent geochemical constraints on the mass accreted by the Earth after the Moon forming event and on the timescale for the growth of Mars, we conclude that the protoplanetary disc at the beginning of the giant impact phase must have had most of its mass in Mars-sized embryos and only a small fraction of the total disc mass in the planetesimal population. From this, we infer that the Moon forming event occurred between $\sim$60 and $\sim$130 My after the formation of the first solids, and was caused most likely by an object with a mass similar to that of Mars.

preprint2014arXiv

Outward migration of Jupiter and Saturn in 3:2 or 2:1 resonance in radiative disks: implications for the Grand Tack and Nice models

Embedded in the gaseous protoplanetary disk, Jupiter and Saturn naturally become trapped in 3:2 resonance and migrate outward. This serves as the basis of the Grand Tack model. However, previous hydrodynamical simulations were restricted to isothermal disks, with moderate aspect ratio and viscosity. Here we simulate the orbital evolution of the gas giants in disks with viscous heating and radiative cooling. We find that Jupiter and Saturn migrate outward in 3:2 resonance in modest-mass ($M_{disk} \approx M_{MMSN}$, where MMSN is the "minimum-mass solar nebula") disks with viscous stress parameter $α$ between $10^{-3}$ and $10^{-2} $. In disks with relatively low-mass ($M_{disk} \lesssim M_{MMSN}$) , Jupiter and Saturn get captured in 2:1 resonance and can even migrate outward in low-viscosity disks ($α\le 10^{-4}$). Such disks have a very small aspect ratio ($h\sim 0.02-0.03$) that favors outward migration after capture in 2:1 resonance, as confirmed by isothermal runs which resulted in a similar outcome for $h \sim 0.02$ and $α\le 10^{-4}$. We also performed N-body runs of the outer Solar System starting from the results of our hydrodynamical simulations and including 2-3 ice giants. After dispersal of the gaseous disk, a Nice model instability starting with Jupiter and Saturn in 2:1 resonance results in good Solar Systems analogs. We conclude that in a cold Solar Nebula, the 2:1 resonance between Jupiter and Saturn can lead to outward migration of the system, and this may represent an alternative scenario for the evolution of the Solar System.

preprint2014arXiv

Separating gas-giant and ice-giant planets by halting pebble accretion

In the Solar System giant planets come in two flavours: 'gas giants' (Jupiter and Saturn) with massive gas envelopes and 'ice giants' (Uranus and Neptune) with much thinner envelopes around their cores. It is poorly understood how these two classes of planets formed. High solid accretion rates, necessary to form the cores of giant planets within the life-time of protoplanetary discs, heat the envelope and prevent rapid gas contraction onto the core, unless accretion is halted. We find that, in fact, accretion of pebbles (~ cm-sized particles) is self-limiting: when a core becomes massive enough it carves a gap in the pebble disc. This halt in pebble accretion subsequently triggers the rapid collapse of the super-critical gas envelope. As opposed to gas giants, ice giants do not reach this threshold mass and can only bind low-mass envelopes that are highly enriched by water vapour from sublimated icy pebbles. This offers an explanation for the compositional difference between gas giants and ice giants in the Solar System. Furthermore, as opposed to planetesimal-driven accretion scenarios, our model allows core formation and envelope attraction within disc life-times, provided that solids in protoplanetary discs are predominantly in pebbles. Our results imply that the outer regions of planetary systems, where the mass required to halt pebble accretion is large, are dominated by ice giants and that gas-giant exoplanets in wide orbits are enriched by more than 50 Earth masses of solids.

preprint2014arXiv

Stellar irradiated discs and implications on migration of embedded planets II: accreting-discs

The strength and direction of migration of embedded low mass planets depends on the disc's structure. It has been shown that, in discs with viscous heating and radiative transport, the migration can be directed outwards. In this paper we investigate the influence of a constant dM/dt-flux through the disc, as well as the influence of the disc's metallicity on the disc's thermodynamics. We focus on dM/dt discs, which have a net mass flux through them. Utilizing the resulting disc structure, we determine the regions of outward migration in the disc. We perform numerical hydrosimulations of dM/dt discs with viscous heating, radiative cooling and stellar irradiation in 2D in the r-z-plane. We use the explicit/implicit hydrodynamical code FARGOCA that includes a full tensor viscosity and stellar irradiation, as well as a two temperature solver that includes radiation transport in the flux-limited diffusion approximation. The migration of embedded planets is studied by using torque formulae. For a disc of gas surface density and viscosity, we find that the discs thermal structure depends on the product Sigma_G νand the amount of heavy elements, while the migration of planets additionally to the mentioned quantities, depends on the amount of viscosity itself. As a result of this, the disc structure can not be approximated by simple power laws. During the lifetime of the disc, the structure of the disc changes significantly in a non-linear way in the inner parts. In the late stages of the discs evolution, outward migration is only possible if the metallicity of the disc is high. For low metallicity, planets would migrate inwards and could potentially be lost. The presented disc structures and migration maps have important consequences on the formation of planets, as they can give hints on the different formation mechanisms for different types of planets as a function of metallicity.

preprint2014arXiv

Stellar irradiated discs and implications on migration of embedded planets III: viscosity transitions

The migration strength and direction of embedded low-mass planets depends on the disc structure. In discs with an efficient radiative transport, the migration can be directed outwards for planets with more than 3-5 Earth masses. This is due to the entropy driven corotation torque, a process that extends the lifetimes of growing planetary embryos. We investigate the influence on the disc structure caused by a jump in the alpha parameter of the viscosity to model a dead-zone structure in the disc. We focus on M-dot discs, which have a constant net mass flux. Using the resulting disc structure, we investigate the consequences for the formation of planetesimals and determine the regions of outward migration for proto-planets. We performed numerical hydrosimulations of M-dot discs in the r-z-plane. We used the explicit/implicit hydrodynamical code FARGOCA that includes a full tensor viscosity and stellar irradiation as well as a two-temperature solver that includes radiation transport in the flux-limited diffusion approximation. Viscosity transitions inside the disc create transitions in density that stop inward migration for small planets through the so-called "planet trap" mechanism. This mechanism also works for planets down to M_P > 0.5M_E. Additionally, the viscosity transitions change the pressure gradient in the disc, which facilitates planetesimal formation via the streaming instability. However, a very steep transition in viscosity is needed to achieve in a pressure bump in the disc. The transition in viscosity facilitates planetesimal formation and can stop the migration of small-mass planets (M_P>0.5M_E), but still does not halt inward migration of smaller planets and planetesimals that are affected by gas drag. A very steep, probably unrealistic viscosity gradient is needed to trap planets of smaller masses and halt gas-drag-driven planetesimal migration at a pressure bump.

preprint2014arXiv

Terrestrial Planet Formation at Home and Abroad

We review the state of the field of terrestrial planet formation with the goal of understanding the formation of the inner Solar System and low-mass exoplanets. We review the dynamics and timescales of accretion from planetesimals to planetary embryos and from embryos to terrestrial planets. We discuss radial mixing and water delivery, planetary spins and the importance of parameters regarding the disk and embryo properties. Next, we connect accretion models to exoplanets. We first explain why the observed hot Super Earths probably formed by in situ accretion or inward migration. We show how terrestrial planet formation is altered in systems with gas giants by the mechanisms of giant planet migration and dynamical instabilities. Standard models of terrestrial accretion fail to reproduce the inner Solar System. The "Grand Tack" model solves this problem using ideas first developed to explain the giant exoplanets. Finally, we discuss whether most terrestrial planet systems form in the same way as ours, and highlight the key ingredients missing in the current generation of simulations.

preprint2014arXiv

Terrestrial Planet Formation in the Presence of Migrating Super-earths

Super-Earths with orbital periods less than 100 days are extremely abundant around Sun-like stars. It is unlikely that these planets formed at their current locations. Rather, they likely formed at large distances from the star and subsequently migrated inward. Here we use N-body simulations to study the effect of super-Earths on the accretion of rocky planets. In our simulations, one or more super-Earths migrates inward through a disk of planetary embryos and planetesimals embedded in a gaseous disk. We tested a wide range of migration speeds and configurations. Fast-migrating super-Earths ($τ_{mig} \sim$0.01-0.1 Myr) only have a modest effect on the protoplanetary embryos and planetesimals. Sufficient material survives to form rocky, Earth-like planets on orbits exterior to the super-Earths'. In contrast, slowly migrating super-Earths shepherd rocky material interior to their orbits and strongly deplete the terrestrial planet-forming zone. In this situation any Earth-sized planets in the habitable zone are extremely volatile-rich and are therefore probably not Earth-like.

preprint2014arXiv

The Absolute Magnitude Distribution of Kuiper Belt Objects

Here we measure the absolute magnitude distributions (H-distribution) of the dynamically excited and quiescent (hot and cold) Kuiper Belt objects (KBOs), and test if they share the same H-distribution as the Jupiter Trojans. From a compilation of all useable ecliptic surveys, we find that the KBO H-distributions are well described by broken power-laws. The cold population has a bright-end slope, $α_{\textrm{1}}=1.5_{-0.2}^{+0.4}$, and break magnitude, $H_{\textrm{B}}=6.9_{-0.2}^{+0.1}$ (r'-band). The hot population has a shallower bright-end slope of, $α_{\textrm{1}}=0.87_{-0.2}^{+0.07}$, and break magnitude $H_{\textrm{B}}=7.7_{-0.5}^{+1.0}$. Both populations share similar faint end slopes of $α_2\sim0.2$. We estimate the masses of the hot and cold populations are $\sim0.01$ and $\sim3\times10^{-4} \mbox{ M$_{\bigoplus}$}$. The broken power-law fit to the Trojan H-distribution has $α_\textrm{1}=1.0\pm0.2$, $α_\textrm{2}=0.36\pm0.01$, and $H_{\textrm{B}}=8.3$. The KS test reveals that the probability that the Trojans and cold KBOs share the same parent H-distribution is less than 1 in 1000. When the bimodal albedo distribution of the hot objects is accounted for, there is no evidence that the H-distributions of the Trojans and hot KBOs differ. Our findings are in agreement with the predictions of the Nice model in terms of both mass and H-distribution of the hot and Trojan populations. Wide field survey data suggest that the brightest few hot objects, with $H_{\textrm{r'}}\lesssim3$, do not fall on the steep power-law slope of fainter hot objects. Under the standard hierarchical model of planetesimal formation, it is difficult to account for the similar break diameters of the hot and cold populations given the low mass of the cold belt.

preprint2014arXiv

The Grand Tack model: a critical review

The `Grand Tack' model proposes that the inner Solar System was sculpted by the giant planets' orbital migration in the gaseous protoplanetary disk. Jupiter first migrated inward then Jupiter and Saturn migrated back outward together. If Jupiter's turnaround or "tack" point was at ~1.5 AU the inner disk of terrestrial building blocks would have been truncated at ~1 AU, naturally producing the terrestrial planets' masses and spacing. During the gas giants' migration the asteroid belt is severely depleted but repopulated by distinct planetesimal reservoirs that can be associated with the present-day S and C types. The giant planets' orbits are consistent with the later evolution of the outer Solar System. Here we confront common criticisms of the Grand Tack model. We show that some uncertainties remain regarding the Tack mechanism itself; the most critical unknown is the timing and rate of gas accretion onto Saturn and Jupiter. Current isotopic and compositional measurements of Solar System bodies -- including the D/H ratios of Saturn's satellites -- do not refute the model. We discuss how alternate models for the formation of the terrestrial planets each suffer from an internal inconsistency and/or place a strong and very specific requirement on the properties of the protoplanetary disk. We conclude that the Grand Tack model remains viable and consistent with our current understanding of planet formation. Nonetheless, we encourage additional tests of the Grand Tack as well as the construction of alternate models.

preprint2014arXiv

The structure of protoplanetary discs around evolving young stars

The formation of planets with gaseous envelopes takes place in protoplanetary accretion discs on time-scales of several millions of years. Small dust particles stick to each other to form pebbles, pebbles concentrate in the turbulent flow to form planetesimals and planetary embryos and grow to planets, which undergo substantial radial migration. All these processes are influenced by the underlying structure of the protoplanetary disc, specifically the profiles of temperature, gas scale height and density. The commonly used disc structure of the Minimum Mass Solar Nebular (MMSN) is a simple power law in all these quantities. However, protoplanetary disc models with both viscous and stellar heating show several bumps and dips in temperature, scale height and density caused by transitions in opacity, which are missing in the MMSN model. These play an important role in the formation of planets, as they can act as sweet spots for the formation of planetesimals via the streaming instability and affect the direction and magnitude of type-I-migration. We present 2D simulations of accretion discs that feature radiative cooling, viscous and stellar heating, and are linked to the observed evolutionary stages of protoplanetary discs and their host stars. These models allow us to identify preferred planetesimal and planet formation regions in the protoplanetary disc as a function of the disc's metallicity, accretion rate and lifetime. We derive simple fitting formulae that feature all structural characteristics of protoplanetary discs during the evolution of several Myr. These fits are straightforward to apply for modelling any growth stage of planets where detailed knowledge of the underlying disc structure is required.

preprint2014arXiv

Water Delivery and Giant Impacts in the 'Grand Tack' Scenario

A new model for terrestrial planet formation (Hansen 2009, Walsh et al. 2011) has explored accretion in a truncated protoplanetary disk, and found that such a configuration is able to reproduce the distribution of mass among the planets in the Solar System, especially the Earth/Mars mass ratio, which earlier simulations have generally not been able to match. Walsh et al. tested a possible mechanism to truncate the disk--a two-stage, inward-then-outward migration of Jupiter and Saturn, as found in numerous hydrodynamical simulations of giant planet formation. In addition to truncating the disk and producing a more realistic Earth/Mars mass ratio, the migration of the giant planets also populates the asteroid belt with two distinct populations of bodies--the inner belt is filled by bodies originating inside of 3 AU, and the outer belt is filled with bodies originating from between and beyond the giant planets (which are hereafter referred to as `primitive' bodies). We find here that the planets will accrete on order 1-2% of their total mass from primitive planetesimals scattered onto planet-crossing orbits during the formation of the planets. For an assumed value of 10% for the water mass fraction of the primitive planetesimals, this model delivers a total amount of water comparable to that estimated to be on the Earth today. While the radial distribution of the planetary masses and the dynamical excitation of their orbits are a good match to the observed system, we find that the last giant impact is typically earlier than 20 Myr, and a substantial amount of mass is accreted after that event. However, 5 of the 27 planets larger than half an Earth mass formed in all simulations do experience large late impacts and subsequent accretion consistent with the dating of the Moon-forming impact and the estimated amount of mass accreted by Earth following that event.

preprint2013arXiv

Analytical Treatment of Planetary Resonances

An ever-growing observational aggregate of extrasolar planets has revealed that systems of planets that reside in or near mean-motion resonances are relatively common. While the origin of such systems is attributed to protoplanetary disk-driven migration, a qualitative description of the dynamical evolution of resonant planets remains largely elusive. Aided by the pioneering works of the last century, we formulate an approximate, integrable theory for first-order resonant motion. We utilize the developed theory to construct an intuitive, geometrical representation of resonances within the context of the unrestricted three-body problem. Moreover, we derive a simple analytical criterion for the appearance of secondary resonances between resonant and secular motion. Subsequently, we demonstrate the onset of rapid chaotic motion as a result of overlap among neighboring first-order mean-motion resonances, as well as the appearance of slow chaos as a result of secular modulation of the planetary orbits. Finally, we take advantage of the integrable theory to analytically show that, in the adiabatic regime, divergent encounters with first-order mean-motion resonances always lead to persistent apsidal anti-alignment.

preprint2013arXiv

Capture of Trojans by Jumping Jupiter

Jupiter Trojans are thought to be survivors of a much larger population of planetesimals that existed in the planetary region when planets formed. They can provide important constraints on the mass and properties of the planetesimal disk, and its dispersal during planet migration. Here we tested a possibility that the Trojans were captured during the early dynamical instability among the outer planets (aka the Nice model), when the semimajor axis of Jupiter was changing as a result of scattering encounters with an ice giant. The capture occurs in this model when Jupiter's orbit and its Lagrange points become radially displaced in a scattering event and fall into a region populated by planetesimals (that previously evolved from their natal transplanetary disk to ~5 AU during the instability). Our numerical simulations of the new capture model, hereafter jump capture, satisfactorily reproduce the orbital distribution of the Trojans and their total mass. The jump capture is potentially capable of explaining the observed asymmetry in the number of leading and trailing Trojans. We find that the capture probability is (6-8)x10^-7 for each particle in the original transplanetary disk, implying that the disk contained (3-4)x10^7 planetesimals with absolute magnitude H<9 (corresponding to diameter D=80 km for a 7% albedo). The disk mass inferred from this work, M_disk=14-28 M_Earth, is consistent with the mass deduced from recent dynamical simulations of the planetary instability.

preprint2013arXiv

Early Dynamical Instabilities in the Giant Planet Systems

The observed wide eccentricity distribution of extrasolar giant planets is thought to be the result of dynamical instabilities and gravitational scattering among planets. Previously, it has been assumed that the orbits in giant planet systems become gravitationally unstable after the gas nebula dispersal. It was not well understood, however, how these unstable conditions were established in the first place. In this work we numerically simulate the evolution of systems of three planets as the planets sequentially grow to Jupiter's mass, and dynamically interact among themselves and with the gas disk. We use the hydro-dynamical code FARGO that we modified by implementing the $N$-body integrator SyMBA. The new code can handle close encounters and collisions between planets. To test their stability, the planetary systems were followed with SyMBA for up to $10^8$ yr after the gas disk dispersal. We find that dynamics of the growing planets is complex, because migration and resonances raise their orbital eccentricities, and cause dynamical instabilities when gas is still around. If the dynamical instabilities occur early, planets can be removed by collisions and ejections, and the system rearranges into a new, more stable configuration. In this case, the planetary systems emerging from the gas disks are expected to be stable, and would need to be destabilized by other means (low-mass planets, planetesimal disks, etc.). Alternatively, for the giant planet system to be intrinsically unstable upon the gas disk dispersal, a special timing would be required with the growth of (at least some of) the giant planets having to occur near the end of the gas disk lifetime.

preprint2013arXiv

Geochemical and planetary dynamical views on the origin of Earth's atmosphere and oceans

Earth's volatile elements (H, C, and N) are essential to maintaining habitable conditions for metazoans and simpler life forms. However, identifying the sources (comets, meteorites, and trapped nebular gas) that supplied volatiles to Earth is not straightforward because secondary processes like mantle degassing, crustal recycling, and escape to space modified the composition of the atmosphere. Here, we review two complementary approaches to investigate the origin of Earth's atmosphere and oceans. The geochemical approach uses volatile element abundances and isotopic compositions to identify the possible contributors to the atmosphere and to disentangle the processes that shaped it. In that respect, noble gases (He, Ne, Ar, Kr, and Xe), elements that are chemically inert and possess several isotopes produced by radioactivity, play a critical role. The dynamical approach uses our knowledge of planetary dynamics to track volatile delivery to the Earth, starting with dust transport in the disk to planet-building processes. The main conclusion is that Earth acquired most of its major volatile elements by accretion of planetesimals or embryos akin to volatile-rich meteorites. At the same time, solar/meteoritic noble gases were captured by embryos and some gases were lost to space, by hydrodynamic escape and large impacts. Comets did not contribute much H, C, and N but may have delivered significant noble gases, which could represent the only fingerprints of the bombardment of our planet with icy bodies. The processes that governed the delivery of volatile elements to the Earth are thought to be relatively common and it is likely that Earth-like planets covered with oceans exist in extra-solar systems.

preprint2013arXiv

The Eos family halo

We study K-type asteroids in the broad surroundings of the Eos family because they seem to be intimately related, according to their colours measured by the Sloan Digital Sky Survey. Such `halos' of asteroid families have been rarely used as constraints for dynamical studies to date. We explain its origin as bodies escaping from the family `core' due to the Yarkovsky semimajor-axis drift and interactions with gravitational resonances, mostly with the 9/4 mean-motion resonance with Jupiter at 3.03 AU. Our N-body dynamical model allows us to independently estimate the age of the family 1.5 to 1.9 Gyr. This is approximately in agreement with the previous age estimate by Vokrouhlický et al. (2006) based on a simplified model (which accounts only for changes of semimajor axis). We can also constrain the geometry of the disruption event which had to occur at the true anomaly f =~ 120 to 180 deg.

preprint2013arXiv

The Long-Term Dynamical Evolution of Planetary Systems

This chapter concerns the long-term dynamical evolution of planetary systems from both theoretical and observational perspectives. We begin by discussing the planet-planet interactions that take place within our own Solar System. We then describe such interactions in more tightly-packed planetary systems. As planet-planet interactions build up, some systems become dynamically unstable, leading to strong encounters and ultimately either ejections or collisions of planets. After discussing the basic physical processes involved, we consider how these interactions apply to extrasolar planetary systems and explore the constraints provided by observed systems. The presence of a residual planetesimal disc can lead to planetary migration and hence cause instabilities induced by resonance crossing; however, such discs can also stabilise planetary systems. The crowded birth environment of a planetary system can have a significant impact: close encounters and binary companions can act to destabilise systems, or sculpt their properties. In the case of binaries, the Kozai mechanism can place planets on extremely eccentric orbits which may later circularise to produce hot Jupiters.

preprint2012arXiv

A Sawtooth-like Timeline for the First Billion Year of Lunar Bombardment

We revisit the early evolution of the Moon's bombardment. Our work combines modeling (based on plausible projectile sources and their dynamical decay rates) with constraints from the lunar crater record, radiometric ages of the youngest lunar basins, and the abundance of highly siderophile elements in the lunar crust and mantle. We deduce that the evolution of the impact flux did not decline exponentially over the first billion years of lunar history, but also there was no prominent and "narrow" impact spike some 3.9 Gy ago, unlike that typically envisioned in the lunar cataclysm scenario. Instead, we show the timeline of the lunar bombardment has a sawtooth-like profile, with an uptick in the impact flux near 4.1 Gy ago. The impact flux at the beginning of this weaker cataclysm was 5-10 times higher than the immediately preceding period. The Nectaris basin should have been one of the first basins formed at the sawtooth. We predict the bombardment rate since about 4.1Gy ago declined slowly and adhered relatively close to classic crater chronology models (Neukum and Ivanov (1994)). Overall we expect that the sawtooth event accounted for about 1/4 of the total bombardment suffered by the Moon since its formation. Consequently, considering that about 12-14 basins formed during the sawtooth event, we expect that the net number of basins formed on the Moon was about 45-50. From our expected bombardment timeline, we derived a new and improved lunar chronology suitable for use on Pre-Nectarian surface units. According to this chronology, a significant portion of the oldest lunar cratered terrains has an age of 4.38-4.42 Gyr. Moreover, the largest lunar basin, South Pole Aitken, is older than 4.3Gy, and therefore was not produced during the lunar cataclysm.

preprint2012arXiv

Building Terrestrial Planets

This paper reviews our current understanding of terrestrial planets formation. The focus is on computer simulations of the dynamical aspects of the accretion process. Throughout the chapter, we combine the results of these theoretical models with geochemical, cosmochemical and chronological constraints, in order to outline a comprehensive scenario of the early evolution of our Solar System. Given that the giant planets formed first in the protoplanetary disk, we stress the sensitive dependence of the terrestrial planet accretion process on the orbital architecture of the giant planets and on their evolution. This suggests a great diversity among the terrestrial planets populations in extrasolar systems. Issues such as the cause for the different masses and accretion timescales between Mars and the Earth and the origin of water (and other volatiles) on our planet are discussed at depth.

preprint2012arXiv

Circum-planetary discs as bottlenecks for gas accretion onto giant planets

With hundreds of exoplanets detected, it is necessary to revisit giant planets accretion models to explain their mass distribution. In particular, formation of sub-jovian planets remains unclear, given the short timescale for the runaway accretion of massive atmospheres. However, gas needs to pass through a circum-planetary disc. If the latter has a low viscosity (as expected if planets form in "dead zones"), it might act as a bottleneck for gas accretion. We investigate what the minimum accretion rate is for a planet under the limit assumption that the circum-planetary disc is totally inviscid, and the transport of angular momentum occurs solely because of the gravitational perturbations from the star. To estimate the accretion rate, we present a steady-state model of an inviscid circum-planetary disc, with vertical gas inflow and external torque from the star. Hydrodynamical simulations of a circum-planetary disc were conducted in 2D, in a planetocentric frame, with the star as an external perturber in order to measure the torque exerted by the star on the disc. The disc shows a two-armed spiral wave caused by stellar tides, propagating all the way in from the outer edge of the disc towards the planet. The stellar torque is small and corresponds to a doubling time for a Jupiter mass planet of the order of 5 Myrs. Given the limit assumptions, this is clearly a lower bound of the real accretion rate. This result shows that gas accretion onto a giant planet can be regulated by circum-planetary discs. This suggests that the diversity of masses of extra-solar planets may be the result of different viscosities in these discs.

preprint2012arXiv

Dissipative Divergence of Resonant Orbits

A considerable fraction of multi-planet systems discovered by the observational surveys of extrasolar planets reside in mild proximity to first-order mean motion resonances. However, the relative remoteness of such systems from nominal resonant period ratios (e.g. 2:1, 3:2, 4:3) has been interpreted as evidence for lack of resonant interactions. Here we show that a slow divergence away from exact commensurability is a natural outcome of dissipative evolution and demonstrate that libration of critical angles can be maintained tens of percent away from nominal resonance. We construct an analytical theory for the long-term dynamical evolution of dissipated resonant planetary pairs and confirm our calculations numerically. Collectively, our results suggest that a significant fraction of the near-commensurate extrasolar planets are in fact resonant and have undergone significant dissipative evolution.

preprint2012arXiv

Dynamics of pebbles in the vicinity of a growing planetary embryo: hydro-dynamical simulations

Understanding the growth of the cores of giant planets is a difficult problem. Recently, Lambrechts and Johansen (2012; LJ12) proposed a new model in which the cores grow by the accretion of pebble-size objects, as the latter drift towards the star due to gas drag. Here, we investigate the dynamics of pebble-size objects in the vicinity of planetary embryos of 1 and 5 Earth masses and the resulting accretion rates. We use hydrodynamical simulations, in which the embryo influences the dynamics of the gas and the pebbles suffer gas drag according to the local gas density and velocities. The pebble dynamics in the vicinity of the planetary embryo is non-trivial, and it changes significantly with the pebble size. Nevertheless, the accretion rate of the embryo that we measure is within an order of magnitude of the rate estimated in LJ12 and tends to their value with increasing pebble-size. We conclude that the model by LJ12 has the potential to explain the rapid growth of giant planet cores. The actual accretion rates however, depend on the surface density of pebble size objects in the disk, which is unknown to date.

preprint2012arXiv

NOTE: Explaining why the Uranian satellites have equatorial prograde orbits despite the large planetary obliquity

We show that the existence of prograde equatorial satellites is consistent with a collisional tilting scenario for Uranus. In fact, if the planet was surrounded by a proto-satellite disk at the time of the tilting and a massive ring of material was temporarily placed inside the Roche radius of the planet by the collision, the proto-satellite disk would have started to precess incoherently around the equator of the planet, up to a distance greater than that of Oberon. Collisional damping would then have collapsed it into a thin equatorial disk, from which the satellites eventually formed. The fact that the orbits of the satellites are prograde requires Uranus to have had a non-negligible initial obliquity (comparable to that of Neptune) before it was finally tilted to 98 degrees.

preprint2012arXiv

On the origin of the Almahata-Sitta meteorite and 2008TC3 asteroid

Asteroid 2008TC3 was a Near Earth Asteroid that impacted the Earth on 2008 October 7. Meteorites were produced by the break-up of 2008TC3 in the high atmosphere and at present, about 600 meteorites - called Almahata Sitta - coming from 2008TC3 have been recovered. A mineralogical study of Almahata Sitta fragments shows that the asteroid 2008TC3 was made of meteorites of different types (ureilites, H, L, and E chondrites). Understanding the origin of this body and how it was put together remain a challenge. Here we perform a detailed spectroscopical and dynamical investigation to show that the most likely source region of 2008TC3 is in the inner Main Belt at low inclination (i<8 degrees). We show that asteroids with spectroscopic classes that can be associated with the different meteorite types of Almahata Sitta are present in the region of the Main Belt that includes the Nysa-Polana family and objects of the Background at low inclination. Searching for a possible scenario of formation for 2008TC3, we show that there is little chance that 2008TC3 was formed by low velocity collisions between asteroids of different mineralogies, in the current asteroid belt. It seems more likely that the heterogeneous composition of 2008TC3 was a inherited from a time when the asteroid belt was in a different dynamical state, most likely in the very early Solar System. Because ureilites are fragments of a large, thermally metamorphosed asteroid, this suggests that the phases of collisional erosion (the break-up of the ureilite parent-body) and collisional accretion (the formation of the parent body of 2008TC3) overlapped for some time in the primordial asteroid belt.

preprint2012arXiv

Statistical Study of the Early Solar System's Instability with 4, 5 and 6 Giant Planets

Several properties of the Solar System, including the wide radial spacing and orbital eccentricities of giant planets, can be explained if the early Solar System evolved through a dynamical instability followed by migration of planets in the planetesimal disk. Here we report the results of a statistical study, in which we performed nearly 10^4 numerical simulations of planetary instability starting from hundreds of different initial conditions. We found that the dynamical evolution is typically too violent, if Jupiter and Saturn start in the 3:2 resonance, leading to ejection of at least one ice giant from the Solar System. Planet ejection can be avoided if the mass of the transplanetary disk of planetesimals was large (M_disk>50 M_Earth), but we found that a massive disk would lead to excessive dynamical damping (e.g., final e_55 < 0.01 compared to present e_55=0.044, where e_55 is the amplitude of the fifth eccentric mode in the Jupiter's orbit), and to smooth migration that violates constraints from the survival of the terrestrial planets. Better results were obtained when the Solar System was assumed to have five giant planets initially and one ice giant, with the mass comparable to that of Uranus and Neptune, was ejected into interstellar space by Jupiter. The best results were obtained when the ejected planet was placed into the external 3:2 or 4:3 resonance with Saturn and M_disk ~ 20 M_Earth. The range of possible outcomes is rather broad in this case, indicating that the present Solar System is neither a typical nor expected result for a given initial state, and occurs, in best cases, with only a ~5% probability (as defined by the success criteria described in the main text). The case with six giant planets shows interesting dynamics but does offer significant advantages relative to the five planet case.

preprint2012arXiv

Stellar irradiated discs and implications on migration of embedded planets I: equilibrium discs

The strength and direction of migration of low mass planets depends on the disc's thermodynamics. In discs where the viscous heating is balanced by radiative transport, the migration can be directed outwards, a process which extends the lifetime of growing planetary embryos. We investigate the influence of opacity and stellar irradiation on the disc thermodynamics. Utilizing the resulting disc structure, we determine the regions of outward migration. We perform two-dimensional numerical simulations of equilibrium discs with viscous heating, radiative cooling and stellar irradiation. We use the hydrodynamical code NIRVANA that includes a full tensor viscosity and stellar irradiation, as well as a two temperature solver that includes radiation transport in the flux-limited diffusion approximation. The migration is studied by using torque formulae. In the constant opacity case, we reproduce the analytical results of a black-body disc: the stellar irradiation dominates in the outer regions -- leading to flaring -- while the viscous heating dominates close to the star. We find that the inner edge of the disc should not be significantly puffed-up by the stellar irradiation. If the opacity depends on the local density and temperature, the structure of the disc is different, and several bumps in the aspect ratio H/r appear, due to transitions between different opacity regimes. The bumps in the disc can shield the outer disc from stellar irradiation. Stellar irradiation is an important factor for determining the disc structure and has dramatic consequences for the migration of embedded planets. Compared to discs with only viscous heating, a stellar irradiated disc features a much smaller region of outward migration for a range of planetary masses. This suggests that the region where the formation of giant planet cores takes place is smaller, which in turn might lead to a shorter growth phase.

preprint2011arXiv

Dynamical evolution of planetary systems

The apparent regularity of the motion of the giant planets of our solar system suggested for decades that said planets formed onto orbits similar to the current ones and that nothing dramatic ever happened during their lifetime. The discovery of extra-solar planets showed astonishingly that the orbital structure of our planetary system is not typical. Many giant extra-solar planets have orbits with semi major axes of $\sim 1$ AU, and some have even smaller orbital radii, sometimes with orbital periods of just a few days. Moreover, most extra-solar planets have large eccentricities, up to values that only comets have in our solar system. Why such a big diversity between our solar system and the extra-solar systems, as well as among the extra-solar systems themselves? This chapter aims to give a partial answer to this fundamental question....

preprint2011arXiv

Formation and Evolution of Planetary Systems in Presence of Highly Inclined Stellar Perturbers

The presence of highly eccentric extrasolar planets in binary stellar systems suggests that the Kozai effect has played an important role in shaping their dynamical architectures. However, the formation of planets in inclined binary systems poses a considerable theoretical challenge, as orbital excitation due to the Kozai resonance implies destructive, high-velocity collisions among planetesimals. To resolve the apparent difficulties posed by Kozai resonance, we seek to identify the primary physical processes responsible for inhibiting the action of Kozai cycles in protoplanetary disks. Subsequently, we seek to understand how newly-formed planetary systems transition to their observed, Kozai-dominated dynamical states. We find that theoretical difficulties in planet formation arising from the presence of a distant companion star, posed by the Kozai effect and other secular perturbations, can be overcome by a proper account of gravitational interactions within the protoplanetary disk. In particular, fast apsidal recession induced by disk self-gravity tends to erase the Kozai effect, and ensure that the disk's unwarped, rigid structure is maintained. Subsequently, once a planetary system has formed, the Kozai effect can continue to be wiped out as a result of apsidal precession, arising from planet-planet interactions. However, if such a system undergoes a dynamical instability, its architecture may change in such a way that the Kozai effect becomes operative. The results presented here suggest that planetary formation in highly inclined binary systems is not stalled by perturbations, arising from the stellar companion. Consequently, planet formation in binary stars is probably no different from that around single stars on a qualitative level. Furthermore, it is likely that systems where the Kozai effect operates, underwent a transient phase of dynamical instability in the past.

preprint2011arXiv

Onset of Secular Chaos in Planetary Systems: Period Doubling & Strange Attractors

As a result of resonance overlap, planetary systems can exhibit chaotic motion. Planetary chaos has been studied extensively in the Hamiltonian framework, however, the presence of chaotic motion in systems where dissipative effects are important, has not been thoroughly investigated. Here, we study the onset of stochastic motion in presence of dissipation, in the context of classical perturbation theory, and show that planetary systems approach chaos via a period-doubling route as dissipation is gradually reduced. Furthermore, we demonstrate that chaotic strange attractors can exist in mildly damped systems. The results presented here are of interest for understanding the early dynamical evolution of chaotic planetary systems, as they may have transitioned to chaos from a quasi-periodic state, dominated by dissipative interactions with the birth nebula.

preprint2011arXiv

The Effect of an Early Planetesimal-Driven Migration of the Giant Planets on Terrestrial Planet Formation

The migration of the giant planets due to the scattering of planetesimals causes powerful resonances to move through the asteroid belt and the terrestrial planet region. Exactly when and how the giant planets migrated is not well known. In this paper we present results of an investigation of the formation of the terrestrial planets during and after the migration of the giant planets. The latter is assumed to have occurred immediately after the dissipation of the nebular disk -- i.e. "early" with respect to the timing of the Late Heavy Bombardment (LHB). The presumed cause of our modeled early migration of the giant planets is angular mometum transfer between the planets and scattered planetesimals.

preprint2010arXiv

Calculation of the enrichment of the giant planet envelopes during the "late heavy bombardment"

The giant planets of our solar system possess envelopes consisting mainly of hydrogen and helium but are also significantly enriched in heavier elements relatively to our Sun. In order to better constrain how these heavy elements have been delivered, we quantify the amount accreted during the so-called "late heavy bombardment", at a time when planets were fully formed and planetesimals could not sink deep into the planets. On the basis of the "Nice model", we obtain accreted masses (in terrestrial units) equal to $0.15\pm0.04 \rm\,M_\oplus$ for Jupiter, and $0.08 \pm 0.01 \rm\,M_\oplus$ for Saturn. For the two other giant planets, the results are found to depend mostly on whether they switched position during the instability phase. For Uranus, the accreted mass is $0.051 \pm 0.003 \rm\,M_\oplus$ with an inversion and $0.030 \pm 0.001 \rm\,M_\oplus$ without an inversion. Neptune accretes $0.048 \pm 0.015 \rm\,M_\oplus$ in models in which it is initially closer to the Sun than Uranus, and $0.066 \pm 0.006 \rm\,M_\oplus$ otherwise. With well-mixed envelopes, this corresponds to an increase in the enrichment over the solar value of $0.033 \pm 0.001$ and $0.074 \pm 0.007$ for Jupiter and Saturn, respectively. For the two other planets, we find the enrichments to be $2.1 \pm 1.4$ (w/ inversion) or $1.2 \pm 0.7$ (w/o inversion) for Uranus, and $2.0 \pm 1.2$ (w/ inversion) or $2.7 \pm 1.6$ (w/o inversion) for Neptune. This is clearly insufficient to explain the inferred enrichments of $\sim 4$ for Jupiter, $\sim 7$ for Saturn and $\sim 45$ for Uranus and Neptune.

preprint2010arXiv

Evidence from the asteroid belt for a violent past evolution of Jupiter's orbit

We use the current orbital structure of large (>50km) asteroids in the main asteroid belt to constrain the evolution of the giant planets when they migrated from their primordial orbits to their current ones. Minton & Malhotra (2009) showed that the orbital distribution of large asteroids in the main belt can be reproduced by an exponentially-decaying migration of the giant planets on a time scale of tau ~ 0.5My. However, self-consistent numerical simulations show that the planetesimal-driven migration of the giant planets is inconsistent with an exponential change in their semi major axes on such a short time scale (Hahn & Malhotra, 1999). In fact, the typical time scale is tau > 5My. When giant planet migration on this time scale is applied to the asteroid belt, the resulting orbital distribution is incompatible with the observed one. However, the planet migration can be significantly sped up by planet-planet encounters. Consider an evolution where both Jupiter and Saturn have close encounters with a Neptune-mass planet (presumably Uranus or Neptune themselves) and where this third planet, after being scattered inwards by Saturn, is scattered outwards by Jupiter. This scenario leads to a very rapid increase in the orbital separation between Jupiter and Saturn that we show here to have only mild effects on the structure of asteroid belt. This type of evolution is called a jumping-Jupiter case. Our results suggest that the total mass and dynamical excitation of the asteroid belt before migration were comparable to those currently observed. Moreover, they imply that, before migration, the orbits of Jupiter and Saturn were much less eccentric than the current ones.

preprint2010arXiv

Review: A Coherent and Comprehensive Model of the Evolution of the Outer Solar System

Since the discovery of the first extra-solar planets, we are confronted with the puzzling diversity of planetary systems. Processes like planet radial migration in gas-disks and planetary orbital instabilities, often invoked to explain the exotic orbits of the extra-solar planets, at first sight do not seem to have played a role in our system. In reality, though, there are several aspects in the structure of our Solar System that cannot be explained in the classic scenario of in-situ formation and smooth evolution of the giant planets. This paper describes a new view of the evolution of the outer Solar System that emerges from the so-called 'Nice model' and its recent extensions. The story provided by this model describes a very "dynamical" Solar System, with giant planets affected by both radial migrations and a temporary orbital instability. Thus, the diversity between our system and those found so far around other stars does not seem to be due to different processes that operated here and elsewhere, but rather stems from the strong sensitivity of chaotic evolutions to small differences in the initial and environmental conditions.

preprint2009arXiv

Asteroids Were Born Big

How big were the first planetesimals? We attempt to answer this question by conducting coagulation simulations in which the planetesimals grow by mutual collisions and form larger bodies and planetary embryos. The size frequency distribution (SFD) of the initial planetesimals is considered a free parameter in these simulations, and we search for the one that produces at the end objects with a SFD that is consistent with asteroid belt constraints. We find that, if the initial planetesimals were small (e.g. km-sized), the final SFD fails to fulfill these constraints. In particular, reproducing the bump observed at diameter D~100km in the current SFD of the asteroids requires that the minimal size of the initial planetesimals was also ~100km. This supports the idea that planetesimals formed big, namely that the size of solids in the proto-planetary disk ``jumped'' from sub-meter scale to multi-kilometer scale, without passing through intermediate values. Moreover, we find evidence that the initial planetesimals had to have sizes ranging from 100 to several 100km, probably even 1,000km, and that their SFD had to have a slope over this interval that was similar to the one characterizing the current asteroids in the same size-range. This result sets a new constraint on planetesimal formation models and opens new perspectives for the investigation of the collisional evolution in the asteroid and Kuiper belts as well as of the accretion of the cores of the giant planets.

preprint2009arXiv

Building the Terrestrial Planets: Constrained Accretion in the Inner Solar System

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. (abridged)

preprint2009arXiv

Constructing the secular architecture of the solar system II: The terrestrial planets

We investigate the dynamical evolution of the terrestrial planets during the planetesimal-driven migration of the giant planets. A basic assumption of this work is that giant planet migration occurred after the completion of terrestrial planet formation, such as in the models that link the former to the origin of the Late Heavy Bombardment. The divergent migration of Jupiter and Saturn causes the g5 eigenfrequency to cross resonances of the form g5=gk with k ranging from 1 to 4. Consequently these secular resonances cause large-amplitude responses in the eccentricities of the terrestrial planets. We show that the resonances g5=g_4 and g5=g3 do not pose a problem if Jupiter and Saturn have a fast approach and departure from their mutual 2:1 mean motion resonance. On the other hand, the resonance crossings g5=g2 and g5=g1 are more of a concern as they tend to yield a terrestrial system incompatible with the current one. We offer two solutions to this problem. The first uses the fact that a secular resonance crossing can also damp the amplitude of a Fourier mode if the latter is large originally. A second scenario involves a 'jumping Jupiter' in which encounters between an ice giant and Jupiter, without ejection of the former, cause the latter to migrate away from Saturn much faster than if migration is driven solely by encounters with planetesimals. In this case, the g5=g2 and g5=g1 resonances can be jumped over, or occur very briefly.

preprint2007arXiv

Planetary embryos and planetesimals residing in thin debris disks

We consider constraints on the planetesimal population residing in the disks of AU Microscopii, Beta Pictoris and Fomalhaut taking into account their observed thicknesses and normal disk opacities. We estimate that bodies of radius 5, 180 and 70 km are responsible for initiating the collisional cascade accounting for the dust production for AU-Mic, Beta-Pic and Fomalhaut's disks, respectively, at break radii from the star where their surface brightness profiles change slope. Larger bodies, of radius 1000km and with surface density of order 0.01 g/cm^2, are required to explain the thickness of these disks assuming that they are heated by gravitational stirring. A comparison between the densities of the two sizes suggests the size distribution in the largest bodies is flatter than that observed in the Kuiper belt. AU Mic's disk requires the shallowest size distribution for bodies with radius greater than 10km suggesting that the disk contains planetary embryos experiencing a stage of runaway growth.