@misc{QuerDonatiKelleretal., author = {Quer, Jannes and Donati, Luca and Keller, Bettina and Weber, Marcus}, title = {An automatic adaptive importance sampling algorithm for molecular dynamics in reaction coordinates}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-62075}, abstract = {In this article we propose an adaptive importance sampling scheme for dynamical quantities of high dimensional complex systems which are metastable. The main idea of this article is to combine a method coming from Molecular Dynamics Simulation, Metadynamics, with a theorem from stochastic analysis, Girsanov's theorem. The proposed algorithm has two advantages compared to a standard estimator of dynamic quantities: firstly, it is possible to produce estimators with a lower variance and, secondly, we can speed up the sampling. One of the main problems for building importance sampling schemes for metastable systems is to find the metastable region in order to manipulate the potential accordingly. Our method circumvents this problem by using an assimilated version of the Metadynamics algorithm and thus creates a non-equilibrium dynamics which is used to sample the equilibrium quantities.}, language = {en} } @article{DonatiWeberKeller, author = {Donati, Luca and Weber, Marcus and Keller, Bettina G.}, title = {A review of Girsanov Reweighting and of Square Root Approximation for building molecular Markov State Models}, series = {Journal of Mathematical Physics}, volume = {63}, journal = {Journal of Mathematical Physics}, number = {12}, publisher = {AIP Publishing}, doi = {10.1063/5.0127227}, pages = {123306-1 -- 123306-21}, abstract = {Dynamical reweighting methods permit to estimate kinetic observables of a stochastic process governed by a target potential U(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.}, language = {en} } @article{DonatiHeidaKelleretal., author = {Donati, Luca and Heida, Martin and Keller, Bettina G. and Weber, Marcus}, title = {Estimation of the infinitesimal generator by square-root approximation}, series = {J. Phys.: Condens. Matter}, volume = {30}, journal = {J. Phys.: Condens. Matter}, number = {42}, doi = {10.1088/1361-648X/aadfc8}, pages = {425201 -- 425201}, abstract = {In recent years, for the analysis of molecular processes, the estimation of time-scales and transition rates has become fundamental. Estimating the transition rates between molecular conformations is—from a mathematical point of view—an invariant subspace projection problem. We present a method to project the infinitesimal generator acting on function space to a low-dimensional rate matrix. This projection can be performed in two steps. First, we discretize the conformational space in a Voronoi tessellation, then the transition rates between adjacent cells is approximated by the geometric average of the Boltzmann weights of the Voronoi cells. This method demonstrates that there is a direct relation between the potential energy surface of molecular structures and the transition rates of conformational changes. We will show also that this approximation is correct and converges to the generator of the Smoluchowski equation in the limit of infinitely small Voronoi cells. We present results for a two dimensional diffusion process and alanine dipeptide as a high-dimensional system.}, language = {en} } @article{DonatiWeberKeller, author = {Donati, Luca and Weber, Marcus and Keller, Bettina G.}, title = {Markov models from the square root approximation of the Fokker-Planck equation: calculating the grid-dependent flux}, series = {Journal of Physics: Condensed Matter}, volume = {33}, journal = {Journal of Physics: Condensed Matter}, number = {11}, doi = {10.1088/1361-648X/abd5f7}, pages = {115902}, abstract = {Molecular dynamics (MD) are extremely complex, yet understanding the slow components of their dynamics is essential to understanding their macroscopic properties. To achieve this, one models the MD as a stochastic process and analyses the dominant eigenfunctions of the associated Fokker-Planck operator, or of closely related transfer operators. So far, the calculation of the discretized operators requires extensive MD simulations. The square-root approximation of the Fokker-Planck equation is a method to calculate transition rates as a ratio of the Boltzmann densities of neighboring grid cells times a flux, and can in principle be calculated without a simulation. In a previous work we still used MD simulations to determine the flux. Here, we propose several methods to calculate the exact or approximate flux for various grid types, and thus estimate the rate matrix without a simulation. Using model potentials we test computational efficiency of the methods, and the accuracy with which they reproduce the dominant eigenfunctions and eigenvalues. For these model potentials, rate matrices with up to \$\mathcal{O}\left(1{0}^{6}\right)\$ states can be obtained within seconds on a single high-performance compute server if regular grids are used.}, language = {en} } @article{QuerDonatiKelleretal., author = {Quer, Jannes and Donati, Luca and Keller, Bettina and Weber, Marcus}, title = {An automatic adaptive importance sampling algorithm for molecular dynamics in reaction coordinates}, series = {SIAM Journal on Scientific Computing}, volume = {40}, journal = {SIAM Journal on Scientific Computing}, number = {2}, doi = {10.1137/17m1124772}, pages = {A653 -- A670}, abstract = {In this article we propose an adaptive importance sampling scheme for dynamical quantities of high dimensional complex systems which are metastable. The main idea of this article is to combine a method coming from Molecular Dynamics Simulation, Metadynamics, with a theorem from stochastic analysis, Girsanov's theorem. The proposed algorithm has two advantages compared to a standard estimator of dynamic quantities: firstly, it is possible to produce estimators with a lower variance and, secondly, we can speed up the sampling. One of the main problems for building importance sampling schemes for metastable systems is to find the metastable region in order to manipulate the potential accordingly. Our method circumvents this problem by using an assimilated version of the Metadynamics algorithm and thus creates a non-equilibrium dynamics which is used to sample the equilibrium quantities.}, language = {en} }