Biased sampling of nonequilibrium trajectories: Can fast switching simulations outperform conventional free energy calculation methods?

Abstract

We have investigated the maximum computational efficiency of reversible work calculations that change control parameters in a finite amount of time. Because relevant nonequilibrium averages are slow to converge, a bias on the sampling of trajectories can be beneficial. Such a bias, however, can also be employed in conventional methods for computing reversible work, such as thermodynamic integration or umbrella sampling. We present numerical results for a simple one-dimensional model and for a Widom insertion in a soft sphere liquid, indicating that, with an appropriately chosen bias, conventional methods are in fact more efficient. We describe an analogy between nonequilibrium dynamics and mappings between equilibrium ensembles, which suggests that the practical inferiority of fast switching is quite general. Finally, we discuss the relevance of adiabatic invariants in slowly driven Hamiltonian systems for the application of Jarzynski's theorem. © 2005 American Chemical Society.