Numerous small bodies inevitably lead to cratering impacts on large planetary bodies during planet formation and evolution. As a consequence of these small impacts, a fraction of the target material escapes from the gravity of the large body, and a fraction of the impactor material accretes onto the target surface, depending on the impact velocities and angles. Here, we study the mass of the high-speed ejecta that escapes from the target gravity by cratering impacts when material strength is neglected. We perform a large number of cratering impact simulations onto a planar rocky target using the smoothed particle hydrodynamics method. We show that the escape mass of the target material obtained from our numerical simulations agrees with the prediction of a scaling law under a point-source assumption when $v_{\rm imp} \gtrsim 12 v_{\rm esc}$, where $v_{\rm imp}$ is the impact velocity and $v_{\rm esc}$ is the escape velocity of the target. However, we find that the point-source scaling law overestimates the escape mass up to a factor of $\sim 70$, depending on the impact angle, when $v_{\rm imp} \lesssim 12 v_{\rm esc}$. Using data obtained from numerical simulations, we derive a new scaling law for the escape mass of the target material for $v_{\rm imp} \lesssim 12 v_{\rm esc}$. We also derive a scaling law that predicts the accretion mass of the impactor material onto the target surface upon cratering impacts by numerically evaluating the escape mass of the impactor material. Our newly derived scaling laws are useful for predicting the escape mass of the target material and the accretion mass of the impactor material for a variety of cratering impacts that would occur on large planetary bodies during planet formation.

Throughout the history of the solar system, Mars has experienced continuous asteroidal impacts. These impacts have produced impact-generated Mars ejecta, and a fraction of this debris is delivered to Earth as Martian meteorites. Another fraction of the ejecta is delivered to the moons of Mars, Phobos and Deimos. Here, we studied the amount and condition of recent delivery of impact ejecta from Mars to its moons. Using state-of-the-art numerical approaches, we report, for the first time, that materials delivered from Mars to its moons are physically and chemically different from the Martian meteorites, which are all igneous rocks with a limited range of ages. We show that Mars ejecta mixed in the regolith of its moons potentially covers all its geological eras and consists of all types of rocks, from sedimentary to igneous. A Martian moons sample-return mission will bring such materials back to Earth, and the samples will provide a wealth of "time-resolved" geochemical information about the evolution of Martian surface environments.

Almost all the planets of our solar system have moons. Each planetary system has however unique characteristics. The Martian system has not one single big moon like the Earth, not tens of moons of various sizes like for the giant planets, but two small moons: Phobos and Deimos. How did form such a system? This question is still being investigated on the basis of the Earth-based and space-borne observations of the Martian moons and of the more modern theories proposed to account for the formation of other moon systems. The most recent scenario of formation of the Martian moons relies on a giant impact occurring at early Mars history and having also formed the so-called hemispheric crustal dichotomy. This scenario accounts for the current orbits of both moons unlike the scenario of capture of small size asteroids. It also predicts a composition of disk material as a mixture of Mars and impactor materials that is in agreement with remote sensing observations of both moon surfaces, which suggests a composition different from Mars. The composition of the Martian moons is however unclear, given the ambiguity on the interpretation of the remote sensing observations. The study of the formation of the Martian moon system has improved our understanding of moon formation of terrestrial planets: The giant collision scenario can have various outcomes and not only a big moon as for the Earth. This scenario finds a natural place in our current vision of the early solar system when conditions were favorable for giant collisions to occur. The next step in exploration of Martian moon is a sample return mission to test the giant collision scenario for their origin, and to provide tests of models of early solar system dynamics since Mars may retain material exchanged between the inner and outer solar system.

Potential microbial contamination of Martian moons, Phobos and Deimos, which can be brought about by transportation of Mars ejecta produced by meteoroid impacts on the Martian surface, has been comprehensively assessed in a statistical approach, based on the most probable history of recent major gigantic meteoroid collisions on the Martian surface. This article is the first part of our study to assess potential microbial density in Mars ejecta departing from the Martian atmosphere, as a source of the second part where statistical analysis of microbial contamination probability is conducted. Potential microbial density on the Martian surface as the source of microorganisms was estimated by analogy to the terrestrial areas having the similar arid and cold environments, from which a probabilistic function was deduced as the asymptotic limit. Microbial survival rate during hypervelocity meteoroid collisions was estimated by numerical analysis of impact phenomena with and without taking internal friction and plastic deformation of the colliding meteoroid and the target ground into consideration. Trajectory calculations of departing ejecta through the Martian atmosphere were conducted with taking account of aerodynamic deceleration and heating by the aid of computational fluid dynamic analysis. It is found that Mars ejecta smaller than 0.03 m in diameter hardly reach the Phobos orbit due to aerodynamic deceleration, or mostly sterilized due to significant aerodynamic heating even though they can reach the Phobos orbit and beyond. Finally, the baseline dataset of microbial density in Mars ejecta departing for Martian moons has been presented for the second part of our study.

This paper presents a case study of microbe transportation in the Mars-satellites system. We examined the spatial distribution of potential impact-transported microbes on the Martian moons using impact physics by following a companion study (Fujita et al.). We used sterilization data from the precede studies. We considered that the microbes came mainly from the Zunil crater on Mars. We found that 70-80% of the microbes are likely to be dispersed all over the moon surface and are rapidly sterilized due to radiation except for those microbes within a thick ejecta deposit produced by meteoroids. The other 20-30% might be shielded from radiation by thick regolith layers that formed at collapsed layers in craters produced by Mars rock impacts. The total number of potentially surviving microbes at the thick ejecta deposits is estimated to be 3-4 orders of magnitude lower than at the Mars rock craters. The microbe concentration is irregular in the horizontal direction and is largely depth-dependent due to the radiation sterilization. The surviving fraction of transported microbes would be only 1 ppm on Phobos and 100 ppm on Deimos, suggesting that the transport processes and radiation severely affect microbe survival. The microbe sampling probability from the Martian moons was also investigated. We suggest that sample return missions from the Martian moons are classified into Unrestricted Earth-Return missions for 30 g samples and 10 cm depth sampling, even in our conservative scenario. We also conducted a full statistical analysis for sampling the regolith of Phobos to include the effects of uncertainties in input parameters on the sampling probability. The most likely probability of microbial contamination for return samples is estimated to be two orders of magnitude lower than the criterion defined by the planetary protection policy of the Committee on Space Research (COSPAR).

It is important to clarify where and when rocky and icy planetesimals are formed in a viscously evolving disk. We wish to understand how local runaway pile-up of solids occurs inside or outside the snow line. We assume an icy pebble contains micron-sized silicate grains that are uniformly mixed with ice and are released during the ice sublimation. Using a local one-dimensional code, we solve the radial drift and the turbulent diffusion of solids and the water vapor, taking account of their sublimation/condensation around the snow line. We systematically investigate effects of back-reactions of the solids to gas on the radial drift and diffusion of solids, scale height evolution of the released silicate particles, and possible difference in effective viscous parameters between that for turbulent diffusion ($\alpha_{\rm tur}$) and that for the gas accretion rate onto the central star ($\alpha_{\rm acc}$). We study the dependence on the ratio of the solid mass flux to the gas ($F_{\rm p/g}$). We show that the favorable locations for the pile-up of silicate grains and icy pebbles are the regions in the proximity of the water snow line inside and outside it, respectively. We found that runaway pile-ups occur when both the back-reactions for radial drift and diffusion are included. In the case with only the back-reaction for the radial drift, no runaway pile-up is found except for extremely high pebble flux, while the condition of streaming instability can be satisfied for relatively large $F_{\rm p/g}$ as found in the past literatures. If the back-reactions for radial diffusion is considered, the runaway pile-up occurs for reasonable value of pebble flux. The runaway pile-up of silicate grains that would lead to formation of rocky planetesimals occurs for $\alpha_{\rm tur} \ll \alpha_{\rm acc}$, while the runaway pile-up of icy pebbles is favored for $\alpha_{\rm tur} \sim \alpha_{\rm acc}$.

Recent studies have revealed that all large (over 1000 km in diameter) trans-Neptunian objects (TNOs) form satellite systems. Although the largest Plutonian satellite, Charon, is thought to be an intact fragment of an impactor directly formed via a giant impact, whether giant impacts can explain the variations in secondary-to-primary mass ratios and spin/orbital periods among all large TNOs remains to be determined. Here we systematically perform hydrodynamic simulations to investigate satellite formation via giant impacts. We find that the simulated secondary-to-primary mass ratio varies over a wide range, which overlaps with observed mass ratios. We also reveal that the satellite systems' current distribution of spin/orbital periods and small eccentricity can be explained only when their spins and orbits tidally evolve: initially as fluid-like bodies, but finally as rigid bodies. These results suggest that all satellites of large TNOs were formed via giant impacts in the early stage of solar system formation, before the outward migration of Neptune, and that they were fully or partially molten during the giant impact era.

Recent works have shown that Martian moons Phobos and Deimos may have accreted within a giant impact-generated disk whose composition is about an equal mixture of Martian material and impactor material. Just after the giant impact, the Martian surface is heated up to $\sim3000-6000$ K and the building blocks of moons, including volatile-rich vapor, are heated up to $\sim2000$ K. In this paper, we investigate the volatile loss from the building blocks of Phobos and Deimos by hydrodynamic escape of vapor and radiation pressure on condensed particles. We show that a non-negligible amount of volatiles ($> 10\%$ of the vapor with temperature $> 1000$ K via hydrodynamic escape, and moderately volatile dusts that condense at $\sim700-2000$ K via radiation pressure) could be removed just after the impact during their first signle orbit from their pericenters to apocenters. Our results indicate that bulk Phobos and Deimos are depleted in volatile elements. Together with future explorations such as JAXA's MMX (Martian Moons eXploration) mission, our results would be used to constrain the origin of Phobos and Deimos.

Rings are ubiquitous around giant planets in our Solar System. They evolve jointly with the nearby satellite system. They could form either during the giant planet formation process or much later, as a result of large scale dynamical instabilities either in the local satellite system, or at the planetary scale. We review here the main characteristics of rings in our solar system, and discuss their main evolution processes and possible origin. We also discuss the recent discovery of rings around small bodies.

Observations and meteorites indicate that the Martian materials are enigmatically distributed within the inner solar system. A mega impact on Mars creating a Martian hemispheric dichotomy and the Martian moons can potentially eject Martian materials. A recent work has shown that the mega-impact-induced debris is potentially captured as the Martian Trojans and implanted in the asteroid belt. However, the amount, distribution, and composition of the debris has not been studied. Here, using hydrodynamic simulations, we report that a large amount of debris ($\sim 1\%$ of Mars' mass), including Martian crust/mantle and the impactor's materials ($\sim 20:80$), are ejected by a dichotomy-forming impact, and distributed between $\sim 0.5-3.0$ astronomical units. Our result indicates that unmelted Martian mantle debris ($\sim 0.02\%$ of Mars' mass) can be the source of Martian Trojans, olivine-rich asteroids in the Hungarian region and the main asteroid belt, and some of which even hit early Earth. A mega impact can naturally implant Martian mantle materials within the inner solar system.

The origin of Phobos and Deimos in a giant impact generated disk is gaining larger attention. Although this scenario has been the subject of many studies, an evaluation of the chemical composition of the Mars' moons in this framework is missing. The chemical composition of Phobos and Deimos is unconstrained. The large uncertainty about the origin of the mid-infrared features, the lack of absorption bands in the visible and near-infrared spectra, and the effects of secondary processes on the moons' surface make the determination of their composition very difficult from remote sensing data. Simulations suggest a formation of a disk made of gas and melt with their composition linked to the nature of the impactor and Mars. Using thermodynamic equilibrium we investigate the composition of dust (condensates from gas) and solids (from a cooling melt) that result from different types of Mars impactors (Mars-, CI-, CV-, EH-, comet-like). Our calculations show a wide range of possible chemical compositions and noticeable differences between dust and solids depending on the considered impactors. Assuming Phobos and Deimos as result of the accretion and mixing of dust and solids, we find that the derived assemblage (dust rich in metallic-iron, sulphides and/or carbon, and quenched solids rich in silicates) can be compatible with the observations. The JAXA's MMX (Martian Moons eXploration) mission will investigate the physical and chemical properties of the Maroons, especially sampling from Phobos, before returning to Earth. Our results could be then used to disentangle the origin and chemical composition of the pristine body that hit Mars and suggest guidelines for helping in the analysis of the returned samples.

Phobos and Deimos are the two small Martian moons, orbiting almost on the equatorial plane of Mars. Recent works have shown that they can accrete within an impact-generated inner dense and outer light disk, and that the same impact potentially forms the Borealis basin, a large northern hemisphere basin on the current Mars. However, there is no a priori reason for the impact to take place close to the north pole (Borealis present location) nor to generate a debris disk in the equatorial plane of Mars (in which Phobos and Deimos orbit). In this paper, we investigate these remaining issues on the giant impact origin of the Martian moons. First, we show that the mass deficit created by the Borealis impact basin induces a global reorientation of the planet to realign its main moment of inertia with the rotation pole (True Polar Wander). This moves the location of the Borealis basin toward its current location. Next, using analytical arguments, we investigate the detailed dynamical evolution of the eccentric inclined disk from the equatorial plane of Mars that is formed by the Martian-moon-forming impact. We find that, as a result of precession of disk particles due to the Martian dynamical flattening $J_{2}$ term of its gravity field and particle-particle inelastic collisions, eccentricity and inclination are damped and an inner dense and outer light equatorial circular disk is eventually formed. Our results strengthen the giant impact origin of Phobos and Deimos that can finally be tested by a future sample return mission such as JAXA's Martian Moons eXploration (MMX) mission.

Phobos and Deimos are the two small moons of Mars. Recent works have shown that they can accrete within an impact-generated disk. However, the detailed structure and initial thermodynamic properties of the disk are poorly understood. In this paper, we perform high-resolution SPH simulations of the Martian moon-forming giant impact that can also form the Borealis basin. This giant impact heats up the disk material (around $\sim 2000$ K in temperature) with an entropy increase of $\sim 1500$ J K$^{-1}$ kg$^{-1}$. Thus, the disk material should be mostly molten, though a tiny fraction of disk material ($< 5\%$) would even experience vaporization. Typically, a piece of molten disk material is estimated to be meter sized due to the fragmentation regulated by their shear velocity and surface tension during the impact process. The disk materials initially have highly eccentric orbits ($e \sim 0.6-0.9$) and successive collisions between meter-sized fragments at high impact velocity ($\sim 3-5$ km s$^{-1}$) can grind them down to $\sim100 \mu$m-sized particles. On the other hand, a tiny amount of vaporized disk material condenses into $\sim 0.1 \mu$m-sized grains. Thus, the building blocks of the Martian moons are expected to be a mixture of these different sized particles from meter-sized down to $\sim 100 \mu$m-sized particles and $\sim 0.1 \mu$m-sized grains. Our simulations also suggest that the building blocks of Phobos and Deimos contain both impactor and Martian materials (at least 35%), most of which come from the Martian mantle (50-150 km in depth; at least 50%). Our results will give useful information for planning a future sample return mission to Martian moons, such as JAXA's MMX (Martian Moons eXploration) mission.

The hypothesis of a recent origin of Saturn's rings and its mid-sized moons is actively debated. It was suggested that a proto-Rhea and a proto-Dione might have collided recently, giving birth to the modern system of mid-sized moons. It is also suggested that the rapid viscous spreading of the debris may have implanted mass inside Saturn's Roche limit, giving birth to the modern Saturn's ring system. However, this scenario has been only investigated in very simplified way for the moment. This paper investigates it in detail to assess its plausibility by using $N$-body simulations and analytical arguments. When the debris disk is dominated by its largest remnant, $N$-body simulations show that the system quickly re-accrete into a single satellite without significant spreading. On the other hand, if the disk is composed of small particles, analytical arguments suggest that the disk experiences dynamical evolutions in three steps. The disk starts significantly excited after the impact and collisional damping dominates over the viscous spreading. After the system flattens, the system can become gravitationally unstable when particles are smaller than $\sim$ 100 m. However, the particles grow faster than spreading. Then, the system becomes gravitationally stable again and accretion continues at a slower pace, but spreading is inhibited. Therefore, the debris is expected to re-accrete into several large bodies. In conclusion, our results show that such a scenario may not form the today's ring system. In contrast, our results suggest that today's mid-sized moons are likely re-accreted from such a catastrophic event.

The origin of rings around giant planets remains elusive. Saturn's rings are massive and made of 90-95% of water ice. In contrast, the much less massive rings of Uranus and Neptune are dark and likely to have higher rock fraction. Here we investigate, for the first time, the tidal disruption of a passing object, including the subsequent formation of planetary rings. First, we perform SPH simulations of the tidal destruction of big differentiated objects ($M_{\rm body}=10^{21-23}$) that experience close encounters with Saturn or Uranus. We find that about $0.1-10$% of the mass of the passing body is gravitationally captured around the planet. However, these fragments are initially big chunks and have highly eccentric orbits around the planet. Then, we perform N-body simulations including the planet's oblateness, starting with data obtained from the SPH simulations. Our N-body simulations show that the chunks are tidally destroyed during their next several orbits. Their individual orbits then start to precess incoherently around the planet's equator, which enhances their encounter velocities on longer-term evolution, resulting in more destructive impacts. These collisions would damp their eccentricities resulting in a progressive collapse of the debris cloud into a thin equatorial and low-eccentricity ring. These high energy impacts are expected to be catastrophic enough to produce small particles. Our numerical results also show that the mass of formed rings is large enough to explain current rings including inner regular satellites around Saturn and Uranus. In the case of Uranus, a body can go deeper inside the planet's Roche limit resulting in a more efficient capture of rocky material compared to Saturn's case in which mostly ice is captured. Thus, our results can naturally explain the compositional difference between the rings of Saturn, Uranus and Neptune.

Centaurs are minor planets orbiting between Jupiter and Neptune that have or had crossing orbits with one or more giant planets. Recent observations and reinterpretation of previous observations have revealed the existence of ring systems around 10199 Chariklo and 2060 Chiron. However, the origin of the ring systems around such a minor planet is still an open question. Here, we propose that the tidal disruption of a differentiated object that experiences a close encounter with a giant planet could naturally form diverse ring-satellite systems around the Centaurs. During the close encounter, the icy mantle of the passing object is preferentially ripped off by the planet's tidal force and the debris is distributed mostly within the Roche limit of the largest remnant body. Assuming the existence of $20-50$wt% silicate core below the icy mantle, a disk of particles is formed when the objects pass within $0.4-0.8$ of the planet's Roche limit with the relative velocity at infinity $3-6$km s$^{-1}$ and 8h initial spin period of the body. The resultant ring mass is $0.1-10$% of the central object's mass. Such particle disks are expected to spread radially, and materials spreading beyond the Roche limit would accrete into satellite(s). Our numerical results suggest that ring formation would be a natural outcome of such extreme close encounters and Centaurs can naturally have such ring systems because they cross the orbits of the giant planets.

Circumplanetary particle disks would be created in the late stage of planetary formation either by impacts of planetary bodies or disruption of satellites or passing bodies, and satellites can be formed by accretion of disk particles spreading across the Roche limit. Previous N-body simulation of lunar accretion focused on the formation of single-satellite systems from disks with large disk-to-planet mass ratios, while recent models of the formation of multiple-satellite systems from disks with smaller mass ratios do not take account of gravitational interaction between formed satellites. In the present work, we investigate satellite accretion from particle disks with various masses, using N-body simulation. In the case of accretion from somewhat less massive disks than the case of lunar accretion, formed satellites are not massive enough to clear out the disk, but can become massive enough to gravitationally shepherd the disk outer edge and start outward migration due to gravitational interaction with the disk. When the radial location of the 2:1 mean motion resonance of the satellite reaches outside the Roche limit, the second satellite can be formed near the disk outer edge, and then the two satellites continue outward migration while being locked in the resonance. Co-orbital satellites are found to be occasionally formed on the orbit of the first satellite. Our simulations also show that stochastic nature involved in gravitational interaction and collision between aggregates in the tidal environment can lead to diversity in the final mass and orbital architecture, which would be expected in satellite systems of exoplanets.