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 - TY - JOUR A1 - Bittracher, Andreas A1 - Moschner, Johann A1 - Koksch, Beate A1 - Netz, Roland A1 - Schütte, Christof T1 - Exploring the locking stage of NFGAILS amyloid fibrillation via transition manifold analysis JF - The European Physical Journal B Y1 - 2021 U6 - https://doi.org/10.1140/epjb/s10051-021-00200-0 VL - 94 ER - TY - JOUR A1 - Lücke, Marvin A1 - Lu, Jingyi A1 - Quevedo, Daniel E. T1 - Coding for secrecy in remote state estimation with an adversary JF - IEEE Transactions on Automatic Control N2 - We study the problem of remote state estimation in the presence of an eavesdropper. A sensor transmits state information over a packet-dropping link to a legitimate user. This information is randomly overheard by an eavesdropper. To reduce information leakage to the eavesdropper, previous studies have shown that by encoding the estimate with the acknowledgments (Acks), perfect secrecy can be achieved. However, this strategy greatly relies on the accuracy of the Acks and may easily fail if the Acks are compromised by cyberattacks. In this article, we tackle this issue by proposing to switch between sending an encoded state and sending the plain state to stay resilient against fake Acks. Specifically, we assume the Acks to be randomly attacked and derive recursive expressions for the minimum-mean-squared error estimates and error covariance matrices at the legitimate user and at the eavesdropper. Based upon this, we propose a transmission policy that depends on the probability of synchronization. We formulate a partially observable Markov decision process to model the evolution of the synchronization status and derive associated optimal transmission policies. Numerical examples are provided to verify the theoretical results. Y1 - 2022 U6 - https://doi.org/10.1109/TAC.2022.3169839 VL - 67 IS - 9 SP - 4955 EP - 4962 ER - TY - JOUR A1 - Chegini, Fatemeh A1 - Kopanicakova, Alena A1 - Krause, Rolf A1 - Weiser, Martin T1 - Efficient Identification of Scars using Heterogeneous Model Hierarchies JF - EP Europace N2 - Aims. Detection and quantification of myocardial scars are helpful both for diagnosis of heart diseases and for building personalized simulation models. Scar tissue is generally charac­terized by a different conduction of electrical excitation. We aim at estimating conductivity-related parameters from endocardial mapping data, in particular the conductivity tensor. Solving this inverse problem requires computationally expensive monodomain simulations on fine discretizations. Therefore, we aim at accelerating the estimation using a multilevel method combining electrophysiology models of different complexity, namely the mono­domain and the eikonal model. Methods. Distributed parameter estimation is performed by minimizing the misfit between simulated and measured electrical activity on the endocardial surface, subject to the mono­domain model and regularization, leading to a constrained optimization problem. We formulate this optimization problem, including the modeling of scar tissue and different regularizations, and design an efficient iterative solver. We consider monodomain grid hierarchies and monodomain-eikonal model hierarchies in a recursive multilevel trust-region method. Results. From several numerical examples, both the efficiency of the method and the estimation quality, depending on the data, are investigated. The multilevel solver is significantly faster than a comparable single level solver. Endocardial mapping data of realistic density appears to be just sufficient to provide quantitatively reasonable estimates of location, size, and shape of scars close to the endocardial surface. Conclusion. In several situations, scar reconstruction based on eikonal and monodomain models differ significantly, suggesting the use of the more accurate but more expensive monodomain model for this purpose. Still, eikonal models can be utilized to accelerate the computations considerably, enabling the use of complex electrophysiology models for estimating myocardial scars from endocardial mapping data. Y1 - 2021 U6 - https://doi.org/10.1093/europace/euaa402 VL - 23 SP - i113 EP - i122 ER - TY - CHAP A1 - Jäger, Klaus A1 - Tillmann, Peter A1 - Katz, Eugene A. A1 - Becker, Christiane T1 - Simulating bifacial perovskite/silicon tandem solar cells in large PV fields T2 - OSA Advanced Photonics Congress Y1 - 2020 U6 - https://doi.org/10.1364/PVLED.2020.PvTh1G.3 SP - PvTh1G.3 ER - TY - JOUR A1 - Jäger, Klaus A1 - Tillmann, Peter A1 - Katz, Eugene A. A1 - Becker, Christiane T1 - Perovskite/silicon tandem solar cells: Effect of luminescent coupling and bifaciality JF - Sol. RRL Y1 - 2021 U6 - https://doi.org/10.1002/solr.202000628 VL - 5 SP - 2000628 ER - TY - JOUR A1 - Helfmann, Luzie A1 - Heitzig, Jobst A1 - Koltai, Péter A1 - Kurths, Jürgen A1 - Schütte, Christof T1 - Statistical analysis of tipping pathways in agent-based models JF - Eur. Phys. J. Spec. Top. N2 - Agent-based models are a natural choice for modeling complex social systems. In such models simple stochastic interaction rules for a large population of individuals on the microscopic scale can lead to emergent dynamics on the macroscopic scale, for instance a sudden shift of majority opinion or behavior. Here we are introducing a methodology for studying noise-induced tipping between relevant subsets of the agent state space representing characteristic configurations. Due to a large number of interacting individuals, agent-based models are high-dimensional, though usually a lower-dimensional structure of the emerging collective behaviour exists. We therefore apply Diffusion Maps, a non-linear dimension reduction technique, to reveal the intrinsic low-dimensional structure. We characterize the tipping behaviour by means of Transition Path Theory, which helps gaining a statistical understanding of the tipping paths such as their distribution, flux and rate. By systematically studying two agent-based models that exhibit a multitude of tipping pathways and cascading effects, we illustrate the practicability of our approach. Y1 - 2021 U6 - https://doi.org/10.1140/epjs/s11734-021-00191-0 VL - 230 SP - 3249 EP - 3271 ER - TY - JOUR A1 - Manley, Phillip A1 - Segantini, Michele A1 - Ahiboz, Doguscan A1 - Hammerschmidt, Martin A1 - Arnaoutakis, Georgios A1 - MacQueen, Rowan W. A1 - Burger, Sven A1 - Becker, Christiane T1 - Double-layer metasurface for enhanced photon up-conversion JF - APL Photon. Y1 - 2021 U6 - https://doi.org/10.1063/5.0040839 VL - 6 SP - 036103 ER - TY - JOUR A1 - Pakhomov, Anton V. A1 - Hammerschmidt, Martin A1 - Burger, Sven A1 - Pertsch, Thomas A1 - Setzpfand, Frank T1 - Modeling of surface-induced second-harmonic generation from multilayer structures by the transfer matrix method JF - Opt. Express Y1 - 2021 U6 - https://doi.org/10.1364/OE.417066 VL - 29 SP - 9098 ER - TY - JOUR A1 - Binkowski, Felix A1 - Wu, Tong A1 - Lalanne, Philippe A1 - Burger, Sven A1 - Govorov, Alexander O. T1 - Hot Electron Generation through Near-Field Excitation of Plasmonic Nanoresonators JF - ACS Photonics Y1 - 2021 U6 - https://doi.org/10.1021/acsphotonics.1c00231 VL - 8 SP - 1243 ER - TY - JOUR A1 - Cvetkovic, Nada A1 - Conrad, Tim A1 - Lie, Han Cheng T1 - A Convergent Discretisation Method for Transition Path Theory for Diffusion Processes JF - Multiscale Modeling & Simulation Y1 - 2021 U6 - https://doi.org/10.1137/20M1329354 VL - 19 IS - 1 SP - 242 EP - 266 PB - Society for Industrial and Applied Mathematics ER - TY - BOOK A1 - Deuflhard, Peter A1 - Weiser, Martin T1 - Numerische Mathematik 3. Adaptive Lösung partieller Differentialgleichungen Y1 - 2020 SN - 978-3-11-069168-9 U6 - https://doi.org/10.1515/9783110689655 PB - de Gruyter ET - 2 ER - TY - JOUR A1 - Sikorski, Alexander A1 - Weber, Marcus A1 - Schütte, Christof T1 - The Augmented Jump Chain JF - Advanced Theory and Simulations N2 - Modern methods of simulating molecular systems are based on the mathematical theory of Markov operators with a focus on autonomous equilibrated systems. However, non-autonomous physical systems or non-autonomous simulation processes are becoming more and more important. A representation of non-autonomous Markov jump processes is presented as autonomous Markov chains on space-time. Augmenting the spatial information of the embedded Markov chain by the temporal information of the associated jump times, the so-called augmented jump chain is derived. The augmented jump chain inherits the sparseness of the infinitesimal generator of the original process and therefore provides a useful tool for studying time-dependent dynamics even in high dimensions. Furthermore, possible generalizations and applications to the computation of committor functions and coherent sets in the non-autonomous setting are discussed. After deriving the theoretical foundations, the concepts with a proof-of-concept Galerkin discretization of the transfer operator of the augmented jump chain applied to simple examples are illustrated. Y1 - 2021 U6 - https://doi.org/10.1002/adts.202000274 VL - 4 IS - 4 PB - Wiley-VCH ER - TY - JOUR A1 - Sechi, Renata A1 - Sikorski, Alexander A1 - Weber, Marcus T1 - Estimation of the Koopman Generator by Newton's Extrapolation JF - Multiscale Modeling and Simulation N2 - This article addresses the problem of estimating the Koopman generator of a Markov process. The direct computation of the infinitesimal generator is not easy because of the discretization of the state space, in particular because of the trade-off inherent in the choice of the best lag time to study the process. Short lag times implies a strong discretization of the state space and a consequent loss of Markovianity. Large lag times bypass events on fast timescales. We propose a method to approximate the generator with the computation of the Newton polynomial extrapolation. This technique is a multistep approach which uses as its input Koopman transfer operators evaluated for a series of lag times. Thus, the estimated infinitesimal generator combines information from different time resolutions and does not bias only fast- or slow-decaying dynamics. We show that the multi-scale Newton method can improve the estimation of the generator in comparison to the computation using finite difference or matrix logarithm methods. Y1 - 2020 U6 - https://doi.org/10.1137/20M1333006 VL - 19 IS - 2 SP - 758 EP - 774 PB - SIAM ER - TY - THES A1 - Venkatareddy, Narendra Lagumaddepalli T1 - Revealing secrets of mussel-glue mimetic peptides - From advanced NMR to computational process modelling Y1 - 2020 ER - TY - THES A1 - Eshtewy, Neveen Ali Salem T1 - Mathematical Modeling of Metabolic-Genetic Networks Y1 - 2020 ER - TY - JOUR A1 - Öttinger, Hans Christian A1 - Montefusco, Alberto A1 - Peletier, Mark A. T1 - A Framework of Nonequilibrium Statistical Mechanics. I. Role and Types of Fluctuations JF - Journal of Non-Equilibrium Thermodynamics N2 - Understanding the fluctuations by which phenomenological evolution equations with thermodynamic structure can be enhanced is the key to a general framework of nonequilibrium statistical mechanics. These fluctuations provide an idealized representation of microscopic details. We consider fluctuation-enhanced equations associated with Markov processes and elaborate the general recipes for evaluating dynamic material properties, which characterize force-flux constitutive laws, by statistical mechanics. Markov processes with continuous trajectories are conveniently characterized by stochastic differential equations and lead to Green–Kubo-type formulas for dynamic material properties. Markov processes with discontinuous jumps include transitions over energy barriers with the rates calculated by Kramers. We describe a unified approach to Markovian fluctuations and demonstrate how the appropriate type of fluctuations (continuous versus discontinuous) is reflected in the mathematical structure of the phenomenological equations. Y1 - 2021 U6 - https://doi.org/10.1515/jnet-2020-0068 VL - 46 IS - 1 SP - 1 EP - 13 PB - De Gruyter ER - TY - JOUR A1 - Montefusco, Alberto A1 - Peletier, Mark A. A1 - Öttinger, Hans Christian T1 - A Framework of Nonequilibrium Statistical Mechanics. II. Coarse-Graining JF - Journal of Non-Equilibrium Thermodynamics N2 - For a given thermodynamic system, and a given choice of coarse-grained state variables, the knowledge of a force-flux constitutive law is the basis for any nonequilibrium modeling. In the first paper of this series we established how, by a generalization of the classical fluctuation-dissipation theorem (FDT), the structure of a constitutive law is directly related to the distribution of the fluctuations of the state variables. When these fluctuations can be expressed in terms of diffusion processes, one may use Green–Kubo-type coarse-graining schemes to find the constitutive laws. In this paper we propose a coarse-graining method that is valid when the fluctuations are described by means of general Markov processes, which include diffusions as a special case. We prove the success of the method by numerically computing the constitutive law for a simple chemical reaction A⇄B. Furthermore, we show that, for such a system, one cannot find a consistent constitutive law by any Green–Kubo-like scheme. Y1 - 2021 U6 - https://doi.org/10.1515/jnet-2020-0069 VL - 46 IS - 1 SP - 15 EP - 33 PB - De Gruyter ER - TY - JOUR A1 - Sunkara, Vikram A1 - Heinz, Gitta A. A1 - Heinrich, Frederik F. A1 - Durek, Pawel A1 - Mobasheri, Ali A1 - Mashreghi, Mir-Farzin A1 - Lang, Annemarie T1 - Combining segmental bulk- and single-cell RNA-sequencing to define the chondrocyte gene expression signature in the murine knee joint JF - bioarxiv (Accepted in Osteoarthr. Cartil.) Y1 - 2020 U6 - https://doi.org/10.1101/2020.06.13.148056 ER - TY - JOUR A1 - Lang, Annemarie A1 - Helfmeier, Sarah A1 - Stefanowski, Jonathan A1 - Kuppe, Aditi A1 - Sunkara, Vikram A1 - Pfeiffenberger, Moritz A1 - Wolter, Angelique A1 - Damerau, Alexandra A1 - Hemmati-Sadeghi, Shabnam A1 - Ringe, Jochen A1 - Haag, Rainer A1 - Hauser, Anja E. A1 - Löhning, Max A1 - Perka, Carsten A1 - Duda, Georg A1 - Hoff, Paula A1 - Schmidt-Bleek, Katharina A1 - Gaber, Timo A1 - Buttgereit, Frank T1 - HIF-stabilization prevents delayed fracture healing JF - bioarxiv Y1 - 2020 U6 - https://doi.org/10.1101/2020.07.02.182832 ER - TY - JOUR A1 - Laydon, Daniel J. A1 - Sunkara, Vikram A1 - Boelen, Lies A1 - Bangham, Charles R. M. A1 - Asquith, Becca T1 - The relative contributions of infectious and mitotic spread to HTLV-1 persistence JF - PLOS Computational Biology Y1 - 2020 U6 - https://doi.org/10.1371/journal.pcbi.1007470 ER - TY - JOUR A1 - Schiela, Anton A1 - Stöcklein, Matthias A1 - Weiser, Martin T1 - A primal dual projection algorithm for efficient constraint preconditioning JF - SIAM Journal on Scientific Computing N2 - We consider a linear iterative solver for large scale linearly constrained quadratic minimization problems that arise, for example, in optimization with PDEs. By a primal-dual projection (PDP) iteration, which can be interpreted and analysed as a gradient method on a quotient space, the given problem can be solved by computing sulutions for a sequence of constrained surrogate problems, projections onto the feasible subspaces, and Lagrange multiplier updates. As a major application we consider a class of optimization problems with PDEs, where PDP can be applied together with a projected cg method using a block triangular constraint preconditioner. Numerical experiments show reliable and competitive performance for an optimal control problem in elasticity. Y1 - 2021 U6 - https://doi.org/10.1137/20M1380739 VL - 43 IS - 6 SP - A4095 EP - A4120 ER - TY - JOUR A1 - Mielke, Alexander A1 - Montefusco, Alberto A1 - Peletier, Mark A. T1 - Exploring families of energy-dissipation landscapes via tilting: three types of EDP convergence JF - Continuum Mechanics and Thermodynamics N2 - This paper revolves around a subtle distinction between two concepts: passing to the limit in a family of gradient systems, on one hand, and deriving effective kinetic relations on the other. The two concepts are strongly related, and in many examples they even appear to be the same. Our main contributions are to show that they are different, to show that well-known techniques developed for the former may give incorrect results for the latter, and to introduce new tools to remedy this. The approach is based on the Energy-Dissipation Principle that provides a variational formulation to gradient-flow equations that allows one to apply techniques from Γ-convergence of functional on states and functionals on trajectories. Y1 - 2021 U6 - https://doi.org/10.1007/s00161-020-00932-x VL - 33 SP - 611 EP - 637 PB - Springer ER - TY - THES A1 - Santiago, Xavier Garcia T1 - Numerical methods for shape optimization of photonic nanostructures Y1 - 2021 U6 - https://doi.org/10.5445/IR/1000131006 ER - TY - JOUR A1 - Raharinirina, N. Alexia A1 - Acevedo-Trejos, Esteban A1 - Merico, Agostino T1 - Modelling the acclimation capacity of coral reefs to a warming ocean JF - PLOS COMPUTATIONAL BIOLOGY N2 - The symbiotic relationship between corals and photosynthetic algae is the foundation of coral reef ecosystems. This relationship breaks down, leading to coral death, when sea temperature exceeds the thermal tolerance of the coral-algae complex. While acclimation via phenotypic plasticity at the organismal level is an important mechanism for corals to cope with global warming, community-based shifts in response to acclimating capacities may give valuable indications about the future of corals at a regional scale. Reliable regional-scale predictions, however, are hampered by uncertainties on the speed with which coral communities will be able to acclimate. Here we present a trait-based, acclimation dynamics model, which we use in combination with observational data, to provide a first, crude estimate of the speed of coral acclimation at the community level and to investigate the effects of different global warming scenarios on three iconic reef ecosystems of the tropics: Great Barrier Reef, South East Asia, and Caribbean. The model predicts that coral acclimation may confer some level of protection by delaying the decline of some reefs such as the Great Barrier Reef. However, the current rates of acclimation will not be sufficient to rescue corals from global warming. Based on our estimates of coral acclimation capacities, the model results suggest substantial declines in coral abundances in all three regions, ranging from 12% to 55%, depending on the region and on the climate change scenario considered. Our results highlight the importance and urgency of precise assessments and quantitative estimates, for example through laboratory experiments, of the natural acclimation capacity of corals and of the speed with which corals may be able to acclimate to global warming. Y1 - 2022 U6 - https://doi.org/10.1371/journal.pcbi.1010099 ER - TY - GEN A1 - Schütte, Christof A1 - Klus, Stefan A1 - Hartmann, Carsten T1 - Overcoming the Timescale Barrier in Molecular Dynamics: Transfer Operators, Variational Principles, and Machine Learning N2 - One of the main challenges in molecular dynamics is overcoming the “timescale barrier”, a phrase used to describe that in many realistic molecular systems, biologically important rare transitions occur on timescales that are not accessible to direct numerical simulation, not even on the largest or specifically dedicated supercomputers. This article discusses how to circumvent the timescale barrier by a collection of transfer operator-based techniques that have emerged from dynamical systems theory, numerical mathematics, and machine learning over the last two decades. We will focus on how transfer operators can be used to approximate the dynamical behavior on long timescales, review the introduction of this approach into molecular dynamics, and outline the respective theory as well as the algorithmic development from the early numerics-based methods, via variational reformulations, to modern data-based techniques utilizing and improving concepts from machine learning. Furthermore, its relation to rare event simulation techniques will be explained, revealing a broad equivalence of variational principles for long-time quantities in MD. The article will mainly take a mathematical perspective and will leave the application to real-world molecular systems to the more than 1000 research articles already written on this subject. T3 - ZIB-Report - 22-25 Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-88637 SN - 1438-0064 ER - TY - JOUR A1 - Ramos Leite, Tiago A1 - Zschiedrich, Lin A1 - Kizilkaya, Orhan A1 - McPeak, Kevin M. T1 - Resonant Plasmonic–Biomolecular Chiral Interactions in the Far-Ultraviolet: Enantiomeric Discrimination of sub-10 nm Amino Acid Films JF - Nano Lett. Y1 - 2022 UR - https://digitalcommons.lsu.edu/chem_engineering_pubs/34/ U6 - https://doi.org/10.1021/acs.nanolett.2c01724 VL - 22 SP - 7343 ER - TY - JOUR A1 - Rej, Sourav A1 - Santiago, Eva Yazmin A1 - Baturina, Olga A1 - Zhang, Yu A1 - Burger, Sven A1 - Kment, Stepan A1 - Govorov, Alexander O. A1 - Naldoni, Alberto T1 - Colloidal Titanium Nitride Nanobars for Broadband Inexpensive Plasmonics and Photochemistry from Visible to Mid-IR Wavelengths JF - Nano Energy Y1 - 2022 U6 - https://doi.org/10.1016/j.nanoen.2022.107989 VL - 104 SP - 107989 ER - TY - CHAP A1 - Berner, Julius A1 - Richter, Lorenz A1 - Ullrich, Karen T1 - An optimal control perspective on diffusion-based generative modeling T2 - Transactions on Machine Learning Research N2 - We establish a connection between stochastic optimal control and generative models based on stochastic differential equations (SDEs) such as recently developed diffusion probabilistic models. In particular, we derive a Hamilton-Jacobi-Bellman equation that governs the evolution of the log-densities of the underlying SDE marginals. This perspective allows to transfer methods from optimal control theory to generative modeling. First, we show that the evidence lower bound is a direct consequence of the well-known verification theorem from control theory. Further, we develop a novel diffusion-based method for sampling from unnormalized densities -- a problem frequently occurring in statistics and computational sciences. Y1 - 2024 UR - https://openreview.net/forum?id=oYIjw37pTP ER - TY - JOUR A1 - Nüsken, Nikolas A1 - Richter, Lorenz T1 - Interpolating between BSDEs and PINNs: deep learning for elliptic and parabolic boundary value problems JF - Journal of Machine Learning N2 - Solving high-dimensional partial differential equations is a recurrent challenge in economics, science and engineering. In recent years, a great number of computational approaches have been developed, most of them relying on a combination of Monte Carlo sampling and deep learning based approximation. For elliptic and parabolic problems, existing methods can broadly be classified into those resting on reformulations in terms of backward stochastic differential equations (BSDEs) and those aiming to minimize a regression-type L2-error (physics-informed neural networks, PINNs). In this paper, we review the literature and suggest a methodology based on the novel diffusion loss that interpolates between BSDEs and PINNs. Our contribution opens the door towards a unified understanding of numerical approaches for high-dimensional PDEs, as well as for implementations that combine the strengths of BSDEs and PINNs. The diffusion loss furthermore bears close similarities to (least squares) temporal difference objectives found in reinforcement learning. We also discuss eigenvalue problems and perform extensive numerical studies, including calculations of the ground state for nonlinear Schr ¨odinger operators and committor functions relevant in molecular dynamics. Y1 - 2023 U6 - https://doi.org/0.4208/jml.220416 VL - 2 IS - 1 SP - 31 EP - 64 ER - TY - JOUR A1 - Hartmann, Carsten A1 - Richter, Lorenz T1 - Nonasymptotic bounds for suboptimal importance sampling JF - SIAM/ASA Journal on Uncertainty Quantification N2 - Importance sampling is a popular variance reduction method for Monte Carlo estimation, where an evident question is how to design good proposal distributions. While in most cases optimal (zero-variance) estimators are theoretically possible, in practice only suboptimal proposal distributions are available and it can often be observed numerically that those can reduce statistical performance significantly, leading to large relative errors and therefore counteracting the original intention. Previous analysis on importance sampling has often focused on asymptotic arguments that work well in a large deviations regime. In this article, we provide lower and upper bounds on the relative error in a nonasymptotic setting. They depend on the deviation of the actual proposal from optimality, and we thus identify potential robustness issues that importance sampling may have, especially in high dimensions. We particularly focus on path sampling problems for diffusion processes with nonvanishing noise, for which generating good proposals comes with additional technical challenges. We provide numerous numerical examples that support our findings and demonstrate the applicability of the derived bounds. Y1 - 2024 U6 - https://doi.org/10.1137/21M1427760 VL - 12 IS - 2 SP - 309 EP - 346 ER - TY - JOUR A1 - Rickert, Lucas A1 - Betz, Fridtjof A1 - Plock, Matthias A1 - Burger, Sven A1 - Heindel, Tobias T1 - High-performance designs for fiber-pigtailed quantum-light sources based on quantum dots in electrically-controlled circular Bragg gratings JF - Opt. Express Y1 - 2023 U6 - https://doi.org/10.1364/OE.486060 VL - 31 SP - 14750 ER - TY - GEN A1 - Rickert, Lucas A1 - Betz, Fridtjof A1 - Plock, Matthias A1 - Burger, Sven A1 - Heindel, Tobias T1 - Data publication for "High-performance designs for fiber-pigtailed quantum-light sources based on quantum dots in electrically-controlled circular Bragg gratings" T2 - Zenodo Y1 - 2022 U6 - https://doi.org/10.5281/zenodo.7360516 SP - 7360516 ER - TY - JOUR A1 - Gorgulla, Christoph A1 - Jayaraj, Abhilash A1 - Fackeldey, Konstantin A1 - Arthanari, Haribabu T1 - Emerging frontiers in virtual drug discovery: From quantum mechanical methods to deep learning approaches JF - Current Opinion in Chemical Biology N2 - Virtual screening-based approaches to discover initial hit and lead compounds have the potential to reduce both the cost and time of early drug discovery stages, as well as to find inhibitors for even challenging target sites such as protein–protein interfaces. Here in this review, we provide an overview of the progress that has been made in virtual screening methodology and technology on multiple fronts in recent years. The advent of ultra-large virtual screens, in which hundreds of millions to billions of compounds are screened, has proven to be a powerful approach to discover highly potent hit compounds. However, these developments are just the tip of the iceberg, with new technologies and methods emerging to propel the field forward. Examples include novel machine-learning approaches, which can reduce the computational costs of virtual screening dramatically, while progress in quantum-mechanical approaches can increase the accuracy of predictions of various small molecule properties. Y1 - 2022 U6 - https://doi.org/10.1016/j.cbpa.2022.102156 VL - 69 SP - 102156 EP - 102156-12 ER - TY - GEN 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 N2 - Our theoretical study concerns an urea-urease-based pH oscillator 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 limit cycle, which controls the dynamics for giant vesicles and dominates the strongly stochastic oscillations in small vesicles of submicrometer size. To this end, we derive reduced models, amenable to analytic treatments, and show that the accuracy of predictions, including the period of oscillations, is highly sensitive to the choice of the 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 accurate description of a single pH oscillator appears crucial for rationalizing experiments and understanding communication of vesicles and synchronization of rhythms. T3 - ZIB-Report - 22-21 Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-88179 SN - 1438-0064 ER - TY - JOUR A1 - Donati, Luca A1 - Weber, Marcus T1 - Assessing transition rates as functions of environmental variables JF - The Journal of Chemical Physics N2 - We present a method to estimate the transition rates of molecular systems under different environmental conditions which cause the formation or the breaking of bonds and require the sampling of the Grand Canonical Ensemble. For this purpose, we model the molecular system in terms of probable "scenarios", governed by different potential energy functions, which are separately sampled by classical MD simulations. Reweighting the canonical distribution of each scenario according to specific environmental variables, we estimate the grand canonical distribution, then we use the Square Root Approximation (SqRA) method to discretize the Fokker-Planck operator into a rate matrix and the robust Perron Cluster Cluster Analysis (PCCA+) method to coarse-grain the kinetic model. This permits to efficiently estimate the transition rates of conformational states as functions of environmental variables, for example, the local pH at a cell membrane. In this work we formalize the theoretical framework of the procedure and we present a numerical experiment comparing the results with those provided by a constant-pH method based on non-equilibrium Molecular Dynamics Monte Carlo simulations. The method is relevant for the development of new drug design strategies which take into account how the cellular environment influences biochemical processes. Y1 - 2022 U6 - https://doi.org/10.1063/5.0109555 VL - 157 IS - 22 SP - 224103-1 EP - 224103-14 PB - AIP Publishing ER - TY - JOUR A1 - Sechi, Renata A1 - Fackeldey, Konstantin A1 - Chewle, Surahit A1 - Weber, Marcus T1 - SepFree NMF: A Toolbox for Analyzing the Kinetics of Sequential Spectroscopic Data JF - Algorithms N2 - This work addresses the problem of determining the number of components from sequential spectroscopic data analyzed by non-negative matrix factorization without separability assumption (SepFree NMF). These data are stored in a matrix M of dimension “measured times” versus “measured wavenumbers” and can be decomposed to obtain the spectral fingerprints of the states and their evolution over time. SepFree NMF assumes a memoryless (Markovian) process to underline the dynamics and decomposes M so that M=WH, with W representing the components’ fingerprints and H their kinetics. However, the rank of this decomposition (i.e., the number of physical states in the process) has to be guessed from pre-existing knowledge on the observed process. We propose a measure for determining the number of components with the computation of the minimal memory effect resulting from the decomposition; by quantifying how much the obtained factorization is deviating from the Markovian property, we are able to score factorizations of a different number of components. In this way, we estimate the number of different entities which contribute to the observed system, and we can extract kinetic information without knowing the characteristic spectra of the single components. This manuscript provides the mathematical background as well as an analysis of computer generated and experimental sequentially measured Raman spectra. Y1 - 2022 U6 - https://doi.org/10.3390/a15090297 VL - 15 IS - 9 SP - 297 ER - TY - JOUR A1 - Djurdjevac Conrad, Natasa A1 - Köppl, Jonas A1 - Djurdjevac, Ana T1 - Feedback Loops in Opinion Dynamics of Agent-Based Models with Multiplicative Noise JF - Entropy N2 - We introduce an agent-based model for co-evolving opinions and social dynamics, under the influence of multiplicative noise. In this model, every agent is characterized by a position in a social space and a continuous opinion state variable. Agents’ movements are governed by the positions and opinions of other agents and similarly, the opinion dynamics are influenced by agents’ spatial proximity and their opinion similarity. Using numerical simulations and formal analyses, we study this feedback loop between opinion dynamics and the mobility of agents in a social space. We investigate the behaviour of this ABM in different regimes and explore the influence of various factors on the appearance of emerging phenomena such as group formation and opinion consensus. We study the empirical distribution, and, in the limit of infinite number of agents, we derive a corresponding reduced model given by a partial differential equation (PDE). Finally, using numerical examples, we show that a resulting PDE model is a good approximation of the original ABM. Y1 - 2022 U6 - https://doi.org/10.3390/e24101352 VL - 24(10) ER - TY - JOUR A1 - Donati, Luca A1 - Weber, Marcus A1 - Keller, Bettina G. T1 - A review of Girsanov Reweighting and of Square Root Approximation for building molecular Markov State Models JF - Journal of Mathematical Physics N2 - Dynamical reweighting methods permit to estimate kinetic observables of a stochastic process governed by a target potential U(x) from trajectories that have been generated at a different potential V(x). In this article, we present Girsanov reweighting and Square Root Approximation (SqRA): the first method reweights path probabilities exploiting the Girsanov theorem and can be applied to Markov State Models (MSMs) to reweight transition probabilities; the second method was originally developed to discretize the Fokker-Planck operator into a transition rate matrix, but here we implement it into a reweighting scheme for transition rates. We begin by reviewing the theoretical background of the methods, then present two applications relevant to Molecular Dynamics (MD), highlighting their strengths and weaknesses. Y1 - 2022 U6 - https://doi.org/10.1063/5.0127227 VL - 63 IS - 12 SP - 123306-1 EP - 123306-21 PB - AIP Publishing ER - TY - JOUR A1 - Quer, Jannes A1 - Ribera Borrell, Enric T1 - Connecting Stochastic Optimal Control and Reinforcement Learning JF - Journal of Mathematical Physics N2 - In this paper the connection between stochastic optimal control and reinforcement learning is investigated. Our main motivation is to apply importance sampling to sampling rare events which can be reformulated as an optimal control problem. By using a parameterised approach the optimal control problem becomes a stochastic optimization problem which still raises some open questions regarding how to tackle the scalability to high-dimensional problems and how to deal with the intrinsic metastability of the system. To explore new methods we link the optimal control problem to reinforcement learning since both share the same underlying framework, namely a Markov Decision Process (MDP). For the optimal control problem we show how the MDP can be formulated. In addition we discuss how the stochastic optimal control problem can be interpreted in the framework of reinforcement learning. At the end of the article we present the application of two different reinforcement learning algorithms to the optimal control problem and a comparison of the advantages and disadvantages of the two algorithms. Y1 - 2024 U6 - https://doi.org/10.1063/5.0140665 VL - 65 IS - 8 ER - TY - JOUR A1 - Ray, Sourav A1 - Fackeldey, Konstantin A1 - Stein, Christoph A1 - Weber, Marcus T1 - Coarse Grained MD Simulations of Opioid interactions with the µ-opioid receptor and the surrounding lipid membrane JF - Biophysica N2 - In our previous studies, a new opioid (NFEPP) was developed to only selectively bind to the 𝜇-opoid receptor (MOR) in inflamed tissue and thus avoid the severe side effects of fentanyl. We know that NFEPP has a reduced binding affinity to MOR in healthy tissue. Inspired by the modelling and simulations performed by Sutcliffe et al., we present our own results of coarse-grained molecular dynamics simulations of fentanyl and NFEPP with regards to their interaction with the 𝜇-opioid receptor embedded within the lipid cell membrane. For technical reasons, we have slightly modified Sutcliffe’s parametrisation of opioids. The pH-dependent opioid simulations are of interest because while fentanyl is protonated at the physiological pH, NFEPP is deprotonated due to its lower pKa value than that of fentanyl. Here, we analyse for the first time whether pH changes have an effect on the dynamical behaviour of NFEPP when it is inside the cell membrane. Besides these changes, our analysis shows a possible alternative interaction of NFEPP at pH 7.4 outside the binding region of the MOR. The interaction potential of NFEPP with MOR is also depicted by analysing the provided statistical molecular dynamics simulations with the aid of an eigenvector analysis of a transition rate matrix. In our modelling, we see differences in the XY-diffusion profiles of NFEPP compared with fentanyl in the cell membrane. Y1 - 2023 U6 - https://doi.org/10.3390/biophysica3020017 VL - 3 IS - 2 SP - 263 EP - 275 ER -