Related papers: Optimized Verlet-like algorithms for molecular dyn…
In ab initio molecular dynamics simulations of real-world problems, the simple Verlet method is still widely used for integrating the equations of motion, while more efficient algorithms are routinely used in classical molecular dynamics.…
We present a new molecular-dynamics algorithm for integrating the equations of motion for a system of particles interacting with mixed continuous/impulsive forces. This method, which we call Impulsive Verlet, is constructed using operator…
The Discrete Element Method is widely employed for simulating granular flows, but conventional integration techniques may produce unphysical results for simulations with static friction when particle size ratios exceed $R \approx 3$. These…
We present a revision to the well known Stormer-Verlet algorithm for simulating second order differential equations. The revision addresses the inclusion of linear friction with associated stochastic noise, and we analytically demonstrate…
A new methodology is developed to integrate numerically the equations of motion for classical many-body systems in molecular dynamics simulations. Its distinguishable feature is the possibility to preserve, independently on the size of the…
Molecular simulations of many particles which move rather according to a brownian than a newtonian type of dynamics, nevertheless, can be performed by means of a "velocity-Verlet-like" algorithm. The derivation of this algorithm requires…
A new scheme for numerical integration of motion for classical systems composed of rigid polyatomic molecules is proposed. The scheme is based on a matrix representation of the rotational degrees of freedom. The equations of motion are…
The Verlet method is still widely used to integrate the equations of motion in ab initio molecular dynamics simulations. We show that the stability limit of the Verlet method may be significantly increased by setting an upper limit on the…
We have presented some practical consequences on the molecular-dynamics simulations arising from the numerical algorithm published recently in paper Int. J. Mod. Phys. C 16, 413 (2005). The algorithm is not a finite-difference method and…
An improved neighbor list algorithm is proposed to reduce unnecessary interatomic distance calculations in molecular simulations. It combines the advantages of Verlet table and cell linked list algorithms by using cell decomposition…
A revised version of the quaternion approach for numerical integration of the equations of motion for rigid polyatomic molecules is proposed. The modified approach is based on a formulation of the quaternion dynamics with constraints. This…
A new algorithm for numerical integration of the rigid-body equations of motion is proposed. The algorithm uses the leapfrog scheme and the quantities involved are angular velocities and orientational variables which can be expressed in…
A new approach is developed to integrate numerically the equations of motion for systems of interacting rigid polyatomic molecules. With the aid of a leapfrog framework, we directly involve principal angular velocities into the integration,…
A recent article in J. Chem. Phys. argues that the two algorithms, the velocity-Verlet, and position-Verlet integrators, commonly used in Molecular Dynamics (MD) simulations, are different \cite{Ni2024}. But not only are the two algorithms…
Simulation of many-particle system evolution by molecular dynamics takes to decrease integration step to provide numerical scheme stability on the sufficiently large time interval. It leads to a significant increase of the volume of…
We present the complete set of stochastic Verlet-type algorithms that can provide correct statistical measures for both configurational and kinetic sampling in discrete-time Langevin systems. The approach is a brute-force general…
Simulations of objects with classical dynamics are in fact a particular version of discrete dynamics since almost all the classical dynamics simulations in natural science are performed with the use of the simple ''Leapfrog" or ''Verlet"…
An algorithm for a family of self-starting high-order implicit time integration schemes with controllable numerical dissipation is proposed for both linear and nonlinear transient problems. This work builds on the previous works of the…
An approach is proposed to improve the efficiency of fourth-order algorithms for numerical integration of the equations of motion in molecular dynamics simulations. The approach is based on an extension of the decomposition scheme by…
A variational formulation is introduced for the Oseen equations written in terms of vor\-ti\-city and Bernoulli pressure. The velocity is fully decoupled using the momentum balance equation, and it is later recovered by a post-process. A…