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

Abstract

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\times 10^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\times 10^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.

arXiv comment10 pages, 10 figures

References in corpus (26)

• Transparent dense sodium
• What do we (not) know theoretically about solar neutrino fluxes?
• 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
• On the constitution of sodium at higher densities
• A multiphase equation of state for carbon addressing high pressures and temperatures
• The importance of finite-temperature exchange-correlation for warm dense matter calculations
• 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
• Metallic Icosahedron Phase of Sodium at Terapascal Pressures
• Extended First-Principles Molecular Dynamics Method From Cold Materials to Hot Dense Plasmas
• Microscopic origins of the anomalous melting behaviour of high-pressure sodium
• Sodium: a charge-transfer insulator at high pressures
• Fast and accurate quantum molecular dynamics of dense plasmas across temperature regimes
• 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
• Characterization of the State of Hydrogen
• Hydrogen-Helium Mixtures at High Pressure
• Plasmons in Sodium under Pressure: Increasing Departure from Nearly-Free-Electron Behavior
• Ab-initio analysis of plasmon dispersion in sodium under pressure

Cited by in corpus (15)

• Nonempirical Semi-local Free-Energy Density Functional for Matter Under Extreme Conditions
• 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
• A Phaseless Auxiliary-Field Quantum Monte Carlo Perspective on the Uniform Electron Gas at Finite Temperatures: Issues, Observations, and Benchmark Study
• Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas
• Equation of state of warm-dense boron nitride combining computation, modeling, and experiment
• Path Integral Monte Carlo and Density Functional Molecular Dynamics Simulations of Warm, Dense MgSiO${}_3$
• First-Principles Simulations of Warm Dense Lithium Fluoride
• Theoretical and experimental investigation of the equation of state of boron plasmas
• D'yakov-Kontorovitch instability of shock waves in hot plasmas
• Magnesium Oxide at Extreme Temperatures and Pressures Studied with First-Principles Simulations
• Benchmarking boron carbide equation of state using computation and experiment
• Thermal and Pressure Ionization in Warm, Dense MgSiO${}_3$ Studied with First-Principles Computer Simulations
• First-principles Equations of State and Structures of Liquid Metals in Multi-megabar Conditions