We discuss major challenges in modeling giant impacts between planetary bodies, focusing on the equations of state (EOS). During the giant impact stage of planet formation, rocky planets are melted and partially vaporized. However, most EOS models fail to reproduce experimental constraints on the thermodynamic properties of the major minerals over the required phase space. Here, we present an updated version of the widely-used ANEOS model that includes a user-defined heat capacity limit in the thermal free energy term. Our revised model for forsterite (Mg$_2$SiO$_4$), a common proxy for the mantles of rocky planets, provides a better fit to material data over most of the phase space of giant impacts. We discuss the limitations of this model and the Tillotson equation of state, a commonly used alternative model.

Terrestrial planets are thought to experience episode(s) of large-scale melting early in their history. Fractionation during magma-ocean freezing leads to unstable stratification within the related cumulate layers due to progressive iron enrichment upward, but the effects of incremental cumulate overturns during MO crystallization remain to be explored. Here, we use geodynamic models with a moving-boundary approach to study convection and mixing within the growing cumulate layer, and thereafter within the fully crystallized mantle. For fractional crystallization, cumulates are efficiently stirred due to subsequent incremental overturns, except for strongly iron-enriched late-stage cumulates, which persist as a stably stratified layer at the base of the mantle for billions of years. Less extreme crystallization scenarios can lead to somewhat more subtle stratification. In any case, the long-term preservation of at least a thin layer of extremely enriched cumulates with Fe#>0.4, as predicted by all our models, is inconsistent with seismic constraints. Based on scaling relationships, however, we infer that final-stage Fe-rich magma-ocean cumulates originally formed near the surface should have overturned as small diapirs, and hence undergone melting and reaction with the host rock during sinking. The resulting moderately iron-enriched metasomatized/hybrid rock assemblages should have accumulated at the base of the mantle, potentially fed an intermittent basal magma ocean, and be preserved through the present-day. Such moderately iron-enriched rock assemblages can reconcile the physical properties of the large low shear-wave velocity provinces in the presentday lower mantle. Thus, we reveal Hadean melting and rock-reaction processes by integrating magmaocean crystallization models with the seismic-tomography snapshot.