@misc{SikorskiSechiHelfmann2021, author = {Sikorski, Alexander and Sechi, Renata and Helfmann, Luzie}, title = {cmdtools}, doi = {https://doi.org/10.5281/zenodo.4749331}, year = {2021}, abstract = {Python implementation of severals tools (PCCA, AJC, SQRA, P/Q estimation) for the analysis of dynamical systems from the transfer operator perspective.}, language = {en} } @article{DonatiWeberKeller2021, author = {Donati, Luca and Weber, Marcus and Keller, Bettina G.}, title = {Markov models from the square root approximation of the Fokker-Planck equation: calculating the grid-dependent flux}, volume = {33}, journal = {Journal of Physics: Condensed Matter}, number = {11}, doi = {10.1088/1361-648X/abd5f7}, pages = {115902}, year = {2021}, abstract = {Molecular dynamics (MD) are extremely complex, yet understanding the slow components of their dynamics is essential to understanding their macroscopic properties. To achieve this, one models the MD as a stochastic process and analyses the dominant eigenfunctions of the associated Fokker-Planck operator, or of closely related transfer operators. So far, the calculation of the discretized operators requires extensive MD simulations. The square-root approximation of the Fokker-Planck equation is a method to calculate transition rates as a ratio of the Boltzmann densities of neighboring grid cells times a flux, and can in principle be calculated without a simulation. In a previous work we still used MD simulations to determine the flux. Here, we propose several methods to calculate the exact or approximate flux for various grid types, and thus estimate the rate matrix without a simulation. Using model potentials we test computational efficiency of the methods, and the accuracy with which they reproduce the dominant eigenfunctions and eigenvalues. For these model potentials, rate matrices with up to \$\mathcal{O}\left(1{0}^{6}\right)\$ states can be obtained within seconds on a single high-performance compute server if regular grids are used.}, language = {en} } @inproceedings{BinkowskiBetzColometal.2021, author = {Binkowski, Felix and Betz, Fridtjof and Colom, Remi and Hammerschmidt, Martin and Schneider, Philipp-Immanuel and Zschiedrich, Lin and Burger, Sven}, title = {Contour integral methods for resonance phenomena in nano-optics}, booktitle = {The 11th International Conference on Metamaterials, Photonic Crystals and Plasmonics (META 2021)}, editor = {Zouhdi, S. and Pawlak, D. and Kudelski, A.}, issn = {2429-1390}, pages = {523}, year = {2021}, language = {en} } @article{DjurdjevacGraeserHerbert2021, author = {Djurdjevac, Ana and Gr{\"a}ser, Carsten and Herbert, Philip J.}, title = {Non-stationary incompressible linear fluid equations in a moving domain}, journal = {arXiv}, arxiv = {http://arxiv.org/abs/2106.08312}, year = {2021}, abstract = {This article considers non-stationary incompressible linear fluid equations in a moving domain. We demonstrate the existence and uniqueness of an appropriate weak formulation of the problem by making use of the theory of time-dependent Bochner spaces. It is not possible to directly apply established evolving Hilbert space theory due to the incompressibility constraint. After we have established the well-posedness, we derive and analyse a time discretisation of the system.}, language = {en} } @article{Djurdjevac2021, author = {Djurdjevac, Ana}, title = {Linear Parabolic Problems in Random Moving Domains}, volume = {9}, number = {2}, doi = {https://doi.org/10.1137/19M1284889}, pages = {848 -- 879}, year = {2021}, abstract = {We consider linear parabolic equations on a random non-cylindrical domain. Utilizing the domain mapping method, we write the problem as a partial differential equation with random coefficients on a cylindrical deterministic domain. Exploiting the deterministic results concerning equations on non-cylindrical domains, we state the necessary assumptions about the velocity filed and in addition, about the flow transformation that this field generates. In this paper we consider both cases, the uniformly bounded with respect to the sample and log-normal type transformation. In addition, we give an explicit example of a log-normal type transformation and prove that it does not satisfy the uniformly bounded condition. We define a general framework for considering linear parabolic problems on random non-cylindrical domains. As the first example, we consider the heat equation on a random tube domain and prove its well-posedness. Moreover, as the other example we consider the parabolic Stokes equation which illustrates the case when it is not enough just to study the plain-back transformation of the function, but instead to consider for example the Piola type transformation, in order to keep the divergence free property.}, language = {en} } @article{BergPloentzkeSiebertetal.2021, author = {Berg, Mascha and Pl{\"o}ntzke, Julia and Siebert, Heike and R{\"o}blitz, Susanna}, title = {Modelling Oscillatory Patterns in the Bovine Estrous Cycle with Boolean Delay Equations}, volume = {83}, journal = {Bulletin of Mathematical Biology}, number = {121}, doi = {10.1007/s11538-021-00942-z}, pages = {1 -- 25}, year = {2021}, abstract = {Boolean delay equations (BDEs), with their relatively simple and intuitive mode of modelling, have been used in many research areas including, for example, climate dynamics and earthquake propagation. Their application to biological systems has been scarce and limited to the molecular level. Here, we derive and present two BDE models. One is directly derived from a previously published ordinary differential equation (ODE) model for the bovine estrous cycle, whereas the second model includes a modification of a particular biological mechanism. We not only compare the simulation results from the BDE models with the trajectories of the ODE model, but also validate the BDE models with two additional numerical experiments. One experiment induces a switch in the oscillatory pattern upon changes in the model parameters, and the other simulates the administration of a hormone that is known to shift the estrous cycle in time. The models presented here are the first BDE models for hormonal oscillators, and the first BDE models for drug administration. Even though automatic parameter estimation still remains challenging, our results support the role of BDEs as a framework for the systematic modelling of complex biological oscillators.}, language = {en} } @misc{Mathew2021, type = {Master Thesis}, author = {Mathew, Maria}, title = {Numerical treatment of radiative heat transfer on cooling process}, pages = {53}, year = {2021}, abstract = {This thesis examines how taking into account surface to surface radiation impacts the cooling process in general. We formulate the non local bound- ary condition after introducing the general setting for the cooling model. In section 3, the mathematical description of the radiative heat transfer is dis- cussed. We cover the implementation of the radiative matrix in section 4, which is followed by a brief explanation of the radiative matrix's structure and several techniques to dealing with the accompanying challenges. We investigate the importance of radiative heat transport by applying the given approach to a two-dimensional geometry and computing the ensuing cooling curves. We compare the findings of our computation to those ac- quired from experiment conducted and find that they are extremely similar. There is a considerable difference (of about 35\%) in the time of cooling of the surface where there is a possibility of influence of radiation from the second surface to that of the surface with no influence at all. Although it is possible to infer that heat convection plays a role in the total result, this has yet to be proved. However, one can clearly see the significance of the surface to surface radiative heat transfer on these parts confirming the research question posed at the begining. The effect of the surface to surface radiative heat transfer has an influence on the resulting cooling time and should be considered in the model.}, language = {en} } @article{BirkRaharinirinaFackeldeyetal.2021, author = {Birk, Ralph and Raharinirina, N. Alexia and Fackeldey, Konstantin and Richter, Tonio Sebastian and Weber, Marcus}, title = {Inferring cultural and social processes based on patterns of statistical relationships between Synodal texts}, year = {2021}, abstract = {In this paper, we explore the relationship patterns between Ancient Egyptian texts of the corpus ``Synodal decrees'', which are originating between 243 and 185 BCE, during the Ptolemaic period. Particularly, we are interested in analyzing the grammatical features of the different texts. Conventional data analysis methods such as correspondence Analysis are very useful to explore the patterns of statistical interdependence between categories of variables. However, it is based on a PCA-like dimension-reduction method and turned out to be unsuitable for our dataset due to the high dimensionality of our data representations. Additionally, the similarity between pairs of texts and pairs of grammatical features is observed through the distance between their representation, but the degree of association between a particular grammatical feature and a text is not. Here, we applied a qualitative Euclidean embedding method that provides a new Euclidean representation of the categories of variables. This new representation of the categories is constructed in such a way that all the patterns of statistical interdependence, similarity, and association, are seen through the Euclidean distance between them. Nevertheless, the PCA-like dimension-reduction method also performed poorly on our new representation. Therefore, we obtained a two-dimensional visualization using non-linear methods such UMAP or t-SNE. Although these dimension-reduction methods reduced the interpretability of interpoint distances, we were still able to identify important similarity patterns between the Synodal text as well as their association patterns with the grammatical features.}, language = {en} } @article{HungenbergSchwedeWulkowetal.2021, author = {Hungenberg, Klaus-Dieter and Schwede, Christian and Wulkow, Michael and Wulkow, Niklas}, title = {Determination of reactivity ratios for acrylic acid and its dimer from classical parameter estimation and Bayesian approach}, journal = {Canadian Journal of Chemical Engineering}, doi = {https://doi.org/10.1002/cjce.24330}, pages = {1 -- 10}, year = {2021}, abstract = {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.}, language = {en} } @article{DriesenElKhatibWulkowetal.2021, author = {Driesen, Joshua and El-Khatib, Ziad and Wulkow, Niklas and Joblin, Mitchell and Vasileva, Iskriyana and Gl{\"u}cker, Andreas and Kruspel, Valentin and Vogel, Catherine}, title = {Data-Powered Positive Deviance during the SARS-CoV-2 Pandemic—An Ecological Pilot Study of German Districts}, volume = {18}, journal = {International Journal of Environmental Research and Public Health}, number = {9765}, doi = {https://doi.org/10.3390/ijerph18189765}, pages = {1 -- 29}, year = {2021}, abstract = {We introduced the mixed-methods Data-Powered Positive Deviance (DPPD) framework as a potential addition to the set of tools used to search for effective response strategies against the SARS-CoV-2 pandemic. For this purpose, we conducted a DPPD study in the context of the early stages of the German SARS-CoV-2 pandemic. We used a framework of scalable quantitative methods to identify positively deviant German districts that is novel in the scientific literature on DPPD, and subsequently employed qualitative methods to identify factors that might have contributed to their comparatively successful reduction of the forward transmission rate. Our qualitative analysis suggests that quick, proactive, decisive, and flexible/pragmatic actions, the willingness to take risks and deviate from standard procedures, good information flows both in terms of data collection and public communication, alongside the utilization of social network effects were deemed highly important by the interviewed districts. Our study design with its small qualitative sample constitutes an exploratory and illustrative effort and hence does not allow for a clear causal link to be established. Thus, the results cannot necessarily be extrapolated to other districts as is. However, the findings indicate areas for further research to assess these strategies' effectiveness in a broader study setting. We conclude by stressing DPPD's strengths regarding replicability, scalability, adaptability, as well asits focus on local solutions, which make it a promising framework to be applied in various contexts,e.g., in the context of the Global South.}, language = {en} } @article{StoutColomBonodetal.2021, author = {Stout, Brian and Colom, Remi and Bonod, Nicolas and McPhedran, Ross C}, title = {Spectral expansions of open and dispersive optical systems: Gaussian regularization and convergence}, volume = {23}, journal = {New J. Phys.}, doi = {10.1088/1367-2630/ac10a6}, pages = {083004}, year = {2021}, language = {en} } @inproceedings{KostreSunkaraDjurdjevacConrad2021, author = {Kostre, Margarita and Sunkara, Vikram and Djurdjevac Conrad, Natasa}, title = {Inference of historical influence networks}, booktitle = {International Conference on Complex Networks \& Their Applications Proceedings}, edition = {International Conference on Complex Networks \& Their Applications}, doi = {10.12752/8558}, pages = {110 -- 113}, year = {2021}, abstract = {We study the romanization process of northern Africa from 50 BC till 300 AD. Our goal is to infer the communication strength between different subregions, based on the evolution of the status of cities. Herefore, we use the general inverse infection model, that infers the weights of a known underlying network, given observations of the spreading on this network. As infection process we choose the SI metapopulation model, where I stands for a city with a Roman status. To solve the minimization problem we use the particle swarm optimization algorithm with a specific choice of parameters.}, language = {de} } @masterthesis{Chemnitz2021, type = {Bachelor Thesis}, author = {Chemnitz, Robin}, title = {Spectral Coarse-Graining of SIS Dynamics on Networks}, year = {2021}, language = {en} } @article{Schweigart2021, author = {Schweigart, Fleur}, title = {KARTHAGO, TUNESIEN Datengesteuerte Modellierung des Romanisierungsprozesses in Nordafrika}, volume = {Fasc. 1}, journal = {e-Forschungsberichte des Deutschen Arch{\"a}ologischen Instituts. Berlin: Deutsches Arch{\"a}ologisches Institut, Zentrale.}, doi = {https://doi.org/10.34780/d469-84h2}, year = {2021}, language = {de} } @inproceedings{RickertGaoUrbanetal.2021, author = {Rickert, Lucas and Gao, Timm and Urban, Felix and Große, Jan and Srocka, Nicole and Rodt, Sven and Musial, Anna and Zolnacz, Kinga and Mergo, Pawel and Dybka, Kamil and Urbanczyk, Waclaw and Sek, Grzegorz and Burger, Sven and Reitzenstein, Stephan and Tobias, Heindel}, title = {A Quantum Key Distribution Testbed Using Plug\&Play Telecom-Wavelength Single-Photons}, booktitle = {Quantum Information and Measurement VI}, doi = {10.1364/QIM.2021.M2C.1}, pages = {M2C.1}, year = {2021}, language = {en} } @misc{RiberaBorrellQuerRichteretal.2021, author = {Ribera Borrell, Enric and Quer, Jannes and Richter, Lorenz and Sch{\"u}tte, Christof}, title = {Improving control based importance sampling strategies for metastable diffusions via adapted metadynamics}, issn = {1438-0064}, year = {2021}, abstract = {Sampling rare events in metastable dynamical systems is often a computationally expensive task and one needs to resort to enhanced sampling methods such as importance sampling. Since we can formulate the problem of finding optimal importance sampling controls as a stochastic optimization problem, this then brings additional numerical challenges and the convergence of corresponding algorithms might as well suffer from metastabilty. In this article we address this issue by combining systematic control approaches with the heuristic adaptive metadynamics method. Crucially, we approximate the importance sampling control by a neural network, which makes the algorithm in principle feasible for high dimensional applications. We can numerically demonstrate in relevant metastable problems that our algorithm is more effective than previous attempts and that only the combination of the two approaches leads to a satisfying convergence and therefore to an efficient sampling in certain metastable settings.}, language = {en} } @article{BauerWeberDiehlWieseneckeretal.2021, author = {Bauer, Wolfgang and Weber, Marcus and Diehl-Wiesenecker, Eva and Galtung, Noa and Prpic, Monika and Somasundaram, Rajan and Tauber, Rudolf and Schwenk, Jochen and Micke, Patrick and Kappert, Kai}, title = {Plasma Proteome Fingerprints Reveal Distinctiveness and Clinical Outcome of SARS-CoV-2 Infection}, volume = {13}, journal = {Viruses}, number = {12}, doi = {10.3390/v13122456}, pages = {2456}, year = {2021}, abstract = {We evaluated how plasma proteomic signatures in patients with suspected COVID-19 can unravel the pathophysiology, and determine kinetics and clinical outcome of the infection. We identified distinct plasma proteins linked to the presence and course of COVID-19. These plasma proteomic findings may translate to a protein fingerprint, helping to assist clinical management decisions.}, language = {en} } @article{FackeldeyGorgullaWeber2021, author = {Fackeldey, Konstantin and Gorgulla, Christoph and Weber, Marcus}, title = {Neue Medikamente dank Supercomputern}, journal = {Spektrum der Wissenschaft}, number = {11}, pages = {40 -- 46}, year = {2021}, abstract = {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.}, language = {de} } @inproceedings{SchneiderGarciaSantiagoWiegandetal.2021, author = {Schneider, Philipp-Immanuel and Garcia-Santiago, Xavier and Wiegand, Benjamin and Anton, Oliver and Krutzik, Markus and Rockstuhl, Carsten and Burger, Sven}, title = {Advances in Bayesian optimization for photonics and quantum atom optics applications}, booktitle = {OSA Advanced Photonics Congress}, doi = {10.1364/IPRSN.2021.JTh1E.2}, pages = {JTh1E.2}, year = {2021}, language = {en} } @article{AhibozAndresenManleyetal.2021, author = {Ahiboz, Doguscan and Andresen, Elina and Manley, Phillip and Resch-Genger, Ute and W{\"u}rth, Christian and Becker, Christiane}, title = {Metasurface-Enhanced Photon Upconversion upon 1550 nm Excitation}, volume = {9}, journal = {Adv. Opt. Mater.}, doi = {10.1002/adom.202101285}, pages = {2101285}, year = {2021}, language = {en} } @inproceedings{AhibozAndresenManleyetal.2021, author = {Ahiboz, Doguscan and Andresen, Elina and Manley, Phillip and Resch Genger, Ute and W{\"u}rth, Christian and Becker, Christiane}, title = {Enhanced Photon Upconversion Using Erbium-Doped Nanoparticles Interacting with Silicon Metasurfaces}, booktitle = {IEEE 48th Photovoltaic Specialists Conference (PVSC)}, doi = {10.1109/PVSC43889.2021.9518495}, pages = {2008 -- 2010}, year = {2021}, language = {en} } @inproceedings{PerevoznikBoseBurgeretal.2021, author = {Perevoznik, Dmitrii and Bose, Surajit and Burger, Sven and Demircan, Ayhan and Morgner, Uwe}, title = {Photonic components in polymers made by femtosecond pulses}, booktitle = {Conference on Lasers and Electro-Optics Europe \& European Quantum Electronics Conference (CLEO/Europe-EQEC)}, doi = {10.1109/CLEO/Europe-EQEC52157.2021.9541906}, pages = {1}, year = {2021}, language = {en} } @article{StraubeWinkelmannSchuetteetal.2021, author = {Straube, Arthur and Winkelmann, Stefanie and Sch{\"u}tte, Christof and H{\"o}fling, Felix}, title = {Stochastic pH oscillations in a model of the urea-urease reaction confined to lipid vesicles}, volume = {12}, journal = {J. Phys. Chem. Lett.}, doi = {10.1021/acs.jpclett.1c03016}, pages = {9888 -- 9893}, year = {2021}, abstract = {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.}, language = {en} } @article{SechiWeberHeyne2021, author = {Sechi, Renata and Weber, Marcus and Heyne, Karsten}, title = {MSM building and projection for the analysis of time-resolved spectra}, journal = {Proceedings in Applied Mathematics and Mechanics}, number = {21}, doi = {https://doi.org/10.1002/pamm.202100102}, url = {http://nbn-resolving.de/https://onlinelibrary.wiley.com/doi/10.1002/pamm.202100102}, pages = {e202100102}, year = {2021}, abstract = {Understanding the kinetics between the components of time-resolved spectra is a crucial step in the study of photo-activatedprocesses. However, modeling the kinetics requires usually some a priori knowledge about the system. In our approach, webuild a Markov State Model (MSM) from the spectral data, and obtain a Koopman transition matrix K(t). With genPCCA,an invariant subspace projection, we project the process into its metastable components. The result of the application of gen-PCCA is a transition matrix Kc(t), from which we can read the transition probability between the metastable components of the reaction. We discuss the application of this analysis method to the transient absorption spectrum of brominated Al-corrole}, language = {en} } @phdthesis{Moldenhauer2021, author = {Moldenhauer, Marian}, title = {Adaptive Algorithms in Optimization under PDE Constraints}, pages = {106}, year = {2021}, abstract = {In this thesis, adaptive algorithms in optimization under PDE constraints have been inves- tigated. In its application, the aim of optimization is to increase the longevity of implants, namely the hip joint implant, and in doing so to minimize stress shielding and simultaneously minimize the influence of locally high stresses, that, above a threshold value, are malign to the bone structure. Under the constraint of the equilibrium of forces, describing an elastodynamic setup, coupled with a contact inequality condition, a computationally expensive problem formulation is given. The first step to make the solution of the given problem possible and efficient was to change over to the spatial equilibrium equation, thus rendering an elastostatic setup. Subsequently the intrinsically dynamic motions - trajectories in the load domain - were converted to the static setup. Thus, the trajectories are marginalized to the load domain and characterized with probability distributions. Therefore the solving of the PDE constraint, the contact problem, is simplified. Yet in the whole optimization process, the solving of the PDE, the spatial equilibrium equation together with the contact condition has the most expensive contribution still and hence needed further reduction. This was achieved by application of Kriging interpolation to the load responses of the integrated distribution of stress difference and the maximum stresses. The interpolation of the two response surfaces only needs comparatively few PDE solves to set up the models. Moreover, the Kriging models can be adaptively extended by sequentially adding sample-response pairs. For this the Kriging inherent variance is used to estimate ideal new sample locations with maximum variance values. In doing so, the overall interpolation variance and therefore the interpolation error is reduced. For the integration of the integrated stress differences and penalty values on the relative high dimensional load domain Monte Carlo integration was implemented, averting the curse of dimension. Here, the motion's probability distribution combined with patient specific data of motion frequencies is taken advantage of, making obsolete the use of the otherwise necessary importance sampling. Throughout the optimization, the FE-discretization error and the subsequently attached errors entering the solution process via PDE discretization and approximative solving of the PDE, Kriging interpolation and Monte Carlo integration need to decrease. While the FE-discretization error and the solution of the elastostatic contact problem were assumed precise enough, numerics showed, that the interpolation and integration errors can be controlled by adaptive refinement of the respective methods. For this purpose comparable error quantities for the particular algorithms were introduced and effectively put to use. For the implant position's optimization, the derivative of the objective function was derived using the implicit function theorem. As the FE-discretization changes with implant position modifications big enough, a special line search had to be used to deal with the discontinuities in the objective function. The interplay and performance of the subalgorithms was demonstrated numerically on a reduced 2D setup of a hip joint with and without the implant. Consequently the load domain and the control variable were also limited to the 2D case.}, language = {en} } @inproceedings{TillmannJaegerKatzetal.2021, author = {Tillmann, Peter and J{\"a}ger, Klaus and Katz, Eugene and Becker, Christiane}, title = {Relaxed current-matching constraints by bifacial operation and luminescent coupling in perovskite/silicon tandem solar cells}, booktitle = {2021 IEEE 48th Photovoltaic Specialists Conference (PVSC)}, doi = {10.1109/PVSC43889.2021.9518891}, pages = {1086}, year = {2021}, language = {en} } @misc{NiemannSchuetteKlus2021, author = {Niemann, Jan-Hendrik and Sch{\"u}tte, Christof and Klus, Stefan}, title = {Simulation data: Data-driven model reduction of agent-based systems using the Koopman generator}, volume = {16}, journal = {PLOS ONE}, number = {5}, doi = {http://doi.org/10.5281/zenodo.4522119}, year = {2021}, language = {en} } @article{TiernoJohansenStraube2021, author = {Tierno, Pietro and Johansen, Tom H. and Straube, Arthur}, title = {Thermally active nanoparticle clusters enslaved by engineered domain wall traps}, volume = {12}, journal = {Nature Commun.}, doi = {10.1038/s41467-021-25931-7}, pages = {5813}, year = {2021}, abstract = {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.}, language = {en} } @article{AndrleHoenickeGwaltetal.2021, author = {Andrle, Anna and H{\"o}nicke, Philipp and Gwalt, Grzegorz and Schneider, Philipp-Immanuel and Kayser, Yves and Siewert, Frank and Soltwisch, Victor}, title = {Shape- and Element-Sensitive Reconstruction of Periodic Nanostructures with Grazing Incidence X-ray Fluorescence Analysis and Machine Learning}, volume = {11}, journal = {Nanomaterials}, doi = {10.3390/nano11071647}, pages = {1647}, year = {2021}, language = {en} } @article{NitzkeFackeldeyVrabec2021, author = {Nitzke, Isabel and Fackeldey, Konstantin and Vrabec, Jadran}, title = {Long range corrections for inhomogeneous fluids containing a droplet or a bubble}, journal = {Molecular Simulation}, doi = {https://doi.org/10.1080/08927022.2021.1954639}, pages = {1 -- 14}, year = {2021}, abstract = {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.}, language = {en} } @inproceedings{PerevoznikBoseBurgeretal.2021, author = {Perevoznik, Dmitrii and Bose, Surajit and Burger, Sven and Demircan, Ayhan and Morgner, Uwe}, title = {Writing Photonic Components in Polymers Using Femtosecond Pulses}, booktitle = {Conference on Lasers and Electro-Optics (CLEO)}, doi = {10.1364/CLEO_AT.2021.ATh4P.4}, pages = {ATh4P.4}, year = {2021}, language = {en} } @inproceedings{SchnauberSchallBounouaretal.2021, author = {Schnauber, Peter and Schall, Johannes and Bounouar, Samir and Srinivasan, Kartik and Davanco, Marcelo and Song, Jin-Dong and Burger, Sven and Rodt, Sven and Reitzenstein, Stephan}, title = {Deterministically fabricated quantum dot - waveguide systems for on-chip quantum optics}, volume = {11806}, booktitle = {Proc. SPIE}, doi = {10.1117/12.2596042}, pages = {118060G}, year = {2021}, language = {en} } @article{WangCockayneChkrebtiietal.2021, author = {Wang, Junyang and Cockayne, Jon and Chkrebtii, Oksana and Sullivan, T. J. and Oates, Chris}, title = {Bayesian numerical methods for nonlinear partial differential equations}, volume = {31}, journal = {Statistics and Computing}, number = {5}, doi = {10.1007/s11222-021-10030-w}, year = {2021}, language = {en} } @article{BetzBinkowskiBurger2021, author = {Betz, Fridtjof and Binkowski, Felix and Burger, Sven}, title = {RPExpand: Software for Riesz projection expansion of resonance phenomena}, volume = {15}, journal = {SoftwareX}, doi = {10.1016/j.softx.2021.100763}, pages = {100763}, year = {2021}, language = {en} } @article{WuArrivaultDurufleetal.2021, author = {Wu, Tong and Arrivault, Denis and Durufl{\´e}, Marc and Gras, Alexandre and Binkowski, Felix and Burger, Sven and Yan, Wei and Lalanne, Philippe}, title = {Efficient hybrid method for the modal analysis of optical microcavities and nanoresonators}, volume = {38}, journal = {J. Opt. Soc. Am. A}, doi = {10.1364/JOSAA.428224}, pages = {1224}, year = {2021}, language = {en} } @misc{BetzBinkowskiBurger2021, author = {Betz, Fridtjof and Binkowski, Felix and Burger, Sven}, title = {RPExpand: Software for Riesz projection expansion of resonance phenomena [Source Code]}, doi = {10.24433/CO.2372084.v1}, year = {2021}, language = {de} } @phdthesis{Pakhomov2021, author = {Pakhomov, Anton}, title = {Efficient modeling and optimization of surface second-harmonic generation from nanophotonic components}, doi = {10.22032/dbt.49369}, year = {2021}, language = {en} } @misc{RayThiesSunkaraetal.2021, author = {Ray, Sourav and Thies, Arne and Sunkara, Vikram and Wulkow, Hanna and Celik, {\"O}zg{\"u}r and Yerg{\"o}z, Fatih and Sch{\"u}tte, Christof and Stein, Christoph and Weber, Marcus and Winkelmann, Stefanie}, title = {Modelling altered signalling of G-protein coupled receptors in inflamed environment to advance drug design}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-82797}, year = {2021}, abstract = {Initiated by mathematical modelling of extracellular interactions between G-protein coupled receptors (GPCRs) and ligands in normal versus diseased (inflamed) environments, we previously reported the successful design, synthesis and testing of the prototype opioid painkiller NFEPP that does not elicit adverse side effects. Uniquely, this design recognised that GPCRs function differently under pathological versus healthy conditions. We now present a novel stochastic model of GPCR function that includes intracellular dissociation of G-protein subunits and modulation of plasma membrane calcium channels associated with parameters of inflamed tissue (pH, radicals). By means of molecular dynamics simulations, we also assessed qualitative changes of the reaction rates due to additional disulfide bridges inside the GPCR binding pocket and used these rates for stochastic simulations of the corresponding reaction jump process. The modelling results were validated with in vitro experiments measuring calcium currents and G-protein activation. We found markedly reduced G-protein dissociation and calcium channel inhibition induced by NFEPP at normal pH, and enhanced constitutive G-protein activation but lower probability of ligand binding with increasing radical concentrations. These results suggest that, compared to radicals, low pH is a more important determinant of overall GPCR function in an inflamed environment. Future drug design efforts should take this into account.}, language = {en} } @misc{ErnstSchuetteSigristetal.2021, author = {Ernst, Ariane and Sch{\"u}tte, Christof and Sigrist, Stephan and Winkelmann, Stefanie}, title = {Variance of filtered signals: Characterization for linear reaction networks and application to neurotransmission dynamics}, issn = {1438-0064}, doi = {10.1016/j.mbs.2021.108760}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-82674}, year = {2021}, abstract = {Neurotransmission at chemical synapses relies on the calcium-induced fusion of synaptic vesicles with the presynaptic membrane. The distance to the calcium channels determines the release probability and thereby the postsynaptic signal. Suitable models of the process need to capture both the mean and the variance observed in electrophysiological measurements of the postsynaptic current. In this work, we propose a method to directly compute the exact first- and second-order moments for signals generated by a linear reaction network under convolution with an impulse response function, rendering computationally expensive numerical simulations of the underlying stochastic counting process obsolete. We show that the autocorrelation of the process is central for the calculation of the filtered signal's second-order moments, and derive a system of PDEs for the cross-correlation functions (including the autocorrelations) of linear reaction networks with time-dependent rates. Finally, we employ our method to efficiently compare different spatial coarse graining approaches for a specific model of synaptic vesicle fusion. Beyond the application to neurotransmission processes, the developed theory can be applied to any linear reaction system that produces a filtered stochastic signal.}, language = {en} } @article{PlockBurgerSchneider2021, author = {Plock, Matthias and Burger, Sven and Schneider, Philipp-Immanuel}, title = {Recent advances in Bayesian optimization with applications to parameter reconstruction in optical nano-metrology}, volume = {11783}, journal = {Proc. SPIE}, arxiv = {http://arxiv.org/abs/2107.05499}, doi = {10.1117/12.2592266}, pages = {117830J}, year = {2021}, language = {en} } @misc{JaegerSutterSchneideretal.2021, author = {J{\"a}ger, Klaus and Sutter, Johannes and Schneider, Philipp-Immanuel and Hammerschmidt, Martin and Becker, Christiane}, title = {Optical simulations of nanotextured perovskite/silicon tandem solar cell}, journal = {HZB Data Service}, doi = {10.5442/ND000005}, year = {2021}, language = {en} } @misc{RaharinirinaWeberBirketal.2021, author = {Raharinirina, N. Alexia and Weber, Marcus and Birk, Ralph and Fackeldey, Konstantin and Klasse, Sarah M. and Richter, Tonio Sebastian}, title = {Different Tools and Results for Correspondence Analysis}, doi = {10.12752/8257}, year = {2021}, abstract = {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.}, language = {en} } @article{TorunSchneiderHammerschmidtetal.2021, author = {Torun, Cem G{\"u}ney and Schneider, Philipp-Immanuel and Hammerschmidt, Martin and Burger, Sven and Munns, Joseph H.D. and Schr{\"o}der, Tim}, title = {Optimized diamond inverted nanocones for enhanced color center to fiber coupling}, volume = {118}, journal = {Appl. Phys. Lett.}, arxiv = {http://arxiv.org/abs/2105.13333}, doi = {10.1063/5.0050338}, pages = {234002}, year = {2021}, language = {en} } @article{TillmannBlaesiBurgeretal.2021, author = {Tillmann, Peter and Bl{\"a}si, Benedikt and Burger, Sven and Hammerschmidt, Martin and H{\"o}hn, Oliver and Becker, Christiane and J{\"a}ger, Klaus}, title = {Optimizing metal grating back reflectors for III-V-on-silicon multijunction solar cells}, volume = {29}, journal = {Opt. Express}, doi = {10.1364/OE.426761}, pages = {22517}, year = {2021}, language = {en} } @article{GorgullaCınaroğluFischeretal.2021, author = {Gorgulla, Christoph and {\c{C}}{\i}naroğlu, S{\"u}leyman and Fischer, Patrick D. and Fackeldey, Konstantin and Wagner, Gerhard and Arthanari, Haribabu}, title = {VirtualFlow Ants—Ultra-Large Virtual Screenings with Artificial Intelligence Driven Docking Algorithm Based on Ant Colony Optimization}, volume = {22}, journal = {Special Issue Artificial Intelligence \& Deep Learning Approaches for Structural Bioinformatics}, number = {11}, doi = {https://doi.org/10.3390/ijms22115807}, pages = {5807}, year = {2021}, abstract = {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.}, language = {en} } @article{BoethBiesenHollnageletal.2021, author = {Boeth, Heide and Biesen, Robert and Hollnagel, Jeremias and Herrmann, Sandra and Ehrig, Rainald and Pelli, Lorenzo and Taylor, William R. and Duda, Georg and Buttgereit, Frank}, title = {Quantification of morning stiffness to assess disease activity and treatment effects in rheumatoid arthritis}, volume = {60}, journal = {Rheumatology}, number = {11}, doi = {10.1093/rheumatology/keab323}, pages = {5282 -- 5291}, year = {2021}, abstract = {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.}, language = {en} } @article{BittracherMoschnerKokschetal.2021, author = {Bittracher, Andreas and Moschner, Johann and Koksch, Beate and Netz, Roland and Sch{\"u}tte, Christof}, title = {Exploring the locking stage of NFGAILS amyloid fibrillation via transition manifold analysis}, volume = {94}, journal = {The European Physical Journal B}, doi = {10.1140/epjb/s10051-021-00200-0}, year = {2021}, language = {en} } @article{LiuZhouWangetal.2021, author = {Liu, Yongtao and Zhou, Zhiguang and Wang, Fan and Kewes, G{\"u}nter and Wen, Shihui and Burger, Sven and Wakiani, Majid Ebrahimi and Xi, Peng and Yang, Jiong and Yang, Xusan and Benson, Oliver and Jin, Dayong}, title = {Axial localization and tracking of self-interference nanoparticles by lateral point spread functions}, volume = {12}, journal = {Nat. Commun.}, doi = {10.1038/s41467-021-22283-0}, pages = {2019}, year = {2021}, language = {en} } @article{JaegerSutterHammerschmidtetal.2021, author = {J{\"a}ger, Klaus and Sutter, Johannes and Hammerschmidt, Martin and Schneider, Philipp-Immanuel and Becker, Christiane}, title = {Prospects of light management in perovskite/silicon tandem solar cells}, volume = {10}, journal = {Nanophotonics}, doi = {10.1515/nanoph-2020-0674}, pages = {1991}, year = {2021}, language = {en} } @article{GelssKlusSchusteretal.2021, author = {Gelss, Patrick and Klus, Stefan and Schuster, Ingmar and Sch{\"u}tte, Christof}, title = {Feature space approximation for kernel-based supervised learning}, volume = {221}, journal = {Knowledge-Based Sytems}, publisher = {Elsevier}, doi = {https://doi.org/10.1016/j.knosys.2021.106935}, year = {2021}, language = {en} } @article{delRazoDibakSchuetteetal.2021, author = {del Razo, Mauricio J. and Dibak, Manuel and Sch{\"u}tte, Christof and No{\´e}, Frank}, title = {Multiscale molecular kinetics by coupling Markov state models and reaction-diffusion dynamics}, volume = {155}, journal = {The Journal of Chemical Physics}, number = {12}, doi = {10.1063/5.0060314}, year = {2021}, language = {en} } @article{WulkowTelgmannHungenbergetal.2021, author = {Wulkow, Niklas and Telgmann, Regina and Hungenberg, Klaus-Dieter and Sch{\"u}tte, Christof and Wulkow, Michael}, title = {Deterministic and Stochastic Parameter Estimation for Polymer Reaction Kinetics I: Theory and Simple Examples}, volume = {30}, journal = {Macromolecular Theory and Simulations}, doi = {10.1002/mats.202100017}, year = {2021}, abstract = {Two different approaches to parameter estimation (PE) in the context of polymerization are introduced, refined, combined, and applied. The first is classical PE where one is interested in finding parameters which minimize the distance between the output of a chemical model and experimental data. The second is Bayesian PE allowing for quantifying parameter uncertainty caused by experimental measurement error and model imperfection. Based on detailed descriptions of motivation, theoretical background, and methodological aspects for both approaches, their relation are outlined. The main aim of this article is to show how the two approaches complement each other and can be used together to generate strong information gain regarding the model and its parameters. Both approaches and their interplay in application to polymerization reaction systems are illustrated. This is the first part in a two-article series on parameter estimation for polymer reaction kinetics with a focus on theory and methodology while in the second part a more complex example will be considered.}, language = {en} } @article{BinkowskiWuLalanneetal.2021, author = {Binkowski, Felix and Wu, Tong and Lalanne, Philippe and Burger, Sven and Govorov, Alexander O.}, title = {Hot Electron Generation through Near-Field Excitation of Plasmonic Nanoresonators}, volume = {8}, journal = {ACS Photonics}, arxiv = {http://arxiv.org/abs/2103.06652}, doi = {10.1021/acsphotonics.1c00231}, pages = {1243}, year = {2021}, language = {en} } @article{PakhomovHammerschmidtBurgeretal.2021, author = {Pakhomov, Anton V. and Hammerschmidt, Martin and Burger, Sven and Pertsch, Thomas and Setzpfand, Frank}, title = {Modeling of surface-induced second-harmonic generation from multilayer structures by the transfer matrix method}, volume = {29}, journal = {Opt. Express}, doi = {10.1364/OE.417066}, pages = {9098}, year = {2021}, language = {en} } @article{ManleySegantiniAhibozetal.2021, author = {Manley, Phillip and Segantini, Michele and Ahiboz, Doguscan and Hammerschmidt, Martin and Arnaoutakis, Georgios and MacQueen, Rowan W. and Burger, Sven and Becker, Christiane}, title = {Double-layer metasurface for enhanced photon up-conversion}, volume = {6}, journal = {APL Photon.}, arxiv = {http://arxiv.org/abs/2012.03587}, doi = {10.1063/5.0040839}, pages = {036103}, year = {2021}, language = {en} } @article{HelfmannHeitzigKoltaietal.2021, author = {Helfmann, Luzie and Heitzig, Jobst and Koltai, P{\´e}ter and Kurths, J{\"u}rgen and Sch{\"u}tte, Christof}, title = {Statistical analysis of tipping pathways in agent-based models}, volume = {230}, journal = {Eur. Phys. J. Spec. Top.}, arxiv = {http://arxiv.org/abs/2103.02883}, doi = {10.1140/epjs/s11734-021-00191-0}, pages = {3249 -- 3271}, year = {2021}, abstract = {Agent-based models are a natural choice for modeling complex social systems. In such models simple stochastic interaction rules for a large population of individuals on the microscopic scale can lead to emergent dynamics on the macroscopic scale, for instance a sudden shift of majority opinion or behavior. Here we are introducing a methodology for studying noise-induced tipping between relevant subsets of the agent state space representing characteristic configurations. Due to a large number of interacting individuals, agent-based models are high-dimensional, though usually a lower-dimensional structure of the emerging collective behaviour exists. We therefore apply Diffusion Maps, a non-linear dimension reduction technique, to reveal the intrinsic low-dimensional structure. We characterize the tipping behaviour by means of Transition Path Theory, which helps gaining a statistical understanding of the tipping paths such as their distribution, flux and rate. By systematically studying two agent-based models that exhibit a multitude of tipping pathways and cascading effects, we illustrate the practicability of our approach.}, language = {en} } @article{AndrleHoenickeGwaltetal.2021, author = {Andrle, Anna and H{\"o}nicke, Philipp and Gwalt, Grzegorz and Schneider, Philipp-Immanuel and Kayser, Yves and Siewert, Frank and Soltwisch, Victor}, title = {GIXRF and machine learning as metrological tools for shape and element sensitive reconstructions of periodic nanostructures}, volume = {11611}, journal = {Proc. SPIE}, doi = {10.1117/12.2586082}, pages = {116110R}, year = {2021}, language = {en} } @article{WinkelmannZonkerSchuetteetal.2021, author = {Winkelmann, Stefanie and Zonker, Johannes and Sch{\"u}tte, Christof and Djurdjevac Conrad, Natasa}, title = {Mathematical modeling of spatio-temporal population dynamics and application to epidemic spreading}, volume = {336}, journal = {Mathematical Biosciences}, publisher = {Elsevier}, arxiv = {http://arxiv.org/abs/2205.05000}, doi = {10.1016/j.mbs.2021.108619}, year = {2021}, abstract = {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.}, language = {en} } @article{CardereraPokuttaSchuetteetal.2021, author = {Carderera, Alejandro and Pokutta, Sebastian and Sch{\"u}tte, Christof and Weiser, Martin}, title = {An efficient first-order conditional gradient algorithm in data-driven sparse identification of nonlinear dynamics to solve sparse recovery problems under noise}, journal = {Journal of Computational and Applied Mathematics}, arxiv = {http://arxiv.org/abs/2101.02630}, year = {2021}, abstract = {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.}, language = {en} } @article{BoltzKurchanLiu2021, author = {Boltz, Horst-Holger and Kurchan, Jorge and Liu, Andrea J.}, title = {Fluctuation distributions of energy minima in complex landscapes}, volume = {3}, journal = {Physical Review Research}, number = {1}, doi = {10.1103/PhysRevResearch.3.013061}, pages = {013061}, year = {2021}, abstract = {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.}, language = {en} } @article{SchneckWeiserWende2021, author = {Schneck, Jakob and Weiser, Martin and Wende, Florian}, title = {Impact of mixed precision and storage layout on additive Schwarz smoothers}, volume = {28}, journal = {Numerical Linear Algebra with Applications}, number = {4}, doi = {10.1002/nla.2366}, year = {2021}, abstract = {The growing discrepancy between CPU computing power and memory bandwidth drives more and more numerical algorithms into a bandwidth-bound regime. One example is the overlapping Schwarz smoother, a highly effective building block for iterative multigrid solution of elliptic equations with higher order finite elements. Two options of reducing the required memory bandwidth are sparsity exploiting storage layouts and representing matrix entries with reduced precision in floating point or fixed point format. We investigate the impact of several options on storage demand and contraction rate, both analytically in the context of subspace correction methods and numerically at an example of solid mechanics. Both perspectives agree on the favourite scheme: fixed point representation of Cholesky factors in nested dissection storage.}, language = {en} } @misc{BeckerJaegerBurger2021, author = {Becker, Christiane and J{\"a}ger, Klaus and Burger, Sven}, title = {Nanophotonics for Solar Energy}, volume = {221}, journal = {Sol. Energy Mater. Sol. Cells}, doi = {10.1016/j.solmat.2020.110916}, pages = {110916}, year = {2021}, language = {en} } @article{NiemannKlusSchuette2021, author = {Niemann, Jan-Hendrik and Klus, Stefan and Sch{\"u}tte, Christof}, title = {Data-driven model reduction of agent-based systems using the Koopman generator}, volume = {16}, journal = {PLOS ONE}, number = {5}, arxiv = {http://arxiv.org/abs/2012.07718}, doi = {10.1371/journal.pone.0250970}, year = {2021}, abstract = {The dynamical behavior of social systems can be described by agent-based models. Although single agents follow easily explainable rules, complex time-evolving patterns emerge due to their interaction. The simulation and analysis of such agent-based models, however, is often prohibitively time-consuming if the number of agents is large. In this paper, we show how Koopman operator theory can be used to derive reduced models of agent-based systems using only simulation or real-world data. Our goal is to learn coarse-grained models and to represent the reduced dynamics by ordinary or stochastic differential equations. The new variables are, for instance, aggregated state variables of the agent-based model, modeling the collective behavior of larger groups or the entire population. Using benchmark problems with known coarse-grained models, we demonstrate that the obtained reduced systems are in good agreement with the analytical results, provided that the numbers of agents is sufficiently large.}, language = {en} } @article{WulkowConradDjurdjevacConradetal.2021, author = {Wulkow, Hanna and Conrad, Tim and Djurdjevac Conrad, Natasa and M{\"u}ller, Sebastian A. and Nagel, Kai and Sch{\"u}tte, Christof}, title = {Prediction of Covid-19 spreading and optimal coordination of counter-measures: From microscopic to macroscopic models to Pareto fronts}, volume = {16}, journal = {PLOS One}, number = {4}, publisher = {Public Library of Science}, doi = {10.1371/journal.pone.0249676}, year = {2021}, language = {en} } @article{BorndoerferDaneckerWeiser2021, author = {Bornd{\"o}rfer, Ralf and Danecker, Fabian and Weiser, Martin}, title = {A Discrete-Continuous Algorithm for Free Flight Planning}, volume = {14}, journal = {Algorithms}, number = {1}, publisher = {MDPI}, issn = {1438-0064}, doi = {10.3390/a14010004}, pages = {4}, year = {2021}, abstract = {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.}, language = {en} } @phdthesis{Santiago2021, author = {Santiago, Xavier Garcia}, title = {Numerical methods for shape optimization of photonic nanostructures}, organization = {Karlsruher Institut f{\"u}r Technologie}, doi = {10.5445/IR/1000131006}, year = {2021}, language = {en} } @article{MielkeMontefuscoPeletier2021, author = {Mielke, Alexander and Montefusco, Alberto and Peletier, Mark A.}, title = {Exploring families of energy-dissipation landscapes via tilting: three types of EDP convergence}, volume = {33}, journal = {Continuum Mechanics and Thermodynamics}, publisher = {Springer}, doi = {10.1007/s00161-020-00932-x}, pages = {611 -- 637}, year = {2021}, abstract = {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.}, language = {en} } @article{SchielaStoeckleinWeiser2021, author = {Schiela, Anton and St{\"o}cklein, Matthias and Weiser, Martin}, title = {A primal dual projection algorithm for efficient constraint preconditioning}, volume = {43}, journal = {SIAM Journal on Scientific Computing}, number = {6}, doi = {10.1137/20M1380739}, pages = {A4095 -- A4120}, year = {2021}, abstract = {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.}, language = {en} } @article{CvetkovicConradLie2021, author = {Cvetkovic, Nada and Conrad, Tim and Lie, Han Cheng}, title = {A Convergent Discretisation Method for Transition Path Theory for Diffusion Processes}, volume = {19}, journal = {Multiscale Modeling \& Simulation}, number = {1}, publisher = {Society for Industrial and Applied Mathematics}, doi = {10.1137/20M1329354}, pages = {242 -- 266}, year = {2021}, language = {en} } @article{JaegerTillmannKatzetal.2021, author = {J{\"a}ger, Klaus and Tillmann, Peter and Katz, Eugene A. and Becker, Christiane}, title = {Perovskite/silicon tandem solar cells: Effect of luminescent coupling and bifaciality}, volume = {5}, journal = {Sol. RRL}, arxiv = {http://arxiv.org/abs/2008.13696}, doi = {10.1002/solr.202000628}, pages = {2000628}, year = {2021}, language = {en} } @article{CheginiKopanicakovaKrauseetal.2021, author = {Chegini, Fatemeh and Kopanicakova, Alena and Krause, Rolf and Weiser, Martin}, title = {Efficient Identification of Scars using Heterogeneous Model Hierarchies}, volume = {23}, journal = {EP Europace}, doi = {10.1093/europace/euaa402}, pages = {i113 -- i122}, year = {2021}, abstract = {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.}, language = {en} } @article{MontefuscoPeletierOettinger2021, author = {Montefusco, Alberto and Peletier, Mark A. and {\"O}ttinger, Hans Christian}, title = {A Framework of Nonequilibrium Statistical Mechanics. II. Coarse-Graining}, volume = {46}, journal = {Journal of Non-Equilibrium Thermodynamics}, number = {1}, publisher = {De Gruyter}, doi = {10.1515/jnet-2020-0069}, pages = {15 -- 33}, year = {2021}, abstract = {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.}, language = {en} } @article{OettingerMontefuscoPeletier2021, author = {{\"O}ttinger, Hans Christian and Montefusco, Alberto and Peletier, Mark A.}, title = {A Framework of Nonequilibrium Statistical Mechanics. I. Role and Types of Fluctuations}, volume = {46}, journal = {Journal of Non-Equilibrium Thermodynamics}, number = {1}, publisher = {De Gruyter}, doi = {10.1515/jnet-2020-0068}, pages = {1 -- 13}, year = {2021}, abstract = {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.}, language = {en} } @article{SechiSikorskiWeber2021, author = {Sechi, Renata and Sikorski, Alexander and Weber, Marcus}, title = {Estimation of the Koopman Generator by Newton's Extrapolation}, volume = {19}, journal = {Multiscale Modeling and Simulation}, number = {2}, publisher = {SIAM}, doi = {10.1137/20M1333006}, pages = {758 -- 774}, year = {2021}, abstract = {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.}, language = {en} } @article{SikorskiWeberSchuette2021, author = {Sikorski, Alexander and Weber, Marcus and Sch{\"u}tte, Christof}, title = {The Augmented Jump Chain}, volume = {4}, journal = {Advanced Theory and Simulations}, number = {4}, publisher = {Wiley-VCH}, arxiv = {http://arxiv.org/abs/2008.04624}, doi = {10.1002/adts.202000274}, year = {2021}, abstract = {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.}, language = {en} } @article{MartensBinkowskiNguyenetal.2021, author = {Martens, Kevin and Binkowski, Felix and Nguyen, Linh and Hu, Li and Govorov, Alexander O. and Burger, Sven and Liedl, Tim}, title = {Long- and Short-Ranged Chiral Interactions in DNA Assembled Plasmonic Chains}, volume = {12}, journal = {Nat. Commun.}, arxiv = {http://arxiv.org/abs/2010.12072}, doi = {10.1038/s41467-021-22289-8}, pages = {2025}, year = {2021}, language = {en} } @article{SharmaZhang2021, author = {Sharma, Upanshu and Zhang, Wei}, title = {Non-reversible sampling schemes on submanifolds}, volume = {59}, journal = {SIAM Journal on Numerical Analysis}, number = {6}, arxiv = {http://arxiv.org/abs/2011.02835}, doi = {10.1137/20M1378752}, pages = {2989 -- 3031}, year = {2021}, abstract = {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.}, language = {en} } @article{FischerEhrigSchaeferetal.2021, author = {Fischer, Sophie and Ehrig, Rainald and Sch{\"a}fer, Stefan and Tronci, Enrico and Mancini, Toni and Egli, Marcel and Ille, Fabian and Kr{\"u}ger, Tillmann H. C. and Leeners, Brigitte and R{\"o}blitz, Susanna}, title = {Mathematical Modelling and Simulation Provides Evidence for New Strategies of Ovarian Stimulation}, volume = {12}, journal = {Frontiers in Endocrinology}, doi = {10.3389/fendo.2021.613048}, year = {2021}, abstract = {New approaches to ovarian stimulation protocols, such as luteal start, random start or double stimulation, allow for flexibility in ovarian stimulation at different phases of the menstrual cycle which is especially useful when time for assisted reproductive technology is limited, e.g. for emergency fertility preservation in cancer patients. It has been proposed that the success of these methods is based on the continuous growth of multiple cohorts ("waves") of follicles throughout the menstrual cycle which leads to the availability of ovarian follicles for ovarian controlled stimulation at several time points. Though several preliminary studies have been published, their scientific evidence has not been considered as being strong enough to integrate these results into routine clinical practice. This work aims at adding further scientific evidence about the efficiency of variable-start protocols and underpinning the theory of follicular waves by using mathematical modelling and numerical simulations. For this purpose, we have modified and coupled two previously published models, one describing the time course of hormones and one describing competitive follicular growth in a normal menstrual cycle. The coupled model is used to test stimulation protocols in silico. Simulation results show the occurrence of follicles in a wave-like manner during a normal menstrual cycle and qualitatively predict the outcome of ovarian stimulation initiated at different time points of the menstrual cycle.}, language = {en} } @article{MironBeronVeraHelfmannetal.2021, author = {Miron, Philippe and Beron-Vera, Francisco and Helfmann, Luzie and Koltai, P{\´e}ter}, title = {Transition paths of marine debris and the stability of the garbage patches}, volume = {31}, journal = {Chaos: An Interdisciplinary Journal of Nonlinear Science}, number = {3}, arxiv = {http://arxiv.org/abs/2009.11234}, doi = {https://doi.org/10.1063/5.0030535}, year = {2021}, abstract = {We used transition path theory (TPT) to infer "reactive" pathways of floating marine debris trajectories. The TPT analysis was applied on a pollution-aware time-homogeneous Markov chain model constructed from trajectories produced by satellite-tracked undrogued buoys from the NOAA Global Drifter Program. The latter involved coping with the openness of the system in physical space, which further required an adaptation of the standard TPT setting. Directly connecting pollution sources along coastlines with garbage patches of varied strengths, the unveiled reactive pollution routes represent alternative targets for ocean cleanup efforts. Among our specific findings we highlight: constraining a highly probable pollution source for the Great Pacific Garbage Patch; characterizing the weakness of the Indian Ocean gyre as a trap for plastic waste; and unveiling a tendency of the subtropical gyres to export garbage toward the coastlines rather than to other gyres in the event of anomalously intense winds.}, language = {en} } @article{KlebanovSprungkSullivan2021, author = {Klebanov, Ilja and Sprungk, Bj{\"o}rn and Sullivan, T. J.}, title = {The linear conditional expectation in Hilbert space}, volume = {27}, journal = {Bernoulli}, number = {4}, arxiv = {http://arxiv.org/abs/2008.12070}, doi = {10.3150/20-BEJ1308}, pages = {2299 -- 2299}, year = {2021}, language = {en} } @article{BoltzSirbuStelzeretal.2021, author = {Boltz, Horst-Holger and Sirbu, Alexei and Stelzer, Nina and Lohse, Martin J. and Sch{\"u}tte, Christof and Annibale, Paolo}, title = {Quantitative spectroscopy of single molecule interaction times}, volume = {46}, journal = {Optic Letters}, number = {7}, doi = {10.1364/OL.413030}, pages = {1538 -- 1541}, year = {2021}, abstract = {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.}, language = {en} } @article{SantiagoBurgerRockstuhletal.2021, author = {Santiago, Xavier Garcia and Burger, Sven and Rockstuhl, Carsten and Schneider, Philipp-Immanuel}, title = {Bayesian optimization with improved scalability and derivative information for efficient design of nanophotonic structures}, volume = {39}, journal = {J. Light. Technol.}, arxiv = {http://arxiv.org/abs/2101.02972}, doi = {10.1109/JLT.2020.3023450}, pages = {167}, year = {2021}, language = {en} } @article{AlhaddadFoerstnerGrothetal.2021, author = {Alhaddad, Samer and F{\"o}rstner, Jens and Groth, Stefan and Gr{\"u}newald, Daniel and Grynko, Yevgen and Hannig, Frank and Kenter, Tobias and Pfreundt, Franz-Josef and Plessl, Christian and Schotte, Merlind and Steinke, Thomas and Teich, J{\"u}rgen and Weiser, Martin and Wende, Florian}, title = {HighPerMeshes - A Domain-Specific Language for Numerical Algorithms on Unstructured Grids}, journal = {Euro-Par 2020: Parallel Processing Workshops.}, publisher = {Springer}, doi = {10.1007/978-3-030-71593-9_15}, pages = {185 -- 196}, year = {2021}, abstract = {Solving partial differential equations on unstructured grids is a cornerstone of engineering and scientific computing. Nowadays, heterogeneous parallel platforms with CPUs, GPUs, and FPGAs enable energy-efficient and computationally demanding simulations. We developed the HighPerMeshes C++-embedded Domain-Specific Language (DSL) for bridging the abstraction gap between the mathematical and algorithmic formulation of mesh-based algorithms for PDE problems on the one hand and an increasing number of heterogeneous platforms with their different parallel programming and runtime models on the other hand. Thus, the HighPerMeshes DSL aims at higher productivity in the code development process for multiple target platforms. We introduce the concepts as well as the basic structure of the HighPer-Meshes DSL, and demonstrate its usage with three examples, a Poisson and monodomain problem, respectively, solved by the continuous finite element method, and the discontinuous Galerkin method for Maxwell's equation. The mapping of the abstract algorithmic description onto parallel hardware, including distributed memory compute clusters is presented. Finally, the achievable performance and scalability are demonstrated for a typical example problem on a multi-core CPU cluster.}, language = {en} } @article{GorgullaDasLeighetal.2021, author = {Gorgulla, Christoph and Das, Krishna M. Padmanabha and Leigh, Kendra E and Cespugli, Marco and Fischer, Patrick D. and Wang, Zi-Fu and Tesseyre, Guilhem and Pandita, Shreya and Shnapir, Alex and Calderaio, Anthony and Hutcheson, Colin and Gechev, Minko and Rose, Alexander and Lewis, Noam and Yaffe, Erez and Luxenburg, Roni and Herce, Henry D. and Durmaz, Vedat and Halazonetis, Thanos D. and Fackeldey, Konstantin and Patten, Justin J. and Chuprina, Alexander and Dziuba, Igor and Plekhova, Alla and Moroz, Yurii and Radchenko, Dmytro and Tarkhanova, Olga and Yavnyuk, Irina and Gruber, Christian C. and Yust, Ryan and Payne, Dave and N{\"a}{\"a}r, Anders M. and Namchuk, Mark N. and Davey, Robert A. and Wagner, Gerhard and Kinney, Jamie and Arthanari, Haribabu}, title = {A Multi-Pronged Approach Targeting SARS-CoV-2 Proteins Using Ultra-Large Virtual Screening}, volume = {24}, journal = {iScience}, number = {2}, publisher = {CellPress}, doi = {10.26434/chemrxiv.12682316}, pages = {102021}, year = {2021}, abstract = {Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), previously known as 2019 novel coronavirus (2019-nCoV), has spread rapidly across the globe, creating an unparalleled global health burden and spurring a deepening economic crisis. As of July 7th, 2020, almost seven months into the outbreak, there are no approved vaccines and few treatments available. Developing drugs that target multiple points in the viral life cycle could serve as a strategy to tackle the current as well as future coronavirus pandemics. Here we leverage the power of our recently developed in silico screening platform, VirtualFlow, to identify inhibitors that target SARS-CoV-2. VirtualFlow is able to efficiently harness the power of computing clusters and cloud-based computing platforms to carry out ultra-large scale virtual screens. In this unprecedented structure-based multi-target virtual screening campaign, we have used VirtualFlow to screen an average of approximately 1 billion molecules against each of 40 different target sites on 17 different potential viral and host targets in the cloud. In addition to targeting the active sites of viral enzymes, we also target critical auxiliary sites such as functionally important protein-protein interaction interfaces. This multi-target approach not only increases the likelihood of finding a potent inhibitor, but could also help identify a collection of anti-coronavirus drugs that would retain efficacy in the face of viral mutation. Drugs belonging to different regimen classes could be combined to develop possible combination therapies, and top hits that bind at highly conserved sites would be potential candidates for further development as coronavirus drugs. Here, we present the top 200 in silico hits for each target site. While in-house experimental validation of some of these compounds is currently underway, we want to make this array of potential inhibitor candidates available to researchers worldwide in consideration of the pressing need for fast-tracked drug development.}, language = {en} } @article{RoehlWeberFackeldey2021, author = {R{\"o}hl, Susanne and Weber, Marcus and Fackeldey, Konstantin}, title = {Computing the minimal rebinding effect for non-reversible processes}, volume = {19}, journal = {Multiscale Modeling and Simulation}, number = {1}, arxiv = {http://arxiv.org/abs/2007.08403}, doi = {https://doi.org/10.1137/20M1334966}, pages = {460 -- 477}, year = {2021}, abstract = {The aim of this paper is to investigate the rebinding effect, a phenomenon describing a "short-time memory" which can occur when projecting a Markov process onto a smaller state space. For guaranteeing a correct mapping by the Markov State Model, we assume a fuzzy clustering in terms of membership functions, assigning degrees of membership to each state. The macro states are represented by the membership functions and may be overlapping. The magnitude of this overlap is a measure for the strength of the rebinding effect, caused by the projection and stabilizing the system. A minimal bound for the rebinding effect included in a given system is computed as the solution of an optimization problem. Based on membership functions chosen as a linear combination of Schur vectors, this generalized approach includes reversible as well as non-reversible processes.}, language = {en} } @article{FackeldeyRoehmNiknejadetal.2021, author = {Fackeldey, Konstantin and R{\"o}hm, Jonas and Niknejad, Amir and Chewle, Surahit and Weber, Marcus}, title = {Analyzing Raman Spectral Data without Separabiliy Assumption}, volume = {3}, journal = {Journal of Mathematical Chemistry}, number = {59}, publisher = {Springer}, arxiv = {http://arxiv.org/abs/2007.06428}, doi = {10.1007/s10910-020-01201-7}, pages = {575 -- 596}, year = {2021}, abstract = {Raman spectroscopy is a well established tool for the analysis of vibration spectra, which then allow for the determination of individual substances in a chemical sample, or for their phase transitions. In the Time-Resolved-Raman-Sprectroscopy the vibration spectra of a chemical sample are recorded sequentially over a time interval, such that conclusions for intermediate products (transients) can be drawn within a chemical process. The observed data-matrix M from a Raman spectroscopy can be regarded as a matrix product of two unknown matrices W and H, where the first is representing the contribution of the spectra and the latter represents the chemical spectra. One approach for obtaining W and H is the non-negative matrix factorization. We propose a novel approach, which does not need the commonly used separability assumption. The performance of this approach is shown on a real world chemical example.}, language = {en} } @article{WulkowKoltaiSchuette2021, author = {Wulkow, Niklas and Koltai, P{\´e}ter and Sch{\"u}tte, Christof}, title = {Memory-Based Reduced Modelling and Data-Based Estimation of Opinion Spreading}, volume = {31}, journal = {Journal of Nonlinear Science}, doi = {10.1007/s00332-020-09673-2}, year = {2021}, abstract = {We investigate opinion dynamics based on an agent-based model and are interested in predicting the evolution of the percentages of the entire agent population that share an opinion. Since these opinion percentages can be seen as an aggregated observation of the full system state, the individual opinions of each agent, we view this in the framework of the Mori-Zwanzig projection formalism. More specifically, we show how to estimate a nonlinear autoregressive model (NAR) with memory from data given by a time series of opinion percentages, and discuss its prediction capacities for various specific topologies of the agent interaction network. We demonstrate that the inclusion of memory terms significantly improves the prediction quality on examples with different network topologies.}, language = {en} } @article{Zhang2021, author = {Zhang, Wei}, title = {Some new results on relative entropy production, time reversal, and optimal control of time-inhomogeneous diffusion processes}, volume = {62}, journal = {Journal of Mathematical Physics}, number = {4}, arxiv = {http://arxiv.org/abs/2006.11212}, doi = {10.1063/5.0038740}, pages = {26}, year = {2021}, abstract = {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.}, language = {en} } @article{KostreSchuetteNoeetal.2021, author = {Kostre, Margarita and Sch{\"u}tte, Christof and No{\´e}, Frank and del Razo Sarmina, Mauricio}, title = {Coupling Particle-Based Reaction-Diffusion Simulations with Reservoirs Mediated by Reaction-Diffusion PDEs}, volume = {19}, journal = {Multiscale Modeling \& Simulation}, number = {4}, publisher = {Society for Industrial and Applied Mathematics}, arxiv = {http://arxiv.org/abs/2006.00003}, doi = {10.1137/20M1352739}, pages = {1659 -- 1683}, year = {2021}, abstract = {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.}, language = {en} } @inproceedings{KrauseWeiser2021, author = {Krause, Rolf and Weiser, Martin}, title = {Multilevel augmented Lagrangian solvers for overconstrained contact formulations}, volume = {71}, booktitle = {ESAIM: ProcS}, doi = {10.1051/proc/202171175}, pages = {175 -- 184}, year = {2021}, abstract = {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.}, language = {en} } @article{NiemannWinkelmannWolfetal.2021, author = {Niemann, Jan-Hendrik and Winkelmann, Stefanie and Wolf, Sarah and Sch{\"u}tte, Christof}, title = {Agent-based modeling: Population limits and large timescales}, volume = {31}, journal = {Chaos: An Interdisciplinary Journal of Nonlinear Science}, number = {3}, issn = {1438-0064}, doi = {10.1063/5.0031373}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-77309}, year = {2021}, abstract = {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.}, language = {en} } @article{GoetschelSchielaWeiser2021, author = {G{\"o}tschel, Sebastian and Schiela, Anton and Weiser, Martin}, title = {Kaskade 7 - a Flexible Finite Element Toolbox}, volume = {81}, journal = {Computers and Mathematics with Applications}, doi = {10.1016/j.camwa.2020.02.011}, pages = {444 -- 458}, year = {2021}, abstract = {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.}, language = {en} } @article{LieSullivanTeckentrup2021, author = {Lie, Han Cheng and Sullivan, T. J. and Teckentrup, Aretha}, title = {Error bounds for some approximate posterior measures in Bayesian inference}, journal = {Numerical Mathematics and Advanced Applications ENUMATH 2019}, publisher = {Springer}, arxiv = {http://arxiv.org/abs/1911.05669}, doi = {10.1007/978-3-030-55874-1_26}, pages = {275 -- 283}, year = {2021}, language = {en} } @article{HelfmannDjurdjevacConradDjurdjevacetal.2021, author = {Helfmann, Luzie and Djurdjevac Conrad, Natasa and Djurdjevac, Ana and Winkelmann, Stefanie and Sch{\"u}tte, Christof}, title = {From interacting agents to density-based modeling with stochastic PDEs}, volume = {16}, journal = {Communications in Applied Mathematics and Computational Science}, number = {1}, arxiv = {http://arxiv.org/abs/1905.13525}, doi = {10.2140/camcos.2021.16.1}, pages = {1 -- 32}, year = {2021}, abstract = {Many real-world processes can naturally be modeled as systems of interacting agents. However, the long-term simulation of such agent-based models is often intractable when the system becomes too large. In this paper, starting from a stochastic spatio-temporal agent-based model (ABM), we present a reduced model in terms of stochastic PDEs that describes the evolution of agent number densities for large populations. We discuss the algorithmic details of both approaches; regarding the SPDE model, we apply Finite Element discretization in space which not only ensures efficient simulation but also serves as a regularization of the SPDE. Illustrative examples for the spreading of an innovation among agents are given and used for comparing ABM and SPDE models.}, language = {en} } @article{KlebanovSikorskiSchuetteetal.2021, author = {Klebanov, Ilja and Sikorski, Alexander and Sch{\"u}tte, Christof and R{\"o}blitz, Susanna}, title = {Objective priors in the empirical Bayes framework}, volume = {48}, journal = {Scandinavian Journal of Statistics}, number = {4}, publisher = {Wiley Online Library}, doi = {10.1111/sjos.12485}, pages = {1212 -- 1233}, year = {2021}, abstract = {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.}, language = {en} }