Qualitative study of the symplectic Störmer–Verlet integrator

Hardy, David J.; Okunbor, Daniel I.
June 1995
Journal of Chemical Physics;6/8/1995, Vol. 102 Issue 22, p8978
Academic Journal
Symplectic numerical integrators, such as the Störmer–Verlet method, are useful in preserving properties that are not preserved by conventional numerical integrators. This paper analyzes the Störmer–Verlet method as applied to the simple harmonic model, whose generalization is an important model for molecular dynamics simulations. Restricting our attention to the one -dimensional case, both the exact solution and the Störmer–Verlet solution to this model are expressed as functions of the number of time steps taken, and then both of these functions are interpreted geometrically. The paper shows the existence of an upper bound on the error from the Störmer–Verlet method, and then an example is worked to demonstrate the closeness of this bound. © 1995 American Institute of Physics.


Related Articles

  • Vibrational energy relaxation of HOD in liquid D2O. Rey, Rossend; Hynes, James T. // Journal of Chemical Physics;2/8/1996, Vol. 104 Issue 6, p2356 

    Molecular Dynamics simulation is used to study the vibrational relaxation of the first excited state of the O–H stretch for HOD dissolved in D2O. The technique applied is based on a Landau–Teller type formula, in which the solvent contribution is computed classically, while the...

  • Large ionic clusters in concentrated aqueous NaCl solution. Degreve, Leo; da Silva, Fernando Luis B. // Journal of Chemical Physics;9/15/1999, Vol. 111 Issue 11, p5150 

    Investigates the stability of the local structure of aqueous sodium chloride solution by molecular dynamics simulations. Mean and maximum life-time of the ion pairs; Stability of non-neutral ion pairs; Finding that non-neutral ion pairs are stabilized by at least one counter-ion forming large...

  • Molecular-dynamics computer simulation of an aqueous NaCl solution: Structure. Zhu, Sheng-Bai; Robinson, G. Wilse // Journal of Chemical Physics;9/15/1992, Vol. 97 Issue 6, p4336 

    A classical molecular-dynamics simulation of a 1.791 molal aqueous NaCl solution is performed using a flexible/polarizable five-site water model. Through an investigation of the ion-water pair-correlation functions and the relative orientation of the waters of hydration, we are able to study the...

  • Highly concentrated salt solutions: Molecular dynamics simulations of structure and transport. Payne, V. A.; Forsyth, M.; Ratner, M. A.; Shriver, D. F.; de Leeuw, S. W. // Journal of Chemical Physics;4/1/1994, Vol. 100 Issue 7, p5201 

    Molecular dynamics (MD) simulations in NaI solutions, where the solvent has been represented by the Stockmayer fluid, were performed as a function of temperature, salt concentration, and solvent dipole strength. At higher temperatures contact ion pairs become more prevalent, regardless of...

  • Solute molecular dynamics in a mesoscale solvent. Malevanets, Anatoly; Kapral, Raymond // Journal of Chemical Physics;4/22/2000, Vol. 112 Issue 16 

    A hybrid molecular dynamics (MD) algorithm which combines a full MD description of solute-solute and solute-solvent interactions with a mesoscale treatment of solvent-solvent interactions is developed. The solvent dynamics is modeled on a mesoscale level by coarse graining the system into cells...

  • A density functional/molecular dynamics study of the structure of liquid nitromethane. Seminario, Jorge M.; Concha, Monica C.; Politzer, Peter // Journal of Chemical Physics;5/22/1995, Vol. 102 Issue 20, p8281 

    A molecular dynamics simulation of liquid nitromethane has been carried out, in which a density functional (DF) procedure was used to compute the properties needed to establish the intermolecular force field. The latter was subsequently updated by further DF calculations. The first and second...

  • A molecular dynamics study of an aqueous LiI solution between Lennard-Jones walls. Spohr, E.; Heinzinger, K. // Journal of Chemical Physics;2/15/1986, Vol. 84 Issue 4, p2304 

    An MD simulation of a 2.2 m aqueous LiI solution between 12-6 Lennard-Jones walls has been performed which extended over 19 ps at an average temperature of 283 K. The basic tetragonal cell contained 200 ST2 water molecules and eight ions of each kind described as Lennard-Jones spheres with an...

  • Influence of the solvent size on the behavior of model-polymers in solution. Komorowski, Joachim; Bruns, Wolfgang // Journal of Chemical Physics;10/1/1995, Vol. 103 Issue 13, p5756 

    Two-dimensional chains with C=4 up to C=20 segments are investigated by hard disks molecular dynamics simulations. Solvent particles are explicitly taken into account. In addition, the ratio S/σC of the diameters of the solvent molecules S and the chain segments C is varied between...

  • Theory of activated rate processes for arbitrary frequency dependent friction: Solution of the turnover problem. Pollak, Eli; Grabert, Hermann; Hänggi, Peter // Journal of Chemical Physics;10/1/1989, Vol. 91 Issue 7, p4073 

    An analytical theory is formulated for the thermal (classical mechanical) rate of escape from a metastable state coupled to a dissipative thermal environment. The working expressions are given solely in terms of the quantities entering the generalized Langevin equation for the particle dynamics....


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics