*Article* **Structural, Optical, Magnetic and Electrical Properties of Sputtered ZnO and ZnO:Fe Thin Films: The Role of Deposition Power**

**Ahmed Faramawy 1,2,\* , Hamada Elsayed 3,4 , Carlo Scian <sup>1</sup> and Giovanni Mattei <sup>1</sup>**


**\*** Correspondence: ahmed\_faramawy@sci.asu.edu.eg

**Abstract:** Structural, optical, magnetic, and electrical properties of zinc oxide (henceforth, ZO) and iron doped zinc oxide (henceforth, ZOFe) films deposited by sputtering technique are described by means of Rutherford backscattering spectrometry, grazing incidence X-ray diffraction, scanning electron microscope (SEM), UV–Vis spectrometer, vibrating sample magnetometer, and room temperature electrical conductivity, respectively. GIXRD analysis revealed that the films were polycrystalline with a hexagonal phase, and all films had a preferred (002) c-axis orientation. The lattice parameters *a* and c of the wurtzite structure were calculated for all films. The *a* parameter remains almost the same (around 3 Å), while c parameter varies slightly with increasing Fe content from 5.18 to 5.31 Å throughout the co-deposition process. The optical gap for undoped and doped ZO was obtained from different numerical methods based on the experimental data and it was increased with the increment of the concentration of Fe dopant from 3.26 eV to 3.35 eV. The highest magnetization (4.26 <sup>×</sup> <sup>10</sup>−<sup>4</sup> emu/g) and lowest resistivity (4.6 <sup>×</sup> <sup>10</sup><sup>7</sup> <sup>Ω</sup>·cm) values of the ZO films were found to be at an Fe content of 5% at. %. An explanation for the dependence of the optical, magnetic, and electrical properties of the samples on the Fe concentrations is also given. The enhanced magnetic properties such as saturated magnetization and coercivity with optical properties reveal that Fe doped ZO thin films are suitable for magneto-optoelectronic (optoelectronic and spintronics) device applications.

**Keywords:** diluted magnetic semiconductor; TM doping; sputtering; magnetic properties; optical band gap; resistivity

### **1. Introduction**

Nowadays, studies on ZO, which is one of the most important 3rd generation semiconductors, support a new stage of comprehensive use due to its multi-functional characteristics. It plays an important role in a very wide range of applications in fields such as: biosensors, biomedical and antibacterial active media, in biological fields [1,2]; UV light-emitting devices, optical waveguides, solar cells in electro-optical fields [3]; transparent high power electronics, varistors, thin film transistors in electronics; gas sensors, water purification and solar photocatalytics in the environmental field [1,4,5]. Moreover, since 2001, many studies have been carried out in the field of Diluted Magnetic Semiconducting Oxides (DMSO) which combine two interesting properties—semiconducting and magnetic—whether to look for a new compound which might be ferromagnetic at room temperature or to find materials which might have a large magnetic moment [6]. The interest in ZO is fueled and fanned by its prospects in optoelectronics applications, having the necessary excellent conditions of ultraviolet and blue emissions because of its unique of chemical and physical properties of a large direct energy gap (e.g., ~3.37 eV at 300 K) and a large exciton binding energy (60 meV) at room temperature which makes it being invaluable compared to numerous other nanomaterials [7]. It is also a promising candidate

**Citation:** Faramawy, A.; Elsayed, H.; Scian, C.; Mattei, G. Structural, Optical, Magnetic and Electrical Properties of Sputtered ZnO and ZnO:Fe Thin Films: The Role of Deposition Power. *Ceramics* **2022**, *5*, 1128–1153. https://doi.org/10.3390/ ceramics5040080

Academic Editors: Amirhossein Pakseresht and Kamalan Kirubaharan Amirtharaj Mosas

Received: 14 October 2022 Accepted: 28 November 2022 Published: 1 December 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/).

for semiconductor spintronics applications; Dietl et al. [8] predicted a Curie temperature of >300 for Mn-doped ZO. In addition, n-type doping in Fe-, Co-, or Ni-alloyed ZO was also predicted to stabilize high—Curie—temperature ferromagnetism. Moreover, the wellknown room temperature ferromagnetism (RTFM) is one of the more studied phenomena in pure ZO nanostructures doped with magnetic ions. This phenomenon has been found at low temperatures as well, and it has been studied both experimentally and theoretically [9].

Furthermore, the remarkable morphological features of ZO make it a striking material that is widely used as a counterpart of a noteworthy scope of applications. As for ZO's morphological perspectives, synthesis and fabrication procedures also play an important role. The various parameters used (namely surfactant, temperature, concentration, and time, etc.) play a crucial role in the growth of various forms of synthetic processes. An enormous number of different methods for the synthesis of surface structured ZO as nanopowder, composites, and films have been published worldwide. There are widely noted techniques for synthesizing different ZO nanomaterials, composites and doped heterostructures, including precipitation, wet-chemical technique, sol-gel procedure, hydrothermal method, solvothermal procedure, sputtering, coating and microwave techniques, etc. [1].

For several of the above-mentioned applications, a stable, high, and reproducible p-doping is mandatory. Though progress has been made, this aspect still represents a major problem. The contribution of the current ZO research is not only on the same topics as earlier, but also includes nanostructures and nanotechnologies, new growth and doping techniques, and focuses more on the application-related aspects of daily life [9,10].

In this study we planned to control and enhance the optical band gap, and the magnetic properties of Zinc Oxide (ZO) thin films doped by a transition element (TM).

Fe is a magnetic transition metal that has several oxidation states, an electronegativity of 1.83 and ionic radii that range from 0.61 to 0.78 Å (depending on the magnetic and oxidation states) [10]. Because the properties of iron are similar to those of aluminum or gallium, it is considered a potentially attractive alternative for the doping of ZO thin films [11]. However, the performance of iron as a dopant is complicated by the existence of 3d electrons in its electronic configuration. Currently, most studies on Fe-doped ZO thin films focus on their wide applications. As we are aware, fewer reports have investigated Fe-doped ZO thin films as transparent conducting oxides [12–14]. Due to the complex electronic structure of iron, it has been difficult to control the different properties of these films. Moreover, iron-doped ZO thin films have been synthesized by several techniques including electrodeposition, spray pyrolysis, spin coating, ion implantation, and the sol–gel method. Sputtering is the most widely used technique for the preparation of Fe-doped ZO thin films [15] because it provides the ability to control the properties of the films through several deposition parameters, such as the sputtering mode, power, sputtering and reactive gas pressures, substrate temperature, and target composition and allows deposition of different materials [16]. Three sputtering techniques have been used to synthesize Fe-doped ZO thin films:


Several of the studies using RF-sputtering for Fe-doped ZO thin film have investigated the effect of doping concentration on the properties of the films [17]. Others have studied the influence of the deposition parameters on the properties of films with fixed dopant concentrations [19].

In this present study, a thin film of iron (Fe) doped zinc oxide (ZO), with various Fe-concentrations, is investigated using RF-magnetron sputtering on silica and silicon substrates. We report the influence of the Fe doping on the microstructure, optical, magnetic and electrical resistivity of the ZO films. The results indicated that ZO film with 5.0 at. % doping-concentration can attain a higher energy gap and magnetic properties, and better electrical resistivity.

#### **2. Materials and Methods**

Undoped and doped ZO films with Fe were deposited on single-crystal n-type Si (100) wafer and silica (HSQ 100 from Heraeus) glass preliminarily washed in piranha solution for magnetron sputtering. The ZO target (99.95% purity, 100-mm diameter) is placed in a radio-frequency driven (RF) source while setting the power at 100 W and Fe target in a direct current (DC) source providing powers of 3, 5, 10, and 15 W to produce (1, 2.4, 4.8, 5) at. % of Fe, respectively. Moreover, in order to control the sputtering yield of Fe and instead of reducing the DC power (which produced instabilities of the DC source), we performed the DC sputtering through a grid with a different mesh size to block part of the Fe atoms. During the deposition, the substrate was rotated at a low speed to enhance the thickness uniformity of deposited films. The two sources were tilted of an angle 30◦ with respect to the axis normal to the sample surface allowing a simultaneous deposition of the two targets and the source-target distance was about 15 cm. The thickness and the composition of the films were properly tuned through an optimized combination of the electrical power of the source and the opening time of the shutters placed in front of the targets. The thickness of the films was 150 nm, which was controlled by the deposition time. Other detailed conditions about the deposition process of the films are summarized in Table 1.

**Table 1.** Typical deposition parameters for sputtering deposition.


Furthermore, the film thickness was measured using a thickness profile-meter and by Rutherford back scattering (RBS), performed with a HVEC 2.5 MeV Van de Graaff accelerator at Legnaro National Laboratory (LNL)—INFN, using 2 MeV α -particles at normal incidence and a scattering angle of 20◦ , which was also used to determine the contents by atomic percentage of Zn, O and Fe ions. Furthermore, the hexagonal phase wurtzite structure of samples was confirmed by grazing incidence X-ray diffraction (GIXRD) (Philips X'Pert PRO MRD triple axis diffractometer (angular resolution 0.01◦ , accuracy 0.001◦ )) and the scanning electron microscope (SEM), Zeiss Field-emission SEM (FE-SEM), was used to determine the surface morphology. The optical properties of the films were measured using UV-Vis optical spectroscopy, Jasco V-670 UV-Vis-NIR Spectrophotometer with working range 190–2500 nm, in order to determine the optical parameters such as band gap, Urbach energy and refractive index for all investigated samples.

Likewise, the magnetic hysteresis loops for the doped and undoped films were investigated by a vibrating sample magnetometer (VSM) (LakeShore model no. 7410 USA) in a maximum applied field of 20 kOe. Finally, the IV characteristics and electric resistivity was also measured for undoped and doped films by a two probe method with a KEITHLEY 6514 system electrometer. For simplicity, throughout the present work, the "balanced stoichiometry" for the Fe doped ZnO will be named ZOFe. The DC power values are indicated as follows: ZOFe/3 W, ZOFe/5 W, ZOFe/10 W and ZOFe/15 W for 3, 5, 10 and 15 Watt for DC deposited power of the Fe source target.

#### **3. Results**

#### *3.1. Elemental Analysis*

The elemental composition of ZO and ZOFe films was measured by Rutherford backscattering spectrometry (RBS). The RBS data analysis is reported in Figure 1 where the experimental data for all investigated films is superimposed to the simulation. The quality of the fit is the direct proof of the sample homogeneity.

**Figure 1.** RBS random spectra for the ZO and ZOFe thin films deposited on Si.

**Figure 1.** RBS random spectra for the ZO and ZOFe thin films deposited on Si. One can notice that the surface edge for the different elements (Zn, O, Si and Fe) can be observed clearly from Figure 1. The large peak on the right-hand side located in the high energy range of the spectrum is due to the zinc atoms whereas the low energy peak signal is related to the oxygen, which is superimposed over the continuous spectrum of a Si substrate. Moreover, small peaks were observed on right-hand side of the Zn peak which indicate the Fe. The yield account of these peaks increased with increase in the DC power value in the co-deposited process from 3 W to 15 W.

RBS is not equally sensitive to various elements. The height of the RBS signal depends on the scattering cross section which is proportional to the square of the atomic number of the element. Thus, for light elements, the RBS signal is much lower than for heavy ones [20]. Therefore, it is quite difficult to separate the signal from such elements as nitrogen, carbon or oxygen, from the high background of much heavier substrates as Si. Consequently, it is not easy to evaluate the contribution of light elements. Zinc, oxygen and iron contents in the ZO and ZOFe thin films were evaluated from the RBS random spectra. It has been found that the value O: Zn content is approximately equal 1 and the results are summarized in Table 2.


**Table 2.** RBS results for the estimation of Zn and Fe content (at. %) in the investigated films.

#### *3.2. GIXRD Analysis*

Figure 2 depicts the GIXRD pattern of the ZO thin films. One can see that it exhibits a hexagonal crystal phase, which was confirmed with a standard JCPDS card (PDF No. 36–1451). In the GIXRD pattern, the two dominant peaks, (002) and (103) were obtained for all films without any secondary phase.

All films exhibited preferential (002) c-axis orientation, which grew perpendicular to the substrate surface. This can interpreted as follows: the sputtered ZO and ZOFe thin films with a c-axis have a preferential orientation due to the lowest surface free energy of (002) plane in ZO. In the equilibrium state, the films grow with the plane of the lowest surface free energy parallel to the surface if there is no effect from the substrate. Additionally, the sp<sup>3</sup> hybridized orbit in zinc oxide forms a tetrahedral coordination. Because each apex is parallel to the c-axis in the wurtzite structure, ZO films prefer to grow in the (002) direction [21].

Table 3 represents some significant data estimated from the diffracted peaks such as 2θ, d-spacing, (hkl), peak intensity, crystallite size, and the texture coefficient for all deposited films.

The texture coefficient *TC(hkl)* is another perspective for preferred orientation and/or preferred growth. *TC(hkl)* can provide quantitative information on the preferred crystal orientation which can be calculated from the expression [22,23]

$$TC\_{\left(hkl\right)} = \frac{I\_{\left(hkl\right)} / I\_{o\left(hkl\right)}}{\left(1/n\right) \ \sum\_{n} I\_{\left(hkl\right)} / I\_{o\left(hkl\right)}}\tag{1}$$

where *I(hkl)* is the measured intensity and *Io(hkl)* is the standard JCPDS intensity and *n* is the number of diffraction peaks. If *TC(hkl)*∼1 for all the (*hkl*) planes considered, then the films have a randomly oriented crystallite similar to the JCPDS standard, while values higher than 1 indicate the abundance of grains in a given (*hkl*) direction.

for all films without any secondary phase.

Figure 2 depicts the GIXRD pattern of the ZO thin films. One can see that it exhibits a hexagonal crystal phase, which was confirmed with a standard JCPDS card (PDF No. 36–1451). In the GIXRD pattern, the two dominant peaks, (002) and (103) were obtained

All films exhibited preferential (002) c-axis orientation, which grew perpendicular to the substrate surface. This can interpreted as follows: the sputtered ZO and ZOFe thin films with a c-axis have a preferential orientation due to the lowest surface free energy of (002) plane in ZO. In the equilibrium state, the films grow with the plane of the lowest surface free energy parallel to the surface if there is no effect from the substrate. Additionally, the sp3 hybridized orbit in zinc oxide forms a tetrahedral coordination. Because each apex is parallel to the c-axis in the wurtzite structure, ZO films prefer to grow in the

*3.2. GIXRD Analysis* 

(002) direction [21].

**Figure 2.** GIXRD patterns of ZO and ZOFe thin films deposited on SiO2. **Figure 2.** GIXRD patterns of ZO and ZOFe thin films deposited on SiO<sup>2</sup> .


Table 3 represents some significant data estimated from the diffracted peaks such as 2θ, d-spacing, (hkl), peak intensity, crystallite size, and the texture coefficient for all de-**Table 3.** Planes angle 2θ, d-spacing, intensity I, crystallite size D, and texture coefficient TC for different DC depositing powers of sputtered ZO and ZOFe from GIXRD.

Values 0 < *TC(hkl)* < 1 indicate the lack of grains oriented in that direction. As *TC(hkl)* increases, the preferential growth of the crystallites in the direction perpendicular to the *hkl* plane is the greater. Since two diffraction peaks were used ((002) and (103)), the maximum value *TC(hkl)* possible is two. increases, the preferential growth of the crystallites in the direction perpendicular to the *hkl* plane is the greater. Since two diffraction peaks were used ((002) and (103)), the maximum value *TC(hkl)*possible is two. The calculated texture coefficient results are displayed in Table 3 and were found to

Values 0 < *TC(hkl)* < 1 indicate the lack of grains oriented in that direction. As *TC(hkl)*

The texture coefficient *TC(hkl)* is another perspective for preferred orientation and/or preferred growth. *TC(hkl)* can provide quantitative information on the preferred crystal

(ଵ ) ⁄ ∑ ூ(ೖ) ூ ⁄ (ೖ)

(1)

() <sup>=</sup>ூ(ೖ) ூ <sup>⁄</sup> (ೖ)

where *I(hkl)* is the measured intensity and *Io(hkl)* is the standard JCPDS intensity and *n* is the number of diffraction peaks. If *TC(hkl)*~1 for all the (*hkl*) planes considered, then the films have a randomly oriented crystallite similar to the JCPDS standard, while values higher

*Ceramics* **2022**, *5,* FOR PEER REVIEW 7

orientation which can be calculated from the expression [22,23]

than 1 indicate the abundance of grains in a given (*hkl*) direction.

The calculated texture coefficient results are displayed in Table 3 and were found to be maximum for the (002) plane and increased by increasing iron content, as shown in Figure 3. Moreover, the (002) polar facet of Fe doped ZO films increased significantly compared to the (103) facet. This also indicates the preferential growth along the c-axis orientation. be maximum for the (002) plane and increased by increasing iron content, as shown in Figure 3. Moreover, the (002) polar facet of Fe doped ZO films increased significantly compared to the (103) facet. This also indicates the preferential growth along the c-axis orientation.

**Figure 3.** Texture coefficient, TC, for (002) and (103) plans of ZO, Fe-doped ZO deposited at different DC power. **Figure 3.** Texture coefficient, TC, for (002) and (103) plans of ZO, Fe-doped ZO deposited at different DC power.

On other side, the crystallite size may be estimated from the full-width at half-maximum (FWHM) of (002) and (103) diffraction peak using the Debye–Scherrer formula [24] which is expressed as: On other side, the crystallite size may be estimated from the full-width at halfmaximum (FWHM) of (002) and (103) diffraction peak using the Debye–Scherrer formula [24] which is expressed as:

$$D = \frac{0.9\ \lambda}{\beta \cos \theta\_{hkl}} \tag{2}$$

where *D* is the average crystallite size, *λ* is the wavelength of the X-ray and takes 1.54 Å for Cu Kα, *β* is the peak width of two peaks at half maximum height (FWHM) in radians, and *θhkl* is the diffraction angle of the crystal plane (*hkl*). The values of the average crytallite size are given in Table 3. where *D* is the average crystallite size, *λ* is the wavelength of the X-ray and takes 1.54 Å for Cu Kα, *β* is the peak width of two peaks at half maximum height (FWHM) in radians, and *θhkl* is the diffraction angle of the crystal plane (*hkl*). The values of the average crytallite size are given in Table 3.

In addition, using the interplanar spacing values (d-spacing) and the related hkl parameters for hexagonal systems, the lattice parameters (*a*, *b*, and *c* where *a* = *b*) of the wurtzite cell of ZO and ZOFe thin films have been calculated by using following formula [22,25]:

$$\frac{1}{d^2} = \frac{4}{3} \left( \frac{h^2 + hk + k^2}{a^2} \right) + \left( \frac{l^2}{c^2} \right) \tag{3}$$

Interplanar d-spacing is calculated according to Bragg's law:

$$\mathbf{n}\lambda = \mathcal{Z}\,\mathrm{d}\_{\mathrm{(hkl)}}\mathrm{Sin}\,\theta\,\tag{4}$$

where d is lattice spacing, a and c are lattice constants, h, k, l are miller indices, θ is the angle of corresponding peak and λ is the wavelength of X-ray used (1.5402 Å). Considering Bragg's law, it is possible to rewrite Equation (4) as follows:

$$\frac{4\sin^2\theta}{\lambda^2} = \frac{4}{3} \left( \frac{h^2 + hk + k^2}{a^2} \right) + \left( \frac{l^2}{c^2} \right) \tag{5}$$

There are two unknowns to be calculated in the formula above. Because of this, while calculating lattice constant *a*, a peak in the form of (h00) should be used to eliminate *c* from the equation. As an alternative, the (00l) peak should be used to get an equation with just one unknown for calculating the *c* constant. Following correct peak selection, the following equations for *a* and *c* constants are obtained [25,26].

$$a = \frac{\lambda \cdot \sqrt{h^2 + hk + k^2}}{\sqrt{3}\sin\theta} \tag{6}$$

$$\mathcal{L} = \frac{\lambda \cdot l}{2 \sin \theta} \tag{7}$$

The calculated lattice parameters for all films are listed in Table 4. As can be seen from Table 4, the *a* parameter remains almost the same, while the *c* parameter varies slightly due to an increase in Fe content throughout the co-deposition process. As we noticed from XRD, no obvious peaks corresponding to iron metal or iron oxide phase were observed. Furthermore, the wurtzite structure of ZO remained unchanged after iron modification. There is a slight left shift in the (002) and (103) peaks of the ZOFe films which could be interpreted as follows. Fe exists stably in both valence states, i.e., Fe2+ and Fe3+ state. The peaks would shift towards a lower angle if Fe exists in Fe2+ (0.78 Å) state due to its larger ionic radius than Zn2+ (0.74 Å), as seen in Figure 4. This confirms that iron has been doped into the ZO lattice. However, if it exists in the Fe3+ (0.68 Å) state, the peak shift will occur towards the higher angle side. As a result, we may conclude that the iron in our samples is predominantly in the 2+ state. *Ceramics* **2022**, *5,* FOR PEER REVIEW 9 shift will occur towards the higher angle side. As a result, we may conclude that the iron in our samples is predominantly in the 2+ state.

**Figure 4.** The peak shifts towards lower angle for all deposited films in planes (**a**) (002) (**b**) (103). **Figure 4.** The peak shifts towards lower angle for all deposited films in planes (**a**) (002) (**b**) (103).

As well, the volume *V*, the internal parameter *U*, and the bond length *L* of wurtzite unit cell for hexagonal ZO and ZOFe thin films have been calculated using the following

The obtained values are tabulated in Table 4 and the variation in , and *L* after Fe incorporation may be attributed to the differences between the atomic radii of the host

Moreover, the average uniform strain, *ezz*, in the lattice along the c-axis in the randomly oriented ZO and ZOFe films has been estimated from the lattice parameters as

> ௭௭ <sup>=</sup> ିబ బ

Where is the lattice parameter of the ZO film calculated from the (002) peak of XRD pattern and the *c*0 is the lattice parameter for the ZO bulk. For hexagonal crystals, the stress (σ) in the plane of the film can be calculated using the biaxial strain model [28,29]:

= ቂ2ଵଷ <sup>−</sup> (భభାభమ)యయ

where *Cij* are elastic stiffness constants for ZO, (*C*11 = 2.1 × 1011 N/m2, *C*33= 2.1 × 1011 N/m2, *C*12 = 2.1 × 1011 N/m2, and *C*13 = 2.1 × 1011 N/m2). This yields the following numerical rela-

The calculated values of stress (*σ*) in the undoped and doped films are listed in Table 4. The negative sign indicates that the films are in a state of compressive stress which not seen in pure ZO film. The total stress in the film typically consists of two components: intrinsic stress which introduced by impurities, defects and lattice distortions in the crystal, and the extrinsic stress familiarized by the lattice mismatch and thermal expan-

tion for the stress derived from the change in the 'c' lattice parameter:

భయ

(ିଶ) = −4.5 ൈ 10ଵଵ௭௭ (11)

<sup>ସ</sup> ; =  c<sup>ଶ</sup> <sup>ቀ</sup>

ଵ <sup>ଶ</sup>  −ቁ

ଶ <sup>+</sup> <sup>మ</sup> ଷ ൨ ଵ/ଶ

(9)

ቃ ௭௭ (10)

(8)

<sup>ଶ</sup> ଶ; = <sup>ଵ</sup>

ଷ ቀ మ <sup>మ</sup><sup>ቁ</sup> <sup>+</sup> <sup>ଵ</sup>

= √ଷ

material (ZO) and the dopant (Fe) [27].

expressions [22,25,27]:

follow:


**Table 4.** Values of lattice constants, cell volume, internal parameter, bond length, strain and the compressive stress for the wurtzite hexagonal ZO and ZOFe films deposited at different DC powers.

As well, the volume *V*, the internal parameter *U*, and the bond length *L* of wurtzite unit cell for hexagonal ZO and ZOFe thin films have been calculated using the following expressions [22,25,27]:

$$V = \frac{\sqrt{3}}{2}a^2 \text{c; } \mathcal{U} = \frac{1}{3} \left(\frac{a^2}{c^2}\right) + \frac{1}{4}; \ L = \left[c^2 \left(\frac{1}{2} - \mathcal{U}\right)^2 + \frac{a^2}{3}\right]^{1/2} \tag{8}$$

The obtained values are tabulated in Table 4 and the variation in *a*, *c* and *L* after Fe incorporation may be attributed to the differences between the atomic radii of the host material (ZO) and the dopant (Fe) [27].

Moreover, the average uniform strain, *ezz*, in the lattice along the c-axis in the randomly oriented ZO and ZOFe films has been estimated from the lattice parameters as follow:

$$e\_{zz} = \frac{c - c\_0}{c\_0} \tag{9}$$

where *c* is the lattice parameter of the ZO film calculated from the (002) peak of XRD pattern and the *c*<sup>0</sup> is the lattice parameter for the ZO bulk. For hexagonal crystals, the stress (σ) in the plane of the film can be calculated using the biaxial strain model [28,29]:

$$
\sigma = \left[ 2\mathbf{C}\_{13} - \frac{(\mathbf{C}\_{11} + \mathbf{C}\_{12})\mathbf{C}\_{33}}{\mathbf{C}\_{13}} \right] e\_{zz} \tag{10}
$$

where *<sup>C</sup>ij* are elastic stiffness constants for ZO, (*C*<sup>11</sup> = 2.1 <sup>×</sup> <sup>10</sup><sup>11</sup> N/m<sup>2</sup> , *<sup>C</sup>*33= 2.1 <sup>×</sup> <sup>10</sup><sup>11</sup> N/m<sup>2</sup> , *<sup>C</sup>*<sup>12</sup> = 2.1 <sup>×</sup> <sup>10</sup><sup>11</sup> N/m<sup>2</sup> , and *<sup>C</sup>*<sup>13</sup> = 2.1 <sup>×</sup> <sup>10</sup><sup>11</sup> N/m<sup>2</sup> ). This yields the following numerical relation for the stress derived from the change in the 'c' lattice parameter:

$$
\sigma(\text{Nm}^{-2}) = -4.5 \times 10^{11} e\_{zz} \tag{11}
$$

The calculated values of stress (*σ*) in the undoped and doped films are listed in Table 4. The negative sign indicates that the films are in a state of compressive stress which not seen in pure ZO film. The total stress in the film typically consists of two components: intrinsic stress which introduced by impurities, defects and lattice distortions in the crystal, and the extrinsic stress familiarized by the lattice mismatch and thermal expansion coefficient mismatch between the film and substrate. When the thickness of a thin film is large, the extrinsic stress in the thin film normally relaxes. All of the films in our current study have a thickness around 150 nanometers. As a result, extrinsic stress will be absent, and the overall predicted stress values will be dominated by intrinsic stress. Furthermore, in the current scenario of sputtering deposition with rising DC power, this intrinsic stress arises as a result of the bombardment of energetic species. Strong compressive stress is present in films deposited at lower temperatures. In the sputtering process, the highly energetic species could be implanted below the ZO surface, causing high intrinsic stresses by creating the defects [29,30]. The existence of stress in the ZO films was reported by Gupta and Mansingh, who ascribed it to oxygen interstitial defects [31]. It is well known that, because of the higher formation energy, the oxygen interstitials are less expected in ZO films. In the

present experiment all the films were deposited in argon atmosphere. So, the possibility of the formation of oxygen interstitials is ruled out and the formation of Zn interstitials and the oxygen vacancies are expected. Thus, in our case, the intrinsic stress in the deposited ZO and ZOFe films seems to arise due to presence of Zn interstitials.

### *3.3. Microstructural Analysis*

Figure 5 shows the morphology of the ZO and Fe doped ZO thin films; the pure zinc oxide film surface seems to be agglomerated, spherical in shape as shown Figure 5a. The surface of the films consists of tightly packed grains forming a smooth surface without any voids and cracks whereas the Fe doped zinc oxide films in Figure 5b,c look like black and white islands of spherical shape. Still, by increasing the DC deposited power for Fe with ZO deposition, the black spots decreased. A SEM image of the sectional view of the samples could help in discovering the shape of the spherical grains and their interface with the substrate as in Figure 5. Moreover, the shape of the grains would be hard to identify because the growth could not be described in a reliable way. Further, the distribution and size in grains is inhomogeneous. Figure 6 shows the EDX spectrum of the ZO and Fe doped ZO thin films. The observed peaks are attributed to zinc, oxygen and iron as principal elemental components with no other impurities being obtained. The Zn:O:Fe atomic percentage ratios of the ZO thin film are listed. *Ceramics* **2022**, *5,* FOR PEER REVIEW 11

**Figure 5.** Cross sectional view of ZO thin film deposited in SiO2 and SEM photos for the (**a**) undoped ZO and Fe doped ZO at DC deposited power (**b**) 3 W and (**c**) 5 W. **Figure 5.** Cross sectional view of ZO thin film deposited in SiO<sup>2</sup> and SEM photos for the (**a**) undoped ZO and Fe doped ZO at DC deposited power (**b**) 3 W and (**c**) 5 W.

**Figure 6.** EDX spectrum of (**a**) undoped ZO and Fe doped ZO at DC deposited power (**b**) 3 W and (**c**) 10 W. **Figure 6.** EDX spectrum of (**a**) undoped ZO and Fe doped ZO at DC deposited power (**b**) 3 W and (**c**) 10 W.

#### *3.4. Optical Properties 3.4. Optical Properties*

#### 3.4.1. Transmittance Spectra 3.4.1. Transmittance Spectra

Zinc oxide and iron doped zinc oxide samples grown on SiO2 are transparent to visible and NIR light; thus it was possible to obtain transmittance spectra and make observations on the material properties. The spectra were taken in the range between 300 and 900 nm, as shown in Figure 7. The optical transmittance spectra of the films show a good transmission in the visible region and a sharp fall in the UV region which corresponds to the band gap. The decrease of the transmittance is due to the interaction of the incident long-wavelength radiation with the free electrons in the films [32]. Moreover, the films have good transmission in the visible region of the spectrum terminated at shorter wavelengths while, in the UV region, the obvious absorption could be observed; it increased with increasing Fe content. Zinc oxide and iron doped zinc oxide samples grown on SiO<sup>2</sup> are transparent to visible and NIR light; thus it was possible to obtain transmittance spectra and make observations on the material properties. The spectra were taken in the range between 300 and 900 nm, as shown in Figure 7. The optical transmittance spectra of the films show a good transmission in the visible region and a sharp fall in the UV region which corresponds to the band gap. The decrease of the transmittance is due to the interaction of the incident longwavelength radiation with the free electrons in the films [32]. Moreover, the films have good transmission in the visible region of the spectrum terminated at shorter wavelengths while, in the UV region, the obvious absorption could be observed; it increased with increasing Fe content.

**Figure 7.** Transmittance of ZO and ZOFe thin films. **Figure 7.** Transmittance of ZO and ZOFe thin films.

#### 3.4.2. Optical Band Gap 3.4.2. Optical Band Gap

The energy gap and refractive index of semiconductors are two essential physical parameters that define their optical and electronic properties. The optical energy gap (*Eg*) determines the threshold for absorption of photons in a semiconductor. In order to determine the optical energy gap (*Eg*), there are several numerical methods based on the experimental data. We are used different methods and compared them to obtain (*Eg*) values which have a good agreement with the literature. The energy gap and refractive index of semiconductors are two essential physical parameters that define their optical and electronic properties. The optical energy gap (*Eg*) determines the threshold for absorption of photons in a semiconductor. In order to determine the optical energy gap (*Eg*), there are several numerical methods based on the experimental data. We are used different methods and compared them to obtain (*Eg*) values which have a good agreement with the literature.

#### (A) Tauc method: (A) Tauc method:

In most of the literature, the energy gap in the crystalline and amorphous semiconductors is determined by the Tauc method [33]. This method is easy to apply due to fitting of the linear function to the graph for parameter: m = {0.5, 1.5, 2, 3} which represents direct allowed, direct forbidden, indirect allowed and indirect forbidden optical transition, respectively. The energy gap determined using the Tauc method is often subject to high uncertainty. In some cases, it is possible to fit a straight line for all values of parameter m. The energy gap of all the films is determined from the absorption coefficient (α) which can be calculated from the transmittance (*T*) of the undoped and doped ZO thin films. The absorption coefficient (α) is calculated from Equation: In most of the literature, the energy gap in the crystalline and amorphous semiconductors is determined by the Tauc method [33]. This method is easy to apply due to fitting of the linear function to the graph for parameter: m = {0.5, 1.5, 2, 3} which represents direct allowed, direct forbidden, indirect allowed and indirect forbidden optical transition, respectively. The energy gap determined using the Tauc method is often subject to high uncertainty. In some cases, it is possible to fit a straight line for all values of parameter m. The energy gap of all the films is determined from the absorption coefficient (α) which can be calculated from the transmittance (*T*) of the undoped and doped ZO thin films. The absorption coefficient (α) is calculated from Equation:

$$a = \left(\frac{1}{d}\right) \ln\left(\frac{1}{T}\right) \tag{12}$$

where (*d*) is film thickness and (*T*) is the transmittance. The absorption edge was analyzed by the following Equation: where (*d*) is film thickness and (*T*) is the transmittance. The absorption edge was analyzed by the following Equation:

$$\left(\omega h\nu\right)^{2} = A\left(h\nu - E\_{\mathcal{S}}\right) \tag{13}$$

where *A* is a constant, *h* is the incident photon energy and *Eg* is the optical energy band gap. Based on Equation (15), the plots of (*αhν*)2 as a function of incident photon energy (*hν*) were obtained for the undoped and doped ZO thin films and are shown in Figure 8. The linear portions of these plots are extrapolated to meet the energy axis and the energy value at (*αhν*)2 = 0 gives the values of the energy gap *Eg* for the thin films. The *Eg* values are listed in Table 5. where *A* is a constant, *hν* is the incident photon energy and *E<sup>g</sup>* is the optical energy band gap. Based on Equation (15), the plots of (*αhν*) <sup>2</sup> as a function of incident photon energy (*hν*) were obtained for the undoped and doped ZO thin films and are shown in Figure 8. The linear portions of these plots are extrapolated to meet the energy axis and the energy value at (*αhν*) <sup>2</sup> = 0 gives the values of the energy gap *E<sup>g</sup>* for the thin films. The *E<sup>g</sup>* values are listed in Table 5.

**Figure 8.** Plot of (*αhν*) <sup>2</sup> vs. (*hν*) for the undoped and Fe doped ZO thin films.

**Figure 8.** Plot of (*αhν*)2 vs. (*hν*) for the undoped and Fe doped ZO thin films. (B)The derivative of transmittance (T) against photon energy (*hν*) **Table 5.** Energy gap values *E<sup>g</sup>* obtained from different methods for undoped ZO and Fe doped ZO thin films.


(B) The derivative of transmittance (T) against photon energy (*hν*)

This method is applicable for transitions between direct band gaps [34]. The transmission through a film may be approximated as:

$$T(E) \approx \left[1 - \mathcal{R}(E)\right]^2 e^{\mathcal{Q}(E)d} \tag{14}$$

where *E* is the energy of the incident light, *R* is the reflectance, and *d* is film thickness. where *C* is a constant and η*r*(*E*) is the energy-dependent index of refraction. The behavior of the reflectance at the band gap energy must be considered while calculating the

α() =

*Ceramics* **2022**, *5,* FOR PEER REVIEW 15

$$\mathfrak{a}(E) = \frac{\mathbb{C}}{E\mathfrak{n}\_r(E)}\sqrt{E - E\_\mathcal{g}} \tag{15}$$

ா(ா) ඥ− (15)

where *C* is a constant and η*r*(*E*) is the energy-dependent index of refraction. The behavior of the reflectance at the band gap energy must be considered while calculating the transmission derivative. It can be shown that as *E*→*Eg*, the quantities (1 − *R(E)*) and *dR(E)*/*dE* remain finite at energies in the environs of *E<sup>g</sup>* due to the presence of a band tail that prohibits singularity- type behavior. The reflectance and its derivative are wellbehaved values near the band gap energy, according to experimental studies concerning ZO thin film properties [35]. As a result, we can write the transmission through a direct gap semiconductor as that prohibits singularity- type behavior. The reflectance and its derivative are well-behaved values near the band gap energy, according to experimental studies concerning ZO thin film properties [35]. As a result, we can write the transmission through a direct gap semiconductor as () = ௗ൫ଵ/ா(ா)൯ඥாିா (16) After taking the first derivative of *T(E)* with respect to energy and then at the limit *E*→*Eg*, this results in a spike towards negative infinity at *E* = *Eg*. Due to band tail states, it

$$T(E) = e^{\mathbb{C}d(1/E\eta\_r(E))\sqrt{E - E\_\mathcal{S}}} \tag{16}$$

After taking the first derivative of *T(E)* with respect to energy and then at the limit *E*→*Eg*, this results in a spike towards negative infinity at *E* = *Eg*. Due to band tail states, it was assumed that *d*η*r*(*E*)/*dE* is continuous around *Eg*. As a result, when plotting *dT*/*dE* versus *E*, a significant singularity will appear at the band gap energy. In realistic instances, absorption tails exist, which soften the divergence and create well-defined peaks around the gap energy, as also can be seen in Figure 9 for ZO and ZOFe/10 W thin films. The values of the band gap of all investigated samples obtained by the derivative method are listed also in Table 5. These values of band gap are very close to the values obtained by Tauc method. Discrepancy between values of band gap obtained by both methods can be explained: the fitting of a line in the linear region of (*αhν*) 2 is quite challenging which adds errors in determination of the band gap. Hence, the derivative method is considered more appropriate when compared to Tauc method. versus *E*, a significant singularity will appear at the band gap energy. In realistic instances, absorption tails exist, which soften the divergence and create well-defined peaks around the gap energy, as also can be seen in Figure 9 for ZO and ZOFe/10 W thin films. The values of the band gap of all investigated samples obtained by the derivative method are listed also in Table 5. These values of band gap are very close to the values obtained by Tauc method. Discrepancy between values of band gap obtained by both methods can be explained: the fitting of a line in the linear region of (*αhν*)2 is quite challenging which adds errors in determination of the band gap. Hence, the derivative method is considered more appropriate when compared to Tauc method. lim ா→ ா ௗ் ௗா = − ቀ−0 − 0 + <sup>ଵ</sup> ቁ → − ∞ (17)

$$\lim\_{E \to E\_{\mathcal{K}}} \frac{dT}{dE} = -\mathcal{C} \left( -0 - 0 + \frac{1}{0} \right) \to -\infty \tag{17}$$

**Figure 9.** The plots of the derivative of the transmission with respect to energy for the undoped and ZOFe/10 W thin films. **Figure 9.** The plots of the derivative of the transmission with respect to energy for the undoped and ZOFe/10 W thin films.

(C) The derivative of transmittance (T) against wavelength (λ)

#### (C) The derivative of transmittance (T) against wavelength (λ) As we noticed in the transmittance curve of ZO films, there are transmittance tails at the high energy side of the absorption edge. Because a sharp absorption edge is generally

*Ceramics* **2022**, *5,* FOR PEER REVIEW 16

As we noticed in the transmittance curve of ZO films, there are transmittance tails at the high energy side of the absorption edge. Because a sharp absorption edge is generally observed in the transmittance spectra of direct band gap semiconductor films, we computed the derivative of transmittance (T) against wavelength (λ) in an effort to accurately determine the optical band gap of undoped and doped ZO films. The sharp absorption edge in the transmittance spectra may result in a sharp peak in a plot of dT/dλ vs. *hν* from which the optical band gap can be successfully determined. Figure 10 shows plots of dT/dλ vs. *hν* for the ZO and ZOFe/10 W thin films. One can see in Figure 10 a broad peak centered around 3.3 eV, corresponding to a band gap of all investigated films, also listed in Table 5. observed in the transmittance spectra of direct band gap semiconductor films, we computed the derivative of transmittance (T) against wavelength (λ) in an effort to accurately determine the optical band gap of undoped and doped ZO films. The sharp absorption edge in the transmittance spectra may result in a sharp peak in a plot of dT/dλ vs. *hν* from which the optical band gap can be successfully determined. Figure 10 shows plots of dT/dλ vs. *hν* for the ZO and ZOFe/10 W thin films. One can see in Figure 10 a broad peak centered around 3.3 eV, corresponding to a band gap of all investigated films, also listed in Table 5.

**Figure 10.** dT/dλ vs. photon energy for the undoped and undoped and ZOFe/10 W thin films. **Figure 10.** dT/dλ vs. photon energy for the undoped and undoped and ZOFe/10 W thin films.

(D) Logarithmic derivative (LD) method:

(D) Logarithmic derivative (LD) method: This approach combines the advantages of both the Tauc [33] and McLean [36] methods: the fit is linear, and the parameter m is calculated as a consequence of the fit-This approach combines the advantages of both the Tauc [33] and McLean [36] methods: the fit is linear, and the parameter m is calculated as a consequence of the fitting. Significantly, the results of all three approaches are consistent in simple cases.

ting. Significantly, the results of all three approaches are consistent in simple cases. As we referred above, the Tauc method is based on linear fit for *αhν*1/m data for selected *m* = {0.5, 1.5, 2, 3} values. The zero point of the fitted linear function is then the energy gap *Eg*. In the McLean method, one can determine m, *Eg* and parameter A in Equation (13) by fitting a power function to *αhν* data. LD transformation starts from the Tauc equation. For the sake of calculating the natural logarithm, let us suppose that all the As we referred above, the Tauc method is based on linear fit for *αhν* 1/m data for selected *m* = {0.5, 1.5, 2, 3} values. The zero point of the fitted linear function is then the energy gap *Eg*. In the McLean method, one can determine m, *E<sup>g</sup>* and parameter A in Equation (13) by fitting a power function to *αhν* data. LD transformation starts from the Tauc equation. For the sake of calculating the natural logarithm, let us suppose that all the quantities in Equation (13) are unitless. Then the natural logarithm will be:

$$
\ln(\text{ad}\mathbf{v}) = m\ln(A) + m\ln\left(\text{hv} - E\_{\mathcal{S}}\right) \tag{18}
$$

(αℎν) =  () +  ൫ℎν − ൯ (18) By differentiation of Equation (18) with respect to *hν*

$$\frac{d\ln(\text{alv})}{d(h\text{v})} = m \left(\frac{1}{h\text{v} - \text{E}\_{\text{g}}}\right) \tag{19}$$

ௗ(ఔ)

ିா The left side of Equation (19) is calculated based on the experimental data. This numerical derivative should be calculated as the difference of transformed measurement The left side of Equation (19) is calculated based on the experimental data. This numerical derivative should be calculated as the difference of transformed measurement data (*hν*, *αhν*) from subsequent measurement steps.

Now, we can determine the precise value of the optical band gap from the ௗ (ఈఔ)

data (*hν*, *αhν*) from subsequent measurement steps.

Now, we can determine the precise value of the optical band gap from the *d ln*(*αhν*) *d*(*hν*) vs. *hν* curve as in Figure 11. We can see that a peak in the curve appeared at the band gap energy (*Eg*); i.e., *hν* = *Eg*. The peak at a particular energy value gives the approximate optical band gap, *E<sup>g</sup>* [37]. energy (*Eg*); i.e., *hν* = *Eg*. The peak at a particular energy value gives the approximate optical band gap, *Eg* [37]. *Ceramics* **2022**, *5,* FOR PEER REVIEW 17 energy (*Eg*); i.e., *hν* = *Eg*. The peak at a particular energy value gives the approximate op-

*Ceramics* **2022**, *5,* FOR PEER REVIEW 17

**Figure 11.** LD of *αhν* as a function of photon energy *hν* with fitted straight lines for the undoped and Fe/10 w doped ZO thin films. **Figure 11.** LD of *αhν* as a function of photon energy *hν* with fitted straight lines for the undoped and Fe/10 w doped ZO thin films. and Fe/10 w doped ZO thin films. On other hand, the type of the optical transition can be estimated from the slope (m)

On other hand, the type of the optical transition can be estimated from the slope (m) of the curve of ln (*αhν*) vs. ln (*hν* − *Eg*) using the *Eg* value of the highest Fe concentration doped ZO (15 at. %), as an example, and to determine the value of (m) which was found about 1/2 from the slope as shown in Figure 12. This indicates that the optical transition through the base glass and also the doped samples is a directly allowed transition which On other hand, the type of the optical transition can be estimated from the slope (m) of the curve of ln (*αhν*) vs. ln (*hν* − *Eg*) using the *E<sup>g</sup>* value of the highest Fe concentration doped ZO (15 at. %), as an example, and to determine the value of (m) which was found about 1/2 from the slope as shown in Figure 12. This indicates that the optical transition through the base glass and also the doped samples is a directly allowed transition which we selected to obtain *E<sup>g</sup>* values by using the Tauc relation. of the curve of ln (*αhν*) vs. ln (*hν* − *Eg*) using the *Eg* value of the highest Fe concentration doped ZO (15 at. %), as an example, and to determine the value of (m) which was found about 1/2 from the slope as shown in Figure 12. This indicates that the optical transition through the base glass and also the doped samples is a directly allowed transition which we selected to obtain *Eg* values by using the Tauc relation. Thus, the estimated values of *Eg* for ZO doped and undoped by Fe are listed in Table 5 to compare with the other values obtained from different approaches.

**Figure 12.** Type of the optical transition used in the undoped and Fe doped ZO thin films.

Thus, the estimated values of *E<sup>g</sup>* for ZO doped and undoped by Fe are listed in Table 5 to compare with the other values obtained from different approaches.

One can see that in Table 5, the obtained values of *E<sup>g</sup>* determined from the Tauc method are close to those determined by dT/dE and dT/dλ and LD methods. Furthermore, the average values of the *E<sup>g</sup>* are 3.26, 3.27, 3.29, 3.33 and 3.39 eV for the undoped and 3 W, 5 W, 10 W and 15 W of Fe-doped ZO thin films, respectively.

It is worthwhile to notice that two points may be inferred from these energy gap values. The first is that from the undoped to doped thin film, the *E<sup>g</sup>* values increased. Furthermore, as the doping concentration of Fe increases, the *E<sup>g</sup>* value increases.

The first point can be interpreted as follows. As we know, the incidence of a photon with energy of *hν* on a semiconductor leads to a transition between the highest occupied state of valence band and the lowest unoccupied state of the conduction band. This phenomenon is applied to calculate the optical band gap energy of ZO thin films as we referred before. In Fe-doped ZO thin films, the optical band gap energy mainly depends on the valence state of Fe ions [38,39]. With Fe-doped ZO thin films, there is a large variability in the optical behavior which leads to inconsistent conclusions. The valence state of Fe is an important factor in the optical properties of ZO thin films. Fe dopants in ZO thin films can be in the two forms of Fe2+ or Fe3+, or even coexistence of these two forms is seen. By addition of Fe into ZO films, the maximum of the valence band is increased and the minimum of the conduction band is decreased which leads to the reduction of band gap. However, substituting Fe3+ ions into Zn2+ leads to providing extra free carrier concentrations. As a result, the Fermi level moves toward the conduction band and the band gap increases. This makes the optical band gap of Fe-doped ZO thin films tunable.

Moreover, it is known that doping significantly alters the absorption spectrum of the pure semiconductor, resulting in degenerate energy levels which cause the Fermi level to push above the conduction band edge. This shift is known as doping induced band-filling called a Moss–Burstein shift and leads to an increase in the band gap from the undoped to doped thin film [40]. In addition, with Fe incorporation, there are initially sp-d exchange interactions that do not strongly affect the band gap. After a further doping percentage, the Moss–Burstein effect dominates the normal sp—d exchange interactions. In addition, ZO is normally n-type material stable at room temperature in which the Fermi level lies in the conduction band when doped with greater Fe [41]; as a result, there is an increase in the band gap with more Fe doping. Since no empty state is available inside the conduction band, so the absorption edge shifts towards the higher energy region, thereby giving rise to an increase in the band gap. With increasing Fe doping the d—d transition of Fe ions also increases and leads to increase in the band gap. According to Parra-Palomino et al. [38], increasing Fe doping concentration to ZO nanocrystals increased the band gap significantly. Chen et al. [42] and Wang et al. [43] found that raising the Fe doping concentration in ZO films can minimize the band gap.

#### 3.4.3. Urbach Energy

One of the main properties which measure the degree of structural disorder in the films is the Urbach energy (*EU*). The imperfection in the structurally disordered film leads to broadening the bands of localized states. As we mentioned, near the band edge, the absorption coefficient has an exponential dependence on photon energy. This dependence is given as follows:

$$\alpha = \alpha\_o exp\left(\frac{h\mathbf{v}}{E\_{ll}}\right) \tag{20}$$

where α*<sup>o</sup>* is the band tailing parameter that can be obtained by:

$$\alpha\_o = \sqrt{\frac{\sigma\_o \left(\frac{4\pi}{c}\right)}{x \Delta E}}\tag{21}$$

where σ*o*, *c* and ∆*E* are the electrical conductivity at absolute zero, the velocity of light and the width of the tail of localized states in the energy gap, respectively, and x = 0.5 for the directly allowed transitions [44]. Then the Urbach energy can emerge from Equation (22) as: where σ*o*, *c* and Δ*E* are the electrical conductivity at absolute zero, the velocity of light and the width of the tail of localized states in the energy gap, respectively, and x = 0.5 for the directly allowed transitions [44]. Then the Urbach energy can emerge from Equation (22) as:

absorption coefficient has an exponential dependence on photon energy. This depend-

α=α ቀ

α <sup>=</sup> ටቀ

 ாೆ

రಘ ቁ

*Ceramics* **2022**, *5,* FOR PEER REVIEW 19

where α*o* is the band tailing parameter that can be obtained by:

ence is given as follows:

$$E\_{II} = \left[\frac{d\left(\ln(\alpha)\right)}{d\,h\mathbf{v}}\right]^{-1} \tag{22}$$

ቁ (20)

௫ா (21)

As shown in Figure 13, the values of the *E<sup>U</sup>* were obtained from the slopes of the linear fitting of the plots of ln (α) versus (*hν*) and are listed in Table 6. E<sup>U</sup> increases from 225 to 403 meV as the Fe concentration increased. The dopant may contribute significantly to the width of localized states within the ZO's optical band. In other words, the *E<sup>U</sup>* is responsible for the valence and conduction bands' tails [45]. Moreover, *E<sup>U</sup>* values confirm the GIXRD analysis which shows that Fe doping is followed by an increase of the strain which can also create structural disorder. As shown in Figure 13, the values of the *EU* were obtained from the slopes of the linear fitting of the plots of ln (α) versus (*hν*) and are listed in Table 6. EU increases from 225 to 403 meV as the Fe concentration increased. The dopant may contribute significantly to the width of localized states within the ZO's optical band. In other words, the *EU* is responsible for the valence and conduction bands' tails [45]. Moreover, *EU* values confirm the GIXRD analysis which shows that Fe doping is followed by an increase of the strain which can also create structural disorder.

**Figure 13.** ln (α) vs. (*hν*) to determine the Urbach energy of ZO and Fe doped ZO thin films. **Figure 13.** ln (α) vs. (*hν*) to determine the Urbach energy of ZO and Fe doped ZO thin films.

3.4.4. Refractive Index **Table 6.** Optical band gap and refractive index values for undoped ZO and Fe doped ZO thin films.


#### 3.4.4. Refractive Index

The two most interesting optical properties of a semiconductor are the absorption edge, or optical energy gap, and the refractive index. The relation between the optical band gap *E<sup>g</sup>* and the refractive index n<sup>o</sup> was proposed for first time by Moss in 1950 on the very general grounds that all energy levels in a solid are scaled down by a factor 1/ε 2 *opt*, where ε*opt* = n<sup>2</sup> is the optical dielectric constant; it was described in detail elsewhere [46,47]. Moss showed that the energy levels of the semiconducting materials can be scaled down by the term 1/n<sup>4</sup> as follows [46]:

$$n\_o^4 E\_\mathcal{S} = 9 \text{ eV} \tag{23}$$

Then Ravindra and Srivastava [48] modified Equation (25) to be:

$$m\_o^4 E\_\% = 108\,\text{eV} \tag{24}$$

The main difference between the Moss and Ravindra relations originated from the estimating of the reflection loss by Ravindra. Moreover, Equation (26) covers most of the semiconducting materials used for solar cell, photocatalytic and sensing applications. Dimitrov and Sakka [49] have deduced the correlation between the refractive index and the optical band gap from the Lorentz- Lorentz (L.L.) Equation:

$$\frac{n\_o^2 - 1}{n\_o^2 + 2} = 1 - \sqrt{\frac{E\_g}{20}}\tag{25}$$

The estimated refractive index based on the optical band gap is shown in Table 6. The refractive index of ZO was compared with the optical data from references [22,49]. As we can see, by increasing the Fe content, the refractive index decreased and as we know the refractive index is correlated to the optical band gap. Subsequently, from Equations (25) and (26), the relation between them is inverse.

#### 3.4.5. The Extinction Coefficient

The extinction coefficient (*k*) reflects the absorption of electromagnetic waves in the semiconductor due to inelastic scattering events and it can be calculated from the transmittance spectra using the relation

$$k = \frac{\alpha \lambda}{4\pi} \tag{26}$$

The dependence of *k* on the wavelength for all investigated thin films is shown in Figure 14. In the UV region (*λ* < 400 nm), the *k* value increases from 0.36 to 0.42 with increasing Fe concentration. Furthermore, the *k* values are very small and relatively constant in most of the visible region. This indicates the smoothness on the surface and homogeneity of the films and enhances its uses in blocking UV radiation [40]. As well, increase in *k* value by doping Fe can be ascribed to increasing topical donor levels formed within the energy levels that lead to an increased extinction coefficient, showing that the electronic transitions occur directly. *Ceramics* **2022**, *5,* FOR PEER REVIEW 21

**Figure 14.** The extinction coefficient (K) as a function of the wavelength (λ) for the undoped and Fe doped ZO films. **Figure 15.** M-H curve showing ferromagnetic loops for the undoped and Fe doped ZO thin films. **Figure 14.** The extinction coefficient (K) as a function of the wavelength (λ) for the undoped and Fe doped ZO films.

temperature ferromagnetism (RTFM) for all the doping concentrations of Fe; they were obtained from room temperature vibrating sample magnetometer (VSM) measurements. It is clear from Figure 15 that undoped ZO and Fe doped ZO films show well-defined hysteresis loops with noticeable coercivity. The loops are S-shaped, which indicates the ferromagnetic nature of the samples. One can notice that there is non-saturation of the M-H loop at higher applied fields, which is a commonly observed phenomenon in semiconductor/metal nanocomposites. Rumpf et al. also observed such non-saturation behavior for ferromagnetic semiconductors and attributed it to the lattice defects and

*3.5. Magnetic Properties* 

magnetic anisotropy [50].

3.5.1. Magnetic Behavior of the Films

doped ZO films.

*3.5. Magnetic Properties* 

#### *3.5. Magnetic Properties* 3.5.1. Magnetic Behavior of the Films

3.5.1. Magnetic Behavior of the Films The hysteresis loops of ZO and Fe doped ZO thin films in Figure 15 show clear room

**Figure 14.** The extinction coefficient (K) as a function of the wavelength (λ) for the undoped and Fe

*Ceramics* **2022**, *5,* FOR PEER REVIEW 21

The hysteresis loops of ZO and Fe doped ZO thin films in Figure 15 show clear room temperature ferromagnetism (RTFM) for all the doping concentrations of Fe; they were obtained from room temperature vibrating sample magnetometer (VSM) measurements. It is clear from Figure 15 that undoped ZO and Fe doped ZO films show well-defined hysteresis loops with noticeable coercivity. The loops are S-shaped, which indicates the ferromagnetic nature of the samples. One can notice that there is non-saturation of the M-H loop at higher applied fields, which is a commonly observed phenomenon in semiconductor/metal nanocomposites. Rumpf et al. also observed such non-saturation behavior for ferromagnetic semiconductors and attributed it to the lattice defects and magnetic anisotropy [50]. temperature ferromagnetism (RTFM) for all the doping concentrations of Fe; they were obtained from room temperature vibrating sample magnetometer (VSM) measurements. It is clear from Figure 15 that undoped ZO and Fe doped ZO films show well-defined hysteresis loops with noticeable coercivity. The loops are S-shaped, which indicates the ferromagnetic nature of the samples. One can notice that there is non-saturation of the M-H loop at higher applied fields, which is a commonly observed phenomenon in semiconductor/metal nanocomposites. Rumpf et al. also observed such non-saturation behavior for ferromagnetic semiconductors and attributed it to the lattice defects and magnetic anisotropy [50].

**Figure 15.** M-H curve showing ferromagnetic loops for the undoped and Fe doped ZO thin films. **Figure 15.** M-H curve showing ferromagnetic loops for the undoped and Fe doped ZO thin films.

The mechanism of RTFM in the ZO based system was still unclear for researchers. The dopant related secondary phase, magnetic cluster, metal precipitation formation will lead to magnetic behavior in the doped ZO. In a doped ZO system, the carrier exchange interactions such as RKKY, indirect double exchange interactions and super exchange interactions induced a magnetic nature [51–53]. In addition, the doped and co-doped ZO samples reveal ferromagnetism at RT due to the presence of intrinsic lattice defects. The presence of lattice defects as strain and texture, and concentration of Fe doping ions may also have a major impact on the room temperature magnetic properties [54,55]. Generally, RT ferromagnetism in TM doped ZO films could be explained on the basis of intrinsic and extrinsic effects associated with the local magnetic moments. If the spin of Fe ions has an exchange interaction with local magnetic moments of ZO then it is intrinsic effect. When there is a direct interaction between the local magnetic moments of the Fe clusters, it is known to be an extrinsic effect.

Accordingly, in Figure 15, the M-H curves of the pure ZO exhibit ferromagnetic (FM) behavior and Fe doped ZO thin films, with Fe concentration 2.4, 4.8 and 5 at. % of Fe dopants, exhibit diamagnetic behavior at a high magnetic field and FM behavior in the lowest field regime. There are several literature reports that pure ZO exhibited ferromagnetic nature at room temperature due to the exchange interaction between localized electron spin moments resulting from oxygen vacancies [56].

For more details about the diamagnetic behavior shown in the pure ZO films, a fitting with Langevin curve was performed, as shown in Figure 16; this type of analysis has been conducted on ZnO doped with transition metals [9]. ting with Langevin curve was performed, as shown in Figure 16; this type of analysis has been conducted on ZnO doped with transition metals [9].

For more details about the diamagnetic behavior shown in the pure ZO films, a fit-

Accordingly, in Figure 15, the M-H curves of the pure ZO exhibit ferromagnetic (FM) behavior and Fe doped ZO thin films, with Fe concentration 2.4, 4.8 and 5 at. % of Fe dopants, exhibit diamagnetic behavior at a high magnetic field and FM behavior in the lowest field regime. There are several literature reports that pure ZO exhibited ferromagnetic nature at room temperature due to the exchange interaction between localized

The mechanism of RTFM in the ZO based system was still unclear for researchers. The dopant related secondary phase, magnetic cluster, metal precipitation formation will lead to magnetic behavior in the doped ZO. In a doped ZO system, the carrier exchange interactions such as RKKY, indirect double exchange interactions and super exchange interactions induced a magnetic nature [51–53]. In addition, the doped and co-doped ZO samples reveal ferromagnetism at RT due to the presence of intrinsic lattice defects. The presence of lattice defects as strain and texture, and concentration of Fe doping ions may also have a major impact on the room temperature magnetic properties [54,55]. Generally, RT ferromagnetism in TM doped ZO films could be explained on the basis of intrinsic and extrinsic effects associated with the local magnetic moments. If the spin of Fe ions has an exchange interaction with local magnetic moments of ZO then it is intrinsic effect. When there is a direct interaction between the local magnetic moments of the Fe

*Ceramics* **2022**, *5,* FOR PEER REVIEW 22

clusters, it is known to be an extrinsic effect.

electron spin moments resulting from oxygen vacancies [56].

**Figure 16.** Fitting of the Langevin function with experimental M (H) curves. **Figure 16.** Fitting of the Langevin function with experimental M (H) curves.

Further, the ferromagnetic property for DMSs can be contributed from the secondary phase related to metal oxides or metal clusters with a ferromagnetic property. DMSs could have an intrinsic effect as well [57]. However, our GIXRD patterns clearly indicate the absence of clustering of metallic iron phase or iron oxide cluster in the films. Moreover, the possibility of attributing the observed ferromagnetic order of Fe doped ZO thin films to the carrier-mediated mechanism or RKKY exchange interaction is also ruled out because the DMSs usually show very high resistance. Therefore, RTFM of our Fe doped ZO films are not due to the carrier-mediated mechanism. A pertinent and appropriate explanation for the observed RTFM [58] in our case is the ferromagnetic exchange Further, the ferromagnetic property for DMSs can be contributed from the secondary phase related to metal oxides or metal clusters with a ferromagnetic property. DMSs could have an intrinsic effect as well [57]. However, our GIXRD patterns clearly indicate the absence of clustering of metallic iron phase or iron oxide cluster in the films. Moreover, the possibility of attributing the observed ferromagnetic order of Fe doped ZO thin films to the carrier-mediated mechanism or RKKY exchange interaction is also ruled out because the DMSs usually show very high resistance. Therefore, RTFM of our Fe doped ZO films are not due to the carrier-mediated mechanism. A pertinent and appropriate explanation for the observed RTFM [58] in our case is the ferromagnetic exchange mechanism involving oxygen vacancies model. This could be explained in terms of two reasons: (1) the number of oxygen vacancies (OV) and zinc interstitials (Zn<sup>i</sup> ) and (2) the exchange interaction between doped transition metal ion and O ion spins.

Further, Kittilstved et al. [59] stated that the ferromagnetism in TM doped ZO is due to the effective dopant defect hybridization arising from the energetic alignment of the shallow donors, Zn interstitial (Zn<sup>i</sup> ) and oxygen vacancy (OV), with TM+/2+ level. Thus, the magnetic coupling between Fe ions with ZO is FM mediated by (Zn<sup>i</sup> ) and (OV) and this may account for the observed RTFM and indicates that defects play a role to gain a FM coupling. In addition, in our case, for low magnetic fields, the ferromagnetic behavior can be attributed to the presence of small magnetic dipoles located at the surface of thin film nanocrystals, which interact with their nearest neighbors inside the crystal involving oxygen vacancies [53].

#### 3.5.2. The Magnetic Parameters

Table <sup>7</sup> shows the saturation magnetization M<sup>s</sup> of undoped ZO~9 <sup>×</sup> <sup>10</sup>−<sup>4</sup> emu/g at 2000 G. Vettumperumal et al. [60] have observed RTFM of ZO thin film with M<sup>s</sup> value ~ 6.9 <sup>×</sup> <sup>10</sup>−<sup>5</sup> emu/g emu/g at 300 K while Ariyakkani et al. [61] have obtained M<sup>s</sup> value ~ 5.6 <sup>×</sup> <sup>10</sup>−<sup>5</sup> emu/g at 300 K for ZO thin film. Therefore, the magnetization value obtained in our work is higher as compared to previous reported values of undoped ZO. Furthermore, it is obvious, from Table 7, that the saturation magnetization of the films is increased by increasing the Fe concentration. As we mentioned before, the origin of ferromagnetism in transition metal doped ZO until now is still not clear. A number of studies indicate that ferromagnetism in transition metal doped ZO may come from secondary magnetic phases [62]. However, others reported that the RTFM in Fe-doped ZO could be due to the oxidation of Fe in mixed valence (Fe2+ and Fe3+) states [63]. The exchange interaction between conductive electrons and local spin polarized electrons is responsible for the magnetic behavior of Fe-doped ZO. Likewise, the occurrence of Fe-Fe super exchange interaction could be one reasons for increased M<sup>s</sup> values with increasing Fe concentration; it is found to be dominant at higher Fe doping concentrations [64,65].


**Table 7.** Magnetization, coercivity and remanence values of undoped and Fe doped ZO thin films.

The coercivity of Fe doped films is lower than the undoped film at any particular Fe concentration, as shown in Table 7. There are several factors influencing the reduction in the coercivity such as magnetic domain size, micro-strain, magnetocrystalline anisotropy, and shape anisotropy. However, such decrease in coercivity could be related to a decrease in the anisotropy field, which reduces the thickness of the domain wall [66]. However, the coercivity (Hc) is low (~95 Oe) for the doped films and it indicates that all Fe doped ZO thin film samples are ferromagnetically soft compared to ZO. The magnetization measurements by Kumar et al. [63] on 1% Fe doped ZO samples shows soft ferromagnetic behavior at room temperature with coercivity (Hc) about 27 Oe. Further, low coercivity and high saturation magnetization are required for spintronics memory devices and this is obtained for the ZO: Fe (5% at 15 W) thin films compared to the other thin films. It can be seen that the coercivities and saturation magnetizations were enhanced with the increase of Fe content. This phenomenon should be investigated further to understand it. Therefore, to enhance ferromagnetic order, the appropriate doping concentration of Fe in ZO wurtzite structure plays a vital role in obtaining an optimum number of charge carriers and oxygen vacancies, which may be useful for spintronics based solid state devices.

Thus, the remanence (Mr) is found to slightly increase with the increase in the high Fe concentration which is related to the saturation magnetization Ms. Table 7 shows the magnetic saturation (Ms), remanence magnetization (Mr) and coercivity (Hc) for all the thin films.

#### *3.6. Electrical Properties*

As already noted, ZO is a large band gap semiconductor and its electric conductivity is very low. Low electric conductivity is one of the most important reasons for limited performance in the devices based on ZO film. Therefore, to improve electron transport rates in ZO film would be a significant way to achieve high performance devices. It has been suggested by many researchers to substitute ZO films with TM [67].

The electrical resistivity of ZO thin films undoped and doped by Fe with different concentrations was investigated by a two probe method. Figure 17a shows the I-V characteristics of the thin film samples having Ohmic behavior. The lowest resistivity of these films (4.6 <sup>×</sup> <sup>10</sup><sup>7</sup> <sup>Ω</sup>·cm) was measured for the films doped with 5% Fe corresponding to the highest DC power of the thin films. Otherwise, the resistivity of the co-sputtered films showed an initial increase for a DC power of 3 W and also for 5 W. This was followed by a

continuous reduction in resistivity as the DC power increased beyond 10 W, reaching the lowest value of (4.6 <sup>×</sup> <sup>10</sup><sup>7</sup> <sup>Ω</sup>·cm) in the films deposited under a DC power of 15 W, i.e., two orders of magnitude lower than the lowest resistivity obtained in the films prepared from the doped targets. These results manifest the pivotal role played by the method of incorporation of the dopant. *Ceramics* **2022**, *5,* FOR PEER REVIEW 25

**Figure 17.** (**a**) I-V characteristics of undoped and Fe doped ZO thin films and (**b**) variations of the electrical resistivity ρDC with Fe concentration of all investigated thin films. **Figure 17.** (**a**) I-V characteristics of undoped and Fe doped ZO thin films and (**b**) variations of the electrical resistivity ρDC with Fe concentration of all investigated thin films.

**4. Conclusions**  ZO and ZOFe films are deposited using RF magnetron sputtering with thicknesses 150 nm. The effect of Fe dopant on the structure, morphology, optical, magnetic and electric properties in Fe-doped ZO films was investigated. Curiously, we found that, both The electrical resistivity of undoped ZO thin films prepared by RF sputtering has been reported to be in the range of (104–10<sup>8</sup> Ω.cm), depending on the substrate temperature [68]. Further, a previous study on co-sputtered Fe—doped ZO thin films reported a decrease of resistivity by four orders of magnitude compared to the undoped film [16].

ZO and ZOFe films show a high crystalline quality with sharp band edge emission. As an interesting point, the energy band gap (*Eg*) was obtained by different numerical methods based on the experimental data in order to determine the most precise values. In addition, it has been pointed out that the refractive index (n) decreased with increasing Fe content which correlated to the optical band gap behavior. Likewise, an increasing of the saturation magnetization (Ms) values was noted with increasing Fe doped concentration while the coercivity (Hc) decreased. Finally, the resistivity of ZO is reduced as the concentration of Fe in ZO is increased slightly or, in other words, the conductivity is increased. The high saturation magnetization and conductivity values as well as the tunable optical properties of the investigated doped thin films make them promising multifunctional nanomaterials for a wide range of applications as As we pointed above, the resistances of the films increase with increase in low dopant concentration of Fe then decrease with higher concentrations, as shown in Figure 17b. The increasing trend of resistivity at low dopant concentration was reported by Cherif et al. [69]. The probable causes of increase in resistivity from undoped to Fe dopants in low concentration could be due to carrier traps formation at the film surface causing impediment to charge carriers and due to Fe incorporation, with the formation of interstitial metal atoms causing a decrease in the oxygen vacancies [70]. Meanwhile, increasing the concentration of Fe led to an increase the number of surface defects in the ZO matrix. Hence, the electrical resistivity of ZO is reduced as the concentration of Fe in ZO is increased slightly or, in other words, the conductivity increased. This makes our samples of thin films a promising material in very wide range of applications.

#### spintronics memory and magneto-optoelectronic devices. **4. Conclusions**

corresponding author.

**Author Contributions:** Conceptualization, A.F. and G.M.; Data curation, H.E. and C.S.; Formal analysis, A.F. and C.S.; Investigation, A.F. and C.S.; Methodology, A.F. and C.S.; Resources, G.M.; Supervision, G.M.; Validation, A.F., H.E. and G.M.; Visualization, A.F., H.E. and G.M.; Writing—original draft, A.F.; Writing—review and editing, A.F., H.E. and G.M. All authors have read and agreed to the published version of the manuscript. **Funding:** This research received no external funding. **Institutional Review Board Statement:** Not applicable. ZO and ZOFe films are deposited using RF magnetron sputtering with thicknesses 150 nm. The effect of Fe dopant on the structure, morphology, optical, magnetic and electric properties in Fe-doped ZO films was investigated. Curiously, we found that, both ZO and ZOFe films show a high crystalline quality with sharp band edge emission. As an interesting point, the energy band gap (*Eg*) was obtained by different numerical methods based on the experimental data in order to determine the most precise values. In addition, it has been pointed out that the refractive index (n) decreased with increasing Fe content which correlated to the optical band gap behavior.

**Data Availability Statement:** The data presented in this study are available on request from the

**Informed Consent Statement:** Not applicable.

Likewise, an increasing of the saturation magnetization (Ms) values was noted with increasing Fe doped concentration while the coercivity (Hc) decreased. Finally, the resistivity of ZO is reduced as the concentration of Fe in ZO is increased slightly or, in other words, the conductivity is increased. The high saturation magnetization and conductivity values as well as the tunable optical properties of the investigated doped thin films make them promising multifunctional nanomaterials for a wide range of applications as spintronics memory and magneto-optoelectronic devices.

**Author Contributions:** Conceptualization, A.F. and G.M.; Data curation, H.E. and C.S.; Formal analysis, A.F. and C.S.; Investigation, A.F. and C.S.; Methodology, A.F. and C.S.; Resources, G.M.; Supervision, G.M.; Validation, A.F., H.E. and G.M.; Visualization, A.F., H.E. and G.M.; Writing—original draft, A.F.; Writing—review and editing, A.F., H.E. and G.M. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **References**

