Thermoelectric and structural properties of transparent Sb-Doped ZnO thin films sputtered in a confocal geometry

Citation: Faria, H.F.; Ribeiro, J.M.; Boll, T.; Tavares, C.J. Thermoelectric and Structural Properties of Transparent Sb-Doped ZnO Thin Films Sputtered in a Confocal Geometry. Coatings 2023,13, 735. https://doi.org/10.3390/ coatings13040735 Academic Editor: Emerson Coy Received: 8 March 2023 Revised: 29 March 2023 Accepted: 1 April 2023 Published: 4 April 2023 Copyright: © 2023 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/). coatings Article Thermoelectric and Structural Properties of Transparent Sb-Doped ZnO Thin Films Sputtered in a Confocal Geometry Helder Filipe Faria 1, Joana Margarida Ribeiro 1, Torben Boll 2and Carlos JoséTavares 1,* 1Physics Centre of Minho and Porto Universities (CF-UM-PT), University of Minho, 4835-386 Guimarães, Portugal 2Institute for Applied Materials (IAM-WK), Institute for Nanotechnology (INT), Karlsruhe Nano Micro Facility (KNMFi), Karlsruhe Institute of Technology (KIT), D-76344 Karlsruhe, Germany *Correspondence: ctavar[email protected]; Tel.: +351-253510474 Abstract: This study focuses on understanding the influence of low Sb doping on ZnO’s electrical, optical, and thermoelectrical properties, while also studying its structural and morphological parameters. For this, several ZnO films with varying Sb target current densities, in the range of 0–0.27 mA/cm 2 , were produced by DC magnetron sputtering in a confocal geometry. As a result, thin ZnO:Sb films with an average transparency in the visible region greater than 80% are obtained, revealing for optimized conditions an absolute Seebeck coefficient of 100 µ V/K and a respective power factor of 1.1 mW · m −1· K −2 at 300 K, effectively modifying the electrical, optical, and thermoelectrical properties of the material and ensuring its suitability for heat harvesting applications. From atom probe tomography experiments, a larger Zn content is registered at triple junctions of the grain boundary, which matches the approximately 25 nm crystallite grain size derived from the X-ray diffraction analysis. Keywords: ZnO; antimony; Seebeck coefficient; thermoelectric; thin films; magnetron sputtering; atom probe tomography 1. Introduction Energy is and will continue to be a crucial part of the everyday life of our society. According to the Annual Energy Outlook 2021 made by the EIA [ 1 ], by 2050, global energy use is expected to increase by nearly 50% compared to 2020. Similarly, the World Energy Outlook 2021 [ 2 ] found evidence supporting the energy demand growth and its effect on carbon emissions. As renewable energies still lack the potential to sustain the demand, other types of solutions have gained importance, namely heat. Heat is present in almost every system, and it has been shown to be responsible for a large portion of all wasted energy [3]. There are several types of heat-related systems, some of which utilize reutilized heat that is dissipated during industrial processes, and some that rely on direct electrical conversion, such as thermoelectric systems. Thermoelectric materials transform heat into electricity (or vice versa) and have drawn much attention as heat-harvesting sources. This technology is based on the Seebeck effect [ 4 ], and for a material to be considered a good thermoelectric, it has to have large Seebeck coefficient for maximum electricity conversion, large electrical conductivity, and small thermal conductivity. These are all factors that influence the performance of the material, and are equated by the figure of merit (ZT), which is a dimensionless quantity that represents the thermoelectric quality of the material, and is given by the following [4]: ZT = S2σ kT(1) Coatings 2023,13, 735. https://doi.org/10.3390/coatings13040735 https://www.mdpi.com/journal/coatings Coatings 2023,13, 735 2 of 21 where Tis the operating temperature in kelvin (K), Sis the Seebeck coefficient (V/K), σ is the electrical conductivity ( Ω−1· cm −1 , or S · cm −1 ), and kis the thermal conductivity (W·m−1·K−1). Alternatively, the thermoelectric power factor (PF, in W · m −1· K −2 ) [ 4 ] can also be used to quantify the material in cases where the thermal conductivity has not been measured or calculated, as is the present case. PF =S2σ(2) Several types of materials, including organic and inorganic, have been used in thin film technology as thermoelectric materials, as they have significant benefits when compared to their less efficient and more expensive bulk forms [ 5 – 7 ]. The low dimensional material design was proved to improve thermoelectric efficiency through more effective phonon scattering, contributing to a reduced lattice thermal conductivity, and, at the same time, increasing the PF [ 4 ], popularizing thin films and nanostructures. Studies on bismuth telluride, a popular thermoelectric material, achieved a PF of 2.5 mW · m −1· K −2 (at 473 K) and 1.2 mW · m −1· K −2 (at room temperature, RT) in lead-doped and flexible undoped Bi 2 Te 3 thin films, respectively [ 6 , 7 ]. However, although being very efficient as thermoelectrics, these materials are opaque and cannot be applied in optically transparent devices. Metal oxides have potential as transparent thermoelectrics. Several experimental studies on TiO 2 - and ZnO-based thermoelectric films have reported interesting thermoelectric performance, albeit with smaller power factors compared with the aforementioned bismuth telluride. For example, TiO 2 :Nb 150 nm thick films with a PF of 500 µ W · m −1· K −2 at 300 K (ZT = 0.18) [ 8 ]; ZnO:Ga,Bi and ZnO:Al,Bi 300 nm thick films with a PF in the range of 70–80 µ W · m −1· K −2 at 300 K [ 9 ]; ZnO:Al 150 nm thick films with a PF of ~200 µ W · m −1· K −2 at 510 K [ 10 ]; ZnO:Ti 100 nm thick films with a PF of ~1 µ W · m −1· K −2 at 473 K [ 11 ]; ZnO:S films with a PF of ~5 µ W · m −1· K −2 at 310 K [ 12 ]; ZnO:Sb 600 nm thick films with a PF of ~400 µ W · m −1· K −2 at 700 K [ 13 ]; p-type ZnO:Sb 500 nm thick films with a PF of ~50 µW·m−1·K−2at 773 K. In addition to electrical conductivity, some materials also possess high transparency in the visible range. These materials are denominated as transparent conducting oxides (TCOs) and have been investigated for a wide range of applications in optoelectronic devices including photovoltaics, displays, transparent electrodes in touch panels, optoelectronic interfaces, and window glass technologies [14,15]. The design of a TCO material with thermoelectric properties would enable the production of electrical energy through heat harvesting from the environment in all kinds of transparent devices. Even though a wide range of thermoelectric materials have been studied and explored, a large ZT value remains connected to bulk materials [ 16 ]. The problem is that, because of their small band gap, these types of materials tend to be optically opaque, not constituting a valid solution to the intended applications [ 17 ]. Therefore, a material that provides all the necessary characteristics of a TCO with good thermoelectric properties is still necessary to allow for practical applications. Zinc oxide (ZnO) is a well-known inorganic compound quite common in modern electronics, mostly due to its abundance and versatility. It has a wide direct band gap of 3.44 eV at low temperatures and 3.37 eV at room temperature [ 18 ]. It crystallizes in the hexagonal wurtzite structure, displaying an intrinsic n-type conductivity, although in recent years it has become a target of interest due to its semiconductor properties as researchers try to control the unintentional n-type conductivity and achieve p-type conductivity [ 19 ]. Furthermore, added to the possible conductivity characteristics, the optical properties of ZnO make it a viable solution for TCO technology, but its relatively large electrical resistivity is a constraint that is usually countered by doping. Usually, n-type doping of ZnO is achieved with Group-III elements, such as aluminum and gallium, that can act as substitutional elements for Zn and group-VII elements, such as chlorine and iodine, that can serve as substitutional elements for oxygen [ 20 – 23 ]. As for p-type doping, is not as straightforward as n-type doping because it is very difficult to obtain p-type doping in wide band gap semiconductors. Coatings 2023,13, 735 3 of 21 Other Group-V elements have been studied, and a promising yet also commonly researched element with good potential is antimony. Several different studies focused on and effectively achieved p-type conductivity, with small electrical resistivity, in Sbdoped ZnO thin films, through a variety of deposition methods, including molecular-beam epitaxy, pulsed laser deposition, RF magnetron sputtering, spray pyrolysis, and oxidative evaporation [ 24 – 28 ]. Still, the most successful study, which used an approach based on sol–gel and hydrothermal methods, achieved an electrical resistivity of 2.1 m Ω· cm [ 29 ]. These studies effectively show the electrical potential of ZnO:Sb thin films, but studies on their thermoelectric properties are still scarce. The novelty of this study is in the fabrication of thermoelectric ZnO:Sb thin films by confocal sputtering, with TCO characteristics and the complete evaluation of their optical, electrical, and thermoelectrical properties. 2. Materials and Methods Undoped and Sb-doped ZnO thin films were prepared by a custom-made DC magnetron sputtering system conceived at the Centre of Physics of the University of Minho (see Figure A1 for the setup’s schematic). A confocal geometry is associated with two circular magnetrons, each with a diameter of 10 cm. One had a target of ZnO with 99.99% purity fabricated by FHR Anlagenbau GmbH (Saxony, Germany), and the other had an Sb target also with 99.99% purity supplied by Photon Export. A series of depositions were executed but, before each one, the chamber was evacuated with a primary rotary pump and a turbomolecular pump to achieve a base pressure of around 10 −4 Pa. Each set of films, for a specific level of Sb doping, were deposited on 76 × 26 × 1 mm 3 glass substrates from Normax Lda (Leiria, Portugal), along with four 10 × 10 × 0.5 mm 3 Si substrates, cut from P/B doped Si-wafer <100>, from Siegert Wafers GmbH (Aachen, Germany), following the schematic in Figure A2. Before assembling the substrates inside the chamber, they were cleaned with isopropyl alcohol (2-propanol) and acetone in an ultrasonic bath, both for 15 min, removing contaminants on the surface of the substrates that may arise during storage and handling, in order maximize the adhesion at the film/substrate interface. The substrate holders were mounted onto a support controlled by a motor that operated in continuous rotation during deposition to ensure uniformity. The distance between the targets and the substrate was kept at 8 cm, and the substrate holder was heated before and during depositions. Research-grade Ar with 99.999% purity was used as a working gas to enable plasma formation. Ahead of every deposition, 3 min of ion etching at 500 V was performed in an Ar atmosphere at a pressure of 1.6 Pa. This process is used to remove oxides and impurities that possibly accumulated and were retained on the surface of the substrates, which could create atomic defects in the substrate, enabling better nucleation of the film during deposition. The series of depositions started with a dummy deposition of 2 h for cleaning (burn) of the ZnO target, which was new. This was followed by a control deposition where the Sb target had no current and was covered with an Al sheet to prevent unwanted interactions with the plasma. Then, a series of depositions, of 35 min each, were performed where the current applied in the Sb target was gradually increased with steps of 2 mA, starting from 0 mA in the first deposition, 1 mA for the second, and ending with 21 mA (in the range of 0–0.27 mA/cm 2 ). This was the only parameter that suffered variations, as all others were kept the same between depositions, intending to effectively study how it would change doping levels and how that change would affect the thin film properties. Already tested deposition process parameters used in other studies in the same setup were followed. Target current density, bias voltage, voltage applied in the ZnO target, and Ar flow were fixed at 10 mA · cm −2 , − 60 V, 400 V, and 40 sccm, respectively, maintaining a working pressure between 0.18 and 0.21 Pa. Optically transparent and conducting films were effectively fabricated, and their properties were accordingly studied. An in-depth study of the electrical properties of the ZnO:Sb thin films was carried out based on the Hall effect, utilizing an Ecopia AMP55T HMS-5000 Hall effect measuring Coatings 2023,13, 735 4 of 21 system, with a DC four-point probe apparatus in the Van der Pauw configuration. The Hall voltage (VH) was determined, and the respective carrier concentration (n) and mobility ( µ ) were calculated with the following Equations (3) and (4) [30]: n= |B|·I e·|VH|·t(3) µ= 1 e·ρ·n(4) where B is the magnetic field, which was kept at a constant 0.560 T, t is the film thickness, I is the injected current, and ρ is the electrical resistivity. Ten measurements were performed for each sample deposited on Si substrates, and an average was computed for each of the significant properties referred above, in an attempt to obtain a more accurate result. Values for optical transmittance and reflectance were acquired through measurements with a Shimadzu UV-2501PC UV–Vis spectrophotometer. The samples deposited on glass substrates were analyzed in the 300 to 900 nm wavelength range and the average for each quantity was estimated between 400 and 700 nm for posterior calculations and analysis. For both studies, a baseline was established using two samples of glass without any film deposited on them, and then one was kept as a reference when studying each sample. Finally, with the transmittance spectra and custom-made software called dCalc, the thickness and optical properties of each sample were calculated. The software utilizes interference patterns to determine the refractive index of the film and, consequently, its thickness, and follows the Swanepoel method to determine the optical properties, namely the absorption index and band gap [31]. A study on the morphology and cross-section of the thin films was performed with scanning electron microscopy (SEM), performed with an FEI NOVA NanoSEM 200. To perform cross-section imaging, the samples were cut with a diamond tip and fixed in the sample table with double-sided carbon tape. Atomic force microscopy (AFM) was also used to study the surface properties of the films, namely roughness, on a CSI Instruments Nano-Observer atomic force microscope. For each sample, two different zones of the film were analyzed, typically one near the center and one near the edge. A mean of both values was obtained afterwards for a more accurate evaluation. The surface chemistry of the produced ZnO:Sb films was studied with an XPS spectrometer (Kratos Axis-Supra instrument, at 3 Bs Group, University of Minho) equipped with a monochromatic Al-K α X-ray radiation source (1486.6 eV) operated at an X-ray power of 225 W. The photoelectron spectra were collected at the take-off angle of 90 ◦ with the sample surface by means of a hemispherical electron energy analyzer operated in the constant analyzer energy lens mode (CAE). The pass energies of 160 eV and 40 eV were used for the survey and high-resolution spectra, respectively. The binding energy was referenced by setting the binding energy of the C1s hydrocarbon peak (the most intensive component of the C1s spectrum) at 284.8 eV. An electron flood gun was used to compensate for surface charging effects. Furthermore, to improve the composition study, energy-dispersive X-ray spectroscopy (EDX) was performed in an EDAXPegasus X4M. X-ray diffraction (XRD) analyses were performed to study the crystallographic structure of the thin films using a Bruker AXS D8 Discover diffractometer operated in θ –2 θ geometry with CuK α radiation. XRD patterns were obtained with a step size of 0.02 ◦ and an integration time of 1.5 s. Afterwards, with the help of the software Fityk, the XRD diffraction patterns were analyzed. Atom probe tomography (APT) was employed to study the morphology of the grain boundaries of the thin films. The experiments were performed at the KNMFi Laboratory for Microscopy and Spectroscopy at the Karlsruhe Institute of Technology, Germany. The thin films were deposited onto a specialized silicon coupon with prefabricated microtips from CAMECA, and subsequently prepared for APT analysis using a Zeiss Auriga 60 Dual Coatings 2023,13, 735 5 of 21 Beam Focused Ion Beam (FIB) in order to obtain a small tip with an apex diameter of 60 nm, retaining a roughly 300 nm thick layer on top of the Si substrate. The tips were measured in a local electrode atom probe (LEAP 4000X HR) in UV-laser mode (30 pJ), with a pulse frequency of 125 kHz, at a temperature of 50 K and a detection rate of 0.4%. The data was reconstructed using the Ivas 3.6.14 software from CAMECA. Finally, the Seebeck coefficient was evaluated using custom-made equipment that consisted of three main parts: the chamber, the pump, and a custom-made controlling software and hardware system. The pump ensured a low vacuum of approximately 1 Pa inside the chamber during the experiments. Inside the chamber were placed the holding system and the temperature control setup. Two Peltier devices (Quick-Ohm Kupper & Co. GmbH, Wuppertal, Germany) operate as a heater and cooler for the edges of the sample, applying their effect in a 25 × 75 mm 2 area. These devices ensure the difference in temperature necessary for calculating the Seebeck coefficient, and the potential difference is measured using a two-probe contact geometry, proven efficient even for higher temperatures [ 32 ]. For each sample, the thermal gradients applied on the films deposited on glass substrates varied between 30 and 55 ◦ C with a step of 5 ◦ C, resulting in six measurements for increasing ∆ T. For each measurement, 10 min of thermal stabilization was allowed to obtain a good and coherent value. These values are controlled and registered by the third part of the system, the hardware and software, where a custom-made interface is responsible for controlling the applied temperature in each Peltier device, controlling its evolution throughout the analysis, and registering each respective potential difference. After the experiment, and with all the values registered, the obtained ∆ V values were plotted as a function of ∆ T and linearly fitted to determine the slope, and subsequently the Seebeck coefficient, similar to other applications of this same system [9]. 3. Results and Discussion In total, 13 samples were deposited for this study. The first one, an undoped ZnO thin film (ZnOX) served as a control sample, as it was deposited with the Sb target covered, so it would not affect the normal deposition of ZnO. Then, the rest of the samples were deposited with optimized conditions and an increasing target current was applied to the Sb target, in order to produce ZnO:Sb films with varying doping levels. These depositions allowed the study of the influence of Sb doping on the structural, optical, and electrical properties of the ZnO thin films. For the sake of simplicity, the title of each film corresponds to the current applied in the Sb target during deposition (ZnO:SbXX = the film deposited with XX current applied to the Sb target, in mA). 3.1. Optical Properties As previously mentioned, a study on the optical transmittance and reflectance of each sample was performed between 300 and 900 nm (see Figure A3 for the full spectra). Then, an average value between 400 and 700 nm for each measurement is calculated and plotted in Figure 1. The results show that the increase in Sb doping, caused by the larger applied current to the Sb target, generates a slight decrease in the average transmittance, but the smallest value remains larger than 80%, so all samples have good optical transparency. As for reflectance curves, there is no relevant change, as it seems that it is independent of the doping level. Hence, even with the presence of Sb, the films preserve the natural good transparency of ZnO in the Vis and UV, which represents one of the key factors of a good TCO material. Coatings 2023,13, 735 6 of 21 Coatings 2023, 13, x FOR PEER REVIEW 6 of 20 Figure 1. Average transmittance and reflectance as a function of the current applied to the Sb target. 3.2. Compositional and Structural Analysis The control (undoped) sample ZnOX and the doped samples ZnO:Sb00, ZnO:Sb01, ZnO:Sb11, and ZnO:Sb21 were subjected to SEM analysis. As result, several micrographs of the cross-sections and of the surfaces of each sample are obtained (see Figure 2) and, with the data, an average thickness of the studied films is calculated. Afterwards, with the transmittance spectra, the thickness of each film is calculated with the software dCalc. Both cases showed that a greater current density in the Sb target resulted in films with a larger thickness, as it is expected. Furthermore, both studies coincided in the values they presented, with a small difference in only a few samples that never reached 100 nm. The thinnest is the control sample that showed 578 nm of thickness, and the thickest is the ZnO:Sb21 sample with 790 nm of thickness. The dCalc software also produced the Tauc plot of each sample, estimating the direct energy band gap. The effect of the Sb doping is reflected in this analysis, as a decrease from 3.37 eV for the undoped ZnO to 3.25 eV for the Sb-doped ZnO is registered, without any significant change between doped samples. (a) (b) (c) (d) 1 3 5 7 9 111315171921 0 10 20 30 40 50 60 70 80 90 100 Sb Target Current (mA) Average Transmitance (%) 0 10 20 30 40 50 60 70 80 90 100 Average Reflectance (%) Figure 1. Average transmittance and reflectance as a function of the current applied to the Sb target. 3.2. Compositional and Structural Analysis The control (undoped) sample ZnOX and the doped samples ZnO:Sb00, ZnO:Sb01, ZnO:Sb11, and ZnO:Sb21 were subjected to SEM analysis. As result, several micrographs of the cross-sections and of the surfaces of each sample are obtained (see Figure 2) and, with the data, an average thickness of the studied films is calculated. Afterwards, with the transmittance spectra, the thickness of each film is calculated with the software dCalc. Both cases showed that a greater current density in the Sb target resulted in films with a larger thickness, as it is expected. Furthermore, both studies coincided in the values they presented, with a small difference in only a few samples that never reached 100 nm. The thinnest is the control sample that showed 578 nm of thickness, and the thickest is the ZnO:Sb21 sample with 790 nm of thickness. Coatings 2023, 13, x FOR PEER REVIEW 6 of 20 Figure 1. Average transmittance and reflectance as a function of the current applied to the Sb target. 3.2. Compositional and Structural Analysis The control (undoped) sample ZnOX and the doped samples ZnO:Sb00, ZnO:Sb01, ZnO:Sb11, and ZnO:Sb21 were subjected to SEM analysis. As result, several micrographs of the cross-sections and of the surfaces of each sample are obtained (see Figure 2) and, with the data, an average thickness of the studied films is calculated. Afterwards, with the transmittance spectra, the thickness of each film is calculated with the software dCalc. Both cases showed that a greater current density in the Sb target resulted in films with a larger thickness, as it is expected. Furthermore, both studies coincided in the values they presented, with a small difference in only a few samples that never reached 100 nm. The thinnest is the control sample that showed 578 nm of thickness, and the thickest is the ZnO:Sb21 sample with 790 nm of thickness. The dCalc software also produced the Tauc plot of each sample, estimating the direct energy band gap. The effect of the Sb doping is reflected in this analysis, as a decrease from 3.37 eV for the undoped ZnO to 3.25 eV for the Sb-doped ZnO is registered, without any significant change between doped samples. (a) (b) (c) (d) 1 3 5 7 9 111315171921 0 10 20 30 40 50 60 70 80 90 100 Sb Target Current (mA) Average Transmitance (%) 0 10 20 30 40 50 60 70 80 90 100 Average Reflectance (%) Figure 2. Cont. Coatings 2023,13, 735 7 of 21 Coatings 2023, 13, x FOR PEER REVIEW 7 of 20 (e) (f) (g) (h) (i) (j) Figure 2. SEM micrographs of the cross-sections and surfaces of the (a,b) ZnOX, (c,d) ZnOSb:00, (e,f) ZnO:Sb01, (g,h) ZnO:Sb11, and (i,j) ZnO:Sb21 samples. The scale bar is 500 nm and 1 μm for the micrographs on the left and right, respectively. Then, the same samples were subjected to EDX analysis. We expected to see an increase in the atomic concentration of Sb with the increase in the current applied to the Sb target, meaning that the doping level is larger, but this did not happen. The concentration of Sb in all doped films is approximately 0.5 at.%. This must be related to the fact the limit of solubility of Sb in the wurtzite structure is around ~0.5 at.%, and a further increase in the Sb target current density in confocal geometry does not overcome this threshold. The absence of Sb in the control sample verifies that covering the target prevented impurities. Next, an XPS analysis is performed, and the fitted spectra for the Zn 2p (doublet), O 1s (singlet), and Sb 3d (doublet) core line levels for an undoped sample (ZnOX) and a ZnO:Sb sample (ZnO:Sb21) are studied (Figures 3 and 4). The composition is approximately stoichiometric for the undoped sample, with a slight deficiency in Zn (48.7 at.%). It should be noted that XPS is only sensible to the topmost atomic layers (~7 nm below the surface); hence, this slight zinc deficiency may be attributed to some evaporation of this material, which competes with sputtering due to its low melting temperature. For the case of the Sb-doped sample, the sub-stoichiometry is reversed, with the film being deficient in oxygen (48.8 at.%), possibly due to the promotion of oxygen vacancies upon Sb doping (1.0 at.%) in the wurtzite structure; further photoluminescence experiments needed to be performed to confirm these vacancies. The latter value of Sb doping is larger than that is measured by EDX (0.34 at.%, in Table A1), since EDX averages the film bulk, while XPS only considers the topmost atomic layers; hence the registered increase in Sb at the surface Figure 2. SEM micrographs of the cross-sections and surfaces of the ( a , b ) ZnOX, ( c , d ) ZnOSb:00, ( e , f ) ZnO:Sb01, ( g , h ) ZnO:Sb11, and ( i , j ) ZnO:Sb21 samples. The scale bar is 500 nm and 1 µ m for the micrographs on the left and right, respectively. The dCalc software also produced the Tauc plot of each sample, estimating the direct energy band gap. The effect of the Sb doping is reflected in this analysis, as a decrease from 3.37 eV for the undoped ZnO to 3.25 eV for the Sb-doped ZnO is registered, without any significant change between doped samples. Then, the same samples were subjected to EDX analysis. We expected to see an increase in the atomic concentration of Sb with the increase in the current applied to the Sb target, meaning that the doping level is larger, but this did not happen. The concentration of Sb in all doped films is approximately 0.5 at.%. This must be related to the fact the limit of solubility of Sb in the wurtzite structure is around ~0.5 at.%, and a further increase in the Sb target current density in confocal geometry does not overcome this threshold. The absence of Sb in the control sample verifies that covering the target prevented impurities. Next, an XPS analysis is performed, and the fitted spectra for the Zn 2p (doublet), O 1s (singlet), and Sb 3d (doublet) core line levels for an undoped sample (ZnOX) and a ZnO:Sb sample (ZnO:Sb21) are studied (Figures 3and 4). The composition is approximately stoichiometric for the undoped sample, with a slight deficiency in Zn (48.7 at.%). It should be noted that XPS is only sensible to the topmost atomic layers (~7 nm below the surface); hence, this slight zinc deficiency may be attributed to some evaporation of this material, which competes with sputtering due to its low melting temperature. For the case of the Sb-doped sample, the sub-stoichiometry is reversed, with the film being deficient in oxygen (48.8 at.%), possibly due to the promotion of oxygen vacancies upon Sb Coatings 2023,13, 735 8 of 21 doping (1.0 at.%) in the wurtzite structure; further photoluminescence experiments needed to be performed to confirm these vacancies. The latter value of Sb doping is larger than that is measured by EDX (0.34 at.%, in Table A1), since EDX averages the film bulk, while XPS only considers the topmost atomic layers; hence the registered increase in Sb at the surface is possibly caused by diffusion processes. From Figures 3a and 4a, it can be seen that the shape and position of the Zn 2p 3/2 peak (1021.3 eV) is practically unaltered with doping, with the doublet separation being the same ( ∆ E = 23.1 eV), and the peaks full width at half maximum (FWHM, β ) are within the range of 1.5–1.7 eV. However, the impact of Sb doping is readily viewed for the O 1s core line (Figures 3b and 4b). Since the O 1s core line overlaps with the Sb 3d 5/2 main doublet peak, the fitting of these two core lines presents a challenge. The used strategy is to first fit the Sb 3d doublet (Figure 4c) and, from the position of the lower intensity doublet peak (Sb 3d 3/2 ), the relative areas between the 3d 5/2 :3d 3/2 (3:2) doublet peaks, the determined spin-orbit splitting (9.6 eV) of these peaks, the deduced position of Sb 3d 5/2 (530.2 eV), and its FWHM (1.2 eV), it is then possible to fit the O 1s core line (Figure 4b) by forcing the insertion of the 3d 5/2 in its envelope structure (component CSb in Table 1and Figure 4c). From the NIST X-ray Photoelectron Spectroscopy Database [ 33 ] it can be derived that the average binding energies for Sb metal, Sb 3+ in Sb 2 O 3 , and Sb 5+ in Sb 2 O 5 , are, respectively 528.04 eV, 530.13 eV, and 531.45 eV. Since from Table 1and Figure 4the binding energy ascribed to the Sb 3d5/2 is 530.2 eV, it can be concluded that the valence state of Sb in the ZnO wurtzite structure is 3+. For both samples, the O 1s core line is fitted with two contributions as follows. For the smaller binding energy (529.9–530.0 eV) component C1 is ascribed to Zn-O bonds and Sb-O bonds (doped sample), and its position only varies marginally (0.1 eV), with unaltered FWHM. Conversely, component C2 (531.6–532.0 eV) is attributed to adsorbed oxygen and OH groups, mostly from surface contamination, and C2 is not taken into account to determine film composition. The area of C1 is slightly larger (66%) for the undoped sample, in comparison to the Sb-doped sample (63%). Coatings 2023, 13, x FOR PEER REVIEW 9 of 20 1016 1018 1020 1022 1024 1026 2p 3/2 Intensity (arb. units) Binding energy (eV) exp. backgr. fit Zn 2p 527 528 529 530 531 532 533 534 Intensity (arb. units) Binding energy (eV) exp. C1 C2 backgr. fit O 1s (a) (b) Figure 3. (a) Zn 2p and (b) O1s XPS core lines and respective fits for the undoped ZnO film (NOx). 1016 1018 1020 1022 1024 1026 2p 3/2 Intensity (arb. units) Binding energy (eV) exp. backgr. fit Zn 2p 528 529 530 531 532 533 534 Intensity (arb. units) Binding energy (eV) exp. C1 CSb C3 backgr. fit O 1s (a) (b) 524 528 532 536 540 544 3d 5/2 3d 3/2 Intensity (arb. units) Binding energy (eV) exp. backgr. fit Sb 3d (c) Figure 4. (a) Zn 2p, (b) O1s, and (c) Sb 3d XPS core lines and respective fits for the Sb-doped ZnO film (ZnO:Sb21). Figure 3. (a) Zn 2p and (b) O1s XPS core lines and respective fits for the undoped ZnO film (NOx). Coatings 2023,13, 735 9 of 21 Coatings 2023, 13, x FOR PEER REVIEW 9 of 20 1016 1018 1020 1022 1024 1026 2p 3/2 Intensity (arb. units) Binding energy (eV) exp. backgr. fit Zn 2p 527 528 529 530 531 532 533 534 Intensity (arb. units) Binding energy (eV) exp. C1 C2 backgr. fit O 1s (a) (b) Figure 3. (a) Zn 2p and (b) O1s XPS core lines and respective fits for the undoped ZnO film (NOx). 1016 1018 1020 1022 1024 1026 2p 3/2 Intensity (arb. units) Binding energy (eV) exp. backgr. fit Zn 2p 528 529 530 531 532 533 534 Intensity (arb. units) Binding energy (eV) exp. C1 CSb C3 backgr. fit O 1s (a) (b) 524 528 532 536 540 544 3d 5/2 3d 3/2 Intensity (arb. units) Binding energy (eV) exp. backgr. fit Sb 3d (c) Figure 4. (a) Zn 2p, (b) O1s, and (c) Sb 3d XPS core lines and respective fits for the Sb-doped ZnO film (ZnO:Sb21). Figure 4. ( a ) Zn 2p, ( b ) O1s, and ( c ) Sb 3d XPS core lines and respective fits for the Sb-doped ZnO film (ZnO:Sb21). Table 1. Composition and core line binding energy positions derived from the XPS fits to the ZnO:Sb films data. ∆ E is the spin-orbit separation of the Zn 2p and Sb 3d doublets. β is the FWHM of the fitted peaks. A relative area percentage (%) is given for the two O 1s contributions (C1 and C2) with respective FWHM ( β 1 and β 2). CSb in the O1s spectra is attributed to the Sb 3d5/2 (530.2 eV) doublet peak. Sample Composition (at.%) Zn 2p O 1s Sb 3d5/2 Position C1 (β1) % Position /∆E/ CSb (β) % /∆E/ (β) (eV) C2 (β2) % (β) (eV) (eV) ZnOX Zn 49.7 529.9 - (1.1) 66% O 50.3 1021.3 23.1 (1.7) - Sb -531.6 (1.7) 34% ZnO:Sb21 Zn 50.2 530.0 (1.1) 63% O 48.8 1021.3 23.1 (1.5) 530.2 (1.2) 13% 530.2 9.6 (1.2) Sb 1.0 532.0 (1.5) 24% Coatings 2023,13, 735 16 of 21 in practical applications as optically transparent thermoelectric thin films. From these appealing results, it would be possible to implement these ZnO:Sb films in electronic devices and, perhaps one day, to use this material in display screens and windows of skyscrapers for thermal heat harvesting, positively contributing to solving the emerging global energy problem. Author Contributions: Conceptualization, C.J.T.; methodology, C.J.T.; software, C.J.T., J.M.R. and H.F.F.; validation, C.J.T. and T.B.; investigation, C.J.T., J.M.R. and H.F.F.; resources, C.J.T. and T.B.; data curation, J.M.R. and H.F.F.; writing—original draft preparation, H.F.F.; writing—review and editing, C.J.T. and J.M.R.; supervision, C.J.T.; project administration, C.J.T. and T.B.; funding acquisition, C.J.T. and T.B. All authors have read and agreed to the published version of the manuscript. Funding: This research is funded by FCT/PIDDAC through the Strategic Funds project reference UIDB/04650/2020-2023. Joana M. Ribeiro is grateful to the Fundação para a Ciência e Tecnologia (FCT, Portugal) for the Ph.D grant SFRH/BD/147221/2019. This work (proposal ID 2021-025-030112) was carried out with the support of the Karlsruhe Nano Micro Facility (KNMFi, www.knmf.kit.edu, accessed on 8 March 2023), a Helmholtz Research Infrastructure at Karlsruhe Institute of Technology (KIT, www.kit.edu, accessed on 8 March 2023). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: Data can be shared upon formal request. Conflicts of Interest: The authors declare no conflict of interest. Appendix A Coatings 2023, 13, x FOR PEER REVIEW 16 of 20 and promising solution for TCO-based applications. From atom probe tomography experiments, a more significant Zn content is registered on triple junctions at the grain boundary. Finally, considering all the properties previously mentioned, several ZnO:Sb films with good characteristics are developed, and the results testify the potential of this material in practical applications as optically transparent thermoelectric thin films. From these appealing results, it would be possible to implement these ZnO:Sb films in electronic devices and, perhaps one day, to use this material in display screens and windows of skyscrapers for thermal heat harvesting, positively contributing to solving the emerging global energy problem. Author Contributions: Conceptualization, C.J.T.; methodology, C.J.T.; software, C.J.T., J.M.R. and H.F.F.; validation, C.J.T. and T.B.; investigation, C.J.T., J.M.R. and H.F.F.; resources, C.J.T. and T.B.; data curation, J.M.R. and H.F.F.; writing—original draft preparation, H.F.F.; writing—review and editing, C.J.T. and J.M.R.; supervision, C.J.T.; project administration, C.J.T. and T.B.; funding acquisition, C.J.T. and T.B. All authors have read and agreed to the published version of the manuscript. Funding: This research is funded by FCT/PIDDAC through the Strategic Funds project reference UIDB/04650/2020-2023. Joana M. Ribeiro is grateful to the Fundação para a Ciência e Tecnologia (FCT, Portugal) for the Ph.D grant SFRH/BD/147221/2019. This work (proposal ID 2021-025-030112) was carried out with the support of the Karlsruhe Nano Micro Facility (KNMFi, www.knmf.kit.edu, accessed on 8 March 2023), a Helmholtz Research Infrastructure at Karlsruhe Institute of Technology (KIT, www.kit.edu, accessed on 8 March 2023). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: Data can be shared upon formal request. Conflicts of Interest: The authors declare no conflict of interest. Appendix A Figure A1. Schematic of the sputtering deposition chamber. Figure A1. Schematic of the sputtering deposition chamber. Coatings 2023,13, 735 17 of 21 Coatings 2023, 13, x FOR PEER REVIEW 17 of 20 Figure A2. Layout of the substrate holder. Appendix B (a) (b) Figure A3. Optical (a) transmittance and (b) reflectance of each sample. 300 400 500 600 700 800 900 0 10 20 30 40 50 60 70 80 90 100 Transmitance (%) Wavelength (nm) ZnO:Sb00 ZnO:Sb01 ZnO:Sb03 ZnO:Sb05 ZnO:Sb07 ZnO:Sb09 ZnO:Sb11 ZnO:Sb13 ZnO:Sb15 ZnO:Sb17 ZnO:Sb19 ZnO:Sb21 300 400 500 600 700 800 900 0 10 20 30 40 50 60 70 80 90 100 Reflectance (%) Wavelength (nm) ZnO:Sb00 ZnO:Sb01 ZnO:Sb03 ZnO:Sb05 ZnO:Sb07 ZnO:Sb09 ZnO:Sb11 ZnO:Sb13 ZnO:Sb15 ZnO:Sb17 ZnO:Sb19 ZnO:Sb21 Figure A2. Layout of the substrate holder. Appendix B Coatings 2023, 13, x FOR PEER REVIEW 17 of 20 Figure A2. Layout of the substrate holder. Appendix B (a) (b) Figure A3. Optical (a) transmittance and (b) reflectance of each sample. 300 400 500 600 700 800 900 0 10 20 30 40 50 60 70 80 90 100 Transmitance (%) Wavelength (nm) ZnO:Sb00 ZnO:Sb01 ZnO:Sb03 ZnO:Sb05 ZnO:Sb07 ZnO:Sb09 ZnO:Sb11 ZnO:Sb13 ZnO:Sb15 ZnO:Sb17 ZnO:Sb19 ZnO:Sb21 300 400 500 600 700 800 900 0 10 20 30 40 50 60 70 80 90 100 Reflectance (%) Wavelength (nm) ZnO:Sb00 ZnO:Sb01 ZnO:Sb03 ZnO:Sb05 ZnO:Sb07 ZnO:Sb09 ZnO:Sb11 ZnO:Sb13 ZnO:Sb15 ZnO:Sb17 ZnO:Sb19 ZnO:Sb21 Figure A3. Optical (a) transmittance and (b) reflectance of each sample. Coatings 2023,13, 735 18 of 21 Coatings 2023, 13, x FOR PEER REVIEW 18 of 20 Figure A4. Zoomed in XRD patterns and respective line fits for ZnOX, ZnO:Sb01, ZnO:Sb03, ZnO:Sb07, ZnO:Sb13, and ZnO:Sb21 films to highlight the small peak at around 36° that is ascribed to the Sb (002) atomic planes, which is not easily discerned in the full XRD patterns of Figure 5. Appendix C Table A1. Samples and respective important properties, absolute Seebeck coefficient |S|, electrical resistivity (ρ), carrier concentration (𝑛), carrier mobility (μ), lattice parameters (a, c), grain size (g) mean, and quadratic roughness (Ra, Rq). Thickness calculated and measured by SEM analysis (t(SEM)), Sb content measured by EDX analysis (Sb(EDX)), average transmittance and reflectance (T ,R ), and energy band gap (Eg). Sample |S| (μV/K) PF (mW∙m−1∙K−2) ρ (Ω∙cm) n (×1018cm−3) μ (cm2/V∙s) a (Å) c (Å) g (nm) Ra (nm) Rq (nm) t(SEM) (nm) Sb(EDX) (at.%) 𝐓  * (%) 𝐑  * (%) Eg (eV) ZnOX 2.8 0.2×10−4 4.46 −4.5 2.6 3.27 5.25 21 4.97 6.35 578.47 - 87.41 16.74 3.26 Zno:Sb00 6.5 0.2×10−3 2.60 −2.9 3.46 - - - - - 627.03 - 87.50 16.54 3.26 ZnoZnO:S b01 27.1 0.4×10−2 0.19 −12.6 1.22 3.28 5.25 25 4.82 6.35 702.30 0.51 87.01 15.21 3.25 ZnO:Sb03 95.9 1.1 0.08 −23.1 1.4 3.28 5.24 25 9.52 11.84 719.01 - 85.83 16.54 3.25 ZnO:Sb05 69.2 0.5 0.09 −53.3 1.4 - - - - - 740.74 - 85.00 16.45 3.25 ZnO:Sb07 74.7 0.4 0.13 −49.3 1.3 3.29 5.24 26 5.23 6.59 757.60 0.50 85.35 16.69 3.24 ZnO:Sb09 93.0 0.8 0.11 −20.0 1.8 - - - - - 771.38 - 85.03 16.42 3.26 ZnO:Sb11 74.4 0.5 0.11 −15.9 1.9 - - - - - 778.07 0.12 84.83 16.33 3.25 ZnO:Sb13 100.4 0.5 0.18 −104.1 0.8 3.27 5.24 25 8.93 11.21 782.57 - 84.59 16.69 3.25 ZnO:Sb15 96.2 0.6 0.16 −43.3 0.9 - - - - - 785.20 - 81.80 15.65 3.24 30 31 32 33 34 35 36 37 38 39 40 0 53 106 159 0 53 106 159 0 53 106 159 0 53 106 159 0 53 106 159 30 31 32 33 34 35 36 37 38 39 40 0 53 106 159 2θ (º) ZnOX Fit ZnO:Sb01 Fit Intensity (a.u.) ZnO:Sb03 Fit ZnO:Sb07 Fit (002) ZnO:Sb13 Fit ZnO:Sb21 Fit Figure A4. Zoomed in XRD patterns and respective line fits for ZnOX, ZnO:Sb01, ZnO:Sb03, ZnO:Sb07, ZnO:Sb13, and ZnO:Sb21 films to highlight the small peak at around 36 ◦ that is ascribed to the Sb (002) atomic planes, which is not easily discerned in the full XRD patterns of Figure 5. Coatings 2023,13, 735 19 of 21 Appendix C Table A1. Samples and respective important properties, absolute Seebeck coefficient |S| , electrical resistivity ( ρ ), carrier concentration ( n ), carrier mobility ( µ ), lattice parameters (a, c), grain size (g) mean, and quadratic roughness (Ra, Rq). Thickness calculated and measured by SEM analysis (t(SEM)), Sb content measured by EDX analysis (Sb(EDX)), average transmittance and reflectance (− T, − R), and energy band gap (Eg). Sample |S| (µV/K) PF (mW·m−1·K−2) ρ (Ω·cm) n (×1018cm−3) µ (cm2/V·s) a (Å) c (Å) g (nm) Ra (nm) Rq (nm) t(SEM) (nm) Sb(EDX) (at.%) − T * (%) − R * (%) Eg (eV) ZnOX 2.8 0.2 ×10−44.46 −4.5 2.6 3.27 5.25 21 4.97 6.35 578.47 - 87.41 16.74 3.26 Zno:Sb00 6.5 0.2 ×10−32.60 −2.9 3.46 - - - - - 627.03 - 87.50 16.54 3.26 ZnoZnO:Sb01 27.1 0.4 ×10−20.19 −12.6 1.22 3.28 5.25 25 4.82 6.35 702.30 0.51 87.01 15.21 3.25 ZnO:Sb03 95.9 1.1 0.08 −23.1 1.4 3.28 5.24 25 9.52 11.84 719.01 - 85.83 16.54 3.25 ZnO:Sb05 69.2 0.5 0.09 −53.3 1.4 - - - - - 740.74 - 85.00 16.45 3.25 ZnO:Sb07 74.7 0.4 0.13 −49.3 1.3 3.29 5.24 26 5.23 6.59 757.60 0.50 85.35 16.69 3.24 ZnO:Sb09 93.0 0.8 0.11 −20.0 1.8 - - - - - 771.38 - 85.03 16.42 3.26 ZnO:Sb11 74.4 0.5 0.11 −15.9 1.9 - - - - - 778.07 0.12 84.83 16.33 3.25 ZnO:Sb13 100.4 0.5 0.18 −104.1 0.8 3.27 5.24 25 8.93 11.21 782.57 - 84.59 16.69 3.25 ZnO:Sb15 96.2 0.6 0.16 −43.3 0.9 - - - - - 785.20 - 81.80 15.65 3.24 ZnO:Sb17 94.5 0.4 0.22 −13.3 1.3 - - - - - 787.07 - 83.69 15.99 3.25 ZnO:Sb19 86.1 0.4 0.20 −13.5 1.1 - - - - - 788.52 - 84.66 16.39 3.25 ZnO:Sb21 96.9 0.3 0.28 −32.7 0.6 3.29 5.24 24 9.84 12.10 789.70 0.34 81.91 15.93 3.24 * Average optical transmittance and reflectance calculated between 400 and 700 nm. Coatings 2023,13, 735 20 of 21 References 1. US Energy Information Administration, IEO 2021 Narrative, n.d. Available online: www.eia.gov (accessed on 19 January 2023). 2. International Energy Agency. World Energy Outlook 2021. 2021. Available online: www.iea.org/weo (accessed on 19 January 2023). 3. Jouhara, H.; Khordehgah, N.; Almahmoud, S.; Delpech, B.; Chauhan, A.; Tassou, S.A. Waste heat recovery technologies and applications. Therm. Sci. Eng. Prog. 2018,6, 268–289. [CrossRef] 4. Beretta, D.; Neophytou, N.; Hodges, J.M.; Kanatzidis, M.G.; Narducci, D.; Gonzalez, M.M.; Beekman, M.; Balke, B.; Cerretti, G.; Tremel, W.; et al. Thermoelectrics: From history, a window to the future. Mater. Sci. Eng. R Rep. 2018,138, 100501. [CrossRef] 5. Dresselhaus, M.S.; Chen, G.; Tang, M.Y.; Yang, R.G.; Lee, H.; Wang, D.Z.; Ren, Z.F.; Fleurial, J.; Gogna, P. New Directions for Low–Dimensional Thermoelectric Materials. Adv. Mater. 2007,19, 1043–1053. [CrossRef] 6. Zhou, Y.; Li, L.; Tan, Q.; Li, J.-F. Thermoelectric properties of Pb-doped bismuth telluride thin films deposited by magnetron sputtering. J. Alloy. Compd. 2014,590, 362–367. [CrossRef] 7. Nuthongkum, P.; Sakdanuphab, R.; Horprathum, M.; Sakulkalavek, A. [Bi]:[Te] Control, Structural and Thermoelectric Properties of Flexible Bi x Te y Thin Films Prepared by RF Magnetron Sputtering at Different Sputtering Pressures. J. Electron. Mater. 2017 , 46, 6444–6450. [CrossRef] 8. Ribeiro, J.; Correia, F.; Kuzmin, A.; Jonane, I.; Kong, M.; Goñi, A.; Reparaz, J.; Kalinko, A.; Welter, E.; Tavares, C. Influence of Nb-doping on the local structure and thermoelectric properties of transparent TiO 2 :Nb thin films. J. Alloy. Compd. 2020 , 838, 155561. [CrossRef] 9. Correia, F.; Salvador, P.; Ribeiro, J.; Mendes, A.; Tavares, C. Effect on the electrical and morphological properties of Bi incorporation into ZnO:Ga and ZnO:Al thin films deposited by confocal magnetron sputtering. Vacuum 2018,152, 252–260. [CrossRef] 10. Liu, S.; Li, G.; Lan, M.; Piao, Y.; Miyazaki, K.; Wang, Q. Effect of growth modes on electrical and thermal transport of thermoelectric ZnO:Al films. Acta Crystallogr. Sect. B Struct. Sci. Cryst. Eng. Mater. 2020,76, 259–266. [CrossRef] 11. Chao-Moo, W.; Vora-Ud, A.; Thaowankaew, S.; Muthitamongkol, P.; Seetawan, T. Effect of annealing treatment on thermoelectric properties of Ti-doped ZnO thin film. AIP Conf. Proc. 2018,2010, 020013. [CrossRef] 12. Gather, F.; Kronenberger, A.; Hartung, D.; Becker, M.; Polity, A.; Klar, P.J.; Meyer, B.K. Possibility of enhancing the thermoelectric figure of merit of ZnO by sulfur incorporation. Appl. Phys. Lett. 2013,103, 082115. [CrossRef] 13. Pham, A.T.T.; Le, O.K.T.; Phan, T.T.T.; Van Hoang, D.; Nguyen, T.H.; Le, N.D.; Phan, T.B.; Tran, V.C. Enhancing transparent thermoelectric properties of Sb-doped ZnO thin films via controlled deposition temperature. Vacuum 2022 ,202, 111137. [CrossRef] 14. Safeen, K.; Micheli, V.; Bartali, R.; Gottardi, G.; Safeen, A.; Ullah, H.; Laidani, N. Influence of intrinsic defects on the electrical and optical properties of TiO2:Nb films sputtered at room temperature. Thin Solid Films 2018,645, 173–179. [CrossRef] 15. Tao, J.; Pan, H.; Wong, L.M.; Wong, T.I.; Chai, J.W.; Pan, J.; Wang, S.J. Mechanism of insulator-to-metal transition in heavily Nb doped anatase TiO2.Mater. Res. Express 2014,1, 015911. [CrossRef] 16. Ding, L.; Akbarzadeh, A.; Tan, L. A review of power generation with thermoelectric system and its alternative with solar ponds. Renew. Sustain. Energy Rev. 2018,81, 799–812. [CrossRef] 17. Faustino, B.M.M.; Gomes, D.; Faria, J.; Juntunen, T.; Gaspar, G.; Bianchi, C.; Almeida, A.; Marques, A.; Tittonen, I.; Ferreira, I. CuI p-type thin films for highly transparent thermoelectric p-n modules. Sci. Rep. 2018,8, 6867. [CrossRef] [PubMed] 18. Mang, A.; Reimann, K.; Rübenacke, S. Band gaps, crystal-field splitting, spin-orbit coupling, and exciton binding energies in ZnO under hydrostatic pressure. Solid State Commun. 1995,94, 251–254. [CrossRef] 19. Janotti, A.; Van de Walle, C. Fundamentals of zinc oxide as a semiconductor. Rep. Prog. Phys. 2009,72, 126501. [CrossRef] 20. Ozgür, Ü.; Alivov, Y.I.; Liu, C.; Teke, A.; Reshchikov, M.A.; Do˘gan, S.; Avrutin, V.; Cho, S.-J.; Morkoç, H. A comprehensive review of ZnO materials and devices. J. Appl. Phys. 2005,98, 041301. [CrossRef] 21. Castro, M.; Cerqueira, M.; Rebouta, L.; Alpuim, P.; Garcia, C.; Júnior, G.; Tavares, C. Influence of hydrogen plasma thermal treatment on the properties of ZnO:Al thin films prepared by dc magnetron sputtering. Vacuum 2014,107, 145–154. [CrossRef] 22. Wang, X.; Huang, X.; Wong, Z.M.; Suwardi, A.; Zheng, Y.; Wei, F.; Wang, S.; Tan, T.L.; Wu, G.; Zhu, Q.; et al. Gallium-Doped Zinc Oxide Nanostructures for Tunable Transparent Thermoelectric Films. ACS Appl. Nano Mater. 2022,5, 8631–8639. [CrossRef] 23. Singh, G.; Shrivastava, S.B.; Jain, D.; Pandya, S.; Shripathi, T.; Ganesan, V. Effect of indium doping on zinc oxide films prepared by chemical spray pyrolysis technique. Bull. Mater. Sci. 2010,33, 581–587. [CrossRef] 24. Xiu, F.X.; Yang, Z.; Mandalapu, L.J.; Zhao, D.T.; Liu, J.L.; Beyermann, W.P. High-mobility Sb-doped p-type ZnO by molecular-beam epitaxy. Appl. Phys. Lett. 2005,87, 152101. [CrossRef] 25. Pan, X.; Ye, Z.; Li, J.; Gu, X.; Zeng, Y.; He, H.; Zhu, L.; Che, Y. Fabrication of Sb-doped p-type ZnO thin films by pulsed laser deposition. Appl. Surf. Sci. 2007,253, 5067–5069. [CrossRef] 26. De Campos, J.A.; Viseu, T.M.; Rolo, A.; Barradas, N.P.; Alves, E.; de Lacerda-Arôso, T.; Cerqueira, M.F. Electrical and Raman Scattering Studies of ZnO:P and ZnO:Sb Thin Films. J. Nanosci. Nanotechnol. 2010,10, 2620–2623. [CrossRef] 27. Kumar, N.S.; Bangera, K.V.; Shivakumar, G.K. Properties of antimony doped ZnO thin films deposited by spray pyrolysis technique. Semiconductors 2015,49, 899–904. [CrossRef] 28. Li, G.; Lin, X.; Liu, S.; Jia, B.; Wang, Q. Sb content dependent thermoelectric properties of the p-type ZnO:Sb films fabricated by oxidation method. Appl. Surf. Sci. 2018,439, 82–87. [CrossRef] Coatings 2023,13, 735 21 of 21 29. Fathima, N.; Pradeep, N.; Balakrishnan, J. Investigations of the effects of electrode geometry and mechanical stress on Antimony doped Zinc Oxide nanostructures based MSM UV photodetectors fabricated on flexible substrates. Sol. Energy Mater. Sol. Cells 2019,194, 207–214. [CrossRef] 30. van der Pauw, L.J. A method of measuring specific resistivity and hall effect of discs of arbitrary shape. In Semiconductor Devices Pioneering Papers; World Scientific Publishing Co Pte Ltd.: Singapore, 1991; pp. 174–182. [CrossRef] 31. Sánchez-González, J.; Díaz-Parralejo, A.; Ortiz, A.; Guiberteau, F. Determination of optical properties in nanostructured thin films using the Swanepoel method. Appl. Surf. Sci. 2006,252, 6013–6017. [CrossRef] 32. Martin, J.; Wong-Ng, W.; Green, M.L. Seebeck Coefficient Metrology: Do Contemporary Protocols Measure Up? J. Electron. Mater. 2015,44, 1998–2006. [CrossRef] 33. NIST X-ray Photoelectron Spectroscopy (XPS) Database, Version 3.5, (n.d.). Available online: http://srdata.nist.gov/xps/ (accessed on 22 March 2003). 34. Correia, F.C.; Ribeiro, J.M.; Ferreira, A.; Reparaz, J.S.; Goñi, A.R.; Boll, T.; Mendes, A.; Tavares, C.J. The effect of Bi doping on the thermal conductivity of ZnO and ZnO:Al thin films. Vacuum 2023,207, 111572. [CrossRef] 35. Ribeiro, J.M.; Rodrigues, F.J.; Correia, F.C.; Pudza, I.; Kuzmin, A.; Kalinko, A.; Welter, E.; Barradas, N.P.; Alves, E.; LaGrow, A.P.; et al. The influence of Sb doping on the local structure and disorder in thermoelectric ZnO:Sb thin films. J. Alloy. Compd. 2023 , 939, 168751. [CrossRef] Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.