On the Application of Widom's Test Particle Method to Homogeneous and Inhomogeneous Fluids

Abstract

Chemical potentials of a homogeneous and an inhomogeneous Lennard-Jones fluid have been determined by molecular dynamics simulations on the vector computer CYBER 205 by applying essentially the fictitious test particle method of Widom. For the homogeneous fluid we find, contrary to the previous result of Guillot and Guissani, that the simulated chemical potential is independent of the particle number. The crucial point, however, is a sufficiently large cut-off radius in the evaluation of the Boltzmann factor. Comparing with our WCA-type perturbation theory, we get agreement in the chemical potentials within 0.1 kT up to the density n[sgrave]³ = 0.80 and a difference of 0.2 kT at n[sgrave]³ = 0.85. For the inhomogeneous case we consider a fluid in a cylindrical pore and integrate Widom's equation over a certain probe volume as suggested earlier by us. Chemical potentials are then calculated independently in five different probe volumes, which are cylindrical shells. The results agree well from the second to the fourth shell. Inaccuracies in the innermost cylinder can be easily explained by bad statistics. In the shell close to the wall the extremely high local density is responsible for the inaccuracies. Extending the probe volume over all cylindrical shells besides the one closest to the wall is thought to yield rather reliable results for the chemical potential. As a by-product of the simulations we also obtained diffusion coefficients, which are given in an appendix.