Next Article in Journal
Diterpenoid Lactones with Anti-Inflammatory Effects from the Aerial Parts of Andrographis paniculata
Next Article in Special Issue
Non-Isothermal Crystallization Kinetics of Poly(4-Hydroxybutyrate) Biopolymer
Previous Article in Journal
Discovery of (5-Phenylfuran-2-yl)methanamine Derivatives as New Human Sirtuin 2 Inhibitors
Previous Article in Special Issue
Thermal Decomposition of Maya Blue: Extraction of Indigo Thermal Decomposition Steps from a Multistep Heterogeneous Reaction Using a Kinetic Deconvolution Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kinetic Processes in Amorphous Materials Revealed by Thermal Analysis: Application to Glassy Selenium

Department of Physical Chemistry, Faculty of Chemical Technology, University of Pardubice, Studentská 95, 53210 Pardubice, Czech Republic
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(15), 2725; https://doi.org/10.3390/molecules24152725
Submission received: 1 July 2019 / Revised: 20 July 2019 / Accepted: 22 July 2019 / Published: 26 July 2019
(This article belongs to the Special Issue Thermal Analysis Kinetics for Understanding Materials Behavior)

Abstract

:
It is expected that viscous flow is affecting the kinetic processes in a supercooled liquid, such as the structural relaxation and the crystallization kinetics. These processes significantly influence the behavior of glass being prepared by quenching. In this paper, the activation energy of viscous flow is discussed with respect to the activation energy of crystal growth and the structural relaxation of glassy selenium. Differential scanning calorimetry (DSC), thermomechanical analysis (TMA) and hot-stage infrared microscopy were used. It is shown that the activation energy of structural relaxation corresponds to that of the viscous flow at the lowest value of the glass transition temperature obtained within the commonly achievable time scale. The temperature-dependent activation energy of crystal growth, data obtained by isothermal and non-isothermal DSC and TMA experiments, as well as direct microscopic measurements, follows nearly the same dependence as the activation energy of viscous flow, taking into account viscosity and crystal growth rate decoupling due to the departure from Stokes–Einstein behavior.

Graphical Abstract

1. Introduction

Glasses are amorphous materials lacking the periodic atomic arrangement typical for crystalline substances. Structurally, they resemble supercooled liquids but behave mechanically like solids. Figure 1 shows the specific volume or enthalpy as a function of temperature for a typical glass-forming liquid. Upon slow cooling from high temperatures, a liquid may crystallize at Tm forming a stable crystalline material. However, if the cooling through this temperature range is fast enough to avoid nucleation and subsequent crystal growth, a metastable supercooled liquid state is attained.
When a supercooled liquid is cooled by a cooling rate of q+1 to lower temperatures, the internal molecular motion slows down and its viscosity significantly increases. At the glass transition temperature (Tg), the time needed for molecular rearrangement becomes comparable to the experimental time scale. At a lower cooling rate (q+2 < q+1) the supercooled liquid stays in metastable equilibrium until lower temperatures. Therefore, the glass transition is not a true thermodynamic phase transition depending on procedural variables such as cooling rate. As a consequence, there is not a single glassy state, and the properties of the glass depend upon how it was obtained [1].
Clear evidence of the non-existence of a single glassy state is a slow, gradual approach of volume or enthalpy towards the extrapolated supercooled liquid equilibrium line that has been called “structural relaxation”. This process, associated with a slow molecular rearrangement, is experimentally observable in the glass transition range (see Figure 1). It seems that structural relaxation is strongly affected by the supercooled liquid dynamics. An understanding of clues between long time-scale structural relaxation and short time-scale molecular dynamics of corresponding supercooled liquids is of fundamental importance for glass science.
On reheating of a glass, the structural relaxation peak is usually observed just above Tg. At higher temperatures the crystallization process takes place [2]. These processes can be followed by thermal analysis methods such as differential thermal analysis (DTA) or differential scanning calorimetry (DSC). Both these methods are quite frequently used to study the structural relaxation and crystallization behavior in glasses. Usually, the kinetic parameters of such activation energies are extracted from experiments taken at different heating rates. In this paper, the physical meaning of these parameters determined for structural relaxation and crystallization of glassy selenium is analyzed and discussed with respect to the viscosity behavior of supercooled selenium.

2. Results

This section provides a concise description of the structural relaxation experiments as well as the crystallization of glassy selenium obtained by differential scanning calorimetry (DSC) and dilatometry.

2.1. Volume and Enthalpy Relaxation

The isothermal volume relaxation is typically studied by a temperature down-jump experiment. In this experimental set-up the classical mercury-filled dilatometer containing a selenium glass is first equilibrated in a thermal bath at temperature T0. Then it is quickly transferred to another temperature bath at temperature T < T0 and the time-dependent volume contraction is recorded. The temperature up-jump experiment starts by equilibration of a dilatometer at a temperature T0. Then it is quickly transferred to another thermal bath at temperature T > T0 and the time dependent volume expansion is recorded. Isothermal volume relaxation can be expressed as a relative departure of an actual volume V from the equilibrium volume V:
δ = V V V .
The relaxation response immediately after the temperature jump from temperature T0 to T is defined by the following equation:
δ 0 = Δ α T 0 T ,
where Δα corresponds to the difference between the thermal expansion coefficient in a selenium supercooled liquid and selenium glass [3]. Figure 2 shows the isothermal volume relaxation response of selenium glass subjected to the temperature down-jump (T0 = 37 °C, T = 32 °C) and up-jump (T0 = 27 °C, TB = 32 °C). Open circles represent experimental data. The logarithmic time axis is normalized to the initial time that corresponds to the thermal equilibration of dilatometer after the temperature jump (ti = 70 s) [3].
It is clearly seen that both relaxation responses after temperature down-jump and up-jump are non-exponential and non-linear. Such behavior can be described by a non-exponential decay function including the reduced time integral [4]:
δ t = δ 0 exp 0 t d t τ ( T , T f ) β ,
where β is the non-exponentiality parameter (0 < β ≤ 1), inversely proportional to the width of the spectrum of relaxation times. It is assumed that τ depends on temperature T, as well as on the instantaneous structure of amorphous material characterized by the fictive temperature Tf [5]. The relaxation time can be expressed by the following equation [6]:
τ ( T , T f ) = A r e l · exp x E r e l R T + ( 1 x ) E r e l R T f ,
where Arel is the pre-exponential constant, x is the non-linearity parameter (0 < x ≤ 1) and Erel is the effective activation energy of the relaxation process. The time-dependent fictive temperature for this Tool–Narayanaswamy–Moynihan (TNM) model can be expressed as
T f = δ ( t ) δ 0 T 0 T + T .
Full lines in Figure 2 were calculated by using Equations (3)–(5) for the following set of parameters: β = 0.58 ± 0.05, x = 0.42 ± 0.05, ln (Arel/s) = −133.0 ± 0.5 and Erel = 355 ± 2 kJ·mol−1 [3].
Enthalpy relaxation cannot be measured directly in a similar way. However, the DSC heating scans exhibit a typical relaxation overshot just above the glass-transition temperature. It can be shown that normalized heat capacity (CpN) measured by DSC can be expressed [6] as the first derivative of the fictive temperature dTf/dT:
C p N = C p C p g C p l C p g = d T f d T ,
where Cp is measured heat capacity, Cpg is heat capacity of a glass and Cpl is heat capacity of a supercooled liquid.
The temperature dependent plots of CpN are shown in Figure 3. Open circles represent experimental data for heating scans at q+ = 10 K·min−1 taken immediately after the cooling scans performed at the indicated rates. Full lines in Figure 3 were calculated by using Equations (3)–(6) for the following set of parameters: β = 0.65 ± 0.05, x = 0.52 ± 0.05, ln (Arel/s) = −133.0 ± 0.5 and Erel = 355 ± 2 kJ·mol−1 [3]. The parameters ln (Arel/s) and Erel are identical as for the volume relaxation. However, the values of non-exponentiality and non-linearity parameters are higher than those obtained for volume relaxation.
Very similar results are also obtained for the thermal expansion coefficient measured by thermomechanical analysis for As2Se3 glass [7,8]. The agreement between the experiment and the TNM model is very good below Tg. Nevertheless, some deviations are observed in the supercooled liquid above Tg due to viscous flow deformation of the sample. Such effect is not relevant in classical mercury dilatometry or DSC experiments described above.

2.2. Nucleation and Crystal Growth

On further reheating of the supercooled liquid above Tg, the nucleation process takes place being followed by crystal growth. In a selenium supercooled liquid, crystals grow from a relatively low-density nuclei population. Well defined and compact spherulitic structures grow from these centers. It seems that the nucleation has negligible effects during the crystal growth. The crystal growth is usually visible on a microscopic level and therefore it can directly be observed by microscopy methods using hot stage.
Temperature-dependent data for crystal growth velocity u [9,10] and viscosity η [10,11,12] in a selenium supercooled liquid are shown in Figure 4. These two kinetic processes are closely bound together as will be discussed later.
We can assume that in a narrow temperature range the crystal growth velocity can be described by a simple exponential dependence on reciprocal temperature u ~ exp(−EG/RT), that should be linear on a logarithmic scale. The activation energy of crystal growth EG can then be obtained from the slope of such linear dependence,
d log u d 1 / T = E G 2.303 · R ,
where the coefficient 2.303 in Equation (7) comes from the conversion from the natural to the decimal logarithm.
From the crystal growth rate data shown in Figure 4 it is clearly seen that the activation energy is gradually decreasing from a relatively high value just above Tg (≅250 kJ·mol−1) to a considerably lower value (≅40 kJ·mol−1) just below the maximum growth rate.
The heat evolved during the crystal growth can easily be recorded by DSC. Figure 5 shows such heat flow at different scanning rates ranging from 1 to 30 K·min−1.
The crystallization peaks shown in Figure 5 involve the whole crystallization process including nucleation and crystal growth. More detailed calorimetric studies [13,14,15,16] indicate a complex behavior involving both bulk and surface crystal growth with nucleation possibly affected by internal stresses. Nevertheless, nucleation and surface crystal growth are negligible for the bulk selenium sample.
The DSC curves can easily be converted to kinetic data. It is assumed that the fraction crystallized, α, can be obtained by partial integration of non-isothermal heat flow, ϕ, after baseline subtraction:
α = 1 Δ H q + T o n T ϕ · d T ,
where ΔHc correspond to enthalpic change of crystallization, q+ is the heating rate and Ton is the starting point of baseline approximation. The heat flow due to the crystallization process can then be written as
ϕ = Δ H c A · exp E c / R T · f α ,
where A is the preexponential factor and Ec is the activation energy of the crystallization process. The f(α) function corresponds to the Johnson–Mehl–Avrami–Kolmogorov (commonly abbreviated as JMA) model of nucleation-growth process:
f α = m 1 α ln 1 α 1 1 / m .
The development of this equation is described in [8] and references quoted in.
Analysis of DSC curves shown in Figure 5 is complicated due to the strong temperature dependence of heat flow as follows from Equation (9). Nevertheless, it has been shown [17,18] that if the measured heat flow is multiplied by T2 and plotted as a function of α, all data taken at different heating rates collapse to one master curve defined as
z α = f α 0 α d α f α ϕ · T 2 .
This function can be expressed for the JMA model as follows:
z α = m 1 α ln 1 α .
Figure 6 shows the z(α) function obtained by transformation of all crystallization peaks from Figure 5 by Equations (8) and (11). These plots are scaled within the 0 < z(α) ≤ 1 range for easier comparison of different data sets (points). The z(α) function (full line) calculated by Equation (12) fits the experimental data quite well. This confirms the applicability of the JMA model for the description of non-isothermal crystallization kinetics in selenium glass. Another method to test the validity of this model is the shape analysis of DSC curve [19]. Slight data scatter might indicate variation in thermal contacts between the sample and DSC sensor.
The activation energy of the non-isothermal crystallization process Ec can be determined by the Kissinger method [20] from the shift of the maximum of the DSC peak Tp with heating rate q+:
d ln q + / T p 2 d 1 / T p = E c R .
In case of the isothermal data, the Friedman method [21] is usually used, where (dα/dt)α and Tα are the conversion rate and temperature corresponding to arbitrarily chosen values of conversion α:
ln d α / d t α = E R T α + c o n s t . .
The Kissinger and Friedman (for α = 0.50) plots are shown in Figure 7 for all the analyzed non-isothermal and isothermal DSC data, respectively. Similarly, as shown above for crystal growth experiments, also here it is seen that the activation energy significantly changes from relatively low values at a higher temperature (≅40 kJ·mol−1) to higher values at a lower temperature (≅117 kJ·mol−1).
Such important variation of the activation energy makes difficult further analysis of experimental data. The next step of kinetic analysis should be an estimation of the kinetic exponent m in Equation (10). Equation (9) can be rewritten in a somewhat different form,
y α = ϕ · exp E c / R T = Δ H c A · f α .
Figure 8 shows the y(α) function obtained by transformation of all crystallization peaks from Figure 5 by Equation (15) for Ec = 105 kJ·mol−1. These plots are scaled within the 0 < y(α) ≤ 1 range for easier comparison of different data sets (points). Full lines shown in Figure 8 were calculated by Equation (10) for two different values of kinetic exponent (m = 1.5 and 2.5). It is clearly visible that a higher value of the kinetic exponent better fits data taken at lower heating rates, and lower values of m does the same for data taken at higher heating rates. However, this effect is artificial being just a consequence of important changes in activation energy discussed above.
Crystallization in glassy selenium was also measured by means of thermomechanical analysis (TMA). The example curve obtained at 0.2 °C·min−1 is shown in Figure 9 as the temperature dependence of sample height decreased. As was shown in [22], one of most reliable and reproducible characteristic temperatures associated with the TMA crystallization measurements is the extrapolated endset temperature Te—its evaluation is suggested in the figure. The minimum achieved sample height is the point at which the sample deformation caused by viscous flow is ceased by the rigid crystalline structure formed within the sample. The occasional increase of sample height occurring during further heating is the consequence of the ongoing outwards surface crystal growth building up on the stiffened sample profile. The extrapolated endset temperatures can be utilized [23] in a similar way as the characteristic temperatures obtained via calorimetry, i.e., e.g., using the Kissinger equation (Equation (13)). The resulting dependence is for the glassy selenium depicted by the red data and axes in Figure 9. Note the curvature of the Kissinger plot data, similar to the one observed for the DSC crystallization data in Figure 7.

3. Discussion

Thermo-kinetic data described in the previous section imply that the activation energy of most kinetic processes observable in glassy materials by thermal analysis varies with temperature. This gives an interesting opportunity to compare the apparent activation energies among the particular processes for the wide range of experimental conditions. In such a situation, the equilibrium viscosity can be seen as the overarching quantity offering the most sensible comparison for all other data. The activation energy of viscous flow Eη can be obtained from the tangent slope for the viscosity data depicted in Figure 4:
d log η d 1 / T = E η 2.303 · R .
The temperature dependence of Eη is shown in Figure 10 (solid line).
Starting with the structural relaxation process, it is somewhat surprising that the activation energy of this process remains constant throughout the whole measured temperature range (9–45 °C). Note that structural relaxation is generally considered to be very closely interlinked with viscous flow. Nonetheless, the constant value of Erel in the above-given temperature range was unambiguously confirmed from the dilatometric data [3] (see, e.g., Figure 2), where Erel was evaluated by non-linear optimization as well as by the linearization method [23], and also from the calorimetric data [3,24] (see, e.g., Figure 3), where the non-linear optimization was complemented by the newly developed methodology [25] based on the shift of the relaxation peak. Values of Erel are shown in Figure 10 for dilatometry (dark blue) and calorimetry (light blue). It is well apparent that the relaxation data exhibit a constant value of Erel, not only in the glass transition temperature range commonly revealed via the non-isothermal measurement techniques (indicated by the red dashed lines), but also well below these temperatures. In the temperature window where the relaxation times considerably exceed the time-scale characteristic for the non-isothermal measurements, long-term isothermal annealing experiments need to be performed. Interestingly, the intersection of the Erel and Eη well corresponds to the lowest value of Tg obtained during cooling within the commonly achievable time scale. The constancy of Erel is a priori given by the definition of the TNM model, where the potential temperature-dependent component of Erel is replaced by the Tf-based term on the right-hand side of Equation (4). Direct incorporation of the Erel(T) dependence into Equation (4) might, however, be the way towards solving the occasionally raised questions [26] associated with the universality of TNM formalism—as such it is certainly worth of further exploration.
The second process occurring during further heating of the glassy materials is crystal growth, observable either microscopically (see the crystal growth data in Figure 4) or macroscopically, usually via calorimetric methods (see, e.g., Figure 5). Although the two approaches and the temperature ranges of their applicability differ to a great extent, the observed process is essentially the same and the corresponding activation energies should exhibit unified temperature dependence. In order to verify this hypothesis, EG was determined from the crystal growth data depicted in Figure 4 by using Equation (7) (via the direct tangential approach), and Ec was determined from both non-isothermal and isothermal DSC measurements by using the Kissinger (Equation (13)) and Friedman (Equation (14)) methods, respectively (see Figure 7 for the two overall dependences). Again, the temperature-resolved tangential approach to the determination of Ec was adopted. The values of EG and Ec are then compared in Figure 10, showing a very good agreement and confirming the universal nature of the activation energy for the crystal growth process. The non-isothermal DSC crystallization measurements were performed (in addition to the bulk samples) also for a finely powdered glassy selenium. Despite the different crystallization mechanism occurring in case of the fine Se powders [13,14], the corresponding Ec(T) dependence also falls on the crystal growth master-curve depicted in Figure 10, further confirming the universality of this aspect of the crystallization process kinetics.
With regard to the relation between the crystal growth rate and viscosity, it results from the Turnbull–Cohen formula where u is inversely proportional to η. However, it has been shown [27] that for a number of materials the so-called decoupling of these two quantities occurs, breaking the Stokes–Einstein formalism [28]. Formally, the decoupling is described by the apparent decoupling parameter ξa:
ξ a = d log u d log η E G E η .
This rather empirical expression can be further corrected by accounting only for the rate at which the structural entities (atoms, molecules) present in the liquid phase are attached to the growth liquid-crystal interface. This correction is based on the elimination of the term fp from the expression for the crystal growth rate u, where fp is the probability of the structural entity, newly attached to the crystal growth interface, remaining within the crystalline phase:
ξ = d log u f p d log η = d log u log 1 exp Δ G l c / R T d log η = d log u k i n d log η E G , k i n E η .
ΔGlc is the difference between the Gibbs energies of supercooled liquid and crystalline phases. Both forms of the decoupling parameter were essentially calculated from the ratio of the respective activation energies; ΔGlc was calculated based on the standard thermodynamic expression [29] from the selenium enthalpy and entropy of fusion and the heat capacity data published in [30].
Temperature dependences of both forms of the decoupling parameter are shown in Figure 11. Values of ξa and ξ indicate that at low temperatures/growth rates almost no decoupling between u and η can be found (in fact, slight negative decoupling occurs below 340 K). However, as the temperature increases the decoupling becomes more prominent and very well recognizable at temperatures above 370 K. Increasing temperature also results in a rising difference between the two forms of decoupling parameters ξ and ξa. This difference is negligible up to ~365 K but then rapidly increases with the exponential increase of the correction term fp. Note that the melting entropy of crystalline selenium is ΔSm ≈ 1.5R [10]. Such a relatively small value in combination with low supercooling causes relatively large difference between ξ and ξa at higher temperatures. A similar effect may also bring N-type heat capacity of supercooled selenium [30], which effectively lowers the fp contribution. Interestingly, both decoupling dependences exhibit a sudden step-like decrease of the decoupling parameter at 400 K. This is the consequence of the change in preferential morphology of the growing crystallites (transition from spherulitic form B to spherulitic form A, as described in [9]). The step-like change very well agrees with the break on the temperature dependences of the corresponding integral data utilized for calculation of EG and Ec, as depicted in Figure 11, where the dashed vertical line indicates the break point. It is noteworthy that the corrected rate of growth ukin for the high-temperature spherulitic form A does not show further increase in decoupling (which would be represented by the decrease of ξ) and is close to ξ = 0.62.
The last data set depicted in Figure 10 is the one corresponding to the TMA crystallization measurements. Evaluation of activation energy from the TMA data was done similarly as in the case of non-isothermal DSC, i.e., by the tangential approach of the Kissinger dependence (Equation (13)), where the extrapolated endset Te (see Figure 9) was used as the characteristic temperature. It was shown for several other chalcogenide glassy systems [31,32] that the activation energy evaluated in this way from the TMA measurements is in a good agreement with Ec from DSC, but the data are (due to the choice of the extrapolated characteristic temperature) shifted to lower temperatures; that is, the crystal growth process is seemingly observed “in advance” and the EA–T dependence gets shifted to lower temperatures. From this point of view, it is the high-temperature points shown in Figure 10 that in this dependence represent the commonly observed behavior. Contrary to what was observed for most other studied chalcogenides, EA rapidly increases at lower temperatures. The exact position of Te depends on many factors (including applied force, sample geometry, nucleation density, location and morphology of forming crystallites, etc.) but essentially can be understood as the competition between the viscous flow and crystal growth rate (which also depends on viscosity). In case of amorphous selenium, the position of Te appears to be driven more by viscous flow (in combination with surface tension) at low temperatures. Thus, the activation energy determined from the TMA measurements at low q+ gets closer to Eη.
There are of course many consequences associated with the temperature variation of crystallization activation energies. In the last part of the Discussion section we will focus on the model-based master-plot evaluation method utilizing the characteristic kinetic functions z(α) and y(α). The former can be expressed as the product f(α)·g(α), see Equation (11), and as such is invariable with EA. On the other hand, function y(α) is proportional to f(α), which utilizes EA during the transformation of experimental data, see Equation (15). It naturally suggests itself to use the full Ec(T) dependence (see Figure 10) in Equation (15). However, the y(α) function is too sensitive to the value of activation energy, and the large variation of Ec throughout each non-isothermal DSC measurement effectively results in y(α) distortions reminiscent of output obtained for the JMA formalism with the sub-unity kinetic exponent m. It is therefore reasonable to replace the full Ec(T) dependence by constant values of Ec, selected for each measurement individually based on the arbitrarily determined characteristic temperature point. If the value of Ec for the temperature corresponding to the maximum transformation rate (maximum of the DSC peak) is used, the resulting values of JMA kinetic exponents are too large to be physically meaningful. Nevertheless, if we consider that the dimensionality of the formed crystallites is being already set at the start of the crystallization process, and, correspondingly, we utilize Ec values corresponding to α = 0.10, the data depicted in Figure 12 are obtained.
As can be seen, most datapoints can be reasonably described by using the JMA kinetic exponent m = 3, which may correspond to the assumption that for bulk selenium samples the crystal growth starts dominantly via formation of three-dimensional volume-located crystallites.

4. Materials and Methods

Glassy selenium was prepared from pure elements (5N, Sigma-Aldrich, Prague, Czech Republic) by melt-quenching. Elementary Se was melted in an evacuated fused silica ampoule, which was then let to cool in air. The glassy material was crushed in agate mortar and sieved through defined mesh so that the various particle size fractions were obtained. Powder DSC crystallization data reported in this paper were obtained for the 20–50 µm fraction [13,14]. Pieces of glass with a diameter larger than 1 mm were used for the DSC relaxation and bulk crystallization measurements [3,13]. Melt-quench in thin ampoules was utilized to prepare cylindrical Se samples. The ampoules were quenched vertically to obtain a glassy ingot, which was then sawed into samples with the following diameters (dm) and heights (hm): For viscosity measurements dm = 6 mm and hm = 2.5 mm [12], for microscopic crystal growth measurements dm = 4 mm and hm = 2 mm [10] and for TMA crystallization measurements dm = 4 mm and hm = 1 mm.
Experimental setups and details of most measurements were already published in the respective papers: Structural relaxation by DSC and dilatometry in [3,24], viscosity in [12], microscopic crystal growth in [10] and non-isothermal powder and bulk crystallization by DSC in [13,14]. The new, previously unpublished data are those for isothermal bulk crystallization measured by DSC and non-isothermal crystallization measured by TMA.
The isothermal DSC data were obtained using a Q2000 DSC (TA Instruments, Prague, Czech Republic) equipped with a cooling accessory, autolid, autosampler and T-zero Technology. Dry nitrogen was used as the purge gas at a rate of 50 cm3·min−1. The calorimeter was calibrated using In, Zn and H2O. The stability of the DSC signal was checked daily. Open T-zero low-mass pans were used. Regarding the applied temperature program, the sample (8–10 mg) was first subjected to a 5 min isotherm at 45 °C and then heated at 100 °C·min−1 to a selected temperature Ti, where the sample was allowed to isothermally crystallize until the crystallization process was complete. The isothermal crystallization temperatures utilized in the case of each particle size fraction were 100, 105, 110, 115, 120, 125, 130, 135, 140, 145, 150, 155 and 160 °C. In order to obtain a baseline for the isothermal measurement, each DSC pan with a crystalline sample of glass was kept in the DSC cell and the above-described temperature procedure was repeated (in this way the data subtracted from the isothermal crystallization signal truly surrogated the presence of an inert material with similar heat capacity, mass, grain size and positioning in the DSC pan/cell). Perfect flatness of the baseline and reproducibility of the crystallization measurements were confirmed. For the extensive testing of the suitability and repeatability of the initial 100 °C·min−1 heating ramp, see [15].
The thermomechanical measurements were realized by using a TMA Q400EM (TA Instruments), where the cylindrical samples were compressed in-between two alumina plates, and the force applied to the sample was 30 mN. A linear heating rate was applied to study the effect of crystal growth suppressing the decrease of sample height caused by viscous flow. The following heating rates were applied between 35 and 170 °C: 0.2, 0.5, 1, 2, 3, 4 and 5 °C·min−1.

Author Contributions

Conceptualization, J.M.; methodology, J.M.; validation, R.S. and J.M.; formal analysis, J.M. and R.S.; investigation, R.S.; resources, J.M. and R.S.; data curation, R.S.; writing—original draft preparation, J.M. and R.S.; supervision, J.M.; project administration, R.S.; and funding acquisition, R.S.

Funding

This research was funded by Czech Science Foundation grant number 17-11753S.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ediger, M.D.; Angell, C.A.; Nagel, S.R. Supercooled liquids and glasses. J. Phys. Chem. 1996, 100, 13200–13212. [Google Scholar] [CrossRef]
  2. Debenedetti, P.G. Metastable Liquids. Concepts and Principles; Princeton University Press: Princeton, NJ, USA, 1996; ISBN 9780691085951. [Google Scholar]
  3. Málek, J.; Svoboda, R.; Pustková, P.; Čičmanec, P. Volume and enthalpy relaxation of a-Se in the glass transition region. J. Non Cryst. Solids 2009, 355, 264–272. [Google Scholar] [CrossRef]
  4. Naraynaswamy, O.S. A Model of Structural Relaxation in Glass. J. Am. Ceram. Soc. 1971, 54, 491–498. [Google Scholar] [CrossRef]
  5. Tool, A.Q. Relation between inelastic deformability and thermal expansion of glass in its annealing range. J. Am. Ceram. Soc. 1946, 29, 240–253. [Google Scholar] [CrossRef]
  6. DeBolt, M.A.; Easteal, A.J.; Macedo, P.B.; Moynihan, C.T. Analysis of Structural Relaxation in Glass Using Rate Heating Data. J. Am. Ceram. Soc. 1976, 59, 16–21. [Google Scholar] [CrossRef]
  7. Málek, J.; Svoboda, R. Structural Relaxation and Viscosity Behavior in Supercooled Liquids at the Glass Transition. In Thermal Analysis of Micro, Nano- and Non-Crystalline Materials; Šesták, J., Šimon, P., Eds.; Hot Topics in Thermal Analysis and Calorimetry; Springer Science+Business Media: Dordrecht, The Netherlands, 2013; Volume 9, pp. 147–173. ISBN 978-90-481-3149-5. [Google Scholar]
  8. Málek, J.; Shánělová, J. Crystallization Kinetics in Amorphous and Glassy Materials. In Thermal Analysis of Micro, Nano- and Non-Crystalline Materials; Šesták, J., Šimon, P., Eds.; Hot Topics in Thermal Analysis and Calorimetry; Springer Science+Business Media: Dordrecht, The Netherlands, 2013; Volume 9, pp. 291–324. ISBN 978-90-481-3149-5. [Google Scholar]
  9. Ryschenkow, G.; Faivre, G. Bulk crystallization of liquid selenium—Primary nucleation, growth kinetics and modes of crystallization. J. Cryst. Growth 1988, 87, 221–235. [Google Scholar] [CrossRef]
  10. Málek, J.; Barták, J.; Shánělová, J. Spherulitic Crystal Growth Velocity in Selenium Supercooled Liquid. Cryst. Growth Des. 2016, 16, 5811–5821. [Google Scholar] [CrossRef]
  11. Bernatz, K.; Echeveria, I.; Simon, S.; Plazek, D. Characterization of the molecular structure of amorphous selenium using recoverable creep compliance measurements. J. Non Cryst. Solids 2002, 307, 790–801. [Google Scholar] [CrossRef]
  12. Koštál, P.; Málek, J. Viscosity of selenium melt. J. Non Cryst. Solids 2010, 356, 2803–2806. [Google Scholar] [CrossRef]
  13. Svoboda, R.; Málek, J. Crystallization kinetics of a-Se, part 1: Interpretation of kinetic functions. J. Therm. Anal. Calorim. 2013, 114, 473–482. [Google Scholar] [CrossRef]
  14. Svoboda, R.; Málek, J. Crystallization kinetics of a-Se, part 2: Deconvolution of a complex process—The final answer. J. Therm. Anal. Calorim. 2014, 115, 81–91. [Google Scholar] [CrossRef]
  15. Svoboda, R.; Málek, J. Crystallization kinetics of a-Se, part 3: Isothermal data. J. Therm. Anal. Calorim. 2015, 119, 1363–1372. [Google Scholar] [CrossRef]
  16. Svoboda, R.; Gutwirth, J.; Málek, J. Crystallization kinetics of a-Se, part 4: Thin films. Philos. Mag. 2014, 94, 3036–3051. [Google Scholar] [CrossRef]
  17. Málek, J. The Kinetic Analysis of Non-Isothermal Data. Thermochim. Acta 1992, 200, 257–269. [Google Scholar] [CrossRef]
  18. Málek, J. The applicability of Johnson-Mehl-Avrami model in the thermal analysis of the crystallization kinetics of glasses. Thermochim. Acta 1995, 267, 61–73. [Google Scholar] [CrossRef]
  19. Málek, J.; Criado, J.M. The Shape of a Thermoanalytical Curve and Its Kinetic Information Content. Thermochim. Acta 1990, 164, 199–209. [Google Scholar] [CrossRef]
  20. Kissinger, H.E. Reaction kinetics in differential thermal analysis. Anal. Chem. 1957, 29, 1702–1706. [Google Scholar] [CrossRef]
  21. Friedman, H.L. Kinetics of thermal degradation of char-forming plastics from thermogravimetry. Application to a phenolic plastic. J. Polym. Sci. Part C 1964, 6, 183–195. [Google Scholar] [CrossRef]
  22. Zmrhalová, Z.; Pilný, P.; Svoboda, R.; Shánělová, J.; Málek, J. Thermal properties and viscous flow behavior of As2Se3 glass. J. Alloys Compd. 2016, 655, 220–228. [Google Scholar] [CrossRef]
  23. Málek, J. Rate-determining factors for structural relaxation in non-crystalline materials II. Normalized volume and enthalpy relaxation rate. Thermochim. Acta 1998, 313, 181–190. [Google Scholar] [CrossRef]
  24. Svoboda, R.; Pustková, P.; Málek, J. Relaxation behavior of glassy selenium. J. Phys. Chem. Sol. 2007, 68, 850–854. [Google Scholar] [CrossRef]
  25. Svoboda, R. Novel equation to determine activation energy of enthalpy relaxation. J. Therm. Anal. Calorim. 2015, 121, 895–899. [Google Scholar] [CrossRef]
  26. Kovacs, A.J. Transition vitreuse dans les polymères amorphes. Etude phénoménologique. Fortschr. Hochpolym. Forsch. 1963, 3, 394–507. [Google Scholar] [CrossRef]
  27. Ediger, M.D.; Harrowell, P.; Yu, L. Crystal growth kinetics, exhibit a fragility-dependent decoupling from viscosity. J. Chem. Phys. 2008, 128, 034709. [Google Scholar] [CrossRef] [PubMed]
  28. Gutzow, I.S.; Schmelzer, J.W.P. The Vitreous State: Thermodynamics, Structure, Rheology, and Crystallization; Springer: Berlin/Heidelberg, Germany, 2013; ISBN 978-3-642-34632-3. [Google Scholar]
  29. Busch, R.; Kim, Y.J.; Johnson, W.L. Thermodynamics and kinetics of the undercooled liquid and the glass transition of the Zr41.5Ti13.8Cu12.5Ni10.0Be22.5 alloy. J. Appl. Phys. 1995, 77, 4039–4043. [Google Scholar] [CrossRef]
  30. Svoboda, R.; Málek, J. Thermal behavior in Se-Te chalcogenide system: Interplay of thermodynamics and kinetics. J. Chem. Phys. 2014, 141, 224507. [Google Scholar] [CrossRef] [PubMed]
  31. Svoboda, R.; Brandová, D.; Chromčíková, M.; Setnička, M.; Chovanec, J.; Černá, A.; Liška, M.; Málek, J. Se-doped GeTe4 glasses for far-infrared optical fibers. J. Alloys Compd. 2017, 695, 2434–2443. [Google Scholar] [CrossRef]
  32. Svoboda, R.; Brandová, D.; Chromčíková, M.; Liška, M. Thermokinetic behavior of Ga-doped GeTe4 glasses. J. Non Cryst. Solids 2019, 512, 7–14. [Google Scholar] [CrossRef]
Sample Availability: Samples of the amorphous selenium are not available from the authors.
Figure 1. Schematic diagram showing the change in enthalpy and volume during glass formation, structural relaxation and crystallization.
Figure 1. Schematic diagram showing the change in enthalpy and volume during glass formation, structural relaxation and crystallization.
Molecules 24 02725 g001
Figure 2. The volume change of Se glass subjected to a temperature down-jump and up-jump ±5 °C. Points correspond to experimental data obtained by dilatometry. Full lines were calculated for the Tool–Narayanaswamy–Moynihan (TNM) model (parameters in text).
Figure 2. The volume change of Se glass subjected to a temperature down-jump and up-jump ±5 °C. Points correspond to experimental data obtained by dilatometry. Full lines were calculated for the Tool–Narayanaswamy–Moynihan (TNM) model (parameters in text).
Molecules 24 02725 g002
Figure 3. The normalized heat capacity of Se glass reflecting the structural relaxation in the glass transition range. Points correspond to selected experimental data-curves obtained by differential scanning calorimetry (DSC). Full lines were calculated for the TNM model (parameters in text).
Figure 3. The normalized heat capacity of Se glass reflecting the structural relaxation in the glass transition range. Points correspond to selected experimental data-curves obtained by differential scanning calorimetry (DSC). Full lines were calculated for the TNM model (parameters in text).
Molecules 24 02725 g003
Figure 4. The crystal growth rate of isothermally grown spherulitic crystals and viscosity of a selenium supercooled liquid as a function of reciprocal temperature.
Figure 4. The crystal growth rate of isothermally grown spherulitic crystals and viscosity of a selenium supercooled liquid as a function of reciprocal temperature.
Molecules 24 02725 g004
Figure 5. The DSC measurement of bulk selenium glass at different heating rates.
Figure 5. The DSC measurement of bulk selenium glass at different heating rates.
Molecules 24 02725 g005
Figure 6. The z(α) function obtained by transformation of DSC data shown in Figure 5. Points were calculated by using Equation (11) (numbers indicate heating rate). The full line was calculated by Equation (12).
Figure 6. The z(α) function obtained by transformation of DSC data shown in Figure 5. Points were calculated by using Equation (11) (numbers indicate heating rate). The full line was calculated by Equation (12).
Molecules 24 02725 g006
Figure 7. Determination of the crystallization activation energy for selenium glass by the Kissinger method (non-isothermal data) and the Friedman method (isothermal data).
Figure 7. Determination of the crystallization activation energy for selenium glass by the Kissinger method (non-isothermal data) and the Friedman method (isothermal data).
Molecules 24 02725 g007
Figure 8. The y(α) function obtained by transformation of DSC data shown in Figure 5. Points were calculated by using Equation (15) (numbers indicate heating rate). Full lines were calculated by Equation (10).
Figure 8. The y(α) function obtained by transformation of DSC data shown in Figure 5. Points were calculated by using Equation (15) (numbers indicate heating rate). Full lines were calculated by Equation (10).
Molecules 24 02725 g008
Figure 9. Black-based data and axes show an example thermomechanical analysis (TMA) crystallization curve (obtained at 0.2 °C·min−1); evaluation of the extrapolated endset is illustrated. Red-based data and axes show the Kissinger plot for the TMA crystallization data.
Figure 9. Black-based data and axes show an example thermomechanical analysis (TMA) crystallization curve (obtained at 0.2 °C·min−1); evaluation of the extrapolated endset is illustrated. Red-based data and axes show the Kissinger plot for the TMA crystallization data.
Molecules 24 02725 g009
Figure 10. Activation energies of viscous flow (solid line), microscopic crystal growth, macroscopic crystallization observed by DSC (two data series, for bulk and powdered samples) and TMA, as well as structural relaxation observed by DSC and dilatometry.
Figure 10. Activation energies of viscous flow (solid line), microscopic crystal growth, macroscopic crystallization observed by DSC (two data series, for bulk and powdered samples) and TMA, as well as structural relaxation observed by DSC and dilatometry.
Molecules 24 02725 g010
Figure 11. Temperature dependences of decoupling parameters ξa and ξ (right-hand axis), crystal growth rate u (black outer left-hand axis) and logarithm of heat flow corresponding to α = 0.50 obtained during isothermal DSC experiments (red inner left-hand axis). The vertical dashed line indicates the transition between the spherulitic B and A crystallite forms.
Figure 11. Temperature dependences of decoupling parameters ξa and ξ (right-hand axis), crystal growth rate u (black outer left-hand axis) and logarithm of heat flow corresponding to α = 0.50 obtained during isothermal DSC experiments (red inner left-hand axis). The vertical dashed line indicates the transition between the spherulitic B and A crystallite forms.
Molecules 24 02725 g011
Figure 12. The y(α) functions obtained by transformation of DSC data shown in Figure 5. Points were calculated by using Equation (15) (numbers indicate heating rate); Ec utilized in the calculations were determined from the corresponding dependence shown in Figure 10 for the temperatures matching α = 0.1 of each respective DSC data curve. The solid red line was calculated by Equation (10) for m = 3.
Figure 12. The y(α) functions obtained by transformation of DSC data shown in Figure 5. Points were calculated by using Equation (15) (numbers indicate heating rate); Ec utilized in the calculations were determined from the corresponding dependence shown in Figure 10 for the temperatures matching α = 0.1 of each respective DSC data curve. The solid red line was calculated by Equation (10) for m = 3.
Molecules 24 02725 g012

Share and Cite

MDPI and ACS Style

Málek, J.; Svoboda, R. Kinetic Processes in Amorphous Materials Revealed by Thermal Analysis: Application to Glassy Selenium. Molecules 2019, 24, 2725. https://doi.org/10.3390/molecules24152725

AMA Style

Málek J, Svoboda R. Kinetic Processes in Amorphous Materials Revealed by Thermal Analysis: Application to Glassy Selenium. Molecules. 2019; 24(15):2725. https://doi.org/10.3390/molecules24152725

Chicago/Turabian Style

Málek, Jiří, and Roman Svoboda. 2019. "Kinetic Processes in Amorphous Materials Revealed by Thermal Analysis: Application to Glassy Selenium" Molecules 24, no. 15: 2725. https://doi.org/10.3390/molecules24152725

APA Style

Málek, J., & Svoboda, R. (2019). Kinetic Processes in Amorphous Materials Revealed by Thermal Analysis: Application to Glassy Selenium. Molecules, 24(15), 2725. https://doi.org/10.3390/molecules24152725

Article Metrics

Back to TopTop