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.

High-speed impact ejecta at velocities comparable to the impact velocity are expected to contribute to material transport between planetary bodies and deposition of ejecta far from the impact crater. We investigated the behavior of high-speed ejecta produced at angles of 45 and 90 degrees, using both experimental and numerical methods. The experimental system developed at the Planetary Exploration Research Center of Chiba Institute of Technology (Japan) allowed us to observe the initial growth of the ejecta. We succeeded in imaging high-speed ejecta at 0.2 $\mathrm{{\mu}}$s intervals for impacts of polycarbonate projectiles of 4.8 mm diameter onto a polycarbonate plate at an impact velocity of ~4 km s$^{-1}$. Smoothed particle hydrodynamics (SPH) simulations of various numerical resolutions were conducted for the same impact conditions as pertaining to the experiments. We compared the morphology and velocities of the ejecta for the experiments and simulations, and we confirmed a close match for high-resolution simulations (with $\geq10^6$ SPH particles representing the projectile). According to the ejecta velocity distributions obtained from our high-resolution simulations, the ejection velocities of the high-speed ejecta for oblique impacts are much greater than those for vertical impacts. The translational motion of penetrating projectiles parallel to the target surface in oblique impacts could cause long-term, sustained acceleration at the root of the ejecta.

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.

Shock-induced metamorphism in meteorites informs us about the collisional environment and history of our solar system. Recently the importance of material strength in impact heating was reported from head-on impact simulations. Here, we perform three-dimensional oblique impact simulations, and confirm the additional heating due to material strength for oblique impacts. Despite a large difference in the peak pressure at the impact point at a given impact velocity, we find that the heated mass for an oblique impact is nearly the same as that for a head-on impact. Thus, our results differ from the previous finding that the heated mass decreases as the impact becomes more oblique, and show that the additional shear heating is more effective for oblique impacts than for head-on impacts. This also indicates that material ejected during oblique impact tends to experience lower shock pressures but higher temperatures.

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.

When planets receive insolation above a certain critical value called the runaway threshold, liquid surface water vaporizes completely, which forms the inner edge of the habitable zone. Because land planets can emit a large amount of radiation from the dry tropics, they have a higher runaway threshold than aqua planets do. Here we systematically investigated the runaway threshold for various surface water distributions using a three-dimensional dynamic atmosphere model. The runaway threshold for the meridionally uniform surface water distribution increases from the typical value for the aqua-planet regime (~130% S0) to one for the land-planet regime (~155% S0) as the dry surface area increases, where S0 is the present Earth's insolation. Although this result is similar to the previous work considering zonally uniform surface water distributions, the runaway threshold for the land-planet regime is quite low compared to that of the previous work. This is because a part of the tropical atmosphere is always wet for the meridionally uniform case. We also considered the surface water distributions determined by the Earth's, Mars' and Venus' topographies. We found that their runaway thresholds are close to that for the meridionally uniform cases, and the amount of water at the boundary between an aqua- and land-planet regime is around 10% of the Earth's ocean. This clearly shows that the runaway threshold is not determined uniquely by the luminosity of the central star, but it has a wide range caused by the surface water distribution of the terrestrial water planet itself.

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

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.

The abundance of highly siderophile elements (HSEs) inferred for Mars' mantle from martian meteorites implies a Late Veneer (LV) mass addition of ~0.8 wt% with broadly chondritic composition. Late accretion to Mars by a differentiated Ceres-sized (~1000 km diameter) object can account for part of the requisite LV mass, and geochronological constraints suggests that this must have occurred no later than ca. 4480 Ma. Here, we analyze the outcome of the hypothetical LV giant impact to Mars with smoothed particle hydrodynamics simulations together with analytical theory. Results show that, in general about 50% of the impactor's metallic core shatters into ~10m fragments that subsequently fragment into sub-mm metallic hail at re-accretion. This returns a promising delivery of HSEs into martian mantle compared to either a head-on and hit-and-run collision; in both cases,<10% of impactor's core materials are fragmented and finally embedded in the martian mantle. Isotopic evidence from martian meteorites, and interpretations from atmospheric mapping data show that a global surface water reservoir could be present during the early Noachian (before ca. 4100 Ma). The millimeter-sized metal hail could thus react with a martian hydrosphere to generate ~3 bars of H2, which is adequate to act as a greenhouse and keep early Mars warm. Yet, we also find that this atmosphere is transient. It typically survives shorter than 3 Myr based on the expected extreme ultraviolet (EUV) flux of the early Sun; if the Sun was a slow rotator an accordingly weaker EUV flux could extend this lifetime to >10 Myr. A dense pre-Noachian CO2 atmosphere should lower the escape efficiency of hydrogen by IR emission. A more detailed hydrodynamic atmospheric model of this early hydrogen atmosphere is warranted to examine its effect on pre-Noachian Mars.

Hydrous minerals are found on the surfaces of asteroids, but their origin is not clear. If their origin is endogenic, the hydrous minerals that were formed in the inner part of a planetesimal (or parent body) should come out on to the surface without dehydration. If their origin is exogenic, the source of hydrous minerals accreting onto asteroids is needed. Collisions in the asteroid belt would be related to both origins because collisions excavate the surface and eject the materials. However, the fate of hydrous minerals in large planetesimals during the collisional process has not been well investigated. Here, we explore planetesimal collisions by using the iSALE-2D code, and investigate the effect of an impact for the target planetesimal containing hydrous minerals. Our numerical results for the fiducial case (5 km/s of the impact velocity) show that hydrous minerals are slightly heated during the collisions. This moderate heating indicates that they can avoid the dehydration reaction and keep their original composition. Some hydrous minerals have larger velocity than the escape velocity of the collision system. This means that hydrous minerals can escape from the planetesimal and support the theory of exogenic origin for the hydrous minerals on asteroids. Meanwhile, the velocity of other hydrous minerals is smaller than the escape velocity of the system. This also indicates the possibility of an endogenic origin for the hydrous minerals on asteroids. Our results suggest that hydrous minerals on asteroids can be provided by planetesimal collisions.

Since its inception six decades ago, astrobiology has diversified immensely to encompass several scientific questions including the origin and evolution of Terran life, the organic chemical composition of extraterrestrial objects, and the concept of habitability, among others. The detection of life beyond Earth forms the main goal of astrobiology, and a significant one for space exploration in general. This goal has galvanized and connected with other critical areas of investigation such as the analysis of meteorites and early Earth geological and biological systems, materials gathered by sample-return space missions, laboratory and computer simulations of extraterrestrial and early Earth environmental chemistry, astronomical remote sensing, and in-situ space exploration missions. Lately, scattered efforts are being undertaken towards the R&D of the novel and as-yet-space-unproven life-detection technologies capable of obtaining unambiguous evidence of extraterrestrial life, even if it is significantly different from Terran life. As the suite of space-proven payloads improves in breadth and sensitivity, this is an apt time to examine the progress and future of life-detection technologies.

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.

In most of the previous studies related to collisional disruption of planetesimals in the gravity regime, Smoothed Particle Hydrodynamics (SPH) simulations have been used. On the other hand, impact simulations using grid-based hydrodynamic code have not been sufficiently performed. In the present study, we execute impact simulations in the gravity regime using the shock-physics code iSALE, which is a grid-based Eulerian hydrocode. We examine the dependence of the critical specific impact energy Q_RD* on impact conditions for a wide range of specific impact energy (Q_R) from disruptive collisions to erosive collisions, and compare our results with previous studies. We find collision outcomes of the iSALE simulation agree well with those of the SPH simulation. Detailed analysis mainly gives three results. (1) The value of Q_RD* depends on numerical resolution, and is close to convergence with increasing numerical resolution. The difference in converged value of Q_RD* between the iSALE code and the SPH code is within 30%. (2) Ejected mass normalized by total mass (M_ej/M_tot) generally depends on various impact conditions. However, when Q_R is normalized by Q_RD* that is calculated for each impact simulation, M_ej/M_tot can be scaled by Q_R/Q_RD*, and is independent of numerical resolution, impact velocity and target size. (3) This similarity law for Q_R/Q_RD* is confirmed for a wide range of specific impact energy. We also derive a semi-analytic formula for Q_RD* based on the similarity law and the crater scaling law. We find that the semi-analytic formula for the case with a non-porous object is consistent with numerical results.

The atmospheres of the terrestrial planets are known to have been modified as a consequence of the impact degassing and atmospheric erosion during the late accretion. Despite the commonality of these processes, there are distinct gaps -- roughly two orders of magnitude -- between the abundances of noble gases and nitrogen in the present-day atmospheres on Venus, Earth, and Mars. The element partitioning on planetary surfaces is thought to be significantly different between the three planets ~4 Ga: the runaway greenhouse on Venus, the carbon-silicate cycle and ocean formation on Earth, and the CO2-ice and H2O-ice formation on Mars. Consequences of element partitioning for the atmospheric evolution during the late accretion onto Venus, Earth, and Mars are investigated with a numerical model. We set upper limits to the partial pressures of CO2 and H2O on Earth and Mars, which corresponds to the state of phase equilibrium and carbon-silicate cycle. The final N2 mass shrinks by ~40% and ~15% for Earth and Mars, respectively. The effect of element partitioning is found to be insufficient to reproduce the gaps. For Venus, the survival of the primordial atmosphere through the late accretion may partially account for the present-day atmosphere. Whereas on Mars, the atmospheric escape due to solar extreme UV and wind may have also influenced the atmospheric evolution.

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.

Liquid water is one of the most important materials affecting the climate and habitability of a terrestrial planet. Liquid water vaporizes entirely when planets receive insolation above a certain value, which is called the runaway greenhouse threshold. This threshold forms the inner most limit of the habitable zone. Here, we investigate the effect of the distribution of surface water on the runaway greenhouse threshold for Earth-sized planets using a three-dimensional dynamic atmosphere model. We considered a 1-bar atmosphere whose composition is similar to the current Earth's atmosphere with a zonally uniform distribution of surface water. As previous studies have already shown, we also recognized two climate regimes: the land planet regime, which has dry low latitude and wet high latitude region, and the aqua planet regime, which is globally wet. We showed that each regime is controlled by the width of the Hadley circulation, the amount of surface water, and the planetary topography. We found that the runaway greenhouse threshold varies continuously with the surface water distribution from about 130% (an aqua planet) to 180% (the extreme case of a land planet) of the present insolation at Earth's orbit. Our results indicate that the inner edge of the habitable zone is not a single sharp boundary, but a border whose location varies depending on planetary surface condition, such as the amount of surface water. Since land planets have wider habitable zones and less cloud cover, land planets would be good targets for future observations investigating planetary habitability.

Hypervelocity impacts cause significant heating of planetary bodies. Such events are recorded by a reset of 40Ar-36Ar ages and/or impact melts. Here, we investigate the influence of friction and plastic deformation in shock-generated comminuted rocks on the degree of impact heating using the iSALE shock-physics code. We demonstrate that conversion from kinetic to internal energy in the targets with strength occurs during pressure release, and additional heating becomes significant for low-velocity impacts (<10 km/s). This additional heat reduces the impact-velocity thresholds required to heat the targets with the 0.1 projectile mass to temperatures for the onset of Ar loss and melting from 8 and 10 km/s, respectively, for strengthless rocks to 2 and 6 km/s for typical rocks. Our results suggest that the impact conditions required to produce the unique features caused by impact heating span a much wider range than previously thought.

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.

Overabundances in highly siderophile elements (HSEs) of Earth's mantle can be explained by conveyance from a singular, immense (3000 km in a diameter) "Late Veneer" impactor of chondritic composition, subsequent to lunar formation and terrestrial core-closure. Such rocky objects of approximately lunar mass (about 0.01 M_E) ought to be differentiated, such that nearly all of their HSE payload is sequestered into iron cores. Here, we analyze the mechanical and chemical fate of the core of such a Late Veneer impactor, and trace how its HSEs are suspended - and thus pollute - the mantle. For the statistically most-likely oblique collision (about 45degree), the impactor's core elongates and thereafter disintegrates into a metallic hail of small particles (about 10 m). Some strike the orbiting Moon as sesquinary impactors, but most re-accrete to Earth as secondaries with further fragmentation. We show that a single oblique impactor provides an adequate amount of HSEs to the primordial terrestrial silicate reservoirs via oxidation of (<m-sized) metal particles with a hydrous, pre-impact, early Hadean Earth.

Hypervelocity ejection of material by impact spallation is considered a plausible mechanism for material exchange between two planetary bodies. We have modeled the spallation process during vertical impacts over a range of impact velocities from 6 to 21 km/s using both grid- and particle-based hydrocode models. The Tillotson equations of state, which are able to treat the nonlinear dependence of density on pressure and thermal pressure in the strongly shocked matter, were used to study the hydrodynamic and thermodynamic response after impacts. The effects of material strength and gravitational acceleration were not considered. A two-dimensional time-dependent pressure field within a 1.5-fold projectile radius from the impact point was investigated in cylindrical coordinates to address the generation of spalled material. A resolution test was also performed to reject ejected materials with peak pressures that were too low due to artificial viscosity. The relationship between ejection velocity veject and peak pressure Ppeak was also derived. Our approach shows that late stage acceleration in an ejecta curtain occurs due to the compressible nature of the ejecta, resulting in an ejection velocity that can be higher than the ideal maximum of the resultant particle velocity after passage of a shock wave. We also calculate the ejecta mass that can escape from a planet like Mars (i.e., veject higher than 5 km/s) that matches the petrographic constraints from Martian meteorites, and which occurs when Ppeak from 30-50 GPa. Although the mass of such ejecta is limited to from 0.1-1 percent of the projectile mass in vertical impacts, this is sufficient for spallation to have been a plausible mechanism for the ejection of Martian meteorites. Finally, we propose that impact spallation is a plausible mechanism for the generation of tektites.

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 Earth was born in violence. Many giant collisions of protoplanets are thought to have occurred during the terrestrial planet formation. Here we investigated the giant impact stage by using a hybrid code that consistently deals with the orbital evolution of protoplanets around the Sun and the details of processes during giant impacts between two protoplanets. A significant amount of materials (up to several tens of percent of the total mass of the protoplanets) is ejected by giant impacts. We call these ejected fragments the giant-impact fragments (GIFs). In some of the erosive hit-and-run and high-velocity collisions, metallic iron is also ejected, which comes from the colliding protoplanets' cores. From ten numerical simulations for the giant impact stage, we found that the mass fraction of metallic iron in GIFs ranges from ~ 1wt% to ~ 25wt%. We also discussed the effects of the GIFs on the dynamical and geochemical characteristics of formed terrestrial planets. We found that the GIFs have the potential to solve the following dynamical and geochemical conflicts: (1) The Earth, currently in a near circular orbit, is likely to have had a highly eccentric orbit during the giant impact stage. The GIFs are large enough in total mass to lower the eccentricity of the Earth to its current value via their dynamical friction. (2) The concentrations of highly siderophile elements (HSEs) in the Earth's mantle are greater than what was predicted experimentally. Re-accretion of the iron-bearing GIFs onto the Earth can contribute to the excess of HSEs. In addition, the estimated amount of iron-bearing GIFs provides significant reducing agent that could transform primitive CO2-H2O atmosphere and ocean into more reducing H2-bearing atmosphere. Thus, GIFs are important for the origin of Earth's life and its early evolution.

Disruptive collisions have been regarded as an important process for planet formation, while non-disruptive, small-scale collisions (hereafter called erosive collisions) have been underestimated or neglected by many studies. However, recent studies have suggested that erosive collisions are also important to the growth of planets, because they are much more frequent than disruptive collisions. Although the thresholds of the specific impact energy for disruptive collisions (Q_RD^*) have been investigated well, there is no reliable model for erosive collisions. In this study, we systematically carried out impact simulations of gravity-dominated planetesimals for a wide range of specific impact energy (Q_R) from disruptive collisions (Q_R ~ Q_RD^*) to erosive ones (Q_R << Q_RD^*) using the smoothed particle hydrodynamics method. We found that the ejected mass normalized by the total mass (M_ej/M_tot) depends on the numerical resolution, the target radius (R_tar) and the impact velocity (v_imp), as well as on Q_R, but that it can be nicely scaled by Q_RD^* for the parameter ranges investigated (R_tar = 30-300 km, v_imp = 2-5 km/s). This means that M_ej/M_tot depends only on Q_R/Q_RD^* in these parameter ranges. We confirmed that the collision outcomes for much less erosive collisions (Q_R < 0.01 Q_RD^*) converge to the results of an impact onto a planar target for various impact angles and that M_ej/M_tot = C * QR/QRD* holds. For disruptive collisions (Q_R ~ Q_RD^*), the curvature of the target has a significant effect on Mej/Mtot. We also examined the angle-averaged value of M_ej/M_tot and found that the numerically obtained relation between angle-averaged M_ej/M_tot and Q_R/Q_RD^* is very similar to the cases for 45-degree impacts. We proposed a new erosion model based on our numerical simulations for future research on planet formation with collisional erosion.

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.

At present, the giant impact (GI) is the most widely accepted model for the origin of the Moon. Most of the numerical simulations of GI have been carried out with the smoothed particle hydrodynamics (SPH) method. Recently, however, it has been pointed out that standard formulation of SPH (SSPH) has difficulties in the treatment of a contact discontinuity such as a core-mantle boundary and a free surface such as a planetary surface. This difficulty comes from the assumption of differentiability of density in SSPH. We have developed an alternative formulation of SPH, density independent SPH (DISPH), which is based on differentiability of pressure instead of density to solve the problem of a contact discontinuity. In this paper, we report the results of the GI simulations with DISPH and compare them with those obtained with SSPH. We found that the disk properties, such as mass and angular momentum produced by DISPH is different from that of SSPH. In general, the disks formed by DISPH are more compact: while formation of a smaller mass moon for low-oblique impacts is expected with DISPH, inhibition of ejection would promote formation of a larger mass moon for high-oblique impacts. Since only the improvement of core-mantle boundary significantly affects the properties of circumplanetary disks generated by GI and DISPH has not been significantly improved from SSPH for a free surface, we should be very careful when some conclusions are drawn from the numerical simulations for GI. And it is necessary to develop the numerical hydrodynamical scheme for GI that can properly treat the free surface as well as the contact discontinuity.

Two habitable planetary states are proposed: an aqua planet like the Earth and a land planet that has a small amount of water. Land planets keep liquid water under larger solar radiation compared to aqua planets. Water loss may change an aqua planet into a land planet, and the planet can remain habitable for a longer time than had it stayed an aqua planet. We calculate planetary evolution with hydrogen escape for different initial water inventories and different distances from the central star. We find that there are two conditions to evolve an aqua planet into a land planet: the critical amount of water on the surface M_ml consistent with a planet being a land planet, and the critical amount of water vapor in the atmosphere M_cv that defines the onset of the runaway greenhouse state. We find that Earth-size aqua planets with initial oceans < 10 % of the Earth's can evolve into land planets if M_cv = 3 m in precipitable water and M_ml = 5 % of the Earth's ocean mass. Such planets can keep liquid water on their surface for another 2 Gyrs. The initial amount of water and M_cv are shown to be important dividing parameters of the planetary evolution path. Our results indicate that massive hydrogen escape could give a fresh start as another kind of habitable planet rather than the end of its habitability.

Collisions are a fundamental process in planet formation. If colliding objects simply merge, a planetary object can grow. However, if the collision is disruptive, planetary growth is prevented. Therefore, the impact conditions under which collisions are destructive are important in understanding planet formation. So far, the critical specific impact energy for a disruptive collision Q_D^* has been investigated for various types of collisions between objects ranging in scale from centimeters to thousands of kilometers. Although the values of Q_D^* have been calculated numerically while taking into consideration various physical properties such as self-gravity, material strength, and porosity, the dependence of Q_D^* on numerical resolution has not been sufficiently investigated. In this paper, using the smoothed particle hydrodynamics (SPH) method, we performed numerical simulations of collisions between planetesimals at various numerical resolutions (from 2 x 10^4 to 5 x 10^6 SPH particles) and investigated the resulting variation in Q_D^*. The value of Q_D^* is shown to decrease as the number of SPH particles increases, and the difference between the Q_D^* values for the lowest and highest investigated resolutions is approximately a factor of two. Although the results for 5 x 10^6 SPH particles do not fully converge, higher-resolution simulations near the impact site show that the value of Q_D^* for the case with 5 x 10^6 SPH particles is close to the expected converged value. Although Q_D^* depends on impact parameters and material parameters, our results indicate that at least 5 x 10^6 SPH particles are required for numerical simulations in disruptive collisions to obtain the value of Q_D^* within 20% error.

In our solar system, Mars-sized protoplanets frequently collided with each other during the last stage of terrestrial planet formation called the giant impact stage. Giant impacts eject a large amount of material from the colliding protoplanets into the terrestrial planet region, which may form debris disks with observable infrared excesses. Indeed, tens of warm debris disks around young solar-type stars have been observed. Here, we quantitatively estimate the total mass of ejected materials during the giant impact stages. We found that $\sim$0.4 times the Earth's mass is ejected in total throughout the giant impact stage. Ejected materials are ground down by collisional cascade until micron-sized grains are blown out by radiation pressure. The depletion timescale of these ejected materials is determined primarily by the mass of the largest body among them. We conducted high-resolution simulations of giant impacts to accurately obtain the mass of the largest ejected body. We then calculated the evolution of the debris disks produced by a series of giant impacts and depleted by collisional cascades to obtain the infrared excess evolution of the debris disks. We found that the infrared excess is almost always higher than the stellar infrared flux throughout the giant impact stage ($\sim$100 Myr) and is sometimes $\sim$10 times higher immediately after a giant impact. Therefore, giant impact stages would explain the infrared excess from most observed warm debris disks. The observed fraction of stars with warm debris disks indicates that the formation probability of our solar system-like terrestrial planets is approximately 10%.

Although the two moons of Mars, Phobos and Deimos, have long been thought to be captured asteroids, recent observations of their compositions and orbits suggest that they may have formed from debris generated by one or more giant impacts of bodies with ~ 0.01 x target mass. Recent studies have both analytically estimated debris produced by giant impacts on Mars and numerically examined the evolution of circum-Mars debris disks. We perform a numerical study (Smoothed Particle Hydrodynamics simulation) of debris retention from giant impacts onto Mars, particularly in relation to a Borealis-scale giant impact (E ~ 3 x 10^29 J) capable of producing the Borealis basin. We find that a Borealis-scale impact is capable of producing a disk of mass ~ 5 x 10^20 kg (~ 1 - 4 % of the impactor mass), sufficient debris to form at least one of the martian moons according to recent numerical studies of martian debris disk evolution. While a Borealis-scale impact may generate sufficient debris to form both Phobos and Deimos, further studies of the debris disk evolution are necessary. Our results can serve as inputs for future studies of martian debris disk evolution.

At the final stage of terrestrial planet formation, known as the giant impact stage, a few tens of Mars-sized protoplanets collide with one another to form terrestrial planets. Almost all previous studies on the orbital and accretional evolution of protoplanets in this stage have been based on the assumption of perfect accretion, where two colliding protoplanets always merge. However, recent impact simulations have shown that collisions among protoplanets are not always merging events, that is, two colliding protoplanets sometimes move apart after the collision (hit-and-run collision). As a first step towards studying the effects of such imperfect accretion of protoplanets on terrestrial planet formation, we investigated the merging criteria for collisions of rocky protoplanets. Using the smoothed particle hydrodynamic (SPH) method, we performed more than 1000 simulations of giant impacts with various parameter sets, such as the mass ratio of protoplanets, $\gamma$, the total mass of two protoplanets, $M_{\rm T}$, the impact angle, $\theta$, and the impact velocity, $v_{\rm imp}$. We investigated the critical impact velocity, $v_{\rm cr}$, at the transition between merging and hit-and-run collisions. We found that the normalized critical impact velocity, $v_{\rm cr}/v_{\rm esc}$, depends on $\gamma$ and $\theta$, but does not depend on $M_{\rm T}$, where $v_{\rm esc}$ is the two-body escape velocity. We derived a simple formula for $v_{\rm cr}/v_{\rm esc}$ as a function of $\gamma$ and $\theta$, and applied it to the giant impact events obtained by \textit{N}-body calculations in the previous studies. We found that 40% of these events should not be merging events.

It is widely accepted that Titan and the mid-sized regular satellites around Saturn were formed in the circum-Saturn disk. Thus, if these mid-sized satellites were simply accreted by collisions of similar ice-rock satellitesimals in the disk, the observed wide diversity in density (i.e., the rock fraction) of the Saturnian mid-sized satellites is enigmatic. A recent circumplanetary disk model suggests satellite growth in an actively supplied circumplanetary disk, in which Titan-sized satellites migrate inward by interaction with the gas and are eventually lost to the gas planet. Here we report numerical simulations of giant impacts between Titan-sized migrating satellites and smaller satellites in the inner region of the Saturnian disk. Our results suggest that in a giant impact with impact velocity > 1.4 times the escape velocity and impact angle of ~45 degree, a smaller satellite is destroyed, forming multiple mid-sized satellites with a very wide diversity in satellite density (the rock fraction = 0-92 wt%). Our results of the relationship between the mass and rock fraction of the satellites resulting from giant impacts reproduce the observations of the Saturnian mid-sized satellites. Giant impacts also lead to internal melting of the formed mid-sized satellites, which would initiate strong tidal dissipation and geological activity, such as those observed on Enceladus today and Tethys in the past. Our findings also imply that giant impacts might have affected the fundamental physical property of the Saturnian mid-sized satellites as well as those of the terrestrial planets in the solar system and beyond.

The final stage of terrestrial planet formation is known as the giant impact stage where protoplanets collide with one another to form planets. So far this stage has been mainly investigated by N-body simulations with an assumption of perfect accretion in which all collisions lead to accretion. However, this assumption breaks for collisions with high velocity and/or a large impact parameter. We derive an accretion condition for protoplanet collisions in terms of impact velocity and angle and masses of colliding bodies, from the results of numerical collision experiments. For the first time, we adopt this realistic accretion condition in N-body simulations of terrestrial planet formation from protoplanets. We compare the results with those with perfect accretion and show how the accretion condition affects terrestrial planet formation. We find that in the realistic accretion model, about half of collisions do not lead to accretion. However, the final number, mass, orbital elements, and even growth timescale of planets are barely affected by the accretion condition. For the standard protoplanetary disk model, typically two Earth-sized planets form in the terrestrial planet region over about 100 M years in both realistic and perfect accretion models. We also find that for the realistic accretion model, the spin angular velocity is about 30% smaller than that for the perfect accretion model that is as large as the critical spin angular velocity for rotational instability. The spin angular velocity and obliquity obey Gaussian and isotropic distributions, respectively, independently of the accretion condition.

The origin of the Earth's ocean has been discussed on the basis of deuterium/hydrogen ratios (D/H) of several sources of water in the solar system. The average D/H of carbonaceous chondrites (CC's) is known to be close to the current D/H of the Earth's ocean, while those of comets and the solar nebula are larger by about a factor of two and smaller by about a factor of seven, respectively, than that of the Earth's ocean. Thus, the main source of the Earth's ocean has been thought to be CC's or adequate mixing of comets and the solar nebula. However, those conclusions are correct only if D/H of water on the Earth has remained unchanged for the past 4.5 Gyr. In this paper, we investigate evolution of D/H in the ocean in the case that the early Earth had a hydrogen-rich atmosphere, the existence of which is predicted by recent theories of planet formation no matter whether the nebula remains or not. Then we show that D/H in the ocean increases by a factor of 2-9, which is caused by the mass fractionation during atmospheric hydrogen loss, followed by deuterium exchange between hydrogen gas and water vapor during ocean formation. This result suggests that the apparent similarity in D/H of water between CC's and the current Earth's ocean does not necessarily support the CC's origin of water and that the apparent discrepancy in D/H is not a good reason for excluding the nebular origin of water.

The high density of the close-in extrasolar planet HD149026b suggests the presence of a huge core in the planet, which challenges planet formation theory. We first derive constraints on the amount of heavy elements and hydrogen/helium present in the planet: We find that preferred values of the core mass are between 50 and 80 M_E. We then investigate the possibility of subcritical core accretion as envisioned for Uranus and Neptune and find that the subcritical accretion scenario is unlikely in the case of HD149026b for at least two reasons: (i) Subcritical planets are such that the ratio of their core mass to their total mass is above ~0.7, in contradiction with constraints for all but the most extreme interior models of HD149026b; (ii) High accretion rates and large isolation mass required for the formation of a subcritical core of 30 M_E are possible only at specific orbital distances in a disk with a surface density of dust equal to at least 10 times that of the minimum mass solar nebula. This value climbs to 30 when considering a 50 M_E core. These facts point toward two main routes for the formation of this planet: (i) Gas accretion that is limited by a slow viscous inflow of gas in an evaporating disk; (ii) A significant modification of the composition of the planet after as accretion has stopped. These two routes are not mutually exclusive. Illustrating the second route, we show that for a wide range of impact parameters, giant impacts lead to a loss of the gas component of the planet and thus may lead to planets that are highly enriched in heavy elements. In the giant impact scenario, we expect an outer giant planet to be present. Observational studies by imaging, astrometry and long term interferometry of this system are needed to better narrow down the ensemble of possibilities.

From an astrobiological point of view, special attention has been paid to the probability of habitable planets in extrasolar systems. The purpose of this study is to constrain a possible range of the mass of a terrestrial planet that can get water. We focus on the process of water production through oxidation of the atmospheric hydrogen--the nebular gas having been attracted gravitationally--by oxide available at the planetary surface. For the water production to work well on a planet, a sufficient amount of hydrogen and enough high temperature to melt the planetary surface are needed. We have simulated the structure of the atmosphere that connects with the protoplanetary nebula for wide ranges of heat flux, opacity, and density of the nebular gas. We have found both requirements are fulfilled for an Earth-mass planet for wide ranges of the parameters. We have also found the surface temperature of planets of <= 0.3 Earth masses is lower than the melting temperature of silicate (~ 1500K). On the other hand, a planet of more than several Earth masses becomes a gas giant planet through runaway accretion of the nebular gas.