Alternative ensemble averages in molecular dynamics simulation of hard spheres

Abstract

ABSTRACTWe consider application to the hard sphere (HS) model of the mapped-averaging framework for generating alternative ensemble averages for thermodynamic properties. Specifically, we develop and examine new formulas for the pressure, the singlet and pair densities, and the cavity-correlation function inside the HS core; the pressure formula in particular is constructed such that it gives an ensemble average that exactly corrects the second-order virial equation of state. The force plays a central role in mapped-averaging expressions, and we write them in a way that accounts for the impulsive, event-driven nature of the HS dynamics. Comparison between results obtained conventionally versus mapped averaging finds that the latter has some advantage at low density, while both perform equally well (in terms of uncertainties for a given amount of sampling) at higher densities.