Abstract
A new version of the test particle method for determining the chemical potential by Monte Carlo simulations is proposed. The method, applicable to any fluid at any density, combines the Widom's test particle insertion method with the ideas of the scaled particle theory, gradual insertion method and multistage sampling. Its applicability is exemplified by evaluating the chemical potential of the hard sphere fluid at a very high density in semi-grand-canonical and grand-canonical ensembles. A theory estimating the efficiency (i.e. statistical errors) of the method is proposed and the results are compared with the Widom's and gradual insertion methods, and the analytic results.