Isobaric molecular dynamics simulations of hard sphere systems.
We describe an implementation of the Andersen algorithm for simulating the molecular dynamics in the isobaric isoenthalpic (NPH) ensemble for the hard sphere potential. The work is based on the adaptation of the Andersen algorithm to hard spheres by de Smedt et al. For a hard sphere system in the NPH ensemble, the particle velocities are not constant between collisions and we describe an efficient method for handling this part of the dynamics. The method is extended to give an NPT ensemble simulation of hard sphere systems by applying an ad hoc rescaling of the velocities. The accuracy of the algorithms is tested by comparison with traditional NVE simulation results for the structural, thermodynamic, and transport properties.
Gruhn
T
T
Department of Chemical Engineering, University of Massachusetts, Amherst, Massachusetts 01003, USA.
Monson
P A
PA