Liouville-von Neumann molecular dynamics.

Area of Science:

• Quantum Chemistry
• Computational Physics
• Materials Science

Background:

• Accurate simulation of molecular dynamics is crucial for understanding chemical reactions and material properties.
• Existing methods like Born-Oppenheimer molecular dynamics face limitations in stability and accuracy for certain systems, particularly those with fractional electron occupations.

Purpose of the Study:

• To present a novel first principles molecular dynamics scheme, Liouville-von Neumann molecular dynamics (LvN MD).
• To combine accurate quantum propagation of electrons with classical nuclei dynamics.
• To overcome limitations of conventional molecular dynamics methods.

Main Methods:

• Utilizes the Liouville-von Neumann equation for density matrix propagation.
• Employs Magnus expansion for the time-evolution operator.
• Involves iterative formation of the Fock operator and integration of the von Neumann equation within each time step.

Main Results:

• The LvN MD scheme is free of constraints and fictitious parameters.
• It avoids the computationally expensive diagonalization of the Fock operator.
• Demonstrates high stability and excellent energy conservation, even for challenging systems.
• Successfully applied to simulations of fullerene formation and the retinal system.

Conclusions:

• LvN MD offers a formally accurate and computationally efficient alternative for quantum molecular dynamics.
• The method's stability and accuracy make it suitable for systems with fractional electron occupations, such as metallic systems.
• This novel approach advances the simulation capabilities in quantum chemistry and condensed matter physics.