Supervised Learning and the Finite-Temperature String Method for Computing Committor Functions and Reaction Rates

TL;DR

This paper enhances algorithms for computing committor functions and reaction rates in rare event simulations by integrating supervised learning and the finite-temperature string method, improving accuracy and efficiency.

Contribution

It introduces modifications combining supervised learning and the finite-temperature string method into existing algorithms for better accuracy in rare event analysis.

Findings

Supervised learning improves committor prediction accuracy.

Finite-temperature string method enables homogeneous sampling.

Accurate reaction rates obtained with fewer samples.

Abstract

A central object in the computational studies of rare events is the committor function. Though costly to compute, the committor function encodes complete mechanistic information of the processes involving rare events, including reaction rates and transition-state ensembles. Under the framework of transition path theory (TPT), recent work [1] proposes an algorithm where a feedback loop couples a neural network that models the committor function with importance sampling, mainly umbrella sampling, which collects data needed for adaptive training. In this work, we show additional modifications are needed to improve the accuracy of the algorithm. The first modification adds elements of supervised learning, which allows the neural network to improve its prediction by fitting to sample-mean estimates of committor values obtained from short molecular dynamics trajectories. The second modification replaces the committor-based umbrella sampling with the finite-temperature string (FTS) method, which enables homogeneous sampling in regions where transition pathways are located. We test our modifications on low-dimensional systems with non-convex potential energy where reference solutions can be found via analytical or the finite element methods, and show how combining supervised learning and the FTS method yields accurate computation of committor functions and reaction rates. We also provide an error analysis for algorithms that use the FTS method, using which reaction rates can be accurately estimated during training with a small number of samples. The methods are then applied to a molecular system in which no reference solution is known, where accurate computations of committor functions and reaction rates can still be obtained.