Search Results
Beyond just “flattening the curve”: Optimal control of epidemics with purely non-pharmaceutical interventions
When effective medical treatment and vaccination are not available, non-pharmaceutical interventions such as social distancing, home quarantine and far-reaching shutdown of public life are the only available strategies to prevent the spread of epidemics. Based on an extended SEIR (susceptible-exposed-infectious-recovered) model and continuous-time optimal control theory, we compute the optimal non-pharmaceutical intervention strategy for the case that a vaccine is never found and complete containment (eradication of the epidemic) is impossible. In this case, the optimal control must meet competing requirements: First, the minimization of disease-related deaths, and, second, the establishment of a sufficient degree of natural immunity at the end of the measures, in order to exclude a second wave. Moreover, the socio-economic costs of the intervention shall be kept at a minimum. The numerically computed optimal control strategy is a single-intervention scenario that goes beyond heuristically motivated interventions and simple “flattening of the curve”. Careful analysis of the computed control strategy reveals, however, that the obtained solution is in fact a tightrope walk close to the stability boundary of the system, where socio-economic costs and the risk of a new outbreak must be constantly balanced against one another. The model system is calibrated to reproduce the initial exponential growth phase of the COVID-19 pandemic in Germany. © 2020, The Author(s).
Efficient Current Injection Into Single Quantum Dots Through Oxide-Confined p-n-Diodes
Current injection into single quantum dots embedded in vertical p-n-diodes featuring oxide apertures is analyzed in the low-injection regime suitable for single-photon emitters. The experimental and theoretical evidence is found for a rapid lateral spreading of the carriers after passing the oxide aperture in the conventional p-i-n-design. By an alternative design employing p-doping up to the oxide aperture, the current spreading can be suppressed resulting in an enhanced current confinement and increased injection efficiencies, both, in the continuous wave and under pulsed excitation.
Numerical simulation of carrier transport in semiconductor devices at cryogenic temperatures
At cryogenic temperatures the electron-hole plasma in semiconductor materials becomes strongly degenerate, leading to very sharp internal layers, extreme depletion in intrinsic domains and strong nonlinear diffusion. As a result, the numerical simulation of the drift-diffusion system suffers from serious convergence issues using standard methods. We consider a one-dimensional p-i-n diode to illustrate these problems and present a simple temperature-embedding scheme to enable the numerical simulation at cryogenic temperatures. The method is suitable for forward-biased devices as they appear e.g. in optoelectronic applications.
Multi-dimensional modeling and simulation of semiconductor nanophotonic devices
Self-consistent modeling and multi-dimensional simulation of semiconductor nanophotonic devices is an important tool in the development of future integrated light sources and quantum devices. Simulations can guide important technological decisions by revealing performance bottlenecks in new device concepts, contribute to their understanding and help to theoretically explore their optimization potential. The efficient implementation of multi-dimensional numerical simulations for computer-aided design tasks requires sophisticated numerical methods and modeling techniques. We review recent advances in device-scale modeling of quantum dot based single-photon sources and laser diodes by self-consistently coupling the optical Maxwell equations with semiclassical carrier transport models using semi-classical and fully quantum mechanical descriptions of the optically active region, respectively. For the simulation of realistic devices with complex, multi-dimensional geometries, we have developed a novel hp-adaptive finite element approach for the optical Maxwell equations, using mixed meshes adapted to the multi-scale properties of the photonic structures. For electrically driven devices, we introduced novel discretization and parameter-embedding techniques to solve the drift-diffusion system for strongly degenerate semiconductors at cryogenic temperature. Our methodical advances are demonstrated on various applications, including vertical-cavity surface-emitting lasers, grating couplers and single-photon sources.
Consistency and convergence for a family of finite volume discretizations of the Fokker–Planck operator
We introduce a family of various finite volume discretization schemes for the Fokker–Planck operator, which are characterized by different Stolarsky weight functions on the edges. This family particularly includes the well-established Scharfetter–Gummel discretization as well as the recently developed square-root approximation (SQRA) scheme. We motivate this family of discretizations both from the numerical and the modeling point of view and provide a uniform consistency and error analysis. Our main results state that the convergence order primarily depends on the quality of the mesh and in second place on the choice of the Stolarsky weights. We show that the Scharfetter–Gummel scheme has the analytically best convergence properties but also that there exists a whole branch of Stolarsky means with the same convergence quality. We show by numerical experiments that for small convection the choice of the optimal representative of the discretization family is highly non-trivial, while for large gradients the Scharfetter–Gummel scheme stands out compared to the others.
Hybrid quantum-classical modeling of quantum dot devices
The design of electrically driven quantum dot devices for quantum optical applications asks for modeling approaches combining classical device physics with quantum mechanics. We connect the well-established fields of semi-classical semiconductor transport theory and the theory of open quantum systems to meet this requirement. By coupling the van Roosbroeck system with a quantum master equation in Lindblad form, we obtain a new hybrid quantum-classical modeling approach, which enables a comprehensive description of quantum dot devices on multiple scales: It allows the calculation of quantum optical figures of merit and the spatially resolved simulation of the current flow in realistic semiconductor device geometries in a unified way. We construct the interface between both theories in such a way, that the resulting hybrid system obeys the fundamental axioms of (non-)equilibrium thermodynamics. We show that our approach guarantees the conservation of charge, consistency with the thermodynamic equilibrium and the second law of thermodynamics. The feasibility of the approach is demonstrated by numerical simulations of an electrically driven single-photon source based on a single quantum dot in the stationary and transient operation regime.
Generalized Scharfetter--Gummel schemes for electro-thermal transport in degenerate semiconductors using the Kelvin formula for the Seebeck coefficient
Many challenges faced in today's semiconductor devices are related to self-heating phenomena. The optimization of device designs can be assisted by numerical simulations using the non-isothermal drift-diffusion system, where the magnitude of the thermoelectric cross effects is controlled by the Seebeck coefficient. We show that the model equations take a remarkably simple form when assuming the so-called Kelvin formula for the Seebeck coefficient. The corresponding heat generation rate involves exactly the three classically known self-heating effects, namely Joule, recombination and Thomson--Peltier heating, without any further (transient) contributions. Moreover, the thermal driving force in the electrical current density expressions can be entirely absorbed in the (nonlinear) diffusion coefficient via a generalized Einstein relation. The efficient numerical simulation relies on an accurate and robust discretization technique for the fluxes (finite volume Scharfetter--Gummel method), which allows to cope with the typically stiff solutions of the semiconductor device equations. We derive two non-isothermal generalizations of the Scharfetter--Gummel scheme for degenerate semiconductors (Fermi--Dirac statistics) obeying the Kelvin formula. The approaches differ in the treatment of degeneration effects: The first is based on an approximation of the discrete generalized Einstein relation implying a specifically modified thermal voltage, whereas the second scheme follows the conventionally used approach employing a modified electric field. We present a detailed analysis and comparison of both schemes, indicating a superior performance of the modified thermal voltage scheme.
Semiconductor laser linewidth theory revisited
More and more applications require semiconductor lasers distinguished not only by large modulation bandwidths or high output powers, but also by small spectral linewidths. The theoretical understanding of the root causes limiting the linewidth is therefore of great practical relevance. In this paper, we derive a general expression for the calculation of the spectral linewidth step by step in a self-contained manner. We build on the linewidth theory developed in the 1980s and 1990s but look from a modern perspective, in the sense that we choose as our starting points the time-dependent coupled-wave equations for the forward and backward propagating fields and an expansion of the fields in terms of the stationary longitudinal modes of the open cavity. As a result, we obtain rather general expressions for the longitudinal excess factor of spontaneous emission (K-factor) and the effective α-factor including the effects of nonlinear gain (gain compression) and refractive index (Kerr effect), gain dispersion, and longitudinal spatial hole burning in multi-section cavity structures. The effect of linewidth narrowing due to feedback from an external cavity often described by the so-called chirp reduction factor is also automatically included. We propose a new analytical formula for the dependence of the spontaneous emission on the carrier density avoiding the use of the population inversion factor. The presented theoretical framework is applied to a numerical study of a two-section distributed Bragg reflector laser.
Delay-induced patterns in a two-dimensional lattice of coupled oscillators
We show how a variety of stable spatio-temporal periodic patterns can be created in 2D-lattices of coupled oscillators with non-homogeneous coupling delays. The results are illustrated using the FitzHugh-Nagumo coupled neurons as well as coupled limit cycle (Stuart-Landau) oscillators. A "hybrid dispersion relation" is introduced, which describes the stability of the patterns in spatially extended systems with large time-delay.
Consistency and convergence for a family of finite volume discretizations of the Fokker--Planck operator
We introduce a family of various finite volume discretization schemes for the Fokker--Planck operator, which are characterized by different weight functions on the edges. This family particularly includes the well-established Scharfetter--Gummel discretization as well as the recently developed square-root approximation (SQRA) scheme. We motivate this family of discretizations both from the numerical and the modeling point of view and provide a uniform consistency and error analysis. Our main results state that the convergence order primarily depends on the quality of the mesh and in second place on the quality of the weights. We show by numerical experiments that for small gradients the choice of the optimal representative of the discretization family is highly non-trivial while for large gradients the Scharfetter--Gummel scheme stands out compared to the others.