Efficient computation, sensitivity, and error analysis of committor probabilities for complex dynamical processes

In many fields of physics, chemistry, and biology, the characterization of rates and pathways between certain states or species is of fundamental interest. The central mathematical object in such situations is the committor probability—a generalized reaction coordinate that measures the progress of...

Full description

Saved in:
Bibliographic Details
Main Authors: Prinz, Jan-Hendrik (Author) , Held, Martin (Author) , Smith, Jeremy C. (Author) , Noé, Frank (Author)
Format: Article (Journal)
Language:English
Published: June 24, 2011
In: Multiscale modeling & simulation
Year: 2011, Volume: 9, Issue: 2, Pages: 545-567
ISSN:1540-3467
DOI:10.1137/100789191
Online Access:Verlag, lizenzpflichtig, Volltext: https://doi.org/10.1137/100789191
Verlag, lizenzpflichtig, Volltext: https://epubs.siam.org/doi/10.1137/100789191
Get full text
Author Notes:Jan-Hendrik Prinz, Martin Held, Jeremy C. Smith, and Frank Noé
Description
Summary:In many fields of physics, chemistry, and biology, the characterization of rates and pathways between certain states or species is of fundamental interest. The central mathematical object in such situations is the committor probability—a generalized reaction coordinate that measures the progress of the process as the probability of proceeding to the target state rather than relapsing to the source state. Here, we conduct a numerical analysis of the committor. First, it is shown that committors can be expressed by the stationary eigenfunctions of a modified dynamical operator, thus relating the committors to the dominant eigenfunctions of the original operator. Based on this reformulation, committors can be efficiently computed for systems with large state spaces. Moreover, a sensitivity analysis of the committor is conducted, which allows its statistical uncertainty from estimation to be quantified within a Bayesian framework. The methods are illustrated on two examples of diffusive dynamics: a two-dimensional model potential with three minima, and a three-dimensional model representing protein-ligand binding.
Item Description:Gesehen am 17.01.2023
Physical Description:Online Resource
ISSN:1540-3467
DOI:10.1137/100789191