# Modeling of Defect Related Reliability Phenomena in SiC Power-MOSFETs

#### 4.5 Efficient Framework for MOS Gate Leakage Currents

As stated in Chapter 1, several leakage current mechanisms through the dielectric in a gate stack are possible. In Figure 4.7, the most common tunneling mechanisms are schematically shown, with direct and Fowler-Nordheim Tunneling currents through an energetic barrier and trap assisted currents via one or multiple defects. In this section, a standard model for calculating currents through energetic barriers together with a derivation for the charge hopping mechanism following  based on a generalized Ramo-Shockley theorem  will be re-framed to be applicable for a MOS gate stack.

##### 4.5.1 Tsu-Esaki Model

The tunneling current through an energetic barrier that separates two electrodes can be computed with the so called Tsu-Esaki formalism  , in which the current density at the semiconductor channel for electrons and holes is expressed as with the effective electron density of states (hole) mass in the semiconductor in the plane parallel to the interface , the elementary charge , the Planck constant , a tunneling probability and the so-called supply function . Assuming Fermi-Dirac statistics at both contacts, the supply function for electron tunneling is given by In order to efficiently calculate the tunneling probability, the WKB approximation can be applied for thin and energetically high barriers as with denoting the effective electron mass in the dielectric.

Thereby, the integral is carried out from position to over the energetic barrier with shape . In the case of a triangular barrier, as shown in Figure 4.8 (left), (4.50) evaluates to the commonly known FN formula for electron tunneling with the potential barrier , the electron energy and the electrical oxide field strength . Also, for a trapezoidal barrier (Direct Tunneling (DT)) as shown in Figure 4.8 (right) with the barrier heights an analytic expression evaluates to From (4.48) and (4.50) it can be seen that the material parameters, i.e. the energy barrier and effective tunnel masses, within the exponent in the WKB factor strongly influence the current density in the inversion and accumulation regime, in which the supply function (4.49) is non-zero. Together with the exact knowledge of the oxide thickness and the dielectric constant , these parameters determine the accuracy of the Tsu-Esaki computation.

Upon careful parameter calibration, the capability of the model to calculate accurate tunneling currents through defect free dielectrics has been demonstrated. By considering this best case scenario, the suitability of dielectrics for new material combinations has been demonstrated with our implementation in Comphy for a large number of compounds considered as insulators for two-dimensional channel materials in [CSJ7]. It was thereby found that hexagonal boron-nitride (hBN), which was often proposed as dielectric for 2D materials, does not meet the requirements of low gate leakage current for pMOS fabrication even in the defect-free case [CSJ7]. As for real devices the leakage current is eventually further enhanced due to trap-assisted tunneling, the next section outlines the calculation of such a hopping current.

##### 4.5.2 Charge Hopping Model

As for a charge hopping current the local current density does not fulfill the continuity condition, i.e. charges “vanish" at one defect site and “pop up" at another instantaneously. In this case, the current at the gate contact can be written as with the integral taken over the device surface at the gate contact area . By selecting a test-function as the solution of the Laplace equation at the gate contact and for all other contacts and by using the divergence theorem, the integral over the whole device volume reads The left term in the right hand side integral thereby vanishes due to the choice of and a conduction current of at the contact. Due to the local violation of the continuity condition in the Pauli-Master equation , the current is not divergence free and is given on a one-dimensional grid between points and as by using (4.10) with being the occupations of the two trapping sites. For a bias difference of between gate and channel interface, the test function evaluates to with the oxide thickness . Inserting (4.55) into (4.54) the current at the gate contact evaluates to the generalized Shockley-Ramo theorem [231, 226], which yields the displacement currents observed at the gate due to each charge movement within , to with the oxide capacitance . The rates in (4.56) consist of both reservoir to defect interactions for the TAT current and charge trapping current contributions and defect to defect rates for the multi-TAT current contribution. This charge hopping current model has been implemented in Comphy, as will be discussed in the next section.