@misc{GutscheLaeuterSchmidt, author = {Gutsche, Philipp and L{\"a}uter, Matthias and Schmidt, Frank}, title = {Parameter-dependent Parallel Block Sparse Arnoldi and D{\"o}hler Algorithms on Distributed Systems}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-58202}, abstract = {We summarize the basics and first results of the analyses within our ZIB Bridge Project and give an outlook on further studies broadening the usage of hardware acceleration within the Finite Element Method (FEM) based solution of Maxwell's equations.}, language = {en} } @misc{HammerschmidtPomplunBurgeretal., author = {Hammerschmidt, Martin and Pomplun, Jan and Burger, Sven and Schmidt, Frank}, title = {Adaptive sampling strategies for efficient parameter scans in nano-photonic device simulations}, issn = {1438-0064}, doi = {10.1117/12.2036363}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-50395}, abstract = {Rigorous optical simulations are an important tool in optimizing scattering properties of nano-photonic devices and are used, for example, in solar cell optimization. The finite element method (FEM) yields rigorous, time-harmonic, high accuracy solutions of the full 3D vectorial Maxwell's equations [1] and furthermore allows for great flexibility and accuracy in the geometrical modeling of these often complex shaped 3D nano-structures. A major drawback of frequency domain methods is the limitation of single frequency evaluations. For example the accurate computation of the short circuit current density of an amorphous silicon / micro-crystalline multi-junction thin film solar cell may require the solution of Maxwell's equations for over a hundred different wavelengths if an equidistant sampling strategy is employed. Also in optical metrology, wavelength scans are frequently used to reconstruct unknown geometrical and material properties of optical systems numerically from measured scatterometric data. In our contribution we present several adaptive numerical integration and sampling routines and study their efficiency in the context of the determination of generation rate profiles of solar cells. We show that these strategies lead to a reduction in the computational effort without loss of accuracy. We discuss the employment of tangential information in a Hermite interpolation scheme to achieve similar accuracy on coarser grids. We explore the usability of these strategies for scatterometry and solar cell simulations.}, language = {en} } @misc{HammerschmidtLockauZschiedrichetal., author = {Hammerschmidt, Martin and Lockau, Daniel and Zschiedrich, Lin and Schmidt, Frank}, title = {Optical modelling of incoherent substrate light-trapping in silicon thin film multi-junction solar cells with finite elements and domain decomposition}, issn = {1438-0064}, doi = {10.1117/12.2036346}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-50410}, abstract = {In many experimentally realized applications, e.g. photonic crystals, solar cells and light-emitting diodes, nano-photonic systems are coupled to a thick substrate layer, which in certain cases has to be included as a part of the optical system. The finite element method (FEM) yields rigorous, high accuracy solutions of full 3D vectorial Maxwell's equations [1] and allows for great flexibility and accuracy in the geometrical modelling. Time-harmonic FEM solvers have been combined with Fourier methods in domain decomposition algorithms to compute coherent solutions of these coupled system. [2,3] The basic idea of a domain decomposition approach lies in a decomposition of the domain into smaller subdomains, separate calculations of the solutions and coupling of these solutions on adjacent subdomains. In experiments light sources are often not perfectly monochromatic and hence a comparision to simulation results might only be justified if the simulation results, which include interference patterns in the substrate, are spectrally averaged. In this contribution we present a scattering matrix domain decomposition algorithm for Maxwell's equations based on FEM. We study its convergence and advantages in the context of optical simulations of silicon thin film multi-junction solar cells. This allows for substrate light-trapping to be included in optical simulations and leads to a more realistic estimation of light path enhancement factors in thin-film devices near the band edge.}, language = {en} } @misc{HammerschmidtHerrmannPomplunetal., author = {Hammerschmidt, Martin and Herrmann, Sven and Pomplun, Jan and Zschiedrich, Lin and Burger, Sven and Schmidt, Frank}, title = {Reduced basis method for Maxwell's equations with resonance phenomena}, series = {Proc. SPIE}, volume = {9630}, journal = {Proc. SPIE}, issn = {1438-0064}, doi = {10.1117/12.2190425}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-55687}, pages = {96300R}, abstract = {Rigorous optical simulations of 3-dimensional nano-photonic structures are an important tool in the analysis and optimization of scattering properties of nano-photonic devices or parameter reconstruction. To construct geometrically accurate models of complex structured nano-photonic devices the finite element method (FEM) is ideally suited due to its flexibility in the geometrical modeling and superior convergence properties. Reduced order models such as the reduced basis method (RBM) allow to construct self-adaptive, error-controlled, very low dimensional approximations for input-output relationships which can be evaluated orders of magnitude faster than the full model. This is advantageous in applications requiring the solution of Maxwell's equations for multiple parameters or a single parameter but in real time. We present a reduced basis method for 3D Maxwell's equations based on the finite element method which allows variations of geometric as well as material and frequency parameters. We demonstrate accuracy and efficiency of the method for a light scattering problem exhibiting a resonance in the electric field.}, language = {en} } @misc{Gutsche2014, type = {Master Thesis}, author = {Gutsche, Philipp}, title = {Convergence Study of the Fourier Modal Method for Nano-optical Scattering Problems in Comparison with the Finite Element Method}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-56084}, pages = {87}, year = {2014}, abstract = {Nano-optical scattering problems play an important role in our modern, technologically driven society. Computers, smartphones and all kinds of electronic devices are manufactured by the semiconductor industry which relies on production using photomasks as well as optical process control. The digital world, e.g. the world wide web, is based on optical interconnects and so-called quantum computers based on optics are supposed to be next generation computers. Moreover, global economic progress demands new and sustainable energy resources and one option is to make use of the power stored in optical radiation from the sun. Additionally, understanding fundamental physics such as the optical properties of asymmetric, or chiral, structures could promote future innovations in engineering. In order to understand and manipulate these kinds of processes, physics provides a well established model: the so-called Maxwell's equations. Stated by James Clerk Maxwell in 1862, this description of the interaction of light and matter still provides a profound basis for the analysis of electromagnetic phenomena. However, real world problems cannot be calculated using simple mathematics. Rather, computer simulations are needed to obtain solutions of the physical model. Finding suitable methods to solve these problems opens up a wide variety of possibilities. On the one hand, there are methods which require long computing times. On the other hand, some algorithms depend on high memory usage. That is why the field of numerics deals with the question which method is optimally suited for specific problems. The aim of this work is to investigate the applicability of the so-called Fourier Modal Method (FMM) to nano-optical scattering problems in general. Since simple analytical solutions are non-existent for most recent physical problems, we use the Finite Element Method (FEM) to double-check performance of the FMM. Mathematics provide reliable procedures to control the errors of numerics using the FEM. Yet up to now it has not been possible to rigorously classify the quality of the Fourier Modal Method's results. It is not fully understood whether the process of investing more and more computing resources yields more accurate results. So, we have to ask ourselves: does the numerical method invariably converge? In spite of this uncertainty when using the FMM, it is a well established method dating back to the 1980s. This numerical method has recently been used to optimize performance of solar cells [19] as well as to improve the optical properties of so-called single-photon sources [41] which are essential for quantum cryptography. The latter is a promising candidate to increase digital security and revolutionise cryptography techniques. Furthermore, with the help of the Fourier Modal Method an important issue in optics has been partly resolved: angular filtering of light was made possible by using a mirror which becomes transparent at a certain viewing angle [77]. In addition, an improved numerical technique to design so-called Photonic Crystal waveguides based on the FMM was developed recently [15]. Photonic Crystals are used in the fields of optical bio-sensing and for the construction of novel semiconductor devices. Moreover, approaches to link the FMM and the FEM try to combine advantages of both methods to obtain fast and accurate results [81]. These ideas are closely linked to the well-known concept of Domain Decomposition within the FEM [88]. Here, one possibility to couple domains is to use the scattering matrix formalism as it is done in the FMM. In the scope of this convergence study, we state Maxwell's equations, particularly for periodic geometries. We describe two physical phenomena of nano-optics, namely chirality and opto-electrical coupling, and define the errors of our simulations. Afterwards, the two investigated methods are analysed with respect to their general properties and a way to unify modelling physics when using both algorithms is presented. With the help of various numerical experiments, we explore convergence characteristics of the FMM and draw conclusions about the ability of this approach to provide accurate results and, consequently, its potential for research on technological innovations.}, language = {en} } @misc{HammerschmidtHerrmannPomplunetal., author = {Hammerschmidt, Martin and Herrmann, Sven and Pomplun, Jan and Burger, Sven and Schmidt, Frank}, title = {Model order reduction for the time-harmonic Maxwell equation applied to complex nanostructures}, series = {Proc. SPIE}, volume = {9742}, journal = {Proc. SPIE}, issn = {1438-0064}, doi = {10.1117/12.2212367}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57239}, pages = {97420M}, abstract = {Fields such as optical metrology and computational lithography require fast and efficient methods for solving the time-harmonic Maxwell's equation. Highly accurate geometrical modeling and numerical accuracy atcomputational costs are a prerequisite for any simulation study of complex nano-structured photonic devices. We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem based on the hp-adaptive finite element solver JCMsuite capable of handling geometric and non-geometric parameter dependencies allowing for online evaluations in milliseconds. We apply the RBM to compute light-scatteringoptical wavelengths off periodic arrays of fin field-effect transistors (FinFETs) where geometrical properties such as the width and height of the fin and gate can vary in a large range.}, language = {en} } @misc{HammerschmidtBarthPomplunetal., author = {Hammerschmidt, Martin and Barth, Carlo and Pomplun, Jan and Burger, Sven and Becker, Christiane and Schmidt, Frank}, title = {Reconstruction of photonic crystal geometries using a reduced basis method for nonlinear outputs}, issn = {1438-0064}, doi = {10.1117/12.2212482}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57249}, abstract = {Maxwell solvers based on the hp-adaptive finite element method allow for accurate geometrical modeling and high numerical accuracy. These features are indispensable for the optimization of optical properties or reconstruction of parameters through inverse processes. High computational complexity prohibits the evaluation of the solution for many parameters. We present a reduced basis method (RBM) for the time-harmonic electromagnetic scattering problem allowing to compute solutions for a parameter configuration orders of magnitude faster. The RBM allows to evaluate linear and nonlinear outputs of interest like Fourier transform or the enhancement of the electromagnetic field in milliseconds. We apply the RBM to compute light-scattering off two dimensional photonic crystal structures made of silicon and reconstruct geometrical parameters.}, language = {en} } @misc{HammerschmidtHerrmannBurgeretal., author = {Hammerschmidt, Martin and Herrmann, Sven and Burger, Sven and Pomplun, Jan and Schmidt, Frank}, title = {Reduced basis method for the optimization of nano-photonic devices}, issn = {1438-0064}, doi = {10.1007/s11082-016-0530-1}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-57556}, abstract = {Optical 3D simulations in many-query and real-time contexts require new solution strategies. We study an adaptive, error controlled reduced basis method for solving parametrized time-harmonic optical scattering problems. Application fields are, among others, design and optimization problems of nano-optical devices as well as inverse problems for parameter reconstructions occuring e. g. in optical metrology. The reduced basis method presented here relies on a finite element modeling of the scattering problem with parametrization of materials, geometries and sources.}, language = {en} } @misc{HammerschmidtWeiserSantiagoetal., author = {Hammerschmidt, Martin and Weiser, Martin and Santiago, Xavier Garcia and Zschiedrich, Lin and Bodermann, Bernd and Burger, Sven}, title = {Quantifying parameter uncertainties in optical scatterometry using Bayesian inversion}, issn = {1438-0064}, url = {http://nbn-resolving.de/urn:nbn:de:0297-zib-64704}, abstract = {We present a Newton-like method to solve inverse problems and to quantify parameter uncertainties. We apply the method to parameter reconstruction in optical scatterometry, where we take into account a priori information and measurement uncertainties using a Bayesian approach. Further, we discuss the influence of numerical accuracy on the reconstruction result.}, language = {en} } @article{Le2019, author = {Le, Thi Thai}, title = {Kelvin-Helmholtz instability in a shallow-water flow with a finite width}, series = {Journal of Mathematical Physics}, volume = {60}, journal = {Journal of Mathematical Physics}, number = {1}, publisher = {AIP Publishing}, doi = {10.1063/1.5126321}, pages = {123101}, year = {2019}, abstract = {We examine an effect of side walls on the linear stability of an interface of tangential-velocity discontinuity in shallow-water flow. The flow is pure horizontal in the plane xy, and the fluid is bounded in a finite width 2d in the y- direction. In region 0 < y < d, the fluid is moving with uniform velocity U but is at rest for -d < y < 0. Without side walls, the flow is unstable for a velocity difference U<√8c U < √8 c, with c being the velocity of gravity waves. In this work, we show that if the velocity difference U is smaller than 2c, the interface is always destabilized, also known as the flow is unstable. The unstable region of an infinite width model is shrunken by the effects of side walls in the case of narrow width, while there is no range for the Froude number for stabilization in the case of large width. These results play an important role in predicting the wave propagations and have a wide application in the fields of industry. As a result of the interaction of waves and the mean flow boundary, the flow is unstable, which is caused by a decrease in the kinetic energy of disturbance.}, language = {en} } @article{JinLeFukumoto, author = {Jin, Liangbing and Le, Thi Thai and Fukumoto, Yasuhide}, title = {Frictional effect on stability of discontinuity interface in tangential velocity of a shallow-water flow}, series = {Physics Letters A}, volume = {383}, journal = {Physics Letters A}, number = {26}, doi = {10.1016/j.physleta.2019.125839}, pages = {125839}, abstract = {We examine a frictional effect on the linear stability of an interface of discontinuity in tangential velocity. The fluid is moving with uniform velocity U in a region but is at rest in the other, and the bottom surface is assumed to exert drag force, quadratic in velocity, on the thin fluid layer. In the absence of the drag, the instability of the Kelvin-Helmholtz type is suppressed for U>√8 c, with c being the propagating speed of the gravity wave. We find by asymptotic analyses for both small and large values of the drag strength that the drag, regardless of its strength, makes the flow unstable for the whole range of the Froude number U/c.}, language = {en} }