    Previous: 8. Finite Difference Wigner Up: 8. Finite Difference Wigner Next: 8.2 Discrete Wigner Transform

# 8.1 Wigner Equation

The transient Wigner equation with relaxation time scattering is (8.1)

Here is the 1-D Wigner function at position , wavenumber and position . The non-local Wigner potential is calculated from the real potential energy by (8.2)

In Equation 8.1 denotes the term for relaxation time scattering, that is (8.3)

Here denotes the relaxation time, denotes carrier density (8.4)

and denotes the equilibrium solution. Assuming that the doping is the same in the left and the right electrode we define as the solution of Equation 8.1 under zero bias (and relaxation time scattering excluded).

Note the factor which comes from the Wigner transformation of the density matrix trace operation.

The Wigner equation has to be supplemented with suitable boundary conditions. It is of first order in , hence a possibility is to specify the solution on one side of the simulation domain. We used inflow boundary conditions: these are of Dirichlet type in Wigner phase space and specify the incoming flow on the left and on the right contacts. Hence half of the boundary conditions is given on the left, half is given on the right boundary. We assume that even under bias, the distribution in the leads is in equilibrium and the incoming part of the distribution at the contacts is given by this equilibrium distribution. The equilibrium distribution is given by a Fermi-Dirac distribution which poses the boundary conditions: (8.5)

The Fermi level is correct when the carrier concentration equals the doping density.

We allow the effective mass to be a space dependent parameter. This is only an approximation to the Wigner transformation of the term in the Schrödinger equation and is known as Frensley's model'' or as classical approximation''. The correct term as obtained by Wigner transformation leads to a convolution integral and is treated in [TOM91], [MH94].