Relative Entropy Methods for the Approximation of Reactive Trajectories
This work addresses challenges in molecular simulation for chemical reaction modeling, but it is incremental as it builds on existing methods for approximating reactive trajectories.
The paper tackles the problem of approximating reactive trajectories in molecular simulation by replacing the exact singular drift with an approximate committor function, resulting in a derived formula for relative entropy and a stochastic gradient descent method for training the approximation.
Motivated by challenges arising in molecular simulation, we study reactive trajectories of the overdamped Langevin dynamics, i.e. trajectories observed as they pass from a set A corresponding to the reagents of a chemical reaction to a set B corresponding to the products. Reactive trajectories are known to have the same distribution as trajectories of the overdamped Langevin dynamics biased by a singular drift related to the committor function. In this work, we assess the effect of replacing the exact singular drift with an approximation based on an approximate committor function. We derive a convenient formula for the relative entropy between the distributions of exact and approximate reactive trajectories, and we propose a stochastic gradient descent method for minimizing the entropy to train an approximate committor function on the fly while computing reactive trajectories. We also devise a model assessment procedure for comparing the qualities of different approximations to the committor function based on the relative entropy.