TY - JOUR A1 - Götschel, Sebastian A1 - Schiela, Anton A1 - Weiser, Martin T1 - Kaskade 7 - a Flexible Finite Element Toolbox JF - Computers and Mathematics with Applications N2 - Kaskade 7 is a finite element toolbox for the solution of stationary or transient systems of partial differential equations, aimed at supporting application-oriented research in numerical analysis and scientific computing. The library is written in C++ and is based on the \textsc{Dune} interface. The code is independent of spatial dimension and works with different grid managers. An important feature is the mix-and-match approach to discretizing systems of PDEs with different ansatz and test spaces for all variables. We describe the mathematical concepts behind the library as well as its structure, illustrating its use at several examples on the way. Y1 - 2021 U6 - https://doi.org/10.1016/j.camwa.2020.02.011 VL - 81 SP - 444 EP - 458 ER - TY - JOUR A1 - Jäger, Klaus A1 - Tillmann, Peter A1 - Becker, Christiane T1 - Detailed illumination model for bifacial solar cells JF - Opt. Express Y1 - 2020 U6 - https://doi.org/10.1364/OE.383570 VL - 28 SP - 4751 ER - TY - JOUR A1 - Feis, Joshua A1 - Beutel, Dominik A1 - Köpfler, Julian A1 - Santiago, Xavier Garcia A1 - Rockstuhl, Carsten A1 - Wegener, Martin A1 - Fernandez-Corbaton, Ivan T1 - Helicity-Preserving Optical Cavity Modes for Enhanced Sensing of Chiral Molecules JF - Phys. Rev. Lett. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevLett.124.033201 VL - 124 SP - 033201 ER - TY - JOUR A1 - Werdehausen, Daniel A1 - Burger, Sven A1 - Staude, Isabelle A1 - Pertsch, Thomas A1 - Decker, Manuel T1 - Flat optics in high numerical aperture broadband imaging systems JF - J. Opt. Y1 - 2020 U6 - https://doi.org/10.1088/2040-8986/ab8ea2 VL - 22 SP - 065607 ER - TY - JOUR A1 - Farchmin, Nando A1 - Hammerschmidt, Martin A1 - Schneider, Philipp-Immanuel A1 - Wurm, Matthias A1 - Bodermann, Bernd A1 - Bär, Markus A1 - Heidenreich, Sebastian T1 - Efficient Bayesian inversion for shape reconstruction of lithography masks JF - J. Micro/Nanolith. MEMS MOEMS Y1 - 2020 U6 - https://doi.org/10.1117/1.JMM.19.2.024001 VL - 19 SP - 024001 ER - TY - JOUR A1 - Binkowski, Felix A1 - Betz, Fridtjof A1 - Colom, Rémi A1 - Hammerschmidt, Martin A1 - Zschiedrich, Lin A1 - Burger, Sven T1 - Quasinormal mode expansion of optical far-field quantities JF - Phys. Rev. B Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevB.102.035432 VL - 102 SP - 035432 ER - TY - JOUR A1 - Musiał, Anna A1 - Żołnacz, Kinga A1 - Srocka, Nicole A1 - Kravets, Oleh A1 - Große, Jan A1 - Olszewski, Jacek A1 - Poturaj, Krzysztof A1 - Wojcik, Grzegorz A1 - Mergo, Paweł A1 - Kamil, Dybka A1 - Dyrkacz, Mariusz A1 - Dlubek, Michal A1 - Lauritsen, Kristian A1 - Bülter, Andreas A1 - Schneider, Philipp-Immanuel A1 - Zschiedrich, Lin A1 - Burger, Sven A1 - Rodt, Sven A1 - Urbanczyk, Wacław A1 - Sek, Grzegorz A1 - Reitzenstein, Stephan T1 - Plug&play fibre-coupled 73 kHz single-photon source operating in the telecom O-band JF - Adv. Quantum Technol. Y1 - 2020 U6 - https://doi.org/10.1002/qute.202000018 VL - 3 SP - 2000018 ER - TY - JOUR A1 - Colom, Rémi A1 - Stout, Brian A1 - Bonod, Nicolas T1 - Modal Expansion of the Impulse Response Function for Predicting the Time Dynamics of Open Optical Cavities JF - arXiv Y1 - 2020 ER - TY - BOOK A1 - Reuter, Bernhard T1 - Generalisierte Markov-Modellierung N2 - Markov State Models (MSM) sind der Goldstandard zur Modellierung biomolekularer Dynamik, da sie die Identifizierung und Analyse metastabiler Zustände ermöglichen. Die robuste Perron-Cluster-Cluster-Analyse (PCCA+) ist ein verbreiteter Spectral-Clustering-Algorithmus, der für das Clustering hochdimensionaler MSM verwendet wird. Da die PCCA+ auf reversible Prozesse beschränkt ist, wird sie zur Generalisierten PCCA+ (G-PCCA) verallgemeinert, die geeignet ist, nichtreversible Prozesse aufzuklären. Bernhard Reuter untersucht hier mittels G-PCCA die nichtthermischen Auswirkungen von Mikrowellen auf die Proteindynamik. Dazu führt er molekulardynamische Nichtgleichgewichtssimulationen des Amyloid-β-(1–40)-Peptids durch und modelliert diese. Y1 - 2020 SN - 978-3-658-29711-4 U6 - https://doi.org/10.1007/978-3-658-29712-1 PB - Springer Spektrum, Wiesbaden ER - TY - THES A1 - Wulkow, Hanna T1 - Regularization of Elliptic Partial Differential Equations using Neural Networks N2 - This thesis presents a method for interpolating data using a neural network. The data is sparse and perturbed and is used as training data for a small neural network. For severely perturbed data, the network does not manage to find a smooth interpolation. But as the data resembles the solution to the one-dimensional and time-independent heat equation, the weak form of this PDE and subsequently its functional can be written down. If the functional is minimized, a solution to the weak form of the heat equation is found. The functional is now added to the traditional loss function of a neural network, the mean squared error between the network prediction and the given data, in order to smooth out fluctuations and interpolate between distanced grid points. This way, the network minimizes both the mean squared error and the functional, resulting in a smoother curve that can be used to predict u(x) for any grid point x. KW - elliptic partial differential equations KW - neural networks Y1 - 2020 ER - TY - GEN A1 - Bittracher, Andreas A1 - Schütte, Christof ED - Junge, Oliver ED - Schütze, O. ED - Froyland, Gary ED - Ober-Blobaum, S. ED - Padberg-Gehle, E. T1 - A weak characterization of slow variables in stochastic dynamical systems T2 - Advances in Dynamics, Optimization and Computation. Series: Studies in Systems, Decision and Control. A volume dedicated to Michael Dellnitz on the occasion of his 60th birthday Y1 - 2020 SN - 978-3-030-51264-4 U6 - https://doi.org/10.1007/978-3-030-51264-4_6 VL - 304 SP - 132 EP - 150 PB - Springer International ER - TY - GEN A1 - Mollenhauer, Mattes A1 - Schuster, Ingmar A1 - Klus, Stefan A1 - Schütte, Christof ED - Junge, Oliver ED - Schütze, O. ED - Froyland, Gary ED - Ober-Blobaum, S. ED - Padberg-Gehle, K. T1 - Singular Value Decomposition of Operators on Reproducing Kernel Hilbert Spaces T2 - Advances om Dynamics, Optimization and Computation. Series: Studies in Systems, Decision and Control. A volume dedicated to Michael Dellnitz on his 60th birthday Y1 - 2020 SN - 978-3-030-51264-4 U6 - https://doi.org/10.1007/978-3-030-51264-4_5 VL - 304 SP - 109 EP - 131 PB - Springer International ER - TY - CHAP A1 - Rodt, Sven A1 - Schneider, Philipp-Immanuel A1 - Zschiedrich, Lin A1 - Heindel, Tobias A1 - Bounouar, Samir A1 - Kantner, Markus A1 - Koprucki, Thomas A1 - Bandelow, Uwe A1 - Burger, Sven A1 - Reitzenstein, Stephan ED - Kneissl, Michael ED - Knorr, Andreas ED - Reitzenstein, Stephan ED - Hoffmann, Axel T1 - Deterministic Quantum Devices for Optical Quantum Communication T2 - Semiconductor Nanophotonics Y1 - 2020 U6 - https://doi.org/10.1007/978-3-030-35656-9_8 VL - 194 SP - 285 EP - 359 PB - Springer CY - Cham ER - TY - JOUR A1 - Gorgulla, Christoph A1 - Boeszoermnyi, Andras A1 - Wang, Zi-Fu A1 - Fischer, Patrick D. A1 - Coote, Paul A1 - Das, Krishna M. Padmanabha A1 - Malets, Yehor S. A1 - Radchenko, Dmytro S. A1 - Moroz, Yurii A1 - Scott, David A. A1 - Fackeldey, Konstantin A1 - Hoffmann, Moritz A1 - Iavniuk, Iryna A1 - Wagner, Gerhard A1 - Arthanari, Haribabu T1 - An open-source drug discovery platform enables ultra-large virtual screens JF - Nature N2 - On average, an approved drug today costs $2-3 billion and takes over ten years to develop1. In part, this is due to expensive and time-consuming wet-lab experiments, poor initial hit compounds, and the high attrition rates in the (pre-)clinical phases. Structure-based virtual screening (SBVS) has the potential to mitigate these problems. With SBVS, the quality of the hits improves with the number of compounds screened2. However, despite the fact that large compound databases exist, the ability to carry out large-scale SBVSs on computer clusters in an accessible, efficient, and flexible manner has remained elusive. Here we designed VirtualFlow, a highly automated and versatile open-source platform with perfect scaling behaviour that is able to prepare and efficiently screen ultra-large ligand libraries of compounds. VirtualFlow is able to use a variety of the most powerful docking programs. Using VirtualFlow, we have prepared the largest and freely available ready-to-dock ligand library available, with over 1.4 billion commercially available molecules. To demonstrate the power of VirtualFlow, we screened over 1 billion compounds and discovered a small molecule inhibitor (iKeap1) that engages KEAP1 with nanomolar affinity (Kd = 114 nM) and disrupts the interaction between KEAP1 and the transcription factor NRF2. We also identified a set of structurally diverse molecules that bind to KEAP1 with submicromolar affinity. This illustrates the potential of VirtualFlow to access vast regions of the chemical space and identify binders with high affinity for target proteins. Y1 - 2020 U6 - https://doi.org/https://doi.org/10.1038/s41586-020-2117-z VL - 580 SP - 663 EP - 668 PB - Springer Nature ER - TY - JOUR A1 - Helfmann, Luzie A1 - Ribera Borrell, Enric A1 - Schütte, Christof A1 - Koltai, Peter T1 - Extending Transition Path Theory: Periodically Driven and Finite-Time Dynamics JF - Journal of Nonlinear Science Y1 - 2020 U6 - https://doi.org/https://doi.org/10.1007/s00332-020-09652-7 VL - 30 SP - 3321 EP - 3366 ER - TY - JOUR A1 - Niemann, Jan-Hendrik A1 - Winkelmann, Stefanie A1 - Wolf, Sarah A1 - Schütte, Christof T1 - Agent-based modeling: Population limits and large timescales JF - Chaos: An Interdisciplinary Journal of Nonlinear Science N2 - Modeling, simulation and analysis of interacting agent systems is a broad field of research, with existing approaches reaching from informal descriptions of interaction dynamics to more formal, mathematical models. In this paper, we study agent-based models (ABMs) given as continuous-time stochastic processes and their pathwise approximation by ordinary and stochastic differential equations (ODEs and SDEs, respectively) for medium to large populations. By means of an appropriately adapted transfer operator approach we study the behavior of the ABM process on long time scales. We show that, under certain conditions, the transfer operator approach allows to bridge the gap between the pathwise results for large populations on finite timescales, i.e., the SDE limit model, and approaches built to study dynamical behavior on long time scales like large deviation theory. The latter provides a rigorous analysis of rare events including the associated asymptotic rates on timescales that scale exponentially with the population size. We demonstrate that it is possible to reveal metastable structures and timescales of rare events of the ABM process by finite-length trajectories of the SDE process for large enough populations. This approach has the potential to drastically reduce computational effort for the analysis of ABMs. Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-77309 SN - 1438-0064 VL - 31 IS - 3 ER - TY - THES A1 - Baumann, Felix T1 - Impact of Contact Surfaces on the Death Estimation N2 - The determination of time of death is one of the central tasks in forensic medicine. A standard method of time of death estimation elies on matching temperature measurements of the corpse with a post-mortem cooling model. In addition to widely used empirical post-mortem models, modelling based on a precise mathematical simulation of the cooling process have been gaining popularity. The simulation based cooling models and the resulting time of death estimates dependon a large variety of parameters. These include hermal properties for different body tissue types, environmental conditions such as temperature and air flow, and the presence of clothing and coverings. In this thesis we focus on a specific arameter - the contact between corpse and underground - and investigate its influence on the time of death estimation. Resulting we aim to answer the question whether it is necessary to consider contact mechanics in the underlying mathematical cooling model. KW - time of death estimation KW - contact mechanics KW - heat equation KW - inverse problems Y1 - 2020 ER - TY - JOUR A1 - Temp, Julia A1 - Labuz, Dominika A1 - Negrete, Roger A1 - Sunkara, Vikram A1 - Machelska, Halina T1 - Pain and knee damage in male and female mice in the medial meniscal transection-induced osteoarthritis JF - Osteoarthritis and Cartilage Y1 - 2019 U6 - https://doi.org/10.1016/j.joca.2019.11.003 ER - TY - THES A1 - Semler, Phillip T1 - Mixed finite elements in cloth simulation N2 - This master thesis investigates the use and behaviour of a mixed finite element formulation for the simulation of garments. The garment is modelled as an isotropic shell and is related to its mid-surface by energetic degeneration. Based on this, an energy functional is constructed, which contains the deformation and the mid-surface vector as degree of freedom. It is then shown why this problem does not correspond to a saddle point problem, but to a non-convex energy minimization. The implementation of the energy minimization takes place with the ZIB-internal FE framework Kaskade7.4, whereby a geometric linear and different geometric non-linear problems are examined, whereby for a selected, non-linear example a comparison is made with an existing implementation on basis of Morley elements. The further evaluations include the analysis of the quantitative and qualitative results, the used solution method, the behaviour of the system energy as well as the used CPU time. KW - cloth simulation KW - shell elements KW - mixed finite elements Y1 - 2020 ER - TY - CHAP A1 - Krause, Rolf A1 - Weiser, Martin T1 - Multilevel augmented Lagrangian solvers for overconstrained contact formulations T2 - ESAIM: ProcS N2 - Multigrid methods for two-body contact problems are mostly based on special mortar discretizations, nonlinear Gauss-Seidel solvers, and solution-adapted coarse grid spaces. Their high computational efficiency comes at the cost of a complex implementation and a nonsymmetric master-slave discretization of the nonpenetration condition. Here we investigate an alternative symmetric and overconstrained segment-to-segment contact formulation that allows for a simple implementation based on standard multigrid and a symmetric treatment of contact boundaries, but leads to nonunique multipliers. For the solution of the arising quadratic programs, we propose augmented Lagrangian multigrid with overlapping block Gauss-Seidel smoothers. Approximation and convergence properties are studied numerically at standard test problems. Y1 - 2021 U6 - https://doi.org/10.1051/proc/202171175 VL - 71 SP - 175 EP - 184 ER - TY - GEN A1 - Musiał, Anna A1 - Żołnacz, Kinga A1 - Srocka, Nicole A1 - Kravets, Oleh A1 - Große, Jan A1 - Olszewski, Jacek A1 - Poturaj, Krzysztof A1 - Wojcik, Grzegorz A1 - Mergo, Paweł A1 - Dybka, Kamil A1 - Dyrkacz, Mariusz A1 - Dlubek, Michal A1 - Lauritsen, Kristian A1 - Bülter, Andreas A1 - Schneider, Philipp-Immanuel A1 - Zschiedrich, Lin A1 - Burger, Sven A1 - Rodt, Sven A1 - Urbanczyk, Wacław A1 - Sek, Grzegorz A1 - Reitzenstein, Stephan T1 - Front Cover: Plug&Play Fiber‐Coupled 73 kHz Single‐Photon Source Operating in the Telecom O‐Band (Adv. Quantum Technol. 6/2020) T2 - Adv. Quantum Technol. Y1 - 2020 U6 - https://doi.org/10.1002/qute.202070061 VL - 3 SP - 2070061 ER - TY - JOUR A1 - Martens, Kevin A1 - Binkowski, Felix A1 - Nguyen, Linh A1 - Hu, Li A1 - Govorov, Alexander O. A1 - Burger, Sven A1 - Liedl, Tim T1 - Long- and Short-Ranged Chiral Interactions in DNA Assembled Plasmonic Chains JF - Nat. Commun. Y1 - 2021 U6 - https://doi.org/10.1038/s41467-021-22289-8 VL - 12 SP - 2025 ER - TY - JOUR A1 - Gorgulla, Christoph A1 - Fackeldey, Konstantin A1 - Wagner, Gerhard A1 - Arthanari, Haribabu T1 - Accounting of Receptor Flexibility in Ultra-Large Virtual Screens with VirtualFlow Using a Grey Wolf Optimization Method JF - Supercomputing Frontiers and Innovations N2 - Structure-based virtual screening approaches have the ability to dramatically reduce the time and costs associated to the discovery of new drug candidates. Studies have shown that the true hit rate of virtual screenings improves with the scale of the screened ligand libraries. Therefore, we have recently developed an open source drug discovery platform (VirtualFlow), which is able to routinely carry out ultra-large virtual screenings. One of the primary challenges of molecular docking is the circumstance when the protein is highly dynamic or when the structure of the protein cannot be captured by a static pose. To accommodate protein dynamics, we report the extension of VirtualFlow to allow the docking of ligands using a grey wolf optimization algorithm using the docking program GWOVina, which substantially improves the quality and efficiency of flexible receptor docking compared to AutoDock Vina. We demonstrate the linear scaling behavior of VirtualFlow utilizing GWOVina up to 128 000 CPUs. The newly supported docking method will be valuable for drug discovery projects in which protein dynamics and flexibility play a significant role. Y1 - 2020 U6 - https://doi.org/10.14529/jsfi200301 VL - 7 IS - 3 SP - 4 EP - 12 ER - TY - JOUR A1 - Kostre, Margarita A1 - Schütte, Christof A1 - Noé, Frank A1 - del Razo Sarmina, Mauricio T1 - Coupling Particle-Based Reaction-Diffusion Simulations with Reservoirs Mediated by Reaction-Diffusion PDEs JF - Multiscale Modeling & Simulation N2 - Open biochemical systems of interacting molecules are ubiquitous in life-related processes. However, established computational methodologies, like molecular dynamics, are still mostly constrained to closed systems and timescales too small to be relevant for life processes. Alternatively, particle-based reaction-diffusion models are currently the most accurate and computationally feasible approach at these scales. Their efficiency lies in modeling entire molecules as particles that can diffuse and interact with each other. In this work, we develop modeling and numerical schemes for particle-based reaction-diffusion in an open setting, where the reservoirs are mediated by reaction-diffusion PDEs. We derive two important theoretical results. The first one is the mean-field for open systems of diffusing particles; the second one is the mean-field for a particle-based reaction-diffusion system with second-order reactions. We employ these two results to develop a numerical scheme that consistently couples particle-based reaction-diffusion processes with reaction-diffusion PDEs. This allows modeling open biochemical systems in contact with reservoirs that are time-dependent and spatially inhomogeneous, as in many relevant real-world applications. Y1 - 2021 U6 - https://doi.org/10.1137/20M1352739 VL - 19 IS - 4 SP - 1659 EP - 1683 PB - Society for Industrial and Applied Mathematics ER - TY - GEN A1 - Bittracher, Andreas A1 - Schütte, Christof T1 - A probabilistic algorithm for aggregating vastly undersampled large Markov chains N2 - Model reduction of large Markov chains is an essential step in a wide array of techniques for understanding complex systems and for efficiently learning structures from high-dimensional data. We present a novel aggregation algorithm for compressing such chains that exploits a specific low-rank structure in the transition matrix which, e.g., is present in metastable systems, among others. It enables the recovery of the aggregates from a vastly undersampled transition matrix which in practical applications may gain a speedup of several orders of mag- nitude over methods that require the full transition matrix. Moreover, we show that the new technique is robust under perturbation of the transition matrix. The practical applicability of the new method is demonstrated by identifying a reduced model for the large-scale traffic flow patterns from real-world taxi trip data. T3 - ZIB-Report - 20-21 Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-78688 UR - https://opus4.kobv.de/opus4-zib/frontdoor/index/index/docId/7587 SN - 1438-0064 ER - TY - JOUR A1 - Pakhomov, Anton A1 - Löchner, Franz A1 - Zschiedrich, Lin A1 - Saravi, Sina A1 - Hammerschmidt, Martin A1 - Burger, Sven A1 - Pertsch, Thomas A1 - Setzpfand, Frank T1 - Far-field polarization signatures of surface optical nonlinearity in noncentrosymmetric semiconductors JF - Sci. Rep. Y1 - 2020 U6 - https://doi.org/10.1038/s41598-020-67186-0 VL - 10 SP - 10545 ER - TY - JOUR A1 - Werdehausen, Daniel A1 - Santiago, Xavier Garcia A1 - Burger, Sven A1 - Staude, Isabelle A1 - Pertsch, Thomas A1 - Rockstuhl, Carsten A1 - Decker, Manuel T1 - Modeling Optical Materials at the Single Scatterer Level: The Transition from Homogeneous to Heterogeneous Materials JF - Adv. Theory Simul. Y1 - 2020 U6 - https://doi.org/10.1002/adts.202000192 VL - 3 SP - 2000192 ER - TY - JOUR A1 - Boltz, Horst-Holger A1 - Sirbu, Alexei A1 - Stelzer, Nina A1 - Lohse, Martin J. A1 - Schütte, Christof A1 - Annibale, Paolo T1 - Quantitative spectroscopy of single molecule interaction times JF - Optic Letters N2 - Single molecule fluorescence tracking provides information at nanometer-scale and millisecond-temporal resolution about the dynamics and interaction of individual molecules in a biological environment. While the dynamic behavior of isolated molecules can be characterized well, the quantitative insight is more limited when interactions between two indistinguishable molecules occur. We address this aspect by developing a theoretical foundation for a spectroscopy of interaction times, i.e., the inference of interaction from imaging data. A non-trivial crossover between a power law to an exponential behavior of the distribution of the interaction times is highlighted, together with the dependence of the exponential term upon the microscopic reaction affinity. Our approach is validated with simulated and experimental datasets. Y1 - 2021 U6 - https://doi.org/10.1364/OL.413030 VL - 46 IS - 7 SP - 1538 EP - 1541 ER - TY - JOUR A1 - Villatoro, José A1 - Zühlke, Martin A1 - Riebe, Daniel A1 - Beitz, Toralf A1 - Weber, Marcus A1 - Löhmannsröben, Hans-Gerd T1 - Sub-ambient pressure IR-MALDI ion mobility spectrometer for the determination of low and high field mobilities JF - Analytical and Bioanalytical Chemistry N2 - A new ion mobility (IM) spectrometer, enabling mobility measurements in the pressure range between 5 and 500 mbar and in the reduced field strength range E/N of 5–90 Td, was developed and characterized. Reduced mobility (K0) values were studied under low E/N (constant value) as well as high E/N (deviation from low field K0) for a series of molecular ions in nitrogen. Infrared matrix-assisted laser desorption ionization (IR-MALDI) was used in two configurations: a source working at atmospheric pressure (AP) and, for the first time, an IR-MALDI source working with a liquid (aqueous) matrix at sub-ambient/reduced pressure (RP). The influence of RP on IR-MALDI was examined and new insights into the dispersion process were gained. This enabled the optimization of the IM spectrometer for best analytical performance. While ion desolvation is less efficient at RP, the transport of ions is more efficient, leading to intensity enhancement and an increased number of oligomer ions. When deciding between AP and RP IR-MALDI, a trade-off between intensity and resolving power has to be considered. Here, the low field mobility of peptide ions was first measured and compared with reference values from ESI-IM spectrometry (at AP) as well as collision cross sections obtained from molecular dynamics simulations. The second application was the determination of the reduced mobility of various substituted ammonium ions as a function of E/N in nitrogen. The mobility is constant up to a threshold at high E/N. Beyond this threshold, mobility increases were observed. This behavior can be explained by the loss of hydrated water molecules. Y1 - 2020 U6 - https://doi.org/10.1007/s00216-020-02735-0 VL - 412 SP - 5247 EP - 5260 ER - TY - JOUR A1 - Weber, Marcus A1 - Weitere Autoren, ED - Reichardt, Christine T1 - DIN SPEC 2343: Übertragung von sprachbasierten Daten zwischen Künstlichen Intelligenzen - Festlegung von Parametern und Formaten JF - Beuth Verlag N2 - Dieses Dokument legt Parameter und Formate für die Übertragung sprachbasierter Daten zwischen verschiedenen KI-Ökosystemen fest. Y1 - 2020 ER - TY - JOUR A1 - Rabben, Robert Julian A1 - Ray, Sourav A1 - Weber, Marcus T1 - ISOKANN: Invariant subspaces of Koopman operators learned by a neural network JF - The Journal of Chemical Physics N2 - 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. Y1 - 2020 U6 - https://doi.org/10.1063/5.0015132 VL - 153 IS - 11 SP - 114109 ER - TY - GEN A1 - Sunkara, Vikram A1 - Raharinirina, N. Alexia A1 - Peppert, Felix A1 - von Kleist, Max A1 - Schütte, Christof T1 - Inferring Gene Regulatory Networks from Single Cell RNA-seq Temporal Snapshot Data Requires Higher Order Moments N2 - 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. T3 - ZIB-Report - 20-25 Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-79664 SN - 1438-0064 ER - TY - JOUR A1 - Klebanov, Ilja A1 - Sprungk, Björn A1 - Sullivan, T. J. T1 - The linear conditional expectation in Hilbert space JF - Bernoulli Y1 - 2021 U6 - https://doi.org/10.3150/20-BEJ1308 VL - 27 IS - 4 SP - 2299 EP - 2299 ER - TY - JOUR A1 - Miron, Philippe A1 - Beron-Vera, Francisco A1 - Helfmann, Luzie A1 - Koltai, Péter T1 - Transition paths of marine debris and the stability of the garbage patches JF - Chaos: An Interdisciplinary Journal of Nonlinear Science N2 - 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. Y1 - 2021 U6 - https://doi.org/https://doi.org/10.1063/5.0030535 VL - 31 IS - 3 ER - TY - JOUR A1 - Mrowiński, Paweł A1 - Schnauber, Peter A1 - Kaganskiy, Arsenty A1 - Schall, Johannes A1 - Burger, Sven A1 - Rodt, Sven A1 - Reitzenstein, Stephan T1 - Directional single-photon emission from deterministic quantum dot waveguide structures JF - Phys. Status Solidi RRL Y1 - 2020 U6 - https://doi.org/10.1002/pssr.202000115 VL - 14 SP - 2000115 ER - TY - GEN A1 - Ray, Sourav A1 - Sunkara, Vikram A1 - Schütte, Christof A1 - Weber, Marcus T1 - How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs N2 - 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. T3 - ZIB-Report - 20-18 KW - Opioid, Ligand-Receptor Interaction, Binding Kinetics, Molecular Dynamics, Metadynamics, SQRA Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-78437 SN - 1438-0064 ER - TY - JOUR A1 - Shi, Liping A1 - Evlyukhin, Andrey B. A1 - Reinhardt, Carsten A1 - Babushkin, Ihar A1 - Zenin, Vladimir A. A1 - Burger, Sven A1 - Malureanu, Radu A1 - Chichkov, Boris N. A1 - Morgner, Uwe A1 - Kovacev, Milutin T1 - Progressive Self-Boosting Anapole-Enhanced Deep-Ultraviolet Third Harmonic During Few-Cycle Laser Radiation JF - ACS Photonics Y1 - 2020 U6 - https://doi.org/10.1021/acsphotonics.0c00753 VL - 7 SP - 1655 ER - TY - JOUR A1 - Fackeldey, Konstantin A1 - Röhm, Jonas A1 - Niknejad, Amir A1 - Chewle, Surahit A1 - Weber, Marcus T1 - Analyzing Raman Spectral Data without Separabiliy Assumption JF - Journal of Mathematical Chemistry N2 - 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. Y1 - 2021 U6 - https://doi.org/10.1007/s10910-020-01201-7 VL - 3 IS - 59 SP - 575 EP - 596 PB - Springer ER - TY - JOUR A1 - Siefke, Thomas A1 - Hurtado, Carol A1 - Dickmann, Johannes A1 - Dickmann, Walter A1 - Käseberg, Tim A1 - Meyer, Jan A1 - Burger, Sven A1 - Zeitner, Uwe A1 - Bodermann, Bernd A1 - Kroker, Stefanie T1 - Quasi-bound states in the continuum for deep subwavelength structural information retrieval for DUV nano-optical polarizers JF - Opt. Express Y1 - 2020 U6 - https://doi.org/10.1364/OE.396044 VL - 28 SP - 23122 ER - TY - JOUR A1 - Röhl, Susanne A1 - Weber, Marcus A1 - Fackeldey, Konstantin T1 - Computing the minimal rebinding effect for non-reversible processes JF - Multiscale Modeling and Simulation N2 - 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. Y1 - 2021 U6 - https://doi.org/https://doi.org/10.1137/20M1334966 VL - 19 IS - 1 SP - 460 EP - 477 ER - TY - JOUR A1 - Tockhorn, Philipp A1 - Sutter, Johannes A1 - Colom, Rémi A1 - Kegelmann, Lukas A1 - Al-Ashouri, Amran A1 - Roß, Marcel A1 - Jäger, Klaus A1 - Unold, Thomas A1 - Burger, Sven A1 - Albrecht, Steve A1 - Becker, Christiane T1 - Improved Quantum Efficiency by Advanced Light Management in Nanotextured Solution-Processed Perovskite Solar Cells JF - ACS Photonics Y1 - 2020 U6 - https://doi.org/10.1021/acsphotonics.0c00935 VL - 7 SP - 2589 ER - TY - JOUR A1 - Gorgulla, Christoph A1 - Das, Krishna M. Padmanabha A1 - Leigh, Kendra E A1 - Cespugli, Marco A1 - Fischer, Patrick D. A1 - Wang, Zi-Fu A1 - Tesseyre, Guilhem A1 - Pandita, Shreya A1 - Shnapir, Alex A1 - Calderaio, Anthony A1 - Hutcheson, Colin A1 - Gechev, Minko A1 - Rose, Alexander A1 - Lewis, Noam A1 - Yaffe, Erez A1 - Luxenburg, Roni A1 - Herce, Henry D. A1 - Durmaz, Vedat A1 - Halazonetis, Thanos D. A1 - Fackeldey, Konstantin A1 - Patten, Justin J. A1 - Chuprina, Alexander A1 - Dziuba, Igor A1 - Plekhova, Alla A1 - Moroz, Yurii A1 - Radchenko, Dmytro A1 - Tarkhanova, Olga A1 - Yavnyuk, Irina A1 - Gruber, Christian C. A1 - Yust, Ryan A1 - Payne, Dave A1 - Näär, Anders M. A1 - Namchuk, Mark N. A1 - Davey, Robert A. A1 - Wagner, Gerhard A1 - Kinney, Jamie A1 - Arthanari, Haribabu T1 - A Multi-Pronged Approach Targeting SARS-CoV-2 Proteins Using Ultra-Large Virtual Screening JF - iScience N2 - 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. Y1 - 2021 U6 - https://doi.org/10.26434/chemrxiv.12682316 VL - 24 IS - 2 SP - 102021 PB - CellPress ER - TY - JOUR A1 - Bonnet, Luc A1 - Akian, Jean-Luc A1 - Savin, Éric A1 - Sullivan, T. J. T1 - Adaptive reconstruction of imperfectly-observed monotone functions, with applications to uncertainty quantification JF - Algorithms N2 - 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. Y1 - 2020 U6 - https://doi.org/10.3390/a13080196 VL - 13 IS - 8 SP - 196 ER - TY - JOUR A1 - Würth, Christian A1 - Manley, Phillip A1 - Voigt, Robert A1 - Ahiboz, Doguscan A1 - Becker, Christiane A1 - Resch-Genger, Ute T1 - Metasurface enhanced sensitized photon upconversion: towards highly efficient low power upconversion applications and nano-scale E-field sensors JF - Nano Lett. Y1 - 2020 U6 - https://doi.org/10.1021/acs.nanolett.0c02548 VL - 20 SP - 6682 ER - TY - JOUR A1 - Alhaddad, Samer A1 - Förstner, Jens A1 - Groth, Stefan A1 - Grünewald, Daniel A1 - Grynko, Yevgen A1 - Hannig, Frank A1 - Kenter, Tobias A1 - Pfreundt, Franz-Josef A1 - Plessl, Christian A1 - Schotte, Merlind A1 - Steinke, Thomas A1 - Teich, Jürgen A1 - Weiser, Martin A1 - Wende, Florian T1 - HighPerMeshes - A Domain-Specific Language for Numerical Algorithms on Unstructured Grids JF - Euro-Par 2020: Parallel Processing Workshops. N2 - 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. Y1 - 2021 U6 - https://doi.org/10.1007/978-3-030-71593-9_15 SP - 185 EP - 196 PB - Springer ER - TY - CHAP A1 - Binkowski, Felix A1 - Betz, Fridtjof A1 - Colom, Rémi A1 - Hammerschmidt, Martin A1 - Zschiedrich, Lin A1 - Burger, Sven T1 - Modal expansion of optical far-field quantities using quasinormal modes T2 - EPJ Web Conf. Y1 - 2020 U6 - https://doi.org/10.1051/epjconf/202023805007 VL - 238 SP - 05007 ER - TY - CHAP A1 - Colom, Rémi A1 - Binkowski, Felix A1 - Betz, Fridtjof A1 - Hammerschmidt, Martin A1 - Zschiedrich, Lin A1 - Burger, Sven T1 - Quasi-normal mode expansion as a tool for the design of nanophotonic devices T2 - EPJ Web Conf. Y1 - 2020 U6 - https://doi.org/10.1051/epjconf/202023805008 VL - 238 SP - 05008 ER - TY - JOUR A1 - Ray, Sourav A1 - Sunkara, Vikram A1 - Schütte, Christof A1 - Weber, Marcus T1 - How to calculate pH-dependent binding rates for receptor-ligand systems based on thermodynamic simulations with different binding motifs JF - Molecular Simulation N2 - 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. Y1 - 2020 U6 - https://doi.org/10.1080/08927022.2020.1839660 VL - 46 IS - 18 SP - 1443 EP - 1452 PB - Taylor and Francis ER - TY - THES A1 - Betz, Fridtjof T1 - Investigation and application of contour integral methods for resonance expansion in nanophotonics Y1 - 2020 ER - TY - JOUR A1 - Santiago, Xavier Garcia A1 - Burger, Sven A1 - Rockstuhl, Carsten A1 - Schneider, Philipp-Immanuel T1 - Bayesian optimization with improved scalability and derivative information for efficient design of nanophotonic structures JF - J. Light. Technol. Y1 - 2021 U6 - https://doi.org/10.1109/JLT.2020.3023450 VL - 39 SP - 167 ER - TY - JOUR A1 - Weber, Marie-Christin A1 - Fischer, Lisa A1 - Damerau, Alexandra A1 - Ponomarev, Igor A1 - Pfeiffenberger, Moritz A1 - Gaber, Timo A1 - Götschel, Sebastian A1 - Lang, Jens A1 - Röblitz, Susanna A1 - Buttgereit, Frank A1 - Ehrig, Rainald A1 - Lang, Annemarie T1 - Macroscale mesenchymal condensation to study cytokine-driven cellular and matrix-related changes during cartilage degradation JF - Biofabrication N2 - 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. Y1 - 2020 U6 - https://doi.org/10.1088/1758-5090/aba08f VL - 12 IS - 4 ER - TY - JOUR A1 - Alves, Sónia A1 - Ehrig, Rainald A1 - Raffalt, Peter C. A1 - Bender, Alwina A1 - Duda, Georg A1 - Agres, Alison N. T1 - Quantifying Asymmetry in Gait: A Weighted Universal Symmetry Index to Evaluate 3D Ground Reaction Forces JF - Frontiers in Bioengineering and Biotechnology N2 - 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. Y1 - 2020 ER - TY - JOUR A1 - Wulkow, Niklas A1 - Koltai, Péter A1 - Schütte, Christof T1 - Memory-Based Reduced Modelling and Data-Based Estimation of Opinion Spreading JF - Journal of Nonlinear Science N2 - 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. Y1 - 2021 U6 - https://doi.org/10.1007/s00332-020-09673-2 VL - 31 ER - TY - JOUR A1 - Straube, Arthur A1 - Kowalik, Bartosz G. A1 - Netz, Roland R. A1 - Höfling, Felix T1 - Rapid onset of molecular friction in liquids bridging between the atomistic and hydrodynamic pictures JF - Commun. Phys. N2 - 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. Y1 - 2020 U6 - https://doi.org/10.1038/s42005-020-0389-0 VL - 3 SP - 126 PB - Nature ER - TY - JOUR A1 - Fischer, Sophie A1 - Ehrig, Rainald A1 - Schäfer, Stefan A1 - Tronci, Enrico A1 - Mancini, Toni A1 - Egli, Marcel A1 - Ille, Fabian A1 - Krüger, Tillmann H. C. A1 - Leeners, Brigitte A1 - Röblitz, Susanna T1 - Mathematical Modelling and Simulation Provides Evidence for New Strategies of Ovarian Stimulation JF - Frontiers in Endocrinology N2 - 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. Y1 - 2021 U6 - https://doi.org/10.3389/fendo.2021.613048 VL - 12 ER - TY - JOUR A1 - Fackeldey, Konstantin A1 - Oster, Mathias A1 - Sallandt, Leon A1 - Schneider, Reinhold T1 - Approximative Policy Iteration for Exit Time Feedback Control Problems driven by Stochastic Differential Equations using Tensor Train format JF - SIAM Journal on Multiscale Modeling and Simulation N2 - 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. Y1 - 2022 U6 - https://doi.org/10.1137/20M1372500 VL - 20 IS - 1 SP - 379 EP - 403 ER - TY - JOUR A1 - Helfmann, Luzie A1 - Djurdjevac Conrad, Natasa A1 - Djurdjevac, Ana A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - From interacting agents to density-based modeling with stochastic PDEs JF - Communications in Applied Mathematics and Computational Science N2 - 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. Y1 - 2021 U6 - https://doi.org/10.2140/camcos.2021.16.1 VL - 16 IS - 1 SP - 1 EP - 32 ER - TY - JOUR A1 - Chewle, Surahit A1 - Emmerling, Franziska A1 - Weber, Marcus T1 - Effect of choice of solvent on crystallization pathway of Paracetamol: An experimental and theoretical case study JF - Crystals N2 - 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. Y1 - 2020 U6 - https://doi.org/10.3390/cryst10121107 VL - 10 IS - 12 SP - 1107 ER - TY - JOUR A1 - Niemann, Jan-Hendrik A1 - Klus, Stefan A1 - Schütte, Christof T1 - Data-driven model reduction of agent-based systems using the Koopman generator JF - PLOS ONE N2 - 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. Y1 - 2021 U6 - https://doi.org/10.1371/journal.pone.0250970 VL - 16 IS - 5 ER - TY - GEN A1 - Becker, Christiane A1 - Jäger, Klaus A1 - Burger, Sven T1 - Nanophotonics for Solar Energy T2 - Sol. Energy Mater. Sol. Cells Y1 - 2021 U6 - https://doi.org/10.1016/j.solmat.2020.110916 VL - 221 SP - 110916 ER - TY - JOUR A1 - Klebanov, Ilja A1 - Schuster, Ingmar T1 - Markov Chain Importance Sampling - a highly efficient estimator for MCMC JF - Journal of Computational and Graphical Statistics N2 - 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. Y1 - 2020 U6 - https://doi.org/10.1080/10618600.2020.1826953 ER - TY - JOUR A1 - Scott, Philip A1 - Santiago, Xavier Garcia A1 - Beutel, Dominik A1 - Rockstuhl, Carsten A1 - Wegener, Martin A1 - Fernandez-Corbaton, Ivan T1 - On enhanced sensing of chiral molecules in optical cavities JF - Appl. Phys. Rev. Y1 - 2020 U6 - https://doi.org/10.1063/5.0025006 VL - 7 SP - 041413 ER - TY - JOUR A1 - Schneck, Jakob A1 - Weiser, Martin A1 - Wende, Florian T1 - Impact of mixed precision and storage layout on additive Schwarz smoothers JF - Numerical Linear Algebra with Applications N2 - 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. Y1 - 2021 U6 - https://doi.org/10.1002/nla.2366 VL - 28 IS - 4 ER - TY - JOUR A1 - Wulkow, Hanna A1 - Conrad, Tim A1 - Djurdjevac Conrad, Natasa A1 - Müller, Sebastian A. A1 - Nagel, Kai A1 - Schütte, Christof T1 - Prediction of Covid-19 spreading and optimal coordination of counter-measures: From microscopic to macroscopic models to Pareto fronts JF - PLOS One Y1 - 2021 U6 - https://doi.org/10.1371/journal.pone.0249676 VL - 16 IS - 4 PB - Public Library of Science ER - TY - JOUR A1 - Zimper, Sebastian A1 - Djurdjevac, Ana A1 - Hartmann, Carsten A1 - Schütte, Christof A1 - Conrad, Natasa Djurdjevac T1 - Mean-field optimal control with stochastic leaders N2 - 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ô-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. Y1 - 2025 ER - TY - JOUR A1 - Hartmann, Carsten A1 - Jöster, Annika A1 - Schütte, Christof A1 - Sikorski, Alexander A1 - Weber, Marcus T1 - Importance sampling of unbounded random stopping times: computing committor functions and exit rates without reweighting N2 - 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. Y1 - 2026 ER - TY - JOUR A1 - Garcia Santiago, Xavier A1 - Hammerschmidt, Martin A1 - Sachs, Johannes A1 - Burger, Sven A1 - Kwon, Hyunah A1 - Knöller, Marvin A1 - Arens, Tilo A1 - Fischer, Peer A1 - Fernandez-Corbaton, Ivan A1 - Rockstuhl, Carsten T1 - Toward maximally electromagnetically chiral scatterers at optical frequencies JF - ACS Photonics Y1 - 2022 U6 - https://doi.org/10.1021/acsphotonics.1c01887 VL - 9 SP - 1954 ER - TY - GEN A1 - Ernst, Ariane A1 - Schütte, Christof A1 - Sigrist, Stephan A1 - Winkelmann, Stefanie T1 - Variance of filtered signals: Characterization for linear reaction networks and application to neurotransmission dynamics N2 - 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. T3 - ZIB-Report - 21-15 KW - linear reaction networks KW - cross-correlation KW - neurotransmission Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-82674 SN - 1438-0064 ER - TY - JOUR A1 - Ernst, Ariane A1 - Schütte, Christof A1 - Sigrist, Stephan A1 - Winkelmann, Stefanie T1 - Variance of filtered signals: Characterization for linear reaction networks and application to neurotransmission dynamics JF - Mathematical Biosciences N2 - 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. Y1 - 2022 U6 - https://doi.org/10.1016/j.mbs.2021.108760 VL - 343 ER - TY - JOUR A1 - Thies, Arne A1 - Sunkara, Vikram A1 - Ray, Sourav A1 - Wulkow, Hanna A1 - Celik, M. Özgür A1 - Yergöz, Fatih A1 - Schütte, Christof A1 - Stein, Christoph A1 - Weber, Marcus A1 - Winkelmann, Stefanie T1 - Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design JF - Scientific Reports N2 - 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. Y1 - 2023 U6 - https://doi.org/10.1038/s41598-023-27699-w VL - 13 IS - 607 ER - TY - GEN A1 - Ray, Sourav A1 - Thies, Arne A1 - Sunkara, Vikram A1 - Wulkow, Hanna A1 - Celik, Özgür A1 - Yergöz, Fatih A1 - Schütte, Christof A1 - Stein, Christoph A1 - Weber, Marcus A1 - Winkelmann, Stefanie T1 - Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design N2 - 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. T3 - ZIB-Report - 21-19 Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-82797 SN - 1438-0064 ER - TY - THES A1 - Pakhomov, Anton T1 - Efficient modeling and optimization of surface second-harmonic generation from nanophotonic components Y1 - 2021 U6 - https://doi.org/10.22032/dbt.49369 ER - TY - JOUR A1 - Prasad, Anup K. A1 - Tiwari, Chandni A1 - Ray, Sourav A1 - Holden, Stephanie A1 - Armstrong, David A. A1 - Rosengren, K. Johan A1 - Rodger, Alison A1 - Panwar, Ajay S. A1 - Martin, Lisandra L. T1 - Secondary Structure Transitions for a Family of Amyloidogenic, Antimircobial Uperin 3 Peptides in Contact with Sodium Dodecyl Sulfate JF - ChemPlusChem Y1 - 2022 U6 - https://doi.org/10.1002/cplu.202100408 VL - 87 SP - e202100408 ER - TY - JOUR A1 - Berg, Mascha A1 - Plöntzke, Julia A1 - Siebert, Heike A1 - Röblitz, Susanna T1 - Modelling Oscillatory Patterns in the Bovine Estrous Cycle with Boolean Delay Equations JF - Bulletin of Mathematical Biology N2 - 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. Y1 - 2021 U6 - https://doi.org/10.1007/s11538-021-00942-z VL - 83 IS - 121 SP - 1 EP - 25 ER - TY - JOUR A1 - Djurdjevac, Ana T1 - Linear Parabolic Problems in Random Moving Domains N2 - 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. Y1 - 2021 U6 - https://doi.org/https://doi.org/10.1137/19M1284889 VL - 9 IS - 2 SP - 848 EP - 879 ER - TY - JOUR A1 - Djurdjevac, Ana A1 - Gräser, Carsten A1 - Herbert, Philip J. T1 - Non-stationary incompressible linear fluid equations in a moving domain JF - arXiv N2 - 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. Y1 - 2021 ER - TY - JOUR A1 - Colom, Remi A1 - Mikheeva, Elena A1 - Achouri, Karim A1 - Zuniga-Perez, Jesus A1 - Bonod, Nicolas A1 - Martin, Olivier J. F. A1 - Burger, Sven A1 - Genevet, Patrice T1 - Crossing of the branch cut: the topological origin of a universal 2π-phase retardation in non-Hermitian metasurfaces JF - Laser Photonics Rev. Y1 - 2023 U6 - https://doi.org/10.1002/lpor.202200976 VL - 17 SP - 2200976 ER - TY - JOUR A1 - Montefusco, Alberto A1 - Schütte, Christof A1 - Winkelmann, Stefanie T1 - A route to the hydrodynamic limit of a reaction-diffusion master equation using gradient structures JF - SIAM Journal on Applied Mathematics N2 - 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. Y1 - 2023 U6 - https://doi.org/10.1137/22M1488831 VL - 83 IS - 2 SP - 837 EP - 861 ER - TY - JOUR A1 - Dellnitz, Michael A1 - Hüllermeier, Eyke A1 - Lücke, Marvin A1 - Ober-Blöbaum, Sina A1 - Offen, Christian A1 - Peitz, Sebastian A1 - Pfannschmidt, Karlson T1 - Efficient Time-Stepping for Numerical Integration Using Reinforcement Learning JF - SIAM Journal on Scientific Computing Y1 - 2023 U6 - https://doi.org/10.1137/21M1412682 VL - 45 IS - 2 SP - A579 EP - A595 ER - TY - JOUR A1 - Peppert, Felix A1 - von Kleist, Max A1 - Schütte, Christof A1 - Sunkara, Vikram T1 - On the Sufficient Condition for Solving the Gap-Filling Problem Using Deep Convolutional Neural Networks JF - IEEE Transactions on Neural Networks and Learning Systems N2 - Deep convolutional neural networks (DCNNs) are routinely used for image segmentation of biomedical data sets to obtain quantitative measurements of cellular structures like tissues. These cellular structures often contain gaps in their boundaries, leading to poor segmentation performance when using DCNNs like the U-Net. The gaps can usually be corrected by post-hoc computer vision (CV) steps, which are specific to the data set and require a disproportionate amount of work. As DCNNs are Universal Function Approximators, it is conceivable that the corrections should be obsolete by selecting the appropriate architecture for the DCNN. In this article, we present a novel theoretical framework for the gap-filling problem in DCNNs that allows the selection of architecture to circumvent the CV steps. Combining information-theoretic measures of the data set with a fundamental property of DCNNs, the size of their receptive field, allows us to formulate statements about the solvability of the gap-filling problem independent of the specifics of model training. In particular, we obtain mathematical proof showing that the maximum proficiency of filling a gap by a DCNN is achieved if its receptive field is larger than the gap length. We then demonstrate the consequence of this result using numerical experiments on a synthetic and real data set and compare the gap-filling ability of the ubiquitous U-Net architecture with variable depths. Our code is available at https://github.com/ai-biology/dcnn-gap-filling. Y1 - 2022 U6 - https://doi.org/10.1109/TNNLS.2021.3072746 VL - 33 IS - 11 SP - 6194 EP - 6205 ER - TY - JOUR A1 - Boeth, Heide A1 - Biesen, Robert A1 - Hollnagel, Jeremias A1 - Herrmann, Sandra A1 - Ehrig, Rainald A1 - Pelli, Lorenzo A1 - Taylor, William R. A1 - Duda, Georg A1 - Buttgereit, Frank T1 - Quantification of morning stiffness to assess disease activity and treatment effects in rheumatoid arthritis JF - Rheumatology N2 - Objectives The clinical parameter “morning stiffness” is widely used to assess the status of rheumatoid arthritis (RA), but its accurate quantitative assessment in a clinical setting has not yet been successful. This lack of individual quantification limits both personalized medication and efficacy evaluation in the treatment of RA. Methods We have developed a novel technology to assess passive resistance of the metacarpophalangeal (MCP) III joint (stiffness) and its Passive Range of Motion (PRoM). Within this pilot study, nineteen female postmenopausal RA patients and nine healthy controls were examined in the evening as well as in the morning of the following day. To verify the specificity of the biomechanical quantification, eleven patients with RA were assessed both prior to and ∼3 h after glucocorticoid therapy. Results While the healthy controls showed only minor changes between afternoon and morning, in RA patients mean±SD PRoM decreased significantly by 18 ± 22% and stiffness increased significantly by 20 ± 18% in the morning compared with the previous afternoon. We found a significant positive correlation between RA activity and biomechanical measures. Glucocorticoids significantly increased mean PRoM by 16 ± 11% and reduced mean stiffness by 23 ± 22%. Conclusion This technology allowed mechanical stiffness to be quantified in MCP joints, and has demonstrated high sensitivity in respect to disease status as well as medication effect in RA patients. Such non-invasive, low risk, and rapid assessment of biomechanical joint stiffness opens a novel avenue for judging therapy efficacy in patients with RA, and potentially also in other non-RA inflammatory joint diseases. Y1 - 2021 U6 - https://doi.org/10.1093/rheumatology/keab323 VL - 60 IS - 11 SP - 5282 EP - 5291 ER - TY - GEN A1 - Pflüger, Mika A1 - Kline, R Joseph A1 - Fernández Herrero, Analía A1 - Hammerschmidt, Martin A1 - Soltwisch, Victor A1 - Krumrey, Michael T1 - Extracting dimensional parameters of gratings produced with self-aligned multiple patterning using grazing-incidence small-angle x-ray scattering [Source Code] Y1 - 2020 U6 - https://doi.org/10.24433/CO.0953516.v2 ER - TY - JOUR A1 - Gao, Timm A1 - Rickert, Lucas A1 - Urban, Felix A1 - Große, Jan A1 - Srocka, Nicole A1 - Rodt, Sven A1 - Musial, Anna A1 - Zolnacz, Kinga A1 - Mergo, Pawel A1 - Dybka, Kamil A1 - Urbanczyk, Waclaw A1 - Sek, Grzegorz A1 - Burger, Sven A1 - Reitzenstein, Stephan A1 - Heindel, Tobias T1 - A quantum key distribution testbed using a plug&play telecom-wavelength single-photon source JF - Appl. Phys. Rev. Y1 - 2022 U6 - https://doi.org/10.1063/5.0070966 VL - 9 SP - 011412 ER - TY - JOUR A1 - Torun, Cem Güney A1 - Schneider, Philipp-Immanuel A1 - Hammerschmidt, Martin A1 - Burger, Sven A1 - Munns, Joseph H.D. A1 - Schröder, Tim T1 - Optimized diamond inverted nanocones for enhanced color center to fiber coupling JF - Appl. Phys. Lett. Y1 - 2021 U6 - https://doi.org/10.1063/5.0050338 VL - 118 SP - 234002 ER - TY - JOUR A1 - Gorgulla, Christoph A1 - Çınaroğlu, Süleyman A1 - Fischer, Patrick D. A1 - Fackeldey, Konstantin A1 - Wagner, Gerhard A1 - Arthanari, Haribabu T1 - VirtualFlow Ants—Ultra-Large Virtual Screenings with Artificial Intelligence Driven Docking Algorithm Based on Ant Colony Optimization JF - Special Issue Artificial Intelligence & Deep Learning Approaches for Structural Bioinformatics N2 - The docking program PLANTS, which is based on ant colony optimization (ACO) algorithm, has many advanced features for molecular docking. Among them are multiple scoring functions, the possibility to model explicit displaceable water molecules, and the inclusion of experimental constraints. Here, we add support of PLANTS to VirtualFlow (VirtualFlow Ants), which adds a valuable method for primary virtual screenings and rescoring procedures. Furthermore, we have added support of ligand libraries in the MOL2 format, as well as on the fly conversion of ligand libraries which are in the PDBQT format to the MOL2 format to endow VirtualFlow Ants with an increased flexibility regarding the ligand libraries. The on the fly conversion is carried out with Open Babel and the program SPORES. We applied VirtualFlow Ants to a test system involving KEAP1 on the Google Cloud up to 128,000 CPUs, and the observed scaling behavior is approximately linear. Furthermore, we have adjusted several central docking parameters of PLANTS (such as the speed parameter or the number of ants) and screened 10 million compounds for each of the 10 resulting docking scenarios. We analyzed their docking scores and average docking times, which are key factors in virtual screenings. The possibility of carrying out ultra-large virtual screening with PLANTS via VirtualFlow Ants opens new avenues in computational drug discovery. Y1 - 2021 U6 - https://doi.org/https://doi.org/10.3390/ijms22115807 VL - 22 IS - 11 SP - 5807 ER - TY - JOUR A1 - Ahiboz, Doguscan A1 - Andresen, Elina A1 - Manley, Phillip A1 - Resch-Genger, Ute A1 - Würth, Christian A1 - Becker, Christiane T1 - Metasurface-Enhanced Photon Upconversion upon 1550 nm Excitation JF - Adv. Opt. Mater. Y1 - 2021 U6 - https://doi.org/10.1002/adom.202101285 VL - 9 SP - 2101285 ER - TY - JOUR A1 - Martens, Kevin A1 - Funck, Timon A1 - Santiago, Eva Y. A1 - Govorov, Alexander O. A1 - Burger, Sven A1 - Liedl, Tim T1 - Onset of Chirality in Plasmonic Meta-Molecules and Dielectric Coupling JF - ACS Nano Y1 - 2022 U6 - https://doi.org/10.1021/acsnano.2c04729 VL - 16 SP - 16143 ER - TY - CHAP A1 - Perevoznik, Dmitrii A1 - Bose, Surajit A1 - Burger, Sven A1 - Demircan, Ayhan A1 - Morgner, Uwe T1 - Photonic components in polymers made by femtosecond pulses T2 - Conference on Lasers and Electro-Optics Europe & European Quantum Electronics Conference (CLEO/Europe-EQEC) Y1 - 2021 U6 - https://doi.org/10.1109/CLEO/Europe-EQEC52157.2021.9541906 SP - 1 ER - TY - CHAP A1 - Ahiboz, Doguscan A1 - Andresen, Elina A1 - Manley, Phillip A1 - Resch Genger, Ute A1 - Würth, Christian A1 - Becker, Christiane T1 - Enhanced Photon Upconversion Using Erbium-Doped Nanoparticles Interacting with Silicon Metasurfaces T2 - IEEE 48th Photovoltaic Specialists Conference (PVSC) Y1 - 2021 U6 - https://doi.org/10.1109/PVSC43889.2021.9518495 SP - 2008 EP - 2010 ER - TY - JOUR A1 - Bittracher, Andreas A1 - Mollenhauer, Mattes A1 - Koltai, Péter A1 - Schütte, Christof T1 - Optimal Reaction Coordinates: Variational Characterization and Sparse Computation JF - Multiscale Modelling & Simulation N2 - Reaction coordinates (RCs) are indicators of hidden, low-dimensional mechanisms that govern the long-term behavior of high-dimensional stochastic processes. We present a novel and general variational characterization of optimal RCs and provide conditions for their existence. Optimal RCs are minimizers of a certain loss function, and reduced models based on them guarantee a good approximation of the statistical long-term properties of the original high-dimensional process. We show that for slow-fast systems, metastable systems, and other systems with known good RCs, the novel theory reproduces previous insight. Remarkably, for reversible systems, the numerical effort required to evaluate the loss function scales only with the variability of the underlying, low-dimensional mechanism, and not with that of the full system. The theory provided lays the foundation for an efficient and data-sparse computation of RCs via modern machine learning techniques. Y1 - 2023 U6 - https://doi.org/10.1137/21M1448367 VL - 21 IS - 2 SP - 449 EP - 488 ER - TY - JOUR A1 - Mollenhauer, Mattes A1 - Klus, Stefan A1 - Schütte, Christof A1 - Koltai, Péter T1 - Kernel Autocovariance Operators of Stationary Processes: Estimation and Convergence JF - Journal of Machine Learning Research N2 - We consider autocovariance operators of a stationary stochastic process on a Polish space that is embedded into a reproducing kernel Hilbert space. We investigate how empirical estimates of these operators converge along realizations of the process under various conditions. In particular, we examine ergodic and strongly mixing processes and obtain several asymptotic results as well as finite sample error bounds. We provide applications of our theory in terms of consistency results for kernel PCA with dependent data and the conditional mean embedding of transition probabilities. Finally, we use our approach to examine the nonparametric estimation of Markov transition operators and highlight how our theory can give a consistency analysis for a large family of spectral analysis methods including kernel-based dynamic mode decomposition. Y1 - 2022 UR - https://jmlr.org/papers/v23/20-442.html VL - 23 IS - 327 SP - 1 EP - 34 ER - TY - CHAP A1 - Rickert, Lucas A1 - Gao, Timm A1 - Urban, Felix A1 - Große, Jan A1 - Srocka, Nicole A1 - Rodt, Sven A1 - Musial, Anna A1 - Zolnacz, Kinga A1 - Mergo, Pawel A1 - Dybka, Kamil A1 - Urbanczyk, Waclaw A1 - Sek, Grzegorz A1 - Burger, Sven A1 - Reitzenstein, Stephan A1 - Tobias, Heindel T1 - A Quantum Key Distribution Testbed Using Plug&Play Telecom-Wavelength Single-Photons T2 - Quantum Information and Measurement VI Y1 - 2021 U6 - https://doi.org/10.1364/QIM.2021.M2C.1 SP - M2C.1 ER - TY - JOUR A1 - Heida, Martin A1 - Sikorski, Alexander A1 - Weber, Marcus T1 - Consistency and order 1 convergence of cell-centered finite volume discretizations of degenerate elliptic problems in any space dimension JF - SIAM Journal on Numerical Analysis N2 - We study consistency of cell-centered finite difference methods for elliptic equations with degenerate coefficients in any space dimension $d \geq 2$. This results in order of convergence estimates in the natural weighted energy norm and in the weighted discrete $L^2$-norm on admissible meshes. The cells of meshes under consideration may be very irregular in size. We particularly allow the size of certain cells to remain bounded from below even in the asymptotic limit. For uniform meshes we show that the order of convergence is at least 1 in the energy semi-norm, provided the discrete and continuous solutions exist and the continuous solution has $H^2$ regularity. Y1 - 2022 U6 - https://doi.org/10.20347/WIAS.PREPRINT.2913 ER - TY - CHAP A1 - Schneider, Philipp-Immanuel A1 - Garcia-Santiago, Xavier A1 - Wiegand, Benjamin A1 - Anton, Oliver A1 - Krutzik, Markus A1 - Rockstuhl, Carsten A1 - Burger, Sven T1 - Advances in Bayesian optimization for photonics and quantum atom optics applications T2 - OSA Advanced Photonics Congress Y1 - 2021 U6 - https://doi.org/10.1364/IPRSN.2021.JTh1E.2 SP - JTh1E.2 ER - TY - JOUR A1 - Colom, Remi A1 - Binkowski, Felix A1 - Betz, Fridtjof A1 - Kivshar, Yuri A1 - Burger, Sven T1 - Enhanced Purcell factor for nanoantennas supporting interfering resonances JF - Phys. Rev. Research Y1 - 2022 U6 - https://doi.org/10.1103/PhysRevResearch.4.023189 VL - 4 SP - 023189 ER - TY - JOUR A1 - Lücke, Marvin A1 - Nüske, Feliks T1 - tgEDMD: Approximation of the Kolmogorov Operator in Tensor Train Format JF - Journal of Nonlinear Science N2 - Extracting information about dynamical systems from models learned off simulation data has become an increasingly important research topic in the natural and engineering sciences. Modeling the Koopman operator semigroup has played a central role in this context. As the approximation quality of any such model critically depends on the basis set, recent work has focused on deriving data-efficient representations of the Koopman operator in low-rank tensor formats, enabling the use of powerful model classes while avoiding over-fitting. On the other hand, detailed information about the system at hand can be extracted from models for the infinitesimal generator, also called Kolmogorov backward operator for stochastic differential equations. In this work, we present a data-driven method to efficiently approximate the generator using the tensor train (TT) format. The centerpiece of the method is a TT representation of the tensor of generator evaluations at all data sites. We analyze consistency and complexity of the method, present extensions to practically relevant settings, and demonstrate its applicability to benchmark numerical examples. Y1 - 2022 U6 - https://doi.org/10.1007/s00332-022-09801-0 VL - 32 ER - TY - CHAP A1 - Musial, Anna A1 - Zolnacz, Kinga A1 - Srocka, Nicole A1 - Kravets, Oleh A1 - Große, Jan A1 - Schneider, Philipp-Immanuel A1 - Olszewski, Jacek A1 - Poturaj, Krzystof A1 - Wojcik, Grzegorz A1 - Mergo, Pawel A1 - Dybka, Kamil A1 - Dyrkacz, Mariusz A1 - Dlubek, Michal A1 - Rodt, Sven A1 - Burger, Sven A1 - Zschiedrich, Lin A1 - Urbanczyk, Waclaw A1 - Sek, Grzegorz A1 - Reitzenstein, Stephan T1 - Stand-alone quantum dot-based single-photon source operating at telecommunication wavelengths T2 - 10th International Conference on Spontaneous Coherence in Excitonic Systems ICSCE10 Y1 - 2020 SP - 39 ER - TY - JOUR A1 - Fackeldey, Konstantin A1 - Gorgulla, Christoph A1 - Weber, Marcus T1 - Neue Medikamente dank Supercomputern JF - Spektrum der Wissenschaft N2 - Die aktuelle Pandemie verdeutlicht, wie wichtig es ist, rasch geeignete Arzneimittel zu finden. In Computer­simulationen gelingt das erheblich schneller als im Labor. Gegen das Coronavirus ließen sich auf diese Weise bereits Wirkstoffkandidaten identifizieren. Y1 - 2021 IS - 11 SP - 40 EP - 46 ER - TY - JOUR A1 - Stout, Brian A1 - Colom, Remi A1 - Bonod, Nicolas A1 - McPhedran, Ross C T1 - Spectral expansions of open and dispersive optical systems: Gaussian regularization and convergence JF - New J. Phys. Y1 - 2021 U6 - https://doi.org/10.1088/1367-2630/ac10a6 VL - 23 SP - 083004 ER - TY - JOUR A1 - Hungenberg, Klaus-Dieter A1 - Schwede, Christian A1 - Wulkow, Michael A1 - Wulkow, Niklas T1 - Determination of reactivity ratios for acrylic acid and its dimer from classical parameter estimation and Bayesian approach JF - Canadian Journal of Chemical Engineering N2 - The reactivity ratios of acrylic acid (AA, M1) and its dimer beta-acroyloxypropionic acid (diAA, M2) are determined from cumulative copolymerization data by two different methods: classical parameter estimation (PE) by minimizing the objective function and a Bayesian analysis. Classical PE gives r1 =0.74 and r2 = 1.23 at the minimum of the residual. From the Bayesian analysis, the probability distribution of the parameter sets is obtained, revealing the existence of parameter sets with rather the same probability. The influence of the number of data and the size of the measurement error are discussed. Y1 - 2021 U6 - https://doi.org/https://doi.org/10.1002/cjce.24330 SP - 1 EP - 10 ER -