Generalized Path Reweighting and History-Dependent Free Energies

Abstract

Transition interface sampling (TIS) and replica exchange TIS (RETIS) are powerful methods for computing rates of rare events inaccessible to straightforward molecular dynamics (MD) simulations. Path reweighting extends their output, enabling the evaluation of diverse thermodynamic and kinetic quantities, including reaction prediction metrics, activation barriers, committor functions, and free energies. The recently developed Infinity-RETIS algorithm boosts parallel efficiency through asynchronous replica exchanges in the infinite-swap limit, eliminating the wall-time bottlenecks of conventional RETIS. This approach introduces fractional samples and biased sampling distributions, requiring a generalized path reweighting framework, for which we derive expressions demonstrating how exact dynamic and thermodynamic variables can be computed. We then focus on a special class of free energy surfaces defined by history-dependent conditions, whose values are influenced by kinetic factors such as particle mass and friction, unlike standard unconditional free energy surfaces. Even with suboptimal reaction coordinates, these conditional free energies can reveal kinetically relevant barriers that may be misrepresented by standard unconditional free energies, thereby providing a rigorous and versatile tool for characterizing complex molecular transitions.

15 pages, 5 figures

References in corpus (11)

• A Novel Path Sampling Method for the Calculation of Rate Constants
• Sampling rare switching events in biochemical networks
• Rate constants for diffusive processes by partial path sampling
• Reaction rate calculation by parallel path swapping
• Free Energy Calculations of Membrane Permeation: Challenges due to Strong Headgroup-Solute Interactions
• Efficiency analysis of reaction rate calculation methods using analytical models I: The 2D sharp barrier
• Contour forward flux sampling: Sampling rare events along multiple collective variables
• Exchanging replicas with unequal cost, infinitely and permanently
• Enhanced path sampling using subtrajectory Monte Carlo moves
• Crystal Nucleation Kinetics and Mechanism: Influence of Interaction Potential
• Combining multiple interface set path ensembles with MBAR reweighting