Search

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 file

While no known asteroid poses a threat to Earth for at least the next century, the catalog of near-Earth asteroids is incomplete for objects whose impacts would produce regional devastation. Several approaches have been proposed to potentially prevent an asteroid impact with Earth by deflecting or disrupting an asteroid. A test of kinetic impact technology was identified as the highest priority space mission related to asteroid mitigation. NASA's Double Asteroid Redirection Test (DART) mission is the first full-scale test of kinetic impact technology. The mission's target asteroid was Dimorphos, the secondary member of the S-type binary near-Earth asteroid (65803) Didymos. This binary asteroid system was chosen to enable ground-based telescopes to quantify the asteroid deflection caused by DART's impact. While past missions have utilized impactors to investigate the properties of small bodies those earlier missions were not intended to deflect their targets and did not achieve measurable deflections. Here we report the DART spacecraft's autonomous kinetic impact into Dimorphos and reconstruct the impact event, including the timeline leading to impact, the location and nature of the DART impact site, and the size and shape of Dimorphos. The successful impact of the DART spacecraft with Dimorphos and the resulting change in Dimorphos's orbit demonstrates that kinetic impactor technology is a viable technique to potentially defend Earth if necessary., Comment: Accepted by Nature

Some active asteroids have been proposed to be the result of impact events. Because active asteroids are generally discovered serendipitously only after their tail formation, the process of the impact ejecta evolving into a tail has never been directly observed. NASA's Double Asteroid Redirection Test (DART) mission, apart from having successfully changed the orbital period of Dimorphos, demonstrated the activation process of an asteroid from an impact under precisely known impact conditions. Here we report the observations of the DART impact ejecta with the Hubble Space Telescope (HST) from impact time T+15 minutes to T+18.5 days at spatial resolutions of ~2.1 km per pixel. Our observations reveal a complex evolution of ejecta, which is first dominated by the gravitational interaction between the Didymos binary system and the ejected dust and later by solar radiation pressure. The lowest-speed ejecta dispersed via a sustained tail that displayed a consistent morphology with previously observed asteroid tails thought to be produced by impact. The ejecta evolution following DART's controlled impact experiment thus provides a framework for understanding the fundamental mechanisms acting on asteroids disrupted by natural impact., Comment: accepted by Nature

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 $\beta$, showing that a particular direction-specific $\beta$ will be directly determined by the DART results, and that a related direction-specific $\beta$ is a figure of merit for a kinetic impact mission. The DART $\beta$ 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., Comment: 30 pages, 7 figures. Planetary Science Journal, in press, accepted 2022 September 22

The Double Asteroid Redirection Test (DART) spacecraft will impact into the asteroid Dimorphos on September 26, 2022 as a test of the kinetic impactor technique for planetary defense. The efficiency of the deflection following a kinetic impactor can be represented using the momentum enhancement factor, Beta, which is dependent on factors such as impact geometry and the specific target material properties. Currently, very little is known about Dimorphos and its material properties that introduces uncertainty in the results of the deflection efficiency observables, including crater formation, ejecta distribution, and Beta. The DART Impact Modeling Working Group (IWG) is responsible for using impact simulations to better understand the results of the DART impact. Pre-impact simulation studies also provide considerable insight into how different properties and impact scenarios affect momentum enhancement following a kinetic impact. This insight provides a basis for predicting the effects of the DART impact and the first understanding of how to interpret results following the encounter. Following the DART impact, the knowledge gained from these studies will inform the initial simulations that will recreate the impact conditions, including providing estimates for potential material properties of Dimorphos and Beta resulting from DARTs impact. This paper summarizes, at a high level, what has been learned from the IWG simulations and experiments in preparation for the DART impact. While unknown, estimates for reasonable potential material properties of Dimorphos provide predictions for Beta of 1-5, depending on end-member cases in the strength regime., Comment: Accepted to PSJ Didymos-DART Focus Issue

NASA's Double Asteroid Redirection Test (DART) spacecraft is planned to impact the natural satellite of (65803) Didymos, Dimorphos, around 23:14 UTC on 26 September 2022, causing a reduction in its orbital period that will be measurable with ground-based observations. This test of kinetic impactor technology will provide the first estimate of the momentum transfer enhancement factor $\beta$ at a realistic scale, wherein ejecta from the impact provides an additional deflection to the target. Earth-based observations, the LICIACube spacecraft (to be detached from DART prior to impact), and ESA's follow-up Hera mission to launch in 2024, will provide additional characterization of the deflection test. Together Hera and DART comprise the Asteroid Impact and Deflection Assessment (AIDA) cooperation between NASA and ESA. Here the predicted dynamical states of the binary system upon arrival and after impact are presented. The assumed dynamically relaxed state of the system will be excited by the impact, leading to an increase in eccentricity and slight tilt of the orbit together with enhanced libration of Dimorphos with amplitude dependent on the currently poorly known target shape. Free rotation around the moon's long axis may also be triggered and the orbital period will experience variations from seconds to minutes over timescales of days to months. Shape change of either body due to cratering or mass wasting triggered by crater formation and ejecta may affect $\beta$ but can be constrained through additional measurements. Both BYORP and gravity tides may cause measurable orbital changes on the timescale of Hera's rendezvous., Comment: 23 pages, 13 figures, published in PSJ

In the canonical model of Moon formation, a Mars-sized protoplanet "Theia" collides with proto-Earth at close to their mutual escape velocity $v_{\rm esc}$ and a common impact angle 45{\deg}. The "graze-and-merge" collision strands a fraction of Theia's mantle into orbit, while Earth accretes most of Theia and its momentum. Simulations show that this produces a hot, high angular momentum, silicate-dominated protolunar system, in substantial agreement with lunar geology, geochemistry, and dynamics. However, a Moon that derives mostly from Theia's mantle, as angular momentum dictates, is challenged by the fact that O, Ti, Cr, radiogenic W, and other elements are indistinguishable in Earth and lunar rocks. Moreover, the model requires an improbably low initial velocity. Here we develop a scenario for Moon formation that begins with a somewhat faster collision, when proto-Theia impacts proto-Earth at ~1.2 $v_{\rm esc}$, also around 45{\deg}. Instead of merging, the bodies come into violent contact for a half-hour and their major components escape, a "hit-and-run collision." N-body evolutions show that the "runner" often returns ~0.1-1 Myr later for a second giant impact, closer to $v_{\rm esc}$; this produces a postimpact disk of ~2-3 lunar masses in smoothed particle hydrodynamics simulations, with angular momentum comparable to canonical scenarios. The disk ends up substantially inclined, in most cases, because the terminal collision is randomly oriented to the first. Proto-Earth contributions to the silicate disk are enhanced by the compounded mixing and greater energy of a collision chain., Comment: 27 pages, 4 figures, 2 tables. Published in PSJ, article available at https://iopscience.iop.org/article/10.3847/PSJ/ac19b2

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 $\gtrsim$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., Comment: 18 pages, 9 figures, 4 tables. Published in PSJ, article available at https://iopscience.iop.org/article/10.3847/PSJ/ac19b1

We compare low velocity impacts that ricochet with the same impact velocity and impact angle into granular media with similar bulk density, porosity and friction coefficient but different mean grain size. The ratio of projectile diameter to mean grain length ranges from 4 in our coarsest medium to 50 in our finest sand. Using high speed video and fluorescent markers, we measure the ratio of pre- to post-impact horizontal and vertical velocity components, which we refer to as coefficients of restitution, and the angle of deflection caused by the impact in the horizontal plane. Coefficients of restitution are sensitive to mean grain size with the ratio associated with the horizontal velocity component about twice as large for our coarsest gravel as that for our finest sand. This implies that coefficients for hydro-static-like, drag-like and lift-like forces, used in empirical force laws, are sensitive to mean grain size. The coefficient that is most strongly sensitive to grain size is the lift coefficient which decreases by a factor of 3 between our coarsest and finest media. The deflection angles are largest in the coarser media and their size approximately depends on grain size to the 3/2 power. This scaling is matched with a model where momentum transfer takes place via collisions with individual grains. The dependence of impact mechanics on substrate size distribution should be considered in future models for populations of objects that impact granular asteroid surfaces., Comment: Link to the high speed videos for the 4 experiments tracked and discussed: https://youtube.com/playlist?list=PLYnpHw29go3pCWwm554fwR6jNFvB5qVQn Submitted to the journal Icarus on September 12, 2021

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/abf0ad

The Discrete Element Method (DEM) is frequently used to model complex granular systems and to augment the knowledge that we obtain through theory, experimentation, and real-world observations. Numerical simulations are a particularly powerful tool for studying the regolith-covered surfaces of asteroids, comets, and small moons, where reduced-gravity environments produce ill-defined flow behaviors. In this work, we present a method for validating soft-sphere DEM codes for both terrestrial and small-body granular environments. The open-source code Chrono is modified and evaluated first with a series of simple two-body-collision tests, and then, with a set of piling and tumbler tests. In the piling tests, we vary the coefficient of rolling friction to calibrate the simulations against experiments with 1 mm glass beads. Then, we use the friction coefficient to model the flow of 1 mm glass beads in a rotating drum, using a drum configuration from a previous experimental study. We measure the dynamic angle of repose, the flowing layer thickness, and the flowing layer velocity for tests with different particle sizes, contact force models, coefficients of rolling friction, cohesion levels, drum rotation speeds and gravity levels. The tests show that the same flow patterns can be observed at Earth and reduced-gravity levels if the drum rotation speed and the gravity-level are set according to the dimensionless parameter known as the Froude number. Chrono is successfully validated against known flow behaviors at different gravity and cohesion levels, and will be used to study small-body regolith dynamics in future works.

Spin off events and impacts can eject boulders from an asteroid surface and rubble pile asteroids can accumulate from debris following a collision between large asteroids. These processes produce a population of gravitational bound objects in orbit that can impact an asteroid surface at low velocity and with a distribution of impact angles. We present laboratory experiments of low velocity spherical projectiles into a fine granular medium, sand. We delineate velocity and impact angles giving ricochets, those giving projectiles that roll-out from the impact crater and those that stop within their impact crater. With high speed camera images and fluorescent markers on the projectiles we track spin and projectile trajectories during impact. We find that the projectile only reaches a rolling without slipping condition well after the marble has reached peak penetration depth. The required friction coefficient during the penetration phase of impact is 4-5 times lower than that of the sand suggesting that the sand is fluidized near the projectile surface during penetration. We find that the critical grazing impact critical angle dividing ricochets from roll-outs, increases with increasing impact velocity. The critical angles for ricochet and for roll-out as a function of velocity can be matched by an empirical model during the rebound phase that balances a lift force against gravity. We estimate constraints on projectile radius, velocity and impact angle that would allow projectiles on asteroids to ricochet or roll away from impact, finally coming to rest distant from their initial impact sites., Comment: Link to the high speed videos for the 3 experiments tracked and discussed: https://www.youtube.com/playlist?list=PLYnpHw29go3rKYOvM4WEQkHWYYIGgRC-E Preliminary results presented at American Geophysical Union (AGU) 2019 Fall Meeting (P43D-3497). Submitted to the journal Icarus on February 3, 2020

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/collresolve

Near-Earth small-body mission targets 162173 Ryugu, 101955 Bennu, and 25143 Itokawa produce gravity fields around 4 orders of magnitude below that of Earth and their irregular shapes, combined with rotational effects produce varying surface potentials. Still, we observe familiar geologic textures and landforms that are the result of the granular physical processes that take place on their surfaces. The nature of these landforms, however, their origins, and how these surfaces react to interrogation by probes, landers, rovers, and penetrators, remain largely unknown, and therefore landing on an asteroid and manipulating its surface material remains a daunting challenge. The AOSAT+ design is a 12U CubeSat that will be in Low-Earth Orbit (LEO) and that will operate as a spinning on-orbit centrifuge. Part of this 12U CubeSat will contain a laboratory that will recreate asteroid surface conditions using crushed meteorite as a regolith proxy. The spinning of the laboratory will simulate the surface gravity of asteroids 2 km and smaller. The result is a bed of realistic regolith, the environment that landers and diggers and maybe astronauts will interact with. A crucial component of this mission involves the reproduction of the experimental results in numerical simulation in order to extract the material parameters of the regolith and its behavior in a sustained, very low but nonzero-gravity environment., Comment: 7 pages, 2 figures, International Astronautical Congress 2019

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., Comment: 8 pages, 5 figures, published in the Astrophysical Journal Letters, 874, 1, L10

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 such investigations allow us to better understand the final configuration and geochemistry of the terrestrial planets, as well as to place our solar system in the context of other exosolar systems. These simulations, however, classically prescribe collisions to result in perfect mergers, but 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, streamlining simulations into a classifier of collision outcomes and a regressor of accretion efficiency. The classifier maps a 4-Dimensional 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 provides insights about the outcome regimes. The regressor is a neural network which is trained to closely mimic the functional relationship between the 4-D 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, 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., Comment: 18 pages, 7 figures, 3 tables. Accepted for publication on ApJ

We consider an impact on an asteroid that is energetic enough to cause resurfacing by seismic reverberation and just below the catastrophic disruption threshold, assuming that seismic waves are not rapidly attenuated. In asteroids with diameter less than 1 km we identify a regime where rare energetic impactors can excite seismic waves with frequencies near those of the asteroid's slowest normal modes. In this regime, the distribution of seismic reverberation is not evenly distributed across the body surface. With mass-spring model elastic simulations, we model impact excitation of seismic waves with a force pulse exerted on the surface and using three different asteroid shape models. The simulations exhibit antipodal focusing and normal mode excitation. If the impulse excited vibrational energy is long lasting, vibrations are highest at impact point, its antipode and at high surface elevations such as an equatorial ridge. A near equatorial impact launches a seismic impulse on a non-spherical body that can be focused on two additional points on an the equatorial ridge. We explore simple flow models for the morphology of vibration induced surface slumping. We find that the initial seismic pulse is unlikely to cause large shape changes. Long lasting seismic reverberation on Bennu caused by a near equatorial impact could have raised the height of its equatorial ridge by a few meters and raised two peaks on it, one near impact site and the other near its antipode., Comment: submitted to Icarus

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., Comment: 8 pages, 7 figures, 2 tables, accepted for publication in MNRAS

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., Comment: 66 pages, 24 figures, submitted to Icarus on 25/Aug/2016

CubeSats and small satellites are emerging as low-cost tools for performing science and exploration in deep space. These new classes of satellite exploit the latest advancement in miniaturization of electronics, power systems, and communication technologies to promise reduced launch cost and development cadence. JPL's MarCO CubeSats, part of the Mars Insight mission, will head on an Earth escape trajectory to Mars in 2018 and serve as communication relays for the Mars Insight Lander during Entry, Descent and Landing. Incremental advancements to the MarCO CubeSats, particularly in propulsion and GNC, could enable these spacecraft to get to another planet or to Near Earth Objects. This can have substantial science return with the right science instrument. We have developed an interplanetary CubeSat concept that includes onboard green monopropellant propulsion system and that can get into a capture orbit around a neighboring planet or chase a small-body. One such candidate is the Martian moon Phobos. Because of the limits of current CubeSat hardware and lack of an accurate ephemeris of Phobos, there will be a 2 to 5 km uncertainty in distance between the spacecraft and Phobos. This presents a major GNC challenge when the CubeSat first attempts to get into visual range of the moon. One solution to this challenge is to develop optical navigation technology that enables the CubeSat to take epicyclic orbits around the most probable location of the target, autonomously search and home-in on the target body. In worst-case scenarios, the technology would narrow down the uncertainty of the small-body location and then use optical flow, a computer vision algorithm to track movement of objects in the field of view. A dimly lit small-body can be detected by the occlusion of one or more surrounding stars. Our studies present preliminary simulations that support the concept., Comment: in Proceedings of the 40th Annual AAS Guidance, Navigation and Control Conference 2017

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., Comment: 50 pages, 9 figures

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., Comment: 19 pages, 11 figures, accepted for publication in Advances in Space Research, Special Issue: Asteroids & Space Debris

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., Comment: Accepted for publication in the proceedings of the IAUS 318 - Asteroids: New Observations, New Models, held at the IAU General Assembly in Honolulu, Hawaii, USA 2015

The regolith-covered surfaces of asteroids preserve records of geophysical processes that have occurred both at their surfaces and sometimes also in their interiors. As a result of the unique micro-gravity environment that these bodies posses, a complex and varied geophysics has given birth to fascinating features that we are just now beginning to understand. The processes that formed such features were first hypothesised through detailed spacecraft observations and have been further studied using theoretical, numerical and experimental methods that often combine several scientific disciplines. These multiple approaches are now merging towards a further understanding of the geophysical states of the surfaces of asteroids. In this chapter we provide a concise summary of what the scientific community has learned so far about the surfaces of these small planetary bodies and the processes that have shaped them. We also discuss the state of the art in terms of experimental techniques and numerical simulations that are currently being used to investigate regolith processes occurring on small-body surfaces and that are contributing to the interpretation of observations and the design of future space missions., Comment: 28 pages, 8 figures. Chapter to appear in the book ASTEROIDS IV, (University of Arizona Press) Space Science Series, edited by P. Michel, F. DeMeo and W. Bottke

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., Comment: 27 pages, 6 figures, Accepted for publication in Planetary and Space Science. arXiv admin note: text overlap with arXiv:1406.5228

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., Comment: 15 pages, 1 table, 14 figures, accepted for publication in MNRAS

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., Comment: 13 pages, 11 figures, 2 tables, Accepted for publication in Astrophysical Journal