TY - JOUR A1 - Meerbach, E. A1 - Schütte, Christof A1 - Fischer, Alexander T1 - Eigenvalue Bounds on Restrictions of Reversible Nearly Uncoupled Markov Chains JF - Lin. Alg. Appl. Y1 - 2005 UR - http://publications.imp.fu-berlin.de/54/ U6 - https://doi.org/10.1016/j.laa.2004.10.018 VL - 398 SP - 141 EP - 160 ER - TY - JOUR A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - A Geometric Approach to Constrained Molecular Dynamics and Free Energy JF - Comm. Math. Sci. Y1 - 2005 UR - http://publications.imp.fu-berlin.de/48/ VL - 3 IS - 1 SP - 1 EP - 20 ER - TY - JOUR A1 - Antony, Jens A1 - Schmidt, Burkhard A1 - Schütte, Christof T1 - Nonadiabatic Effects on Peptide Vibrational Dynamics Induced by Conformational Changes JF - J. Chem. Phys. Y1 - 2005 UR - http://publications.imp.fu-berlin.de/57/ U6 - https://doi.org/10.1063/1.1829057 VL - 122 IS - 1 SP - 014309 ER - TY - JOUR A1 - Horenko, Illia A1 - Dittmer, E. A1 - Schütte, Christof T1 - Reduced Stochastic Models for Complex Molecular Systems JF - Comp. Vis. Sci. Y1 - 2005 UR - http://publications.imp.fu-berlin.de/59/ U6 - https://doi.org/10.1007/s00791-006-0021-1 VL - 9 IS - 2 SP - 89 EP - 102 ER - TY - JOUR A1 - Dellnitz, M. A1 - Neumann, M. A1 - Schütte, Christof T1 - Special Issue on Matrices and Mathematical Biology JF - Lin. Alg. Appl. Y1 - 2005 UR - http://publications.imp.fu-berlin.de/60/ U6 - https://doi.org/10.1016/j.laa.2004.12.006 VL - 398 SP - 1 EP - 245 ER - TY - CHAP A1 - Schütte, Christof A1 - Forster, R. A1 - Meerbach, E. A1 - Fischer, Alexander ED - Kornhuber, Ralf ED - Hoppe, Ronald H. W. ED - PÃ\copyrightriaux, J. ED - Pironneau, O. ED - Widlund, Olof ED - Xu, J. T1 - Uncoupling-Coupling Techniques for Metastable Dynamical Systems T2 - Domain Decomposition Methods in Science and Engineering Y1 - 2005 UR - http://publications.imp.fu-berlin.de/124/ VL - 40 SP - 115 EP - 129 PB - Springer ER - TY - JOUR A1 - Walter, J. A1 - Schütte, Christof A1 - Pavliotis, G. A1 - Stuart, Andrew T1 - Averaging of Stochastic Differential Equations JF - preprint Y1 - 2004 UR - http://publications.imp.fu-berlin.de/63/ ER - TY - JOUR A1 - Schütte, Christof A1 - Walter, J. A1 - Hartmann, Carsten A1 - Huisinga, Wilhelm T1 - An Averaging Principle for Fast Degrees of Freedom Exhibiting Long-Term Correlations JF - Multiscale Model. Simul. Y1 - 2004 UR - http://publications.imp.fu-berlin.de/62/ U6 - https://doi.org/10.1137/030600308 VL - 2 IS - 3 SP - 501 EP - 526 ER - TY - JOUR A1 - Preis, R. A1 - Dellnitz, M. A1 - Hessel, M. A1 - Schütte, Christof A1 - Meerbach, E. T1 - Dominant Paths Between Almost Invariant Sets of Dynamical Systems JF - preprint Y1 - 2004 UR - http://publications.imp.fu-berlin.de/64/ ER - TY - JOUR A1 - Horenko, Illia A1 - Weiser, Martin A1 - Schmidt, Burkhard A1 - Schütte, Christof T1 - Fully Adaptive Propagation of the Quantum-Classical Liouville Equation JF - J. Chem. Phys. Y1 - 2004 UR - http://publications.imp.fu-berlin.de/65/ U6 - https://doi.org/10.1063/1.1691015 VL - 120 IS - 19 SP - 8913 EP - 8923 ER - TY - CHAP A1 - Deuflhard, Peter A1 - Schütte, Christof ED - Hill, James ED - Moore, Ross T1 - Molecular Conformation Dynamics and Computational Drug Design T2 - Applied Mathematics Entering the 21st Century Y1 - 2004 UR - http://publications.imp.fu-berlin.de/125/ VL - 116 SP - 91 EP - 119 PB - SIAM ER - TY - JOUR A1 - Metzner, Ph. A1 - Dittmer, E. A1 - Jahnke, Tobias A1 - Schütte, Christof T1 - Generator Estimation of Markov Jump Processes JF - J. Comp. Phys. Y1 - 2007 UR - http://publications.imp.fu-berlin.de/35/ U6 - https://doi.org/10.1016/j.jcp.2007.07.032 VL - 227 IS - 1 SP - 353 EP - 375 ER - TY - JOUR A1 - Noé, Frank A1 - Horenko, Illia A1 - Schütte, Christof A1 - Smith, J. T1 - Hierarchical Analysis of Conformational Dynamics in Biomolecules JF - J. Chem. Phys. Y1 - 2007 UR - http://publications.imp.fu-berlin.de/37/ U6 - https://doi.org/10.1063/1.2714539 VL - 126 IS - 15 SP - 155102 ER - TY - JOUR A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Comment on Two Distinct Notions of Free Energy JF - Physica D Y1 - 2007 UR - http://publications.imp.fu-berlin.de/30/ U6 - https://doi.org/10.1016/j.physd.2007.02.006 VL - 228 IS - 1 SP - 59 EP - 63 ER - TY - JOUR A1 - Lorenz, S. A1 - Diederichs, Elmar A1 - Telgmann, R. A1 - Schütte, Christof T1 - Discrimination of Dynamical System Models for Biological and Chemical Processes JF - J. Comp. Chem. Y1 - 2007 UR - http://publications.imp.fu-berlin.de/33/ U6 - https://doi.org/10.1002/jcc.20674 VL - 28 IS - 8 SP - 1384 EP - 1399 ER - TY - CHAP A1 - Held, M. A1 - Meerbach, E. A1 - Hinderlich, S. A1 - Reutter, W. A1 - Schütte, Christof ED - Hansmann, U. ED - Meinke, J. ED - Mohanty, S. ED - Zimmermann, O. T1 - Conformational Studies of UDP-GlcNAc in Environments of Increasing Complexity T2 - From Computational Biophysics to Systems Biology Y1 - 2007 UR - http://publications.imp.fu-berlin.de/113/ VL - 36 SP - 145 EP - 148 PB - John von Neumann Institut for Computing ER - TY - JOUR A1 - Horenko, Illia A1 - Hartmann, Carsten A1 - Schütte, Christof A1 - Noé, Frank T1 - Data-based Parameter Estimation of Generalized Multidimensional Langevin Processes JF - Phys. Rev. E Y1 - 2007 UR - http://publications.imp.fu-berlin.de/31/ U6 - https://doi.org/10.1103/PhysRevE.76.016706 VL - 76 IS - 01 SP - 016706 ER - TY - JOUR A1 - Fischer, Alexander A1 - Waldhausen, S. A1 - Horenko, Illia A1 - Meerbach, E. A1 - Schütte, Christof T1 - Identification of Biomolecular Conformations from Incomplete Torsion Angle Observations by Hidden Markov Models JF - J. Comp. Chem. Y1 - 2007 UR - http://publications.imp.fu-berlin.de/38/ U6 - https://doi.org/10.1002/jcc.20692 VL - 28 IS - 15 SP - 2453 EP - 2464 ER - TY - CHAP A1 - Meerbach, E. A1 - Schütte, Christof A1 - Horenko, Illia A1 - Schmidt, Burkhard ED - Kühn, O. ED - Wöste, L. T1 - Metastable Conformational Structure and Dynamics T2 - Analysis and Control of Ultrafast Photoinduced Reactions Y1 - 2007 UR - http://publications.imp.fu-berlin.de/114/ U6 - https://doi.org/10.1007/978-3-540-68038-3_9 VL - 87 SP - 796 EP - 806 PB - Springer CY - Berlin ER - TY - JOUR A1 - Noé, Frank A1 - Smith, J. A1 - Schütte, Christof T1 - A network-based approach to biomolecular dynamics JF - From Computational Biophysics to Systems Biology (CBSB07). Editors Y1 - 2007 UR - http://publications.imp.fu-berlin.de/319/ VL - NIC Series 36 PB - John von Neumann Institute for Computing, Jülich ER - TY - JOUR A1 - Horenko, Illia A1 - Dittmer, E. A1 - Fischer, Alexander A1 - Schütte, Christof T1 - Automated Model Reduction for Complex Systems exhibiting Metastability JF - Mult. Mod. Sim. Y1 - 2006 UR - http://publications.imp.fu-berlin.de/44/ U6 - https://doi.org/10.1137/050623310 VL - 5 IS - 3 SP - 802 EP - 827 ER - TY - CHAP A1 - Schütte, Christof A1 - Nettesheim, Peter ED - Keil, F. ED - Mackens, W. ED - Voss, H. ED - Werther, J. T1 - Nonadiabatic Effects in Quantum-Classical Molecular Dynamics T2 - Scientific Computing in Chemical Engineering II Y1 - 1999 UR - http://publications.imp.fu-berlin.de/136/ SP - 42 EP - 56 PB - Springer ER - TY - CHAP A1 - Nettesheim, Peter A1 - Schütte, Christof ED - Deuflhard, Peter ED - Hermans, J. ED - Leimkuhler, Benedict ED - Marks, A. ED - Reich, Sebastian ED - Skeel, R. T1 - Numerical Integrators for Quantum-Classical Molecular Dynamics T2 - Computational Molecular Dynamics Y1 - 1999 UR - http://publications.imp.fu-berlin.de/137/ VL - 4 SP - 396 EP - 411 PB - Springer ER - TY - JOUR A1 - Bornemann, Folkmar A. A1 - Schütte, Christof T1 - On the Singular Limit of the Quantum-Classical Molecular Dynamics Model JF - J. Appl. Math. Y1 - 1999 UR - http://publications.imp.fu-berlin.de/94/ U6 - https://doi.org/10.1137/S0036139997318834 VL - 59 IS - 4 SP - 1208 EP - 1224 ER - TY - THES A1 - Schütte, Christof T1 - Partial Wigner transforms and the quantum-classical Liouville equation Y1 - 1999 UR - http://publications.imp.fu-berlin.de/95/ ER - TY - JOUR A1 - Fischer, Alexander A1 - Cordes, Frank A1 - Schütte, Christof T1 - Hybrid Monte Carlo with adaptive temperature in mixed-canonical ensemble JF - J. Comp. Chem. Y1 - 1998 UR - http://publications.imp.fu-berlin.de/97/ U6 - https://doi.org/10.1002/(SICI)1096-987X(19981130)19:15<1689::AID-JCC2>3.0.CO;2-J VL - 19 IS - 15 SP - 1689 EP - 1697 ER - TY - JOUR A1 - Bornemann, Folkmar A. A1 - Schütte, Christof T1 - A mathematical investigation of the Car-Parrinello Method JF - Num. Math. Y1 - 1998 UR - http://publications.imp.fu-berlin.de/96/ U6 - https://doi.org/10.1007/s002110050316 VL - 78 IS - 3 SP - 359 EP - 376 ER - TY - JOUR A1 - Schütte, Christof A1 - Bornemann, Folkmar A. T1 - Homogenization Approach to Smoothed Molecular Dynamics JF - Nonlinear Analysis Y1 - 1997 UR - http://publications.imp.fu-berlin.de/99/ U6 - https://doi.org/10.1016/S0362-546X(97)00216-2 VL - 30 IS - 3 SP - 1805 EP - 1814 ER - TY - JOUR A1 - Bornemann, Folkmar A. A1 - Schütte, Christof T1 - Homogenization of Hamiltonian Systems with a Strong Constraining Potential JF - Physica D Y1 - 1997 UR - http://publications.imp.fu-berlin.de/100/ U6 - https://doi.org/10.1016/S0167-2789(96)00245-X VL - 102 IS - 1-2 SP - 57 EP - 77 ER - TY - JOUR A1 - Nettesheim, Peter A1 - Huisinga, Wilhelm A1 - Schütte, Christof T1 - Chebyshev-Approximation for Wavepacket-Dynamics JF - preprint Y1 - 1996 UR - http://publications.imp.fu-berlin.de/102/ ER - TY - JOUR A1 - Nettesheim, Peter A1 - Bornemann, Folkmar A. A1 - Schmidt, Burkhard A1 - Schütte, Christof T1 - An Explicit and Symplectic Integrator for Quantum-Classical Molecular Dynamics JF - Chem. Phys. Lett. Y1 - 1996 UR - http://publications.imp.fu-berlin.de/101/ U6 - https://doi.org/10.1016/0009-2614(96)00471-X VL - 256 IS - 6 SP - 581 EP - 588 ER - TY - JOUR A1 - Huisinga, Wilhelm A1 - Meyn, S. A1 - Schütte, Christof T1 - Phase Transitions and Metastability in Markovian and Molecular Systems JF - Ann. Appl. Prob. Y1 - 2004 UR - http://publications.imp.fu-berlin.de/67/ U6 - https://doi.org/10.1214/aoap/1075828057 VL - 14 IS - 1 SP - 419 EP - 458 ER - TY - CHAP A1 - Schütte, Christof A1 - Huisinga, Wilhelm ED - Ciarlet, P. ED - Le Bris, Claude T1 - Biomolecular Conformations can be Identified as Metastable Sets of Molecular Dynamics T2 - Special Volume Y1 - 2003 UR - http://publications.imp.fu-berlin.de/126/ VL - X SP - 699 EP - 744 PB - Elsevier ER - TY - JOUR A1 - Huisinga, Wilhelm A1 - Schütte, Christof A1 - Stuart, Andrew T1 - Extracting Macroscopic Stochastic Dynamics JF - Comm. Pure Appl. Math. Y1 - 2003 UR - http://publications.imp.fu-berlin.de/72/ U6 - https://doi.org/10.1002/cpa.10057 VL - 56 IS - 2 SP - 234 EP - 269 ER - TY - CHAP A1 - Schütte, Christof A1 - Huisinga, Wilhelm A1 - Meyn, S. ED - Namachchivaya, N. ED - Lin, Y. T1 - Metastability of Diffusion Processes T2 - Nonlinear Stochastic Dynamics Y1 - 2003 UR - http://publications.imp.fu-berlin.de/127/ VL - 110 SP - 71 EP - 81 PB - Springer ER - TY - CHAP A1 - Fischer, Alexander A1 - Schütte, Christof A1 - Deuflhard, Peter A1 - Cordes, Frank ED - Schlick, T. ED - Gan, H. T1 - Hierarchical Uncoupling-Coupling of Metastable Conformations T2 - Computational Methods for Macromolecules Y1 - 2002 UR - http://publications.imp.fu-berlin.de/128/ VL - 24 SP - 235 EP - 259 PB - Springer ER - TY - JOUR A1 - Horenko, Illia A1 - Schmidt, Burkhard A1 - Schütte, Christof T1 - Multidimensional Classical Liouville Dynamics with Quantum Initial Conditions JF - J. Chem. Phys. Y1 - 2002 UR - http://publications.imp.fu-berlin.de/73/ U6 - https://doi.org/10.1063/1.1498467 VL - 117 IS - 10 SP - 4643 EP - 4650 ER - TY - JOUR A1 - Horenko, Illia A1 - Salzmann, Ch. A1 - Schmidt, Burkhard A1 - Schütte, Christof T1 - Quantum-Classical Liouville Approach to Molecular Dynamics JF - J. Chem. Phys. Y1 - 2002 UR - http://publications.imp.fu-berlin.de/75/ U6 - https://doi.org/10.1063/1.1522712 VL - 117 IS - 24 SP - 11075 EP - 11088 ER - TY - JOUR A1 - Horenko, Illia A1 - Schmidt, Burkhard A1 - Schütte, Christof T1 - A Theoretical Model for Molecules Interacting with Intense Laser Pulses JF - J. Chem. Phys. Y1 - 2001 UR - http://publications.imp.fu-berlin.de/79/ U6 - https://doi.org/10.1063/1.1398577 VL - 115 IS - 13 SP - 5733 EP - 5743 ER - TY - CHAP A1 - Schütte, Christof A1 - Huisinga, Wilhelm A1 - Deuflhard, Peter ED - Fiedler, B. T1 - Transfer Operator Approach to Conformational Dynamics in Biomolecular Systems T2 - Ergodic Theory, Analysis, and Efficient Simulation of Dynamical Systems Y1 - 2001 UR - http://publications.imp.fu-berlin.de/129/ SP - 191 EP - 223 PB - Springer ER - TY - CHAP A1 - Schütte, Christof A1 - Huisinga, Wilhelm ED - Sreenivas, R. ED - Jones, D. T1 - Biomolecular Conformations as Metastable Sets of Markov Chains T2 - Proceedings of the 38th Annual Allerton Conference on Communication, Control, and Computing, Monticello, Illinoins/USA Y1 - 2000 UR - http://publications.imp.fu-berlin.de/130/ SP - 1106 EP - 1115 PB - University of Illinois at Urbana-Champaign ER - TY - GEN A1 - Zhang, Wei A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Effective Dynamics Along Given Reaction Coordinates, and Reaction Rate Theory N2 - In molecular dynamics and related fields one considers dynamical descriptions of complex systems in full (atomic) detail. In order to reduce the overwhelming complexity of realistic systems (high dimension, large timescale spread, limited computational resources) the projection of the full dynamics onto some reaction coordinates is examined in order to extract statistical information like free energies or reaction rates. In this context, the effective dynamics that is induced by the full dynamics on the reaction coordinate space has attracted considerable attention in the literature. In this article, we contribute to this discussion: We first show that if we start with an ergodic diffusion processes whose invariant measure is unique then these properties are inherited by the effective dynamics. Then, we give equations for the effective dynamics, discuss whether the dominant timescales and reaction rates inferred from the effective dynamics are accurate approximations of such quantities for the full dynamics, and compare our findings to results from approaches like Zwanzig-Mori, averaging, or homogenization. Finally, by discussing the algorithmic realization of the effective dynamics, we demonstrate that recent algorithmic techniques like the ”equation-free” approach and the ”heterogeneous multiscale method” can be seen as special cases of our approach. T3 - ZIB-Report - 16-35 KW - Ergodic diffusion KW - reaction coordinate KW - effective dynamics KW - model reduction KW - equation-free KW - heterogeneous multiscale method Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-59706 SN - 1438-0064 ER - TY - GEN A1 - Ciccotti, Giovanni A1 - Ferrario, Mauro A1 - Schütte, Christof T1 - Molecular Dynamics vs. Stochastic Processes: Are We Heading Anywhere? T2 - Special Issue: Understanding Molecular Dynamics via Stochastic Processes, Entropy Y1 - 2018 U6 - https://doi.org/10.3390/e20050348 VL - 20 IS - 5 ER - TY - JOUR A1 - Zhang, Wei A1 - Schütte, Christof T1 - Reliable approximation of long relaxation timescales in molecular dynamics JF - Entropy Y1 - 2017 U6 - https://doi.org/10.3390/e19070367 VL - 19 IS - 7 ER - TY - JOUR A1 - Zhang, Wei A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Effective dynamics along given reaction coordinates, and reaction rate theory JF - Faraday Discussions Y1 - 2016 U6 - https://doi.org/10.1039/C6FD00147E IS - 195 SP - 365 EP - 394 ER - TY - JOUR A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Objective priors in the empirical Bayes framework JF - Scandinavian Journal of Statistics N2 - When dealing with Bayesian inference the choice of the prior often remains a debatable question. Empirical Bayes methods offer a data-driven solution to this problem by estimating the prior itself from an ensemble of data. In the nonparametric case, the maximum likelihood estimate is known to overfit the data, an issue that is commonly tackled by regularization. However, the majority of regularizations are ad hoc choices which lack invariance under reparametrization of the model and result in inconsistent estimates for equivalent models. We introduce a nonparametric, transformation-invariant estimator for the prior distribution. Being defined in terms of the missing information similar to the reference prior, it can be seen as an extension of the latter to the data-driven setting. This implies a natural interpretation as a trade-off between choosing the least informative prior and incorporating the information provided by the data, a symbiosis between the objective and empirical Bayes methodologies. Y1 - 2021 U6 - https://doi.org/10.1111/sjos.12485 VL - 48 IS - 4 SP - 1212 EP - 1233 PB - Wiley Online Library ER - TY - JOUR A1 - Hartmann, Carsten A1 - Richter, Lorenz A1 - Schütte, Christof A1 - Zhang, Wei T1 - Variational characterization of free energy: theory and algorithms JF - Entropy Y1 - 2017 U6 - https://doi.org/10.3390/e19110626 VL - 19 IS - 11 ER - TY - JOUR A1 - Dibak, Manuel A1 - J. del Razo, Mauricio A1 - De Sancho, David A1 - Schütte, Christof A1 - Noé, Frank T1 - MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations JF - Journal of Chemical Physics N2 - 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. Y1 - 2018 U6 - https://doi.org/10.1063/1.5020294 VL - 148 IS - 214107 ER - TY - JOUR A1 - Djurdjevac Conrad, Natasa A1 - Helfmann, Luzie A1 - Zonker, Johannes A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - Human mobility and innovation spreading in ancient times: a stochastic agent-based simulation approach JF - EPJ Data Science N2 - 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. Y1 - 2018 U6 - https://doi.org/10.1140/epjds/s13688-018-0153-9 VL - 7 IS - 1 PB - EPJ Data Science ET - EPJ Data Science ER - TY - JOUR A1 - Möller, Jan A1 - Isbilir, Ali A1 - Sungkaworn, Titiwat A1 - Osberg, Brenda A1 - Karathanasis, Christos A1 - Sunkara, Vikram A1 - Grushevsky, Eugene O A1 - Bock, Andreas A1 - Annibale, Paolo A1 - Heilemann, Mike A1 - Schütte, Christof A1 - Lohse, Martin J. T1 - Single molecule mu-opioid receptor membrane-dynamics reveal agonist-specific dimer formation with super-resolved precision JF - Nature Chemical Biology Y1 - 2020 U6 - https://doi.org/10.1038/s41589-020-0566-1 VL - 16 SP - 946 EP - 954 ER - TY - JOUR A1 - Bittracher, Andreas A1 - Schütte, Christof T1 - A probabilistic algorithm for aggregating vastly undersampled large Markov chains JF - Physica D: Nonlinear Phenomena Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-75874 VL - 416 ER - TY - JOUR A1 - Bittracher, Andreas A1 - Klus, Stefan A1 - Hamzi, Boumediene A1 - Schütte, Christof T1 - Dimensionality Reduction of Complex Metastable Systems via Kernel Embeddings of Transition Manifolds JF - Journal of Nonlinear Science N2 - We present a novel kernel-based machine learning algorithm for identifying the low-dimensional geometry of the effective dynamics of high-dimensional multiscale stochastic systems. Recently, the authors developed a mathematical framework for the computation of optimal reaction coordinates of such systems that is based on learning a parameterization of a low-dimensional transition manifold in a certain function space. In this article, we enhance this approach by embedding and learning this transition manifold in a reproducing kernel Hilbert space, exploiting the favorable properties of kernel embeddings. Under mild assumptions on the kernel, the manifold structure is shown to be preserved under the embedding, and distortion bounds can be derived. This leads to a more robust and more efficient algorithm compared to the previous parameterization approaches. Y1 - 2021 U6 - https://doi.org/10.1007/s00332-020-09668-z VL - 31 ER - TY - JOUR A1 - Klus, Stefan A1 - Nüske, Feliks A1 - Peitz, Sebastian A1 - Niemann, Jan-Hendrik A1 - Clementi, Cecilia A1 - Schütte, Christof T1 - Data-driven approximation of the Koopman generator: Model reduction, system identification, and control JF - Physica D: Nonlinear Phenomena Y1 - 2020 U6 - https://doi.org/10.1016/j.physd.2020.132416 VL - 406 ER - TY - BOOK A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - Stochastic Dynamics in Computational Biology T3 - Frontiers in Applied Dynamical Systems: Reviews and Tutorials Y1 - 2020 SN - 978-3-030-62386-9 U6 - https://doi.org/10.1007/978-3-030-62387-6 VL - 8 PB - Springer International Publishing ER - TY - JOUR A1 - Schulz, Robert A1 - Yamamoto, Kenji A1 - Klossek, André A1 - Rancan, Fiorenza A1 - Vogt, Annika A1 - Schütte, Christof A1 - Rühl, Eckart A1 - Netz, Roland R. T1 - Modeling of Drug Diffusion Based on Concentration Profiles in Healthy and Damaged Human Skin JF - Biophysical Journal N2 - Based on experimental drug concentration profiles in healthy as well as tape-stripped ex vivo human skin, we model the penetration of the antiinflammatory drug dexamethasone into the skin layers by the one-dimensional generalized diffusion equation. We estimate the position-dependent free-energy and diffusivity profiles by solving the conjugated minimization problem, in which the only inputs are concentration profiles of dexamethasone in skin at three consecutive penetration times. The resulting free-energy profiles for damaged and healthy skin show only minor differences. In contrast, the drug diffusivity in the first 10 μm of the upper skin layer of damaged skin is 200-fold increased compared to healthy skin, which reflects the corrupted barrier function of tape-stripped skin. For the case of healthy skin, we examine the robustness of our method by analyzing the behavior of the extracted skin parameters when the number of input and output parameters are reduced. We also discuss techniques for the regularization of our parameter extraction method. Y1 - 2019 U6 - https://doi.org/10.1016/j.bpj.2019.07.027 VL - 117 IS - 5 SP - 998 EP - 1008 ER - TY - JOUR A1 - Becker, Fabian A1 - Djurdjevac Conrad, Natasa A1 - Eser, Raphael A. A1 - Helfmann, Luzie A1 - Schütt, Brigitta A1 - Schütte, Christof A1 - Zonker, Johannes T1 - The Furnace and the Goat—A spatio-temporal model of the fuelwood requirement for iron metallurgy on Elba Island, 4th century BCE to 2nd century CE JF - PLOS ONE Y1 - 2020 UR - https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0241133 U6 - https://doi.org/10.1371/journal.pone.0241133 VL - 15 SP - 1 EP - 37 PB - Public Library of Science ER - TY - GEN A1 - Bittracher, Andreas A1 - Schütte, Christof ED - Junge, Oliver ED - Schütze, O. ED - Froyland, Gary ED - Ober-Blobaum, S. ED - Padberg-Gehle, E. T1 - A weak characterization of slow variables in stochastic dynamical systems T2 - Advances in Dynamics, Optimization and Computation. Series: Studies in Systems, Decision and Control. A volume dedicated to Michael Dellnitz on the occasion of his 60th birthday Y1 - 2020 SN - 978-3-030-51264-4 U6 - https://doi.org/10.1007/978-3-030-51264-4_6 VL - 304 SP - 132 EP - 150 PB - Springer International ER - TY - GEN A1 - Mollenhauer, Mattes A1 - Schuster, Ingmar A1 - Klus, Stefan A1 - Schütte, Christof ED - Junge, Oliver ED - Schütze, O. ED - Froyland, Gary ED - Ober-Blobaum, S. ED - Padberg-Gehle, K. T1 - Singular Value Decomposition of Operators on Reproducing Kernel Hilbert Spaces T2 - Advances om Dynamics, Optimization and Computation. Series: Studies in Systems, Decision and Control. A volume dedicated to Michael Dellnitz on his 60th birthday Y1 - 2020 SN - 978-3-030-51264-4 U6 - https://doi.org/10.1007/978-3-030-51264-4_5 VL - 304 SP - 109 EP - 131 PB - Springer International ER - TY - JOUR A1 - Helfmann, Luzie A1 - Ribera Borrell, Enric A1 - Schütte, Christof A1 - Koltai, Peter T1 - Extending Transition Path Theory: Periodically Driven and Finite-Time Dynamics JF - Journal of Nonlinear Science Y1 - 2020 U6 - https://doi.org/https://doi.org/10.1007/s00332-020-09652-7 VL - 30 SP - 3321 EP - 3366 ER - TY - JOUR A1 - Niemann, Jan-Hendrik A1 - Winkelmann, Stefanie A1 - Wolf, Sarah A1 - Schütte, Christof T1 - Agent-based modeling: Population limits and large timescales JF - Chaos: An Interdisciplinary Journal of Nonlinear Science N2 - Modeling, simulation and analysis of interacting agent systems is a broad field of research, with existing approaches reaching from informal descriptions of interaction dynamics to more formal, mathematical models. In this paper, we study agent-based models (ABMs) given as continuous-time stochastic processes and their pathwise approximation by ordinary and stochastic differential equations (ODEs and SDEs, respectively) for medium to large populations. By means of an appropriately adapted transfer operator approach we study the behavior of the ABM process on long time scales. We show that, under certain conditions, the transfer operator approach allows to bridge the gap between the pathwise results for large populations on finite timescales, i.e., the SDE limit model, and approaches built to study dynamical behavior on long time scales like large deviation theory. The latter provides a rigorous analysis of rare events including the associated asymptotic rates on timescales that scale exponentially with the population size. We demonstrate that it is possible to reveal metastable structures and timescales of rare events of the ABM process by finite-length trajectories of the SDE process for large enough populations. This approach has the potential to drastically reduce computational effort for the analysis of ABMs. Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-77309 SN - 1438-0064 VL - 31 IS - 3 ER - TY - JOUR A1 - Kostre, Margarita A1 - Schütte, Christof A1 - Noé, Frank A1 - del Razo Sarmina, Mauricio T1 - Coupling Particle-Based Reaction-Diffusion Simulations with Reservoirs Mediated by Reaction-Diffusion PDEs JF - Multiscale Modeling & Simulation N2 - Open biochemical systems of interacting molecules are ubiquitous in life-related processes. However, established computational methodologies, like molecular dynamics, are still mostly constrained to closed systems and timescales too small to be relevant for life processes. Alternatively, particle-based reaction-diffusion models are currently the most accurate and computationally feasible approach at these scales. Their efficiency lies in modeling entire molecules as particles that can diffuse and interact with each other. In this work, we develop modeling and numerical schemes for particle-based reaction-diffusion in an open setting, where the reservoirs are mediated by reaction-diffusion PDEs. We derive two important theoretical results. The first one is the mean-field for open systems of diffusing particles; the second one is the mean-field for a particle-based reaction-diffusion system with second-order reactions. We employ these two results to develop a numerical scheme that consistently couples particle-based reaction-diffusion processes with reaction-diffusion PDEs. This allows modeling open biochemical systems in contact with reservoirs that are time-dependent and spatially inhomogeneous, as in many relevant real-world applications. Y1 - 2021 U6 - https://doi.org/10.1137/20M1352739 VL - 19 IS - 4 SP - 1659 EP - 1683 PB - Society for Industrial and Applied Mathematics ER - TY - GEN A1 - Bittracher, Andreas A1 - Schütte, Christof T1 - A probabilistic algorithm for aggregating vastly undersampled large Markov chains N2 - Model reduction of large Markov chains is an essential step in a wide array of techniques for understanding complex systems and for efficiently learning structures from high-dimensional data. We present a novel aggregation algorithm for compressing such chains that exploits a specific low-rank structure in the transition matrix which, e.g., is present in metastable systems, among others. It enables the recovery of the aggregates from a vastly undersampled transition matrix which in practical applications may gain a speedup of several orders of mag- nitude over methods that require the full transition matrix. Moreover, we show that the new technique is robust under perturbation of the transition matrix. The practical applicability of the new method is demonstrated by identifying a reduced model for the large-scale traffic flow patterns from real-world taxi trip data. T3 - ZIB-Report - 20-21 Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-78688 UR - https://opus4.kobv.de/opus4-zib/frontdoor/index/index/docId/7587 SN - 1438-0064 ER - TY - JOUR A1 - Boltz, Horst-Holger A1 - Sirbu, Alexei A1 - Stelzer, Nina A1 - Lohse, Martin J. A1 - Schütte, Christof A1 - Annibale, Paolo T1 - Quantitative spectroscopy of single molecule interaction times JF - Optic Letters N2 - Single molecule fluorescence tracking provides information at nanometer-scale and millisecond-temporal resolution about the dynamics and interaction of individual molecules in a biological environment. While the dynamic behavior of isolated molecules can be characterized well, the quantitative insight is more limited when interactions between two indistinguishable molecules occur. We address this aspect by developing a theoretical foundation for a spectroscopy of interaction times, i.e., the inference of interaction from imaging data. A non-trivial crossover between a power law to an exponential behavior of the distribution of the interaction times is highlighted, together with the dependence of the exponential term upon the microscopic reaction affinity. Our approach is validated with simulated and experimental datasets. Y1 - 2021 U6 - https://doi.org/10.1364/OL.413030 VL - 46 IS - 7 SP - 1538 EP - 1541 ER - TY - GEN A1 - Sunkara, Vikram A1 - Raharinirina, N. Alexia A1 - Peppert, Felix A1 - von Kleist, Max A1 - Schütte, Christof T1 - Inferring Gene Regulatory Networks from Single Cell RNA-seq Temporal Snapshot Data Requires Higher Order Moments N2 - Due to the increase in accessibility and robustness of sequencing technology, single cell RNA-seq (scRNA-seq) data has become abundant. The technology has made significant contributions to discovering novel phenotypes and heterogeneities of cells. Recently, there has been a push for using single-- or multiple scRNA-seq snapshots to infer the underlying gene regulatory networks (GRNs) steering the cells' biological functions. To date, this aspiration remains unrealised. In this paper, we took a bottom-up approach and curated a stochastic two gene interaction model capturing the dynamics of a complete system of genes, mRNAs, and proteins. In the model, the regulation was placed upstream from the mRNA on the gene level. We then inferred the underlying regulatory interactions from only the observation of the mRNA population through~time. We could detect signatures of the regulation by combining information of the mean, covariance, and the skewness of the mRNA counts through time. We also saw that reordering the observations using pseudo-time did not conserve the covariance and skewness of the true time course. The underlying GRN could be captured consistently when we fitted the moments up to degree three; however, this required a computationally expensive non-linear least squares minimisation solver. There are still major numerical challenges to overcome for inference of GRNs from scRNA-seq data. These challenges entail finding informative summary statistics of the data which capture the critical regulatory information. Furthermore, the statistics have to evolve linearly or piece-wise linearly through time to achieve computational feasibility and scalability. T3 - ZIB-Report - 20-25 Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-79664 SN - 1438-0064 ER - TY - GEN A1 - Ray, Sourav A1 - Sunkara, Vikram A1 - Schütte, Christof A1 - Weber, Marcus T1 - How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs N2 - Molecular simulations of ligand-receptor interactions are a computational challenge, especially when their association- (``on''-rate) and dissociation- (``off''-rate) mechanisms are working on vastly differing timescales. In addition, the timescale of the simulations themselves is, in practice, orders of magnitudes smaller than that of the mechanisms; which further adds to the complexity of observing these mechanisms, and of drawing meaningful and significant biological insights from the simulation. One way of tackling this multiscale problem is to compute the free-energy landscapes, where molecular dynamics (MD) trajectories are used to only produce certain statistical ensembles. The approach allows for deriving the transition rates between energy states as a function of the height of the activation-energy barriers. In this article, we derive the association rates of the opioids fentanyl and N-(3-fluoro-1-phenethylpiperidin-4-yl)- N-phenyl propionamide (NFEPP) in a $\mu$-opioid receptor by combining the free-energy landscape approach with the square-root-approximation method (SQRA), which is a particularly robust version of Markov modelling. The novelty of this work is that we derive the association rates as a function of the pH level using only an ensemble of MD simulations. We also verify our MD-derived insights by reproducing the in vitro study performed by the Stein Lab, who investigated the influence of pH on the inhibitory constant of fentanyl and NFEPP (Spahn et al. 2017). MD simulations are far more accessible and cost-effective than in vitro and in vivo studies. Especially in the context of the current opioid crisis, MD simulations can aid in unravelling molecular functionality and assist in clinical decision-making; the approaches presented in this paper are a pertinent step forward in this direction. T3 - ZIB-Report - 20-18 KW - Opioid, Ligand-Receptor Interaction, Binding Kinetics, Molecular Dynamics, Metadynamics, SQRA Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-78437 SN - 1438-0064 ER - TY - JOUR A1 - Ray, Sourav A1 - Sunkara, Vikram A1 - Schütte, Christof A1 - Weber, Marcus T1 - How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs JF - Molecular Simulation N2 - Molecular simulations of ligand–receptor interactions are a computational challenge, especially when their association- (‘on’-rate) and dissociation- (‘off’-rate) mechanisms are working on vastly differing timescales. One way of tackling this multiscale problem is to compute the free-energy landscapes, where molecular dynamics (MD) trajectories are used to only produce certain statistical ensembles. The approach allows for deriving the transition rates between energy states as a function of the height of the activation-energy barriers. In this article, we derive the association rates of the opioids fentanyl and N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenyl propionamide (NFEPP) in a μ-opioid receptor by combining the free-energy landscape approach with the square-root-approximation method (SQRA), which is a particularly robust version of Markov modelling. The novelty of this work is that we derive the association rates as a function of the pH level using only an ensemble of MD simulations. We also verify our MD-derived insights by reproducing the in vitro study performed by the Stein Lab. Y1 - 2020 U6 - https://doi.org/10.1080/08927022.2020.1839660 VL - 46 IS - 18 SP - 1443 EP - 1452 PB - Taylor and Francis ER - TY - JOUR A1 - Wulkow, Niklas A1 - Koltai, Péter A1 - Schütte, Christof T1 - Memory-Based Reduced Modelling and Data-Based Estimation of Opinion Spreading JF - Journal of Nonlinear Science N2 - We investigate opinion dynamics based on an agent-based model and are interested in predicting the evolution of the percentages of the entire agent population that share an opinion. Since these opinion percentages can be seen as an aggregated observation of the full system state, the individual opinions of each agent, we view this in the framework of the Mori–Zwanzig projection formalism. More specifically, we show how to estimate a nonlinear autoregressive model (NAR) with memory from data given by a time series of opinion percentages, and discuss its prediction capacities for various specific topologies of the agent interaction network. We demonstrate that the inclusion of memory terms significantly improves the prediction quality on examples with different network topologies. Y1 - 2021 U6 - https://doi.org/10.1007/s00332-020-09673-2 VL - 31 ER - TY - JOUR A1 - Bittracher, Andreas A1 - Koltai, Péter A1 - Klus, Stefan A1 - Banisch, Ralf A1 - Dellnitz, Michael A1 - Schütte, Christof T1 - Transition manifolds of complex metastable systems: Theory and data-driven computation of effective dynamics JF - Jounal of Nonlinear Science N2 - 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. Y1 - 2018 U6 - https://doi.org/10.1007/s00332-017-9415-0 VL - 28 IS - 2 SP - 471 EP - 512 ER - TY - JOUR A1 - Koltai, Peter A1 - Wu, Hao A1 - Noé, Frank A1 - Schütte, Christof T1 - Optimal data-driven estimation of generalized Markov state models for non-equilibrium dynamics JF - Computation Y1 - 2018 U6 - https://doi.org/10.3390/computation6010022 VL - 6 IS - 1 PB - MDPI CY - Basel, Switzerland ER - TY - JOUR A1 - Dibak, Manuel A1 - del Razo, Mauricio J. A1 - de Sancho, David A1 - Schütte, Christof A1 - Noé, Frank T1 - MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations JF - Journal of Chemical Physics N2 - 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. Y1 - 2018 U6 - https://doi.org/10.1063/1.5020294 VL - 148 IS - 21 ER - TY - JOUR A1 - Djurdjevac Conrad, Natasa A1 - Fuerstenau, Daniel A1 - Grabundzija, Ana A1 - Helfmann, Luzie A1 - Park, Martin A1 - Schier, Wolfram A1 - Schütt, Brigitta A1 - Schütte, Christof A1 - Weber, Marcus A1 - Wulkow, Niklas A1 - Zonker, Johannes T1 - Mathematical modeling of the spreading of innovations in the ancient world JF - eTopoi. Journal for Ancient Studies Y1 - 2018 U6 - https://doi.org/10.17171/4-7-1 SN - ISSN 2192-2608 VL - 7 ER - TY - GEN A1 - Helfmann, Luzie A1 - Djurdjevac Conrad, Natasa A1 - Djurdjevac, Ana A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - From interacting agents to density-based modeling with stochastic PDEs N2 - Many real-world processes can naturally be modeled as systems of interacting agents. However, the long-term simulation of such agent-based models is often intractable when the system becomes too large. In this paper, starting from a stochastic spatio-temporal agent-based model (ABM), we present a reduced model in terms of stochastic PDEs that describes the evolution of agent number densities for large populations. We discuss the algorithmic details of both approaches; regarding the SPDE model, we apply Finite Element discretization in space which not only ensures efficient simulation but also serves as a regularization of the SPDE. Illustrative examples for the spreading of an innovation among agents are given and used for comparing ABM and SPDE models. T3 - ZIB-Report - 19-21 Y1 - 2019 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-73456 SN - 1438-0064 ER - TY - JOUR A1 - Helfmann, Luzie A1 - Djurdjevac Conrad, Natasa A1 - Djurdjevac, Ana A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - From interacting agents to density-based modeling with stochastic PDEs JF - Communications in Applied Mathematics and Computational Science N2 - Many real-world processes can naturally be modeled as systems of interacting agents. However, the long-term simulation of such agent-based models is often intractable when the system becomes too large. In this paper, starting from a stochastic spatio-temporal agent-based model (ABM), we present a reduced model in terms of stochastic PDEs that describes the evolution of agent number densities for large populations. We discuss the algorithmic details of both approaches; regarding the SPDE model, we apply Finite Element discretization in space which not only ensures efficient simulation but also serves as a regularization of the SPDE. Illustrative examples for the spreading of an innovation among agents are given and used for comparing ABM and SPDE models. Y1 - 2021 U6 - https://doi.org/10.2140/camcos.2021.16.1 VL - 16 IS - 1 SP - 1 EP - 32 ER - TY - JOUR A1 - Niemann, Jan-Hendrik A1 - Klus, Stefan A1 - Schütte, Christof T1 - Data-driven model reduction of agent-based systems using the Koopman generator JF - PLOS ONE N2 - The dynamical behavior of social systems can be described by agent-based models. Although single agents follow easily explainable rules, complex time-evolving patterns emerge due to their interaction. The simulation and analysis of such agent-based models, however, is often prohibitively time-consuming if the number of agents is large. In this paper, we show how Koopman operator theory can be used to derive reduced models of agent-based systems using only simulation or real-world data. Our goal is to learn coarse-grained models and to represent the reduced dynamics by ordinary or stochastic differential equations. The new variables are, for instance, aggregated state variables of the agent-based model, modeling the collective behavior of larger groups or the entire population. Using benchmark problems with known coarse-grained models, we demonstrate that the obtained reduced systems are in good agreement with the analytical results, provided that the numbers of agents is sufficiently large. Y1 - 2021 U6 - https://doi.org/10.1371/journal.pone.0250970 VL - 16 IS - 5 ER - TY - JOUR A1 - Wulkow, Hanna A1 - Conrad, Tim A1 - Djurdjevac Conrad, Natasa A1 - Müller, Sebastian A. A1 - Nagel, Kai A1 - Schütte, Christof T1 - Prediction of Covid-19 spreading and optimal coordination of counter-measures: From microscopic to macroscopic models to Pareto fronts JF - PLOS One Y1 - 2021 U6 - https://doi.org/10.1371/journal.pone.0249676 VL - 16 IS - 4 PB - Public Library of Science ER - TY - GEN A1 - Bornemann, Folkmar A. A1 - Schütte, Christof T1 - Adaptive Accuracy Control for Microcanonical Car-Parrinello Simulations N2 - The Car-Parrinello (CP) approach to ab initio molecular dynamics serves as an approximation to time-dependent Born-Oppenheimer (BO) calculations. It replaces the explicit minimization of the energy functional by a fictitious Newtonian dynamics and therefore introduces an artificial mass parameter $\mu$ which controls the electronic motion. A recent theoretical investigation shows that the CP-error, i.e., the deviation of the CP--solution from the BO-solution {\em decreases} like $\mu^{1/2}$ asymptotically. Since the computational effort {\em increases} like $\mu^{-1/2}$, the choice of $\mu$ has to find a compromise between efficiency and accuracy. The asymptotical result is used in this paper to construct an easily implemented algorithm which automatically controls $\mu$: the parameter $\mu$ is repeatedly adapted during the simulation by choosing $\mu$ as large as possible while pushing an error measure below a user-given tolerance. The performance and reliability of the algorithm is illustrated by a typical example. T3 - ZIB-Report - SC-96-20 Y1 - 1996 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-2311 ER - TY - GEN A1 - Bornemann, Folkmar A. A1 - Schütte, Christof T1 - Homogenization of Highly Oscillatory Hamiltonian Systems N2 - The paper studies Hamiltonian systems with a strong potential forcing the solutions to oscillate on a very small time scale. In particular, we are interested in the limit situation where the size $\epsilon$ of this small time scale tends to zero but the velocity components remain oscillating with an amplitude variation of order ${\rm O}(1)$. The process of establishing an effective initial value problem for the limit positions will be called {\em homogenization} of the Hamiltonian system. This problem occurs in mechanics as the problem of realization of holonomic constraints, in plasma physics as the problem of guiding center motion, in the simulation of biomolecules as the so called smoothing problem. We suggest the systematic use of the notion of {\em weak convergence} in order to approach this problem. This methodology helps to establish unified and short proofs of the known results which throw light on the inherent structure of the problem. Moreover, we give a careful and critical review of the literature. T3 - ZIB-Report - SC-95-39 Y1 - 1995 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-2050 ER - TY - GEN A1 - Bornemann, Folkmar A. A1 - Schütte, Christof T1 - On the Singular Limit of the Quantum-Classical Molecular Dynamics Model N2 - \noindent In molecular dynamics applications there is a growing interest in so-called {\em mixed quantum-classical} models. These models describe most atoms of the molecular system by the means of classical mechanics but an important, small portion of the system by the means of quantum mechanics. A particularly extensively used model, the QCMD model, consists of a {\em singularly perturbed}\/ Schrödinger equation nonlinearly coupled to a classical Newtonian equation of motion. This paper studies the singular limit of the QCMD model for finite dimensional Hilbert spaces. The main result states that this limit is given by the time-dependent Born-Oppenheimer model of quantum theory---provided the Hamiltonian under consideration has a smooth spectral decomposition. This result is strongly related to the {\em quantum adiabatic theorem}. The proof uses the method of {\em weak convergence} by directly discussing the density matrix instead of the wave functions. This technique avoids the discussion of highly oscillatory phases. On the other hand, the limit of the QCMD model is of a different nature if the spectral decomposition of the Hamiltonian happens not to be smooth. We will present a generic example for which the limit set is not a unique trajectory of a limit dynamical system but rather a {\em funnel} consisting of infinitely many trajectories. T3 - ZIB-Report - SC-97-07 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-2761 ER - TY - GEN A1 - Deuflhard, Peter A1 - Dellnitz, Michael A1 - Junge, Oliver A1 - Schütte, Christof T1 - Computation of Essential Molecular Dynamics by Subdivision Techniques I: Basic Concept N2 - The paper presents the concept of a new type of algorithm for the numerical computation of what the authors call the {\em essential dynamics\/} of molecular systems. Mathematically speaking, such systems are described by Hamiltonian differential equations. In the bulk of applications, individual trajectories are of no specific interest. Rather, time averages of physical observables or relaxation times of conformational changes need to be actually computed. In the language of dynamical systems, such information is contained in the natural invariant measure (infinite relaxation time) or in almost invariant sets ("large" finite relaxation times). The paper suggests the direct computation of these objects via eigenmodes of the associated Frobenius-Perron operator by means of a multilevel subdivision algorithm. The advocated approach is different to both Monte-Carlo techniques on the one hand and long term trajectory simulation on the other hand: in our setup long term trajectories are replaced by short term sub-trajectories, Monte-Carlo techniques are just structurally connected via the underlying Frobenius-Perron theory. Numerical experiments with a first version of our suggested algorithm are included to illustrate certain distinguishing properties. A more advanced version of the algorithm will be presented in a second part of this paper. T3 - ZIB-Report - SC-96-45 Y1 - 1996 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-2553 ER - TY - GEN A1 - Nettesheim, Peter A1 - Huisinga, Wilhelm A1 - Schütte, Christof T1 - Chebyshev-Approximation for Wavepacket-Dynamics: better than expected N2 - The aim of this work is to study the accuracy and stability of the Chebyshev--approximation method as a time--discretization for wavepacket dynamics. For this frequently used discretization we introduce estimates of the approximation and round--off error. These estimates mathematically confirm the stability of the Chebyshev--approximation with respect to round--off errors, especially for very large stepsizes. But the results also disclose threads to the stability due to large spatial dimensions. All theoretical statements are illustrated by numerical simulations of an analytically solvable example, the harmonic quantum oszillator. T3 - ZIB-Report - SC-96-47 Y1 - 1996 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-2575 ER - TY - GEN A1 - Bornemann, Folkmar A. A1 - Schütte, Christof T1 - A Mathematical Investigation of the Car-Parrinello Method N2 - The Car-Parrinello method for ab-initio molecular dynamics avoids the explicit minimization of energy functionals given by functional density theory in the context of the quantum adiabatic approximation (time-dependent Born-Oppenheimer approximation). Instead, it introduces a fictitious classical dynamics for the electronic orbitals. For many realistic systems this concept allowed first-principle computer simulations for the first time. In this paper we study the {\em quantitative} influence of the involved parameter $\mu$, the fictitious electronic mass of the method. In particular, we prove by use of a carefully chosen two-time-scale asymptotics that the deviation of the Car-Parrinello method from the adiabatic model is of order ${\rm O}(\mu^{1/2})$ --- provided one starts in the ground state of the electronic system and the electronic excitation spectrum satisfies a certain non-degeneracy condition. Analyzing a two-level model problem we prove that our result cannot be improved in general. Finally, we show how to use the gained quantitative insight for an automatic control of the unphysical ``fake'' kinetic energy of the method. T3 - ZIB-Report - SC-96-19 Y1 - 1996 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-2302 ER - TY - GEN A1 - Schütte, Christof A1 - Bornemann, Folkmar A. T1 - Homogenization Approach to Smoothed Molecular Dynamics N2 - {\footnotesize In classical Molecular Dynamics a molecular system is modelled by classical Hamiltonian equations of motion. The potential part of the corresponding energy function of the system includes contributions of several types of atomic interaction. Among these, some interactions represent the bond structure of the molecule. Particularly these interactions lead to extremely stiff potentials which force the solution of the equations of motion to oscillate on a very small time scale. There is a strong need for eliminating the smallest time scales because they are a severe restriction for numerical long-term simulations of macromolecules. This leads to the idea of just freezing the high frequency degrees of freedom (bond stretching and bond angles) via increasing the stiffness of the strong part of the potential to infinity. However, the naive way of doing this via holonomic constraints mistakenly ignores the energy contribution of the fast oscillations. The paper presents a mathematically rigorous discussion of the limit situation of infinite stiffness. It is demonstrated that the average of the limit solution indeed obeys a constrained Hamiltonian system but with a {\em corrected soft potential}. An explicit formula for the additive potential correction is given via a careful inspection of the limit energy of the fast oscillations. Unfortunately, the theory is valid only as long as the system does not run into certain resonances of the fast motions. Behind those resonances, there is no unique limit solution but a kind of choatic scenario for which the notion ``Takens chaos'' was coined. For demonstrating the relevance of this observation for MD, the theory is applied to a realistic, but still simple system: a single butan molecule. The appearance of ``Takens chaos'' in smoothed MD is illustrated and the consequences are discussed.} T3 - ZIB-Report - SC-96-31 Y1 - 1996 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-2410 ER - TY - GEN A1 - Bornemann, Folkmar A. A1 - Nettesheim, Peter A1 - Schütte, Christof T1 - Quantum-Classical Molecular Dynamics as an Approximation to Full Quantum Dynamics N2 - This paper presents a mathematical derivation of a model for quantum-classical molecular dynamics (QCMD) as a {\em partial} classical limit of the full Schrödinger equation. This limit is achieved in two steps: separation of the full wavefunction and short wave asymptotics for its ``classical'' part. Both steps can be rigorously justified under certain smallness assumptions. Moreover, the results imply that neither the time-dependent self-consistent field method nor mixed quantum-semi-classical models lead to better approximations than QCMD since they depend on the separation step, too. On the other hand, the theory leads to a characterization of the critical situations in which the models are in danger of largely deviating from the solution of the full Schrödinger equation. These critical situations are exemplified in an illustrative numerical simulation: the collinear collision of an Argon atom with a harmonic quantum oscillator. T3 - ZIB-Report - SC-95-26 Y1 - 1995 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1922 ER - TY - GEN A1 - Bornemann, Folkmar A. A1 - Schütte, Christof T1 - A Mathematical Approach to Smoothed Molecular Dynamics: Correcting Potentials for Freezing Bond Angles N2 - The interaction potential of molecular systems which are typically used in molecular dynamics can be split into two parts of essentially different stiffness. The strong part of the potential forces the solution of the equations of motion to oscillate on a very small time scale. There is a strong need for eliminating the smallest time scales because they are a severe restriction for numerical long-term simulations of macromolecules. This leads to the idea of just freezing the high frequency degrees of freedom (bond stretching and bond angles). However, the naive way of doing this via holonomic constraints is bound to produce incorrect results. The paper presents a mathematically rigorous discussion of the limit situation in which the stiffness of the strong part of the potential is increased to infinity. It is demonstrated that the average of the limit solution indeed obeys a constrained Hamiltonian system but with a {\em corrected soft potential}. An explicit formula for the additive potential correction is given and its significant contribution is demonstrated in an illustrative example. It appears that this correcting potential is definitely not identical with the Fixman-potential as was repeatedly assumed in the literature. T3 - ZIB-Report - SC-95-30 Y1 - 1995 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1960 ER - TY - GEN A1 - Schütte, Christof A1 - Dinand, Manfred A1 - Zumbusch, Gerhard A1 - Brinkmann, Ralf T1 - Dynamics of Erbium-doped Waveguide Lasers: Modelling, Reliable Simulation, and Comparison with Experiments N2 - A theoretical investigation of the dynamic properties of integrated optical Er--doped waveguide lasers is presented. It includes the construction of a physical model and of numerical techniques which allow reliable simulations of the dynamical behaviour of the laser signal depending on essential parameters of the laser device and on its external, time--dependent pump radiation. Therefore, a physical theory is developed which describes the propagation of light and its interaction with the active substrate in the laser cavity. This is realized in two steps. First, a {\em fundamental model} based on Maxwell's equations and on rate equations for the transitions in the active medium is constructed. Since this turns out to prohibit reliable simulations, it is, in a second step, reformulated via averaging in time and space which suppresses the fluctuations on the fastest time scales but represents them correctly. For this {\em reduced model} reliable and efficient simulation techniques using adaptive control schemes are designed and implemented. We apply the linear--implicit Euler discretization with extrapolation in time and a multilevel quadrature scheme in space. Finally, the model is justified in comparison with experimental observations in four cases of technological relevance. T3 - ZIB-Report - SC-95-19 Y1 - 1995 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1856 ER - TY - GEN A1 - Schütte, Christof T1 - Smoothed Molecular Dynamics for Thermally Embedded Systems N2 - This paper makes use of statistical mechanics in order to construct effective potentials for Molecular Dynamics for systems with nonstationary thermal embedding. The usual approach requires the computation of a statistical ensemble of trajectories. In the context of the new model the evaluation of only one single trajectory is sufficient for the determination of all interesting quantities, which leads to an enormous reduction of computational effort. This single trajectory is the solution to a corrected Hamiltonian system with a new potential $\tilde{V}$. It turns out that $\tilde{V}$ can be defined as spatial average of the original potential $V$. Therefore, the Hamiltonian dynamics defined by $\tilde{V}$ is smoother than that effected by $V$, i.e. a numerical integration of its evolution in time allows larger stepsizes. Thus, the presented approach introduces a Molecular Dynamics with smoothed trajectories originating from spatial averaging. This is deeply connected to time--averaging in Molecular Dynamics. These two types of {\em smoothed Molecular Dynamics} share advantages (gain in efficiency, reduction of error amplification, increased stability) and problems (necessity of closing relations and adaptive control schemes) which will be explained in detail. T3 - ZIB-Report - SC-95-14 Y1 - 1995 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1809 ER - TY - GEN A1 - Schütte, Christof T1 - A quasiresonant smoothing algorithm for the fast analysis of selective vibrational excitation. N2 - One key problem in modern chemistry is the simulation of the dynamical reaction of a molecule subjected to external radiation. This is described by the Schrödinger equation, which, after eigenfunction expansion, can be written in form of a system of ordinary differential equations, whose solutions show a highly oscillatory behaviour. The oscillations with high frequencies and small amplitudes confine the stepsizes of any numerical integrator -- an effect, which, in turn, blows up the simulation time. Larger stepsizes can be expected by averaging these fast oscillations, thus smoothing the trajectories. This idea leads to the construction of a quasiresonant smoothing algorithm (QRS). In QRS, a natural and computationally available splitting parameter $\delta$ controls the smoothing properties. The performance of QRS is demonstrated in two applications treating the selective excitation of vibrational states by picosecond laser pulses. In comparison with standard methods a speedup factor of 60--100 is observed. A closer look to purely physically motivated quasiresonant approximations such as WFQRA shows some additional advantages of the above smoothing idea. Among these the possibility of an adaptive formulation of QRS via the parameter $\delta$ is of particular importance. T3 - ZIB-Report - SC-93-10 Y1 - 1993 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1066 ER - TY - GEN A1 - Schütte, Christof A1 - Hohmann, Andreas A1 - Dinand, Manfred T1 - Numerical Simulation of Relaxation Oscillations of Waveguide-Lasers. N2 - An analysis of relaxation oscillations in local Er-doped optically pumped lasers is reported. It is based on a time dependent rate equation model for a quasi-two-level-system with wavelength dependent emission- and absorption cross-sections. For the first time a numerically reliable simulation of the characteristic laser behaviour was possible: the onset and decay of the oscillations, the time-dependent repetition period and the steady state signal output power. The characteristic waveguide parameters, as the erbium-concentration profile, the polarization dependent pump- and signal mode intensity profiles, the scattering losses, the cavity length and the front and rear reflectivities were all taken into account. The basic formulas are general and can also be used for Er-doped fiber lasers. Mathematically the problem can be characterized as a large boundary value problem, which can approximately be replaced by a stiff initial value problem of ordinary differential equations. The used algorithmic replacement procedure is motivated and discussed in detail. Here, pump- and signal evolution versus time are presented for an planar Er-diffused $\rm Ti$:$\rm LiNbO_{3}$ waveguide laser. The numerically obtained results show a nearly quantitative agreement with experimental investigations. Simultanously they supply knowledge about non-measureable (space-dependent population dynamic of the Er-atoms) and till today not measured data (dynamical response of the laser by a sharp peak in the external pump). T3 - ZIB-Report - SC-93-14 Y1 - 1993 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-1106 ER - TY - GEN A1 - Schütte, Christof A1 - Wulkow, Michael T1 - Quantum Theory with Discrete Spectra and Countable Systems of Differential Equations - A Numerical Treatment of RamanSpectroscopy. N2 - Models for occupation dynamics in discrete quantum systems lead to large or even infinite systems of ordinary differential equations. Some new mathematical techniques, developed for the simulation of chemical processes, make a numerical solution of countable systems of ordinary differential equations possible. Both, a basic physical concept for the construction of such systems and the structure of the numerical tools for solving them are presented. These conceptual aspects are illustrated by a simulation of an occupation process from spectroscopy. In this example the structures of rotation spectra observed in infrared spectroscopy are explained and some possibilities for an extension of the model are shown. T3 - ZIB-Report - SC-92-07 Y1 - 1992 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-778 ER - TY - GEN A1 - Schütte, Christof A1 - Nettesheim, Peter T1 - Non-Adiabatic Effects in Quantum-Classical Molecular Dynamics N2 - In molecular dynamics applications there is a growing interest in mixed quantum-classical models. The article is concerned with the so-called QCMD model. This model describes most atoms of the molecular system by the means of classical mechanics but an important, small portion of the system by the means of a wavefunction. We review the conditions under which the QCMD model is known to approximate the full quantum dynamical evolution of the system. In most quantum-classical simulations the {\em Born-Oppenheimer model} (BO) is used. In this model, the wavefunction is adiabatically coupled to the classical motion which leads to serious approximation deficiencies with respect to non-adiabatic effects in the fully quantum dynamical description of the system. In contrast to the BO model, the QCMD model does include non-adiabatic processes, e.g., transitions between the energy levels of the quantum system. It is demonstrated that, in mildly non-adiabatic scenarios, so-called {\em surface hopping} extensions of QCMD simulations yield good approximations of the non-adiabatic effects in full quantum dynamics. The algorithmic strategy of such extensions of QCMD is explained and the crucial steps of its realization are discussed with special emphasis on the numerical problems caused by highly oscillatory phase effects. T3 - ZIB-Report - SC-98-38 KW - quantum-classical molecular dynamics KW - non-adiabatic processes KW - Schrödinger equation KW - highly oscillatory phase KW - adiabatic limit KW - quantum adiabati Y1 - 1998 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3817 ER - TY - GEN A1 - Schütte, Christof A1 - Fischer, Alexander A1 - Huisinga, Wilhelm A1 - Deuflhard, Peter T1 - A Direct Approach to Conformational Dynamics based on Hybrid Monte Carlo N2 - Recently, a novel concept for the computation of essential features of the dynamics of Hamiltonian systems (such as molecular dynamics) has been proposed. The realization of this concept had been based on subdivision techniques applied to the Frobenius--Perron operator for the dynamical system. The present paper suggests an alternative but related concept that merges the conceptual advantages of the dynamical systems approach with the appropriate statistical physics framework. This approach allows to define the phrase ``conformation'' in terms of the dynamical behavior of the molecular system and to characterize the dynamical stability of conformations. In a first step, the frequency of conformational changes is characterized in statistical terms leading to the definition of some Markov operator $T$ that describes the corresponding transition probabilities within the canonical ensemble. In a second step, a discretization of $T$ via specific hybrid Monte Carlo techniques is shown to lead to a stochastic matrix $P$. With these theoretical preparations, an identification algorithm for conformations is applicable. It is demonstrated that the discretization of $T$ can be restricted to few essential degrees of freedom so that the combinatorial explosion of discretization boxes is prevented and biomolecular systems can be attacked. Numerical results for the n-pentane molecule and the triribonucleotide adenylyl\emph{(3'-5')}cytidylyl\emph{(3'-5')}cytidin are given and interpreted. T3 - ZIB-Report - SC-98-45 KW - conformation KW - conformational dynamics KW - hybrid Monte Carlo KW - reweighting KW - essential degrees of freedom KW - transition probabilities KW - Markov operator Y1 - 1999 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3889 ER - TY - GEN A1 - Nettesheim, Peter A1 - Schütte, Christof T1 - Numerical Integrators for Quantum-Classical Molecular Dynamics N2 - It was revealed that the QCMD model is of canonical Hamiltonian form with symplectic structure, which implies the conservation of energy. An efficient and reliable integrator for transfering these properties to the discrete solution is the symplectic and explicit {\sc Pickaback} algorithm. The only drawback of this kind of integrator is the small stepsize in time induced by the splitting techniques used to discretize the quantum evolution operator. Recent investigations concerning Krylov iteration techniques result in alternative approaches which overcome this difficulty for a wide range of problems. By using iterative methods in the evaluation of the quantum time propagator, these techniques allow for the stepsize to adapt to the coupling between the classical and the quantum mechanical subsystem. This yields a drastic reduction of the numerical effort. The pros and cons of both approaches as well as the suitable applications are discussed in the last part. T3 - ZIB-Report - SC-97-42 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3111 ER - TY - GEN A1 - Schütte, Christof A1 - Bornemann, Folkmar A. T1 - Approximation Properties and Limits of the Quantum-Classical Molecular Dynamics Model N2 - In molecular dynamics applications there is a growing interest in including quantum effects for simulations of larger molecules. This paper is concerned with {\em mixed quantum-classical} models which are currently discussed: the so-called QCMD model with variants and the time-dependent Born-Oppenheimer approximation. All these models are known to approximate the full quantum dynamical evolution---under different assumptions, however. We review the meaning of these assumptions and the scope of the approximation. In particular, we characterize those typical problematic situations where a mixed model might largely deviate from the full quantum evolution. One such situation of specific interest, a non-adiabatic excitation at certain energy level crossings, can promisingly be dealt with by a modification of the QCMD model that we suggest. T3 - ZIB-Report - SC-97-41 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3102 ER - TY - GEN A1 - Fischer, Alexander A1 - Cordes, Frank A1 - Schütte, Christof T1 - Hybrid Monte Carlo with Adaptive Temperature in a Mixed-Canonical Ensemble: Efficient Conformational Analysis of RNA N2 - A hybrid Monte Carlo method with adaptive temperature choice is presented, which exactly generates the distribution of a mixed-canonical ensemble composed of two canonical ensembles at low and high temperature. The analysis of resulting Markov chains with the reweighting technique shows an efficient sampling of the canonical distribution at low temperature, whereas the high temperature component facilitates conformational transitions, which allows shorter simulation times. \\The algorithm was tested by comparing analytical and numerical results for the small n-butane molecule before simulations were performed for a triribonucleotide. Sampling the complex multi-minima energy landscape of these small RNA segments, we observed enforced crossing of energy barriers. T3 - ZIB-Report - SC-97-67 Y1 - 1997 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3364 ER - TY - GEN A1 - Huisinga, Wilhelm A1 - Best, Christoph A1 - Cordes, Frank A1 - Roitzsch, Rainer A1 - Schütte, Christof T1 - From Simulation Data to Conformational Ensembles: Structure and Dynamics based Methods N2 - Statistical methods for analyzing large data sets of molecular configurations within the chemical concept of molecular conformations are described. The strategies are based on dependencies between configurations of a molecular ensemble; the article concentrates on dependencies induces by a) correlations between the molecular degrees of freedom, b) geometrical similarities of configurations, and c) dynamical relations between subsets of configurations. The statistical technique realizing aspect a) is based on an approach suggested by {\sc Amadei et al.} (Proteins, 17 (1993)). It allows to identify essential degrees of freedom of a molecular system and is extended in order to determine single configurations as representatives for the crucial features related to these essential degrees of freedom. Aspects b) and c) are based on statistical cluster methods. They lead to a decomposition of the available simulation data into {\em conformational ensembles} or {\em subsets} with the property that all configurations in one of these subsets share a common chemical property. In contrast to the restriction to single representative conformations, conformational ensembles include information about, e.g., structural flexibility or dynamical connectivity. The conceptual similarities and differences of the three approaches are discussed in detail and are illustrated by application to simulation data originating from a hybrid Monte Carlo sampling of a triribonucleotide. T3 - ZIB-Report - SC-98-36 KW - conformational ensemble KW - cluster method KW - structural and dynamical similarity KW - representative KW - conformation KW - essential degrees of freedom KW - transi Y1 - 1998 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3797 ER - TY - JOUR A1 - Zimper, Sebastian A1 - Djurdjevac, Ana A1 - Hartmann, Carsten A1 - Schütte, Christof A1 - Conrad, Natasa Djurdjevac T1 - Mean-field optimal control with stochastic leaders N2 - We consider interacting agent systems with a large number of stochastic agents (or particles) influenced by a fixed number of external stochastic lead agents. Such examples arise, for example in models of opinion dynamics, where a small number of leaders (influencers) can steer the behaviour of a large population of followers. In this context, we study a partial mean-field limit where the number of followers tends to infinity, while the number of leaders stays constant. The partial mean-field limit dynamics is then given by a McKean-Vlasov stochastic differential equation (SDE) for the followers, coupled to a controlled Itô-SDE governing the dynamics of the lead agents. For a given cost functional that the lead agents seek to minimise, we show that the unique optimal control of the finite agent system convergences to the optimal control of the limiting system. This establishes that the low-dimensional control of the partial (mean-field) system provides an effective approximation for controlling the high-dimensional finite agent system. In addition, we propose a stochastic gradient descent algorithm that can efficiently approximate the mean-field control. Our theoretical results are illustrated on opinion dynamics model with lead agents, where the control objective is to drive the followers to reach consensus in finite time. Y1 - 2025 ER - TY - JOUR A1 - Hartmann, Carsten A1 - Jöster, Annika A1 - Schütte, Christof A1 - Sikorski, Alexander A1 - Weber, Marcus T1 - Importance sampling of unbounded random stopping times: computing committor functions and exit rates without reweighting N2 - Rare events in molecular dynamics are often related to noise-induced transitions between different macroscopic states (e.g., in protein folding). A common feature of these rare transitions is that they happen on timescales that are on average exponentially long compared to the characteristic timescale of the system, with waiting time distributions that have (sub)exponential tails and infinite support. As a result, sampling such rare events can lead to trajectories that can be become arbitrarily long, with not too low probability, which makes the reweighting of such trajectories a real challenge. Here, we discuss rare event simulation by importance sampling from a variational perspective, with a focus on applications in molecular dynamics, in particular the computation of committor functions. The idea is to design importance sampling schemes that (a) reduce the variance of a rare event estimator while controlling the average length of the trajectories and (b) that do not require the reweighting of possibly very long trajectories. In doing so, we study different stochastic control formulations for committor and mean first exit times, which we compare both from a theoretical and a computational point of view, including numerical studies of some benchmark examples. Y1 - 2026 ER - TY - GEN A1 - Deuflhard, Peter A1 - Schütte, Christof T1 - Molecular Conformation Dynamics and Computational Drug Design N2 - The paper surveys recent progress in the mathematical modelling and simulation of essential molecular dynamics. Particular emphasis is put on computational drug design wherein time scales of $msec$ up to $min$ play the dominant role. Classical long-term molecular dynamics computations, however, would run into ill-conditioned initial value problems already after time spans of only $psec=10^{-12} sec$. Therefore, in order to obtain results for times of pharmaceutical interest, a combined deterministic-stochastic model is needed. The concept advocated in this paper is the direct identification of metastable conformations together with their life times and their transition patterns. It can be interpreted as a {\em transfer operator} approach corresponding to some underlying hybrid Monte Carlo process, wherein short-term trajectories enter. Once this operator has been discretized, which is a hard problem of its own, a stochastic matrix arises. This matrix is then treated by {\em Perron cluster analysis}, a recently developed cluster analysis method involving the numerical solution of an eigenproblem for a Perron cluster of eigenvalues. In order to avoid the 'curse of dimension', the construction of appropriate boxes for the spatial discretization of the Markov operator requires careful consideration. As a biomolecular example we present a rather recent SARS protease inhibitor. T3 - ZIB-Report - 03-20 KW - conformation dynamics KW - Monte Carlo methods KW - transfer operations KW - Hamiltonian dynamics KW - Smoluchowski dynamics KW - metastable sets KW - Perron cluster an Y1 - 2003 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-7427 ER - TY - THES A1 - Schütte, Christof T1 - A Quasiresonant Smoothing Algorithm for Solving Large Highly Differential Equations from Quantum Chemistry. N2 - In Quantum Chemistry the field of Laser--Assisted Molecular Control'' has received a considerable amount of attention recently. One key problem in this new field is the simulation of the dynamical reaction of a molecule subjected to external radiation. This problem is described by the Schrödinger equation, which, after eigenfunction expansion, can be written in the form of a large system of ordinary differential equations, the solutions of which show a highly oscillatory behaviour. The oscillations with high frequencies and small amplitudes confine the stepsizes of any numerical integrator -- an effect, which, in turn, blows up the simulation time. Larger stepsizes can be expected by averaging these fast oscillations, thus smoothing the trajectories. Standard smoothing techniques (averaging, filtering) would kill the whole process and thus, lead to wrong numerical results. To avoid this unwanted effect and nevertheless speed up computations, this paper presents a quasiresonant smoothing algorithm (QRS). In QRS, a natural splitting parameter $\delta$ controls the smoothing properties. An adaptive QRS--version (AQRS) is presented which includes an error estimation scheme for choosing this parameter $\delta$ in order to meet a given accuracy requirement. In AQRS $\delta$ is permanently adapted to the solution properties for computing the chemically necessary information'' only. The performance of AQRS is demonstrated in several test problems from the field Laser--Assisted Selective Excitation of Molecules'' in which the external radiation is a picosecond laser pulse. In comparison with standard methods speedup factors of the order of $10^2$ are observed. T3 - ZIB-Report - TR-94-04 Y1 - 1994 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-5090 ER - TY - GEN A1 - Schütte, Christof A1 - Huisinga, Wilhelm T1 - On Conformational Dynamics induced by Langevin Processes N2 - The function of many important biomolecules is related to their dynamic properties and their ability to switch between different {\em conformations}, which are understood as {\em almost invariant} or {\em metastable} subsets of the positional state space of the system. Recently, the present authors and their coworkers presented a novel algorithmic scheme for the direct numerical determination of such metastable subsets and the transition probability between them. Although being different in most aspects, this method exploits the same basic idea as {\sc Dellnitz} and {\sc Junge} in their approach to almost invariance in discrete dynamical systems: the almost invariant sets are computed via certain eigenvectors of the Markov operators associated with the dynamical behavior. In the present article we analyze the application of this approach to (high--friction) Langevin models describing the dynamical behavior of molecular systems coupled to a heat bath. We will see that this can be related to theoretical results for (symmetric) semigroups of Markov operators going back to {\sc Davies}. We concentrate on a comparison of our approach in respect to random perturbations of dynamical systems. T3 - ZIB-Report - SC-99-25 KW - Smoluchowski equation KW - Fokker--Planck equation KW - semigroup of Markov operators KW - canonical ensemble KW - small noise KW - first exit time KW - half time perio Y1 - 1999 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-4130 ER - TY - GEN A1 - Schütte, Christof T1 - Partial Wigner Transforms and the Quantum--Classical Liouville Equation N2 - In molecular dynamics applications there is a growing interest in mixed quantum-classical models. The {\em quantum-classical Liouville equation} (QCL) describes most atoms of the molecular system under consideration by means of classical phase space density but an important, small portion of the system by means of quantum mechanics. The QCL is derived from the full quantum dynamical (QD) description by applying the Wigner transform to the classical part'' of the system only. We discuss the conditions under which the QCL model approximates the full QD evolution of the system. First, analysis of the asymptotic properties of the Wigner transform shows that solving the QCL yields a first order approximation of full quantum dynamics. Second, we discuss the adiabatic limit of the QCL. This discussion shows that the QCL solutions may be interpretated as classical phase space densities, at least near the adiabatic limit. Third, it is demonstrated that the QCL yields good approximations of {\em non-adiabatic quantum effects,} especially near so-called {\em avoided crossings} where most quantum-classical models fail. T3 - ZIB-Report - SC-99-10 KW - QCMD KW - quantum-classical Liouville equation KW - surface hopping KW - Wigner transform KW - asymptotic expansion KW - nonadiabatic effects Y1 - 1999 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-3983 ER -