The transient Wigner equation with relaxation time scattering is
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
Here denotes the relaxation time, denotes carrier density
(8.4) |
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) |
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].
Previous: 8. Finite Difference Wigner Up: 8. Finite Difference Wigner Next: 8.2 Discrete Wigner Transform