A Parallel-pulling Protocol for Free-Energy Evaluation
arXiv:1105.4840 · doi:10.1103/PhysRevE.85.036702
Abstract
Jarzynski's equality (JE) allows us to compute free energy differences from distributions of work. In molecular dynamics simulations, the traditional way of constructing work distributions is to perform as many pulling simulations as possible. But reliable work distributions are not always produced in a finite number of simulations. The computational cost of using JE is not less than other commonly used methods such as Thermodynamic Integration and Umbrella Sampling methods. Here we first show a different proof of JE based on the idea of step-wise pulling procedures that is efficient in computing free energies by using JE. The key point in our proof is that the processes of turning-on/off a harmonic potential to perform work are described by double Heaviside functions of time. We then show that the distributions of work performed by the potential can be easily generated from the distributions of a reaction coordinate along a pathway. Based on the proof, we propose sequential and parallel step-wise pulling protocols for generating work distributions that require suitable relaxation time at each pulling step. The criterion for reliable work distributions is that there must be sufficient mutual overlaps between the adjacent distributions of the reaction coordinate along the pathway. We arrive at an alternative formula (besides JE) to compute free energy differences from the averaged values of the reaction coordinate. The combination of JE and the alternative formula provides a viable way to determine the accuracy of computed free energy differences. For the stretching of a deca-alanine molecule, our approach requires 21 parallel simulations and relaxation time as small as 0.4 ns for each simulation to estimate free energy differences with an uncertainty of about 13%.
13 pages, 6 figures, 3 appendices
References in corpus (9)
- Verification of the Crooks fluctuation theorem and recovery of RNA folding free energies
- An experimental test of the Jarzynski equality in a mechanical experiment
- Comparison of far-from-equilibrium work relations
- A "fast growth" method of computing free energy differences
- Optimized free energies from bidirectional single-molecule force spectroscopy
- Failure of the work-Hamiltonian connection for free energy calculations
- Note on the Jarzynski Equality
- Comment on: Failure of the Work-Hamiltonian Connection for Free-Energy Calculations [Phys Rev Lett 100, 020601 (2008), arXiv:0704.0761]
- Comment on "Failure of the work-Hamiltonian connection for free-energy calculations" by Jose M. G. Vilar and J. Miguel Rubi