Calculation of free energy through successive umbrella sampling.
We consider an implementation of umbrella sampling in which the pertinent range of states is subdivided into small windows that are sampled consecutively and linked together. This allows us to simulate without a weight function or to extrapolate the results to the neighboring window in order to estimate a weight function. Additionally, we present a detailed error analysis in which we demonstrate that the error in umbrella sampling is controlled and, in the absence of sampling difficulties, independent of the window sizes. In this case, the efficiency of our implementation is comparable to a multicanonical simulation with a very good weight function, which in our scheme does not need to be known ahead of time. The analysis also allows us to detect sampling difficulties such as correlations between adjacent windows and provides a test of equilibration. We exemplify the scheme by simulating the liquid-vapor coexistence in a Lennard-Jones system.
Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139-4307, USA. virnau@mit.edu
Peter
P
Virnau
Marcus
M
Muller