Select a project:

Project Summary

Molecular dynamics and related computational methods enable the description of bio­logical systems with all-atom detail. However, these approaches are limited regarding simulation times and system sizes. A systematic way to bridge the micro-macro scale range between molecular dynamics and experiments is to apply coarse-graining (CG) techniques. The basic idea of (CG) is to replace the high-dimensional all-atom descrip­tion of the system by a reduced representation that preserves a suf.cient accuracy of the properties of interest. Obviously, numerical evaluation of the coarse-grained system would require less resources so that an increase of orders of magnitude in the simu­lated time and length scales can be achieved in this way. Many different CG approaches have been introduced over the years. Most of them are validated by means of numerical experiments only, while reliable theoretical insight into their approximation properties is missing. The ambitious goal of this project is to put CG approximations on a solid mathematical footing.

To this end, we will start with the three main challenges every approach to CG is facing: (1) How to identify a suitable reduced state space, i.e., a lower-dimensional sub­space of the all-atom state space, onto which the full-atom dynamics can be projected without destroying its essential properties (by essential properties of the dynamics we mean their longest timescales and the transport and kinetic properties associated with them), (2) how to find a closed representation of the projected dynamics in terms of the result­ing coarse graining coordinates (CGC), and (3) how to provide an ef.cient numerical realization. Our approach to (1) will rely on novel CG error estimators exploiting ideas from dimension-adaptive sparse quadrature. CG is then understood as a projection of the transfer operator, re.ecting the full-atom dynamics, onto the reduced state space. Concerning (2) we will investigate existing, prevalent methods, including very recent approaches to the extraction of CG dynamics for deterministic, discrete dynamical sys­tems, with respect to their applicability in the context of projected transfer operators. Basic multigrid ideas of scale separation and localization will be used in (3) to derive and analyze a multiscale discretization of the transfer operator together with a multi­level strategy integrating the adaptive concepts from (1) and (2). In order to validate our approach and compare it with existing, heuristic CG strategies, we will .rst apply this multilevel strategy to systems with fast-slow scale separation, where analytical results are available. We will also investigate whether existing strategies can be justi.ed by our abstract results. The theoretical considerations will be complemented by applications to small but realistic molecular systems. The long-term goal of the project is a multilevel coarse grained description of supramolecular aggregation processes related to neurode­generative diseases, with simulations complemented and validated by experiments.

Project publications

Dellnitz, M. and Klus, S. and Ziessler, A. (2017) A Set-Oriented Numerical Approach for Dynamical Systems with Parameter Uncertainty. SIAM Journal on Applied Dynamical Systems, 16 (1). pp. 120-138. ISSN 1536-0040

Klus, S. and Nüske, F. and Koltai, Péter and Wu, H. and Kevrekidis, Ioannis and Schütte, Ch. and Noé, F. (2017) Data-driven model reduction and transfer operator approximation. J. Nonlin. Sci. . (Submitted)

Klus, S. and Schütte, Ch. (2016) Towards tensor-based methods for the numerical approximation of the Perron-Frobenius and Koopman operator. Journal of Computational Dynamics . ISSN 2158-2491

Klus, S. and Koltai, Péter and Schütte, Ch. (2016) On the numerical approximation of the Perron-Frobenius and Koopman operator. Journal of Computational Dynamics, 3 (1). pp. 51-79. ISSN 2158-2491

Klus, S. and Gelß, P. and Peitz, S. and Schütte, Ch. (2016) Tensor-based dynamic mode decomposition. SIAM Journal on Scientific Computing . ISSN ISSN 1064-8275 (print); 1095-7197 (electronic) (Submitted)

Kornhuber, R. and Yserentant, H. (2016) Numerical Homogenization of Elliptic Multiscale Problems by Subspace Decomposition. Multiscale Model. Simul., 14 (3). pp. 1017-1036. ISSN print: 1540-3459; online: 1540-3467

Enciso, M. and Schütte, Ch. and Delle Site, L. (2015) Influence of pH and sequence in peptide aggregation via molecular simulation. Journal of Chemical Physics, 143 (24). p. 243130. ISSN 0021-9606

Schütte, Ch. and Sarich, M. (2015) A Critical Appraisal of Markov State Models. SFB 1114 Preprint, 15-18 . ISSN 1438-0064

Project Details

  • Status:  In Progress
  • Head(s):  Prof. Schütte, Prof. Kornhuber, Prof. Koksch
  • Project member(s):  Dr. Stefan Klus Dr. Johann Moschner, Dr. Marco Sarich, Maren-Wanda Wolf, Niklas Wulkow, Han Cheng Lie, Patrick Gelß, Dr. Andreas Bittracher
  • Participating Institution(s):  FU Berlin
  • Area:  B: Uniform meso scale behavior in scaling cascades
  • Positions available?: 
Back to top