This investigation addresses the discovery of a proposed impact crater located off the coast of Miami, FL under the North Atlantic Ocean. A preliminary analysis of bathymetry data obtained from the National Oceanic and Atmospheric Administration implied a morphology consistent with a complex crater produced by a hypervelocity impact event of extraterrestrial origin. The proposed impact features include a central peak, concentric rings, and an ejecta field to the northwest. Analysis of geological data from the US Geological Survey places the strata overlying the proposed impact site as Miami Limestone, Pleistocene, accumulated during Marine Isotope Stage 5e, thereby placing the maximum age of the proposed impact crater at 80 ka to 130 ka. Three other competing hypotheses for the formation of the structure, namely a controlled maritime explosion, radial lava flow from volcano, or a depressed bioherm, doline, or karst (solutional depression) were explored throughout the investigation. To confirm the proposed structure as an impact crater, an in-situ underwater expedition was organized by Planetary Sciences, Inc. specifically to ascertain whether planar formations, shatter cones, and shock metamorphic and or other meteoritic properties were present. After analyzing the geological samples collected at the proposed impact crater, examining the morphology of analogous geologic structures, and evaluating competing hypotheses, we conclude that the structure is a solutional doline formed by the uneven dissolution of the Miami Limestone, and, accordingly, do not recommend that the structure be indexed in the Earth Impact Database

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.

Analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 AU, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partial-to-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions) while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature, even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

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.

We consider how an impact generated seismic pulse affects the surface of an asteroid distant from the impact site. With laboratory experiments on dry polydisperse gravel mixtures, we track the trajectories of particles ejected from the surface by a single strong upward propagating pressure pulse. High speed video images show that ejecta trajectories are independent of particle size, and collisions primarily take place upon landing. When they land particles are ballistically sorted, as proposed by Shinbrot et al. (2017), leaving larger particles on the surface and smaller particles more widely dispersed. A single strong pulse can leave previously buried boulders stranded on the surface. Boulder stranding due to an impact excited seismic pulse is an additional mechanism that could leave large boulders present on the surface of rubble asteroids such as 162173 Ryugu, 101955 Bennu and 25143 Itokawa.

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.

A statistical analysis of 15,210 electron velocity distribution function (VDF) fits, observed within $\pm$2 hours of 52 interplanetary (IP) shocks by the $Wind$ spacecraft near 1 AU, is presented. This is the second in a three-part series on electron VDFs near IP shocks. The electron velocity moment statistics for the dense, low energy core, tenuous, hot halo, and field-aligned beam/strahl are a statistically significant list of values illustrated with both histograms and tabular lists for reference and baselines in future work. The beam/strahl fit results in the upstream are currently the closest thing to a proper parameterization of the beam/strahl electron velocity moments in the ambient solar wind. This work will also serve as a 1 AU baseline and reference for missions like $Parker \ Solar \ Probe$ and $Solar \ Orbiter$. The median density, temperature, beta, and temperature anisotropy values for the core(halo)[beam/strahl] components, with subscripts $ec$($eh$)[$eb$], of all fit results respectively are $n{\scriptstyle_{ec(h)[b]}}$ $\sim$ 11.3(0.36)[0.17] $cm^{-3}$, $T{\scriptstyle_{ec(h)[b], tot}}$ $\sim$ 14.6(48.4)[40.2] $eV$, $\beta{\scriptstyle_{ec(h)[b], tot}}$ $\sim$ 0.93(0.11)[0.05], and $\mathcal{A}{\scriptstyle_{ec(h)[b]}}$ $\sim$ 0.98(1.03)[0.93]. The nuanced details of the fitting method and data product description were published in Paper I and the detailed analysis of the results will be shown in Paper III.

Asteroids and other Small Solar System Bodies (SSSBs) are of high general and scientific interest in many aspects. The origin, formation, and evolution of our Solar System (and other planetary systems) can be better understood by analysing the constitution and physical properties of small bodies in the Solar System. Currently, two space missions (Hayabusa2, OSIRIS-REx) have recently arrived at their respective targets and will bring a sample of the asteroids back to Earth. Other small body missions have also been selected by, or proposed to, space agencies. The threat posed to our planet by near-Earth objects (NEOs) is also considered at the international level, and this has prompted dedicated research on possible mitigation techniques. The DART mission, for example, will test the kinetic impact technique. Even ideas for industrial exploitation have risen during the last years. Lastly, the origin of water and life on Earth appears to be connected to asteroids. Hence, future space mission projects will undoubtedly target some asteroids or other SSSBs. In all these cases and research topics, specific knowledge of the structure and mechanical behaviour of the surface as well as the bulk of those celestial bodies is crucial. In contrast to large telluric planets and dwarf planets, a large proportion of such small bodies is believed to consist of gravitational aggregates ('rubble piles') with no -- or low -- internal cohesion, with varying macro-porosity and surface properties (from smooth regolith covered terrain, to very rough collection of boulders), and varying topography (craters, depressions, ridges) [...].

Analysis of 15314 electron velocity distribution functions (VDFs) within $\pm$2 hours of 52 interplanetary (IP) shocks observed by the \emph{Wind} spacecraft near 1 AU are introduced. The electron VDFs are fit to the sum of three model functions for the cold dense core, hot tenuous halo, and field-aligned beam/strahl component. The best results were found by modeling the core as either a bi-kappa or a symmetric (or asymmetric) bi-self-similar velocity distribution function, while both the halo and beam/strahl components were best fit to bi-kappa velocity distribution function. This is the first statistical study to show that the core electron distribution is better fit to a self-similar velocity distribution function than a bi-Maxwellian under all conditions. The self-similar distribution deviation from a Maxwellian is a measure of inelasticity in particle scattering from waves and/or turbulence. The range of values defined by the lower and upper quartiles for the kappa exponents are $\kappa{\scriptstyle_{ec}}$ $\sim$ 5.40--10.2 for the core, $\kappa{\scriptstyle_{eh}}$ $\sim$ 3.58--5.34 for the halo, and $\kappa{\scriptstyle_{eb}}$ $\sim$ 3.40--5.16 for the beam/strahl. The lower-to-upper quartile range of symmetric bi-self-similar core exponents are $s{\scriptstyle_{ec}}$ $\sim$ 2.00--2.04, and asymmetric bi-self-similar core exponents are $p{\scriptstyle_{ec}}$ $\sim$ 2.20--4.00 for the parallel exponent, and $q{\scriptstyle_{ec}}$ $\sim$ 2.00--2.46 for the perpendicular exponent. The nuanced details of the fit procedure and description of resulting data product are also presented. The statistics and detailed analysis of the results are presented in Paper II and Paper III of this three-part study.

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.

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.

The science and origins of asteroids is deemed high priority in the Planetary Science Decadal Survey. Major scientific goals for the study of planetesimals are to decipher geological processes in SSSBs not determinable from investigation via in-situ experimentation, and to understand how planetesimals contribute to the formation of planets. Ground based observations are not sufficient to examine SSSBs, as they are only able to measure what is on the surface of the body; however, in-situ analysis allows for further, close up investigation as to the surface characteristics and the inner composure of the body. To this end, the Asteroid Mobile Imager and Geologic Observer (AMIGO) an autonomous semi-inflatable robot will operate in a swarm to efficiently characterize the surface of an asteroid. The stowed package is 10x10x10 cm (equivalent to a 1U CubeSat) that deploys an inflatable sphere of ~1m in diameter. Three mobility modes are identified and designed: ballistic hopping, rotation during hops, and up-righting maneuvers. Ballistic hops provide the AMIGO robot the ability to explore a larger portion of the asteroid's surface to sample a larger area than a stationary lander. Rotation during the hop entails attitude control of the robot, utilizing propulsion and reaction wheel actuation. In the event of the robot tipping or not landing up-right, a combination of thrusters and reaction wheels will correct the robot's attitude. The AMIGO propulsion system utilizes sublimate-based micro-electromechanical systems (MEMS) technology as a means of lightweight, low-thrust ballistic hopping and coarse attitude control. Each deployed AMIGO will hop across the surface of the asteroid multiple times.

There are thousands of asteroids in near-Earth space and millions expected in the Main Belt. They are diverse in their physical properties and compositions. They are also time capsules of the early Solar System making them valuable for planetary science, and are strategic for resource mining, planetary defense/security and as interplanetary depots. But we lack direct knowledge of the geophysical behavior of an asteroid surface under milligravity conditions, and therefore landing on an asteroid and manipulating its surface material remains a daunting challenge. Towards this goal we are putting forth plans for a 12U CubeSat that will be in Low Earth Orbit and that will operate as a spinning centrifuge on-orbit. In this paper, we will present an overview of the systems engineering and instrumentation design on the spacecraft. Parts of this 12U CubeSat will contain a laboratory that will recreate asteroid surface conditions by containing crushed meteorite. The laboratory will spin at 1 to 2 RPM during the primary mission to simulate surface conditions of asteroids 2 km and smaller, followed by an extended mission where the spacecraft will spin at even higher RPM. The result is a bed of realistic regolith, the environment that landers and diggers and maybe astronauts will interact with. The CubeSat is configured with cameras, lasers, actuators and small mechanical instruments to both observe and manipulate the regolith at low simulated gravity conditions. A series of experiments will measure the general behavior, internal friction, adhesion, dilatancy, coefficients of restitution and other parameters that can feed into asteroid surface dynamics simulations. Effective gravity can be varied, and external mechanical forces can be applied.

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.

There are thousands of asteroids in near-Earth space and millions in the Main Belt. They are diverse in physical properties and composition and are time capsules of the early solar system. This makes them strategic locations for planetary science, resource mining, planetary defense/security and as interplanetary depots and communication relays. Landing on a small asteroid and manipulating its surface materials remains a major unsolved challenge fraught with high risk. The asteroid surface may contain everything from hard boulders to soft regolith loosely held by cohesion and very low-gravity. Upcoming missions Hayabusa II and OSIRIS-REx will perform touch and go operations to mitigate the risks of landing on an asteroid. This limits the contact time and requires fuel expenditure for hovering. An important unknown is the problem of getting stuck or making a hard impact with the surface. The Spacecraft Penetrator for Increasing Knowledge of NEOs (SPIKE) mission concept will utilize a small-satellite bus that is propelled using a xenon-fueled ion engine and will contain an extendable, low-mass, high-strength boom with a tip containing force-moment sensors. SPIKE will enable contact with the asteroid surface, where it will perform detailed regolith analysis and seismology as well as penetrometry, while keeping the main spacecraft bus at a safe distance. Using one or more long stilts frees the spacecraft from having to hover above the asteroid and thus substantially reduces or eliminates fuel use when doing science operations. This enables much longer missions that include a series of hops to multiple locations on the small-body surface.

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.

We advocate a low-cost strategy for long-duration research into the 'milligravity' environment of asteroids, comets and small moons, where surface gravity is a vector field typically less than 1/1000 the gravity of Earth. Unlike the microgravity environment of space, there is a directionality that gives rise, over time, to strangely familiar geologic textures and landforms. In addition to advancing planetary science, and furthering technologies for hazardous asteroid mitigation and in-situ resource utilization, simplified access to long-duration milligravity offers significant potential for advancing human spaceflight, biomedicine and manufacturing. We show that a commodity 3U ($10\times10\times34$ cm$^3$) cubesat containing a laboratory of loose materials can be spun to 1 rpm = $2\pi/60$ s$^{-1}$ on its long axis, creating a centrifugal force equivalent to the surface gravity of a kilometer-sized asteroid. We describe the first flight demonstration, where small meteorite fragments will pile up to create a patch of real regolith under realistic asteroid conditions, paving the way for subsequent missions where landing and mobility technology can be flight-proven in the operational environment, in Low-Earth Orbit (LEO). The 3U design can be adapted for use onboard the International Space Station (ISS) to allow for variable gravity experiments under ambient temperature and pressure for a broader range of experiments.

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.

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.

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.

We investigate properties of a solar wind-like plasma including a secondary alpha particle population exhibiting a parallel temperature anisotropy with respect to the background magnetic field, using linear and quasi-linear predictions and by means of one-dimensional hybrid simulations. We show that anisotropic alpha particles can drive a parallel fire hose instability analogous to that generated by protons, but that, remarkably, the instability can be triggered also when the parallel plasma beta of alpha particles is below unity. The wave activity generated by the alpha anisotropy affects the evolution of the more abundant protons, leading to their anisotropic heating. When both ion species have sufficient parallel anisotropies both of them can drive the instability, and we observe generation of two distinct peaks in the spectra of the fluctuations, with longer wavelengths associated to alphas and shorter ones to protons. If a non-zero relative drift is present, the unstable modes propagate preferentially in the direction of the drift associated with the unstable species. The generated waves scatter particles and reduce their temperature anisotropy to marginally stable state, and, moreover, they significantly reduce the relative drift between the two ion populations. The coexistence of modes excited by both species leads to saturation of the plasma in distinct regions of the beta/anisotropy parameter space for protons and alpha particles, in good agreement with in situ solar wind observations. Our results confirm that fire hose instabilities are likely at work in the solar wind and limit the anisotropy of different ion species in the plasma.

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.

We investigate properties of the plasma fluid motion in the large amplitude low frequency fluctuations of highly Alfv\'enic fast solar wind. We show that protons locally conserve total kinetic energy when observed from an effective frame of reference comoving with the fluctuations. For typical properties of the fast wind, this frame can be reasonably identified by alpha particles, which, owing to their drift with respect to protons at about the Alfv\'en speed along the magnetic field, do not partake in the fluid low frequency fluctuations. Using their velocity to transform proton velocity into the frame of Alfv\'enic turbulence, we demonstrate that the resulting plasma motion is characterized by a constant absolute value of the velocity, zero electric fields, and aligned velocity and magnetic field vectors as expected for unidirectional Alfv\'enic fluctuations in equilibrium. We propose that this constraint, via the correlation between velocity and magnetic field in Alfv\'enic turbulence, is at the origin of the observed constancy of the magnetic field: while the constant velocity corresponding to constant energy can be only observed in the frame of the fluctuations, the correspondingly constant total magnetic field, invariant for Galilean transformations, remains the observational signature, in the spacecraft frame, of the constant total energy in the Alfv\'en turbulence frame.

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.

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.

Classical novae are expected to contribute to the 1809-keV Galactic $\gamma$-ray emission by producing its precursor $^{26}$Al, but the yield depends on the thermonuclear rate of the unmeasured $^{25}$Al($p,\gamma$)$^{26}$Si reaction. Using the $\beta$ decay of $^{26}$P to populate the key $J^{\pi}=3^+$ resonance in this reaction, we report the first evidence for the observation of its exit channel via a $1741.6 \pm 0.6 (\textrm{stat}) \pm 0.3 (\textrm{syst})$ keV primary $\gamma$ ray, where the uncertainties are statistical and systematic, respectively. By combining the measured $\gamma$-ray energy and intensity with other experimental data on $^{26}$Si, we find the center-of-mass energy and strength of the resonance to be $E_r = 414.9 \pm 0.6(\textrm{stat}) \pm 0.3 (\textrm{syst}) \pm 0.6(\textrm{lit.})$ keV and $\omega\gamma = 23 \pm 6 (\textrm{stat})^{+11}_{-10}(\textrm{lit.})$ meV, respectively, where the last uncertainties are from adopted literature data. We use hydrodynamic nova simulations to model $^{26}$Al production showing that these measurements effectively eliminate the dominant experimental nuclear-physics uncertainty and we estimate that novae may contribute up to 30% of the Galactic $^{26}$Al.

Collisionless shocks are loosely defined as shocks where the transition between pre-and post-shock states happens on a length scale much shorter than the collisional mean free path. In the absence of collision to enforce thermal equilibrium post-shock, electrons and ions need not have the same temperatures. While the acceleration of electrons for injection into shock acceleration processes to produce cosmic rays has received considerable attention, the related problem of the shock heating of quasi-thermal electrons has been relatively neglected. In this paper we review that state of our knowledge of electron heating in astrophysical shocks, mainly associated with supernova remnants (SNRs), shocks in the solar wind associated with the terrestrial and Saturnian bowshocks, and galaxy cluster shocks. The solar wind and SNR samples indicate that the ratio of electron temperature to ion temperature declines with increasing shock speed or Alfvenic Mach number. We discuss the extent to which such behavior can be understood via cosmic ray-generated waves in a shock precursor, which then subsequently damp by heating electrons. Finally, we speculate that a similar mechanism may be at work for both solar wind and SNR shocks.

Under sufficiently high electric field gradients, electron behaviour within exactly perpendicular shocks is unstable to the so-called trajectory instability. We extend previous work paying special attention to shortiscale, high amplitude structures as observed within the electric field profile. Via test particle simulations, we show that such structures can cause the electron distribution to heat in a manner that violates conservation of the first adiabatic invariant. This is the case even if the overall shock width is larger than the upstream electron gyroradius. The spatial distance over which these structures occur therefore constitutes a new scale length relevant to the shock heating problem. Furthermore, we find that the spatial location of the short-scale structure is important in determining the total effect of non-adiabatic behaviour - a result that has not been previously noted.

Electrons can be accelerated to ultrarelativistic energies at strong (high-Mach number) collisionless shock waves that form when stellar debris rapidly expands after a supernova. Collisionless shock waves also form in the flow of particles from the Sun (the solar wind), and extensive spacecraft observations have established that electron acceleration at these shocks is effectively absent whenever the upstream magnetic field is roughly parallel to the shock surface normal (quasi-parallel conditions). However, it is unclear whether this magnetic dependence of electron acceleration also applies to the far stronger shocks around young supernova remnants, where local magnetic conditions are poorly understood. Here we present Cassini spacecraft observations of an unusually strong solar system shock wave (Saturn's bow shock) where significant local electron acceleration has been confirmed under quasi-parallel magnetic conditions for the first time, contradicting the established magnetic dependence of electron acceleration at solar system shocks. Furthermore, the acceleration led to electrons at relativistic energies (~MeV), comparable to the highest energies ever attributed to shock-acceleration in the solar wind. These observations demonstrate that at high-Mach numbers, like those of young supernova remnant shocks, quasi-parallel shocks become considerably more effective electron accelerators.

We report here on serendipitous observations of the intense gamma-ray flare from SGR 1806-20 that occured on 27 December 2004. Unique data from the Cluster and Double Star-2 satellites, designed to study the Earth's magnetosphere, provide the first observational evidence of three separate timescales within the early (first 100ms) phases of this class of events. These observations reveal that, in addition to the initial very steep (<0.25ms) X-ray onset, there is firstly a 4.9ms exponential rise timescale followed by a continued exponential rise in intensity on a timescale of 70ms. These three timescales are a prominent feature of current theoretical models including the timescale (several ms) for fracture propagation in the crust of the neutron star.