Modeling of Defect Related Reliability Phenomena
in SiC Power-MOSFETs
4.3 Non-Radiative Multi-Phonon Theory
As discussed in the previous section, the SRH model fails to describe charge transfer to and from oxide defects due to its purely elastic tunneling character that results in a “tunneling front" from the gate channel into the oxide and
underestimation of defect charge transition time constants in thin oxides [73]. The missing model component within the SRH picture is the structural
transformation of the defect upon charge capture. In order to compute exact transition rates for both radiative and non-radiative transitions, a full quantum mechanical treatment of the defect and the surrounding phonon bath is
required [209]. However, such an approach is computationally unfeasible when applied to a large ensemble of defects with a wide variation of
structural properties, as is the case in large area power MOSFETs with amorphous insulators. In this section, the detailed derivation of classical NMP transition rates from the exact quantum mechanical solution as given by
Waldhoer [210] will be iterated. In such a form, the NMP model allows for transient computations for a large defect ensemble as is
required to study charge trapping on a device level. It should be noted that these classical approximations of the transfer rates have shown to be able to capture the main features of charge trapping related effects in a large number
of Si-based [211], as-well as novel 2D [140] and wide-band
gap devices [158, 212]. Unlike in radiative transitions,
in which the excess energy to overcome an energetic barrier is induced by a photon, both the energetic activation of the carrier and the dissipation of the excess energy upon relaxation to the final state’s equilibrium configuration is
due to interaction with multiple phonons. Therefore, for the derivation of the non-radiative multi-phonon transition rates, both electronic states and vibrational states
as well as their mutual interaction have to be considered. By using the simplifications induced by the Born-Oppenheimer approximation [213], the transition rates from an initial vibronic state
to a final state
can be calculated with Fermi’s Golden Rule [214]
with the state eigenenergies ,
and the matrix element
Within the Born-Oppenheimer approximation, the first order perturbation Hamiltonian can be separated in an electronic
and vibrational
part. This approximation is justified due to the large difference in timescales, at which the electrons and nuclei respond to changes of their coordinates and momenta within the system and the matrix element
becomes [210]
Due to the orthogonality of the different electronic wave-functions, the second term in (4.18) vanishes and the matrix element is determined by the overlap integral of initial and final vibrational states multiplied by the electronic matrix element, as stated by the Franck-Condon principle [215, 216].
As for typical operating conditions a large number of different vibrational modes will contribute to the overall rate, the thermal average over all partial rates (4.16) needs to be considered. With the initial state assumed to be in thermal equilibrium this results in the expression of the rate by a product
with being the electronic matrix element describing the electronic coupling between the defect wavefunction and the charge reservoir. The line-shape function
describes the vibrational interactions in the classical high temperature limit, in which the difference between the vibrational energies
is much smaller then the thermal energy
and the classical barrier
defined by the crossing point between the Potential Energy Surfaces (PESs) of the initial and final state (c.f. Figure 4.3), i.e.
, and can be approximated by [210]
The exponential term typically dominates , whereas the prefactor
, which depends on the particular PES shape, can be neglected [210].
Up to this point, the NMP formalism is fairly general. In order to reduce the complexity introduced by considering a full PES, the harmonic approximation of the one-dimensional Potential Energy Curve (PEC) of the neutral and negative state at their energetic minimum is used, as shown in Figure 4.3 (left). While, the PES shape is in general arbitrarily complex, the assumption of parabolic PECs has been shown to deliver a reasonably accurate approximation for important defect candidates in Si MOSFETs [210, 217]. For the calculation of charge transfer in MOSFETs, the rates (4.19) need to be calculated for the interaction with a carrier reservoir, which can be the channel conduction (valence) band or the gate contact. As the most relevant degradation mechanism for power switch applications is PBTI in a nMOSFET, the rates are derived for charge transfer between the channel conduction band and a defect, as shown in Figure 4.3 (right).
Figure 4.3: Harmonic approximations of the Potential Energy Curves (PECs) of two defect states along the dominant reaction coordinate for charge transfer between the channel conduction band and the defect. In the classical high temperature limit the barriers
and
are determined by the intersection point. These barriers are uniquely defined by the relaxation energy
, the trap level
relative to the conduction band edge
and the curvature ratio
. Note that the minima of the parabolas shift relative to each other with different field strength across the oxide (c.f dashed lines).
The system’s total energy in the neutral defect state is given by
with the curvature of the neutral PEC and the reaction coordinate difference between the minima
. The energy of the electron
in the channel conduction band is included in the first term
and as indicated in Figure 4.3 each band state the electron can dwell in results in a separate PEC for the neutral state. When an electron
with energy
is captured from the conduction band, the defect becomes negatively charged and the total energy reads [210]
As only energy differences are relevant for calculating the transition rates, cancels out. With the additional definition of the thermodynamic trap level
the transition rate (4.19) becomes a function of
and ET as [210]
For a continuous energy state reservoir an integration over the energy space is required yielding [210]
Thereby denotes the density of states in the conduction band and the electron (hole) occupation
(
) of the reservoir needs to be weighed by Fermi-Dirac statistics [210]
Finally, to solve the integral (4.24) over the band states, further approximations are necessary, in order to reduce the numerical expense. First, by using the bandedge approximation, which assumes that the carriers within a semiconductor are predominantly located around the bandedges, the matrix elements and line-shape function can be factored out of the integral [210]
The left over terms in the integral simply yield the electron concentration
Furthermore, the electronic matrix element is typically approximated by a simple tunneling factor [74], which is reasonable given the
strong localization of the defect wave functions [210]. This approximation yields [72]
with the electron thermal velocity , the capture cross section
and the tunneling probability, which is typically calculated by using a Wentzel-Kramers-Brillouin (WKB) approximation [218], which gives analytical expressions for trapezoidal and triangular barriers. Together with the line-shape function in the classical limit, we
arrive at the analytical expression for the electron capture rate from the conduction band edge [210]
Analogously and by using the relation the electron emission rate is given by [210]
In the same fashion analytical rates for electron capture (emission) can be calculated for the interaction with the channel valence band (
) and the gate contact
(
). The total capture and emission rates for solving the Master equation (4.4) are then given by the sum of all capture
and the sum of all emission partial rates. As can be seen from the equations (4.29) and (4.30), the rates are dominated by the barriers
, which can be analytically calculated from the intersection point of the PECs in the harmonic approximation (c.f. Figure 4.3) by [73]
with the relaxation energy and the square root of curvature ratio
. The reverse barrier
can be calculated by
with the energetic difference between the PECs minima
and
the electric field strength within the oxide. The left panel of Figure 4.3 thereby illustrates how the minima of the PECs are shifted
relative to each other by applying
across the oxide in the gate stack, which imposes the gate bias dependence of the energy barriers. Note, that this is the main difference of the Grasser two-stage NMP model, compared to the Kirton and Uren
model [67], which extended the SRH model by a bias independent Boltzmann factor, which is only a good
approximation in the case of interface defects [73].