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.

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.