Collision Chains among the Terrestrial Planets. II. An Asymmetry between Earth and Venus

During the late stage of terrestrial planet formation, hit-and-run collisions are about as common as accretionary mergers, for expected velocities and angles of giant impacts. Average hit-and-runs leave two major remnants plus debris: the target and impactor, somewhat modified through erosion, escaping at lower relative velocity. Here we continue our study of the dynamical effects of such collisions. We compare the dynamical fates of intact runners that start from hit-and-runs with proto-Venus at 0.7 au and proto-Earth at 1.0 au. We follow the orbital evolutions of the runners, including the other terrestrial planets, Jupiter, and Saturn, in an N-body code. We find that the accretion of these runners can take ≳10 Myr (depending on the egress velocity of the first collision) and can involve successive collisions with the original target planet or with other planets. We treat successive collisions that the runner experiences using surrogate models from machine learning, as in previous work, and evolve subsequent hit-and-runs in a similar fashion. We identify asymmetries in the capture, loss, and interchange of runners in the growth of Venus and Earth. Hit-and-run is a more probable outcome at proto-Venus, being smaller and faster orbiting than proto-Earth. But Venus acts as a sink, eventually accreting most of its runners, assuming typical events, whereas proto-Earth loses about half, many of those continuing to Venus. This leads to a disparity in the style of late-stage accretion that could have led to significant differences in geology, composition, and satellite formation at Earth and Venus

Realistic On-the-fly Outcomes of Planetary Collisions: Machine Learning Applied to Simulations of Giant Impacts

Planet formation simulations are capable of directly integrating the evolution of hundreds to thousands of planetary embryos and planetesimals as they accrete pairwise to become planets. In principle, these investigations allow us to better understand the final configuration and geochemistry of the terrestrial planets, and also to place our solar system in the context of other exosolar systems. While these simulations classically prescribe collisions to result in perfect mergers, recent computational advances have begun to allow for more complex outcomes to be implemented. Here we apply machine learning to a large but sparse database of giant impact studies, which allows us to streamline the simulations into a classifier of collision outcomes and a regressor of accretion efficiency. The classifier maps a four-dimensional (4D) parameter space (target mass, projectile-to-target mass ratio, impact velocity, impact angle) into the four major collision types: merger, graze-and-merge, hit-and-run, and disruption. The definition of the four regimes and their boundary is fully data-driven. The results do not suffer from any model assumption in the fitting. The classifier maps the structure of the parameter space and it provides insights into the outcome regimes. The regressor is a neural network that is trained to closely mimic the functional relationship between the 4D space of collision parameters, and a real-variable outcome, the mass of the largest remnant. This work is a prototype of a more complete surrogate model, that will be based on extended sets of simulations (big data), that will quickly and reliably predict specific collision outcomes for use in realistic N-body dynamical studies of planetary formation.NASA Planetary Science Division; University of ArizonaThis item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]

A New Database of Giant Impacts over a Wide Range of Masses and with Material Strength: A First Analysis of Outcomes

In the late stage of terrestrial planet formation, planets are predicted to undergo pairwise collisions known as giant impacts. Here we present a high-resolution database of giant impacts for differentiated colliding bodies of iron-silicate composition, with target masses ranging from 10^-4 M_Earth up to super-Earths (5 M_Earth). We vary impactor-to-target mass ratio, core-mantle (iron-silicate) fraction, impact velocity, and impact angle. Strength in the form of friction is included in all simulations. We find that due to strength, collisions with bodies smaller than about 2*10^-3 M_Earth can result in irregular shapes, compound core structures, and captured binaries. We observe that the characteristic escaping velocity of smaller remnants (debris) is approximately half of the impact velocity, significantly faster than currently assumed in N-body simulations of planet formation. Incorporating these results in N-body planet formation studies would provide more realistic debris-debris and debris-planet interactions.Comment: Accepted for publication in PSJ; Table 2 is available in full in an ancillary fil

The Effect of Inefficient Accretion on Planetary Differentiation

Pairwise collisions between terrestrial embryos are the dominant means of accretion during the last stage of planet formation. Hence, their realistic treatment in N-body studies is critical to accurately model the formation of terrestrial planets and to develop interpretations of telescopic and spacecraft observations. In this work, we compare the effects of two collision prescriptions on the core-mantle differentiation of terrestrial planets: a model in which collisions are always completely accretionary (``perfect merging'') and a more realistic model based on neural networks that has been trained on hydrodynamical simulations of giant impacts. The latter model is able to predict the loss of mass due to imperfect accretion and the evolution of non-accreted projectiles in hit-and-run collisions. We find that the results of the neural-network model feature a wider range of final core mass fractions and metal-silicate equilibration pressures, temperatures, and oxygen fugacities than the assumption of perfect merging. When used to model collisions in N-body studies of terrestrial planet formation, the two models provide similar answers for planets more massive than 0.1 Earth's masses. For less massive final bodies, however, the inefficient-accretion model predicts a higher degree of compositional diversity. This phenomenon is not reflected in planet formation models of the solar system that use perfect merging to determine collisional outcomes. Our findings confirm the role of giant impacts as important drivers of planetary diversity and encourage a realistic implementation of inefficient accretion in future accretion studies.Comment: 21 pages, 2 tables, 7 figures. Published open access on PSJ: https://iopscience.iop.org/article/10.3847/PSJ/abf0a

Realistic On-the-fly Outcomes of Planetary Collisions II: Bringing Machine Learning to N-body Simulations

Terrestrial planet formation theory is at a bottleneck, with the growing realization that pairwise collisions are treated far too simply. Here, and in our companion paper (Cambioni et al. 2019) that introduces the training methodology, we demonstrate the first application of machine learning to more realistically model the late stage of planet formation by giant impacts. We present surrogate models that give fast, reliable answers for the masses and velocities of the two largest remnants of a giant impact, as a function of the colliding masses and their impact velocity and angle, with the caveat that our training data do not yet include pre-impact rotation or variable thermal conditions. We compare canonical N-body scenarios of terrestrial planet formation assuming perfect merger (Chambers 2001) with our more realistic treatment that includes inefficient accretions and hit-and-run collisions. The result is a protracted tail of final events lasting ~200 Myr, and the conversion of about half the mass of the initial population to debris. We obtain profoundly different solar system architectures, featuring a much wider range of terrestrial planet masses and enhanced compositional diversity.Comment: 20 pages, 10 figures, 3 tables; accepted for publication in ApJ; Table 1 is available in full in an ancillary file; the code developed in this work is available at https://github.com/aemsenhuber/collresolv

After DART: Using the First Full-scale Test of a Kinetic Impactor to Inform a Future Planetary Defense Mission

NASA’s Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ∼10% and characterization of the Didymos binary system. But DART is a single experiment; how could these results be used in a future planetary defense necessity involving a different asteroid? We examine what aspects of Dimorphos’s response to kinetic impact will be constrained by DART results; how these constraints will help refine knowledge of the physical properties of asteroidal materials and predictive power of impact simulations; what information about a potential Earth impactor could be acquired before a deflection effort; and how design of a deflection mission should be informed by this understanding. We generalize the momentum enhancement factor β, showing that a particular direction-specific β will be directly determined by the DART results, and that a related direction-specific β is a figure of merit for a kinetic impact mission. The DART β determination constrains the ejecta momentum vector, which, with hydrodynamic simulations, constrains the physical properties of Dimorphos’s near-surface. In a hypothetical planetary defense exigency, extrapolating these constraints to a newly discovered asteroid will require Earth-based observations and benefit from in situ reconnaissance. We show representative predictions for momentum transfer based on different levels of reconnaissance and discuss strategic targeting to optimize the deflection and reduce the risk of a counterproductive deflection in the wrong direction

After DART: Using the First Full-scale Test of a Kinetic Impactor to Inform a Future Planetary Defense Mission

NASA’s Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ∼10% and characterization of the Didymos binary system. But DART is a single experiment; how could these results be used in a future planetary defense necessity involving a different asteroid? We examine what aspects of Dimorphos’s response to kinetic impact will be constrained by DART results; how these constraints will help refine knowledge of the physical properties of asteroidal materials and predictive power of impact simulations; what information about a potential Earth impactor could be acquired before a deflection effort; and how design of a deflection mission should be informed by this understanding. We generalize the momentum enhancement factor β, showing that a particular direction-specific β will be directly determined by the DART results, and that a related direction-specific β is a figure of merit for a kinetic impact mission. The DART β determination constrains the ejecta momentum vector, which, with hydrodynamic simulations, constrains the physical properties of Dimorphos’s near-surface. In a hypothetical planetary defense exigency, extrapolating these constraints to a newly discovered asteroid will require Earth-based observations and benefit from in situ reconnaissance. We show representative predictions for momentum transfer based on different levels of reconnaissance and discuss strategic targeting to optimize the deflection and reduce the risk of a counterproductive deflection in the wrong direction

Fast boulder fracturing by thermal fatigue detected on stony asteroids.

Spacecraft observations revealed that rocks on carbonaceous asteroids, which constitute the most numerous class by composition, can develop millimeter-to-meter-scale fractures due to thermal stresses. However, signatures of this process on the second-most populous group of asteroids, the S-complex, have been poorly constrained. Here, we report observations of boulders' fractures on Dimorphos, which is the moonlet of the S-complex asteroid (65803) Didymos, the target of NASA's Double Asteroid Redirection Test (DART) planetary defense mission. We show that the size-frequency distribution and orientation of the mapped fractures are consistent with formation through thermal fatigue. The fractures' preferential orientation supports that these have originated in situ on Dimorphos boulders and not on Didymos boulders later transferred to Dimorphos. Based on our model of the fracture propagation, we propose that thermal fatigue on rocks exposed on the surface of S-type asteroids can form shallow, horizontally propagating fractures in much shorter timescales (100 kyr) than in the direction normal to the boulder surface (order of Myrs). The presence of boulder fields affected by thermal fracturing on near-Earth asteroid surfaces may contribute to an enhancement in the ejected mass and momentum from kinetic impactors when deflecting asteroids

Asteroid (101955) Bennu’s weak boulders and thermally anomalous equator

Thermal inertia and surface roughness are proxies for the physical characteristics of planetary surfaces. Global maps of these two properties distinguish the boulder population on near-Earth asteroid (NEA) (101955) Bennu into two types that differ in strength, and both have lower thermal inertia than expected for boulders and meteorites. Neither has strongly temperature-dependent thermal properties. The weaker boulder type probably would not survive atmospheric entry and thus may not be represented in the meteorite collection. The maps also show a high–thermal inertia band at Bennu’s equator, which might be explained by processes such as compaction or strength sorting during mass movement, but these explanations are not wholly consistent with other data. Our findings imply that other C-complex NEAs likely have boulders similar to those on Bennu rather than finer-particulate regoliths. A tentative correlation between albedo and thermal inertia of C-complex NEAs may be due to relative abundances of boulder types

Evidence for multi-fragmentation and mass shedding of boulders on rubble-pile binary asteroid system (65803) Didymos.

Asteroids smaller than 10 km are thought to be rubble piles formed from the reaccumulation of fragments produced in the catastrophic disruption of parent bodies. Ground-based observations reveal that some of these asteroids are today binary systems, in which a smaller secondary orbits a larger primary asteroid. However, how these asteroids became binary systems remains unclear. Here, we report the analysis of boulders on the surface of the stony asteroid (65803) Didymos and its moonlet, Dimorphos, from data collected by the NASA DART mission. The size-frequency distribution of boulders larger than 5 m on Dimorphos and larger than 22.8 m on Didymos confirms that both asteroids are piles of fragments produced in the catastrophic disruption of their progenitors. Dimorphos boulders smaller than 5 m have size best-fit by a Weibull distribution, which we attribute to a multi-phase fragmentation process either occurring during coalescence or during surface evolution. The density per km2 of Dimorphos boulders ≥1 m is 2.3x with respect to the one obtained for (101955) Bennu, while it is 3.0x with respect to (162173) Ryugu. Such values increase once Dimorphos boulders ≥5 m are compared with Bennu (3.5x), Ryugu (3.9x) and (25143) Itokawa (5.1x). This is of interest in the context of asteroid studies because it means that contrarily to the single bodies visited so far, binary systems might be affected by subsequential fragmentation processes that largely increase their block density per km2. Direct comparison between the surface distribution and shapes of the boulders on Didymos and Dimorphos suggest that the latter inherited its material from the former. This finding supports the hypothesis that some asteroid binary systems form through the spin up and mass shedding of a fraction of the primary asteroid