The current properties of small bodies provide important clues to their origin and history. However, how much small bodies were processed by past collisions and to what extent they retain a record of processes that took place during the formation and early evolution of the Solar System is still poorly understood. Here we study the degree of collisional heating and compaction by analysing the large set of previous simulations of small body break-ups by Jutzi et al. (2019), which used porous targets of 50 - 400 km in diameter and investigated a large range of impact velocities, angles as well as energies. We find that the degree of impact processing is generally larger than found in previous studies which considered smaller objects (e.g. Jutzi et al., 2017; Schwartz et al., 2018). However, there is a clear dichotomy in terms of impact processing: the escaping material always experiences stronger heating than the material bound to the largest remnant. Assuming they originate from the same parent body, some of the observed differences between the recently visited asteroids Ryugu and Bennu may be explained by a different location of the material eventually forming these asteroids in the original parent body. Our results also provide constraints on the initial size of cometary nuclei.

The overall shape, internal structure and surface morphology of small bodies such as asteroids and comets are determined to a large degree by the last global-scale impact or disruption event. Depending on the specific energy, impacts lead to a large spectrum of outcomes. Sub-catastrophic disruptions take place in an energy range between cratering impacts and catastrophic disruptions. Although less energetic than catastrophic events, they can still significantly alter the overall shape and structure of the target body. This has been demonstrated recently in the case of bi-lobe cometary nuclei (Jutzi and Benz, 2017). Here we present results of a subsequent study on the shapes of asteroids resulting from such collisions. Sizes ranging from a few hundred meters to a few kilometers are considered. We show that impacts on elongated rotating asteroids frequently lead to the formation of contact binaries. Our results confirm that this mechanism is robust and works for a large range of asteroid sizes and impact velocities. Scaling-laws for the prediction of the size and velocity dependent specific energies required for successful bi-lobe formation are presented. Based on these scaling laws, the expected frequency of such sub-catastrophic impacts is calculated and is compared to the one of catastrophic disruptions, which require much higher specific energies and are more rare. Our analysis suggest that the shapes and structures of a large fraction of small asteroids as observed today may be the result of the last major sub-catastrophic impact.

About 10% of the observed asteroids have rotational periods lower than P = 3 h and they seem to be relatively close to the spin barrier. Yet, the rotation has often been neglected in simulations of asteroid collisions. To determine the effect of rotation, we perform a large number of SPH/N-body impact simulations with rotating targets. We developed a new unified SPH/N-body code with self-gravity, suitable for simulations of both fragmentation phase and gravitational reaccumulation. The code has been verified against previous ones (Benz and Asphaug 1994), but we also tested new features, e.g. rotational stability, tensile stability, etc. Using the new code, we ran simulations with D_pb = 10 km and 100 km monolithic targets and compared synthetic asteroid families created by these impacts with families corresponding to non-rotating targets. The rotation affects mostly cratering events at oblique impact angles. The total mass ejected by these collision can be up to five times larger for rotating targets. We further compute the transfer of the angular momentum and determine conditions under which impacts accelerate or decelerate the target. While individual cratering collisions can cause both acceleration and deceleration, the deceleration prevails on average, collisions thus cause a systematic spin-down of asteroid population.

Understanding the collisional fragmentation and subsequent reaccumulation of fragments is crucial for studies of the formation and evolution of the small-body populations. Using an SPH / N-body approach, we investigate the size-frequency distributions (SFDs) resulting from the disruption of 100 km-diameter targets consisting of porous material, including the effects of pore-crushing as well as friction. Overall, the porous targets have a significantly higher impact strength (Q*D) than the rubble-pile parent bodies investigated previously (Benavidez et al., 2012) and show a behavior more similar to non-porous monolithic targets (Durda et al., 2007). Our results also confirm that for a given specific impact energy, the SFDs resulting from a parent body disruption are strongly dependent on the size scale.

The Cassini spacecraft revealed the spectacular, highly irregular shapes of the small inner moons of Saturn, ranging from the unique "ravioli-like" forms of Pan and Atlas to the highly elongated structure of Prometheus. Closest to Saturn, these bodies provide important clues regarding the formation process of small moons in close orbits around their host planet, but their range of irregular shapes has not been explained yet. Here we show that the spectrum of shapes among Saturn's small moons is a natural outcome of merging collisions among similar-sized moonlets possessing physical properties and orbits that are consistent with those of the current moons. A significant fraction of such merging collisions take place either at the first encounter or after 1-2 hit-and-run events, with impact velocities in the range of 1-5 times the mutual escape velocity. Close to head-on mergers result in flattened objects with large equatorial ridges, as observed on Atlas and Pan. With slightly more oblique impact angles, collisions lead to elongated, Prometheus-like shapes. These results suggest that the current forms of the small moons provide direct evidence of the processes at the final stages of their formation, involving pairwise encounters of moonlets of comparable size. Finally, we show that this mechanism may also explain the formation of Iapetus' equatorial ridge, as well as its oblate shape.

For decades the limited thickness of Vesta's basaltic crust, revealed by the link between the asteroid and the howardite-eucrite-diogenite family of meteorites, and its survival to collisional erosion offered an important constraint for the study of the early evolution of the Solar System. Some results of the Dawn mission, however, cast doubts on our understanding of Vesta's interior composition and of the characteristics of its basaltic crust, weakening this classical constraint. In this work we investigate the late accretion and erosion experienced by Vesta's crust after its differentiation and recorded in the composition of eucrites and diogenites and show that it offers an astrochemical window into the earliest evolution of the Solar System. In our proof-of-concept case study focusing on the late accretion and erosion of Vesta's crust during the growth and migration of Jupiter, the water enrichment of eucrites appears to be a sensitive function of Jupiter's migration while the enrichment in highly-siderophile elements of diogenites appears to be particularly sensitive to the size-frequency distribution of the planetesimals. The picture depicted by the enrichments created by late accretion in eucrites and diogenites is not qualitatively affected by the uncertainty on the primordial mass of Vesta. Crustal erosion, instead, is more significantly affected by said uncertainty and Vesta's crust survival appears to be mainly useful to study violent collisional scenarios where highly energetic impacts can strip significant amounts of vestan material while limitedly contributing to Vesta's late accretion. Our results suggest that the astrochemical record of the late accretion and erosion of Vesta's crust provided by eucrites and diogenites can be used as a tool to investigate any process or scenario associated to the evolution of primordial Vesta and of the early Solar System.

Asteroid (162173) Ryugu is the target object of Hayabusa2, an asteroid exploration and sample return mission led by Japan Aerospace Exploration Agency (JAXA). Ground-based observations indicate that Ryugu is a C-type near-Earth asteroid with a diameter of less than 1 km, but the knowledge of its detailed properties is still very limited. This paper summarizes our best understanding of the physical and dynamical properties of Ryugu based on remote sensing and theoretical modeling. This information is used to construct a design reference model of the asteroid that is used for formulation of mission operations plans in advance of asteroid arrival. Particular attention is given to the surface properties of Ryugu that are relevant to sample acquisition. This reference model helps readers to appropriately interpret the data that will be directly obtained by Hayabusa2 and promotes scientific studies not only for Ryugu itself and other small bodies but also for the Solar System evolution that small bodies shed light on.

We present a study of an impacting descent probe that increases the science return of spacecraft orbiting or passing an atmosphere-less planetary body of the solar system, such as the Galilean moons of Jupiter. The descent probe is a carry-on small spacecraft (< 100 kg), to be deployed by the mother spacecraft, that brings itself onto a collisional trajectory with the targeted planetary body in a simple manner. A possible science payload includes instruments for surface imaging, characterisation of the neutral exosphere, and magnetic field and plasma measurement near the target body down to very low-altitudes (~1 km), during the probe's fast (~km/s) descent to the surface until impact. The science goals and the concept of operation are discussed with particular reference to Europa, including options for flying through water plumes and after-impact retrieval of very-low altitude science data. All in all, it is demonstrated how the descent probe has the potential to provide a high science return to a mission at a low extra level of complexity, engineering effort, and risk. This study builds upon earlier studies for a Callisto Descent Probe (CDP) for the former Europa-Jupiter System Mission (EJSM) of ESA and NASA, and extends them with a detailed assessment of a descent probe designed to be an additional science payload for the NASA Europa Mission.

Giant impacts have been suggested to explain various characteristics of terrestrial planets and their moons. However, so far in most models only the immediate effects of the collisions have been considered, while the long-term interior evolution of the impacted planets was not studied. Here we present a new approach, combining 3-D shock physics collision calculations with 3-D thermochemical interior evolution models. We apply the combined methods to a demonstration example of a giant impact on a Mars-sized body, using typical collisional parameters from previous studies. While the material parameters (equation of state, rheology model) used in the impact simulations can have some effect on the long-term evolution, we find that the impact angle is the most crucial parameter for the resulting spatial distribution of the newly formed crust. The results indicate that a dichotomous crustal pattern can form after a head-on collision, while this is not the case when considering a more likely grazing collision. Our results underline that end-to-end 3-D calculations of the entire process are required to study in the future the effects of large-scale impacts on the evolution of planetary interiors.

We model large-scale ($\approx$2000km) impacts on a Mars-like planet using a Smoothed Particle Hydrodynamics code. The effects of material strength and of using different Equations of State on the post-impact material and temperature distributions are investigated. The properties of the ejected material in terms of escaping and disc mass are analysed as well. We also study potential numerical effects in the context of density discontinuities and rigid body rotation. We find that in the large-scale collision regime considered here (with impact velocities of 4km/s), the effect of material strength is substantial for the post-impact distribution of the temperature and the impactor material, while the influence of the Equation of State is more subtle and present only at very high temperatures.

We examine the relevance of tidal heating for large Trans-Neptunian Objects, with a focus on its potential to melt and maintain layers of subsurface liquid water. Depending on their past orbital evolution, tidal heating may be an important part of the heat budget for a number of discovered and hypothetical TNO systems and may enable formation of, and increased access to, subsurface liquid water. Tidal heating induced by the process of despinning is found to be particularly able to compete with heating due to radionuclide decay in a number of different scenarios. In cases where radiogenic heating alone may establish subsurface conditions for liquid water, we focus on the extent by which tidal activity lifts the depth of such conditions closer to the surface. While it is common for strong tidal heating and long lived tides to be mutually exclusive, we find this is not always the case, and highlight when these two traits occur together.

The origin of the particular shape of a small body like comet 67P/Churyumov-Gerasimenko (67P/C-G) is a topic of active research. How and when it acquired its peculiar characteristics has distinct implications on the origin of the solar system and its dynamics. We investigate how shapes like the one of comet 67P/C-G can result from a new type of low-energy, sub-catastrophic impacts involving elongated, rotating bodies. We focus on parameters potentially leading to bi-lobed structures. We also estimate the probability for such structures to survive subsequent impacts. We use a smooth particle hydrodynamics (SPH) shock physics code to model the impacts, the subsequent reaccumulation of material and the reconfiguration into a stable final shape. The energy increase as well as the degree of compaction of the resulting bodies are tracked in the simulations. Our modelling results suggest that the formation of bi-lobed structures like 67P/C-G is a natural outcome of the low energy, sub-catastrophic collisions considered here. Sub-catastrophic impacts have the potential to alter the shape of a small body significantly, without leading to major heating or compaction. The currently observed shapes of cometary nuclei, such as 67P/C-G, maybe a result of such a last major shape forming impact.

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.

We define a set of procedures to numerically study the fate of ejecta produced by the impact of an artificial projectile with the aim of deflecting an asteroid. Here we develop a simplified, idealized model of impact conditions that can be adapted to fit the details of specific deflection-test scenarios, such as what is being proposed for the AIDA project. Ongoing studies based upon the methodology described here can be used to inform observational strategies and safety conditions for an observing spacecraft. To account for ejecta evolution, the numerical strategies we are employing are varied and include a large N-Body component, a smoothed-particle hydrodynamics (SPH) component, and an application of impactor scaling laws. Simulations that use SPH-derived initial conditions show high-speed ejecta escaping at low angles of inclination, and very slowly moving ejecta lofting off the surface at higher inclination angles, some of which re-impacts the small-body surface. We are currently investigating the realism of this and other models' behaviors. Next steps will include the addition of solar perturbations to the model and applying the protocol developed here directly to specific potential mission concepts such as the proposed AIDA scenario.

Global scale impacts modify the physical or thermal state of a substantial fraction of a target asteroid. Specific effects include accretion, family formation, reshaping, mixing and layering, shock and frictional heating, fragmentation, material compaction, dilatation, stripping of mantle and crust, and seismic degradation. Deciphering the complicated record of global scale impacts, in asteroids and meteorites, will lead us to understand the original planet-forming process and its resultant populations, and their evolution in time as collisions became faster and fewer. We provide a brief overview of these ideas, and an introduction to models.

As a complement to experimental and theoretical approaches, numerical modeling has become an important component to study asteroid collisions and impact processes. In the last decade, there have been significant advances in both computational resources and numerical methods. We discuss the present state-of-the-art numerical methods and material models used in "shock physics codes" to simulate impacts and collisions and give some examples of those codes. Finally, recent modeling studies are presented, focussing on the effects of various material properties and target structures on the outcome of a collision.

In the last decade, thanks to the development of sophisticated numerical codes, major breakthroughs have been achieved in our understanding of the formation of asteroid families by catastrophic disruption of large parent bodies. In this review, we describe numerical simulations of asteroid collisions that reproduced the main properties of families, accounting for both the fragmentation of an asteroid at the time of impact and the subsequent gravitational interactions of the generated fragments. The simulations demonstrate that the catastrophic disruption of bodies larger than a few hundred meters in diameter leads to the formation of large aggregates due to gravitational reaccumulation of smaller fragments, which helps explain the presence of large members within asteroid families. Thus, for the first time, numerical simulations successfully reproduced the sizes and ejection velocities of members of representative families. Moreover, the simulations provide constraints on the family dynamical histories and on the possible internal structure of family members and their parent bodies.

We present recent improvements of the modeling of the disruption of strength dominated bodies using the Smooth Particle Hydrodynamics (SPH) technique. The improvements include an updated strength model and a friction model, which are successfully tested by a comparison with laboratory experiments. In the modeling of catastrophic disruptions of asteroids, a comparison between old and new strength models shows no significant deviation in the case of targets which are initially non-porous, fully intact and have a homogeneous structure (such as the targets used in the study by Benz&Asphaug (1999). However, for many cases (e.g. initially partly or fully damaged targets, rubble-pile structures, etc.) we find that it is crucial that friction is taken into account and the material has a pressure dependent shear strength. Our investigations of the catastrophic disruption threshold $Q^*_{D}$ as a function of target properties and target sizes up to a few 100 km show that a fully damaged target modeled without friction has a $Q^*_{D}$ which is significantly (5-10 times) smaller than in the case where friction is included. When the effect of the energy dissipation due to compaction (pore crushing) is taken into account as well, the targets become even stronger ($Q^*_{D}$ is increased by a factor of 2-3). On the other hand, cohesion is found to have an negligible effect at large scales and is only important at scales $\lesssim$ 1km. Our results show the relative effects of strength, friction and porosity on the outcome of collisions among small ($\lesssim$ 1000 km) bodies. These results will be used in a future study to improve existing scaling laws for the outcome of collisions (e.g. Leinhardt&Stewart, 2012).

Numerical simulations of asteroid break-ups, including both the fragmentation of the parent body and the gravitational interactions between the fragments, have allowed us to reproduce successfully the main properties of asteroid families formed in different regimes of impact energy, starting from a non-porous parent body. In this paper, using the same approach, we concentrate on a single regime of impact energy, the so-called catastrophic threshold usually designated by Q*D, which results in the escape of half of the target's mass. Thanks to our recent implementation of a model of fragmentation of porous materials, we can characterize Q*D for both porous and non-porous targets with a wide range of diameters. We can then analyze the potential influence of porosity on the value of Q*D, and by computing the gravitational phase of the collision in the gravity regime, we can characterize the collisional outcome in terms of the fragment size and ejection speed distributions, which are the main outcome properties used by collisional models to study the evolutions of the different populations of small bodies. We also check the dependency of Q*D on the impact speed of the projectile. In the strength regime, which corresponds to target sizes below a few hundreds of meters, we find that porous targets are more difficult to disrupt than non-porous ones. In the gravity regime, the outcome is controlled purely by gravity and porosity in the case of porous targets. In the case of non-porous targets, the outcome also depends on strength. We then propose some power-law relationships between Q*D and both target's size and impact speed that can be used in collisional evolution models.

In this paper, we extend our Smooth Particle Hydrodynamics (SPH) impact code to include the effect of porosity at a sub-resolution scale by adapting the so-called $P-alpha$ model. Many small bodies in the different populations of asteroids and comets are believed to contain a high degree of porosity and the determination of both their collisional evolution and the outcome of their disruption requires that the effect of porosity is taken into account in the computation of those processes. Here, we present our model and show how porosity interfaces with the elastic-perfectly plastic material description and the brittle fracture model generally used to simulate the fragmentation of non-porous rocky bodies. We investigate various compaction models and discuss their suitability to simulate the compaction of (highly) porous material. Then, we perform simple test cases where we compare results of the simulations to the theoretical solutions. We also present a Deep Impact-like simulation to show the effect of porosity on the outcome of an impact. Detailed validation tests will be presented in a next paper by comparison with high-velocity laboratory experiments on porous materials (Jutzi et al., in preparation). Once validated at small scales, our new impact code can then be used at larger scales to study impacts and collisions involving brittle solids including porosity, such as the parent bodies of C-type asteroid families or cometary materials, both in the strength- and in the gravity-dominated regime.