ترغب بنشر مسار تعليمي؟ اضغط هنا

Path reweighting is a principally exact method to estimate dynamic properties from biased simulations - provided that the path probability ratio matches the stochastic integrator used in the simulation. Previously reported path probability ratios mat ch the Euler-Maruyama scheme for overdamped Langevin dynamics. Since MD simulations use Langevin dynamics rather than overdamped Langevin dynamics, this severely impedes the application of path reweighting methods. Here, we derive the path probability ratio $M_L$ for Langevin dynamics propagated by a variant of the Langevin Leapfrog integrator. This new path probability ratio allows for exact reweighting of Langevin dynamics propagated by this integrator. We also show that a previously derived approximate path probability ratio $M_{mathrm{approx}}$ differs from the exact $M_L$ only by $mathcal{O}(xi^4Delta t^4)$, and thus yields highly accurate dynamic reweighting results. ($Delta t$ is the integration time step, $xi$ is the collision rate.) The results are tested and the efficiency of path-reweighting is explored using butane as an example.
Markov State Models (MSM) are widely used to elucidate dynamic properties of molecular systems from unbiased Molecular Dynamics (MD). However, the implementation of reweighting schemes for MSMs to analyze biased simulations, for example produced by e nhanced sampling techniques, is still at an early stage of development. Several dynamical reweighing approaches have been proposed, which can be classified as approaches based on (i) Kramers rate theory, (ii) rescaling of the probability density flux, (iii) reweighting by formulating a likelihood function, (iv) path reweighting. We present the state-of-the-art and discuss the methodological differences of these methods, their limitations and recent applications.
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا