Abstract

This paper proposes a poststress time evolution model for sub-10-nm thick SiO2 films for degradation prediction and the extraction of trap-related parameters. The model is based on the understanding that the degradation in thin SiO2 films continues within the unstressed interval. The phenomenon is captured by an analytical expression that indicates that the time evolution of SiO2 film degradation roughly consists of two stages and that the degradation is more likely to occur if water molecules are present. It is demonstrated that the simple analytical model successfully reproduces measured results. It is also suggested that the degradation process considered here is related to oxygen diffusion in the resistive transition process.

1. Introduction

It is well known that stress-induced leakage current (SILC) is observed in sub-10-nm thick films after the application of constant-current stress (CCS) or constant-voltage stress (CVS) [1, 2]. In this thickness range (4~8 nm), oxide breakdown and oxide degradation are among the most important factors influencing reliability. SILC observations suggest that traps are created in the film over time [1, 2]. The SILC of tunnel oxide films is, at present, the major factor impeding the downscaling of nonvolatile memory devices because a high oxide field is needed if the film is to permit data to be written and read reliably [3]. According to the ITRS Roadmap, the tunnel oxide film thickness () should be held to around 6 to 7 nm [4]. However, it is inevitable that sub-6-nm thick films will be needed given the continuing decrease in supply voltage. Furthermore, recent works have reported that SILC can be related to the increase in neutral electron traps generated by hot electrons [5], and dielectric breakdown occurs as soon as a certain critical density of neutral electron traps in the oxide is reached [6, 7]. Although an electrical stress experiment has already been performed to investigate the time evolution of SILC events, the importance of the unstressed interval has not been well considered [5].

Recently, the author proposed a physical model of the time evolution of SILC events within the unstressed interval by characterizing the post-SILC degradation of films [8]; the diffusion of and molecules induces many Si–OH bonds (neutral center) which influence electron conduction, after the unstressed time interval (), through the neutral electron traps, such as Si– and SiO–, which are generated by the electrical stress [8]. However, the model was rather complex since it took account of the many chemical processes possible.

Silicate films and transition metal-oxide films are also attracting attention due to the resistive transition phenomenon [911] which should yield the creation of new nonvolatile memory structures [12, 13]. Therefore, a full characterization of silicon oxide film is very important in deeply understanding the mechanism of resistive transition and the resulting performance degradation.

This paper develops a feasible time evolution model to predict the breakdown due to the stress-relaxation process during the unstressed period by analyzing the degradation of metal-oxide-semiconductor (MOS) capacitors with 5.2-nm thick films. The paper addresses extra active trap creation within the unstressed interval. Possible theoretical models are proposed and compared to measured results; they suggest that degradation proceeds in two stages, and several important physical parameters are extracted from the experiment’s results. Finally, we address oxygen diffusion in the resistive switching process of film [9, 10].

2. Device Fabrication and Stress Method

The MOS capacitor (.2 nm) examined here has an -type poly-Si//Si substrate stack structure; its -type (001) silicon substrate has resistivity of 4 Ωcm. The film is fabricated at 950 C by the rapid thermal oxidation technique, followed by the chemical vapor deposition of phosphorus-doped amorphous silicon [14]. Wet etching is used to pattern the poly-Si film into a gate electrode (150 μm × 200 μm) to minimize processing damage.

In the present experiment, a short electrical stress was imposed on the MOS capacitor as shown in Figure 1. The CCS condition was 500 mA/c and was applied to the film until a fixed number of electrons were injected. The density of electronic charge () injected during one CCS cycle was either 6 C/c or 12 C/c. After one CCS cycle, the gate electrode was held open for various periods () at room temperature. The current-voltage (-) characteristics of the post-CCS film were at first evaluated just after the CCS cycle in order to observe how the SILC develops; was evaluated, where denotes the stress time. After the unstressed interval , - characteristics were again measured in order to confirm subsequent SILC development; was evaluated. Plotting the curve of as a function of gate voltage () yields a specific peak at of 3.6 V (not shown here) [8, 15]. Accordingly, we use this value of as representative of SILC hereafter.

In the following section, the time-to-breakdown () is taken to be the integration of stress time and the hard breakdown (HBD) event is defined as the time at which the gate voltage stepped down by half during the CCS test. In order to determine the Weibull slope, , at least 50 capacitors were measured for each case and the variation, taken to be the difference between the highest and the lowest confidence bounds of the value, is plotted as a function of sample size by using the maximum likelihood estimation at the 95% confidence level.

3. Results and Discussion

3.1. Statistical Consideration of

In conventional CCS tests, the electrical stress is continuously imposed until hard breakdown occurs. In the present experiment, however, the electrical stress is periodically imposed on the MOS capacitor as shown in Figure 1. After one cycle of CCS, the gate electrode is held open for one week at room temperature. This sequence was repeated until the MOS capacitor exhibited HBD; hereafter we call this “periodic stress method with unstressed interval ”. The proposed method is different from the conventional unipolar-pulsed stress [16, 17]; the duty cycle is 50% in the conventional experiment, while it is much smaller than 50% and the frequency is much less than 1 Hz in the present one.

Weibull distributions of are shown in Figure 2. Three different current stress categories are compared; the first one is the conventional CCS without any unstressed interval, the second one is the PSMI with of 6 C/c and one week unstressed interval, and the third one is the PSMI with of 12 C/c and one week unstressed interval. In Figure 2, the PSMI without net is used to provide a baseline for the basic periodic CCS method. In the PSMI method without any unstressed interval, the CCS ( = 6 C/c or 12 C/c) was repeated without any unstressed interval. It should also be noted that all plots are identical to each other in Figure 2. Consequently, it can be concluded that the conventional CCS method and the PSMI without net exhibit no substantial difference.

Next, we investigated how the finite value of of PSMI influenced oxide breakdown. Weibull distributions of under PSMI with of 12 C/c and one-week unstressed interval are shown in Figure 3(a), where distributions with and without net are shown for comparison. The two Weibull distributions of in Figure 3(a) cross in the center of the distribution, and the slopes are very different. The Weibull slope for the PSMI with the one-week unstressed interval is steeper than that without the net interval. This suggests that uniform degradation occurs inside the film over the one-week unstressed interval because it is expected that the initial total energy accumulated at the molecule by the electrical stress yields many traps during the unstressed time interval as discussed in the following section; this degradation process should accelerate the breakdown and then the Weibull slope is steeper than that in the conventional CCS. It can be considered that the traps created in the active CCS period () generate some extra traps around themselves in the . This would be possible when defects created by the CCS relax their structural strain. Therefore, note that the poststress spatial extension of CCS-induced traps is a secondary result of failure events around initial trap sites [18]. Although the reason, why the average value for the PMSI experiment is larger than that for the conventional CCS experiment is not clear, I guess that the PSMI may lead to a certain degradation phase different from the conventional degradation that yields the SILC [8].

Weibull distributions of under PSMI with of 6 C/c and the one-week unstressed interval are shown in Figure 3(b), where distributions with and without net are shown for comparison. The two Weibull distributions of in Figure 3(b) are parallel to each other; the slopes are almost identical to those in Figure 2; the Weibull distribution for PSMI with the one-week unstressed interval is shifted to the left. This suggests that the one-week unstressed interval advances the percolation process of traps inside the film. In other words, it can be considered that the area of trap states created in the active CCS period () expands due to some additional mechanisms [18] and generates some extra traps around themselves within . It is anticipated that the spatial expansion of trap states reduced the average value of for the PSMI method because it should accelerate the percolation process of traps. This is possible when the structural strain induced by defects created by CCS relaxes; this is also a secondary result of failure events around initial trap sites [18].

As a result, the above results strongly suggest that the generation of extra traps expands in three dimensions. Therefore, it is anticipated that such phenomena observed in this paper are related to those observed in [18]. However, we cannot conclude that the degradation observed in this study is the protrusion in the film[19].

3.2. Time Evolution of SILC Component

We evaluated the time evolution of the increase in gate current after one CCS cycle [8]; the increase is defined by , where is 12 s for = 6 C/c or 24 s for = 12 C/c. In order to determine the normalized increase, , at least 5 capacitors were measured using various values of . In each case, the experimental error was 4.6% for = 6 C/c and 3.5% for = 12 C/c; these error values are within the allowable margin. It has already been described that the increase in SILC is directly related to the increase in the number of traps yielding SILC under the stress applied [5]. Actually, the number of neutral traps generated is proportional to the injected charge, and the increase in SILC is proportional to the number of neutral traps (not shown here). dependence on is shown in Figure 4. It can be seen that increases and then saturates as increases. This suggests that extra active traps contributing to SILC are created after CCS [8].

We analyzed the time evolution of . We assumed that initial traps created by CCS have a certain, large, internal energy and that the energy is dissipated as structural relaxation occurs around the initial defects. Thus extra trap creation should stop when the internal energy of the initial defects falls below the threshold value of trap creation; in other words, should saturate. Since the structural relaxation of the defect region induces additional traps surrounding the initial defect region, stemming from the extra traps can be characterized by the effective volume of structural relaxation region .

3.3. Possible Physical Models

In this paper, we reconsider the above empirical understanding of the degradation. Under electrical stress, some of the electrons injected from the electrode lose energy through inelastic scattering processes. Molecules accumulate energy by these events (detailed in the appendix). After the electrical stress is stopped, molecules having any extra energy transfer the energy to surrounding molecules, which yields electronic traps at the surrounding molecules.

It can be anticipated that the differential equation describing this phenomenon is given by where is the effective volume enclosing the created traps (c in unit), is the energy per volume accumulated at the specific molecule by the electrical stress integration (eVsec/c in unit) (see appendix), is the initial total energy accumulated at the molecule by the electrical stress (eV in unit), is the energy dissipated by a single new trap creation (eV in unit), and is the trap density created (c in unit). When it is assumed that has the initial value of ~0 c, the solution of (1) is obtained as where is the time constant characterizing the saturation of SILC increment. The left-hand side of (2) means the total trap number at time . Equation (2) indicates that the total trap number saturates as the unstressed period () increases. Assuming that the increment of SILC component is proportional to according to the consideration described before (1), calculation results are shown by lines in Figure 4 together with measured results, where the parameters of best fit are used. For Figure 4 it is assumed that where is a constant. Calculated results successfully reproduce the increment in SILC component in the early stage of the time evolution; the early stage suggests that the reaction process is determined by just the time constant ; consider  sec for 2 C/c and  sec for  C/c. The time constant means the averaged “total time” in which each trap is entirely created by structural relaxation. However, the calculated curves do not reproduce all of the measured results, so a second process must be posited. In the previous paper [8], the author demonstrated that the stage of medium degradation is ruled by the diffusion processes of molecule species and [8]. Therefore, it can be considered that is not proportional to just extra trap density and as shown in Figure 4. This means that the extra active trap creation process must involve another differential equation because the degradation mechanism is not a simple reaction process.

We investigated how to reproduce the measured results by applying an analytical model. Following the previous consideration, it is anticipated that the normalized degradation of SILC is proportional to the normalized conductance degradation as shown by [5] where is the conductance of the poststressed film and is the incremental conductance of the film. When the effective radius of the initial defect is given as , the defect volume is given by . It follows that the assumption that the incremental SILC is ruled primarily by the in-depth expansion of the extra trap region [20] yields the following relation: On the other hand, the assumption that the incremental SILC is ruled primarily by the in-plane lateral expansion [21] of the extra trap region yields the following relation:

(i) Possibility of One-Dimensional Percolation Process of Trap Creation. Combining (2), (5), and (6), we have where is a constant, is the time constant characterizing the saturation of the incremental SILC, and is the energy per volume accumulated by the electrical stress integration (eVsec/c in unit) assuming the one-dimensional percolation process of trap creation.

Calculated curves are shown in Figure 5 together with the measured results; the value of is about  sec for the injected charge of 6 C/c and about  sec for the injected charge of 12 C/c. The proposed model roughly traces the second stage demonstrated by all measured results. As a consequence, the relation can be used for . This suggests that the increase in SILC is related to the one-dimensional percolation process of trap creation [20]. However, if we attempt to assume that the second stage is primarily ruled by a mechanism independent of that of the first stage, we must note that (8) does not comprehensively explain all measured results by itself. In this paper, we assume that the second stage is composed of two different mechanisms (separate from the mechanism appearing in the first stage). Therefore, we can consider that the second stage of degradation is primarily ruled by the two-dimensional expansion of extra trap creation.

(ii) Possibility of In-Plane Diffusion Process of Trap Creation [21]. Combining (2), (5), and (7) yields where is a constant, is the time constant characterizing the saturation of the incremental SILC, and is the energy per volume accumulated by the electrical stress integration (eVsec/c in unit) assuming the in-plane diffusion process of trap creation.

In order to examine whether the two mechanisms modeled by (2) and (10) rule the process, the following equation is assumed: Here, it is assumed that for simplicity. Calculated curves are shown in Figure 6 together with measured results. The value of is about  sec for the injected charge of 6 C/c and about  sec for the injected charge of 12 C/c. The proposed model (12) reproduces the early-to-final stage of all measured results. As a consequence, the relation can be used to explain the primary mechanism of the second stage. This strongly suggests that the increase in SILC is related to the two-dimensional percolation process of extra trap creation, in other words, the diffusion process of some specific material as suggested in [8, 22].

3.4. Extraction of Physical Parameters

(i) Using (A.2) and (A.3), we can estimate the values of and . They are given as = 5.5 ×  [J] (for = 6 C/c); = 1.1 ×  [J] (for = 12 C/c); = 5.2 × 109 [Js/c] (for = 6 C/c); 2.1 × 1010 [Js/c] (for = 12 C/c).Although the estimated values of seem large, they lie in an acceptable range when we estimate the number of electrons scattered by atoms under electrical stress. Values of suggest that the energy accumulated by one-second of electrical stress is ~1 J/atom, and we can consider that they are also in a reasonable range.

Next we estimate the trap density created () using (A.4). When we use the values of extracted from experimental results ( s for = 6 C/c and  s for = 12 C/c), we have = 1.7 × 103 [J/c] (for = 6 C/c), = 1.4 × 104 [J/c] (for = 12 C/c).Using estimated values, we have > 3.1 × 1016 [/c] (for = 6 C/c), > 1.3 × 1017 [/c] (for = 12 C/c).These values mean that the density of traps created is much larger than 1010 c. These values of physical parameters are acceptable based on various SILC-related experimental results in past articles [13, 58, 17, 2023].

(ii) When it is assumed that some material is diffusing within the film, the diffusion constant () can be estimated by , where is the radius of the in-plane expansion of the extra trap creation. This yields the following important relation between the diffusion constant () and the trap density ():We have already found that the expansion of extra trap creation is three-dimensional [23]; this suggests , where is film thickness. Thus, we haveSince  s to  s, the estimated diffusion constant of damaged film examined here has a value of to  c/s. In unstressed film, it is known that the diffusion constant of molecules is about  c/s [22]. As suggested in [8], the diffusion constant of molecules in a poststressed film in a certain unstressed interval should be about ten times larger than that in an unstressed film [8]. Therefore, the analysis proposed here is reasonable and very useful for extracting several parameters from measured results.

According to related studies [8, 22], molecules are the most likely material at room temperature. It is already known that many neutral traps are created just after electron injection, so it is anticipated that those neutral traps are activated by being terminated with H or OH species; that is, O–Si–H and HO–Si–O structures are created [8]. From this discussion, it can be considered that all traps are created as neutral traps just after stress termination and that molecules activate the neutral traps by bonding with them [8]. We consider that Figure 6 demonstrates the veracity of this simplified model in assuming the diffusion of molecules through the film at room temperature.

In the case of the resistive transition phenomenon, the electric-field-induced diffusion of oxygen ions is usually assumed [8, 24, 25]. The above consideration strongly suggests that the oxygen diffusion is supported by hydrogen molecules [8] and simple ionic transport.

4. Summary

In this paper, a poststress time evolution model for sub-10-nm thick films was proposed for degradation prediction and the extraction of trap-related parameters. The model is based on the fact that the degradation in thin films continues within the unstressed interval. The phenomenon was represented by an analytical expression that reflects the fact that film degradation is more likely to occur in the presence of water molecules. The proposed model successfully predicted the increase in the diffusion constant of water molecules, and it also successfully reproduced the measured results. Therefore, the proposed model is useful for the purpose of extracting degradation parameters from measured results.

Appendix

On the Accumulation Energy per Volume around the SiO2 Molecule

When it is assumed that is the oxide layer thickness, is the electric field, is the stress time, is the lattice constant, is the stress current, and is the gate electrode area, we have

On the other hand, relates the total energy accumulated at the specific molecule (, , or ) as Combining (A.1) and (A.2), we have From (3), (9), and (11), we have

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The author expresses his thanks to Mr. A. Masuo, graduate student (presently, Panasonic Corp., Japan), for his technical support in experiments.