MISCIBILITY CALCULATIONS FOR WATER AND HYDROGEN IN GIANT PLANETS

Mixture of hydrogen and methane under planetary interior conditions

We employ first-principles molecular dynamics simulations to provide equation-of-state data, pair distribution functions (PDFs), diffusion coefficients, and band gaps of a mixture of hydrogen and methane under planetary interior conditions as relevant for Uranus, Neptune, and similar icy exoplanets. We test the linear mixing approximation, which is fulfilled within a few percent for the chosen P-T conditions. Evaluation of the PDFs reveals that methane molecules dissociate into carbon clusters and free hydrogen atoms at temperatures greater than 3000 K. At high temperatures, the clusters are found to be short-lived. Furthermore, we calculate the electrical conductivity from which we derive the non-metal-to-metal transition region of the mixture. We also calculate the electrical conductivity along the P-T profile of Uranus [N. Nettelmann et al., Planet. Space Sci., 2013, 77, 143-151] and observe the transition of the mixture from a molecular to an atomic fluid as a function of the radius of the planet. The density and temperature ranges chosen in our study can be achieved using dynamic shock compression experiments and seek to aid such future experiments. Our work also provides a relevant data set for a better understanding of the interior, evolution, luminosity, and magnetic field of the ice giants in our solar system and beyond.

Miscibility of rock and ice in the interiors of water worlds

Super-Earths and sub-Neptunes are the most common planet types in our galaxy. A subset of these planets is predicted to be water worlds, bodies that are rich in water and poor in hydrogen gas. The interior structures of water worlds have been assumed to consist of water surrounding a rocky mantle and iron core. In small planets, water and rock form distinct layers with limited incorporation of water into silicate phases, but these materials may interact differently during the growth and evolution of water worlds due to greater interior pressures and temperatures. Here, we use density functional molecular dynamics (DFT-MD) simulations to study the miscibility and interactions of enstatite (MgSiO₃), a major end-member silicate phase, and water (H₂O) at extreme conditions in water world interiors. We explore pressures ranging from 30 to 120 GPa and temperatures from 500 to 8000 K. Our results demonstrate that enstatite and water are miscible in all proportions if the temperature exceeds the melting point of MgSiO₃. Furthermore, we performed smoothed particle hydrodynamics simulations to demonstrate that the conditions necessary for rock-water miscibility are reached during giant impacts between water-rich bodies of 0.7–4.7 Earth masses. Our simulations lead to water worlds that include a mixed layer of rock and water.

Gibbs-ensemble Monte Carlo simulation of H2-H2O mixtures

The miscibility gap in hydrogen-water mixtures is investigated by conducting Gibbs-ensemble Monte Carlo simulations with analytical two-body interaction potentials between the molecular species. We calculate several demixing curves at pressures below 150 kbar and temperatures of 1000 K ≤T≤ 2000 K. Despite the approximations introduced by the two-body interaction potentials, our results predict a large miscibility gap in hydrogen-water mixtures at similar conditions as found in experiments. Our findings are in contrast to those from ab initio simulations and provide a renewed indication that hydrogen-water immiscibility regions may have a significant impact on the structure and evolution of ice giant planets like Uranus and Neptune.

Ab initio Gibbs ensemble Monte Carlo simulations of the liquid-vapor equilibrium and the critical point of sodium

The ab initio (ai) Gibbs ensemble (GE) Monte Carlo (MC) method coupled with Kohn-Sham density functional theory is successful in predicting the liquid-vapour equilibrium of insulating systems. Here we show that the aiGEMC method can be used to study also metallic systems, where the excited electronic states play an important role and cannot be neglected. For this we include the electronic free energy in the formulation of the effective energy of the system to be used in the acceptance criteria for the MC moves. The application of this aiGEMC method to sodium yields a good agreement with available experimental data on the liquid-vapour equilibrium densities. We predict a critical point for sodium at 2338 ± 108 K and 0.24 ± 0.03 g cm-3. The liquid structure stemming from aiGEMC simulations is very similar to the one from ab initio molecular dynamics. Since this method can determine phase transition without computing the Gibbs free energy, it may offer a new possibility to study other materials with a reasonable computational cost.

The interiors of Uranus and Neptune: current understanding and open questions

Uranus and Neptune form a distinct class of planets in our Solar System. Given this fact, and ubiquity of similar-mass planets in other planetary systems, it is essential to understand their interior structure and composition. However, there are more open questions regarding these planets than answers. In this review, we concentrate on the things we do not know about the interiors of Uranus and Neptune with a focus on why the planets may be different, rather than the same. We next summarize the knowledge about the planets' internal structure and evolution. Finally, we identify the topics that should be investigated further on the theoretical front as well as required observations from space missions. This article is part of a discussion meeting issue 'Future exploration of ice giant systems'.

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20 000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T=10000  K and P=370  GPa, which is higher in pressure than was inferred with quasiharmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. This may render it difficult to observe this phase with shock experiments because of kinetic effects.

Stability of H3O at extreme conditions and implications for the magnetic fields of Uranus and Neptune

The anomalous nondipolar and nonaxisymmetric magnetic fields of Uranus and Neptune have long challenged conventional views of planetary dynamos. A thin-shell dynamo conjecture captures the observed phenomena but leaves unexplained the fundamental material basis and underlying mechanism. Here we report extensive quantum-mechanical calculations of polymorphism in the hydrogen-oxygen system at the pressures and temperatures of the deep interiors of these ice giant planets (to >600 GPa and 7,000 K). The results reveal the surprising stability of solid and fluid trihydrogen oxide (H₃O) at these extreme conditions. Fluid H₃O is metallic and calculated to be stable near the cores of Uranus and Neptune. As a convecting fluid, the material could give rise to the magnetic field consistent with the thin-shell dynamo model proposed for these planets. H₃O could also be a major component in both solid and superionic forms in other (e.g., nonconvecting) layers. The results thus provide a materials basis for understanding the enigmatic magnetic-field anomalies and other aspects of the interiors of Uranus and Neptune. These findings have direct implications for the internal structure, composition, and dynamos of related exoplanets.

Growth model interpretation of planet size distribution

The radii and orbital periods of 4,000+ confirmed/candidate exoplanets have been precisely measured by the Kepler mission. The radii show a bimodal distribution, with two peaks corresponding to smaller planets (likely rocky) and larger intermediate-size planets, respectively. While only the masses of the planets orbiting the brightest stars can be determined by ground-based spectroscopic observations, these observations allow calculation of their average densities placing constraints on the bulk compositions and internal structures. However, an important question about the composition of planets ranging from 2 to 4 Earth radii (R⊕) still remains. They may either have a rocky core enveloped in a H2-He gaseous envelope (gas dwarfs) or contain a significant amount of multicomponent, H2O-dominated ices/fluids (water worlds). Planets in the mass range of 10-15 M⊕, if half-ice and half-rock by mass, have radii of 2.5 R⊕, which exactly match the second peak of the exoplanet radius bimodal distribution. Any planet in the 2- to 4-R⊕ range requires a gas envelope of at most a few mass percentage points, regardless of the core composition. To resolve the ambiguity of internal compositions, we use a growth model and conduct Monte Carlo simulations to demonstrate that many intermediate-size planets are "water worlds."

First-principles equation of state and shock compression predictions of warm dense hydrocarbons

We use path integral Monte Carlo and density functional molecular dynamics to construct a coherent set of equations of state (EOS) for a series of hydrocarbon materials with various C:H ratios (2:1, 1:1, 2:3, 1:2, and 1:4) over the range of 0.07-22.4gcm^{-3} and 6.7×10^{3}-1.29×10^{8}K. The shock Hugoniot curve derived for each material displays a single compression maximum corresponding to K-shell ionization. For C:H = 1:1, the compression maximum occurs at 4.7-fold of the initial density and we show radiation effects significantly increase the shock compression ratio above 2 Gbar, surpassing relativistic effects. The single-peaked structure of the Hugoniot curves contrasts with previous work on higher-Z plasmas, which exhibit a two-peak structure corresponding to both K- and L-shell ionization. Analysis of the electronic density of states reveals that the change in Hugoniot structure is due to merging of the L-shell eigenstates in carbon, while they remain distinct for higher-Z elements. Finally, we show that the isobaric-isothermal linear mixing rule for carbon and hydrogen EOS is a reasonable approximation with errors better than 1% for stellar-core conditions.