@article{BittracherKoltaiKlusetal., author = {Bittracher, Andreas and Koltai, P{\´e}ter and Klus, Stefan and Banisch, Ralf and Dellnitz, Michael and Sch{\"u}tte, Christof}, title = {Transition manifolds of complex metastable systems: Theory and data-driven computation of effective dynamics}, series = {Jounal of Nonlinear Science}, volume = {28}, journal = {Jounal of Nonlinear Science}, number = {2}, doi = {10.1007/s00332-017-9415-0}, pages = {471 -- 512}, abstract = {We consider complex dynamical systems showing metastable behavior but no local separation of fast and slow time scales. The article raises the question of whether such systems exhibit a low-dimensional manifold supporting its effective dynamics. For answering this question, we aim at finding nonlinear coordinates, called reaction coordinates, such that the projection of the dynamics onto these coordinates preserves the dominant time scales of the dynamics. We show that, based on a specific reducibility property, the existence of good low-dimensional reaction coordinates preserving the dominant time scales is guaranteed. Based on this theoretical framework, we develop and test a novel numerical approach for computing good reaction coordinates. The proposed algorithmic approach is fully local and thus not prone to the curse of dimension with respect to the state space of the dynamics. Hence, it is a promising method for data-based model reduction of complex dynamical systems such as molecular dynamics.}, language = {en} } @article{KoltaiWuNoeetal., author = {Koltai, Peter and Wu, Hao and No{\´e}, Frank and Sch{\"u}tte, Christof}, title = {Optimal data-driven estimation of generalized Markov state models for non-equilibrium dynamics}, series = {Computation}, volume = {6}, journal = {Computation}, number = {1}, publisher = {MDPI}, address = {Basel, Switzerland}, doi = {10.3390/computation6010022}, pages = {22}, language = {en} } @article{DibakdelRazodeSanchoetal., author = {Dibak, Manuel and del Razo, Mauricio J. and de Sancho, David and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations}, series = {Journal of Chemical Physics}, volume = {148}, journal = {Journal of Chemical Physics}, number = {21}, doi = {10.1063/1.5020294}, abstract = {Molecular dynamics (MD) simulations can model the interactions between macromolecules with high spatiotemporal resolution but at a high computational cost. By combining high-throughput MD with Markov state models (MSMs), it is now possible to obtain long time-scale behavior of small to intermediate biomolecules and complexes. To model the interactions of many molecules at large length scales, particle-based reaction-diffusion (RD) simulations are more suitable but lack molecular detail. Thus, coupling MSMs and RD simulations (MSM/RD) would be highly desirable, as they could efficiently produce simulations at large time and length scales, while still conserving the characteristic features of the interactions observed at atomic detail. While such a coupling seems straightforward, fundamental questions are still open: Which definition of MSM states is suitable? Which protocol to merge and split RD particles in an association/dissociation reaction will conserve the correct bimolecular kinetics and thermodynamics? In this paper, we make the first step toward MSM/RD by laying out a general theory of coupling and proposing a first implementation for association/dissociation of a protein with a small ligand (A + B ⇌ C). Applications on a toy model and CO diffusion into the heme cavity of myoglobin are reported.}, language = {en} } @article{vonKleistSchuetteZhang, author = {von Kleist, Max and Sch{\"u}tte, Christof and Zhang, Wei}, title = {Statistical analysis of the first passage path ensemble of jump processes}, series = {Journal of Statistical Physics}, volume = {170}, journal = {Journal of Statistical Physics}, doi = {10.1007/s10955-017-1949-x}, pages = {809 -- 843}, abstract = {The transition mechanism of jump processes between two different subsets in state space reveals important dynamical information of the processes and therefore has attracted considerable attention in the past years. In this paper, we study the first passage path ensemble of both discrete-time and continuous-time jump processes on a finite state space. The main approach is to divide each first passage path into nonreactive and reactive segments and to study them separately. The analysis can be applied to jump processes which are non-ergodic, as well as continuous-time jump processes where the waiting time distributions are non-exponential. In the particular case that the jump processes are both Markovian and ergodic, our analysis elucidates the relations between the study of the first passage paths and the study of the transition paths in transition path theory. We provide algorithms to numerically compute statistics of the first passage path ensemble. The computational complexity of these algorithms scales with the complexity of solving a linear system, for which efficient methods are available. Several examples demonstrate the wide applicability of the derived results across research areas.}, language = {en} } @article{KlusNueskeKoltaietal., author = {Klus, Stefan and N{\"u}ske, Feliks and Koltai, Peter and Wu, Hao and Kevrekidis, Ioannis and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {Data-driven model reduction and transfer operator approximation}, series = {Journal of Nonlinear Science}, volume = {28}, journal = {Journal of Nonlinear Science}, number = {3}, doi = {10.1007/s00332-017-9437-7}, pages = {985 -- 1010}, language = {en} } @article{DibakJdelRazoDeSanchoetal., author = {Dibak, Manuel and J. del Razo, Mauricio and De Sancho, David and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations}, series = {Journal of Chemical Physics}, volume = {148}, journal = {Journal of Chemical Physics}, number = {214107}, doi = {10.1063/1.5020294}, abstract = {Molecular dynamics (MD) simulations can model the interactions between macromolecules with high spatiotemporal resolution but at a high computational cost. By combining high-throughput MD with Markov state models (MSMs), it is now possible to obtain long-timescale behavior of small to intermediate biomolecules and complexes. To model the interactions of many molecules at large lengthscales, particle-based reaction-diffusion (RD) simulations are more suitable but lack molecular detail. Thus, coupling MSMs and RD simulations (MSM/RD) would be highly desirable, as they could efficiently produce simulations at large time- and lengthscales, while still conserving the characteristic features of the interactions observed at atomic detail. While such a coupling seems straightforward, fundamental questions are still open: Which definition of MSM states is suitable? Which protocol to merge and split RD particles in an association/dissociation reaction will conserve the correct bimolecular kinetics and thermodynamics? In this paper, we make the first step towards MSM/RD by laying out a general theory of coupling and proposing a first implementation for association/dissociation of a protein with a small ligand (A + B <--> C). Applications on a toy model and CO diffusion into the heme cavity of myoglobin are reported.}, language = {en} } @article{DjurdjevacConradHelfmannZonkeretal., author = {Djurdjevac Conrad, Natasa and Helfmann, Luzie and Zonker, Johannes and Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {Human mobility and innovation spreading in ancient times: a stochastic agent-based simulation approach}, series = {EPJ Data Science}, volume = {7}, journal = {EPJ Data Science}, number = {1}, edition = {EPJ Data Science}, publisher = {EPJ Data Science}, doi = {10.1140/epjds/s13688-018-0153-9}, pages = {24}, abstract = {Human mobility always had a great influence on the spreading of cultural, social and technological ideas. Developing realistic models that allow for a better understanding, prediction and control of such coupled processes has gained a lot of attention in recent years. However, the modeling of spreading processes that happened in ancient times faces the additional challenge that available knowledge and data is often limited and sparse. In this paper, we present a new agent-based model for the spreading of innovations in the ancient world that is governed by human movements. Our model considers the diffusion of innovations on a spatial network that is changing in time, as the agents are changing their positions. Additionally, we propose a novel stochastic simulation approach to produce spatio-temporal realizations of the spreading process that are instructive for studying its dynamical properties and exploring how different influences affect its speed and spatial evolution.}, language = {en} } @misc{CiccottiFerrarioSchuette, author = {Ciccotti, Giovanni and Ferrario, Mauro and Sch{\"u}tte, Christof}, title = {Molecular Dynamics vs. Stochastic Processes: Are We Heading Anywhere?}, series = {Special Issue: Understanding Molecular Dynamics via Stochastic Processes, Entropy}, volume = {20}, journal = {Special Issue: Understanding Molecular Dynamics via Stochastic Processes, Entropy}, number = {5}, doi = {10.3390/e20050348}, language = {en} } @article{KlusBittracherSchusteretal., author = {Klus, Stefan and Bittracher, Andreas and Schuster, Ingmar and Sch{\"u}tte, Christof}, title = {A kernel-based approach to molecular conformation analysis}, series = {Journal of Chemical Physics}, volume = {149}, journal = {Journal of Chemical Physics}, number = {24}, doi = {10.1063/1.5063533}, abstract = {We present a novel machine learning approach to understanding conformation dynamics of biomolecules. The approach combines kernel-based techniques that are popular in the machine learning community with transfer operator theory for analyzing dynamical systems in order to identify conformation dynamics based on molecular dynamics simulation data. We show that many of the prominent methods like Markov State Models, EDMD, and TICA can be regarded as special cases of this approach and that new efficient algorithms can be constructed based on this derivation. The results of these new powerful methods will be illustrated with several examples, in particular the alanine dipeptide and the protein NTL9.}, language = {en} } @article{KoltaiSchuette, author = {Koltai, P{\´e}ter and Sch{\"u}tte, Christof}, title = {A multiscale perturbation expansion approach for Markov state modeling of non-stationary molecular dynamics}, series = {SIAM J. Multiscale Model. Simul.}, volume = {16}, journal = {SIAM J. Multiscale Model. Simul.}, number = {4}, publisher = {SIAM}, doi = {10.1137/17M1146403}, pages = {1455 -- 1485}, abstract = {We investigate metastable dynamical systems subject to non-stationary forcing as they appear in molecular dynamics for systems driven by external fields. We show, that if the strength of the forcing is inversely proportional to the length of the slow metastable time scales of the unforced system, then the effective behavior of the forced system on slow time scales can be described by a low-dimensional reduced master equation. Our construction is explicit and uses the multiscale perturbation expansion method called two-timing, or method of multiple scales. The reduced master equation—a Markov state model—can be assembled by constructing two equilibrium Markov state models; one for the unforced system, and one for a slightly perturbed one.}, language = {en} } @article{KlusGelssPeitzetal., author = {Klus, Stefan and Gelß, Patrick and Peitz, Sebastian and Sch{\"u}tte, Christof}, title = {Tensor-based dynamic mode decomposition}, series = {Nonlinearity}, volume = {31}, journal = {Nonlinearity}, number = {7}, publisher = {IOP Publishing Ltd \& London Mathematical Society}, doi = {10.1088/1361-6544/aabc8f}, language = {en} } @article{DjurdjevacConradFuerstenauGrabundzijaetal., author = {Djurdjevac Conrad, Natasa and Fuerstenau, Daniel and Grabundzija, Ana and Helfmann, Luzie and Park, Martin and Schier, Wolfram and Sch{\"u}tt, Brigitta and Sch{\"u}tte, Christof and Weber, Marcus and Wulkow, Niklas and Zonker, Johannes}, title = {Mathematical modeling of the spreading of innovations in the ancient world}, series = {eTopoi. Journal for Ancient Studies}, volume = {7}, journal = {eTopoi. Journal for Ancient Studies}, issn = {ISSN 2192-2608}, doi = {10.17171/4-7-1}, language = {en} }