4.2 Basis Functions and Matrix Representation

For the purpose of discretization one can expand the HAMILTONian, self-energies, and GREEN's functions in terms of some basis functions to obtain the corresponding matrices. In the tight-binding method one can take the basis functions to be any set of localized functions, such as atomic s- and p-orbitals [116,122], WANNIER functions [242], and so forth. A common approximation used to describe the HAMILTONian of layered structures consists of non-vanishing interactions only between nearest neighbor layers. That is, each layer $ i$ interacts only with itself and its nearest neighbor layers $ i-1$ and $ i+1$. Then the single particle HAMILTONian of the layered structure is a block tri-diagonal matrix, where diagonal blocks $ \ensuremath{{\underline{H}}}_i$ represent the Hamiltonian of layer $ i$ and off-diagonal blocks $ \ensuremath{{\underline{t}}}_{i,i+1}$ represent interaction between layers $ i$ and $ i+1$

$\displaystyle \ensuremath{{\underline{H}}} \ = \ \left[ \begin{array}{ccccccccc...
... \bullet & \\ & & & & & & & \bullet & \bullet & \bullet \end{array} \right] \ ,$ (4.10)

where $ \ensuremath{{\underline{t}}}_{i+1,i}=\ensuremath{{\underline{t}}}^\dagger_{i,i+1}$. The matrix representation of the kinetic equations (3.71) and (3.72) are

\begin{displaymath}\begin{array}{l} [E \ensuremath{{\underline{I}}} - \ensuremat...
...}^\mathrm{r} \ = \ \ensuremath{{\underline{I}}} \ , \end{array}\end{displaymath} (4.11)


\begin{displaymath}\begin{array}{l} \ensuremath{{\underline{G}}}^{\gtrless} \ = ...
...Scat} \ \ensuremath{{\underline{G}}}^\mathrm{a} \ , \end{array}\end{displaymath} (4.12)

where $ \ensuremath{{\underline{\Sigma}}}_\mathrm{Scat}$ is the self-energy due to scattering processes and $ \ensuremath{{\underline{G}}}^\mathrm{a} ={[\ensuremath{{\underline{G}}}^\mathrm{r}]}^\dagger$ [60].
Figure 4.5: Partitioning of the simulation domain into device region and left and right contacts. Each point corresponds to a layer.
\includegraphics[width=0.48\textheight]{figures/Layered.eps}
One can partition the layered structure into left contact with index $ L$, device region with index $ D$, and right contact with index $ R$ (Fig. 4.5). The device corresponds to the region where one solves the transport equations and the contacts are the highly conducting regions connected to the device. While the device region consists of only $ N$ layers, the matrix equation corresponding to (4.11) is infinitely dimensional due to the semi-infinite contacts. As shown in Appendix G the influence of the semi-infinite contacts can be folded into the device region by adding a self-energy to the device region. This can be viewed as an additional self-energy due to the transitions between the device and the contacts. In the next sections the matrix representation of the HAMILTONian and self-energies are discussed in detail.

M. Pourfath: Numerical Study of Quantum Transport in Carbon Nanotube-Based Transistors