Modeling of Defect Related Reliability Phenomena
in SiC Power-MOSFETs
4.6 Compact Physics Framework (Comphy)
The development of the Compact Physics Framework (Comphy) [203, 211] aimed at reducing the computational effort needed to extract the large number of parameters as required for detailed BTI models, i.e. the four-state NMP or the GSHR model which are partly implemented in commercial TCAD frameworks [232, 233, 234]. However, it has been demonstrated that the mean BTI degradation observed in large area MOSFETs can be explained physically with a two-state NMP model, which requires only a reduced defect parameter set. Comphy and its implementation of the two-state model have been initially applied to reproduce the degradation for a large number of process splits with gate stacks employing SiO2 and HfO2 including different gate contact materials . As a large ensemble of defects can be handled efficiently, it makes the framework particularly suitable for calculating BTI in SiC-MOSFETs, allowing for the extraction of physical defect parameters by using established material parameters. In the following, the main physical models for providing the electrostatic quantities required to calculate the charge transfer kinetics and transient ∆Vth will be outlined. Also, the extension of the framework for efficiently calculating gate leakage currents employing the models presented in the previous section will be presented.
4.6.1 Electrostatic Quantities
In order to calculate the effective trap level of a defect relative to the channel carrier reservoirs, the relation of the applied gate bias to the surface potential at the channel / oxide interface needs to be known. Therefore, by assuming a uniform doping concentration and charge neutrality exists deep in the bulk semiconductor (far away from the interface), the approximation of the surface charge as a function of the surface potential can be used, which reads 
Thereby, denotes the electron Debye length, the semiconductor permittivity and the carrier concentrations at thermal equilibrium, which are derived from the doping concentrations and the band gap . In the case of a charge free insulator, the potential drop across the gate stack is given by
with the work-function difference between the gate and channel materials. As the inverse relation is required to calculate the channel electrostatic quantities like Fermi-Level and carrier densities with employing the Joyce-Dixon approximation , (4.58) is solved numerically for employing an iterative scheme. Figure 4.9 shows a comparison of the approximations used in Comphy to a full numerical solution of the Poisson equation across a poly-Si/SiO2/SiC stack with a Finite-Volume method and Fermi-Dirac statistics.
4.6.2 Threshold Voltage Shift
With the surface potential derived in the previous section, the Fermi-Levels in the channel and gate are known and can be used to calculate the transition rates (4.29) and (4.30) for each input tuple . These allow the computation of the transient defect occupancy using (4.9). By using a simple charge sheet approximation [237, 69], the perturbation of the potential due to defects reads
with being the distance of the defect from the channel/oxide interface.
Often a solution of the occupation for AC signals is sought in order to calculate ∆Vth within circuit simulations. Therefore, an efficient numerical solution for the occupations given by a two-state Master equation for arbitrary shaped periodic two-level bias signals has been proposed by Giering et. al . For the special case of a periodic digital gate AC signal with period at biases at high level and at low level an analytical expression, evaluated after periods, can be found by 
The availability of an analytic expression to calculate the defect occupation for periodic AC signals allows for efficient extrapolation of the device degradation ∆Vth and therefore under operating conditions as will be shown in Section 5.1.6.
4.6.3 Gate Leakage Current Computation
Irrespective of the fact that the Poisson equation is only solved in one dimension within the Comphy framework, the derivation of the correct spatial tunneling distance requires a three dimensional defect distribution. Therefore, based on the input quantities defect density NT and average defect number that should be sampled, a volume with is defined that fulfills . The dielectric thickness is fixed by the gate stack input quantities, and a square area. Within this volume, the defects are sampled using a Monte-Carlo scheme according to their uniform spatial and normally distributed energetic parameters. To avoid low defect densities, in particular those which lead to below 1 nm, defects that result in smaller to any neighboring defect are removed from the drawn sample and re-drawn, i.e. rejection sampling is implemented. The actual number of defects within a slab is further sampled based on a Poisson distribution, as such is observed in small area devices [239, 152], with probability of defects within the slab given by
For the energetic defect parameters (ET ,ER ) samples are drawn based on a normal distribution with ( , ). The expectation values and variations of the current density and threshold voltage shifts can then be calculated from repeating the simulation on slabs, as schematically shown for the simulation space in Figure 4.10.
The computation of the currents in steady state is based on solving the coupled system of equations (4.10) for zero derivative, which reads
with the total in and out rates and at each defect , which consist of reservoir and defect interaction rates. Therefore, for each step with the input tuple ( , , ) all rates as described in Section 4.4 are computed and (4.64) is solved with a Newton scheme. For the transient case, the equation is discretized with an implicit Euler scheme, resulting in the occupations as
The discrete equation can then again be solved with a Newton scheme and the resulting occupations together with the individual rate contributions. With the occupations, the TAT current (4.56) can be evaluated ant the total gate current sums up from the contributions of the TAT and band-to-band current (Tsu-Esaki) (4.48) to a total gate current
With the full modeling framework at hand to compute both TAT currents via single and multiple steps, the significance of a charge transition between two contacts via multiple defects in terms of its contribution to a total leakage current will be evaluated in the following section.