As stated in the beginning of this chapter, BTI is, as its name suggests, a gate voltage, device temperature and time dependent threshold voltage shift ΔVth, which indicates a charge buildup in the oxide of the investigated MOS structures. As the magnitude of ΔVth over stress time ts can be expressed by a power law [130, 131], this spawned a debate over the correct power law exponent. At this stage of research recovery after bias temperature stress has only been superficially investigated and the reaction-diffusion (RD) models [130, 132] could explain the published measurement data. All flavours of the RD models assume a charge buildup at the semiconductor-oxide-interface through hydrogen-based reaction and diffusion. It was soon noticed that the measurement techniques and all time transients differ from publication to publication and that time delays and transients need to be well defined. Additionally, measurement and ΔVth extraction techniques for BTI have been investigated [123]. Utilizing well defined experiments the quantitative relations between stress time, recovery time, device temperature and gate bias or oxide field can be identified. As such it was, for example, found that ΔVth depends quadratically on the magnitude of the oxide field. However, the authors of [115] had to introduce a permanent component of BTI, due to the limited recovery time in the experiments. This permanent component is still under debate and is necessary to reproduce recovery traces, due to the stark asymmetry of stress and recovery over time as can also be seen in Figure 6.1. By asymmetry we mean that the recovery of 1mV of ΔVth after BTI stress takes longer than the buildup of 1mV of ΔVth during stress, provided the device temperature is kept constant. Thus when applying a model to experimentally recorded ΔVth data, the data outside the measurement window is accounted for by a permanent non-recoverable component although a plateau of recovery has to date never been, reproducibly, observed [133]. With the published amount of data, it was finally conclusively found that any reaction diffusion model, particularily on the microscopic scale, cannot possibly explain all the observed BTI characteristics [112]. Finally, with the emergence of time dependent defect spectroscopy (TDDS), it was possible to investigate the charge capture and emission of single defects in small-area MOSFETs and to identify a multiphonon process as the charge exchange process between the substrate and the defects in the oxide.
In the course of this thesis three different models for BTI have been implemented into the drift diffusion simulator MinimosNT. First a simple two well model, second the Two-Stage-Model [134] and last a four state non-radiative multiphonon (NMP) model [4]. In contrast to the four-state NMP model alone, which can reproduce the recoverable component, the Two-Stage-Model predicts both the recoverable and the permanent component of BTI. Thus, the four-state NMP model is often combined with a two well model, since the two well model is used to reproduce the permanent component of BTI. The permanent degradation is usually modelled as a buildup of interface states, whereas the recoverable component is modelled as formation and annealing of oxide traps. Thus the two well model describes the time evolution of Nit, whilst the four-state NMP model describes the time evolution of Not. The Two-Stage model, in contrast to the other two models, includes a description of Not (recoverable) and Nit (permanent) over time. Since the four state NMP model combined with the two well model can explain more BTI characteristics than the other two models alone, this is the model which will be explained in detail.
In order to model the weak temperature and the quadratic oxide field dependence of the charge capture process causing the threshold voltage shifts many phenomenological models have been developed. In addition to that many authors, e.g. [135], tried to explain the charge capture process by extending the successful SRH [136] trapping model with a tunneling coefficient to account for the fact that a charge trapping defect is not necessarily located directly at the oxide interface. However, it was shown that these models cannot sufficiently explain the physics involved in BTI [113], especially the capture and emission time constants found using TDDS. Nevertheless, a double-well model serves as a good example how interface defects caused by BTI can be modelled.
Starting from the SRH model introduced in Section 2.2.2, one can explain the formation of new defects causing a shift in threshold voltage by introducing a model for the trap concentration Nit. To capture the oxide field dependence and the temperature dependence of the formation and annealing of defects an Arrhenius-type law is often employed. In a two-well model for Nit, a defect can be in two states. It can either be active, where trapping and detrapping of charge carriers is governed by SRH statistics (charge trapping), or inactive, where transitions are modelled by defect reaction rates. The barrier height Ea determining the defect reactions rates is assumed to be oxide field dependent and temperature activated. For the model in Figure 6.7 the defect reaction rates are
k12 | = ν0 exp, | (6.2) |
k21 | = ν0 exp, | (6.3) |
Another often employed approach is to explain the defect buildup by initially electrical inactive traps and hole trapping instead. To this end the trap is assumed to reside in the oxide, a distance xT away from the semiconductor-oxide interface and the trapping kinetics are described by the SRH model and elastic tunneling [137]. The modified SRH rates per trap read
k12 = k12SRHλ(x T,Eox) and k21 = k21SRHλ(x T,Eox), | (6.4) |
The first researchers to extend the SRH theory [136] formulated a charge trapping theory which was loosely based on non-radiative multiphonon transitions (NMP) [141]. These models considered charge carrier tunneling and a thermally activated, oxide field dependent process as demonstrated in the previous section [142]. In this model [142] the capture and emission time constants, which are obtained from the reaction rates using Markov-Chain theory [143], read
τc | = τ0 λ(xT)exp, | (6.5) |
τe | = τ0 λ(xT) exp, | (6.6) |
The non-radiative multiphonon theory is based on the possible scattering between multiple phonons with a single electron, which eventually gets trapped in the process. It is theorized that this particular electron can, depending on the configuration and bonds of the surrounding atoms, not only emit photons to loose this gained energy, but now has sufficent energy to occupy a previously unoccupied state (charge trapping). What the atomic configuration at the defect site is, is still not completely clear [112]. In order to derive the reaction rates necessary to describe the charge trapping process lots of approximations to the full Schrödinger equation (cf. Equation (2.1)) describing the many particle problem are required. First the Born-Approximation is applied to separate the equation into two loosely coupled equations: one for the system of electrons and one for the system of nuclei. Per electron i and nucleus j the system of equations is
ψi(r,R) | = V i(R)ψi(r,R), | (6.7) |
ϕiJ(R) | = EiJϕiJ(R), | (6.8) |
| (6.9) |
where
| (6.10) |
is the electronic matrix element, V ′ is the adiabatic perturbation operator, ⟨ψi| is the electronic wave function corresponding to state i. In Equation (6.9) fij is the so-called lineshape function which formally equates to,
| (6.11) |
where the operator avg denotes the thermal average operator over all initial states of the nuclei I. The lineshape function describes the likelihood of a transition from the nuclei state I to J and the electronic matrix element describes the electronic transition probability. In the NMP theory both need to be non-vanishing in order to have a defect reaction. For the lineshape function to be non-vanishing the energies of the final and initial state need to be very close. In Equation (6.11) this has been approximated by the Dirac-Delta. Thus the Franck-Condon factor |⟨ϕiI|ϕjJ⟩|2 is interpreted as a transition probability, determined by the overlap of the two vibrational wave functions ϕiI and ϕjJ (cf. Figure 6.8).
The lineshape function, its formal derivation, numeric evaluation using data from density functional theory (DFT) calculations and subsequent approximations to Equation (6.11) have been extensively discussed previously [112]. For the following discussion the classic approximation of the line shape function is briefly summarized. Parameterizing the adiabatic potentials (cf. Figure 6.8) using the reaction coordinate concept [147] yields
V i(q) | = ci(q - qi)2 + d i, | (6.12) |
V j(q) | = cj(q - qj)2 + d j, | (6.13) |
| (6.14) |
where q1 and q2 are the intersection points between V i(q) and V j(q). Depending on the actual parameters of the adiabatic potentials V i and V j, multiple further approximations are possible to obtain simpler formulae for the lineshape function. The barrier heights from the parameterized potentials, which have been used in the previous equation are
ΔV ij | = 2 | (6.15) |
= 2, | (6.16) |
Depending on the location of the intersection between V i and V j two cases of phonon coupling can be distinguished. When the dominant intersection point q1 (cf. Figure 6.8) is located between the minima of the potentials (qi < q1 < qj) the process is referred to as strong phonon coupling, otherwise it is referred to as weak phonon coupling [113].
To obtain useful transition rates for semiconductors, where a continuum of energy levels must be considered, an integration over the electronic energy is necessary. Employing the definition of the density of states (cf. Equation (2.7)) the rates for capturing a hole (or emitting an electron) from/to the valence or conduction band reads
k12 = | νp ∫ -∞Ev g(E)(1 - fp(x,E,t))A 12(E,xT)f12(ΔE)dE | ||
+ νn ∫ Ec∞g(E)(1 - fn(x,E,t))A 12(E,xT)f12(ΔE)dE. | (6.17) |
k21 = | νp ∫ -∞Ev g(E)fp(x,E,t)A 21(E,xT)f21(ΔE)dE | ||
+ νn ∫ Ec∞g(E)fn(x,E,t)A 21(E,xT)f21(ΔE)dE. | (6.18) |
ΔV 12 | ≈ + + . | (6.19) |
k12 = | σ0 expvthppexp | ||
+ σ0 expvthnnexpexp | (6.20) | ||
k21 = | σ0 expvthppexpexp | ||
+ σ0 expvthnnexp | (6.21) |
In [121] the strong temperature dependence of BTI was explicitly shown. Also in TDDS measurements an oxide field independent, but temperature dependent regime was found [105]. To account for these findings it was theorized that this is due to a structural relaxation of the nuclei forming the defect site after the non-radiative multiphonon transition took place. Due to structural relaxation the effective energy barriers for a subsequent NMP transition change. Such a mechanism can also explain the strongly temperature dependent recovery, where a high temperature accelerates the recovery.
It was theorized and investigated by [148] that a dimer configuration or a hydrogen bridge can be broken by capturing a hole (or emitting an electron) by an NMP transition and that the resulting configuration can structurally relax such that the inverse NMP transition becomes less likely (cf. Figure 6.9). It is assumed that the process of structural relaxation is solely dependent on the temperature. Thus the transition rates are modelled using an Arrhenius law,
k2′2 = ν0 exp and k22′ = ν0 exp, | (6.22) |
NMP as a charge capture and emisson process as well as structural relaxation have been combined by [151] into a four state model (cf. Figure 6.11 and Figure 6.10), which was shown to successfully describe the capture and emission time constants found in TDDS experiments. The model is constructed such that a defect can in each state either exchange charges with the substrate or undergo structural relaxation. However, it is noteworthy that this model does not attempt to explain the step-heights of single defects. Instead the ΔVth contribution of each defect is attributed to the interaction between the charged defect and the random discrete dopants in the channel underneath the defect (cf. Chapter 4).
In order to determine whether or not a defect is charged, a system of equations describing the state transitions per defect are required. To this end, the framework of first-order continuous-time Markov-Chains [143] can be directly applied [29]. A single defect can only be in one state at a time. More precisely Xi(t) is the random variable for state i at time t, which is exactly 1 if the defect is in the ith state and 0 otherwise. The condition that any defect has to be in any of its N states, where N is a finite integer number, can be expressed by
| (6.23) |
Next, the probabilities for state transitions need to be defined. Since a transition from i to j depends on the occupation probabilities Xi(t) and Xj(t), where Xi(t) = 1 and Xj(t) = 0, one can only define conditional probabilities. Formally, for an infinitesimal small time step dt the conditional transition probability from state i to j reads,
| (6.24) |
In the above equation kij are the transition rates, which have been covered in the preceding sections. Inserting Equation (6.24) into Equation (6.23) and defining that
| (6.25) |
is the probability for the defect to stay in a certain state, one obtains
| (6.26) |
Equation (6.26) can now be written as
| (6.27) |
In addition Equation (6.23) can be reexpressed to
| (6.28) |
Equation (6.27) and Equation (6.28) yield N + 1 equations for N unknowns pi(t). Since the system is overdetermined one can ommit a single equation from the equation system. The pi(t) are occupancy probabilities, which determine the probability of the defect to be in state i. Thus the Xi(t) are at a single instance in time either 1 or 0 with the probability pi(t). Sometimes one is only interested in the expecation values of the occupancies. Consequently, applying the expectation operator E{} to Equation (6.27) and Equation (6.28) one obtains
∂tpj(t) = kijpi(t) - kjipj(t) | ⇒ ∂tfj(t) = kijfi(t) - kjifj(t), | (6.29) |
E | ⇒∑ iNf i(t) = 1, | (6.30) |
Putting all of the above together, one is able to write down the system of equations for the four state NMP model from Figure 6.11 and Figure 6.10. The full, time averaged system for the drift diffusion model reads,
ν0 ≈ 1013∕s, λ(xT) denotes the tunneling coefficient and superscripts v and c have been added to account for electron and hole dependent parameters, respectively. This system of equations has more than ten parameters per band edge. One parameter set is the position of the defect in the oxide, where only the trap depth xT directly enters the equations, which together with the dopant positions determines the step-height by electrostatic interaction. The parameters mainly determining the capture and emission times are the energy barriers (cf. Figure 6.11), the trap depth xT, the Huang-Rhys factors Sℏω, S′ℏω′ as well as the curvature ratios R and R′. All other quantities in the equation system are either fixed values, such as the band weights Nc and Nv, or directly obtained from the transport model.
The trapping model is also tightly coupled to the transport model via the trapped charge and charge recombination. The trapped charge Qt is computed via
Qt = ±∥q∥(X2′ + X2)orE{Qt} = ±∥q∥(f2′ + f2), | (6.42) |
Γn{fν n,fν p } | = ∑ defects ∑ ν ∫ k12′p1(1 - p2′)(1 - fν n)gn - k2′1p2′(1 - p1)fν ngnd3k, | (6.43) |
Γp{fν n,fν p } | = ∑ defects ∑ ν ∫ k2′1p2′(1 - p1)(1 - fν p )gp - k12′p1(1 - p2′)fν p gpd3k, | (6.44) |
Rn | = ∑ defectsk12′p1(1 - p2′)n - k2′1p2′(1 - p1)n, | (6.45) |
Rp | = ∑ defectsk2′1p2′(1 - p1)p - k12′p1(1 - p2′)p. | (6.46) |
For comparison with experiment, the capture and emission times τc and τe need to be calculated. Since the four state NMP model has been derived using first-order continuous-time Markov-Chain theory one can use the concept of first passage times from the theory. In the four state model, there are two pathways for charge capture and emission. Charge capture either proceeds via state 1′ or 2′. Thus, following [113], we define the capture and emission times as
τc | = , | τe | = , | (6.47) |
τc(1′) | = , | τc(2′) | = , | (6.48) |
τe(1′) | = , | τe(2′) | = . | (6.49) |