Some computational aspects of Jarzynki's non-equilbrium work theorem
Prof. Christoph Dellago
Université de Vienne
Mon, May. 19th 2008, 14:15
Salle Claude Itzykson, Bât. 774, Orme des Merisiers
Jarzynski's recent theorem relates equilibrium free energies to the statistics of work accumulated in irreversible transformations between ensembles. Fast switching simulations based on this theorem can be carried out with arbitrary switching rate promising inexpensive free energy estimation algorithms. Straightforward application of the fast switching method, however, is efficient only if the system is driven away from equilibrium only mildly thus requiring the generation of long non-equilibrium trajectories. The reason for this limitation lies in the work distribution at fast switching rates: the values of work contributing most to the free energy are rarely sampled causing large statistical errors. Here, I will discuss several computational strategies, including large time step integration and bias path sampling, to circumvent this problem. The algorithms will be illustrated using the calculation of the chemical potential of simple fluid and the stretching of a small polypeptide as examples.