Ramon Brasser

How primordial is the structure of comet 67P? - Combined collisional and dynamical models suggest a late formation

Context. 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. Aims. We investigate the effects of collisions on a comet with a structure like 67P/Churyumov-Gerasimenko (67P). We develop scaling laws for the critical specific impact energies Q reshape 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. Although the focus of this work is on a structure of this kind, the analysis is also performed for more generic bi-lobe shapes, for which we define the critical specific energy Q bil . The simulation outcomes are also analyzed in terms of impact heating and the evolution of the porosity. Methods. The effects of impacts on comet 67P are studied using a state-of-the-art smooth particle hydrodynamics shock physics code. In the 3D simulations, a publicly available shape model of 67P is applied and a range of impact conditions and material properties are investigated. The resulting critical specific impact energy Q reshape (as well as Q bil for generic bi-lobe shapes) 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. Results. We find that in any scenario considered here, comet 67P 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. Conclusions. Our study indicates that the observed bi-lobe structure of comet 67P 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.

Late veneer and late accretion to the terrestrial planets

Abstract 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 ⊕ Myr − 1 ) to avoid wholesale melting of Earths 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 ⊕ , 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.

A Strategy for Origins of Life Research

Contents 1. Introduction 1.1. A workshop and this document 1.2. Framing origins of life science 1.2.1. What do we mean by the origins of life (OoL)? 1.2.2. Defining life 1.2.3. How should we characterize approaches to OoL science? 1.2.4. One path to life or many? 2. A Strategy for Origins of Life Research 2.1. Outcomes—key questions and investigations 2.1.1. Domain 1: Theory 2.1.2. Domain 2: Practice 2.1.3. Domain 3: Process 2.1.4. Domain 4: Future studies 2.2. EON Roadmap 2.3. Relationship to NASA Astrobiology Roadmap and Strategy documents and the European AstRoMap  Appendix I  Appendix II  Supplementary Materials  References

An updated estimate of the number of Jupiter-family comets using a simple fading law

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.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.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 language of exoplanet ranking metrics needs to change

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.

TILTING SATURN WITHOUT TILTING JUPITER: CONSTRAINTS ON GIANT PLANET MIGRATION

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 Saturns 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 todays 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 Neptunes 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%).

N-body simulations of planet formation via pebble accretion - I. First results

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.

Asteroid Family Associations of Active Asteroids

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.

The Structure of the Distant Kuiper Belt in a Nice Model Scenario

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.

The terrestrial late veneer from core disruption of a lunar-sized impactor

Abstract Overabundances in highly siderophile elements (HSEs) of Earths mantle can be explained by conveyance from a singular, immense ( D ∼ 3000 km ) “Late Veneer” impactor of chondritic composition, subsequent to lunar formation and terrestrial core-closure. Such rocky objects of approximately lunar mass (∼0.01 M⊕) 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 (∼45°), the impactors core elongates and thereafter disintegrates into a metallic hail of small particles (∼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 (