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.

In two recent papers published in MNRAS, Namouni and Morais (2018, 2020) claimed evidence for the interstellar origin of some small Solar System bodies, including i) objects in retrograde co-orbital motion with the giant planets, and ii) the highly-inclined Centaurs. Here, we discuss the flaws of those papers that invalidate the authors' conclusions. Numerical simulations backwards in time are not representative of the past evolution of real bodies. Instead, these simulations are only useful as a means to quantify the short dynamical lifetime of the considered bodies and the fast decay of their population. In light of this fast decay, if the observed bodies were the survivors of populations of objects captured from interstellar space in the early Solar System, these populations should have been implausibly large (e.g. about 10 times the current main asteroid belt population for the retrograde coorbital of Jupiter). More likely, the observed objects are just transient members of a population that is maintained in quasi-steady state by a continuous flux of objects from some parent reservoir in the distant Solar System. We identify in the Halley type comets and the Oort cloud the most likely sources of retrograde coorbitals and highly-inclined Centaurs.

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.

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.

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 \pi (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$.

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.

The Al/Si and Mg/Si ratios in non-carbonaceous chondrites are lower than the solar (i.e., CI-chondritic) values, in sharp contrast to the non-CI carbonaceous meteorites and the Earth, which are enriched in refractory elements and have Mg/Si ratios that are solar or larger. We show that the formation of a first generation of planetesimals during the condensation of refractory elements implies the subsequent formation of residual condensates with strongly sub-solar Al/Si and Mg/Si ratios. The mixing of residual condensates with different amounts of material with solar refractory/Si element ratios explains the Al/Si and Mg/Si values of non-carbonaceous chondrites. To match quantitatively the observed ratios, we find that the first-planetesimals should have accreted when the disk temperature was ~1,330-1,400 K depending on pressure and assuming a solar C/O ratio of the disk. We discuss how this model relates to our current understanding of disk evolution, grain dynamics, and planetesimal formation. We also extend the discussion to moderately volatile elements (e.g., Na), explaining how it may be possible that the depletion of these elements in non-carbonaceous chondrites is correlated with the depletion of refractory elements (e.g., Al). Extending the analysis to Cr, we find evidence for a higher than solar C/O ratio in the protosolar disk's gas when/where condensation from a fractionated gas occurred. Finally, we discuss the possibility that the supra-solar Al/Si and Mg/Si ratios of the Earth are due to the accretion of ~40% of the mass of our planet from the first-generation of refractory-rich planetesimals.

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.

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.

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.

The dynamical structure of the Solar System can be explained by a period of orbital instability experienced by the giant planets. While a late instability was originally proposed to explain the Late Heavy Bombardment, recent work favors an early instability. We model the early dynamical evolution of the outer Solar System to self-consistently constrain the most likely timing of the instability. We first simulate the dynamical sculpting of the primordial outer planetesimal disk during the accretion of Uranus and Neptune from migrating planetary embryos during the gas disk phase, and determine the separation between Neptune and the inner edge of the planetesimal disk. We performed simulations with a range of migration histories for Jupiter. We find that, unless Jupiter migrated inwards by 10 AU or more, the instability almost certainly happened within 100 Myr of the start of Solar System formation. There are two distinct possible instability triggers. The first is an instability that is triggered by the planets themselves, with no appreciable influence from the planetesimal disk. Of those, the median instability time is $\sim4$Myr. Among self-stable systems -- where the planets are locked in a resonant chain that remains stable in the absence of a planetesimal's disk-- our self-consistently sculpted planetesimal disks nonetheless trigger a giant planet instability with a median instability time of 37-62 Myr for a reasonable range of migration histories of Jupiter. The simulations that give the latest instability times are those that invoked long-range inward migration of Jupiter from 15 AU or beyond; however these simulations over-excited the inclinations of Kuiper belt objects and are inconsistent with the present-day Solar System. We conclude on dynamical grounds that the giant planet instability is likely to have occurred early in Solar System history.

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.

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.

This paper reviews habitability conditions for a terrestrial planet from the point of view of geosciences. It addresses how interactions between the interior of a planet or a moon and its atmosphere and surface (including hydrosphere and biosphere) can affect habitability of the celestial body. It does not consider in detail the role of the central star but focusses more on surface conditions capable of sustaining life. We deal with fundamental issues of planetary habitability, i.e. the environmental conditions capable of sustaining life, and the above-mentioned interactions can affect the habitability of the celestial body. We address some hotly debated questions including: - How do core and mantle affect the evolution and habitability of planets? - What are the consequences of mantle overturn on the evolution of the interior and atmosphere? - What is the role of the global carbon and water cycles? - What influence do comet and asteroid impacts exert on the evolution of the planet? - How does life interact with the evolution of the Earth's geosphere and atmosphere? - How can knowledge of the solar system geophysics and habitability be applied to exoplanets? In addition, we address the identification of preserved life tracers in the context of the interaction of life with planetary evolution.

Gas-giant planets, like Jupiter and Saturn, acquire massive gaseous envelopes during the approximately 3 Myr-long lifetimes of protoplanetary discs. In the core accretion scenario, the formation of a solid core of around 10 Earth masses triggers a phase of rapid gas accretion. Previous 3D grid-based hydrodynamical simulations found runaway gas accretion rates corresponding to approximately 10 to 100 Jupiter masses per Myr. Such high accretion rates would result in all planets with larger-than-10-Earth-mass cores forming Jupiter-like planets, in clear contrast to the ice giants in the Solar System and the observed exoplanet population. In this work, we use 3D hydrodynamical simulations, that include radiative transfer, to model the growth of the envelope on planets with different masses. We find that gas flows rapidly through the outer part of the envelope, but this flow does not drive accretion. Instead, gas accretion is the result of quasi-static contraction of the inner envelope, which can be orders of magnitude smaller than the mass flow through the outer atmosphere. For planets smaller than Saturn, we measure moderate gas accretion rates that are below 1 Jupiter mass per Myr. Higher mass planets, however, accrete up to 10 times faster and do not reveal a self-driven mechanism that can halt gas accretion. Therefore, the reason for the final masses of Saturn and Jupiter remains difficult to understand, unless their completion coincided with the dissipation of the Solar Nebula.

The evolution of the meteorite flux to the Earth can be studied by determining the terrestrial ages of meteorite collected in hot deserts. We have measured the terrestrial ages of 54 stony meteorites from the El M\'edano area, in the Atacama Desert, using the cosmogenic nuclide chlorine 36. With an average age of 710 ka, this collection is the oldest collection of non fossil meteorites at the Earth's surface. This allows both determining the average meteorite flux intensity over the last 2 Myr (222 meteorites larger than 10 g per km2 per Myr) and discussing its possible compositional variability over the Quaternary period. A change in the flux composition, with more abundant H chondrites, occurred between 0.5 and 1 Ma, possibly due to the direct delivery to Earth of a meteoroid swarm from the asteroid belt.

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.

Early dynamical evolution of close-in planetary systems is shaped by an intricate combination of planetary gravitational interactions, orbital migration, and dissipative effects. While the process of convergent orbital migration is expected to routinely yield resonant planetary systems, previous analyses have shown that the semi-major axes of initially resonant pairs of planets will gradually diverge under the influence of long-term energy damping, producing an overabundance of planetary period ratios in slight excess of exact commensurability. While this feature is clearly evident in the orbital distribution of close-in extrasolar planets, the existing theoretical picture is limited to the specific case of the planetary three-body problem. In this study, we generalise the framework of dissipative divergence of resonant orbits to multi-resonant chains, and apply our results to the current observational census of well-characterised three-planet systems. Focusing on the 2:1 and 3:2 commensurabilities, we identify three 3-planet systems, whose current orbital architecture is consistent with an evolutionary history wherein convergent migration first locks the planets into a multi-resonant configuration and subsequent dissipation repels the orbits away from exact commensurability. Nevertheless, we find that the architecture of the overall sample of multi planetary systems is incompatible with this simple scenario, suggesting that additional physical mechanisms must play a dominant role during the early stages of planetary systems' dynamical evolution.

Giant planets migrate though the protoplanetary disc as they grow. We investigate how the formation of planetary systems depends on the radial flux of pebbles through the protoplanetary disc and on the planet migration rate. Our N-body simulations confirm previous findings that Jupiter-like planets in orbits outside the water ice line originate from embryos starting out at 20-40 AU when using nominal type-I and type-II migration rates and a pebble flux of 100-200 Earth masses per million years, enough to grow Jupiter within the lifetime of the solar nebula. The planetary embryos placed up to 30AU migrate into the inner system (r<1AU) and form super-Earths or hot and warm gas giants, producing systems that are inconsistent with the configuration of the solar system, but consistent with some exoplanetary systems. We also explore slower migration rates which allow the formation of gas giants from embryos originating from the 5-10AU region, which are stranded exterior to 1 AU at the end of the gas-disc phase. We identify a pebble flux threshold below which migration dominates and moves the planetary core to the inner disc, where the pebble isolation mass is too low for the planet to accrete gas efficiently. Giant planet growth requires a sufficiently-high pebble flux to enable growth to out-compete migration. Even higher pebble fluxes produce systems with multiple gas giants. We show that planetary embryos starting interior to 5AU do not grow into gas giants, even if migration is slow and the pebble flux is large. Instead they grow to the mass regime of super-Earths. This stunted growth is caused by the low pebble isolation mass in the inner disc and is independent of the pebble flux. Additionally we show that the long term evolution of our formed planetary systems can produce systems with hot super-Earths and outer gas giants as well as systems of giants on eccentric orbits (abridged).

At least 30% of main sequence stars host planets with sizes between 1 and 4 Earth radii and orbital periods of less than 100 days. We use N-body simulations including a model for gas-assisted pebble accretion and disk-planet tidal interaction to study the formation of super-Earth systems. We show that the integrated pebble mass reservoir creates a bifurcation between hot super-Earths or hot-Neptunes ($\lesssim15M_{\oplus}$) and super-massive planetary cores potentially able to become gas giant planets ($\gtrsim15M_{\oplus}$). Simulations with moderate pebble fluxes grow multiple super-Earth-mass planets that migrate inwards and pile up at the disk's inner edge forming long resonant chains. We follow the long-term dynamical evolution of these systems and use the period ratio distribution of observed planet-pairs to constrain our model. Up to $\sim$95% of resonant chains become dynamically unstable after the gas disk dispersal, leading to a phase of late collisions that breaks the resonant configuration. Our simulations match observations if we combine a dominant fraction ($\gtrsim95\%$) of unstable systems with a sprinkling ($\lesssim5\%$) of stable resonant chains (the Trappist-1 system represents one such example). Our results demonstrate that super-Earth systems are inherently multiple (${\rm N\geq2}$) and that the observed excess of single-planet transits is a consequence of the mutual inclinations excited by the planet-planet instability. In simulations in which planetary seeds are initially distributed in the inner and outer disk, close-in super-Earths are systematically ice-rich. This contrasts with the interpretation that most super-Earths are rocky based on bulk density measurements of super-Earths and photo-evaporation modeling of their bimodal radius distribution. We investigate the conditions needed to form rocky super-Earths. The formation of rocky super-Earths (abridged)

Super-Earths are found in tighter orbits than the Earth's around more than one third of main sequence stars. It has been proposed that super-Earths are scaled-up terrestrial planets that formed similarly, through mutual accretion of planetary embryos, but in discs much denser than the solar protoplanetary disc. We argue instead that terrestrial planets and super-Earths have two distinct formation pathways that are regulated by the disc's pebble reservoir. Through numerical integrations, which combine pebble accretion and N-body gravity between embryos, we show that a difference of a factor of two in the pebble mass-flux is enough to change the evolution from the terrestrial to the super-Earth growth mode. If the pebble mass-flux is small, then the initial embryos within the ice line grow slowly and do not migrate substantially, resulting in a widely spaced population of Mars-mass embryos when the gas disc dissipates. Without gas being present, the embryos become unstable and a small number of terrestrial planets are formed by mutual collisions. The final terrestrial planets are at most 5 Earth masses. Instead, if the pebble mass-flux is high, then the initial embryos within the ice line rapidly become sufficiently massive to migrate through the gas disc. Embryos concentrate at the inner edge of the disc and growth accelerates through mutual merging. This leads to the formation of a system of closely spaced super-Earths in the 5 to 20 Earth-mass range, bounded by the pebble isolation mass. Generally, instabilities of these super-Earth systems after the disappearance of the gas disc trigger additional merging events and dislodge the system from resonant chains. The pebble flux - which controls the transition between the two growth modes - may be regulated by the initial reservoir of solids in the disc or the presence of more distant giant planets that can halt the radial flow of pebbles.

The near-Earth asteroid population suggests the existence of an inner Main Belt source of asteroids that belongs to the spectroscopic X-complex and has moderate albedos. The identification of such a source has been lacking so far. We argue that the most probable source is one or more collisional asteroid families that escaped discovery up to now. We apply a novel method to search for asteroid families in the inner Main Belt population of asteroids belonging to the X-complex with moderate albedo. Instead of searching for asteroid clusters in orbital elements space, which could be severely dispersed when older than some billions of years, our method looks for correlations between the orbital semimajor axis and the inverse size of asteroids. This correlation is the signature of members of collisional families, which drifted from a common centre under the effect of the Yarkovsky thermal effect. We identify two previously unknown families in the inner Main Belt among the moderate-albedo X-complex asteroids. One of them, whose lowest numbered asteroid is (161) Author, is ~3 Gyrs-old, whereas the second one, whose lowest numbered object is (689) Zita, can be as old as the Solar System. Members of this latter family have orbital eccentricities and inclinations that spread them over the entire inner Main Belt, which is an indication that this family could be primordial, i.e. it formed before the giant planet orbital instability. The vast majority of moderate-albedo X-complex asteroids of the inner-Main Belt are genetically related, as they can be included into few asteroid families. Only nine X-complex asteroids with moderate albedo of the inner Main Belt cannot be included in asteroid families. We suggest that these bodies formed by direct accretion of the solids in the protoplanetary disk, and are thus surviving planetesimals.

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.

In order to fully understand the shapes of asteroids families in the 3-dimensional space of the proper elements $(a_{\rm p}, e_{\rm p}, \sin I_{\rm p})$ it is necessary to compare observed asteroids with N-body simulations. To this point, we describe a rigorous yet simple method which allows for a selection of the observed asteroids, assures the same size-frequency distribution of synthetic asteroids, accounts for a background population, and computes a $\chi^2$ metric. We study the Eos family as an example, and we are able to fully explain its non-isotropic features, including the distribution of pole latitudes $\beta$. We confirm its age $t = (1.3\pm0.3)\,{\rm Gyr}$; while this value still scales with the bulk density, it is verified by a Monte-Carlo collisional model. The method can be applied to other populous families (Flora, Eunomia, Hygiea , Koronis, Themis, Vesta, etc.).

When and how planets form in protoplanetary disks is still a topic of discussion. Exoplanet detection surveys and protoplanetary disk surveys are now providing results that allow us to have new insights. We collect the masses of confirmed exoplanets and compare their dependence with stellar mass with the same dependence for protoplanetary disk masses measured in ~1-3 Myr old star-forming regions. The latter are recalculated by us using the new estimates of their distances derived from Gaia DR2 parallaxes. We note that single and multiple exoplanetary systems form two different populations, probably pointing to a different formation mechanism for massive giant planets around very low mass stars. While expecting that the mass in exoplanetary systems is much lower than the measured disk masses, we instead find that exoplanetary systems masses are comparable or higher than the most massive disks. This same result is found also by converting the measured planet masses into heavy-element content (core masses for the giant planets and full masses for the super-Earth systems) and by comparing this value with the disk dust masses. Unless disk dust masses are heavily underestimated, this is a big conundrum. An extremely efficient recycling of dust particles in the disk cannot solve this conundrum. This implies that either the cores of planets have formed very rapidly (<0.1-1 Myr) and large amount of gas is expelled on the same timescales from the disk, or that disks are continuously replenished of fresh planet-forming material from the environment. These hypotheses can be tested by measuring disk masses in even younger targets and by better understanding if and how the disks are replenished by their surroundings.

Massive planets form within the lifetime of protoplanetary disks and undergo orbital migration due to planet-disk interactions. When the first planet reaches the inner edge of the disk its migration stops and the second planet is locked in resonance. We detail the resonant trapping comparing semi-analytical formulae and numerical simulations in the case of two equal-mass coplanar planets trapped in first order resonances. We describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using expanded and non-expanded Hamiltonians. We show that during convergent migration the planets evolve along these families of equilibria. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that the resonant configurations achieved by migration can become unstable after the removal of the disk. We probe the stability of the resonances as a function of planetary mass. To do this, we fictitiously increase the masses of the planets, adiabatically maintaining the low-amplitude libration regime until instability. We discuss two hypotheses for instability: low-order secondary resonance of the libration frequency with a fast synodic frequency, and minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low-amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter still play the decisive role in the destabilisation of resonant pairs. We show that, as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.

The link between the dynamical evolution of the giant planets and the Kuiper Belt orbital structure can provide clues and insight about the dynamical history of the Solar System. The classical region of the Kuiper Belt has two populations (the cold and hot populations) with completely different physical and dynamical properties. These properties have been explained in the framework of a subset of the simulations of the Nice Model, in which Neptune remained on a low-eccentricity orbit (Neptune's eccentricity is never larger than 0.1) throughout the giant planet instability. However, recent simulations have showed that the remaining Nice model simulations, in which Neptune temporarily acquires a large-eccentricity orbit (larger than 0.1), are also consistent with the preservation of the cold population (inclination smaller than 4 degrees), if the latter formed in situ. However, the resulting a cold population showed in many of the simulations eccentricities larger than those observed for the real population. We focus on a short period of time which is characterized by Neptune's large eccentricity and a slow precession of Neptune's perihelion. We show that if self-gravity is considered in the disk, the precession rate of the particles longitude of perihelion is slowed down. This, combined with the effect of mutual scattering among the bodies, which spreads all orbital elements, allows some objects to return to low eccentricities. However, we show that if the cold population originally had a small total mass, this effect is negligible. Thus, we conclude that the only possibilities to keep at low eccentricity some cold-population objects during a high-eccentricity phase of Neptune are that (i) either Neptune's precession was rapid, as suggested by Batygin et al. (2011) or (ii) Neptune's slow precession phase was long enough to allow some particles to experience a full secular cycle.

Context. Giant planets open gaps in their protoplanetary and subsequently suffer so-called type II migration. Schematically, planets are thought to be tightly locked within their surrounding disks, and forced to follow the viscous advection of gas onto the central star. This fundamental principle however has recently been questioned, as migrating planets were shown to decouple from the gas' radial drift. Aims. In this framework, we question whether the traditionally used linear scaling of migration rate of a giant planet with the disk's viscosity still holds. Additionally, we assess the role of orbit-crossing material as part of the decoupling mechanism. Methods. We have performed 2D (r, {\theta}) numerical simulations of point-mass planets embedded in locally isothermal {\alpha}-disks in steady-state accretion, with various values of {\alpha}. Arbitrary planetary accretion rates were used as a means to diminish or nullify orbit-crossing flows. Results. We confirm that the migration rate of a gap-opening planet is indeed proportional to the disk's viscosity, but is not equal to the gas drift speed in the unperturbed disk. We show that the role of gap-crossing flows is in fact negligible. Conclusions. From these observations, we propose a new paradigm for type II migration : a giant planet feels a torque from the disk that promotes its migration, while the gap profile relative to the planet is restored on a viscous timescale, thus limiting the planet migration rate to be proportional to the disk's viscosity. Hence, in disks with low viscosity in the planet region, type II migration should still be very slow. Key words. protoplanetary disks; planet-disk interactions; planets and satellites: formation

The main asteroid belt (MB) is low in mass but dynamically excited. Here we propose a new mechanism to excite the MB during the giant planet ('Nice model') instability, which is expected to have featured repeated close encounters between Jupiter and one or more ice giants ('Jumping Jupiter' -- JJ). We show that, when Jupiter temporarily reaches a high enough level of excitation, both in eccentricity and inclination it induces strong forced vectors of eccentricity and inclination across the MB region. Because during the JJ instability Jupiter's orbit `jumps' around, the forced vectors keep changing both in magnitude and phase throughout the whole MB region. The entire cold primordial MB is thus excited as a natural outcome of the JJ instability. The level of such an excitation, however, is typically larger than the current orbital excitation observed in the MB. We show that the subsequent evolution of the Solar System is capable of reshaping the resultant over-excited MB to its present day orbital state, and that a strong mass depletion ($\sim$90$\%$) is associated to the JJ instability phase and its subsequent evolution throughout the age of the Solar System.

A number of multiplanet systems are observed to contain planets very close to mean motion resonances, although there is no significant pileup of precise resonance pairs. We present theoretical and numerical studies on the outcome of capture into first-order mean motion resonances (MMRs) using a parametrized planet migration model that takes into account nonlinear eccentricity damping due to planet-disk interaction. This parametrization is based on numerical hydrodynamical simulations and is more realistic than the simple linear parametrization widely used in previous analytic studies. We find that nonlinear eccentricity damping can significantly influence the stability and outcome of resonance capture. In particular, the equilibrium eccentricity of the planet captured into MMRs become larger, and the captured MMR state tends to be more stable compared to the prediction based on the simple migration model. In addition, when the migration is sufficiently fast or/and the planet mass ratio is sufficiently small, we observe a novel phenomenon of eccentricity overshoot, where the planet's eccentricity becomes very large before settling down to the lower equilibrium value. This can lead to the ejection of the smaller planet if its eccentricity approaches unity during the overshoot. This may help explain the lack of low-mass planet companion of hot Jupiters when compared to warm Jupiters.

The debiased absolute-magnitude and orbit distributions as well as source regions for near-Earth objects (NEOs) provide a fundamental frame of reference for studies of individual NEOs and more complex population-level questions. We present a new four-dimensional model of the NEO population that describes debiased steady-state distributions of semimajor axis, eccentricity, inclination, and absolute magnitude $H$ in the range $17<H<25$. The modeling approach improves upon the methodology originally developed by Bottke et al. (2000; Science 288, 2190-2194) in that it is, for example, based on more realistic orbit distributions and uses source-specific absolute-magnitude distributions that allow for a power-law slope that varies with $H$. We divide the main asteroid belt into six different entrance routes or regions (ER) to the NEO region: the $\nu_6$, 3:1J, 5:2J and 2:1J resonance complexes as well as Hungarias and Phocaeas. In addition we include the Jupiter-family comets as the primary cometary source of NEOs. We calibrate the model against NEO detections by Catalina Sky Surveys' stations 703 and G96 during 2005-2012, and utilize the complementary nature of these two systems to quantify the systematic uncertainties associated to the resulting model. The most important ERs are the $\nu_6$ and 3:1J resonance complexes with JFCs contributing a few percent of NEOs on average. We predict that there are $962^{+52}_{-56}$ ($802^{+48}_{-42}\times10^3$) NEOs with $H<17.75$ ($H<25$) and these numbers are in agreement with the most recent estimates found in the literature (the uncertainty estimates only account for the random component).

No planets exist inside the orbit of Mercury and the terrestrial planets of the solar system exhibit a localized configuration. According to thermal structure calculation of protoplanetary disks, a silicate condensation line (~ 1300 K) is located around 0.1 au from the Sun except for the early phase of disk evolution, and planetesimals could have formed inside the orbit of Mercury. A recent study of disk evolution that includes magnetically driven disk winds showed that the gas disk obtains a positive surface density slope inside ~ 1 au from the central star. In a region with positive midplane pressure gradient, planetesimals undergo outward radial drift. We investigate the radial drift of planetesimals and type I migration of planetary embryos in a disk that viscously evolves with magnetically driven disk winds. We show a case in which no planets remain in the close-in region. Radial drifts of planetesimals are simulated using a recent disk evolution model that includes effects of disk winds. The late stage of planet formation is also examined by performing N-body simulations of planetary embryos. We demonstrate that in the middle stage of disk evolution, planetesimals can undergo convergent radial drift in a magnetorotational instability (MRI)-inactive disk, in which the pressure maximum is created, and accumulate in a narrow ring-like region with an inner edge at ~ 0.7 au from the Sun. We also show that planetary embryos that may grow from the narrow planetesimal ring do not exhibit significant type I migration in the late stage of disk evolution. The origin of the localized configuration of the terrestrial planets of the solar system, in particular the deficit of close-in planets, can be explained by the convergent radial drift of planetesimals in disks with a positive pressure gradient in the close-in region.

Planets with masses larger than about 0.1 Earth-masses undergo rapid inward migration (type I migration) in a standard protoplanetary disk. Recent magnetohydrodynamical simulations revealed the presence of magnetically driven disk winds, which would alter the disk profile and the type I migration in the close-in region. We investigate orbital evolution of planetary embryos in disks that viscously evolve under the effects of disk winds. The aim is to discuss effects of altered disk profiles on type I migration. In addition, we aim to examine whether observed distributions of close-in super-Earths can be reproduced by simulations that include effects of disk winds. We perform N-body simulations of super-Earth formation from planetary embryos, in which a recent model for disk evolution is used. We explore a wide range of parameters and draw general trends. We also carry out N-body simulations of close-in super-Earth formation from embryos in such disks under various conditions. We find that the type I migration is significantly suppressed in many cases. Even in cases in which inward migration occurs, the migration timescale is lengthened to 1 Myr, which mitigates the type I migration problem. This is because the gas surface density is decreased and has a flatter profile in the close-in region due to disk winds. We find that when the type I migration is significantly suppressed, planets undergo late orbital instability during the gas depletion, leading to a non-resonant configuration. We also find that observed distributions of close-in super-Earths (e.g., period ratio, mass ratio) can be reproduced. In addition, we show that in some results of simulations, systems with a chain of resonant planets, like the TRAPPIST-1 system, form.

Planetary systems can evolve dynamically even after the full growth of the planets themselves. There is actually circumstantial evidence that most planetary systems become unstable after the disappearance of gas from the protoplanetary disk. These instabilities can be due to the original system being too crowded and too closely packed or to external perturbations such as tides, planetesimal scattering, or torques from distant stellar companions. The Solar System was not exceptional in this sense. In its inner part, a crowded system of planetary embryos became unstable, leading to a series of mutual impacts that built the terrestrial planets on a timescale of ~100 My. In its outer part, the giant planets became temporarily unstable and their orbital configuration expanded under the effect of mutual encounters. A planet might have been ejected in this phase. Thus, the orbital distributions of planetary systems that we observe today, both solar and extrasolar ones, can be different from the those emerging from the formation process and it is important to consider possible long-term evolutionary effects to connect the two.

In planetary science, accretion is the process in which solids agglomerate to form larger and larger objects and eventually planets are produced. The initial conditions are a disc of gas and microscopic solid particles, with a total mass of about 1% of the gas mass. Accretion has to be effective and fast. Effective, because the original total mass in solids in the solar protoplanetary disk was probably of the order of ~ 300 Earth masses, and the mass incorporated into the planets is ~100 Earth masses. Fast, because the cores of the giant planets had to grow to tens of Earth masses in order to capture massive doses of hydrogen and helium from the disc before the dispersal of the latter, i.e. in a few millions of years. There is probably not one accretion process but several, depending on the scale at which accretion operates. A first process is the sticking of microscopic dust into larger grains and pebbles. A second process is the formation of an intermediate class of objects called planetesimals. A third accretion process has to lead from planetesimals to planets. Actually, several processes can be involved in this step, from collisional coagulation among planetesimals to the accretion of small particles under the effect of gas drag, to giant impacts between protoplanets. This chapter will detail all these processes, adopting a historical perspective: i.e. from the classic processes investigated in the past decades to those unveiled in the last years.

The thermal properties of the surfaces of asteroids determine the magnitude of the drift rate cause by the Yarkovsky force. In the general case of Main Belt asteroids, the Yarkovsky force is indirectly proportional to the thermal inertia, $\Gamma$. Following the proposed relationship between $\Gamma$ and asteroid diameter $D$, we find that asteroids' Yarkovsky drift rates might have a more complex size dependence than previous thought, leading to a curved family V-shape boundary in semi-major axis, a, vs. 1/$D$ space. This implies that asteroids are drifting faster at larger sizes than previously considered decreasing on average the known ages of asteroid families. The V-Shape curvature is determined for $>$25 families located throughout the Main Belt to quantify the Yarkovsky size-dependent drift rate. We find that there is no correlation between family age and V-shape curvature. In addition, the V-shape curvature decreases for asteroid families with larger heliocentric distances suggesting that the relationship between $\Gamma$ and $D$ is weaker in the outer MB possibly due to homogenous surface roughness among family members.

The growth of a planetary core by pebble accretion stops at the so called pebble isolation mass, when the core generates a pressure bump that traps drifting pebbles outside its orbit. If the isolation mass is very small, then gas accretion is protracted and the planet remains at a few Earth masses with a mainly solid composition. For larger values of the pebble isolation mass, the planet might be able to accrete gas from the protoplanetary disc and grow into a gas giant. Previous works have determined a scaling of the pebble isolation mass with cube of the disc aspect ratio. Here we expand on previous measurements and explore the dependency of the pebble isolation mass on all relevant parameters of the protoplanetary disc. We use 3D hydrodynamical simulations to measure the pebble isolation mass and derive a simple scaling law that captures the dependence on the local disc structure and the turbulent viscosity parameter $\alpha$. We find that small pebbles, coupled to the gas, with Stokes number $\tau_{\rm f}<0.005$ can drift through the partial gap at pebble isolation mass. However, as the planetary mass increases, particles must be decreasingly smaller to penetrate through the pressure bump. Turbulent diffusion of particles, however, can lead to an increase of the pebble isolation mass by a factor of two, depending on the strength of the background viscosity and on the pebble size. We finally explore the implications of the new scaling law of the pebble isolation mass on the formation of planetary systems by numerically integrating the growth and migration pathways of planets in evolving protoplanetary discs. Compared to models neglecting the dependence of the pebble isolation mass on the $\alpha$-viscosity, our models including this effect result in larger core masses for giant planets. These larger core masses are more akin to the core masses of the giant planets in the Solar System.

The timeline of the lunar bombardment in the first Gy of the Solar System remains unclear. Some basin-forming impacts occurred 3.9-3.7Gy ago. Many other basins formed before, but their exact ages are not precisely known. There are two possible interpretations of the data: in the cataclysm scenario there was a surge in the impact rate approximately 3.9Gy ago, while in the accretion tail scenario the lunar bombardment declined since the era of planet formation and the latest basins formed in its tail-end. Here, we revisit the work of Morbidelli et al.(2012) that examined which scenario could be compatible with both the lunar crater record in the 3-4Gy period and the abundance of highly siderophile elements (HSE) in the lunar mantle. We use updated numerical simulations of the fluxes of impactors. Under the traditional assumption that the HSEs track the total amount of material accreted by the Moon since its formation, we conclude that only the cataclysm scenario can explain the data. The cataclysm should have started ~3.95Gy ago. However we show that HSEs could have been sequestered from the lunar mantle due to iron sulfide exsolution during magma ocean crystallization, followed by mantle overturn. Based on the hypothesis that the lunar magma ocean crystallized about 100-150My after Moon formation, and therefore that HSEs accumulated in the lunar mantle only after this time, we show that the bombardment in the 3-4Gy period can be explained in the accretion tail scenario. This hypothesis would also explain why the Moon appears so depleted in HSEs relative to the Earth. We also extend our analysis of the cataclysm and accretion tail scenarios to the case of Mars. The accretion tail scenario requires a global resurfacing event on Mars ~4.4Gy ago, possibly associated with the formation of the Borealis basin, and it is consistent with the HSE budget of the planet.

The origin of the orbital structure of the cold component of the Kuiper belt is still a hot subject of investigation. Several features of the solar system suggest that the giant planets underwent a phase of global dynamical instability, but the actual dynamical evolution of the planets during the instability is still debated. To explain the structure of the cold Kuiper belt, Nesvorny (2015, AJ 150,68) argued for a soft instability, during which Neptune never achieved a very eccentric orbit. Here we investigate the possibility of a more violent instability, from an initially more compact fully resonant configuration of 5 giant planets. We show that the orbital structure of the cold Kuiper belt can be reproduced quite well provided that the cold population formed in situ, with an outer edge between 44 and 45 au and never had a large mass.

The observational census of trans-Neptunian objects with semi-major axes greater than ~250 AU exhibits unexpected orbital structure that is most readily attributed to gravitational perturbations induced by a yet-undetected, massive planet. Although the capacity of this planet to (i) reproduce the observed clustering of distant orbits in physical space, (ii) facilitate dynamical detachment of their perihelia from Neptune, and (iii) excite a population of long-period centaurs to extreme inclinations is well established through numerical experiments, a coherent theoretical description of the dynamical mechanisms responsible for these effects remains elusive. In this work, we characterize the dynamical processes at play, from semi-analytic grounds. We begin by considering a purely secular model of orbital evolution induced by Planet Nine, and show that it is at odds with the ensuing stability of distant objects. Instead, the long-term survival of the clustered population of long-period KBOs is enabled by a web of mean-motion resonances driven by Planet Nine. Then, by taking a compact-form approach to perturbation theory, we show that it is the secular dynamics embedded within these resonances that regulates the orbital confinement and perihelion detachment of distant Kuiper belt objects. Finally, we demonstrate that the onset of large-amplitude oscillations of orbital inclinations is accomplished through capture of low-inclination objects into a high-order secular resonance and identify the specific harmonic that drives the evolution. In light of the developed qualitative understanding of the governing dynamics, we offer an updated interpretation of the current observational dataset within the broader theoretical framework of the Planet Nine hypothesis.

Earth possesses a persistent, internally-generated magnetic field, whereas no trace of a dynamo has been detected on Venus, at present or in the past, although a high surface temperature and recent resurfacing events may have removed paleomagnetic evidence. Whether or not a terrestrial body can sustain an internally generated magnetic field by convection inside its metallic fluid core is determined in part by its initial thermodynamic state and its compositional structure, both of which are in turn set by the processes of accretion and differentiation. Here we show that the cores of Earth- and Venus-like planets should grow with stable compositional stratification unless disturbed by late energetic impacts. They do so because higher abundances of light elements are incorporated into the liquid metal that sinks to form the core as the temperatures and pressures of metal-silicate equilibration increase during accretion. We model this process and determine that this establishes a stable stratification that resists convection and inhibits the onset of a geodynamo. However, if a late energetic impact occurs, it could mechanically stir the core creating a single homogenous region within which a long-lasting geodynamo would operate. While Earth's accretion has been punctuated by a late giant impact with likely enough energy to mix the core (e.g. the impact that formed the Moon), we hypothesize that the accretion of Venus is characterized by the absence of such energetic giant impacts and the preservation of its primordial stratifications.

Magnetically-driven disk winds would alter the surface density slope of gas in the inner region of a protoplanetary disk $(r \lesssim 1 {\rm au})$. This in turn affects planet formation. Recently, the effect of disk wind torque has been considered with the suggestion that it would carve out the surface density of the disk from inside and would induce global gas flows (wind-driven accretion). We aim to investigate effects of global gas flows on type I migration and also examine planet formation. A simplified approach was taken to address this issue, and N-body simulations with isolation-mass planets were also performed. In previous studies, the effect of gas flow induced by turbulence-driven accretion has been taken into account for its desaturation effect of the corotation torque. If more rapid gas flows (e.g., wind-driven accretion) are considered, the desaturation effect can be modified. In MRI-inactive disks, in which the wind-driven accretion dominates the disk evolution, the gas flow at the midplane plays an important role. If this flow is fast, the corotation torque is efficiently desaturated. Then, the fact that the surface density slope can be positive in the inner region due to the wind torque can generate an outward migration region extended to super-Earth mass planets. In this case, we observe that no planets fall onto the central star in N-body simulations with migration forces imposed to reproduce such migration pattern. We also see that super-Earth mass planets can undergo outward migration. Relatively rapid gas flows affects type I migration and thus the formation of close-in planets.

Seven of the nine known Mars Trojan asteroids belong to an orbital cluster named after its largest member 5261 Eureka. Eureka is likely the progenitor of the whole cluster, which formed at least 1 Gyr ago. It was suggested that the thermal YORP effect spun-up Eureka resulting with fragments being ejected by the rotational-fission mechanism. Eureka's spectrum exhibits a broad and deep absorption band around 1 {\mu}m, indicating an olivine-rich composition. Here we show evidence that the Trojan Eureka cluster progenitor could have originated as impact debris excavated from the Martian mantle. We present new near-infrared observations of two Trojans (311999 2007 NS2 and 385250 2001 DH47) and find that both exhibit an olivine-rich reflectance spectrum similar to Eureka's. These measurements confirm that the progenitor of the cluster has an achondritic composition. Olivine-rich reflectance spectra are rare amongst asteroids but are seen around the largest basins on Mars. They are also consistent with some Martian meteorites (e.g. Chassigny), and with the material comprising much of the Martian mantle. Using numerical simulations, we show that the Mars Trojans are more likely to be impact ejecta from Mars than captured olivine-rich asteroids transported from the main belt. This result links directly specific asteroids to debris from the forming planets.

Ejection velocity fields of asteroid families are largely unconstrained due to the fact that members disperse relatively quickly on Myr time-scales by secular resonances and the Yarkovsky effect. The spreading of fragments in $a$ by the Yarkovsky effect is indistinguishable from the spreading caused by the initial ejection of fragments. By examining families $<$20 Myrs-old, we can use the V-shape identification technique to separate family shapes that are due to the initial ejection velocity field and those that are due to the Yarkovsky effect. $<$20 Myr-old asteroid families provide an opportunity to study the velocity field of family fragments before they become too dispersed. Only the Karin family's initial velocity field has been determined and scales inversely with diameter, $D^{-1}$. We have applied the V-shape identification technique to constrain young families' initial ejection velocity fields by measuring the curvature of their fragments' V-shape correlation in semi-major axis, $a$, vs. $D^{-1}$ space. Curvature from a straight line implies a deviation from a scaling of $D^{-1}$. We measure the V-shape curvature of 11 young asteroid families including the \FYnospace, Aeolia, Brangane, Brasilia, Clarissa, Iannini, Karin, Konig, Koronis(2), Theobalda and Veritas asteroid families. We find that the majority of asteroid families have initial ejection velocity fields consistent with $\sim D^{-1}$ supporting laboratory impact experiments and computer simulations of disrupting asteroid parent bodies.

Comets are icy objects that orbitally evolve from the trans-Neptunian region (the Kuiper belt and beyond) into the inner Solar System, where they are heated by solar radiation and become active due to sublimation of water ice. Here we perform end-to-end simulations in which cometary reservoirs are formed in the early Solar System and evolved over 4.5 Gyr. The gravitational effects of Planet 9 (P9), hypothesized to circle the Sun on a wide orbit, are included in some of our simulations. Different models are considered for comets to be active, including a simple assumption that comets remain active for Np(q) perihelion passages with perihelion distance q<2.5 au. The orbital distribution and number of active comets produced in our model is compared to observations. The orbital distribution of ecliptic comets (ECs) is well reproduced in models with Np(2.5)=500 and without P9. With P9, the inclination distribution of model ECs is wider than the observed one. We find that the known Halley-type comets (HTCs) have a nearly isotropic inclination distribution (with only a slight preference for prograde orbits). In our model, the HTCs appear to be an extension of the population of returning Oort-cloud comets (OCCs) to shorter orbital periods. The inclination distribution of model HTCs becomes broader with increasing Np, but the existing observational data are not good enough to constrain Np from orbital fits. Np(2.5)>1000 is required to obtain a steady-state population of large active HTCs that is consistent with observations. To fit the ratio of the returning-to-new OCCs, by contrast, our model implies that Np(2.5)<10, possibly because the detected long-period comets are smaller and much easier to disrupt than observed HTCs.

In this work we aim to constrain the slope of the size distribution of main-belt asteroids, at their primordial state. To do so we turn out attention to the part of the main asteroid belt between 2.82 and 2.96~AU, the so-called "pristine zone", which has a low number density of asteroids and few, well separated asteroid families. Exploiting these unique characteristics, and using a modified version of the hierarchical clustering method we are able to remove the majority of asteroid family members from the region. The remaining, background asteroids should be of primordial origin, as the strong 5/2 and 7/3 mean-motion resonances with Jupiter inhibit transfer of asteroids to and from the neighboring regions. The size-frequency distribution of asteroids in the size range $17<D(\rm{km})<70$ has a slope $q\simeq-1$. Using Monte-Carlo methods, we are able to simulate, and compensate for the collisional and dynamical evolution of the asteroid population, and get an upper bound for its size distribution slope $q=-1.43$. In addition, applying the same 'family extraction' method to the neighboring regions, i.e. the middle and outer belts, and comparing the size distributions of the respective background populations, we find statistical evidence that no large asteroid families of primordial origin had formed in the middle or pristine zones.

"Hot super-Earths" (or "Mini-Neptunes") between 1 and 4 times Earth's size with period shorter than 100 days orbit 30-50\% of Sun-like type stars. Their orbital configuration -- measured as the period ratio distribution of adjacent planets in multi-planet systems -- is a strong constraint for formation models. Here we use N-body simulations with synthetic forces from an underlying evolving gaseous disk to model the formation and long-term dynamical evolution of super-Earth systems. While the gas disk is present, planetary embryos grow and migrate inward to form a resonant chain anchored at the inner edge of the disk. These resonant chains are far more compact than the observed super-Earth systems. Once the gas dissipates resonant chains may become dynamically unstable. They undergo a phase of giant impacts that spreads the systems out. Disk turbulence has no measurable effect on the outcome. Our simulations match observations if a small fraction of resonant chains remain stable, while most super-Earths undergo a late dynamical instability. Our statistical analysis restricts the contribution of stable systems to less than $25\%$. Our results also suggest that the large fraction of observed single planet systems does not necessarily imply any dichotomy in the architecture of planetary systems. Finally, we use the low abundance of resonances in Kepler data to argue that, in reality, the survival of resonant chains happens likely only in $\sim 5\%$ of the cases. This leads to a mystery: in our simulations only 50-60\% of resonant chains became unstable whereas at least 75\% (and probably 90-95\%) must be unstable to match observations.

It is well known that asteroids and comets fall into the Sun. Metal pollution of white dwarfs and transient spectroscopic signatures of young stars like $\beta$-Pic provide growing evidence that extra solar planetesimals can attain extreme orbital eccentricities and fall onto their parent stars. We aim to develop a general, practically implementable, semi-analytical theory of secular eccentricity excitation of small bodies in mean motion resonances with an eccentric planet valid for arbitrary values of the eccentricities and including the short-range force due to General Relativity. Our semi-analytic model for the restricted planar three-body problem does not make use of any series expansion and therefore is valid for any values of eccentricities and semi-major axes ratios. The model is based on the application of the adiabatic principle, which is valid when the precession period of the longitude of pericenter of the planetesimal is much longer than the libration period in the mean motion resonance. This holds down to vanishingly small eccentricities in resonances of order larger than 1. We provide a Mathematica notebook with the implementation of the model allowing direct use to the interested reader. We confirm that the 4:1 mean motion resonance with a moderately eccentric planet is the most powerful one to lift the eccentricity of planetesimals from nearly circular orbits to star-grazing ones. However, if the planet is too eccentric, we find that this resonances becomes unable to pump the planetesimal's eccentricity to very high value. The inclusion of the General Relativity effect imposes a condition on the mass of the planet to drive the planetsimals into star-grazing orbits. For a planetesimal at $\sim$1 AU around a solar-mass star (or white dwarf), we find a threshold planetary mass of about 17 Earth masses. We finally derive an analytical formula for this critical mass.

Recent works on planetary migration show that the orbital structure of the Kuiper belt can be very well reproduced if before the onset of the planetary instability Neptune underwent a long-range planetesimal-driven migration up to $\sim$28 au. However, considering that all giant planets should have been captured in mean motion resonances among themselves during the gas-disk phase, it is not clear whether such a very specific evolution for Neptune is possible, nor whether the instability could have happened at late times. Here, we first investigate which initial resonant configuration of the giant planets can be compatible with Neptune being extracted from the resonant chain and migrating to $\sim$28 au before that the planetary instability happened. We address the late instability issue by investigating the conditions where the planets can stay in resonance for about 400 My. Our results indicate that this can happen only in the case where the planetesimal disk is beyond a specific minimum distance $\delta_{stab}$ from Neptune. Then, if there is a sufficient amount of dust produced in the planetesimal disk, that drifts inwards, Neptune can enter in a slow dust-driven migration phase for hundreds of Mys until it reaches a critical distance $\delta_{mig}$ from the disk. From that point, faster planetesimal-driven migration takes over and Neptune continues migrating outward until the instability happens. We conclude that, although an early instability reproduces more easily the evolution of Neptune required to explain the structure of the Kuiper belt, such evolution is also compatible with a late instability.

The Asteroid Belt is characterized by the radial mixing of bodies with different physical properties, a very low mass compared to Minimum Mass Solar Nebula expectations and has an excited orbital distribution. Models of the evolution of the Asteroid Belt show that the origin of its structure is strongly linked to the process of terrestrial planet formation. The Grand Tack model presents a possible solution to the conundrum of reconciling the small mass of Mars with the properties of the Asteroid Belt, including the mass depletion, radial mixing and orbital excitation. However, while the inclination distribution produced in the Grand Tack model is in good agreement with the one observed, the eccentricity distribution is skewed towards values larger than those found today. Here, we evaluate the evolution of the orbital properties of the Asteroid Belt from the end of the Grand Tack model (at the end of the gas nebula phase when planets emerge from the dispersing gas disk), throughout the subsequent evolution of the Solar System including an instability of the Giant Planets approximately 400 My later. Before the instability, the terrestrial planets were modeled on dynamically cold orbits with Jupiter and Saturn locked in a 3:2 mean motion resonance. The model continues for an additional 4.1 Gy after the giant planet instability. Our results show that the eccentricity distribution obtained in the Grand Tack model evolves towards one very similar to that currently observed, and the semimajor axis distribution does the same. The inclination distribution remains nearly unchanged with a slight preference for depletion at low inclination; this leads to the conclusion that the inclination distribution at the end of the Grand Tack is a bit over-excited.

There is an active debate about whether the properties of comets as observed today are primordial or, alternatively, if they are a result of collisional evolution or other processes. We investigate the effects of collisions on a comet with a structure like 67P/C-G. We develop scaling laws for the critical specific impact energies required for a significant shape alteration. These are then used in simulations of the combined dynamical and collisional evolution of comets in order to study the survival probability of a primordially formed object with a shape like 67P/C-G. The effects of impacts on comet 67P/C-G are studied using a SPH shock physics code. The resulting critical specific impact energy defines a minimal projectile size which is used to compute the number of shape-changing collisions in a set of dynamical simulations. These simulations follow the dispersion of the trans-Neptunian disk during the giant planet instability, the formation of a scattered disk, and produce 87 objects that penetrate into the inner solar system with orbits consistent with the observed JFC population. The collisional evolution before the giant planet instability is not considered here. Hence, our study is conservative in its estimation of the number of collisions. We find that in any scenario considered here, comet 67P/C-G would have experienced a significant number of shape-changing collisions, if it formed primordially. This is also the case for generic bi-lobe shapes. Our study also shows that impact heating is very localized and that collisionally processed bodies can still have a high porosity. Our study indicates that the observed bi-lobe structure of comet 67P/C-G may not be primordial, but might have originated in a rather recent event, possibly within the last 1 Gy. This may be the case for any kilometer-sized two-component cometary nuclei.

Over the past two decades, large strides have been made in the field of planet formation. Yet fundamental questions remain. Here we review our state of understanding of five fundamental bottlenecks in planet formation. These are: 1) the structure and evolution of protoplanetary disks; 2) the growth of the first planetesimals; 3) orbital migration driven by interactions between proto-planets and gaseous disk; 4) the origin of the Solar System's orbital architecture; and 5) the relationship between observed super-Earths and our own terrestrial planets. Given our lack of understanding of these issues, even the most successful formation models remain on shaky ground.

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.

There are only a few known main belt (MB) asteroid families with ages greater than 2 Gyr (Bro\v{z} 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.

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.

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.

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.

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.

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.

We present three-dimensional simulations with nested meshes of the dynamics of the gas around a Jupiter mass planet with the JUPITER and FARGOCA codes. We implemented a radiative transfer module into the JUPITER code to account for realistic heating and cooling of the gas. We focus on the circumplanetary gas flow, determining its characteristics at very high resolution ($80\%$ of Jupiter's diameter). In our nominal simulation where the temperature evolves freely by the radiative module and reaches 13000 K at the planet, a circumplanetary envelope was formed filling the entire Roche-lobe. Because of our equation of state is simplified and probably overestimates the temperature, we also performed simulations with limited maximal temperatures in the planet region (1000 K, 1500 K, and 2000 K). In these fixed temperature cases circumplanetary disks (CPDs) were formed. This suggests that the capability to form a circumplanetary disk is not simply linked to the mass of the planet and its ability to open a gap. Instead, the gas temperature at the planet's location, which depends on its accretion history, plays also fundamental role. The CPDs in the simulations are hot and cooling very slowly, they have very steep temperature and density profiles, and are strongly sub-Keplerian. Moreover, the CPDs are fed by a strong vertical influx, which shocks on the CPD surfaces creating a hot and luminous shock-front. In contrast, the pressure supported circumplanetary envelope is characterized by internal convection and almost stalled rotation.

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.

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.

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.

Earth and Moon are shown here to be composed of oxygen isotope reservoirs that are indistinguishable, with a difference in {\Delta}"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 {\Delta}"17O within approximately 1 per mil of the terrestrial value, suggesting that these impactors were water rich.

The terrestrial planets and the asteroids dominant in the inner asteroid belt are water poor. However, in the protoplanetary disk the temperature should have decreased below water condensation level well before the disk was photoevaporated. Thus, the global water depletion of the inner Solar System is puzling. We show that, even if the inner disk becomes cold, there cannot be direct condensation of water. This is because the snowline moves towards the Sun more slowly than the gas itself. The appearance of ice in a range of heliocentric distances swept by the snowline can only be due to the radial drift of icy particles from the outer disk. However, if a sufficiently massive planet is present, the radial drift of particles is interrupted, because the disk acquires a superKeplerian rotation just outside of the planetary orbit. From this result, we propose that the precursor of Jupiter achieved about 20 Earth masses when the snowline was still around 3 AU. This effectively fossilized the snowline at that location. Although cooling, the disk inside of the Jovian orbit remained ice-depleted because the flow of icy particles from the outer system was intercepted by the planet. This scenario predicts that planetary systems without giant planets should be much more rich in water in their inner regions than our system. We also show that the inner edge of the planetesimal disk at 0.7AU, required in terrestrial planet formation models to explain the small mass of Mercury and the absence of planets inside of its orbit, could be due to the silicate condensation line, fossilized at the end of the phase of streaming instability that generated the planetesimal seeds. Thus, when the disk cooled, silicate particles started to drift inwards of 0.7AU without being sublimated, but they could not be accreted by any pre-existing planetesimals.

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.

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.

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.

For the very first time we present 3D simulations of planets embedded in stellar irradiated discs. It is well known that thermal effects could reverse the direction of planetary migration from inwards to outwards, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. When considering stellar irradiation in addition to viscous friction as a source of heating, the outer disc changes from a shadowed to a flared structure. Using a suited analytical formula it has been shown that in the flared part of the disc the migration is inwards; planets can migrate outwards only in shadowed regions of the disc, { because the radial gradient of entropy is stronger there}. In order to confirm this result numerically, we have computed the total torque acting on planets held on fixed orbits embedded in stellar irradiated 3D discs using the hydrodynamical code FARGOCA. We find qualitatively good agreement between the total torque obtained with numerical simulations and the one predicted by the analytical formula. For large masses (>20 Earth masses) we find quantitative agreement, and we obtain outwards migration regions for planets up to 60 Earth masses in the early stages of accretional discs. We find nevertheless that the agreement with the analytic formula is quite fortuitous because the formula underestimates the size of the horseshoe region; this error is compensated by imperfect estimates of other terms, most likely the cooling rate and saturation.

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.

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.

The basic structure of the solar system is set by the presence of low-mass terrestrial planets in its inner part and giant planets in its outer part. This is the result of the formation of a system of multiple embryos with approximately the mass of Mars in the inner disk and of a few multi-Earth-mass cores in the outer disk, within the lifetime of the gaseous component of the protoplanetary disk. What was the origin of this dichotomy in the mass distribution of embryos/cores? We show in this paper that the classic processes of runaway and oligarchic growth from a disk of planetesimals cannot explain this dichotomy, even if the original surface density of solids increased at the snowline. Instead, the accretion of drifting pebbles by embryos and cores can explain the dichotomy, provided that some assumptions hold true. We propose that the mass-flow of pebbles is two-times lower and the characteristic size of the pebbles is approximately ten times smaller within the snowline than beyond the snowline (respectively at heliocentric distance $r<r_{ice}$ and $r>r_{ice}$, where $r_{ice}$ is the snowline heliocentric distance), due to ice sublimation and the splitting of icy pebbles into a collection of chondrule-size silicate grains. In this case, objects of original sub-lunar mass would grow at drastically different rates in the two regions of the disk. Within the snowline these bodies would reach approximately the mass of Mars while beyond the snowline they would grow to $\sim 20$ Earth masses. The results may change quantitatively with changes to the assumed parameters, but the establishment of a clear dichotomy in the mass distribution of protoplanets appears robust, provided that there is enough turbulence in the disk to prevent the sedimentation of the silicate grains into a very thin layer.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 $\alpha$ 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 ($\alpha \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 $\alpha \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.

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 ($\tau_{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.

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.

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.

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.

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.

Vesta has a complex cratering history, with ancient terrains as well as recent large impacts that have led to regional resurfacing. Crater counts can help constrain the relative ages of different units on Vesta's surface, but converting those crater counts to absolute ages requires a chronology function. We present a cratering chronology based on the best current models for the dynamical evolution of asteroid belt, and calibrate it to Vesta using the record of large craters on its surface. While uncertainties remain, our chronology function is broadly consistent with an ancient surface of Vesta as well as other constraints such as the bombardment history of the rest of the inner Solar System and the Ar-Ar age distribution of howardite, eucrite and diogenite (HED) meteorites from Vesta.

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.

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.

In this paper, we address the migration of small mass planets in 3D radiative disks. Indeed, migration of small planets is known to be too fast inwards in locally isothermal conditions. However, thermal effects could reverse its direction, potentially saving planets in the inner, optically thick parts of the protoplanetary disc. This effect has been seen for masses larger than 5 Earth masses, but the minimum mass for this to happen has never been probed numerically, although it is of crucial importance for planet formation scenarios. We have extended the hydro-dynamical code FARGO to 3D, with thermal diffusion. With this code, we perform simulations of embedded planets down to 2 Earth masses. For a set of discs parameters for which outward migration has been shown in the range of $[5, 35]$ Earth masses, we find that the transition to inward migration occurs for masses in the range $[3, 5]$ Earth masses. The transition appears to be due to an unexpected phenomenon: the formation of an asymmetric cold and dense finger of gas driven by circulation and libration streamlines. We recover this phenomenon in 2D simulations where we control the cooling effects of the gas through a simple modeling of the energy equation.

We examine the gas circulation near a gap opened by a giant planet in a protoplanetary disk. We show with high resolution 3D simulations that the gas flows into the gap at high altitude over the mid-plane, at a rate dependent on viscosity. We explain this observation with a simple conceptual model. From this model we derive an estimate of the amount of gas flowing into a gap opened by a planet with Hill radius comparable to the scale-height of a layered disk (i. e. a disk with viscous upper layer and inviscid midplane). Our estimate agrees with modern MRI simulations(Gressel et al., 2013). We conclude that gap opening in a layered disk can not slow down significantly the runaway gas accretion of Saturn to Jupiter-mass planets.

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 \nu 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.

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, $\alpha_{\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, $\alpha_{\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 $\alpha_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 $\alpha_\textrm{1}=1.0\pm0.2$, $\alpha_\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.

In the core-accretion model the nominal runaway gas-accretion phase brings most planets to multiple Jupiter masses. However, known giant planets are predominantly Jupiter-mass bodies. Obtaining longer timescales for gas accretion may require using realistic equations of states, or accounting for the dynamics of the circumplanetary disk (CPD) in low-viscosity regime, or both. Here we explore the second way using global, three-dimensional isothermal hydrodynamical simulations with 8 levels of nested grids around the planet. In our simulations the vertical inflow from the circumstellar disk (CSD) to the CPD determines the shape of the CPD and its accretion rate. Even without prescribed viscosity Jupiter's mass-doubling time is $\sim 10^4$ years, assuming the planet at 5.2 AU and a Minimum Mass Solar Nebula. However, we show that this high accretion rate is due to resolution-dependent numerical viscosity. Furthermore, we consider the scenario of a layered CSD, viscous only in its surface layer, and an inviscid CPD. We identify two planet-accretion mechanisms that are independent of the viscosity in the CPD: (i) the polar inflow -- defined as a part of the vertical inflow with a centrifugal radius smaller than 2 Jupiter-radii and (ii) the torque exerted by the star on the CPD. In the limit of zero effective viscosity, these two mechanisms would produce an accretion rate 40 times smaller than in the simulation.

Dawson and Murray-Clay (2012) pointed out that the inner part of the cold population in the Kuiper belt (that with semi major axis a<43.5 AU) has orbital eccentricities significantly smaller than the limit imposed by stability constraints. Here, we confirm their result by looking at the orbital distribution and stability properties in proper element space. We show that the observed distribution could have been produced by the slow sweeping of the 4/7 mean motion resonance with Neptune that accompanied the end of Neptune's migration process. The orbital distribution of the hot Kuiper belt is not significantly affected in this process, for the reasons discussed in the main text. Therefore, the peculiar eccentricity distribution of the inner cold population can not be unequivocally interpreted as evidence that the cold population formed in-situ and was only moderately excited in eccentricity; it can simply be the signature of Neptune's radial motion, starting from a moderately eccentric orbit. We discuss how this agrees with a scenario of giant planet evolution following a dynamical instability and, possibly, with the radial transport of the cold population.

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.

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.

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.

Debris discs are a consequence of the planet formation process and constitute the fingerprints of planetesimal systems. Their solar system's counterparts are the asteroid and Edgeworth-Kuiper belts. The DUNES survey aims at detecting extra-solar analogues to the Edgeworth-Kuiper belt around solar-type stars, putting in this way the solar system into context. The survey allows us to address some questions related to the prevalence and properties of planetesimal systems. We used {\it Herschel}/PACS to observe a sample of nearby FGK stars. Data at 100 and 160 $\mu$m were obtained, complemented in some cases with observations at 70 $\mu$m, and at 250, 350 and 500 $\mu$m using SPIRE. The observing strategy was to integrate as deep as possible at 100 $\mu$m to detect the stellar photosphere. Debris discs have been detected at a fractional luminosity level down to several times that of the Edgeworth-Kuiper belt. The incidence rate of discs around the DUNES stars is increased from a rate of $\sim$ 12.1% $\pm$ 5% before \emph{Herschel} to $\sim$ 20.2% $\pm$ 2%. A significant fraction ($\sim$ 52%) of the discs are resolved, which represents an enormous step ahead from the previously known resolved discs. Some stars are associated with faint far-IR excesses attributed to a new class of cold discs. Although it cannot be excluded that these excesses are produced by coincidental alignment of background galaxies, statistical arguments suggest that at least some of them are true debris discs. Some discs display peculiar SEDs with spectral indexes in the 70-160$\mu$m range steeper than the Rayleigh-Jeans one. An analysis of the debris disc parameters suggests that a decrease might exist of the mean black body radius from the F-type to the K-type stars. In addition, a weak trend is suggested for a correlation of disc sizes and an anticorrelation of disc temperatures with the stellar age.

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.

One of the outstanding problems of the dynamical evolution of the outer solar system concerns the observed population ratio between the Oort Cloud (OC) and the Scattered Disc (SD): observations suggest that this ratio lies between 100 and 1000 but simulations that produce these two reservoirs simultaneously consistently yield a value of the order of 10. Here we stress that the populations in the OC and SD are inferred from the observed fluxes of new Long Period Comets (LPCs) and Jupiter-family comets (JFCs), brighter than some reference total magnitude. However, the population ratio estimated in the simulations of formation of the SD and OC refers to objects bigger than a given size. There are multiple indications that LPCs are intrinsically brighter than JFCs, i.e. an LPC is smaller than a JFC with the same total absolute magnitude. When taking this into account we revise the SD/JFC population ratio from our simulations relative to Duncan and Levison (1997), and then deduce from the observations that the size-limited population ratio between the OC and the SD is 44 (-34)(+54). From simulations we obtain 12 +- 1 but the agreement cannot be rejected by the null hypothesis.

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.

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\'y 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.

In the Nice model, the late heavy bombardment (LHB) is related to an orbital instability of giant planets which causes a fast dynamical dispersion of a transneptunian cometary disk. We study effects produced by these hypothetical cometary projectiles on main-belt asteroids. In particular, we want to check whether the observed collisional families provide a lower or an upper limit for the cometary flux during the LHB. We present an updated list of observed asteroid families as identified in the space of synthetic proper elements by the hierarchical clustering method, colour data, albedo data and dynamical considerations and we estimate their physical parameters. We selected 12 families which may be related to the LHB according to their dynamical ages. We then used collisional models and N-body orbital simulations to gain insight into the long-term dynamical evolution of synthetic LHB families over 4 Gyr. We account for the mutual collisions, the physical disruptions of comets, the Yarkovsky/YORP drift, chaotic diffusion, or possible perturbations by the giant-planet migration. Assuming a "standard" size-frequency distribution of primordial comets, we predict the number of families with parent-body sizes D_PB >= 200 km which seems consistent with observations. However, more than 100 asteroid families with D_PB >= 100 km should be created at the same time which are not observed. This discrepancy can be nevertheless explained by the following processes: i) asteroid families are efficiently destroyed by comminution (via collisional cascade), ii) disruptions of comets below some critical perihelion distance (q <~ 1.5 AU) are common. Given the freedom in the cometary-disruption law, we cannot provide stringent limits on the cometary flux, but we can conclude that the observed distribution of asteroid families does not contradict with a cometary LHB.

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.

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.

This paper explores the possibility that the progenitors of the small satellites of Pluto got captured in the Pluto-Charon system from the massive heliocentric planetesimal disk in which Pluto was originally embedded into. We find that, if the dynamical excitation of the disk is small, temporary capture in the Pluto-Charon system can occur with non-negligible probability, due to the dynamical perturbations exerted by the binary nature of the Pluto-Charon pair. However, the captured objects remain on very elliptic orbits and the typical capture time is only 100 years. In order to explain the origin of the small satellites of Pluto, we conjecture that some of these objects got disrupted during their Pluto-bound phase by a collision with a planetesimal of the disk. This could have generated a debris disk, which damped under internal collisional evolution, until turning itself into an accretional disk that could form small satellites on circular orbits, co-planar with Charon. Unfortunately, we find that objects large enough to carry a sufficient amount of mass to generate the small satellites of Pluto have collisional lifetimes orders of magnitude longer than the capture time. Thus, this scenario cannot explain the origin of the small satellites of Pluto, which remains elusive.

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.

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.

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.

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.

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.

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.

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.

Observations indicate that the gaseous circumstellar disks around young stars vary significantly in size, ranging from tens to thousands of AU. Models of planet formation depend critically upon the properties of these primordial disks, yet in general it is impossible to connect an existing planetary system with an observed disk. We present a method by which we can constrain the size of our own protosolar nebula using the properties of the small body reservoirs in the solar system. In standard planet formation theory, after Jupiter and Saturn formed they scattered a significant number of remnant planetesimals into highly eccentric orbits. In this paper, we show that if there had been a massive, extended protoplanetary disk at that time, then the disk would have excited Kozai oscillations in some of the scattered objects, driving them into high-inclination (i > 50 deg), low-eccentricity orbits (q > 30 AU). The dissipation of the gaseous disk would strand a subset of objects in these high-inclination orbits; orbits that are stable on Gyr time scales. To date, surveys have not detected any Kuiper Belt Objects with orbits consistent with this dynamical mechanism. Using these non-detections by the Deep Ecliptic Survey (DES) and the Palomar Distant Solar System Survey we are able to rule out an extended gaseous protoplanetary disk (R_D > 80 AU) in our solar system at the time of Jupiter's formation. Future deep all sky surveys such as the Large Synoptic Survey Telescope (LSST) will all us to further constrain the size of the protoplanetary disk.

Modeling the formation of the ice giants Uranus and Neptune is a long-lasting problem in planetary science. Due to gas-drag, collisional damping, and resonant shepherding, the planetary embryos repel the planetesimals away from their reach and thus they stop growing (Levison et al. 2010). This problem persists independently of whether the accretion took place at the current locations of the ice giants or closer to the Sun. Instead of trying to push the runaway/oligarchic growth of planetary embryos up to 10-15 Earth masses, we envision the possibility that the planetesimal disk could generate a system of planetary embryos of only 1-3 Earth masses. Then we investigate whether these embryos could have collided with each other and grown enough to reach the masses of current Uranus and Neptune. Our results point to two major problems. First, there is typically a large difference in mass between the first and the second most massive core formed and retained beyond Saturn. Second, in many simulations the final planetary system has more than two objects beyond Saturn. The growth of a major planet from a system of embryos requires strong damping of eccentricities and inclinations from the disk of gas. But strong damping also favors embryos and cores to find a stable resonant configuration, so that systems with more than two surviving objects are found. In addition to these problems, in order to have substantial mutual accretion among embryos, it is necessary to assume that the surface density of the gas was several times higher than that of the minimum-mass solar nebula. However this contrasts with the common idea that Uranus and Neptune formed in a gas-starving disk, which is suggested by the relatively small amount of hydrogen and helium contained in the atmospheres of these planets. Only one of our simulations "by chance" successfully reproduced the structure of the outer Solar System.

Jupiter and Saturn formed in a few million years (Haisch et al. 2001) from a gas-dominated protoplanetary disk, and were susceptible to gas-driven migration of their orbits on timescales of only ~100,000 years (Armitage 2007). Hydrodynamic simulations show that these giant planets can undergo a two-stage, inward-then-outward, migration (Masset & Snellgrove 2001, Morbidelli & Crida 2007, Pierens & Nelson 2008). The terrestrial planets finished accreting much later (Klein et al. 2009), and their characteristics, including Mars' small mass, are best reproduced by starting from a planetesimal disk with an outer edge at about one astronomical unit from the Sun (Wetherill 1978, Hansen 2009) (1 AU is the Earth-Sun distance). Here we report simulations of the early Solar System that show how the inward migration of Jupiter to 1.5 AU, and its subsequent outward migration, lead to a planetesimal disk truncated at 1 AU; the terrestrial planets then form from this disk over the next 30-50 million years, with an Earth/Mars mass ratio consistent with observations. Scattering by Jupiter initially empties but then repopulates the asteroid belt, with inner-belt bodies originating between 1 and 3 AU and outer-belt bodies originating between and beyond the giant planets. This explains the significant compositional differences across the asteroid belt. The key aspect missing from previous models of terrestrial planet formation is the substantial radial migration of the giant planets, which suggests that their behaviour is more similar to that inferred for extrasolar planets than previously thought.

The European Space Agency's Rosetta spacecraft passed by the main belt asteroid (21) Lutetia the 10th July 2010. With its ~100km size, Lutetia is one of the largest asteroids ever imaged by a spacecraft. During the flyby, the on-board OSIRIS imaging system acquired spectacular images of Lutetia's northern hemisphere revealing a complex surface scarred by numerous impact craters, reaching the maximum dimension of about 55km. In this paper, we assess the cratering history of the asteroid. For this purpose, we apply current models describing the formation and evolution of main belt asteroids, that provide the rate and velocity distributions of impactors. These models, coupled with appropriate crater scaling laws, allow us to interpret the observed crater size-frequency distribution (SFD) and constrain the cratering history. Thanks to this approach, we derive the crater retention age of several regions on Lutetia, namely the time lapsed since their formation or global surface reset. We also investigate the influence of various factors -like Lutetia's bulk structure and crater obliteration- on the observed crater SFDs and the estimated surface ages. From our analysis, it emerges that Lutetia underwent a complex collisional evolution, involving major local resurfacing events till recent times. The difference in crater density between the youngest and oldest recognized units implies a difference in age of more than a factor of 10. The youngest unit (Beatica) has an estimated age of tens to hundreds of Myr, while the oldest one (Achaia) formed during a period when the bombardment of asteroids was more intense than the current one, presumably around 3.6Gyr ago or older.

We model the long-term evolution of the Hilda collisional family located in the 3/2 mean-motion resonance with Jupiter. Its eccentricity distribution evolves mostly due to the Yarkovsky/YORP effect and assuming that: (i) impact disruption was isotropic, and (ii) albedo distribution of small asteroids is the same as for large ones, we can estimate the age of the Hilda family to be $4_{-1}^{+0}\,{\rm Gyr}$. We also calculate collisional activity in the J3/2 region. Our results indicate that current collisional rates are very low for a 200\,km parent body such that the number of expected events over Gyrs is much smaller than one. The large age and the low probability of the collisional disruption lead us to the conclusion that the Hilda family might have been created during the Late Heavy Bombardment when the collisions were much more frequent. The Hilda family may thus serve as a test of orbital behavior of planets during the LHB. We tested the influence of the giant-planet migration on the distribution of the family members. The scenarios that are consistent with the observed Hilda family are those with fast migration time scales $\simeq 0.3\,{\rm Myr}$ to $3\,{\rm Myr}$, because longer time scales produce a family that is depleted and too much spread in eccentricity. Moreover, there is an indication that Jupiter and Saturn were no longer in a compact configuration (with period ratio $P_{\rm S}/P_{\rm J} > 2.09$) at the time when the Hilda family was created.

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....

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.

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.

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.

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.

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.

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.

We present the first far-IR observations of the solar-type stars delta Pav, HR 8501, 51 Peg and zeta^2 Ret, taken within the context of the DUNES Herschel Open Time Key Programme (OTKP). This project uses the PACS and SPIRE instruments with the objective of studying infrared excesses due to exo-Kuiper belts around nearby solar-type stars. The observed 100 um fluxes from delta Pav, HR 8501, and 51 Peg agree with the predicted photospheric fluxes, excluding debris disks brighter than Ldust/Lstar ~ 5 x 10^-7 (1 sigma level) around those stars. A flattened, disk-like structure with a semi-major axis of ~ 100 AU in size is detected around zeta^2 Ret. The resolved structure suggests the presence of an eccentric dust ring, which we interpret as an exo-Kuiper belt with Ldust/Lstar ~ 10^-5.

The habitability of planets is strongly affected by impacts from comets and asteroids. Indications from the ages of Moon rocks suggest that the inner Solar System experienced an increased rate of impacts roughly 3.8 Gya known as the Late Heavy Bombardment (LHB). Here we develop a model of how the Solar System would have appeared to a distant observer during its history based on the Nice model of Gomes et al. (2005). We compare our results with observed debris discs. We show that the Solar System would have been amongst the brightest of these systems before the LHB. Comparison with the statistics of debris disc evolution shows that such heavy bombardment events must be rare occurring around less than 12% of Sun-like stars.

Recent observations of exoplanets by direct imaging, reveal that giant planets orbit at a few dozens to more than a hundred of AU from their central star. The question of the origin of these planets challenges the standard theories of planet formation. We propose a new way of obtaining such far planets, by outward migration of a pair of planets formed in the 10 AU region. Two giant planets in mean motion resonance in a common gap in the protoplanetary disk migrate outwards, if the inner one is significantly more massive than the outer one. Using hydrodynamical simulations, we show that their semi major axes can increase by almost one order of magnitude. In a flared disk, the pair of planets should reach an asymptotic radius. This mechanism could account for the presence of Fomalhaut b ; then, a second, more massive planet, should be orbiting Fomalhaut at about 75 AU.

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.

It is known that near-Earth objects (NEOs) during their orbital evolution may often undergo close approaches to the Sun. Indeed it is estimated that up to ~70% of them end their orbital evolution colliding with the Sun. Starting from the present orbital properties, it is possible to compute the most likely past evolution for every NEO, and to trace its distance from the Sun. We find that a large fraction of the population may have experienced in the past frequent close approaches, and thus, as a consequence, a considerable Sun-driven heating, not trivially correlated to the present orbits. The detailed dynamical behaviour, the rotational and the thermal properties of NEOs determine the exact amount of the resulting heating due to the Sun. In the present paper we discuss the general features of the process, providing estimates of the surface temperature reached by NEOs during their evolution. Moreover, we investigate the effects of this process on meteor-size bodies, analyzing possible differences with the NEO population. We also discuss some possible effects of the heating which can be observed through remote sensing by ground-based surveys or space missions.

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.

We review the problem of the formation of terrestrial planets, with particular emphasis on the interaction of dynamical and geochemical models. The lifetime of gas around stars in the process of formation is limited to a few million years based on astronomical observations, while isotopic dating of meteorites and the Earth-Moon system suggest that perhaps 50-100 million years were required for the assembly of the Earth. Therefore, much of the growth of the terrestrial planets in our own system is presumed to have taken place under largely gas-free conditions, and the physics of terrestrial planet formation is dominated by gravitational interactions and collisions. The earliest phase of terrestrial-planet formation involve the growth of km-sized or larger planetesimals from dust grains, followed by the accumulations of these planetesimals into ~100 lunar- to Mars-mass bodies that are initially gravitationally isolated from one-another in a swarm of smaller planetesimals, but eventually grow to the point of significantly perturbing one-another. The mutual perturbations between the embryos, combined with gravitational stirring by Jupiter, lead to orbital crossings and collisions that drive the growth to Earth-sized planets on a timescale of 10-100 million years. Numerical treatment of this process has focussed on the use of symplectic integrators which can rapidly integrate the thousands of gravitationally-interacting bodies necessary to accurately model planetary growth. While the general nature of the terrestrial planets--their sizes and orbital parameters--seem to be broadly reproduced by the models, there are still some outstanding dynamical issues. One of these is the presence of an embryo-sized body, Mars, in our system in place of the more massive objects that simulations tend to yield. [Abridged]

The Nice model of Gomes et al. (2005) suggests that the migration of the giant planets caused a planetesimal clearing event which led to the Late Heavy Bombardment (LHB) at 880 Myr. Here we investigate the IR emission from the Kuiper belt during the history of the Solar System as described by the Nice model. We describe a method for easily converting the results of n-body planetesimal simulations into observational properties (assuming black-body grains and a single size distribution) and further modify this method to improve its realism (using realistic grain properties and a three-phase size distribution). We compare our results with observed debris discs and evaluate the plausibility of detecting an LHB-like process in extrasolar systems. Recent surveys have shown that 4% of stars exhibit 24 um excess and 16% exhibit 70 um excess. We show that the Solar System would have been amongst the brightest of these systems before the LHB at both 24 and 70 um. We find a significant increase in 24 um emission during the LHB, which rapidly drops off and becomes undetectable within 30 Myr, whereas the 70 um emission remains detectable until 360 Myr after the LHB. Comparison with the statistics of debris disc evolution shows that such depletion events must be rare occurring around less than 12% of Sun-like stars and with this level of incidence we would expect approximately 1 of the 413 Sun-like, field stars so far detected to have a 24 um excess to be currently going through an LHB. We also find that collisional processes are important in the Solar System before the LHB and that parameters for weak Kuiper belt objects are inconsistent with the Nice model interpretation of the LHB.

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)

The origin of Saturn\' s massive ring system is still unknown. Two popular scenarios - the tidal splitting of passing comets and the collisional destruction of a satellite - rely on a high cometary flux in the past. In the present paper we attempt to quantify the cometary flux during the Late Heavy Bombardment (LHB) to assess the likelihood of both scenarios. Our analysis relies on the so-called Nice model of the origin of the LHB (Tsiganis et al., 2005; Morbidelli et al., 2005; Gomes et al., 2005) and on the size distribution of the primordial trans-Neptunian planetesimals constrained in Charnoz & Morbidelli (2007). We find that the cometary flux on Saturn during the LHB was so high that both scenarios for the formation of Saturn rings are viable in principle. However, a more detailed study shows that the comet tidal disruption scenario implies that all four giant planets should have comparable ring systems whereas the destroyed satellite scenario would work only for Saturn, and perhaps Jupiter. This is because in Saturn\'s system, the synchronous orbit is interior to the Roche Limit, which is a necessary condition for maintaining a satellite in the Roche zone up to the time of the LHB. We also discuss the apparent elimination of silicates from the ring parent body implied by the purity of the ice in Saturn \' s rings. The LHB has also strong implications for the survival of the Saturnian satellites: all satellites smaller than Mimas would have been destroyed during the LHB, whereas Enceladus would have had from 40% to 70% chance of survival depending on the disruption model. In conclusion, these results suggest that the LHB is the sweet moment for the formation of a massive ring system around Saturn.

We study the conditions for collisions between planetesimals to be accretional or disruptive in turbulent disks, through analytical arguments based on fluid dynamical simulations and orbital integrations. In turbulent disks, the velocity dispersion of planetesimals is pumped up by random gravitational perturbations from density fluctuations of the disk gas. When the velocity dispersion is larger than the planetesimals' surface escape velocity, collisions between planetesimals do not result in accretion, and may even lead to their destruction. In disks with a surface density equal to that of the ``minimum mass solar nebula'' and with nominal MRI turbulence, we find that accretion proceeds only for planetesimals with sizes above $\sim 300$ km at 1AU and $\sim 1000$ km at 5AU. We find that accretion is facilitated in disks with smaller masses. However, at 5AU and for nominal turbulence strength, km-sized planetesimals are in a highly erosive regime even for a disk mass as small as a fraction of the mass of Jupiter. The existence of giant planets implies that either turbulence was weaker than calculated by standard MRI models or some mechanism was capable of producing Ceres-mass planetesimals in very short timescales. In any case, our results show that in the presence of turbulence planetesimal accretion is most difficult in massive disks and at large orbital distances.

We explore the origin and orbital evolution of the Kuiper belt in the framework of a recent model of the dynamical evolution of the giant planets, sometimes known as the Nice model. This model is characterized by a short, but violent, instability phase, during which the planets were on large eccentricity orbits. One characteristic of this model is that the proto-planetary disk must have been truncated at roughly 30 to 35 AU so that Neptune would stop migrating at its currently observed location. As a result, the Kuiper belt would have initially been empty. In this paper we present a new dynamical mechanism which can deliver objects from the region interior to ~35 AU to the Kuiper belt without excessive inclination excitation. Assuming that the last encounter with Uranus delivered Neptune onto a low-inclination orbit with a semi-major axis of ~27 AU and an eccentricity of ~0.3, and that subsequently Neptune's eccentricity damped in ~1 My, our simulations reproduce the main observed properties of the Kuiper belt at an unprecedented level.

A well-known bottleneck for the core-accretion model of giant-planet formation is the loss of the cores into the star by Type-I migration, due to the tidal interactions with the gas disk. It has been shown that a steep surface-density gradient in the disk, such as the one expected at the boundary between an active and a dead zone, acts as a planet trap and prevents isolated cores from migrating down to the central star. We study the relevance of the planet trap concept for the accretion and evolution of systems of multiple planetary embryos/cores. We performed hydrodynamical simulations of the evolution of systems of multiple massive objects in the vicinity of a planet trap. The planetary embryos evolve in 3 dimensions, whereas the disk is modeled with a 2D grid. Synthetic forces are applied onto the embryos to mimic the damping effect that the disk has on their inclinations. Systems with two embryos tend to acquire stable, separated and non-migrating orbits, with the more massive embryo placed at the planet trap and the lighter one farther out in the disk. Systems of multiple embryos are intrinsically unstable. Consequently, a long phase of mutual scattering can lead to accreting collisions among embryos; some embryos are injected into the inner part of the disk, where they can be evacuated into the star by Type I migration. The system can resume a stable, non-migrating configuration only when the number of surviving embryos decreases to a small value (~2-4). This can explain the limited number of giant planets in our solar system. These results should apply in general to any case in which the Type-I migration of the inner embryo is prevented by some mechanism, and not solely to the planet trap scenario.

We study the orbital evolution of the 4 giant planets of our solar system in a gas disk. Our investigation extends the previous works by Masset and Snellgrove (2001) and Morbidelli and Crida (2007, MC07), which focussed on the dynamics of the Jupiter-Saturn system. The only systems that we found to reach a steady state are those in which the planets are locked in a quadruple mean motion resonance (i.e. each planet is in resonance with its neighbor). In total we found 6 such configurations. For the gas disk parameters found in MC07, these configurations are characterized by a negligible migration rate. After the disappearance of the gas, and in absence of planetesimals, only two of these six configurations (the least compact ones) are stable for a time of hundreds of millions of years or more. The others become unstable on a timescale of a few My. Our preliminary simulations show that, when a planetesimal disk is added beyond the orbit of the outermost planet, the planets can evolve from the most stable of these configurations to their current orbits in a fashion qualitatively similar to that described in Tsiganis et al. (2005).

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.

This chapter discusses the dynamical properties of the Kuiper belt population. Then, it focuses on the characteristics of the Kuiper belt that cannot be explained by its evolution in the framework of the current solar system. We review models of primordial solar system evolution that have been proposed to reproduce the Kuiper belt features, outlining advantages and problems of each of them.

We study the effect of a Jovian planet on the gas distribution of a protoplanetary disc, using a new numerical scheme that allows us to take into consideration the global evolution of the disc, down to an arbitrarily small inner physical radius. We find that Jovian planets do not open cavities in the inner part of the disc (i.e. interior to their orbits) unless (a) the inner physical edge of the disc is close to the planet's location or (b) the planet is much more massive than the disc. In all other cases the planet simply opens a gap in the gas density distribution, whose global profile is essentially unchanged relative to the one that it would have if the planet were absent. We recognize, though, that the dust distribution can be significantly different from the gas distribution and that dust cavities might be opened in some situations, even if the gas is still present in the inner part of the disc. Concerning the migration of the planet, we find that classical Type-II migration (with speed proportional to the viscosity of the disc) occurs only if the gap opened by the planet is deep and clean. If there is still a significant amount of gas in the gap, the migration of the planet is generally slower than the theoretical Type-II migration rate. In some situations, migration can be stopped or even reversed. We develop a simple model that reproduces satisfactorily the migration rate observed in the simulations, for a wide range of disc viscosities and planet masses and locations relative to the inner disc edge. Our results are relevant for extra-solar planetary systems, as they explain (a) why some hot Jupiters did not migrate all the way down to their parent stars and (b) why the outermost of a pair of resonant planets is typically the most massive one.

Chiang et al. 2006, hereafter C06 have recently proposed that the observed structure of the Kuiper belt could be the result of a dynamical instability of a system of ~5 primordial ice giant planets in the outer Solar System. According to this scenario, before the instability occurred, these giants were growing in a highly collisionally damped environment according to the arguments in Goldreich et al. (2004a,b, hereafter G04). Here we test this hypothesis with a series of numerical simulations using a new code designed to incorporate the dynamical effects of collisions. We find that we cannot reproduce the observed Solar System. In particular, G04 and C06 argue that during the instability, all but two of the ice giants would be ejected from the Solar System by Jupiter and Saturn, leaving Uranus and Neptune behind. We find that ejections are actually rare and that instead the systems spread outward. This always leads to a configuration with too many planets that are too far from the Sun. Thus, we conclude that both G04's scheme for the formation of Uranus and Neptune and C06's Kuiper belt formation scenario are not viable in their current forms.

We have investigated the final accretion stage of terrestrial planets from Mars-mass protoplanets that formed through oligarchic growth in a disk comparable to the minimum mass solar nebula (MMSN), through N-body simulation including random torques exerted by disk turbulence due to Magneto-Rotational-Instability. For the torques, we used the semi-analytical formula developed by Laughlin et al.(2004). The damping of orbital eccentricities (in all runs) and type-I migration (in some runs) due to the tidal interactions with disk gas are also included. We found that the orbital eccentricities pumped up by the turbulent torques and associated random walks in semimajor axes tend to delay isolation of planets, resulting in more coagulation of planets than in the case without turbulence. The eccentricities are still damped after planets become isolated. As a result, the number of final planets decreases with increase in strength of the turbulence, while Earth-mass planets with small eccentricities are still formed. In the case of relatively strong turbulence, the number of final planets are 4-5 at 0.5-2AU, which is consistent with Solar system, for relatively wide range of disk surface density (~10^{-4}-10^{-2} times MMSN).

Numerical simulations of planet-disk interactions are usually performed with hydro-codes that -- because they consider only an annulus of the disk, over a 2D grid -- can not take into account the global evolution of the disk. However, the latter governs planetary migration of type II, so that the accuracy of the planetary evolution can be questioned. To develop an algorithm that models the local planet-disk interactions together with the global viscous evolution of the disk, we surround the usual 2D grid with a 1D grid ranging over the real extension of the disk. The 1D and 2D grids are coupled at their common boundaries via ghost rings, paying particular attention to the fluxes at the interface, especially the flux of angular momentum carried by waves. The computation is done in the frame centered on the center of mass to ensure angular momentum conservation. The global evolution of the disk and the local planet-disk interactions are both well described and the feedback of one on the other can be studied with this algorithm, for a negligible additional computing cost with respect to usual algorithms.

We estimate the total number and the slope of the size frequency distribution (SFD) of dormant Jupiter Family Comets (JFCs) by fitting a one-parameter model to the known population. We first select 61 Near Earth Objects (NEOs) that are likely to be dormant JFCs because their orbits are dynamically coupled to Jupiter (Bottke et al, 2002). Then, from the numerical simulations of Levison & Duncan (1997), we construct an orbit distribution model for JFCs in the NEO orbital element space. We assume an orbit independent SFD for all JFCs, the slope of which is our unique free parameter. Finally, we compute observational biases for dormant JFCs using a calibrated NEO survey simulator (Jedicke et al. 2003). By fitting the biased model to the data, we estimate that there are ~75 dormant JFCs with H<18 in the NEO region and that the slope of their cumulative SFD is -1.5 +/- 0.3. Our slope for the SFD of dormant JFCs is very close to that of active JFCs as determined by Weissman and Lowry (2003). Thus, we argue that when JFCs fade they are likely to become dormant rather than to disrupt and that the fate of faded comets is size independent. Our results imply that the size distribution of the JFC progenitors - the scattered disk trans-Neptunian population - either (i) has a similar and shallow SFD or (i') is slightly steeper and physical processes acting on the comets in a size-dependent manner creates the shallower active comet SFD. Our measured slope, typical of collisionally evolved populations with a size dependent impact strength (Benz and Asphaug 1999), suggests that scattered disk bodies reached collisional equilibrium inside the proto-planetary disk prior to their removal from the planetary region.

[Abridged] The tidal torque exerted by a protoplanetary disk with power law surface density and temperature profiles onto an embedded protoplanetary embryo is generally a negative quantity that leads to the embryo inwards migration. Here we investigate how the tidal torque balance is affected at a disk surface density radial jump. The jump has two consequences : - it affects the differential Lindblad torque. In particular if the disk is merely empty on the inner side, the differential Lindblad torque almost amounts to the large negative outer Lindblad torque. - It affects the corotation torque, which is a quantity very sensitive to the local gradient of the disk surface density. In particular if the disk is depleted on the inside and if the jump occurs radially over a few pressure scale-heights, the corotation torque is a positive quantity that is much larger than in a power-law disk. We show by means of customized numerical simulations of low mass planets embedded in protoplanetary nebulae with a surface density jump that the second effect is dominant, that is that the corotation torque largely dominates the differential Lindblad torque on the edge of a central depletion. As a consequence the type I migration of low mass objects reaching the jump should be halted, and all these objects should be trapped there.

This text was originally written to accompany a series of lectures that I gave at the `35th Saas-Fee advanced course' in Switzerland and at the Institute for Astronomy of the University of Hawaii. It reviews my current understanding of the dynamics of comets and of the origin and primordial sculpting of their reservoirs. It starts discussing the structure of the Kuiper belt and the current dynamics of Kuiper belt objects, including scattered disk objects. Then it discusses the dynamical evolution of Jupiter family comets from the trans-Neptunian region, and of long period comets from the Oort cloud. The formation of the Oort cloud is then reviewed, as well as the primordial sculpting of the Kuiper belt. Finally, these issues are revisited in the light of a new model of giant planets evolution that has been developed to explain the origin of the late heavy bombardment of the terrestrial planets.

Although it is well known that a massive planet opens a gap in a proto-planetary gaseous disk, there is no analytic description of the surface density profile in and near the gap. The simplest approach, which is based upon the balance between the torques due to the viscosity and the gravity of the planet and assumes local damping, leads to gaps with overestimated width, especially at low viscosity. Here, we take into account the fraction of the gravity torque that is evacuated by pressure supported waves. With a novel approach, which consists of following the fluid elements along their trajectories, we show that the flux of angular momentum carried by the waves corresponds to a pressure torque. The equilibrium profile of the disk is then set by the balance between gravity, viscous and pressure torques. We check that this balance is satisfied in numerical simulations, with a planet on a fixed circular orbit. We then use a reference numerical simulation to get an ansatz for the pressure torque, that yields gap profiles for any value of the disk viscosity, pressure scale height and planet to primary mass ratio. Those are in good agreement with profiles obtained in numerical simulations over a wide range of parameters. Finally, we provide a gap opening criterion that simultaneously involves the planet mass, the disk viscosity and the aspect ratio.

The Hannay angle has been previously studied for a celestial circular restricted three-body system by means of an adiabatic approach. In the present work, three main results are obtained. Firstly, a formal connection between perturbation theory and the Hamiltonian adiabatic approach shows that both lead to the Hannay angle; it is thus emphasised that this effect is already contained in classical celestial mechanics, although not yet defined nor evaluated separately. Secondly, a more general expression of the Hannay angle, valid for an action-dependent potential is given; such a generalised expression takes into account that the restricted three-body problem is a time-dependent, two degrees of freedom problem even when restricted to the circular motion of the test body. Consequently, (some of) the eccentricity terms cannot be neglected {\it a priori}. Thirdly, we present a new numerical estimate for the Earth adiabatically driven by Jupiter. We also point out errors in a previous derivation of the Hannay angle for the circular restricted three-body problem, with an action-independent potential.

Explaining the origin of the orbit of 2000 CR105 (a ~ 230AU, q ~ 45AU) is a major test for our understanding of the primordial evolution of the outer Solar System. Gladman et al. (2001) showed that this objects could not have been a normal member of the scattered disk that had its perihelion distance increased by chaotic diffusion. In this paper we explore four seemingly promising mechanisms for explaining the origin of the orbit of this peculiar object: (i) the passage of Neptune through a high-eccentricity phase, (ii) the past existence of massive planetary embryos in the Kuiper belt or the scattered disk, (iii) the presence of a massive trans-Neptunian disk at early epochs which exerted tides on scattered disk objects, and (iv) encounters with other stars. Of all these mechanisms, the only one giving satisfactory results is the passage of a star. Indeed, our simulations show that the passage of a solar mass star at about 800 AU only perturbs objects with semi-major axes larger than roughly 200 AU to large perihelion distances. This is in good agreement with the fact that 2000 CR105 has a semi-major axis of 230AU and no other bodies with similar perihelion distances but smaller semi-major axes have yet been discovered. The discovery of 2003 VB12, (a=450AU, q=75AU) announced a few days before the submission of this paper, strengthen our conclusions.