TY - CHAP A1 - Binkowski, Felix A1 - Betz, Fridtjof A1 - Colom, Remi A1 - Hammerschmidt, Martin A1 - Schneider, Philipp-Immanuel A1 - Zschiedrich, Lin A1 - Burger, Sven ED - Zouhdi, S. ED - Pawlak, D. ED - Kudelski, A. T1 - Contour integral methods for resonance phenomena in nano-optics T2 - The 11th International Conference on Metamaterials, Photonic Crystals and Plasmonics (META 2021) Y1 - 2021 SN - 2429-1390 SP - 523 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - JOUR A1 - Straube, Arthur A1 - Winkelmann, Stefanie A1 - Schütte, Christof A1 - Höfling, Felix T1 - Stochastic pH oscillations in a model of the urea–urease reaction confined to lipid vesicles JF - J. Phys. Chem. Lett. N2 - The urea-urease clock reaction is a pH switch from acid to basic that can turn into a pH oscillator if it occurs inside a suitable open reactor. We numerically study the confinement of the reaction to lipid vesicles, which permit the exchange with an external reservoir by differential transport, enabling the recovery of the pH level and yielding a constant supply of urea molecules. For microscopically small vesicles, the discreteness of the number of molecules requires a stochastic treatment of the reaction dynamics. Our analysis shows that intrinsic noise induces a significant statistical variation of the oscillation period, which increases as the vesicles become smaller. The mean period, however, is found to be remarkably robust for vesicle sizes down to approximately 200 nm, but the periodicity of the rhythm is gradually destroyed for smaller vesicles. The observed oscillations are explained as a canard-like limit cycle that differs from the wide class of conventional feedback oscillators. Y1 - 2021 U6 - https://doi.org/10.1021/acs.jpclett.1c03016 VL - 12 SP - 9888 EP - 9893 ER - TY - CHAP A1 - Tillmann, Peter A1 - Jäger, Klaus A1 - Katz, Eugene A1 - Becker, Christiane T1 - Relaxed current-matching constraints by bifacial operation and luminescent coupling in perovskite/silicon tandem solar cells T2 - 2021 IEEE 48th Photovoltaic Specialists Conference (PVSC) Y1 - 2021 U6 - https://doi.org/10.1109/PVSC43889.2021.9518891 SP - 1086 ER - TY - JOUR A1 - Tierno, Pietro A1 - Johansen, Tom H. A1 - Straube, Arthur T1 - Thermally active nanoparticle clusters enslaved by engineered domain wall traps JF - Nature Commun. N2 - The stable assembly of fluctuating nanoparticle clusters on a surface represents a technological challenge of widespread interest for both fundamental and applied research. Here we demonstrate a technique to stably confine in two dimensions clusters of interacting nanoparticles via size-tunable, virtual magnetic traps. We use cylindrical Bloch walls arranged to form a triangular lattice of ferromagnetic domains within an epitaxially grown ferrite garnet film. At each domain, the magnetic stray field generates an effective harmonic potential with a field tunable stifness. The experiments are combined with theory to show that the magnetic confinement is effectively harmonic and pairwise interactions are of dipolar nature, leading to central, strictly repulsive forces. For clusters of magnetic nanoparticles, the stationary collective states arise from the competition between repulsion, confinement and the tendency to fill the central potential well. Using a numerical simulation model as a quantitative map between the experiment and theory we explore the field-induced crystallization process for larger clusters and unveil the existence of three different dynamical regimes. The present method provides a model platform for investigations of the collective phenomena emerging when strongly confined nanoparticle clusters are forced to move in an idealized, harmonic-like potential. Y1 - 2021 U6 - https://doi.org/10.1038/s41467-021-25931-7 VL - 12 SP - 5813 ER - TY - JOUR A1 - Andrle, Anna A1 - Hönicke, Philipp A1 - Gwalt, Grzegorz A1 - Schneider, Philipp-Immanuel A1 - Kayser, Yves A1 - Siewert, Frank A1 - Soltwisch, Victor T1 - Shape- and Element-Sensitive Reconstruction of Periodic Nanostructures with Grazing Incidence X-ray Fluorescence Analysis and Machine Learning JF - Nanomaterials Y1 - 2021 U6 - https://doi.org/10.3390/nano11071647 VL - 11 SP - 1647 ER - TY - JOUR A1 - Nitzke, Isabel A1 - Fackeldey, Konstantin A1 - Vrabec, Jadran T1 - Long range corrections for inhomogeneous fluids containing a droplet or a bubble JF - Molecular Simulation N2 - Long range corrections for molecular simulations of inhomogeneous fluids with a spherical interface are presented. Correction terms for potential energy, force and virial are derived for the monatomic Lennard–Jones fluid. The method is generalised to the Mie potential and arbitrary molecular structures, employing a numerically efficient centre of mass cut-off scheme. The results are validated by a series of droplet simulations for one-centre and two-centre Lennard–Jones fluids with different cut-off radii rc. Systems with rc=8σ provide a check of self-consistence. Further, a system containing a bubble is investigated for the one-centre Lennard–Jones fluid. The equilibrium properties are almost completely independent on the cut-off radius. In comparison with vapour–liquid equilibrium data for systems without a curved interface, all properties show the expected behaviour. Simulation data are used to approximate the surface tension, which is in good agreement with the findings for planar interfaces, thus verifying the present corrections. Y1 - 2021 U6 - https://doi.org/https://doi.org/10.1080/08927022.2021.1954639 SP - 1 EP - 14 ER - TY - CHAP A1 - Perevoznik, Dmitrii A1 - Bose, Surajit A1 - Burger, Sven A1 - Demircan, Ayhan A1 - Morgner, Uwe T1 - Writing Photonic Components in Polymers Using Femtosecond Pulses T2 - Conference on Lasers and Electro-Optics (CLEO) Y1 - 2021 U6 - https://doi.org/10.1364/CLEO_AT.2021.ATh4P.4 SP - ATh4P.4 ER - TY - CHAP A1 - Schnauber, Peter A1 - Schall, Johannes A1 - Bounouar, Samir A1 - Srinivasan, Kartik A1 - Davanco, Marcelo A1 - Song, Jin-Dong A1 - Burger, Sven A1 - Rodt, Sven A1 - Reitzenstein, Stephan T1 - Deterministically fabricated quantum dot – waveguide systems for on-chip quantum optics T2 - Proc. SPIE Y1 - 2021 U6 - https://doi.org/10.1117/12.2596042 VL - 11806 SP - 118060G ER - TY - JOUR A1 - Betz, Fridtjof A1 - Binkowski, Felix A1 - Burger, Sven T1 - RPExpand: Software for Riesz projection expansion of resonance phenomena JF - SoftwareX Y1 - 2021 U6 - https://doi.org/10.1016/j.softx.2021.100763 VL - 15 SP - 100763 ER - TY - JOUR A1 - Wu, Tong A1 - Arrivault, Denis A1 - Duruflé, Marc A1 - Gras, Alexandre A1 - Binkowski, Felix A1 - Burger, Sven A1 - Yan, Wei A1 - Lalanne, Philippe T1 - Efficient hybrid method for the modal analysis of optical microcavities and nanoresonators JF - J. Opt. Soc. Am. A Y1 - 2021 U6 - https://doi.org/10.1364/JOSAA.428224 VL - 38 SP - 1224 ER - TY - GEN A1 - Betz, Fridtjof A1 - Binkowski, Felix A1 - Burger, Sven T1 - RPExpand: Software for Riesz projection expansion of resonance phenomena [Source Code] Y1 - 2021 U6 - https://doi.org/10.24433/CO.2372084.v1 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 - 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 - 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 - Plock, Matthias A1 - Burger, Sven A1 - Schneider, Philipp-Immanuel T1 - Recent advances in Bayesian optimization with applications to parameter reconstruction in optical nano-metrology JF - Proc. SPIE Y1 - 2021 U6 - https://doi.org/10.1117/12.2592266 VL - 11783 SP - 117830J ER - TY - GEN A1 - Raharinirina, N. Alexia A1 - Weber, Marcus A1 - Birk, Ralph A1 - Fackeldey, Konstantin A1 - Klasse, Sarah M. A1 - Richter, Tonio Sebastian T1 - Different Tools and Results for Correspondence Analysis N2 - This is a list of codes generated from ancient egyptian texts. The codes are used for a correspondence analysis (CA). Codes and CA software are available from the linked webpage. Y1 - 2021 U6 - https://doi.org/10.12752/8257 N1 - A detailed description of the software can be found in the code repository at https://github.com/AlexiaNomena/Correspondence_Analysis_User_Friendly (repository version of CA software might include updates). 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 - Tillmann, Peter A1 - Bläsi, Benedikt A1 - Burger, Sven A1 - Hammerschmidt, Martin A1 - Höhn, Oliver A1 - Becker, Christiane A1 - Jäger, Klaus T1 - Optimizing metal grating back reflectors for III-V-on-silicon multijunction solar cells JF - Opt. Express Y1 - 2021 U6 - https://doi.org/10.1364/OE.426761 VL - 29 SP - 22517 ER - TY - JOUR A1 - Liu, Yongtao A1 - Zhou, Zhiguang A1 - Wang, Fan A1 - Kewes, Günter A1 - Wen, Shihui A1 - Burger, Sven A1 - Wakiani, Majid Ebrahimi A1 - Xi, Peng A1 - Yang, Jiong A1 - Yang, Xusan A1 - Benson, Oliver A1 - Jin, Dayong T1 - Axial localization and tracking of self-interference nanoparticles by lateral point spread functions JF - Nat. Commun. Y1 - 2021 U6 - https://doi.org/10.1038/s41467-021-22283-0 VL - 12 SP - 2019 ER - TY - JOUR A1 - Jäger, Klaus A1 - Sutter, Johannes A1 - Hammerschmidt, Martin A1 - Schneider, Philipp-Immanuel A1 - Becker, Christiane T1 - Prospects of light management in perovskite/silicon tandem solar cells JF - Nanophotonics Y1 - 2021 U6 - https://doi.org/10.1515/nanoph-2020-0674 VL - 10 SP - 1991 ER - TY - JOUR A1 - Binkowski, Felix A1 - Wu, Tong A1 - Lalanne, Philippe A1 - Burger, Sven A1 - Govorov, Alexander O. T1 - Hot Electron Generation through Near-Field Excitation of Plasmonic Nanoresonators JF - ACS Photonics Y1 - 2021 U6 - https://doi.org/10.1021/acsphotonics.1c00231 VL - 8 SP - 1243 ER - TY - JOUR A1 - Pakhomov, Anton V. A1 - Hammerschmidt, Martin A1 - Burger, Sven A1 - Pertsch, Thomas A1 - Setzpfand, Frank T1 - Modeling of surface-induced second-harmonic generation from multilayer structures by the transfer matrix method JF - Opt. Express Y1 - 2021 U6 - https://doi.org/10.1364/OE.417066 VL - 29 SP - 9098 ER - TY - JOUR A1 - Manley, Phillip A1 - Segantini, Michele A1 - Ahiboz, Doguscan A1 - Hammerschmidt, Martin A1 - Arnaoutakis, Georgios A1 - MacQueen, Rowan W. A1 - Burger, Sven A1 - Becker, Christiane T1 - Double-layer metasurface for enhanced photon up-conversion JF - APL Photon. Y1 - 2021 U6 - https://doi.org/10.1063/5.0040839 VL - 6 SP - 036103 ER - TY - JOUR A1 - Andrle, Anna A1 - Hönicke, Philipp A1 - Gwalt, Grzegorz A1 - Schneider, Philipp-Immanuel A1 - Kayser, Yves A1 - Siewert, Frank A1 - Soltwisch, Victor T1 - GIXRF and machine learning as metrological tools for shape and element sensitive reconstructions of periodic nanostructures JF - Proc. SPIE Y1 - 2021 U6 - https://doi.org/10.1117/12.2586082 VL - 11611 SP - 116110R ER - TY - JOUR A1 - Winkelmann, Stefanie A1 - Zonker, Johannes A1 - Schütte, Christof A1 - Djurdjevac Conrad, Natasa T1 - Mathematical modeling of spatio-temporal population dynamics and application to epidemic spreading JF - Mathematical Biosciences N2 - Agent based models (ABMs) are a useful tool for modeling spatio-temporal population dynamics, where many details can be included in the model description. Their computational cost though is very high and for stochastic ABMs a lot of individual simulations are required to sample quantities of interest. Especially, large numbers of agents render the sampling infeasible. Model reduction to a metapopulation model leads to a significant gain in computational efficiency, while preserving important dynamical properties. Based on a precise mathematical description of spatio-temporal ABMs, we present two different metapopulation approaches (stochastic and piecewise deterministic) and discuss the approximation steps between the different models within this framework. Especially, we show how the stochastic metapopulation model results from a Galerkin projection of the underlying ABM onto a finite-dimensional ansatz space. Finally, we utilize our modeling framework to provide a conceptual model for the spreading of COVID-19 that can be scaled to real-world scenarios. Y1 - 2021 U6 - https://doi.org/10.1016/j.mbs.2021.108619 VL - 336 PB - Elsevier ER - TY - JOUR A1 - Carderera, Alejandro A1 - Pokutta, Sebastian A1 - Schütte, Christof A1 - Weiser, Martin T1 - An efficient first-order conditional gradient algorithm in data-driven sparse identification of nonlinear dynamics to solve sparse recovery problems under noise JF - Journal of Computational and Applied Mathematics N2 - Governing equations are essential to the study of nonlinear dynamics, often enabling the prediction of previously unseen behaviors as well as the inclusion into control strategies. The discovery of governing equations from data thus has the potential to transform data-rich fields where well-established dynamical models remain unknown. This work contributes to the recent trend in data-driven sparse identification of nonlinear dynamics of finding the best sparse fit to observational data in a large library of potential nonlinear models. We propose an efficient first-order Conditional Gradient algorithm for solving the underlying optimization problem. In comparison to the most prominent alternative algorithms, the new algorithm shows significantly improved performance on several essential issues like sparsity-induction, structure-preservation, noise robustness, and sample efficiency. We demonstrate these advantages on several dynamics from the field of synchronization, particle dynamics, and enzyme chemistry. Y1 - 2021 ER - TY - JOUR A1 - Boltz, Horst-Holger A1 - Kurchan, Jorge A1 - Liu, Andrea J. T1 - Fluctuation distributions of energy minima in complex landscapes JF - Physical Review Research N2 - We discuss the properties of the distributions of energies of minima obtained by gradient descent in complex energy landscapes. We find strikingly similar phenomenology across several prototypical models. We particularly focus on the distribution of energies of minima in the analytically well-understood p-spin-interaction spin-glass model. We numerically find non-Gaussian distributions that resemble the Tracy-Widom distributions often found in problems of random correlated variables, and nontrivial finite-size scaling. Based on this, we propose a picture of gradient-descent dynamics that highlights the importance of a first-passage process in the eigenvalues of the Hessian. This picture provides a concrete link to problems in which the Tracy-Widom distribution is established. Aspects of this first-passage view of gradient-descent dynamics are generic for nonconvex complex landscapes, rationalizing the commonality that we find across models. Y1 - 2020 U6 - https://doi.org/10.1103/PhysRevResearch.3.013061 VL - 3 IS - 1 SP - 013061 ER - TY - JOUR A1 - 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 - 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 - 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 - 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 - Borndörfer, Ralf A1 - Danecker, Fabian A1 - Weiser, Martin T1 - A Discrete-Continuous Algorithm for Free Flight Planning JF - Algorithms N2 - We propose a hybrid discrete-continuous algorithm for flight planning in free flight airspaces. In a first step, our DisCOptER method discrete-continuous optimization for enhanced resolution) computes a globally optimal approximate flight path on a discretization of the problem using the A* method. This route initializes a Newton method that converges rapidly to the smooth optimum in a second step. The correctness, accuracy, and complexity of the method are goverened by the choice of the crossover point that determines the coarseness of the discretization. We analyze the optimal choice of the crossover point and demonstrate the asymtotic superority of DisCOptER over a purely discrete approach. KW - shortest path KW - flight planning KW - free flight KW - discrete-continuous algorithm KW - optimal control KW - discrete optimization Y1 - 2020 U6 - https://doi.org/10.3390/a14010004 SN - 1438-0064 VL - 14 IS - 1 SP - 4 PB - MDPI ER - TY - THES A1 - Santiago, Xavier Garcia T1 - Numerical methods for shape optimization of photonic nanostructures Y1 - 2021 U6 - https://doi.org/10.5445/IR/1000131006 ER - TY - JOUR A1 - Mielke, Alexander A1 - Montefusco, Alberto A1 - Peletier, Mark A. T1 - Exploring families of energy-dissipation landscapes via tilting: three types of EDP convergence JF - Continuum Mechanics and Thermodynamics N2 - This paper revolves around a subtle distinction between two concepts: passing to the limit in a family of gradient systems, on one hand, and deriving effective kinetic relations on the other. The two concepts are strongly related, and in many examples they even appear to be the same. Our main contributions are to show that they are different, to show that well-known techniques developed for the former may give incorrect results for the latter, and to introduce new tools to remedy this. The approach is based on the Energy-Dissipation Principle that provides a variational formulation to gradient-flow equations that allows one to apply techniques from Γ-convergence of functional on states and functionals on trajectories. Y1 - 2021 U6 - https://doi.org/10.1007/s00161-020-00932-x VL - 33 SP - 611 EP - 637 PB - Springer ER - TY - JOUR A1 - Schiela, Anton A1 - Stöcklein, Matthias A1 - Weiser, Martin T1 - A primal dual projection algorithm for efficient constraint preconditioning JF - SIAM Journal on Scientific Computing N2 - We consider a linear iterative solver for large scale linearly constrained quadratic minimization problems that arise, for example, in optimization with PDEs. By a primal-dual projection (PDP) iteration, which can be interpreted and analysed as a gradient method on a quotient space, the given problem can be solved by computing sulutions for a sequence of constrained surrogate problems, projections onto the feasible subspaces, and Lagrange multiplier updates. As a major application we consider a class of optimization problems with PDEs, where PDP can be applied together with a projected cg method using a block triangular constraint preconditioner. Numerical experiments show reliable and competitive performance for an optimal control problem in elasticity. Y1 - 2021 U6 - https://doi.org/10.1137/20M1380739 VL - 43 IS - 6 SP - A4095 EP - A4120 ER - TY - JOUR A1 - Cvetkovic, Nada A1 - Conrad, Tim A1 - Lie, Han Cheng T1 - A Convergent Discretisation Method for Transition Path Theory for Diffusion Processes JF - Multiscale Modeling & Simulation Y1 - 2021 U6 - https://doi.org/10.1137/20M1329354 VL - 19 IS - 1 SP - 242 EP - 266 PB - Society for Industrial and Applied Mathematics ER - TY - JOUR A1 - Jäger, Klaus A1 - Tillmann, Peter A1 - Katz, Eugene A. A1 - Becker, Christiane T1 - Perovskite/silicon tandem solar cells: Effect of luminescent coupling and bifaciality JF - Sol. RRL Y1 - 2021 U6 - https://doi.org/10.1002/solr.202000628 VL - 5 SP - 2000628 ER - TY - JOUR A1 - Chegini, Fatemeh A1 - Kopanicakova, Alena A1 - Krause, Rolf A1 - Weiser, Martin T1 - Efficient Identification of Scars using Heterogeneous Model Hierarchies JF - EP Europace N2 - Aims. Detection and quantification of myocardial scars are helpful both for diagnosis of heart diseases and for building personalized simulation models. Scar tissue is generally charac­terized by a different conduction of electrical excitation. We aim at estimating conductivity-related parameters from endocardial mapping data, in particular the conductivity tensor. Solving this inverse problem requires computationally expensive monodomain simulations on fine discretizations. Therefore, we aim at accelerating the estimation using a multilevel method combining electrophysiology models of different complexity, namely the mono­domain and the eikonal model. Methods. Distributed parameter estimation is performed by minimizing the misfit between simulated and measured electrical activity on the endocardial surface, subject to the mono­domain model and regularization, leading to a constrained optimization problem. We formulate this optimization problem, including the modeling of scar tissue and different regularizations, and design an efficient iterative solver. We consider monodomain grid hierarchies and monodomain-eikonal model hierarchies in a recursive multilevel trust-region method. Results. From several numerical examples, both the efficiency of the method and the estimation quality, depending on the data, are investigated. The multilevel solver is significantly faster than a comparable single level solver. Endocardial mapping data of realistic density appears to be just sufficient to provide quantitatively reasonable estimates of location, size, and shape of scars close to the endocardial surface. Conclusion. In several situations, scar reconstruction based on eikonal and monodomain models differ significantly, suggesting the use of the more accurate but more expensive monodomain model for this purpose. Still, eikonal models can be utilized to accelerate the computations considerably, enabling the use of complex electrophysiology models for estimating myocardial scars from endocardial mapping data. Y1 - 2021 U6 - https://doi.org/10.1093/europace/euaa402 VL - 23 SP - i113 EP - i122 ER - TY - JOUR A1 - Montefusco, Alberto A1 - Peletier, Mark A. A1 - Öttinger, Hans Christian T1 - A Framework of Nonequilibrium Statistical Mechanics. II. Coarse-Graining JF - Journal of Non-Equilibrium Thermodynamics N2 - For a given thermodynamic system, and a given choice of coarse-grained state variables, the knowledge of a force-flux constitutive law is the basis for any nonequilibrium modeling. In the first paper of this series we established how, by a generalization of the classical fluctuation-dissipation theorem (FDT), the structure of a constitutive law is directly related to the distribution of the fluctuations of the state variables. When these fluctuations can be expressed in terms of diffusion processes, one may use Green–Kubo-type coarse-graining schemes to find the constitutive laws. In this paper we propose a coarse-graining method that is valid when the fluctuations are described by means of general Markov processes, which include diffusions as a special case. We prove the success of the method by numerically computing the constitutive law for a simple chemical reaction A⇄B. Furthermore, we show that, for such a system, one cannot find a consistent constitutive law by any Green–Kubo-like scheme. Y1 - 2021 U6 - https://doi.org/10.1515/jnet-2020-0069 VL - 46 IS - 1 SP - 15 EP - 33 PB - De Gruyter ER - TY - JOUR A1 - Öttinger, Hans Christian A1 - Montefusco, Alberto A1 - Peletier, Mark A. T1 - A Framework of Nonequilibrium Statistical Mechanics. I. Role and Types of Fluctuations JF - Journal of Non-Equilibrium Thermodynamics N2 - Understanding the fluctuations by which phenomenological evolution equations with thermodynamic structure can be enhanced is the key to a general framework of nonequilibrium statistical mechanics. These fluctuations provide an idealized representation of microscopic details. We consider fluctuation-enhanced equations associated with Markov processes and elaborate the general recipes for evaluating dynamic material properties, which characterize force-flux constitutive laws, by statistical mechanics. Markov processes with continuous trajectories are conveniently characterized by stochastic differential equations and lead to Green–Kubo-type formulas for dynamic material properties. Markov processes with discontinuous jumps include transitions over energy barriers with the rates calculated by Kramers. We describe a unified approach to Markovian fluctuations and demonstrate how the appropriate type of fluctuations (continuous versus discontinuous) is reflected in the mathematical structure of the phenomenological equations. Y1 - 2021 U6 - https://doi.org/10.1515/jnet-2020-0068 VL - 46 IS - 1 SP - 1 EP - 13 PB - De Gruyter ER - TY - JOUR A1 - Sechi, Renata A1 - Sikorski, Alexander A1 - Weber, Marcus T1 - Estimation of the Koopman Generator by Newton's Extrapolation JF - Multiscale Modeling and Simulation N2 - This article addresses the problem of estimating the Koopman generator of a Markov process. The direct computation of the infinitesimal generator is not easy because of the discretization of the state space, in particular because of the trade-off inherent in the choice of the best lag time to study the process. Short lag times implies a strong discretization of the state space and a consequent loss of Markovianity. Large lag times bypass events on fast timescales. We propose a method to approximate the generator with the computation of the Newton polynomial extrapolation. This technique is a multistep approach which uses as its input Koopman transfer operators evaluated for a series of lag times. Thus, the estimated infinitesimal generator combines information from different time resolutions and does not bias only fast- or slow-decaying dynamics. We show that the multi-scale Newton method can improve the estimation of the generator in comparison to the computation using finite difference or matrix logarithm methods. Y1 - 2020 U6 - https://doi.org/10.1137/20M1333006 VL - 19 IS - 2 SP - 758 EP - 774 PB - SIAM ER - TY - JOUR A1 - Sikorski, Alexander A1 - Weber, Marcus A1 - Schütte, Christof T1 - The Augmented Jump Chain JF - Advanced Theory and Simulations N2 - Modern methods of simulating molecular systems are based on the mathematical theory of Markov operators with a focus on autonomous equilibrated systems. However, non-autonomous physical systems or non-autonomous simulation processes are becoming more and more important. A representation of non-autonomous Markov jump processes is presented as autonomous Markov chains on space-time. Augmenting the spatial information of the embedded Markov chain by the temporal information of the associated jump times, the so-called augmented jump chain is derived. The augmented jump chain inherits the sparseness of the infinitesimal generator of the original process and therefore provides a useful tool for studying time-dependent dynamics even in high dimensions. Furthermore, possible generalizations and applications to the computation of committor functions and coherent sets in the non-autonomous setting are discussed. After deriving the theoretical foundations, the concepts with a proof-of-concept Galerkin discretization of the transfer operator of the augmented jump chain applied to simple examples are illustrated. Y1 - 2021 U6 - https://doi.org/10.1002/adts.202000274 VL - 4 IS - 4 PB - Wiley-VCH ER - TY - JOUR A1 - 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 - Sharma, Upanshu A1 - Zhang, Wei T1 - Non-reversible sampling schemes on submanifolds JF - SIAM Journal on Numerical Analysis N2 - Calculating averages with respect to probability measures on submanifolds is often necessary in various application areas such as molecular dynamics, computational statistical mechanics and Bayesian statistics. In recent years, various numerical schemes have been proposed in the literature to study this problem based on appropriate reversible constrained stochastic dynamics. In this paper we present and analyse a non-reversible generalisation of the projection-based scheme developed by one of the authors [ESAIM: M2AN, 54 (2020), pp. 391-430]. This scheme consists of two steps - starting from a state on the submanifold, we first update the state using a non-reversible stochastic differential equation which takes the state away from the submanifold, and in the second step we project the state back onto the manifold using the long-time limit of a ordinary differential equation. We prove the consistency of this numerical scheme and provide quantitative error estimates for estimators based on finite-time running averages. Furthermore, we present theoretical analysis which shows that this scheme outperforms its reversible counterpart in terms of asymptotic variance. We demonstrate our findings on an illustrative test example. Y1 - 2020 U6 - https://doi.org/10.1137/20M1378752 VL - 59 IS - 6 SP - 2989 EP - 3031 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - 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 - Zhang, Wei T1 - Some new results on relative entropy production, time reversal, and optimal control of time-inhomogeneous diffusion processes JF - Journal of Mathematical Physics N2 - This paper studies time-inhomogeneous nonequilibrium diffusion processes, including both Brownian dynamics and Langevin dynamics. We derive upper bounds of the relative entropy production of the time-inhomogeneous process with respect to the transient invariant probability measures. We also study the time reversal of the reverse process in Crooks' fluctuation theorem. We show that the time reversal of the reverse process coincides with the optimally controlled forward process that leads to zero variance importance sampling estimator based on Jarzynski's equality. Y1 - 2021 U6 - https://doi.org/10.1063/5.0038740 VL - 62 IS - 4 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 - 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 - 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 - 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 - Lie, Han Cheng A1 - Sullivan, T. J. A1 - Teckentrup, Aretha T1 - Error bounds for some approximate posterior measures in Bayesian inference JF - Numerical Mathematics and Advanced Applications ENUMATH 2019 Y1 - 2020 U6 - https://doi.org/10.1007/978-3-030-55874-1_26 SP - 275 EP - 283 PB - Springer ER - TY - JOUR A1 - Bittracher, Andreas A1 - Klus, Stefan A1 - Hamzi, Boumediene A1 - Schütte, Christof T1 - Dimensionality Reduction of Complex Metastable Systems via Kernel Embeddings of Transition Manifolds JF - Journal of Nonlinear Science N2 - We present a novel kernel-based machine learning algorithm for identifying the low-dimensional geometry of the effective dynamics of high-dimensional multiscale stochastic systems. Recently, the authors developed a mathematical framework for the computation of optimal reaction coordinates of such systems that is based on learning a parameterization of a low-dimensional transition manifold in a certain function space. In this article, we enhance this approach by embedding and learning this transition manifold in a reproducing kernel Hilbert space, exploiting the favorable properties of kernel embeddings. Under mild assumptions on the kernel, the manifold structure is shown to be preserved under the embedding, and distortion bounds can be derived. This leads to a more robust and more efficient algorithm compared to the previous parameterization approaches. Y1 - 2021 U6 - https://doi.org/10.1007/s00332-020-09668-z VL - 31 ER - TY - JOUR A1 - Bittracher, Andreas A1 - Schütte, Christof T1 - A probabilistic algorithm for aggregating vastly undersampled large Markov chains JF - Physica D: Nonlinear Phenomena Y1 - 2021 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-75874 VL - 416 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 - Schäfer, Florian A1 - Sullivan, T. J. A1 - Owhadi, Houman T1 - Compression, inversion, and approximate PCA of dense kernel matrices at near-linear computational complexity JF - Multiscale Modeling and Simulation Y1 - 2021 U6 - https://doi.org/10.1137/19M129526X VL - 19 IS - 2 SP - 688 EP - 730 ER - TY - JOUR A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Objective priors in the empirical Bayes framework JF - Scandinavian Journal of Statistics N2 - When dealing with Bayesian inference the choice of the prior often remains a debatable question. Empirical Bayes methods offer a data-driven solution to this problem by estimating the prior itself from an ensemble of data. In the nonparametric case, the maximum likelihood estimate is known to overfit the data, an issue that is commonly tackled by regularization. However, the majority of regularizations are ad hoc choices which lack invariance under reparametrization of the model and result in inconsistent estimates for equivalent models. We introduce a nonparametric, transformation-invariant estimator for the prior distribution. Being defined in terms of the missing information similar to the reference prior, it can be seen as an extension of the latter to the data-driven setting. This implies a natural interpretation as a trade-off between choosing the least informative prior and incorporating the information provided by the data, a symbiosis between the objective and empirical Bayes methodologies. Y1 - 2021 U6 - https://doi.org/10.1111/sjos.12485 VL - 48 IS - 4 SP - 1212 EP - 1233 PB - Wiley Online Library ER -