Portal de Periódicos da CAPES

An Efficient Way to Conserve the Total Energy in Molecular Dynamics Simulations; Boundary Effects on Energy Conservation and Dynamic Properties

1988; Taylor & Francis; Volume: 1; Issue: 5 Linguagem: Inglês

10.1080/08927028808080956

1029-0435

Olle Teleman,

Spectroscopy and Quantum Chemical Studies

Abstract Molecular Dynamics simulations of pure water are reported, where the so-called jump correction i.e., inclusion of the Dirac delta force corresponding to the potential discontinuity at cutoff, ensures that the force is the negative gradient of the potential everywhere. Results are compared to those obtained for potentials truncated in the ordinary way and for potentials where the discontinuity is removed by means of a switching function. The jump correction and the switching functions provide excellent conservation of the total energy, and remove the need for temperature scaling at equilibrium. Radial distributions are insensitive to boundary conditions. Translational and rotational diffusion vary considerably with boundary conditions and also with cutoff distance.