First-principles Equation of State and Shock Compression Predictions of Warm Dense Hydrocarbons

Abstract

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\times 10^3-1.29\times 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 EOSs is a reasonable approximation with errors better than 1% for stellar-core conditions.

arXiv comment7 pages, 4 figures. Accepted by Physical Review E

References in corpus (20)

• The Interiors of Giant Planets: Models and Outstanding Questions
• Path Integral Monte Carlo Calculation of the Deuterium Hugoniot
• Molecular Dissociation in Hot, Dense Hydrogen
• FPEOS: A First-Principles Equation of State Table of Deuterium for Inertial Confinement Fusion Applications
• A multiphase equation of state for carbon addressing high pressures and temperatures
• All-Electron Path Integral Monte Carlo Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas
• Path Integral Monte Carlo Simulation of the Low-Density Hydrogen Plasma
• First Principles Calculations of Shock Compressed Fluid Helium
• Path Integral Monte Carlo and Density Functional Molecular Dynamics Simulations of Hot, Dense Helium
• Calculation of a Deuterium Double Shock Hugoniot from Ab initio Simulations
• First-Principles Equation of State Calculations of Warm Dense Nitrogen
• First-Principles Equation of State and Electronic Properties of Warm Dense Oxygen
• Simulations of Dense Atomic Hydrogen in the Wigner Crystal Phase
• First-Principles Prediction of the Softening of the Silicon Shock Hugoniot Curve
• Equation of state for partially ionized carbon at high temperatures
• Characterization of the State of Hydrogen
• Thermophysical properties for shock compressed polystyrene
• Hydrogen-Helium Mixtures at High Pressure
• The properties of heavy elements in giant planet envelopes
• Ab initio calculation of thermodynamic, transport, and optical properties of CH${}_2$ plastics

Cited by in corpus (19)

• First-Principles Equation of State Database for Warm Dense Matter Computation
• First principles simulations of dense hydrogen
• SQDFT: Spectral Quadrature method for large-scale parallel $\mathcal{O}(N)$ Kohn-Sham calculations at high temperature
• Ab initio path integral Monte Carlo approach to the momentum distribution of the uniform electron gas at finite temperature without fixed nodes
• Ionization potential depression and Pauli blocking in degenerate plasmas at extreme densities
• Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas
• Equation of state of warm-dense boron nitride combining computation, modeling, and experiment
• Requirements for very high temperature Kohn-Sham density functional simulations and how to bypass them
• Path Integral Monte Carlo and Density Functional Molecular Dynamics Simulations of Warm, Dense MgSiO${}_3$
• Theoretical and experimental investigation of the equation of state of boron plasmas
• Finite-temperature coupled cluster: Efficient implementation and application to prototypical systems
• Magnesium Oxide at Extreme Temperatures and Pressures Studied with First-Principles Simulations
• Model of Ramp Compression of Diamond from Ab Initio Simulations
• Benchmarking boron carbide equation of state using computation and experiment
• Plane-Wave-Based Stochastic-Deterministic Density Functional Theory for Extended Systems
• Nature of the bonded-to-atomic transition in liquid silica to TPa pressures
• First-principles Equations of State and Structures of Liquid Metals in Multi-megabar Conditions
• Pauli Blocking in Degenerate Plasmas and the Separable Potential Approach
• Using density matrix quantum Monte Carlo for calculating exact-on-average energies for ab-initio Hamiltonians in a finite basis set