Abstract

A review is given on the application of the reflectance ellipsometry for optical characterization of bulk materials and thin films with thickness between /20 and 2 (at  nm). The knowledge of the optical constants (refractive index, , and extinction coefficient, ) of thin films is of a great importance from the point of view of modelling and controlling the manufacture of various optical elements, such as waveguides, diffraction gratings, and microlenses. The presented results concern the optical properties of thin films from multicomponent chalcogenide glasses on the base of As2S3 and GeS2 determined by multiple-angle-of-incidence ellipsometry and regarded as a function of the composition and thickness. The homogeneity of the films is verified by applying single-angle calculations at different angles. Due to decomposition of the bulk glass during thermal evaporation, an optical inhomogeneity of the thin As (Ge)-S-Bi(Tl) films is observed. The profile of in depth of thin As-S-Tl (Bi) films was investigated by evaporation of discrete layers. It is demonstrated that homogenous layers from the previous compounds with controlled composition can be deposited by coevaporation of As2S3 and metals or their compounds (Bi, Tl, In2S3).

1. Introduction

The optical constants of the thin films are commonly determined by methods such as spectrophotometry, prism-coupling technique, and ellipsometry. The ellipsometric methods are fast, nondestructive, highly sensitive to the presence of surface inclusions and provide higher accuracy for determination of the refractive index and thickness of thin layers in comparison with the spectrophotometric methods [14]. Currently various forms of ellipsometry exist and are applied for investigation of the physical properties of thin films [3]. The most common among these methods are the spectral ellipsometry and multiple-angle-of-incidence ellipsometry [5]. Spectral ellipsometry includes measurement and analysis of the complex reflection ratio in wide range of wavelengths at a fixed angle of light incidence. The method provides information for the optical response functions (complex refractive index, , or complex permittivity, ) [6]. The main problem of the spectral forms of ellipsometry is the selection of an appropriate dispersion model for description of the optical constants [6]. Single-wavelength ellipsometry offers an advantage in this respect since it does not require a dispersion law and enables the in situ control of the process of thin film deposition [5].

Chalcogenide glasses exhibit transmittance in wide spectral range (1–20 μm), low phonon energy and bandgap, high linear and nonlinear refractive index, and high photo-sensitivity when illuminated by near bandgap light [79]. Mixing chalcogenide glasses with a relevant oxide allows extended variation of the refractive index, for example, in GeO2-GeS2 system [10]. The chalcogenide glasses can be deposited as thin films by a variety of techniques as thermal and flash evaporation [11, 12], spin-coating [13] and pulsed laser deposition [14]. As a result, they have found practical applications in manufacturing of various elements for optics in visible and infrared spectral region [15], photonics [16], and all optical processing [17] such as waveguides [18], optics amplifiers [19], diffraction gratings [20], microlenses [21], and phase-change recording media [22]. To realize all these devices and applications, one should know and be able to precisely control the optical parameters of the materials. Ellipsometric measurements have been applied for investigation of the optical properties of chalcogenide films [2335] and of different photo-induced effects as photocrystaleization [29, 31], volume changes [26, 34], and photodissolution of metals [35]. The comparison of dispersion models for analysis of the measured spectra of ellipsometric angles (Ψ and Δ) of thin As2S3 was made in [36]. The authors [37] found that the model of Urbach-Cody-Lorentz and parameterization of the density of the electronic states model give the best fit of spectrophoptometric (transmittance and reflectance) and/or ellipsometric experimental data for thin As-S films.

This is a review paper on the optical properties and their changes after illumination of thin chalcogenide films from As-S(Se), As-S-Ge(Se), and As(Ge)-S-Tl(Bi, In) as studied by means of multiple-angle-of-incidence ellipsometry. We discuss the applicability of multiple-angle-of-incidence ellipsometry (MAI-ellipsometry) for characterization of chalcogenide layers from the above systems exclusively when irregular gradient of the refractive index is involved, resulting from compositional variation in the depth of the layers. The possibility to deposit homogeneous multicomponent chalcogenide films by coevaporation is demonstrated.

2. Experimental Details

Bulk glasses from the system AsxS100−x (), AsxSe100−x (), GexS100−x (), As40S60−xSex (), and (As2S3)100−xBi(Tl)x () were synthesized in a quartz ampoule from elements of purity 99,999% by the method of melt quenching [25, 26, 38]. Thin films with thickness between 30 and 1000 nm were deposited from a Mo or Ta boat with a preliminarily calculated amount of bulk materials for each particular thickness. The layers were deposited by thermal evaporation at deposition rates of 0.4–0.7 nm/s on optical glass BK-7 substrates. The films were exposed to white light (Mercury lamp, 20 mW/cm2). The composition of the bulk samples and thin films deposited on graphite substrates was determined by scanning electron microscopy with an X-ray microanalyzer (Jeol Superprobe 733, Japan) [39]. The composition of the films from binary and ternary As-S-Se and As-S-Ge systems [25, 26] is close to that of the bulk samples, while it is observed that the amount of Bi or Tl in the thin As-S-Bi(Tl) layers is smaller compared to that in the bulk samples [39].

The ellipsometric measurements were carried out at three angles of light incidence: 45, 50, and 55° using LEM-3M1 ellipsometer equipped with low-intensity He-Ne laser ( nm). The measurements were performed at angles of incidence below the Brewster’s angle in order to remove reflection from the back side of the substrate. The accuracy of determination of and was 0.01° and 0.02°, respectively. The systematic errors were eliminated through averaging of the four zones measurements. The duration of each ellipsometric measurement was about 5 min, and all measurements were carried out within 24 hours after the films deposition. Thus, the fast relaxation in selenium rich films [29] and oxidation of germanium containing films [40] would have not produced an artefact.

3. Theory and Modelling

Reflective ellipsometry is applied for quantitative estimation of an optical system by the so-called relative coefficient of reflectance, [1]. The basic equation of the ellipsometry is where and are complex reflection coefficients for - and -polarization, respectively. is the ratio amplitude of the reflection coefficients and and is the difference between their phases: where and are the modulus of and and and are their phases. In the case of a three-phase system (ambient, thin film, substrate) (Figure 1(a)), (1) can be written as: where , , , and are the Fresnel coefficients for ambient, thin layer and thin layer, substrate boundaries, respectively, and is the phase thickness: where is the complex refractive index of the thin layer with thickness . and are the angle of incidence and wavelength of the light, respectively.

A schematic model of an inhomogeneous layer is presented in Figure 1(b). In this case the Fresnel coefficients in (3) for ambient, thin layer and thin layer, substrate boundaries will include different refractive indices for the thin film, and , respectively [41].

In (4) for the phase thickness, , the refractive index and angle of incidence are generalized by their averaged values, the product of which satisfies the following integral: Two-layer models are usually used for analysis of roughness and surface phenomena as oxidation. Rough-surface models involving ridges, pyramids, or other similar geometric forms have been used to simulate surface roughness [42]. For the selection of the most appropriate model, where no data is available for the profile, best-fit result out of the experimental set of values of and should be used.

The dependence of and on the system parameters in (3) for a three-phase system can be expressed by the following analytical function: where is the refractive index of ambient (in case of air, = 1), and and are the refractive index and extinction coefficient of the substrate, respectively. An essential problem of all forms of ellipsometry is the impossibility of inverting (6) in order to obtain the complex refractive index parameters. The methods applied for solving the system of (6) using multiangle measurements were described in [1]. In many cases these parameters are varied in order to minimize some function. One of the most often used estimators is where is the number of the experimental points (angles of incidence in our case), and are the theoretical values of the ellipsometric angles at a determined group of optical parameters for the applied model, and and are the measured values. The fitting process includes a choice of suitable initial values for the unknown parameters and the minimization of the function using the fitting procedure based on the Levenberg-Marquardt algorithm. The initial estimates of the parameters were simultaneously varied until a fit was obtained between the experimental data and the values calculated from the thin film formula (3). The good coincidence between the measured and theoretical values of the ellipsometric angles within the defined accuracy is an evidence of felicitous choice of a model of the investigated system [5, 43]. The error is calculated from the elements of the correlation matrix which the Levenberg-Marquardt algorithm produces during the minimization procedure. The calculation procedure is described in [1, 5].

4. Results and Discussion

Knowledge of the optical parameters of the single layers is important for successful engineering of multilayer structures. Firstly, we determined the optical parameters of the thin films using the single homogenous absorbing layer model and applying multiangle measurements (45, 50, and 55°). In this case the refractive index, , extinction coefficient, , and thickness, , are simultaneously calculated (three-parameter calculation). We found that the maximal experimental error of was [26]. The basic advantage of the multiangle measurements consists in their independence on the half-wave point of light incidence angle [44]. Further, we used the calculated optical parameters from multiangle measurement as an initial estimate of the unknown parameters for single-angle calculations in order to determine the degree of inhomogeneity of the thin layers. We searched for two optical parameters of the thin film and the third one has been iterated with small steps around its approximate value until minimum value of was found [26]. The single-angle calculation allows the film’s thickness to be calculated with an accuracy better than 0.3%. We showed that in case of homogenous thin films, both single- and multiple-angle methods gave the same results.

In Table 1 a comparison of the thickness and refractive index values of thin films calculated from ellipsometric measurements and/or from their transmission spectra using Swanepoel’s method [45] is presented. The measurements were performed on one spot of each sample. The areas of the incident probe beams of the spectrophotometer and the ellipsometer were 50 and 1 mm2, respectively.

For thin As-S films the values for the refractive index within the limits of the experimental error obtained through ellipsometry are in a good agreement with those calculated from spectrophotometric measurements. The single-angle calculations at different angles show that thin films from As-S(Se), Ge-S, As-Ge-S systems are homogenous or their inhomogeneity cannot be revealed since the discrepancies are within the experimental error. It is seen that the optical constants and thickness for layers from Ge-S-Ga(In, Tl, Bi) and As-S-Bi(Tl) systems as derived by ellipsometry differ from those calculated from spectrophotometric measurements. At the same time the difference in the values of independently determined at 45°, 50°, and 55° is higher than the experimental error . It is seen from (4) and (5) that in case of inhomogeneous coatings, due to different paths of the laser beam of the ellipsometer, different values for would be observed. These considerations suggest inhomogeneity of the thin films deposited form bulk glasses from As-S-Bi(Tl) systems. The derived values of , thus, can be regarded just as effective values. We used the method of thermal coevaporation to deposit thin films from ternary As2S3-In2S3 system [46] and thin films from multicomponent GeSe2-Sb2Se3-AgI system [47]. In Table 1 the refractive index of thin As2S3/Bi, As2S3/Tl, and As2S3/In2S3 films deposited by coevaporation is given. It is seen from single-angle calculations that the co-evaporation from two sources makes possible the deposition of homogenous coatings from the above multicomponent glasses.

Our results for the optical constants and thickness of thin films from As-S-Ge, AsxSe100−x, and GexS100−x systems and their changes after illumination by halogen lamp are presented in Table 2. It is seen that the refractive index of the chalcogenide films varies in wide diapason, that is, from 2.05 to 3.15 for GeS2 and As2Se3, respectively. The results for and of the films from As-S=(Se) systems are in agreement with the data reported in the literature [29, 4850], With the substitution of Se for S in thin As-S-Se layers, the refractive index of the layers increases and the highest value is reached for As2Se3 film. The changes in the optical parameters of thin As-Se films under illumination are more pronounced compared to those in As-S layers. Maximum changes in the refractive index of the As-Se films were observed for the thin film with composition As50Se50 (). The refractive index of the thin films from Ge-S system decreases after exposure. The greatest changes were observed for thin film with stoichiometric composition GeS2.

The photo-induced volume effects in thin chalcogenide films are connected with changes of their thickness. From practical point of view, volume changes make chalcogenide layers perspective materials for fabrication of microlenses [21]. An effect of photocontraction was found after exposure to light of thin films from As-S(Se) (see Table 2). In Figure 2 the relative changes of the thickness, as a function of the composition for As-S, Ge-S, As40S60−xSex, and As-S-Ge films are presented. It is seen that the layer from As-S and As40S60−xSex systems exhibited effect of photocontraction. The changes in the thickness increase with sulphur contain in layers from binary As-S system. A pure effect of photoexpansion is observed for thin Ge-S layers after illumination, while for As-S-Ge system, both effects of photocontraction and photoexpansion are presented. We have obtained the maximal value of photoexpansion for thin film with composition As10Ge30S60 (). We have obtained the maximal value of photoexpansion for thin film with composition As10Ge30S60 (). This result is in agreement with data published by other authors [51].

To understand the origin of inhomogeneity of the thin films, we investigated the thickness dependence of the refractive index. For the most photonic devices, it is necessary that the thin layers are with an exact thickness. The accumulated knowledge of the refractive index and its photo-induced changes would support the modelling and manufacturing of 1D photonic crystal [52]. The results from our investigations on the thickness dependence of the optical properties of thin films from As-Se, Ge-S, As-Ge-S, and As-S-Se systems [2527] showed that the refractive index is independent of for  nm. For acquiring of information about the possible changes of the optical constants during the process of measurement with He-Ne laser, we compared the values of , , and for thin films obtained by MAI-ellipsometry with those determined by two and three spectrophotometric measurements [53] (Table 3).

In Figure 3(a) a theoretical, as well as, experimental - dependence at room temperature is given for thin As2S3 films ( nm) with refractive index at = 632.8 nm and angle of light incidence . It is seen that the experimental results follow the theoretical curve for  nm. The thickness dependence of the ellipsometric angles - of thin films with compositions As40S30Se30 and As40Se60 is presented in Figures 3(b) and 3(c), respectively. The parameters of theoretical curves are and or and at  nm and angle of light incidence for thin films with composition As40S30Se30 and As40Se60, respectively.

It is observed that is independent of the thickness for  nm. We believe that the optical properties of the films thinner than 50 nm are strongly influenced by the substrate and the process of the thin films growth. The single homogenous absorbing layer model is only an approximation, since it does not consider either the existence of the substrate/film interface, the surface roughness, or oxidation. We can expect that the error of determination of increases extremely for  nm [54, 55]. The relative error in the thickness determination is independent of for layers with  nm (), and downwards it increases up to 1% for layers with d30 nm. The error of the refractive index determination is weakly affected () by the layer’s thickness [24].

It is seen from Table 1 that the composition of thin films from As-S-Bi is different from that of the respective bulk sample. We made the composition analysis of the residue in the boat that was found to be As0.5S59.3Bi40.2 after evaporation of thin film from bulk As38S56Bi6 glass. That is why firstly we analyzed the thin As-S-Bi films using the two-layer model considering an absorbing sublayer on the upper part of the coating (closer to air) when the whole quantity in the boat was evaporated. This model gives the best fit of the experimental results (Figure 4). The results show a higher extinction coefficient, , of the upper sublayer. The thicknesses of the layers for thin films with compositions As39S58Bi3 and As34S54Bi10 system were 39 and 70 nm, respectively. The X-ray microanalyses at different acceleration voltage [39] showed that concentration of Bi or Tl increases in direction from the substrate to the interface of the layer with the ambient. In [56] it was found by mass-spectroscopy analysis that the vapour of As2S3, formed in the process of thermal evaporation, consists of different sulphur- or arsenic-rich fractions. These fractions polymerize on the substrate forming the glass network of the thin film. This finding suggests that the melt in the boat should become poorer of light fractions during the thermal evaporation, and therefore the upper part of the thin films should be enriched with structural units formed from heavier fractions. Most probably this is the reason why in some cases inhomogeneous layers are produced by thermal evaporation. This phenomenon is more pronounced for thin films from the ternary As-S-Tl(Bi) systems [39, 57].

Usually inhomogeneous thin films are presented as multilayer coatings where the refractive index is assumed to be variable under some mathematical function—linear, exponential and Gaussian [58]. In many cases, it is not possible to determine the manner in which the refractive index of the film varies. For estimation of the optical homogeneity in the volume of the thin layers, thin films (about 80–100 nm thick) were deposited consecutively on different substrates until the whole quantity in the boat was evaporated. We expect that the dependence of the refractive index of the discrete layers on the order of their deposition would be the same as the profile of in the one micron integral coating. Figure 5 represents the distribution of the refractive index, , on the total thickness of the sublayers for unexposed As40S60 and As38S56Tl6 layers. It is seen that the refractive index for as-deposited As40S60 films does not change with the stack order, while considerably increases for the last two numbers of the system As38S56Tl6. It can be concluded therefore that the As2S3 thin films are homogeneous in composition in depth of the films, while during the process of vacuum deposition of As-S-Tl coatings, a process of decomposition takes place. This fact suggests that we can consider As38S56Tl6 film to be composed of two parts: homogenous sublayer with low refractive index and inhomogeneous sublayer on the upper surface of the coating. The X-ray microanalysis of the composition of the sublayers revealed variation between 0.1–0.6 at % of the thallium content in the samples from 1 to 6, while the compositions for the sublayers 7 and 8 were As32.6S62.3Tl8.1 and As5.6S59.5Tl34.9, respectively. It is known that in case linear profile of the average values for refractive index, , is equal to the refractive index at the thickness level 0.5  [59]. The dependence presented in Figure 5 can be extrapolated by exponential function of the type , where is the value of the refractive index of the film close to substrate. Having determined the optical parameters and thickness of the single layers, we calculated the theoretical ellipsometric angles Ψ and Δ of the complex structure [60]. Further to finding the effective value of we used these values to calculate refractive index of such coating as the integral film. We found that the effective value of was 2.679 that is close to the one determined for one-micron layer evaporated from bulk glass (see Table 1).

The thin films deposited from bulk As38S56Bi6 glass showed similar dependence of on their deposition number (Figure 6). This dependence suggests that metal content increases extremely in the last two samples. The calculated values for for is lower than those calculated from two-layer model (see Figure 4(b)) since in the former case it is averaged over higher thickness ( nm).

In the existing methods the profile is determined using virtual slicing of the thin film into sublayers, suggestion of some function of the dependence of the refractive index in the depth of the coating (linear, exponential, or polynomial) and minimization of the discrepancy between theoretically calculated and measured values of and [61]. In contrast, the previously suggested method requires no mathematical function and provides reliable information for the refractive index profile in the depth of thin films. The gradient of due to variation of composition in the depth of the layer could replace the design of column structures in the construction of multilayer structure [62, 63].

5. Conclusion

The surface roughness plays a key role in the formation of the interface between layers in multilayer structures. The effects of surface oxidation observed in thin films [40] could be determining in the formation of this interface and their presence indirectly demonstrated. In [24] it was shown that the introducing of a surface layer with considerably improves the accuracy of the fit of the experimental data for Ψ and Δ for As2S3 film. The thickness of the surface overlayer is in the range of 2–15 nm [26, 28]. Our investigation on surface roughness by White Light Interferometric Profiler [64] shows that the root mean square (rms) is about 1.0–8.6 nm for As-S-Se/Ag layer with thickness about 700 nm. In spite of the great number of papers dealing with surface phenomena in chalcogenide thin films, the investigations by ellipsometric measurements are still in the beginning.

In the present paper the optical properties of thin films on the base of As2S3 and GeS2 as a function of the composition and thickness were investigated by MAI ellipsometry. The optical parameters of the thin films from the previous systems were determined with an accuracy of 0.2-0.3% for thickness in the range from /20 to 2. A good agreement was demonstrated between the values of the optical constants as obtained by ellipsometric and spectrophotometric measurements in the case of homogeneous thin films from As-S(Se) and As-S-Ge systems. The possibility is shown for the refractive index to be varied in wide range from 2.05 to 3.15. The inhomogeneity of thin As-S-Tl(Bi) films leads to discrepancy of the results derived from both types of measurements. The refractive index profile of thin As38S56Tl6 and As38S56Bi6 films is investigated. It was found that the variation of the refractive index with the distance from the substrate to the surface was about 0.8 and 0.1, respectively. This result presents an opportunity for utilization of inhomogeneous layers from As-S-Tl(Bi) system in photonic multilayer structures. The possibility for deposition of homogenous As-S-Tl(Bi, In) layers by coevaporation of As2S3 and metal compounds is demonstrated.

Acknowledgment

The authors would like to acknowledge the financial support from Project BG05PO001-3.3-005/0001 “Science & Business” branch, operational program “Human Resource Development.”