First-Principles Prediction of the Softening of the Silicon Shock Hugoniot Curve

Abstract

Shock compression of silicon (Si) under extremely high pressures (>100 Mbar) was investigated by using two first-principles methods of orbital-free molecular dynamics (OFMD) and path integral Monte Carlo (PIMC). While pressures from the two methods agree very well, PIMC predicts a second compression maximum because of 1s electron ionization that is absent in OFMD calculations since Thomas-Fermi-based theories lack shell structure. The Kohn-Sham density functional theory is used to calculate the equation of state (EOS) of warm dense silicon for low-pressure loadings (P < 100 Mbar). Combining these first-principles EOS results, the principal shock Hugoniot curve of silicon for pressures varying from 1 Mbar to above 10 Gbar was derived. We find that silicon is 20% or more softer than what was predicted by widely-used EOS models. Existing high-pressure experimental data (P = 1 - 2 Mbar) seem to indicate this softening behavior of Si, which calls for future strong-shock experiments (P > 10 Mbar) to benchmark our results.

30 pages, 6 figures

References in corpus (2)

• All-Electron Path Integral Monte Carlo Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas
• First Principles Calculations of Shock Compressed Fluid Helium

Cited by in corpus (7)

• First-Principles Equation of State Database for Warm Dense Matter Computation
• First-principles Equation of State and Shock Compression Predictions of Warm Dense Hydrocarbons
• Equation of state of warm-dense boron nitride combining computation, modeling, and experiment
• Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas
• Path Integral Monte Carlo and Density Functional Molecular Dynamics Simulations of Warm, Dense MgSiO${}_3$
• Magnesium Oxide at Extreme Temperatures and Pressures Studied with First-Principles Simulations
• Atom-in-jellium equations of state in the high energy density regime