Temperature-Dependent Spectroscopic Ellipsometry of Thin Polymer Films

Thin polymer films have found many important applications in organic electronics, such as active layers, protective layers, or antistatic layers. Among the various experimental methods suitable for studying the thermo-optical properties of thin polymer films, temperature-dependent spectroscopic ellipsometry plays a special role as a nondestructive and very sensitive optical technique. In this Review Article, issues related to the physical origin of the dependence of ellipsometric angles on temperature are surveyed. In addition, the Review Article discusses the use of temperature-dependent spectroscopic ellipsometry for studying phase transitions in thin polymer films. The benefits of studying thermal transitions using different cooling/heating speeds are also discussed. Furthermore, it is shown how the analysis and modeling of raw ellipsometric data can be used to determine the thermal properties of thin polymer films.


INTRODUCTION
Thin polymer films have found many important applications in organic electronics, such as active layers, protective layers, or antistatic layers. A good, specific example of this is the wide range of applications for organic semiconductor poly(3,4ethylene dioxythiophene) (PEDOT) doped with poly(4styrenesulfonate) (PSS), which include energy conversion applications (solar cells), antistatic and conductive coatings, capacitors, touch panels, organic light emitting diodes, and printed organic electronics. 1−9 Other examples of widely studied thin organic layers are thin layers of poly(3hexylthiophene-2,5-diyl) (P3HT) and (6,6)-phenyl-C61-butyric (PCBM) explored as active layers in bulk organic solar cells. 10−16 Optimization of their effectiveness has become an active branch of research and can be carried out using various physical methods, 17−21 e.g., by thermal or vapor annealing, 19 by dopant addition, 18,22 or by introducing auxiliary layers to photovoltaic structures. 23 To properly design organic optoelectronic devices, it is necessary to learn the thermo-optical properties of the thin polymeric layers used in their construction. These materials must ensure adequate time stability of many important parameters in addition to the specific physical properties characteristic of the planned device. In particular, it is extremely important to determine the optimal operating temperature range, including peak operating temperatures. Therefore, the glass transition temperature (T g ) 24−28 of polymeric materials and their thin films used in organic electronics is one of the most important physical parameters routinely taken into account at the device planning stage. Its value determines the stability limits of the material microstructure and the stability of material parameters such as thermal expansion or stiffness modulus. For these reasons, the glass transition is the subject of intensive studies in the material science, chemistry, and physics of condensed matter. The results of these studies have been presented in many scientific articles, including both theoretical and experimental studies, see refs 29−39 and the references therein.
Among the various experimental methods suitable for studying the thermo-optical properties of thin polymer films, temperature-dependent spectroscopic ellipsometry plays a special role. 40 A recent review by Erber et al. 40 concerning the determination of the glass transition of polymers in nanoscopic films takes into account the role of the substrate, interfaces, and their dimensional confinement. The authors discussed the influence of polymer film thickness, i.e., the confinement effect, on T g . They pointed out the unique role of temperature-dependent ellipsometry as a nondestructive and very sensitive optical technique. The use of spectroscopic ellipsometry as a highly sensitive and noninvasive method to obtain fundamental information about conjugated polymer films has been overviewed by Campoy-Quiles et al. 41 They discussed the wide possibilities of using in situ ellipsometry to characterize phase behavior in quasi-isothermal experiments. Their work also includes a brief discussion of information that can be derived directly from the raw ellipsometric angles. This discussion is based on the empirical observation that the ellipsometric parameters usually change approximately linearly with the thickness of the layer.
However, issues related to the physical origin of the dependence of ellipsometric angles on temperature have not been adequately discussed in the literature so far. This review article is attempting to fill this gap. In addition, this article discusses the use of temperature-dependent spectroscopic ellipsometry for studying phase transitions in thin polymer films. The benefits of studying thermal transitions using different cooling/heating speeds and a modulated-temperature technique 35 are also discussed. Furthermore, it shows how the analysis and modeling of raw ellipsometric data can be used to determine the thermal properties of thin polymer films. This work also discusses the latest, interesting results regarding thermo-optical properties of thin organic layers irreducibly adsorbed on the substrate. 42

STRUCTURAL ORDER WITHIN POLYMER FILMS AND PHASE TRANSITIONS OCCURRING IN THEM
It is not an exaggeration to say that we live in an era of plastics. Plastics successfully replace metal, and their important applications in various industries are so numerous that it is difficult to count them. Their success results mainly from their ease of production and processing, as well as from the diversity of their properties. All plastics are polymeric materials, but not vice versa. Their properties are widely discussed in the scientific literature, e.g. 45−49 These properties, which can be examined ellipsometrically, will be briefly discussed here. Our description mainly concerns the features of thin polymer films associated with their microstructure and morphology. In general, long-range three-dimensional order at the atomic scale occurs in the crystalline state of matter. On the other hand, the lack of such a structural order is characteristic of the amorphous phase. 50−53 In nature, the crystalline state is more energetically favorable because the energy of an ordered system of atoms or molecules is lower than the energy of an unstructured one.
Usually, the polymers contain crystalline and amorphous regions. A helpful parameter in determining whether a polymer is partially crystalline is the degree of crystallinity. The two-phase microstructure of a polymer can be related to the degree of crystallinity, C, as follows: 54 where ρ c is the density of completely crystalline polymer, ρ s is the density of the sample under consideration, and ρ a is the density of completely amorphous polymer. Crystalline phase formation depends on the architecture of the polymer chain, the polymerization method, and the method of crystallization (which is related to the temperature and the rate of cooling). An example of an almost completely amorphous polymer for which crystallinity degree is approximately zero is atactic polypropylene. Examples of partially crystalline polymers include poly(ethylene terephthalate) (PET) and polyethylene (PE). PET has a degree of crystallinity of 30−40%, and low-density polyethylene (LDPE) has C in the range of 45−55%, while high-density polyethylene (HDPE) has C possible in the range 70−80%. 61 Most commercial polymeric materials are partially crystalline polymers with different degrees of crystallinity. 59 In amorphous polymers, it is only possible to distinguish the local order within a diameter of about 1 nm. The crystalline regions may be as small as about 2 nm in one or two crystallographic directions, and usually less than 50 nm in the remaining. 62,63 The crystallization kinetics of polymeric materials is frequently described by the Avrami equation, 64−67 which represents changes in the volume fraction of the crystallized phase as a function of time.
The thermal transition with the characteristic features of a second-order phase transition (continuous, not stepwise change according to Landau's theory) that can be observed in amorphous polymers, and therefore also in partially crystalline polymers, is the glass transition in which a polymer melt changes on cooling to a polymer glass or a polymer glass changes on heating to a polymer melt. 68 The glass transition is not an equilibrium transformation, unlike other phase transitions. 68−70 This means that during cooling, the supercooled polymer melt undergoes kinetic slowdown rather than the rapid loss of entropy. While cooling, at the glass transition the viscosity of the material increases rapidly and is accompanied by a change in the heat capacity, a decrease in the coefficient of expansion and a significant increase in the relaxation time, among others. Fakhraai and Forrest 71 probed slow dynamics in supported thin polymer films using variable cooling rate ellipsometric measurements of polystyrene. By relating the cooling rate to the relaxation time, they showed that the relaxation dynamics of thin films appear to follow the Arrhenius equation within an activation energy that decreases with decreasing film thicknesses.
Thermal transitions occurring in semicrystalline polymers can be described by appropriate characteristic temperatures on a typical DSC curve, as shown in Figure 1. The first of these is the glass transition occurring at T g . The following shows the cold crystallization at T cc , crystallization at T c , as shown in Figure 2, and melting at T m . At higher temperatures, the processes of oxidation and degradation of polymers are also visible.

THEORETICAL BACKGROUND OF
TEMPERATURE-DEPENDENT ELLIPSOMETRY 3.1. Basics of Ellipsometry. Ellipsometry is an optical technique that enables accurate, indirect optical measurements of film thickness and the dielectric properties of materials. Essentially, an ellipsometer is an optical device that measures the ratio of the intensity of orthogonally polarized light components reaching its detector. The very high sensitivity of this optical technique allows a researcher to follow changes in the polarization of light caused by its interaction with matter.
The basics of ellipsometry are widely described in many books 72−79 and scientific articles. 80−83 For the purpose of this Review Article, a brief introduction on reflection ellipsometry is presented here. Ellipsometry is based on an optical phenomenon in which linearly polarized light is reflected (or transmitted) at the boundary of two media, changing the polarization to elliptical. In the case of optically isotropic materials, this property of reflected light can be described in the following matrix equation: where E ν a denotes the complex amplitude of the electric field component of the incident light beam, a = i, or the reflected light beam, a = r, which in addition can be polarized in two orthogonal directions, either ν = p or ν = s; = ν ν ν r E E / r i def are polarization-dependent Fresnel reflection coefficients. The diagonal form of the reflection matrix in eq 2 indicates the absence of mixing p and s polarized components of light waves upon reflection in the group of materials considered. Standard ellipsometry measures the complex reflectance ratio ρ: 81 where Ψ and Δ are ellipsometric angles. This pair of angles fully describes the polarization ellipse, in a plane perpendicular to the direction of light propagation, and through eq 3 their relationship with the measured quantity ρ. In turn, this ellipse can equally well be described by another pair of angles, namely, ellipticity and azimuth, which better reflect the name of this experimental technique as ellipsometry. For the simplest optical system, consisting of two isotropic infinite media, designated by indexes 0 and 1, with one boundary plane in which reflection and transmission of the incident light beam occurs, which is shown in Figure 3, the explicit form of reflection coefficients r p = r 01 p and r s = r 01 s is as follows: 79 Here, ñ0 and ñ1 are the complex refractive indices for these two media, and corresponding angles of the light propagation θ 0 and θ 1 within the specific medium are related to each other by Snell's law: It follows thus from eqs 3, 4, and 5 that the experimentally measured complex quantity ρ for the two-phase optical system is a function of n0, n1, and θ 0 . Of course, this list should be supplemented with a wavelength of light λ, which has not yet been explicitly considered, because it manifests itself implicitly through the spectral dependence of the relevant refractive indexes. For completeness, it is worth adding that in the complex refractive index̃= + n n ik def the real part is the refractive index n, which describes the change in the speed of light propagation in a given medium, and its imaginary part is   the extinction coefficient k describing the absorption properties of this medium; for a transparent medium, k = 0. In the real world, it often happens that the boundary between two media is in the form of a thin layer with different optical properties. This can be caused, for example, inevitably by the presence of surface roughness or intentionally by deposition of a thin film on the substrate. In such cases, the three-phase optical system, sketched in Figure 3, is a much more correct description. Due to the fact that ellipsometry is an experimental technique showing high sensitivity to the thickness of such a layer, we will also present formulas for the complex reflection coefficients r p and r s for a three-phase optical system. To this end, it can be seen that such an optical system consists of two parallel boundary planes spaced apart by the layer thickness. Therefore, the total amplitude of the reflected wave, in addition to the contribution from the reflection of the incident beam from the upper surface, also includes the contribution derived from all secondary reflections in the layer. Each time the secondary light wave travels the distance between these planes, it obtains an additional phase shift. All this can be expressed mathematically as the following infinite geometric series. 80 Here, t 01 ν and t 10 ν are the Fresnel transmission coefficients of the secondary waves entering and leaving the layer while the change in phase b depending on the wavelength λ, the thickness of the layer h, the refractive index of the layer n1, and cos(θ 1 ) is as follows: Finally, the complex reflection coefficients r p and r s for a three-phase optical system are calculated on the basis of eq 2 for E ν r with ν = p, s using eqs 6 and 7. It is worth noting that these equations depend on the product hn1, by the quantity b. As such, this relationship causes difficulties in accurately determining the thickness of the layer for thin layers, for example, thinner than 20 nm; see, for example, ref 81 and references therein. On the other hand, the use of variable angle spectroscopic ellipsometry 79 can significantly increase the accuracy of determining both the refractive index and the thickness of the layer. Of course, the variable angle in this technique means changing the angle, θ 0 , of incidence of the light beam. The hierarchy of optical models in terms of their complexity extends from the simplest two-phase models through three-phase models to multiphase or multilayer models. We postpone the discussion of the basic methods for determining the thickness of the layer and/or its refractive index using ellipsometry to Section 5. Nevertheless, it should be mentioned here that, in ellipsometry, often such quantities as complex refractive index and complex dielectric function, ε, are used interchangeably because of their mutual relationship ε̃= ñ2. For optically isotropic materials, these dielectric functions have no directional dependence. Uniaxial anisotropic materials with their optical axis oriented perpendicular to the surface, such as, e.g., spin-coated PEDOT:PSS or ordered P3HT phase, are still described by eq 2. 2,82 However, their principal dielectric functions are no longer scalar values, but ε x = ε y ≠ ε z , where ε τ , τ = x, y, z denote τ-component of the dielectric function of a given material. In the case of materials with uniaxial optical anisotropy whose optical axis is oriented in the z direction, their component of the dielectric function in the plane, ε ε ε = = x y o def , is described as ordinary and that off-plane, ε ε = z e def , respectively, as extraordinary. At this stage, the question of why structural phase transitions can be investigated by means of spectroscopic ellipsometry should be asked. To answer this question, the dependence of the refractive index on temperature should be considered.
3.2. Temperature-Dependent Spectroscopic Ellipsometry. Adding the possibility of heating or cooling the sample in a controlled manner during ellipsometric measurements enables thermal analysis of the tested sample. This opens a method to study the dielectric properties of thin polymer films during glass transition or other structural phase changes. Transitions, such as crystallization or melting, as well as the entire spectrum of thermal properties of different materials, can be examined by temperature-dependent ellipsometry, just like in thermal analysis. Optical properties and film thickness are influenced by the thermal expansion. The thickness of the film usually increases, and the refractive index decreases with a temperature increase; see, e.g., refs 83−95. This can be seen by writing an expression for the thermally induced increment dh of the film thickness h upon a temperature increase by dT: where α h is the linear thermal expansion coefficient (TEC), explicitly defined as The expression for an increase in the refractive index, dn, can be obtained from the Lorentz−Lorenz equation, 96−98 which relates n with the mass density, δ, molecular weight, M, and molecular polarizability, γ, which in turn also depends on δ, by the following relationship: 10) where N A is the Avogadro number. Therefore, three different thermal effects δ(T), γ(δ(T)), and γ(T) can contribute to the dependence of the refractive index, n, on the temperature. So, the thermo-optic coefficient (TOC), β, defined as can be written as 99 is the polarity coefficient, α V is the thermal expansion coefficient (causing a change in the mass density), and the subscript 0 indicates that the values of the respective quantities are set at a certain reference temperature, T 0 . It should be noted that it is the experimental study of the linear relationship β(α V ) that can be used to determine the values of both quantities Λ 0 and δ ∂ ∂ ( ) n T . Namely, (i) the presence of the temperature-dependent density δ T ( ) in eq 10 leads to the appearance of α V in the first term of eq 12, which corresponds to the number 1 in the first parentheses; (ii) the nonzero value of γ(δ(T)) causes a deviation from 1 by Λ 0 , which is proportional to α V ; and (iii) the direct dependence of γ on T causes a β shift which is independent of α V . We will also notice that the minus sign preceding the first term in eq 12 takes into account the fact that, as the temperature rises, the density decreases. The coefficients of thermal expansion in eqs 8 and 12 are not the same and are, therefore, marked differently. In eq 8, α h is the linear coefficient of thermal expansion, whereas in eq 12, α V is the volume coefficient of thermal expansion. 100 For isotropic bulk materials, a simple relationship, α h ≈ (1/3)α V , is fulfilled. 101 However, for thin polymer films, volume thermal expansion may be highly anisotropic and dependent on the polymer−substrate interaction character. 99,102 Wide areas of polymer applications have been stimulating intensive research into polymer thermo-optical properties for many years. The thermo-optic coefficients of 10 polymers with diverse chemical structures were analyzed by Zhang et al. 99 The authors found a linear relationship between TOCs and corresponding TECs, with a directional coefficient of −0.56 and an intercept δ ∂ ∂ ( ) n T of −3.7 × 10 −6 which directly relates to the linear dependence β(α V ) described by eq 12. We show this relationship in Figure 4, because it allows easy estimation of the value of one of the two parameters involved, as well as comparison of the determined values with the published data.
Ellipsometry is sensitive both to the thickness of the film and its refractive index, i.e., ρ(n, h). 103 Therefore, four thermal effects contribute to the dependence of ρ(Ψ, Δ) on temperature. Three of them, mentioned above, contribute to the dependence of the refractive index n on temperature, which within the linear approximation and definition β given by eq 11 can be written as follows: However, the fourth effect results from the dependence of the film thickness, h, on the temperature. Let us explain in more detail the origin of the fourth effect. Basically, the coefficient of thermal expansion is an anisotropic quantity. Given a thin layer with uniaxial symmetry, its average volume TEC can be written as α V ≈ α x + α y + α z , where α i , i = x, y, and z, are linear TEC in x, y, and z directions, and α x = α y ≠ α z , respectively. Now it is easy to see, from eqs 12 and 13, that the thermo-optical properties (n(T), β 0 ) of the layer are affected by α V with a nontrivial volumetric effect. However, h(T) will still change only under the influence of α z . Campoy-Quiles et al. 41 gave a brief overview of information that can be obtained directly from the raw ellipsometric angles. Their discussion is based on the empirical observation that the ellipsometric parameters usually change approximately linearly with the thickness of the layer. However, on the basis of the above discussion, we can determine the range of its application, namely, to isotropic films, i.e., those for which the following equality holds: Summarizing the information provided here, we will list the quantities necessary to determine the thermo-optical properties of thin polymer films, which can be determined using temperature-dependent ellipsometry, namely, n 0 , β 0 , and α h . Of course, the use of spectroscopic ellipsometry allows us to determine the spectral dispersion n 0 and β 0 . However, if we additionally know α V , we can also determine Λ 0 , taking into for polymers is a small quantity. 99 It should be remembered that the assumption α V ≈ 3α h can often be made, especially to estimate the value of Λ 0 for relatively thick films.
According to the literature discussed below, it is often possible to quite accurately determine the characteristic temperature (temperatures) of the thermal transition based on the analysis of raw ellipsometric data. The relevant examples are presented and discussed in Section 4. However, studies of thermal transitions are most often based on determining the appropriate changes in physical quantities, such as thickness h or refractive index n, using ellipsometric data modeling. The issues related to the modeling of ellipsometric data are discussed in more detail in Section 5.

DETERMINING THE TEMPERATURE OF THERMAL
TRANSITIONS BASED ON ANALYSIS OF RAW ELLIPSOMETRIC DATA The issues related to determining glass transition temperatures in thin polymer films based on the analysis of raw ellipsometric data will be discussed first. The observed size of changes in ellipsometric angles caused by a 100°C sample temperature change does not exceed a few degrees at best. This is directly related to the low numerical value of material parameters such as TEC and TOC. Representative data for polymers are shown in Figure 4. However, one can expect that the linear approximation to the ellipsometric angles, discussed in the previous section, should be sufficiently well met in a wide temperature range on both sides of the thermal transition. This observation also indicates a way to detect thermal transitions. Namely, straight lines should be drawn in the low-and high-temperature ranges for the considered ellipsometric angle, and the point of their intersection will determine the temperature of a given thermal transition.
Clough et al. 110 provide an example in which the glass transition temperature was determined based on the analysis of raw ellipsometric data. They studied thin polystyrene films coated on SiO x /Si(100) substrates and performed temperature scans using a single-wavelength Stokes ellipsometer. They annealed the samples at 125°C for 1 h before measurement. The measurements were taken during the cooling cycle in the temperature range 140−36°C at 1°C/ min. In Figure 5, the results for a 30 nm thick PS layer are shown. In Figure 5a, the dependence of Ψ vs temperature is presented. The good quality of the experimental data is clearly visible; however, the contrast between the lower and higher border temperatures is not very high in this thin film. Therefore, in order to increase the sensitivity, the authors determined T g using the temperature dependence of the Ψ derivative as shown in Figure 5b. In this approach, T g is defined as the midpoint of the higher and lower border temperatures.
One of the oldest examples of the use of one-wavelength ellipsometry is the work of Beaucage et al., 91 in which thin layers of polystyrene were studied. This paper presents the dependence of the refractive index as a function of temperature inverse, on the basis of which the glass transition temperature of the tested material was determined.
Another example in which the glass transition temperature was precisely determined on the basis of the analysis of raw ellipsometric data is the work of Fakhraai and Forrest. 71 They studied thin polystyrene (PS) films spin-coated on a Pt(50 nm)/Ti/SiN(50 nm)/Si substrate and performed the temperature scans using an EXACTA 2000 fast nulling ellipsometer. The authors determined T g from the temperature dependence of the analyzer angle derivative at a cooling rate of 6 K/min. In their approach, the T g was defined as the midpoint of the high border temperature T + and the low border temperature T − .
A similar method for determining the glass transition temperature can be found in the work of Wang et al. 111 These authors investigated the effect of film thickness on its molecular structure and T g for poly[N-9′-heptadecan-2,7carbazol-alt-5,5-(4′,7′-di-2-thienyl-2′,1′,3′-benzothia-diazol)] (PCDTBT). Figure 6b, taken from their work, shows Ψ as a function of T at λ = 650 nm for three films that are 102, 62, and 30 nm thick. It is easy to see that their glass transition temperatures, which are 128, 120, and 110°C, respectively, are closely related to the thickness of the samples. So, this example shows that T g decreases as the layers become thinner.
The glass transition temperature can be determined from the temperature dependencies of both ellipsometric angles, Ψ and Δ. Ouakili et al. 113 describe the multiple glass transition phenomena of polystyrene and poly(methyl methacrylate) thin films with thicknesses in range 5−200 nm. Figure 7 presents an example of the dependence Δ(T) for two thin films of polystyrene with a thickness of 9 and 70 nm. The glass transition T g of these samples is placed at about 80 and 95°C, respectively.
Similar curves for thin PS films are presented in work of Chandran and Basu. 114 They studied the variation of thin films' glass transition temperatures as well as interface and surface effects for polymer−nanoparticle blend systems by  Raw ellipsometric data often cannot be fitted well linearly. Sometimes the calculation of the derivative dΔ/dT allows the characteristic temperatures to be determined with greater accuracy. For example, this is clearly seen in Figure 8a, taken from work of Geng and Tsui, 115 where the temperature dependence of Δ(T) is of this nature, while the analysis of dΔ(T)/dT allows a more accurate T g determination. These tests have been carried out on PMMA films with the same molecular mass deposited on silicon substrates for multiple film thicknesses.
It should be emphasized also that Δ(T), in Figure 8a, changes by a few degrees as T increases, but Ψ changes by only a few tenths of a degree, so the noise level should be taken into account. Therefore, the wavelength and angle of incidence at which the best sensitivity can be achieved should also be considered. A good example of this can be found in the work of Lee et al., 142 from which we present Figure 9, showing the spectral dependences of Ψ(λ) and Δ(λ) at 30, 80, 130, and 170°C, for 110 nm thick PS film deposited on the grafted PS layer. In order to determine the glass transition temperature, two spectral widows have been chosen in this example, namely, 553−563 nm (A) and 734−744 nm (B).
These spectral windows were selected to provide good sensitivity to temperature changes at lower wavelengths and to avoid interference caused by noise at higher wavelengths.
The other characteristic temperatures of thermodynamic phase transitions occurring in thin polymer films were also determined on the basis of the analysis of raw ellipsometric data. The phenomena observed were related to the cold crystallization, melting on the heating cycle, or crystallization on the cooling cycle, in accordance with the thermal transitions sketched in Figure 2. They are usually accompanied by rapid changes in the value of ellipsometric angles as opposed to the glass transition case discussed above. 88−117 In Figure 10, the temperature dependence of the ellipsometric angles Ψ(T) and Δ(T), as determined by Campoy-Quiles et al., 41 illustrate these relationships. The authors studied samples of thin films of regioregular (with side chains in head−tail order) P3HT. The kinks observed in Figure 10a,b just below 100°C during the heating cycle are due to the glass transition of the amorphous phase. The significant changes around 200°C correspond to the melting of the     Figure  11a, where the temperatures of these thermal transitions were marked. In turn, Figure 11b shows how these transitions are affected by film thickness. The thickness of the samples was in the range 17−59 nm, and both transitions were visible only for films whose thickness was greater than 40 nm. The dipper insight into the origin of these thermal transitions can be provided by ellipsometric modeling; see the discussion at the beginning of Section 5.4. We will only mention here that the interpretation regarding the effect of layer thickness on the characteristic temperatures of thermal transitions is confirmed by AFM. 118 It is also worth noting that, as in the case of DSC, subsequent heating and cooling cycles can be used in temperature-dependent ellipsometry to monitor structural changes in the polymer layer under investigation.
Up to now, in this section, we have focused on issues related to determining the glass transition temperatures of thin polymer films based on the analysis of raw ellipsometric data, presenting representative examples taken from a rich collection of published works. We will now discuss cases in which such analysis provides equivocal results. The problem is the accuracy of determining the temperature of a given thermal transition, which can be affected by several factors, such as those related to the quality of the sample, i.e., surface roughness, and thickness unevenness, 40,80 as well as the presence of residual solvent in the samples, 119−121 error measurement, noise, a large dispersion of data, and a small range of linear dependence of data values. Problems related to sample imperfections can and should be identified by measuring the degree of depolarization. Also, there are methods, such as optical models, that can take into account the physical imperfection of samples. 80 This always happens at the expense of using more complex optical models containing a greater number of specific parameters. Therefore, the best practice is accurate sample selection and the use of only nondepolarizing samples in tests by temperature-dependent spectroscopic ellipsometry.
Experimental data may also be dominated by the thermal properties of the substrates on which the samples are deposited 122 and/or may depend on sample/substrate interaction. 102,123,124 For this reason, detailed knowledge about the temperature-dependent dielectric properties of the substrate is required, and appropriate substrate selection can greatly facilitate the analysis of ellipsometric data. A good example illustrating the influence of the substrate on ellipsometric results is provided by Glor and Fakhraai. 148 Separate physical processes that can affect the measurement data can be related to the sorption or desorption of atmospheric gases on the surface of the sample. In addition, water vapor may also condense on the surface of the sample. Efremov et al. 125 illustrate well how a vacuum and the presence of residual gases affect the results of temperaturedependent ellipsometric measurements. Figure 12    , and this indicates that the expansion coefficient has a negative value in this temperature range. The presence of residual gases, whose main ingredient is water, appearing due to desorption from the surface of the sample, causes the appearance of negative thermal expansion coefficient values.
We finish this chapter with a brief overview of information that can be obtained directly from raw ellipsometric data in temperature-dependent ellipsometry. Namely, studies of the literature show that it is common practice to determine the glass transition temperature based on raw ellipsometric data, e.g., refs 99142. Representative examples from the literature on the use of raw ellipsometric data obtained from temperature-dependent ellipsometry are shown in Table 1. Some of these studies are limited to this method of determining the temperature, T g . It should be remembered, however, that raw temperature-dependent ellipsometric data reflect changes in both layer thickness and refractive index caused by temperature changes. Because these contributions can be competitive, thermal transitions in raw data may also be poorly visible in such cases. For this reason, there is a large group of works in which, apart from analyzing raw ellipsometric data, appropriate ellipsometric modeling is also performed.

TEMPERATURE-DEPENDENT ELLIPSOMETRIC
MODELING In the previous section, the determination of thermal transition temperatures based on analysis of raw ellipsometric data was discussed. Issues regarding the determination of characteristic temperatures using ellipsometric data modeling based on more advanced numerical techniques will be discussed in more detail here. Let us briefly discuss the basic approaches used to determine optical properties by ellipsometry. The simplest one, called the direct inversion method, is just the mathematical operation performed in order to calculate a complex refractive index from measured ellipsometric angles Ψ and Δ. For example, for the simplest optical model that consists of two phases, the refractive index n1 of medium 1 can be calculated using eqs 3, 4, and 5 as follows: 80 Therefore, eq 14 can be used to calculate n1(λ), point by point, from the known dependence of n0, Ψ, and Δ on the wavelength, λ. It should be emphasized that direct inversion does not introduce any parametrization on the spectral dependence of n1. In contrast, more advanced approaches use such a parametrization. They are intended for a specific group of materials and for a specific range of the light spectrum. 72,79,80 A useful and most commonly used parametrization is known as the Cauchy optical model. 72−80 The Cauchy dispersion describes the dependence of the refractive index and extinction coefficient on the wavelength according to the following expressions: 132 Here, n(λ, T) is the refractive index, and k(λ, T) is the extinction coefficient. Temperature-dependent parameters n i and k i (where i = 0, 1 and 2) are the Cauchy model parameters. A clear advantage of this optical model is that for transparent materials the value of the extinction coefficient can be set to zero, i.e., k(λ, T) = 0, in their spectral transparency range. Moreover, this model also provides a simple description of materials with weak spectral dependence at long wavelengths. 72,79,80 Still, a simple optical system consisting of a thin polymer film deposited on a thick substrate, as shown in Figure 3, contains three optical phases: ambient medium, film, and substrate. For this system, ρ = ρ(n0, n1, h, n2, θ, λ) is a function of physical parameters such as the thickness of the polymer film h, the wavelength λ of the incident light at an angle θ, and complex refractive indexes, which in turn may also depend on λ and additional parameters describing their dispersion. In principle, the unknown parameters are determined by fitting the modeled ρ m (Ψ m , Δ m ) to the ellipsometric data ρ(Ψ, Δ) using the following relationships: where the dependence on θ, λ, and T is not explicitly written. Most often in the experimental implementation, ñ0 is 1, because it is atmospheric air, and the method of calculating ñ2 has been described above, so there are only two quantities in eq 17 to determine, namely, ñ1 and h. It still can be considered as a simple case when one of these quantities is known. Because then the second can be accurately determined by adjusting the parameters describing the optical model of the sample to the ellipsometric data. The cases where both the thickness of the film and its refractive index are unknown are more complex. However, it often happens that polymer films are transparent in a fairly wide spectral region. In such cases, it is best to place the film on the light-absorbing substrate and solve the ellipsometric equation for a three-phase system. In The Journal of Physical Chemistry B pubs.acs.org/JPCB Review Article the considered cases, the refractive index and film thickness can still be relatively easily determined, using the combined Cauchy model in the transparent area, followed by direct inversion for the entire spectrum; for more details, see, e.g., ref 80 and references therein. However, in practice it looks a bit different; namely, the commercial software provided with the ellipsometer uses robust algorithms based on linear regression, 72 which are much more universal and can also be used in much more complex optical systems, e.g., for multilayer systems. Ellipsometric modeling is widely described in the literature, including textbooks, e.g., 72,79 so we will shorten the discussion on this topic to the necessary minimum. Analysis of ellipsometry data is usually performed using linear regression analysis, and the optical constants and film structures, including layer thickness, but not only this, are determined by minimizing fit errors calculated on the basis of the mean squared error function (MSE). 72 During modeling, ellipsometric angles Ψ m and Δ m are calculated for specific values of their parameters and compared with corresponding Ψ and Δ determined experimentally. All this is done as part of one error function that can simultaneously contain data for different incidence angles and experiments with multiple samples, e.g., samples with different layer thicknesses. This procedure is repeated iteratively for changed parameter values in such a way as to minimize the error function. 72 At the end of this discussion, we emphasize that, in ellipsometric modeling, verification of results, including analysis of statistical errors and correlations between parameters, is a very important step. A solid theoretical basis for determining temperaturedependent optical properties is provided by the linear n(T) analysis. However, this is despite the fact that the unified linear analysis of reflection and transmission ellipsometry was developed by Azzam et al. 143 already in 1975. The corresponding extension of this linear analysis to temperature-dependent ellipsometry has not yet been reported. Therefore, the techniques discussed so far are used, explicitly taking into account the dependence of the parameters of the optical model and the thickness of the layer on the temperature. However, some tips on temperature-dependent ellipsometry can be given. Namely, an important requirement for reliable modeling of ellipsometric data is that the number of unknown parameters of the model should not be greater than the number of equations used to determine their   values. 10,104,105 For this reason, ellipsometric studies are often carried out in the spectrally transparent range of the sample being tested, using the direct inversion method or the Cauchy optical model. 107 Another important requirement is the validation of the optical model, e.g., by its use to describe the dielectric properties of the tested material at a reference temperature, usually room temperature. In addition, for supported films, knowledge about temperature-dependent dielectric properties of the substrate is also required. 148 Usually, the assumption that the refractive index of ambient atmosphere is temperature-independent is sufficiently good (the refractive index of air is approximately constant and near 1). More practical, interesting, and very valuable information on the modeling of temperature-dependent ellipsometric data is provided, for example, by Glor and Fakhraai in ref 148. Also note that when discussing the published results below, we will denote TEC and layer thickness using the symbols from the original articles. 5.1. Determination of Characteristic Temperatures of Thermal Transitions Using Data Modeling. The direct inversion method was used, for example, by Lee et al. 142 These authors studied thin layers of PS terminated with hydroxyl groups (PSOH) and PS. Ellipsometric angles Ψ and Δ were measured at an incidence angle of 70°in a very narrow 10 nm spectral range (553−563 nm and 734−744 nm). Figure 13 shows the temperature dependence Ψ, Δ, and thickness, and the appropriate glass transition temperature is marked on all graphs. In addition, the coefficient of thermal expansion of the material before and after the vitreous transition was also marked. 142 More reliable results of ellipsometric modeling can be obtained if a wider spectrum of wavelengths is used. This can be achieved by parametrizing the spectral dependence of the complex refractive index. This approach has been successfully used in temperature-dependent spectroscopic ellipsometry to determine h(T) and n(λ, T), for example, in the work of Hajduk et al., 132 where phase transitions occurring in thin films of P3HT, PCBM, and their blends were studied. Figure 14 shows, for example, the dependence of film thickness on temperature for thin layers of these materials.
The characteristic temperatures visible on these diagrams for pure materials are attributed to glass transitions at 10°C for P3HT and 110°C for PCBM, respectively, while cold crystallization occurs at higher temperatures. However, in the case of P3HT:PCBM, Figure 14c, the presence of such a large number of thermal transitions indicates the phase separation in this material; for more details, see 132.
A similar way of using the Cauchy model can be found in Park et al. 145 The authors investigated the dependence of T g on the thickness and composition of thin poly(costyrenemethacrylate) P(S-co-PMMA) and poly(2-vinylpyridine-costyrene) layers deposited on Si substrates. The glass transition temperature was determined as the abscissa of the intersection point of two linear approximations to the curve on both sides of a given thermal transition. Another example of the use of ellipsometric modeling can be found in the work of Sharp and Forrest,146 who studied the segmental dynamics in thin layers of isotactic poly(methyl methacrylate). More examples of determining T g based on h(T) and n(T) curves from ellipsometric modeling can be found in. 30,83,84,86,88,93,95 In work of E. Bittrich et al. 92 thin films of polyimide were studied. They characterized glass transition of these films and analyzed their thermal expansion. Representative results from their work are shown in Figure 15. Although the glass transition temperature has not been determined directly in the plots, it can be seen that the thermal transition is influenced by the film thickness.
At this point, it should be said that the refractive index of the polymer layers can sometimes increase, and their thickness decreases with increasing temperature. An example can be found, e.g., in the work of Jaglarz et al. 147 in which thin layers of poly(3-hexylthiophene) P3HT and poly(3-octylthiophene) P3OT were tested. The above behavior was observed at elevated temperatures exceeding 200°C, and the authors attributed this to structural changes in the films. 147 The effect of film shrinkage during their heating can also be caused by the loss of residual solvent and is known in the literature. 121,182 A good example of this can be found in the article by Baker et al., 183 presenting methods of studying structural relaxation in polymer films using ellipsometry.
The necessity of standardizing the method of determining the glass transition temperature based on ellipsometric measurements was indicated by Glor and Fakhraai. 148 They proposed an appropriate protocol using the variable cooling The Journal of Physical Chemistry B pubs.acs.org/JPCB Review Article rate technique. This proposal was based on methodology already explored, inter alia, to study thin polystyrene films in order to show the correlation between the average dynamics of the film and the dynamics of its free surface. 149−154 The proposed protocol consists of several steps, starting from rigorous recipes for preparing thin polymer films and determining their thickness. Then, the glass transition temperature can be determined, depending on the cooling rate. The authors indicated the necessary elements of the analysis of the average dynamics of the studied layers. The cooling-rate-dependent T g measurement (CR-T g ) protocol was included in Glor et al. 155 The authors examined thin films of poly(2-vinylpyridine) (P2VP) (with thickness 217 nm), deposited on silicon substrates. They recorded thickness changes during temperature modulation. Relevant plots of the dependence of thickness versus time are shown in Figure  16. 155 This example also clearly demonstrates how important it is to accurately account for the thermal expansion coefficient of the substrate. In many works aimed at determining T g , the behavior of the thermal expansion coefficient α h was analyzed as a function of temperature in accordance with eq 9; see, e.g., refs 36, 156−159. However, in practice, instead of eq 9, finite element approximation is used to numerically differentiate layer thickness as a function of temperature. The explicit form for α h (T) in this approximation is as follows: where ΔT is the temperature step. Kawana (19) where w is the temperature width and G, A, M, and B are the parameters describing linear behavior in the vitreous and rubbery states. The changes in α h (T) according to eq 19 are shown in Figure 17. Quite recently, Erber at al. 40 used the analytically integrated form of eq 19 to describe the dependence of the layer thickness on temperature: We conclude this section with a brief overview regarding the application of ellipsometric modeling to data obtained from temperature-dependent ellipsometry ( Table 2). In the ellipsometric modeling, independent of the temperature, the measured ellipsometric data is used to determine the values of physical quantities, such as h, n, and k. For optical systems with a large number of parameters, the dependence of ellipsometric data on the angle of incidence is also used. However, in temperature-dependent ellipsometry, the amount of recorded data increases so much that these measurements are made practically only for one fixed value of the angle of incidence of light. Therefore, in this experimental method two physical quantities, usually h(T) and n(T), are determined on the basis of measured ρ(Ψ, Δ, T). Among the various experimental techniques used to determine these quantities, temperature-modulated ellipsometry also gives the possibility to distinguish thermal processes that are reversible from those that are irreversible. 148

Influence of Thin Polymer Film Thickness on T g .
The reduction of the glass transition temperature when the thickness of the film decreases is reported in many scientific articles; see, e.g., the review by Erber et al. 40

and ref 42.
However, there is evidence that very thin polymer films can also have an increased T g value than the corresponding bulk materials due to interaction with substrate. 42,102 Over 20 years of research on the effects of geometric confinement of polymeric films did not lead to unequivocal results. 40,42,91,92,157 There are a large number of scientific works reporting that the glass transition temperature of thin a Poly(n-butyl methacrylate). b Styrene/n-butyl methacrylate copolymers.
The Journal of Physical Chemistry B pubs.acs.org/JPCB Review Article polymer films differs significantly from the bulk value, T g,bulk . 29,34,35,[42][43][44][84][85][86][87]126,137,158,168,169 However, the interpretation of these results is not unambiguous and is often controversial. 91,92 Undoubtedly, it is widely accepted that surface effects play a decisive role here. Namely, the physical conditions at the interface of the polymer film with the air and/or substrate can be significantly different from those at the center of the film. Striking examples are works devoted to irreversibly adsorbed ultrathin polymer layers, whose thermooptical properties are precisely the result of a strong interaction with the substrate. 42−44 Therefore, taking into account that a glassy transition occurs at temperatures where independent movement of polymer chain segments with lengths of just a few mer units becomes possible. 184 It is clear that the controversies mentioned above relate to the thickness of the polymer film, for which the influence of surface layers begins to be significant. Erber et al. 40 noticed that most robust results suggest that T g deviations occur below 15 nm film thicknesses and can be assigned to interfacial interactions. Figure 18 presents results from refs 40, 71, and 116, showing the T g deviation from the bulk value, T g,bulk , for thin PS films, deposited on different substrates. In this figure the horizontal line at ΔT = −5°distinguishes somehow arbitrary data with a visible deviation from T g,bulk . Moreover, the substrate type had no effect on the glass transition temperature for layers with a thickness above 40 nm, while previous reports significantly moved the limit of this effect, even to 100 nm. 170−172 Kim et al. 87 showed the dependence of the thermal expansion coefficient on temperature for two PS layers with thicknesses of 475 and 23 nm. In the case of a very thin layer, the points are very scattered, which makes it difficult to determine the glass transition temperature. The authors used the fast Fourier transform method over 60 data points to obtain a smooth curve and determine T g more accurately. To explain the effect of scattered points, a three-layer model, first proposed in DeMaggio et al., 172 can be used, where the supported film consists of three layers: the free-surface layer, the middle layer, and the substrate layer (adsorbed layer). Confirmation can be found in Ellison and Torkelson's work; 173 they reported that the value of averaged T g in thin polymer films depends on their thickness, on the basis of fluorescence measurements. On the other hand, the T g value in the surface layer of the polymer decreases with respect to the bulk T g , unlike in the substrate layer. Each component layer consists of cooperatively rearranging regions, in which length is related quantitatively with the length of nanoconfinement (its confinement in nanoscale). The thickness of the layer in which there is a deviation in T g from T g,bulk depends on the interaction strength at the interface. When there is a strong interaction between polymer film and substrate, the thickness of such a layer may exceed 100 nm. Further explanations can be found in Kawana and Jones, 159 who presented an analysis of ellipsometric data collected on thin PS films (with thickness in the 10−170 nm range) deposited on Si substrates. Similarly to ref 87, the authors observed a significant dispersion of points in the α(T) dependence for layers thinner than 30 nm. An example can be found in Figure  19, which shows the TEC of 29 and 164 nm thick polystyrene films as a function of temperature.
This was explained by the presence of a thin surface layer with liquid properties. In the study, 158 the temperature dependence of the TEC α is normalized to the averaged value of the molten (liquid) phase α liq , for films with different thicknesses. To rationalize results, it was assumed that α(T) of Here, the labels "liq" and "glass" refer to the quantities in the liquid and glass phases, respectively. Note that this approximation accounts for the existence of the two separated phases in layers with temperature-dependent thicknesses h liq (T) and h − h liq (T), respectively. Figure 20 presents this thickness as a function of temperature. The 10 nm thick liquid-like layer exists even at temperatures well below glass transition.

Dependence of the Coefficient of Thermal Expansion on the Thickness of Irreversibly Adsorbed
Layers. Spontaneous interactions that can arise between the rigid surface of the substrate and the soft surface of the polymer cause polymer chains to be more adsorbed to the substrate. Such adsorption may be irreversible due to van der Waals forces. One of the key papers describing the impact of interactions between the layer and the substrate on the glass transition temperature was the work of Wallace et al. 102 The tested thin polystyrene layers were deposited on hydrogenated silicon substrates. The thickness of tested films was from 7.5 to 198.8 nm. In this work, the values of the thermal expansion coefficient were measured by means of X-ray reflection. It was shown that thin films less than 40 nm thick did not show a glass transition until the temperature value was at least 60°C higher than the T g of the bulk. The coefficient of thermal expansion of these films was described as follows: 102 where α n glass is the thermal expansion coefficient of the layer constrained within the substrate in the glassy state and α n melt is the thermal expansion coefficient of the molten layer, A is the constant characteristic thickness, and h is the thickness of the layer we are interested in. The value of A indicates approximately the thickness for which the change in glass transition temperature is not noticeable. This is the distance over which the substrate affects the polymer film and its T g .
For films thicker than 40 nm, the chains on the free surface have high mobility, which is due to the fact that the glass transition temperature of the surface layer is lower than that in the bulk. However, this relationship does not work for very thin layers (below 40 nm), when the substrate plays a major role and affects the thermal expansion of the material, regardless of whether the chains on the surface have greater mobility. In such cases, we can talk about the compensation of the effect of the free surface of the layer on the temperature T g through its interaction with the substrate. An interesting example from the literature in which this effect was studied is the work of Perez-de-Eulate et al. 174 describing how irreversible adsorption reduces the effect of the free surface on the T g of thin layers of poly(4-tert-butylstyrene) (PTBS). In particular, the authors determined the dependence of the thickness of an irreversible adsorbed layer on the annealing time of a thin layer of PTBS. Thin PTBS films were deposited by spin coating from polymer solutions on silicon wafers coated with a native oxide layer. The films were deposited at a constant thickness of 200 nm and annealed for various times at 180°C. Chains that were not adsorbed were removed with the same solvent, and the sample after 10 min of exposure under atmospheric conditions was measured by spectroscopic ellipsometry. Ellipsometric angles Ψ and Δ were fitted using Cauchy's dispersion, and the multilayer optical model consisted of the following layers: air/PTBS/SiOx/Si. The relation of Dalnoki -Veres 175 is used for describing the effect of annealing on the thickness dependence of the T g : where ξ 0 is the altered segmental mobility and R s is the rheological temperature at the free interface. This quantity is derived from parametrization of the dynamics and is indicating the temperature at which the molten bulk layer has the same mobility like a free surface and ξ 0 has the length scale associated with perturbation of the glass transition. Figure 21 presents, for PTBS films, the comparison between the time evolution of ΔT g with the recovery of bulk T g upon annealing and with the kinetics of irreversible adsorption. The obtained results suggest that irreversible adsorption at the polymer/substrate interface limits or totally erases the impact of the free surface on the thermal T g . The free surface does not affect T g in very thin layers, because the increased mobility of the chains on this surface is balanced by the reduced mobility of the chains that are irreducibly adsorbed on the substrate. Therefore, such very thin layers change into ultraviscous plates, whose flow properties are smaller compared to those of the bulk material. 176,177 The mechanism of irreversible adsorption was studied, e.g., in work of Davis et al., 178 where the pure polymers and copolymers of PS and PMMA deposited on silicon wafers were investigated. The difference between adsorption of pure polymers and copolymer is shown in Figure 22.
Polystyrene films show looser surface adsorption than PMMA films, while the P(S-r-MMA) copolymer has more substrate-bound chains than PS and does not show large loops. It is worth noting that the copolymers show a dependence of T g on the composition of homopolymers, which allows controlling the glass transition temperature by changing their composition accordingly. It is these irreversible adsorbed layers with heterogeneous composition that behave differently. The relationship between the increase in thickness of such a layer and the annealing time at constant temperature (T g,bulk + 45 K) was measured by ellipsometry. It was shown that the rate of increase in the thickness of the adsorbed film was initially logarithmic (for the first 3 h), and then these changes became almost constant and linear (for the next 32 h). The trend of linear growth rates indicates that thermal chain fluctuations play a smaller role in reducing adsorption than the interaction energy with the substrate, which is associated with so-called packing frustration. A brief overview of the impact of irreversible adsorption on the mechanisms of thermal expansion of confined 1D polymer layers was given in the paper of Braatz et al. 42 In these investigations, thin films of polystyrene were deposited on aluminum and silicon oxide substrates. The thermal expansion coefficient in the direction normal to the substrate was determined using temperature-dependent spectroscopic ellispometry. The ellipsometric angles were fitted to the multilayer model. In particular, thin supported films have been described using a three-layer model to explicitly take into account the presence of a surface layer and one that comes into contact with the substrate. The latter layer is about 4 nm thick and is also known as the dead layer 43,44 in which there is no molecular movement; its coefficient of thermal expansion is approximately 0. The structure of these layers is densely packed with loose adsorbed chains. 179 It has been shown in a study 43 that dead layers have a thermal expansion coefficient greater than zero under the influence of packing frustration. 180,181 The dependence of TEC coefficients (here the thermal expansion of irreversible adsorbed layer in glass and melted state) on the thickness of the adsorbed layer, h ads , is shown in eq 23 below. Thermal expansion of the interfacial layers, in the melt and glassy states, is enhanced with respect to the bulk values. The film thickness has a similar impact on both thermal expansion coefficients, and the final expression for this dependence can be written as follows: 42 This equation is valid for h ads > δ(h ads ), where δ is the thickness of the dead layer. Thorough analysis revealed that excessive thermal expansion occurs near the adsorption layer. This is a characteristic feature of very thin polymer layers, which results from packaging frustration and is independent of the chain length. For the formation, dead layers are responsible the forces immobilizing the chains on the surface, which are also subjected to the packing frustration. It can be concluded that, in fact, the thickness of the dead layer increases for films that are thick enough not to undergo excessive thermal expansion at the joint. In addition, the increase in thermal expansion coefficients TEC indicates that the adsorbed layers are not clearly associated with an increase in T g . 43 5.4. Surface Effects and Thermal Transition Depth Profiles in Thin Polymer Films. We have already discussed in Section 4 that the thermal transitions occurring in partially crystalline polymers are connected with microstructural reorganization of polymeric material and can be described by appropriate characteristic temperatures. Therefore, thermodynamic phase transitions, such as cold crystallization, crystallization, or melting of crystals, can be successfully investigated using the techniques discussed so far. Changes in the raw ellipsometric angles associated with thermodynamic phase transitions are often more pronounced. An example of this is the stepwise crystallization shown in Figure 11, discussed in Section 4, as observed in thin PET films by Xu et al. 118 However, significantly deeper insight into the origin of these thermal transitions and their nature can be obtained using ellipsometric modeling. Thus, in the mentioned example, a multilayer optical model was used, which included the thickness of the layer adsorbed directly on the substrate,   the thickness of the layer in which the interactions of the substrate disappear, followed by the bulk type layer, and finally the thin surface layer in which the increased mobility of polymer chains occurs. The use of this optical model has allowed for a much more comprehensive interpretation of the experimental results and, for example, the identification of a smaller peak as derived from thermally induced structural changes in the surface layer; for more details, see ref 118. It is worth noting here that surface effects regarding the difference in the glass transition temperature between the surface layer and the polystyrene core layer are described in refs 159−161. Very interesting research in this topic is the work on determining the depth profiles of thermal transitions in thin polymer films. Such studies were carried out by Muller et al. 135 They tested four different organic semiconductors (among others, the well-known and widely used P3HT), which are used as active layers in optoelectronic devices. Transition temperature depth profiles calculated by the authors included the following: liquid-crystalline melting temperature ( -T m l c ), crystalline melting temperature (T m ), and recrystallization temperature (T c ). The authors showed on the example of the 70 nm thick APFO3 film that the average temperature of liquid-crystalline melting ( -T m l c ) is dependent on the wavelength of the applied light. More specifically, the estimated average melting temperature at a wavelength of 800 nm was approximately 185°C, while the reading at λ = 550 nm had a value close to 215°C. In order to determine the depth profile of -T m l c , they proposed a model that is a combination of two phases, where separate sets of the optical parameters describe the liquid-crystalline phase and the isotropic liquid phase, respectively. Next, they assumed that below the − T m l c temperature, the liquid crystal layer has no sublayers, and the isotropic liquid is at temperatures greater than − T m l c . In addition, in the case of an intermediate temperature range, the authors examined several models of intermediate film configurations. A two-layer structure consisting of a liquid crystal layer on top of the isotropic liquid layer provides the lowest standard deviation when fitting the model parameters to the ellipsometric spectra over the entire temperature range. 135 An example of thermal transition temperature depth profiles for five thin polymer films as determined by Muller et al. 135 is shown in Figure 23.
As can be seen in this figure, the surface effect is clearly visible for all layers of tested polymer systems. It can be distinguished by rapid changes in transition temperature in a thin surface layer.
The examples discussed in this section show that temperature-dependent spectroscopic ellipsometry also makes it possible to distinguish local thermal transformations (such as glass transition T g , crystallization temperature T c , or melting point T m ) on the surface (T t surface ) from thermal transitions in the bulk of the sample (T t bulk ).

OUTLOOK AND ENDING REMARKS
Temperature-dependent spectroscopic ellipsometry is a nondestructive optical technique suitable for experimental polymer science. Frequently, thermal transitions occurring in thin polymer films are clearly visible, so it is possible to determine characteristic temperatures directly from raw ellipsometric data, while detailed analysis of ellipsometric data can provide accurate quantitative information on the dependence of the thickness of polymer films and their optical parameters on temperature. In this work, the methods used to analyze and model raw ellipsometric data have been discussed. Since ellipsometry is sensitive both to the thickness of the film and its refractive index, four thermal effects contribute to the dependence of ellipsometric angles on temperature. Three processes are derived from the dependence of the refractive index n on temperature. The fourth effect results from the dependence of the film thickness on the temperature. The raw ellipsometric data analysis is used to determine the characteristic temperature(s) of the thermal process(es). This analysis works well under the following conditions: (i) the thermooptical coefficient and thermal expansion coefficient are very weakly dependent on temperature, so the linear relationship for ρ is satisfied; (ii) the data is not noisy and highly dispersed; and (iii) the substrate, on which the film is deposited, has been selected appropriately. Temperaturedependent ellipsometric modeling is more complex and timeconsuming because it requires the use of more advanced numerical techniques to analyze large data sets. On the other hand, the advantage of this approach is the ability to determine the temperature dependence of the thermal expansion coefficient, TEC, and the spectral dependence of the thermo-optic coefficient, TOC. The highest level of modeling for temperature-dependent spectroscopic ellipsom- The Journal of Physical Chemistry B pubs.acs.org/JPCB Review Article etry data is to determine the thermal transition depth profiles. It is a spectroscopic technique that not only allows the determination of the characteristic temperatures of thermal transitions of thin polymer films but also indirectly allows a probe of their thermal relaxation times, e.g., by using the variable rate cooling technique. Most results published on this topic concern the reduction of the glass transition temperature of thin polymer films due to the reduction of their thickness. However, recent, reliable results indicate that the method of sample preparation has a big impact. There is now strong experimental evidence that a thin layer of molten polymer, the thickness of which depends on temperature, is formed on the free surface of the polymeric material. In addition, the glass transition temperature depends on the thickness of the thin layer and the nature of its interaction with the substrate. Temperature-dependent spectroscopic ellipsometry has great potential for more frequent use in determining the thermo-optical properties of thin polymer layers. The development of experimental techniques undoubtedly contributes to this. For example, ellipsometers are currently available that collect data from the entire spectrum in a split second. The examples of application of this experimental technique discussed here clearly show the high accuracy of ellipsometry in determining the TOC and TEC values of thin polymer layers. The attractiveness of temperature-dependent ellipsometry is also increased due to the possibility of filtering thermal processes that are reversible from those that are irreversible, which are offered by the technique of temperature modulation. 35,148 We hope that our work will contribute to a significant increase in the number of articles presenting a more complete analysis of the thermo-optical properties of polymer layers, especially those that determine the spectral dispersion of the thermo-optic coefficient and which are clearly lacking in the existing literature. It should also be mentioned that there are visible deficiencies in the theoretical foundations of temperature-dependent ellipsometry, e.g., a lack of linear analysis for ellipsometric angles. Therefore, we trust that our article will give the right impulse and initiate appropriate theoretical research in this direction. It can also be assumed that interesting experimental work will further focus on explaining the relationship between T g reduction in irreversibly adsorbed layers, the role of the substrate, packing frustration, and free volume.