François Soubiran

Understanding Jupiter's interior

This article provides an overview of how models of giant planet interiors are constructed. We review measurements from past space missions that provide constraints for the interior structure of Jupiter. We discuss typical three-layer interior models that consist of a dense central core and an inner metallic and an outer molecular hydrogen-helium layer. These models rely heavily on experiments, analytical theory, and first-principle computer simulations of hydrogen and helium to understand their behavior up to the extreme pressures ∼10 Mbar and temperatures ∼10 000 K. We review the various equations of state used in Jupiter models and compare them with shock wave experiments. We discuss the possibility of helium rain, core erosion and double diffusive convection may have important consequences for the structure and evolution of giant planets. In July 2016 the Juno spacecraft entered orbit around Jupiter, promising high-precision measurements of the gravitational field that will allow us to test our understanding of gas giant interiors better than ever before.

MISCIBILITY CALCULATIONS FOR WATER AND HYDROGEN IN GIANT PLANETS

We present results from ab initio simulations of liquid water?hydrogen mixtures in the range from 2 to 70 GPa and from 1000 to 6000 K, covering conditions in the interiors of ice giant planets and parts of the outer envelope of gas giant planets. In addition to computing the pressure and the internal energy, we derive the Gibbs free energy by performing a thermodynamic integration. For all conditions under consideration, our simulations predict hydrogen and water to mix in all proportions. The thermodynamic behavior of the mixture can be well described with an ideal mixing approximation. We suggest that a substantial fraction of water and hydrogen in giant planets may occur in homogeneously mixed form rather than in separate layers. The extent?of mixing depends on the planet?s interior dynamics and its conditions of formation, in particular on how much hydrogen was present when icy planetesimals were delivered. Based on our results, we do not predict water?hydrogen mixtures to phase separate during any stage of the evolution of giant planets. We also show that the hydrogen content of an exoplanet is much higher if the mixed interior is assumed.

First-principles equation of state and electronic properties of warm dense oxygen

We perform all-electron path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of oxygen. Our simulations cover a wide density-temperature range of 1-100 g cm(-3) and 10(4)-10(9) K. By combining results from PIMC and DFT-MD, we are able to compute pressures and internal energies from first-principles at all temperatures and provide a coherent equation of state. We compare our first-principles calculations with analytic equations of state, which tend to agree for temperatures above 8 × 10(6) K. Pair-correlation functions and the electronic density of states reveal an evolving plasma structure and ionization process that is driven by temperature and density. As we increase the density at constant temperature, we find that the ionization fraction of the 1s state decreases while the other electronic states move towards the continuum. Finally, the computed shock Hugoniot curves show an increase in compression as the first and second shells are ionized.

Equation of state and shock compression of warm dense sodium—A first-principles study

As one of the simple alkali metals, sodium has been of fundamental interest for shock physics experiments, but knowledge of its equation of state (EOS) in hot, dense regimes is not well known. By combining path integral Monte Carlo (PIMC) results for partially ionized states [B. Militzer and K. P. Driver, Phys. Rev. Lett. 115, 176403 (2015)] at high temperatures and density functional theory molecular dynamics (DFT-MD) results at lower temperatures, we have constructed a coherent equation of state for sodium over a wide density-temperature range of 1.93-11.60 g/cm3 and 103-1.29×108 K. We find that a localized, Hartree-Fock nodal structure in PIMC yields pressures and internal energies that are consistent with DFT-MD at intermediate temperatures of 2×106 K. Since PIMC and DFT-MD provide a first-principles treatment of electron shell and excitation effects, we are able to identify two compression maxima in the shock Hugoniot curve corresponding to K-shell and L-shell ionization. Our Hugoniot curves provide a benchmark for widely used EOS models: SESAME, LEOS, and Purgatorio. Due to the low ambient density, sodium has an unusually high first compression maximum along the shock Hugoniot curve. At beyond 107 K, we show that the radiation effect leads to very high compression along the Hugoniot curve, surpassing relativistic corrections, and observe an increasing deviation of the shock and particle velocities from a linear relation. We also compute the temperature-density dependence of thermal and pressure ionization processes.

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.

THE PROPERTIES OF HEAVY ELEMENTS IN GIANT PLANET ENVELOPES

The core accretion model for giant planet formation suggests a two layer picture for the initial structure of Jovian planets, with heavy elements in a dense core and a thick H-He envelope. Late planetesimal accretion and core erosion could potentially enrich the H-He envelope in heavy elements, which is supported by the three-fold solar metallicity that was measured in Jupiters atmosphere by the Galileo entry probe. In order to reproduce the observed gravitational moments of Jupiter and Saturn, models for their interiors include heavy elements,

Properties of hydrogen, helium, and silicon dioxide mixtures in giant planet interiors

Z

Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas

, in various proportions. However, their effect on the equation of state of the hydrogen-helium mixtures has not been investigated beyond the ideal mixing approximation. In this article, we report results from \textit{ab initio} simulations of fully interacting H-He-

Hydrogen–water mixtures in giant planet interiors studied with ab initio simulations

Z

Electrical conductivity and magnetic dynamos in magma oceans of Super-Earths

mixtures in order to characterize their equation of state and to analyze possible consequences for the interior structure and evolution of giant planets. Considering C, N, O, Si, Fe, MgO and SiO