This section is organized as follows. First, Subsection 3.1, the structural, surface morphology and composition of interested samples used this study are display. The THz optical and electrical properties of those samples are display in Subsection 3.2. Finally, we briefly study the optical properties of those samples in visible light spectrum for a comparison in the Subsection 3.3.
3.1. Surface Morphology, Composition and Structral Properties
In
Figure 1a, we present the SEM images of the top view of the as-deposited ITO thin film, revealing a smooth surface morphology. Additionally,
Figure 1b displays the cross-sectional view of the as-deposited ITO thin film, indicating a thickness of 275±6 nm.
Figure 1c illustrates the surface morphology of the ITO thin film treated with RTA at 400°C, showing a distinguishable structural transition from that of the as-deposited film while maintaining a flat surface. The thickness of ITO thin film treated with RTA at 400°C is not investigated (as specified by N/A in
Figure 1d). On the other hand, the surface morphology of the ITO thin film treated with slightly higher temperature of 600°C is showing a similar that that of RTA at 400°C. This is shown in
Figure 1e. We note further the thickness of the RTA-treated at 600°C ITO thin film is slightly reduced to 273±3 nm, as confirmed by the cross-sectional SEM image in
Figure 1f. This reduction in thickness may be attributed to the minimization of defects or an increase in density of surface area to volume ratio [
39]. Notably, a remarkable change in surface morphology of the ITO thin film treated with RTA at 800°C is observed, as depicted in
Figure 1g, where prominent but non-uniform particles are evident. From image analysis of the SEM pictures (ImageJ, University of Wisconsin-Madison, USA), we estimate that the average grain sizes are 80 nm
2, 175 nm
2 and 694 nm
2 for those annealed at 400°C, 600°C and 800°C, respectively. The results indicate a significant increase in grain sizes of the ITO thin films as the RTA temperature increases from 400 to 800°C, a phenomenon also observed in furnace annealing at 550°C for 1 hour [
40]. Moreover, the thickness of the ITO thin film treated with RTA at 800°C is approximately 273±2 nm, similar to that of the film treated at 600°C, as confirmed by the cross-sectional image shown in
Figure 1h. While we did not perform X-Ray Diffraction (XRD) measurements, it is plausible that our samples, as-deposited is amorphous and annealed are polycrystalline.
Next, AFM was employed to determine the surface roughness and particle size of the ITO thin films.
Figure 2a–d display the surface topography of the as-deposited, RTA-treated at 400°C, RTA-treated at 600°C, and RTA-treated at 800°C ITO thin films, respectively. The root mean square (RMS) roughnesses of the as-deposited ITO thin film is 1.5 nm, whereas RTA-treated at 400°C, 600°C, and 800°C ITO thin films are 3.4, 3.5 nm, and 8.3 nm, respectively. The increasing surface roughness with increasing temperature is suspected to be due to structural changes as the ITO thin films become polycrystalline and the grain sizes increase with heat treatment. This increasing trend is in agreement with RTA treated ITO thin film reported by Song et al [
25]. It is worth to noting that the RMS roughness nearly same for RTA treated at 400°C and 600°C ITO thin film, showing again similar trend in agrement with Song et al. [
25]. To confirm this hypothesis on grain sizes increase with heat treatment, the particle size of the films was examined using AFM, as shown in
Figure 5e–h. The calculated particle sizes of the as-deposited, RTA-treated at 400°C, RTA-treated at 600°C, and RTA-treated at 800°C ITO thin films are 37 nm, 62 nm, 85 nm, and 243 nm, respectively. These are in good agreement with the high-resolution SEM study (see
Figure 1). The trend of particle size increasing as the increased of annealing temperature was also found in previous literature [
41].
Energy-dispersive X-ray spectroscopy (EDS) was utilized for the elemental analysis of the samples of interest.
Figure 3a–c show the EDS spectra confirming the presence of all related elements, i.e., indium (In), tin (Sn), and oxygen (O). Furthermore, the weight and atomic percentages of all elements recorded by EDS for these samples are displayed in
Table 1. The EDS spectrscopy studiy on RTA-treated at 400°C ITO thin film is not conducted.
We note a trend of increasing weight and atomic percentage of oxygen (O) from the as-deposited (weight: 22.78%, atomic: 68.01%) to RTA-treated at 600°C (weight: 22.98%, atomic: 68.23%) to RTA-treated at 800°C (weight: 24.71%, atomic: 70.27%) ITO samples. This trend is suspected to be due to the increasing oxygen partial pressure in the environment. Similar phenomena were reported in spinel-type nickel manganite oxide (NiMn2O4-δ) polycrystalline samples treated at high temperatures of 800°C [
42]. Conversely, indium (In L) recorded higher values in both weight (69.48%) and atomic (28.75%) percentages in RTA-treated 600°C ITO samples compared to the as-deposited ones (weight: 67.44%, atomic: 28.05%). In contrast, tin (Sn L) exhibited lower values in both weight (7.54%) and atomic (3.02%) percentages in RTA-treated 600°C ITO samples compared to the as-deposited ITO thin film (weight: 9.77%, atomic: 3.93%). In a related study by Thirumoorthi et al. [
43], and Seki et al. [
44], revealed that the resistivity is decreased, carrier concentration is increased and mobility is increased as increasing Sn L concentrations due to the valence difference between Sn L and In L ion ions. Thus, we expect an improvement in THz electrical properties with RTA-treated at 600°C compared to the as-deposited and RTA-treated at 800°C ITO thin film.
3.2. THz Optical and Electrical Properties
The THz frequency-dependent optical properties of as-deposited and RTA-treated ITO thin films coated on HR silicon and the reference (bare HR silicon) were studied by using a home-made transmission-type THz-TDS system. The THz pulse train transmitted through the reference and ITO/HR silicon are shown in
Figure 4a. A small-time delay due to discernible phase shift of the THz pulse transmitted through the samples under study can readily be seen. Next, a Fast Fourier Transform (FFT) algorithm was applied to convert the THz time-domain signal to its spectral amplitudes. These are displayed in
Figure 4b. The data, i.e., time-delay or phase shift and attenuation imply significant interaction between the THz wave and materials investigated. The inset of
Figure 4b shows an example of a nearly perfect phase linearity of transmitted THz field through all the samples in the frequency range from 0.2 to 1.2 THz. This is a testimony of the reliability of the THz-TDS spectrometer.
For extraction of the optical constants, we employ the same approach as in our previous works on various thin films and nanostructures [
33,
45,
46,
47,
48,
49,
50]. We reiterate here for the benefits of readers: Consider the ITO layer (subscripted as 2) to be sandwiched between air (subscripted as 1) and a substrate (subscripted as 3) with complex refractive indices n
2, n
1 and n
3, respectively. The THz wave is assumed to be incident on this thin layer with a thickness of d from the air toward the substrate. Considering multiple reflections [
51], the transmitted THz electric field through the ITO-coated substrate can be written as
where E
0, t
12, t
23, t
31, r
21, and r
23 are the incident electric field, transmission, reflection coefficients between various interfaces; ω is the angular frequency; c is the speed of the light in vacuum; D is the thickness of substrate. Similarly, a reference THz field defined as the THz field transmitted through the bare substrate can be written as
The transfer function of the THz wave through the ITO layer then can be written as
Experimentally, t
exp(ω, n
2) is obtained from ratios of the frequency components of the Fourier transforms of the measured THz waveforms through the as-deposited and annealed ITO/HR-Si and the bare substrate. An error function
Error(ω, n
2) is defined as follows:
By minimizing the error function, the real and imaginary parts of the complex refractive index n2 of the thin film layer can be deduced.
Figure 5 shows the real (n) and imaginary (κ) part of refractive indices of as-deposited and RTA-treated ITO films at 400°C, 600° and 800°C as a function of frequency. In the frequency range of 0.2-1.0 THz, as-deposited ITO thin film exhibits a decreasing trend in n from ~ 65 to ~ 33, and κ from ~ 66 to ~ 27. Such decreasing trend and values of optical constants for as-deposited ITO thin film are very close to those of the DC-magnetron-sputtered and electron-beam-deposited ITO thin films reported previously by us [
45,
46,
47,
48,
51] as well as pulsed laser deposited ITO thin film reported by Wang et al. [
52]. This implies that the THz optical constants of ITO films deposited by various methods are very similar. The THz optical constants of RTA-treated ITO thin films also exhibit analogous decreasing trend with increasing frequency. Note, however, the n and κ of ITO films RTA-treated at 400°C and 600°C are almost identical, from ~ 143 to ~ 72 and ~ 131 to ~ 67, respectively over 0.2 – 1.0 THz, but much higher in values in comparison to those of the as-deposited ITO thin film. In contrast, the n values of RTA-treated at 800°C ITO thin film revert back to about the same values as those of the as-deposited sample while the extinction coefficients, κ, are somewhat lower, ranging from ~ 45 to ~ 19, as compared to ~ 66 to ~ 27 of the as-deposited ITO film in the same THz frequency band. These trends correlate with the electrical characteristics of the films to be discussed in a latter paragraph of this section.
Figure 5.
(a) The real (n) and (b) imaginary (κ) part of refractive indices of ITO films as a function of frequency. Black open square: as-deposited, red circles: RTA-treated at 400°C, green open triangles: RTA-treated at 600°C and blue open triangles: RTA-treated at 800°C.
Figure 5.
(a) The real (n) and (b) imaginary (κ) part of refractive indices of ITO films as a function of frequency. Black open square: as-deposited, red circles: RTA-treated at 400°C, green open triangles: RTA-treated at 600°C and blue open triangles: RTA-treated at 800°C.
We are not able to verify the THz refractive indices of those RTA-treated ITO thin films due to lack of published literature at the present time. Yet, the increasing of THz refractive index with RTA-treated at 400°C and 600°C is supported with the experimental measurements as well as the numerical estimate in reflectance of DC-magnetron sputtered ITO thin films deposited on glass substrate and post-annealed in the tubular furnace up to 300°C [
35]. From our SEM and AFM results, we assumed the increase of refractive indices for RTA-treated ITO films at 400°C and 600°C compared to those of as-deposited samples might be attributed to the increasing of grain size and an improvement in crystalline structure. The behavior of ITO films RTA-treated at 800°C is unique. Structural and morphology studies indicate that the grain sizes and surface roughness for these samples are much larger than those heat-treated at lower temperatures. This could be correlated with the observed trend in this work. We note that almost all of reported studies on the effects of annealing on ITO films were conducted at temperatures below 600°C [
53]. Further, these works focus on the optical properties of ITO films, in the visible region (400-700nm). Therefore, we also studied “optical “properties of these ITO thin film in the visible window and present them in the next sub-section.
We next extracted the transmittance of those samples using the ratio of frequency domain THz field amplitudes through the ITO/HR silicon to the frequency domain THz field through the reference (bare HR silicon), see
Figure 6a. The THz transmittance of as-deposited and RTA-treated ITO thin films are slow-varying in the band of 0.2 – 1.0 THz. The average transmittance of as-deposited, RTA-treated at 400°C, RTA-treated at 600°C and RTA-treated at 800°C are 27 %, 8%, 8 % and 39 %, respectively. We note the average THz transmittance of as-deposited ITO thin film is somewhat lower than those thinner (~ 100 nm) as-deposited ITO thin films (T ~ 40%) used in our earlier works [
54]. This is primarily attributed to the thickness dependent free carrier absorption of ITO thin film. The THz transmittance of ITO thin film treated with RTA at 400°C and 600°C are as low 8 %. This is consistent with higher extinction coefficients of RTA-annealed ITO films at either 400°C or 600°C shown in
Figure 5b. Our results indicate that annealing temperature is a critical parameter affecting transmittance and reflectance of the ITO film.
In
Figure 6b and c, it can be seen that the absorbance and reflectance of ITO thin films depend strongly on the annealing temperature. Here, absorbance, A, is calculated by A=αd/2.303, where the absorption coefficient α is related to the extinction coefficient κ by 𝛼=4𝜋𝜅/𝜆𝐴. In turn, reflectance of the samples is determined by R =1-T-A. We note that the decreasing trend of reflectance for increasing temperature is consistent with the square root dependence of conductivity according to the Hagen-Rubens approximation [
55],
, by using the conductivities σ listed in
Table 2 of this work. Both absorbance and reflectance are found to be much higher for ITO films heat treated by RTA at 400°C and 600°C than that of as-deposited ones. We speculate that RTA-treated ITO thin films at either 400°C or 600°C could potentially be useful as dichroic mirrors or absorbers for the THz band [
35,
56]. On the other hand, ITO thin films RTA-treated beyond 800°C could be an effective THz TCE material for THz EO devices. We note the transmittance of ITO-coated device can be further enhanced by employing a grating-like, i.e., wire-grid structure, as we have demonstrated earlier [
54].
Electrical characteristics of the ITO films can be extracted from the optical constants determined above. The complex conductivity, σ* is related to the dielectric function,
Where ε
0 is the free-space permittivity, ε
∞ ~ 4 is the high-frequency dielectric constant of ITO, consisting of contributions from bound electron. Therefore, we can determine THz conductivities of ITO samples from their complex refractive index, i.e.,
In degenerate semiconductors such as ITO, the Drude free-electron model has been most widely used to describe its electrical characteristics. That is,
[
57], where ω
p and τ are plasma frequency and scattering time of carriers in the material. The real part of conductivity, Re{σ}, exhibits a maximum at zero (DC) frequency, and its value decreases with increasing frequency. For Im{σ} [
46,
47,
48], it can only be positive and approaches a maximum as the frequency approaches that of the inverse of the scattering rate. We have found, however, that the THz conductivities of ITO films and nanostructures exhibit non-Drude-like behavior, e.g., depressed values of DC conductivity and negative values for Im{σ} [
46,
47,
48]. Therefore, the Drude-Smith model, which takes into account the carrier localization effect, is chosen to fit the experimentally deduced conductivity of the ITO material, In this model [
58],
where c is the persistence of velocity and that varies from -1 to 0. The value of c ( -1 < c < 0) can be associated with the degree of backscattering that a carrier suffers after a collision. When parameter c = 0, the Drude model is recovered. Carriers undergo complete backscattering or localization for c = -1. Although it is a rough approximation, this simple empirical model works remarkably well for a variety of materials. The SEM data show that our ITO samples consists of grains tens to hundreds nanometer in dimension. Thus the non-Drude-like behavior can be contributed by backscattering of carriers at grain boundaries [
59,
60]. On the other hand, the dopants in ITO will become positive ions after providing free electrons [
61,
62]. Due to the Coulomb interaction between dopant ions and free electrons, the former will also become scattering centers hindering movement of the electrons as well. As we shall see in
Figure 7, the Drude-Smith model fits the conductivity data better than the simple Drude model. All of the fitting parameters are summarized in
Table 2. Meanwhile, the extracted motilities (µ) are 29 cm
2/V∙s, 110 cm
2/V∙s, 226 cm
2/V∙s and 70 cm
2/V∙s; plasma frequency (ω
p) are 3315 rad.THz, 3711 rad.THz, 2580 rad.THz and 1648 rad.THz; scattering times (τ) are 6 fs, 12 fs, 21 fs and 20 fs; the fitting parameters (c) are -0.83, -0.68, -0.63 and -0.87; carrier concentration (N) are 2.17×10
20 cm
-3, 2.72×10
20 cm
-3, 1.31×10
20 cm
-3 and 5.58×10
19 cm
-3; conductivities (σ) are 1019 Ω
-1‧cm
-1, 4897 Ω
-1‧cm
-1, 4881 Ω
-1‧cm
-1and 628 Ω
-1‧cm
-1; and resistivity (ρ) are 9.8×10
-4 Ω·cm, 2.0×10
-4 Ω·cm, 2.0×10
-4 Ω·cm and 15.9×10
-4 Ω·cm for as-deposited, RTA-treated at 400°C, RTA-treated at 600°C and RTA-treated at 800°C, respectively.
3.3. Annealing Effects on UV-VIS-NIR Optical Properties of ITO Films
It is interesting to compare the optical properties of samples studied in this work in the ultraviolet (UV), visible (VIS) and near infrared (NIR) band as a bearing on their far infrared (FIR) or THz frequency range. We have experimentally measured transmittance and reflectance of as-deposited and RTA-treated ITO thin films coated on fused silica as well as that of the substrate itself in the energy range of 0.6 eV – 6 eV, or the wavelength band of 200 – 2000 nm. These are shown in
Figure 8. The average transmittance in the VIS-NIR region (400-800nm) band for as-deposited ITO film was Tav ~ 72%. After annealing at 400°C, 600°C, and 800°C, Tav ~ 68 %, ~ 80 % and 83%, respectively. Considering the contribution by the substrate, In, all samples can be considered highly transparent in this band. Higher transmittance of samples after heat treatment by annealing can be attributed to increasing sizes of crystallites after annealing as shown in the SEM data (see
Figure 1). It is well-known that improved crystallinity leads to reduced scattering of incident light, which enhanced their transmission through the film. Transmittance drops off sharply in the NIR and UV band for all samples. The absorption edge shifted to shorter wavelength for samples annealed at increasingly higher temperatures. This point will be addressed further later on.
The average reflectance of as-deposited ITO thin film in the UV-VIS-NIR region, Rav ~ 14 %. Whereas, the reflectance of ITO thin film with RTA-treated at 400°C, 600°C and 800°C are Rav ~ 10 %, Rav ~12% and Rav ~ 13%, respectively in VIS spectrum wavelengths. In addition, the reflectance curve shows small but yet significant oscillations due to multiple reflections of light within the sample. This also reasonably good agreement with the previous finding by others [ref]. Further, the results show a general increasing trend in reflectance in NIR region with temperature. The trend is quite similar to the literatures reported previous [ref]. It is also observed that the reflectivity increased to 70 % as the annealing temperature was increased to 600°C. This is tentatively attributed to increase of carrier concentration after annealing. Yet, the reflectivity is shown slightly lower for ITO thin film RTA-treated at 800°C than those of other RTA-treated samples. This is probably impacted by surface roughness effects [ref]. Thus, our results indicating annealing can be used to tune the optical transmittance and reflectivity of ITO thin film over a broad band from UV-VIS-NIR to the THz region by changing structural, surface morphology and electrical properties.
We now focus on spectral regions of ITO films that are weakly absorbing and also shown prominent maximum (TM) and minimum (Tm) transmission spectrum, i.e., the visible window of 400-800 nm (see
Figure 8a). Following Swanepoel [
63], the refractive index of ITO the film in its weakly absorbing VIS-NIR (400 to 1000 nm) band can be estimated as
, where
, T
M and T
m are envelopes of maxima and minima of the transmission curve and
is the refractive index of fused silica substrate and is determined using its transmission spectrum (T) to be ~ 1.43, and found to be fairly independent of wavelength. The calculated refractive index n of as-deposited and those RTA-treated ITO thin film in wavelength range from 400 to 1000 nm are shown in
Figure 9a. Note that the refractive index of 1.74 in this band is almost wavelength-independent and approximately the same for all types of ITO thin film studied in this work. Note a small dispersion-like feature of the refractive index was observed near 500 nm. It can be seen that the refractive indices of annealed ITO films are lower than those of the as-deposited one in this dispersion region. In the literature, the reported refractive indices for as deposited as well as annealed ITO thin films vary from 1.6 to 2 in this band [
64].
The thickness d of those ITO thin films can be determined using the envelopes of maxima and minima of the transmission curve in the visible spectra,
, where n
1 and n
2 are the refractive indices at two adjacent maxima or minima, T
M or T
m , corresponding to wavelengths at λ
1 and λ
2, respectively. The estimated thickness of as-deposited and RTA-treated ITO thin films at 400°C, 600°C and 800°C are ~ 470 nm, ~ 410 nm, ~ 483 nm and ~ 468 nm, respectively. Knowing the thickness of ITO thin films, the absorption coefficient α of these samples are extracted using the relationship between transmittance (T), reflectance (R) in VIS spectrum range and thickness d, given:
. Next, the imaginary part of refractive index κ as the function of wavelength is evaluated using another simple formula shown as:
, plotted in
Figure 9b. We found that κ in the wavelength range of 400 to 1000 nm is slow varying and approximately 0.04, 0.18, 0.13 and 0.7 for as-deposited and RTA-treated ITO thin films at 400°C, 600°C and 800°C, respectively. Such estimated values are also in agreement with those reported in previous works [
64].
Crystalline ITO has been shown to be a direct gap semiconductor [
65]. From the absorption coefficient α as plotted in
Figure 10a, the optical bandgap Eg can be extracted using the extrapolation method proposed by Tauc: (αhʋ)2 = A (hʋ-Eg) [
66]; Where A is a constant and h = 6.626×10-34 J.s is the Planck constant. These are shown for the samples under study in
Figure 9b. The optical band gap Eg of as-deposited ITO thin film was found to be 3.6 eV. On the other, the optical band gap Eg for RTA-treated ITO thin films at 400°C, 600°C and 800°C are estimated to be approximately 3.97 eV, 4.10 eV and 4.11 eV, respectively. The increasing Eg as the annealing temperature increases ascribed primary to the improvement of the crystallinity and reducing the defects in thin film. Nevertheless, these values are good agreement with typical value reported previously. For example, Abd-Elnaiem et al. [
67], extracted Eg for as-deposited and for annealed ITO thin film with thickness of 80 nm at 300°C to be ~ 3.2 eV and~ 3.8 eV, respectively. Meanwhile, Eg for a 200 nm-thick ITO films heat-treated with RTA at 600°C was determined to be as high as 4.04 eV by Song et al. [
25].