A review of Girsanov Reweighting and of Square Root Approximation for building molecular Markov State Models
Authors:
- Luca DonatiLDLuca Donati
- Marcus WeberMWMarcus Weber
- Bettina G. KellerBKBettina G. Keller
Preprint
Abstract
Dynamical reweighting methods permit to estimate kinetic observables of a stochastic process governed by a target potential $\tilde{V}(x)$ from trajectories that have been generated at a different potential $V(x)$. In this article, we present Girsanov reweighting and Square Root Approximation (SqRA): the first method reweights path probabilities exploiting the Girsanov theorem and can be applied to Markov State Models (MSMs) to reweight transition probabilities; the second method was originally developed to discretize the Fokker-Planck operator into a transition rate matrix, but here we implement it into a reweighting scheme for transition rates. We begin by reviewing the theoretical background of the methods, then present two applications relevant to Molecular Dynamics (MD), highlighting their strengths and weaknesses.
Keywords
All Related Versions
This publication has 35 references indexed in Scilit:
- Dynamical reweighting methods for Markov modelsCurrent Opinion in Structural Biology, 2020
- Protein Conformational Populations and Functionally Relevant SubstatesAccounts of Chemical Research, 2013
- OpenMM 4: A Reusable, Extensible, Hardware Independent Library for High Performance Molecular SimulationJournal of Chemical Theory and Computation, 2012
- An Analysis of the Validity of Markov State Models for Emulating the Dynamics of Classical Molecular Systems and EnsemblesJournal of Chemical Theory and Computation, 2011
- Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculationsComputer Physics Communications, 2001
- Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscapeChemical Physics Letters, 2000
- A Direct Approach to Conformational Dynamics Based on Hybrid Monte CarloJournal of Computational Physics, 1999
- Chaos, Fractals, and Noise1994
- Location of saddle points and minimum energy paths by a constrained simplex optimization procedureTheoretical Chemistry Accounts, 1979
- The Onsager-Machlup function as Lagrangian for the most probable path of a diffusion processCommunications in Mathematical Physics, 1978