Semi-Empirical Born–Oppenheimer Molecular Dynamics (SEBOMD) within the Amber Biomolecular Package

Semi-empirical quantum methods from the neglect of differential diatomic overlap (NDDO) family such as MNDO, AM1, or PM3 are fast albeit approximate quantum methods. By combining them with linear scaling methods like the divide & conquer (D&C) method, it is possible to quickly evaluate the e...

Full description

Saved in:
Bibliographic Details
Published in:Journal of chemical information and modeling 2019-01, Vol.59 (1), p.206-214
Main Authors: Marion, Antoine, Gokcan, Hatice, Monard, Gerald
Format: Article
Language:English
Subjects:
Approximation
Boundary conditions
Chemical Sciences
Convergence
Embedding
Energy conservation
Molecular dynamics
or physical chemistry
Potential energy
Quantum mechanics
Self consistent fields
Simulation
Theoretical and
Water
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:Semi-empirical quantum methods from the neglect of differential diatomic overlap (NDDO) family such as MNDO, AM1, or PM3 are fast albeit approximate quantum methods. By combining them with linear scaling methods like the divide & conquer (D&C) method, it is possible to quickly evaluate the energy of systems containing hundreds to thousands of atoms. We here present our implementation in the Amber biomolecular package of a SEBOMD module that provides a way to run semi-empirical Born–Oppenheimer molecular dynamics. At each step of a SEBOMD, a fully converged self-consistent field (SCF) calculation is performed to obtain the semiempirical quantum potential energy of a molecular system encaged or not in periodic boundary conditions. We describe the implementation and the features of our SEBOMD implementation. We show the requirements to conserve the total energy in NVE simulations, and how to accelerate SCF convergence through density matrix extrapolation. Specific ways of handling periodic boundary conditions using mechanical embedding or electrostatic embedding through a tailored quantum Ewald summation is developed. The parallel performance of SEBOMD simulations using the D&C scheme are presented for liquid water systems of various sizes, and a comparison between the traditional full diagonalization scheme and the D&C approach for the reproduction of the structure of liquid water illustrates the potentiality of SEBOMD to simulate molecular systems containing several hundreds of atoms for hundreds of picoseconds with a quantum mechanical potential in a reasonable amount of CPU time.
ISSN:1549-9596
1549-960X
DOI:10.1021/acs.jcim.8b00605