@misc{FackeldeyRoeblitzScharkoietal.2011, author = {Fackeldey, Konstantin and R{\"o}blitz, Susanna and Scharkoi, Olga and Weber, Marcus}, title = {Soft Versus Hard Metastable Conformations in Molecular Simulations}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-13189}, number = {11-27}, year = {2011}, abstract = {Particle methods have become indispensible in conformation dynamics to compute transition rates in protein folding, binding processes and molecular design, to mention a few. Conformation dynamics requires at a decomposition of a molecule's position space into metastable conformations. In this paper, we show how this decomposition can be obtained via the design of either ``soft'' or ``hard'' molecular conformations. We show, that the soft approach results in a larger metastabilitiy of the decomposition and is thus more advantegous. This is illustrated by a simulation of Alanine Dipeptide.}, language = {en} } @misc{Fackeldey2010, author = {Fackeldey, Konstantin}, title = {Coupling Meshbased and Meshfree Methods by a Transfer Operator Approach}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-11755}, number = {10-12}, year = {2010}, abstract = {In contrast to the well known meshbased methods like the finite element method, meshfree methods do not rely on a mesh. However besides their great applicability, meshfree methods are rather time consuming. Thus, it seems favorable to combine both methods, by using meshfree methods only in a small part of the domain, where a mesh is disadvantageous, and a meshbased method for the rest of the domain. We motivate, that this coupling between the two simulation techniques can be considered as saddle point problem and show the stability of this coupling. Thereby a novel transfer operator is introduced, which interacts in the transition zone, where both methods coexist.}, language = {en} } @misc{WeberFackeldey2013, author = {Weber, Marcus and Fackeldey, Konstantin}, title = {Computing the Minimal Rebinding Effect Included in a Given Kinetics}, issn = {1438-0064}, doi = {10.1137/13091124X}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-17796}, year = {2013}, abstract = {The rebinding effect is a phenomenon which occurs when observing a ligand-receptor binding process. On the macro scale this process comprises the Markov property. This Makovian view is spoiled when switching to the atomistic scale of a binding process. We therefore suggest a model which accurately describes the rebinding effect on the atomistic scale by allowing ''intermediate'' bound states. This allows us to define an indicator for the magnitude of rebinding and to formulate an optimization problem. The results form our examples show good agreement with data form laboratory.}, language = {en} } @misc{LieFackeldeyWeber2013, author = {Lie, Han Cheng and Fackeldey, Konstantin and Weber, Marcus}, title = {A square root approximation of transition rates for a Markov State Model}, issn = {1438-0064}, doi = {10.1137/120899959}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-42195}, year = {2013}, abstract = {Trajectory- or mesh-based methods for analyzing the dynamical behavior of large molecules tend to be impractical due to the curse of dimensionality - their computational cost increases exponentially with the size of the molecule. We propose a method to break the curse by a novel square root approximation of transition rates, Monte Carlo quadrature and a discretization approach based on solving linear programs. With randomly sampled points on the molecular energy landscape and randomly generated discretizations of the molecular configuration space as our initial data, we construct a matrix describing the transition rates between adjacent discretization regions. This transition rate matrix yields a Markov State Model of the molecular dynamics. We use Perron cluster analysis and coarse-graining techniques in order to identify metastable sets in configuration space and approximate the transition rates between the metastable sets. Application of our method to a simple energy landscape on a two-dimensional configuration space provides proof of concept and an example for which we compare the performance of different discretizations. We show that the computational cost of our method grows only polynomially with the size of the molecule. However, finding discretizations of higher-dimensional configuration spaces in which metastable sets can be identified remains a challenge.}, language = {en} } @misc{NielsenFackeldeyWeber2013, author = {Nielsen, Adam and Fackeldey, Konstantin and Weber, Marcus}, title = {On a Generalized Transfer Operator}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-43162}, year = {2013}, abstract = {We introduce a generalized operator for arbitrary stochastic processes by using a pre-kernel, which is a generalization of the Markov kernel. For deterministic processes, such an operator is already known as the Frobenius-Perron operator, which is defined for a large class of measures. For Markov processes, there exists transfer operators being only well defined for stationary measures in \$L^2\$. Our novel generalized transfer operator is well defined for arbitrary stochastic processes, in particular also for deterministic ones. We can show that this operator is acting on \$L^1\$. For stationary measures, this operator is also an endomorphism of \$L^2\$ and, therefore, allows for a mathematical analysis in Hilbert spaces.}, language = {en} } @article{FackeldeyKoltaiNeviretal.2019, author = {Fackeldey, Konstantin and Koltai, Peter and Nevir, Peter and Rust, Henning and Schild, Axel and Weber, Marcus}, title = {From metastable to coherent sets - Time-discretization schemes}, volume = {29}, journal = {Chaos: An Interdisciplinary Journal of Nonlinear Science}, doi = {10.1063/1.5058128}, pages = {012101 -- 012101}, year = {2019}, abstract = {In this article, we show that these well-established spectral algorithms (like PCCA+, Perron Cluster Cluster Analysis) 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-time-discretization scheme. The article gives an overview about different time-discretization schemes and shows their applicability in two different fields of application.}, language = {en} } @article{ReuterFackeldeyWeber2019, author = {Reuter, Bernhard and Fackeldey, Konstantin and Weber, Marcus}, title = {Generalized Markov modeling of nonreversible molecular kinetics}, volume = {17}, journal = {The Journal of Chemical Physics}, number = {150}, doi = {10.1063/1.5064530}, pages = {174103}, year = {2019}, abstract = {Markov state models are to date the gold standard for modeling molecular kinetics since they enable the identification and analysis of metastable states and related kinetics in a very instructive manner. The state-of-the-art Markov state modeling methods and tools are very well developed for the modeling of reversible processes in closed equilibrium systems. On the contrary, they are largely not well suited to deal with nonreversible or even nonautonomous processes of nonequilibrium systems. Thus, we generalized the common Robust Perron Cluster Cluster Analysis (PCCA+) method to enable straightforward modeling of nonequilibrium systems as well. The resulting Generalized PCCA (G-PCCA) method readily handles equilibrium as well as nonequilibrium data by utilizing real Schur vectors instead of eigenvectors. This is implemented in the G-PCCA algorithm that enables the semiautomatic coarse graining of molecular kinetics. G-PCCA is not limited to the detection of metastable states but also enables the identification and modeling of cyclic processes. This is demonstrated by three typical examples of nonreversible systems.}, language = {en} } @article{ErnstFackeldeyVolkameretal.2019, author = {Ernst, Natalia and Fackeldey, Konstantin and Volkamer, Andrea and Opatz, Oliver and Weber, Marcus}, title = {Computation of temperature-dependent dissociation rates of metastable protein-ligand complexes}, volume = {45}, journal = {Molecular Simulation}, number = {11}, doi = {10.1080/08927022.2019.1610949}, pages = {904 -- 911}, year = {2019}, abstract = {Molecular simulations are often used to analyse the stability of protein-ligand complexes. The stability can be characterised by exit rates or using the exit time approach, i.e. by computing the expected holding time of the complex before its dissociation. However determining exit rates by straightforward molecular dynamics methods can be challenging for stochastic processes in which the exit event occurs very rarely. Finding a low variance procedure for collecting rare event statistics is still an open problem. In this work we discuss a novel method for computing exit rates which uses results of Robust Perron Cluster Analysis (PCCA+). This clustering method gives the possibility to define a fuzzy set by a membership function, which provides additional information of the kind 'the process is being about to leave the set'. Thus, the derived approach is not based on the exit event occurrence and, therefore, is also applicable in case of rare events. The novel method can be used to analyse the temperature effect of protein-ligand systems through the differences in exit rates, and, thus, open up new drug design strategies and therapeutic applications.}, language = {en} } @article{LieFackeldeyWeber2013, author = {Lie, Han Cheng and Fackeldey, Konstantin and Weber, Marcus}, title = {A Square Root Approximation of Transition Rates for a Markov State Model}, volume = {34}, journal = {SIAM. J. Matrix Anal. Appl.}, number = {2}, doi = {10.1137/120899959}, pages = {738 -- 756}, year = {2013}, language = {en} } @article{ScharkoiFackeldeyMerkulowetal.2013, author = {Scharkoi, Olga and Fackeldey, Konstantin and Merkulow, Igor and Andrae, Karsten and Weber, Marcus and Nehls, Irene}, title = {Conformational Analysis of Alternariol on the Quantum Level}, volume = {19}, journal = {J. Mol. Model.}, number = {6}, doi = {10.1007/s00894-013-1803-2}, pages = {2567 -- 2572}, year = {2013}, language = {en} } @misc{BujotzekSchuettNielsenetal.2013, author = {Bujotzek, Alexander and Sch{\"u}tt, Ole and Nielsen, Adam and Fackeldey, Konstantin and Weber, Marcus}, title = {Efficient Conformational Analysis by Partition-of-Unity Coupling}, journal = {Math Chem}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-42570}, year = {2013}, abstract = {Obtaining a sufficient sampling of conformational space is a common problem in molecular simulation. We present the implementation of an umbrella-like adaptive sampling approach based on function-based meshless discretization of conformational space that is compatible with state of the art molecular dynamics code and that integrates an eigenvector-based clustering approach for conformational analysis and the computation of inter-conformational transition rates. The approach is applied to three example systems, namely n-pentane, alanine dipeptide, and a small synthetic host-guest system, the latter two including explicitly modeled solvent.}, language = {en} } @incollection{DurmazFackeldeyWeber2011, author = {Durmaz, Vedat and Fackeldey, Konstantin and Weber, Marcus}, title = {A rapidly Mixing Monte Carlo Method for the Simulation of Slow Molecular Processes}, booktitle = {Applications of Monte Carlo Methods in Biology, Medicine and Other Fields of Science}, editor = {Mode, Ch.}, publisher = {InTech}, year = {2011}, language = {en} } @inproceedings{FackeldeyKrauseKrause2008, author = {Fackeldey, Konstantin and Krause, Dorian and Krause, Rolf}, title = {Quadrature and Implementation of the Weak Coupling Method}, booktitle = {MMM 2008: Proceedings of the 4th International Conference on Multiscale Materials Modelling, 27-31 October 2008, Tallahassee, Fla.}, pages = {62 -- 65}, year = {2008}, language = {en} } @inproceedings{FackeldeyKrauseKrause2009, author = {Fackeldey, Konstantin and Krause, Dorian and Krause, Rolf}, title = {Weak Coupling Algorithms in Multiscale Simulations}, booktitle = {Proceedings of the third Inter. Confer. on Comp. Meth. for Coupled Problems in Science and Engrg}, pages = {1023pp}, year = {2009}, language = {en} } @inproceedings{FackeldeyKrauseKrause2010, author = {Fackeldey, Konstantin and Krause, Dorian and Krause, Rolf}, title = {Numerical Validation of Constraints Based Multiscale Methods}, booktitle = {Lecture Notes in Computational Science and Engineering}, editor = {Griebel, Michael and Schweitzer, M.}, pages = {141 -- 154}, year = {2010}, language = {en} } @article{FackeldeyKrauseKrauseetal.2011, author = {Fackeldey, Konstantin and Krause, Dorian and Krause, Rolf and Lenzen, Christoph}, title = {Coupling Molecular Dynamics and Continua with Weak Constraints}, volume = {9}, journal = {SIAM Multiscale Model. Simul.}, doi = {/10.1137/100782097}, pages = {1459 -- 1494}, year = {2011}, language = {en} } @inproceedings{FackeldeyKrauseSchweitzer2008, author = {Fackeldey, Konstantin and Krause, Rolf and Schweitzer, M.}, title = {Stability of Energy Transfer in the Weak Coupling Method}, volume = {65}, booktitle = {Meshfree Methods for Partial Differential Equations IV}, editor = {Griebel, Michael and Schweitzer, M.}, pages = {111 -- 121}, year = {2008}, language = {en} } @article{Fackeldey2011, author = {Fackeldey, Konstantin}, title = {Coupling Meshbased and Meshfree Methods by a Transfer Operator Approach}, volume = {69}, journal = {IJPAM}, number = {2}, pages = {185 -- 202}, year = {2011}, language = {en} } @article{FackeldeyKrause2009, author = {Fackeldey, Konstantin and Krause, Rolf}, title = {Multiscale Coupling in Function Space - Weak Coupling between Molecular Dynamics and Continuum Mechanics}, volume = {79}, journal = {International Journal for Numerical Methods in Engineering}, number = {12}, doi = {10.1002/nme.2626}, pages = {1517 -- 1535}, year = {2009}, language = {en} } @phdthesis{Fackeldey2009, author = {Fackeldey, Konstantin}, title = {The Weak Coupling Method for Coupling Continuum Mechanics with Molecular Dynamics}, year = {2009}, language = {en} } @inproceedings{FackeldeyKrause2007, author = {Fackeldey, Konstantin and Krause, Rolf}, title = {Solving Frictional Contact Problems with Multigrid Efficiency}, volume = {50}, booktitle = {Proc.of the 16th International Conference on Domain Decomposition Methods}, editor = {Widlund, Olof}, pages = {547 -- 554}, year = {2007}, language = {en} } @inproceedings{FackeldeyKrause2007, author = {Fackeldey, Konstantin and Krause, Rolf}, title = {Weak coupling in function space}, volume = {7}, booktitle = {Proc. Appl. Math. Mech.}, number = {1}, pages = {2020113pp}, year = {2007}, language = {en} } @inproceedings{FackeldeyKrause2008, author = {Fackeldey, Konstantin and Krause, Rolf}, title = {CM/MD Coupling - A Function Space Oriented Multiscale-Coupling Approach}, volume = {8}, booktitle = {Proc. Appl. Math. Mech.}, number = {1}, pages = {10495pp}, year = {2008}, language = {en} } @inproceedings{FackeldeyBujotzekWeber2012, author = {Fackeldey, Konstantin and Bujotzek, Alexander and Weber, Marcus}, title = {A meshless discretization method for Markov state models applied to explicit water peptide folding simulations}, volume = {89}, booktitle = {Meshfree Methods for Partial Differential Equations VI}, publisher = {Springer}, pages = {141 -- 154}, year = {2012}, language = {en} } @inproceedings{Fackeldey2012, author = {Fackeldey, Konstantin}, title = {Multiscale Methods in Time and Space}, volume = {17}, booktitle = {Progress in Industrial Mathematics at ECMI 2010}, pages = {619 -- 626}, year = {2012}, language = {en} } @article{FackeldeyKlimmWeber2012, author = {Fackeldey, Konstantin and Klimm, Martina and Weber, Marcus}, title = {A Coarse Graining Method for the Dimension Reduction of the State Space of Biomolecules}, volume = {5}, journal = {Journal of Mathematical Chemistry}, number = {9}, pages = {2623 -- 2635}, year = {2012}, language = {en} } @article{KrauseFackeldeyKrause2014, author = {Krause, Dorian and Fackeldey, Konstantin and Krause, Rolf}, title = {A Parallel Multiscale Simulation Toolbox for Coupling Molecular Dynamics and Finite Elements}, journal = {Singular Phenomena and Scaling in Mathematical Models}, editor = {Griebel, Michael}, publisher = {Springer International Publishing}, doi = {10.1007/978-3-319-00786-1_14}, pages = {327 -- 346}, year = {2014}, language = {en} } @article{HaackFackeldeyRoeblitzetal.2013, author = {Haack, Fiete and Fackeldey, Konstantin and R{\"o}blitz, Susanna and Scharkoi, Olga and Weber, Marcus and Schmidt, Burkhard}, title = {Adaptive spectral clustering with application to tripeptide conformation analysis}, volume = {139}, journal = {The Journal of Chemical Physics}, doi = {10.1063/1.4830409}, pages = {110 -- 194}, year = {2013}, language = {en} } @article{ReuterWeberFackeldeyetal.2018, author = {Reuter, Bernhard and Weber, Marcus and Fackeldey, Konstantin and R{\"o}blitz, Susanna and Garcia, Martin E.}, title = {Generalized Markov State Modeling Method for Nonequilibrium Biomolecular Dynamics: Exemplified on Amyloid β Conformational Dynamics Driven by an Oscillating Electric Field}, volume = {14}, journal = {Journal of Chemical Theory and Computation}, number = {7}, doi = {10.1021/acs.jctc.8b00079}, pages = {3579 -- 3594}, year = {2018}, abstract = {Markov state models (MSMs) have received an unabated increase in popularity in recent years, as they are very well suited for the identification and analysis of metastable states and related kinetics. However, the state-of-the-art Markov state modeling methods and tools enforce the fulfillment of a detailed balance condition, restricting their applicability to equilibrium MSMs. To date, they are unsuitable to deal with general dominant data structures including cyclic processes, which are essentially associated with nonequilibrium systems. To overcome this limitation, we developed a generalization of the common robust Perron Cluster Cluster Analysis (PCCA+) method, termed generalized PCCA (G-PCCA). This method handles equilibrium and nonequilibrium simulation data, utilizing Schur vectors instead of eigenvectors. G-PCCA is not limited to the detection of metastable states but enables the identification of dominant structures in a general sense, unraveling cyclic processes. This is exemplified by application of G-PCCA on nonequilibrium molecular dynamics data of the Amyloid β (1-40) peptide, periodically driven by an oscillating electric field.}, language = {en} } @misc{WeberFackeldey2015, author = {Weber, Marcus and Fackeldey, Konstantin}, title = {G-PCCA: Spectral Clustering for Non-reversible Markov Chains}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55505}, year = {2015}, abstract = {Spectral clustering methods are based on solving eigenvalue problems for the identification of clusters, e.g., the identification of metastable subsets of a Markov chain. Usually, real-valued eigenvectors are mandatory for this type of algorithms. The Perron Cluster Analysis (PCCA+) is a well-known spectral clustering method of Markov chains. It is applicable for reversible Markov chains, because reversibility implies a real-valued spectrum. We extend this spectral clustering method also to non-reversible Markov chains and give some illustrative examples. The main idea is to replace the eigenvalue problem by a real-valued Schur decomposition. By this extension, non-reversible Markov chains can be analyzed. Furthermore, the chains need not have a positive stationary distribution. And additionally to metastabilities, dominant cycles and sinks can be identified, too.}, language = {en} } @article{AndraeMerkelDurmazetal.2014, author = {Andrae, Karsten and Merkel, Stefan and Durmaz, Vedat and Fackeldey, Konstantin and K{\"o}ppen, Robert and Weber, Marcus and Koch, Matthias}, title = {Investigation of the Ergopeptide Epimerization Process}, volume = {2}, journal = {Computation}, number = {3}, doi = {10.3390/computation2030102}, pages = {102 -- 111}, year = {2014}, abstract = {Ergopeptides, like ergocornine and a-ergocryptine, exist in an S- and in an R-configuration. Kinetic experiments imply that certain configurations are preferred depending on the solvent. The experimental methods are explained in this article. Furthermore, computational methods are used to understand this configurational preference. Standard quantum chemical methods can predict the favored configurations by using minimum energy calculations on the potential energy landscape. However, the explicit role of the solvent is not revealed by this type of methods. In order to better understand its influence, classical mechanical molecular simulations are applied. It appears from our research that "folding" the ergopeptide molecules into an intermediate state (between the S- and the R-configuration) is mechanically hindered for the preferred configurations.}, language = {en} } @article{WeberFackeldey2014, author = {Weber, Marcus and Fackeldey, Konstantin}, title = {Computing the Minimal Rebinding Effect Included in a Given Kinetics}, volume = {12}, journal = {Multiscale Model. Simul.}, number = {1}, doi = {10.1137/13091124X}, pages = {318 -- 334}, year = {2014}, abstract = {The rebinding effect is a phenomenon which occurs when observing a ligand-receptor binding process. On the macro scale this process comprises the Markov property. This Makovian view is spoiled when switching to the atomistic scale of a binding process. We therefore suggest a model which accurately describes the rebinding effect on the atomistic scale by allowing ''intermediate'' bound states. This allows us to define an indicator for the magnitude of rebinding and to formulate an optimization problem. The results form our examples show good agreement with data form laboratory.}, 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{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} } @inproceedings{FackeldeyBujotzek2014, author = {Fackeldey, Konstantin and Bujotzek, Alexander}, title = {Local Quantum-Like Updates in Classical Molecular Simulation Realized Within an Uncoupling-Coupling Approach}, volume = {19}, booktitle = {Progress in Industrial Mathematics at ECMI 2012}, doi = {10.1007/978-3-319-05365-3_42}, pages = {309 -- 313}, year = {2014}, abstract = {In this article a method to improve the precision of the classical molecular dynamics force field by solving an approximation problem with scattered quantum mechanical data is presented. This novel technique is based on two steps. In the first step a partition of unity scheme is used for partitioning the state space by meshfree basis functions. As a consequence the potential can be localized for each basis function. In a second step, for one state in each meshfree basis function, the precise QM-based charges are computed. These local QM-based charges are then used, to optimize the local potential function. The performance of this method is shown for the alanine tripeptide.}, 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} } @article{WeberFackeldey2014, author = {Weber, Marcus and Fackeldey, Konstantin}, title = {Local Refinements in Classical Molecular Dynamics Simulations}, volume = {490}, journal = {J. Phys. Conf. Ser.}, pages = {012016}, year = {2014}, language = {en} } @article{AndraeDurmazFackeldeyetal.2013, author = {Andrae, Karsten and Durmaz, Vedat and Fackeldey, Konstantin and Scharkoi, Olga and Weber, Marcus}, title = {Medizin aus dem Computer}, volume = {62}, journal = {Der Anaesthesist}, number = {7}, publisher = {Springer}, doi = {10.1007/s00101-013-2202-x}, pages = {561 -- 557}, year = {2013}, language = {de} } @misc{Fackeldey2010, author = {Fackeldey, Konstantin}, title = {Challenges in Atomistic-to-Continuum Coupling}, volume = {2010}, journal = {ZIB Report}, edition = {10-12}, doi = {10.1155/2015/834517}, year = {2010}, abstract = {This paper is concerned with the design, analysis, and implementation of concurrent coupling approaches where different (atomic and continuous) models are used simultaneously within a single simulation process. Thereby, several problems or pitfalls can happen, for example, the reflection of molecular movements at the "boundary" between the atomic and continuum regions which leads to an unphysical increase in energy in the atomic model. We investigate the problems with the aim of giving an introduction into this field and preventing errors for scientists starting their research towards multiscale methods.}, language = {en} } @phdthesis{Fackeldey2015, author = {Fackeldey, Konstantin}, title = {Crossing the Scales in Structural Mechanics and Molecular Research}, year = {2015}, language = {en} } @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} } @misc{WitzigBeckenbachEifleretal.2016, author = {Witzig, Jakob and Beckenbach, Isabel and Eifler, Leon and Fackeldey, Konstantin and Gleixner, Ambros and Grever, Andreas and Weber, Marcus}, title = {Mixed-Integer Programming for Cycle Detection in Non-reversible Markov Processes}, issn = {1438-0064}, doi = {10.1137/16M1091162}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-60353}, year = {2016}, abstract = {In this paper, we present a new, optimization-based method to exhibit cyclic behavior in non-reversible stochastic processes. While our method is general, it is strongly motivated by discrete simulations of ordinary differential equations representing non-reversible biological processes, in particular molecular simulations. Here, the discrete time steps of the simulation are often very small compared to the time scale of interest, i.e., of the whole process. In this setting, the detection of a global cyclic behavior of the process becomes difficult because transitions between individual states may appear almost reversible on the small time scale of the simulation. We address this difficulty using a mixed-integer programming model that allows us to compute a cycle of clusters with maximum net flow, i.e., large forward and small backward probability. For a synthetic genetic regulatory network consisting of a ring-oscillator with three genes, we show that this approach can detect the most productive overall cycle, outperforming classical spectral analysis methods. Our method applies to general non-equilibrium steady state systems such as catalytic reactions, for which the objective value computes the effectiveness of the catalyst.}, language = {en} } @article{WitzigBeckenbachEifleretal.2018, author = {Witzig, Jakob and Beckenbach, Isabel and Eifler, Leon and Fackeldey, Konstantin and Gleixner, Ambros and Grever, Andreas and Weber, Marcus}, title = {Mixed-Integer Programming for Cycle Detection in Non-reversible Markov Processes}, volume = {16}, journal = {Multiscale Modeling and Simulation}, number = {1}, issn = {1438-0064}, doi = {10.1137/16M1091162}, pages = {248 -- 265}, year = {2018}, abstract = {In this paper, we present a new, optimization-based method to exhibit cyclic behavior in non-reversible stochastic processes. While our method is general, it is strongly motivated by discrete simulations of ordinary differential equations representing non-reversible biological processes, in particular molecular simulations. Here, the discrete time steps of the simulation are often very small compared to the time scale of interest, i.e., of the whole process. In this setting, the detection of a global cyclic behavior of the process becomes difficult because transitions between individual states may appear almost reversible on the small time scale of the simulation. We address this difficulty using a mixed-integer programming model that allows us to compute a cycle of clusters with maximum net flow, i.e., large forward and small backward probability. For a synthetic genetic regulatory network consisting of a ring-oscillator with three genes, we show that this approach can detect the most productive overall cycle, outperforming classical spectral analysis methods. Our method applies to general non-equilibrium steady state systems such as catalytic reactions, for which the objective value computes the effectiveness of the catalyst.}, language = {en} }