MOLECULAR-DYNAMICS IMPLEMENTATION OF THE GIBBS ENSEMBLE CALCULATION

A molecular dynamics version of the Gibbs ensemble calculation is proposed. This calculation is based on an extended Hamiltonian formalism that treats the temperature, volume, and the coupling of a single particle to the rest of the system as continuous dynamical degrees of freedom with their own equations of motion. Calculations on the truncated and shifted Lennard-Jones fluid are presented and compared to the Gibbs ensemble Monte Carlo results of Smit. Quantitative agreement is found between the molecular dynamics and Monte Carlo calculations. © 1994 American Institute of Physics.