Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Abstract

Accurately modeling dense plasmas over wide ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/cm${}^3$ to 16 g/cm${}^3$ and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature (SQ) method. We find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the SESAME database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

References in corpus (17)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Advanced capabilities for materials modelling with Quantum ESPRESSO
• Nearsightedness of Electronic Matter
• All-Electron Path Integral Monte Carlo Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas
• First-Principles Equation of State Database for Warm Dense Matter Computation
• Review of the First Charged-Particle Transport Coefficient Comparison Workshop
• Extended First-Principles Molecular Dynamics Method From Cold Materials to Hot Dense Plasmas
• Carbon ionization at Gbar pressures: an ab initio perspective on astrophysical high-density plasmas
• Fast and Universal Kohn Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma
• Machine Learning Interatomic Potential for Simulations of Carbon at Extreme Conditions
• Pseudoatom molecular dynamics
• Equation of state of warm-dense boron nitride combining computation, modeling, and experiment
• Equation of state for partially ionized carbon at high temperatures
• Pair-distribution functions of two-temperature two-mass systems: Comparison of MD, HNC, CHNC, QMC and Kohn-Sham calculations for dense hydrogen
• Requirements for very high temperature Kohn-Sham density functional simulations and how to bypass them
• Real-space formulation of the stress tensor for $\mathcal{O}(N)$ density functional theory: application to high temperature calculations
• The structure in warm dense carbon

Cited by in corpus (6)

• Electronic Density Response of Warm Dense Matter
• First principles simulations of dense hydrogen
• Abinit 2025: New Capabilities for the Predictive Modeling of Solids and Nanomaterials
• From Density Response to Energy Functionals and Back: An ab initio perspective on Matter Under Extreme Conditions
• Imposing Correct Jellium Response Is Key to Predict the Density Response by Orbital-Free DFT
• Excited States in Warm and Hot Dense Matter