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.

During planet formation gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration prescriptions proposed thus far describing the planet-disc interaction force and the timescales of orbital change in the supersonic regime are inconsistent with one another. Here we discuss the details of existing prescriptions in the literature and derive a new simple and intuitive formulation for planet-disc interactions based on dynamical friction that can be applied in both supersonic and subsonic cases. While the existing prescriptions assume particular disc models, ours include the explicit dependence on the disc parameters; hence it can be applied to discs with any radial surface density and temperature dependence (except for the local variations with radial scales less than the disc scale height). Our prescription will reduce the uncertainty originating from different literature formulations of planet migration and will be an important tool to study planet accretion processes, especially when studying the formation of close-in low-mass planets that are commonly found in exoplanetary systems.

Subsequent to the Moon's formation, late accretion to the terrestrial planets modified their silicate crusts and mantles. We combine dynamical N-body and Monte Carlo simulations to determine impact probabilities, impact velocities, and expected mass augmentation onto the terrestrial planets from three sources: planetesimals left over from primary accretion, asteroids from the E-belt, and comets arriving from the outer Solar System. We present estimates of the amount of cometary material striking the terrestrial planets in an early episode of planetesimal-driven giant planet migration (Mojzsis et al., 2019). The Moon and Mars suffer proportionally higher cometary accretion than Venus and Earth. The background mass addition from small leftover planetesimals to Earth and Mars is far less than independent estimates based on their mantle abundances of highly-siderophile elements. This supports the theory that both planets were struck by single large bodies that delivered most of their mass augmentation since primary accretion. We present fits to the impact chronologies valid from 4500 Ma to ca. 3700 Ma. The lunar timeline thus obtained is at odds with both the calibrated Neukum (Neukum et al., 2001) and Werner (Werner et al., 2014; Werner, 2019) chronologies. For Mars, the match with its calibrated Werner chronology is no better; by increasing the mass of the E-belt by a factor of four the dynamical lunar and martian chronologies are in line with that of Werner (2019) and match constraints from the current population of Hungaria asteroids. Yet, neither of our dynamical timelines fit well with that of Neukum. The dynamical lunar and martian chronologies are also different from each other. Consequently, the usual extrapolation of such chronologies from one planetary body to the other is technically inappropriate.

The widespread presence of ribonucleic acid (RNA) catalysts and cofactors in Earth's biosphere today suggests that RNA was the first biopolymer to support Darwinian evolution. However, most "path-hypotheses" to generate building blocks for RNA require reduced nitrogen-containing compounds not made in useful amounts in the CO2-N2-H2O atmospheres of the Hadean. We review models for Earth's impact history that invoke a single ~10^23 kg impactor (Moneta) to account for measured amounts of platinum, gold, and other siderophilic ("iron-loving") elements on the Earth and Moon. If it were the last sterilizing impactor, Moneta would have reduced the atmosphere but not its mantle, opening a "window of opportunity" for RNA synthesis, a period when RNA precursors rained from the atmosphere to land holding oxidized minerals that stabilize advanced RNA precursors and RNA. Surprisingly, this combination of physics, geology, and chemistry suggests a time when RNA formation was most probable, ~120 +/- 100 million years after Moneta's impact, or ~4.36 +/- 0.1 billion years ago. Uncertainties in this time are driven by uncertainties in rates of productive atmosphere loss and amounts of sub-aerial land.

The abundance of highly siderophile elements (HSEs) inferred for Mars' mantle from martian meteorites implies a Late Veneer (LV) mass addition of ~0.8 wt% with broadly chondritic composition. Late accretion to Mars by a differentiated Ceres-sized (~1000 km diameter) object can account for part of the requisite LV mass, and geochronological constraints suggests that this must have occurred no later than ca. 4480 Ma. Here, we analyze the outcome of the hypothetical LV giant impact to Mars with smoothed particle hydrodynamics simulations together with analytical theory. Results show that, in general about 50% of the impactor's metallic core shatters into ~10m fragments that subsequently fragment into sub-mm metallic hail at re-accretion. This returns a promising delivery of HSEs into martian mantle compared to either a head-on and hit-and-run collision; in both cases,<10% of impactor's core materials are fragmented and finally embedded in the martian mantle. Isotopic evidence from martian meteorites, and interpretations from atmospheric mapping data show that a global surface water reservoir could be present during the early Noachian (before ca. 4100 Ma). The millimeter-sized metal hail could thus react with a martian hydrosphere to generate ~3 bars of H2, which is adequate to act as a greenhouse and keep early Mars warm. Yet, we also find that this atmosphere is transient. It typically survives shorter than 3 Myr based on the expected extreme ultraviolet (EUV) flux of the early Sun; if the Sun was a slow rotator an accordingly weaker EUV flux could extend this lifetime to >10 Myr. A dense pre-Noachian CO2 atmosphere should lower the escape efficiency of hydrogen by IR emission. A more detailed hydrodynamic atmospheric model of this early hydrogen atmosphere is warranted to examine its effect on pre-Noachian Mars.

The TRAPPIST-1 planetary system consists of seven planets within 0.05 au of each other, five of which are in a multi-resonant chain. {These resonances suggest the system formed via planet migration; subsequent tidal evolution has damped away most of the initial eccentricities. We used dynamical N-body simulations to estimate how long it takes for the multi-resonant configuration that arises during planet formation to break. From there we use secular theory to pose limits on the tidal parameters of planets b and c. We calibrate our results against multi-layered interior models constructed to fit the masses and radii of the planets, from which the tidal parameters are computed independently.} The dynamical simulations show that the planets typically go unstable 30 Myr after their formation. {Assuming synchronous rotation throughout} we compute $\frac{k_2}{Q} \gtrsim 2\times 10^{-4}$ for planet b and $\frac{k_2}{Q} \gtrsim 10^{-3}$ for planet c. Interior models yield $(0.075-0.37) \times 10^{-4}$ for TRAPPIST-1 b and $(0.4-2)\times 10^{-4}$ for TRAPPIST-1 c. The agreement between the {dynamical and interior} models is not too strong, but is still useful to constrain the dynamical history of the system. We suggest that this two-pronged approach could be of further use in other multi-resonant systems if the planet's orbital and interior parameters are sufficiently well known.

The source of impactors capable of producing the cratering asymmetry on Triton is open to debate. While hyperbolic comets are the likely source of impactors, the theoretically predicted crater distribution due to impacts from comets is incompatible with the observed crater distribution on Triton. In this study, we adopted a dynamical approach to investigate the crater distribution produced by three sources of impactors: Neptune-centric impactors on prograde, low-inclination orbits, Neptune-centric impactors on high-inclination orbits such as its irregular satellites, and heliocentric impactors with hyperbolic orbits. We employ N-body simulations which take into account the dynamics of the impacts. We found that none of the impactor sources alone are able to reconcile with the crater distribution on Triton. The crater distribution produced by a low-inclination Neptune-centric source is narrowly concentrated on Triton's leading hemisphere whereas the heliocentric impactors produce too many craters on the trailing hemisphere compared to the observed crater distribution. All three impactor populations produce craters on the trailing hemisphere of Triton. We derived the impact distribution of a mixed population of impactors using a Monte Carlo approach and found that if the population consists of mainly heliocentric impactors (70%), the resulting impact distribution fits the observed crater distribution on Triton's leading hemisphere well.

Immediately after their formation, the terrestrial planets experienced intense impact bombardment by comets, leftover planetesimals from primary accretion, and asteroids. This temporal interval in solar system evolution, termed late accretion, thermally and chemically modified solid planetary surfaces and may have impeded the emergence of life on the Hadean Earth. The sources and tempo of late accretion are, however, vague. Here, we present a timeline that relates variably retentive radiometric ages from asteroidal meteorites, to new dynamical models of late accretion that invokes giant planet migration. Reconciliation of the geochronological data with dynamical models shows that giant planet migration immediately leads to an intense 30 Myr influx of comets to the entire solar system. The absence of whole-sale crustal reset ages after 4450 Ma for the most resilient chronometers from Earth, Moon, Mars, Vesta and various meteorite parent bodies confines the onset of giant planet migration to no later than ca. 4480 Ma. Waning impacts from planetesimals, asteroids (and a minor cometary component) continue to strike the inner planets through a protracted monotonic decline in impactor flux; this is in agreement with predictions from crater chronology. Amended global 3-D thermal analytical bombardment models derived from our new impact mass-production functions show that persistent niches for prebiotic chemistry on the early Hadean Earth could endure late accretion for at least the last 4400 Myr.

Planetary migration is a major challenge for planet formation theories. The speed of Type I migration is proportional to the mass of a protoplanet, while the final decade of growth of a pebble-accreting planetary core takes place at a rate that scales with the mass to the two-thirds power. This results in planetary growth tracks (i.e., the evolution of a protoplanet's mass versus its distance from the star) that become increasingly horizontal (migration-dominated) with rising mass of the protoplanet. It has been shown recently that the migration torque on a protoplanet is reduced proportional to the relative height of the gas gap carved by the growing planet. Here we show from 1-D simulations of planet-disc interaction that the mass at which a planet carves a 50% gap is approximately 2.3 times the pebble isolation mass. Our measurements of the pebble isolation mass from 1-D simulations match published 3-D results relatively well, except at very low viscosities where the 3-D pebble isolation mass is significantly higher, possibly due to gap edge instabilities not captured in 1-D. The pebble isolation mass demarks the transition from pebble accretion to gas accretion. Gas accretion to form gas-giant planets therefore takes place over a few astronomical units of migration after reaching first the pebble isolation mass and, shortly after, the 50% gap mass. Our results demonstrate how planetary growth can outperform migration, both during core accretion and during gas accretion, even when the Stokes number of the pebbles is small, St~0.01, and the pebble-to-gas flux ratio in the protoplanetary disc is in the nominal range of 0.01-0.02. We find that planetary growth is very rapid in the first million years of the protoplanetary disc and that the probability for forming gas-giant planets increases with the initial size of the protoplanetary disc and with decreasing turbulent diffusion.

Exoplanet interior modelling usually makes the assumption that the elemental abundances of a planet are identical to those of its host star. Host stellar abundances are good proxies of planetary abundances, but only for refractory elements. This is particularly true for terrestrial planets, as evidenced by the relative differences in bulk chemical composition between the Sun and the Earth and other inner solar system bodies. The elemental abundances of a planet host star must therefore be devolatilised in order to correctly represent the bulk chemical composition of its terrestrial planets. Furthermore, nickel and light elements make an important contribution alongside iron to the core of terrestrial planets. We therefore adopt an extended chemical network of the core, constrained by an Fe/Ni ratio of 18 $\pm$ 4 (by number). By applying these constraints to the Sun, our modelling reproduces the composition of the mantle and core, as well as the core mass fraction of the Earth. We also apply our modelling to four exoplanet host stars with precisely measured elemental abundances: Kepler-10, Kepler-20, Kepler-21 and Kepler-100. If these stars would also host terrestrial planets in their habitable zone, we find that such planets orbiting Kepler-21 would be the most Earth-like, while those orbiting Kepler-10 would be the least. To assess the similarity of a rocky exoplanet to the Earth in terms of interior composition and structure, high-precision host stellar abundances are critical. Our modelling implies that abundance uncertainties should be better than $\sim$ 0.04 dex for such an assessment to be made.

The intensity and effects of early impact bombardment on the major satellites of the giant planets during an episode of giant planet migration is still poorly known. We use a combination of dynamical N-body and Monte Carlo simulations to determine impact probabilities, impact velocities, and expected masses that collide with these satellites to determine the chronology of impacts during the migration. Volatile loss through bombardment is typically 20% for Miranda, a few percents for the larger Uranian satellites and negligible for the Galilean satellites. Due to its small size and the high impact velocity there is a >99% chance that Miranda suffered a catastrophic impact that shattered the satellite. Subsequent re-accretion from a circum-Uranian ring could account for its peculiar surface morphology and low density. The probability to destroy Ariel and Umbriel is 15% and 1% for Titania and Oberon. Approximately 90% of the mass in planetesimals that passes through the Jovian and Uranian satellite systems (about $4 {\rm \ M_{\oplus}}$ and $2 {\rm \ M_{\oplus}}$ respectively) does so in about 15 Myr. This extremely rapid and intense bombardment causes repeated local crustal melting on all satellites. The combination of these effects results in an entirely different impact chronology than that of the inner solar system. We conclude that the simple extrapolation of the lunar chronology to the outer solar system satellites is not correct. The tail end (after 25 Myr) of the chronology function has an e-folding time of 100 Myr at Jupiter, but follows a cumulative Weibull distribution at Uranus, making direct comparisons between the gas and ice giant planets difficult. Based on our results the surfaces of the Uranian satellites, Callisto, and possibly Ganymede, are all about the same age, and are roughly 150 Myr younger than the timing of the dynamical instability.

The first detected interstellar object 'Oumuamua that passed within 0.25au of the Sun on 2017 September 9 was presumably ejected from a stellar system. We use its newly determined non-Keplerian trajectory together with the reconstructed Galactic orbits of 7 million stars from Gaia DR2 to identify past close encounters. Such an "encounter" could reveal the home system from which 'Oumuamua was ejected. The closest encounter, at 0.60pc (0.53-0.67pc, 90% confidence interval), was with the M2.5 dwarf HIP 3757 at a relative velocity of 24.7km/s, 1Myr ago. A more distant encounter (1.6pc) but with a lower encounter (ejection) velocity of 10.7km/s was with the G5 dwarf HD 292249, 3.8Myr ago. Two more stars have encounter distances and velocities intermediate to these. The encounter parameters are similar across six different non-gravitational trajectories for 'Oumuamua. Ejection of 'Oumuamua by scattering from a giant planet in one of the systems is plausible, but requires a rather unlikely configuration to achieve the high velocities found. A binary star system is more likely to produce the observed velocities. None of the four home candidates have published exoplanets or are known to be binaries. Given that the 7 million stars in Gaia DR2 with 6D phase space information is just a small fraction of all stars for which we can eventually reconstruct orbits, it is a priori unlikely that our current search would find 'Oumuamua's home star system. As 'Oumuamua is expected to pass within 1pc of about 20 stars and brown dwarfs every Myr, the plausibility of a home system depends also on an appropriate (low) encounter velocity.

The formation and orbital evolution of Saturn's inner mid-sized moons are still debated. The most puzzling aspects are 1) how the Tethys-Dione pair and the Mimas-Enceladus pair passed through their strong 3:2 mean-motion resonances during the tidal orbital evolution, and 2) the current strong heat flow from Enceladus, which is a few orders of magnitude higher than the tidal energy dissipation caused by the present orbital eccentricity of Enceladus. Here we perform N-body simulations of the moons' orbital evolution from various initial conditions -- assuming that the moons were formed from Saturn's hypothetical massive ring -- and investigate possible paths to solve the above difficulties. If the influence of the rings is neglected, we find that the Tethys-Dione pair cannot avoid becoming trapped in the mean-motion resonances as they recede from Saturn, and that the Tethys-Enceladus pair cannot avoid collisions after the resonance trapping, in case Saturn's quality factor is smaller than 15,000. These findings are inconsistent with the current orbital configuration. However, taking into account both the eccentricity excitation and the orbital expansion caused by the ring torque, we find that these resonance captures are avoided. With the high eccentricity pumped up by the torque, Enceladus passes through the resonances with Tethys, and the Dione-Tethys pair passes through their 2:1 and possibly the 3:2 resonance as well. After Enceladus resides beyond the 2:1 resonance with the outer ring edge, the eccentricity can be tidally damped. While this is a promising path of evolution, in most runs, Enceladus collides with Tethys by the excited eccentricity. There is a hint that a ring mass decrease could avoid the collision. The tidal heat due to the eccentricity excitation by the ring torque is stored in the moons and slowly radiated away. It may account for the current high heat flow.

The formation of multiple close-in low-mass exoplanets is still a mystery. The challenge is to build a system wherein the outermost planet is beyond 0.2 AU from the star. Here we investigate how the prescription for type I planet migration affects the ability to trap multiple planets in a resonant chain near the inner edge of the protostellar disc. A sharp edge modelled as a hyperbolic tangent function coupled with supersonic corrections to the classical type I migration torques results in the innermost planets being pushed inside the cavity through resonant interaction with farther planets because migration is starward at slightly supersonic eccentricities. Planets below a few Earth masses are generally trapped in a resonant chain with the outermost planet near the disc edge, but long-term stability is not guaranteed. For more massive planets the migration is so fast that the eccentricity of the innermost resonant pair is excited to highly supersonic levels due to decreased damping on the innermost planet as it is pushed inside the cavity; collisions frequently occur and the system consists one or two intermediate-mass planets residing closer to the star than the disc's inner edge. We found a neat pileup of resonant planets outside the disc edge only if the corotation torque does not rapidly diminish at high eccentricity. We call for detailed studies on planet migration near the disc's inner edge, which is still uncertain, and for an improved understanding of eccentricity damping and disc torques in the supersonic regime.

Dynamical models of planet formation coupled with cosmochemical data from martian meteorites show that Mars' isotopic composition is distinct from that of Earth. Reconciliation of formation models with meteorite data require that Mars grew further from the Sun than its present position. Here, we evaluate this compositional difference in more detail by comparing output from two $N$-body planet formation models. The first of these planet formation models simulates what is termed the "Classical" case wherein Jupiter and Saturn are kept in their current orbits. We compare these results with another model based on the "Grand Tack", in which Jupiter and Saturn migrate through the primordial asteroid belt. Our estimate of the average fraction of chondrite assembled into Earth and Mars assumes that the initial solid disk consists of only sources of enstatite chondrite composition in the inner region, and ordinary chondrite in the outer region. Results of these analyses show that both models tend to yield Earth and Mars analogues whose accretion zones overlap. The Classical case fares better in forming Mars with its documented composition (29% to 68% enstatite chondrite plus 32% to 67% ordinary chondrite) though the Mars analogues are generally too massive. However, if we include the restriction of mass on the Mars analogues, the Classical model does not work better. We also further calculate the isotopic composition of $^{17} \rm O$, $^{50} \rm Ti$, $^{54} \rm Cr$, $^{142} \rm Nd$, $^{64} \rm Ni$, and $^{92} \rm Mo$ in the martian mantle from the Grand Tack simulations. We find that it is possible to match the calculated isotopic composition of all the above elements in Mars' mantle with their measured values, but the resulting uncertainties are too large to place good restriction on the early dynamical evolution and birth place of Mars.

Radiometric dating indicates that Mars accreted in the first ~4 Myr of solar system formation, which coincides with the formation and possible migration of Jupiter. While nebular gas from the protoplanetary disk was still present, Jupiter may have migrated inwards and tacked at 1.5 AU in a 3:2 resonance with Saturn. This migration excited planetary building blocks in the inner solar system, resulting in extensive mixing and planetesimal removal. Here we evaluate the plausible nature of Mars' building blocks, focusing in particular on how its growth was influenced by the formation and migration of Jupiter. We use a combination of dynamical simulations and an isotopic mixing model that traces the accretion of elements with different affinities for metal. Dynamical simulations show that Jupiter's migration causes the late stages of Earth's and Mars' accretion to be dominated by EC-type (enstatite chondrite) material due to the loss of OC (ordinary chondrite) planetesimals. Our analysis of available isotopic data for SNC meteorites shows that Mars consists of approximately 68%$^{+0}_{-39}$ EC + 32%$^{+35}_{-0}$ OC by mass (2{\sigma}). The large uncertainties indicate that isotopic analyses of martian samples are for the most part too imprecise to definitely test model predictions; in particular it remains uncertain whether or not Mars accreted predominantly EC material in the latter stages of its formation history. Dynamical simulations also provide no definitive constraint on Mars' accretion history due to the great variety of dynamical pathways that the martian embryo exhibits. The present work calls for new measurements of isotopic anomalies in SNC meteorites targeting siderophile elements (most notably Ni, Mo and Ru) to constrain Mars' accretion history and its formation location.

We report on the results of a systematic search for associated asteroid families for all active asteroids known to date. We find that 10 out of 12 main-belt comets (MBCs) and 5 out of 7 disrupted asteroids are linked with known or candidate families, rates that have ~0.1% and ~6% probabilities, respectively, of occurring by chance, given an overall family association rate of 37% for asteroids in the inner solar system. We find previously unidentified family associations between 238P/Read and the candidate Gorchakov family, 311P/PANSTARRS and the candidate Behrens family, 324P/La Sagra and the Alauda family, 354P/LINEAR and the Baptistina family, P/2013 R3-B (Catalina-PANSTARRS) and the Mandragora family, P/2015 X6 (PANSTARRS) and the Aeolia family, P/2016 G1 (PANSTARRS) and the Adeona family, and P/2016 J1-A/B (PANSTARRS) and the Theobalda family. All MBCs with family associations belong to families that contain asteroids with primitive taxonomic classifications and low average reported albedos (pV_avg < 0.10), while disrupted asteroids with family associations belong to families that contain asteroids that span wider ranges of taxonomic types and average reported albedos (0.06 < pV_avg < 0.25). These findings are consistent with MBC activity being closely correlated to composition (i.e., whether an object is likely to contain ice), while disrupted asteroid activity is not as sensitive to composition. Given our results, we describe a sequence of processes by which the formation of young asteroid families could lead to the production of present-day MBCs.

Impact bombardment on Mercury in the solar system late accretion phase (ca. 4.4 to 3.8 Ga) caused considerable mechanical, chemical and thermal reworking of its silicate reservoirs (crust and mantle). Depending on the frequency, size and velocity of such impactors, effects included regional and global scale crustal melting, and thermal perturbations of the mercurian mantle. We use a 3D transient heating model to test the effects of two bombardment scenarios on early (pre Tolstojan) Mercury mantle and crust. Results show that rare impacts by the largest (greater than 100 km diameter) bodies deliver sufficient heat to the shallow mercurian mantle producing high-temperature ultra-magnesian (komatiitic s.s.) melts. Impact heating leading to effusive (flood) volcanism can account for eponymous High-Magnesium Region (HMR) observed during the MErcury Surface, Space Environment, GEochemistry Ranging (MESSENGER) mission. We find that late accretion to Mercury induced volumetrically significant crustal melting (less than 58 vol.percent), mantle heating and melt production, which, combined with extensive resurfacing (less than 100 percent), also explains why its oldest cratering record was effectively erased, consistent with crater-counting statistics.

Overabundances in highly siderophile elements (HSEs) of Earth's mantle can be explained by conveyance from a singular, immense (3000 km in a diameter) "Late Veneer" impactor of chondritic composition, subsequent to lunar formation and terrestrial core-closure. Such rocky objects of approximately lunar mass (about 0.01 M_E) ought to be differentiated, such that nearly all of their HSE payload is sequestered into iron cores. Here, we analyze the mechanical and chemical fate of the core of such a Late Veneer impactor, and trace how its HSEs are suspended - and thus pollute - the mantle. For the statistically most-likely oblique collision (about 45degree), the impactor's core elongates and thereafter disintegrates into a metallic hail of small particles (about 10 m). Some strike the orbiting Moon as sesquinary impactors, but most re-accrete to Earth as secondaries with further fragmentation. We show that a single oblique impactor provides an adequate amount of HSEs to the primordial terrestrial silicate reservoirs via oxidation of (<m-sized) metal particles with a hydrous, pre-impact, early Hadean Earth.

We have found many Earth-sized worlds but we have no way of determining if their surfaces are Earth-like. This makes it impossible to quantitatively compare habitability, and pretending we can risks damaging the field.

Context. Planet formation with pebbles has been proposed to solve a couple of long-standing issues in the classical formation model. Some sophisticated simulations have been done to confirm the efficiency of pebble accretion. However, there has not been any global N-body simulations that compare the outcomes of planet formation via pebble accretion with observed extrasolar planetary systems. Aims. In this paper, we study the effects of a range of initial parameters of planet formation via pebble accretion, and present the first results of our simulations. Methods. We incorporate the pebble accretion model by Ida et al. (2016) in the N-body code SyMBA (Duncan et al. 1998), along with the effects of gas accretion, eccentricity and inclination damping and planet migration in the disc. Results. We confirm that pebble accretion leads to a variety of planetary systems, but have difficulty in reproducing observed properties of exoplanetary systems, such as planetary mass, semimajor axis, and eccentricity distributions. The main reason behind this is a too-efficient type I migration, which sensitively depends on the disc model. However, our simulations also lead to a few interesting predictions. First, we find that formation efficiencies of planets depend on the stellar metallicities, not only for giant planets, but also for Earths (Es) and Super-Earths (SEs). The dependency for Es/SEs is subtle. Although higher metallicity environments lead to faster formation of a larger number of Es/SEs, they also tend to be lost later via dynamical instability. Second, our results indicate that a wide range of bulk densities observed for Es and SEs is a natural consequence of dynamical evolution of planetary systems. Third, the ejection trend of our simulations suggest that one free-floating E/SE may be expected for two smaller-mass planets.

Once the terrestrial planets had mostly completed their assembly, bombardment continued by planetesimals left-over from accretion. Highly siderophile element (HSE) abundances in Mars' mantle imply its late accretion supplement was 0.8 wt.%; Earth and the Moon obtained an additional 0.7 wt.% and 0.02 wt.%, respectively. The disproportionately high Earth/Moon accretion ratio is explicable by stochastic addition of a few remaining Ceres-sized bodies that preferentially targeted Earth. Here we show that Mars' late accretion budget also requires a colossal impact, a plausible visible remnant of which is the hemispheric dichotomy. The addition of sufficient HSEs to the martian mantle entails an impactor of at least 1200 km in diameter to have struck Mars before ca. 4430 Ma, by which time crust formation was well underway. Thus, the dichotomy could be one of the oldest geophysical features of the martian crust. Ejected debris could be the source material for its satellites.

Overcoming type I migration and preventing low-mass planets from spiralling into the central star is a long-studied topic. It is well known that outward migration is possible in viscous-heated discs relatively close to the central star because the entropy gradient can be sufficiently steep that the positive corotation torque overcomes the negative Lindblad torque. Yet efficiently trapping planets in this region remains elusive. Here we study disc conditions that yield outward migration for low-mass planets under specific planet migration prescriptions. Using a steady-state disc model with a constant alpha-viscosity, outward migration is only possible when the negative temperature gradient exceeds ~0.87. We derive an implicit relation for the maximum mass at which outward migration is possible as a function of viscosity and disc scale height. We apply these criteria, using a simple power-law disc model, to planets that have reached their pebble isolation mass after an episode of rapid accretion. It is possible to trap planets with the pebble isolation mass farther than the inner edge of the disc provided that alpha_crit >~ 0.004 for discs older than 1 Myr. In very young discs the high temperature causes the planets to grow to masses exceeding the maximum for outward migration. As the disc evolves, these more massive planets often reach the central star, generally only towards the end of the disc's lifetime. Saving super-Earths is therefore a delicate interplay between disc viscosity, the opacity profile and the temperature gradient in the viscously heated inner disc.

With approximately one ninth of Earth's mass, Mars is widely considered to be a stranded planetary embryo that never became a fully-grown planet. A currently popular planet formation theory predicts that Mars formed near Earth and Venus and was subsequently scattered outwards to its present location. In such a scenario, the compositions of the three planets are expected to be similar to each other. However, bulk elemental and isotopic data for martian meteorites demonstrate that key aspects of Mars' composition are markedly different from that of Earth. This suggests that Mars formed outside of the terrestrial feeding zone during primary accretion. It is therefore probable that Mars always remained significantly farther from the Sun than Earth; its growth was stunted early and its mass remained relatively low. Here we identify a potential dynamical pathway that forms Mars in the asteroid belt and keeps it outside of Earth's accretion zone while at the same time accounting for strict age and compositional constraints, as well as mass differences. Our uncommon pathway (approximately 2% probability) is based on the Grand Tack scenario of terrestrial planet formation, in which the radial migration by Jupiter gravitationally sculpts the planetesimal disc at Mars' current location. We conclude that Mars' formation requires a specific dynamical pathway, while this is less valid for Earth and Venus. We further predict that} Mars' volatile budget is most likely different from Earth's and that Venus formed close enough to our planet that it is expected to have a nearly identical composition from common building blocks.

This work explores the orbital distribution of minor bodies in the outer Solar System emplaced as a result of a Nice model migration from the simulations of Brasser & Morbidelli (2013). This planetary migration scatters a planetesimal disk from between 29-34 AU and emplaces a population of objects into the Kuiper belt region. From the 2:1 Neptune resonance and outward, the test particles analyzed populate the outer resonances with orbital distributions consistent with trans-Neptunian object (TNO) detections in semi-major axis, inclination, and eccentricity, while capture into the closest resonances is too efficient. The relative populations of the simulated scattering objects and resonant objects in the 3:1 and 4:1 resonances are also consistent with observed populations based on debiased TNO surveys, but the 5:1 resonance is severely underpopulated compared to population estimates from survey results. Scattering emplacement results in the expected orbital distribution for the majority of the TNO populations, however the origin of the large observed population in the 5:1 resonance remains unexplained.

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.

It is generally accepted that silicate-metal (`rocky') planet formation relies on coagulation from a mixture of sub-Mars sized planetary embryos and (smaller) planetesimals that dynamically emerge from the evolving circum-solar disc in the first few million years of our Solar System. Once the planets have, for the most part, assembled after a giant impact phase, they continue to be bombarded by a multitude of planetesimals left over from accretion. Here we place limits on the mass and evolution of these planetesimals based on constraints from the highly siderophile element (HSE) budget of the Moon. Outcomes from a combination of N-body and Monte Carlo simulations of planet formation lead us to four key conclusions about the nature of this early epoch. First, matching the terrestrial to lunar HSE ratio requires either that the late veneer on Earth consisted of a single lunar-size impactor striking the Earth before 4.45 Ga, or that it originated from the impact that created the Moon. An added complication is that analysis of lunar samples indicates the Moon does not preserve convincing evidence for a late veneer like Earth. Second, the expected chondritic veneer component on Mars is 0.06 weight percent. Third, the flux of terrestrial impactors must have been low ( <=10^(-6) M_earth/Myr) to avoid wholesale melting of Earth's crust after 4.4~Ga, and to simultaneously match the number of observed lunar basins. This conclusion leads to an Hadean eon which is more clement than assumed previously. Last, after the terrestrial planets had fully formed, the mass in remnant planetesimals was ~10^(-3) M_earth, lower by at least an order of magnitude than most previous models suggest. Our dynamically and geochemically self-consistent scenario requires that future N-body simulations of rocky planet formation either directly incorporate collisional grinding or rely on pebble accretion.

The Grand Tack model of terrestrial planet formation has emerged in recent years as the premier scenario used to account for several observed features of the inner solar system. It relies on early migration of the giant planets to gravitationally sculpt and mix the planetesimal disc down to ~1 AU, after which the terrestrial planets accrete from material left in a narrow circum-solar annulus. Here we have investigated how the model fares under a range of initial conditions and migration course-change (`tack') locations. We have run a large number of N-body simulations with a tack location of 1.5 AU and 2 AU and tested initial conditions using equal mass planetary embryos and a semi-analytical approach to oligarchic growth. We make use of a recent model of the protosolar disc that takes account of viscous heating, include the full effect of type 1 migration, and employ a realistic mass-radius relation for the growing terrestrial planets. Results show that the canonical tack location of Jupiter at 1.5 AU is inconsistent with the most massive planet residing at 1 AU at greater than 95% confidence. This favours a tack farther out at 2 AU for the disc model and parameters employed. Of the different initial conditions, we find that the oligarchic case is capable of statistically reproducing the orbital architecture and mass distribution of the terrestrial planets, while the equal mass embryo case is not.

Recent observations started revealing the compositions of protostellar discs and planets beyond the Solar System. In this paper, we explore how the compositions of terrestrial planets are affected by dynamical evolution of giant planets. We estimate the initial compositions of building blocks of these rocky planets by using a simple condensation model, and numerically study the compositions of planets formed in a few different formation models of the Solar System. We find that the abundances of refractory and moderately volatile elements are nearly independent of formation models, and that all the models could reproduce the abundances of these elements of the Earth. The abundances of atmophile elements, on the other hand, depend on the scattering rate of icy planetesimals into the inner disc as well as the mixing rate of the inner planetesimal disc. For the classical formation model, neither of these mechanisms are efficient and the accretion of atmophile elements during the final assembly of terrestrial planets appears to be difficult. For the Grand Tack model, both of these mechanisms are efficient, which leads to a relatively uniform accretion of atmophile elements in the inner disc. It is also possible to have a "hybrid" scenario where the mixing is not very efficient but the scattering is efficient. The abundances of atmophile elements in this case increases with orbital radii. Such a scenario may occur in some of the extrasolar planetary systems which are not accompanied by giant planets or those without strong perturbations from giants. We also confirm that the Grand Tack scenario leads to the distribution of asteroid analogues where rocky planetesimals tend to exist interior to icy ones, and show that their overall compositions are consistent with S-type and C-type chondrites, respectively.

The migration and encounter histories of the giant planets in our Solar System can be constrained by the obliquities of Jupiter and Saturn. We have performed secular simulations with imposed migration and N-body simulations with planetesimals to study the expected obliquity distribution of migrating planets with initial conditions resembling those of the smooth migration model, the resonant Nice model and two models with five giant planets initially in resonance (one compact and one loose configuration). For smooth migration, the secular spin-orbit resonance mechanism can tilt Saturn's spin axis to the current obliquity if the product of the migration time scale and the orbital inclinations is sufficiently large (exceeding 30 Myr deg). For the resonant Nice model with imposed migration, it is difficult to reproduce today's obliquity values, because the compactness of the initial system raises the frequency that tilts Saturn above the spin precession frequency of Jupiter, causing a Jupiter spin-orbit resonance crossing. Migration time scales sufficiently long to tilt Saturn generally suffice to tilt Jupiter more than is observed. The full N-body simulations tell a somewhat different story, with Jupiter generally being tilted as often as Saturn, but on average having a higher obliquity. The main obstacle is the final orbital spacing of the giant planets, coupled with the tail of Neptune's migration. The resonant Nice case is barely able to simultaneously reproduce the {orbital and spin} properties of the giant planets, with a probability ~0.15%. The loose five planet model is unable to match all our constraints (probability <0.08%). The compact five planet model has the highest chance of matching the orbital and obliquity constraints simultaneously (probability ~0.3%).

It has long been hypothesised that the Jupiter-family comets (JFCs) come from the scattered disc, an unstable planetesimal population beyond Neptune. This viewpoint has been widely accepted, but a few issues remain, the most prominent of which are the total number of visible JFCs with a perihelion distance q<2.5 AU and the corresponding number of objects in the scattered disc. In this work we give a robust estimate of the number of visible JFCs with q<2.5 AU and diameter D>2.3 km based on recent observational data. This is combined with numerical simulations that use a simple fading law applied to JFCs that come close to the Sun. For this we numerically evolve thousands of comets from the scattered disc through the realm of the giant planets and keep track of their number of perihelion passages with perihelion distance q<2.5 AU, below which the activity is supposed to increase considerably. We can simultaneously fit the JFC inclination and semi-major axis distribution accurately with a delayed power law fading function and the maximum perihelion distance below which the observational data is complete is q~2.3 AU. From observational data we calculate that a JFC with diameter D=2.3 km has a typical total absolute magnitude H_T=10.8, and the steady-state number of active JFCs with diameter D>2.3 km and q<2.5 AU is of the order of 300 (but with large uncertainties), approximately a factor two higher than earlier estimates. The increased JFC population results in a scattered disc population of 6 billion objects and decreases the observed Oort cloud to scattered disc population ratio to 13, virtually the same as the value of 12 obtained with numerical simulations.

The detached object Sedna is likely at the inner edge of the Oort cloud, more precisely the inner Oort cloud (IOC). Until recently it was the sole member of this population. The recent discovery of the detached object 2012 VP113 has confirmed that there should be more objects in this region. Three additional IOC candidates with orbits much closer to Neptune have been proposed in the past decade since Sedna's discovery: 2000 CR105, 2004 VN112 and 2010 GB174. Sedna and 2012 VP113 have perhelia near 80 AU and semi-major axes over 250 AU. The latter three have perihelia between 44 AU and 50 AU and semi-major axes between 200 AU and 400 AU. Here we determine whether the latter three objects belong to the IOC or are from the Kuiper Belt's Extended Scattered Disc (ESD) using numerical simulations. We assume that the IOC was formed when the Sun was in its birth cluster. We analyse the evolution of the IOC and the Scattered Disc (SD) during an episode of late giant planet migration. We examine the impact of giant planet migration in the context of four and five planets. We report that the detached objects 2004 VN112 and 2010 GB174 are likely members of the IOC that were placed there while the Sun was in its birth cluster or during an episode of Solar migration in the Galaxy. The origin of 2000 CR105 is ambiguous but it is likely it belongs to the ESD. Based on our simulations we find that the maximum perihelion distance of SD objects is 41 AU when the semi-major axis is higher than 250 AU. Objects closer in are subject to mean motion resonances with Neptune that may raise their perihelia. The five planet model yields the same outcome. We impose a conservative limit and state that all objects with perihelion distance q>45 AU and semi-major axis a>250 AU belong to the inner Oort cloud.

HARPS and it Kepler results indicate that half of solar-type stars host planets with periods P<100 d and masses M < 30 M_E. These super Earth systems are compact and dynamically cold. Here we investigate the stability of the super Earth system around the K-dwarf HD40307. It could host up to six planets, with one in the habitable zone. We analyse the system's stability using numerical simulations from initial conditions within the observational uncertainties. The most stable solution deviates 3.1 sigma from the published value, with planets e and f not in resonance and planets b and c apsidally aligned. We study the habitability of the outer planet through the yearly-averaged insolation and black-body temperature at the pole. Both undergo large variations because of its high eccentricity and are much more intense than on Earth. The insolation variations are precession dominated with periods of 40 kyr and 102 kyr for precession and obliquity if the rotation period is 3 d. A rotation period of about 1.5 d could cause extreme obliquity variations because of capture in a Cassini state. For faster rotation rates the periods converge to 10 kyr and 20 kyr. The large uncertainty in the precession period does not change the overall outcome.

The origin of the Halley-type comets (HTCs) is one of the last mysteries of the dynamical evolution of the Solar System. Prior investigation into their origin has focused on two source regions: the Oort cloud and the Scattered Disc. From the former it has been difficult to reproduce the non-isotropic, prograde skew in the inclination distribution of the observed HTCs without invoking a multi-component Oort cloud model and specific fading of the comets. The Scattered Disc origin fares better but suffers from needing an order of magnitude more mass than is currently advocated by theory and observations. Here we revisit the Oort cloud origin and include cometary fading. Our observational sample stems from the JPL catalogue. We only keep comets discovered and observed after 1950 but place no a priori restriction on the maximum perihelion distance of observational completeness. We then numerically evolve half a million comets from the Oort cloud through the realm of the giant planets and keep track of their number of perihelion passages with perihelion distance q<2.5AU, below which the activity is supposed to increase considerably. We can simultaneously fit the HTC inclination and semi-major axis distribution very well with a power law fading function of the form m^-k, where m is the number of perihelion passages with q<2.5 AU and k is the fading index. We match both the inclination and semi-major axis distributions when k~1 and the maximum imposed perihelion distance of the observed sample is q~1.8AU. The value of k is higher than the one obtained for the Long-Period Comets (LPCs), with k~0.7. This increase in k is most likely the result of cometary surface processes. We argue the HTC sample is now most likely complete for q<1.8AU. We calculate that the steady-state number of active HTCs with diameter D>2.3km and q<1.8AU is of the order of 100.

Evidence in the Solar System suggests that the giant planets underwent an epoch of radial migration that was very rapid, with an e-folding timescale shorter than 1~Myr. It is probable that the cause of this migration was that the giant planets experienced an orbital instability that caused them to encounter each other, resulting in radial migration. Several works suggest that this dynamical instability occurred `late', long after all the planets had formed and the solar nebula had dissipated. Assuming that the terrestrial planets had already formed, then their orbits would have been affected by the migration of the giant planets. As a result, how did the orbits of the terrestrial planets change? And can we use this migration to obtain information on the primordial orbits of the terrestrial planets? We directly model a large number of terrestrial planet systems and their response to giant planet migration. We study the change in the Angular Momentum Deficit (AMD) of the terrstrials. We conclude that the primordial AMD should have been lower than ~70\% of the current value, but higher than 10\%. We find that a scenario with five giant planets better satisfies the orbital constraints of the terrestrial planets. We predict that Mars was initially on an eccentric and inclined orbit while the orbits of Mercury, Venus and Earth were more circular and coplanar. The lower primordial dynamical excitement and the peculiar partitioning between planets impose new constraints for terrestrial planet formation simulations.

The number of multi-planet systems known to be orbiting their host stars with orbital periods that place them in mean motion resonances is growing. For the most part, these systems are in first-order resonances and dynamical studies have focused their efforts towards understanding the origin and evolution of such dynamically resonant commensurabilities. We report here the discovery of two super-Earths that are close to a second-order dynamical resonance, orbiting the metal-poor ([Fe/H]=-0.43 dex) and inactive G2V star HD41248. We analysed 62 HARPS archival radial velocities for this star, that until now, had exhibited no evidence for planetary companions. Using our new Bayesian Doppler signal detection algorithm, we find two significant signals in the data, with periods of 18.357 days and 25.648 days, indicating they could be part of a 7:5 second-order mean motion resonance. Both semi-amplitudes are below 3m/s and the minimum masses of the pair are 12.3 and 8.6Mearth, respectively. Our simulations found that apsidal alignment stabilizes the system, and even though libration of the resonant angles was not seen, the system is affected by the presence of the resonance and could yet occupy the 7:5 commensurability, which would be the first planetary configuration in such a dynamical resonance. Given the multitude of low-mass multiplanet systems that will be discovered in the coming years, we expect more of these second-order resonant configurations will emerge from the data, highlighting the need for a better understanding of the dynamical interactions between forming planetesimals.

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.

We investigate the obliquity and spin period of Earth-Moon like systems after 4.5 Gyr of tidal evolution with various satellite masses and initial planetary obliquity and discuss their relations to the habitability of the planet. We find three possible outcomes: either i) the system is still evolving, ii) the system is double synchronous or iii) the satellite has collided with the planet. The transition between case i) and ii) is abrupt and occurs at slightly larger satellite mass ($m_s \sim 0.02m_p$) than the lunar mass. We suggest that cases ii) and iii) are less habitable than case i). Using results from models of giant impacts and satellite accretion, we found that the systems that mimic our own with rotation period $12 < P_p < 48$ h and current planetary obliquity $\varepsilon_p < 40^\circ$ or $\varepsilon_p > 140^\circ$ only represent 14% of the possible outcomes. Elser et al. (2011) conclude that the probability of a terrestrial planet having a heavy satellite is 13%. Combining these results suggests that the probability of ending up with a system such as our own is of the order of 2%.

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.

We analyse the origin of three Centaurs with perihelia in the range 15 AU to 30 AU, inclinations above 70 deg and semi-major axes shorter than 100 AU. Based on long-term numerical simulations we conclude that these objects most likely originate from the Oort cloud rather than the Kuiper Belt or Scattered Disc. We estimate that there are currently between 1 and 200 of these high-inclination, high-perihelion Centaurs with absolute magnitude H<8.

We re-examine the formation of the inner Oort comet cloud while the Sun was in its birth cluster with the aid of numerical simulations. This work is a continuation of an earlier study (Brasser et al., 2006) with several substantial modifications. First, the system consisting of stars, planets and comets is treated self-consistently in our N-body simulations, rather than approximating the stellar encounters with the outer Solar System as hyperbolic fly-bys. Second, we have included the expulsion of the cluster gas, a feature that was absent previously. Third, we have used several models for the initial conditions and density profile of the cluster -- either a Hernquist or Plummer potential -- and chose other parameters based on the latest observations of embedded clusters from the literature. {These other parameters result in the stars being on radial orbits and the cluster collapses.} Similar to previous studies, in our simulations the inner Oort cloud is formed from comets being scattered by Jupiter and Saturn and having their pericentres decoupled from the planets by perturbations from the cluster gas and other stars. We find that all inner Oort clouds formed in these clusters have an inner edge ranging from 100 AU to a few hundred AU, and an outer edge at over 100\,000 AU, with little variation in these values for all clusters. All inner Oort clouds formed are consistent with the existence of (90377) Sedna, an inner Oort cloud dwarf planetoid, at the inner edge of the cloud: Sedna tends to be at the innermost 2% for Plummer models, while it is 5% for Hernquist models. We emphasise that the existence of Sedna is a generic outcome. We define a `concentration radius' for the inner Oort cloud and find that its value increases with increasing number of stars in the cluster, ranging from 600 AU to 1500 AU for Hernquist clusters and from 1500 AU to 4000 AU for Plummer clusters...

Agnor & Hamilton (2006) demonstrated that the disruption of a binary was an effective mechanism to capture Triton. The subsequent evolution of Triton's post-capture orbit could have proceeded through gravitational tides. The study by Agnor & Hamilton (2006) is repeated in the framework of the Nice model to determine the post-capture orbit of Triton. After capture it is then subjected to tidal evolution. The perturbations from the Sun and the figure of Neptune are included. The perturbations from the Sun acting on Triton cause it to spend a long time in its high-eccentricity phase, usually of the order of 10 Myr, while the typical time to circularise to its current orbit is some 200 Myr. The current orbit of Triton is consistent with an origin through binary capture and tidal evolution, even though the model prefers Triton to be closer to Neptune than it is today. The probability of capturing Triton in this manner is approximately 0.7%. Since the capture of Triton was at most a 50% event -- since only Neptune has one, but Uranus does not -- we deduce that in the primordial trans-Neptunian disc there were 100 binaries with at least one Triton-sized member. Morbidelli et al. (2009) concludes there were some 1000 Triton-sized bodies in the trans-Neptunian proto-planetary disc, so the primordial binary fraction with at least one Triton-sized member is 10%. This value is consistent with theoretical predictions, but at the low end. If Triton was captured at the same time as Neptune's irregular satellites, the far majority of these, including Nereid, would be lost. This suggests either that Triton was captured on an orbit with a small semi-major axis a < 50 R_N (a rare event), or that it was captured before the dynamical instability of the Nice model, or that some other mechanism was at play. The issue of keeping the irregular satellites remains unresolved.

We performed numerical simulations of the obliquity evolution of Mars during the Noachian era, at which time the giant planets were on drastically different orbits than today. For the preferred primordial configuration of the planets we find that there are two large zones where the Martian obliquity is stable and oscillates with an amplitude lower than 20$^\circ$. These zones occur at obliquities below 30$^\circ$ and above 60$^\circ$; intermediate values show either resonant or chaotic behaviour depending on the primordial orbits of the terrestrial planets.

It has previously been shown that varying the numerical timestep during a symplectic orbital integration leads to a random walk in energy and angular momentum, destroying the phase space-conserving property of symplectic integrators. Here we show that when altering the timestep symplectic correctors can be used to reduce this error to a negligible level. Furthermore, these correctors can also be employed to avoid a large error introduction when changing the Hamiltonian's partitioning. We have constructed a numerical integrator using this technique that is nearly as accurate as widely used fixed-step routines. In addition, our algorithm is drastically faster for integrations of highly eccentricitic, large semimajor axis orbits, such as those found in the Oort Cloud.

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

Using numerical simulations, we show that smooth migration of the giant planets through a planetesimal disk leads to an orbital architecture that is inconsistent with the current one: the resulting eccentricities and inclinations of their orbits are too small. The crossing of mutual mean motion resonances by the planets would excite their orbital eccentricities but not their orbital inclinations. Moreover, the amplitudes of the eigenmodes characterising the current secular evolution of the eccentricities of Jupiter and Saturn would not be reproduced correctly; only one eigenmode is excited by resonance-crossing. We show that, at the very least, encounters between Saturn and one of the ice giants (Uranus or Neptune) need to have occurred, in order to reproduce the current secular properties of the giant planets, in particular the amplitude of the two strongest eigenmodes in the eccentricities of Jupiter and Saturn.