Abstract:In traditional finite-temperature Kohn–Sham density functional theory (KSDFT), the partial occupation of a large number of high-energy KS eigenstates restricts the use of first-principles molecular dynamics methods at extremely high temperatures. However, stochastic density functional theory (SDFT) can overcome this limitation. Recently, SDFT and the related mixed stochastic–deterministic density functional theory, based on a plane-wave basis set, have been implemented in the first-principles electronic structure software ABACUS [Q. Liu and M. Chen, Phys. Rev. B 106, 125132 (2022)]. In this study, we combine SDFT with the Born–Oppenheimer molecular dynamics method to investigate systems with temperatures ranging from a few tens of eV to 1000 eV. Importantly, we train machine-learning-based interatomic models using the SDFT data and employ these deep potential models to simulate large-scale systems with long trajectories. Subsequently, we compute and analyze the structural properties, dynamic properties, and transport coefficients of warm dense matter.

What problem does this paper attempt to address?

This paper aims to solve the problem of understanding the behavior of materials under extremely high - temperature conditions (such as warm dense matter and hot dense plasma). Specifically, the traditional Kohn - Sham density functional theory (KSDFT) at finite temperature is costly and difficult to implement when dealing with warm dense matter (WDM) because it needs to consider the partial occupation of a large number of high - energy Kohn - Sham (KS) eigenstates. To solve this problem, the paper proposes to combine the Stochastic Density Functional Theory (SDFT) with the Deep Potential Molecular Dynamics (DPMD) method to study systems with a temperature range from tens of electron volts to 1,000 electron volts. ### Main problems 1. **High computational cost**: In warm dense matter, as the temperature rises, the number of electronic states to be considered increases significantly, which makes the traditional KSDFT method infeasible. 2. **Difficult convergence**: For some physical properties, it is difficult to obtain convergent results using a small number of atoms for simulation, unless a large - sized simulation unit and a long - time dynamic trajectory are used. ### Solutions 1. **Combination of SDFT and DPMD**: - **SDFT**: By introducing random orbitals, it overcomes the limitation of partially occupying a large number of high - energy KS eigenstates, thus reducing the computational cost. - **DPMD**: Using the deep potential model generated by machine learning, it can simulate large - scale systems and long - time dynamic trajectories, thus obtaining more accurate physical properties. 2. **Verification and application**: - **Accuracy verification**: By analyzing the statistical error introduced by random orbitals and comparing it with the results of the traditional KSDFT method, the accuracy of SDFT is verified. - **Simulating warm dense boron**: Combining the Born - Oppenheimer Molecular Dynamics (BOMD) method and SDFT to generate training data and construct a deep potential model for simulating warm dense boron. - **Extension of time and space scales**: Through DPMD simulation, the time and space scales of warm dense boron are significantly extended, and convergent physical property data are obtained. ### Experimental results - **Statistical error analysis**: By calculating the root - mean - square error (RMSE) under different numbers of random orbitals and k - points, it is verified that the statistical error of the SDFT method decreases with the increase of the number of random orbitals and k - points. - **Comparison with KSDFT**: Under multiple temperature and density conditions, the differences in total energy, pressure, and ionization degree calculated by the SDFT and KSDFT methods are all less than 0.41%, indicating that SDFT has high precision. - **Force comparison**: The atomic forces predicted by SDFT are highly consistent with the KSDFT results, and the RMSE is less than 0.103 eV/Å, further verifying the reliability of SDFT. - **Density of states (DOS)**: By calculating the density of states of the B and C systems, the accuracy of the SDFT method at different temperatures is verified. ### Conclusion The paper successfully solves the simulation problem of warm dense matter under extremely high - temperature conditions by combining the SDFT and DPMD methods, providing an efficient and accurate computational tool for studying warm dense matter.

Combining stochastic density functional theory with deep potential molecular dynamics to study warm dense matter

Combining stochastic density functional theory with deep potential molecular dynamics to study warm dense matter

Plane-Wave-Based Stochastic-Deterministic Density Functional Theory for Extended Systems

Ab Initio Validation on the Connection Between Atomistic and Hydrodynamic Description to Unravel the Ion Dynamics of Warm Dense Matter

Structure and dynamics of warm dense aluminum: a molecular dynamics study with density functional theory and deep potential

Warm dense matter simulation via electron temperature dependent deep potential molecular dynamics

Stochastic density functional theory combined with Langevin dynamics for warm dense matter

Dynamical structure factors of Warm Dense Matter from Time-Dependent Orbital-Free and Mixed-Stochastic-Deterministic Density Functional Theory

Large Scale and Quantum Accurate Molecular Dynamics Simulation:Liquid Iron under Extreme Condition

Nonlocal Free-Energy Density Functional for a Broad Range of Warm Dense Matter Simulations

Accelerating finite-temperature Kohn-Sham density functional theory with deep neural networks

Accuracy of Kohn-Sham density functional theory for warm- and hot-dense matter equation of state

Deepks: A Comprehensive Data-Driven Approach Toward Chemically Accurate Density Functional Theory

A Time-Dependent Random State Approach for Large-Scale Density Functional Calculations

Accelerating equilibration in first-principles molecular dynamics with orbital-free density functional theory

Nonlocal free-energy density functional for warm dense matter

Ab initio simulation of warm dense matter

Orbital-Free Density Functional Theory: An Attractive Electronic Structure Method for Large-Scale First-Principles Simulations

Understanding intermolecular interactions of large systems in ground state and excited state by using density functional based tight binding methods