Next Article in Journal
Bi/Mn-Doped BiOCl Nanosheets Self-Assembled Microspheres toward Optimized Photocatalytic Performance
Next Article in Special Issue
Si-Based Polarizer and 1-Bit Phase-Controlled Non-Polarizing Beam Splitter-Based Integrated Metasurface for Extended Shortwave Infrared
Previous Article in Journal
Recent Advances on the Design and Applications of Antimicrobial Nanomaterials
Previous Article in Special Issue
Numerical Approach to the Plasmonic Enhancement of Cs2AgBiBr6 Perovskite-Based Solar Cell by Embedding Metallic Nanosphere
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Increasing the Accuracy of the Characterization of a Thin Semiconductor or Dielectric Film on a Substrate from Only One Quasi-Normal Incidence UV/Vis/NIR Reflectance Spectrum of the Sample

1
Scientific Research Section (NIS), Technical University, 1000 Sofia, Bulgaria
2
Department of Microelectronics, Faculty of Electronics Engineering and Technologies, Technical University, 1000 Sofia, Bulgaria
3
Faculty of Science, Department of Condensed-Matter Physics, University of Cadiz, Puerto Real, 11510 Cadiz, Spain
4
Photovoltaic Solar Energy Unit, Centre for Energy, Environmental and Technological Research (CIEMAT), Avenida Complutense 40, 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Nanomaterials 2023, 13(17), 2407; https://doi.org/10.3390/nano13172407
Submission received: 27 July 2023 / Revised: 13 August 2023 / Accepted: 18 August 2023 / Published: 24 August 2023

Abstract

:
OEMT is an existing optimizing envelope method for thin-film characterization that uses only one transmittance spectrum, T(λ), of the film deposited on the substrate. OEMT computes the optimized values of the average thickness, d ¯ , and the thickness non-uniformity, Δd, employing variables for the external smoothing of T(λ), the slit width correction, and the optimized wavelength intervals for the computation of d ¯ and Δd, and taking into account both the finite size and absorption of the substrate. Our group had achieved record low relative errors, <0.1%, in d ¯ of thin semiconductor films via OEMT, whereas the high accuracy of d ¯ and Δd allow for the accurate computation of the complex refractive index, N ˙ (λ), of the film. In this paper is a proposed envelope method, named OEMR, for the characterization of thin dielectric or semiconductor films using only one quasi-normal incidence UV/Vis/NIR reflectance spectrum, R(λ), of the film on the substrate. The features of OEMR are similar to the described above features of OEMT. OEMR and several popular dispersion models are employed for the characterization of two a-Si films, only from R(λ), with computed d ¯ = 674.3 nm and Δd = 11.5 nm for the thinner film. It is demonstrated that the most accurate characterizations of these films over the measured spectrum are based on OEMR.

1. Introduction

Thin semiconductor and dielectric films are used in photonic circuits, solar cells, thin-film transistors, holography, thin-film batteries, etc. [1,2]. Using such a film in a device requires knowledge of its characteristics, which depends on its composition and preparation technology [3,4]. The optical characteristics of a film can be derived via optical characterization, where a narrow light beam irradiates the surface of the film, thus forming a light spot there [5,6]. The main spectral characteristic of a film, which can be determined using such techniques, is its complex refractive index N ˙ ( λ ) = n ( λ ) + i k ( λ ) , where n(λ) is the refractive index, k(λ) is the extinction coefficient, λ is the wavelength of light, and “i” is the imaginary unit [7,8]. The main thickness parameters of a film, measurable using such techniques, are the average thickness, d ¯ , over the light spot and the thickness non-uniformity, ∆d = [max(d) − min(d)]/2 ≥ 0, over the same spot, where d is the thickness in a particular point within the spot. Additionally, ∆d can be described as a larger-scale version of the surface roughness of the film [9].
Film characterizations are facilitated by employing dispersion models (DMs), usually based on electron oscillators, which can be classified as narrow-spectrum DMs and broad-spectrum DMs. Narrow-spectrum DMs usually formulate n(λ) or k(λ) only in the spectral region of weak and medium absorption, where n2(λ) >> k2(λ) and E < Eb > Eg, where E (eV) = 1239.8/λ (nm) is the photon energy, Eg is the optical bandgap, and Eb is the photon energy above which interband transitions prevail. The Wemple–DiDomenico model, which is of this type, representing an undamped single oscillator, is valid for covalent, ionic, glassy, and amorphous semiconductor materials and provides the following [10,11]:
1 n 2 1 E 0 E d E 2 E 0 E d for   E < E b > E g ,
where E0 and Ed are the energy and the strength of the oscillator. Therefore, E0 and Ed can be calculated via the linear regression from a Wemple–DiDomenico plot depicting (n2 − 1)−1 as a function of E2 [11]. It is also seen from Equation (1) that the static refractive index can be expressed as n(0) = n(E = 0) = (1 + Ed/E0)1/2.
In addition, the absorption coefficient α = 4πk/λ of amorphous materials in the region of weak and medium absorption can be approximated using the Urbach rule, corresponding to the structural disorder-generating localized electronic states and the Urbach tail as follows [12,13]:
α α 0 exp ( E E U ) for   E < E b > E g ,
where α0 = α(E→0) > 0, and the Urbach energy EU quantifies the energetic disorder in the band edges. However, Equation (2) is inaccurate for very long wavelengths because it provides α(E→0) ≈ α0 > 0, although α(E = 0) = 4πk(0)/λ ~ 4πk(0)E(0) = 0.
Broad-spectrum DMs formulate the complex dielectric function ε ˙ ( E ) = ε r ( E ) + i ε i ( E ) = N ˙ 2 = n 2 k 2 + i 2 n k or the respective complex refractive index of a material over wider spectra, including the UV/Vis/NIR spectral region with wavelengths λ ≈ [150, 3000] nm [14,15]. Arguably the most popular broad-spectrum DMs are the Tauc–Lorentz DM of Jellison and Modine (TL) [16,17,18], the Campi–Coriasso DM (CC) [19,20], and the new amorphous dispersion formula DM (NA) [21,22]. TL uses the Tauc joint density of states and a Lorentz oscillator representing a damped harmonic oscillator. CC is a version of TL, in which the normalized oscillator strength is shifted by Eg, which is achieved by substituting E with EEg. NA has been established in order to give a Lorenzian shape to the expressions for n(E) and k(E). The main problem concerning the utilization of TL, CC, and NA for the characterization of amorphous materials is that in all of them, k(λ) = 0 is assumed for E < Eg, thus leading to characterization errors in the range λ > λg (nm) = 1239.8/Eg (eV).
Unlike the above three DMs, the Gaussian DM (GA) provides symmetrical εi(E), which makes it applicable for descriptions of absorption with peaks at E < Eg associated with the disorder, dangling bonds, impurities, and vacancies [23,24]. Furthermore, the contributions to the dielectric function caused via absorption with peaks above the measured energy spectral range can be included as a pole (PE) in this DM, i.e., a Lorentz oscillator with vanishing broadening [25,26].
With regard to the above two paragraphs, the number of employed oscillators is printed before the two-letter abbreviation of any DM allowing the usage of multiple oscillators henceforth in this paper. The exemplary calculated real and imaginary parts of the complex dielectric function and its respective normal incidence reflectance, R = [(n − 1)2 + k2]/[(n + 1)2 + k2], at the boundary between air and material, are depicted in Figure 1 for 1TL, 1CC, 1NA, and 1GA.
It is inferred from Figure 1a that max[εr(E)] and max[εi(E)] are positioned at the highest energies for 1NA as well as that the refractive index n{E = [0,max(εr)]} is the largest also for 1NA. Accordingly, the data from Figure 1c indicate that the reflectance of UV/Vis/NIR light is significantly larger, assuming the validity of NA compared to TL, CC, and GA. It is also apparent from Figure 1 that the results regarding TL are the closest to these for CC.
To increase the accuracy of the characterization of amorphous materials, the broad-spectrum single oscillator Tauc–Lorentz–Urbach DM (TLU) has been developed. There are two versions of TLU, and in both of them is the assumed validity of TL for E > Eb and the Urbach tail for EEb. The TLU developed by Foldyna (TLUF) has used the expression εi(E) = const/E × exp(E/EU) in the range EEb, assuming that the Urbach tail is formulated via Equation (2), and the variation in n is negligible in this range [27,28]. However, TLUF is inaccurate for very long wavelengths because the above expression provides εi(E→0) = ∞. This particular problem has been resolved in the TLU of Rodriguez de Marcos (TLUR) [29,30], which utilizes the expression εi(E) = const × E × exp(E/EU) in the range EEb, as it provides εi(E→0) = 0.
Notably, both TLUF and TLUR cannot employ more than one oscillator, which should decrease their accuracy for materials with intricate electronic band structures. The more accurate characterization of such materials could be achieved using DMs with multiple oscillators. In this respect, multiple-oscillator TL, CC, and GA can be included in the universal DM of Franta et al. (UD), which, in principle, allows for film characterization over the entire electromagnetic spectrum [31,32]. In UD, εr(E) and εi(E) are obtained by adding, respectively, parametrized electric susceptibility contributions from the individual oscillators, whereas each contribution to εr(E) is formulated using Kramers–Kronig integration utilizing the corresponding contribution to εi(E). Moreover, UD can readily be applied to amorphous materials because it can include parametrized electric susceptibility contributions to εr(E) and εi(E) from the Urbach tail.
Spectroscopic ellipsometry (SE) and spectrophotometry are primary and complementary techniques for the optical characterization of thin films. SE uses oblique incidence of linearly polarized light, a change in its polarization state due to its reflection from the film, and a selected DM [5,33]. For example, ∆d can be computed via an SE analysis of the depolarization of the reflected light [34,35]. When the film is on a light-transmitting substrate, reflections from the back surface of the substrate are problematic for SE because they create an incoherent additive (without phase information) to the coherent interaction (with phase information) between the light reflected from the two surfaces of the film. This problem can be circumvented by appropriate backside roughening [36] or attaching light scattering refractive index matching material to the back surface of the substrate [37]. Furthermore, the main SE parameters psi-delta have been computed in recent papers [38,39] for a mostly transparent film on a transparent substrate, depending on the spectral bandwidth ∆λ, the thickness non-uniformity ∆d of the film, its surface roughness, and its absorption. Nevertheless, it seems that there are no publications reporting the computation of all film characteristics printed in red in Figure 2, in the case of a weakly absorbing substrate, from experimental SE psi-delta spectra and taking into account ∆λ (or its respective spectral slit width). The most likely reason for this is the complexity of the phase change in the polarization state of light, during the reflection, depending on all of these film characteristics.
In the spectrophotometry of thin films are measured transmittance spectra T(λ) and/or reflectance spectra R(λ), which are most often of a sample consisting of a film deposited on a substrate [40,41]. The spectra T(λ) and R(λ) are scanned using a spectrophotometer, which can commonly operate in the UV/Vis/NIR spectral region [14,15]. T(λ) of such a sample is measured usually at the normal incidence of the light beam to the film and R(λ) at quasi-normal incidence to the film, both of which result in diminishing sideway leaking of light, disregarding the polarization state of light, and reducing mathematical complications [42,43]. Also, it is preferable to perform film characterization only from T(λ) or R(λ) because repositioning the sample for using both of them can lead to not identical light spots on the surface of the film and, therefore, to different averages of film thicknesses over these light spots [40,44].
In this present paper, the optical characteristics of the substrate are designated by the subscript “s”. Moreover, we studied only the quasi-normal incidence reflectance R(λ) of a thin film on a substrate and the normal incidence transmittance T(λ) of such a sample. Furthermore, only the common case of a dominantly coherent interaction between the light reflected from the two surfaces of this film is considered, which is formulated as d ¯ << λ2/(2nλ), as well as incoherent interaction between the light reflected from the two surfaces of the substrate, i.e., ds >> λ2/(2nsλ) [45]. The reflectance R(λ) of a thin film with non-uniform thickness, named non-uniform film, on such a substrate, known as a thick planar substrate, is illustrated in Figure 2.
There are several spectrophotometric methods for thin film characterizations employing R(λ), as shown in Figure 2. For example, n(λ) and k(λ) have been calculated from the measurements of d, R(λ), and T(λ) or the absorbance A(λ) of the sample, assuming that the film has uniform thickness d (such film is named uniform film) and the substrate is infinite [46,47,48,49]. Furthermore, it is well established that λg < 2000 nm (Eg > 0.6 eV) is valid for the vast majority of semiconductors and dielectrics [50,51], which means that they have a significant spectral range of wavelengths λ > λg < 2000 nm with weak absorption and n2(λ) >> k2(λ) << 1 in the UV/Vis/NIR region. The facts from the last two paragraphs indicate that there are apparent minima and maxima in this spectral range of T(λ) and R(λ) of the sample from Figure 2 due to thin-film interference, provided that the semiconductor or dielectric film is not too thin, i.e., d ¯ > λ/(2n) [45]. In such a case, it is possible to draw two envelopes around the apparent maxima and around the apparent minima of T(λ) or R(λ), as well as to apply the envelope method (EM) for film characterization in this spectral range [52,53]. The tangency points between the spectrum T(λ) or R(λ) and its two envelopes can be represented by the respective tangency wavelengths λt(l), where the integer l is the tangency wavelength number counted from the longer wavelengths end of the spectrum. Spectra with apparent minima and maxima due to thin-film interference are referred to as interference spectra.
EM is a spectrophotometric method that utilizes an interference fringes equation, including all λt(l), together with the interference spectrum T(λ) or R(λ) and its two envelopes [52,53]. This additional information makes it possible to achieve a higher accuracy of computation of the d ¯ , ∆d, and n(λ) of the film, in the range of its weak absorption, by using EM in comparison with exclusively DM-based parametrization or SE. In this study, EM for a spectrum T(λ) is abbreviated as EMT, and EM for R(λ) is denoted as EMR, whereas the ultimately computed characteristics of the film are designated by the subscript “c”.
The most popular EMT is the method of Swanepool for uniform film [54] and for non-uniform film [55], both developed assuming a finite transparent thick planar substrate with n(λ) > ns(λ) > 1. Our group enhanced this EMT to account for substrate absorption, as ∆dc and the interval l  [l1c, l2c] of tangency wavelength numbers participating in the characterization were deduced based on observation [56]. Furthermore, we improved EMT by developing an optimizing envelope method for T(λ), abbreviated as OEMT [57,58,59]. OEMT employs external smoothing of the noisy interference spectrum T(λ), which generates externally smoothed interference spectrum Tsm0(λ) touching externally T(λ), instead of common smoothing generating Tsm0(λ) passing closer to internally to T(λ). Such externally smoothed Tsm0(λ) is intended to offset the partial coherence of light in the film due to light scattering mainly from the roughness of the surface film/air, and Tsm0(λ) is slit width corrected thus providing a spectrum Tsm(λ) used thereafter in the film characterization. In essence, OEMT computes optimized values d ¯ c of the average thickness and ∆dc of the thickness non-uniformity of the film over the light spot via the minimization of a selected error metric ERM for wavelengths λ within an optimized interval [λt(l2c), λt(l1c)] corresponding to an optimized interval [l1c, l2c].
The computed results derived from our previous EMT and OEMT characterizations of the thickness parameters of thin films of amorphous semiconductors are displayed in Table 1. The same table also includes the relative error RE( d ¯ c) of the computed average film thickness d ¯ c over the light spot, which corresponds to ERM representing the root mean square error of the estimations of the film thickness calculated for every l from the interval [l1c, l2c].
It is seen from the data about EMT characterizations and RE( d ¯ c) in Table 1 that the accuracy of the computed average film thickness d ¯ c of the uniform film tends to decrease when the substrate is absorbing in comparison with the transparent substrate. Also, for EMT characterizations, the accuracy of d ¯ c apparently decreases even more for the non-uniform film. The data about RE( d ¯ c) from Table 1 also reveal that using OEMT characterization, with fixed l1 = 1, even of a non-uniform film on absorbing substrate leads to the accuracy of d ¯ c commensurate with that of EMT characterization of a uniform film on a transparent substrate.
Furthermore, the data from the fourth and fifth text lines in Table 1 indicate that OEMT provides d ¯ c with significantly smaller RE( d ¯ c) compared to its respective EMT. It is also seen that using OEMT with external smoothing and optimized ∆dc and [l1c, l2c] results in the relative error RE( d ¯ c) < 0.1% in the computation of the average film thickness. This very high accuracy of d ¯ c is due to the very high accuracy of the optimized Δdc and [l1c, l2c] because d ¯ c, ∆dc, and [l1c, l2c] are obtained at the same step of the OEMT algorithm. In addition to this, it is sensible that employing very accurate d ¯ c and ∆dc can lead to the computation of very accurate spectral dependencies nc(λ) and kc(λ) of the film. With respect to the above, it is not surprising that our comparative study [62] showed that OEMT ensures the most accurate characterization of a-Si films, with significantly different average thicknesses, amongst four methods selected as most likely to provide an accurate characterization of such films. Nevertheless, neither external smoothing of interference spectra nor the use of an optimized interval [l1c, l2c] of tangency wavelength numbers are available in commercial software for spectroscopic ellipsometry [63,64] or spectrophotometry [65,66].
EMR for thin film characterization only from R(λ) was proposed in [67] for uniform film on a finite transparent thick planar substrate with n(λ) > ns(λ) > 1, whereas this EMR, also known as the method of Minkov, has been utilized by research groups in several countries [68,69,70,71,72]. Moreover, some modifications of EMR from [67] have been dedicated to the characterization of semiconductor or dielectric thin film from R(λ), assuming that the film has uniform thickness and the thick planar substrate is transparent [73,74,75,76,77]. The non-uniformity Δd of the film has been included in formulae for such R(λ) in [78,79,80], however in none of these papers have been taken into account absorption in the substrate, external smoothing of R(λ), or optimization of d ¯ c, ∆dc, or [l1c, l2c]. As already described, though, accounting for these factors in OEMT results in increasing the accuracy of characterization of a thin film on a thick planar substrate with n(λ) > ns(λ) > 1. Considering the above, the main goal of this study is to develop and use an optimizing envelope method for the characterization of a thin semiconductor or dielectric film on a finite and generally non-transparent thick planar substrate, only from quasi-normal incidence R(λ) of the sample. This method, abbreviated as OEMR, has to resolve the above-mentioned deficiencies of the approaches from [67,68,69,70,71,72,73,74,75,76,77,78,79,80]. Another goal is to determine the most accurate method for the characterization of a-Si films with dissimilar thicknesses, only from UV/Vis/NIR spectrum R(λ), by comparing the accuracies of such characterizations employing different DM-based methods.

2. Materials and Methods

2.1. Theory

A formula about the reflectance Ru(λ) of a uniform thin film on a finite and generally non-transparent (with ks > 0 for some λ) thick planar substrate was derived in [81] via flow graph theory, assuming coherent interaction between light reflected from the two surfaces of the film and incoherent interaction between light reflected from the two surfaces of the substrate. In the case of quasi-normal incidence of light to the film, this formula is rewritten as follows:
R u ( λ ) = ρ a , s 2 + ( n s 2 + k s 2 ) ( τ a , f 2 τ f , s 2 ρ s , a x s ) 2 1 ( ρ s , a ρ s , a x s ) 2 ,
  • whereas t ˙ a , f = τ a , f exp ( i ξ a , f ) = 2 N ˙ a / ( N ˙ a + N ˙ f ) ,   t ˙ f , s = τ f , s exp ( i ξ f , s ) = 2 N ˙ f / ( N ˙ f + N ˙ s ) ,
  • r ˙ a , f = ρ a , f exp ( i ψ a , f ) = ( N ˙ a N ˙ f ) / ( N ˙ a + N ˙ f ) ,   r ˙ f , s = ρ f , s exp ( i ψ f , s ) = ( N ˙ f N ˙ s ) / ( N ˙ f + N ˙ s ) ,
  • r ˙ s , a = ρ s , a exp ( i ψ s , a ) = ( N ˙ s N ˙ a ) / ( N ˙ s + N ˙ a ) ,
  • r ˙ a , s = ρ a , s exp ( i ψ a , s ) = r ˙ a , f + r ˙ f , s ς ˙ 2 1 + r ˙ a , f r ˙ f , s ς ˙ 2 ,   r ˙ s , a = ρ s , a exp ( i ψ s , a ) = r ˙ f , s r ˙ a , f ς ˙ 2 1 + r ˙ a , f r ˙ f , s ς ˙ 2 ,   x S = exp ( 4 π k S d S / λ ) ,
  • φ = 4 π n d / λ ,   α = 4 π k / λ ,   x = exp ( 4 π k d / λ ) = exp ( α d ) ,   x ¯ = exp ( α d ¯ ) ,
  • ς ˙ = exp ( i 2 π N ˙ d / λ ) = exp ( 2 π k d / λ ) exp ( i 2 π n d / λ ) = x exp ( i φ / 2 ) ,
the subscripts “a” and “f” refer to air and film, respectively, and the superscript ’ indicates backpropagation of light. The expression of Ru(λ) in Equation (3) using complex numbers is convenient for development of computer code for optical characterization of thin films. The incoherent interaction between light reflected from the two surfaces of the substrate can be confirmed via absence of interference extrema in reflectance spectrum Rs(λ) of the bare substrate (without film on it).
Furthermore, the reflectance spectrum R(λ) of non-uniform thin film, with uniformly distributed thickness d  [ d ¯ − ∆d, d ¯ + ∆d] over the light spot, on a finite and generally non-transparent thick planar substrate can be obtained by numerical integration of Ru [78,82]:
R ( λ ) = 1 φ 2 φ 1 φ 1 φ 2 R u d φ ,
where φ 1 = 4 π n ( d ¯ Δ d ) / λ , φ 2 = 4 π n ( d ¯ + Δ d ) / λ . Equation (4) takes into account both the substrate absorption and the finite size of the thick substrate, which represents a novelty in spectrophotometry of thin non-uniform films.
On the other hand, envelopes can be computed for every interference reflectance spectrum with at least three apparent maxima and three apparent minima, whereas the envelope along the maxima is designated as R+(λ), and the envelope along the minima is R(λ) [83]. Also, the tangency wavelengths λt between the smoothed and slit width corrected reflectance spectrum Rsm(λ) and its envelopes R+(λ) and R(λ) should satisfy the interference fringes equation. In the region of weak and medium absorption in the film, and n(λ) > ns(λ) > 1, which is commonly valid for a semiconductor or dielectric film, e.g., on a glass substrate [84,85], the interference fringe equation is expressed as follows [78,82]:
2 n ( λ t ) d ¯ = m l ( λ t ) λ t ( l ) where { m l 1 / 2 half - integer   for   all   tangency   wavelengths   λ t ( l )   from   the   envelope R + ( λ ) , m l 1 integer   for   all   tangency   wavelengths   λ t ( l )   from   the   envelope R ( λ ) ,
l = 1, 2, … lM is the tangency wavelength number counted from the longer wavelengths end of R(λ), and m l ( λ t ) is the interference order.
Moreover, the following expression calculates the estimated film thickness:
d 1 ( l ) = d ¯ [ λ t ( l ) ] = λ t ( l ) λ t ( l + 1 ) 4 [ λ t ( l ) n ( l + 1 ) λ t ( l + 1 ) n ( l ) ]
Equation (6) is obtained by rewriting Equation (5) for a pair of adjacent tangency wavelengths λt(l) and λt(l + 1) corresponding to a pair of tangency points Rsm[λt(l)] and Rsm[λt(l + 1)] between Rsm(λ) and its envelopes (one of these points is from the envelope R+(λ), and the other is from the envelope R(λ)) [45,82]. Therefore, the average film thickness d ¯ over the light spot can be approximated using Equations (5) and (6) in the region of weak and medium absorption in the film.
Taking into account Equation (5), the variable φ from Equation (3) can be expressed as follows for the envelope R+(λ):
φ = 4 π n [ d ¯ + ( d d ¯ ) ] / λ from   Equation   ( 5 ) for R + ( λ ) = 2 π . ( half - integer ) + 4 π n ( d d ¯ ) / λ .
Since φ participates in Equations (3) and (4) only via the functions sin( φ ) and cos( φ ) (as seen in the lines of Equation (3), including φ and ς ˙ ), which are periodical functions with a period of 2π, a formula for the upper envelope R+(λ) of Rsm(λ) is derived by substituting Equation (7) in Equation (4) as follows:
R + ( λ ) = 1 φ 2 + φ 1 + φ 1 + φ 2 + R u ( φ + ) d φ + = 1 Δ φ φ 1 + φ 2 + R u ( φ + ) d φ + ,
where
φ + = 4 π n ( d d ¯ ) / λ + π ,   φ 1 + = 4 π n Δ d / λ + π ,   φ 2 + = 4 π n Δ d / λ + π ,   Δ φ = φ 2 + φ 1 + = 8 π n Δ d / λ .
Similarly, φ can be expressed, also from Equation (5), for the envelope R_(λ) as follows:
φ = 4 π n [ d ¯ + ( d d ¯ ) ] / λ from   Equation   ( 5 ) for   R _ ( λ ) = 2 π . ( integer ) + 4 π n ( d d ¯ ) / λ .
Correspondingly, a substitution of Equation (9) in Equation (4) provides a formula for the lower envelope R_(λ) of Rsm(λ) as follows:
R ( λ ) = 1 φ 2 φ 1 φ 1 φ 2 R u ( φ ) d φ = 1 Δ φ φ 1 φ 2 R u ( φ ) d φ ,
where
φ = 4 π n ( d d ¯ ) / λ ,   φ 1 = 4 π n Δ d / λ ,   φ 2 = 4 π n Δ d / λ ,   Δ φ = φ 2 φ 1 = 8 π n Δ d / λ .
Furthermore, the following analytic approximations of the envelopes R+(λ) and R(λ) are obtained by substituting Ru(λ) from Equation (3) into Equations (8) and (10) and solving the respective symbolic integrals as follows:
R 0 ± ( λ ) 1 [ a 3 0 a 2 0 θ a 3 0 2 b 3 0 2 ( n s 2 + k s 2 ) ( τ a , f 0 2 τ f , s 0 2 ρ s , a x ¯ x s ) 2 b 3 0 θ ( a 1 0 b 3 0 a 3 0 b 1 0 ) a 3 0 2 b 3 0 2 ] tan 1 [ a 3 0 b 3 0 a 3 0 ± b 3 0 tan ( θ ) ] + ( n s 2 + k s 2 ) ( τ a , f 0 2 τ f , s 0 2 ρ s , a x ¯ x s ) 2 b 1 0 θ ( a 1 0 b 3 0 a 3 0 b 1 0 ) a 1 0 2 b 1 0 2 tan 1 [ a 1 0 b 1 0 a 1 0 ± b 1 0 tan ( θ ) ] ,
where
a 1 0 = 1 ( ρ a , f 0 ρ s , a x ¯ x S ) 2 + ( ρ a , f 0 ρ f , s 0 x ¯ ) 2 ( ρ f , s 0 ρ s , a x s ) 2 ,   b 1 0 = 2 ρ a , f 0 ρ f , s 0 x ¯ [ 1 ( ρ s , a x S ) 2 ] , a 2 0 = ρ a , f 0 2 + ( ρ f , s 0 x ¯ ) 2 ,   b 3 0 = 2 ρ a , f 0 ρ f , s 0 x ¯ ,   a 3 0 = 1 + ( ρ a , f 0 ρ f , s 0 x ¯ ) 2 ,   θ = 2 π n Δ d / λ = Δ φ / 4 , τ a , f 0 = 2 ( n + 1 ) ,   τ f , s 0 = 2 n ( n + n s ) 2 + k s 2 ,   ρ a , f 0 = n 1 n + 1 ,   ρ f , s 0 = ( n n s ) 2 + k s 2 ( n + n s ) 2 + k s 2 ,
as the upper sign from the “ ± ” and “ ” symbols corresponds to R+(λ), and the lower of these signs refer to R(λ). Equation (11) is derived by replacing the transmittance x = exp ( α d ) in the film (for light passing once between its two surfaces) with its averaging transmittance x ¯ = exp ( α d ¯ ) over the light spot, and ignoring k(λ) in the Fresnel coefficients denoted by “τ” and “ρ” in Equation (3). The novel Equation (11) makes it possible to expand the EMR framework to account for both the substrate absorption and the finite size of the thick planar substrate, unlike the respective expressions used in EMR from [67,68,69,70,71,72,73,74,75,76,77,78].

2.2. Features of Simulated Quasi-Normal Incidence Interference Spectra R(λ) and Their Envelopes

Several model specimens corresponding to a-Si:H film on a thick planar glass substrate are introduced for representing the typical behavior of interference reflectance spectra R(λ) of a thin semiconductor film on a thick planar substrate. The refractive index and the extinction coefficient of the model film are n ( λ ) = 2.6 + 3 × 10 5 / λ 2 and k ( λ ) = ( λ / 4 π ) × 10 1.5 × 10 6 / λ 2 8 , where λ (nm) [54]. The refractive index ns(λ) of the model substrate is that of a standard microscope slide glass substrate G50 of Levenhuk [83,84], and its extinction coefficient is either ks(λ) > 0 of the same substrate of Levenhuk or ks(λ) = 0. Interference reflectance spectra of the model specimens are calculated from Equation (3) for uniform film and from Equation (4) for non-uniform film, and interference transmittance spectrum Tu(λ) of one model specimen with uniform film is calculated as in [59,86]. Such simulated reflectance spectra and transmittance spectrum, their envelopes and corresponding absorbance spectrum of the model specimens, and other optical characteristics of their films and substrates, are presented in Figure 3a,b,d. Figure 3c shows differences ∆R between reflectance spectra as well as differences ∆R+ and ∆R between their respective envelopes, computed either by numerical integration (NI) with Nst integration steps from Equations (4), (8) and (10) or by the analytic approximations from Equation (11).
Analysis of the graphs from Figure 3 reveals the following dependencies:
(2.2.1)
Ru(λ) + Tu(λ) ≈ 1 in the region of weak absorption in the film because the absorbance of the specimen with transparent substrate is Au(λ) = 1 − Ru(λ) + Tu(λ) ≈ 0 in this region, which is illustrated with the dashed magenta line in Figure 3a. Therefore, the interference patterns of R(λ) and T(λ) have very similar features in the region of weak absorption in the film. This fact indicates that advances used for accurate characterization of a thin film from T(λ) (such as those of OEMT), in the region of weak absorption in the film, can be readily employed for accurate characterization of a thin film from R(λ). Moreover, it is seen also from Figure 3a that the reflectance spectrum of the specimen shrinks, and its lower envelope Ru−(λ) drifts above the reflectance Rs(λ) of the bare substrate with decreasing λ in the region of medium absorption in the film, where Au(λ) rises.
(2.2.2)
The lower envelope Ru−(λ) of R(λ) almost coincides with Rs(λ) in the region of weak absorption in uniform film (Figure 3a), and R(λ) is positioned above Rs(λ) in the region of weak absorption in non-uniform film (Figure 3b). However, such a drift of R(λ) above Rs(λ) can be due to both thickness non-uniformity of the film and absorption in the film, taking into account the last comment from (2.2.1). On the other hand, it was indicated in the introduction that OEMT renders optimized values of both the average thickness d ¯ and the thickness non-uniformity ∆d of the film based on analysis for its region of weak absorption. As a result of the above and the principle similarity with OEMT, it is expected that employing OEMR can lead to superior accuracy characterization of semiconductor and dielectric films, only from R(λ), over the entire UV/Vis/NIR R(λ).
(2.2.3)
In the region where the averaging transmittance in the film is x ¯ = exp ( α d ¯ ) ≈ 0, the film is opaque, the envelopes R+(λ) and R(λ) merge (as seen in Figure 3b), the incident light is reflected only at the surface film/air, and R(λ) = ρa,f2(λ) = [(n − 1)2 + k2]/[(n + 1)2 + k2]. Therefore, the interference-free part of R(λ) supplies information about n(λ) and k(λ), which can be used for their accurate determination in this region, unlike the transmittance spectrum because T(λ) ≈ 0 there.
(2.2.4)
Regarding Figure 3c, the absolute values of all differences in reflectance, calculated by NI with Nst = 100 and Nst = 30, do not exceed 10−4. Since 10−4 is quite a small value compared to R(λ) and the maximum difference in thicknesses d over the light spot of this model film is 2∆d = 60 nm, all NI in Equations (4), (8) and (10) are executed using Nst ≥ ∆d(nm) ∗ 30/60 = ∆d(nm)/2 from here on in this paper. Furthermore, the absolute values of the differences ∆R+ and ∆R, calculated by NI with Nst = 100 and the analytical approximation from Equation (11), do not exceed 2 × 10−5 over the region of weak absorption in the film, which exemplifies the accuracy of Equation (11).
(2.2.5)
The lower envelope R(λ) of R(λ) is more dependent on the substrate characteristics than the upper envelope R+(λ), as comprehended from Figure 3d. Correspondingly, the envelope R+(λ) should be more dependent on the film characteristics than the envelope R(λ), which indicates that employing R+(λ) is more likely to provide accurate film characteristics compared to R(λ). The appearances of R(λ) and its envelopes in Figure 3d also clarify that presence of only two apparent extrema of R(λ) at each of its envelopes should not be sufficient for precise calculation of R+(λ) and R(λ).
(2.2.6)
There are notable differences in R(λ) and its two envelopes for ks > 0 and ks = 0, respectively, in the region of significant absorption in the substrate, enclosed by a red colored rectangle in Figure 3b,d. Therefore, taking into account the absorption in the substrate should result in increasing the accuracy of characterization of a thin film only from R(λ).

2.3. The Algorithm of OEMR and Its Details

The OEMR is developed for the common case of quasi-normal incidence of light to a specimen consisting of a thin semiconductor or dielectric film on a thick planar substrate with n(λ) > ns(λ) > 1. The algorithm of OEMR is based on the algorithms of EMR from [78] and OEMT from [57,58], and its steps are shown in Figure 4.
In this present study, at step A2 of the algorithm can be independently employed three different kinds of initial smoothing of R(λ) (needed for removing its false extrema due to noise), providing a smoothed spectrum Rsm0(λ) in the region of the interference pattern. More specifically, “internally smoothed” Rsm0(λ) touches internally the noisy spectrum R(λ), “externally smoothed” Rsm0(λ) touches externally R(λ), and “medium smoothed” Rsm0(λ) passes in the middle of the previous two.
The smoothed spectrum Rsm(λ) is obtained, at step A3, using a slit width correction of Rsm0(λ), based on the fact that a similar correction has been employed for T(λ) in [54,57] as follows:
R sm ( λ t ) = R sm 0 ( λ t ) ± [ R sm 0 ( λ t )   Δ s Δ l ( λ t ) ] 2
where Δs is the spectral slit width, and Δ l ( λ t ) = Δ l [ λ t ( l ) ] = λ t ( l 1 )     λ t ( l + 1 ) is the line width at λt(l), as “+” from the “ ± ” symbol refers to a maximum of Rsm0(λ) and “−” to its minimum. In general, slit width correction is needed for thicker film samples with some linewidths Δl(λt) < 10Δs in their R(λ) [87]. At step A4 of the algorithm, the envelopes R+(λ) and R(λ) of Rsm(λ) are computed using left and right “boundary points”, as described in [88,89] as well as the interpolation procedure from [58], for each of the two envelopes. The envelopes around internally smoothed Rsm(λ) are named “internal envelopes of R(λ)”, the envelopes around externally smoothed Rsm(λ)—“external envelopes of R(λ)”, and the envelopes around medium smoothed Rsm(λ)—“medium envelopes of R(λ)”.
Part B of the algorithm of OEMR, presented in Figure 5, performed optimization, providing optimized results for the interval [l1c, l2c] of tangency wavelength numbers, the lowest interference order m1c = ml[λt(l = 1)] = ml(l = 1), the thickness non-uniformity ∆dc of the film, and its average thickness d ¯ c over the light spot.
With respect to steps B1 and B12 of the algorithm, the performance of seven error metrics (ERM) for employment in OEMT from [58,59] (only for T(λ)) have been investigated in [90]. It has been concluded there that the presented below error metrics ERM1 and ERM2 provide the most accurate OEMT characterization for a wide variety of semiconductor and dielectric thin films on a finite generally non-transparent thick planar substrate:
E R M 1 ( i 1 , l 1 , l 2 ) = 1 l 2 l 1 l = l 1 l 2 1 [ m l ( i 1 , l ) m e ( i 1 , l ) ] 2 l 2 l 1 = 1 N l ( l 1 , l 2 ) 1 l = l 1 l 2 1 [ m l ( i 1 , l ) m e ( i 1 , l ) ] 2 l 2 l 1 = R M S E ( m e ) N l 1 0 , E R M 2 ( i 1 , l 1 , l 2 ) = 1 l 2 l 1 + 1 l = l 1 l 2 [ d ¯ e ( i 1 , l 1 , l 2 ) d 2 ( i 1 , l ) ] 2 l 2 l 1 + 1 = 1 N l ( l 1 , l 2 ) l = l 1 l 2 [ d ¯ e ( i 1 , l 1 , l 2 ) d 2 ( i 1 , l ) ] 2 l 2 l 1 + 1 = σ ( d ¯ e ) N l 0 , E R M i ( i 1 , l 1 , l 2 ) = E R M 1 ( i 1 , l 1 , l 2 )   or   E R M 2 ( i 1 , l 1 , l 2 ) ,
where RMSE means root mean square error, σ designates standard deviation, and N l ( l 1 , l 2 ) is the number of tangency wavelengths in the interval [l1, l2]. The suitability of ERMi (from Equation (13)) for OEMR is justified via the validity of the approximation T(λ) + R(λ) ≈ 1 in the region of weak absorption in the film, as illustrated in Figure 3a and described in paragraph (2.2.1). A principle advantage of ERMi is that it provides a wider interval [l1, l2], over which the computed reflectance spectrum Rc(λ) (obtained via substitution of nc(λ), kc(λ), d ¯ c, and ∆dc in Equation (4)) represents a good approximation of the measured R(λ) compared to, e.g., the error metric σ( d ¯ e) from the numerator of ERM2.
Moreover, the error metric used at step B15 of the algorithm is expressed as follows:
E R M a ( l 1 ) = min { E R M i [ i 1 ( l 1 ) , fixed   l 1 , l 2 ( l 1 ) ] } and   their   respective   provides   i 1 ( l 1 ) , l 1 , l 2 ( l 1 ) m 1 a ( l 1 ) = m e [ i 1 ( l 1 ) , l = 1 ] ,   Δ d a ( l 1 ) = Δ d [ i 1 ( l 1 ) ] ,   d ¯ a ( l 1 ) = d ¯ e [ i 1 ( l 1 ) , fixed   l 1 , l 2 ( l 1 ) ] .
Furthermore, the relative error in the computed average film thickness d ¯ c is expressed by ERMc for ERM2, obtained at step B17, as follows:
R E ( d ¯ c ) = σ ( d ¯ c ) d ¯ c = E R M c ( E R M 2 ) ( l 2 c l 1 c + 1 ) d ¯ c 0 and   its   respective provides   l 1 c , l 2 c ( l 1 c ) [ λ t ( l 1 c ) ,   λ t ( l 2 c ) ] [ l 1 c ,   l 2 c ] , m 1 c = m 1 a ( l = 1 ) ,   Δ d c = Δ d a ( l 1 c ) ,   d ¯ c = d ¯ a ( l 1 c ) .
The optimization in part B of the algorithm of OEMR is achieved by a direct comparison between the error metric values corresponding to all credible sets (l1,l2d) rather than using a standard minimization algorithm, which can terminate in a local minimum of the error metric. This is carried out by including an independent loop for every one of the parameters l1, l2, and Δd, as seen in Figure 5. There is, therefore, only one possible set of thickness parameters ∆dc and d ¯ c of the film, corresponding to the global minimum of ERMi(l1,l2d), as these thickness parameters are computed at the end of part B of the algorithm. Notably, the algorithm of OEMR is validated for the model a-Si:H film with n(λ) and k(λ) introduced in the first paragraph from Section 2.2, d ¯ = 1000 nm and Δd = 30 nm on 3.28 mm thick Borofloat33 substrate, as it turns out that the errors in both d ¯ c and ∆dc do not exceed 0.025%.
Furthermore, the computed spectral dependencies nc(λ) for the refractive index and kc(λ) for the extinction coefficient of the film can be determined, as described in part C of the algorithm of OEMR from Figure 4, similarly to using OEMT [57,58]. Moreover, nc(λ) and kc(λ) can be also derived via some existing broad spectrum dispersion models, however employing the already accurately computed via OEMR thickness parameters d ¯ c and ∆dc of the film. After d ¯ c, ∆dc, nc(λ), kc(λ), and their respective spectrum Rc(λ) from Equation (4) are computed, for all wavelengths from the measured R(λ), the figure of merit is calculated as follows:
F O M = 1000 × j = 1 N j { R [ λ ( j ) ] R c [ λ ( j ) ] } 2 N j = 1000 × R M S E ( R c ) 0
Equation (16) can be utilized as a measure of closeness of Rc(λ) to R(λ), where Nj is the number of all λ  [min(λ), max(λ)]. Therefore, the film characterization providing smallest FOM is considered to be the most accurate characterization of that film.

3. Results

3.1. Summary of the Preparation of Samples and the Measurement of Their Reflectance Spectra

a-Si thin films with dissimilar thicknesses were deposited on non-transparent glass substrates via RF-magnetron sputtering using applied RF power of 525 W, target-to-substrate distance of 6.1 cm, and Ar gas. Sample SP1 consists of a-Si film sputtered at Ar pressure of 4.4 Pa on 0.9 mm thick Corning7059 substrate, and the sample SP2 represents a-Si film sputtered at Ar pressure of 0.13 Pa on 3.28 mm thick Borofloat33 substrate. The deposition rate was 73.8 nm/min for the film SP1 (from the sample SP1) and 90.0 nm/min for the film SP2 (from the sample SP2). The reflectance spectra R(λ) of these samples were measured by a Perkin-Elmer Lambda 1050 UV/visible/NIR double-beam spectrophotometer using specular reflectance accessory B0086703 with Φ = 6° [91,92], 2 nm slit width, 1 nm data collection step, and 10 mm × 3 mm illuminated area of the film. The a-Si thin films SP1 and SP2 are characterized next, only from their corresponding R(λ), using OEMR and several dispersion models.

3.2. Characterizations of the a-Si Films by OEMR

The two a-Si films SP1 and SP2 are characterized via OEMR employing internal envelopes, external envelopes, or medium envelopes of their reflectance spectra R(λ). The results from the execution of parts A and B from the algorithm of OEMR for these films are presented in Table 2. Data from OEMR characterizations using the deconvolution of these R(λ) and their respective deconvolution envelopes are also included in Table 2. The deconvolutions are performed assuming Gaussian distribution, with spectral slit width Δs, of the spectrum of light incident on the film.
It is seen from the last two columns of Table 2 that the most accurate results for d ¯ c are achieved by employing the external envelopes of R(λ), in which case, the differences in d ¯ c and ∆dc obtained via the error metrics ERM1 and ERM2 do not exceed 0.06% and 0.9%, respectively, for both studied films. Therefore, the optimized thickness parameters d ¯ c and ∆dc of a given film are derived by using OEMR with the external envelopes of R(λ), and the values of these parameters are selected to be equal to half of the sums of their corresponding values obtained via ERM1 and ERM2.
Reflectance spectra, including illustrations of the smoothing and the boundaries of the intervals [l1c, l2c], corresponding to the optimized thickness parameters d ¯ c and ∆dc from Table 2, are shown in Figure 6, for both studied samples. The dashed red lines in Figure 6c,d represent the computed reflectance spectra Rc(λ) corresponding to the characterizations of the films SP1 and SP2 via the synthetic dispersion model “OEMR with 3NA” described in Section 3.5.
Taking into account the interference fringes in Equation (5), the significantly larger number of interference extrema of R(λ) in Figure 6d compared to Figure 6c confirms that the film SP2 is much thicker than the film SP1. Furthermore, the larger difference R(λ) − Rs(λ) > 0 from Figure 6d than from Figure 6c reaffirms that the film SP2 has larger thickness non-uniformity ∆d compared to the film SP1, according to the description from paragraph (2.2.2). The fact that l1c > 1, illustrated in both Figure 6c,d, indicates that OEMR tends to disregard a small number of the largest wavelengths extrema of Rsm(λ). This should be related to the unavoidable small inaccuracy of the right boundary points R+[max(λ)] and R[max(λ)] from the envelopes R+(λ) and R(λ).
Images regarding the execution of part C from the OEMR algorithm are presented in Figure 7 for both studied films. Each of the Wemple–DiDomenico plots from Figure 7a,c includes green circles corresponding to all nc(λt) calculated at step C1 from the OEMR algorithm, as well as a representation of their linear regression performed over a selected interval [l1w, l2w] and depicted by a solid green line. The computed energy E0 and strength Ed of the respective undamped single oscillator, as well as the static refractive index n(0), are also printed in these plots, whereas the expansion of the linear regression beyond the interval [l1w, l2w] is illustrated with a dashed magenta line. The computed refractive index nc(λ) is derived via such an expansion of the linear regression, for wavelengths λ within the interval [λt(l2w), max(λ)], as described in step C2 from the OEMR algorithm. Furthermore, the extinction coefficients kc(λ) from Figure 7b,d are computed, for wavelengths λ also within the interval [λt(l2w), max(λ)], using nc(λ) as well as Equations (8) and (10) about the envelopes R+(λ) and R(λ). The figures of merit FOM, printed in Figure 7b,d, are calculated using Equation (16), however applied to the above interval of λ.
It is concluded from Figure 6c,d and Figure 7a,c and Equation (1) that DM of Wemple–DiDomenico, corresponding to the undamped single oscillator, describes well n(λ) in the interval [l1w, l2w], representing the region of weak to medium absorption in the a-Si films. In Figure 7a, the small shifts of the green circles from the dashed magenta line for l < l1w = 3 should be related to the already mentioned small inaccuracy of the envelopes R+(λ) and R(λ) around their right boundary points R+[max(λ)] and R[max(λ)]. Moreover, the drift of the green circles away from the dashed magenta line for l > l2w is due to the inaccuracy of DM of Wemple–DiDomenico, which does not include k(λ), in the region of strong to medium absorption in the film. Furthermore, the smallest FOMs correspond to the green curves from Figure 7b,d, which indicates that k(λ) is computed most accurately from R+(λ) and its respective Equation (8) compared to R(λ) or [R+(λ) + R(λ)]/2. Nevertheless, the inaccuracy of representation of n(λ) via the linear regression for l > l2w from the Wemple–DiDomenico plot indicates that the accuracy of computation of k(λ) from R+(λ) should decrease with decreasing λ in the region of strong to medium absorption in the film. Based on the above comments from this paragraph, it is expected that the accuracy of film characterization over the entire UV/Vis/NIR R(λ) can be increased by combining OEMR with DM describing more accurately N ˙ ( λ ) of the film in its region of medium and strong absorption.
Some amorphous materials can contain voids [93,94], and the voids volume fraction (with respect to the entire volume of the material) has been approximated as follows:
f void     [ 1 + 2 n 2 ( 0 ) ] [ n dense 2 ( 0 ) n 2 ( 0 ) ] 3 n 2 ( 0 ) [ n dense 2 ( 0 ) 1 ]
where ndense(0) is the static refractive index of the same amorphous material but without voids [60,93]. On the other hand, E0 ≈ 4.0 and Ed ≈ 44.4 for crystalline Si [95], whereby n(0) ≈ 3.48 is obtained from Equation (1) applied to this material, while the difference between ndense(0) for a-Si without voids and n(0) for crystalline Si is ≈ 0.28 based on data from [96]. Correspondingly, the static refractive index of a-Si without voids is estimated to be ndense(0) ≈ 3.48 + 0.28 = 3.76. Furthermore, replacing the values of n(0), printed in Figure 7a,c, in Equation (17), provides fvoid ≈ 5.4% for the film SP1 and fvoid ≈ 0.4% for the film SP2. Therefore, the somewhat larger value of kc(λ) in the region of weak absorption in the film SP1 (compared to the film SP2) should be due to light reflections from boundaries of such voids, which create incoherent additive to R(λ), thus slightly smearing its interference pattern. Moreover, effective thickness filled with voids in a film can be expressed as d ¯ void = d ¯ c fvoid, which leads to d ¯ void = 36.4 nm for the film SP1 and d ¯ void = 15.4 nm for the film SP2.

3.3. Characterizations of the Films via Single Electron Oscillator DMs

The a-Si films are characterized using broad spectrum single electron oscillator DMs 1TL, 1NA, TLUR, and TLUF as well as Equation (4), which includes implicitly the film thickness parameters d ¯ and ∆d. The computed results from these characterizations are presented in Table 3, where εr(∞) and nc(∞) are the complex dielectric function and the complex refractive index of the film for photon energy E→∞, assuming that kr(∞) = 0. The first three columns of this Table present the values of the oscillator parameters, as the first of them corresponds to the oscillator amplitude, the second is the central energy, and the third is the broadening term. Each symbol of these parameters contains subscripts, of which the first two letters represent the employed DM, and the number after them is the number of the oscillator, i.e., 1 for the DMs featured in Table 3. The initial approximations of the oscillator parameters are selected to be in the midst of the intervals of recommended values for these parameters prescribed by their corresponding DM. In this paper, all film characterizations employing broad spectrum oscillator DMs are performed via the minimization of FOM from Equation (16), separately utilizing Nelder–Mead simplex direct search algorithm and interior point algorithm [97,98]. The differences are recorded between the values of FOM at every two successive steps of the minimization, and the respective film characterization is considered to be completed when the absolute values of these differences fall below 0.001 for both of these algorithms.
The data from Table 3 indicate that TLUF provides the most accurate characterization of both a-Si films (i.e., characterization with the smallest FOM) amongst the four broad-spectrum single oscillator DMs featured in this Table. In particular, it is seen from FOM in Table 3 that the characterizations of both a-Si films using TLUR are less accurate than those via TLUF. The main reason for this is that TLUR employs the formula εi(E) = const × E × exp(E/EU) for the film in the entire range EEb [29], which, however, is inaccurate in the region of weak absorption from the measured UV/Vis/NIR spectrum R(λ). Instead, in this region εi(E) = 2n(E)k(E) ~ k(E) = αλ/(4π) = 1239.8α/(4πE) const/E × exp(E/EU), taking into account the expression about α from Equation (3) and assuming the existence of Urbach tail in this region represented by Equation (2), whereas exactly the last formula for εi(E) is used in TLUF [27].
It was indicated in the introduction that nc(λ) and kc(λ) are expressed in the framework of every broad-spectrum electron oscillator DM as functions only of the oscillator parameters. Correspondingly, the computed reflectance spectrum Rc(λ) for such DM is obtained via the substitution of the computed film thickness parameters d ¯ c and ∆dc and the oscillators parameters in Equation (4). Therefore, a larger absolute value of the difference Rc(λ) − R(λ) in a particular spectral region indicates a larger inaccuracy of its respective DM in this region. With regard to this, the differences Rc(λ) − R(λ) for the single oscillator DMs 1TL, 1NA, TLUR, and TLUF are computed, using the data from Table 3, and illustrated in Figure 8a,b, for both a-Si films.

3.4. Characterizations of the Films via Multiple Electron Oscillator DMs

Aiming to resolve the above-mentioned problem of the insufficient accuracy of characterization via a single oscillator DM in the region of medium absorption, both films are also characterized based on using three electron oscillator DMs 3TL and 3UD. To additionally reduce FOM, two single electron oscillators, 1PE and 1GA, are added to 3TL, as in [25], whereas the initial approximations of the central energies of 1PE and 1GA are above and below those of 3TL, respectively; the resultant DM being abbreviated as “3TL with 1PE and 1GA”. Moreover, Urbach tail and 1GA are added to 3UD, as in [31], and this DM is abbreviated to “3UD with UT and 1GA”. The results from the characterizations of the a-Si films via 3TL with 1PE and 1GA, as well as via 3UD with UT and 1GA are included in Table 4.
The a-Si films are also characterized via DMs, including two electron oscillators 2TL or 2UD as well as 1GA and 1PE or UT, whereas these DMs are abbreviated to “2TL with 1PE and 1GA” and “2UD with UT and 1GA” (they are simplifications of “3TL with 1PE and 1GA” and “3UD with UT and 1GA”). It turns out, though, that 2TL with 1PE and 1GA, as well as 2UD with UT and 1GA, are also insufficiently accurate in the region of medium absorption, for both a-Si films, similar to the single oscillator DMs 1TL, 1NA, TLUR, and TLUF, whose performance is featured in Table 3. Therefore, the results from the characterization of the a-Si films via 2TL with 1PE and 1GA, as well as via 2UD with UT and 1GA are not included in this paper.

3.5. Characterizations of the Films via a Synthetic DM including OEMR

This novel synthetic DM represents a constrained combination of OEMR for longer wavelengths and an existing electron oscillator DM for shorter wavelengths. It employs d ¯ c and ∆dc computed at step B17 of the algorithm of OEMR, nc(λλw) is calculated from the Wemple–DiDomenico plot at step C2 of the algorithm, where λw(nm) = 1239.8/Ew(eV), and kc(λλc) is computed from the envelope R+(λ) at step C3 of the algorithm, where λc(nm) = 1239.8/Ec(eV). Both nc(λλw) and kc(λλc) are obtained via the electron oscillator DM, whereas their respective nc(λλw) and kc(λλc) are constrained to be identical with those for OEMR. The electron oscillator DM for shorter wavelengths, employed for the characterization of the as-Si films, is selected to be 3TL or 3NA (either of them containing three oscillators), thus having similarity to the already used 3TL with 1PE and 1GA as well as 3UD with UT and 1GA. The corresponding synthetic DMs are abbreviated as “OEMR with 3TL” and “OEMR with 3NA”. Computed parameters from characterizations of the a-Si films via OEMR with 3TL and OEMR with 3NA are also included in Table 4. Notably, in all four DMs featured in Table 4 are used the physically plausible approximations nc(∞) = 1 and kc(∞) = 0, as indicated, e.g., in [25,31].
A comparison of the FOM data from Table 3 and Table 4 shows that the smallest FOMs are achieved in the characterizations via the synthetic dispersion model OEMR with 3NA, for both films SP1 and SP2, thus indicating that these characterizations are the most accurate ones for these films. For this reason, the computed reflectance spectra Rc corresponding to the characterizations of the films SP1 and SP2 via OEMR with 3NA are included as dashed red lines in Figure 6c,d. Also, the fact that FOM is smaller for the characterization of either of the films SP1 and SP2 via OEMR with 3NA compared to OEMR with 3TL (as seen from Table 4) indicates that the appearances of n(λλw) and k(λλc) of these films can be explained using 3NA oscillators rather than 3TL oscillators.
Moreover, in Figure 8c,d are also included the differences Rc(λ) − R(λ) corresponding to the characterizations of the a-Si films via the four DMs featured in Table 4. Furthermore, the extinction coefficient k(λ) of the film can be expressed via its absorption coefficient α(E), whereas k(λ) = λα(λ)/(4π) = 1239.8α(E)/(4πE) from Equation (3), and a logarithm from α(E) would be used for the representation of the Urbach tail according to Equation (2). With respect to the above, the refractive index nc(λ) and log10αc(E), computed using the four DMs featured in Table 4, are shown in Figure 9. Taking into account that the smallest FOMs in the characterizations of the a-Si films are achieved by employing OEMR with 3NA (see Table 3 and Table 4), the most accurate nc(λ) and log10αc(E) are shown in red colors in Figure 9.
A comparison between Figure 6a,c and Figure 9a,c, shows that the refractive index of the film SP1 is smaller than that of the film SP2, in the region of weak absorption. This is undoubtedly due to the presence of voids with a larger volume fraction fvoid ≈ 5.4% for the film SP1 compared to fvoid ≈ 0.4% for the film SP2, as described in the paragraph after Equation (17).
It is noticed from Figure 9a,c that nc(λ) corresponding to the characterizations of the a-Si films via OEMR with 3NA, and depicted by the red colored lines, is significantly larger in the region of strong absorption compared to the characterizations based on 3TA or 3UD in this region. This result is in agreement with the comment after Figure 1 that the refractive index n{E = [0,max(εr)]} is larger for 1NA compared to 1TL, 1CC, and 1GA. Moreover, according to Equation (2), the insufficient linearity of the dependencies log10αc(E), presented in red color in Figure 9b,d, for photon energies E below Eg (taken from the data printed in red in Table 4) indicates that the assumption of Urbach tail is inaccurate, especially for the film SP2. A detailed analysis of the optical and electrical characteristics of such a-Si films, based on similar results for nc(λ) and log10αc(E), is intended to be published elsewhere.

4. Discussion

In the characterization of a thin film on a thick planar substrate only from R(λ) via a broad-spectrum DM, the film thickness parameters are computed together with the oscillator parameters by fitting the computed spectrum Rc(λ) to R(λ), although DM represents an approximation of the true dispersion. In the case of interference spectrum R(λ), however, it is possible to achieve higher accuracy of the computed film thickness parameters d ¯ c and ∆dc by employing EMR because it uses the interference fringes Equation (5) in addition to R(λ).
Furthermore, the proposed OEMR represents an improvement of EMR from [67,78]. The advantages of OEMR (compared to EMR) are due to the following: accounting for the finite size and the absorption in the substrate, the external smoothing of R(λ) to offset the partial coherence of light in the film, the slit width correction, the enhanced drawing of both envelopes based on [58,89], and the optimized computed thickness parameters d ¯ c and ∆dc derived in an optimized interval [λt(l2c), λt(l1c)]. In fact, these advantages of OEMR compared to EMR are identical to the advantages of OEMT compared to EMT. With respect to this, the data and comments regarding Table 1 show that such OEMT provides a more accurate computed average film thickness d ¯ c than EMT or a limited version of OEMT, for a-Si, a-As40S60 and a-As98Te2 films. Therefore, the characterization of these films using OEMR should be more accurate than EMR.
Furthermore, RE( d ¯ c) is calculated from Equation (15) for an EMR characterization of amorphous As33S67 film published in Table 2 from [78] (Δd was taken into account there, unlike assuming Δd = 0 in [67]), whereas the result is RE( d ¯ c) ≈ 0.393%. This can be compared with RE( d ¯ c) ≈ 0.099% for the a-Si film SP1 and RE( d ¯ c) ≈ 0.172% for the a-Si film SP2, taken from Table 2 of this paper. These data indicate that employing OEMR leads to decreasing the error in the computation of the average film thickness d ¯ by more than 2.2 times in comparison with using EMR.
Table 5 includes data about FOM and the thickness parameters d ¯ c and ∆dc of both a-Si films already computed via OEMR with 3NA, TLUF (being the best characterization featured in Table 3), and the best characterizations not employing OEMR from Table 4. Taking into account the content of the above three paragraphs, it is assumed tentatively that the parameters d ¯ c and ∆dc computed via OEMR with 3NA (they are the same as those from OEMR) are identical with their respective true thickness parameters d ¯ and ∆d for each of the two a-Si films. Based on these data, we calculated the relative errors REc( d ¯ c) = ( d ¯ c d ¯ )/ d ¯ and REc(∆dc) = (∆dc − ∆d)/∆d corresponding to the characterizations via TLUF and by the best characterizations not employing OEMR from Table 4, as these errors are also included in Table 5. Therefore, REc( d ¯ c) and REc(∆dc) represent the relative errors in the computation of d ¯ and ∆d via TLUF or the best characterizations not employing OEMR from Table 4.
The results from Table 5 indicate that smaller FOM for a particular film characterization only from UV/Vis/NIR R(λ) is associated with smaller errors in the computation of both d ¯ and ∆d, as well as that the errors in the computation of d ¯ are smaller than those for ∆d. Based on the above comments from this section, it is argued that the most accurate thickness parameters d ¯ and ∆d of a thin semiconductor or dielectric film on a thick planar substrate can be computed, from only one quasi-normal incidence UV/Vis/NIR R(λ), by using the proposed here OEMR. Therefore, OEMR should provide the most accurate n(λ) and k(λ) of such film in its region of weak to medium absorption. Notably, OEMR is applicable for any kind of thick planar substrate, independent of its degree of absorption.
Furthermore, utilizing the already computed via OEMR d ¯ and ∆d together with a given broad spectrum oscillator DM should result in the computation of more accurate parameters of its oscillators, which should lead to the computation of more accurate n(λ) and k(λ) of the film. Aiming to inculcate this idea into DM, the proposed synthetic DM includes OEMR for longer wavelengths and an existing electron oscillator DM for shorter wavelengths. Therefore, it is not surprising that the most accurate characterizations of the a-Si films SP1 and SP2 over the entire measured spectrum are performed by employing such synthetic DM, which turns out to be OEMR with 3NA for these films.
It should be also possible to use OEMR for the characterization of optical metasurfaces, such as those reported in [99,100]. However, this would require development of a model of the reflection from the metasurface and a respective equation for its R(λ).

Author Contributions

Conceptualization, D.M. and E.M.; methodology, D.M.; software, D.M. and D.N.; validation, D.M., E.M., and G.A.; formal analysis, D.M.; investigation, R.R. (Rostislav Rusev) and S.R.; resources, S.R. and R.R. (Rossen Radonov); data curation, D.N.; writing—original draft preparation, D.M.; writing—review and editing, E.M.; visualization, R.R. (Rostislav Rusev); supervision, G.A. and R.R. (Rossen Radonov); project administration, G.A. and R.R. (Rossen Radonov); funding acquisition, G.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Regional Development Fund within the Operational Programme “Bulgarian national recovery and resilience plan”, procedure for direct provision of grants Establishing of a network of research higher education institutions in Bulgaria”, under Project BG-RRP-2.004-0005 “Improving the research capacity and quality to achieve international recognition and resilience of TU-Sofia (IDEAS)”.

Data Availability Statement

Data are available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kumar, S.; Aswal, D.K. Recent Advances in Thin Films (Materials Horizons: From Nature to Nanomaterials), 1st ed.; Springer: Berlin/Heidelberg, Germany, 2020; pp. 211–343. [Google Scholar]
  2. Vanalakar, S.; Patil, S.; Deshmukh, S. Solid State Thin Films Deposition and Their Applications: Thin Film Deposition; Lambert Academic Publishing: Saarbrukken, Germany, 2020; pp. 52–147. [Google Scholar]
  3. Stenzel, O.; Ohlídal, M. Optical Characterization of Thin Solid Films, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2016; pp. 93–99. [Google Scholar]
  4. Stenzel, O. Optical Coatings: Material Aspects in Theory and Practice, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2016; pp. 35–51. [Google Scholar]
  5. Tompkins, H.G.; Hilfiker, J.N. Spectroscopic Ellipsometry: Practical Application to Thin Film Characterization, 1st ed.; Momentum Press: New York, NY, USA, 2015; pp. 84–88. [Google Scholar]
  6. Germer, T.A.; Zwinkels, J.C.; Tsai, B.K. Spectrophotometry: Accurate Measurement of Optical Properties of Materials, 1st ed.; Academic Press: Cambridge, MA, USA, 2014; pp. 75–97. [Google Scholar]
  7. Stenzel, O. The Physics of Thin Film Optical Spectra, 1st ed.; Springer: Berlin/Heidelberg, Germany, 2016; pp. 79–81. [Google Scholar]
  8. Koster, G.; Rijnders, G. In Situ Characterization of Thin Film Growth, 1st ed.; Woodhead Publishing: Sawston, UK, 2011; pp. 64–79. [Google Scholar]
  9. Sun, Z.; To, S.; Zhang, S.; Zhang, G. Theoretical and experimental investigation into non-uniformity of surface generation in micro-milling. Int. J. Mech. Sci. 2018, 140, 313–324. [Google Scholar] [CrossRef]
  10. Wemple, S.H.; DiDomenico, M., Jr. Behavior of the Electronic Dielectric Constant in Covalent and Ionic Materials. Phys. Rev. B Condens. Matter 1971, 3, 1338–1351. [Google Scholar] [CrossRef]
  11. Wemple, S.H. Refractive-index behavior of amorphous semiconductors and glasses. Phys. Rev. B 1973, 7, 3767–3777. [Google Scholar] [CrossRef]
  12. Pan, Y.; Inam, F.; Zhang, M.; Drabold, D.A. Atomistic origin of Urbach tails in amorphous silicon. Phys. Rev. Lett. 2008, 100, 206403. [Google Scholar] [CrossRef]
  13. Zaynobidinov, S.; Ikramov, R.G.; Jalalov, R.M. Urbach energy and the tails of the density of states in amorphous semiconductors. J. Appl. Spectrosc. 2011, 78, 223–227. [Google Scholar] [CrossRef]
  14. Poelman, D.; Smet, P.J. Methods for the determination of the optical constants of thin films from single transmission measurements: A critical review. J. Phys. D 2003, 36, 1850–1857. [Google Scholar] [CrossRef]
  15. Sadek, M.A.H. Principles and Practices of Spectrophotometer; Scholars’ Press: Chisinau, Moldova, 2019; pp. 37–50. [Google Scholar]
  16. Jellison, G.E.; Modine, F.A. Parameterization of the optical functions of amorphous materials in the interband region. Appl. Phys. Lett. 1996, 69, 371–373, Erratum in Appl. Phys. Lett. 1996, 69, 2137. [Google Scholar] [CrossRef]
  17. Blankenhagen, B.; Tonova, D.; Ullmann, J. Application of the Tauc-Lorentz formulation to the interband absorption of optical coating materials. Appl. Opt. 2002, 41, 3137–3141. [Google Scholar] [CrossRef]
  18. Likhachev, D.V. Evaluation of different dispersion models for the thickness-dependent optical properties of thin TiN films in combination with X-ray reflectometry data. Rev. Sci. Instrum. 2019, 90, 023909. [Google Scholar] [CrossRef]
  19. Campi, D.; Coriasso, C. Relationships between optical properties and band parameters in amorphous tetrahedrally bonded materials. Mater. Lett. 1988, 7, 134–137. [Google Scholar]
  20. Franta, D.; Cermák, M.; Vohánka, I.; Ohlídal, I. Dispersion models describing interband electronic transitions combining Tauc’s law and Lorentz model. Thin Solid Films 2017, 631, 12–22. [Google Scholar] [CrossRef]
  21. New Amorphous Dispersion Formula. Available online: https://www.horiba.com/fileadmin/uploads/Scientific/Downloads/OpticalSchool_CN/TN/ellipsometer/New_Amorphous_Dispersion_Formula.pdf (accessed on 17 July 2023).
  22. Dkhilalli, F.; Borchani, S.M.; Rasheed, M.; Barille, R.; Guidara, K.; Megdiche, M. Structural, dielectric, and optical properties of the zinc tungstate ZnWO4 compound. J. Mater. Sci. Mater. Electron. 2018, 29, 6297–6307. [Google Scholar] [CrossRef]
  23. Meneses, D.D.S.; Malki, M.; Echegut, P. Structure and lattice dynamics of binary lead silicate glasses investigated by infrared spectroscopy. J. Non-Cryst. Solids 2006, 352, 769–776. [Google Scholar] [CrossRef]
  24. Fischer, D.; Hertwig, A.; Beck, U.; Negendank, D.; Lohse, V.; Kormunda, M.; Esser, N. Dependence of the optical constants and the performance in the SPREE gas measurement on the thickness of doped tin oxide overcoatings. Appl. Surf. Sci. 2017, 421B, 480–486. [Google Scholar] [CrossRef]
  25. Li, Y.; Kocherga, M.; Park, S.; Lata, M.; McLamb, M.; Boreman, G.; Schmedake, T.A.; Hofmann, T. Optical dielectric function of Si(2,6-bis(benzimidazol-2’-yl)pyridine)2 determined by spectroscopic ellipsometry. Opt. Mater. Express 2019, 9, 3469. [Google Scholar] [CrossRef]
  26. Lorentz Dispersion Model. Available online: https://pdf4pro.com/amp/view/lorentz-dispersion-model-horiba-20da3.html (accessed on 21 June 2023).
  27. Foldyna, M.; Postava, K.; Bouchal, J.; Pitora, J.; Yamaguchi, T. Model dielectric function of amorphous materials including Urbach tail. Proc. SPIE 2003, 5445, 301–305. [Google Scholar]
  28. Márquez, E.; Ruíz-Pérez, J.J.; Ballester, M.; Márquez, A.P.; Blanco, E.; Minkov, D.; Fernandez, S.M.; Saugar, E. Optical Characterization of H-Free a-Si Layers Grown by rf-Magnetron Sputtering by Inverse Synthesis Using Matlab: Tauc–Lorentz–Urbach Parameterization. Coatings 2021, 11, 1324. [Google Scholar] [CrossRef]
  29. Rodruguez-de Marcos, L.; Larruquert, J.I. Analytic optical-constant model derived from Tauc-Lorentz and Urbach tail. Opt. Express 2016, 24, 28561–28572. [Google Scholar] [CrossRef]
  30. Ballester, M.; García, M.; Márquez, A.P.; Blanco, E.; Fernández, S.M.; Minkov, D.; Katsaggelos, A.K.; Cossairt, O.; Willomitzer, F.; Marquez, E. Application of the Holomorphic Tauc-Lorentz-Urbach Function to Extract the Optical Constants of Amorphous Semiconductor Thin Films. Coatings 2022, 12, 1549. [Google Scholar] [CrossRef]
  31. Franta, D.; Necas, D.; Ohlidal, I.; Giglia, A. Dispersion model for optical thin films applicable in wide spectral range. Proc. SPIE 2015, 9628, 342–353. [Google Scholar]
  32. Franta, D.; Necas, D.; Ohlidal, I.; Giglia, A. Optical characterization of SiO2 thin films using universal dispersion model over wide spectral range. Proc. SPIE 2016, 9628, 253–267. [Google Scholar]
  33. Fujiwara, H. Spectroscopic Ellipsometry: Principles and Applications, 1st ed.; Wiley: New York, NY, USA, 2007; pp. 67–96. [Google Scholar]
  34. Nečas, D.; Ohlídal, I.; Franta, D. Variable-angle spectroscopic ellipsometry of considerably non-uniform thin films. J. Opt. 2011, 13, 085705. [Google Scholar] [CrossRef]
  35. Hilfiker, J.N.; Hale, J.S.; Herzinger, C.M.; Tiwald, T.; Hong, N.; Schoche, S.; Arwin, H. Estimating depolarization with the Jones matrix quality factor. Appl. Surf. Sci. 2017, 421, 494–499. [Google Scholar] [CrossRef]
  36. Johnson, B.I.; Hodges, G.; Cushman, C.; Banerjee, J.; Smith, N.J.; Linford, M.R. Appropriate Backside Roughening is Key for Good Spectroscopic Ellipsometry Analysis of Transparent Materials. Vac. Technol. Coat. 2018, 19, 26–33. [Google Scholar]
  37. Synowicky, R.A. Suppression of backside reflections from transparent substrates. Phys. Stat. Sol. 2008, 5, 1085–1088. [Google Scholar] [CrossRef]
  38. Patel, D.I.; Jacobsen, C.; Hilfiker, J.N.; Linford, M.R. Fitting the Spectroscopic Ellipsometry Data from a Rather Thick (Organic?) Film on Fused Silica Part 1-Using the ‘Global Fit’. Vac. Technol. Coat. 2021, 22, 34–38. [Google Scholar]
  39. Hilfiker, J.N.; Linford, M.R. Fitting the Spectroscopic Ellipsometry Data from a Rather Thick (Organic?) Film on Fused Silica Part 2-Thickness non-uniformity, Absorption and Roughness. Vac. Technol. Coat. 2021, 22, 34–39. [Google Scholar]
  40. Kaflé, B.P. Chemical Analysis and Material Characterization by Spectrophotometry, 1st ed.; Kindle: Seattle, WA, USA, 2019; pp. 72–75. [Google Scholar]
  41. LAMBDA 365 UV/Vis Spectrophotometer. Available online: https://www.perkinelmer.com/product/lambda-365-spectrophotometer-uv-express-n4100020 (accessed on 23 May 2023).
  42. Tien, P.K. Light Waves in Thin Films and Integrated Optics. Appl. Opt. 1971, 10, 2395–2413. [Google Scholar] [CrossRef]
  43. Solid Sample Reflectance Measurements. Available online: https://www.shimadzu.com/an/service-support/technical-support/analysis-basics/fundamentals-uv/reflectance_measurements.html (accessed on 14 July 2023).
  44. Noble, S.D.; Crowe, T.G. Sample holder and methodology for measuring the reflectance and transmittance of narrow-leaf samples. Appl. Opt. 2007, 46, 2395–2413. [Google Scholar] [CrossRef]
  45. Minkov, D. DSc Thesis: Characterization of Thin Films and Surface Cracks by Electromagnetic Methods and Technologies; Technical University: Sofia, Bulgaria, 2018; pp. 114–119. [Google Scholar]
  46. Mostafa, N.Y.; Badawi, A.; Ahmed, S.I. Influence of Cu and Ag doping on structure and optical properties of In2O3 thin film prepared by spray pyrolysis. Results Phys. 2018, 10, 126–131. [Google Scholar] [CrossRef]
  47. Kocyigit, A.; Erdal, M.O.; Yıldırım, M. Effect of Indium Doping on Optical Parameter Properties of Sol–Gel-Derived ZnO Thin Films. Z. Naturforsch. A 2019, 74, 915–923. [Google Scholar] [CrossRef]
  48. Alsaad, A.M.; Al-Dairy, A.R.; Ahmad, A.A.; Al-Anbar, A.S.; Al_Bataineh, Q.M. Synthesis and characterization of as-grown doped polymerized (PMMA-PVA)/ZnO NPs hybrid thin films. Polym. Bull. 2021, 79, 2019–2040. [Google Scholar] [CrossRef]
  49. Althubiti, N.A.; Atta, A.; Al-Harbi, N.; Sendi, R.K.; Abdelhamied, M.M. Structural, characterization and linear/nonlinear optical properties of oxygen beam irradiated PEO/NiO composite films. Opt. Quantum Electron. 2023, 55, 348. [Google Scholar] [CrossRef]
  50. Thin-Film Interference. Available online: https://en.wikipedia.org/wiki/Thin-film_interference/ (accessed on 16 June 2021).
  51. Thin-Film Interference. Available online: https://www.britannica.com/science/light/Thin-film-interference/ (accessed on 16 June 2021).
  52. Minkov, D.; Vateva, E.; Skordeva, E.; Arsova, D.; Nikiforova, M. Optical properties of Ge-As-S thin films. J. Non-Cryst. Solids 1987, 90, 481–484. [Google Scholar] [CrossRef]
  53. Wagner, T.; Marquez, E.; Fernandez-Pena, J.; Gonzalez-Leal, J.M.; Ewen, P.J.S.; Kasap, S.O. The kinetics of the photo-induced solid-state chemical reaction in Ag/As33S67 bilayers and its reaction products. Philos. Mag. B 1999, 79, 223–237. [Google Scholar] [CrossRef]
  54. Swanepoel, R. Determination of the thickness and optical constants of amorphous silicon. J. Phys. E 1983, 16, 1214–1222. [Google Scholar] [CrossRef]
  55. Swanepoel, R. Determination of surface roughness and optical constants of inhomogeneous amorphous silicon films. J. Phys. E Sci. Instrum. 1984, 17, 896–903. [Google Scholar] [CrossRef]
  56. Leal, J.M.G.; Alcon, R.P.; Angel, J.A.; Minkov, D.A.; Marquez, E. Influence of substrate absorption on the optical and geometrical characterization of thin dielectric films. Appl. Opt. 2002, 41, 7300–7308. [Google Scholar] [CrossRef]
  57. Minkov, D.A.; Gavrilov, G.M.; Angelov, G.V.; Moreno, G.M.D.; Vazquez, C.G.; Ruano, S.M.F.; Marquez, E. Optimisation of the envelope method for characterisation of optical thin film on substrate specimens from their normal incidence transmittance spectrum. Thin Solid Films 2018, 645, 370–378. [Google Scholar] [CrossRef]
  58. Minkov, D.A.; Angelov, G.V.; Nestorov, R.N.; Marquez, E. Perfecting the dispersion model free characterization of a thin film on a substrate specimen from its normal incidence interference transmittance spectrum. Thin Solid Films 2020, 706, 137984. [Google Scholar] [CrossRef]
  59. Minkov, D.; Marquez, E.; Angelov, G.; Gavrilov, G.; Ruano, S.M.F.; Saugar, E. Further increasing the accuracy of characterization of a thin dielectric or semiconductor film on a substrate from its interference transmittance spectrum. Materials 2021, 14, 4681. [Google Scholar] [CrossRef]
  60. Márquez, E.; Saugar, E.; Díaz, J.M.; García-Vázquez, C.; Fernández-Ruano, S.M.; Blanco, E.; Ruiz-Pérez, J.J.; Minkov, D.A. The influence of Ar pressure on the structure and optical properties of non-hydrogenated a-Si thin films grown by rf magnetron sputtering onto room-temperature glass substrates. J. Non-Cryst. Solids 2019, 517, 32–43. [Google Scholar] [CrossRef]
  61. Minkov, D.A.; Angelov, G.; Nestorov, R.; Nezhdanov, A.; Usanov, D.; Kudryashov, M.; Mashin, A. Optical characterization of AsxTe100-x Films Grown by Plasma Deposition Based on the Advanced Optimizing Envelope Method. Materials 2020, 13, 2981. [Google Scholar] [CrossRef]
  62. Minkov, D.A.; Angelov, G.V.; Nestorov, R.N.; Marquez, E.; Blanco, E.; Ruiz-Perez, J.J. Comparative study of the accuracy of characterization of thin films a-Si on glass substrates from their interference normal incidence transmittance spectrum by the Tauc-Lorentz-Urbach, the Cody-Lorentz-Urbach, the optimized envelopes and the optimized graphical methods. Mater. Res. Express 2019, 6, 036410. [Google Scholar]
  63. Available online: https://qd-uki.co.uk/ellipsometers/j-a-woollam-wvase-ellipsometer-software/ (accessed on 3 May 2023).
  64. Available online: https://www.horiba.com/en_en/products/detail/action/show/Product/deltapsi2-software-1648/ (accessed on 23 August 2022).
  65. Available online: https://www.perkinelmer.com/product/kit-uvwinlab-v7-3-std-software-l6100127 (accessed on 19 April 2023).
  66. Available online: https://www.thermofisher.com/bg/en/home/industrial/spectroscopy-elemental-isotope-analysis/molecular-spectroscopy/uv-vis-spectrophotometry/software.html (accessed on 26 January 2023).
  67. Minkov, D.A. Calculation of the optical constants of a thin layer upon a transparent substrate from the reflection spectrum. J. Phys. D 1989, 22, 1157–1161. [Google Scholar] [CrossRef]
  68. Gonzalez-Leal, J.M.; Marquez, E.; Bernal-Oliva, A.M.; Ruiz-Perez, J.J.; Jimenez-Garay, R. Derivation of the optical constants of thermally-evaporated uniform films of binary chalcogenide glasses using only their reflection spectra. Thin Solid Films 1998, 317, 223–227. [Google Scholar] [CrossRef]
  69. Gungor, T.; Saka, B. Calculation of the optical constants of a thin layer upon a transparent substrate from the reflection spectrum using a genetic algorithm. Thin Solid Films 2004, 467, 319–325. [Google Scholar] [CrossRef]
  70. Tabassum, S.; Dong, L.; Kumar, R. Determination of dynamic variations in the optical properties of graphene oxide in response to gas exposure based on thin-film interference. Opt. Express 2018, 26, 6331–6344. [Google Scholar] [CrossRef]
  71. Hassanien, A.S.; Aly, K.A.; Elsaeedy, H.I.; Alqahtani, A. Optical characterization and dispersion discussions of the novel thermally evaporated thin a-S50-xGe10CdxTe40 films. Appl. Phys. A 2022, 128, 1021. [Google Scholar] [CrossRef]
  72. Dangi, R.; Basnet, B.; Pandey, M.; Bhusal, S.; Budhathoki, B.; Parajuli, K.; Tiwari, S.K.; Kafle, B.P. Effect of Oxygen Vacancy on the Crystallinity and Optical Band Gap in Tin Oxide Thin Film. Energies 2023, 16, 2653. [Google Scholar] [CrossRef]
  73. Rusli, N.; Amaratunga, G.A.J. Determination of the optical constants and thickness of thin films on slightly absorbing substrates. Appl. Opt. 1995, 34, 7914–7924. [Google Scholar] [CrossRef]
  74. Ruiz-Pérez, J.; Márquez, E.; Minkov, D.; Reyes, J.; Ramírez-Malo, J.B.; Villares, P.; Jiménez-Garay, R. Computation of the optical constants of thermally-evaporated thin films of GeSe2 chalcogenide glass from their reflection spectra. Phys. Scr. 1996, 53, 76–82. [Google Scholar] [CrossRef]
  75. Márquez, E.; González-Leal, J.M.; Prieto-Alcón, R.; Vlcek, M.; Stronski, A.; Wagner, T.; Minkov, D. Optical characterization of thermally evaporated thin films of As40S40Se20 chalcogenide glass by reflectance measurements. Appl. Phys. A 1998, 67, 371–378. [Google Scholar]
  76. Dahshan, A.; Sharma, P.; Aly, K.A. Optical constants of Ge-Sb-Se-I chalco-halide glasses using a single reflectance spectrum. Infrared Phys. Technol. 2019, 102, 102997. [Google Scholar] [CrossRef]
  77. Boukhris, I.; Kebaili, I.; Znaidia, S.; Neffati, R.; Hegazy, H.H.; Aly, K.A.; Mehta, N.; Dahshan, A. Optical constants of Sn-doped amorphous Ge-As-Te thin films and their physical characterization. Phys. B Condens. 2020, 583, 412066. [Google Scholar] [CrossRef]
  78. Ruiz-Perez, J.J.; Gonzalez-Leal, J.M.; Minkov, D.; Marquez, E. Method for determining the optical constants of thin dielectric films with variable thickness using only their shrunk reflection spectra. J. Phys. D 2001, 34, 2489–2496. [Google Scholar] [CrossRef]
  79. Pisarkiewicz, T. Reflection spectrum for a thin film with non-uniform thickness. J. Phys. D. 1994, 27, 160–164. [Google Scholar] [CrossRef]
  80. Necas, D.; Ohlıdal, I.; Franta, D. The reflectance of non-uniform thin films. J. Opt. A Pure Appl. Opt. 2009, 11, 045202. [Google Scholar] [CrossRef]
  81. Minkov, D. Flow-graph approach for optical analysis of planar structures. Appl. Opt. 1994, 33, 7698–7703. [Google Scholar] [CrossRef]
  82. Gonzalez-Leal, J.M.; Prieto-Alcon, R.; Stuchlik, M.; Vlcek, M.; Elliott, S.R.; Marquez, E. Determination of the surface roughness and refractive index of amorphous As40S60 films deposited by spin coating. Opt. Mater. 2004, 27, 147–154. [Google Scholar] [CrossRef]
  83. Minkov, D.; Swanepoel, R. Computerization of the optical characterization of a thin dielectric film. Opt. Eng. 1993, 32, 3333–3337. [Google Scholar] [CrossRef]
  84. Levenhuk Zoom & Joy. Available online: https://www.levenhuk.com/catalogue/accessories/levenhuk-n18-ng-prepared-slides-set/#.XsbUvmgzaUk (accessed on 6 June 2020).
  85. Tripathy, S.K. Refractive indices of semiconductors from energy gaps. Opt. Mater. 2015, 46, 240–246. [Google Scholar] [CrossRef]
  86. Minkov, D.; Nedelchev, L.; Angelov, G.; Marquez, E.; Blagoeva, B.; Mateev, G.; Nazarova, D. Hybrid Dispersion Model Characterization of PAZO Azopolymer Thin Films over the Entire Transmittance Spectrum Measured in the UV/VIS/NIR Spectral Region. Materials 2022, 15, 8617. [Google Scholar] [CrossRef]
  87. Principles of UV/Vis Spectroscopy (7) Bandwidth. Available online: https://www.jasco-global.com/principle/principles-of-uv-vis-spectroscopy-7-bandwidth/ (accessed on 24 April 2023).
  88. Gavrilov, G.M.; Minkov, D.A.; Marquez, E.; Ruano, S.M.F. Advanced computer drawing envelopes of transmittance spectra of thin film specimens. Int. J. Adv. Res. Sci. Eng. Technol. 2016, 3, 163–168. [Google Scholar]
  89. Minkov, D.A.; Gavrilov, G.M.; Marquez, E.; Ruano, S.M.F.; Stoynova, A.V. Development of algorithm for computer drawing envelopes of interference reflectance spectra for thin film specimens. Optik 2017, 132, 320–328. [Google Scholar] [CrossRef]
  90. Nestorov, R.N. Selection of error metric for accurate characterization of a thin dielectric or semiconductor film on glass substrate by the optimizing envelope method. Int. J. Adv. Res. Sci. Eng. Technol. 2020, 7, 1–11. [Google Scholar]
  91. Fixed Angle Reflectance Accessory—Fixed 60. Available online: https://www.perkinelmer.com/product/relative-specular-reflectance-accessory-b0086703 (accessed on 18 April 2023).
  92. Relative Specular Reflectance Accessory, 6 Degrees for PerkinElmer Lambda 800/900-B0086703. Available online: https://uvison.com/chromatography-supplies/perkinelmer/perkinelmer-accessories/perkinelmer-relative-specular-reflectance-accessory-6-degrees-for-perkinelmer-lambda-800-900-b0086703 (accessed on 17 May 2023).
  93. Bruggeman, D.A.G. Calculation of various physical constants of heterogeneous substances and dielectric constants and conductivities of mixed bodies of isotropic substances. Ann. Phys. 1935, 416, 636–664. [Google Scholar] [CrossRef]
  94. Guerrero, E.; Strubbe, D.A. Computational generation of voids in a-Si and a-Si:H by cavitation at low density. Phys. Rev. Mater. 2020, 4, 025601. [Google Scholar] [CrossRef]
  95. Chen, H.; Shen, W.Z. Perspectives in the characteristics and applications of Tauc-Lorentz dielectric function model. Eur. Phys. J. B 2005, 43, 503–507. [Google Scholar] [CrossRef]
  96. Dodd, M.J.A.; Polman, A.; Zijlstra, T.; Drift, E.W.J.M. Amorphous silicon waveguides for microphotonics. J. Appl. Phys. 2002, 92, 649–653. [Google Scholar] [CrossRef]
  97. Brent, R.P. Algorithms for Minimization without Derivatives, 1st ed.; Dover Publications: Mineola, NY, USA, 2013; pp. 21–48. [Google Scholar]
  98. Kochenderfer, M.J.; Wheeler, T.A. Algorithms for Optimization, 1st ed.; MIT Press: Cambridge, MA, USA, 2019; pp. 103–154. [Google Scholar]
  99. Floris, F.; Fornasari, L.; Marini, A.; Bellani, V.; Banfi, F.; Roddaro, S.; Ercolani, D.; Rocci, M.; Beltram, F.; Cecchini, M.; et al. Self-Assembled InAs Nanowires as Optical Reflectors. Nanomaterials 2017, 7, 7110400. [Google Scholar] [CrossRef] [PubMed]
  100. Floris, F.; Fornasari, L.; Bellani, V.; Marini, A.; Banfi, F.; Marabelli, F.; Beltram, F.; Ercolani, D.; Battiato, S.; Sorba, L.; et al. Strong Modulations of Optical Reflectance in Tapered Core–Shell Nanowires. Materials 2019, 12, 12213572. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Calculated UV/Vis/NIR spectra for single oscillator DMs 1TL, 1CC, 1NA, and 1GA with identical central energy E0 = 3 eV, broadening term B0 = 2 eV, ε ˙ (0) = ε ˙ (E = 0) = n2(0) = 5 and ε ˙ (∞) = ε ˙ (E→∞) = n2(∞) = 1, as these values are selected to represent a semiconductor or dielectric material. (a) εr(E) and εi(E); (b) εr(λ) and εi(λ); (c) normal incidence R(λ) for the boundary air/material.
Figure 1. Calculated UV/Vis/NIR spectra for single oscillator DMs 1TL, 1CC, 1NA, and 1GA with identical central energy E0 = 3 eV, broadening term B0 = 2 eV, ε ˙ (0) = ε ˙ (E = 0) = n2(0) = 5 and ε ˙ (∞) = ε ˙ (E→∞) = n2(∞) = 1, as these values are selected to represent a semiconductor or dielectric material. (a) εr(E) and εi(E); (b) εr(λ) and εi(λ); (c) normal incidence R(λ) for the boundary air/material.
Nanomaterials 13 02407 g001
Figure 2. A sketch of reflectance R(λ) of a sample consisting of a non-uniform thin film on a thick planar substrate. The unknown optical characteristics of the film are printed in red color. The angle Φ is several degrees for quasi-normal incidence of light.
Figure 2. A sketch of reflectance R(λ) of a sample consisting of a non-uniform thin film on a thick planar substrate. The unknown optical characteristics of the film are printed in red color. The angle Φ is several degrees for quasi-normal incidence of light.
Nanomaterials 13 02407 g002
Figure 3. Simulated interference reflectance and transmittance spectra, their envelopes, and their features, for quasi-normal incidence of light to model specimens corresponding to a thin film a-Si:H on a thick planar glass substrate. (a) For a uniform film on a finite transparent substrate; (b) For two specimens of non-uniform film on finite transparent and non-transparent substrate, respectively; (c) For the non-uniform film on finite non-transparent substrate from (b); (d) For two specimens of much thinner non-uniform film on non-transparent finite and semi-infinite substrate, respectively. The region of significant absorbance in the substrate is exhibited by a dashed red colored rectangle in (b,d).
Figure 3. Simulated interference reflectance and transmittance spectra, their envelopes, and their features, for quasi-normal incidence of light to model specimens corresponding to a thin film a-Si:H on a thick planar glass substrate. (a) For a uniform film on a finite transparent substrate; (b) For two specimens of non-uniform film on finite transparent and non-transparent substrate, respectively; (c) For the non-uniform film on finite non-transparent substrate from (b); (d) For two specimens of much thinner non-uniform film on non-transparent finite and semi-infinite substrate, respectively. The region of significant absorbance in the substrate is exhibited by a dashed red colored rectangle in (b,d).
Nanomaterials 13 02407 g003
Figure 4. The algorithm of OEMR, exhibiting its detailed parts A and C for characterization of a thin semiconductor or dielectric film on a finite and generally non-transparent thick planar substrate.
Figure 4. The algorithm of OEMR, exhibiting its detailed parts A and C for characterization of a thin semiconductor or dielectric film on a finite and generally non-transparent thick planar substrate.
Nanomaterials 13 02407 g004
Figure 5. The optimizing part B of the algorithm of OEMR.
Figure 5. The optimizing part B of the algorithm of OEMR.
Nanomaterials 13 02407 g005
Figure 6. Reflectance spectra of the samples SP1 and SP2, their smoothing, and external envelopes. (a,b) Magnified images of the smoothing around some maxima and minima of R(λ) of SP2, respectively; (c) R(λ) of SP1, its external envelopes and substrate characteristics; (d) R(λ) of SP2, its external envelopes and substrate characteristics.
Figure 6. Reflectance spectra of the samples SP1 and SP2, their smoothing, and external envelopes. (a,b) Magnified images of the smoothing around some maxima and minima of R(λ) of SP2, respectively; (c) R(λ) of SP1, its external envelopes and substrate characteristics; (d) R(λ) of SP2, its external envelopes and substrate characteristics.
Nanomaterials 13 02407 g006
Figure 7. Results regarding the computation of nc(λ) and kc(λ) of both a-Si films via execution of part C from the OEMR algorithm. (a,c) Wemple-DiDomenico plots; (b,d) kc(λ) computed from one of the envelopes R+(λ), R(λ), or their half-sum.
Figure 7. Results regarding the computation of nc(λ) and kc(λ) of both a-Si films via execution of part C from the OEMR algorithm. (a,c) Wemple-DiDomenico plots; (b,d) kc(λ) computed from one of the envelopes R+(λ), R(λ), or their half-sum.
Nanomaterials 13 02407 g007
Figure 8. Differences Rc(λ) − R(λ) between the computed and the measured reflectance spectra for several DMs (a,b) for the single oscillator DMs 1TL, 1NA, TLUR, and TLUF; (c,d) for the DMs featured in Table 4.
Figure 8. Differences Rc(λ) − R(λ) between the computed and the measured reflectance spectra for several DMs (a,b) for the single oscillator DMs 1TL, 1NA, TLUR, and TLUF; (c,d) for the DMs featured in Table 4.
Nanomaterials 13 02407 g008
Figure 9. Spectral dependencies of nc and log10αc of the a-Si films SP1 and SP2, computed over the entire measured UV/Vis/NIR spectrum using each of the four DMs represented in Table 4. (a,b) For SP1; (c,d) For SP2. The most accurate nc(λ) and log10αc(E) are displayed by the red colored lines.
Figure 9. Spectral dependencies of nc and log10αc of the a-Si films SP1 and SP2, computed over the entire measured UV/Vis/NIR spectrum using each of the four DMs represented in Table 4. (a,b) For SP1; (c,d) For SP2. The most accurate nc(λ) and log10αc(E) are displayed by the red colored lines.
Nanomaterials 13 02407 g009
Table 1. Computed data regarding the thickness parameters of amorphous thin films characterized via EMT or OEMT. The text and data in red color represent the most accurate results for these parameters.
Table 1. Computed data regarding the thickness parameters of amorphous thin films characterized via EMT or OEMT. The text and data in red color represent the most accurate results for these parameters.
Attributes of the Sample;
and the Characterization
Film Material, Sample, [Reference] d ¯ c (nm), Δdc (nm),
[l1c, l2c]
RE( d ¯ c) (%)
uniform film, transparent substrate; EMT, common smoothing of T(λ), [l1c, l2c] is not optimizeda-Si, A029, [60]1173, 0, [1, 12]0.215
a-Si, A074, [60]1269, 0, [2, 14]0.213
uniform film, absorbing substrate; EMT, common smoothing of T(λ), [l1c, l2c] is not optimizeda-As40S60, [56]1409, 0, [11, 20]0.261
non-uniform film, absorbing substrate; EMT, common smoothing of T(λ), Δdc and [l1c, l2c] are not optimizeda-As40S60, [56]2689, 50, [3, 24]0.496
non-uniform film, absorbing substrate; OEMT, external smoothing of T(λ), l1c = 1, Δdc and l2c are optimizeda-Si, A041, [58]3929.9, 53.5, [1, 17]0.245
a-As98Te2, [59]1983.8, 22.7, [1, 12]0.133
non-uniform film, absorbing substrate; OEMT, external smoothing of T(λ), Δdc and [l1c, l2c] are optimizeda-Si, A041, [59]3949.2, 53.0, [5, 14]0.090
a-As98Te2, [61]1983.2, 23.9, [3, 9]0.043
Table 2. Results for the two a-Si films, computed by parts A and B from the algorithm of OEMR. The data corresponding to smallest RE( d ¯ c), for each of these films, are printed in red color, as their respective optimized thickness parameters d ¯ c and ∆dc are in bold red.
Table 2. Results for the two a-Si films, computed by parts A and B from the algorithm of OEMR. The data corresponding to smallest RE( d ¯ c), for each of these films, are printed in red color, as their respective optimized thickness parameters d ¯ c and ∆dc are in bold red.
FILM SP1
ERMERMc[l1c, l2c]m1cdc (nm) d ¯ c (nm)RE( d ¯ c) (%)
DECONVOLUTION ENVELOPES OF R(λ)
ERM16.61 × 10−4[2, 7]210.7672.5
ERM20.1263 nm[2, 7]210.7672.50.1127
INTERNAL ENVELOPES OF R(λ)
ERM15.63 × 10−4[2, 7]211.7674.5
ERM20.1144 nm[2, 7]211.8674.70.1017
MEDIUM ENVELOPES OF R(λ)
ERM15.64 × 10−4[2, 7]211.6674.5
ERM20.1127 nm[2, 7]211.7674.50.1003
EXTERNAL ENVELOPES OF R(λ)
ERM14.62 × 10−4[2, 7]211.4674.2
ERM20.1115 nm[2, 7]211.5674.40.0992
Optimized thickness parameters11.5674.3
FILM SP2
ERMERMc[l1c, l2c]m1cdc (nm) d ¯ c (nm)RE( d ¯ c) (%)
DECONVOLUTION ENVELOPES OF R(λ)
ERM11.72 × 10−3[7, 44]1042.53589.9
ERM20.2540 nm[5, 44]1139.83701.60.2745
INTERNAL ENVELOPES OF R(λ)
ERM11.20 × 10−3[8, 40]1140.73713.0
ERM20.2100 nm[5, 41]1140.53704.10.2098
MEDIUM ENVELOPES OF R(λ)
ERM11.15 × 10−3[4, 40]1237.53837.7
ERM20.2133 nm[4, 40]1237.43835.50.2224
EXTERNAL ENVELOPES OF R(λ)
ERM11.02 × 10−3[1, 39]1237.03848.7
ERM20.1782 nm[3, 39]1236.93845.40.1715
Optimized thickness parameters37.03847.1
Table 3. Computed parameters from utilization of four broad spectrum single electron oscillator DMs, employed independently from each other for characterizations of the a-Si films SP1 and SP2. The results regarding the most accurate amongst the respective four characterizations are printed in blue color, for each of these films.
Table 3. Computed parameters from utilization of four broad spectrum single electron oscillator DMs, employed independently from each other for characterizations of the a-Si films SP1 and SP2. The results regarding the most accurate amongst the respective four characterizations are printed in blue color, for each of these films.
FILM SP1
1TL
ATL1 (eV)ETL1 (eV)BTL1 (eV)Eg (eV)εr (∞) d ¯ c (nm)Δdc (nm)FOM
91.253.2811.0641.2562.152742.830.619.16
1NA
fNA1ENA1 (eV)BNA1 (eV)Eg (eV)nc (∞) d ¯ c (nm)Δdc (nm)FOM
0.12223.5430.60260.0002.604712.225.719.30
TLUR
ATL1 (eV)ETL1 (eV)BTL1 (eV)Eg (eV)Eb (eV) d ¯ c (nm)Δdc (nm)FOM
134.03.3131.2691.4502.236729.927.914.65
TLUF
ATL1 (eV)ETL1 (eV)BTL1 (eV)Eg (eV)Eb (eV)εr () d ¯ c (nm)Δdc (nm)FOM
130.23.3261.3041.3942.2111.000718.324.813.66
FILM SP2
1TL
ATL1 (eV)ETL1 (eV)BTL1 (eV)Eg (eV)εr (∞) d ¯ c (nm)Δdc (nm)FOM
46.612.9900.93310.66845.9923846.040.017.35
1NA
fNA1ENA1 (eV)BNA1 (eV)Eg (eV)nc (∞) d ¯ c (nm)Δdc (nm)FOM
0.11513.2240.58850.0003.1973700.535.915.78
TLUR
ATL1 (eV)ETL1 (eV)BTL1 (eV)Eg (eV)Eb (eV) d ¯ c (nm)Δdc (nm)FOM
88.812.6942.1770.69410.83183571.739.751.1
TLUF
ATL1 (eV)ETL1 (eV)BTL1 (eV)Eg (eV)Eb (eV)εr (∞) d ¯ c (nm)Δdc (nm)FOM
96.783.3821.5500.89531.3582.8143701.933.57.65
Table 4. Computed parameters from characterizations of the a-Si films based on four broad spectrum oscillator DMs, each of them including three oscillators from the same type in the spectral region of strong absorption in the film. The data regarding the most accurate characterizations of the films by these DMs are in red color.
Table 4. Computed parameters from characterizations of the a-Si films based on four broad spectrum oscillator DMs, each of them including three oscillators from the same type in the spectral region of strong absorption in the film. The data regarding the most accurate characterizations of the films by these DMs are in red color.
FILM SP1
3TL with 1PE and 1GA
ATL1 (eV)ETL1 (eV)BTL1 (eV)ATL2 (eV)ETL2 (eV)BTL2 (eV)ATL3 (eV)ETL3 (eV)BTL3 (eV)Eg (eV)
10.63 4.1950.507257.523.6730.956615.512.9920.90860.7942
APE1 (eV2)EPE1 (eV)fGA1EGA1 (eV)BGA1 (eV) d ¯ c (nm)Δdc (nm)FOM
0.07330.34930.36190.70700.4239664.70.06646.80
3UD with UT and 1GA
AUD1EUD1 (eV)BUD1 (eV)AUD2EUD2 (eV)BUD2 (eV)AUD3EUD3 (eV)BUD3 (eV)Eg (eV)
8.4075.0240.00429092.223.6000.567918.062.9190.51701.209
NvcEh (eV)fUTEU (eV)fGA1EGA1 (eV)BGA1 (eV) d ¯ c (nm)Δdc (nm)FOM
286.433.110.04288 0.79720.000550.000420.3111670.45.06.64
OEMR with 3TL
ATL1 (eV)ETL1 (eV)BTL1 (eV)ATL2 (eV)ETL2 (eV)BTL2 (eV)ATL3 (eV)ETL3 (eV)BTL3 (eV)Eg (eV)
65.743.7461.19912.733.2840.70218.3992.7900.78750.8749
Ew (eV)Ec (eV)FOM
1.6631.5206.55
OEMR with 3NA
fNA1ENA1 (eV)BNA1 (eV)fNA2ENA2 (eV)BNA2 (eV)fNA3ENA3 (eV)BNA3 (eV)Eg (eV)
1.1474.2681.3090.43971.6860.48830.011641.5450.037181.428
Ew (eV)Ec (eV)FOM
1.7181.5285.46
FILM SP2
3TL with 1PE and 1GA
ATL1 (eV)ETL1 (eV)BTL1 (eV)ATL2 (eV)ETL2 (eV)BTL2 (eV)ATL3 (eV)ETL3 (eV)BTL3 (eV)Eg (eV)
99.563.2361.15832.912.5440.979948.171.6830.76901.424
APE1 (eV2)EPE1 (eV)fGA1EGA1 (eV)BGA1 (eV) d ¯ c (nm)Δdc (nm)FOM
272.010.091.0381.6430.42713758.335.95.49
3UD with UT and 1GA
AUD1EUD1 (eV)BUD1 (eV)AUD2EUD2 (eV)BUD2 (eV)AUD3EUD3 (eV)BUD3 (eV)Eg (eV)
14.563.2710.63238.8821.9911.3880.27741.58480.0013211.593
NvcEh (eV)fUTEU (eV)fGA1EGA1 (eV)BGA1 (eV) d ¯ c (nm)Δdc (nm)FOM
345.741.700.33010.16620.0007161.29890.3746923741.535.25.84
OEMR with 3TL
ATL1 (eV)ETL1 (eV)BTL1 (eV)ATL2 (eV)ETL2 (eV)BTL2 (eV)ATL3 (eV)ETL3 (eV)BTL3 (eV)Eg (eV)
5.5793.3000.63850.24992.8370.39548.6232.6800.41940.00022
Ew (eV)Ec (eV)FOM
2.3381.21916.22
OEMR with 3NA
fNA1ENA1 (eV)BNA1 (eV)fNA2ENA2 (eV)BNA2 (eV)fNA3ENA3 (eV)BNA3 (eV)Eg (eV)
0.70394.9010.93770.017621.5350.11400.86741.2420.33891.232
Ew (eV)Ec (eV)FOM
1.3201.3055.18
Table 5. FOMs and thickness parameters already computed using selected characterizations of the a-Si films, as well as relative errors in the computations of these parameters with respect to their values computed by OEMR.
Table 5. FOMs and thickness parameters already computed using selected characterizations of the a-Si films, as well as relative errors in the computations of these parameters with respect to their values computed by OEMR.
FILM SP1
Characterization methodFOM d ¯ c (nm)REc( d ¯ c) (%)dc (nm)REc(∆dc) (%)
OEMR with 3NA5.46674.3-11.5-
TLUF13.66718.36.524.8115.7
3UD with UT and 1GA6.64670.4−0.65.0−56.5
FILM SP2
Characterization methodFOM d ¯ c (nm)RE( d ¯ c) (%)dc (nm)REc(∆dc) (%)
OEMR with 3NA5.183847.1-37.0-
TLUF7.653701.9−3.833.5−9.5
3TL with 1PE and 1GA5.493758.3−2.335.9−3.0
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.

Share and Cite

MDPI and ACS Style

Minkov, D.; Angelov, G.; Marquez, E.; Radonov, R.; Rusev, R.; Nikolov, D.; Ruano, S. Increasing the Accuracy of the Characterization of a Thin Semiconductor or Dielectric Film on a Substrate from Only One Quasi-Normal Incidence UV/Vis/NIR Reflectance Spectrum of the Sample. Nanomaterials 2023, 13, 2407. https://doi.org/10.3390/nano13172407

AMA Style

Minkov D, Angelov G, Marquez E, Radonov R, Rusev R, Nikolov D, Ruano S. Increasing the Accuracy of the Characterization of a Thin Semiconductor or Dielectric Film on a Substrate from Only One Quasi-Normal Incidence UV/Vis/NIR Reflectance Spectrum of the Sample. Nanomaterials. 2023; 13(17):2407. https://doi.org/10.3390/nano13172407

Chicago/Turabian Style

Minkov, Dorian, George Angelov, Emilio Marquez, Rossen Radonov, Rostislav Rusev, Dimitar Nikolov, and Susana Ruano. 2023. "Increasing the Accuracy of the Characterization of a Thin Semiconductor or Dielectric Film on a Substrate from Only One Quasi-Normal Incidence UV/Vis/NIR Reflectance Spectrum of the Sample" Nanomaterials 13, no. 17: 2407. https://doi.org/10.3390/nano13172407

APA Style

Minkov, D., Angelov, G., Marquez, E., Radonov, R., Rusev, R., Nikolov, D., & Ruano, S. (2023). Increasing the Accuracy of the Characterization of a Thin Semiconductor or Dielectric Film on a Substrate from Only One Quasi-Normal Incidence UV/Vis/NIR Reflectance Spectrum of the Sample. Nanomaterials, 13(17), 2407. https://doi.org/10.3390/nano13172407

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop