** Journal of Crystallization Process and Technology ** Vol. 2 No. 4 (2012) , Article ID: 23436 , 8 pages DOI:10.4236/jcpt.2012.24025

Investigation of Crystallization Kinetics in Glassy Se and Binary Se_{98}M_{2} (M = Ag, Cd, Zn) Alloys Using DSC Technique in Non-Isothermal Mode

^{ }^{ }^{ }^{}

Department of Physics, Glass Science Laboratory, Banaras Hindu University, Varanasi, India.

Email: ^{*}dr_neeraj_mehta@yahoo.co.in

Received July 11^{th}, 2012; revised August 16^{th}, 2012; accepted September 7^{th}, 2012

**Keywords:** Differential Scanning Calorimetry; Crystallization Kinetics; Chalcogenide Glasses

ABSTRACT

The crystallization kinetics of glassy Se and binary Se_{98}M_{2} (M = Ag, Cd, Zn) alloys have been studied at different heating rates (5, 10, 15, 20 K·min^{–}^{1}) using Differential Scanning Calorimetric (DSC) technique. The crystallization temperature (T_{c}) is determined from exothermic peak obtained in DSC scans of present samples. The variation in peak crystallization temperature (T_{c}) with the heating rate (β) has been used to investigate the growth kinetics using Kissinger, Augis-Bennet and Matusita-Sakka models. The activation energy of crystallization (E_{c}) has been found to increase with Ag additive and to decrease with Zn and Cd additive. The value of various kinetic parameters such as rate constant (K_{p}), Avrami index (n), thermal stability (S) and Hruby number (H_{r}) have been calculated under non-isothermal mode. The maximum change in different kinetic parameters has been found after the incorporation of Ag additive.

1. Introduction

Thermal analysis of chalcogenide glasses in term of crystallization kinetics is of particular interest from fundamental and applications point of views. Chalcogenide glasses have good infrared transparency, high refractive index and low optical loss, so the possible application of these glasses are in optical wave guide, infrared optical fibres, chemical sensors, bio-sensor, photoreceptors in copying machines and X-ray imaging plates, I.R. optical lenses, windows and high sensitivity ionic sensors [1-6]. Crystallization kinetics is an important study in the amorphous as well as glassy solids. The thermal stability and the dynamics of crystallization process can be determined by formulism of crystallization kinetics. The crystallization kinetics depend on various parameters; such as the mode of crystallization, the number of quenched nuclei, the activation energy of diffusion and also the difference in free energy between the amorphous and possible crystalline phases [7-10].

The reason of selection of glassy Se based chalcogenide alloys in the present study are the fact that these material exhibits reversible phase transformation which has a wide range of application in memory and data storage devices. The phase transformations of chalcogenide glasses in binary and ternary systems have been widely studied by nonisothermal technique, because of their usage in numerous applications in semiconductor technology [1-6].

The effect of various metallic additives on the different physical properties such as thermally activated d.c. conduction and a.c. conduction, glass transition and crystallization kinetics of glassy Se has been studied by various workers by changing the composition but fewer studies have been reported by changing the additives. The effect of incorporation of Ag, Zn and Cd on the glass forming ability, hardness and kinetics of crystallization of the binary system SSe_{20} has been studied by El-Den [11]. He found that there is an increase in the value of the hardness after incorporation of the added metals. He explained this result in terms of increase in the strength of chemical bonds in ternary alloys. Singh et al. [12] reported the d.c. conductivity measurements at high electric fields in vacuum-evaporated thin films of amorphous Se_{80}Te_{20}, Se_{75}Te_{20}Ge_{5} and Se_{75}Te_{20}Sb_{5} systems. From analysis of the experimental data they confirmed the presence of space charge limited conduction (SCLC) in all the glassy materials. They observed an increase in density of defect states (DOS) in case of Sb and a decrease in (DOS) for Ge. The change in DOS on impurity incorporation is explained in terms of the change in structure of these glasses.

The temperature dependence of direct current (dc) conductivity has been reported by Khan et al. in thin films of glassy Se_{95}M_{5} (where M = Ga, Sb, Bi) in order to identify the conduction mechanism and to observe the doping effect of different metals on amorphous selenium [13]. They found that the conduction in high temperature range (314 - 375 K) is due to thermally activated tunnelling of charge carriers in the band tails of localized states; and in the low temperature range (219 - 314 K) conduction takes place through variable range hopping in the localized states near the Fermi level. They also reported the existence of space charge limited conduction (SCLC) in these glassy alloys. Sharma and Kumar studied the influence of additives (Sb and Ag) on a.c. conductivity and dielectric properties of glassy Se_{70}Te_{30} alloy [14]. Their experimental results indicate that a.c. conductivity and dielectric parameters depend on temperature, frequency and the impurity incorporated in Se-Te glassy system. They found that the Ag and Sb additive have opposite influence in the dielectric parameters when they incorporated at cost of Te. The results obtained by them are being correlated by the nature of covalent character of the studied composition and with the change in density of defect states.

These days, our laboratory is also engaged in this paticular direction to find some fruitful results from the significant and remarkable problem of “Metal-induced Effects”. Recently, Dohare et al. [15] have observed the effect of additives on the crystallization kinetics of Se_{70}Te_{30} with the same composition of additives as a chemical modifier. In the present work, investigation on the crystallization kinetics of the glassy Se and glassy Se_{98}M_{2} (M = Ag, Zn, Cd) alloys were carried out by means of non-isothermal DSC measurements.

2. Material Preparation

Glassy Se and Se_{98}M_{2} (M = Ag, Zn, Cd) alloys were prepared by quenching technique. The accurate proportion of high purity 5N (99.999%) elements, in agreement with their atomic percentages, were weighed using an electronic balance with the least count of 0.1 mg. The materials were then sealed in evacuated (~10^{–}^{5 }Torr) quartz ampoules (length ~ 5 cm and internal diameter ~ 8 mm). Each ampoule was kept inside the furnace at 800˚C (where the temperature was raised at a rate of 3˚C/min - 4˚C/min). During heating, all the ampoules were continuously rocked, by revolving a ceramic rod to which the ampoules were tucked away in the furnace for obtaining homogeneous glassy alloys.

After rocking for about 12 hours, the obtained melts were cooled rapidly by removing the ampoules from the furnace and dipping in ice-cooled water rapidly. The ingots of the samples were then taken out by breaking the quartz ampoules.

3. Experimental

The confirmation of glassy phase of as prepared samples was done by using XRD technique. X-ray diffractometer (Simadzu PW1710, Japan) with CuK_{α} radiation (λ = 1.5405 Å) was used for this purpose. The X-ray tube voltage and current were 30 kV and 20 mA, respectively. The scan range was 5˚ - 80˚ (2θ) and the scan speed was 1˚/min. The XRD pattern of as-prepared sample of glassy Se_{98}Cd_{2} alloy is shown in Figure 1. The glassy nature of other samples (Se, Se_{98}Ag_{2} and Se_{98}Zn_{2}) alloys was also confirmed by XRD technique.

The glasses, thus prepared, were ground to make fine powder for DSC studies (flow type). 10 mg to 20 mg of the powder was heated at constant heating rate and the changes in heat flow with respect to an empty reference pan were measured. The thermal behaviour was investigated using differential scanning calorimeter. The temperature precision of this equipment is ±0.1 K with an average standard error of about ±1 K in the measured values.

4. Results

Measurements were made under almost identical conditions at four different heating rates 5, 10, 15 and 20 Kmin^{–}^{1} so that a comparison of various kinetic parameters of crystallization could be made in order to understand the effect of Ag, Cd and Zn element in binary Se_{98}M_{2} alloy. Figure 2 shows typical DSC scans at heating rate of 10 K/min for glassy Se and Se_{98}M_{2} (M = Ag, Zn, Cd) alloys. Each DSC scan showed a well defined exothermic peak at the crystallization temperature T_{c}_{}. From Figure 2, it is also clear that all the chalcogenide glasses characterized in the present study show a single glass transition endothermic peak and a single exothermic crystallization peak. The single endothermic glass

Figure 1. XRD pattern of glassy Se_{98}Cd_{2} alloy.

(a)(b)(c)(d)

Figure 2. (a)-(d) DSC scans at a heating rate 10 K/min for glassy Se and binary Se_{98}M_{2} (M = Ag, Cd, Zn) alloys.

transition peaks at T_{g} and single exothermic crystallization peaks at T_{c} indicates the homogeneity of these glasses. This is an essential requirement of phase change optical recording materials. Similar scans were obtained for the other heating rates.

From the data of experiments, various kinetic parameters (T_{c}, E_{c}, n, S, H_{r}, H_{s} and K_{P}) have been determined.

4.1. Experimental Evaluation of Crystallization Temperature (T_{c})

In a DSC scan, the crystallization temperature is exhibited as a temperature corresponding to an exothermic peak due to a sudden change in specific heat. The values of T_{c} are given in Table 1. From this table, we notice a systematic shift in T_{c} to higher temperatures with an increase in the heating rate β.

It is clear from Table 1 that T_{c} of binary alloy is decreased after the incorporation of Ag and Zn but increase with Cd additive. It is well known that the crystallization temperature is associated with the nucleation and growth process that dominates the devitrification of most glassy solids. Glassy Se contains long polymeric chains and eight member rings in various proportions as its structural units, which are dissociated by the addition of Ag, Cd and Zn in glassy Se. Thus, the incorporation of Ag, Cd and Zn additives in glassy Se affects the nucleation and growth rate by cross-linking the long polymeric chains and eight member rings. This is probably the reason of change in the crystallization temperature of binary alloys at comparatively lower crystallization temperatures.

4.2. Estimation of Activation Energy of Crystallization (E_{c})

The crystallization kinetics of amorphous alloys has been extensively studied in past using the classical JohnsonMehl-Avrami (JMA) theoretical model [16-18]^{ }in which the crystallization fraction (a) can be described as a function of time (t) according to the formula:

(1)

Here, n is the Avrami index and, K is the rate constant

Table 1. Peak crystallization temperature of glassy Se_{ }and binary Se_{98}M_{2} (M = Ag, Cd, Zn) alloys at different heating rates.

which is given by:

(2)

where, E_{c} is the activation energy of crystallization, R is the universal gas constant and K_{o}_{ }is the pre-exponential factor of rate constant.

According to Kissinger [19], peak temperature of crystallization T_{c}_{ }in terms of the heating rate β, can be expressed as:

(3)

In this equation, the slope (−E_{c}/R) of the straight line for the plot of vs 1000/T_{c} is used to calculate the activation energy of crystallization.

The extent of crystallization (α) at a temperature, T is expressed as follows:

(4)

The above expression is derived by Matusita and Sakka [20,21]^{ }from the classical JMA equation. Since the value of α, is independent of β at T = T_{c} [22], so at T = T_{c}, the Equation (4) takes the form:

(5)

In this equation, the slope (−E_{c}/R) of the straight line for the plot of ln(β) vs 1000/T_{c} is used to calculate the activation energy of crystallization.

The activation energy of crystallization E_{c} can also be determined by an approximation method developed by Augis and Bennett [23]. The relation used by them is of the form:

(6)

The activation energy of crystallization has been evaluated by this equation using the plots of ln(β/T_{c})_{ } against 1000/T_{c}. This method has an extra advantage that the intercept of ln(β/T_{c})_{ }versus 1000/T_{c} gives the value of pre-exponential factor K_{0} (Equation (2)). The activation energy of crystallization for glassy Se and Se_{98}M_{2} (M = Ag, Zn, Cd) alloys has been calculated by MatusitaSakka method, Augis-Bennett method and Kissinger method. For this purpose, the plots ln(β) versus 10^{3}/T_{c},_{ }ln(β/T_{c}) versus 10^{3}/T_{c} and versus 10^{3}/T_{c} are plotted for the present glasses. Such plots for glassy Se and Se_{98}M_{2} (M = Ag, Zn, Cd) alloys are shown in Figure 3. The obtained increasing sequence, (E_{c})_{Cd} < (E_{c})_{Zn} < (E_{c})_{Ag} of activation energy of crystallization are given in Table 2. This sequence can be explained to some extent in terms of average heat of atomization. The average heat of atomization H_{S}, is based on chemical bonding aspects and H_{S} can be given as:

(7)

Here, (H_{S})_{X} and (H_{S})_{Y}_{ }are the heat of atomization of atoms X, and Y respectively and [a + b] is the weight of % of constituent elements. The plot of E_{c} vs H_{S} is shown in Figure 4 for all the three methods. This plot clearly shows that the average heat of atomization of binary Se_{98}M_{2} (M = Ag, Cd, Zn) alloys increases in the sequence (H_{S})_{Cd} < (H_{S})_{Zn} < (H_{S})_{Ag}. The values of average heat of atomization of binary Se_{98}M_{2} (M = Ag, Zn, Cd) alloys are given in Table 3. Thus, one can conclude that higher the average heat of atomization of the binary alloy, greater is the activation energy of crystallization.

4.3. Thermal Stability (S) and Glass Forming Ability (H_{r})

Saad and Poulin proposed a thermal stability parameter [24], is defined as:

(8)

where, T_{c} is the peck crystallization temperature, T_{0} is the onset crystallization temperature and T_{g} is the glass transition temperature. Thermal stability parameter reflects the resistance to devitrification after the formation of

Table 2. Activation energy of crystallization for glassy Se_{ }and binary Se_{98}M_{2} (M = Ag, Cd, Zn) alloys determined by different non-isothermal methods.

Table 3. Average heat of atomization for glassy Se and binary Se_{98}M_{2} (M = Ag, Cd, Zn) alloys.

(a)(b)(c)(d)

Figure 3. (a)-(d) Plots of ln(β) vs 10^{3}/T_{c}, ln(β/T_{c}) vs 10^{3}/T_{c} and ln (β/T_{c}^{2}) vs 10^{3}/T_{c} for glassy Se and binary Se_{98}M_{2} (M = Ag, Zn, Cd) alloys.

glass. The calculated value of stability parameter, S, is given in Table 4.

Hruby [25] introduce a value, H_{r} which is strong indicator of GFA, is defined as:

(9)

where, T_{c} is the peck crystallization temperature, T_{g} is the glass transition temperature and T_{m} is the melting temperature.

Higher the values, H_{r} signify a delay in the nucleation retardation of the growth process. The value of H_{r} is also listed in Table 4. From Table 4, one can see that the optimum values of S and H_{r} parameters are obtained in case of Ag additive.

4.4. Rate of Crystallization (K_{p})

To obtain more information about morphology of crystal growth, the following equation from Gao-Wang model^{ }[26] has been used:

(10)

Table 4. Thermal stability, S Glass forming ability, H_{r} and rate of crystallization, K_{p} for glassy Se_{ }and binary Se_{98}M_{2}_{ }(M = Ag, Cd, Zn) alloys at the heating rate 10 K·min^{-1}.

Figure 4. Plot of E_{c} vs H_{S} for glassy Se and binary Se_{98}M_{2} (M = Ag, Zn, Cd) alloys.

Here, K_{p} represents the value of rate constant at peak crystallization temperature T_{c} for a constant heating rate. The values of K_{p} for glassy Se and binary Se_{98}M_{2} (M = Ag, Zn, Cd) alloys are given in Table 4. The heating rate dependence of K_{p} is shown in Figure 5. From this figure, it is clear that the value of K_{p} is also highest in case of Ag additive.

4.5. Avrami Index (n)

For constant temperature, Equation (4) can be written as:

(11)

From this equation, the value of n can be calculated by plotting vs lnβ curves at different temperatures. The fraction “α” (extent of conversion) crystallized at any temperature T is given as α = A_{T}/A, where A is the total area of exothermic peak between the temperature T_{i} where the peak begins (i.e. the crystallization starts) and the temperature T_{f} where the peak ends (i.e. the crystallization is complete). A_{T} is the partial area of exothermic peak between the temperatures T_{i} and T.

Figure 6 shows the variation of with ln β for glassy Se and Se_{98}M_{2} (M = Ag, Zn, Cd) alloys at three constant temperatures. Using Equation (11), the value of Avrami index, n has been determined from the slope of these curves at three different temperatures. Similar plots were obtained for other glassy alloys. The temperature dependence of Avrami index for present glassy alloys is given in Table 5.

When the liquid is cooled in the glass transition region, the relaxation times for molecular movements become comparable to the experimental time scale. Therefore, the diffusive movements become comparable to the experi-mental time scale. Therefore, the diffusive motion of the liquid is trapped and the system falls out of thermal equilibrium. At this moment, the size of the nuclei does not reach the critical size required to initiate the nuclea tion process and hence the glass is assumed to have no

Figure 5. Heating rate dependence of K_{p} for glassy Se and binary glassy Se_{98}M_{2} (M = Ag, Zn, Cd) alloys.

(a)(b)(c)(d)

Figure 6. Plots of ln [ln (1 − α)^{−}^{1}] with ln β for glassy Se and binary Se_{98}M_{2} (M = Ag, Zn, Cd) alloys at different temperatures.

nuclei (of critical size). According to Matusita et al. [27], when the glass is heated in the DSC furnace, the rate of crystal nucleation in glass reaches the maximum at a temperature somewhat higher than the glass transition temperature and then decreases rapidly with increasing temperature, while the rate of crystal growth reaches the maximum at a temperature much higher than the temperature at which the nucleation rate is highest. When the glass is heated at a constant rate, the crystal nuclei are formed only at lower temperatures and crystal grow in size at higher temperatures without any increase in number.

The values of Avrami exponent n depend on two major factors. First factor is the mechanism of crystal growth which may be mono-, bior tridimensional. Second is the growth limiting factors either diffusion or motion of the outer crystal surface. The first situation corresponds to diffusion controlled growth (DCG) and the second to interface controlled growth (ICG). Finally, nucleation rate, which expresses the rate at which new nuclei are formed, may be zero, corresponding to a constant number of nuclei; it may be constant, increase or decrease versus time. Table 6 summarizes the n values for ICG and DCG mechanisms assuming growth is tridimensional [28]. The same n value may be found for different situations, and additional information is required. The mean value of n is less than 1.5 for Se while it is between 1 and 3 for binary alloys. Thus, one can conclude that the growth in the present glasses is diffusion controlled. It also supports the fact that nucleation rate becomes negligible in the region of crystallization peak.

Table 5 indicates that n decreases with increase in temperature. It is well known that crystallization of chalcogenide glasses is associated with nucleation and growth process and the extent of crystallization α increases with increase in temperature.

In other words it tends to its maximum value 1. The decrease in order parameter with increasing temperature for present glassy alloys suggests that the character of crystallization goes over from nucleation-driven in the beginning to essentially a growth driven regime by the end of crystallization process.

The values of n generally indicate the mode of crystallization. It is suggested that n = 4, 3, 2 correspond to the volume nucleation, three-dimensional, two-dimensional and one-dimensional growth, respectively and n = 1 suggests surface nucleation and one-dimensional growth. From Table 5, it is clear that the average value of n is nearly 1.5 for Se while it is nearly 1 for Se_{98}Ag_{2}, 2 for Se_{98}Zn_{2}, 3 for Se_{98}Cd_{2} alloys.

5. Conclusions

The calorimetric measurements have been performed in glassy Se and binary Se_{98}M_{2} (M = Ag, Zn, Cd) alloy us ing DSC technique. The major conclusions are:

Table 5. Temperature dependence of Avrami index, n for glassy Se_{ }and Se_{98}M_{2} (M = Ag, Zn, Cd) alloys.

Table 6. Possible values of Avrami’s exponent, n for interface-controlled growth (ICG) and diffusion controlled growth (DCG).

1) The addition of Ag, Cd, Zn additives at lower concentration (2 atomic %) in glassy Se significantly affects its thermally activated non-isothermal crystallization.

2) The maximum change in different kinetic parameters is observed for Ag additive. Thus Ag can be used for optimization of different kinetic parameters of crystallization in glassy Se as compared to Cd and Zn additives.

3) The activation energy of crystallization in binary alloys is increased with the increase in average heat of atomization.

4) The crystallization mechanism is changed from onedimensional growth to two dimensional and three dimensional growths after the inclusion of Zn and Cd additives respectively in glassy Se. However there is no effect of Ag incorporation on the crystallization growth.

6. Acknowledgements

NM is thankful to the Board of Research in Nuclear Sciences (BRNS), Mumbai, India for providing financial assistance under DAE Research Award for Young Scientists [Scheme No. 2011/20/37P/02/BRNS]. CD is thankful to UGC, New Delhi, India for providing financial support under RGNF programme to continue his research work. We also wish to express our sincere thanks to Professor O.N. Srivastava, Department of Physics, B.H.U., for undertaking the XRD measurements.

REFERENCES

- T. Ohta, “Phase-Change Optical Memory Promotes the DVD Optical Disc,” Journal of Optoelectronics and Advanced Materials, Vol. 3, No. 3, 2001, pp. 609-626.
- M. Andriesh, M. S. Lovu and S. D. Shutov, “Chalcogenide Non-Crystalline Semiconductors in Optoelectronics,” Journal of Optoelectronics and Advanced Materials, Vol. 4, No. 3, 2002, pp. 631-647.
- J. Keirsse, C. Boussard-Pledel, O. Loreal, O. Sire, B. Bureau, P. Leroyer, B. Turlin and J. Lucas, “Chalcogenide Glass Fibers Used in Bio-Sensors,” Journal of Non-Crystalline Solids, Vol. 326, 2003, pp. 430-433.
- S. Kumar and K. Singh, “Glass Transition Thermal Stability and Glass Forming Tendency of Se
_{90-x}Te_{5}Sn_{5}In_{x}Multi-Component Chalcogenide Glasses,” Thermochimica Acta, Vol. 528, 2012, pp. 32-37. doi:10.1016/j.tca.2011.11.005 - M. Popescu, “Disorder chalcogenide optoelectro-nic materials: Phenomena and Application,” Journal of Optoelectronics and Advanced Materials, Vol. 7, 2005, pp. 2189-2210.
- K. Tanaka, “Chalcogenide Glasses—In Encyclopedia of Materials: Science and Technology,” Elsevier Science Ltd., Oxford, 2001, p. 1123. doi:10.1016/B0-08-043152-6/00210-2
- Ahmad, S. A. Khan, A. A. Al-Ghamdi, F. A. Al-Agel, K. Sinha, M. Zulfequar and M. Husain, “Kinetic of NonIsothermal Crystallization of Ternary Se
_{80-x}Te_{20}Zn_{x},” Journal of Alloys and Compounds, Vol. 497, 2010, pp. 215-220. doi:10.1016/j.jallcom.2010.03.015 - S. R. Ovshinsky, “Reversible Electrical Switching Phenomena in Disordered Structures,” Physical Review Letters, Vol. 21, No. 20, 1968, pp. 1450-1453. doi:10.1103/PhysRevLett.21.1450
- S. D. Kaloshkin and I. A. Tomilin, “The Crystallization Kinetics of Amorphous Alloys,” Thermochimica Acta, Vol. 280-281, 1996, pp. 303-317. doi:10.1016/0040-6031(96)02926-7
- P. Pradeep, N. S. Saxena, M. P. Saxena and A. Kumar, “Isothermal Crystallization Study of Se
_{70}Te_{28}Cd_{2 }Chalcogenide Glasses,” Physica Scripta, Vol. 54, No. 2, 1996, pp. 207-209. doi:10.1088/0031-8949/54/2/016 - M. B. El-Den, “Study of Hardness and Crystallization Kinetics Due to Addition of Metals in SSe
_{20}Chalcogenide Glasses,” Egyptian Journal of Solids, Vol. 24, 2001, pp. 171-179. - S. P Singh, S. Kumar and A. Kumar, “Effect of Impurity (Sb and Ag) Incorporation on the a.c. Conductivity and Dielectric Properties of a-Se
_{70}Te_{30}Glassy Alloy,” Physica B: Condensed Matter, Vol. 407, No. 3, 2012, pp. 457-463. doi:10.1016/j.physb.2011.11.014 - M. A. M. Khan, S. Kumar, M. W. Khan, M. Husain and M. Zulfequar, “Electrical Transport Mechanisms in Glassy Se
_{95}M_{5}(M=Ga, Sb, Bi),” Materials Research Bulletin, Vol. 45, No. 6, 2010, pp. 727-732. doi:10.1016/j.materresbull.2010.02.011 - J. Sharma and A. Kumar, “Effect of Impurity (Sb and Ag) Incorporation on the a.c. Conductivity and Dielectric Properties of Glassy Se
_{70}Te_{30}Glassy Alloy,” Physica B: Condensed Matter, Vol. 407, 2012, pp. 457-463. doi:10.1016/j.physb.2011.11.014 - C. Dohare, N. Mehta and A. Kumar, “Effect of Some Metallic Additives (Ag, Cd, Zn) on the Crystallization Kinetics of Glassy Se
_{70}Te_{30}Alloy,” Materials Chemistry and Physics, Vol. 127, 2011, pp. 208-213 - W. A. Johnson and R. F. Mehl, “Reaction Kinetics in Processes of Nucleation and Growth,” Transactions of the American Institute of Mining and Metallurgical Engineers, Vol. 135, 1939, pp. 416-442.
- M. Avrami, “Kinetics of Phase Change I,” Journal of Physical Chemistry, Vol. 7, 1939, pp. 1103-1112. doi:10.1063/1.1750380
- M. Avrami, “Kinetics of Phase Change II,” Journal of Physical Chemistry, Vol. 8, No. 2, 1940, pp. 212-224. doi:10.1063/1.1750631
- H. E. Kissinger, “Reaction Kinetics in Differential Thermal Analysis,” Analytical Chemistry, Vol. 29, No. 11, 1957, pp. 1702-1706. doi:10.1021/ac60131a045
- K. Matusita and S. Sakka, “Kinetic Study of the Crystallization of Glass by Differential Scanning Calorimetry,” Physics and Chemistry of Glasses, Vol. 20, 1979, pp. 81-84.
- K. Matusita, S. Sakka, “Study of Non-Isothermal Crystallization of Glass by Thermal Analysis,” Bulletin of the Institute for Chemical Research, Vol. 59, No. 3, 1981, pp. 159-171.
- T. Ozawa, “Kinetic Analysis of Derivative Curves in Thermal Analysis,” Journal of Thermal Analysis, Vol. 2, No. 3, 1970, pp. 301-324. doi:10.1007/BF01690134
- J. A. Augis and J. E. Bennett, “Calculation of the Avrami Parameters for Heterogeneous Solid State Reactions Using a Modification of the Kissinger Method,” Journal of Thermal Analysis, Vol. 13, 1978, pp. 283-292.
- M. Saad and M. Poulin, “Glass Forming Ability Criteria,” Materials Science Forum, Vol. 19-20, 1987, pp. 11-18.
- Hruby, “Evaluation of Glass Forming Tendency by Means of DTA,” Czechoslovak Journal of Physics, Section B, Vol. 22, 1972, pp. 1187-1193
- Y. Q. Gao and W. Wang, “On the Activation Energy of Crystallization in Metallic Glasses,” Journal of Non-Crystalline Solids, Vol. 81, No. 1-2, 1986, pp. 129-134.
- K. Matusita, T. Konatsu and R. Yokota, “Kinetics of Non-Isothermal Crystallization Process and Activation Energy for Crystal Growth in Amorphous Materials,” Journal of Materials Science, Vol. 19, 1984, pp. 291-296.
- J. W. Christian, “The Theory of Transformation in Metals and Alloys,” 2nd Edition, Pergamon, New York, 1971.

NOTES

^{*}Corresponding author.