Respa

Reversible Reference System Propagator Algorithms (r-RESPA) 299... 

In many cases the dynamical system consists of fast degrees of freedom, labeled x, and slow degrees of freedom, labeled y. An example is that of a fluid containing polyatomic molecules. The internal vibrations of the molecules are often very fast compared to their translational and orientational motions. Although this and other systems, like proteins, have already been treated using RESPA,[17, 34, 22, 23, 24, 25, 26] another example, and the one we focus on here, is that of a system of very light particles (of mass m) dissolved in a bath of very heavy particles (mass M).[14] The positions of the heavy particles are denoted y and the positions of the light particles rire denoted by X. In this case the total Liouvillian of the system is ... 

Another immediate application of r-RESPA is to the case when the force can be subdivided into a short range part and a long range part. One way for effectuating this break up is to introduce a switching function, s x) that is unity at short inter-particle separations and 0 at large inter-particle separations. We introduced this strategy in our earlier non-reversible RESPA paper [15] where we expressed the total force as. 

Reversible Reference System Propagator Algorithms (r-RESPA) 307 Thus the propagator in Eq. (27) produces the following dynamics algorithm ... 

Fig. 1. CPU times (in hours) for 1 ps MD runs for various proteins using three different methods, direct velocity Verlet with a time-step 0.5 fs, r-RESPA with direct evaluation of electrostatic forces and an overall time-step of 4.0 fs, and r-RESPA/TFMM with an overall time-step 4.0 fs (combination of (2,2,2,2) in force breakup).The energy conservation parameter log AE for the three methods are comparable. The CPU time (hours) is for RISC6000 /MODEL 590 computer.

Since many systems of interest in chemistry have intrinsic multiple time scales it is important to use integrators that deal efficiently with the multiple time scale problem. Since our multiple time step algorithm, the so-called reversible Reference System Propagator Algorithm (r-RESPA) [17, 24, 18, 26] is time reversible and symplectic, they are very useful in combination with HMC for constant temperature simulations of large protein systems. 

This article is organized as follows Sect. 2 explains why it seems important to use symplectic integrators. Sect. 3 describes the Verlet-I/r-RESPA impulse MTS method, Sect. 4 presents the 5 femtosecond time step barrier. Sect. 5 introduce a possible solution termed the mollified impulse method (MOLLY), and Sect. 6 gives the results of preliminary numerical tests with MOLLY. 

The idea is illustrated by Fig. 1. These equations constitute a readily understandable and concise representation of the widely used Verlet-I/r-RESPA impulse MTS method. The method was described first in [8, 9] but tested... 

In the r-RESPA method this operator is decomposed into two or more parts, for example... 

The underlying theory of r-RESPA is somewhat involved, but the final result and const quent implementation is actually rather straightforward, being very closely related to th velocity Verlet integration scheme. For our four-way decomposition the algorithm woul... 

---