Size Evolution of Close-in Super-Earths through Giant Impacts and Photoevaporation

The Kepler transit survey with follow-up spectroscopic observations has discovered numerous super-Earth sized planets and revealed intriguing features of their sizes, orbital periods, and their relations between adjacent planets. For the first time, we investigate the size evolution of planets via both giant impacts and photoevaporation to compare with these observed features. We calculate the size of a protoplanet, which is the sum of its core and envelope sizes, by analytical models. N-body simulations are performed to evolve planet sizes during the giant impact phase with envelope stripping via impact shocks. We consider the initial radial profile of the core mass and the initial envelope mass fractions as parameters. Inner planets can lose their whole envelopes via giant impacts, while outer planets can keep their initial envelopes, because they do not experience giant impacts. Photoevaporation is simulated to evolve planet sizes afterward. Our results suggest that the period-radius distribution of the observed planets would be reproduced if we perform simulations in which the initial radial profile of the core mass follows a wide range of power-law distributions and the initial envelope mass fractions are ∼0.1. Moreover, our model shows that the adjacent planetary pairs have similar sizes and regular spacings, with slight differences from detailed observational results such as the radius gap.

On the fate of water in the formation of rocky planets

<p>Rocky planets go through at least one and likely multiple magma ocean stages, produced by the giant impacts of accretion. Planetary data and models show that giant impacts do not dehydrate either the mantle or the atmosphere of their target planets. The magma ocean liquid consists of melted target material and melted impactor, and so will be dominated by silicate melt, and also contain dissolved volatiles including water, carbon, and sulfur compounds.</p><p>As the magma ocean cools and solidifies, water and other volatiles will be incorporated into the nominally anhydrous mantle phases up to their saturation limits, and will otherwise be enriched in the remaining, evolving magma ocean liquids. The water content of the resulting cumulate mantle is therefore the sum of the traces in the mineral grains, and any water in trapped interstitial liquids. That trapped liquid fraction may in fact be by far the largest contributor to the cumulate water budget.</p><p>The water and other dissolved volatiles in the evolving liquids may quickly reach the saturation limit of magmas near the surface, where pressure is low, but degassing the magma ocean is likely more difficult than has been assumed in some of our models. To degas into the atmosphere, the gases must exsolve from the liquid and form bubbles, and those bubbles must be able to rise quickly enough to avoid being dragged down by convection and re-dissolved at higher pressures. If bubbles are buoyant enough (that is, large enough) to decouple from flow and rise, then they are also dynamically unstable and liable to be torn into smaller bubbles and re-entrained. This conundrum led to the hypothesis that volatiles do not significantly degas until a high level of supersaturation is reached, and the bubbles form a buoyant layer and rise in diapirs in a continuum dynamics sense. This late degassing would have the twin effects of increasing the water content of the cumulates, and of speeding up cooling and solidification of the planet.</p><p>Once the mantle is solidified, the timeclock until the start of plate tectonics begins. Modern plate tectonics is thought to rely on water to lower the viscosity of the asthenosphere, but plate tectonics is also thought to be the process by which water is brought into the mantle. Magma ocean solidification, however, offers two relevant processes. First, following solidification the cumulate mantle is gravitationally unstable and overturns to stability, carrying water-bearing minerals from the upper mantle through the transition zone and into the lower mantle. Upon converting to lower-mantle phases, these minerals will release their excess water, since lower mantle phases have lower saturation limits, thus fluxing the upper mantle with water. Second, the mantle will be near its solidus temperature still, and thus its viscosity will be naturally low. When fluxed with excess water, the upper mantle would be expected to form a low degree melt, which if voluminous enough with rise to help form the earliest crust, and if of very low degree, will further reduce the viscosity of the asthenosphere.</p>

Partial core vaporization during giant impacts inferred from the entropy and the critical point of iron

<h3>The prevailing theory of the origin of the Moon is the giant impact hypothesis, in which a Mars-sized impactor collides with the proto-Earth in the late stage of accretion, and the Moon is subsequently formed from the proto-lunar disk made of the ejected materials. As the laboratory-scale experiments are not able to simulate planetary-scale impacts, our understanding of the giant impact mostly comes from hydrodynamic simulations. However, the results of these simulations heavily depend upon the available equation of state to describe the thermodynamic response of the constituent materials of the proto-Earth and impactor to shock waves.</h3><h3> </h3><h3>Iron as a building block material of the terrestrial planets naturally received significant attention. But the major effort has been put to determine its phase diagram up to the Earth’s core conditions (126-360 GPa and 3000-7000 K) and beyond. The studies of iron at low densities are still scarce and the position of the critical point (CP) is uncertain. As the liquid-vapor dome ends at CP, the position of the latter determines the time evolution of the proto-lunar disk during its condensation.</h3><h3> </h3><h3>In order to assess whether the core of the planets undergoes significant vaporization during a giant impact, we employ <em>ab initio</em> molecular-dynamics simulations to explore iron over a wide density region encompassing the critical point (CP) and the Hugoniot lines of the shocked iron cores. We determine the critical point of iron in the temperature range of 9000-9350 K, and the density range of 1.85-2.40 g/cm<sup>3</sup>, corresponding to a pressure range of 4-7 kbars [1]. This implies that the iron core of the proto-Earth may become supercritical after giant impacts. We show that the iron core of Theia partially vaporized during the Giant Impact. Part of this vapour may have remained in the disk, to eventually participate in the Moon’s small core. Similarly, during the late veneer stage a large fraction of the planetesimals have their cores undergoing partial vaporization. This would help to mix the highly siderophile elements into magma ponds or oceans.</h3><h3> </h3><h3>References:</h3><h3>[1] Z. Li, R. Caracas, F. Soubiran, Partial core vaporization during Giant Impacts inferred from the entropy and the critical point of iron, Earth Planet. Sci. Letters, 2020, https://doi.org/10.1016/j.epsl.2020.116463</h3><h3> </h3>

Planetary embryo collisions and the wiggly nature of extreme debris discs

ABSTRACT In this paper, we present results from a multistage numerical campaign to begin to explain and determine why extreme debris disc detections are rare, what types of impacts will result in extreme debris discs and what we can learn about the parameters of the collision from the extreme debris discs. We begin by simulating many giant impacts using a smoothed particle hydrodynamical code with tabulated equations of state and track the escaping vapour from the collision. Using an N-body code, we simulate the spatial evolution of the vapour generated dust post-impact. We show that impacts release vapour anisotropically not isotropically as has been assumed previously and that the distribution of the resulting generated dust is dependent on the mass ratio and impact angle of the collision. In addition, we show that the anisotropic distribution of post-collision dust can cause the formation or lack of formation of the short-term variation in flux depending on the orientation of the collision with respect to the orbit around the central star. Finally, our results suggest that there is a narrow region of semimajor axis where a vapour generated disc would be observable for any significant amount of time implying that giant impacts where most of the escaping mass is in vapour would not be observed often but this does not mean that the collisions are not occurring.

Crustal porosity reveals the bombardment history of the Moon

Planetary bombardment histories provide critical information regarding the formation and evolution of the Solar System and of the planets within it. These records evidence transient instabilities in the Solar System’s orbital evolution, giant impacts such as the Moon-forming impact, and material redistribution. Such records provide insight into planetary evolution, including the deposition of energy, delivery of materials, and crustal processing, specifically the modification of porosity. Bombardment histories are traditionally constrained from the surface expression of impacts — these records, however, are degraded by various geologic processes. Here we show that the Moon’s porosity contains a more complete record of its bombardment history. We find that the terrestrial planets were subject to double the number of ≥20-km-diameter-crater-forming impacts than are recorded on the lunar highlands, fewer than previously thought to have occurred. We show that crustal porosity doesn’t slowly increase as planets evolve, but instead is generated early in a planet’s evolution when most basins formed and decreases as planets evolve. We show that porosity constrains the relative ages of basins formed early in a planet’s evolution, a timeframe for which little information exists. These findings demonstrate that the Solar System was less violent than previously thought. Fewer volatiles and other materials were delivered to the terrestrial planets, consistent with estimates of the delivery of siderophiles and water to the Moon. High crustal porosity early in the terrestrial planets’ evolution slowed their cooling and enhanced their habitability. Several lunar basins formed early than previously considered, casting doubt on the existence of a late heavy bombardment.

Steam Atmospheres and Magma Oceans on Planets

A magma ocean is a global layer of partially or fully molten rocks. Significant melting of terrestrial planets likely occurs due to heat release during planetary accretion, such as decay heat of short-lived radionuclides, impact energy released by continuous planetesimal accretion, and energetic impacts among planetary-sized bodies (giant impacts). Over a magma ocean, all water, which is released upon impact or degassed from the interior, exists as superheated vapor, forming a water-dominated, steam atmosphere. A magma ocean extending to the surface is expected to interact with the overlying steam atmosphere through material and heat exchange. Impact degassing of water starts when the size of a planetary body becomes larger than Earth’s moon or Mars. The degassed water could build up and form a steam atmosphere on protoplanets growing by planetesimal accretion. The atmosphere has a role in preventing accretion energy supplied by planetesimals from escaping, leading to the formation of a magma ocean. Once a magma ocean forms, part of the steam atmosphere would start to dissolve into the surface magma due to the high solubility of water into silicate melt. Theoretical studies indicated that as long as the magma ocean is present, a negative feedback loop can operate to regulate the amount of the steam atmosphere and to stabilize the surface temperature so that a radiative energy balance is achieved. Protoplanets can also accrete the surrounding H2-rich disk gas. Water could be produced by oxidation of H2 by ferrous iron in the magma. The atmosphere and water on protoplanets could be a mixture of outgassed and disk-gas components. Planets formed by giant impact would experience a global melting on a short timescale. A steam atmosphere could grow by later outgassing from the interior. Its thermal blanketing and greenhouse effects are of great importance in controlling the cooling rate of the magma ocean. Due to the presence of a runaway greenhouse threshold, the crystallization timescale and water budget of terrestrial planets can depend on the orbital distance from the host star. The terrestrial planets in our solar system essentially have no direct record of their earliest history, whereas observations of young terrestrial exoplanets may provide us some insight into what early terrestrial planets and their atmosphere are like. Evolution of protoplanets in the framework of pebble accretion remains unexplored.

Losing oceans: The effects of composition on the thermal component of impact-driven atmospheric loss

ABSTRACT The formation of the Solar system’s terrestrial planets concluded with a period of giant impacts. Previous works examining the volatile loss caused by the impact shock in the moon-forming impact find atmospheric losses of at most 20–30 per cent and essentially no loss of oceans. However, giant impacts also result in thermal heating, which can lead to significant atmospheric escape via a Parker-type wind. Here we show that H2O and other high-mean molecular weight outgassed species can be efficiently lost through this thermal wind if present in a hydrogen-dominated atmosphere, substantially altering the final volatile inventory of terrestrial planets. We demonstrate that a giant impact during terrestrial planet formation can remove several Earth oceans’ worth of H2O, and other heavier volatile species, together with a primordial hydrogen-dominated atmosphere. These results may offer an explanation for the observed depletion in Earth’s light noble gas budget and for its depleted xenon inventory, which suggest that Earth underwent significant atmospheric loss by the end of its accretion. Because planetary embryos are massive enough to accrete primordial hydrogen envelopes and because giant impacts are stochastic and occur concurrently with other early atmospheric evolutionary processes, our results suggest a wide diversity in terrestrial planet volatile budgets.