Diffusion maps, reduction coordinates, and low dimensional representation of stochastic systems

R. R. Coifman, S. Lafon, I. G. Kevrekidis, M. Maggioni, B. Nadler

Research output: Contribution to journalArticlepeer-review

219 Scopus citations

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 languageEnglish (US)
Pages (from-to)842-864
Number of pages23
JournalMultiscale Modeling and Simulation
Volume7
Issue number2
DOIs
StatePublished - 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

Fingerprint

Dive into the research topics of 'Diffusion maps, reduction coordinates, and low dimensional representation of stochastic systems'. Together they form a unique fingerprint.

Cite this