Pressure is a key parameter in the physics and chemistry of planet formation and evolution. Previous studies have erroneously assumed that internal pressures monotonically increase with the mass of a body. Using smoothed particle hydrodynamics and potential field method calculations, we demonstrate that the hot, rapidly-rotating bodies produced by giant impacts can have much lower internal pressures than cool, slowly-rotating planets of the same mass. Pressures subsequently increase due to thermal and rotational evolution of the body. Using the Moon-forming impact as an example, we show that the internal pressures after the collision could have been less than half that in present-day Earth. The current pressure profile was not established until Earth cooled and the Moon receded, a process that may take up to 10s Myr after the last giant impact. Our work defines a new paradigm for pressure evolution during accretion of terrestrial planets: stochastic changes driven by impacts.

Scientific visualization tools are currently not optimized to create cinematic, production-quality representations of numerical data for the purpose of science communication. In our pipeline \texttt{Estra}, we outline a step-by-step process from a raw simulation into a finished render as a way to teach non-experts in the field of visualization how to achieve production-quality outputs on their own. We demonstrate feasibility of using the visual effects software Houdini for cinematic astrophysical data visualization, informed by machine learning clustering algorithms. To demonstrate the capabilities of this pipeline, we used a post-impact, thermally-equilibrated Moon-forming synestia from \cite{Lock18}. Our approach aims to identify "physically interpretable" clusters, where clusters identified in an appropriate phase space (e.g. here we use a temperature-entropy phase-space) correspond to physically meaningful structures within the simulation data. Clustering results can then be used to highlight these structures by informing the color-mapping process in a simplified Houdini software shading network, where dissimilar phase-space clusters are mapped to different color values for easier visual identification. Cluster information can also be used in 3D position space, via Houdini's Scene View, to aid in physical cluster finding, simulation prototyping, and data exploration. Our clustering-based renders are compared to those created by the Advanced Visualization Lab (AVL) team for the full dome show "Imagine the Moon" as proof of concept. With \texttt{Estra}, scientists have a tool to create their own production-quality, data-driven visualizations.

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.

Giant impacts dominate the final stages of terrestrial planet formation and set the configuration and compositions of the final system of planets. A giant impact is believed to be responsible for the formation of Earth's Moon, but the specific impact parameters are under debate. Because the canonical Moon-forming impact is the most intensely studied scenario, it is often considered the archetypal giant impact. However, a wide range of impacts with different outcomes are possible. Here we examine the total energy budgets of giant impacts that form Earth-mass bodies and find that they differ substantially across the wide range of possible Moon-forming events. We show that gravitational potential energy exchange is important, and we determine the regime in which potential energy has a significant effect on the collision outcome. Energy is deposited heterogeneously within the colliding planets, increasing their internal energies, and portions of each body attain sufficient entropy for vaporization. After gravitational re-equilibration, post-impact bodies are strongly thermally stratified, with varying amounts of vaporized and supercritical mantle. The canonical Moon-forming impact is a relatively low energy event and should not be considered the archetype of accretionary giant impacts that form Earth-mass planets. After a giant impact, bodies are significantly inflated in size compared to condensed planets of the same mass, and there are substantial differences in the magnitudes of their potential, kinetic and internal energy components. As a result, the conditions for metal-silicate equilibration and the subsequent evolution of the planet may vary widely between different impact scenarios.

Quantifying the energy budget of Earth in the first few million years following the Moon-forming giant impact is vital to understanding Earth's initial thermal state and the dynamics of lunar tidal evolution. After the impact, the body was substantially vaporized and rotating rapidly, very different from the planet we know today. The subsequent evolution of Earth's energy budget, as the body cooled and angular momentum was transferred during lunar tidal recession, has not been accurately calculated with all relevant energy components included. Here, we use giant impact simulations and planetary structure models to calculate the energy budget at stages in Earth's evolution. We show that the figure and internal structure of Earth changed substantially during its post-impact evolution and that changes in kinetic, potential, and internal energy were all significant. These changes have important implications for the dynamics of tidal recession and the thermal structure of early Earth.

The giant impact hypothesis remains the leading theory for lunar origin. However, current models struggle to explain the Moon's composition and isotopic similarity with Earth. Here we present a new lunar origin model. High-energy, high-angular momentum giant impacts can create a post-impact structure that exceeds the corotation limit (CoRoL), which defines the hottest thermal state and angular momentum possible for a corotating body. In a typical super-CoRoL body, traditional definitions of mantle, atmosphere and disk are not appropriate, and the body forms a new type of planetary structure, named a synestia. Using simulations of cooling synestias combined with dynamic, thermodynamic and geochemical calculations, we show that satellite formation from a synestia can produce the main features of our Moon. We find that cooling drives mixing of the structure, and condensation generates moonlets that orbit within the synestia, surrounded by tens of bars of bulk silicate Earth (BSE) vapor. The moonlets and growing moon are heated by the vapor until the first major element (Si) begins to vaporize and buffer the temperature. Moonlets equilibrate with BSE vapor at the temperature of silicate vaporization and the pressure of the structure, establishing the lunar isotopic composition and pattern of moderately volatile elements. Eventually, the cooling synestia recedes within the lunar orbit, terminating the main stage of lunar accretion. Our model shifts the paradigm for lunar origin from specifying a certain impact scenario to achieving a Moon-forming synestia. Giant impacts that produce potential Moon-forming synestias were common at the end of terrestrial planet formation.

In the giant impact hypothesis for lunar origin, the Moon accreted from an equatorial circum-terrestrial disk; however the current lunar orbital inclination of 5 degrees requires a subsequent dynamical process that is still debated. In addition, the giant impact theory has been challenged by the Moon's unexpectedly Earth-like isotopic composition. Here, we show that tidal dissipation due to lunar obliquity was an important effect during the Moon's tidal evolution, and the past lunar inclination must have been very large, defying theoretical explanations. We present a new tidal evolution model starting with the Moon in an equatorial orbit around an initially fast-spinning, high-obliquity Earth, which is a probable outcome of giant impacts. Using numerical modeling, we show that the solar perturbations on the Moon's orbit naturally induce a large lunar inclination and remove angular momentum from the Earth-Moon system. Our tidal evolution model supports recent high-angular momentum giant impact scenarios to explain the Moon's isotopic composition and provides a new pathway to reach Earth's climatically favorable low obliquity.

During accretion, terrestrial bodies attain a wide range of thermal and rotational states, which are accompanied by significant changes in physical structure (size, shape, pressure and temperature profile, etc.). However, variations in structure have been neglected in most studies of rocky planet formation and evolution. Here, we present a new code, HERCULES, that solves for the equilibrium structure of planets as a series of overlapping constant-density spheroids. Using HERCULES and a smoothed particle hydrodynamics code, we show that Earth-like bodies display a dramatic range of morphologies. For any rotating planetary body, there is a thermal limit beyond which the rotational velocity at the equator intersects the Keplerian orbital velocity. Beyond this corotation limit (CoRoL), a hot planetary body forms a structure, which we name a synestia, with a corotating inner region connected to a disk-like outer region. By analyzing calculations of giant impacts and models of planet formation, we show that typical rocky planets are substantially vaporized multiple times during accretion. For the expected angular momentum of growing planets, a large fraction of post-impact bodies will exceed the CoRoL and form synestias. The common occurrence of hot, rotating states during accretion has major implications for planet formation and the properties of the final planets. In particular, the structure of post-impact bodies influences the physical processes that control accretion, core formation and internal evolution. Synestias also lead to new mechanisms for satellite formation. Finally, the wide variety of possible structures for terrestrial bodies also expands the mass-radius range for rocky exoplanets.