Abstract

Recently, metal oxide nanolayered semiconductors revealed their increasing usefulness as UV detectors and transparent conductors in optoelectronic devices. In the present paper, a simple fabrication process has been carried out to prepare layered TCO compounds with different controlled thicknesses. Conjoint physical investigations allowed discussion of the validity of optimality in terms of geometrical parameters. A synthetic function based on pondered physical parameters was a practical guide to reach optimality.

1. Introduction

Recently, several nanolayered transparent conducting oxides TCO have been investigated [120] for their interesting optical, mechanical, and electrical performance and besides all their high optical transparence in the visible domain. Commonly, these oxides could be synthesized using several methods [21, 22] like reactive evaporation, electron beam evaporation (EBE), pulsed laser deposition (PLD), chemical vapor deposition (CVD), sol-gel coating, and chemical spray pyrolysis.

In this study, TCO have been prepared by spray pyrolysis technique using, as main precursor, a solution of zinc acetate dissolved in deionized water. Several characterization techniques such as XRD, optical spectra, and hardness investigations have been applied to differently thick obtained samples [2326]. The thickness-dependent performance of the different as-grown layers has been additionally investigated in terms of mechanical, optothermal and microhardness behaviours. The existence of an optimal thickness was investigated using the BPES protocol and compared to some recently published results [27].

2. Experiment

2.1. ZnO Layers Elaboration

ZnO thin layers have been prepared by the technique of chemical reactive technique in liquid phase spray. The obtained layers’ structural and morphological properties, as well as synthesising details, have been recently published [2326].

ZnO thin films having different thicknesses were obtained on 20 × 10 × 3 mm3 Pyrex glass substrates by spraying an aqueous solution (mixture of water and Propanol 2) containing zinc acetate (Zn (CH3CO2)2) and a small amount of acetic acid, in order to avoid the precipitation of zinc hydroxide Zn(OH)2. The solution and gas (Nitrogen) flow rates were kept constant at 2 cm3·min−1 and 4.0 l min−1, respectively, corresponding to a minispray pyrolysis. The substrate temperature Ts, of the order of 460°C (optimal value under the given experimental conditions, as already confirmed by A. Amlouk et al. [23, 24]), was used to prepare these films according to in-room precedent successful attempts of deposition by the spray pyrolysis process [2834]. Six samples have been elaborated: with thicknesses 𝑑=0.16, 0.4, 0.58, 0.64, 0.84, and 1.03 μm, respectively.

2.2. Characterization Techniques

Precedent X-ray diffraction analyses of the prepared layers were performed by a copper-source diffractometer (Analytical X Pert PROMPD), with the wavelength (λ = 1.54056 Å). The optical transmission and reflection measures have been achieved using a common spectrophotometer (Shimadzu UV 3100S) equipped with an integrating sphere (LISR 3200). The spectrophotometer consists of double-beam monochromator (UV-2OOS, Shimazdu apparatus) with monochromatization performed by two flat silicon crystals in Laue diffraction on the (111) plane and enough energy to make several types of accurate measures in a wide wave-length range (220–1800 nm). The load-independent hardness measurements have been performed by a standard Microhardness Vickers (Hv) apparatus.

3. Already Achieved Analyses and Discussion

XRD diagram analysis (Figure 1) shows that the thinner layers (𝑑<0.6μm) develop an exclusive preferred orientation of the crystallites with respect to the (002) reflection. Differently, the thicker ones, are identified, beyond the (002) peak, by an extra XRD peak: (101).

These features correspond to the hexagonal würtzite system (JCPDS card, file n° : 361451 𝑎=3.24982,𝑐=5.20661Å), which is generally associated, as stated by Dubey et al. [35], Kumar et al. [34] and Paraguay et al. [36], with the appearance of irregularly spaced pores. This structure results in an increasing level of surface roughness [37, 38]. It seems that a minimal thickness is needed for the stability of a regular hexagonal würtzite system. Beyond this minimal thicknesses, the epitaxial growth seems to occur with no alteration in the lattice parameters, as discussed earlier by Kaeble [38].

The already published reflection R(λ) spectra [27] presented an accurate guide for estimating the thickness of each sample. In fact, the oscillating behavior of the λ-dependent reflection in the visible domain traduces the fact that the reflectance reaches discrete p-indexed maximal values 𝑅(𝑝)max for some particular p-indexed values 𝜆𝑝max of the wavelength [3942].

These p occurrences are traduced by the following equations: 𝑅(𝑝)max𝑛𝜆(𝜆)=𝑝max2𝑛𝑆𝑛𝜆𝑝max2+𝑛𝑆,𝑛𝜆𝑝max=𝑛𝑆1+𝑅max(𝜆)1𝑅max,(𝜆)(1) where 𝑛𝑠 is the substrate optical index (𝑛𝑠1.55) and 𝑛 is the layer λ-dependent optical index deduced from the optical transmission-reflection spectra.

The knowledge of two successive couples of maximal values [40, 41] allows calculating the film’s thickness d: 𝜆𝑑=𝑝max𝜆𝑝+1max2𝜆𝑝max𝑛𝜆𝑝+1max𝜆𝑝+1max𝑛𝜆𝑝max.(2)

4. Additional Mechanical Moh’s Hardness and Optothermal Investigations

Hardness is the resistance offered by a material to localized plastic deformation caused by scratching or indentation. Load-independent Moh’s hardness (𝐻𝑀) estimated from the scratch test is on the Moh’s scale, which could be calculated out from the Microhardness Vickers (𝐻𝑣) using the formula [42]: 𝐻𝑀=0.675×3𝐻𝑉.(3) The set of hardness of measurements has been performed using a common diamond-pyramidal-indenter under a prefixed load. The obtained imprints dimensions have been exploited for yielding the Microhardness Vickers (𝐻𝑣) of each sample as a synthetic value at room temperature using a Vickers diamond pyramidal indenter equipped with squared pyramid with 136° as summit angle (Figure 2). During experiment, it has been verified experimentally that the indenter penetration depth was always less than the layer’s thickness. The experimental values are gathered in (Figure 3).

The differences noticed in (Figure 3) can be explained, for the thinner layers, by the hardness contribution of the Pyrex substrate. The remaining layers, develop a mean hardness varying in a narrow interval. If the contribution of the substrate can be omitted, it seems that there is a maximum hardness zone situated in the interval of thickness [0.60–0.80 μm]. This noticed feature is probably due to the establishment of the hexagonal würtzite system which is harder to demolish than unorganized structures [33, 34, 36, 37].

The effective absorptivity 𝛼, as defined in precedent studies [35, 43], is the mean normalized absorbance weighted by ̃𝐼(𝜆)AM1.5, the solar standard irradiance: 𝛼=10𝐼̃𝜆AM1.5̃𝜆𝑑̃𝜆×𝛼10𝐼̃𝜆)AM1.5𝑑̃𝜆,𝜆𝜆min,𝜆max̃[],𝜆𝜆0,1min=300.0nm;𝜆max=1800.0nm,(4) where ̃𝐼(𝜆)AM1.5 is the reference solar spectral irradiance, fitted using the Boubaker polynomials expansion scheme BPES [4447]: ̃𝐼(𝜆)=[(1/2𝑁0)𝑁0𝑛=1𝜃𝑛𝐵4𝑛(̃𝜆×𝛽𝑛)], where 𝛽𝑛 are the Boubaker polynomials [45, 46] 𝐵4𝑛 minimal positive roots, 𝜃𝑛 are given coefficients, 𝑁0 is a given integer, ̃𝛼(𝜆) is the normalized absorbance spectrum, and ̃𝜆 is the normalised wavelength.

The normalized absorbance spectrum ̃𝛼(𝜆) is deduced from the BPES by establishing a set of 𝑁 experimental measured values of the transmittance-reflectance vector (𝑇𝑖(̃𝜆𝑖);𝑅𝑖(̃𝜆𝑖))|𝑖=1𝑁 versus the normalized wavelength ̃𝜆𝑖|𝑖=1𝑁. Then the system (5) is set: 𝑅̃𝜆=12𝑁0𝑁0𝑛=1𝜉𝑛×𝐵4𝑛̃𝜆×𝛽𝑛,𝑇̃𝜆=12𝑁0𝑁0𝑛=1𝜉𝑛×𝐵4𝑛̃𝜆×𝛽𝑛,(5) where 𝛽𝑛 are the 4n-Boubaker polynomials 𝐵4𝑛 minimal positive roots, 𝑁0 is a given integer and 𝜉𝑛 and 𝜉𝑛 are coefficients determined through the Boubaker polynomials expansion scheme (BPES).

The normalized absorbance spectrum ̃𝛼(𝜆) is deuced from the relation: 𝛼̃𝜆=1𝑑424̃𝜆ln1𝑅𝑇̃𝜆4+̃𝜆2ln1𝑅𝑇̃𝜆4,(6) where 𝑑 is the layer thickness.

The obtained value of normalized absorbance spectrum ̃𝛼(𝜆) is a final guide to the determination of the effective absorptivity 𝛼 through (4).

The Amlouk-Boubaker optothermal expansivity ΨAB is a thermophysical parameter defined in precedent studies [35, 4349], as a 3D expansion velocity of the transmitted heat inside the material. It is expressed in m3 s−1, and calculated by the following: ΨAB=𝐷,𝛼(7) where 𝐷 is the thermal diffusivity and 𝛼 is the effective absorptivity.

This parameter gathers two totally different parameters: the thermal diffusivity which is a material intrinsic property, and the effective absorptivity which links the material response to the solar spectrum. The Amlouk-Boubaker optothermal expansivity ΨAB presents the advantage of describing the conjoint thermal-optical response of the material to the whole spectrum, not only to monochromatic excitation.

As optimality is a complex notion when several physical parameters are involved, a synthetic aggregate is proposed on the bases of both Amlouk-Boubaker optothermal expansivity ΨAB and load-independent Moh’s Hardness 𝐻𝑀. For this purpose, and according to the definition of the Boubaker polynomials expansion scheme (BPES) [4447], the conjoint values of the N (𝑁=6) obtained d-dependent optothermal expansivity ΨAB and Moh’s Hardness 𝐻𝑀 are ordered as ((ΨAB)𝑘,(𝐻𝑀)𝑘)|𝑘=1𝑁 within the ranges [(ΨAB)min,(ΨAB)max] and [(𝐻𝑀)min,(𝐻𝑀)max]. Consecutively, the two normalized variables Ψ𝐴𝐵, 𝑑 and 𝐻𝑀 are defined as follows: Ψ𝐴𝐵=ΨABΨABminΨABmaxΨABmin,𝐻𝑀=𝐻𝑀𝐻𝑀min𝐻𝑀max𝐻𝑀min,𝑑=𝑑𝑑min𝑑max𝑑min.(8) The thickness is hence expressed as follows: 𝑑ΨAB,𝐻𝑀=|||||12𝑁0𝑁0𝑞=1𝜉𝑞𝐵4𝑞Ψ𝐴𝐵𝛼𝑞12𝑁0𝑁0𝑞=1𝜉𝑞𝐵4𝑞𝐻𝑀𝛼𝑞|||||,(9) where 𝛼𝑞 is 4q-Boubaker polynomial minimal root [45], 𝑁0 is a prefixed integer, 𝜉𝑞|𝑞=1𝑁0 and 𝜉𝑞|𝑞=1𝑁0 are unknown coefficients. The final step consists of computing (using Matlab subroutine SSO45) the set of coefficients (𝜉𝑞,opt,𝜉𝑞,opt)|𝑞=1𝑁0 which minimizes the real functional Ψ(𝑥): Ψ𝑁,𝑁0=𝑁𝑘=1|||||12𝑁0𝑁0𝑞=1𝜉𝑞𝐵4𝑞Ψ𝐴𝐵𝛼𝑞12𝑁0𝑁0𝑞=1𝜉𝑞𝐵4𝑞𝐻𝑀𝛼𝑞|||||.(10) The obtained solution 𝑑(ΨAB,𝐻𝑀), for 𝑁=6, 𝑁0=20, is normalized then plotted in Figure 4.

5. Discussion and Perspectives

A study has been recently presented by Lee et al. [48] on the thickness-related optimization inside IZO/ZnO/PET multilayered structures. The authors have prospected the TCO-induced exclusive electrical performance (carriers’ concentration, mobility…) and conjectured the existence of an eventual optimal thickness. Earlier, Shan and Yu [49] attempted to verify optimality for a thickness around 500 nm. These intriguing observations have been confirmed through the features of Figure 4. In fact, the optimality is obtained for the lowest point (𝑑0.61μm). This optimality is deduced by estimating the thickness which corresponds to a maximum of the load-independent Moh’s Hardness conjointly with a minimal value of Amlouk-Boubaker optothermal expansivity ΨAB. For this purpose, an optimisation protocol has been established for calculating a pondered function 𝐹 of Moh’s Hardness and Amlouk-Boubaker optothermal expansivity ΨAB with positive weights for Moh’s Hardness and negative ones for ΨAB. For each value of the measured couple (Moh’s Sardness, ΨAB), the correspondent values of the thickness 𝑑 along with the test function 𝐹 were recorded. Optimality was reached for the value of which corresponded to maximum value of 𝐹.

Moreover, the mechanical study revealed that, in order to reach the intrinsic hardness value of the layer (required for TCO applications in shock circumstances or aggressive mediums…), while avoiding substrate effects, one must either use a load as low as possible, or ensure a thickness exceeding the micrometer. Still, these values can be refined by using the nanohardness measurements to obtain the accurate hardness values of such thin films.

6. Conclusion

In summary, we have discussed the effects of the controlled thickness of TCO thin films prepared using a simplified and already optimized spray pyrolysis setup. Besides commonly considered parameters of the obtained layers, new and original characterization protocols have been carried out. It has been stated that an optimal ZnO thin film thickness is required for TCO purposes via simultaneous morphological-structural-physical investigations.