@misc{WinkelmannSchuettevonKleist2013, author = {Winkelmann, Stefanie and Sch{\"u}tte, Christof and von Kleist, Max}, title = {Markov Control Processes with Rare State Observation: Theory and Application to Treatment Scheduling in HIV-1}, issn = {1438-0064}, doi = {10.4310/CMS.2014.v12.n5.a4}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-41955}, year = {2013}, abstract = {Markov Decision Processes (MDP) or Partially Observable MDPs (POMDP) are used for modelling situations in which the evolution of a process is partly random and partly controllable. These MDP theories allow for computing the optimal control policy for processes that can continuously or frequently be observed, even if only partially. However, they cannot be applied if state observation is very costly and therefore rare (in time). We present a novel MDP theory for rare, costly observations and derive the corresponding Bellman equation. In the new theory, state information can be derived for a particular cost after certain, rather long time intervals. The resulting information costs enter into the total cost and thus into the optimization criterion. This approach applies to many real world problems, particularly in the medical context, where the medical condition is examined rather rarely because examination costs are high. At the same time, the approach allows for efficient numerical realization. We demonstrate the usefulness of the novel theory by determining, from the national economic perspective, optimal therapeutic policies for the treatment of the human immunodefficiency virus (HIV) in resource-rich and resource-poor settings. Based on the developed theory and models, we discover that available drugs may not be utilized efficiently in resource-poor settings due to exorbitant diagnostic costs.}, language = {en} } @misc{WeberFackeldeySchuette2017, author = {Weber, Marcus and Fackeldey, Konstantin and Sch{\"u}tte, Christof}, title = {Set-free Markov State Building}, issn = {1438-0064}, doi = {10.1063/1.4978501}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-62167}, year = {2017}, abstract = {Molecular dynamics (MD) simulations face challenging problems since the timescales of interest often are much longer than what is possible to simulate and even if sufficiently long simulation are possible the complex nature of the resulting simulation data makes interpretation difficult. Markov State Models (MSMs) help to overcome these problems by making experimentally relevant timescales accessible via coarse grained representations that also allows for convenient interpretation. However, standard set-based MSMs exhibit some caveats limiting their approximation quality and statistical significance. One of the main caveats results from the fact that typical MD trajectories repeatedly re-cross the boundary between the sets used to build the MSM which causes statistical bias in estimating the transition probabilities between these sets. In this article, we present a set-free approach to MSM building utilizing smooth overlapping ansatz functions instead of sets and an adaptive refinement approach. This kind of meshless discretization helps to overcome the recrossing problem and yields an adaptive refinement procedure that allows to improve the quality of the model while exploring state space and inserting new ansatz functions into the MSM.}, language = {en} } @misc{WinkelmannSchuette2016, author = {Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {The spatiotemporal master equation: approximation of reaction-diffusion dynamics via Markov state modeling}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-60999}, year = {2016}, abstract = {Accurate modeling and numerical simulation of reaction kinetics is a topic of steady interest. We consider the spatiotemporal chemical master equation (ST-CME) as a model for stochastic reaction-diffusion systems that exhibit properties of metastability. The space of motion is decomposed into metastable compartments and diffusive motion is approximated by jumps between these compartments. Treating these jumps as first-order reactions, simulation of the resulting stochastic system is possible by the Gillespie method. We present the theory of Markov state models (MSM) as a theoretical foundation of this intuitive approach. By means of Markov state modeling, both the number and shape of compartments and the transition rates between them can be determined. We consider the ST-CME for two reaction-diffusion systems and compare it to more detailed models. Moreover, a rigorous formal justification of the ST-CME by Galerkin projection methods is presented.}, language = {en} } @misc{KlebanovSikorskiSchuetteetal.2016, author = {Klebanov, Ilja and Sikorski, Alexander and Sch{\"u}tte, Christof and R{\"o}blitz, Susanna}, title = {Empirical Bayes Methods, Reference Priors, Cross Entropy and the EM Algorithm}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-61230}, year = {2016}, abstract = {When estimating a probability density within the empirical Bayes framework, the non-parametric maximum likelihood estimate (NPMLE) usually tends to overfit the data. This issue is usually taken care of by regularization - a penalization term is subtracted from the marginal log-likelihood before the maximization step, so that the estimate favors smooth solutions, resulting in the so-called maximum penalized likelihood estimation (MPLE). The majority of penalizations currently in use are rather arbitrary brute-force solutions, which lack invariance under transformation of the parameters(reparametrization) and measurements. This contradicts the principle that, if the underlying model has several equivalent formulations, the methods of inductive inference should lead to consistent results. Motivated by this principle and using an information-theoretic point of view, we suggest an entropy-based penalization term that guarantees this kind of invariance. The resulting density estimate can be seen as a generalization of reference priors. Using the reference prior as a hyperprior, on the other hand, is argued to be a poor choice for regularization. We also present an insightful connection between the NPMLE, the cross entropy and the principle of minimum discrimination information suggesting another method of inference that contains the doubly-smoothed maximum likelihood estimation as a special case.}, language = {en} } @misc{KlebanovSikorskiSchuetteetal.2016, author = {Klebanov, Ilja and Sikorski, Alexander and Sch{\"u}tte, Christof and R{\"o}blitz, Susanna}, title = {Empirical Bayes Methods for Prior Estimation in Systems Medicine}, issn = {1438-0064}, arxiv = {http://arxiv.org/abs/1612.01403}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-61307}, year = {2016}, abstract = {One of the main goals of mathematical modelling in systems medicine related to medical applications is to obtain patient-specific parameterizations and model predictions. In clinical practice, however, the number of available measurements for single patients is usually limited due to time and cost restrictions. This hampers the process of making patient-specific predictions about the outcome of a treatment. On the other hand, data are often available for many patients, in particular if extensive clinical studies have been performed. Therefore, before applying Bayes' rule separately to the data of each patient (which is typically performed using a non-informative prior), it is meaningful to use empirical Bayes methods in order to construct an informative prior from all available data. We compare the performance of four priors - a non-informative prior and priors chosen by nonparametric maximum likelihood estimation (NPMLE), by maximum penalized lilelihood estimation (MPLE) and by doubly-smoothed maximum likelihood estimation (DS-MLE) - by applying them to a low-dimensional parameter estimation problem in a toy model as well as to a high-dimensional ODE model of the human menstrual cycle, which represents a typical example from systems biology modelling.}, language = {en} } @misc{ZhangHartmannSchuette2016, author = {Zhang, Wei and Hartmann, Carsten and Sch{\"u}tte, Christof}, title = {Effective Dynamics Along Given Reaction Coordinates, and Reaction Rate Theory}, issn = {1438-0064}, doi = {10.1039/C6FD00147E}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-59706}, year = {2016}, abstract = {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.}, language = {en} } @misc{BittracherSchuette2020, author = {Bittracher, Andreas and Sch{\"u}tte, Christof}, title = {A probabilistic algorithm for aggregating vastly undersampled large Markov chains}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-78688}, year = {2020}, abstract = {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.}, language = {en} } @misc{SunkaraRaharinirinaPeppertetal.2020, author = {Sunkara, Vikram and Raharinirina, N. Alexia and Peppert, Felix and von Kleist, Max and Sch{\"u}tte, Christof}, title = {Inferring Gene Regulatory Networks from Single Cell RNA-seq Temporal Snapshot Data Requires Higher Order Moments}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-79664}, year = {2020}, abstract = {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.}, language = {en} } @misc{RaySunkaraSchuetteetal.2020, author = {Ray, Sourav and Sunkara, Vikram and Sch{\"u}tte, Christof and Weber, Marcus}, title = {How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs}, issn = {1438-0064}, doi = {10.1080/08927022.2020.1839660}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-78437}, year = {2020}, abstract = {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.}, language = {en} } @misc{HelfmannDjurdjevacConradDjurdjevacetal.2019, author = {Helfmann, Luzie and Djurdjevac Conrad, Natasa and Djurdjevac, Ana and Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {From interacting agents to density-based modeling with stochastic PDEs}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-73456}, year = {2019}, abstract = {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.}, language = {en} } @misc{BornemannSchuette1996, author = {Bornemann, Folkmar A. and Sch{\"u}tte, Christof}, title = {Adaptive Accuracy Control for Microcanonical Car-Parrinello Simulations}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2311}, number = {SC-96-20}, year = {1996}, abstract = {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.}, language = {en} } @misc{BornemannSchuette1995, author = {Bornemann, Folkmar A. and Sch{\"u}tte, Christof}, title = {Homogenization of Highly Oscillatory Hamiltonian Systems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2050}, number = {SC-95-39}, year = {1995}, abstract = {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.}, language = {en} } @misc{BornemannSchuette1997, author = {Bornemann, Folkmar A. and Sch{\"u}tte, Christof}, title = {On the Singular Limit of the Quantum-Classical Molecular Dynamics Model}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2761}, number = {SC-97-07}, year = {1997}, abstract = {\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{\"o}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.}, language = {en} } @misc{DeuflhardDellnitzJungeetal.1996, author = {Deuflhard, Peter and Dellnitz, Michael and Junge, Oliver and Sch{\"u}tte, Christof}, title = {Computation of Essential Molecular Dynamics by Subdivision Techniques I: Basic Concept}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2553}, number = {SC-96-45}, year = {1996}, abstract = {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.}, language = {en} } @misc{NettesheimHuisingaSchuette1996, author = {Nettesheim, Peter and Huisinga, Wilhelm and Sch{\"u}tte, Christof}, title = {Chebyshev-Approximation for Wavepacket-Dynamics: better than expected}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2575}, number = {SC-96-47}, year = {1996}, abstract = {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.}, language = {en} } @misc{BornemannSchuette1996, author = {Bornemann, Folkmar A. and Sch{\"u}tte, Christof}, title = {A Mathematical Investigation of the Car-Parrinello Method}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2302}, number = {SC-96-19}, year = {1996}, abstract = {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.}, language = {en} } @misc{SchuetteBornemann1996, author = {Sch{\"u}tte, Christof and Bornemann, Folkmar A.}, title = {Homogenization Approach to Smoothed Molecular Dynamics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-2410}, number = {SC-96-31}, year = {1996}, abstract = {{\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.}}, language = {en} } @misc{BornemannNettesheimSchuette1995, author = {Bornemann, Folkmar A. and Nettesheim, Peter and Sch{\"u}tte, Christof}, title = {Quantum-Classical Molecular Dynamics as an Approximation to Full Quantum Dynamics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-1922}, number = {SC-95-26}, year = {1995}, abstract = {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{\"o}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{\"o}dinger equation. These critical situations are exemplified in an illustrative numerical simulation: the collinear collision of an Argon atom with a harmonic quantum oscillator.}, language = {en} } @misc{BornemannSchuette1995, author = {Bornemann, Folkmar A. and Sch{\"u}tte, Christof}, title = {A Mathematical Approach to Smoothed Molecular Dynamics: Correcting Potentials for Freezing Bond Angles}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-1960}, number = {SC-95-30}, year = {1995}, abstract = {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.}, language = {en} } @misc{SchuetteDinandZumbuschetal.1995, author = {Sch{\"u}tte, Christof and Dinand, Manfred and Zumbusch, Gerhard and Brinkmann, Ralf}, title = {Dynamics of Erbium-doped Waveguide Lasers: Modelling, Reliable Simulation, and Comparison with Experiments}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-1856}, number = {SC-95-19}, year = {1995}, abstract = {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.}, language = {en} } @misc{Schuette1995, author = {Sch{\"u}tte, Christof}, title = {Smoothed Molecular Dynamics for Thermally Embedded Systems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-1809}, number = {SC-95-14}, year = {1995}, abstract = {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.}, language = {en} } @misc{Schuette1993, author = {Sch{\"u}tte, Christof}, title = {A quasiresonant smoothing algorithm for the fast analysis of selective vibrational excitation.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-1066}, number = {SC-93-10}, year = {1993}, abstract = {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{\"o}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.}, language = {en} } @misc{SchuetteHohmannDinand1993, author = {Sch{\"u}tte, Christof and Hohmann, Andreas and Dinand, Manfred}, title = {Numerical Simulation of Relaxation Oscillations of Waveguide-Lasers.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-1106}, number = {SC-93-14}, year = {1993}, abstract = {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).}, language = {en} } @misc{SchuetteWulkow1992, author = {Sch{\"u}tte, Christof and Wulkow, Michael}, title = {Quantum Theory with Discrete Spectra and Countable Systems of Differential Equations - A Numerical Treatment of RamanSpectroscopy.}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-778}, number = {SC-92-07}, year = {1992}, abstract = {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.}, language = {en} } @misc{SchuetteNettesheim1998, author = {Sch{\"u}tte, Christof and Nettesheim, Peter}, title = {Non-Adiabatic Effects in Quantum-Classical Molecular Dynamics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3817}, number = {SC-98-38}, year = {1998}, abstract = {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.}, language = {en} } @misc{SchuetteFischerHuisingaetal.1999, author = {Sch{\"u}tte, Christof and Fischer, Alexander and Huisinga, Wilhelm and Deuflhard, Peter}, title = {A Direct Approach to Conformational Dynamics based on Hybrid Monte Carlo}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3889}, number = {SC-98-45}, year = {1999}, abstract = {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.}, language = {en} } @misc{NettesheimSchuette1997, author = {Nettesheim, Peter and Sch{\"u}tte, Christof}, title = {Numerical Integrators for Quantum-Classical Molecular Dynamics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3111}, number = {SC-97-42}, year = {1997}, abstract = {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.}, language = {en} } @misc{SchuetteBornemann1997, author = {Sch{\"u}tte, Christof and Bornemann, Folkmar A.}, title = {Approximation Properties and Limits of the Quantum-Classical Molecular Dynamics Model}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3102}, number = {SC-97-41}, year = {1997}, abstract = {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.}, language = {en} } @misc{FischerCordesSchuette1997, author = {Fischer, Alexander and Cordes, Frank and Sch{\"u}tte, Christof}, title = {Hybrid Monte Carlo with Adaptive Temperature in a Mixed-Canonical Ensemble: Efficient Conformational Analysis of RNA}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3364}, number = {SC-97-67}, year = {1997}, abstract = {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.}, language = {en} } @misc{HuisingaBestCordesetal.1998, author = {Huisinga, Wilhelm and Best, Christoph and Cordes, Frank and Roitzsch, Rainer and Sch{\"u}tte, Christof}, title = {From Simulation Data to Conformational Ensembles: Structure and Dynamics based Methods}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3797}, number = {SC-98-36}, year = {1998}, abstract = {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.}, language = {en} } @misc{DeuflhardSchuette2003, author = {Deuflhard, Peter and Sch{\"u}tte, Christof}, title = {Molecular Conformation Dynamics and Computational Drug Design}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-7427}, number = {03-20}, year = {2003}, abstract = {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.}, language = {en} } @misc{SchuetteHuisinga1999, author = {Sch{\"u}tte, Christof and Huisinga, Wilhelm}, title = {On Conformational Dynamics induced by Langevin Processes}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4130}, number = {SC-99-25}, year = {1999}, abstract = {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.}, language = {en} } @misc{Schuette1999, author = {Sch{\"u}tte, Christof}, title = {Partial Wigner Transforms and the Quantum--Classical Liouville Equation}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3983}, number = {SC-99-10}, year = {1999}, abstract = {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.}, language = {en} } @misc{SchuetteHuisingaDeuflhard1999, author = {Sch{\"u}tte, Christof and Huisinga, Wilhelm and Deuflhard, Peter}, title = {Transfer Operator Approach to Conformational Dynamics in Biomolecular Systems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-4247}, number = {SC-99-36}, year = {1999}, abstract = {The article surveys the development of novel mathematical concepts and algorithmic approaches based thereon in view of their possible applicability to biomolecular design. Both a first deterministic approach, based on the Frobenius-Perron operator corresponding to the flow of the Hamiltonian dynamics, and later stochastic approaches, based on a spatial Markov operator or on Langevin dynamics, can be subsumed under the unified mathematical roof of the transfer operator approach to effective dynamics of molecular systems. The key idea of constructing specific transfer operators especially taylored for the purpose of conformational dynamics appears as the red line throughout the paper. Different steps of the algorithm are exemplified by a trinucleotide molecular system as a small representative of possible RNA drug molecules.}, language = {en} } @misc{DeuflhardHuisingaFischeretal.1998, author = {Deuflhard, Peter and Huisinga, Wilhelm and Fischer, Alexander and Sch{\"u}tte, Christof}, title = {Identification of Almost Invariant Aggregates in Reversible Nearly Uncoupled Markov Chains}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3469}, number = {SC-98-03}, year = {1998}, abstract = {The topic of the present paper bas been motivated by a recent computational approach to identify chemical conformations and conformational changes within molecular systems. After proper discretization, the conformations show up as almost invariant aggregates in reversible nearly uncoupled Markov chains. Most of the former work on this subject treated the direct problem: given the aggregates, analyze the loose coupling in connection with the computation of the stationary distribution (aggregation/disaggregation techniques). In contrast to that the present paper focuses on the inverse problem: given the system as a whole, identify the almost invariant aggregates together with the associated transition probabilites. A rather simple and robust algorithm is suggested and illustrated by its application to the n-pentane molecule.}, language = {en} } @misc{SchuetteFischerHuisingaetal.1998, author = {Sch{\"u}tte, Christof and Fischer, Alexander and Huisinga, Wilhelm and Deuflhard, Peter}, title = {A Hybrid Monte Carlo Method for Essential Molecular Dynamics}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-3474}, number = {SC-98-04}, year = {1998}, abstract = {Recently, a novel concept for the computation of essential features of Hamiltonian systems (such as those arising in molecular dynamics) has been proposed. The realization of that concept was based on subdivision techniques applied to the Frobenius--Perron operator for the dynamical system. The present paper suggests an alternative but related concept based on statistical mechanics, which allows to attack realistic molecular systems. 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 hybrid Monte Carlo techniques (based on short term subtrajectories only) is shown to lead to a stochastic matrix \$P\$. With these theoretical preparations, an identification algorithm for conformations is applicable (to be presented elsewhere). Numerical results for the n-pentane molecule are given and interpreted.}, language = {en} } @misc{WangHartmannSchuette2013, author = {Wang, Han and Hartmann, Carsten and Sch{\"u}tte, Christof}, title = {Linear response theory and optimal control for a molecular system under nonequilibrium conditions}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-18944}, year = {2013}, abstract = {In this paper, we propose a straightforward generalization of linear response theory to systems in nonequilibrium that are subject to nonequilibrium driving. We briefly revisit the standard linear response result for equilibrium systems, where we consider Langevin dynamics as a special case, and then give an alternative derivation using a change-of-measure argument that does not rely on any stationarity or reversibility assumption. This procedure moreover easily enables us to calculate the second order correction to the linear response formula (which may or may not be useful in practice). Furthermore, we outline how the novel nonequilibirum linear response formula can be used to compute optimal controls of molecular systems for cases in which one wants to steer the system to maximize a certain target expectation value. We illustrate our approach with simple numerical examples.}, language = {en} } @misc{WillenbockelSchuette2015, author = {Willenbockel, Christian Tobias and Sch{\"u}tte, Christof}, title = {A Variational Bayesian Algorithm for Clustering of Large and Complex Networks}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-54588}, year = {2015}, abstract = {We propose the Blockloading algorithm for the clustering of large and complex graphs with tens of thousands of vertices according to a Stochastic Block Model (SBM). Blockloading is based on generalized Variational Bayesian EM (VBEM) schemes and works for weighted and unweighted graphs. Existing Variational (Bayesian) EM methods have to consider each possible number of clusters sepa- rately to determine the optimal number of clusters and are prone to converge to local optima making multiple restarts necessary. These factors impose a severe restriction on the size and complexity of graphs these methods can handle. In con- trast, the Blockloading algorithm restricts restarts to subnetworks in a way that provides error correction of an existing cluster assignment. The number of clusters need not be specified in advance because Blockloading will return it as a result. We show that Blockloading outperforms all other variational methods regarding reliability of the results and computational efficiency.}, language = {en} } @misc{DjurdjevacConradWeberSchuette2015, author = {Djurdjevac Conrad, Natasa and Weber, Marcus and Sch{\"u}tte, Christof}, title = {Finding dominant structures of nonreversible Markov processes}, issn = {1438-0064}, doi = {10.1137/15M1032272}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55739}, year = {2015}, abstract = {Finding metastable sets as dominant structures of Markov processes has been shown to be especially useful in modeling interesting slow dynamics of various real world complex processes. Furthermore, coarse graining of such processes based on their dominant structures leads to better understanding and dimension reduction of observed systems. However, in many cases, e.g. for nonreversible Markov processes, dominant structures are often not formed by metastable sets but by important cycles or mixture of both. This paper aims at understanding and identifying these different types of dominant structures for reversible as well as nonreversible ergodic Markov processes. Our algorithmic approach generalizes spectral based methods for reversible process by using Schur decomposition techniques which can tackle also nonreversible cases. We illustrate the mathematical construction of our new approach by numerical experiments.}, language = {en} } @misc{GelssMateraSchuette2015, author = {Gelß, Patrick and Matera, Sebastian and Sch{\"u}tte, Christof}, title = {Solving the master equation without kinetic Monte Carlo: tensor train approximations for a CO oxidation model}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55743}, year = {2015}, abstract = {In multiscale models of heterogeneous catalysis, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. This usually is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO_2(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased stiffness.}, language = {en} } @misc{EncisoSchuetteDelleSite2015, author = {Enciso, Marta and Sch{\"u}tte, Christof and Delle Site, Luigi}, title = {Influence of pH and sequence in peptide aggregation via molecular simulation}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55805}, year = {2015}, abstract = {We employ a recently developed coarse-grained model for peptides and proteins where the effect of pH is automatically included. We explore the effect of pH in the aggregation process of the amyloidogenic peptide KTVIIE and two related sequences, using three different pH environments. Simulations using large systems (24 peptides chains per box) allow us to correctly account for the formation of realistic peptide aggregates. We evaluate the thermodynamic and kinetic implications of changes in sequence and pH upon peptide aggregation, and we discuss how a minimalistic coarse- grained model can account for these details.}, language = {en} } @misc{SchuetteSarich2015, author = {Sch{\"u}tte, Christof and Sarich, Marco}, title = {A Critical Appraisal of Markov State Models}, issn = {1438-0064}, doi = {10.1140/epjst/e2015-02421-0}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-54218}, year = {2015}, abstract = {Markov State Modelling as a concept for a coarse grained description of the essential kinetics of a molecular system in equilibrium has gained a lot of atten- tion recently. The last 10 years have seen an ever increasing publication activity on how to construct Markov State Models (MSMs) for very different molecular systems ranging from peptides to proteins, from RNA to DNA, and via molecu- lar sensors to molecular aggregation. Simultaneously the accompanying theory behind MSM building and approximation quality has been developed well be- yond the concepts and ideas used in practical applications. This article reviews the main theoretical results, provides links to crucial new developments, outlines the full power of MSM building today, and discusses the essential limitations still to overcome.}, language = {en} } @misc{BanischDjurdjevacConradSchuette2015, author = {Banisch, Ralf and Djurdjevac Conrad, Natasa and Sch{\"u}tte, Christof}, title = {Reactive flows and unproductive cycles for random walks on complex networks}, issn = {1438-0064}, doi = {10.1140/epjst/e2015-02417-8}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-54239}, year = {2015}, abstract = {We present a comprehensive theory for analysis and understanding of transition events between an initial set A and a target set B for general ergodic finite-state space Markov chains or jump processes, including random walks on networks as they occur, e.g., in Markov State Modelling in molecular dynamics. The theory allows us to decompose the probability flow generated by transition events between the sets A and B into the productive part that directly flows from A to B through reaction pathways and the unproductive part that runs in loops and is supported on cycles of the underlying network. It applies to random walks on directed networks and nonreversible Markov processes and can be seen as an extension of Transition Path Theory. Information on reaction pathways and unproductive cycles results from the stochastic cycle decomposition of the underlying network which also allows to compute their corresponding weight, thus characterizing completely which structure is used how often in transition events. The new theory is illustrated by an application to a Markov State Model resulting from weakly damped Langevin dynamics where the unproductive cycles are associated with periodic orbits of the underlying Hamiltonian dynamics.}, language = {en} } @misc{DjurdjevacConradBanischSchuette2014, author = {Djurdjevac Conrad, Natasa and Banisch, Ralf and Sch{\"u}tte, Christof}, title = {Modularity of Directed Networks: Cycle Decomposition Approach}, issn = {1438-0064}, doi = {10.3934/jcd.2015.2.1}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-51166}, year = {2014}, abstract = {The problem of decomposing networks into modules (or clusters) has gained much attention in recent years, as it can account for a coarsegrained description of complex systems, often revealing functional subunits of these systems. A variety of module detection algorithms have been proposed, mostly oriented towards finding hard partitionings of undirected networks. Despite the increasing number of fuzzy clustering methods for directed networks, many of these approaches tend to neglect important directional information. In this paper, we present a novel random walk based approach for finding fuzzy partitions of directed, weighted networks, where edge directions play a crucial role in defining how well nodes in a module are interconnected. We will show that cycle decomposition of a random walk process connects the notion of network modules and information transport in a network, leading to a new, symmetric measure of node communication. Finally, we will use this measure to introduce a communication graph, for which we will show that although being undirected it inherits all necessary information about modular structures from the original network.}, language = {en} } @misc{SarichSchuette2014, author = {Sarich, Marco and Sch{\"u}tte, Christof}, title = {Utilizing hitting times for finding metastable sets in non-reversible Markov chains}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-51209}, year = {2014}, abstract = {Techniques for finding metastable or almost invariant sets have been investigated, e.g., for deterministic dynamical systems in set-oriented numerics, for stochastic processes in molecular dynamics, and for random walks on complex networks. Most prominent algorithms are based on spectral apporaches and identify metastable sets via the doimant eigenvalues of the transfer operator associated with the dynamical system under consideration. These algorithms require the dominant eigenvalues to be real-valued. However, for many types of dynamics, e.g. for non-reversible Markov chains, this condition is not met. In this paper we utilize the hitting time apporach to metastable sets and demonstrate how the wellknown statements about optimal metastable decompositions of reversible chains can be reformulated for non-reversible chains if one switches from a spectral approach to an exit time approach. The performance of the resulting algorithm is illustrated by numerical experiments on random walks on complex networks.}, language = {en} } @misc{WangSchuette2014, author = {Wang, Han and Sch{\"u}tte, Christof}, title = {Building Markov State Models for Periodically Driven Non-Equilibrium Systems}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-53167}, year = {2014}, abstract = {Recent years have seen an increased interest in non-equilibrium molecular dynamics (NEMD) simulations, especially for molecular systems with periodic forcing by external fields, e.g., in the context of studying effects of electromagnetic radiation on the human body tissue. Lately, an NEMD methods with local thermostating has been proposed that allows for studying non-equilibrium processes in a statistically reliable and thermodynamically consistent way. In this article, we demonstrate how to construct Markov State Models (MSMs) for such NEMD simulations. MSM building has been well-established for systems in equilibrium where MSMs with just a few (macro-)states allow for accurate reproduction of the essential kinetics of the molecular system under consideration. Non-equilibrium MSMs have been lacking so far. The article presents how to construct such MSMs and illustrates their validity and usefulness for the case of conformation dynamics of alanine dipeptide in an external electric field.}, language = {en} } @misc{GulSchuetteBernhard2015, author = {Gul, Raheem and Sch{\"u}tte, Christof and Bernhard, Stefan}, title = {Mathematical modeling and sensitivity analysis of arterial anastomosis in arm arteries}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-54339}, year = {2015}, abstract = {Cardiovascular diseases are one of the major problems in medicine today and the number of patients increases worldwide. To find the most efficient treatment, prior knowledge about function and dysfunction of the cardiovas- cular system is required and methods need to be developed that identify the disease in an early stage. Mathematical modeling is a powerful tool for prediction and investigation of cardiovascular diseases. It has been shown that the Windkessel model, being based on an analogy between electrical circuits and fluid flow, is a simple but effective method to model the human cardiovascular system. In this paper, we have applied parametric local sensitivity analysis (LSA) to a linear elastic model of the arm arteries, to find and rank sensitive param- eters that may be helpful in clinical diagnosis. A computational model for end-to-side anastomosis (superior ulnar collateral anastomosis with posterior ulnar recurrent, SUC-PUR) is carried out to study the effects of some clinically relevant haemodynamic parameters like blood flow resistance and terminal re- sistance on pressure and flow at different locations of the arm artery. In this context, we also discuss the spatio-temporal dependency of local sensitivities. The sensitivities with respect to cardiovascular parameters reveal the flow resistance and diameter of the vessels as most sensitive parameters. These parameters play a key role in diagnosis of severe stenosis and aneurysms. In contrast, wall thickness and elastic modulus are found to be less sensitive.}, language = {en} } @misc{GuptaGramatkeEinspanieretal.2017, author = {Gupta, Pooja and Gramatke, Annika and Einspanier, Ralf and Sch{\"u}tte, Christof and von Kleist, Max and Sharbati, Jutta}, title = {In silicio cytotoxicity assessment on cultured rat intestinal cells deduced from cellular impedance measurements}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-62666}, year = {2017}, abstract = {Early and reliable identification of chemical toxicity is of utmost importance. At the same time, reduction of animal testing is paramount. Therefore, methods that improve the interpretability and usability of in vitro assays are essential. xCELLigence's real-time cell analyzer (RTCA) provides a novel, fast and cost effective in vitro method to probe compound toxicity. We developed a simple mathematical framework for the qualitative and quantitative assessment of toxicity for RTCA measurements. Compound toxicity, in terms of its 50\% inhibitory concentration IC_{50} on cell growth, and parameters related to cell turnover were estimated on cultured IEC-6 cells exposed to 10 chemicals at varying concentrations. Our method estimated IC50 values of 113.05, 7.16, 28.69 and 725.15 μM for the apparently toxic compounds 2-acetylamino-fluorene, aflatoxin B1, benzo-[a]-pyrene and chloramphenicol in the tested cell line, in agreement with literature knowledge. IC_{50} values of all apparent in vivo non-toxic compounds were estimated to be non-toxic by our method. Corresponding estimates from RTCA's in-built model gave false positive (toxicity) predictions in 5/10 cases. Taken together, our proposed method reduces false positive predictions and reliably identifies chemical toxicity based on impedance measurements. The source code for the developed method including instructions is available at https://git.zib.de/bzfgupta/toxfit/tree/master.}, language = {en} } @misc{WinkelmannSchuette2017, author = {Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {Hybrid Models for Chemical Reaction Networks: Multiscale Theory and Application to Gene Regulatory Systems}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-64264}, year = {2017}, abstract = {Well-mixed stochastic chemical kinetics are properly modelled by the chemical master equation (CME) and associated Markov jump processes in molecule number space. If the reactants are present in large amounts, however, corresponding simulations of the stochastic dynamics become computationally expensive and model reductions are demanded. The classical model reduction approach uniformly rescales the overall dynamics to obtain deterministic systems characterized by ordinary differential equations, the well-known mass action reaction rate equations. For systems with multiple scales there exist hybrid approaches that keep parts of the system discrete while another part is approximated either using Langevin dynamics or deterministically. This paper aims at giving a coherent overview of the different hybrid approaches, focusing on their basic concepts and the relation between them. We derive a novel general description of such hybrid models that allows to express various forms by one type of equation. We also check in how far the approaches apply to model extensions of the CME for dynamics which do not comply with the central well-mixed condition and require some spatial resolution. A simple but meaningful gene expression system with negative self-regulation is analysed to illustrate the different approximation qualities of some of the hybrid approaches discussed.}, language = {en} } @misc{BittracherKoltaiKlusetal.2017, 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}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-63822}, year = {2017}, 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} }