TY - JOUR A1 - Straube, Arthur A1 - Winkelmann, Stefanie A1 - Höfling, Felix T1 - Accurate reduced models for the pH oscillations in the urea-urease reaction confined to giant lipid vesicles JF - The Journal of Physical Chemistry B N2 - This theoretical study concerns a pH oscillator based on the urea-urease reaction confined to giant lipid vesicles. Under suitable conditions, differential transport of urea and hydrogen ion across the unilamellar vesicle membrane periodically resets the pH clock that switches the system from acid to basic, resulting in self-sustained oscillations. We analyse the structure of the phase flow and of the limit cycle, which controls the dynamics for giant vesicles and dominates the pronouncedly stochastic oscillations in small vesicles of submicrometer size. To this end, we derive reduced models, which are amenable to analytic treatments that are complemented by numerical solutions, and obtain the period and amplitude of the oscillations as well as the parameter domain, where oscillatory behavior persists. We show that the accuracy of these predictions is highly sensitive to the employed reduction scheme. In particular, we suggest an accurate two-variable model and show its equivalence to a three-variable model that admits an interpretation in terms of a chemical reaction network. The faithful modeling of a single pH oscillator appears crucial for rationalizing experiments and understanding communication of vesicles and synchronization of rhythms. Y1 - 2023 U6 - https://doi.org/10.1021/acs.jpcb.2c09092 VL - 127 IS - 13 SP - 2955 EP - 2967 ER - TY - JOUR A1 - Sikorski, Alexander A1 - Ribera Borrell, Enric A1 - Weber, Marcus T1 - Learning Koopman eigenfunctions of stochastic diffusions with optimal importance sampling and ISOKANN JF - Journal of Mathematical Physics N2 - The dominant eigenfunctions of the Koopman operator characterize the metastabilities and slow-timescale dynamics of stochastic diffusion processes. In the context of molecular dynamics and Markov state modeling, they allow for a description of the location and frequencies of rare transitions, which are hard to obtain by direct simulation alone. In this article, we reformulate the eigenproblem in terms of the ISOKANN framework, an iterative algorithm that learns the eigenfunctions by alternating between short burst simulations and a mixture of machine learning and classical numerics, which naturally leads to a proof of convergence. We furthermore show how the intermediate iterates can be used to reduce the sampling variance by importance sampling and optimal control (enhanced sampling), as well as to select locations for further training (adaptive sampling). We demonstrate the usage of our proposed method in experiments, increasing the approximation accuracy by several orders of magnitude. Y1 - 2024 U6 - https://doi.org/10.1063/5.0140764 VL - 65 SP - 013502 ER - TY - JOUR A1 - Huynh, Ngoc A1 - Chegini, Fatemeh A1 - Pavarino, Luca A1 - Weiser, Martin A1 - Scacchi, Simone T1 - Convergence analysis of BDDC preconditioners for hybrid DG discretizations of the cardiac cell-by-cell model JF - SIAM Journal on Scientific Computing N2 - A Balancing Domain Decomposition by Constraints (BDDC) preconditioner is constructed and analyzed for the solution of hybrid Discontinuous Galerkin discretizations of reaction-diffusion systems of ordinary and partial differential equations arising in cardiac cell-by-cell models. The latter are different from the classical Bidomain and Monodomain cardiac models based on homogenized descriptions of the cardiac tissue at the macroscopic level, and therefore they allow the representation of individual cardiac cells, cell aggregates, damaged tissues and nonuniform distributions of ion channels on the cell membrane. The resulting discrete cell-by-cell models have discontinuous global solutions across the cell boundaries, hence the proposed BDDC preconditioner is based on appropriate dual and primal spaces with additional constraints which transfer information between cells (subdomains) without influencing the overall discontinuity of the global solution. A scalable convergence rate bound is proved for the resulting BDDC cell-by-cell preconditioned operator, while numerical tests validate this bound and investigate its dependence on the discretization parameters. Y1 - 2023 VL - 45 IS - 6 SP - A2836 EP - A2857 ER - TY - JOUR A1 - Boltz, Horst-Holger A1 - Kurchan, Jorge A1 - Liu, Andrea J. T1 - Fluctuation distributions of energy minima in complex landscapes JF - Physical Review Research N2 - We discuss the properties of the distributions of energies of minima obtained by gradient descent in complex energy landscapes. We find strikingly similar phenomenology across several prototypical models. We particularly focus on the distribution of energies of minima in the analytically well-understood p-spin-interaction spin-glass model. We numerically find non-Gaussian distributions that resemble the Tracy-Widom distributions often found in problems of random correlated variables, and nontrivial finite-size scaling. Based on this, we propose a picture of gradient-descent dynamics that highlights the importance of a first-passage process in the eigenvalues of the Hessian. This picture provides a concrete link to problems in which the Tracy-Widom distribution is established. Aspects of this first-passage view of gradient-descent dynamics are generic for nonconvex complex landscapes, rationalizing the commonality that we find across models. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevResearch.3.013061 VL - 3 IS - 1 SP - 013061 ER - TY - JOUR A1 - Carderera, Alejandro A1 - Pokutta, Sebastian A1 - Schütte, Christof A1 - Weiser, Martin T1 - An efficient first-order conditional gradient algorithm in data-driven sparse identification of nonlinear dynamics to solve sparse recovery problems under noise JF - Journal of Computational and Applied Mathematics N2 - Governing equations are essential to the study of nonlinear dynamics, often enabling the prediction of previously unseen behaviors as well as the inclusion into control strategies. The discovery of governing equations from data thus has the potential to transform data-rich fields where well-established dynamical models remain unknown. This work contributes to the recent trend in data-driven sparse identification of nonlinear dynamics of finding the best sparse fit to observational data in a large library of potential nonlinear models. We propose an efficient first-order Conditional Gradient algorithm for solving the underlying optimization problem. In comparison to the most prominent alternative algorithms, the new algorithm shows significantly improved performance on several essential issues like sparsity-induction, structure-preservation, noise robustness, and sample efficiency. We demonstrate these advantages on several dynamics from the field of synchronization, particle dynamics, and enzyme chemistry. Y1 - 2021 ER - TY - JOUR A1 - Wulkow, Niklas A1 - Telgmann, Regina A1 - Hungenberg, Klaus-Dieter A1 - Schütte, Christof A1 - Wulkow, Michael T1 - Deterministic and Stochastic Parameter Estimation for Polymer Reaction Kinetics I: Theory and Simple Examples JF - Macromolecular Theory and Simulations N2 - Two different approaches to parameter estimation (PE) in the context of polymerization are introduced, refined, combined, and applied. The first is classical PE where one is interested in finding parameters which minimize the distance between the output of a chemical model and experimental data. The second is Bayesian PE allowing for quantifying parameter uncertainty caused by experimental measurement error and model imperfection. Based on detailed descriptions of motivation, theoretical background, and methodological aspects for both approaches, their relation are outlined. The main aim of this article is to show how the two approaches complement each other and can be used together to generate strong information gain regarding the model and its parameters. Both approaches and their interplay in application to polymerization reaction systems are illustrated. This is the first part in a two-article series on parameter estimation for polymer reaction kinetics with a focus on theory and methodology while in the second part a more complex example will be considered. Y1 - 2021 U6 - https://doi.org/10.1002/mats.202100017 VL - 30 ER - TY - JOUR A1 - del Razo, Mauricio J. A1 - Dibak, Manuel A1 - Schütte, Christof A1 - Noé, Frank T1 - Multiscale molecular kinetics by coupling Markov state models and reaction-diffusion dynamics JF - The Journal of Chemical Physics Y1 - 2021 U6 - https://doi.org/10.1063/5.0060314 VL - 155 IS - 12 ER - TY - JOUR A1 - Gelss, Patrick A1 - Klus, Stefan A1 - Schuster, Ingmar A1 - Schütte, Christof T1 - Feature space approximation for kernel-based supervised learning JF - Knowledge-Based Sytems Y1 - 2021 U6 - https://doi.org/https://doi.org/10.1016/j.knosys.2021.106935 VL - 221 PB - Elsevier ER - TY - JOUR A1 - Jäger, Klaus A1 - Sutter, Johannes A1 - Hammerschmidt, Martin A1 - Schneider, Philipp-Immanuel A1 - Becker, Christiane T1 - Prospects of light management in perovskite/silicon tandem solar cells JF - Nanophotonics Y1 - 2021 U6 - https://doi.org/10.1515/nanoph-2020-0674 VL - 10 SP - 1991 ER - TY - JOUR A1 - Liu, Yongtao A1 - Zhou, Zhiguang A1 - Wang, Fan A1 - Kewes, Günter A1 - Wen, Shihui A1 - Burger, Sven A1 - Wakiani, Majid Ebrahimi A1 - Xi, Peng A1 - Yang, Jiong A1 - Yang, Xusan A1 - Benson, Oliver A1 - Jin, Dayong T1 - Axial localization and tracking of self-interference nanoparticles by lateral point spread functions JF - Nat. Commun. Y1 - 2021 U6 - https://doi.org/10.1038/s41467-021-22283-0 VL - 12 SP - 2019 ER -