Search Results

Now showing 1 - 5 of 5
  • Item
    Existence of bounded discrete steady state solutions of the van Roosbroeck system with monotone Fermi-Dirac statistic functions
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2014) Gärtner, Klaus
    If the statistic function is modified, the equations can be derived by a variational formulation or just using a generalized Einstein relation. In both cases a dissipative generalization of the Scharfetter-Gummel scheme citeSch_Gu, understood as a one-dimensional constant current approximation, is derived for strictly monotone coefficient functions in the elliptic operator $nabla cdot bal ff(v) nabla $, $v$ chemical potential, while the hole density is defined by $p=cal F(v)le e^v.$ A closed form integration of the governing equation would simplify the practical use, but mean value theorem based results are sufficient to prove existence of bounded discrete steady state solutions on any boundary conforming Delaunay grid. These results hold for any piecewise, continuous, and monotone approximation of $bal ff(v)$ and $cal F(v)$.
  • Item
    Generalized Scharfetter--Gummel schemes for electro-thermal transport in degenerate semiconductors using the Kelvin formula for the Seebeck coefficient
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2019) Kantner, Markus
    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.
  • Item
    Computational and analytical comparison of flux discretizations for the semiconductor device equations beyond Boltzmann statistics
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Farrell, Patricio; Koprucki, Thomas; Fuhrmann, Jürgen
    For a Voronoi finite volume discretization of the van Roosbroeck system with general charge carrier statistics we compare three thermodynamically consistent numerical fluxes known in the literature. We discuss an extension of the Scharfetter-Gummel scheme to non-Boltzmann (e.g. Fermi-Dirac) statistics. It is based on the analytical solution of a two-point boundary value problem obtained by projecting the continuous differential equation onto the interval between neighboring collocation points. Hence, it serves as a reference flux. The exact solution of the boundary value problem can be approximated by computationally cheaper fluxes which modify certain physical quantities. One alternative scheme averages the nonlinear diffusion (caused by the non-Boltzmann nature of the problem), another one modifies the effective density of states. To study the differences between these three schemes, we analyze the Taylor expansions, derive an error estimate, visualize the flux error and show how the schemes perform for a carefully designed p-i-n benchmark simulation. We present strong evidence that the flux discretization based on averaging the nonlinear diffusion has an edge over the scheme based on modifying the effective density of states.
  • Item
    Assessing the quality of the excess chemical potential flux scheme for degenerate semiconductor device simulation
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2020) Abdel, Dilara; Farrell, Patricio; Fuhrmann, Jürgen
    The van Roosbroeck system models current flows in (non-)degenerate semiconductor devices. Focusing on the stationary model, we compare the excess chemical potential discretization scheme, a flux approximation which is based on a modification of the drift term in the current densities, with another state-of-the-art Scharfetter-Gummel scheme, namely the diffusion-enhanced scheme. Physically, the diffusion-enhanced scheme can be interpreted as a flux approximation which modifies the thermal voltage. As a reference solution we consider an implicitly defined integral flux, using Blakemore statistics. The integral flux refers to the exact solution of a local two point boundary value problem for the continuous current density and can be interpreted as a generalized Scharfetter-Gummel scheme. All numerical discretization schemes can be used within a Voronoi finite volume method to simulate charge transport in (non-)degenerate semiconductor devices. The investigation includes the analysis of Taylor expansions, a derivation of error estimates and a visualization of errors in local flux approximations to extend previous discussions. Additionally, drift-diffusion simulations of a p-i-n device are performed.
  • Item
    Multi-dimensional modeling and simulation of semiconductor nanophotonic devices
    (Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2019) Kantner, Markus; Höhne, Theresa; Koprucki, Thomas; Burger, Sven; Wünsche, Hans-Jürgen; Schmidt, Frank; Mielke, Alexander; Bandelow, Uwe
    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.