Requirements for very high temperature Kohn-Sham density functional simulations and how to bypass them

Abstract

In high temperature density functional theory simulations (from tens of eV to keV) the total number of Kohn-Sham orbitals is a critical quantity to get accurate results. To establish the relationship between the number of orbitals and the level of occupation of the highest orbital, we derived a model based on the electron gas properties at finite temperature. This model predicts the total number of orbitals required to reach a given level of occupation and thus a stipulated precision. Levels of occupation as low as 10-4, and below, must be considered to get converged results better than 1%, making high temperature simulations very time consuming beyond a few tens of eV. After assessing the predictions of the model against previous results and ABINIT minimizations, we show how the extended FPMD method of Zhang et al. [PoP 23 042707, 2016] allows to bypass these strong constraints on the number of orbitals at high temperature.

arXiv comment7 pages, 7 figures

References in corpus (19)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Nearsightedness of Electronic Matter
• Equation of state of fully ionized electron-ion plasmas
• Issues and Challenges in Orbital-free Density Functional Calculations
• Non-empirical Generalized Gradient Approximation Free Energy Functional for Orbital-free Simulations
• Review of the First Charged-Particle Transport Coefficient Comparison Workshop
• Extended First-Principles Molecular Dynamics Method From Cold Materials to Hot Dense Plasmas
• Stochastic Density Functional Theory at Finite Temperatures
• Fast and accurate quantum molecular dynamics of dense plasmas across temperature regimes
• 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
• SQDFT: Spectral Quadrature method for large-scale parallel $\mathcal{O}(N)$ Kohn-Sham calculations at high temperature
• First-principles Equation of State and Shock Compression Predictions of Warm Dense Hydrocarbons
• Ionic and electronic transport properties in dense plasmas by orbital-free density functional theory
• Towards accurate orbital-free simulations: a generalized gradient approximation for the non-interacting free energy density functional
• Comparison of Density Functional Approximations and the Finite-temperature Hartree-Fock Approximation in Warm Dense Lithium
• Static and dynamic properties of multi-ionic plasma mixtures
• Orbital-free extension to Kohn-Sham density functional theory equation of state calculations: application to silicon dioxide
• Magnesium Oxide at Extreme Temperatures and Pressures Studied with First-Principles Simulations

Cited by in corpus (10)

• First principles simulations of dense hydrogen
• Benchmarking Exchange-Correlation Functionals in the Spin-Polarized Inhomogeneous Electron Gas under Warm Dense Conditions
• Abinit 2025: New Capabilities for the Predictive Modeling of Solids and Nanomaterials
• Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics
• Mixed Stochastic-Deterministic Time-Dependent Density Functional Theory: Application to Stopping Power of Warm Dense Carbon
• Stochastic and Mixed Density Functional Theory within the projector augmented wave formalism for the simulation of warm dense matter
• Plane-Wave-Based Stochastic-Deterministic Density Functional Theory for Extended Systems
• Effect of Ionic Disorder on the Principal Shock Hugoniot
• Nonlocal free-energy density functional for warm dense matter
• Random-walk shielding-potential viscosity model for warm dense metals