ABSTRACT:The split integration symplectic method (SISM) for molecular dynamics (MD) integration using normal mode analysis based on a factorization of the Liouville propagator is presented. This approach is quite distinct from others that use fractional-step methods, owing to the analytical treatment of high-frequency motions. The method involves splitting the total Hamiltonian of the system into a harmonic part and the remaining part. Then the Hamilton equations are solved using a second-order generalized leapfrog integration scheme in which the purely harmonic Hamiltonian (which represents the main contribution of the chemical bonds and angles) is treated analytically, i.e., independent of the step size of integration, by a normal mode analysis that is carried out only once, at the beginning of calculation. The whole integration step combines analytical evolution of the harmonic part of the Hamiltonian with a correction arising from the remaining part. The proposed algorithm requires only one force evaluation per integration step. The algorithm was tested on a simple system of linear chain molecules. Results demonstrate the method makes possible the integration of the MD equations over larger time steps without loss of stability while being economical in computer time.