next up previous contents
Next: 2.4.2 Boltzmann Transport Equation Up: 2.4 Phonon Transport Previous: 2.4 Phonon Transport   Contents

2.4.1 Landauer Formula

In semiconductors, the largest component of the heat current is due to the phonon transport. Using the Landauer theory for phonon transport, the heat flow is proportional to the phonon transmission function $ \overline{T}_\mathrm{ph}(\omega)$ , the Bose-Einstein distribution function $ n(\omega)$ of source and drain contacts, and the phonon energy/frequency $ \omega$ as:

$\displaystyle I_q=\frac{1}{h}\int_{0}^{+\infty} \overline{T}_{\mathrm{ph}}(\ome...
...left [ n_{\mathrm{s}}(\omega)-n_{\mathrm{d}}(\omega) \right ] \ d(\hbar \omega)$ (2.52)

In the linear response regime of the Landauer theory the heat current is proportional to the temperature difference by:

\begin{displaymath}\begin{split}I_q&= \frac{1}{h}\int_{0}^{+\infty} \overline{T}...
...c{\partial n(\omega)}{\partial T} \ d(\hbar \omega) \end{split}\end{displaymath} (2.53)

This can alternatively be written as:

$\displaystyle I_q=-\Delta T \frac{k_{\mathrm{B}}^2T\pi^2}{3h}\int_{0}^{+\infty}...
...erline{T}_{\mathrm{ph}}(\omega) W_{\mathrm{ph}}(\hbar \omega) \ d(\hbar \omega)$ (2.54)

where the phononic window function is given by [57]:

\begin{displaymath}\begin{split}W_{\mathrm{ph}}(\hbar \omega)&=-\frac{3}{\pi^2}\...
...ar \omega}{2 k_{\mathrm{B}}T}\right ) \right ]^{-2} \end{split}\end{displaymath} (2.55)

As $ I_q=-K_\mathrm{ph} \Delta T$ , one can express the lattice contribution to the thermal conductance as [58]:

$\displaystyle K_\mathrm{ph}=\frac{k_{\mathrm{B}}^2 T\pi^2}{3h}\int_{0}^{+\infty}\overline{T}_\mathrm{ph}(\omega) W_{\mathrm{ph}}(\omega) d(\hbar\omega)$ (2.56)

To extract the ballistic phonon transmission $ \overline {T}_{\mathrm{ph}}(E)$ in Eq. 2.56, similar to the electronic quantities (Eqs. 2.39-2.42), the phononic density of modes $ M_{\mathrm{Ph}}(E)$ is calculated using the phononic bandstructures, where $ E=\hbar \omega$ denotes the energy of a phonon with frequency $ \omega$ . In our calculations for phonon density of states $ DOS_{\mathrm{ph}}$ , we broadened the delta function by $ 1~\mathrm{meV}$ . This helps to numerically smoothen the DOS function, without affecting the thermal conductance calculations.

The FCM can also be effectively coupled to NEGF for the investigation of coherent phonon transport in low dimensional systems. The NEGF method has been traditionally employed for electronic transport studies, but has been extended to phonon studies as well [54]. Here, the phononic device Green's function is obtained by

$\displaystyle G(E)=\left (E^2I-D-\Sigma_s-\Sigma_d \right)^{-1}$ (2.57)

where $ D$ is device dynamic matrix and $ E=\hbar \omega$ is the phonon energy [59]. Note that this is a Green's function of a classical wave equation, which is second order in time. Therefore, the square of the energy appears in Eq. 2.57. This is in contrast to the Schrodinger equation for quantum treatment of electrons which is first order in time, such that the eigen energy $ E$ appears linearly in that case.

Finally, similar to the electron transport, the transmission probability of phonons through the rough channel can then be calculated using:

$\displaystyle \overline{T}_{\mathrm{ph}}(E)=\mathrm{Trace}[\Gamma_{\mathrm{s}}G\Gamma_{\mathrm{d}}G^{\dagger}]$ (2.58)

Here, the phononic broadening and contact self-energy matrices are obtained in a similar way as the electronic ones (Eq. 2.51) with the substitutions $ H\rightarrow D$ and $ E\rightarrow \hbar^2\omega^2$ .
next up previous contents
Next: 2.4.2 Boltzmann Transport Equation Up: 2.4 Phonon Transport Previous: 2.4 Phonon Transport   Contents
H. Karamitaheri: Thermal and Thermoelectric Properties of Nanostructures