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.

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.