Growing protoplanets experience a number of impacts during the accretion stage. A large impactor hits the surface of a protoplanet and produces a magma ocean, where the impactor's iron emulsifies and experiences metal-silicate equilibration with the mantle of the protoplanet while it descends towards the base of the magma ocean. This process repeatedly occurs and determines the chemical compositions of both mantle and core. The partitioning is controlled by parameters such as the equilibration pressure and temperature, which are often associated with or assumed to be proportional to the pressure and temperature at the base of the magma ocean. The pressure and temperature depend on both the depth and shape of a magma ocean because a spatially confined magma ocean, namely, a melt pool, can have a larger equilibrium pressure than a radially uniform (global) magma ocean even if their melt volumes are the same. Here, we develop scaling laws for (1) the total internal energy gain due to an impact, and (2) the heat distribution within the mantle based on more than 100 smoothed particle hydrodynamic (SPH) simulations. We use Legendre polynomials to describe these scaling laws and determine their coefficients by linear regression, minimizing the error between our model and SPH simulations. The input parameters are the impact angle, total mass, impact velocity, and impactor-to-total mass ratio. We find that the internal energy gain by a large impact is well characterized by the summation of the kinetic energy and accretional potential energy release as a function of the impact angle. We determine that the equilibrium pressure at the base of a melt pool can be higher (by $10-50 \%$) than those obtained from conventional radially-uniform global magma ocean models. These melt scaling laws are publicly available on GitHub (https://github.com/mikinakajima/MeltScalingLaw).

Impact cratering is likely a primary agent of regolith generation on airless bodies. Regolith production via impact cratering has long been a key topic of study since the Apollo era. The evolution of regolith due to impact cratering, however, is not well understood. A better formulation is needed to help quantify the formation mechanism and timescale of regolith evolution. Here, we propose an analytically derived stochastic model that describes the evolution of regolith generated by small, simple craters. We account for ejecta blanketing as well as regolith infilling of the transient crater cavity. Our results show that the regolith infilling plays a key role in producing regolith. Our model demonstrates that, because of the stochastic nature of impact cratering, the regolith thickness varies laterally, which is consistent with earlier work. We apply this analytical model to the regolith evolution at the Apollo 15 site. The regolith thickness is computed considering the observed crater size-frequency distribution of small, simple lunar craters (< 381 m in radius for ejecta blanketing and < 100 m in radius for the regolith infilling). Allowing for some amount of regolith coming from the outside of the area, our result is consistent with an empirical result from the Apollo 15 seismic experiment. Finally, we find that the timescale of regolith growth is longer than that of crater equilibrium, implying that even if crater equilibrium is observed on a cratered surface, it is likely the regolith thickness is still evolving due to additional impact craters.

By the time icy objects grow to about the mass of Europa and silicate bodies grow to approximately lunar mass, a large high-speed impact can generate an intact melt region that allows the dense material to gravitationally segregate, forming a large density anomaly. If this anomaly generates sufficient differential stress, it rapidly segregates to the object's center, triggering whole body differentiation. We used an impact melting model based on the Hugoniot equations, the linear shock-particle velocity relationship, and the empirical relationship between shock pressure and distance coupled with a Monte Carlo simulation of the late accretion process to determine conditions under which large impacts trigger differentiation in icy bodies. In a gas-free environment, impacts of projectiles in the satellite's accretion zone have a small probability of triggering differentiation in surviving proto-satellites as small as Triton. The probability increases to 100% by the time surviving proto-satellites grow to the mass of Europa. The impact of heliocentric particles captured by the central planet also effectively triggers differentiation in the icy satellites. If the largest projectiles in the heliocentric distribution are comparable to the mass of the growing satellites, they trigger differentiation in all surviving proto-satellites by the time they grow to the mass of Callisto. However, Callisto and Ganymede have different probabilities of impact-induced differentiation if they captured a small fraction (0.01 and 0.5%) of the particles they accreted from a distribution of small (largest particles ~10% of the satellite's mass) heliocentric objects. Impact induced differentiation explains the apparent difference between these two objects in a plausible way.

Because of the high energies involved, giant impacts that occur during planetary accretion cause large degrees of melting. The depth of melting in the target body after each collision determines the pressure and temperature conditions of metal-silicate equilibration and thus geochemical fractionation that results from core-mantle differentiation. The accretional collisions involved in forming the terrestrial planets of the inner Solar System have been calculated by previous studies using N-body accretion simulations. Here we use the output from such simulations to determine the volumes of melt produced and thus the pressure and temperature conditions of metal-silicate equilibration, after each impact, as Earth-like planets accrete. For these calculations a parametrised melting model is used that takes impact velocity, impact angle and the respective masses of the impacting bodies into account. The evolution of metal-silicate equilibration pressures (as defined by evolving magma ocean depths) during Earth's accretion depends strongly on the lifetime of impact-generated magma oceans compared to the time interval between large impacts. In addition, such results depend on starting parameters in the N-body simulations, such as the number and initial mass of embryos. Thus, there is the potential for combining the results, such as those presented here, with multistage core formation models to better constrain the accretional history of the Earth.

Comet 103P/Hartley 2 has diurnally controlled, CO2-driven activity on the tip of the small lobe of its bilobate nucleus. Such activity is unique among the comet nuclei visited by spacecraft, and suggests that CO2 ice is very near the surface, which is inconsistent with our expectations of an object that thermophysically evolved for ~45 million years prior to entering the Jupiter Family of comets. Here we explain this pattern of activity by showing that a very plausible recent episode of rapid rotation (rotation period of ~11 [10-13] hours) would have induced avalanches in Hartley 2's currently active regions that excavated down to CO2-rich ices and activated the small lobe of the nucleus. At Hartley 2's current rate of spindown about its principal axis, the nucleus would have been spinning fast enough to induce avalanches ~3-4 orbits prior to the DIXI flyby (~1984-1991). This coincides with Hartley 2's discovery in 1986, and implies that the initiation of CO2 activity facilitated the comet's discovery. During the avalanches, the sliding material would either be lofted off the surface by gas activity, or possibly gained enough momentum moving downhill (toward the tip of the small lobe) to slide off the nucleus. Much of this material would have failed to reach escape velocity, and would reimpact the nucleus, forming debris deposits. The similar size frequency distribution of the mounds observed on the surface of Hartley 2 and chunks of material in its inner coma suggest that the 20-40 meter mounds observed by the DIXI mission on the surface of Hartley 2 are potentially these fallback debris deposits. As the nucleus spun down from a rotation period of ~11 hours to 18.34 hours at the time of the DIXI flyby, the location of potential minima, where materials preferentially settle, migrated about the surface, allowing us to place relative ages on most of the imaged terrains of the nucleus.

Previously proposed mechanisms have difficulty explaining the disruption of Comet C/2012 S1 (ISON) as it approached the Sun. We describe a novel cometary disruption mechanism whereby comet nuclei fragment and disperse through dynamic sublimation pressure, which induces differential stresses within the interior of the nucleus. When these differential stresses exceed its material strength, the nucleus breaks into fragments. We model the sublimation process thermodynamically and propose that it is responsible for the disruption of Comet ISON. We estimate the bulk unconfined crushing strength of Comet ISON's nucleus and the resulting fragments to be 0.5 Pa and 1-9 Pa respectively, assuming typical Jupiter Family Comet (JFC) albedos. However, if Comet ISON has an albedo similar to Pluto, this strength estimate drops to 0.2 Pa for the intact nucleus and 0.6-4 Pa for its fragments. Regardless of assumed albedo, these are similar to previous strength estimates of JFCs. This suggests that, if Comet ISON is representative of dynamically new comets, then low bulk strength is a primordial property of some comet nuclei, and not due to thermal processing during migration into the Jupiter Family.

Spectral modeling of the large infrared excess in the Spitzer IRS spectra of HD 172555 suggests that there is more than 10^19 kg of sub-micron dust in the system. Using physical arguments and constraints from observations, we rule out the possibility of the infrared excess being created by a magma ocean planet or a circumplanetary disk or torus. We show that the infrared excess is consistent with a circumstellar debris disk or torus, located at approximately 6 AU, that was created by a planetary scale hypervelocity impact. We find that radiation pressure should remove submicron dust from the debris disk in less than one year. However, the system's mid-infrared photometric flux, dominated by submicron grains, has been stable within 4 percent over the last 27 years, from IRAS (1983) to WISE (2010). Our new spectral modeling work and calculations of the radiation pressure on fine dust in HD 172555 provide a self-consistent explanation for this apparent contradiction. We also explore the unconfirmed claim that 10^47 molecules of SiO vapor are needed to explain an emission feature at 8 um in the Spitzer IRS spectrum of HD 172555. We find that unless there are 10^48 atoms or 0.05 Earth masses of atomic Si and O vapor in the system, SiO vapor should be destroyed by photo-dissociation in less than 0.2 years. We argue that a second plausible explanation for the 8 um feature can be emission from solid SiO, which naturally occurs in submicron silicate "smokes" created by quickly condensing vaporized silicate.