We present a new equation of state (EOS) for dense hydrogen/helium mixtures which covers a range of densities from $10^{-8}$ to $10^6$ g.cm$^{-3}$, pressures from $10^{-9}$ to $10^{13}$ GPa and temperatures from $10^{2}$ to $10^{8}$ K. The calculations combine the EOS of Saumon, Chabrier & vanHorn (1995) in the low density, low temperature molecular/atomic domain, the EOS of Chabrier & Potekhin (1998) in the high-density, high-temperature fully ionized domain, the limits of which differ for H and He, and ab initio quantum molecular dynamics (QMD) calculations in the intermediate density and temperature regime, characteristic of pressure dissociation and ionization. The EOS for the H/He mixture is based on the so-called additive volume law and thus does not take into account the interactions between the two species. A major improvement of the present calculations over existing ones is that we calculate the entropy over the entire density-temperature domain, a necessary quantity for stellar or planetary evolution calculations. The EOS results are compared with existing experimental data, namely Hugoniot shock experiments for pure H and He, and with first principle numerical simulations for both the single elements and the mixture. This new EOS covers a wide range of physical and astrophysical conditions, from jovian planets to solar-type stars, and recovers the existing relativistic EOS at very high densities, in the domains of white dwarfs and neutron stars.

Carbon-hydrogen plasmas and hydrocarbon materials are of broad interest to laser shock experimentalists, high energy density physicists, and astrophysicists. Accurate equations of state (EOS) of hydrocarbons are valuable for various studies from inertial confinement fusion (ICF) to planetary science. By combining path integral Monte Carlo (PIMC) results at high temperatures and density functional theory molecular dynamics (DFT-MD) results at lower temperatures, we compute the EOS for hydrocarbons at 1473 separate ($\rho,T$)-points distributed over a range of compositions. These methods accurately treat electronic excitation and many-body interaction effects and thus provide a benchmark-quality EOS that surpasses that of semi-empirical and Thomas-Fermi-based methods in the warm dense matter regime. By comparing our first-principles EOS to the LEOS 5112 model for CH, we validate the specific heat assumptions in this model but suggest that the Grueneisen parameter is too large at low temperature. Based on our first-principles EOS, we predict the Hugoniot curve of polystyrene to be 2-5% softer at maximum compression than that predicted by orbital-free DFT and SESAME 7593. By investigating the atomic structure and chemical bonding, we show a drastic decrease in the lifetime of chemical bonds in the pressure interval of 0.4-4 megabar. We find the assumption of linear mixing to be valid for describing the EOS and the shock Hugoniot curve of the dense, partially ionized hydrocarbons under consideration. We make predictions of the shock compression of glow-discharge polymers and investigate the effects of oxygen content and C:H ratio on their Hugoniot curve. Our full suite of first-principles simulation results may be used to benchmark future theoretical investigations pertaining to hydrocarbon EOS, and should be helpful in guiding the design of future gigabar experiments.

We use path integral Monte Carlo and density functional molecular dynamics to construct a coherent set of equation of state 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.4$ g/cm$^{3}$ and $6.7\times10^3-1.29\times10^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 EOSs is a reasonable approximation with errors better than 1% for stellar-core conditions.

Recent observations of Jupiter and Saturn provided by spacecraft missions, such as Juno and Cassini, compel us to revise and improve our models of giant planet interiors. Even though hydrogen and helium are by far the dominant species in these planets, heavy elements can play a significant role in the structure and evolution of the planet. For instance, giant-planet cores may be eroded by their surrounding fluid, which would result in a significantly increased concentration of heavy elements in the hydrogen-helium envelope. Furthermore, the heavy elements could inhibit convection by creating a stabilizing gradient of composition. In order to explore the effects of core erosion, we performed ab initio simulations to study structural, diffusion and viscosity properties of dense multi-component mixtures of hydrogen, helium, and silicon dioxide at relevant pressure-temperature conditions. We computed radial distribution functions to identify changes in the chemical behavior of the mixture and to reveal dissociation trends with pressure and temperature. The computed diffusion coefficients of the different species as well as the viscosity provide constraints for the time scale of the dynamics of the core erosion and the mixing of its constituents into the envelope, which will help improve planetary models.

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/cm$^{3}$ and $10^3-1.29\times10^8$ 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\times10^6$ 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 10$^7$ 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.

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.

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 Jupiter's 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, $Z$, 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-$Z$ 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$_2$, we show that the behavior of heavy elements in H-He mixtures may still be represented by an ideal mixture if the effective volumes and internal energies are chosen appropriately. In the case of oxygen, we also compute the effect on the entropy. We find the resulting changes in the temperature-pressure profile to be small. A homogeneous distribution of 2\% oxygen by mass changes the temperature in Jupiter's interior by only 80 K.

We study water-hydrogen mixtures under planetary interior conditions using ab initio molecular dynamics simulations. We determine the thermodynamic properties of various water-hydrogen mixing ratios at temperatures of 2000 and 6000 K for pressures of a few tens of GPa. These conditions are relevant for ice giant planets and for the outer envelope of the gas giants. We find that at 2000 K the mixture is in a molecular regime, while at 6000 K the dissociation of hydrogen and water is important and affects the thermodynamic properties. We study the structure of the liquid and analyze the radial distribution function. We provide estimates for the transport properties, diffusion and viscosity, based on autocorrelation functions. We obtained viscosity estimates of the order of a few tenths of mPa.s for the conditions under consideration. These results are relevant for dynamo simulations of ice 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 a substantial fraction of water and hydrogen in giant planets may occur in homogeneously mixed form rather than in separate layers. The extend 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.