Dynamics of the terrestrial planets from a large number of N-body simulations

The agglomeration of planetary embryos and planetesimals was the final stage of terrestrial planet formation. This process is modeled using N-body accretion simulations, whose outcomes are tested by comparing to observed physical and chemical Solar System properties. The outcomes of these simulations are stochastic, leading to a wide range of results, which makes it difficult at times to identify the full range of possible outcomes for a given dynamic environment. We ran fifty high-resolution simulations each with Jupiter and Saturn on circular or eccentric orbits, whereas most previous studies ran an order of magnitude fewer. This allows us to better quantify the probabilities of matching various observables, including low probability events such as Mars formation, and to search for correlations between properties. We produce many good Earth analogues, which provide information about the mass evolution and provenance of the building blocks of the Earth. Most observables are weakly correlated or uncorrelated, implying that individual evolutionary stages may reflect how the system evolved even if models do not reproduce all of the Solar Systemʹs properties at the end. Thus individual N-body simulations may be used to study the chemistry of planetary accretion as particular accretion pathways may be representative of a given dynamic scenario even if that simulation fails to reproduce many of the other observed traits of the Solar System.