Abstract:Accurate and efficient integration of the equations of motion is indispensable for molecular dynamics (MD) simulations. Despite the massive use of the conventional leapfrog (LF) integrator in modern computational tools within the framework of MD propagation, further development for better performance is still possible. The alternative version of LF in the middle thermostat scheme (LFmiddle) achieves a higher order of accuracy and efficiency and maintains stable dynamics even with the integration time stepsize extended by several folds. In this work, we perform a benchmark test of the two integrators (LF and LF-middle) in extensive conventional and enhanced sampling simulations, aiming at quantifying the time-stepsizeinduced variations of global properties (e.g., detailed potential energy terms) as well as of local observables (e.g., free energy changes or bondlengths) in practical simulations of complex systems. The test set is composed of six chemically and biologically relevant systems, including the conformational change of dihedral flipping in the N -methylacetamide and an AT (AdenineThymine) tract, the intra-molecular proton transfer inside malonaldehyde, the binding free energy calculations of benzene and phenol targeting T4 lysozyme L99A, the hydroxyl bond variations in ethaline deep eutectic solvent, and the potential energy of the blue-light using flavin photoreceptor. It is observed that the time-step-induced error is smaller for the LFmiddle scheme. The outperformance of LF-middle over the conventional LF integrator is much more significant for global properties than local observables. Overall, the current work demonstrates that the LF-middle scheme should be preferably applied to obtain accurate thermodynamics in the simulation of practical chemical and biological systems.

Comparison of effective and stable Langevin dynamics integrators

Comparison of effective and stable Langevin dynamics integrators

Comparison of Modern Langevin Integrators for Simulations of Coarse-Grained Polymer Melts

Numerical integrators for confined Langevin dynamics

On the Definition of Velocity in Discrete-Time, Stochastic Langevin Simulations

Stable and fast semi-implicit integration of the stochastic Landau-Lifshitz equation

Contraction Rate Estimates of Stochastic Gradient Kinetic Langevin Integrators

Some robust integrators for large time dynamics

Fast and Accurate Numerical Integration of the Langevin Equation with Multiplicative Gaussian White Noise

Comparison of Very Smooth Cell-Model Trajectories Using Five Symplectic and Two Runge-Kutta Integrators

Statistical Error of Numerical Integrators for Underdamped Langevin Dynamics with Deterministic And Stochastic Gradients

Contraction and Convergence Rates for Discretized Kinetic Langevin Dynamics

Extensive Numerical Tests of Leapfrog Integrator in Middle Thermostat Scheme in Molecular Simulations

Numerical Methods with Coordinate Transforms for Efficient Brownian Dynamics Simulations

Kernels and integration cycles in complex Langevin simulations

Stressed and Unstressed Pronouns: Complementary Preferences

Stationary State Distribution and Efficiency Analysis of the Langevin Equation Via Real or Virtual Dynamics.

The velocity jump Langevin process and its splitting scheme: long time convergence and numerical accuracy

Numerics with coordinate transforms for efficient Brownian dynamics simulations