During close approaches to planets or stars, the morphological and dynamical properties of rubble-pile small bodies can be modified, and some may catastrophically break up. This phenomenon is of particular interest for understanding the evolution and population of small bodies, and for making predictions for future encounters. Previous numerical explorations typically used methods that do not adequately represent the nature of rubble piles. The encounter outcomes and influence factors are still poorly constrained. Based on recent advances in modeling rubble piles, we aim to provide a better understanding of the tidal encounter processes through soft-sphere discrete element modeling (SSDEM) and to establish a database of encounter outcomes and the dependencies on encounter conditions and rubble-pile properties. We performed thousands of numerical simulations using the SSDEM implemented in the N-body code pkdgrav to study the dynamical evolution of rubble piles during close encounters with the Earth. The effects of encounter conditions, material strength, arrangement, and resolution of constituent particles are explored. Three typical tidal encounter outcomes are classified, i.e., deformation, mass shedding, and disruption, ranging from mild modifications to severe damages of the progenitor. The encounter speed and distance required for causing disruption events are much smaller than those predicted by previous studies, indicating a smaller creation rate of tidally disrupted small body populations. Extremely elongated fragments with axis ratios ~1:6 can be formed by moderate encounters. Our analyses of the spin-shape evolution of the largest remnants reveal reshaping mechanisms of rubble piles in response to tidal forces, consistent with stable rubble-pile configurations derived by continuum theory. A case study for SL9 suggests a low bulk density (0.2-0.3 g/cc) for its progenitor.

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.

We extend the most recent orbital and absolute magnitude Near Earth Object (NEO) model (Granvik et al., 2018) to provide a statistical description of NEO geometric albedos. Our model is calibrated on NEOWISE albedo data for the NEO population and reproduces these data very well once a simple model for the NEOWISE observational biases is applied. The results are consistent with previous estimates. There are about 1,000 NEOs with diameter D>1km and the mean albedo to convert absolute magnitude into diameter is 0.147. We do not find any statistically significant evidence that the albedo distribution of NEOs depends on NEO size. Instead, we find evidence that the disruption of NEOs at small perihelion distances found in Granvik et al. (2016) occurs preferentially for dark NEOs. The interval between km-sized bodies striking the Earth should occur on average once every 750,000 years. Low and high albedo NEOs are well mixed in orbital space, but a trend remains with higher albedo objects being at smaller semimajor axes and lower albedo objects more likely found at larger semimajor axes.

With an estimated diameter in the 320 to 350 km range, (704) Interamnia is the fifth largest main belt asteroid and one of the few bodies that fills the gap in size between the four largest bodies with $D$ > 400 km (Ceres, Vesta, Pallas and Hygiea) and the numerous smaller bodies with $D$ $\lesssim$ 200 km. However, despite its large size, little is known about the shape and spin state of Interamnia and, therefore, about its bulk composition and past collisional evolution. We aimed to test at what size and mass the shape of a small body departs from a nearly ellipsoidal equilibrium shape (as observed in the case of the four largest asteroids) to an irregular shape as routinely observed in the case of smaller ($D$ $\lesssim$ 200 km) bodies. We observed Interamnia as part of our ESO VLT/SPHERE large program (ID: 199.C-0074) at thirteen different epochs. In addition, several new optical lightcurves were recorded. These data, along with stellar occultation data from the literature, were fed to the All-Data Asteroid Modeling (ADAM) algorithm to reconstruct the 3D-shape model of Interamnia and to determine its spin state. Interamnia's volume-equivalent diameter of 332 $\pm$ 6 km implies a bulk density of $\rho$=1.98 $\pm$ 0.68 gcm$^{-3}$ , which suggests that Interamnia - like Ceres and Hygiea - contains a high fraction of water ice, consistent with the paucity of apparent craters. Our observations reveal a shape that can be well approximated by an ellipsoid, and that is compatible with a fluid hydrostatic equilibrium at the 2 $\sigma$ level. The rather regular shape of Interamnia implies that the size and mass limit, under which the shapes of minor bodies with a high amount of water ice in the subsurface become irregular, has to be searched among smaller ($D$ $\lesssim$ 300km) less massive ($m$ $\lesssim$ 3x10$^{19}$ kg) bodies.

In keeping with the Luxembourg government's initiative to support the future use of space resources, ASIME 2018 was held in Belval, Luxembourg on April 16-17, 2018. The goal of ASIME 2018: Asteroid Intersections with Mine Engineering, was to focus on asteroid composition for advancing the asteroid in-space resource utilisation domain. What do we know about asteroid composition from remote-sensing observations? What are the potential caveats in the interpretation of Earth-based spectral observations? What are the next steps to improve our knowledge on asteroid composition by means of ground-based and space-based observations and asteroid rendez-vous and sample return missions? How can asteroid mining companies use this knowledge? ASIME 2018 was a two-day workshop of almost 70 scientists and engineers in the context of the engineering needs of space missions with in-space asteroid utilisation. The 21 Questions from the asteroid mining companies were sorted into the four asteroid science themes: 1) Potential Targets, 2) Asteroid-Meteorite Links, 3) In-Situ Measurements and 4) Laboratory Measurements. The Answers to those Questions were provided by the scientists with their conference presentations and collected by A. Graps or edited directly into an open-access collaborative Google document or inserted by A. Graps using additional reference materials. During the ASIME 2018, first day and second day Wrap-Ups, the answers to the questions were discussed further. New readers to the asteroid mining topic may find the Conversation boxes and the Mission Design discussions especially interesting.

162173 Ryugu, the target of Hayabusa2, has a round shape with an equatorial ridge, which is known as a spinning top-shape. A strong centrifugal force is a likely contributor to Ryugu's top-shaped features. Observations by Optical Navigation Camera onboard Hayabusa2 show a unique longitudinal variation in geomorphology; the western side of this asteroid, later called the western bulge, has a smooth surface and a sharp equatorial ridge, compared to the other side. Here, we propose a structural deformation process that generated the western bulge. Applying the mission-derived shape model, we employ a finite element model technique to analyze the locations that experience structural failure within the present shape. Assuming that materials are uniformly distributed, our model shows the longitudinal variation in structurally failed regions when the spin period is shorter than ~3.75 h. Ryugu is structurally intact in the subsurface region of the western bulge while other regions are sensitive to structural failure. We infer that this variation is indicative of the deformation process that occurred in the past, and the western bulge is more relaxed structurally than the other region. Our analysis also shows that this deformation process might occur at a spin period between ~3.5 h and ~3.0 h, providing the cohesive strength ranging between ~4 Pa and ~10 Pa.

Asteroid (7) Iris is an ideal target for disk-resolved imaging owing to its brightness (V$\sim$7-8) and large angular size of 0.33 arcsec during its apparitions. Iris is believed to belong to the category of large unfragmented asteroids that avoided internal differentiation, implying that its current shape and topography may record the first few 100 Myr of the solar system's collisional evolution. We recovered information about the shape and surface topography of Iris from disk-resolved VLT/SPHERE/ZIMPOL images acquired in the frame of our ESO large program. We used the All-Data Asteroid Modeling (ADAM) shape reconstruction algorithm to model the 3D shape of Iris, using optical disk-integrated data and disk-resolved images from SPHERE as inputs. We analyzed the SPHERE images to infer the asteroid's global shape and the morphology of its main craters. We present the volume-equivalent diameter D$_{{\rm eq}}$=214$\pm$5 km, and bulk density $\rho$=2.7$\pm$0.3 g cm$^{-3}$ of Iris. Its shape appears to be consistent with that of an oblate spheroid with a large equatorial excavation. We identified eight putative surface features 20--40 km in diameter detected at several epochs, which we interpret as impact craters. Craters on Iris have depth-to-diameter ratios that are similar to those of analogous 10 km craters on Vesta. The bulk density of Iris is consistent with that of its meteoritic analog, namely LL ordinary chondrites. Considering the absence of a collisional family related to Iris and the number of large craters on its surface, we suggest that its equatorial depression may be the remnant of an ancient (at least 3 Gyr) impact. Iris's shape further opens the possibility that large planetesimals formed as almost perfect oblate spheroids. Finally, we attribute the difference in crater morphology between Iris and Vesta to their different surface gravities.

Context. Over the past decades, several interplanetary missions have studied small bodies in situ, leading to major advances in our understanding of their geological and geophysical properties. These missions, however, have had a limited number of targets. Among them, the NASA Dawn mission has characterised in detail the topography and albedo variegation across the surface of asteroid (4) Vesta down to a spatial resolution of approximately 20 m/pixel scale. Aims. Here our aim was to determine how much topographic and albedo information can be retrieved from the ground with VLT/SPHERE in the case of Vesta, having a former space mission (Dawn) providing us with the ground truth that can be used as a benchmark. Methods. We observed Vesta with VLT/SPHERE/ZIMPOL as part of our ESO large programme (ID 199.C-0074) at six different epochs, and deconvolved the collected images with a parametric point spread function (PSF). We then compared our images with synthetic views of Vesta generated from the 3D shape model of the Dawn mission, on which we projected Vesta's albedo information. Results. We show that the deconvolution of the VLT/SPHERE images with a parametric PSF allows the retrieval of the main topographic and albedo features present across the surface of Vesta down to a spatial resolution of approximately 20-30 km. Contour extraction shows an accuracy of approximately 1 pixel (3.6 mas). The present study provides the very first quantitative estimate of the accuracy of ground-based adaptive-optics imaging observations of asteroid surfaces. Conclusions. In the case of Vesta, the upcoming generation of 30-40m telescopes (ELT, TMT, GMT) should in principle be able to resolve all of the main features present across its surface, including the troughs and the north-south crater dichotomy, provided that they operate at the diffraction limit.

Context. CM-like asteroids (Ch and Cgh classes) are a major population within the broader C-complex, encompassing about 10% of the mass of the main asteroid belt. Their internal structure has been predicted to be homogeneous, based on their compositional similarity as inferred from spectroscopy (Vernazza et al., 2016, AJ 152, 154) and numerical modeling of their early thermal evolution (Bland & Travis, 2017, Sci. Adv. 3, e1602514). Aims. Here we aim to test this hypothesis by deriving the density of the CM-like asteroid (41) Daphne from detailed modeling of its shape and the orbit of its small satellite. Methods. We observed Daphne and its satellite within our imaging survey with the Very Large Telescope extreme adaptive-optics SPHERE/ZIMPOL camera (ID 199.C-0074, PI P. Vernazza) and complemented this data set with earlier Keck/NIRC2 and VLT/NACO observations. We analyzed the dynamics of the satellite with our Genoid meta-heuristic algorithm. Combining our high-angular resolution images with optical lightcurves and stellar occultations, we determine the spin period, orientation, and 3-D shape, using our ADAM shape modeling algorithm. Results. The satellite orbits Daphne on an equatorial, quasi-circular, prograde orbit, like the satellites of many other large main-belt asteroids. The shape model of Daphne reveals several large flat areas that could be large impact craters. The mass determined from this orbit combined with the volume computed from the shape model implies a density for Daphne of 1.77+/-0.26 g/cm3 (3 {\sigma}). This density is consistent with a primordial CM-like homogeneous internal structure with some level of macroporosity (~17%). Conclusions. Based on our analysis of the density of Daphne and 75 other Ch/Cgh-type asteroids gathered from the literature, we conclude that the primordial internal structure of the CM parent bodies was homogeneous.

Asteroid (16) Psyche is the target of the NASA Psyche mission. It is considered one of the few main-belt bodies that could be an exposed proto-planetary metallic core and that would thus be related to iron meteorites. Such an association is however challenged by both its near- and mid-infrared spectral properties and the reported estimates of its density. Here, we aim to refine the density of (16) Psyche to set further constraints on its bulk composition and determine its potential meteoritic analog. We observed (16) Psyche with ESO VLT/SPHERE/ZIMPOL as part of our large program (ID 199.C-0074). We used the high angular resolution of these observations to refine Psyche's three-dimensional (3D) shape model and subsequently its density when combined with the most recent mass estimates. In addition, we searched for potential companions around the asteroid. We derived a bulk density of 3.99\,$\pm$\,0.26\,g$\cdot$cm$^{-3}$ for Psyche. While such density is incompatible at the 3-sigma level with any iron meteorites ($\sim$7.8\,g$\cdot$cm$^{-3}$), it appears fully consistent with that of stony-iron meteorites such as mesosiderites (density $\sim$4.25\,$\cdot$cm$^{-3}$). In addition, we found no satellite in our images and set an upper limit on the diameter of any non-detected satellite of 1460\,$\pm$\,200}\,m at 150\,km from Psyche (0.2\%\,$\times$\,R$_{Hill}$, the Hill radius) and 800\,$\pm$\,200\,m at 2,000\,km (3\%\,$\times$\,$R_{Hill}$). Considering that the visible and near-infrared spectral properties of mesosiderites are similar to those of Psyche, there is merit to a long-published initial hypothesis that Psyche could be a plausible candidate parent body for mesosiderites.

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 debiased absolute-magnitude and orbit distributions as well as source regions for near-Earth objects (NEOs) provide a fundamental frame of reference for studies of individual NEOs and more complex population-level questions. We present a new four-dimensional model of the NEO population that describes debiased steady-state distributions of semimajor axis, eccentricity, inclination, and absolute magnitude $H$ in the range $17<H<25$. The modeling approach improves upon the methodology originally developed by Bottke et al. (2000; Science 288, 2190-2194) in that it is, for example, based on more realistic orbit distributions and uses source-specific absolute-magnitude distributions that allow for a power-law slope that varies with $H$. We divide the main asteroid belt into six different entrance routes or regions (ER) to the NEO region: the $\nu_6$, 3:1J, 5:2J and 2:1J resonance complexes as well as Hungarias and Phocaeas. In addition we include the Jupiter-family comets as the primary cometary source of NEOs. We calibrate the model against NEO detections by Catalina Sky Surveys' stations 703 and G96 during 2005-2012, and utilize the complementary nature of these two systems to quantify the systematic uncertainties associated to the resulting model. The most important ERs are the $\nu_6$ and 3:1J resonance complexes with JFCs contributing a few percent of NEOs on average. We predict that there are $962^{+52}_{-56}$ ($802^{+48}_{-42}\times10^3$) NEOs with $H<17.75$ ($H<25$) and these numbers are in agreement with the most recent estimates found in the literature (the uncertainty estimates only account for the random component).

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.

Binary near-Earth asteroid (65803) Didymos is the target of the proposed NASA Double Asteroid Redirection Test (DART), part of the Asteroid Impact & Deflection Assessment (AIDA) mission concept. In this mission, the DART spacecraft is planned to impact the secondary body of Didymos, perturbing mutual dynamics of the system. The primary body is currently rotating at a spin period close to the spin barrier of asteroids, and materials ejected from the secondary due to the DART impact are likely to reach the primary. These conditions may cause the primary to reshape, due to landslides, or internal deformation, changing the permanent gravity field. Here, we propose that if shape deformation of the primary occurs, the mutual orbit of the system would be perturbed due to a change in the gravity field. We use a numerical simulation technique based on the full two-body problem to investigate the shape effect on the mutual dynamics in Didymos after the DART impact. The results show that under constant volume, shape deformation induces strong perturbation in the mutual motion. We find that the deformation process always causes the orbital period of the system to become shorter. If surface layers with a thickness greater than ~0.4 m on the poles of the primary move down to the equatorial region due to the DART impact, a change in the orbital period of the system and in the spin period of the primary will be detected by ground-based measurement.

Solar system small bodies come in a wide variety of shapes and sizes, which are achieved following very individual evolutional paths through billions of years. This paper focuses on the reshaping process of rubble-pile asteroids driven by meteorite impacts. In our study, numerous possible equilibrium configurations are obtained via Monte Carlo simulation, and the structural stability of these configurations is determined via eigen analysis of the geometric constructions. The eigen decomposition reveals a connection between the cluster's reactions and the types of external disturbance. Numerical simulations are performed to verify the analytical results. The gravitational N-body code pkdgrav is used to mimic the responses of the cluster under intermittent non-dispersive impacts. We statistically confirm that the stability index, the total gravitational potential and the volume of inertia ellipsoid show consistent tendency of variation. A common regime is found in which the clusters tend towards crystallization under intermittent impacts, i.e., only the configurations with high structural stability survive under the external disturbances. The results suggest the trivial non-disruptive impacts might play an important role in the rearrangement of the constituent blocks, which may strengthen these rubble piles and help to build a robust structure under impacts of similar magnitude. The final part of this study consists of systematic simulations over two parameters, the projectile momentum and the rotational speed of the cluster. The results show a critical value exists for the projectile momentum, as predicted by theory, below which all clusters become responseless to external disturbances; and the rotation proves to be significant for it exhibits an "enhancing" effect on loose-packed clusters, which coincides with the observation that several fast-spinning asteroids have low bulk densities.

This paper presents a quantitative study of the evolution of the ejecta cloud released from a hypervelocity impact on a binary asteroid. We performed numerical simulations of the post-impact dynamics of the ejecta cloud in the framework of the current mission scenario of AIDA mission project. A grid search of launching sites of ejecta was defined over the globe of Didymoon, and considering a wide range of possible ejection speeds, we determined the dependency of ejecta fate on launching sites (projectile impact sites) and speeds. This range allows us to track all the complex cases that include different types of dynamical fates. Two major mechanisms are found to be working broadly during the post-ejection evolution of the ejecta cloud: 1) ejecta on mean motion resonance orbits with Didymoon produce long-term quasi-periodic showers onto Didymoon over at least a couple of weeks after the projectile impact, 2) ejecta on non-resonant orbits produce a rapid and high re-accretion flux. This rapid and high flux occurs just once because ejecta on such orbits leave the system unless they experience a collision during their first encounter. For the second part of this study, we performed full-scale simulations of the ejecta cloud released from 6 hypothetical impact sites. We considered two kinds of material composing Didymoon's subsurface and then combined a power-law size distribution of the ejecta with an ejection speed distribution. We find that the ejecta cloud evolution can be divided in two periods. It starts with a first violent period (<10 hr) with fast re-accretion or ejection from the system. A second period is found to be more sensitive to the launching site than the first one. During this second period, ejecta will either re-accrete or being ejected from the system, depending both on their sizes and on their average survival time in close proximity of the binary components.

As the target of the proposed Asteroid Impact & Deflection Assessment (AIDA) mission, the near-Earth binary asteroid 65803 Didymos represents a special class of binary asteroids, those whose primaries are at risk of rotational disruption. To gain a better understanding of these binary systems and to support the AIDA mission, this paper investigates the creep stability of the Didymos primary by representing it as a cohesionless self-gravitating granular aggregate subject to rotational acceleration. To achieve this goal, a soft-sphere discrete element model (SSDEM) capable of simulating granular systems in quasi-static states is implemented and a quasi-static spin-up procedure is carried out. We devise three critical spin limits for the simulated aggregates to indicate their critical states triggered by reshaping and surface shedding, internal structural deformation, and shear failure, respectively. The failure condition and mode, and shear strength of an aggregate can all be inferred from the three critical spin limits. The effects of arrangement and size distribution of constituent particles, bulk density, spin-up path, and interparticle friction are numerically explored. The results show that the shear strength of a spinning self-gravitating aggregate depends strongly on both its internal configuration and material parameters, while its failure mode and mechanism are mainly affected by its internal configuration. Additionally, this study provides some constraints on the possible physical properties of the Didymos primary based on observational data and proposes a plausible formation mechanism for this binary system. With a bulk density consistent with observational uncertainty and close to the maximum density allowed for the asteroid, the Didymos primary in certain configurations can remain geo-statically stable without including cohesion.

The aim of the Asteroid Science Intersections with In-Space Mine Engineering (ASIME) 2016 conference on September 21-22, 2016 in Luxembourg City was to provide an environment for the detailed discussion of the specific properties of asteroids, with the engineering needs of space missions that utilize asteroids. The ASIME 2016 Conference produced a layered record of discussions from the asteroid scientists and the asteroid miners to understand each other's key concerns and to address key scientific questions from the asteroid mining companies: Planetary Resources, Deep Space Industries and TransAstra. These Questions were the focus of the two day conference, were addressed by scientists inside and outside of the ASIME Conference and are the focus of this White Paper. The Questions from the asteroid mining companies have been sorted into the three asteroid science themes: 1) survey, 2) surface and 3) subsurface and 4) Other. The answers to those Questions have been provided by the scientists with their conference presentations or edited directly into an early open-access collaborative Google document (August 2016-October 2016), or inserted by A. Graps using additional reference materials. During the ASIME 2016 last two-hours, the scientists turned the Questions from the Asteroid Miners around by presenting their own key concerns: Questions from the Asteroid Scientists . These answers in this White Paper will point to the Science Knowledge Gaps (SKGs) for advancing the asteroid in-space resource utilisation domain.

The Chelyabinsk meteorite is a highly shocked, low porosity, ordinary chondrite, probably similar to S- or Q-type asteroids. Therefore, nanoindentation experiments on this meteorite allow us to obtain key data to understand the physical properties of near-Earth asteroids. Tests at different length scales provide information about the local mechanical properties of the minerals forming this meteorite: reduced Young's modulus, hardness, elastic recovery, and fracture toughness. Those tests are also useful to understand the potential to deflect threatening asteroids using a kinetic projectile. We found that the differences in mechanical properties between regions of the meteorite, which increase or reduce the efficiency of impacts, are not a result of compositional differences. A low mean particle size, attributed to repetitive shock, can increase hardness, while low porosity promotes a higher momentum multiplication. Momentum multiplication is the ratio between the change in momentum of a target due to an impact, and the momentum of the projectile, and, therefore higher values imply more efficient impacts. In the Chelyabinsk meteorite, the properties of the light-colored lithology materials facilitate obtaining higher momentum multiplication values, compared to the other regions described for this meteorite. Also, we found a low value of fracture toughness in the shock-melt veins of Chelyabinsk, which would promote the ejection of material after an impact and, therefore, increase the momentum multiplication. These results are relevant considering the growing interest in missions to test asteroid deflection, such as the recent collaboration between the European Space Agency and NASA, known as the Asteroid Impact and Deflection Assessment mission.

An understanding of the post-impact dynamics of ejecta clouds are crucial to the planning of a kinetic impact mission to an asteroid, and also has great implications for the history of planetary formation. The purpose of this article to track the evolution of ejecta produced by AIDA mission, which targets for kinetic impact the secondary of near-Earth binary asteroid 65803 Didymos on 2022, and to feedback essential informations to AIDA's ongoing phase-A study. We present a detailed dynamic model for the simulation of an ejecta cloud from a binary asteroid that synthesizes all relevant forces based on a previous analysis of the mechanical environment. We apply our method to gain insight into the expected response of Didymos to the AIDA impact, including the subsequent evolution of debris and dust. The crater scaling relations from laboratory experiments are employed to approximate the distributions of ejecta mass and launching speed. The size composition of fragments is modeled with a power law fitted from observations of real asteroid surface. A full-scale demonstration is simulated using parameters specified by the mission. We report the results of the simulation, which include the computed spread of the ejecta cloud and the recorded history of ejecta accretion and escape. The violent period of the ejecta evolution is found to be short, and is followed by a stage where the remaining ejecta is gradually cleared. Solar radiation pressure proves to be efficient in cleaning dust-size ejecta, and the simulation results after two weeks shows that large debris on polar orbits (perpendicular to the binary orbital plane) has a survival advantage over smaller ejecta and ejecta that keep to low latitudes.

Deflection missions to near-Earth asteroids will encounter non-negligible uncertainties in the physical and orbital parameters of the target object. In order to reliably assess future impact threat mitigation operations such uncertainties have to be quantified and incorporated into the mission design. The implementation of deflection demonstration missions offers the great opportunity to test our current understanding of deflection relevant uncertainties and their consequences, e.g., regarding kinetic impacts on asteroid surfaces. In this contribution, we discuss the role of uncertainties in the NEOTwIST asteroid deflection demonstration concept, a low-cost kinetic impactor design elaborated in the framework of the NEOShield project. The aim of NEOTwIST is to change the spin state of a known and well characterized near-Earth object, in this case the asteroid (25143) Itokawa. Fast events such as the production of the impact crater and ejecta are studied via cube-sat chasers and a flyby vehicle. Long term changes, for instance, in the asteroid's spin and orbit, can be assessed using ground based observations. We find that such a mission can indeed provide valuable constraints on mitigation relevant parameters. Furthermore, the here proposed kinetic impact scenarios can be implemented within the next two decades without threatening Earth's safety.

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.

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.

Our knowledge of the strengths of small bodies in the Solar System is limited by our poor understanding of their internal structures, and this, in turn, clouds our understanding of the formation and evolution of these bodies. Observations of the rotational states of asteroids whose diameters are larger than a few hundreds of meters have revealed that they are dominated by gravity and that most are unlikely to be monoliths; however, there is a wide range of plausible internal structures. Numerical and analytical studies of shape and spin limits of gravitational aggregates and their collisional evolution show a strong dependence on shear strength. In order to study this effect, we carry out a systematic exploration of the dependence of collision outcomes on dissipation and friction parameters of the material components making up the bodies. We simulate the catastrophic disruption (leading to the largest remnant retaining 50% of the original mass) of km-size asteroids modeled as gravitational aggregates using pkdgrav, a cosmology N-body code adapted to collisional problems and recently enhanced with a new soft-sphere collision algorithm that includes more realistic contact forces. We find that for a range of three different materials, higher friction and dissipation values increase the catastrophic disruption threshold by about half a magnitude. Furthermore, we find that pre-impact rotation systematically increases mass loss on average, regardless of the target's internal configuration. Our results have important implications for the efficiency of planet formation via planetesimal growth, and also more generally to estimate the impact energy threshold for catastrophic disruption, as this generally has only been evaluated for non-spinning bodies without detailed consideration of material properties.

The Design Reference Asteroid (DRA) is a compilation of all that is known about the OSIRIS-REx mission target, asteroid (101955) Bennu. It contains our best knowledge of the properties of Bennu based on an extensive observational campaign that began shortly after its discovery, and has been used to inform mission plan development and flight system design. The DRA will also be compared with post-encounter science results to determine the accuracy of our Earth-based characterization efforts. The extensive observations of Bennu in 1999 has made it one of the best-characterized near-Earth asteroids. Many physical parameters are well determined, and span a number of categories: Orbital, Bulk, Rotational, Radar, Photometric, Spectroscopic, Thermal, Surface Analog, and Environment Properties. Some results described in the DRA have been published in peer-reviewed journals while others have been reviewed by OSIRIS-REx Science Team members and/or external reviewers. Some data, such as Surface Analog Properties, are based on our best knowledge of asteroid surfaces, in particular those of asteroids Eros and Itokawa. This public release of the OSIRIS-REx Design Reference Asteroid is a annotated version of the internal OSIRIS-REx document OREX-DOCS-04.00-00002, Rev 9 (accepted by the OSIRIS-REx project on 2014-April-14). The supplemental data products that accompany the official OSIRIS-REx version of the DRA are not included in this release. We are making this document available as a service to future mission planners in the hope that it will inform their efforts.

This study is carried out in the framework of sample-return missions to asteroids that use a low-speed projectile as the primary component of its sampling mechanism (e.g., JAXA's Hayabusa and Hayabusa2 missions). We perform numerical simulations of such impacts into granular materials using different projectile shapes under Earth's gravity. We then compare the amounts of ejected mass obtained in our simulations against what was found in experiments that used similar setups, which allows us to validate our numerical approach. For the targets, we consider 2 different monodisperse grain-diameter sizes: 5 mm and 3 mm. The impact speed of the projectile is 11 m s$^{-1}$ directed downward, perpendicular to the surface of the targets. Using an implementation of the soft-sphere discrete element method (SSDEM) in the $N$-Body gravity tree code PKDGRAV, previously validated in the context of low-speed impacts into sintered glass bead agglomerates, we find a noticeable dependence of the amount of ejected mass on the projectile shape. As found in experiments, in the case of the larger target grain size (5 mm), a conically shaped projectile ejects a greater amount of mass than do projectiles of other shapes, including disks and spheres. We then find that numerically the results are sensitive to the normal coefficient of restitution of the grains, especially for impacts into targets comprised of smaller grains (3 mm). We also find that static friction plays a more important role for impacts into targets comprised of the larger grains. As a preliminary demonstration, one of these considered setups is simulated in a microgravity environment. As expected, a reduction in gravity increases both the amount of ejected mass and the timescale of the impact process. Our methodology is also adaptable to the conditions of sampling mechanisms included in specific mission designs.

Asteroid (99942) Apophis' close approach in 2029 will be one of the most significant small-body encounter events in the near future and offers a good opportunity for in situ exploration to determine the asteroid's surface properties and measure any tidal effects that might alter its regolith configuration. Resurfacing mechanics has become a new focus for asteroid researchers due to its important implications for interpreting surface observations, including space weathering effects. This paper provides a prediction for the tidal effects during the 2029 encounter, with an emphasis on whether surface refreshing due to regolith movement will occur. The potential shape modification of the object due to the tidal encounter is first confirmed to be negligibly small with systematic simulations, thus only the external perturbations are taken into account for this work (despite this, seismic shaking induced by shifting blocks might still play a weak role and we will look into this mechanism in future work). A two-stage approach is developed to model the responses of asteroid surface particles (the regolith) based on the soft-sphere implementation of the parallel N-body gravity tree code pkdgrav. A full-body model of Apophis is sent past the Earth on the predicted trajectory to generate the data of all forces acting at a target point on the surface. A sandpile constructed in the local frame is then used to approximate the regolith materials; all the forces the sandpile feels during the encounter are imposed as external perturbations to mimic the regolith's behavior in the full scenario. The local mechanical environment on the asteroid surface is represented in detail, leading to an estimation of the change in global surface environment due to the encounter. Typical patterns of perturbation are presented that depend on the asteroid orientation and sense of rotation at perigee.

Out of the handful of asteroids that have been imaged, some have distributions of blocks that are not easily explained. In this paper, we investigate the possibility that seismic shaking leads to the size sorting of particles in asteroids. In particular, we focus on the so-called Brazil Nut Effect (BNE) that separates large particles from small ones under vibrations. We study the BNE over a wide range of parameters by using the N-body code PKDGRAV, and find that the effect is largely insensitive to the coefficients of restitution, but sensitive to friction constants and oscillation speeds. Agreeing with the previous results, we find that convection drives the BNE, where the intruder rises to the top of the particle bed. For the wide-cylinder case, we also observe a "whale" effect, where the intruder follows the convective current and does not stay at the surface. We show that the non-dimensional critical conditions for the BNE agree well with previous studies. We also show that the BNE is scalable for low-gravity environments and that the rise speed of an intruder is proportional to the square root of the gravitational acceleration. Finally, we apply the critical conditions to observed asteroids, and find that the critical oscillation speeds are comparable to the seismic oscillation speeds that are expected from non-destructive impacts.

We carry out a systematic exploration of the effect of pre-impact rotation on the outcomes of low-speed collisions between planetesimals modeled as gravitational aggregates. We use pkdgrav, a cosmology code adapted to collisional problems and recently enhanced with a new soft-sphere collision algorithm that includes more realistic contact forces. A rotating body has lower effective surface gravity than a non-rotating one and therefore might suffer more mass loss as the result of a collision. What is less well understood, however, is whether rotation systematically increases mass loss on average regardless of the impact trajectory. This has important implications for the efficiency of planet formation via planetesimal growth, and also more generally for the determination of the impact energy threshold for catastrophic disruption (leading to the largest remnant retaining 50% of the original mass), as this has generally only been evaluated for non-spinning bodies. We find that for most collision scenarios, rotation lowers the threshold energy for catastrophic dispersal. For head-on collisions, we develop a semi-analytic description of the change in the threshold description as a function of the target's pre-impact rotation rate, and find that these results are consistent with the "universal law" of catastrophic disruption developed by Leinhardt & Stewart. Using this approach, we introduce re-scaled catastrophic disruption variables that take into account the interacting mass fraction of the target and the projectile in order to translate oblique impacts into equivalent head-on collisions.

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular "atmosphere" that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.

Despite their very low surface gravities, asteroids exhibit a number of different geological processes involving granular matter. Understanding the response of this granular material subject to external forces in microgravity conditions is vital to the design of a successful asteroid sub-surface sampling mechanism, and in the interpretation of the fascinating geology on an asteroid. We have designed and flown a Taylor-Couette shear cell to investigate granular flow due to rotational shear forces under the conditions of parabolic flight microgravity. The experiments occur under weak compression. First, we present the technical details of the experimental design with particular emphasis on how the equipment has been specifically designed for the parabolic flight environment. Then, we investigate how a steady state granular flow induced by rotational shear forces differs in varying gravitational environments. We find that the effect of constant shearing on the granular material, in a direction perpendicular to the effective acceleration, does not seem to be strongly influenced by gravity. This means that shear bands can form in the presence of a weak gravitational field just as on Earth.

We investigate the role of gravity on convection in a dense granular shear flow. Using a microgravity modified Taylor-Couette shear cell under the conditions of parabolic flight microgravity, we demonstrate experimentally that secondary, convective-like flows in a sheared granular material are close to zero in microgravity and enhanced under high-gravity conditions, though the primary flow fields are unaffected by gravity. We suggest that gravity tunes the frictional particle-particle and particle-wall interactions, which have been proposed to drive the secondary flow. In addition, the degree of plastic deformation increases with increasing gravitational forces, supporting the notion that friction is the ultimate cause.

Surfaces of planets and small bodies of our Solar System are often covered by a layer of granular material that can range from a fine regolith to a gravel-like structure of varying depths. Therefore, the dynamics of granular materials are involved in many events occurring during planetary and small-body evolution thus contributing to their geological properties. We demonstrate that the new adaptation of the parallel N-body hard-sphere code pkdgrav has the capability to model accurately the key features of the collective motion of bidisperse granular materials in a dense regime as a result of shaking. As a stringent test of the numerical code we investigate the complex collective ordering and motion of granular material by direct comparison with laboratory experiments. We demonstrate that, as experimentally observed, the scale of the collective motion increases with increasing small-particle additive concentration. We then extend our investigations to assess how self-gravity and external gravity affect collective motion. In our reduced-gravity simulations both the gravitational conditions and the frequency of the vibrations roughly match the conditions on asteroids subjected to seismic shaking, though real regolith is likely to be much more heterogeneous and less ordered than in our idealised simulations. We also show that collective motion can occur in a granular material under a wide range of inter-particle gravity conditions and in the absence of an external gravitational field. These investigations demonstrate the great interest of being able to simulate conditions that are to relevant planetary science yet unreachable by Earth-based laboratory experiments.

Despite their very low surface gravities, the surfaces of asteroids and comets are covered by granular materials - regolith - that can range from a fine dust to a gravel-like structure of varying depths. Understanding the dynamics of granular materials is, therefore, vital for the interpretation of the surface geology of these small bodies and is also critical for the design and/or operations of any device planned to interact with their surfaces. We present the first measurements of transient weakening of granular material after shear reversal in microgravity as well as a summary of experimental results recently published in other journals, which may have important implications for small-body surfaces. Our results suggest that the force contact network within a granular material may be weaker in microgravity, although the influence of any change in the contact network is felt by the granular material over much larger distances. This could mean that small body surfaces are even more unstable than previously imagined. However, our results also indicate that the consequences of, e.g., a meteorite impact or a spacecraft landing, may be very different depending on the impact angle and location, and depending on the prior history of the small body surface.

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.

Asteroids are leftover pieces from the era of planet formation that help us understand conditions in the early Solar System. Unlike larger planetary bodies that were subject to global thermal modification during and subsequent to their formation, these small bodies have kept at least some unmodified primordial material from the solar nebula. However, the structural properties of asteroids have been modified considerably since their formation. Thus, we can find among them a great variety of physical configurations and dynamical histories. In fact, with only a few possible exceptions, all asteroids have been modified or completely disrupted many times during the age of the Solar System. This picture is supported by data from space mission encounters with asteroids that show much diversity of shape, bulk density, surface morphology, and other features. Moreover, the gravitational attraction of these bodies is so small that some physical processes occur in a manner far removed from our common experience on Earth. Thus, each visit to a small body has generated as many questions as it has answered. In this review we discuss the current state of research into asteroid disruption processes, focusing on collisional and rotational mechanisms. We find that recent advances in modeling catastrophic disruption by collisions have provided important insights into asteroid internal structures and a deeper understanding of asteroid families. Rotational disruption, by tidal encounters or thermal effects, is responsible for altering many smaller asteroids, and is at the origin of many binary asteroids and oddly shaped bodies.

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.

For a long time, gravitational instability in the disk of planetesimals has been suspected to be the main engine responsible for the beginning of dust growth, its advantage being that it provides for rapid growth. Its real importance in planetary formation is still debated, mainly because the potential presence of turbulence can prevent the settling of particles into a gravitationally unstable layer. However, several mechanisms could yield strongly inhomogeneous distributions of solids in the disk: radial drift, trapping in vortices, perturbations by other massive bodies, etc. In this paper we present a numerical study of a gravitationally unstable layer. This allows us to go beyond the classical analytical study of linear perturbations, exploring a highly non-linear regime. A hierarchical growth of structure in the presence of dissipation (gas drag) can yield large, virialized clusters of planetesimals, the first time such clusters have been observed in the context of planetesimal disks.

Until recently, most asteroids were thought to be solid bodies whose shapes were determined largely by collisions with other asteroids. It now seems that many asteroids are little more than rubble piles, held together by self-gravity; this means that their shapes may be strongly distorted by tides during close encounters with planets. Here we report on numerical simulations of encounters between a ellipsoid-shaped rubble-pile asteroid and the Earth. After an encounter, many of the simulated asteroids develop the same rotation rate and distinctive shape (i.e., highly elongated with a single convex side, tapered ends, and small protuberances swept back against the rotation direction) as 1620 Geographos. Since our numerical studies show that these events occur with some frequency, we suggest that Geographos may be a tidally distorted object. In addition, our work shows that 433 Eros, which will be visited by the NEAR spacecraft in 1999, is much like Geographos, which suggests that it too may have been molded by tides in the past.