Predictive and Efficient Modeling of Hot Carrier Degradation with Drift-Diffusion Based Carrier Transport Models
Chapter 8 Role of Electron-Electron Scattering
An important aspect in modeling hot-carrier degradation for short-channel devices is electron-electron scattering (EES) [85, 86, 13, 24]. The steps for calculating the carrier DFs with EES are shown in Figure 8.1. EES populates the high-energy fraction of the carrier ensemble, thereby considerably changing the shape of the carrier DF, which is visible in a characteristic hump pronounced at high energies, see Figure 8.2 [11, 15]. Apparently EES also determines the temperature behavior of HCD in scaled devices [11, 88]. To obtain the carrier DF, the BTE should be solved taking this scattering effect into account. However, it might be possible to approximate the DFs utilizing the knowledge of the effects of EES on the DF. In this section, the drift-diffusion based model for the carrier energy distribution
function, Section 4.3, which was derived to describe hot-carrier degradation in LDMOS transistors, is extended for the case of decananometer devices nMOSFETs with gate lengths of 65, 150, and 300 nm.
Special attention is paid to the effect of EES on the carrier DFs. To approximately consider the important effect of electron-electron scattering on the shape
of the distribution function, a balance equation is solved for the in- and out-scattering rates. The DFs obtained from the extended analytic approach are compared with those calculated with the deterministic Boltzmann transport equation solver ViennaSHE. Both sets of DFs are then used with the hot-carrier degradation model, Section 4.4, to calculate changes in the drain current as a function of stress time. The results of the extended DD-based model are compared with the results obtained using the full solution of the Boltzmann transport equation and with experimental data. The accuracy and limits of the applicability of the DD-based model is also studied on devices with a range of gate lengths. This model allows to avoid the computationally expensive solution of the BTE but is expected to still provide excellent accuracy.
8.1 Extended DD-Based Distribution Function Model
The reference electron DFs have been calculated with ViennaSHE for three different devices of the same architecture with SiON gate dielectric but with different gate lengths of 65, 150, and 300 nm and for three different stress conditions ( = = 1.8, 2.0, 2.2 V). The resulting DFs for the 65 nm device are shown in Figures 8.3 8.5. From these DF curves one
can see that the DF shape substantially changes from the source to the drain. Indeed, DFs evaluated near the source and in the channel have a maximum visible at low and moderate energies with the corresponding energy labeled as in Figure 8.2. As for the drain DFs, they typically feature a Maxwellian tail visible at low energies followed by a plateau. The end of this plateau also corresponds to the knee energy , see Figure 8.2. Another characteristic energy where the DF changes its curvature is related to the onset of the EES hump and is labeled as .
Such sophisticated DFs can be approximated by following the algorithm depicted in Figure 8.1. First, Equation 8.1 allows for a rough evaluation of the electron DF without the effect of electron-electron scattering using the parameters , , and evaluated as earlier in Section 4.3
The parameter assumes a different value after every knee position, see Figure 8.2. In the channel and source regions is assigned a value of -1 before , 1 between and , and 2 after . While for the drain region, has a value of 0.2 before , 1 between and , and 2 after .
The first knee energy in the low energy region is evaluated analytically using the electric field from the DD simulations:
where is the electric field, while , , , and are fitting parameters with values 0.4157 eV, 1.3, 11.04, cmV , respectively. The second knee energy, where EES starts to dominate the high energy tail, is obtained by considering the balance of in- and out-scattering rates. The in-scattering rate corresponds to EES, while the out-scattering includes such mechanisms as interactions of electrons with phonons and ionized impurities [206, 207].
The collision integral, for a nondegenerate semiconductor, is written as [175]:
suggesting that carriers with momentum out-scatter and carriers with in-scatter due to carriers with and respectively. Since there is an equal probability for a carrier being in/out scattered, Equation 8.4 can be written as:
Equilibrium is reached when the in- and out-scattering rates balance.
For the scattering rates in Equation 8.3 the same standard expressions [175, 208] were used in both ViennaSHE and the DD-based model. The scattering rate of carriers by accoustic phonons can be written as [175]
where is the accoustic deformation potential, is the elastic constant ( , is the velocity of sound and is the mass density). For a spherical parabolic band the rate becomes
Note that for a fair comparison the same physical parameters are used in both the SHE code and the analytic approximation.
Optical phonon scattering can be due to absorption or emission of the optical phonons. The corresponding rates are given by Equation 8.8 and 8.9, respectively [175]:
where is the accoustic deformation potential, is the energy of optical phonon, is the density of optical phonons ( ).
The ionized impurity scattering rates are [175]:
where
with being the Debye length given as:
As for electron-electron scattering, we use the following expression [208]:
where is the inverse Debye length. It is important to emphasize that calculation of the rate requires a DF. For energies the DF calculated with Equation 4.22 can be used, while for the DF is perturbed due to EES. Once the knee energies are obtained, the DFs are evaluated by changing the value of the parameter at the knees in Equation 4.22. Same approach is followed for the 150 and 300 nm devices to obtain the DFs in Figure 8.8.