Carrier transport

Drift-Diffusion Model

The continuity equations in the presence of generation G recombination R of electron-hole pairs read

(2.2.20)ent+(ejn(x))=e(G(x)R(x)),ept+ejp(x)=e(G(x)R(x)),

where the current is proportional to the gradient of quasi Fermi levels EF,n/p(x)

(2.2.21)jn(x)=μn(x)n(x)EF,n(x),jp(x)=μp(x)p(x)EF,p(x).

Here the charge current has the unit of (area)1(time)1. μn/p are the mobilities of each carrier. In nextnano++, μn/p are determined using the mobility model specified in the input file under currents{ }.

Hereafter we consider stationary solutions and set n˙=p˙=0. The governing equations then reduce to

(2.2.22)μn(x)n(x)EF,n(x)=(G(x)R(x)),μp(x)p(x)EF,p(x)=G(x)R(x),

which we call current equation.

We can also say that the current equation governs the relationship between the carrier densities n(x), p(x) and quasi Fermi levels EF,n/p(x).

The nextnano++ tool solves this equation and Poisson equation (and also Schrödinger equation) self-consistently.

In their solution, the corresponding calculation of the carrier densities (n(x,ϕ,EF,n),p(x,ϕ,EF,p)) and Poisson equation are firstly iterated for a given quasi-Fermi levels until the carreir densities converge. Then the resulting carrier densities are substituted into the current equation and the quasi-Fermi levels are updated. This whole cycle is iterated until the quasi-Fermi levels satisfies the convergence criteria, which can be tuned by the users from run{ current_poisson{ } } or run{ quantum_current_poisson{ } }.


Last update: 04/12/2024