Verlet Velocity Algorithm

The velocity Verlet algorithm may be derived by considering a standard approximate decomposition of the Liouville operator which preserves reversibility and is symplectic (which implies that volume in phase space is conserved). This approach [47] has had several beneficial consequences. 

A straightforward derivation (not reproduced here) shows that the effect of the diree successive steps embodied in equation (b3.3.7), with the above choice of operators, is precisely the velocity Verlet algorithm. This approach is particularly usefiil for generating multiple time-step methods. 

Each of these operators is unitary U —t) = U t). Updating a time step with the propagator Uf( At)U At)Uf At) yields the velocity-Verlet algorithm. Concatenating the force operator for successive steps yields the leapfrog algorithm ... 

There is a number of algorithms to solve equations (1) and (2) that differ appreciably in their properties which are beyond the scope of the present article. In the discussion below we use the velocity Verlet algorithm. However, better approaches can be employed [2-5]. We define a rule - F X t), At) that modifies X t) to X t + At) and repeat the application of this rule as desired. For example the velocity Verlet algorithm ( rule ) is ... 

Our discussion so far has considered the use of SHAKE with the Verlet algorithm Versions have also been derived for other integration schemes, such as the leap-froj algorithm, the predictor-corrector methods and the velocity Verlet algorithm. In the cast of the velocity Verlet algorithm, the method has been named RATTLE [Anderson 1983]... 

The equations of motion are integrated using a modified velocity Verlet algorithm. The modification is required because the force depends upon the velocity the extra step involves... 

Dissipative particle dynamics (DPD) is a technique for simulating the motion of mesoscale beads. The technique is superficially similar to a Brownian dynamics simulation in that it incorporates equations of motion, a dissipative (random) force, and a viscous drag between moving beads. However, the simulation uses a modified velocity Verlet algorithm to ensure that total momentum and force symmetries are conserved. This results in a simulation that obeys the Navier-Stokes equations and can thus predict flow. In order to set up these equations, there must be parameters to describe the interaction between beads, dissipative force, and drag. 

If the magnitudes of the dissipative force, random noise, or the time step are too large, the modified velocity Verlet algorithm will not correctly integrate the equations of motion and thus give incorrect results. The values that are valid depend on the particle sizes being used. A system of reduced units can be defined in which these limits remain constant. 

As j approaches zero, this method reduces to the velocity Verlet algorithm ... 

An even better handling of the velocities is obtained by another variant of the basic Verlet integrator, known as the velocity Verlet algorithm. This is a Verlet-type algorithm that stores positions, velocities, and accelerations all at the same time t and minimizes roundoff errors [14]. The velocity Verlet algorithm is written... 

Again, elimination of the velocities from these equations recovers the Verlet algorithm. In practice, the velocity Verlet algorithm consists of the following steps ... 

In this expression m5, V and a are functions of x and can computed at f and t+At. The velocity at the half-steps is directly provided by the velocity Verlet algorithm. To summarize, to calculate d(m )/df at time t + At/2 we need to collect data from time steps f - At, f, f + At, and f + 2At. Recall again that these equations are used to calculate the mean force along , not to advance the system in time. 

On-the-fly molecular dynamics have been employed in order to simulate the photochemistry of carbonyl-containing compounds. The on-the-fly mechanism implemented in the MNDO program is the velocity-Verlet algorithm. Here an additional aspect of the usage of a computational cheap semiempirical method is visible. In order to provide realistic relative yields of different photochemical reactions, a large enough sample of trajectories is needed. For these systems, a substantial amount of trajectories (around 100) has been calculated for a relatively long timescale (up to 100 ps). 

A very simple implementation of a molecular dynamics trajectory calculation is achieved by using a velocity Verlet algorithm to calculate the... 

This projected velocity Verlet algorithm has been found to be an efficient and simple minimization algorithm for many of the methods discussed here. 

Our AIMD simulations are all-electron and self-consistent at each 0.4 femtoseconds (fs) time step. Variational fitting ensures accurate forces for any finite orbital or fitting basis sets and any finite numerical grid. These forces are used to propagate the nuclear motion according to the velocity Verlet algorithm [22]. The accuracy of these methods is indicated by the fact that during the 500 time-step simulations of methyl iodide dissociation described below, the center of mass moved by less than 10-6 A. 

The equations of motion of the mers can be integrated using a variety of techniques. In the work presented here, I used a velocity-Verlet algorithm [96]. The time step At depends on the range of interaction, the temperature and the viscous damping. In the purely repulsive case, with no shear and Y=0.5r, one can... 

A full description of the outline and the capabilities of the NEWTON-X package is given elsewhere [38], In brief, the nuclear motion is represented by classical trajectories computed by numerical integration of Newton s equations using the Velocity-Verlet algorithm [39], Temperature influence can be added by means of the Andersen thermostat [40], The molecule is considered to be in some specific... 

A second test is to see if varying the time step At used in the simulation leads to the expected change in A . The velocity Verlet algorithm is a second-order algorithm, and a decrease of the time step by a factor of 2 should lead to a decrease of the A value by a factor of 4.13,127... 

The careful reader should have realized that we choose not to break up this operator with another Trotter factorization, as was done for the extended system case. In practice, one does not multiple-time-step the modified velocity Verlet algorithm because it will, in general, have a unit Jacobian. Thus, one would like the best representation of the operator that can be obtained in closed form. However, even in the case of a modified velocity Verlet operator that has a nonunit Jacobian, multiple-time-stepping this procedure can be costly because of the multiple force evaluations. Generally, if the integrator is stable without multiple-time-step procedures, avoid them. The solution to this first-order inhomogeneous differential equation is standard and can be found in texts on differential equations (see, e.g.. Ref. 53). 

Our primary goal was the simulation of entire atomic systems, thus made of electrons and nuclei. As mentioned earlier (see Sect. 2.3), in a large class of systems (e.g. not too high temperature) one can decouple the motion of nuclei and electrons within the Born-Oppenheimer approximation. The previous section was then devoted to the Density Functional Theory solution of the electronic structure problem at fixed ionic positions. By computing the Hellmann-Feynman forces (11) we can now propagate the dynamics of an ensemble of (classical) nuclei as described in Sect. 2.3, using e.g. the velocity verlet algorithm [117]. 

To integrate these equations of motion [173,174] one can use a velocity verlet algorithm [117] propagating similarly the nuclei positions and the electronic degrees of freedom Cj(G). The orthonormality conditions can be taken into account by the SHAKE and RATTLE method [173, 175, 176]. Constant ionic temperature can also be imposed through the use of thermostats [169,170,173,177]. 

There are various, essentially equivalent, versions of the Verlet algorithm, including the original method employed by Verlet [13,44] in his investigations of the properties of the Leimard-Jones fluid, and a leapfrog form [45]. Here we concentrate on the velocity Verlet algorithm [46], which may be written... 

---