Here, we present a protocol for cooling rate dependent ellipsometry experiments, which can determine the glass transition temperature (Tg), average dynamics, fragility and the expansion coefficient of the super-cooled liquid and glass for a variety of glassy materials.
This report aims to fully describe the experimental technique of using ellipsometry for cooling rate dependent Tg (CR-Tg) experiments. These measurements are simple high-throughput characterization experiments, which can determine the glass transition temperature (Tg), average dynamics, fragility and the expansion coefficient of the super-cooled liquid and glassy states for a variety of glassy materials. This technique allows for these parameters to be measured in a single experiment, while other methods must combine a variety of different techniques to investigate all of these properties. Measurements of dynamics close to Tg are particularly challenging. The advantage of cooling rate dependent Tg measurements over other methods which directly probe bulk and surface relaxation dynamics is that they are relatively quick and simple experiments, which do not utilize fluorophores or other complicated experimental techniques. Furthermore, this technique probes the average dynamics of technologically relevant thin films in temperature and relaxation time (τα) regimes relevant to the glass transition (τα > 100 sec). The limitation to using ellipsometry for cooling rate dependent Tg experiments is that it cannot probe relaxation times relevant to measurements of viscosity (τα << 1 sec). Other cooling rate dependent Tg measurement techniques, however, can extend the CR-Tg method to faster relaxation times. Furthermore, this technique can be used for any glassy system so long as the integrity of the film remains throughout the experiment.
The seminal work of Keddie Jones and Corey1 showed that the glass transition temperature (Tg) of ultra-thin polystyrene films decreases with respect to the bulk value at thicknesses lower than 60 nm. Ever since, many experimental studies2-11 have supported the hypothesis that the observed reductions in Tg are caused by a layer of enhanced mobility near the free surface of these films. However, these experiments are indirect measures of a single relaxation time, and thus there is a debate12-18 centered on a direct correlation between average thin film dynamics and the dynamics at the air/polymer interface.
To answer this debate, many studies have directly measured the dynamics of the free surface (τsurface). Nanoparticle embedding,19,20 nanohole relaxation,21 and fluorescence22 studies show that the air/polymer interface has dynamics orders of magnitude faster than the bulk alpha relaxation time (τα) with a much weaker temperature dependence than that of τα. Because of its weak temperature dependence, the τsurface of these films,19-22 and enhanced dynamics of thin polystyrene films,23,24 intersects the bulk alpha relaxation (τα) at a single point T*, which is a few degrees above Tg, and at a τα of ≈ 1 sec. The presence of T* could explain why experiments which probe relaxation times faster than * fail to see any thickness dependence on the Tg of ultra-thin Polystyrene films.13-18 Lastly, while direct measurements of the enhanced mobile layer show that it has a thickness of 4-8 nm,20-22 there is evidence that the propagation length of the dynamics at the air/polymer interface is much larger than the thickness of the mobile surface layer.5,25,26
This report aims to fully describe a protocol for using ellipsometry for cooling rate dependent Tg (CR-Tg) experiments. CR-Tg have been previously used to describe the average dynamics of ultra-thin films of polystyrene.23,24,27,28 Furthermore, This technique was recently used to show a direct correlation between the average dynamics in ultra-thin polystyrene films, and the dynamics at the free surface.23 The advantage of CR-Tg measurements over other types of measurements such as fluorescence, nanoparticle embedding, nanohole relaxation, nanocalorimetry, dielectric spectroscopy, and Brillouin light scattering, studies is that they are relatively quick and simple experiments that do not utilize fluorophores or other complicated experimental techniques. Recent advances in spectroscopic ellipsometry allow this technique to be used to efficiently determine the optical properties of ultra-thin films of polymers and other types of hybrid materials with exceptional accuracy. As such, this technique probes the average dynamics of technologically applicable thin films in temperature and time regimes relevant to the glass transition (T ≤ Tg, τα ≥ 100 sec). Furthermore, this technique will provide information on the expansion coefficients of the glassy and supper cooled liquid states as well as the fragility of the system, which can then be compared with data for bulk films. Lastly, CR- Tg experiments can be used for any glassy system so long as the integrity of the film remains throughout the experiment.
1. Film Preparation
2. Determining Film Thickness
3. Cooling Rate Dependent Tg Measurements
4. Determining Values of Tg
5. Analyzing Average Thin Film Dynamics
Fitting Raw Ellipsometry Data
Polystyrene films are transparent in the wavelength range of the ellipsometer (500-1,600 nm). Thus a Cauchy model is a good model for describing the index of refraction of polystyrene films. Figure 1A shows an example of Ψ(λ) and Δ(λ) for a thick (274 nm) film of polystyrene, and the resulting fit to the Cauchy model . For films thicker than 10 nm, both the A and B parameters of the Cauchy equation should be fit to accurately model the wavelength dependence of the index of refraction. The Cauchy model is only physical when n is a decreasing function of wavelength, λ. Figure 1B shows an example of a physical index as evident by the always decreasing value of n and k=0. For films thinner than 10 nm, the short path length of light means only the A parameter in the Cauchy equation should be fit. In these extremely thin films, having B as an open fit parameter can drive the ellipsometry fit to an unphysical index, even if the "fit" of Ψ(λ) and Δ(λ) has a small mean squared error (MSE). Such an example can be seen in Figure 2. For some materials it may be necessary to fit higher order terms in the Cauchy model or use a more sophisticated optical model in order to accurately fit optical properties.
Figure 1. Physical Ellipsometry Fit. (A) An example of Ψ(λ) (red solid line) and Δ(λ) (green solid line) of a 110 nm film of polystyrene, and the resulting fit (black dashed line). (B) An example of the physical index n (red line) and k (blue line) produced by the fit in part A. Please click here to view a larger version of this figure.
Figure 2. Unphysical Ellipsometry Fit. (A) An example of Ψ(λ) (red solid line) and Δ(λ) (green solid line) of an 8 nm film of polystyrene, and the resulting fit (black dashed line). (B) An example of the unphysical index n (red line) and k (blue line) produced by the fit in part A. Please click here to view a larger version of this figure.
Fitting Cooling Rate Dependent Tg Experiments
When fitting the thickness of a film throughout the temperature profile, it is important to remember both the polystyrene film and the Si wafer substrate will expand, and their optical properties will change with temperature. Thus, in order to calculate accurate expansion coefficients, the index of the Si substrate must be fit with a temperature dependent model to account for changes in the optical properties of Si. An easy way to check to see if the Si substrate is modeled correctly is to see if the fit's MSE changes significantly with temperature. Figure 3A shows an example of a thickness, temperature, and MSE profile for a fit that models the temperature dependence of the index of Si correctly, while Figure 3B shows the same profiles when the fit that does not correctly account for the changes in the optical properties of Si substrate. Notice that the MSE values in Figure 3B vary greatly with temperature. The MSE decrease in Figure 3A is due to switching from an acquisition time of 1 sec to 3 sec.
Figure 3. Cooling Rate Tg profiles. (A) An example of a typical temperature, thickness, and MSE profile for a single CR-Tg experiment on a 110 nm polystyrene film when correctly accounting for the temperature dependent index of the Si Substrate. (B) An example of a typical temperature, thickness, and MSE profile for a single CR-Tg experiment on the same film when incorrectly accounting for the temperature dependent index of the Si Substrate. Please click here to view a larger version of this figure.
Assigning Tg
The Tg can be calculated from a thickness vs. temperature plot for a given cooling ramp. Figure 4 shows an example of such a curve. The Tg is defined as the temperature where a supercooled liquid falls out of equilibrium upon cooling. In these ellipsometry experiments, the Tg is defined as the temperature at which linear fits to the supercooled liquid and glassy regimes intersect. Figure 4 highlights these regimes as red and blue, respectively. These regimes should be chosen such that the calculated expansion coefficients agree with previous bulk measurements, if available. This method would eliminate subjectivity from the selection process, which could lead to artificially high or low expansion coefficients, and therefore less accurate measures of Tg. Additionally, the expansion coefficients should be independent of film thickness and cooling rate, which can provide guidance in cases where bulk values of expansion coefficient are not available. The expansion coefficients can be calculated by dividing the slopes of the two regimes by the film thickness. Using this method for determining Tg, the Tg for a 110 nm film of polystyrene is measured to be 372 ± 2 K at 10 K/min, and the expansion coefficients of the supper-cooled liquid and glass are 5.7 x 10-4 ± 3 x 10-5 K-1 and 1.5 x 10-4 ± 3 x 10-5 K-1, respectively, which agree well with previously determined values.29 The errors on the values of Tg, and the expansion coefficients are a result of reasonable changes in the selected regions for the super-cooled and glassy regimes.
Figure 4. Assigning Tg. A typical plot of thickness vs. temperature for a 110 nm film of 342 kg/mol PS at a cooling rate of 10 K/min. The shaded parts of the curve represent the super-cooled liquid (red) and glassy (blue) regimes chosen for the purposes of assigning Tg. Tg is defined as the temperature at which the two linear fits intersect. Using this method, the Tg for a 110 nm film of polystyrene is measured to be 372 ± 2 K at 1 K/min and the expansion coefficients of the supper-cooled liquid and glass are 5.7 x 10-4 ± 3 x 10-5 K-1 and 1.5 x 10-4 ± 3 x 10-5 K-1, respectively. Please click here to view a larger version of this figure.
Analyzing Average Film Dynamics
The cooling rate dependent Tg data can be related to the average relaxation time at Tg through the empirical relation that at a cooling rate of 10 K/min, the system falls out of equilibrium when the average relaxation time is equal to 100 sec, i.e., cooling rate x τα ≈ 1000.24 Applying this relation to the data in Figure 5A, a plot of log(cooling rate) vs. 1/Tg (Figure 5B) can be used to evaluate how accurate this relation is for polystyrene, and how well the CR-Tg method describes bulk dynamics for a thick film. The red data in Figure 5B are the bulk dynamics of polystyrene as determined via dielectric spectroscopy.16 While the cooling rate x τα ≈ 1000 relation is purely empirical, and may change slightly based on the experimental technique used to determine bulk dynamics, or the specific glass former being tested,30,31 Figure 5B shows that The cooling rate dependent Tg data for a 110 nm film of polystyrene agrees well with this data. This figure also shows that CR-Tg can be used to extend the dynamical range of the measurements to low temperature, which are usually not accessible by dielectric relaxation measurements. Furthermore, the slope of a linear fit of the Log(CR) vs. 1/Tg data is related to the activation energy of the glass transition. This activation energy relates to the fragility (m) of the glassy film at Tg by the relation;
The second term is only correct if an Arrhenius fit to the data is used as an approximation. Using this method, the fragility for a 110 nm PS film is measured to be 162 ± 21. This value is in good agreement with reported values for bulk polystyrene in the literature (150) from dynamic scanning calorimetry measurements.32
Figure 5. Analyzing Average Film Dynamics via CR-Tg Experiments. (A) A Plot of Tg vs. Cooling Rate for a 110 nm film of polystyrene. (B) Plot of Log (Cooling Rate) vs. 1000/Tg for the same film (black circles). With the relation (cooling rate) x τ = 1000, the results of a CR-Tg experiment on 110 nm PS are plotted alongside direct measures of bulk dynamics of PS, using dielectric relaxation16 with no further shifting factors (red open squares). The red dashed line is a Volgel Fulcher Tammann equation fit to the dielectric relaxation data from reference 16. The resultant fit parameters are τ0 = 1012, B = 13,300 K, and T0 = 332 K. The value of T* from Ref 23 is plotted here as a blue star. From the plot, the fragility is measured to be 162 ± 21. This value is in good agreement with previously reported values in the literature (150).32 Please click here to view a larger version of this figure.
Cooling-Rate dependent Tg measurements are high throughput characterization experiments that can determine the Tg, the expansion coefficient of the glass and the super-cooled liquid, the temperature dependence of the average dynamics, and the fragility of a particular glassy material in a single experiment. Furthermore, unlike fluorescence, embedding, or nanohole relaxation experiments, CR-Tg experiments are relatively quick and simple because they do not utilize fluorophores or other complicated experimental techniques. Due to the sensitivity of ellipsometry, this method can be used on films of thicknesses as thin as a few nanometers and as thick as a few microns, so long as the fitting procedure is correct. This allows for quick and simple analysis of both the temperature dependence and thickness dependence of the average dynamics and fragility.
In order to perform these measurements successfully, there are a few critical steps where extra care must be taken. It is imperative that the ellipsometry fit be correct. As explained previously, it is critical that the temperature dependence of the optical properties of the Si substrate be taken into account. Failing to do this could lead to incorrect values of Tg and incorrect values of the expansion coefficient. Also, it is important to clamp the film firmly to the heating element. This helps ensure good thermal contact, which is imperative for accurately defining Tg values at fast rates. Finally, when assigning values of Tg, the chosen supercooled liquid and glassy regimes must not include the glass transition itself. The glass transition is defined as the part of transition where the slope of the thickness vs. temperature data is changing between the supercooled liquid and glassy regimes. Including this change in slope in either linear fit would artificially change the calculate value of Tg. To remove subjectivity from the selection process, choose super-cooled liquid and glassy regimes which produce expansion coefficients that agree with reported values.
Another advantage of this protocol is that it can be amended to allow for the analysis of any glass-former. The only aspect of this protocol that would need to be modified to test the dynamics of a different glass former is the temperature profile. As long as the bulk Tg of the glass former is known, the maximum and minimum temperatures can be altered to ensure that the film undergoes a glass transition, but also doesn't degrade. The maximum temperature should be approximately Tg + 20 K, and the minimum temperature should be at least Tg – 40 K. Also, the chosen cooling rates can be varied to probe other time scales of interest for a particular type of polymer film.
Despite its advantages, there are limitations to this technique. Because this technique indirectly probes an average relaxation time through the cooling rate of the experiment, the time scales this method probes are limited to the maximum cooling rate available by the method of temperature control. For the ellipsometry procedure presented here, the fastest available cooling rate is 150 K/min, which relates to a relaxation time of τ = 6.66 sec. While this time scale is slow enough to be relevant to the glass transition, it is much slower that the time scales relevant to the viscosity of polymer melts. Such time scales usually are determined via rheology or dielectric spectroscopy, but CR-Tg measurements can probe these time scales if the cooling rate is fast enough. This can easily be achieved using nanocalorimetry or flash DSC.33,34
Because of the high-throughput nature of this technique, it allows for many different kinds of materials to be tested. While this report focused on CR-Tg experiments of polystyrene films, this same method could easily be applied to a range of glassy materials from long chain polymers to small organic molecules used in organic electronic technologies. As long as the integrity of the film holds through the experiment, the temperature dependence and thickness dependence of average dynamics and fragility can be determined.
The authors have nothing to disclose.
The Authors would like to acknowledge James A. Forrest for help in the initial idea for this technique.26 This work was supported by funding from the University of Pennsylvania and was partially supported by the MRSEC program of the National Science Foundation under award no. DMR-11- 20901 at the University of Pennsylvania.
Toluene | Sigma Aldrich | 179418-1L | This can be purchased from any chemical company. |
Atactic Polystyrene | Polymer Source Inc. | P-4092-S | This can be purchased from any chemical company. |
THMS 600 temperature stage | Linkam | THMS 600 | any temperature stage that can be fit to an ellipsometer could be used. |
M2000V Spectroscopic Ellipsometer | J.A. Woollam | M200V | This procedure should be applicable for any spectroscopic ellipsometer. |
Spin Coater | Laurell Technologies | WS-650-23B | This Procedure is possible with any spin coater |
Sample vials | Fisher Scientific | 02-912-379 | Any sample vials will do |
Silicon wafers | Virginia semi conductors | 325S1410694D |