Approximating free energy and committor landscapes in standard transition path sampling using virtual interface exchange.

Abstract

Transition path sampling is a powerful technique for investigating rare transitions, especially when the mechanism is unknown and one does not have access to the reaction coordinate. Straightforward application of transition path sampling does not directly provide the free energy landscape nor the kinetics. This drawback has motivated the development of path sampling extensions able to simultaneously access both kinetics and thermodynamics, such as transition interface sampling, and the reweighted path ensemble. However, performing transition interface sampling is more involved than standard two-state transition path sampling and still requires (some) insight into the reaction to define interfaces. While packages that can efficiently compute path ensembles for transition interface sampling are now available, it would be useful to directly compute the free energy from a single standard transition path sampling simulation. To achieve this, we present here an approximate method, denoted virtual interface exchange transition path sampling, that makes use of the rejected pathways in a form of waste recycling. The method yields an approximate reweighted path ensemble that allows an immediate view of the free energy landscape from a standard single transition path sampling simulation, as well as enables a committor analysis.