Molecular dynamics and related computational methods enable the description of biological 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 description 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 simulated 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 subspace 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 resulting 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 systems, 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 multilevel 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 neurodegenerative diseases, with simulations complemented and validated by experiments.
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