Abstract:Molecular dynamics (MD) simulations employing ab initio quantum mechanical and molecular mechanical (ai-QM/MM) potentials are considered to be the state of the art, but the high computational cost associated with the ai-QM calculations remains a theoretical challenge for their routine application. Here, we present a modified protocol of the multiple time step (MTS) method for accelerating ai-QM/MM MD simulations of condensed-phase reactions. Within a previous MTS protocol [Nam J. Chem. Theory Comput. 2014, 10, 4175], reference forces are evaluated using a low-level (semiempirical QM/MM) Hamiltonian and employed at inner time steps to propagate the nuclear motions. Correction forces, which arise from the force differences between high-level (ai-QM/MM) and low-level Hamiltonians, are applied at outer time steps, where the MTS algorithm allows the time-reversible integration of the correction forces. To increase the outer step size, which is bound by the highest-frequency component in the correction forces, the semiempirical QM Hamiltonian is recalibrated in this work to minimize the magnitude of the correction forces. The remaining high-frequency modes, which are mainly bond stretches involving hydrogen atoms, are then removed from the correction forces. When combined with a Langevin or SIN(R) thermostat, the modified MTS-QM/MM scheme remains robust with an up to 8 (with Langevin) or 10 fs (with SIN(R)) outer time step (with 1 fs inner time steps) for the chorismate mutase system. This leads to an over 5-fold speedup over standard ai-QM/MM simulations, without sacrificing the accuracy in the predicted free energy profile of the reaction.

Smoothed Molecular Dynamics for Large Step Time Integration

Molecular Dynamics‐smoothed Molecular Dynamics (MD‐SMD) Adaptive Coupling Method with Seamless Transition

An Improved Smoothed Molecular Dynamics Method by Alternating with Molecular Dynamics

Parallel Smoothed Molecular Dynamics Method and Coupling with Molecular Dynamics

Adaptive Smoothed Molecular Dynamics for Multiscale Modeling

AN ADAPTIVE COUPLING MOLECULAR DYNAMICS-SMOOTHED MOLECULAR DYNAMICS (MD-SMD) METHOD FOR NANOMECHANICS

An improved smoothed molecular dynamics method with high‐order shape function

Seamless Coupling of Molecular Dynamics and Material Point Method Via Smoothed Molecular Dynamics

A critical assessment and improvements on the efficiency of smoothed molecular dynamics method

Accelerating ab initio QM/MM Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semi-empirical QM/MM Hamiltonian

Large-scale simulation of graphene and structural superlubricity with improved smoothed molecular dynamics method

Accelerating Ab Initio Quantum Mechanical and Molecular Mechanical (QM/MM) Molecular Dynamics Simulations with Multiple Time Step Integration and a Recalibrated Semiempirical QM/MM Hamiltonian.

Multiscale Simulation of the Responses of Discrete Nanostructures to Extreme Loading Conditions Based on the Material Point Method

A Concurrent Multiscale Method Based on Smoothed Molecular Dynamics for Large-Scale Parallel Computation at Finite Temperature

Comparative Study of the Water Response to External Force at Nanoscale and Mesoscale

A New Handshaking of Tight-Binding and Molecular Dynamics in Multi-Scale Simulation

Ab Initio Validation on the Connection Between Atomistic and Hydrodynamic Description to Unravel the Ion Dynamics of Warm Dense Matter

Finite Difference Approach for Multiscale Computations of Atomic Chain at Finite Temperature

A Review of Enhanced Sampling Approaches for Accelerated Molecular Dynamics

Coupling the molecular motion and collision processes in numerical simulations

Capturing nanoscale effects by peridynamics