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

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

In contrast to the water-poor inner solar system planets, stochasticity during planetary formation and order of magnitude deviations in exoplanet volatile contents suggest that rocky worlds engulfed in thick volatile ice layers are the dominant family of terrestrial analogues among the extrasolar planet population. However, the distribution of compositionally Earth-like planets remains insufficiently constrained, and it is not clear whether the solar system is a statistical outlier or can be explained by more general planetary formation processes. Here we employ numerical models of planet formation, evolution, and interior structure, to show that a planet's bulk water fraction and radius are anti-correlated with initial $^{26}$Al levels in the planetesimal-based accretion framework. The heat generated by this short-lived radionuclide rapidly dehydrates planetesimals prior to accretion onto larger protoplanets and yields a system-wide correlation of planet bulk abundances, which, for instance, can explain the lack of a clear orbital trend in the water budgets of the TRAPPIST-1 planets. Qualitatively, our models suggest two main scenarios of planetary systems' formation: high-$^{26}$Al systems, like our solar system, form small, water-depleted planets, whereas those devoid of $^{26}$Al predominantly form ocean worlds, where the mean planet radii between both scenarios deviate by up to about 10%.

Rocky planetesimals in the early solar system melted internally and evolved chemically due to radiogenic heating from Al-26. Here we quantify the parametric controls on magma genesis and transport using a coupled petrological and fluid mechanical model of reactive two-phase flow. We find the mean grain size of silicate minerals to be a key control on magma ascent. For grain sizes larger than $\approx$ 1 mm, melt segregation produces distinct radial structure and chemical stratification. This stratification is most pronounced for bodies formed at around 1 Myr after formation of Ca,Al-rich inclusions. These findings suggest a link between the time and orbital location of planetesimal formation and their subsequent structural and chemical evolution. According to our models, the evolution of partially molten planetesimal interiors falls into two categories. In the magma ocean scenario, the whole interior of a planetesimal experiences nearly complete melting, which would result in turbulent convection and core-mantle differentiation by the rainfall mechanism. In the magma sill scenario, segregating melts gradually deplete the deep interior of the radiogenic heat source. In this case, magma may form melt-rich layers beneath a cool and stable lid, while core formation would proceed by percolation. Our findings suggest that grain sizes prevalent during the internal heating stage governed magma ascent in planetesimals. Regardless of whether evolution progresses toward a magma ocean or magma sill structure, our models predict that temperature inversions due to rapid Al-26 redistribution are limited to bodies formed earlier than $\approx$ 1 Myr after CAIs. We find that if grain size was smaller than $\approx$ 1 mm during peak internal melting, only elevated solid-melt density contrasts (such as found for enstatite chondrite compositions) would allow substantial melt segregation to occur.

The short-lived $^{182}$Hf-$^{182}$W decay system is a powerful chronometer for constraining the timing of metal-silicate separation and core formation in planetesimals and planets. Neutron capture effects on W isotopes, however, significantly hamper the application of this tool. In order to correct for neutron capture effects, Pt isotopes have emerged as a reliable in-situ neutron dosimeter. This study applies this method to IAB iron meteorites, in order to constrain the timing of metal segregation on the IAB parent body. The $\epsilon^{182}$W values obtained for the IAB iron meteorites range from -3.61 $\pm$ 0.10 to -2.73 $\pm$ 0.09. Correlating $\epsilon^{\mathrm{i}}$Pt with $^{182}$W data yields a pre-neutron capture $^{182}$W of -2.90 $\pm$ 0.06. This corresponds to a metal-silicate separation age of 6.0 $\pm$ 0.8 Ma after CAI for the IAB parent body, and is interpreted to represent a body-wide melting event. Later, between 10 and 14 Ma after CAI, an impact led to a catastrophic break-up and subsequent reassembly of the parent body. Thermal models of the interior evolution that are consistent with these estimates suggest that the IAB parent body underwent metal-silicate separation as a result of internal heating by short-lived radionuclides and accreted at around 1.4 $\pm$ 0.1 Ma after CAIs with a radius of greater than 60 km.

Water content and the internal evolution of terrestrial planets and icy bodies are closely linked. The distribution of water in planetary systems is controlled by the temperature structure in the protoplanetary disk and dynamics and migration of planetesimals and planetary embryos. This results in the formation of planetesimals and planetary embryos with a great variety of compositions, water contents and degrees of oxidation. The internal evolution and especially the formation time of planetesimals relative to the timescale of radiogenic heating by short-lived 26Al decay may govern the amount of hydrous silicates and leftover rock-ice mixtures available in the late stages of their evolution. In turn, water content may affect the early internal evolution of the planetesimals and in particular metal-silicate separation processes. Moreover, water content may contribute to an increase of oxygen fugacity and thus affect the concentrations of siderophile elements within the silicate reservoirs of Solar System objects. Finally, the water content strongly influences the differentiation rate of the icy moons, controls their internal evolution and governs the alteration processes occurring in their deep interiors.

Chondrules are the dominant bulk silicate constituent of chondritic meteorites and originate from highly energetic, local processes during the first million years after the birth of the Sun. So far, an astrophysically consistent chondrule formation scenario, explaining major chemical, isotopic and textural features, remains elusive. Here, we examine the prospect of forming chondrules from planetesimal collisions. We show that intensely melted bodies with interior magma oceans became rapidly chemically equilibrated and physically differentiated. Therefore, collisional interactions among such bodies would have resulted in chondrule-like but basaltic spherules, which are not observed in the meteoritic record. This inconsistency with the expected dynamical interactions hints at an incomplete understanding of the planetary growth regime during the protoplanetary disk phase. To resolve this conundrum, we examine how the observed chemical and isotopic features of chondrules constrain the dynamical environment of accreting chondrite parent bodies by interpreting the meteoritic record as an impact-generated proxy of planetesimals that underwent repeated collision and reaccretion cycles. Using a coupled evolution-collision model we demonstrate that the vast majority of collisional debris feeding the asteroid main belt must be derived from planetesimals which were partially molten at maximum. Therefore, the precursors of chondrite parent bodies either formed primarily small, from sub-canonical aluminum-26 reservoirs, or collisional destruction mechanisms were efficient enough to shatter planetesimals before they reached the magma ocean phase. Finally, we outline the window in parameter space for which chondrule formation from planetesimal collisions can be reconciled with the meteoritic record and how our results can be used to further constrain early solar system dynamics.

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

The thermal history and internal structure of chondritic planetesimals, assembled before the giant impact phase of chaotic growth, potentially yield important implications for the final composition and evolution of terrestrial planets. These parameters critically depend on the internal balance of heating versus cooling, which is mostly determined by the presence of short-lived radionuclides (SLRs), such as aluminum-26 and iron-60, as well as the heat conductivity of the material. The heating by SLRs depends on their initial abundances, the formation time of the planetesimal and its size. It has been argued that the cooling history is determined by the porosity of the granular material, which undergoes dramatic changes via compaction processes and tends to decrease with time. In this study we assess the influence of these parameters on the thermo-mechanical evolution of young planetesimals with both 2D and 3D simulations. Using the code family I2ELVIS/I3ELVIS we have run numerous 2D and 3D numerical finite-difference fluid dynamic models with varying planetesimal radius, formation time and initial porosity. Our results indicate that powdery materials lowered the threshold for melting and convection in planetesimals, depending on the amount of SLRs present. A subset of planetesimals retained a powdery surface layer which lowered the thermal conductivity and hindered cooling. The effect of initial porosity was small, however, compared to those of planetesimal size and formation time, which dominated the thermo-mechanical evolution and were the primary factors for the onset of melting and differentiation. We comment on the implications of this work concerning the structure and evolution of these planetesimals, as well as their behavior as possible building blocks of terrestrial planets.

Dauphas and Pourmand (2011) [Nature 473, 489--492] estimated the accretion timescale of Mars to be 1.8 $^{+0.9}_{-1.0}$ Myr from the W isotopes of martian meteorites. This timescale was derived assuming perfect metal-silicate equilibration between the impactor and the target's mantle. However, in the case of a small impactor most likely only a fraction of the target's mantle is involved in the equilibration, while only a small part of the impactor's core equilibrates in the case of a giant impact. We examined the effects of imperfect equilibration using results of high-resolution $N$-body simulations for the oligarchic growth stage. These effects were found to be small as long as a planetary embryo has a deep liquid magma ocean during its accretion. The effect due to partial involvement of the target's mantle in equilibration is small due to the low metal-silicate partition coefficient for W suggested from the low Hf/W ratio of the martian mantle. The effect due to partial involvement of the impactor's core is also small because a large fraction of the embryo mass is delivered from small planetesimals, which are likely to fully equilibrate in the deep magma ocean on the embryo. The accretion timescale of Mars estimated by the Hf-W chronology is shorter than that expected for the minimum mass solar nebula model as long as more than 10% of each impactor's core re-equilibrates with the martian mantle and the final stages of accretion are prolonged. This probably indicates that accretion of Mars proceeded rapidly due to solid and gas surface densities significantly larger than those for the minimum mass solar nebula or due to accretion of small fragments or pebbles.