Logo image
Unbiased estimation of equilibrium, rates, and committors from Markov state model analysis
Preprint   Open access

Unbiased estimation of equilibrium, rates, and committors from Markov state model analysis

John D Russo, Jeremy Copperman, David Aristoff, Gideon Simpson and Daniel M Zuckerman
27 May 2021
url
https://doi.org/10.48550/arxiv.2105.13402View
Preprint (Author's original)arXiv.org - Non-exclusive license to distribute Open

Abstract

Markov state models (MSMs) have been broadly adopted for analyzing molecular dynamics trajectories, but the approximate nature of the models that results from coarse-graining into discrete states is a long-known limitation. We show theoretically that, despite the coarse graining, in principle MSM-like analysis can yield unbiased estimation of key observables. We describe unbiased estimators for equilibrium state populations, for the mean first-passage time (MFPT) of an arbitrary process, and for state committors - i.e., splitting probabilities. Generically, the estimators are only asymptotically unbiased but we describe how extension of a recently proposed reweighting scheme can accelerate relaxation to unbiased values. Exactly accounting for 'sliding window' averaging over finite-length trajectories is a key, novel element of our analysis. In general, our analysis indicates that coarse-grained MSMs are asymptotically unbiased for steady-state properties only when appropriate boundary conditions (e.g., source-sink for MFPT estimation) are applied directly to trajectories, prior to calculation of the appropriate transition matrix.

Metrics

10 Record Views

Details

Logo image