Mail:
Dept. of Chemistry
Ohio State University
100 W. 18th Ave.
Columbus, OH 43210
Office:
412 CBEC
Email:
herbert@
chemistry.ohio-state.edu
Spin density of the aqueous electron, e^{–}(aq), from an AIMD simulation |
Efficient algorithms developed in our group are expanding the frontiers of ab initio molecular dynamics (AIMD) simulations, in which on-the-fly quantum chemistry is used to calculate the forces required for a classical molecular dynamics simulation, thus obviating the need for parametrized force fields. This is especially important for unusual or reactive species (such as the hydrated electron and other aqueous radicals), for which it is extremely difficult to construct accurate force fields. Improvements to the efficiency of AIMD algorithms translate directly into longer simulation time scales and larger system sizes, for a fixed quantity of computer resources. Our Fock matrix extrapolation algorithm is amongst the preferred ways to accelerate AIMD. We have also developed the QM/MM-Ewald method for incorporating periodic boundary conditions into QM/MM calculations, in a manner that is robust in arbitrary basis sets.
Monitoring energy conservation in a fragment-based AIMD simulation. Only with rigorously correct gradients is energy conserved. |
Especially for a system such as liquid water, which would seem to divide up into small fragments in a natural way, there is significant interest in using fragment-based quantum chemistry approaches to accelerate AIMD calculations, and numerous algorithms have been proposed. Many of these represent a form of QM/MM simulation in which fragment-based wave function calculations are used to obtain classical MM force-field parameters "on the fly". This form of classical embedding, however, can significantly complicate the analytic gradient formalism for the underlying QM method, and many incorrect analytic gradients have been reported in the literature. Energy conservation in classical MD simulations is known to be exquisitely sensitive to the accuracy and consistency of the forces, so many of these fragment-based approaches are fundamentally ill-posed for dynamics applications. We have recently reported a variational formalism for fragment-based QM/MM embedding, which greatly simplifies the formulation of analytic gradients, without introducing approximations, providing for energy-conserving simulations based on fragmentation.