Comparison of thermodynamically consistent charge carrier flux discretizations for Fermi-Dirac and Gauss-Fermi statistics

Abstract

We compare three thermodynamically consistent ScharfetterGummel schemes for different distribution functions for the carrier densities, including the FermiDirac integral of order 1/2 and the GaussFermi integral. The most accurate (but unfortunately also most costly) generalized ScharfetterGummel scheme requires the solution of an integral equation. We propose a new method to solve this integral equation numerically based on Gauss quadrature and Newtons method. We discuss the quality of this approximation and plot the resulting currents for FermiDirac and GaussFermi statistics. Finally, by comparing two modified (diffusion-enhanced and inverse activity based) ScharfetterGummel schemes with the more accurate generalized scheme, we show that the diffusion-enhanced ansatz leads to considerably lower flux errors, confirming previous results (J. Comp. Phys. 346:497-513, 2017).