(image) (image) [Previous] [Next]

Charge Trapping and Single-Defect Extraction in Gallium-Nitride Based MIS-HEMTs

Chapter 5 Defect Modeling

In this chapter, the essential models which will be used to describe BTI phenomena throughout this work are presented. After a brief explanation of CET maps being a clever way to visualize experimental eMSM data in Section 5.1, the theoretical foundation of the NMP theory to describe charge trapping is discussed in Section 5.2.1. Some evidence for the importance of phonon-assisted charge transitions in oxide defects and the framework of the well-known NMP four-state model and the derivation of equations for the time-constants is covered in Section 5.2.2. Section 5.2.3 introduces an extension to the NMP theory for charge trapping in semiconductors by also including the interactions with the local bands.

5.1 Capture-Emission Time Maps

When performing eMSM measurements, the recorded stress and recovery traces contain detailed information about the active defects at a given temperature and bias condition [100]. One way to visualize the capture and emission times of the active defects and their impact on (math image) is using CET maps [115, 118, 119]. In the case of large-area devices, the response of the individual active defects per area \( \Delta \tau ^2 \) is grouped together into a density \( g_{ij} \) defined as [110]

(5.1) \begin{equation} g_{ij}=g(\tauc _{,i} ,\taue _{,j}) = \sum _k{\frac {\eta _ka_k}{\Delta \tau ^2}\mathrm {rect}\left (\frac {\tauc _{,k}-\tauc _{,i}}{\Delta \tau ^2}\right )\mathrm {rect}\left (\frac {\taue _{,k}-\taue
_{,i}}{\Delta \tau ^2}\right )}, \eqlabel {model:gij} \end{equation}

with \( \eta _k \) being the individual step height, \( a_k \) the maximum occupancy change of the defect and the rectangle functions

(5.2) \begin{equation} \mathrm {rect}\left (\frac {\tau }{\Delta \tau }\right ) = \begin{cases} 1\quad |\tau -\Delta \tau | \leq 1/2\\ 0\quad \mathrm {otherwise} \end {cases}.   \end{equation}

(image)

Figure 5.1: Instead of considering a large number of defects, similar defects around around \( (\tauc _{,i} ,\taue _{,j}) \) are grouped into \( g_ij \) defined in (5.1). The size of the circles represents the impact of the defect on (math image), which is given by its step height \( \eta _k \) times the maximum occupancy change \( a_k \) (from [110]).

With this definition, (5.1) simply collects all defects within a certain area of (math image) and (math image) into \( g_{ij} \) (see Figure 5.1). The change in the occupancy after stressing for \( t_\mathrm {s} \) seconds and \( t_\mathrm {r} \) seconds of recovery, normalized by the maximum change in occupancy is given by [110]

(5.3) \begin{equation} h(t_\mathrm {s},t_\mathrm {r};\tauc _{,i},\taue _{,j}) = \frac {\Delta f(t_\mathrm {s},t_\mathrm {r};\tauc ,\taue )}{a_k}=(1-\e ^{-t_\mathrm {s}/\tauc _{,i}})\e ^{-t_\mathrm {r}/\taue _{,j}}.   \eqlabel
{model:h} \end{equation}

The total degradation can be calculated by summing over all capture and emission times.

(5.4) \begin{equation} \dVth (t_\mathrm {s},t_\mathrm {r}) \approx \sum _i\sum _j{g(\tauc _{,i} ,\taue _{,j})h(t_\mathrm {s},t_\mathrm {r};\tauc _{,i},\taue _{,j})} \eqlabel {model:dVth} \end{equation}

In the limit \( \Delta \tau \rightarrow 0 \), the rectangle functions in (5.1) can be replaced by Dirac functions. If the exponential terms in (5.3) are replaced by two unit step functions around (math image) and (math image), one obtains a simple relationship between between (math image) and \( g \) [110]:

(5.5) \begin{equation} \dVth (t_\mathrm {s} ,t_\mathrm {r} ) \approx \int _0^{t_\mathrm {s}}{d\tauc } \int _{t_\mathrm {r}}^\infty {d\taue }\ g(\tauc ,\taue ) \eqlabel {model:CETVth} \end{equation}

In other words, the degradation is obtained by summing up all the defects being charged until \( t_\mathrm {s} \), but not yet discharged after \( t_\mathrm {r} \). This means that the \( g \) can easily be obtained directly from the measurements by calculating the mixed partial derivative of the measured stress and recovery traces [110] as

(5.6) \begin{equation} g(\tauc ,\taue ) \approx -\frac {\partial ^2\dVth (\tauc ,\taue )}{\partial \tauc \partial \taue } \end{equation}

Thus, the density \( g \) in the CET maps stands for the capture and emission times of the defects weighted by their individual impact on (math image). It should be noted that the data in a CET map in general is only valid for a particular set of stress and recovery biases at a specific temperature. In large-area devices where a lot of defects contribute to (math image), CET maps are usually used to visualize the average response of a large number of defects for a certain technology at a certain stress condition. An example for a CET map obtained from eMSM measurements on GaN MIS-HEMTs is given in Figure 5.2.

(image)

(image)

Figure 5.2: The recovery data recorded in eMSM measurements (left) can be directly used to calculate the CET map (right). With equation (5.5), the density plot in the CET map can be used to calculate the (math image) drift for an arbitrary stress and recovery time (from [75]).