Abstract:First-principles molecular dynamics simulations, which can reveal the microscopic mechanisms and the macroscopic properties of molecules and condensed systems, have become a significant driving force in the development of physics. In the past decades, the ground state electronic properties of various materials can be accurately described using the firstprinciples method based on the Born-Oppenheimer approximation. However, in nonadiabatic processes, such as the ultrafast excited-state dynamics where laser interacts with molecules and condensed matters, the Born-Oppenheimer approximation evolving nuclear wave function on the ground state potential energy plane is no longer valid because of the nonnegligible nonadiabatic couplings between electrons and nuclei. In order to investigate the novel physical phenomena in nonadiabatic processes, many nonadiabatic molecular dynamics methods have been proposed, for instance, full quantum dynamics and mixed quantum-classical dynamics. When the nuclear quantum effects are trivial, the mixed quantum-classical dynamics is effective, which solves the time-dependent electron Schrodinger equation and describes the motion of the nuclei in the form of Newtonian mechanics. Two methods based on the mixed quantum-classical dynamics, the fewest-switches trajectory surface-hopping method (FSSH) and Ehrenfest dynamics, have accepted much attention since they could be easily implemented and combined with real-time time-dependent density functional theory (rt-TDDFT) for high-precision calculations. The combination of quantum and classical dynamics reduces the computational cost and allows the applications for real materials. At low temperatures, the nuclear quantum effects, including quantum tunneling and zero-point vibrations, cannot be ignored. Ring-polymer molecular dynamics (RPMD) based on the imaginary-time path integral is widely used to consider the nuclear quantum effects, which describes the quantum behavior of atomic nuclei by sampling a larger number of quantum configurations and paths. By combining RPMD and mixed quantum-classical nonadiabatic dynamics, one can describe quantum motions of nuclei and the excited electrons in quantum nonadiabatic dynamics simulations. Here, we summarize the recent progress of the rt-TDDFT methods based on numerical atomic orbital basis and planewave basis set, as well as highly efficient RP-TDAP method combining rt-TDDFT and RPMD towards a quantum description of electronic-nuclear dynamics. Finally, we show several representative applications employing these methods, including high harmonic generation modulated by two-color light, charge density wave dynamics, photocatalytic water decomposition as well as the quantum evolution of nuclear wave packets. These applications are in good agreement with experimental measurements, which demonstrates the reliability of these methods and their advantages in the corresponding research areas. These developments and applications are in a significant step from ground state methods to full quantum simulation of the coupled motion of electrons and nuclei, providing a new perspective for understanding and predicting the quantum interactions and dynamical behavior of condensed materials in the degree of atomic spatial scales and attosecond time scales.

Extended Application of Kohn-Sham First-Principles Molecular Dynamics Method with Plane Wave Approximation at High Energy—from Cold Materials to Hot Dense Plasmas

Extended First-Principles Molecular Dynamics Method From Cold Materials to Hot Dense Plasmas

Kinetic Investigations of Nonlinear Electrostatic Excitations in Quantum Plasmas.

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

Link between K-absorption edges and thermodynamic properties of warm-dense plasmas established by improved first-principles method

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

Link Between K Absorption Edges And Thermodynamic Properties Of Warm Dense Plasmas Established By An Improved First-Principles Method

Equation of state of dense plasmas by ab initio simulations: Bridging the gap between quantum molecular dynamics and orbital-free molecular dynamics at high temperature

Development of Various Methods to the Investigation of the Spectral Properties and Collision Dynamics of H‐like Ions Taking Place in Dense and Hot Plasma Environments

Particle-in-cell Simulation Method for Macroscopic Degenerate Plasmas.

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

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

Influence of Plasma Fluctuation on the Electronic Structures and Electron Collisional Excitation and Ionization Processes

Large-scale, Long-Term Nonadiabatic Electron Molecular Dynamics for Describing Material Properties and Phenomena in Extreme Environments.

High Performance Computing for First-Principles Kohn-Sham Density Functional Theory Towards Exascale Supercomputers

Extension of the SpK Atomic Physics Code to Generate Global Equation of State Data

A molecular dynamics framework coupled with smoothed particle hydrodynamics for quantum plasma simulations

Advances in Methods and Applications of Nonadiabatic Quantum Dynamics Simulation of Condensed Matters

Ab Initio Study of Thermodynamically Consistent Equation of State of Warm Dense Aluminum Plasma

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