Abstract:An important challenge for modern molecular simulations is to develop robust algorithms for molecular dynamics (MD), which employ a time interval (i.e., step size) of evolution as large as possible to achieve high efficiency while maintaining accuracy for sampling the distribution of configuration and momentum, e.g., for the canonical ensemble where the number of particles (N), the volume (V), and the temperature (T) are constant. Starting from the velocity-Verlet algorithm, in the recent paper we have proposed a unified scheme—the middle scheme for constructing efficient thermostatting integrators for configurational sampling of the canonical ensemble. In the present paper we extend the middle scheme to derive new thermostatting integrators based on the leap-frog algorithm. The leap-frog middle scheme is able to recover the exact marginal distribution of configuration and that of momentum in the harmonic limit, regardless of the time interval as long as the evolution of the MD trajectory is stable. It is also proved that the leap-frog middle scheme leads to the same marginal distribution of configuration as the velocity- Verlet middle scheme does even for general anharmonic systems. As illustrated in numerical examples (such as a 1-dimensional quartic potential, liquid water at room temperature with or without bond-length constraints, etc.), while both middle schemes perform better than conventional approaches (e.g., the side scheme and the default Langevin dynamics method in AMBER) and can calculate configuration-dependent properties with the same accuracy by using a larger time interval for from at least 4 times to even an order of magnitude, the leap-frog middle scheme produces a more accurate marginal distribution of momentum than that the velocity-Verlet middle scheme does. The sampling efficiency analysis by employing the characteristic correlation time step number (of the potential or of the Hamiltonian) suggests that both middle schemes share the same efficiency in sampling the configuration space and lead to similar efficiency in sampling the whole phase space (i.e., both configuration and momentum). When both only-configuration- dependent and only-momentum-dependent physical observables are considered, we suggest that the leap-frog middle scheme should in principle offer a better approach. When we study properties involving cross-terms of configuration and momentum, e.g., the fluctuation of the Hamiltonian (related to the heat capacity) and the energy current and its fluctuation (related to the thermal conductivity), analytic results in the harmonic limit and numerical data for anharmonic systems indicate that which of the two middle schemes performs better depends on the system and on the specific property. Note that the leap-frog middle scheme is of first-order accuracy for sampling the whole phase space while the velocity-Verlet middle scheme is of second-order accuracy. It is reasonable to choose the velocity-Verlet middle scheme to perform simulations for such properties for general molecular systems. Because either the leap-frog or the velocity-Verlet algorithm is the default integrator in molecular dynamics simulation software packages, the two unified middle schemes can offer novel, simple, accurate, and efficient thermostatting integrators for performing MD (or path integral MD (PIMD)) for the classical (or quantum) canonical ensemble.

A Reversible Unwrapping Algorithm for Constant Pressure Molecular Dynamics Simulations

Unwrapping NPT simulations to calculate diffusion coefficients

Periodic Boundary Conditions for Long-time Nonequilibrium Molecular Dynamics Simulations of Incompressible Flows

A Method for Molecular Dynamics on Curved Surfaces

Constant-temperature molecular-dynamics algorithms for mixed hard-core/continuous potentials

Effective conservation of energy and momentum algorithm using switching potentials suitable for molecular dynamics simulation of thermodynamical systems

Path Integral Molecular Dynamics within the Grand Canonical-like Adaptive Resolution Technique: Simulation of Liquid Water

Adaptive Resolution Force Probe Simulations: Coarse Graining in the Ideal Gas Approximation

Constant pressure and temperature discrete-time Langevin molecular dynamics

Asymmetric periodic boundary conditions for molecular dynamics and coarse-grained simulations of nucleic acids

Discontinuous Molecular Dynamics for Rigid Bodies: Applications

Interpolating many-body wave functions for accelerated molecular dynamics on the near-exact electronic surface

Virtual Molecular Dynamics

One-dimensional moving window atomistic framework to model long-time shock wave propagation

A Leap-Frog Algorithm-Based Efficient Unified Thermostat Scheme for Molecular Dynamics

Fluctuating hydrodynamics for multiscale simulation of inhomogeneous fluids: mapping all-atom molecular dynamics to capillary waves

Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques

A local resampling trick for focused molecular dynamics

A Scalable, Linear-Time Dynamic Cutoff Algorithm for Molecular Simulations of Interfacial Systems

Smoothed Molecular Dynamics for Large Step Time Integration

Hybrid Particle-Field Molecular Dynamics Under Constant Pressure