TY - CHAP A1 - Lassila, Antti A1 - Shpak, Maksim A1 - Vaigu, A. A1 - Smid, M. A1 - Velasco, A. A1 - Pokatilov, A. A1 - Fatadin, I. A1 - Ferguson, R. A1 - Robinson, D. A1 - Walbaum, T. A1 - Burger, Sven A1 - Haapalainen, Mertsi A1 - Vahimaa, P. A1 - Lukner, H. A1 - Pernice, W. A1 - Pitwon, R. A1 - Deninger, A. A1 - Morel, Jacques A1 - Castagna, Natascia A1 - Peters, O. A1 - Husu, H. T1 - Joint research for photonics fibres, waveguides and interconnects T2 - Proc. Europ. Opt. Soc. Ann. Meet. 2016 (EOSAM) Y1 - 2016 SN - 978-1-5108-4796-5 SP - 25 ER - TY - CHAP A1 - Zschiedrich, Lin A1 - Schmidt, Frank A1 - Pomplun, Jan A1 - Burger, Sven T1 - Finite elements for nano-optics applications T2 - Proc. Europ. Opt. Soc. Ann. Meet. 2016 (EOSAM) Y1 - 2016 SN - 978-1-5108-4796-5 SP - 178 ER - TY - CHAP A1 - Barth, Carlo A1 - Roder, Sebastian A1 - Brodoceanu, Daniel A1 - Kraus, Tobias A1 - Burger, Sven A1 - Hammerschmidt, Martin A1 - Schmidt, Frank A1 - Becker, Christiane T1 - Increased fluorescence of PbS quantum dots on photonic crystal slab structures T2 - Proc. Europ. Opt. Soc. Ann. Meet. 2016 (EOSAM) Y1 - 2016 SN - 978-1-5108-4796-5 SP - 181 ER - TY - CHAP A1 - Hammerschmidt, Martin A1 - Barth, Carlo A1 - Burger, Sven A1 - Becker, Christiane A1 - Schmidt, Frank T1 - Determining 2D photonic crystal geometries from reflectance spectra with a reduced basis method T2 - Proc. Europ. Opt. Soc. Ann. Meet. 2016 (EOSAM) Y1 - 2016 SN - 978-1-5108-4796-5 SP - 281 ER - TY - GEN A1 - Winkelmann, Stefanie T1 - Markov Control with Rare State Observation: Average Optimality N2 - This paper investigates the criterion of long-term average costs for a Markov decision process (MDP) which is not permanently observable. Each observation of the process produces a fixed amount of \textit{information costs} which enter the considered performance criterion and preclude from arbitrarily frequent state testing. Choosing the \textit{rare} observation times is part of the control procedure. In contrast to the theory of partially observable Markov decision processes, we consider an arbitrary continuous-time Markov process on a finite state space without further restrictions on the dynamics or the type of interaction. Based on the original Markov control theory, we redefine the control model and the average cost criterion for the setting of information costs. We analyze the constant of average costs for the case of ergodic dynamics and present an optimality equation which characterizes the optimal choice of control actions and observation times. For this purpose, we construct an equivalent freely observable MDP and translate the well-known results from the original theory to the new setting. T3 - ZIB-Report - 16-59 KW - Markov decision process KW - partial observability KW - average optimality KW - information costs Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-60981 SN - 1438-0064 ER - TY - GEN A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - The spatiotemporal master equation: approximation of reaction-diffusion dynamics via Markov state modeling N2 - Accurate modeling and numerical simulation of reaction kinetics is a topic of steady interest. We consider the spatiotemporal chemical master equation (ST-CME) as a model for stochastic reaction-diffusion systems that exhibit properties of metastability. The space of motion is decomposed into metastable compartments and diffusive motion is approximated by jumps between these compartments. Treating these jumps as first-order reactions, simulation of the resulting stochastic system is possible by the Gillespie method. We present the theory of Markov state models (MSM) as a theoretical foundation of this intuitive approach. By means of Markov state modeling, both the number and shape of compartments and the transition rates between them can be determined. We consider the ST-CME for two reaction-diffusion systems and compare it to more detailed models. Moreover, a rigorous formal justification of the ST-CME by Galerkin projection methods is presented. T3 - ZIB-Report - 16-60 KW - reaction-diffusion KW - stochastic chemical kinetics KW - chemical master equation Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-60999 SN - 1438-0064 ER - TY - JOUR A1 - Gutsche, Philipp A1 - Mäusle, Raquel A1 - Burger, Sven T1 - Locally Enhanced and Tunable Optical Chirality in Helical Metamaterials JF - Photonics Y1 - 2016 U6 - https://doi.org/10.3390/photonics3040060 VL - 3 SP - 60 ER - TY - CHAP A1 - Gutsche, Philipp A1 - Mäusle, Raquel A1 - Burger, Sven T1 - Tailoring local optical chirality in helical metamaterials T2 - 2016 10th International Congress on Advanced Electromagnetic Materials in Microwaves and Optics Y1 - 2016 U6 - https://doi.org/10.1109/MetaMaterials.2016.7746440 SP - 73 EP - 75 ER - TY - GEN A1 - Wilson, Jon A1 - Herrmann, Sven A1 - Gutsche, Philipp A1 - Burger, Sven A1 - McPeak, Kevin T1 - The Chiral Coefficient: Rapid Optimization of Broadband Plasmonic Chirality T2 - 2016 MRS Fall Meeting & Exhibit Y1 - 2016 ER - TY - GEN A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Empirical Bayes Methods, Reference Priors, Cross Entropy and the EM Algorithm N2 - When estimating a probability density within the empirical Bayes framework, the non-parametric maximum likelihood estimate (NPMLE) usually tends to overfit the data. This issue is usually taken care of by regularization - a penalization term is subtracted from the marginal log-likelihood before the maximization step, so that the estimate favors smooth solutions, resulting in the so-called maximum penalized likelihood estimation (MPLE). The majority of penalizations currently in use are rather arbitrary brute-force solutions, which lack invariance under transformation of the parameters(reparametrization) and measurements. This contradicts the principle that, if the underlying model has several equivalent formulations, the methods of inductive inference should lead to consistent results. Motivated by this principle and using an information-theoretic point of view, we suggest an entropy-based penalization term that guarantees this kind of invariance. The resulting density estimate can be seen as a generalization of reference priors. Using the reference prior as a hyperprior, on the other hand, is argued to be a poor choice for regularization. We also present an insightful connection between the NPMLE, the cross entropy and the principle of minimum discrimination information suggesting another method of inference that contains the doubly-smoothed maximum likelihood estimation as a special case. T3 - ZIB-Report - 16-56 KW - parameter estimation KW - Bayesian inference KW - Bayesian hierarchical modeling KW - hyperparameter KW - hyperprior KW - EM algorithm KW - NPMLE KW - MPLE KW - DS-MLE KW - principle of maximum entropy KW - cross entropy KW - minimum discrimination information KW - reference prior KW - Jeffreys prior Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-61230 SN - 1438-0064 ER - TY - GEN A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Empirical Bayes Methods for Prior Estimation in Systems Medicine N2 - One of the main goals of mathematical modelling in systems medicine related to medical applications is to obtain patient-specific parameterizations and model predictions. In clinical practice, however, the number of available measurements for single patients is usually limited due to time and cost restrictions. This hampers the process of making patient-specific predictions about the outcome of a treatment. On the other hand, data are often available for many patients, in particular if extensive clinical studies have been performed. Therefore, before applying Bayes’ rule separately to the data of each patient (which is typically performed using a non-informative prior), it is meaningful to use empirical Bayes methods in order to construct an informative prior from all available data. We compare the performance of four priors - a non-informative prior and priors chosen by nonparametric maximum likelihood estimation (NPMLE), by maximum penalized lilelihood estimation (MPLE) and by doubly-smoothed maximum likelihood estimation (DS-MLE) - by applying them to a low-dimensional parameter estimation problem in a toy model as well as to a high-dimensional ODE model of the human menstrual cycle, which represents a typical example from systems biology modelling. T3 - ZIB-Report - 16-57 KW - Parameter estimation KW - Bayesian inference KW - Bayesian hierarchical modelling KW - NPMLE KW - MPLE KW - DS-MLE KW - EM algorithm KW - Jeffreys prior KW - reference prior KW - hyperparameter KW - hyperprior KW - principle of maximum entropy Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-61307 SN - 1438-0064 ER - TY - CHAP A1 - Jäger, Klaus A1 - Hammerschmidt, Martin A1 - Köppel, Grit A1 - Burger, Sven A1 - Becker, Christiane T1 - On Accurate Simulations of Thin-Film Solar Cells With a Thick Glass Superstrate T2 - Light, Energy and the Environment 2016 Y1 - 2016 U6 - https://doi.org/10.1364/PV.2016.PM3B.5 SP - PM3B.5 ER - TY - CHAP A1 - Gutsche, Philipp A1 - Mäusle, Raquel A1 - Burger, Sven T1 - Circular polarization phenomena in chiral nano-optical devices T2 - Light, Energy and the Environment 2016 Y1 - 2016 U6 - https://doi.org/10.1364/FTS.2016.JW4A.4 SP - JW4A.4 ER - TY - GEN A1 - Sali, Adrian T1 - Coupling of Monodomain and Eikonal Models for Cardiac Electrophysiology N2 - The primary goal of this paper is to study the coupling of monodomain and eikonal models for the numerical simulation of cardiac electrophysiology. Eikonal models are nonlinear elliptic equations describing the excitation time of the cardiac tissue. They are often used as very fast approximations for monodomain or bidomain models - parabolic reaction-diffusion systems describing the excitation wavefront in terms of ionic currents. The excitation front is a thin region with high gradients, whereas excitation times vary over larger domains. Hence, eikonal equations can be solved on much coarser grids than monodomain equations. Moreover, as eikonal models are not time-dependent, no time integration is needed. Eikonal models are derived from monodomain models making additional assumptions and using certain approximations. While generally the approximation is rather good, several specific situations are not well captured by eikonal models. We consider coupling the two models, i.e. using the monodomain model in regions where more accurate results or the shape of the wavefront are needed, and the eikonal model in the remaining parts of the domain, where the excitation time is sufficient. Restricting the monodomain simulation to a small subdomain reduces the computational effort considerably. Numerical methods for the simulation of the individual models are presented, with the finite element method as the main ingredient. Coupling conditions as well as algorithms for implementing the coupling are explained. The approximation quality and efficiency of the coupled model is illustrated on simple geometries using an Aliev-Panfilov membrane model. T3 - ZIB-Report - 16-50 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-60508 SN - 1438-0064 ER - TY - THES A1 - Sali, Adrian T1 - Coupling of Monodomain and Eikonal Models for Cardiac Electrophysiology N2 - The primary goal of this paper is to study the coupling of monodomain and eikonal models for the numerical simulation of cardiac electrophysiology. Eikonal models are nonlinear elliptic equations describing the excitation time of the cardiac tissue. They are often used as very fast approximations for monodomain or bidomain models - parabolic reaction-diffusion systems describing the excitation wavefront in terms of ionic currents. The excitation front is a thin region with high gradients, whereas excitation times vary over larger domains. Hence, eikonal equations can be solved on much coarser grids than monodomain equations. Moreover, as eikonal models are not time-dependent, no time integration is needed. Eikonal models are derived from monodomain models making additional assumptions and using certain approximations. While generally the approximation is rather good, several specific situations are not well captured by eikonal models. We consider coupling the two models, i.e. using the monodomain model in regions where more accurate results or the shape of the wavefront are needed, and the eikonal model in the remaining parts of the domain, where the excitation time is sufficient. Restricting the monodomain simulation to a small subdomain reduces the computational effort considerably. Numerical methods for the simulation of the individual models are presented, with the finite element method as the main ingredient. Coupling conditions as well as algorithms for implementing the coupling are explained. The approximation quality and efficiency of the coupled model is illustrated on simple geometries using an Aliev-Panfilov membrane model. Y1 - 2016 ER - TY - GEN A1 - de Lasson, Jakob Rosenkrantz A1 - Frandsen, Lars Hagedorn A1 - Burger, Sven A1 - Gutsche, Philipp A1 - Kim, Oleksiy S. A1 - Breinbjerg, Olav A1 - Sigmund, Ole A1 - Mork, Jesper A1 - Gregersen, Niels T1 - Comparison of four computational methods for computing Q factors and resonance wavelengths in photonic crystal membrane cavities T2 - Proceedings of the 7th International Conference on Metamaterials, Photonic Crystals and Plasmonics Y1 - 2016 UR - http://jakobrdl.dk/presentationfiles/2016/META16Poster_Niels_July2016.pdf ER - TY - JOUR A1 - Villatoro, José A1 - Zühlke, Martin A1 - Riebe, Daniel A1 - Beitz, Toralf A1 - Weber, Marcus A1 - Riedel, Jens A1 - Löhmannsröben, Hans-Gerd T1 - IR-MALDI ion mobility spectrometry: physical source characterization and application as HPLC detector JF - International Journal for Ion Mobility Spectrometry Y1 - 2016 U6 - https://doi.org/10.1007/s12127-016-0208-1 IS - 19/4 SP - 197 EP - 297 ER - TY - JOUR A1 - Lange, Alexander T1 - Reconstruction of disease transmission rates: applications to measles, dengue, and influenza JF - Journal of Theoretical Biology Y1 - 2016 U6 - https://doi.org/10.1016/j.jtbi.2016.04.017 VL - 400 SP - 138 EP - 153 ER - TY - CHAP A1 - Lange, Alexander A1 - Plöntzke, Julia A1 - Schäfer, Stefan A1 - Röblitz, Susanna T1 - Follicular maturation in cows: mathematical models and data T2 - 10. European Conference on Mathematical and Theoretical Biology Y1 - 2016 ER - TY - GEN A1 - Omari, Mohamed A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Röblitz, Susanna T1 - Modelling of glucose-insulin metabolism and its effect on the estrous cycle in bovine T2 - Annual Meeting of the Population Approach Group in Europe Y1 - 2016 ER - TY - JOUR A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - The Spatiotemporal Master Equation: Approximation of Reaction-Diffusion Dynamics via Markov State Modeling JF - Journal of Chemical Physics N2 - Accurate modeling and numerical simulation of reaction kinetics is a topic of steady interest.We consider the spatiotemporal chemical master equation (ST-CME) as a model for stochastic reaction-diffusion systems that exhibit properties of metastability. The space of motion is decomposed into metastable compartments and diffusive motion is approximated by jumps between these compartments. Treating these jumps as first-order reactions, simulation of the resulting stochastic system is possible by the Gillespie method. We present the theory of Markov state models (MSM) as a theoretical foundation of this intuitive approach. By means of Markov state modeling, both the number and shape of compartments and the transition rates between them can be determined. We consider the ST-CME for two reaction-diffusion systems and compare it to more detailed models. Moreover, a rigorous formal justification of the ST-CME by Galerkin projection methods is presented. Y1 - 2016 U6 - https://doi.org/10.1063/1.4971163 VL - 145 IS - 21 ER - TY - JOUR A1 - Abass, Aimi A1 - Gutsche, Philipp A1 - Maes, Bjorn A1 - Rockstuhl, Carsten A1 - Martins, Emiliano R T1 - Insights into directional scattering: from coupled dipoles to asymmetric dimer nanoantennas JF - Opt. Express N2 - Strong and directionally specific forward scattering from optical nanoantennas is of utmost importance for various applications in the broader context of photovoltaics and integrated light sources. Here, we outline a simple yet powerful design principle to perceive a nanoantenna that provides directional scattering into a higher index substrate based on the interference of multiple electric dipoles. A structural implementation of the electric dipole distribution is possible using plasmonic nanoparticles with a fairly simple geometry, i.e. two coupled rectangular nanoparticles, forming a dimer, on top of a substrate. The key to achieve directionality is to choose a sufficiently large size for the nanoparticles. This promotes the excitation of vertical electric dipole moments due to the bi-anisotropy of the nanoantenna. In turn, asymmetric scattering is obtained by ensuring the appropriate phase relation between the vertical electric dipole moments. The scattering strength and angular spread for an optimized nanoantenna can be shown to be broadband and robust against changes in the incidence angle. The scattering directionality is maintained even for an array configuration of the dimer. It only requires the preferred scattering direction of the isolated nanoantenna not to be prohibited by interference. Y1 - 2016 U6 - https://doi.org/10.1364/OE.24.019638 VL - 24 IS - 17 SP - 19638 EP - 19650 ER - TY - GEN A1 - Lie, Han Cheng A1 - Sullivan, T. J. T1 - Cameron--Martin theorems for sequences of Cauchy-distributed random variables N2 - Given a sequence of Cauchy-distributed random variables defined by a sequence of location parameters and a sequence of scale parameters, we consider another sequence of random variables that is obtained by perturbing the location or scale parameter sequences. Using a result of Kakutani on equivalence of infinite product measures, we provide sufficient conditions for the equivalence of laws of the two sequences. T3 - ZIB-Report - 16-40 Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-60230 SN - 1438-0064 ER - TY - GEN A1 - Sullivan, T. J. T1 - Well-posed Bayesian inverse problems and heavy-tailed stable Banach space priors N2 - This article extends the framework of Bayesian inverse problems in infinite-dimensional parameter spaces, as advocated by Stuart (Acta Numer. 19:451–559, 2010) and others, to the case of a heavy-tailed prior measure in the family of stable distributions, such as an infinite-dimensional Cauchy distribution, for which polynomial moments are infinite or undefined. It is shown that analogues of the Karhunen–Loève expansion for square-integrable random variables can be used to sample such measures. Furthermore, under weaker regularity assumptions than those used to date, the Bayesian posterior measure is shown to depend Lipschitz continuously in the Hellinger metric upon perturbations of the misfit function and observed data. T3 - ZIB-Report - 16-30 KW - Bayesian inverse problems KW - heavy-tailed distribution KW - Karhunen–Loève expansion KW - stable distribution KW - uncertainty quantification KW - well-posedness Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-59422 SN - 1438-0064 ER - TY - GEN A1 - Zhang, Wei A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Effective Dynamics Along Given Reaction Coordinates, and Reaction Rate Theory N2 - In molecular dynamics and related fields one considers dynamical descriptions of complex systems in full (atomic) detail. In order to reduce the overwhelming complexity of realistic systems (high dimension, large timescale spread, limited computational resources) the projection of the full dynamics onto some reaction coordinates is examined in order to extract statistical information like free energies or reaction rates. In this context, the effective dynamics that is induced by the full dynamics on the reaction coordinate space has attracted considerable attention in the literature. In this article, we contribute to this discussion: We first show that if we start with an ergodic diffusion processes whose invariant measure is unique then these properties are inherited by the effective dynamics. Then, we give equations for the effective dynamics, discuss whether the dominant timescales and reaction rates inferred from the effective dynamics are accurate approximations of such quantities for the full dynamics, and compare our findings to results from approaches like Zwanzig-Mori, averaging, or homogenization. Finally, by discussing the algorithmic realization of the effective dynamics, we demonstrate that recent algorithmic techniques like the ”equation-free” approach and the ”heterogeneous multiscale method” can be seen as special cases of our approach. T3 - ZIB-Report - 16-35 KW - Ergodic diffusion KW - reaction coordinate KW - effective dynamics KW - model reduction KW - equation-free KW - heterogeneous multiscale method Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-59706 SN - 1438-0064 ER - TY - CHAP A1 - Schäfer, Stefan A1 - Plöntzke, Julia A1 - Röblitz, Susanna T1 - Mathematical Modelling of Follicular Maturation in Cows and Women T2 - 49. Jahrestagung der Physiologie und Pathologie der Fortpflanzung und gleichzeitig 41. Veterinär-Humanmedizinische Gemeinschaftstagung 2016, Leipzig. Y1 - 2016 ER - TY - GEN A1 - Cockayne, Jon A1 - Oates, Chris A1 - Sullivan, T. J. A1 - Girolami, Mark T1 - Probabilistic Meshless Methods for Partial Differential Equations and Bayesian Inverse Problems N2 - This paper develops a class of meshless methods that are well-suited to statistical inverse problems involving partial differential equations (PDEs). The methods discussed in this paper view the forcing term in the PDE as a random field that induces a probability distribution over the residual error of a symmetric collocation method. This construction enables the solution of challenging inverse problems while accounting, in a rigorous way, for the impact of the discretisation of the forward problem. In particular, this confers robustness to failure of meshless methods, with statistical inferences driven to be more conservative in the presence of significant solver error. In addition, (i) a principled learning-theoretic approach to minimise the impact of solver error is developed, and (ii) the challenging setting of inverse problems with a non-linear forward model is considered. The method is applied to parameter inference problems in which non-negligible solver error must be accounted for in order to draw valid statistical conclusions. T3 - ZIB-Report - 16-31 KW - Probabilistic Numerics KW - Partial Differential Equations KW - Inverse Problems KW - Meshless Methods KW - Gaussian Processes KW - Pseudo-Marginal MCMC Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-59513 SN - 1438-0064 ER - TY - CHAP A1 - Jäger, Klaus A1 - Barth, Carlo A1 - Hammerschmidt, Martin A1 - Herrmann, Sven A1 - Burger, Sven A1 - Schmidt, Frank A1 - Becker, Christiane ED - European Society for Quantum Solar Energy Conversion, T1 - Sinusoidal Nanotextures for Enhanced Light Management in Thin-Film Solar Cells T2 - 28th Workshop on Quantum Solar Energy Conversion - (QUANTSOL) Y1 - 2016 ER - TY - JOUR A1 - Barth, Carlo A1 - Burger, Sven A1 - Becker, Christiane T1 - Symmetry-Dependency of Anticrossing Phenomena in Slab-Type Photonic Crystals JF - Opt. Express Y1 - 2016 U6 - https://doi.org/10.1364/OE.24.010931 VL - 24 SP - 10931 ER - TY - JOUR A1 - Barth, Carlo A1 - Probst, Jürgen A1 - Herrmann, Sven A1 - Hammerschmidt, Martin A1 - Becker, Christiane T1 - Numerical characterization of symmetry properties for photonic crystals with hexagonal lattice JF - Proc. SPIE N2 - We present a numerical method to characterize the symmetry properties of photonic crystal (PhC) modes based on field distributions, which themselves can be obtained numerically. These properties can be used to forecast specific features of the optical response of such systems, e.g. which modes are allowed to couple to external radiation fields. We use 2D PhCs with a hexagonal lattice of holes in dielectric as an example and apply our technique to reproduce results from analytical considerations. Further, the method is extended to fully vectorial problems in view of 3D PhCs and PhC slabs, its functionality is demonstrated using test cases and, finally, we provide an efficient implementation. The technique can thus readily be applied to output data of all band structure computation methods or even be embedded – gaining additional information about the mode symmetry. Y1 - 2016 U6 - https://doi.org/10.1117/12.2227094 VL - 9885 SP - 988506 ER - TY - JOUR A1 - Jäger, Klaus A1 - Köppel, Grit A1 - Barth, Carlo A1 - Hammerschmidt, Martin A1 - Herrmann, Sven A1 - Burger, Sven A1 - Schmidt, Frank A1 - Becker, Christiane T1 - Sinusoidal gratings for optimized light management in c-Si thin-film solar cells JF - Proc. SPIE Y1 - 2016 U6 - https://doi.org/10.1117/12.2225459 VL - 9898 SP - 989808 ER - TY - JOUR A1 - Soltwisch, Victor A1 - Haase, Anton A1 - Wernecke, Jan A1 - Probst, Jürgen A1 - Schoengen, Max A1 - Burger, Sven A1 - Krumrey, Michael A1 - Scholze, Frank T1 - Correlated diffuse x-ray scattering from periodically nanostructured surfaces JF - Phys. Rev. B Y1 - 2016 U6 - https://doi.org/10.1103/PhysRevB.94.035419 VL - 94 SP - 035419 ER - TY - THES A1 - Hammerschmidt, Martin T1 - Optical simulation of complex nanostructured solar cells with a reduced basis method N2 - Simulations of optical processes and complex nanostructured devices have become omnipresent in recent years in several fields of current research and industrial applications, not limited to the field of photovoltaics. Devices or processes are optimized with respect to a certain objective where the underlying physical processes are described by partial differential equations. In photovoltaics and photonics electromagnetic fields are investigated which are governed by Maxwell’s equations. In this thesis a reduced basis method for the solution of the parameter dependent electromagnetic scattering problem with arbitrary parameters is developed. The method is developed with the specific challenges arising in optical simulations of thin-film silicon solar cells in mind. These are large in domain size and have a complex three-dimensional structure, making optimization tasks infeasible if high-accuracy of the electromagnetic field solution is required. The application of the empirical interpolation methods allows to expand an arbitrary parameter dependence affinely. Thus not only geometries, but also material tensors and source fields can be parameterized. Additionally, the required non-linear post-processing steps of the electromagnetic field to derive energy fluxes or volume absorption are addressed. The reduced basis method allows to reduce the computational costs by orders of magnitude compared to efficient finite element solvers. In addition, an efficient tailored domain decomposition algorithm is presented to model incoherent layers or illuminations in optical systems efficiently. This is of particular interest for solar cells in superstrate configuration where the absorber is illuminated through a glass substrate. The developed methods are employed in application examples taken from collaborations with experimentalists active in the joint lab “BerOSE” (Berlin Joint Lab for Optical Simulations for Energy Research). The optical model of a thin-film silicon multi-junction with incoherent light-trapping is characterized in great detail. The computational gains through hybrid, hp adaptive finite elements are studied and the incoherent domain decomposition algorithm is applied to model a more realistic light-trapping by the glass substrate. The numerical examples of a hexagonal nano-hole array and multi-junction silicon solar cell with a tunable intermediate reflector layer show that the reduced basis method is well suited as a forward solver for modeling and optimization tasks arising in photovoltaics and photonics. Reduced models for illumination and geometric parameters are built providing up to five orders of magnitude savings in computational costs. Resonance phenomena present in the nano-hole array example are detected and the model adapts itself automatically. KW - reduced basis method; model reduction; electromagnetics; photonics; solar cells; nanooptics; Maxwells equations; incoherence Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:188-fudissthesis000000102429-1 UR - http://www.diss.fu-berlin.de/diss/receive/FUDISS_thesis_000000102429 ER - TY - GEN A1 - Nielsen, Adam T1 - The Monte Carlo Computation Error of Transition Probabilities N2 - In many applications one is interested to compute transition probabilities of a Markov chain. This can be achieved by using Monte Carlo methods with local or global sampling points. In this article, we analyze the error by the difference in the $L^2$ norm between the true transition probabilities and the approximation achieved through a Monte Carlo method. We give a formula for the error for Markov chains with locally computed sampling points. Further, in the case of reversible Markov chains, we will deduce a formula for the error when sampling points are computed globally. We will see that in both cases the error itself can be approximated with Monte Carlo methods. As a consequence of the result, we will derive surprising properties of reversible Markov chains. T3 - ZIB-Report - 16-37 KW - Reversible Markov chain KW - Monte Carlo methods KW - Computation error KW - Measurable state space KW - Markov operator Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-59933 SN - 1438-0064 ER - TY - JOUR A1 - Nielsen, Adam T1 - The Monte Carlo Computation Error of Transition Probabilities JF - Statistics & Probability Letters N2 - In many applications one is interested to compute transition probabilities of a Markov chain. This can be achieved by using Monte Carlo methods with local or global sampling points. In this article, we analyze the error by the difference in the $L^2$ norm between the true transition probabilities and the approximation achieved through a Monte Carlo method. We give a formula for the error for Markov chains with locally computed sampling points. Further, in the case of reversible Markov chains, we will deduce a formula for the error when sampling points are computed globally. We will see that in both cases the error itself can be approximated with Monte Carlo methods. As a consequence of the result, we will derive surprising properties of reversible Markov chains. Y1 - 2016 U6 - https://doi.org/10.1016/j.spl.2016.06.011 VL - 118 SP - 163 EP - 170 PB - Elsevier ER - TY - JOUR A1 - Lie, Han Cheng A1 - Sullivan, T. J. T1 - Cameron--Martin theorems for sequences of Cauchy-distributed random variables JF - arXiv N2 - Given a sequence of Cauchy-distributed random variables defined by a sequence of location parameters and a sequence of scale parameters, we consider another sequence of random variables that is obtained by perturbing the location or scale parameter sequences. Using a result of Kakutani on equivalence of infinite product measures, we provide sufficient conditions for the equivalence of laws of the two sequences. Y1 - 2016 SP - 1608.03784 ER - TY - JOUR A1 - Zhang, Wei A1 - Hartmann, Carsten A1 - Schütte, Christof T1 - Effective dynamics along given reaction coordinates, and reaction rate theory JF - Faraday Discussions Y1 - 2016 U6 - https://doi.org/10.1039/C6FD00147E IS - 195 SP - 365 EP - 394 ER - TY - THES A1 - Durmaz, Vedat T1 - Atomistic Binding Free Energy Estimations for Biological Host–Guest Systems N2 - Accurate quantifications of protein–ligand binding affinities by means of in silico methods increasingly gain importance in various scientific branches including toxicology and pharmacology. In silico techniques not only are generally less demanding than laboratory experiments regarding time as well as cost, in particular, if binding assays or synthesis protocols need to be developed in advance. At times, they also provide the only access to risk assessments on novel chemical compounds arising from biotic or abiotic degradation of anthropogenic substances. However, despite the continuous technological and algorithmic progress over the past decades, binding free energy estimations through molecular dynamics simulations still pose an enormous computational challenge owed to the mathematical complexity of solvated macromolecular systems often consisting of hundreds of thousands of atoms. The goals of this thesis can roughly be divided into two categories dealing with different aspects of host–guest binding quantification. On the one side algorithmic strategies for a comprehensive exploration and decomposition of conformational space in conjunction with an automated selection of representative molecular geometries and binding poses have been elaborated providing initial structures for free energy calculations. In light of the dreaded trapping problem typically associated with molecular dynamics simulations, the focus was laid on a particularly systematic generation of representatives covering a broad range of physically accessible molecular conformations and interaction modes. On the other side and ensuing from these input geometries, binding affinity models based on the linear interaction energy (LIE) method have been developed for a couple of (bio)molecular systems. The applications included a successful prediction of the liquid-chromatographic elution order as well as retention times of highly similar hexabromocyclododecane (HBCD) stereoisomers, a novel empirical LIE–QSAR hybrid binding affinity model related to the human estrogen receptor α (ERα), and, finally, the (eco)toxicological prioritization of transformation products originating from the antibiotic sulfamethoxazole with respect to their binding affinities to the bacterial enzyme dihydropteroate synthase. Altogether, a fully automated approach to binding mode and affinity estimation has been presented that is content with an arbitrary geometry of a small molecule under observation and a spatial vector specifying the binding site of a potential target molecule. According to our studies, it is superior to conventional docking and thermodynamic average methods and primarily suggesting binding free energy calculation on the basis of several heavily distinct complex geometries. Both chromatographic retention times of HBCD and binding affinities to ERα yielded squared coefficients of correlation with experimental results significantly higher than 0.8. Approximately 85 % (100 %) of predicted receptor–ligand binding modes deviated less than 1.53 Å (2.05 Å) from available crystallographic structures. KW - molecular modelling KW - molecular dynamics KW - molecular simulation KW - binding affinity KW - prediction KW - free energy KW - receptor ligand KW - retention time KW - liquid chromatography Y1 - 2016 UR - http://www.diss.fu-berlin.de/diss/receive/FUDISS_thesis_000000103765 PB - FU Dissertationen Online ER - TY - CHAP A1 - Ehrig, Rainald A1 - Dierkes, Thomas A1 - Schäfer, Stefan A1 - Röblitz, Susanna A1 - Tronci, Enrico A1 - Mancini, Toni A1 - Salvo, Ivano A1 - Alimguzhin, Vadim A1 - Mari, Federico A1 - Melatti, Igor A1 - Massini, Annalisa A1 - Krüger, Tillmann H. C. A1 - Egli, Marcel A1 - Ille, Fabian A1 - Leeners, Brigitte T1 - An integrative approach for model driven computation of treatments in reproductive medicine T2 - BIOMAT - Proceedings of the 15th International Symposium on Mathematical and Computational Biology, Rorkee, India Y1 - 2016 ER - TY - GEN A1 - Weiser, Martin A1 - Ghosh, Sunayana T1 - Theoretically optimal inexact SDC methods N2 - In several inital value problems with particularly expensive right hand side computation, there is a trade-off between accuracy and computational effort in evaluating the right hand sides. We consider inexact spectral deferred correction (SDC) methods for solving such non-stiff initial value problems. SDC methods are interpreted as fixed point iterations and, due to their corrective iterative nature, allow to exploit the accuracy-work-tradeoff for a reduction of the total computational effort. On one hand we derive an error model bounding the total error in terms of the right hand side evaluation errors. On the other hand, we define work models describing the computational effort in terms of the evaluation accuracy. Combining both, a theoretically optimal tolerance selection is worked out by minimizing the total work subject to achieving the requested tolerance. T3 - ZIB-Report - 16-52 KW - Spectral deferred correction, initial value problems, error propagation, adaptive control of tolerances Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-53140 SN - 1438-0064 ER - TY - THES A1 - Burtscher, Felicia T1 - Algorithmic approaches to the change-point problem for ovulation detection in the human menstrual cycle Y1 - 2016 ER - TY - JOUR A1 - Klus, Stefan A1 - Schütte, Christof T1 - Towards tensor-based methods for the numerical approximation of the Perron-Frobenius and Koopman operator JF - Journal of Computational Dynamics N2 - The global behavior of dynamical systems can be studied by analyzing the eigenvalues and corresponding eigenfunctions of linear operators associated with the system. Two important operators which are frequently used to gain insight into the system's behavior are the Perron-Frobenius operator and the Koopman operator. Due to the curse of dimensionality, computing the eigenfunctions of high-dimensional systems is in general infeasible. We will propose a tensor-based reformulation of two numerical methods for computing finite-dimensional approximations of the aforementioned infinite-dimensional operators, namely Ulam's method and Extended Dynamic Mode Decomposition (EDMD). The aim of the tensor formulation is to approximate the eigenfunctions by low-rank tensors, potentially resulting in a significant reduction of the time and memory required to solve the resulting eigenvalue problems, provided that such a low-rank tensor decomposition exists. Typically, not all variables of a high-dimensional dynamical system contribute equally to the system's behavior, often the dynamics can be decomposed into slow and fast processes, which is also reflected in the eigenfunctions. Thus, the weak coupling between different variables might be approximated by low-rank tensor cores. We will illustrate the efficiency of the tensor-based formulation of Ulam's method and EDMD using simple stochastic differential equations. Y1 - 2016 U6 - https://doi.org/10.3934/jcd.2016007 VL - 3 IS - 2 SP - 139 EP - 161 ER - TY - JOUR A1 - Hartmann, Carsten A1 - Schütte, Christof A1 - Zhang, Wei T1 - Model reduction algorithms for optimal control and importance sampling of diffusions JF - Nonlinearity N2 - We propose numerical algorithms for solving optimal control and importance sampling problems based on simplified models. The algorithms combine model reduction techniques for multiscale diffusions and stochastic optimization tools, with the aim of reducing the original, possibly high-dimensional problem to a lower dimensional representation of the dynamics, in which only a few relevant degrees of freedom are controlled or biased. Specifically, we study situations in which either a reaction coordinate onto which the dynamics can be projected is known, or situations in which the dynamics shows strongly localized behavior in the small noise regime. No explicit assumptions about small parameters or scale separation have to be made. We illustrate the approach with simple, but paradigmatic numerical examples. Y1 - 2016 U6 - https://doi.org/10.1088/0951-7715/29/8/2298 VL - 29 IS - 8 SP - 2298 EP - 2326 ER - TY - JOUR A1 - Klus, Stefan A1 - Koltai, Peter A1 - Schütte, Christof T1 - On the numerical approximation of the Perron-Frobenius and Koopman operator JF - Journal of Computational Dynamics N2 - Information about the behavior of dynamical systems can often be obtained by analyzing the eigenvalues and corresponding eigenfunctions of linear operators associated with a dynamical system. Examples of such operators are the Perron-Frobenius and the Koopman operator. In this paper, we will review di� fferent methods that have been developed over the last decades to compute � infinite-dimensional approximations of these in� finite-dimensional operators - in particular Ulam's method and Extended Dynamic Mode Decomposition (EDMD) - and highlight the similarities and di� fferences between these approaches. The results will be illustrated using simple stochastic di� fferential equations and molecular dynamics examples. Y1 - 2016 U6 - https://doi.org/10.3934/jcd.2016003 VL - 3 IS - 1 SP - 51 EP - 77 ER - TY - CHAP A1 - Götschel, Sebastian A1 - Höhne, Christian A1 - Kolkoori, Sanjeevareddy A1 - Mitzscherling, Steffen A1 - Prager, Jens A1 - Weiser, Martin T1 - Ray Tracing Boundary Value Problems: Simulation and SAFT Reconstruction for Ultrasonic Testing T2 - Proceedings 19th World Conference on Non-Destructive Testing (WCNDT 2016) Y1 - 2016 UR - https://www.wcndt2016.com/portals/wcndt/bb/Fr1H4.pdf ER - TY - CHAP A1 - Götschel, Sebastian A1 - Maierhofer, Christiane A1 - Müller, Jan A1 - Rothbart, Nick A1 - Weiser, Martin T1 - Quantitative Defect Reconstruction in Active Thermography for Fiber-Reinforced Composites T2 - Proceedings 19th World Conference on Non-Destructive Testing (WCNDT 2016) Y1 - 2016 UR - https://www.wcndt2016.com/portals/wcndt/bb/Th4C4.pdf ER - TY - JOUR A1 - Koltai, Peter A1 - Ciccotti, Giovanni A1 - Schütte, Christof T1 - On Markov state models for non-equilibrium molecular dynamics JF - The Journal of Chemical Physics Y1 - 2016 U6 - https://doi.org/10.1063/1.4966157 N1 - 2016 Editor's Choice of The Journal of Chemical Physics VL - 145 IS - 174103 ER - TY - THES A1 - Schuh, Hannah T1 - Identification of Anisotropies in Welding Seams Using Ultrasound Measurements Y1 - 2016 ER - TY - THES A1 - Manley, Phillip T1 - Simulation of Plasmonic Nanoparticles in Thin Film Solar Cells Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:188-fudissthesis000000102531-2 UR - http://www.diss.fu-berlin.de/diss/receive/FUDISS_thesis_000000102531 ER - TY - JOUR A1 - Eisenhauer, David A1 - Jäger, Klaus A1 - Köppel, Grit A1 - Rech, Bernd A1 - Becker, Christiane T1 - Optical Properties of Smooth Anti-reflective Three-dimensional Textures for Silicon Thin-film Solar Cells JF - Energy Procedia Y1 - 2016 U6 - https://doi.org/10.1016/j.egypro.2016.11.314 VL - 102 SP - 27 EP - 35 ER - TY - JOUR A1 - Gelß, Patrick A1 - Matera, Sebastian A1 - Schütte, Christof T1 - Solving the master equation without kinetic Monte Carlo: Tensor train approximations for a CO oxidation model JF - Journal of Computational Physics N2 - In multiscale modeling of heterogeneous catalytic processes, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. Usually, this is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO2(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased Y1 - 2016 U6 - https://doi.org/10.1016/j.jcp.2016.03.025 VL - 314 SP - 489 EP - 502 ER - TY - JOUR A1 - Vega, Iliusi A1 - Schütte, Christof A1 - Conrad, Tim T1 - Finding metastable states in real-world time series with recurrence networks JF - Physica A: Statistical Mechanics and its Applications N2 - In the framework of time series analysis with recurrence networks, we introduce a self-adaptive method that determines the elusive recurrence threshold and identifies metastable states in complex real-world time series. As initial step, we introduce a way to set the embedding parameters used to reconstruct the state space from the time series. We set them as the ones giving the maximum Shannon entropy of the diagonal line length distribution for the first simultaneous minima of recurrence rate and Shannon entropy. To identify metastable states, as well as the transitions between them, we use a soft partitioning algorithm for module finding which is specifically developed for the case in which a system shows metastability. We illustrate our method with a complex time series example. Finally, we show the robustness of our method for identifying metastable states. Our results suggest that our method is robust for identifying metastable states in complex time series, even when introducing considerable levels of noise and missing data points. Y1 - 2016 U6 - https://doi.org/10.1016/j.physa.2015.10.041 VL - 445 SP - 1 EP - 17 ER -