The formation of the giant planets in our solar system, and likely a majority of giant exoplanets, is commonly explained by the accretion of nebular hydrogen and helium onto a large core of terrestrial-like composition. The fate of this core has important consequences for the evolution of the interior structure of the planet. It has recently been shown that H2O, MgO and SiO2 dissolve in liquid metallic hydrogen at high temperature and pressure. In this study, we perform ab initio calculations to study the solubility of an innermost metallic core. We find dissolution of iron to be strongly favored above 2000 K over the entire pressure range (0.4-4 TPa) considered. We compare with and summarize the results for solubilities on other probable core constituents. The calculations imply that giant planet cores are in thermodynamic disequilibrium with surrounding layers, promoting erosion and redistribution of heavy elements. Differences in solubility behavior between iron and rock may influence evolution of interiors, particularly for Saturn-mass planets. Understanding the distribution of iron and other heavy elements in gas giants may be relevant in understanding mass-radius relationships, as well as deviations in transport properties from pure hydrogen-helium mixtures.

Using density functional molecular dynamics free energy calculations, we show that the body-centered-cubic phase of superionic ice previously believed to be the only phase is in fact thermodynamically unstable compared to a novel phase with oxygen positions in fcc lattice sites. The novel phase has a lower proton mobility than the bc phase and may exhibit a higher melting temperature. We predict a transition between the two phases at a pressure of 1 +/- 0.5 Mbar, with potential consequences for the interiors of ice giants such as Uranus and Neptune.

Water and hydrogen at high pressure make up a substantial fraction of the interiors of giant planets. Using ab initio random structure search methods we investigate the ground-state crystal structures of water, hydrogen, and hydrogen-oxygen compounds. We find that, at pressures beyond 14 Mbar, excess hydrogen is incorporated into the ice phase to form a novel structure with H4O stoichiometry. We also predict two new ground state structures, P2_1/m and I4/mmm, for post-C2/m water ice.

Using ab initio simulations we investigate whether water ice is stable in the cores of giant planets, or whether it dissolves into the layer of metallic hydrogen above. By Gibbs free energy calculations we find that for pressures between 10 and 40 Mbar the ice-hydrogen interface is unstable at temperatures above approximately 3000 K, far below the temperature of the core-mantle boundaries in Jupiter and Saturn that are of the order of 10000 K. This implies that the cores of solar and extrasolar giant planets are at least partially eroded.

Gas giants are believed to form by the accretion of hydrogen-helium gas around an initial protocore of rock and ice. The question of whether the rocky parts of the core dissolve into the fluid H-He layers following formation has significant implications for planetary structure and evolution. Here we use ab initio calculations to study rock solubility in fluid hydrogen, choosing MgO as a representative example of planetary rocky materials, and find MgO to be highly soluble in H for temperatures in excess of approximately 10000 K, implying significant redistribution of rocky core material in Jupiter and larger exoplanets.

The Galileo probe showed that Jupiter's atmosphere is severely depleted in neon compared to protosolar values. We show, via ab initio simulations of the partitioning of neon between hydrogen and helium phases, that the observed depletion can be explained by the sequestration of neon into helium-rich droplets within the postulated hydrogen-helium immiscibility layer of the planet's interior. We also demonstrate that this mechanism will not affect argon, explaining the observed lack of depletion of this gas. This provides strong indirect evidence for hydrogen-helium immiscibility in Jupiter.