A Versatile Multiple Time Step Scheme for Efficient ab Initio Molecular Dynamics Simulations.

Area of Science:

• Computational Chemistry
• Molecular Dynamics
• Quantum Mechanics

Background:

• Accurate molecular dynamics simulations are crucial for understanding chemical and biological processes.
• Standard methods often face computational limitations, restricting the timescale and system size that can be studied.
• Hybrid Quantum Mechanics/Molecular Mechanics (QM/MM) methods offer a balance but can still be computationally demanding.

Purpose of the Study:

• To implement and validate a time-reversible, multiple time step (MTS) method for Born-Oppenheimer molecular dynamics.
• To enable the use of combined low and high-level electronic structure methods within the MTS framework.
• To achieve significant computational speedup without compromising simulation accuracy.

Main Methods:

• Developed a flexible MTS implementation for full Quantum Mechanics (QM) and QM/MM simulations.
• Utilized a combination of different electronic structure methods (e.g., Density Functional Theory, wave function-based methods) for force calculations.
• Employed a generalized Langevin stochastic thermostat for enhanced stability and accuracy.

Main Results:

• Demonstrated the successful combination of diverse electronic structure methods within the MTS scheme.
• Achieved stable and accurate molecular dynamics trajectories with time steps of several femtoseconds.
• Obtained a 5- to 6-fold overall speedup compared to standard Velocity Verlet integration at an unaltered accuracy level.

Conclusions:

• The presented MTS method significantly reduces the computational cost of accurate ab initio molecular dynamics simulations.
• This approach allows for longer simulation timescales, comparable to classical molecular dynamics.
• The method provides a powerful tool for studying complex systems with high fidelity.