Efficient and Accurate Born-Oppenheimer Molecular Dynamics for Large Molecular Systems

Area of Science:

• Computational Chemistry
• Molecular Dynamics
• Spectroscopy

Background:

• Born-Oppenheimer molecular dynamics (BOMD) is crucial for simulating molecular behavior.
• Accurate and efficient methods are needed for large molecular systems and complex simulations.
• Calculating vibrational spectra and thermodynamics requires robust computational tools.

Purpose of the Study:

• To introduce an efficient computational scheme for Born-Oppenheimer molecular dynamics (BOMD) simulations.
• To enable high-quality vibrational spectra and thermodynamic calculations for large molecular systems.
• To explore the parallel performance and scalability of the implemented method.

Main Methods:

• Combines the corrected small basis set Hartree-Fock (HF-3c) method with extended Lagrangian BOMD (XL-BOMD).
• Utilizes graphics processing units (GPUs) for accelerated calculation of two-electron integrals.
• Implements a strong scaling parallelization strategy for performance analysis.

Main Results:

• Successfully computed high-quality vibrational spectra for beta-carotene, paclitaxel, and liquid water (up to 500 atoms).
• Demonstrated the method's capability for large molecular systems, including those with explicit solvent molecules.
• Validated the computational efficiency and reliability of the new BOMD scheme.

Conclusions:

• The presented BOMD scheme is a cost-efficient and reliable tool for vibrational spectra and thermodynamics of large molecular systems.
• The method is suitable for systems exceeding 500 atoms, including explicit solvent.
• The computational approach offers significant advantages for complex molecular simulations.