Abstract
The concise representation of complex high dimensional stochastic systems via a few reduced coordinates is an important problem in computational physics, chemistry, and biology. In this paper we use the first few eigenfunctions of the backward Fokker.Planck diffusion operator as a coarse-grained low dimensional representation for the long-term evolution of astochastic system and show that they are optimal under a certain mean squared error criterion. We denote the mapping from physical space to these eigenfunctions as the diffusion map. While in high dimensional systems these eigenfunctions are difficult to compute numerically by conventional methods such as finite differences or finite elements, we describe a simple computational data-driven method to approximate them from a large set of simulated data. Our method is based on defining an appropriatelyweighted graph on the set of simulated data and computing the first few eigenvectors and eigenvalues of the corresponding random walk matrix on this graph. Thus, our algorithm incorporates the local geometry and densityat each point into a global picture that merges data from different simulationruns in a natural way. Furthermore, we describe lifting and restriction operators between the diffusion map space and the original space. These operators facilitate the description of the coarse-grained dynamics, possibly in the form of a low dimensional effective free energy surface parameterized by the diffusion map reduction coordinates. They also enable a systematic exploration of such effective free energy surfaces through the design of additional intelligently biased computational experiments. Weconclude by demonstrating our method in a few examples.
Original language | English (US) |
---|---|
Pages (from-to) | 842-864 |
Number of pages | 23 |
Journal | Multiscale Modeling and Simulation |
Volume | 7 |
Issue number | 2 |
DOIs | |
State | Published - 2008 |
All Science Journal Classification (ASJC) codes
- General Chemistry
- Modeling and Simulation
- Ecological Modeling
- General Physics and Astronomy
- Computer Science Applications
Keywords
- Diffusion maps
- Dimensional reduction
- Fokker Planck operator
- Metastable states
- Normalized graph Laplacian
- Stochastic dynamical systems