Nonlinear diffusion, boundary layers and nonsmoothness: Analysis of challenges in drift–diffusion semiconductor simulations

Abstract

We study different discretizations of the van Roosbroeck system for charge transport in bulk semiconductor devices that can handle nonlinear diffusion. Three common challenges corrupting the precision of numerical solutions will be discussed: boundary layers, discontinuities in the doping profile, and corner singularities in L-shaped domains. We analyze and benchmark the error and the convergence order of finite difference, finite-element as well as advanced Scharfetter–Gummel type finite-volume discretization schemes. The most problematic of these challenges are boundary layers in the quasi-Fermi potentials near ohmic contacts, which can have a drastic impact on the convergence order. Using a novel formal asymptotic expansion, our theoretical analysis reveals that these boundary layers are logarithmic and significantly shorter than the Debye length.