@article{SchuetteNielsenWeber2015, author = {Sch{\"u}tte, Christof and Nielsen, Adam and Weber, Marcus}, title = {Markov State Models and Molecular Alchemy}, volume = {113}, journal = {Molecular Physics}, number = {1}, doi = {10.1080/00268976.2014.944597}, pages = {69 -- 78}, year = {2015}, abstract = {In recent years Markov State Models (MSMs) have attracted a consid- erable amount of attention with regard to modelling conformation changes and associated function of biomolecular systems. They have been used successfully, e.g., for peptides including time-resolved spectroscopic experiments, protein function and protein folding , DNA and RNA, and ligand-receptor interaction in drug design and more complicated multivalent scenarios. In this article a novel reweighting scheme is introduced that allows to construct an MSM for certain molecular system out of an MSM for a similar system. This permits studying how molecular properties on long timescales differ between similar molecular systems without performing full molecular dynamics simulations for each system under con- sideration. The performance of the reweighting scheme is illustrated for simple test cases including one where the main wells of the respective energy landscapes are located differently and an alchemical transformation of butane to pentane where the dimension of the state space is changed.}, language = {en} } @misc{Lie2012, type = {Master Thesis}, author = {Lie, Han Cheng}, title = {A Meshless Method for Conformational Analysis}, year = {2012}, 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} } @article{NielsenWeber2015, author = {Nielsen, Adam and Weber, Marcus}, title = {Computing the nearest reversible Markov chain}, volume = {22}, journal = {Numerical Linear Algebra with Applications}, number = {3}, doi = {10.1002/nla.1967}, pages = {483 -- 499}, year = {2015}, abstract = {Reversible Markov chains are the basis of many applications. However, computing transition probabilities by a finite sampling of a Markov chain can lead to truncation errors. Even if the original Markov chain is reversible, the approximated Markov chain might be non-reversible and will lose important properties, like the real valued spectrum. In this paper, we show how to find the closest reversible Markov chain to a given transition matrix. It turns out that this matrix can be computed by solving a convex minimization problem.}, language = {en} } @book{OPUS4-5339, title = {Molecular Dynamics Simulation}, journal = {Entropy (special issue)}, editor = {Ciccotti, Giovanni and Ferrario, Mauro and Sch{\"u}tte, Christof}, publisher = {Multidisciplinary Digital Publishing Institute (MDPI)}, year = {2014}, language = {en} } @article{HartmannSchuetteZhang2016, author = {Hartmann, Carsten and Sch{\"u}tte, Christof and Zhang, Wei}, title = {Model reduction algorithms for optimal control and importance sampling of diffusions}, volume = {29}, journal = {Nonlinearity}, number = {8}, doi = {10.1088/0951-7715/29/8/2298}, pages = {2298 -- 2326}, year = {2016}, abstract = {We propose numerical algorithms for solving optimal control and importance sampling problems based on simplified models. The algorithms combine model reduction techniques for multiscale diffusions and stochastic optimization tools, with the aim of reducing the original, possibly high-dimensional problem to a lower dimensional representation of the dynamics, in which only a few relevant degrees of freedom are controlled or biased. Specifically, we study situations in which either a reaction coordinate onto which the dynamics can be projected is known, or situations in which the dynamics shows strongly localized behavior in the small noise regime. No explicit assumptions about small parameters or scale separation have to be made. We illustrate the approach with simple, but paradigmatic numerical examples.}, language = {en} } @article{KlusKoltaiSchuette2016, author = {Klus, Stefan and Koltai, Peter and Sch{\"u}tte, Christof}, title = {On the numerical approximation of the Perron-Frobenius and Koopman operator}, volume = {3}, journal = {Journal of Computational Dynamics}, number = {1}, doi = {10.3934/jcd.2016003}, pages = {51 -- 77}, year = {2016}, abstract = {Information about the behavior of dynamical systems can often be obtained by analyzing the eigenvalues and corresponding eigenfunctions of linear operators associated with a dynamical system. Examples of such operators are the Perron-Frobenius and the Koopman operator. In this paper, we will review di� fferent methods that have been developed over the last decades to compute � infinite-dimensional approximations of these in� finite-dimensional operators - in particular Ulam's method and Extended Dynamic Mode Decomposition (EDMD) - and highlight the similarities and di� fferences between these approaches. The results will be illustrated using simple stochastic di� fferential equations and molecular dynamics examples.}, language = {en} } @article{KlusNueskeKoltaietal.2018, author = {Klus, Stefan and N{\"u}ske, Feliks and Koltai, Peter and Wu, Hao and Kevrekidis, Ioannis and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {Data-driven model reduction and transfer operator approximation}, volume = {28}, journal = {Journal of Nonlinear Science}, number = {3}, doi = {10.1007/s00332-017-9437-7}, pages = {985 -- 1010}, year = {2018}, language = {en} } @article{RuedrichSarichSchuette2017, author = {R{\"u}drich, S. and Sarich, Marco and Sch{\"u}tte, Christof}, title = {Utilizing hitting times for finding metastable sets in non-reversible Markov chains}, journal = {Journal of Comp. Dynamics}, year = {2017}, language = {en} } @article{WeberFackeldeySchuette2017, author = {Weber, Marcus and Fackeldey, Konstantin and Sch{\"u}tte, Christof}, title = {Set-Free Markov State Model Building}, volume = {146}, journal = {Journal of Chemical Physics}, number = {12}, doi = {10.1063/1.4978501}, year = {2017}, language = {en} } @article{QuerLie2017, author = {Quer, Jannes and Lie, Han Cheng}, title = {Some connections between importance sampling and enhanced sampling methods in molecular dynamics}, journal = {Journal of Chemical Physics}, year = {2017}, abstract = {Enhanced sampling methods play an important role in molecular dynamics, because they enable the collection of better statistics of rare events that are important in many physical phenomena. We show that many enhanced sampling methods can be viewed as methods for performing importance sampling, by identifying important correspondences between the language of molecular dynamics and the language of probability theory. We illustrate these connections by highlighting the similarities between the rare event simulation method of Hartmann and Sch{\"u}tte (J. Stat. Mech. Theor. Exp., 2012), and the enhanced sampling method of Valsson and Parrinello (Phys. Rev. Lett. 113, 090601). We show that the idea of changing a probability measure is fundamental to both enhanced sampling and importance sampling.}, language = {en} } @misc{QuerLie2017, author = {Quer, Jannes and Lie, Han Cheng}, title = {Some connections between importance sampling and enhanced sampling methods in molecular dynamics}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-64289}, year = {2017}, abstract = {Enhanced sampling methods play an important role in molecular dynamics, because they enable the collection of better statistics of rare events that are important in many physical phenomena. We show that many enhanced sampling methods can be viewed as methods for performing importance sampling, by identifying important correspondences between the language of molecular dynamics and the language of probability theory. We illustrate these connections by highlighting the similarities between the rare event simulation method of Hartmann and Sch{\"u}tte (J. Stat. Mech. Theor. Exp., 2012), and the enhanced sampling method of Valsson and Parrinello (Phys. Rev. Lett. 113, 090601). We show that the idea of changing a probability measure is fundamental to both enhanced sampling and importance sampling.}, language = {en} } @article{FackeldeyNiknejadWeber2017, author = {Fackeldey, Konstantin and Niknejad, Amir and Weber, Marcus}, title = {Finding Metastabilities in Reversible Markov Chains based on Incomplete Sampling: Case of Molecular Simulation}, journal = {Special Matrices}, number = {5/1}, doi = {10.1515/spma-2017-0006}, pages = {73 -- 81}, year = {2017}, language = {en} } @article{FackeldeyWeber2017, author = {Fackeldey, Konstantin and Weber, Marcus}, title = {GenPCCA -- Markov State Models for Non-Equilibrium Steady States}, journal = {Big data clustering: Data preprocessing, variable selection, and dimension reduction. WIAS Report No. 29}, doi = {10.20347/WIAS.REPORT.29}, pages = {70 -- 80}, year = {2017}, language = {en} } @misc{ZhangWangHartmannetal.2014, author = {Zhang, Wei and Wang, Han and Hartmann, Carsten and Weber, Marcus and Sch{\"u}tte, Christof}, title = {Applications of the cross-entropy method to importance sampling and optimal control of diffusions}, issn = {1438-0064}, doi = {10.1137/14096493X}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-49720}, year = {2014}, abstract = {We study the cross-entropy method for diffusions. One of the results is a versatile cross-entropy algorithm that can be used to design efficient importance sampling strategies for rare events or to solve optimal control problems. The approach is based on the minimization of a suitable cross-entropy functional, with a parametric family of exponentially tilted probability distributions. We illustrate the new algorithm with several numerical examples and discuss algorithmic issues and possible extensions of the method.}, language = {en} } @misc{WeberQuer2015, author = {Weber, Marcus and Quer, Jannes}, title = {Estimating exit rates in rare event dynamical systems via extrapolation}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-56622}, year = {2015}, abstract = {In this article we present a new idea for approximating exit rates for diffusion processes living in a craggy landscape. We are especially interested in the exit rates of a process living in a metastable regions. Due to the fact that Monte Carlo simulations perform quite poor and are very computational expensive in this setting we create several similar situations with a smoothed potential. For this we introduce a new parameter \$\lambda \in [0,1]\$ (\$\lambda = 1\$ very smoothed potential, \$\lambda=0\$ original potential) into the potential which controls the influence the smoothing. We then sample the exit rate for different parameters \$\lambda\$ the exit rate from a given region. Due to the fact that \$\lambda\$ is connected to the exit rate we can use this dependency to approximate the real exit rate. The method can be seen as something between hyperdynamics and temperature accelerated MC.}, language = {en} } @article{BujotzekSchuettNielsenetal.2014, author = {Bujotzek, Alexander and Sch{\"u}tt, Ole and Nielsen, Adam and Fackeldey, Konstantin and Weber, Marcus}, title = {ZIBgridfree: Efficient Conformational Analysis by Partition-of-Unity Coupling}, volume = {52}, journal = {Journal of Mathematical Chemistry}, number = {3}, doi = {10.1007/s10910-013-0265-1}, pages = {781 -- 804}, year = {2014}, language = {de} } @misc{FackeldeyKoltaiNeviretal.2017, author = {Fackeldey, Konstantin and Koltai, P{\´e}ter and N{\´e}vir, Peter and Rust, Henning and Schild, Axel and Weber, Marcus}, title = {From Metastable to Coherent Sets - time-discretization schemes}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-66074}, year = {2017}, abstract = {Given a time-dependent stochastic process with trajectories x(t) in a space \$\Omega\$, there may be sets such that the corresponding trajectories only very rarely cross the boundaries of these sets. We can analyze such a process in terms of metastability or coherence. Metastable sets M are defined in space \$M\subset\Omega\$, coherent sets \$M(t)\subset\Omega\$ are defined in space and time. Hence, if we extend the space by the time-variable t, coherent sets are metastable sets in \$\Omega\times[0,\infty]\$. This relation can be exploited, because there already exist spectral algorithms for the identification of metastable sets. In this article we show that these well-established spectral algorithms (like PCCA+) also identify coherent sets of non-autonomous dynamical systems. For the identification of coherent sets, one has to compute a discretization (a matrix T) of the transfer operator of the process using a space-timediscretization scheme. The article gives an overview about different time-discretization schemes and shows their applicability in two different fields of application.}, language = {en} } @article{LorenzWolfBraunetal.2017, author = {Lorenz, Philipp and Wolf, Frederik and Braun, Jonas and Djurdjevac Conrad, Natasa and H{\"o}vel, Philipp}, title = {Capturing the Dynamics of Hashtag-Communities}, journal = {International Workshop on Complex Networks and their Applications}, doi = {10.1007/978-3-319-72150-7_33}, pages = {401 -- 413}, year = {2017}, language = {en} } @article{HeidaSikorskiWeber2022, author = {Heida, Martin and Sikorski, Alexander and Weber, Marcus}, title = {Consistency and order 1 convergence of cell-centered finite volume discretizations of degenerate elliptic problems in any space dimension}, journal = {SIAM Journal on Numerical Analysis}, doi = {10.20347/WIAS.PREPRINT.2913}, year = {2022}, abstract = {We study consistency of cell-centered finite difference methods for elliptic equations with degenerate coefficients in any space dimension \$d \geq 2\$. This results in order of convergence estimates in the natural weighted energy norm and in the weighted discrete \$L^2\$-norm on admissible meshes. The cells of meshes under consideration may be very irregular in size. We particularly allow the size of certain cells to remain bounded from below even in the asymptotic limit. For uniform meshes we show that the order of convergence is at least 1 in the energy semi-norm, provided the discrete and continuous solutions exist and the continuous solution has \$H^2\$ regularity.}, language = {en} } @article{SechiWeberHeyne2021, author = {Sechi, Renata and Weber, Marcus and Heyne, Karsten}, title = {MSM building and projection for the analysis of time-resolved spectra}, journal = {Proceedings in Applied Mathematics and Mechanics}, number = {21}, doi = {https://doi.org/10.1002/pamm.202100102}, url = {http://nbn-resolving.de/https://onlinelibrary.wiley.com/doi/10.1002/pamm.202100102}, pages = {e202100102}, year = {2021}, abstract = {Understanding the kinetics between the components of time-resolved spectra is a crucial step in the study of photo-activatedprocesses. However, modeling the kinetics requires usually some a priori knowledge about the system. In our approach, webuild a Markov State Model (MSM) from the spectral data, and obtain a Koopman transition matrix K(t). With genPCCA,an invariant subspace projection, we project the process into its metastable components. The result of the application of gen-PCCA is a transition matrix Kc(t), from which we can read the transition probability between the metastable components of the reaction. We discuss the application of this analysis method to the transient absorption spectrum of brominated Al-corrole}, language = {en} } @misc{Sikorski2024, author = {Sikorski, Alexander}, title = {ISOKANN.jl}, doi = {10.5281/zenodo.11519359}, year = {2024}, abstract = {Julia pacarkge implementing the ISOKANN algorithm for the computation of invariant subspaces of Koopman operators}, language = {en} }