Data ‐ Driven Modeling of Low Frequency Noise Using Capture ‐ Emission Energy Maps

: A new approach for modeling low frequency noise is presented to enable the predictions of noise behavior from negative bias temperature instability (NBTI). The noise model is based on a capture ‐ emission energy (CEE) map describing the probability density function of widely distrib ‐ uted defect capture ‐ emission activation energies. To enlarge the capture ‐ emission energy window and to perform the accurate estimation of the recoverable component of CEE, the Gaussian mixture model (GMM) is applied to the CEE map. This approach provides an efficient identification of noise sources and an in ‐ depth noise analysis under both stationary and cyclo ‐ stationary conditions.


Introduction
In nanoscale CMOS devices, the stochastic behavior of a few defects can have a significant impact on the reliability of devices and their circuits [1][2][3][4][5][6][7]. The negative bias temperature instability (NBTI) is among the most critical reliability issues and can be described by capture-emission energy (CEE) maps [3][4][5][8][9][10][11][12][13][14]. Accurately modeling the NBTI relies on a combined methodology comprising thorough CEE maps and efficient models. However, CEE maps are limited by the stress-recovery measurement window [8,9,[13][14][15][16]. An accurate estimation of the parameters of CEE distribution is crucial to enable full lifetime calculation of NBTI without measurement time extrapolation. The kinetics of defects responsible for NBTI has been explained by the random capture-emission process of charge carriers caused by oxide traps [8][9][10][11][12][13][14]17,18]. Moreover, it has been found that the recoverable component of NBTI and random telegraph signal (RTS) noise are due to same defects [18][19][20][21][22][23][24]. The discrete fluctuations of RTSs with widely distributed time or energy scales are the main source of low frequency (1⁄ ) noise in nanoscale MOSFETs [19,20,[25][26][27][28]. For the statistical analysis and rigorous modeling of noise behavior, it is crucial to incorporate the distribution of defects based on CEE maps, including large signal AC operations [26][27][28][29]. In particular, since the static random-access memory (SRAM) cell is a highly scaled device, the threshold voltage variation due to RTS is a major source of degradation in SRAM characteristics [29]. In the CEE map framework, the time-dependent variability of degradation such as NBTI and noise is reproduced under different bias and duty cycle conditions. In this study, CEE-based noise models are newly presented to allow a quantitative prediction for RTS and 1⁄ noise. The models utilize a Gaussian mixture model (GMM) for enlarging the capture-emission energy window and clustering the recoverable component from the fitted CEE map.

Theoretical Model
The capture-emission process of charge carriers is consistent with non-radiative multi-phonon theory based on linear electron-phonon coupling [1]. In this theory, the capture-emission times are expressed as , exp , ⁄ where is the temperature-independent pre-factor, is the Boltzmann constant, is the temperature, and and are the capture and emission activation energy, respectively. The NBTI-induced threshold voltage shifts ∆ is related to the capture-emission processes of individual defects with a broad distribution of characteristic times. The step height observed in ∆ for small-area devices is due to individual defects, while for large-area devices similar defects are grouped together using a suitable distribution of defects [1,18]. This distribution is visualized in the CEE map, which can be directly obtained from experimental recovery curves of ∆ for various stress conditions by using , ∆ , ⁄ [1]. The defect density map is represented in agreement with bivariate Gaussian distribution, describing the permanent and recoverable components. In the activation energy space, the density , is constructed by the joint probability density function [8,9,13,14] , where ⁄ is a voltage-dependent amplitude with the stress voltage and the constants and . , and , are the mean values and standard deviations of , , respectively. A correlation coefficient is defined as covariance , / . For defects being charged up to the stress time and not yet being discharged at the recovery time , the threshold voltage shift is calculated by integrating over the activation energy map as follows [8,9,13,14] ∆ , , , where , and , are the recoverable and permanent component, respectively, and , , , is the occupancy probability map for the applied stress and recovery waveform. The occupancy probability after and is described by [1] , , , 1 For AC stress and recovery patterns with the signal duty cycle , and the signal frequency f, the probability to be occupied after the high gate-source voltage will be written as [29] , , , Trap-related processes are the most likely cause of the fluctuation in the number of charge carriers. The contribution of multiple-trap levels leads to a noise power spectrum expressed by the superposition of a single Lorentzian noise as follows [26] , 4 ∆ where ∆ , , , , ⁄ is the variance of the fluctuation ∆ in the number of carriers , is the amplitude, is the number of traps following a Poisson distri-bution, and is the angular frequency. With the assistance of capture-emission activation energies, a change in charge state can occur even if the energy levels of the traps are higher or lower than the Fermi level compared to a few , resulting in 1/f noise [30][31][32][33][34]. Assuming that is random variables identically distributed and the effective times constants are statistically independent, the noise power spectrum of the occupancy , , , of the traps in the energy element ∆ and ∆ is given by [33,34] where * is the trap density per unit area. Given that the trap distribution is described by the normalized density * , of the recoverable component, the average noise power spectrum for a device with a width and length is converted to an integral expression , , , where 1 ⁄ 1 ⁄ is the effective time. Note that the term , is similar to the occupancy of traps 1 described by a Fermi function , reflecting the fact that traps with higher occupancy probability contribute to noise power. The equivalent gate voltage noise power spectrum of MOSFETs in the linear region is given by where is the gate oxide capacitance per unit area. Under cyclo-stationary square wave excitation, the noise behavior is governed by the time during the on-state, , and , , and the times during the off-state, , and For an accurate estimation of CEE distribution, a GMM of the distributions appears to be a natural choice for two components of Gaussian distributions. Given a two-dimensional random variable , , with mean and covariance Σ , the GM density function of two Gaussian components is given by [35][36][37][38] | , Σ where is the weight of the i-th mixture component, subject to 0 1 and 1. The i-th independent normal distribution is The expectation maximization (EM) algorithm is used to compute the parameters of each Gaussian component [35][36][37]. It is an iterative algorithm that estimates a maximum likelihood solution of the parameters, based on an expectation (E) step and a maximization (M) step. The GMM is used for clustering Gaussian components, by realizing that two bivariate Gaussian distributions of fitted CEE map can represent clusters. The soft clustering method is applied to estimate cluster membership posterior probabilities [38]. The method assigns a cluster membership score to each data point and ranks the points by their cluster membership score, allowing some data points to belong to multiple clusters.

Results and Discussion
The devices used for measurements were fabricated using a manufacturable remote plasma nitride oxide (RPNO) process. The gate oxide of devices used for measurements was grown with both furnace and rapid-thermal based processes. Following N2-plasma nitridation and post-nitridation anneal of the gate oxide, an undoped poly-Si gate layer was deposited. For p+ gates, B was implanted, and rapid-thermal annealing was then performed in N2 for different durations. The threshold voltage shifts are obtained on pMOSFETs (W = 10 m , L = 0.14 m , 2.2nm), performing the ultrafast measure-stress-measurements [8,9]. With a measuring delay time of 1 μs , the temperature-accelerated measurement scheme was used to extend to the experimental stress-recovery window of 10 s. Noise measurements are conducted using an HP 3582A spectrum analyzer (1Hz-25kHz) and a Brookdeal 5004 low noise voltage amplifier connected at the output terminal. A bias resistance at the output terminal is chosen to be 10 times larger than the small signal channel resistance. The noise power spectrum is adjusted to subtract the noise due to the measurement system, which consists of the noise due to the amplifier and the bias resistance, and the thermal noise of device under test. Figure 1a shows experimental recovery data of ∆ following stress time with numerical data for 2.0V and 125 ℃. Numerical data of ∆ are obtained by integrating the CEE distribution in (2), which is calculated by taking the derivatives of measured data, as shown in Figure 1b. To investigate how the trap distribution influences RTS noise, the integrand of for the recoverable component in (7) is numerically calculated and visualized as a function of and under stationary DC condition at a particular frequency. It is shown in Figure 3a that a peak is distinct at 1.10eV and 0.65eV for 1Hz, and moves to relatively low activation energies 0.75eV and 0.45eV as the frequency increases to 1kHz. As illustrated in the inset of Figure 3a, these results arise from different RTSs, indicating that faster components of shallow trap states for 1kHz are situated close to lower activation energies [39]. To examine the behavior of cyclo-stationary RTS, the times are modeled as , • , and , , ⁄ , where the factor is assumed to be greater than one [27,39]. In the off state, and are multiplied and divided by , as they increase and decrease with a decrease in the gate voltage, respectively [27,39]. For 50, a peak moves to higher emission energy , as shown in Figure 3b.
1⁄ noise power is calculated by using (7) with the recoverable component obtained from GMM with soft clustering. It is observed in Figure 4a that the recoverable component of CEE map gives a 1⁄ spectrum form with a slope of −1. For asymmetric RTSs ( ,   ,   ⁄ , cyclo-stationary condition may lead to an increase in noise power [27]. Figure 4b shows a linear relationship between and ∆ after a stress time 1000s. The same strong correlation is also observed for the device with SiON/poly-Si gate [19].

Conclusions
In conclusion, with a recoverable component of NBTI stress-recovery, a new noise model is formulated based on a fitted CEE map with soft clustering of GMM. The proposed model offers the following unique features: 1) The model is based on the occupancy of defects in the capture-emission energy map. Therefore, it allows modeling of the noise under both stationary as well as cyclo-stationary excitation conditions. 2) The GMM with soft clustering enlarges the CEE window and performs independent estimations of two distinct components of CEE map, enabling the accurate calculation of noise characteristics. 3) New theoretical expressions can capture the noise behavior in terms of activation energies and provide the identification of noise sources and the reliable prediction of 1⁄ noise from NBTI.

Data Availability Statement:
No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest:
The author declares no conflict of interest.