Optimal Rejection-Free Path Sampling

Abstract: We propose an efficient novel path sampling-based framework designed to accelerate the investigation of rare events in complex molecular systems. A key innovation is the shift from sampling restricted path ensemble distributions, as in transition path sampling, to directly sampling the distribution of shooting points. This allows for a rejection-free algorithm that samples the entire path ensemble efficiently. Optimal sampling is achieved by applying a selection bias that is the inverse of the free energy along a reaction coordinate. The optimal reaction coordinate, the committor, is iteratively constructed as a neural network using AI for Molecular Mechanism Discovery (AIMMD), concurrently with the free energy profile, which is obtained through reweighting the sampled path ensembles. We showcase our algorithm on theoretical and molecular bechnmarks, and demonstrate how it provides at the same time molecular mechanism, free energy, and rates at a moderate computational cost.