@article{WeberWeitereAutoren2020, author = {Weber, Marcus and Weitere Autoren,}, title = {DIN SPEC 2343: {\"U}bertragung von sprachbasierten Daten zwischen K{\"u}nstlichen Intelligenzen - Festlegung von Parametern und Formaten}, journal = {Beuth Verlag}, editor = {Reichardt, Christine}, year = {2020}, abstract = {Dieses Dokument legt Parameter und Formate f{\"u}r die {\"U}bertragung sprachbasierter Daten zwischen verschiedenen KI-{\"O}kosystemen fest.}, language = {de} } @article{RabbenRayWeber2020, author = {Rabben, Robert Julian and Ray, Sourav and Weber, Marcus}, title = {ISOKANN: Invariant subspaces of Koopman operators learned by a neural network}, volume = {153}, journal = {The Journal of Chemical Physics}, number = {11}, doi = {10.1063/5.0015132}, pages = {114109}, year = {2020}, abstract = {The problem of determining the rate of rare events in dynamical systems is quite well-known but still difficult to solve. Recent attempts to overcome this problem exploit the fact that dynamic systems can be represented by a linear operator, such as the Koopman operator. Mathematically, the rare event problem comes down to the difficulty in finding invariant subspaces of these Koopman operators K. In this article, we describe a method to learn basis functions of invariant subspaces using an artificial neural Network.}, language = {en} } @misc{SunkaraRaharinirinaPeppertetal.2020, author = {Sunkara, Vikram and Raharinirina, N. Alexia and Peppert, Felix and von Kleist, Max and Sch{\"u}tte, Christof}, title = {Inferring Gene Regulatory Networks from Single Cell RNA-seq Temporal Snapshot Data Requires Higher Order Moments}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-79664}, year = {2020}, abstract = {Due to the increase in accessibility and robustness of sequencing technology, single cell RNA-seq (scRNA-seq) data has become abundant. The technology has made significant contributions to discovering novel phenotypes and heterogeneities of cells. Recently, there has been a push for using single-- or multiple scRNA-seq snapshots to infer the underlying gene regulatory networks (GRNs) steering the cells' biological functions. To date, this aspiration remains unrealised. In this paper, we took a bottom-up approach and curated a stochastic two gene interaction model capturing the dynamics of a complete system of genes, mRNAs, and proteins. In the model, the regulation was placed upstream from the mRNA on the gene level. We then inferred the underlying regulatory interactions from only the observation of the mRNA population through~time. We could detect signatures of the regulation by combining information of the mean, covariance, and the skewness of the mRNA counts through time. We also saw that reordering the observations using pseudo-time did not conserve the covariance and skewness of the true time course. The underlying GRN could be captured consistently when we fitted the moments up to degree three; however, this required a computationally expensive non-linear least squares minimisation solver. There are still major numerical challenges to overcome for inference of GRNs from scRNA-seq data. These challenges entail finding informative summary statistics of the data which capture the critical regulatory information. Furthermore, the statistics have to evolve linearly or piece-wise linearly through time to achieve computational feasibility and scalability.}, language = {en} } @article{KlebanovSprungkSullivan2021, author = {Klebanov, Ilja and Sprungk, Bj{\"o}rn and Sullivan, T. J.}, title = {The linear conditional expectation in Hilbert space}, volume = {27}, journal = {Bernoulli}, number = {4}, arxiv = {http://arxiv.org/abs/2008.12070}, doi = {10.3150/20-BEJ1308}, pages = {2299 -- 2299}, year = {2021}, language = {en} } @article{MironBeronVeraHelfmannetal.2021, author = {Miron, Philippe and Beron-Vera, Francisco and Helfmann, Luzie and Koltai, P{\´e}ter}, title = {Transition paths of marine debris and the stability of the garbage patches}, volume = {31}, journal = {Chaos: An Interdisciplinary Journal of Nonlinear Science}, number = {3}, arxiv = {http://arxiv.org/abs/2009.11234}, doi = {https://doi.org/10.1063/5.0030535}, year = {2021}, abstract = {We used transition path theory (TPT) to infer "reactive" pathways of floating marine debris trajectories. The TPT analysis was applied on a pollution-aware time-homogeneous Markov chain model constructed from trajectories produced by satellite-tracked undrogued buoys from the NOAA Global Drifter Program. The latter involved coping with the openness of the system in physical space, which further required an adaptation of the standard TPT setting. Directly connecting pollution sources along coastlines with garbage patches of varied strengths, the unveiled reactive pollution routes represent alternative targets for ocean cleanup efforts. Among our specific findings we highlight: constraining a highly probable pollution source for the Great Pacific Garbage Patch; characterizing the weakness of the Indian Ocean gyre as a trap for plastic waste; and unveiling a tendency of the subtropical gyres to export garbage toward the coastlines rather than to other gyres in the event of anomalously intense winds.}, language = {en} } @article{MrowińskiSchnauberKaganskiyetal.2020, author = {Mrowiński, Paweł and Schnauber, Peter and Kaganskiy, Arsenty and Schall, Johannes and Burger, Sven and Rodt, Sven and Reitzenstein, Stephan}, title = {Directional single-photon emission from deterministic quantum dot waveguide structures}, volume = {14}, journal = {Phys. Status Solidi RRL}, arxiv = {http://arxiv.org/abs/2005.09322}, doi = {10.1002/pssr.202000115}, pages = {2000115}, year = {2020}, language = {en} } @misc{RaySunkaraSchuetteetal.2020, author = {Ray, Sourav and Sunkara, Vikram and Sch{\"u}tte, Christof and Weber, Marcus}, title = {How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs}, issn = {1438-0064}, doi = {10.1080/08927022.2020.1839660}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-78437}, year = {2020}, abstract = {Molecular simulations of ligand-receptor interactions are a computational challenge, especially when their association- (``on''-rate) and dissociation- (``off''-rate) mechanisms are working on vastly differing timescales. In addition, the timescale of the simulations themselves is, in practice, orders of magnitudes smaller than that of the mechanisms; which further adds to the complexity of observing these mechanisms, and of drawing meaningful and significant biological insights from the simulation. One way of tackling this multiscale problem is to compute the free-energy landscapes, where molecular dynamics (MD) trajectories are used to only produce certain statistical ensembles. The approach allows for deriving the transition rates between energy states as a function of the height of the activation-energy barriers. In this article, we derive the association rates of the opioids fentanyl and N-(3-fluoro-1-phenethylpiperidin-4-yl)- N-phenyl propionamide (NFEPP) in a \$\mu\$-opioid receptor by combining the free-energy landscape approach with the square-root-approximation method (SQRA), which is a particularly robust version of Markov modelling. The novelty of this work is that we derive the association rates as a function of the pH level using only an ensemble of MD simulations. We also verify our MD-derived insights by reproducing the in vitro study performed by the Stein Lab, who investigated the influence of pH on the inhibitory constant of fentanyl and NFEPP (Spahn et al. 2017). MD simulations are far more accessible and cost-effective than in vitro and in vivo studies. Especially in the context of the current opioid crisis, MD simulations can aid in unravelling molecular functionality and assist in clinical decision-making; the approaches presented in this paper are a pertinent step forward in this direction.}, language = {en} } @article{ShiEvlyukhinReinhardtetal.2020, author = {Shi, Liping and Evlyukhin, Andrey B. and Reinhardt, Carsten and Babushkin, Ihar and Zenin, Vladimir A. and Burger, Sven and Malureanu, Radu and Chichkov, Boris N. and Morgner, Uwe and Kovacev, Milutin}, title = {Progressive Self-Boosting Anapole-Enhanced Deep-Ultraviolet Third Harmonic During Few-Cycle Laser Radiation}, volume = {7}, journal = {ACS Photonics}, doi = {10.1021/acsphotonics.0c00753}, pages = {1655}, year = {2020}, language = {en} } @article{FackeldeyRoehmNiknejadetal.2021, author = {Fackeldey, Konstantin and R{\"o}hm, Jonas and Niknejad, Amir and Chewle, Surahit and Weber, Marcus}, title = {Analyzing Raman Spectral Data without Separabiliy Assumption}, volume = {3}, journal = {Journal of Mathematical Chemistry}, number = {59}, publisher = {Springer}, arxiv = {http://arxiv.org/abs/2007.06428}, doi = {10.1007/s10910-020-01201-7}, pages = {575 -- 596}, year = {2021}, abstract = {Raman spectroscopy is a well established tool for the analysis of vibration spectra, which then allow for the determination of individual substances in a chemical sample, or for their phase transitions. In the Time-Resolved-Raman-Sprectroscopy the vibration spectra of a chemical sample are recorded sequentially over a time interval, such that conclusions for intermediate products (transients) can be drawn within a chemical process. The observed data-matrix M from a Raman spectroscopy can be regarded as a matrix product of two unknown matrices W and H, where the first is representing the contribution of the spectra and the latter represents the chemical spectra. One approach for obtaining W and H is the non-negative matrix factorization. We propose a novel approach, which does not need the commonly used separability assumption. The performance of this approach is shown on a real world chemical example.}, language = {en} } @article{SiefkeHurtadoDickmannetal.2020, author = {Siefke, Thomas and Hurtado, Carol and Dickmann, Johannes and Dickmann, Walter and K{\"a}seberg, Tim and Meyer, Jan and Burger, Sven and Zeitner, Uwe and Bodermann, Bernd and Kroker, Stefanie}, title = {Quasi-bound states in the continuum for deep subwavelength structural information retrieval for DUV nano-optical polarizers}, volume = {28}, journal = {Opt. Express}, doi = {10.1364/OE.396044}, pages = {23122}, year = {2020}, language = {en} } @article{RoehlWeberFackeldey2021, author = {R{\"o}hl, Susanne and Weber, Marcus and Fackeldey, Konstantin}, title = {Computing the minimal rebinding effect for non-reversible processes}, volume = {19}, journal = {Multiscale Modeling and Simulation}, number = {1}, arxiv = {http://arxiv.org/abs/2007.08403}, doi = {https://doi.org/10.1137/20M1334966}, pages = {460 -- 477}, year = {2021}, abstract = {The aim of this paper is to investigate the rebinding effect, a phenomenon describing a "short-time memory" which can occur when projecting a Markov process onto a smaller state space. For guaranteeing a correct mapping by the Markov State Model, we assume a fuzzy clustering in terms of membership functions, assigning degrees of membership to each state. The macro states are represented by the membership functions and may be overlapping. The magnitude of this overlap is a measure for the strength of the rebinding effect, caused by the projection and stabilizing the system. A minimal bound for the rebinding effect included in a given system is computed as the solution of an optimization problem. Based on membership functions chosen as a linear combination of Schur vectors, this generalized approach includes reversible as well as non-reversible processes.}, language = {en} } @article{TockhornSutterColometal.2020, author = {Tockhorn, Philipp and Sutter, Johannes and Colom, R{\´e}mi and Kegelmann, Lukas and Al-Ashouri, Amran and Roß, Marcel and J{\"a}ger, Klaus and Unold, Thomas and Burger, Sven and Albrecht, Steve and Becker, Christiane}, title = {Improved Quantum Efficiency by Advanced Light Management in Nanotextured Solution-Processed Perovskite Solar Cells}, volume = {7}, journal = {ACS Photonics}, doi = {10.1021/acsphotonics.0c00935}, pages = {2589}, year = {2020}, language = {en} } @article{GorgullaDasLeighetal.2021, author = {Gorgulla, Christoph and Das, Krishna M. Padmanabha and Leigh, Kendra E and Cespugli, Marco and Fischer, Patrick D. and Wang, Zi-Fu and Tesseyre, Guilhem and Pandita, Shreya and Shnapir, Alex and Calderaio, Anthony and Hutcheson, Colin and Gechev, Minko and Rose, Alexander and Lewis, Noam and Yaffe, Erez and Luxenburg, Roni and Herce, Henry D. and Durmaz, Vedat and Halazonetis, Thanos D. and Fackeldey, Konstantin and Patten, Justin J. and Chuprina, Alexander and Dziuba, Igor and Plekhova, Alla and Moroz, Yurii and Radchenko, Dmytro and Tarkhanova, Olga and Yavnyuk, Irina and Gruber, Christian C. and Yust, Ryan and Payne, Dave and N{\"a}{\"a}r, Anders M. and Namchuk, Mark N. and Davey, Robert A. and Wagner, Gerhard and Kinney, Jamie and Arthanari, Haribabu}, title = {A Multi-Pronged Approach Targeting SARS-CoV-2 Proteins Using Ultra-Large Virtual Screening}, volume = {24}, journal = {iScience}, number = {2}, publisher = {CellPress}, doi = {10.26434/chemrxiv.12682316}, pages = {102021}, year = {2021}, abstract = {Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), previously known as 2019 novel coronavirus (2019-nCoV), has spread rapidly across the globe, creating an unparalleled global health burden and spurring a deepening economic crisis. As of July 7th, 2020, almost seven months into the outbreak, there are no approved vaccines and few treatments available. Developing drugs that target multiple points in the viral life cycle could serve as a strategy to tackle the current as well as future coronavirus pandemics. Here we leverage the power of our recently developed in silico screening platform, VirtualFlow, to identify inhibitors that target SARS-CoV-2. VirtualFlow is able to efficiently harness the power of computing clusters and cloud-based computing platforms to carry out ultra-large scale virtual screens. In this unprecedented structure-based multi-target virtual screening campaign, we have used VirtualFlow to screen an average of approximately 1 billion molecules against each of 40 different target sites on 17 different potential viral and host targets in the cloud. In addition to targeting the active sites of viral enzymes, we also target critical auxiliary sites such as functionally important protein-protein interaction interfaces. This multi-target approach not only increases the likelihood of finding a potent inhibitor, but could also help identify a collection of anti-coronavirus drugs that would retain efficacy in the face of viral mutation. Drugs belonging to different regimen classes could be combined to develop possible combination therapies, and top hits that bind at highly conserved sites would be potential candidates for further development as coronavirus drugs. Here, we present the top 200 in silico hits for each target site. While in-house experimental validation of some of these compounds is currently underway, we want to make this array of potential inhibitor candidates available to researchers worldwide in consideration of the pressing need for fast-tracked drug development.}, language = {en} } @article{BonnetAkianSavinetal.2020, author = {Bonnet, Luc and Akian, Jean-Luc and Savin, {\´E}ric and Sullivan, T. J.}, title = {Adaptive reconstruction of imperfectly-observed monotone functions, with applications to uncertainty quantification}, volume = {13}, journal = {Algorithms}, number = {8}, arxiv = {http://arxiv.org/abs/2007.05236}, doi = {10.3390/a13080196}, pages = {196}, year = {2020}, abstract = {Motivated by the desire to numerically calculate rigorous upper and lower bounds on deviation probabilities over large classes of probability distributions, we present an adaptive algorithm for the reconstruction of increasing real-valued functions. While this problem is similar to the classical statistical problem of isotonic regression, the optimisation setting alters several characteristics of the problem and opens natural algorithmic possibilities. We present our algorithm, establish sufficient conditions for convergence of the reconstruction to the ground truth, and apply the method to synthetic test cases and a real-world example of uncertainty quantification for aerodynamic design.}, language = {en} } @article{WuerthManleyVoigtetal.2020, author = {W{\"u}rth, Christian and Manley, Phillip and Voigt, Robert and Ahiboz, Doguscan and Becker, Christiane and Resch-Genger, Ute}, title = {Metasurface enhanced sensitized photon upconversion: towards highly efficient low power upconversion applications and nano-scale E-field sensors}, volume = {20}, journal = {Nano Lett.}, doi = {10.1021/acs.nanolett.0c02548}, pages = {6682}, year = {2020}, language = {en} } @article{AlhaddadFoerstnerGrothetal.2021, author = {Alhaddad, Samer and F{\"o}rstner, Jens and Groth, Stefan and Gr{\"u}newald, Daniel and Grynko, Yevgen and Hannig, Frank and Kenter, Tobias and Pfreundt, Franz-Josef and Plessl, Christian and Schotte, Merlind and Steinke, Thomas and Teich, J{\"u}rgen and Weiser, Martin and Wende, Florian}, title = {HighPerMeshes - A Domain-Specific Language for Numerical Algorithms on Unstructured Grids}, journal = {Euro-Par 2020: Parallel Processing Workshops.}, publisher = {Springer}, doi = {10.1007/978-3-030-71593-9_15}, pages = {185 -- 196}, year = {2021}, abstract = {Solving partial differential equations on unstructured grids is a cornerstone of engineering and scientific computing. Nowadays, heterogeneous parallel platforms with CPUs, GPUs, and FPGAs enable energy-efficient and computationally demanding simulations. We developed the HighPerMeshes C++-embedded Domain-Specific Language (DSL) for bridging the abstraction gap between the mathematical and algorithmic formulation of mesh-based algorithms for PDE problems on the one hand and an increasing number of heterogeneous platforms with their different parallel programming and runtime models on the other hand. Thus, the HighPerMeshes DSL aims at higher productivity in the code development process for multiple target platforms. We introduce the concepts as well as the basic structure of the HighPer-Meshes DSL, and demonstrate its usage with three examples, a Poisson and monodomain problem, respectively, solved by the continuous finite element method, and the discontinuous Galerkin method for Maxwell's equation. The mapping of the abstract algorithmic description onto parallel hardware, including distributed memory compute clusters is presented. Finally, the achievable performance and scalability are demonstrated for a typical example problem on a multi-core CPU cluster.}, language = {en} } @inproceedings{BinkowskiBetzColometal.2020, author = {Binkowski, Felix and Betz, Fridtjof and Colom, R{\´e}mi and Hammerschmidt, Martin and Zschiedrich, Lin and Burger, Sven}, title = {Modal expansion of optical far-field quantities using quasinormal modes}, volume = {238}, booktitle = {EPJ Web Conf.}, doi = {10.1051/epjconf/202023805007}, pages = {05007}, year = {2020}, language = {en} } @inproceedings{ColomBinkowskiBetzetal.2020, author = {Colom, R{\´e}mi and Binkowski, Felix and Betz, Fridtjof and Hammerschmidt, Martin and Zschiedrich, Lin and Burger, Sven}, title = {Quasi-normal mode expansion as a tool for the design of nanophotonic devices}, volume = {238}, booktitle = {EPJ Web Conf.}, doi = {10.1051/epjconf/202023805008}, pages = {05008}, year = {2020}, language = {en} } @article{RaySunkaraSchuetteetal.2020, author = {Ray, Sourav and Sunkara, Vikram and Sch{\"u}tte, Christof and Weber, Marcus}, title = {How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs}, volume = {46}, journal = {Molecular Simulation}, number = {18}, publisher = {Taylor and Francis}, doi = {10.1080/08927022.2020.1839660}, pages = {1443 -- 1452}, year = {2020}, abstract = {Molecular simulations of ligand-receptor interactions are a computational challenge, especially when their association- ('on'-rate) and dissociation- ('off'-rate) mechanisms are working on vastly differing timescales. One way of tackling this multiscale problem is to compute the free-energy landscapes, where molecular dynamics (MD) trajectories are used to only produce certain statistical ensembles. The approach allows for deriving the transition rates between energy states as a function of the height of the activation-energy barriers. In this article, we derive the association rates of the opioids fentanyl and N-(3-fluoro-1-phenethylpiperidin-4-yl)-N-phenyl propionamide (NFEPP) in a μ-opioid receptor by combining the free-energy landscape approach with the square-root-approximation method (SQRA), which is a particularly robust version of Markov modelling. The novelty of this work is that we derive the association rates as a function of the pH level using only an ensemble of MD simulations. We also verify our MD-derived insights by reproducing the in vitro study performed by the Stein Lab.}, language = {en} } @misc{Betz2020, type = {Master Thesis}, author = {Betz, Fridtjof}, title = {Investigation and application of contour integral methods for resonance expansion in nanophotonics}, year = {2020}, language = {en} } @article{SantiagoBurgerRockstuhletal.2021, author = {Santiago, Xavier Garcia and Burger, Sven and Rockstuhl, Carsten and Schneider, Philipp-Immanuel}, title = {Bayesian optimization with improved scalability and derivative information for efficient design of nanophotonic structures}, volume = {39}, journal = {J. Light. Technol.}, arxiv = {http://arxiv.org/abs/2101.02972}, doi = {10.1109/JLT.2020.3023450}, pages = {167}, year = {2021}, language = {en} } @article{WeberFischerDamerauetal.2020, author = {Weber, Marie-Christin and Fischer, Lisa and Damerau, Alexandra and Ponomarev, Igor and Pfeiffenberger, Moritz and Gaber, Timo and G{\"o}tschel, Sebastian and Lang, Jens and R{\"o}blitz, Susanna and Buttgereit, Frank and Ehrig, Rainald and Lang, Annemarie}, title = {Macroscale mesenchymal condensation to study cytokine-driven cellular and matrix-related changes during cartilage degradation}, volume = {12}, journal = {Biofabrication}, number = {4}, doi = {10.1088/1758-5090/aba08f}, year = {2020}, abstract = {Understanding the pathophysiological processes of cartilage degradation requires adequate model systems to develop therapeutic strategies towards osteoarthritis (OA). Although different in vitro or in vivo models have been described, further comprehensive approaches are needed to study specific disease aspects. This study aimed to combine in vitro and in silico modeling based on a tissue-engineering approach using mesenchymal condensation to mimic cytokine-induced cellular and matrix-related changes during cartilage degradation. Thus, scaffold-free cartilage-like constructs (SFCCs) were produced based on self-organization of mesenchymal stromal cells (mesenchymal condensation) and i) characterized regarding their cellular and matrix composition or secondly ii) treated with interleukin-1β (IL-1β) and tumor necrosis factor α (TNFα) for 3 weeks to simulate OA-related matrix degradation. In addition, an existing mathematical model based on partial differential equations was optimized and transferred to the underlying settings to simulate distribution of IL-1β, type II collagen degradation and cell number reduction. By combining in vitro and in silico methods, we aim to develop a valid, efficient alternative approach to examine and predict disease progression and effects of new therapeutics.}, language = {en} } @article{AlvesEhrigRaffaltetal.2020, author = {Alves, S{\´o}nia and Ehrig, Rainald and Raffalt, Peter C. and Bender, Alwina and Duda, Georg and Agres, Alison N.}, title = {Quantifying Asymmetry in Gait: A Weighted Universal Symmetry Index to Evaluate 3D Ground Reaction Forces}, journal = {Frontiers in Bioengineering and Biotechnology}, year = {2020}, abstract = {Though gait asymmetry is used as a metric of functional recovery in clinical rehabilitation, there is no consensus on an ideal method for its evaluation. Various methods have been proposed but are limited in scope, as they can often use only positive signals or discrete values extracted from time-scale data as input. By defining five symmetry axioms, a framework for benchmarking existing methods was established and a new method was described here for the first time: the weighted universal symmetry index (wUSI), which overcomes limitations of other methods. Both existing methods and the wUSI were mathematically compared to each other and in respect to their ability to fulfill the proposed symmetry axioms. Eligible methods that fulfilled these axioms were then applied using both discrete and continuous approaches to ground reaction force (GRF) data collected from healthy gait, both with and without artificially induced asymmetry using a single instrumented elbow crutch. The wUSI with a continuous approach was the only symmetry method capable of determining GRF asymmetries in different walking conditions in all three planes of motion. When used with a continuous approach, the wUSI method was able to detect asymmetries while avoiding artificial inflation, a common problem reported in other methods. In conclusion, the wUSI is proposed as a universal method to quantify three-dimensional GRF asymmetries, which may also be expanded to other biomechanical signals.}, language = {en} } @article{WulkowKoltaiSchuette2021, author = {Wulkow, Niklas and Koltai, P{\´e}ter and Sch{\"u}tte, Christof}, title = {Memory-Based Reduced Modelling and Data-Based Estimation of Opinion Spreading}, volume = {31}, journal = {Journal of Nonlinear Science}, doi = {10.1007/s00332-020-09673-2}, year = {2021}, abstract = {We investigate opinion dynamics based on an agent-based model and are interested in predicting the evolution of the percentages of the entire agent population that share an opinion. Since these opinion percentages can be seen as an aggregated observation of the full system state, the individual opinions of each agent, we view this in the framework of the Mori-Zwanzig projection formalism. More specifically, we show how to estimate a nonlinear autoregressive model (NAR) with memory from data given by a time series of opinion percentages, and discuss its prediction capacities for various specific topologies of the agent interaction network. We demonstrate that the inclusion of memory terms significantly improves the prediction quality on examples with different network topologies.}, language = {en} } @article{StraubeKowalikNetzetal.2020, author = {Straube, Arthur and Kowalik, Bartosz G. and Netz, Roland R. and H{\"o}fling, Felix}, title = {Rapid onset of molecular friction in liquids bridging between the atomistic and hydrodynamic pictures}, volume = {3}, journal = {Commun. Phys.}, publisher = {Nature}, doi = {10.1038/s42005-020-0389-0}, pages = {126}, year = {2020}, abstract = {Friction in liquids arises from conservative forces between molecules and atoms. Although the hydrodynamics at the nanoscale is subject of intense research and despite the enormous interest in the non-Markovian dynamics of single molecules and solutes, the onset of friction from the atomistic scale so far could not be demonstrated. Here, we fill this gap based on frequency-resolved friction data from high-precision simulations of three prototypical liquids, including water. Combining with theory, we show that friction in liquids emerges abruptly at a characteristic frequency, beyond which viscous liquids appear as non-dissipative, elastic solids. Concomitantly, the molecules experience Brownian forces that display persistent correlations. A critical test of the generalised Stokes-Einstein relation, mapping the friction of single molecules to the visco-elastic response of the macroscopic sample, disproves the relation for Newtonian fluids, but substantiates it exemplarily for water and a moderately supercooled liquid. The employed approach is suitable to yield insights into vitrification mechanisms and the intriguing mechanical properties of soft materials.}, language = {en} } @article{FischerEhrigSchaeferetal.2021, author = {Fischer, Sophie and Ehrig, Rainald and Sch{\"a}fer, Stefan and Tronci, Enrico and Mancini, Toni and Egli, Marcel and Ille, Fabian and Kr{\"u}ger, Tillmann H. C. and Leeners, Brigitte and R{\"o}blitz, Susanna}, title = {Mathematical Modelling and Simulation Provides Evidence for New Strategies of Ovarian Stimulation}, volume = {12}, journal = {Frontiers in Endocrinology}, doi = {10.3389/fendo.2021.613048}, year = {2021}, abstract = {New approaches to ovarian stimulation protocols, such as luteal start, random start or double stimulation, allow for flexibility in ovarian stimulation at different phases of the menstrual cycle which is especially useful when time for assisted reproductive technology is limited, e.g. for emergency fertility preservation in cancer patients. It has been proposed that the success of these methods is based on the continuous growth of multiple cohorts ("waves") of follicles throughout the menstrual cycle which leads to the availability of ovarian follicles for ovarian controlled stimulation at several time points. Though several preliminary studies have been published, their scientific evidence has not been considered as being strong enough to integrate these results into routine clinical practice. This work aims at adding further scientific evidence about the efficiency of variable-start protocols and underpinning the theory of follicular waves by using mathematical modelling and numerical simulations. For this purpose, we have modified and coupled two previously published models, one describing the time course of hormones and one describing competitive follicular growth in a normal menstrual cycle. The coupled model is used to test stimulation protocols in silico. Simulation results show the occurrence of follicles in a wave-like manner during a normal menstrual cycle and qualitatively predict the outcome of ovarian stimulation initiated at different time points of the menstrual cycle.}, language = {en} } @article{FackeldeyOsterSallandtetal.2022, author = {Fackeldey, Konstantin and Oster, Mathias and Sallandt, Leon and Schneider, Reinhold}, title = {Approximative Policy Iteration for Exit Time Feedback Control Problems driven by Stochastic Differential Equations using Tensor Train format}, volume = {20}, journal = {SIAM Journal on Multiscale Modeling and Simulation}, number = {1}, arxiv = {http://arxiv.org/abs/2010.04465}, doi = {10.1137/20M1372500}, pages = {379 -- 403}, year = {2022}, abstract = {We consider a stochastic optimal exit time feedback control problem. The Bellman equation is solved approximatively via the Policy Iteration algorithm on a polynomial ansatz space by a sequence of linear equations. As high degree multi-polynomials are needed, the corresponding equations suffer from the curse of dimensionality even in moderate dimensions. We employ tensor-train methods to account for this problem. The approximation process within the Policy Iteration is done via a Least-Squares ansatz and the integration is done via Monte-Carlo methods. Numerical evidences are given for the (multi dimensional) double well potential and a three-hole potential.}, language = {en} } @article{HelfmannDjurdjevacConradDjurdjevacetal.2021, author = {Helfmann, Luzie and Djurdjevac Conrad, Natasa and Djurdjevac, Ana and Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {From interacting agents to density-based modeling with stochastic PDEs}, volume = {16}, journal = {Communications in Applied Mathematics and Computational Science}, number = {1}, arxiv = {http://arxiv.org/abs/1905.13525}, doi = {10.2140/camcos.2021.16.1}, pages = {1 -- 32}, year = {2021}, abstract = {Many real-world processes can naturally be modeled as systems of interacting agents. However, the long-term simulation of such agent-based models is often intractable when the system becomes too large. In this paper, starting from a stochastic spatio-temporal agent-based model (ABM), we present a reduced model in terms of stochastic PDEs that describes the evolution of agent number densities for large populations. We discuss the algorithmic details of both approaches; regarding the SPDE model, we apply Finite Element discretization in space which not only ensures efficient simulation but also serves as a regularization of the SPDE. Illustrative examples for the spreading of an innovation among agents are given and used for comparing ABM and SPDE models.}, language = {en} } @article{ChewleEmmerlingWeber2020, author = {Chewle, Surahit and Emmerling, Franziska and Weber, Marcus}, title = {Effect of choice of solvent on crystallization pathway of Paracetamol: An experimental and theoretical case study}, volume = {10}, journal = {Crystals}, number = {12}, doi = {10.3390/cryst10121107}, pages = {1107}, year = {2020}, abstract = {The choice of solvents influences crystalline solid formed during the crystallization of active pharmaceutical ingredients (API). The underlying effects are not always well understood because of the complexity of the systems. Theoretical models are often insufficient to describe this phenomenon. In this study, the crystallization behavior of the model drug paracetamol in different solvents was studied based on experimental and molecular dynamics data. The crystallization process was followed in situ using time-resolved Raman spectroscopy. Molecular dynamics with simulated annealing algorithm was used for an atomistic understanding of the underlying processes. The experimental and theoretical data indicate that paracetamol molecules adopt a particular geometry in a given solvent predefining the crystallization of certain polymorphs.}, language = {en} } @article{NiemannKlusSchuette2021, author = {Niemann, Jan-Hendrik and Klus, Stefan and Sch{\"u}tte, Christof}, title = {Data-driven model reduction of agent-based systems using the Koopman generator}, volume = {16}, journal = {PLOS ONE}, number = {5}, arxiv = {http://arxiv.org/abs/2012.07718}, doi = {10.1371/journal.pone.0250970}, year = {2021}, abstract = {The dynamical behavior of social systems can be described by agent-based models. Although single agents follow easily explainable rules, complex time-evolving patterns emerge due to their interaction. The simulation and analysis of such agent-based models, however, is often prohibitively time-consuming if the number of agents is large. In this paper, we show how Koopman operator theory can be used to derive reduced models of agent-based systems using only simulation or real-world data. Our goal is to learn coarse-grained models and to represent the reduced dynamics by ordinary or stochastic differential equations. The new variables are, for instance, aggregated state variables of the agent-based model, modeling the collective behavior of larger groups or the entire population. Using benchmark problems with known coarse-grained models, we demonstrate that the obtained reduced systems are in good agreement with the analytical results, provided that the numbers of agents is sufficiently large.}, language = {en} } @misc{BeckerJaegerBurger2021, author = {Becker, Christiane and J{\"a}ger, Klaus and Burger, Sven}, title = {Nanophotonics for Solar Energy}, volume = {221}, journal = {Sol. Energy Mater. Sol. Cells}, doi = {10.1016/j.solmat.2020.110916}, pages = {110916}, year = {2021}, language = {en} } @article{KlebanovSchuster2020, author = {Klebanov, Ilja and Schuster, Ingmar}, title = {Markov Chain Importance Sampling - a highly efficient estimator for MCMC}, journal = {Journal of Computational and Graphical Statistics}, doi = {10.1080/10618600.2020.1826953}, year = {2020}, abstract = {Markov chain (MC) algorithms are ubiquitous in machine learning and statistics and many other disciplines. Typically, these algorithms can be formulated as acceptance rejection methods. In this work we present a novel estimator applicable to these methods, dubbed Markov chain importance sampling (MCIS), which efficiently makes use of rejected proposals. For the unadjusted Langevin algorithm, it provides a novel way of correcting the discretization error. Our estimator satisfies a central limit theorem and improves on error per CPU cycle, often to a large extent. As a by-product it enables estimating the normalizing constant, an important quantity in Bayesian machine learning and statistics.}, language = {en} } @article{ScottSantiagoBeuteletal.2020, author = {Scott, Philip and Santiago, Xavier Garcia and Beutel, Dominik and Rockstuhl, Carsten and Wegener, Martin and Fernandez-Corbaton, Ivan}, title = {On enhanced sensing of chiral molecules in optical cavities}, volume = {7}, journal = {Appl. Phys. Rev.}, doi = {10.1063/5.0025006}, pages = {041413}, year = {2020}, language = {en} } @article{SchneckWeiserWende2021, author = {Schneck, Jakob and Weiser, Martin and Wende, Florian}, title = {Impact of mixed precision and storage layout on additive Schwarz smoothers}, volume = {28}, journal = {Numerical Linear Algebra with Applications}, number = {4}, doi = {10.1002/nla.2366}, year = {2021}, abstract = {The growing discrepancy between CPU computing power and memory bandwidth drives more and more numerical algorithms into a bandwidth-bound regime. One example is the overlapping Schwarz smoother, a highly effective building block for iterative multigrid solution of elliptic equations with higher order finite elements. Two options of reducing the required memory bandwidth are sparsity exploiting storage layouts and representing matrix entries with reduced precision in floating point or fixed point format. We investigate the impact of several options on storage demand and contraction rate, both analytically in the context of subspace correction methods and numerically at an example of solid mechanics. Both perspectives agree on the favourite scheme: fixed point representation of Cholesky factors in nested dissection storage.}, language = {en} } @article{WulkowConradDjurdjevacConradetal.2021, author = {Wulkow, Hanna and Conrad, Tim and Djurdjevac Conrad, Natasa and M{\"u}ller, Sebastian A. and Nagel, Kai and Sch{\"u}tte, Christof}, title = {Prediction of Covid-19 spreading and optimal coordination of counter-measures: From microscopic to macroscopic models to Pareto fronts}, volume = {16}, journal = {PLOS One}, number = {4}, publisher = {Public Library of Science}, doi = {10.1371/journal.pone.0249676}, year = {2021}, language = {en} } @article{ZimperDjurdjevacHartmannetal.2025, author = {Zimper, Sebastian and Djurdjevac, Ana and Hartmann, Carsten and Sch{\"u}tte, Christof and Conrad, Natasa Djurdjevac}, title = {Mean-field optimal control with stochastic leaders}, arxiv = {http://arxiv.org/abs/2512.19201}, year = {2025}, abstract = {We consider interacting agent systems with a large number of stochastic agents (or particles) influenced by a fixed number of external stochastic lead agents. Such examples arise, for example in models of opinion dynamics, where a small number of leaders (influencers) can steer the behaviour of a large population of followers. In this context, we study a partial mean-field limit where the number of followers tends to infinity, while the number of leaders stays constant. The partial mean-field limit dynamics is then given by a McKean-Vlasov stochastic differential equation (SDE) for the followers, coupled to a controlled It{\^o}-SDE governing the dynamics of the lead agents. For a given cost functional that the lead agents seek to minimise, we show that the unique optimal control of the finite agent system convergences to the optimal control of the limiting system. This establishes that the low-dimensional control of the partial (mean-field) system provides an effective approximation for controlling the high-dimensional finite agent system. In addition, we propose a stochastic gradient descent algorithm that can efficiently approximate the mean-field control. Our theoretical results are illustrated on opinion dynamics model with lead agents, where the control objective is to drive the followers to reach consensus in finite time.}, language = {en} } @article{HartmannJoesterSchuetteetal.2026, author = {Hartmann, Carsten and J{\"o}ster, Annika and Sch{\"u}tte, Christof and Sikorski, Alexander and Weber, Marcus}, title = {Importance sampling of unbounded random stopping times: computing committor functions and exit rates without reweighting}, arxiv = {http://arxiv.org/abs/2601.01489}, year = {2026}, abstract = {Rare events in molecular dynamics are often related to noise-induced transitions between different macroscopic states (e.g., in protein folding). A common feature of these rare transitions is that they happen on timescales that are on average exponentially long compared to the characteristic timescale of the system, with waiting time distributions that have (sub)exponential tails and infinite support. As a result, sampling such rare events can lead to trajectories that can be become arbitrarily long, with not too low probability, which makes the reweighting of such trajectories a real challenge. Here, we discuss rare event simulation by importance sampling from a variational perspective, with a focus on applications in molecular dynamics, in particular the computation of committor functions. The idea is to design importance sampling schemes that (a) reduce the variance of a rare event estimator while controlling the average length of the trajectories and (b) that do not require the reweighting of possibly very long trajectories. In doing so, we study different stochastic control formulations for committor and mean first exit times, which we compare both from a theoretical and a computational point of view, including numerical studies of some benchmark examples.}, language = {en} } @article{GarciaSantiagoHammerschmidtSachsetal.2022, author = {Garcia Santiago, Xavier and Hammerschmidt, Martin and Sachs, Johannes and Burger, Sven and Kwon, Hyunah and Kn{\"o}ller, Marvin and Arens, Tilo and Fischer, Peer and Fernandez-Corbaton, Ivan and Rockstuhl, Carsten}, title = {Toward maximally electromagnetically chiral scatterers at optical frequencies}, volume = {9}, journal = {ACS Photonics}, arxiv = {http://arxiv.org/abs/2112.04422}, doi = {10.1021/acsphotonics.1c01887}, pages = {1954}, year = {2022}, language = {en} } @misc{ErnstSchuetteSigristetal.2021, author = {Ernst, Ariane and Sch{\"u}tte, Christof and Sigrist, Stephan and Winkelmann, Stefanie}, title = {Variance of filtered signals: Characterization for linear reaction networks and application to neurotransmission dynamics}, issn = {1438-0064}, doi = {10.1016/j.mbs.2021.108760}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-82674}, year = {2021}, abstract = {Neurotransmission at chemical synapses relies on the calcium-induced fusion of synaptic vesicles with the presynaptic membrane. The distance to the calcium channels determines the release probability and thereby the postsynaptic signal. Suitable models of the process need to capture both the mean and the variance observed in electrophysiological measurements of the postsynaptic current. In this work, we propose a method to directly compute the exact first- and second-order moments for signals generated by a linear reaction network under convolution with an impulse response function, rendering computationally expensive numerical simulations of the underlying stochastic counting process obsolete. We show that the autocorrelation of the process is central for the calculation of the filtered signal's second-order moments, and derive a system of PDEs for the cross-correlation functions (including the autocorrelations) of linear reaction networks with time-dependent rates. Finally, we employ our method to efficiently compare different spatial coarse graining approaches for a specific model of synaptic vesicle fusion. Beyond the application to neurotransmission processes, the developed theory can be applied to any linear reaction system that produces a filtered stochastic signal.}, language = {en} } @article{ErnstSchuetteSigristetal.2022, author = {Ernst, Ariane and Sch{\"u}tte, Christof and Sigrist, Stephan and Winkelmann, Stefanie}, title = {Variance of filtered signals: Characterization for linear reaction networks and application to neurotransmission dynamics}, volume = {343}, journal = {Mathematical Biosciences}, doi = {10.1016/j.mbs.2021.108760}, year = {2022}, abstract = {Neurotransmission at chemical synapses relies on the calcium-induced fusion of synaptic vesicles with the presynaptic membrane. The distance to the calcium channels determines the release probability and thereby the postsynaptic signal. Suitable models of the process need to capture both the mean and the variance observed in electrophysiological measurements of the postsynaptic current. In this work, we propose a method to directly compute the exact first- and second-order moments for signals generated by a linear reaction network under convolution with an impulse response function, rendering computationally expensive numerical simulations of the underlying stochastic counting process obsolete. We show that the autocorrelation of the process is central for the calculation of the filtered signal's second-order moments, and derive a system of PDEs for the cross-correlation functions (including the autocorrelations) of linear reaction networks with time-dependent rates. Finally, we employ our method to efficiently compare different spatial coarse graining approaches for a specific model of synaptic vesicle fusion. Beyond the application to neurotransmission processes, the developed theory can be applied to any linear reaction system that produces a filtered stochastic signal.}, language = {en} } @article{ThiesSunkaraRayetal.2023, author = {Thies, Arne and Sunkara, Vikram and Ray, Sourav and Wulkow, Hanna and Celik, M. {\"O}zg{\"u}r and Yerg{\"o}z, Fatih and Sch{\"u}tte, Christof and Stein, Christoph and Weber, Marcus and Winkelmann, Stefanie}, title = {Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design}, volume = {13}, journal = {Scientific Reports}, number = {607}, doi = {10.1038/s41598-023-27699-w}, year = {2023}, abstract = {We previously reported the successful design, synthesis and testing of the prototype opioid painkiller NFEPP that does not elicit adverse side effects. The design process of NFEPP was based on mathematical modelling of extracellular interactions between G-protein coupled receptors (GPCRs) and ligands, recognizing that GPCRs function differently under pathological versus healthy conditions. We now present an additional and novel stochastic model of GPCR function that includes intracellular dissociation of G-protein subunits and modulation of plasma membrane calcium channels and their dependence on parameters of inflamed and healthy tissue (pH, radicals). The model is validated against in vitro experimental data for the ligands NFEPP and fentanyl at different pH values and radical concentrations. We observe markedly reduced binding affinity and calcium channel inhibition for NFEPP at normal pH compared to lower pH, in contrast to the effect of fentanyl. For increasing radical concentrations, we find enhanced constitutive G-protein activation but reduced ligand binding affinity. Assessing the different effects, the results suggest that, compared to radicals, low pH is a more important determinant of overall GPCR function in an inflamed environment. Future drug design efforts should take this into account.}, language = {en} } @misc{RayThiesSunkaraetal.2021, author = {Ray, Sourav and Thies, Arne and Sunkara, Vikram and Wulkow, Hanna and Celik, {\"O}zg{\"u}r and Yerg{\"o}z, Fatih and Sch{\"u}tte, Christof and Stein, Christoph and Weber, Marcus and Winkelmann, Stefanie}, title = {Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-82797}, year = {2021}, abstract = {Initiated by mathematical modelling of extracellular interactions between G-protein coupled receptors (GPCRs) and ligands in normal versus diseased (inflamed) environments, we previously reported the successful design, synthesis and testing of the prototype opioid painkiller NFEPP that does not elicit adverse side effects. Uniquely, this design recognised that GPCRs function differently under pathological versus healthy conditions. We now present a novel stochastic model of GPCR function that includes intracellular dissociation of G-protein subunits and modulation of plasma membrane calcium channels associated with parameters of inflamed tissue (pH, radicals). By means of molecular dynamics simulations, we also assessed qualitative changes of the reaction rates due to additional disulfide bridges inside the GPCR binding pocket and used these rates for stochastic simulations of the corresponding reaction jump process. The modelling results were validated with in vitro experiments measuring calcium currents and G-protein activation. We found markedly reduced G-protein dissociation and calcium channel inhibition induced by NFEPP at normal pH, and enhanced constitutive G-protein activation but lower probability of ligand binding with increasing radical concentrations. These results suggest that, compared to radicals, low pH is a more important determinant of overall GPCR function in an inflamed environment. Future drug design efforts should take this into account.}, language = {en} } @phdthesis{Pakhomov2021, author = {Pakhomov, Anton}, title = {Efficient modeling and optimization of surface second-harmonic generation from nanophotonic components}, doi = {10.22032/dbt.49369}, year = {2021}, language = {en} } @article{PrasadTiwariRayetal.2022, author = {Prasad, Anup K. and Tiwari, Chandni and Ray, Sourav and Holden, Stephanie and Armstrong, David A. and Rosengren, K. Johan and Rodger, Alison and Panwar, Ajay S. and Martin, Lisandra L.}, title = {Secondary Structure Transitions for a Family of Amyloidogenic, Antimircobial Uperin 3 Peptides in Contact with Sodium Dodecyl Sulfate}, volume = {87}, journal = {ChemPlusChem}, doi = {10.1002/cplu.202100408}, pages = {e202100408}, year = {2022}, language = {en} } @article{BergPloentzkeSiebertetal.2021, author = {Berg, Mascha and Pl{\"o}ntzke, Julia and Siebert, Heike and R{\"o}blitz, Susanna}, title = {Modelling Oscillatory Patterns in the Bovine Estrous Cycle with Boolean Delay Equations}, volume = {83}, journal = {Bulletin of Mathematical Biology}, number = {121}, doi = {10.1007/s11538-021-00942-z}, pages = {1 -- 25}, year = {2021}, abstract = {Boolean delay equations (BDEs), with their relatively simple and intuitive mode of modelling, have been used in many research areas including, for example, climate dynamics and earthquake propagation. Their application to biological systems has been scarce and limited to the molecular level. Here, we derive and present two BDE models. One is directly derived from a previously published ordinary differential equation (ODE) model for the bovine estrous cycle, whereas the second model includes a modification of a particular biological mechanism. We not only compare the simulation results from the BDE models with the trajectories of the ODE model, but also validate the BDE models with two additional numerical experiments. One experiment induces a switch in the oscillatory pattern upon changes in the model parameters, and the other simulates the administration of a hormone that is known to shift the estrous cycle in time. The models presented here are the first BDE models for hormonal oscillators, and the first BDE models for drug administration. Even though automatic parameter estimation still remains challenging, our results support the role of BDEs as a framework for the systematic modelling of complex biological oscillators.}, language = {en} } @article{Djurdjevac2021, author = {Djurdjevac, Ana}, title = {Linear Parabolic Problems in Random Moving Domains}, volume = {9}, number = {2}, doi = {https://doi.org/10.1137/19M1284889}, pages = {848 -- 879}, year = {2021}, abstract = {We consider linear parabolic equations on a random non-cylindrical domain. Utilizing the domain mapping method, we write the problem as a partial differential equation with random coefficients on a cylindrical deterministic domain. Exploiting the deterministic results concerning equations on non-cylindrical domains, we state the necessary assumptions about the velocity filed and in addition, about the flow transformation that this field generates. In this paper we consider both cases, the uniformly bounded with respect to the sample and log-normal type transformation. In addition, we give an explicit example of a log-normal type transformation and prove that it does not satisfy the uniformly bounded condition. We define a general framework for considering linear parabolic problems on random non-cylindrical domains. As the first example, we consider the heat equation on a random tube domain and prove its well-posedness. Moreover, as the other example we consider the parabolic Stokes equation which illustrates the case when it is not enough just to study the plain-back transformation of the function, but instead to consider for example the Piola type transformation, in order to keep the divergence free property.}, language = {en} } @article{DjurdjevacGraeserHerbert2021, author = {Djurdjevac, Ana and Gr{\"a}ser, Carsten and Herbert, Philip J.}, title = {Non-stationary incompressible linear fluid equations in a moving domain}, journal = {arXiv}, arxiv = {http://arxiv.org/abs/2106.08312}, year = {2021}, abstract = {This article considers non-stationary incompressible linear fluid equations in a moving domain. We demonstrate the existence and uniqueness of an appropriate weak formulation of the problem by making use of the theory of time-dependent Bochner spaces. It is not possible to directly apply established evolving Hilbert space theory due to the incompressibility constraint. After we have established the well-posedness, we derive and analyse a time discretisation of the system.}, language = {en} } @article{ColomMikheevaAchourietal.2023, author = {Colom, Remi and Mikheeva, Elena and Achouri, Karim and Zuniga-Perez, Jesus and Bonod, Nicolas and Martin, Olivier J. F. and Burger, Sven and Genevet, Patrice}, title = {Crossing of the branch cut: the topological origin of a universal 2π-phase retardation in non-Hermitian metasurfaces}, volume = {17}, journal = {Laser Photonics Rev.}, arxiv = {http://arxiv.org/abs/2202.05632}, doi = {10.1002/lpor.202200976}, pages = {2200976}, year = {2023}, language = {en} } @article{MontefuscoSchuetteWinkelmann2023, author = {Montefusco, Alberto and Sch{\"u}tte, Christof and Winkelmann, Stefanie}, title = {A route to the hydrodynamic limit of a reaction-diffusion master equation using gradient structures}, volume = {83}, journal = {SIAM Journal on Applied Mathematics}, number = {2}, arxiv = {http://arxiv.org/abs/2201.02613}, doi = {10.1137/22M1488831}, pages = {837 -- 861}, year = {2023}, abstract = {The reaction-diffusion master equation (RDME) is a lattice-based stochastic model for spatially resolved cellular processes. It is often interpreted as an approximation to spatially continuous reaction-diffusion models, which, in the limit of an infinitely large population, may be described by means of reaction-diffusion partial differential equations. Analyzing and understanding the relation between different mathematical models for reaction-diffusion dynamics is a research topic of steady interest. In this work, we explore a route to the hydrodynamic limit of the RDME which uses gradient structures. Specifically, we elaborate on a method introduced in [J. Maas and A. Mielke, J. Stat. Phys., 181 (2020), pp. 2257-2303] in the context of well-mixed reaction networks by showing that, once it is complemented with an appropriate limit procedure, it can be applied to spatially extended systems with diffusion. Under the assumption of detailed balance, we write down a gradient structure for the RDME and use the method in order to produce a gradient structure for its hydrodynamic limit, namely, for the corresponding RDPDE.}, language = {en} } @article{DellnitzHuellermeierLueckeetal.2023, author = {Dellnitz, Michael and H{\"u}llermeier, Eyke and L{\"u}cke, Marvin and Ober-Bl{\"o}baum, Sina and Offen, Christian and Peitz, Sebastian and Pfannschmidt, Karlson}, title = {Efficient Time-Stepping for Numerical Integration Using Reinforcement Learning}, volume = {45}, journal = {SIAM Journal on Scientific Computing}, number = {2}, doi = {10.1137/21M1412682}, pages = {A579 -- A595}, year = {2023}, language = {en} }