We want to understand the global observable consequences of different physical processes and the initial system properties on the demographics of the planetary population. We select distributions of initial conditions that are representative of known protoplanetary discs and we use the Bern model to synthesise five planetary populations around 1 Sun-mass stars. Each population has a different initial number of Moon-mass embryos per disc to test for convergence: 1, 10, 20, 50, and 100. The last one is our nominal planetary population. The properties of giant planets do not change much starting with 10 embryos per system. For terrestrial planets inside 1 au, only the 100-embryos population is able to attain the giant-impact stage. In that population, each system contains 8 planets larger than 1 $M_\oplus$ on average. The fraction of systems with giants planets at all orbital distances is 18%, but only 1.6% are at > 10 au. Systems with giants contain on average 1.6 such planets. Inside of 1 au, the planet type with the highest occurrence rate and multiplicity are super-Earths (2.4 and 3.7), followed by terrestrial planets (1.6 and 2.8). Their frequency peaks at a stellar [Fe/H] of 0.0 and -0.2 respectively. The frequency of more massive planet types (Neptunian, giants) increases monotonically with [Fe/H]. Low-mass planets in the habitable zone are found around 44% of the stars with a multiplicity of 1.3, with a mean [Fe/H] of -0.11. This is one of the most comprehensive simulations of planetary system formation and evolution to date. For observations, the syntheses provide a large data set to search for comparison synthetic systems that show how they formed. For theory, they provide the framework to observationally test the consequences of specific physical processes. This is a important ingredient towards the development of a standard model of planetary formation and evolution.

Aims. Comparing theoretical models of exoplanet formation with observations allows one to make key step forward towards an understanding of planetary systems. It however requires a model able to (i) predict all the necessary observable quantities (not only masses and orbits, but also radii, luminosities, magnitudes, or evaporation rates) and (ii) address the large range in relevant planetary masses (from Mars mass to super-Jupiters) and distances (from stellar-grazing to wide orbits). Methods. We have developed a combined global planetary formation and evolution model that solves as directly as possible the underlying differential equations for the structure and evolution of the gas disc, the dynamical state of the planetesimals, the internal structure of the planets yielding their gas accretion rate and internal structure, the accretion rate of planetesimals, disc-driven orbital migration, and the gravitational interaction of concurrently forming planets via a full N-body calculation. Importantly, the model now also follows the long-term evolution of the planets on Gigayear timescales after formation including the effects of cooling and contraction, atmospheric escape, bloating, and stellar tides. Results. To test the model, we compared it with classical scenarios of Solar System formation. For the terrestrial planets, we find that we obtain a giant impact phase provided enough embryos (~100) are initially emplaced in the disc. For the giant planets, we find that Jupiter-mass planets must accrete their core shortly before the dispersal of the gas disc to prevent strong inward migration that would bring them to the inner edge of the disc. Conclusions. he model can form planetary systems with a wide range of properties. We find that giant-planet bearing systems are more diverse than the ones with only terrestrial planets.

This paper examines the effectiveness of an asteroid exploration architecture comprised of multiple nanosatellite sized spacecraft deployed from a single mother ship into a heliocentric orbit in the main asteroid belt where the mothership is ideally located in region of high density. Basic mission requirements associated with a Mother-Daughter architecture are established utilizing a relatively large number (10-20) daughter spacecraft distributed from a mothership within the asteroid belt for the purpose of executing sample and return missions. A number of trade analyses are performed to establish system performance to changes in initial orbit, delta-V capability and maximum small spacecraft flight time. The balance between the initial delta-V burn and asteroid velocity matching are also examined, with a goal of minimizing the amount of fuel needed in the small spacecraft. Preliminary requirements for the system are established using these results, and a conceptual design is presented for comparison to other asteroid exploration techniques. Preliminary results indicate that the aforementioned concept of a mothership with small spacecraft is viable and should be considered as an alternative approach to first order surveying of the asteroid belt.

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.

Graze-and-merge collisions (GMCs) are common multi-step mergers occurring in low-velocity off-axis impacts between similar sized planetary bodies. The first impact happens at somewhat faster than the mutual escape velocity; for typical impact angles this does not result in immediate accretion, but the smaller body is slowed down so that it loops back around and collides again, ultimately accreting. The scenario changes in the presence of a third major body, i.e. planets accreting around a star, or satellites around a planet. We find that when the loop-back orbit remains inside roughly 1/3 of the Hill radius from the target, then the overall process is not strongly affected. As the loop-back orbit increases in radius, the return velocity and angle of the second collision become increasingly random, with no record of the first collision's orientation. When the loop-back orbit gets to about 3/4 of the Hill radius, the path of smaller body is disturbed up to the point that it will usually escape the target.

The Lunar Gateway is expected to be positioned on-orbit around the Moon or in a Halo orbit at the L2 Lagrange point. The proposed Lunar Gateway is a game-changer for enabling new, high-priority lunar science utilizing Cu-beSats and presents a refreshing new opportunity for utilization of these small spacecraft as explorers. In context, CubeSats are being stretched to their limits as interplanetary explorers. The main technological hurdles include high-bandwidth communications and reliable high delta-v propulsion. Advances in deep-space attitude determination and control has been made possible from the recent NASA JPL MarCO missions. Due to these limitations, CubeSats are primarily designed to be dropped-off from a larger mission. The limited mass and volume have required compromises of the onboard science instruments, longer wait times to send back science data to Earth, shorter mission durations or higher accepted risk. With the Lunar Gateway being planned to be closer to the Moon, it will provide significant savings for a propulsion system and provide a primary relay for communication apart from the DSN and enable tele-operated command/control. These three factors can simplify the mission enabling routine deployment of CubeSats into lunar orbit and enable surface missions. In this paper, we present preliminary designs of 2 CubeSat lunar landers that will explore the lunar pits, Mare Tranquilitatis and the remnant magnetic fields Reiner Gamma.

In similar-sized planetary collisions, a significant part of the impactor often misses the target and continues downrange. We follow the dynamical evolution of "runners" from giant impacts to determine their ultimate fate. Surprisingly, runners re-impact their target planets only about half of the time, for realistic collisional and dynamical scenarios. Otherwise they remain in orbit for tens of millions of years (the limit of our N-body calculations) and longer, or sometimes collide with a different planet than the first one. When the runner does return to collide again with the same arget planet, its impact velocity is mainly constrained by the outcome of the prior collision. Impact angle and orientation, however, are unconstrained by the prior collision.

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.

We present a new method for the determination of the surface properties of airless bodies from measurements of the emitted infrared flux. Our approach uses machine learning techniques to train, validate, and test a neural network representation of the thermophysical behavior of the atmosphereless body given shape model, illumination and observational geometry of the remote sensors. The networks are trained on a dataset of thermal simulations of the emitted infrared flux for different values of surface rock abundance, roughness, and values of the thermal inertia of the regolith and of the rock components. These surrogate models are then employed to retrieve the surface thermal properties by Markov Chain Monte Carlo Bayesian inversion of observed infrared fluxes. We apply the method to the inversion of simulated infrared fluxes of asteroid (101195) Bennu -- according to a geometry of observations similar to those planned for NASA's OSIRIS-REx mission -- and infrared observations of asteroid (25143) Itokawa. In both cases, the surface properties of the asteroid -- such as surface roughness, thermal inertia of the regolith and rock component, and relative rock abundance -- are retrieved; the contribution from the regolith and rock components are well separated. For the case of Itokawa, we retrieve a rock abundance of about 85% for pebbles larger than the diurnal skin depth, which is about 2 cm. The thermal inertia of the rock is found to be lower than the expected value for LL chondrites, indicating that the rocks on Itokawa could be fractured. The average thermal inertia of the surface is around 750 $J s^{-1/2} K^{-1} m^{-2}$ and the measurement of thermal inertia of the regolith corresponds to an average regolith particle diameter of about 10 mm, consistently with in situ measurements as well as results from previous studies.

There are more than 750,000 asteroids identified in the main belt. These asteroids are diverse in composition and size. Some of these asteroids can be traced back to the early solar system and can provide insight into the origins of the so-lar system, origins of Earth and origins of life. Apart from being important tar-gets for science exploration, asteroids are strategically placed due to their low-gravity well, making it low-cost to transport material onto and way from them. They hold valuable resources such as water, carbon, metals including iron, nickel and platinum to name a few. These resources maybe used in refueling depots for interplanetary spacecraft and construction material for future space colonies, communication relays and space telescopes. The costs of getting to the main asteroid belt, combined with large numbers of objects to be explored encourage the application of small spacecraft swarms. The size and capability of the result-ing nano-spacecraft can make detection from Earth difficult. This paper dis-cusses a method by which a spacecraft can establish ephemeris autonomously using line of sight measurements to nearby asteroids with Extended Kalman Filtering techniques, without knowing accurate ephemeris of either the asteroids or the spacecraft initially. A description of the filter implementation is followed by examples of resultant performance.

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.

Exploration of small bodies, namely comets and asteroids remain a challenging endeavor due to their low gravity. The risk is so high that missions such as Hayabusa II and OSIRIS-REx will be performing touch and go missions to obtain samples. The next logical step is to perform longer-term mobility on the surface of these asteroid. This can be accomplished by sending small landers of a 1 kg or less with miniature propulsion systems that can just offset the force of asteroid gravity. Such a propulsion system would ideally be used to hop on the surface of the asteroid. Hopping has been found to be most efficient form of mobility on low-gravity. Use of wheels for rolling presents substantial challenges as the wheel can't gain traction to roll. The Asteroid Mobile Imager and Geologic Observer (AMIGO) utilizes 1 kg landers that are stowed in a 1U CubeSat configuration and deployed, releasing an inflatable that is 1 m in diameter. The inflatable is attached to the top of the 1U lander, enabling high speed communications and a means of easily tracking lander from an overhead mothership. Milligravity propulsion is required for the AMIGO landers to perform ballistic hops on the asteroid surface. The propulsion system is used to navigate the lander across the surface of the asteroid under the extremely low gravity while taking care to not exceed escape velocity.Although the concept for AMIGO missions is to use multiple landers, the more surface area evaluated by each lander the better. Without a propulsion system, each AMIGO will have a limited range of observable area. The propulsion system also serves as a rough attitude control system (ACS), as it enables pointing and regulation over where the lander is positioned via an array of MEMS thrusters.

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.

Terrestrial planets grew in a series of similar-sized collisions that swept up most of the next-largest bodies. Theia was accreted by the Earth to form the Moon according to the theory. Planetesimals likewise may have finished their accretion in a sequence of 'junior giant impacts', scaled down in size and velocity. This chapter considers the complicated physics of pairwise accretion, as planetesimals grow to planetary scales, and considers how the inefficiency of that process influences the origin of planetesimals and the diversity of meteorites and primary asteroids.

The wide availability of Commercial Off-The-Shelf (COTS) electronics that can withstand Low Earth Orbit conditions has opened avenue for wide deployment of CubeSats and small-satellites. CubeSats thanks to their low developmental and launch costs offer new opportunities for rapidly demonstrating on-orbit surveillance capabilities. In our earlier work, we proposed development of SWIMSat (Space based Wide-angle Imaging of Meteors) a 3U CubeSat demonstrator that is designed to observe illuminated objects entering the Earth's atmosphere. The spacecraft would operate autonomously using a smart camera with vision algorithms to detect, track and report of objects. Several CubeSats can track an object in a coordinated fashion to pinpoint an object's trajectory. An extension of this smart camera capability is to track unilluminated objects utilizing capabilities we have been developing to track and navigate to Near Earth Objects (NEOs). This extension enables detecting and tracking objects that can't readily be detected by humans. The system maintains a dense star map of the night sky and performs round the clock observations. Standard optical flow algorithms are used to obtain trajectories of all moving objects in the camera field of view. Through a process of elimination, certain stars maybe occluded by a transiting unilluminated object which is then used to first detect and obtain a trajectory of the object. Using multiple cameras observing the event from different points of view, it may be possible then to triangulate the position of the object in space and obtain its orbital trajectory. In this work, the performance of our space object detection algorithm coupled with a spacecraft guidance, navigation, and control system is demonstrated.

The geophysics of extrasolar planets is a scientific topic often regarded as standing largely beyond the reach of near-term observations. This reality in no way diminishes the central role of geophysical phenomena in shaping planetary outcomes, from formation, to thermal and chemical evolution, to numerous issues of surface and near-surface habitability. We emphasize that for a balanced understanding of extrasolar planets, it is important to look beyond the natural biases of current observing tools, and actively seek unique pathways to understand exoplanet interiors as best as possible during the long interim prior to a time when internal components are more directly accessible. Such pathways include but are not limited to: (a) enhanced theoretical and numerical modeling, (b) laboratory research on critical material properties, (c) measurement of geophysical properties by indirect inference from imprints left on atmospheric and orbital properties, and (d) the purpose-driven use of Solar System object exploration expressly for its value in comparative planetology toward exoplanet-analogs. Breaking down barriers that envision local Solar System exploration, including the study of Earth's own deep interior, as separate from and in financial competition with extrasolar planet research, may greatly improve the rate of needed scientific progress for exoplanet geophysics. As the number of known rocky and icy exoplanets grows in the years ahead, we expect demand for expertise in 'exogeoscience' will expand at a commensurately intense pace. We highlight key topics, including: how water oceans below ice shells may dominate the total habitability of our galaxy by volume, how free-floating nomad planets may often attain habitable subsurface oceans supported by radionuclide decay, and how deep interiors may critically interact with atmospheric mass loss via dynamo-driven magnetic fields.

Anorthosites that comprise the bulk of the lunar crust are believed to have formed during solidification of a Lunar Magma Ocean (LMO) in which these rocks would have floated to the surface. This early flotation crust would have formed a thermal blanket over the remaining LMO, prolonging solidification. Geochronology of lunar anorthosites indicates a long timescale of LMO cooling, or re-melting and re-crystallization in one or more late events. To better interpret this geochronology, we model LMO solidification in a scenario where the Moon is being continuously bombarded by returning projectiles released from the Moon-forming giant impact. More than one lunar mass of material escaped the Earth-Moon system onto heliocentric orbits following the giant impact, much of it to come back on returning orbits for a period of 100 Myr. If large enough, these projectiles would have punctured holes in the nascent floatation crust of the Moon, exposing the LMO to space and causing more rapid cooling. We model these scenarios using a thermal evolution model of the Moon that allows for production (by cratering) and evolution (solidification and infill) of holes in the flotation crust that insulates the LMO. For effective hole production, solidification of the magma ocean can be significantly expedited, decreasing the cooling time by more than a factor of 5. If hole production is inefficient, but shock conversion of projectile kinetic energy to thermal energy is efficient, then LMO solidification can be somewhat prolonged, lengthening the cooling time by 50% or more.

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.

We provide a fast method for computing constraints on impactor pre-impact orbits, applying this to the late giant impacts in the Solar System. These constraints can be used to make quick, broad comparisons of different collision scenarios, identifying some immediately as low-probability events, and narrowing the parameter space in which to target follow-up studies with expensive N-body simulations. We benchmark our parameter space predictions, finding good agreement with existing N-body studies for the Moon. We suggest that high-velocity impact scenarios in the inner Solar System, including all currently proposed single impact scenarios for the formation of Mercury, should be disfavoured. This leaves a multiple hit-and-run scenario as the most probable currently proposed for the formation of Mercury.

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

We 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.

Near Earth Asteroids (NEAs) are thought to contain a wealth of resources, including water, iron, titanium, nickel, platinum and silicates. Future space missions that can exploit these resources by performing In-Situ Resource Utilization (ISRU) gain substantial benefit in terms of range, payload capacity and mission flexibility. Compared to the Moon or Mars, the milligravity on some asteroids demands a fraction of the energy for digging and accessing hydrated regolith just below the surface. However, asteroids and small-bodies, because of their low gravity present a major challenge in landing, surface excavation and resource capture. These challenges have resulted in adoption of a "touch and go techniques", like the upcoming Osiris-rex sample-return mission. Previous asteroid excavation efforts have focused on discrete capture events (an extension of sampling technology) or whole-asteroid capture and processing. This paper analyzes the control of a bucket-wheel design for asteroid or small-body excavation. Our study focuses on system design of two counter rotating bucket-wheels that are attached to a hovering spacecraft. Regolith is excavated and heated to 1000 C to extract water. The water in turn is electrolyzed to produce hydrogen and oxygen for rocket fuel. We analyze control techniques to maximize traction of the bucket-wheels on the asteroid surface and minimize lift-off the surface, together with methods to dig deeper into the asteroid surface. Our studies combine analytical models, with simulation and hardware testing. For initial evaluation of material-spacecraft dynamics and mechanics, we assume lunar-like regolith for bulk density, particle size and cohesion. Our early studies point towards a promising pathway towards refinement of this technology for demonstration aboard a future space mission.

Exploration of asteroids and small-bodies can provide valuable insight into the origins of the solar system, into the origins of Earth and the origins of the building blocks of life. However, the low-gravity and unknown surface conditions of asteroids presents a daunting challenge for surface exploration, manipulation and for resource processing. This has resulted in the loss of several landers or shortened missions. Fundamental studies are required to obtain better readings of the material surface properties and physical models of these small bodies. The Asteroid Origins Satellite 1 (AOSAT 1) is a CubeSat centrifuge laboratory that spins at up to 4 rpm to simulate the milligravity conditions of sub 1 km asteroids. Such a laboratory will help to de-risk development and testing of landing and resource processing technology for asteroids. Inside the laboratory are crushed meteorites, the remains of asteroids. The laboratory is equipped with cameras and actuators to perform a series of science experiments to better understand material properties and asteroid surface physics. These results will help to improve our physics models of asteroids. The CubeSat has been designed to be low-cost and contains 3-axis magnetorquers and a single reaction-wheel to induce spin. In our work, we first analyze how the attitude control system will de-tumble the spacecraft after deployment. Further analysis has been conducted to analyze the impact and stability of the attitude control system to shifting mass (crushed meteorites) inside the spacecraft as its spinning in its centrifuge mode. AOSAT 1 will be the first in a series of low-cost CubeSat centrifuges that will be launched setting the stage for a larger, permanent, on-orbit centrifuge laboratory for experiments in planetary science, life sciences and manufacturing.

Current spacecraft need to launch with all of their required fuel for travel. This limits the system performance, payload capacity, and mission flexibility. One compelling alternative is to perform In-Situ Resource Utilization (ISRU) by extracting fuel from small bodies in local space such as asteroids or small satellites. Compared to the Moon or Mars, the microgravity on an asteroid demands a fraction of the energy for digging and accessing hydrated regolith just below the surface. Previous asteroid excavation efforts have focused on discrete capture events (an extension of sampling technology) or whole-asteroid capture and processing. This paper proposes an optimized bucket wheel design for surface excavation of an asteroid or small-body. Asteroid regolith is excavated and water extracted for use as rocket propellant. Our initial study focuses on system design, bucket wheel mechanisms, and capture dynamics applied to ponded materials known to exist on asteroids like Itokawa and Eros and small satellites like Phobos and Deimos. For initial evaluation of material-spacecraft dynamics and mechanics, we assume lunar-like regolith for bulk density, particle size and cohesion. We shall present our estimates for the energy balance of excavation and processing versus fuel gained. Conventional electrolysis of water is used to produce hydrogen and oxygen. It is compared with steam for propulsion and both show significant delta-v. We show that a return trip from Deimos to Earth is possible for a 12 kg craft using ISRU processed fuel.

Many asteroids are likely rubble-piles that are a collection of smaller objects held together by gravity and possibly cohesion. These asteroids are seismically shaken by impacts, which leads to excitation of their constituent particles. As a result it has been suggested that their surfaces and sub-surface interiors may be governed by a size sorting mechanism known as the Brazil Nut Effect. We study the behavior of a model asteroid that is a spherical, self-gravitating aggregate with a binary size-distribution of particles under the action of applied seismic shaking. We find that above a seismic threshold, larger particles rise to the surface when friction is present, in agreement with previous studies that focussed on cylindrical and rectangular box configurations. Unlike previous works we also find that size sorting takes place even with zero friction, though the presence of friction does aid the sorting process above the seismic threshold. Additionally we find that while strong size sorting can take place near the surface, the innermost regions remain unsorted under even the most vigorous shaking.

In this work, we analyze Jovian Trojans in the space of suitable resonant elements and we identify clusters of possible collisional origin by two independent methods: the hierarchical clustering and a so-called "randombox". Compared to our previous work (Bro\v{z} and Rozehnal 2011), we study a twice larger sample. Apart from Eurybates, Ennomos and $1996\,\rm RJ$ families, we have found three more clusters --- namely families around asteroids (20961)~Arkesilaos, (624)~Hektor in the $L_4$ libration zone and (247341)~$2001\,\rm UV_{209}$ in $L_5$. The families fulfill our stringent criteria, i.e. a high statistical significance, an albedo homogeneity and a steeper size-frequency distribution than that of background. In order to understand their nature, we simulate their long term collisional evolution with the Boulder code (Morbidelli et al. 2009) and dynamical evolution using a modified SWIFT integrator (Levison and Duncan, 1994). Within the framework of our evolutionary model, we were able to constrain the the age of the Hektor family to be either 1 to 4 Gyr or, less likely, 0.1 to 2.5 Gyr, depending on initial impact geometry. Since (624) Hektor itself seems to be a bilobed--shape body with a satellite (Marchis et al. 2014), i.e. an exceptional object, we address its association with the D--type family and we demonstrate that the moon and family could be created during a single impact event. We simulated the cratering event using a Smoothed Particle Hydrodynamics (SPH, Benz and Asphaug, 1994). This is also the first case of a family associated with a D--type parent body.

Mini-Neptunes and volatile-poor super-Earths coexist on adjacent orbits in proximity to host stars such as Kepler-36 and Kepler-11. Several post-formation processes have been proposed for explaining the origin of the compositional diversity: the mass loss via stellar XUV irradiation, degassing of accreted material, and in-situ accumulation of the disk gas. Close-in planets are also likely to experience giant impacts during the advanced stage of planet formation. This study examines the possibility of transforming volatile-rich super-Earths / mini-Neptunes into volatile-depleted super-Earths through giant impacts. We present the results of three-dimensional giant impact simulations in the accretionary and disruptive regimes. Target planets are modeled with a three-layered structure composed of an iron core, silicate mantle and hydrogen/helium envelope. In the disruptive case, the giant impact can remove most of the H/He atmosphere immediately and homogenize the refractory material in the planetary interior. In the accretionary case, the planet can retain more than half of the gaseous envelope, while a compositional gradient suppresses efficient heat transfer as its interior undergoes double-diffusive convection. After the giant impact, a hot and inflated planet cools and contracts slowly. The extended atmosphere enhances the mass loss via both a Parker wind induced by thermal pressure and hydrodynamic escape driven by the stellar XUV irradiation. As a result, the entire gaseous envelope is expected to be lost due to the combination of those processes in both cases. We propose that Kepler-36b may have been significantly devolatilized by giant impacts, while a substantial fraction of Kepler-36c's atmosphere may remain intact. Furthermore, the stochastic nature of giant impacts may account for the large dispersion in the mass--radius relationship of close-in super-Earths and mini-Neptunes.

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

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

We advance the modeling of rubble-pile solid bodies by re-examining the tidal breakup of comet Shoemaker-Levy 9, an event that occurred during a 1.33 Jupiter radii encounter with that planet in July 1992. Tidal disruption of the comet nucleus led to a chain of sub-nuclei about 100-1000 m in diameter; these went on to collide with the planet two years later (Chodas & Yeomans 1996). They were intensively studied prior to and during the collisions, making SL9 the best natural benchmark for physical models of small body disruption. For the first time in the study of this event, we use numerical codes treating rubble-piles as collections of polyhedra (Korycansky & Asphaug 2009). This introduces forces of dilatation and friction, and inelastic response. As in our previous studies (Asphaug & Benz 1994,1996) we conclude that the progenitor must have been a rubble-pile, and we obtain approximately the same pre-breakup diameter (about 1.5 km) in our best fits to the data. We find that the inclusion of realistic fragment shapes leads to grain locking and dilatancy, so that even in the absence of friction or other dissipation we find that disruption is overall more difficult than in our spheres-based simulations. We constrain the comet's bulk density at about 300-400 kg/m^3, half that of our spheres-based predictions and consistent with recent estimates derived from spacecraft observations.

This document covers five broad topics from theory of small bodies: planetesimal formation, cosmochemistry, thermal evolution, collisions and dynamics. Each of these topics is described in a separate section, where we prioritize selected main issues over completeness. The text points out the principal unresolved problems in each area, and suggests ways how progress can be made. This includes support for new code development, observations and experimental work that can be used to constrain theory, new research directions, and studies of cross-over regimes where the system's behavior is determined by several competing processes. The suggested development areas are placed in the context of NASA space exploration.

We use a smooth particle hydrodynamics method (SPH) to simulate colliding rocky and icy bodies from cm-scale to hundreds of km in diameter, in an effort to define self-consistently the threshold for catastrophic disruption. Unlike previous efforts, this analysis incorporates the combined effects of material strength (using a brittle fragmentation model) and self-gravitation, thereby providing results in the ``strength regime'' and the ``gravity regime'', and in between. In each case, the structural properties of the largest remnant are examined.