An exact multiple-time-step variational formulation for the committor and the transition rate
This addresses a methodological bottleneck for researchers studying rare events in molecular dynamics or stochastic systems, though it appears incremental as it builds on existing variational formulations.
The paper tackled the bias in estimating the committor and transition rate from trajectory data due to lag time sensitivity, introducing an alternative expression that yields exact committor estimates at any lag time, with numerical tests showing significantly reduced sensitivity to lag time choice.
For a transition between two stable states, the committor is the probability that the dynamics leads to one stable state before the other. It can be estimated from trajectory data by minimizing an expression for the transition rate that depends on a lag time. We show that an existing such expression is minimized by the exact committor only when the lag time is a single time step, resulting in a biased estimate in practical applications. We introduce an alternative expression that is minimized by the exact committor at any lag time. Numerical tests on benchmark systems demonstrate that our committor and resulting transition rate estimates are much less sensitive to the choice of lag time. We derive an additional expression for the transition rate, relate the transition rate expression to a variational approach for kinetic statistics based on the mean-squared residual, and discuss further numerical considerations with the aid of a decomposition of the error into dynamic modes.