This article is concerned with the averaging principle and its extensions
for stochastic dynamical systems with fast and slow degrees of
freedom. It is demonstrated how the \conventional" averaging principle
results from asymptotic multiscale analysis, how one can construct
an indicator for its (in-)appropriateness, and how, if inappropriate, it
may be extended into an improved approximation. The conventional
scheme contains averages over the entire accessible state space of the
fast degrees of freedom and may thus fail if these fast degrees of freedom
exhibit long-term (auto-)correlations. In contrast, the improved
scheme combines several conditional averages with a Markov jump process
that is designed to represent the
ipping process between the
conditional averages and thus incorporates the important long-term
correlations. All important steps of the derivation are illustrated by
numerical experiments. Application to problems from molecular dynamics
is discussed.
We consider a particle constrained to a submanifold ? of the
configuration space Rm. Using that the notion of holonomic constraints coincides
with integrability of the corresponding vector field, we show how this property
naturally determines local coordinates on ?. We give a rigorous justification for
the calculation of the mean force along a constrained coordinate, and we provide
a concise geometrical interpretation of the different contributions to the mean
force in terms of the unconstrained vector field and extrinsic curvature properties
of ? in Rm. Our approach gives rise to a Hybrid Monte-Carlo based algorithm
that can be used to compute the mean force acting on selected coordinates in the
context of thermodynamic free energy statistics.
We consider the problem of automatically extracting simplified models out of complex high--dimensional and t
ime--dependent data. The simplified model is given by a linear Langevin equation with time--varying coeffici
ents. The reduced model may still be high--dimensional, but it is physically intuitive and much easier to in
terpret than the original data. In particular we can distinguish whether dynamical effects are influenced b
y friction, noise, or deterministic motion. The parameters for the reduced model are obtained by a robust an
d efficient numerical predictor--corrector scheme which relies on analytical solutions to a maximum-likeliho
od problem provided the time steps between successive observations are not too large. If the data set is ver
y heterogeneous the time series is better described not by a single model, but by a collection of reduced mo
dels. This scenario is accounted for by embedding the parameter estimation procedure into the framework of h
idden Markov models, i.e., we decompose the data into several subsets, each of which gives rise to an approp
riate linear Langevin model. The switching between the local model is done by a Markov jump process. The opt
imal decomposition into submodels can then be regarded as one global Langevin model with piecewise constant
coefficients. We illustrate the performance of the algorithm by means of several examples. Especially we foc
us on the numerical error as a function of the time step of the observation sequence.
We consider the problem of automatically extracting simplified models out of complex high-dimensional and time-dependent data. The simplified model is given by a linear Langevin equation with time-varying coefficients. The reduced model may still be high-dimensional, but it is physically intuitive and much easier to interpret than the original data. In particular we can distinguish whether certain dynamical effects are influenced by friction, noise, or systematic drift.
The parameters for the reduced model are obtained by a robust and efficient numerical predictor-corrector scheme which relies on analytical solutions to a maximum-likelihood problem provided the time steps between successive observations are not too large. Our approach emphasizes the specific hypoelliptic structure of the Langevin equation given high-dimensional observation data, and therefore can be considered as complemetary to the procedure recently proposed in \emph{Horenko et al. (submitted SIAM MMS, 2007)} by one of the authors, or to the problem of incomplete (one-dimensional) observations \emph{Pokern et al. (submitted to JRSSB, 2007)}. If the data set is very heterogeneous the time series is better described not by a single model, but by a collection of reduced models. This scenario is accounted for by embedding the parameter estimation procedure into the framework of hidden Markov models which it is particularly suited to treat high-dimensional data. That is, we decompose the data into several subsets, each of which gives rise to an appropriate linear Langevin model, where the switching between the local model is done by a Markov jump process. The optimal decomposition into submodels can then be regarded as one global Langevin model with piecewise constant coefficients. We illustrate the performance of the algorithm by means of several examples. Especially we focus on the numerical error as a function of the time step of the observation sequence.
We comment on two different notions of the thermodynamical free energy
that are used in Hamiltonian molecular dynamics. Both concepts have
different scopes of applications as was pointed out recently in the
context of high--friction Langevin dynamics. We show that problems
that rely on either definition can be treated in a uniform way using
constrained molecular dynamics. Not only proves this useful in
designing algorithms that sample the free energy landscape, but it
also clarifies the relation between seemingly contradictory results
that are present in the literature.
We present a formal procedure for structure-preserving model reduction of linear second-order control problems that appear in a variety of physical contexts, e.g., vibromechanical systems or electrical circuit design. Typical balanced truncation methods that project onto the subspace of the largest Hankel singular values fail to preserve the problem's physical structure and may suffer from lack of stability. In this paper we adopt the framework of generalized Hamiltonian systems that covers the class of relevant problems and that allows for a generalization of balanced truncation to second-order problems.
It turns out that the Hamiltonian structure, stability and passivity are preserved if the truncation is done by imposing a holonomic constraint on the system rather than standard Galerkin projection.
We derive a Crooks-Jarzynski-type identity for computing free energy differences between metastable states that is based on nonequilibrium diffusion processes. Furthermore we outline a brief derivation of an infinite-dimensional stochastic partial differential equation that can be used to efficiently generate the ensemble of trajectories connecting the metastable states.
We study balanced model reduction for stable bilinear systems in the limit of partly vanishing Hankel singular values. We show that the dynamics admit a splitting into fast and slow subspaces and prove an averaging principle for the slow dynamics. We illustrate our method with an example from stochastic control (density evolution of a dragged Brownian particle) and discuss issues of structure preservation and positivity.
We study balanced model reduction of partially-observed linear stochastic differential equa- tions of Langevin type. Balancing the equations of motion gives rise to a singularly perturbed system of equations with slow and fast degrees of freedom, and we prove that in the limit of the fast variables becoming infinitely fast, the solutions converge to the solution of a reduced-order Langevin equation. We illustrate the method with several numerical examples and discuss the relation to model reduction of deterministic control systems that have an underlying Hamiltonian structure.
We propose a nonequilibrium sampling method for computing free energy profiles along a given reaction coordinate. The method consists of two parts: a controlled Langevin sampler that generates nonequilibrium bridge paths conditioned by the reaction coordinate, and Jarzynski’s formula for reweighting the paths. Our derivation of the equa- tions of motion of the sampler is based on stochastic perturbation of a controlled dissipative Hamiltonian system, for which we prove Jarzynski’s identity as a special case of the Feynman-Kac formula. We illustrate our method by means of a suitable numerical example and briefly discuss issues of optimally choosing the control protocol for the reaction coordinate.