Abstract

We report the investigations of the ion migration polarization in the yttria stabilized zirconia (YSZ) thin films in the Metal-Oxide-Metal (MOM) and Metal-Oxide-Semiconductor (MOS) stacks due to the drift of the oxygen vacancies under the external bias voltage applied between the electrodes. The parameters characterizing the drift of the oxygen vacancies in YSZ such as the ion drift activation energy, mobile ion concentration, and the drift mobility have been determined in the temperature range 300–500 K. These data are important for deeper understanding of the fundamental mechanisms of the electroforming and resistive switching in the YSZ-based MOM and MOS stacks, which are promising for the Resistive Random Access Memory (RRAM) and other memristor device applications.

1. Introduction

The phenomenon of resistive switching in thin dielectric films has attracted much attention in last decade [1]. This phenomenon consists in a reversible persistent change of the resistance of thin (several tens nanometers thick) dielectric films sandwiched between two conductive electrodes under the external bias voltage applied between the electrodes. The electron devices utilizing the resistance switching effects are called memristors [2]. The potential fields of application of the memristors include the development of the nonvolatile computer memory of new generation (Resistive Random Access Memory, RRAM) [3, 4], neuromorphic electron devices [57] and so on.

Today’s understanding of the origin of the resistive switching in the metal oxides is based on the concept of the migration (the drift and/or diffusion) of the oxygen ions in the external electric field [8]. This process can be also described in terms of the migration of the charged oxygen vacancies. The oxygen vacancies in the electric fields between the electrodes of the memristor structure tend to arrange in the conductive filaments, which shortcut the electrodes (so-called forming process) [9]. As a result, the memristor switches from the initial high resistance state (HRS) to the low resistance state (LRS). The HRS can be restored by applying a voltage pulse of the opposite polarity destroying the filaments. Further operation of the memristor is provided by the restoring and rapture of the filaments resulting in switching the memristor from HRS to LRS and back (referred to as the SET process and the RESET one, resp.). It is obvious from the above that the oxygen ion transport plays a key role in the mechanism of resistive switching in the metal oxides. For deeper understanding of the resistive switching mechanism in these materials, it is important to know the parameters of the oxygen ion transport, such as the activation energies for the oxygen ion migration, the mobile ion concentration, and mobility. Knowing these parameters is essential to find the ways to improve the performance (first of all—the reproducibility and the reliability) of the memristor devices based on the metal oxides.

Yttria stabilized zirconia (YSZ) ZrO2(Y) is a promising material for the memristor devices [1013]. It is featured by a high concentration of oxygen vacancies that provides the high anion mobility [14]. It is clear from the above, that the availability of the oxygen vacancies is essential for the manifestation of resistive switching in the metal oxides since the oxygen ion transport in these materials goes through the jumps of the oxygen ions into the adjacent oxygen vacancy sites. Usually, the necessary oxygen vacancy concentration in the metal oxide films for the memristor devices is achieved either by deposition of the nonstoichiometric oxide films or by the postdeposition vacuum annealing. In contrary, in YSZ the necessary oxygen vacancy concentration is provided by doping by the stabilizing oxide Y2O3 [15]. This intrinsic property of YSZ opens a way to a precise control of the oxygen vacancy concentration by varying the molar fraction of Y2O3 (the vacancy concentration is 1/2 that of the Y atoms). To date, the prototype memristor devices based on the YSZ films have been fabricated and demonstrated a good durability (up to 107 SET/RESET cycles) [12, 13].

In the present study, we have measured the oxygen ion transport parameters in the thin YSZ (≈12% mol. Y2O3) films deposited by Magnetron Sputtering by studying the ion migration polarization [16] in the Metal-Oxide-Metal (MOM) and Metal-Oxide-Semiconductor (MOS) stacks based on these films. Several methods have been applied to identify the polarization mechanism and to determine the ion transport parameters: measuring the capacitance-voltage (C-V) curves, the high-frequency (HF) conductivity (G-V curves), the depolarization current [17], and the quasistatic capacitance of the MOS and MOM stacks. Besides, a novel method based on measuring the temperature dependence of the hysteresis in the cyclic current-voltage (I-V) curves has been proposed in the present study. This method is a variant of the classical one of measuring the dynamic I-V curves [18].

2. Materials and Methods

The MOM stacks were deposited onto the Si substrates oxidized and covered by the Ti (25 nm thick) and TiN (25 nm thick) metallization layers. The MOS stacks were deposited on the - and -Si(001) substrates with the specific resistivity of 4.5 Ω·cm. The YSZ films were deposited at the substrate temperature °C by Radio Frequency Magnetron Sputtering of the ZrO2 + 12% Y2O3 target in the Ar-O2 (50 : 50% mol.) gas mixture. The YSZ film thicknesses were 10, 15, and 40 nm. The Au gates of 40 nm in thickness with Zr or Ta sublayers were deposited onto the YSZ films at °C through a shadow mask. The values of the sublayer thickness were 3, 5, and 15 nm; the ones of the gate area were ≈10−3 and ≈10−2 cm2. To measure the temperature dependencies of the C-V, G-V, and the dynamic I-V curves, the samples were placed into a vacuum-pumped thermostat providing the variation of the sample temperature within the range 77–300 K. The measurements in the temperature range 300–510 K were performed in ambient air using either the isothermic regime or the rapid heating/cooling one. The uncertainty of the temperature measurements was ≈1°C.

The measurements of the C-V, G-V, and I-V curves were performed using Agilent® В1500А semiconductor device analyzer. The HF C-V and G-V curves were measured at the frequency  kHz with the ramp voltage sweep within the range −3–+3 V with a constant sweep rate . The dynamic I-V curves were measured at several values of ranging from 0.002 to 0.02 V/s. The ranges of the ramp voltage sweep were 0-1 V and 0–2 V. These measurements were performed in the isothermic regime at the temperature varied within the range 300–510 K.

In order to measure the depolarization current, the samples have been heated up to 510 K under the polarizing voltage  V and then were quenched down by the immersion into a Dewar flask filled with liquid nitrogen in order to freeze the polarized state. Then the samples were heated up with the rate ≈0.4°C/s at , and the depolarization current was measured. The quasistatic capacitance and its temperature dependencies were studied by measuring the discharge current of the MOS and MOM stacks after the charging by a constant voltage ranging from 0.1 to 2 V. It should be stressed here that all measurements were performed avoiding the forming and resistive switching; that is, the measured oxygen ion transport parameters were the intrinsic ones of the as-deposited YSZ films.

3. Results and Discussion

Typical HF C-V and G-V curves for the MOS stack with the YSZ insulating layer and Zr underlying one measured at various temperatures are presented in Figure 1. The values of for the HF G-V curves were normalized to the circular measurement frequency . All measurement results are given for the gate area  cm2; the positive values of the bias voltage correspond to the positive potential at the metal gates with respect to the substrates. An anomalous hysteresis was observed in the C-V and G-V curves at  K. This kind of hysteresis can be attributed to the ion migration in YSZ [19], whereas the normal hysteresis typical for the C-V and G-V curves of the MOS stacks originates from the exchange by electrons or holes between the semiconductor substrate and the surface states at the semiconductor/insulator interface [18].

The values of area within the hysteresis loops increased with increasing temperature in the temperature range 350–400 K. At higher temperatures, the hysteresis loop width tended to saturate and even decreased with increasing temperature. The C-V and G-V curves with such behavior of the temperature dependence of the hysteresis loop width were observed in the MOS stacks with relatively low density of surface states at the YSZ/Si interfaces. The values of the surface state density were obtained from the frequency and temperature dependencies of the admittance of the MOS stacks using the model of single level surface states [20]. The ionization energy  eV and the density of the surface states  cm−2 were found. The specific density of the surface states at ( is the flat band voltage) was estimated using the differential method [20] to be <5 · 1011 cm−2·eV−1.

In the samples with the Ta sublayers, the values of were higher than in the samples with the Zr ones (≈4 · 1012 cm−2·eV−1). The anomalous hysteresis in the C-V and G-V curves disappeared completely at higher temperatures (Figure 2). This result confirms the conclusion drawn above that the decrease of the anomalous hysteresis in the C-V and G-V curves of the YSZ-based MOS stacks with increasing temperature originated from the recharging of the surface states at the YSZ/Si interface.

Figure 3 presents the temperature dependencies of the flat band voltage shift at the flat band capacitance level in the anomalous hysteresis loops of the MOS stacks with different . For all MOS stacks with relatively low , the temperature intervals of growing with increasing temperature were observed. Within these intervals, the temperature dependencies of obeyed Arrhenius law with the activation energy  eV (Figure 3, curves (1) and (2)).

At higher temperatures, the saturation and even some decay of with increasing temperature was observed. Within this temperature interval, the highest concentration of the mobile ions responsible for the anomalous hysteresis was reached. The reduced surface density of the mobile ions could be estimated using the following formula [20]:where is the capacitance of the insulator determined from the saturated capacitance in the C-V curves (Figures 1 and 2) and is the elementary charge. As follows from (1), is proportional to . As the temperature dependence of obeys Arrhenius law in certain temperature range, the same should hold for . The origin of such temperature dependence of is the activation nature of the ion conductivity in YSZ as well as in other solid electrolytes [15].

The value of determined from curve (4) in Figure 3 was ≈2 · 1012 cm−2. Assuming the mobile ions to be distributed uniformly inside the YSZ film with the thickness  nm, the mobile ion concentration can be estimated as  cm−3. This value is an underestimate obviously since typical value of the oxygen vacancy concentration in single crystal YSZ is ~1021 cm−3 [21]. This disagreement could be attributed to a nonuniform distribution of the electric field inside the YSZ film. Indeed, once the oxygen vacancies are charged positively, the maximum electric field strength should take place at the YSZ/gate interfaces. It may seem reasonable to assume that the ion drift takes place primarily in these thin YSZ layers adjacent to the metal gates, and the remaining volume of the YSZ film does not contribute into the ion transport. However, this hypothesis cannot explain too low values of obtained for the MOS stacks. If it were true, the thickness of the layer containing the mobile ions could be estimated as ~10−4 nm for d = 40 nm. This value is much less than the lattice constant of YSZ (≈0.63 nm).

The underestimate of obtained from the anomalous hysteresis in the C-V curves of the MOS stacks was attributed to the simultaneous coexistence of the electron transport (the recharging of the surface states) along with the oxygen ion transport in YSZ. Both effects result in the hysteresis in the C-V and G-V curves. However, the contributions of these two effects into the resulting values of have the opposite signs, and these ones compensate each other. The effect of compensation is more pronounced at higher temperatures because the higher the temperature, the more the impact of the surface states on the hysteresis in the C-V and G-V curves. Also, the impact of the surface states should increase with increasing . This effect was observed in the experiment indeed (cf. Figures 1 and 2) which confirms the explanation given above.

In order to obtain more reliable data on the mobile oxygen ion concentration in the YSZ films under study, we performed the investigations of the ion migration polarization in the MOM stacks based on these films. It is well known that studying the ion transport in the dielectrics meets a number of problems often (the difficulty to find out the relation of the electric current in the external circuit to the ion transport inside the insulator, the electrochemical phenomena at the electrodes, etc. [22]). Therefore, we have applied a number of alternative measurement methods. Typical dynamic I-V curves of an YSZ-based MOM stack measured at various temperatures are presented in Figure 4. The dynamic I-V curves did not depend on the polarity of the bias voltage . The direction of the hysteresis loop was always the same (as shown in Figure 4).

The hysteresis observed in the dynamic I-V curves could be attributed either to the ion migration polarization of the insulator or to a trap one. After sweeping the bias voltage from 0 V up to its maximum value and back down to 0 V, a nonzero negative residual current was observed at . So far, the observed hysteresis in the dynamic I-V curves can be attributed to the electret effect [23, 24]. The magnitude of the effect increased with increasing temperature. It should be noted that the electret effect was most pronounced in the same temperature range, where the ion polarization-related anomalous hysteresis was manifested in the C-V and G-V curves of the YSZ-based MOS stacks (cf. Figures 1, 2, and 4). So far, the hysteresis observed in the dynamic I-V curves of the YSZ film-based MOM stacks can be related to the inertia of the oxygen ion drift when increasing and decreasing the bias voltage . Namely, the ion transport in the YSZ films can affect the electric field strength distribution inside these ones, which, in turn, results in the transient electric currents in the external electric circuit reflecting the ion drift inside the YSZ films. The hysteresis loop related to the ion migration polarization in the higher temperature range did not depend on the ramp voltage sweep rate (varied from 0.002 up to 0.02 V/s). So far, one can exclude the effect of the displacement current. The charging and discharging of the MOM capacitors cannot take effect as well since typical time constant of the current relaxation was <10−6 s as determined from the measurements of the small signal parameters of the equivalent scheme.

Since the magnitude of the hysteresis loop did not depend on the sign of the bias voltage applied to the metal electrodes at given measurement temperature, one can conclude the charges of the ions and traps in the investigated samples to be distributed uniformly within the YSZ film before the polarization [23]. The net charge can be determined from the area of the hysteresis loop using the following relation [18]:

Figure 5 presents the temperature dependencies of (curves (1) and (2)), of the forward current at zero reverse one (see Figure 4 for explanation) (curve (3)), and of the residual current (curve (4)). These temperature dependencies can be approximated by the combinations of two straight lines in the Arrhenius axes with the activation energies ≈0.55 eV (at  K) and ≈0.2 eV (at  K). The former activation energy can be related to the activation of the ion transport. This process was also manifested in the MOS stacks with relatively low density of surface states as well as in the MOM ones with the YSZ layers of various thickness (15 nm and 40 nm) regardless to the gate sublayer material (Zr or Та). The latter activation process can be related to the traps in YSZ. This effect was observed in the YSZ films deposited by Magnetron Sputtering only and was not observed in single crystal YSZ [25].

Assuming the uniform distribution of the charge in the YSZ films before the polarization, the surface density of the charge was estimated according to the formula

The values of the reduced surface concentration of the mobile ions in the YSZ films obtained from the dynamic I-V curves of the MOS stacks were not saturated with increasing temperature (unlike the ones obtained from the C-V curves of the MOS stacks) and can be estimated as ≈5 · 1014 cm−2 at  K. Correspondingly, the reduced surface density of the traps in the films of ≈40 nm in thickness was estimated to be ≈1.2 · 1012 cm−2. On the other hand, as it has been already mentioned above, the volume concentration of the oxygen vacancies in the single crystal YSZ with 12% mol. of Y2O3 is ~1021 cm−3 which corresponds to the surface vacancy concentration of ~1014 cm−2. So far, the oxygen vacancies in the investigated YSZ films were concentrated in a layer at least 5 times thinner than the YSZ film itself.

Also, the temperature dependencies of the quasistatic capacitance of the MOM stacks with the YSZ film thicknesses  nm (curve (6)) and ≈10 nm (curve (7)) are presented in Figure 5. These data support the estimates of the activation energies given above. Starting from  K corresponding to the beginning of manifestation of the ion migration polarization, the values of begin to grow with increasing temperature for all values of . At higher temperatures, this growth obeyed Arrhenius law with the activation energies ranging from 0.50 to 0.53 eV. These values are close to the ones for the ion migration obtained from the measurement of the C-V curves of the MOS stacks and of the dynamic I-V curves of the MOM ones. Assuming the dielectric permittivity of YSZ to be constant, one can estimate the thickness of the conductive layer in the YSZ film from the capacitance ratio at the maximum capacitance and temperature as ≈1.6 nm for  nm and ≈6.9 nm for  nm. At the same time, the quasistatic capacitance was not saturated with increasing temperature. This fact points to an incomplete ion polarization. The nonlinear electrical properties of the YSZ-based MOS and MOM stacks evidence the presence of the ions in the insulating layers, which become mobile in strong enough electric field (~105 V/cm) at higher temperatures and concentrated in a thin space charge layers adjacent to the metal gates.

Typical kinetics of and of after the polarization of the MOM stack with the YSZ film thickness  nm are presented in Figure 6. The sample was polarized at constant  V and  K and was cooled down to 300 K rapidly. Then the sample was heated up again at . The maximum absolute value of took place near  K. The analysis of the depolarization curve has shown the temperature dependence of to obey Arrhenius law when increasing . The activation energy was ≈0.55 eV. This value is close to the ones obtained from the C-V measurements of the MOS stacks as well as from the measurements of the dynamic I-V curves and of the quasistatic capacitance of the MOM stacks. One can determine the charge of the polarizing ions from the area under the depolarization curve:

The resulting values almost match the ones determined from the dynamic I-V curves. The activation energy of the ion conductivity for the oxygen vacancies in single crystal zirconia stabilized in the cubic phase is known to fall into the range 1.0–1.2 eV [14, 21]. As it has been shown in [14], this activation energy consists of two nearly equal components: (i) the energy of detachment of the ion from the lattice and (ii) the energy of its further motion over the potential barriers as it takes place in the classic mechanism of the ion migration polarization [23]. In the samples studied in the present work, the activation energies scattered from 0.5 to 0.6 eV. This could probably be related to the polycrystalline (nanocrystalline) structure of the investigated YSZ films so that the oxygen ions can move along the grain boundaries [26, 27]. Earlier, High Resolution Transmission Electron Microscopy studies demonstrated the YSZ films in the investigated stacks to have a nanocrystalline structure with the average grain sizes ≈40 nm (in the 40 nm thick YSZ films); the grain boundaries were perpendicular to the substrate [12]. The theory of the ion transport along the grain boundaries is absent to date. Probably, the activation energy of the ion transport along the grain boundaries is less than the one for the ion transport in single crystals [26].

One can estimate the ion mobility from the dynamic I-V curves presented in Figure 4. Assuming a uniform distribution of the ions in the YSZ films before the polarization, one getsThe film conductivity can be estimated from the slope of the linear parts of the reverse branches of the dynamic I-V curves shown in Figure 4 (when sweeping back from 1 V down to 0 V). This implies that if the polarization originates from the ion migration in the insulator, then the variation of the electric current in the external circuit equals the one of the ion current, that is, the oxygen vacancies arrive at the cathode and contribute to the electric current in the external circuit, but do not capture the electrons. In this case, one can estimate the ion mobility from the following relations for the ion conductivity :

Let us assume the concentration of the mobile ions  cm−3 to be constant, and the number of the polarized ions to grow with increasing temperature because of the increasing of the probability of the ion jumps over the potential barrier.

The temperature dependencies of and are presented in Figure 7. The values of increase with increasing temperature from  cm2/V·s at 425 K up to  cm2/V·s at 502 K. This result can be considered as an evidence of the polarization mechanism related to increasing the ion mobility with increasing temperature. In turn, this fact evidences the barrier mechanism of the ion transport in the YSZ films [14, 28]. The increase of the ion mobility with increasing temperature obeys Arrhenius law with the activation energy ranging from 0.50 to 0.55 eV corresponding to the height of the potential barrier between the ions. This value coincides with the ones obtained in the present study from the dynamic I-V curves, from the depolarization current measurements, and from the quasistatic capacitance of the MOM stacks as well as from the C-V measurements of the MOS stacks.

Also, one can estimate the drift mobility of the ions from the time-of-flight of the ions between the electrodes of the MOM stacks determined from the kinetics of charging due to the ion migration polarization when applying a bias voltage between the electrodes [17]. Assuming the absence of diffusion as well as of the electrochemical reactions at the electrodes (so that the ion transport is characterized by the ion drift mobility only) [29, 30], in the case of normal (Gaussian) transport one has That is, is proportional to the ratio where is the electric field strength in the insulator. This estimate gives the value of ion mobility  cm2/V·s at 500 K. This value is close to the one obtained from the dynamic I-V curves in spite of neglecting the phenomena at the electrodes as well as the electrostatic interaction between the ions, which may be essential at high ion concentrations in the investigated YSZ films (~1019–1020 cm−3). The parameters of the oxygen ion transport in the YSZ films of different thicknesses determined for various materials and thicknesses of the sublayers under the Au gates by various methods are summarized in Table 1.

As follows from Table 1, the parameters of the ion transport in the MOS and MOM stacks measured by various methods did not depend on the materials and thicknesses of the gate sublayers. Also, it follows from Table 1 that the C-V measurements of the MOS stacks gave too low values of the mobile ion concentration as compared to the ones obtained by other methods. As it has been already pointed out above, this discrepancy could be ascribed to the impact of the surface states at the YSZ/Si interfaces in the MOS stacks.

4. Conclusion

In the present work, we studied the ion migration polarization in the MOM and MOS stacks utilizing the YSZ films deposited by Magnetron Sputtering as the gate insulators. The values of the activation energy for the oxygen ion polarization in YSZ were measured using several different methods. The measurements of the dynamic I-V curves, of the depolarization current, and of the quasistatic capacitance of the MOM stacks gave close values of the activation energy (≈0.5 eV) regardless to the materials of the gates (Au/Zr, Au/Ta) and substrates (ТiN/Ti, ТiN) as well as to the thickness of the YSZ films (ranging from 10 to 40 nm). This result evidences the absence of the discharging of the ions at the electrodes and of the diffusion mechanism of the oxygen ion transport in YSZ. Also, the values of the mobile ion concentration and mobility in the YSZ films were determined for the temperature range 300–510 K. The data on the oxygen ion transport parameters in the YSZ films deposited by Magnetron Sputtering can be useful in modeling the processes of forming and resistive switching in the memristor devices based on the YSZ films, wherein the formation of the conductive filaments goes through the arrangement of the oxygen vacancies in the electric field between the electrodes.

Also, it was found that the C-V measurements of the MOS stacks with relatively high density of the surface states at the YSZ/Si interface gave the underestimated values of the mobile oxygen ion concentration. The origin of the underestimate was related to the impact of the surface states at the YSZ/Si interface. Namely, the anomalous hysteresis in the C-V and G-V curves of the MOS stacks related to the ion migration polarization of YSZ was compensated by the normal hysteresis originating from to the trapping of the electrons or holes onto the surface states.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors gratefully acknowledge the financial support from Ministry of Education and Science, Russian Federation (Project no. RFMEFI58717X0042).