# Dielectric Relaxation of Lanthanide-Based Ternary Oxides: Physical and Mathematical Models

Authors: Chun Zhao; C. Z. Zhao; Jing Tao; M. Werner; S. Taylor; P. R. Chalker
Journal: Journal of Nanomaterials
 (2012)

---

## Abstract

Cerium-doped hafnium oxides (CexHf1−xO2) and lanthanum-doped zirconium oxides (LaxZr1−xO2) were investigated. The highest dielectric constants, k, were obtained from lightly doped oxides with an La content of x=0.09 and a Ce content of x=0.1, for which k-values of 33~40 were obtained. The dielectric relaxation appears to be related to the size of crystal grains formed during annealing, which was dependent on the doping level. The physical and mathematical models were used to analyze the relationship between k-values and frequencies. The variations in the k-values up to megahertz frequencies for both CexHf1−xO2 and LaxZr1−xO2 are simulated based on the Curie-von Schweidler (CS) or Havriliak-Negami (HN) relationships. Concerning the lightly doped CexHf1−xO2 and LaxZr1−xO2, the data extracted are best modeled by the HN law, while LaxZr1−xO2 with doping level from x=0.22 to 0.63 are best modelled based on the CS law.


---

## Body

## 1. Introduction

As the thickness of SiON gate dielectric thin films used in complementary metal oxide semiconductor (CMOS) devices is reduced toward 1 nm, the gate leakage current level becomes unacceptable [1]. To overcome this leakage problem, extensive efforts have been focused on finding alternative gate dielectrics for the 45 nm node and beyond technologies [2]. ZrO2 and HfO2 are two of the most promising high dielectric constant (high-k) materials for the replacement of SiO2 in MOSFETs, but it has been found that crystallization occurs at ~500°C for pure hafnia [3]. Doping with lanthanum (La) or cerium (Ce) increases the crystallization temperature [4]. Doping hafnia and zirconia thin films with rare earth elements can also stabilize the metastable tetragonal or cubic phase following annealing which enhances the dielectric constant [5]. Since dielectric relaxation and the associated losses can impair MOSFET performance, the larger dielectric relaxation of most high-k dielectrics compared to SiO2 is a significant issue for their use [6].Capacitance-voltage (C-V) measurements are a fundamental characterization technique for MOS devices. C-V measurements have been widely used to extract the dielectric constant (k-value) and dielectric loss of high-k materials using the capacitance measured in strong accumulation. Several models and analytical formulae have been thoroughly investigated for correcting the data from frequency dispersion [7–11]. When the above effects were taken into account, frequency dispersion in the accumulation region of C-V curves is still observed from high-k dielectric thin films due to the decrease of k-value with frequency (f), which is identified as a dielectric relaxation effect [12, 13].However, there is insufficient information about dielectric relaxation of high-k thin films, which prompts us to investigate the phenomenon and the underlying mechanism. In this paper, the dielectric relaxation and frequency dependence of the k-value of CexHf1-xO2 and LaxZr1-xO2 dielectrics were examined. The physical and mathematical models used were according to recent research for fitting the measured k-f curves. Variations in the k-value up to megahertz frequencies are simulated based on the Curie-von Schweidler (CS) or Havriliak-Negami (HN) relationships.


## 2. Device Preparation and Measurements

AmorphousLaxZr1-xO2 thin films (x=0.09, 0.22, 0.35, 0.63) and CexHf1-xO2 thin films (x=0.1) were deposited on n-type Si(100) substrates using liquid injection atomic layer deposition (ALD), carried out on an Aixtron AIX 200FE AVD reactor fitted with the “Trijet” liquid injector system. The doping level was varied up to a concentration level of 63%, that is, x=0.63. The interfacial layer between high-k thin film and silicon substrate is a ~1 nm native SiO2 determined by cross-section transmission electron microscopy (XTEM). Samples were then annealed at 900°C  for 15 minutes in an N2 ambient to crystallize the thin films.All the MOS capacitors were fabricated by thermal evaporation of Au gates through a shadow mask with an effective area of4.9×10-4cm2. The backside contact of the Si wafers was cleaned with a buffer HF solution and subsequently a 200 nm thickness of Al was deposited by thermal evaporation to minimize the effect of series resistance. The physical properties of the thin films were studied using auger electron spectroscopy, X-ray diffraction, medium energy ion scattering, XTEM, and atomic force microscopy. Their electrical properties were investigated by current-voltage (I-V), high-frequency capacitance-voltage (C-V), and capacitance-frequency (C-f) measurements.


## 3. Experimental Results and Discussion

TypicalC-Vresults of SiO2, HfO2, and Ce0.1Hf0.9O2 dielectrics are shown in Figure 1 and Figure 2, respectively. The C-V results shown in Figure 1 are taken from a reference SiO2 sample (Figure 1(a)) and from an HfO2 sample (Figure 1(b)). No significant variation in C-V response was observed in the reference SiO2 and HfO2 dielectrics, which confirms that the C-V measurement system does not introduce any frequency dispersion. Furthermore, Figure 1(b) shows that the interfacial layer (~1 nm native SiO2) between the HfO2 and the Si substrate also introduces negligible frequency dispersion for the samples used here [9, 11].C-V measurements from SiO2 (a) and HfO2 (b) at different frequencies. No  significant frequency response was observed in the reference SiO2 and HfO2 dielectrics.
(a)(b)Figure 2
C-V measurements from a CexHf1-xO2 thin film at different frequencies (from 100 Hz to 200 kHz).  Frequency dispersion was obtained regardless of the interfacial layer of MOS structures and measurement system, which was identified due to the frequency dependence of k-value (termed dielectric relaxation).However, for the various frequencies ranging from 100 Hz to 200 kHz applied on the annealedx=0.1 sample of CexHf1-xO2 thin films, different C-V results are exhibited, especially obtained in strong accumulation (unlike the SiO2 and HfO2 in Figure 1). Figure 2 shows that substantial frequency dispersion was observed during C-V measurements on the CexHf1-xO2 sample. The results of Figure 2 demonstrate that frequency dispersion can occur occasionally regardless of the interfacial layer of MOS structures and measurement system. In this case, the frequency dependence of the k-value is due to dielectric relaxation. The dielectric relaxation is modelled in more detail in follow, using the CS and Kohlrausch-Williams-Watts (KWW) or HN law to accurately simulate the measured phenomena.The dielectric relaxation in the time domain can be described by a power-law time dependence (CS law),t-n, or a stretched exponential time dependence (KWW law), exp[-(t/t0)m], where n and m are parameters ranging between 0 and 1 and t0 is a characteristic relaxation time [14–17]. There exist two alternative physical approaches to the interpretation of dielectric relaxation: the parallel and series models [18].The parallel model represents the classical relaxation of a large assembly of individual relaxing entities such as dipoles, each of which relaxes with an exponential probability but has a different relaxation timetk [17, 18]. The total relaxation process corresponds to a summation over the available modes k, given a frequency domain response function. Parallel processes are possible among independent noninteracting entities. The parallel model may be approximated by the HN relationship. The alternative approach is the series model, whereby each successive stage of relaxation is conditioned by the immediate past history. The series model may be used to describe briefly the origins of the CS law (the t-n behavior).Dielectric relaxation can be modeled by either the CS law or the HN law (KWW law) for high-k thin films. However, in some complex condensed systems, neither the pure parallel nor the pure series approach is accepted and instead observed results interpolate smoothly between these extremes [19]. The CS behavior is verified to be faster than the HN function at short times and slower than the HN function at long times.The general type of dielectric relaxation can be described from the CS law and the KWW law which are formulated as follows:(1)dPCSdt∝t-nwith0≤n≤1,PKWW∝exp⁡[-(tτ)β0]with0≤β0≤1,
where PCS is the CS polarization and the exponent n indicates the degree of dielectric relaxation. PKWW is the KWW polarization, τ is the KWW relaxation time, and β0 is the PKWW parameter. The complex susceptibilities,χCS and χKWW (frequency domain), related to the CS law and KWW law, are simply the Fourier transforms of the above time-domain responses [20]: (2)χCS(ω)∝F(dPCSdt)∝∫0∞t-nexp⁡(-iωt)dt∝(iω)n-1
or    (3)χCS=A(iω)n-1,
wherenandA are the relaxation parameters, and(4)χKWW(ω)∝F(dPKWWdt)∝∫0∞β0τ(tτ)β0-1×exp⁡[-(tτ)β0]exp⁡(-iωt)dt.
In the frequency domain, after a Fourier transform, the corresponding dielectric response (3) can be approximated by a simple HN relationship:(5)χHN(ω)=(εs-ε∞)(1+(iωτ)1-α)β1,
where εs and ε∞ are the static and high-frequency limit permittivities, respectively, ω=2πf is the angular frequency, the τ is the HN relaxation time, and a and β1 are the HN relaxation parameters. The parameter α is a distribution parameter, which accounts for the decrease of the peak loss and broadening of the distribution. This was also used in the Cole-Cole equation [18].The reason for replacing the KWW law using the HN relationship is to avoid the complicated numerical calculation implicit in (4). It was reported by Bokov and Ye that any Fourier transform of the KWW function into the frequency domain can be approximated by an HN function, but not vice versa [20].Then the intrinsic effect of frequency dispersion caused by the dielectric relaxation can be described by the CS law (3) and/or the HN relationship (5) [21, 22]. Separating the real and imaginary parts of the CS equation (3) gives the following equation for ε′ and ε′′ [23]: (6)ε′(ω)=ε∞+Bωn-1,ε′′(ω)=Cωn-1,
where B and C are constants. The real and imaginary parts of HN equation (5) may be separated which gives the following equation for ε′and ε′′[24]: (7)ε′(ω)=ε∞+(εs-ε∞)×cos⁡β1φ[1+2(ωτo)1-αsin⁡⁡(1/2)απ+(ωτo)2(1-α)]β1/2,(8)ε′′(ω)=ε∞+(εs-ε∞)×sin⁡β1φ[1+2(ωτo)1-αsin⁡⁡(1/2)απ+(ωτo)2(1-α)]β1/2,
where(9)φ=arctg[(ωτo)1-αcos⁡⁡(1/2)πα1+2(ωτo)1-αsin⁡⁡(1/2)απ].
The measured relationship between the k-value and frequency can be extracted from Figure 2, as shown in Figure 3 (cross-symbols). The experimentally observed dielectric relaxation for CexHf1-xO2 was modeled using the HN relationship (7), as shown in Figure 3 (solid line); that is, the cross-symbols in Figure 3 are extracted from the measured data and the solid line is from the model using the HN relationship  (7). For the sample used (CexHf1-xO2, the composition of Ce is x=0.10), the relaxation parameter α,β1, and the relaxation time τ are 0.01, 0.0535, and 0.102 s, respectively. The fitted result, indicated by the solid line, is very close to the measured data.Figure 3
Dielectric relaxation ofCexHf1-xO2 observed in C-Vmeasurements (Figure 2) can be fitted by the HN relationship. The cross-symbols are the measured data and the solid line is from HN relationship given by  (7).For comparison to theCexHf1-xO2 results shown in Figure 3, thek-f data from the LaxZr1-xO2 thin films are shown in Figure 4. The square symbols are the measured data from La0.09Zr0.91O2 where the solid line is modeled from the HN law. The diamond symbols, triangle symbols, and circle symbols are experimental data from the LaxZr1-xO2 thin films with doping levels of 22%, 35%, and 63%, where the solid lines are all simulated using the CS law. The parameters (α, τ,andn) used are also listed in the figure.Figure 4
k-f data taken from the measurement C-f curves of LaxZr1-xO2 thin films. The square symbols are measured from La0.09Zr0.91O2, and the solid line is modeled according to the HN law. The diamond symbols, triangle symbols, and circle symbols are tested, respectively, from the LaxZr1-xO2 with doping levels of 22%, 35%, and 63%, and solid lines in these cases are all simulated using the CS law. The parameters of CS and HN laws (α, τ, and n) are shown in brackets.Figure4 summarizes the frequency dependence of k-values for different La concentrations x. The zirconia thin film with an La concentration of x=0.35 had a relatively flat frequency response and a k-value of ~17 (on average). In the same way the La-doped zirconium oxides with doping level of 22% and 63% showed average values of 14 and 11, respectively. In contrast the lightly doped sample (9%) had a substantially increased k-value but suffered from a more severe dielectric relaxation. A k-value of 35 was obtained at 1 kHz, but this value was found to reduce to a k-value of 25 at 1 MHz.Thek-f relationship of the LaxZr1-xO2 with the doping level of 9% is better modeled by the HN law (rather than the CS law), and the relaxation parameter α,β1, and the relaxation time τ are 0.67, 0.53, and 2.56×10-5s, respectively. The k-values of the LaxZr1-xO2 dielectrics (x=0.22, 0.35, and 0.63) clearly show power-law dependence on frequency (f) and may be modeled using the CS law, k∝fn-1 (0 ≤ n ≤ 1), where the value of the exponent (n) indicates the degree of dielectric relaxation [6]. n values of 0.981, 0.98, and 0.985 are obtained for La compositions (x) of 0.22, 0.35, and 0.63, respectively. Comparing n values to k-values, it seems that the larger k-value suffers from more severe dielectric relaxation (smaller n value). The largest k-value appears in La0.09Zr0.91O2, but the most severe dielectric relaxation occurs, and this could be directly related to the size of crystal grains formed during annealing [25].Considering the lightly doped high-k thin films (CexHf1-xO2 and LaxZr1-xO2), the extracted data are both best modeled by the HN law, which show with significant permittivity enhancement accompanied by a more pronounced dielectric relaxation.For theCexHf1-xO2 and LaxZr1-xO2 dielectrics possible causes for the observed dielectric relaxation have been suggested [26]. It has been reported that a decrease in crystal grain size can cause an increase in the dielectric relaxation in ferroelectric relaxor ceramics [27, 28]. In addition to the doping level affecting the phase that the thin films crystallize during annealing, the doping level also affects the size of the crystal grains formed. It is possible therefore that the dielectric relaxation behavior observed in Figures 3 and 4 is due to the level of mechanical stress in the crystalline grains, which depends on the grain size, as has been found in the case of ferroelectric ceramics [23].


## 4. Conclusion

Doping hafnia and zirconia thin films with rare earth elements can stabilize the metastable tetragonal or cubic phase following annealing which enhances the dielectric constant. The level of enhancement is closely related to the doping level. Dielectric constants of 39 and 33 were obtained for La0.09Zr0.91O2 and Ce0.1Hf0.9O2 thin films, respectively. After taking into account the effect of the measurement system and the interfacial layer between the high-k dielectrics and silicon substrate, the change of the real permittivity with frequency has been modeled by either the CS or HN relationships. Comprehensive physical and mathematical models are provided and discussed for modelling the dielectric relaxation behaviour of high-k thin films under various conditions. For the lightly doped high-k thin films (Ce0.1Hf0.9O2 and La0.09Zr0.91O2), the measured data are both best modeled by the HN expression with significant permittivity enhancement and serious dielectric relaxation. The lanthanum-doped zirconium oxides (LaxZr1-xO2) with doping level from 22%, 35%, to 63% are best modeled based on the CS law in comparison with lightly doped LaxZr1-xO2. Films with larger k-values suffer from more severe dielectric relaxation, which is probably related to the size of crystal grains formed during annealing.



---

*Source: 241470-2012-01-18.xml*