Measuring Coexisting Densities from a Two-Phase Molecular Dynamics Simulation by Voronoi Tessellations

A new algorithm is presented that allows for the determination of bulk liquid and vapor densities from a two-phase Molecular Dynamics (2φMD) simulation. This new method does not use any arbitrary cutoffs for phase definitions; rather it uses single-phase simulations as a self-consistency check. The method does not use any spatial bins for generating histograms of local properties, thereby avoiding the statistical issues associated with bins. Finally, it allows one to approach very close to the critical point. The new method utilizes Voronoi tessellations to determine the molecular volume of every point at every instance in a molecular dynamics simulation. Since the molecular volume is calculated throughout the simulation, statistical parameters such as the average molecular volume and average molecular variance are easy to obtain. To define the phases, the normalized variance of the molecular volume from 1φMD and 2φMD is used as a self-consistency check. The new method gives new insight into the nature of the near-subcritical fluid. The critical properties from this analysis are T_c = 1.293 and ρ_c = 0.313. Direct simulation of the two-phase system was performed up to a temperature of 1.292. The results show excellent agreement to experimental results and Gibbs Ensemble Monte Carlo for coexisting densities. We see that well below the critical temperature, some particles are neither liquid nor vapor. These interfacial particles are primarily, but not exclusively, concentrated at the bulk interface. However, as we approach the critical point, some particles are considered both liquid and vapor. These interfacial particles are distributed through the system.

Copyright (2007) American Institute of Physics. This article may be downloaded for personal use only. Any other use requires prior permission of the author and the American Institute of Physics.