1. Introduction
The universal dispersion model (UDM) describes individual elementary electronic and phonon excitations in materials as separate contributions [
1]. This dispersion model satisfies three fundamental conditions: time reversal symmetry, Kramers–Kronig consistency and conformity with the sum rules. The UDM is a suitable dispersion model for a very precise description of the response function of optical materials in a wide spectral range [
2,
3,
4,
5,
6,
7].
Optical materials used in thin-film interference optics often have a polycrystalline structure [
8,
9,
10,
11]. For optical purposes, it is necessary to prepare polycrystalline thin films with a very fine structure, where the mean size of the grains is much smaller than the wavelength of light for which the system is designed. However, this is difficult to achieve in practice. The size of the grains accompanied by surface roughness is sensitive to deposition conditions and often changes during the growth [
12,
13,
14]. In the first phase of deposition, polycrystalline material with very fine grains grows. The size of the grains then gradually increases with the thickness of the film, until the film becomes optically unusable. The internal scattering in the film (extinction) increases with the size of the grains and the scattering losses at the interfaces also increase due to the roughness [
15,
16]. In addition, films with large grains show a porosity that negatively affects the extinction coefficient in the transparent region due to the presence of localized states [
17,
18,
19]. The localized states originate in the surface states of the grains as well as in the adsorbed substances in the pores, which constitute unstable components depending on the environmental conditions. All these defects negatively affect the properties of the films and layered systems containing these films. This is why only limited work has been performed so far and the optical constants of many fluoride films are thus known only in the transparent region with large uncertainty [
20].
This article will focus on GdF3 material that is used as a high-index material for multilayer coatings in the UV region due to its high band-gap energy. It will be shown that GdF3 films exhibit the defects mentioned above and that the UDM is an effective tool by which to describe the optical properties of these films.
3. Experimental Arrangement
Experimental data were acquired using two ellipsometers and three spectrophotometers: Woollam IR-VASE, Horiba Jobin Yvon UVISEL, Bruker Vertex 80v, Perkin Elmer Lambda 1050 and McPherson VUVAS 1000. These data covered a wide spectral range of 0.00868–10.3 eV (70–83,300
; 0.12–143
). The measurements obtained using individual instruments in individual spectral ranges resulted in 16 data sets for each sample. Thus, a total of 80 experimental data sets were used in the optical characterization (for further details see the
Supplementary Materials).
The ellipsometry was represented by the vector of associated ellipsometric parameters
[
21]. The associated ellipsometric parameters correspond to the three independent elements of the normalized Mueller matrix of isotropic systems:
where
R is the average reflectance (
T instead of
R appears for the transmitted light). Ellipsometric data were measured from both sides of the sample in the infrared region, where the substrate is transparent. The difference ellipsometry data
were used in addition to the ellipsometric data
measured from the side with the GdF
3 film (front side) and the data
measured from the opposite side (back side). Although these data are not independent, the inclusion of difference ellipsometry in the data processing helps to reduce the influence of the substrate and systematic errors [
2,
5]. In the UV and VIS spectral range, where the silicon substrate is not transparent, the ellipsometry was also measured from the back side in order to determine the thickness of the native oxide layer. For this purpose it was sufficient to make measurements for a relatively small number of spectral points at one angle of incidence.
As for spectrophotometry, both absolute quantities (reflectance and transmittance) and relative reflectance were processed. The relative reflectance was calculated as the ratio of the reflectance
from the front side and the reflectance
from the back side of the sample:
According to our experience, similar to difference ellipsometry, the relative reflectance helps to reduce the correlations between the sought parameters by compensating for systematic errors in measurement [
2,
5].
5. Structural Model of Samples
The samples consisted of double-sided polished silicon single-crystal substrates with the GdF
3 films deposited on one side (see
Figure 1). Since the substrates were exposed to air, it was assumed that the back sides were covered by thin native oxide layers (NOLs), which were modeled as thin homogeneous films with thickness
. Since no special effort was made to remove the NOLs from the surfaces of the silicon substrates before the deposition of the films, it is possible that parts of the NOLs remained below the GdF
3 films. Moreover, the treating of substrates by the oxygen plasma prior to the deposition of the films could also be a contributing factor to the NOLs. In the structural model, this is represented by transition layers between the GdF
3 films and the silicon substrates, which were modeled as thin homogeneous films with thickness
and with optical constants identical to the NOLs. The optical constants of float-zone silicon substrates and the NOLs were modeled by UDM with parameters fixed in values corresponding to our previous study of crystalline silicon wafers [
6,
26] (for optical constants see the
Supplementary Materials). In general, the optical response of the NOLs is similar to amorphous
. The exact values of the refractive index play only a minor role because of the correlation with the NOL thickness.
The GdF
3 films with nominal thickness above 20 nm were modeled using inhomogeneous films with randomly rough upper boundaries (see
Figure 1). Since the films were deposited on the polished substrate, it was assumed that the lower boundaries of the GdF
3 films as well as the boundaries of the transition layers were smooth. In the structural model, the inhomogeneous film with thickness
is approximated by a stack of
thin homogeneous films. The mean thicknesses of the films in this stack are equal to
with the exception of the first and last films, which have a mean thickness
. The value of
m was chosen differently for each film such that the ratio of nominal thickness and
m was 20 nm. The thinnest films with nominal thickness 20 nm were assumed to be homogeneous. Moreover, because it was found that the roughness of the thinnest films was very small (see
Section 7,
Table 1), they were not considered in the structural model.
The optical constants
of the films in the stack were chosen as
where the index
numbers the films and
represents dependence of the dielectric function on the coordinate
z perpendicular to the surface of the substrate. The arrangement with the first and last films in the stack having half the thickness of the other films is advantageous, since the convergence to the exact result corresponding to the films with a continuous profile of the optical constants (i.e., to the limit
) is faster than if all the films had the same thickness [
27]. All the boundaries in this stack were assumed to be rough, except for the boundary at
, which separates the GdF
3 film from the transition layer. The irregularities in the roughness were identical in shape, but their heights were scaled proportionally to the distance of the mean planes from the boundary at
. Therefore the root-mean-square (RMS) value of the heights of irregularities of the
j-th boundary is
where
denotes the positions of the mean planes and the value
corresponds to the roughness on the surface of the GdF
3 film.
The Fresnel coefficients of the layered systems with randomly rough boundaries were calculated using the Rayleigh–Rice theory (RRT). The reflection coefficients are expressed as
where
distinguishes the p and s polarized waves, the symbol
denotes the reflection coefficient calculated as if all the boundaries were smooth and the symbol
represents the correction calculated using the RRT. It can be written in the following form [
28,
29]
where the summation is over the rough boundaries, the variables
and
represent the spatial frequencies, the symbol
is the size of the ambient wavevector and
denotes the incidence angle. The symbol
denotes a complicated function, which apart from the indicated dependence on
and
depends on the wavelength of light, incidence angle, optical constants and mean thicknesses of the films, but does not depend on roughness. The information about roughness is expressed using the power spectral density function PSDF
, which, in our model, is given as
where
is the autocorrelation length. The transmission coefficient can be expressed by formulae analogous to (
6) and (
7). Details concerning the calculation of the reflection and transmission coefficients of multilayer systems with rough boundaries using the RRT can be found in [
29].
The reflection and transmission coefficients of the native oxide layer on the back side can be expressed using one of the standard methods. The optical quantities describing the whole sample are then calculated from the reflection and transmission coefficients calculated for layered systems on the front and back sides of the substrate. With the exception of the spectrophotometry in the FIR region measured with high spectral resolution, the thickness of the substrate is much larger than the coherence length. Therefore, the influence of the substrate was included using the method working with incoherent light described in [
21], which utilizes the Mueller matrices. A formalism working with coherent light had to be used in the case of FIR measurements where interference in the substrates was observed.
6. Dispersion Models
Isotropic media without spatial dispersion were assumed. The universal dispersion model (UDM) [
2,
3,
4,
5,
6,
7] was used for all the media in the system, i.e., for GdF
3 films, crystalline silicon substrate c-Si and for the native oxide layer (NOL). One of the advantages of the UDM is that the optical constants of various materials can be easily shared between different studies using a set of parameters instead of describing the optical constants using tables.
The dispersion model of GdF
3 films consists of three contributions representing interband electronic transitions, absorption involving localized states and phonon excitations:
6.1. Interband Electronic Transitions
Most of the interband electronic transitions from the occupied valence band to the unoccupied conduction band occur outside the experimental spectral range because GdF3 has a large band-gap energy .
The interband electronic transitions were modeled using two contributions
The first contribution describes a broad absorption band that has only two parameters, i.e., band gap energy
and transition strength
. The imaginary part of the response function is calculated as:
where
is the Heaviside step function. This function has quadratic dependence above the band gap energy and classical
asymptotic behavior for large values of photon energy. The real part of the response function is expressed on the basis of the Kramers–Kronig integral as follows:
The second contribution describes a relatively narrow absorption excitonic structure near the band gap. This contribution was modeled using the Campi–Coriasso model [
1,
30,
31,
32]:
where parameters
and
determine the width and central energy of the excitonic structure. The model contains two other parameters, i.e., band gap energy
and transition strength
. The band gap energy parameter
is shared with the first contribution. The real part of the response function is also expressed on the basis of the Kramers–Kronig integral and it can be written using the closed form expression [
32].
It is practical to control the strength of the interband transitions between the valence and conduction bands with only one density parameter
, which is achieved using the relations
The parameters
and
are the intrinsic parameters of the dispersion model, while the parameters
and
are the parameters used in the fitting of the experimental data. This trick is important for modeling films with refractive index profiles, because while the density parameter
depends on the coordinate perpendicular to the film surface, the parameter
is held constant in the profile. This reduces the number of sought parameters. In this work, the linear profile of the density parameter was assumed
where
is the contribution of interband transitions to the susceptibility normalized by the sum rule integral such that:
This contribution introduces
z-dependence of the dielectric function needed for expressing the optical constants
of the thin films in the stack representing the rough inhomogeneous GdF
3 film (
4).
6.2. Absorption Involving Localized States
The adsorbed components in the pores give rise to many absorption centers resulting in absorption bands at energies lying below the band gap energy of the GdF
3 host material. This absorption can be interpreted as excitations of the localized states of electrons with energies lying inside the band of forbidden energies. The response function can be effectively described as the sum of Gaussian broadened discrete transitions (five transitions were used in the case of the GdF
3 films) as follows [
1,
33,
34,
35,
36]
where
,
and
are the transition strengths of excitations, mean energies of the excitations and the Gaussian broadening parameters (RMS value). The symbol
denotes the Dawson function (integral) [
37].
Relation (
17) includes only the excitation of electrons from occupied localized electron states to unoccupied localized states, but localized states also contribute to the absorption involving delocalized (extended) states, i.e., to the so-called Urbach tail. The Urbach tail can be explained as the transitions of the electrons from the localized valence states to the extended unoccupied conduction states and the transitions from the extended valence states to the localized unoccupied states [
38]. In other words, the Urbach tail is the transition region between the absorption on localized states and the region of interband transitions. The Fermi energy is assumed to be at the center of the band of forbidden energies. Thus, the minimum excitation energy, which is the half of the band gap energy, corresponds to excitations from the occupied states at the Fermi level to the bottom of the conduction band or from the top of the valence band to the unoccupied states at the Fermi level. In our model, no absorption is assumed below
, then the exponential is used to model the absorption up to the band gap energy
and, finally, the region above
is modeled by a rational function giving classical asymptotic behavior at large energies [
7]
where
is the Urbach energy. The parameters
,
and
must be chosen so the function
is continuous up to the second derivative at
and simultaneously the sum rule integral is equal to the transition strength
, which together with
and
constitute parameters of this model. The real part of the response function is calculated from the imaginary part using the Kramers–Kronig integral [
7].
The parameters expressing the transition strengths on localized states
and the transition strength of the Urbach tail
are internal parameters linked to the fitted parameter
in the following way:
where parameters
are fitted parameters expressing relative strengths. The total transition strength associated with the localized states is then
In the model of inhomogeneous GdF
3 films, it was assumed that the bottom parts of the films are so dense that there are no adsorbed components. The profile of the density of localized states was thought to be linear, thus, the response function was calculated as
where the transition strength parameter
is proportional to the density of the localized states on the top of the film and
is the normalized contribution defined in same way as
in (
16).
6.3. Phonon Excitations
All phonon excitations are modeled by a discrete spectrum broadened by the Voigt broadening function [
1,
6]. Thus, the contribution of the phonon absorption peaks to susceptibility can be written using the complex Faddeeva function
[
39,
40,
41,
42,
43] in the following form
where
and
are the root-mean-square (RMS) width of the Gaussian part and the full-width at half-maximum (FWHM) of the Lorentzian part of the Voigt profile. These parameters are calculated from the FWHM of the Voigt profile
using the following approximate equations
The parameter
takes values between 0 (fully Gaussian model) and 1 (fully Lorentzian model). Because Equation (
22) cannot be used with
, an alternative form based on the damped harmonic oscillator (Lorentz) model must be used for fully Lorentzian contributions
Each absorption peak depends on four parameters: the transition strength
, the phonon excitation energy
, the FWHM of the peak
and the weight between the Gaussian and Lorentzian part
.
The first seven phonon excitations represent one-phonon absorption processes in polycrystalline GdF
3 host material, i.e., the most pronounced absorption processes in the films characterized by sharp structures. In the same region, two broad phonon absorption bands were added for modeling one-phonon processes that correspond to a response from the GdF
3 disordered structure. Two additional peaks in the spectral range above the one-phonon absorption were used for the weak structure modeling multi-phonon processes. The other 19 phonon excitations describe the vibration spectra of adsorbed components. For the description of the one-phonon absorption processes the asymmetric peak approximation was used; thus, these seven absorption structures were modeled by Equation (
22) with the following substitution
where the parameters
determine the asymmetry of the absorption peaks and
is the normalization constant. To ensure the physical correctness of this model, it is necessary to assume that the parameters
are not independent. If
is chosen as the dependent parameter, its value must be calculated as
The normalization constant is calculated as follows
These conditions ensure the validity of the sum rule in the following form
Note that the asymmetric peak approximation may give a nonphysical response function with negative values of
for large values of the
parameters. It is difficult to formulate strict criteria for
parameters that ensure the non-negative values of
. Therefore, this problem was solved by introducing a penalization function that increased the residual sum of squares if
was negative in some part of the spectra. Details can be found in the
Supplementary Materials.
In a similar way as for the electronic excitations, the profile of the phonon response function was assumed. In the case of the sharp seven one-phonon excitations, the linear profile of the dispersion parameters
,
and
was introduced
Due to the dependency on the parameters
and
, it is not possible to define a normalized phonon response function similar to the electronic case. In other words, the distributions of the phonon excitations are
z-dependent. The remaining phonon excitations’ only transition strength parameters are
z-dependent, except for very weak multi-phonon excitations, which are assumed to be homogeneous.
Because the positions and widths of phonon excitations are usually characterized using wavenumbers, the parameters and specified in were used instead of the parameters and specified in eV.
7. Results
In the following subsections the results of the optical characterization are presented. Only the plots of experimental data relevant for the presented discussion are shown here; however, the plots of all the experimental data and their fits can be found in the
Supplementary Materials. It should also be noted that the uncertainty of the determined parameters is shown using concise bracket notation. For example, the thickness of the thickest film is determined in value of 611.75(6) nm, which represents the uncertainty interval
with 68% confidence. This uncertainty includes only random measurement errors and does not include an estimate of systematic errors.
7.1. Film Structure
In
Table 1, the structural parameters of the GdF
3 films are listed. The film thickness was monitored using the QCM method, which measures the mass of the deposited material. Therefore, the small differences between the nominal thickness and the thickness determined in the optical characterization can be explained by the density profile of the film.
The roughness determined in the optical characterization was verified by atomic force microscopy (AFM). The RMS values of the heights of irregularities (roughness) and the autocorrelation lengths obtained by the statistical analysis of the AFM scans are also introduced in
Table 1. The statistical analysis was performed using Gwyddion open-source software [
44].
It is evident that the roughness parameters determined in the optical characterization correspond well with those determined in the AFM study. The AFM study justifies the use of the same autocorrelation lengths for samples #1–#4 and the assumption of the smoothness of sample #5 in the optical study. The autocorrelation lengths determined in the AFM study exhibit slow growth with film thickness. However, because the autocorrelation length is correlated with the RMS value of the heights in the optical study, only one autocorrelation length is used for all the samples.
7.2. Electronic Excitations
The density profile of the porous GdF
3 film is in principle proportional to the transition strength profile of interband excitations
.
Figure 2 shows the transition strength profiles of all five samples calculated on the basis of the values listed in
Table 1 and
Table 2. The rough surfaces are indicated by zigzag lines with a range corresponding to
in this figure. It is evident that for a given value of
z, the transition strength (density) of the thicker films is larger than for the thinner films. This can be explained by assuming that the material is not deposited only on the tops of the films but also partially into the pores during the growth of these films. The estimates of the packing densities on the tops of the films with respect to the bottoms of the films calculated as the ratio
are shown in
Table 2. These estimates show a decrease with film thickness. This is not surprising because the size of crystals in this polycrystalline film grows with thickness, which results in larger pores.
The spectral dependencies of the dielectric function on the top and bottom of the thickest inhomogeneous GdF
3 film are shown in
Figure 3. The excitonic structure centered around 12.2 eV, which is described by the parameters listed in
Table 2, lies beyond the spectral range of the experimental data. Moreover, only reflectance data measured by the VUV spectrophotometer were available above 6.5 eV. Therefore, the credibility of the optical constants in this spectral region, and especially above 10 eV where the spectral range of the VUV spectrophotometer ends, is limited.
The reflectance spectra measured at near normal incidence using the VUV spectrophotometer under vacuum and using the standard instrument in air are also shown in
Figure 3. Evidently, these data do not coincide very well, which is caused by the volatility of the components adsorbed in the pores. The adsorbed components are strongly influenced by the surrounding conditions. Even the optical measurements performed under the same experimental conditions on the same sample can differ in a significant way, because these components are not stable in time. This is evident from the two ellipsometric measurements taken on the same sample six weeks apart, which are shown in
Figure 4. The amount of adsorbed material is almost impossible to control; thus, the concentration of the volatile components was modeled separately for each experimental data point. Only one parameter, the transition strength of excitations involving the localized states
, was used for this purpose. As mentioned in the part devoted to the dispersion models, this parameter exhibits a linear profile, with zero at the bottom of the films. This means that while the dielectric function on the bottoms of the films is the same for all the samples, the dielectric function on the tops of the films is different not only for different samples, but it is slightly different also for different measurements (the one in
Figure 3 corresponds to the thickest film and measurement by the VUV spectrophotometer). The values of
are listed in
Table 3; the other dispersion parameters describing localized states, which are common to all samples and measurements, are listed in
Table 4. From the dates of measurements shown in
Table 3, it is evident that the amount of adsorbed material increased with passing time.
Strictly speaking, the Urbach energy should be used for a weak absorption below the band gap energy of solid materials (such as crystalline or amorphous semiconductors, glasses, etc.). The Urbach energy is in the order of tens of meV. In this work, the Urbach energy is also used more loosely for absorption on adsorbed components. From this point of view, the much higher value of the Urbach energy should not be surprising.
7.3. Phonon Excitations
The phonon absorption spectra of GdF
3 are located in the far-infrared region as can be seen in
Figure 5, where the high-resolution relative reflectance spectra of three selected samples are plotted. The spectral resolution of the FIR spectrophotometry data is 1 cm
−1; thus, the interference patterns originating in the substrate were measured. Therefore, the influence of the substrate must be taken into account by a partially coherent model. Note that the interference observed in the substrate allowed us to determine its thickness with high precision (see
Table 1). The evident evolution of the response function with sample thickness can be seen in the bottom part of the
Figure 5, which shows the imaginary parts of the dielectric functions in the FIR region. The shape of the phonon structure is influenced by the structural changes in the GdF
3 host material. It should be noted that in the spectral region below 250
, only the spectrophotometric experimental data were measured, which influences the accuracy of the determined dielectric function in this region. From the relative reflectance data for the thinnest GdF
3 film (sample #5), it can be seen that the one-phonon structures are at the level of the measurement noise. Therefore, it was assumed that the the phonon response function at the bottom of this film is the same as on the top of sample #4. That is why the response functions of the two thinnest films in
Figure 5 overlap. The dispersion parameters of the phonon contributions of the GdF
3 host material are listed in
Table 5 and
Table 6. The last lines of these tables show the total transition strengths of the crystalline and disordered phases and their sum. One can see the increase in the crystalline phase and the decrease in the disordered phase with increasing thickness of the films. The total transition strengths of both phases decrease with film thickness. This corresponds to the decrease in the packing density calculated as the ratio of the total transition strengths of phonon excitations of the given sample and that corresponding to sample 5, which is assumed to be completely dense (see
Table 6). The reason for the discrepancy between the packing densities calculated from the transition strengths of electronic excitations and phonon excitations could be that it was not possible to correctly determine the transition strength of the thinnest sample (as explained above), which is used as a reference with 100% packing density.
It should be noted that there are also two weak absorption peaks at wavenumbers 418 cm−1 and 473 cm−1 representing multi-phonon excitations.
In the mid-infrared region one can see the vibration spectra of adsorbed components. The spectral resolution 8 cm
−1 of the MIR spectrophotometry is much lower than that of the FIR spectrophotometry, and the interference patterns in the substrate are not observed in this case. The top panel of
Figure 6 shows the relative reflectance spectra for the three thickest films. The characteristic fingerprints of OH (water) and CH (organics) vibrations in the region from 2800 to 3600 cm
−1 are clearly visible in these spectra. Thus, it is apparent that the adsorbed components consist of water and organic molecules. In the region below 1800 cm
−1, not only the vibration spectra of the adsorbed components, but also the multi-phonon absorption in the silicon substrate can be seen. In the bottom panel of
Figure 6, the imaginary parts of the dielectric functions of the GdF
3 films are plotted. It is easy to identify the characteristic absorption peaks of water and organic components even in the region where they overlap with the multi-phonon absorption spectra of silicon. Since the focus of this work is on the optical characterization of GdF
3 host material and not on the exact characterization of the adsorbed components, only rough identification of the vibration spectra in the MIR region was performed.
The absorption structures of adsorbed components were modeled using 19 symmetric Voigt peaks with different transition strengths used for individual data sets. The strengths of the absorption structures were fitted independently for each measurement covering the MIR spectral region. In all other cases, these strengths were fixed in values determined for the MIR transmittance measurement. The other parameters of these peaks, i.e., the position, width and shape parameters, are common for all the experimental data.
7.4. Optical Constants
The optical constants of the host GdF
3 material, i.e., those corresponding to the bottoms of the films where they are assumed to be dense, are shown in
Figure 7. The tabulated data can be found in the
Supplementary Materials. The refractive index at 193 nm is 1.688, which is in good agreement with the values 1.68–1.69 published in [
20].
The tabulated data for the optical constants of float-zone silicon substrates and the NOLs can be also found in the
Supplementary Materials.
8. Conclusions
The optical characterization of gadolinium fluoride (GdF
3) films was performed using the universal dispersion model [
1] implemented in newAD2 software [
22]. In the framework of this dispersion model, the asymmetric peak approximation with the Voigt broadening function was used to describe the one-phonon absorption. The ellipsometric and spectrophotometric experimental data measured for five samples in a wide spectral range using several instruments were processed simultaneously. The nominal thicknesses of the films ranged from 20 to 600 nm. The films exhibited several imperfections, which were found to be more severe for thicker films; namely, surface roughness and refractive index profile. It was shown that the optical properties of the thicker films change over time due to the volatile adsorbed components. The adsorbed components were identified as water and organics from the IR absorption spectra. All these defects suggest that the films have a porous structure. For the precise characterization of the GdF
3 host material, it was necessary to include this instability of the optical properties together with the imperfection of these films in the model used in the optical characterization.
Although the films are far from ideal, the thinner films are not affected by these defects as much as the thicker ones; therefore, the optical properties obtained in this work can be used to describe these films in thin-film systems designed for the ultraviolet spectral region.
The presented method of optical characterization can be applied to a wide range of dielectric materials exhibiting porosity, surface roughness and refractive index profile, i.e., especially for other fluoride films used in multilayer coatings.