Refine
Year of publication
- 2013 (7) (remove)
Document Type
- Article (4)
- ZIB-Report (3)
Language
- English (7)
Is part of the Bibliography
- no (7)
Keywords
- rare events (2)
- stochastic control (2)
- Girsanov transformation (1)
- Markov state models (1)
- change of measure (1)
- cumulant generating function (1)
- dynamic programming (1)
- linear response (1)
- long timescales (1)
- moleculare dynamics (1)
Institute
A good deal of molecular dynamics simulations aims at predicting and quantifying rare events, such as the folding of a protein or a phase transition. Simulating rare events is often prohibitive, especially if the equations of motion are high-dimensional, as is the case in molecular dynamics. Various algorithms have been proposed for efficiently computing mean first passage times, transition rates or reaction pathways. This article surveys and discusses recent developments in the field of rare event simulation and outlines a new approach that combines ideas from optimal control and statistical mechanics. The optimal control approach described in detail resembles the use of Jarzynski's equality for free energy calculations, but with an optimized protocol that speeds up the sampling, while (theoretically) giving variance-free estimators of the rare events statistics. We illustrate the new approach with two numerical examples and discuss its relation to existing methods.
Rare but important transition events between long lived states are a key feature of many molecular systems. In many cases the computation of rare event statistics by direct molecular dynamics (MD) simulations is infeasible even on the most powerful computers because of the immensely long simulation timescales needed. Recently a technique for spatial discretization of the molecular state space designed to help overcome such problems, so-called Markov State Models (MSMs), has attracted a lot of attention. We review the theoretical background and algorithmic realization of MSMs and illustrate their use by some numerical examples. Furthermore we introduce a novel approach to using MSMs for the efficient solution of optimal control problems that appear in applications where one desires to optimize molecular properties by means of external controls.
Linear response theory and optimal control for a molecular system under nonequilibrium conditions
(2013)
In this paper, we propose a straightforward generalization of linear
response theory to systems in nonequilibrium that are subject to
nonequilibrium driving. We briefly revisit the standard linear response
result for equilibrium systems, where we consider Langevin dynamics
as a special case, and then give an alternative derivation using a
change-of-measure argument that does not rely on any stationarity or
reversibility assumption. This procedure moreover easily enables us
to calculate the second order correction to the linear response formula
(which may or may not be useful in practice). Furthermore, we outline
how the novel nonequilibirum linear response formula can be used to
compute optimal controls of molecular systems for cases in which one
wants to steer the system to maximize a certain target expectation
value. We illustrate our approach with simple numerical examples.
Grand-canonical-like molecular-dynamics simulations by using an adaptive-resolution technique
(2013)