3. Discrete Electronic Perturbation Is Unstable
To carefully investigate a discrete electronic perturbation, a reliable oxide film is convenient for excluding non-targeted fluctuation factors, as far as possible. In [
8], we prepared such a convenient sample using a 130 nm standard complementary metal-oxide-semiconductor (CMOS) process.
Figure 3 contains an illustration and transmission electron microscope (TEM) images describing the sample geometry. The left-hand side (a) is a cross-sectional view along the cigar-shaped polysilicon above the substrate. The substrate surface is divided into two parts, labeled “S/D” and “G” by using the shallow trench isolation (STI). C1 is the capacitance between the cigar-shaped polysilicon and the substrate surface G. C2 is the capacitance between the cigar-shaped polysilicon and the substrate surface S/D. If the width of S/D is larger than that of G, then C2 is greater than C1, then the following capacitance coupling ratio (
) becomes greater than 0.5.
The center (b) is a bird’s eye view of this sample. Along the perpendicular to (a), we have the source and drain diffusion layers, labeled “S” and “D”, respectively, on the substrate surface S/D. These diffusion layers sandwich the channel. We have two diffusion layers labeled “G” on the substrate surface G. The diffusion layers, G, can connect to each other by the back-end of line (BEOL). The gate lengths on both substrate surfaces are the same as the width of the cigar-shaped polysilicon (130 nm). Thus, we can easily design by tuning the gate width of both substrate surfaces (S/D and G). In this geometry, the gate oxide on the substrate surface G can serve as a tunnel oxide, while . If a positive voltage is applied to the diffusion layers S/D, then the vertical electric field is concentrated above the substrate surface G. When the voltage is large enough, an inversion layer is generated between the two diffusion layers G. Electrons coming from the diffusion layers G to the inversion layer can be injected into the cigar-shaped polysilicon by Fowler–Nordheim tunneling (FNT).
The righthand side image (c) is a TEM image of a cross-section in the plane perpendicular to the cigar-shaped polysilicon. It is the same as the cross-sectional view of the metal-oxide-semiconductor field effect transistor (MOSFET) of the 130 nm standard CMOS process, except for having no contact with the cigar-shaped polysilicon. Since the cigar-shaped polysilicon can in this way serve as the floating gate (FG), this sample is useful for various experiments of embedded type non-volatile memories (e.g., a battery-less timer [
8]). Since the standard CMOS process is used to fabricate it, the quality of oxide and the substrate interface is within the scope of mass-production. Furthermore, we can obtain a sufficient number of sample cells in a mass-production line in which the cell-to-cell variation is under control.
We investigated the current flow of electrons from the diffusion layers, G to FG. The measurement scheme is as follows. We chose a sample cell with the gate lengths of 0.52 µm and 0.4 µm in the substrate surface S/D and the substrate surface G, respectively. This means that
. We applied a voltage to the source and drain diffusion layers S/D, while the substrate and the diffusion layers G were grounded. The substrate voltage was actually applied to the well contact, which is far from the channel between the diffusion layers S/D. Thus, we observed that
approximates the FG potential to the inversion layer being generated between the diffusion layers G. While sweeping the voltage from zero to a voltage at which an abrupt increase appears (i.e., stopped by the ampere limiter), we measured the current flowing to the diffusion layers G by subtracting the substrate current from it. This procedure was repeated multiple times in succession for each cell. However, between any two sequential measurements of each cell, the voltage became zero once. The breakdown voltage of a 130 nm node is about 7.2 V [
9]. The repeated applications of nearly
to the tunnel oxide may generate an electronic perturbation. The oxide thickness is 3 nm or less in a 130 nm node [
9]. This is a little bit thinner than necessary to generate the FNT at 10 MV/cm, which could cause the FNT in flash memories. The FNT is not caused by damage to the oxides. The tunnel gate area was about 0.052 µm
2 (=400 nm
130 nm). The number of repetitions of voltage application was 30 for the measurement of each cell.
In
Figure 4a, we show the results of the first measurement of a chosen cell (along the vertical axis) over the electric field until 11.3 MV/cm (along the horizontal axis). The signal of the measured current was lower than the measurement limit. This limit was due to the equipment noise, and not of interest in this work. The leakage, if it exists, can be regarded as direct tunneling with no trap-assistance, that is, tunneling not assisted by traps (i.e., direct tunneling), as illustrated in
Figure 4b. Since we were interested in a discrete electronic perturbation in the oxide, we further continued the sequential measurement of this cell.
In
Figure 5a, we show the results of the 2nd to 18th measurements of the same cell as the first measurement. The vertical axis is for current and the horizontal is for electric field. A moderate increase of the current occurs from 5 MV/cm. This moderate increase may be attributable to 1-trap process in a local trap having been generated in the previous measurement (stressing by more than 10 MV/cm), as illustrated in
Figure 5b. This is called the stress-induced leakage current (SILC) [
10]. The leakage current is sensitive to the location and the energy level of the traps. Therefore, if the location and level of traps vary over the sequential measurement of the same cell, then the leakage current may vary in proportion to the fluctuation of the trap location and level. In (c), we consider the appearance of a second trap. If the location of the second trap is in a different tunnel path from that of the first trap, then it can increase the leakage current by about two-fold. Thus, the current of the 2nd to 18th measurements makes a bundle, as in
Figure 5a. We, furthermore, continued the sequential measurement of the same cell.
In
Figure 6a, we show the results of the 19th and 20th measurements of the same cell (along the vertical axis) over the electric field (along the horizontal axis). There was a discontinuous increase from the moderate increase by ten-fold in the current level in the 19th and 20th measurements of the same cell from 2 MV/cm. Ielmini et al. found that it is hard to explain this by using two traps which are in different tunnel paths [
11]. Following them, we can regard this phenomenon as related to the two-trap process (i.e., two traps existing in a tunnel path [
11]), as illustrated in
Figure 6b. Accordingly, a second trap may be generated due to the stressing in the 18th measurement of the same cell. This implies that a local trap is a discrete perturbation. In (c), we can consider that there is another trap in a different tunnel path. However, the two-trap process may dominate it. We, furthermore, continued the sequential measurement of the same cell.
In
Figure 7a, we show the results of the 19th and 20th measurements of the same cell along the vertical axis over the electric field along the horizontal axis. What is surprising is that the current became lower than the measuring limit in the 21st to 30th measurements of the same cell. It is possible that the traps (i.e., discrete electronic perturbations) became inactivated or vanished. Therefore, we can regard the current in the 21st to 30th measurements as tunneling not assisted by traps (i.e., returning back to #1), as illustrated in
Figure 7b.
6. Power Spectrum Analysis of Discrete Electronic Perturbations
As we can understand from the discussion of
Figure 8 and
Figure 9, the emission from the cathode or the IL is the key aspect that we have to analyze. The emission rate must be dependent of the electric field across the oxides, however, in this chapter, we a-priori assume that the emission rate is not explicitly dependent on time. The time dependence of the leakage current is assumed to be attributable only to the time-varying charge in the IL. In other words, we assume that the emission process is a stationary process. The autocorrelation of transient leakage current,
, can be thus regarded as independent of time (
), where
is the lag. We can depict it as
. We apply the Wiener–Khinchin theorem (6a) to analyze the power spectrum of the measured transient current,
, where
is the frequency,
depicts the expectation value of
[
18,
19,
20]. In (6b),
is proportional to
, where
is the spectral index, and ranges from 0.8 to 1.45 [
21,
22]. It was observed that the spectral index is sensitive to the location of electronic perturbations; at the interface or inside oxides [
21,
23].
When
, it is the
-fluctuation, i.e., Flicker noise [
24,
25,
26,
27,
28]. When
, it is the
-fluctuation, i.e., Brownian noise [
29]. In general, the transient current varies over the samples to be measured. Therefore, it is indispensable to exclude the sample variations from the measured transient leakage current [
30]. In
Figure 10a, there is an illustration of the measured sample of a metal/high-K stack/metal (MIM) capacitor. The top and bottom electrodes are composed of metal (TiN). The top electrode is applied with a gate voltage (
) to measure the current, while the bottom electrode is grounded. The high-K oxide stack is composed of a 5 nm TiO
2 layer and 10 nm ZrO
2 layer from the bottom. There may be an interface between the 10 nm ZrO
2 and the 5 nm TiO
2, which corresponds to the IL in
Figure 8 and
Figure 9. The gate area is 2500 µm
2. The open probe fluctuation is about
A, which is equivalent to the measuring limit in
Figure 4,
Figure 5,
Figure 6 and
Figure 7. Both the 10 nm ZrO
2 and the 5 nm TiO
2 are polycrystalline, as shown in the TEM image of
Figure 10b. It was found that ZrO
2 is cubic or tetragonal and TiO
2 is rutile (tetragonal). Some of the grains inside the polycrystalline may face the IL. The distributions of grain orientations, grain sizes, and grain boundary geometry inside the oxide stack may cause the sample variations to be an obstacle for the power spectrum analysis.
In
Figure 11a,b, we plotted the transient current measured using nine samples (labeled #1–#9) having the same geometry as
Figure 11a. The vertical axis is the measured current density and the horizontal is time. In (a), the transient current of the samples was measured under a positive fixed electric field (+0.5 MV/cm). In (b), the transient current of the samples was measured under a negative fixed electric field (−0.5 MV/cm). First of all, the plot shows that the sample-to-sample variation (named, the sample variation) was quite large. Next, we can see an uneven fluctuation (discrete electronic perturbations) in the data of the sample #1 under −0.5 MV/cm and so forth.
To analyze the electronic perturbation appropriately, we have to exclude the sample variation. In
Figure 12, we illustrate the band diagram of
Figure 10 with a moderate positive
applied, which corresponds to +0.5 MV/cm. The horizontal axis is the
z-axis (depth direction) and the vertical axis is the conduction band edge of the stacking film sandwiched by the electrodes. From the bottom, there are three interface layers (IL-1, IL-2, and IL-3) at the interfaces between the bottom electrode (TiN) and the 5 nm TiO
2, between the 5 nm TiO
2 and the 10 nm ZrO
2, and between the 10 nm ZrO
2 and the top electrode (TiN), respectively. In this figure, we assumed no electrons were stored in the interface layers (i.e., the IL-1, the IL-2, and the IL-3). The barrier height of the 5 nm TiO
2 is 0.35 eV from the bottom electrode [
31,
32,
33,
34]. The barrier height of the 10 nm ZrO
2 from the top electrode is 0.94 eV [
31,
32,
33,
34]. In general, the dielectric constant decreases with the band gap (or tunnel barrier). Therefore, the electric field across the 10 nm ZrO
2 is higher than the 5 nm TiO
2, while we can ignore the stored electrons in the interface layers.
In a usual analysis, the IL-1 and IL-3 may be ignored by supposing that the electrons come from an electrode to the other electrode passing through homogenous oxides (having no sample variations). However, both the oxide stack and the interface layers must have sample variations due to grain-related distributions. In other words, for +0.5 MV/cm, an electron can be emitted from the IL-1 to the IL-2 through the 5 nm TiO
2 with sample variations and then be emitted to the IL-3 through the 10 nm ZrO
2 with sample variations. Due to the higher barrier and the thicker width, the 10 nm ZrO
2 dominates the serial resistance of this oxide stack. Therefore, for a positive electric field, to analyze the transient current, we must consider the emission of electrons from the IL-2 (i.e., the cathode interface) to the IL-3 (i.e., the anode interface). For a negative electric field, we must consider the emission of electrons from the IL-3 (i.e., the cathode interface) to the IL-2 (i.e., the anode interface). That is, the IL-2 and IL-3 are the cathode interfaces for a positive and negative electric field, respectively. The IL-2 and IL-3 are the anode interfaces for a negative and positive electric field, respectively. The DT may likely occur between the IL-2 and the IL-3 due to the high tunnel barrier under a moderate electric field. Thus, we have the sample variations (see
Figure 11) that may be attributable to the distribution of grains inside the 10 nm ZrO
2 sandwiched by the IL-2 and IL-3. If we have a sufficient number of stored electrons in the IL-2, the band diagram is modulated to
Figure 13. This corresponds to
Figure 9b.
The 10 nm ZrO
2 is a polycrystalline having variations of grain orientation, grain sizes, and the grain boundary patterns inside. The IL-2 and IL-3 may also be affected by the orientations, the sizes, and the boundary patterns of the grains facing the interfaces. This may be the source of the sample variations. In [
30], we regulated these complex sources of sample variation and then obtained a formula to analyze the transient DT current density, as follows: where
is the stationary DT current density between the IL-2 and the IL-3.
We, a-priori, ignore the influence of the anode interface on the sample variations and then consider the grain-related variations in the cathode interface and the 10 nm ZrO2. However, is the average of surface roughness, which varies over grains facing the cathode interface, and is the average of the longest dwell time of electrons which are emitted from the cathode interface over the grains. The is always positive and less than unity. As is close to unity, the surface roughness reduces at the cathode interface.
Subsequently, suppose there is a trapping site inside the 10 nm ZrO
2, as illustrated in
Figure 14. The horizontal axis is the
z-axis (depth direction) and the vertical axis is the conduction band edge of the stacking film sandwiched by the electrodes. If the energy levels of the trapping sites are low enough, then the TAT may likely occur from the cathode interface (IL-2) to the anode interface (IL-3) via the trapping site. This may dominate the direct tunneling between the IL-2 and the IL-3. In a similar way to the derivation of (7), we obtained the formula to analyze the transient TAT current density, as follows: where
is the stationary TAT current density,
is the dwell time of electrons emitted from the trapping sites for the TAT process, and
is a positive number.
However, the trapping sites may also vary in depth in the energy diagram between samples. In one case, the electrons are stuck in the trapping sites. This cannot contribute to the measured leakage current. If is smaller than unity, the TAT dominates. Otherwise, the electron is stuck as a fixed charge there.
The parameters (
,
,
,
,
,
) characterize the sample variations and not the discrete electronic perturbations. By tuning these parameters, we obtained an excellent agreement with the measured data, without uneven fluctuations (i.e., discrete electronic perturbations). An example is sample No. 6 (#6), shown in
Figure 15. In (a), the vertical axis is the current density and the horizontal one is time. The measured data under the electric fields (0.5 MV/cm, 0.7 MV/cm, and 1 MV/cm) are depicted by black open circles. The lowest current can be excellently reproduced only by using (7) (i.e., DT only). The others can be excellently reproduced by using both (7) and (8) (i.e., both DT and TAT). The TAT can become involved as the electric field increases. In (b), the vertical axis is for the power spectrum densities and the horizontal one is for frequency. The black, red and blue open circles depict the power spectrum densities converted from the measured transient current densities at +0.5 MV/cm, +0.7 MV/cm, and +1.0 MV/cm using (6a), respectively. The calibrated
is 0.224, 0.321, and 0.408 at +0.5 MV/cm, +0.7 MV/cm, and +1.0 MV/cm using (6b), respectively. It increases as TAT becomes involved. The calibrated parameters are shown in
Table 1. In
Figure 16a, the measured transient current of samples No. 2 (#2), No. 6 (#6), and No. 8 (#8) under +0.5 MV/cm are depicted by black open circles. The vertical axis is the current density and the horizontal one is time. We obtained an excellent agreement by using only (7) (i.e., DT). In (b), the vertical axis is the power spectrum densities and the horizontal one is frequency. The red, black, and blue open circles depict the power spectrum densities converted from the measured transient current densities of samples #2, #6, and #8 at +0.5 MV/cm using (6a), respectively. The black open circles are the same as in
Figure 15b. In addition, we observed that
and
for #2 and #8 at +0.5 MV/cm using (6b), respectively. The calibrated parameters are also shown in
Table 1.
In
Figure 17, the vertical axis is
and the horizontal is the surface roughness index (
). The open dots depict those calibrated using sample #6 at +0.5 MV/cm, +0.7 MV/cm, and +1.0 MV/cm in
Figure 16a,b, respectively. Since these transient currents can only be fit using (7) (i.e., DT only), the number of fitting parameters is three (
,
, and
). See
Table 1. Let us choose
pA/cm
2 and
s from the sample #6 at +0.5 MV/cm in
Table 1 (fit by DT). Regarding
as a variable, we can calculate the transient DT current using (7). Substituting it for (6a), we can calibrate
using (6b) for a given
, as plotted using the line in
Figure 17. It is thus found that
increases as
decreases from unity (i.e., no surface roughness). We can thus find that the surface roughness violates the power law of the Flicker noise (
). See
Table 1 again. In sample #2 at +0.5 MV/cm, the transient current was reproduced using only (7) (i.e., only DT) and
, which is very high. However,
, which is very small compared with the others. This means that #2 has a larger surface roughness. According to the trend in
Figure 17, we can assume that
would reduce from 0.813, as the surface roughness could have been suppressed in sample #2. If we can reduce the surface roughness, then
may also have decreased in samples #6 and #8 at +0.5 MV/cm from 0.224 and 0.108, respectively. Thus, we can observe that
(i.e., the spectral index is 1) when the electrons are directly emitted from the cathode to the anode interfaces (i.e., DT). The DT turns out the Flicker noise and the surface roughness violates the power of the Flicker noise. In sample #6, furthermore,
increases due to the TAT as the electric field increases. The TAT also violates the power law of the Flicker noise.
As mentioned above, we can observe a fluctuation which may be discrete and subject to neither (7) nor (8). For example, there is an uneven transient current for sample No. 1 (#1) under −0.5 MV/cm in
Figure 11b. This measured data is replotted using lines and dots in
Figure 18a. The vertical axis is the current density and the horizontal one is time. There are three discrete perturbations. In (b) of this figure, by tuning the fitting parameters of (7) and (8), we can obtain the continuous transient current depicted by the red dash line, with the calibrated fitting parameters shown in
Table 2. By subtracting it from the measured transient current (depicted dot and line), we obtained the red line. It appears as the summation of a stationary current (
A) and the uneven component of the fluctuating current with no continuous decay. The calibration of the parameters was carried out to exclude the continuous decay from the red line. That is, we can deduce the uneven component from the measured data. We considered the uneven components as random-telegraph noise (RTN) in the transient leakage current [
35]. Machlup assumed bi-states (i.e., discrete) and used (6a) and (6b) to study RTN [
18]. Its origin has been considered to be related to number fluctuation [
24], to mobility fluctuation [
36,
37,
38], to both number and mobility fluctuations [
39,
40,
41], and phonon scattering [
42].
Let us remember that multiple grains of the 10 nm ZrO
2 may face either of the cathode or anode interfaces. Accordingly, there may be grain boundary patterns at the interfaces. Some traps are on the grain boundaries and others are out of the grain boundaries at the interfaces. We can suppose that the details of trap states are more complicated on the grain boundaries. This may cause the barrier height fluctuations (
and
) in the cathode and anode interfaces, respectively. The electrons captured by the trap sites on the grain boundaries may change the trap states by emitting and absorbing the phonon energy (
). For the phonon absorption,
is positive. For the phonon emission,
is negative. Thus, the dwell time of electrons stored in the cathode interface may fluctuate, as follows: where the suffixes of
depict the positive and negative electric fields, respectively, and
is the dwell time without the uneven component.
The denominators of the terms in the exponent of (9) are written as follows: where
is the effective tunnel mass and assumed to be half the rest of the electron mass,
is the thickness of the insulator (10 nm in this case),
and
are the conduction band edges in the cathode and anode interfaces, respectively, and
is the compensation of
.
If
, then
and vice versa. In the bottom plot of
Figure 18c, by using these parameters related to phonon scattering, we can extract
in the vertical axis over time in the horizontal axis. The upper is the replot of the uneven component depicted by the red line in (b). However, for the calibration of parameters in (9) and (10), we synchronized the occurrence of phonon scattering in our model and that of discrete perturbations in the extracted uneven component (i.e., in the measurement). In (d) of this figure, we obtained an excellent agreement with the measurement by using (7)–(10) with the phonon scattering (i.e., RTN), as depicted by the blue line. The calibrated parameters for phonon scattering are shown in
Table 3 (see #1 therein). In (e), furthermore, we can convert this measurement to the power spectrum density in the frequency domain using (6a). We can thus extract
using (6b).
We also observed an uneven transient current with samples #3 at −0.5 MV/cm, #4 at +1.0 MV/cm, #7 at +0.7 MV/cm, and #8 at −0.7 MV/cm. In a similar manner, using (7)–(10), we obtained an excellent agreement for the measured transient current densities, as shown in
Figure 19. The vertical axis is current density and the horizontal one is time. In (a), we plot the fitting results of #1 at −0.5 MV/cm (replot of
Figure 18d), #3 at −0.5 MV/cm, and #4 at +1.0 MV/cm. In (b), we plot the fitting result of #7 at +0.7 MV/cm and #8 at −0.7 MV/cm. The calibrated parameters are shown in
Table 2 and
Table 3. Subsequently, using the same method to obtain
Figure 18c, we extracted phonon energies using these samples, as shown in
Figure 20. In the upper plots, we show the discrete perturbation (i.e., uneven component) over time. In the bottom plot, we show the extracted
over time. However, for the calibration of parameters in (9) and (10), we synchronized the occurrence of discrete electronic perturbations and phonon scattering. The
extracted using the same method to obtain
Figure 18e was around 1.0, as shown in
Table 2. Thus, we can find that phonon scattering generates randomness due to the
-fluctuation, that is, Brownian noise.
By comparing with the theoretical infrared phonon scattering [
43], we can suppose that if
then grains are cubic or tetragonal, which is consistent with the 10 nm ZrO
2 of the samples used in the present measurement. See
Figure 10b. Otherwise, they are monoclinic grains. In this way, we can deduce some physical intuitions from the phonon scattering analysis. By repeating the measurement of the same samples with the same conditions, sometimes we could obtain such uneven fluctuations. On other occasions, we could not obtain it since the discrete electronic perturbations were unstable, as we have shown in
Figure 4,
Figure 5,
Figure 6 and
Figure 7.
However, it may be still an open problem how phonons scatter with a discrete local trap. Grasser et al. reported a time-dependent analysis of metastable defect states [
44] and Ielmini et al. performed a structure relaxation due to trapped holes [
45]. It would be interesting to know how or if such details are related to phonon scattering, as we could extract some more detailed physical intuitions from phonon scattering energies.
Nowadays, there are many kinds of semiconductor memory. As conventional ones with MIM capacitors, we have static random-access memory (SRAM), dynamic random-access memory (DRAM), NOR flash, and NAND flash. As emerging types, we have phase-change random access memory (PCRAM), magnetic random access memory (MRAM), resistivity random access memory (RRAM), and so forth. The analysis of spectral index will be critically important for investigating various kinds of reliability issues, not only in the conventional ones, but also in the emerging ones. For example, a complex RTN was reported in RRAM [
23,
46].