Open Access is an initiative that aims to make scientific research freely available to all. To date our community has made over 100 million downloads. It’s based on principles of collaboration, unobstructed discovery, and, most importantly, scientific progression. As PhD students, we found it difficult to access the research we needed, so we decided to create a new Open Access publisher that levels the playing field for scientists across the world. How? By making research easy to access, and puts the academic needs of the researchers before the business interests of publishers.
We are a community of more than 103,000 authors and editors from 3,291 institutions spanning 160 countries, including Nobel Prize winners and some of the world’s most-cited researchers. Publishing on IntechOpen allows authors to earn citations and find new collaborators, meaning more people see your work not only from your own field of study, but from other related fields too.
To purchase hard copies of this book, please contact the representative in India:
CBS Publishers & Distributors Pvt. Ltd.
www.cbspd.com
|
customercare@cbspd.com
Guilin University of Electronic Technology, P. R. China
Huaiying Zhou
Guilin University of Electronic Technology, P. R. China
*Address all correspondence to:
1. Introduction
Bulk metallic glasses (BMGs) have received a great deal of attention due to scientific and technological interest ever since the first successful synthesis of an amorphous phase in the Au–Si system in 1960 (Klement et al., 1960). There has been a lot of interest to identify parameters to assess the glass forming ability (GFA) of various alloy systems and compositions. A great deal of scientific efforts for quantification of GFA of alloys has been devoted to investigation of the GFA of alloys. There have been a lot of parameters to assess the glass forming ability (GFA) of various alloy systems and compositions. As a result, many criteria, including the confusion rule and the deep eutectic rule, for evaluating the glass forming ability (GFA) of an amorphous alloy have been proposed. Among them, the criteria used usually are the supercooled liquid region △Tx(=Tx–Tg, where Tg and Tx are the glass transition temperature and the crystallization temperature, respectively) (Inoue et al., 1993), the reduced glass transition temperature Trg(=Tg/Tl, where Tl is the liquidus temperature) (Turnbull, 1969) and the recently defined parameters γ(=Tx/(Tg+Tl)) (Lu & Liu, 2002), δ(=Τx/(Tl-Tg)) (Chen, et al., 2005), β[=TxTg/(Tl+Tx)2] (Yuan, et al., 2008), ϕ(=△Trg(Tx/Tg)0.143) (Fan, et al. 2007), ω[=Tl(Tl+Tx)/(Tx(Tl–Tx))] (Ji & Pan, 2009), γc[=(3Tx–2Tg)/Tl) (Guo, 2010), and so on. These criteria have generally proved useful parameters for evaluating the GFA of an amorphous alloy. In order to guide the design of alloy compositions with high GFA, Inoue et al. (Inoue et al., 1998) and Johnson (Johnson, 1999) have proposed the following empirical rules: (I) multicomponent systems, (II) significant atomic size ratios above 12%, (III) negative heat of mixing and (IV) deep eutectic rule based on the Trg criterion. However, Al-based metallic glasses with rare earth metal additions (Guo et al., 2000), rare earth (RE) based glasses and some binary BMGs such as Zr-Cu, Ni-Nb binary alloy (Xia et al., 2006), provide important exception from this generality, because most of above mentioned GFA parameters and rules capable of searching metallic glasses with high GFA are not applicable to these Al–based and RE-based amorphous systems. Furthermore, all the above parameters need the alloy to be first prepared in glassy form to be able to measure the crystallization temperature Tx, the liquidus temperature Tl, and/or the glass transition temperature Tg. Hence, the above parameters are not predictive in nature, as they cannot predict a good glass forming composition without actually making that alloy and rapidly solidifying it into the glassy state. It is well known that crystallization is the only event that prevents the formation of an amorphous phase. Metallic glass formation is always a competing process between the undercooled melt and the resulting crystalline phases. The GFA of a melt is thus virtually determined by the stability of the undercooled melt and the competing crystalline phases. Thermodynamic analysis could be useful in evaluating the stability of the undercooled melt and the formation enthalpies of crystalline phases. So far, several attempts have been made successfully to investigate the GFA and predict glass forming range (GFR) in several binary and ternary amorphous alloy systems, using a pure thermodynamic approach or a combined thermodynamics and kinetics approach.
From a thermodynamic point of view, there are generally following methods for calculating the GFA and predicting glass forming range (GFR) of an alloy system. The first approach is based on the T0 curve, which has been used to predict the GFR on several binary and some ternary systems. The quality of these predictions depends critically on the accuracy of the thermodynamic description. The second method is based on the semi-empirical Miedema’s model, which has been successfully applied to calculate and predict the glass forming range of some binary or ternary systems. The third consideration is directly employed on the calculation of the driving forces of crystalline phases (minimum driving force criterion) in a supercooled melt using calculation of phase diagram (CALPHAD) database. By employing driving force criterion with the obtained thermodynamic description for the investigated system, the GFA and predicted GFR of an alloy system were determined by comparing the driving force of crystalline phases precipitated from an undercooled melt. This evaluation has been successfully used to evaluate the GFA of several binary or ternary systems. Especially, it can be used to analyze the GFA of some alloy systems with unique glass forming ability, such as Al-based system. The other thermodynamic considerations, such as suppression of the formation of intermetallic phases, have been introduced.
From a combined thermodynamics and kinetics approach, the GFA of the alloys were evaluated by introducing thermodynamic quantities obtained from CALPHAD method into Davies–Uhlmann kinetic formulations. In this evaluation, by assuming homogeneous nucleation without pre-existing nuclei and following the simplest treatment based on Johnson-Mehl-Avrami’s isothermal transformation kinetics, the time–temperature-transformation (TTT) curves were obtained, which are a measure of the time t for formation of the phase Φ with a minimum detectable mass of crystal as a function of temperature. The critical cooling rates (Rc) for the glass formation calculated on the basis of the TTT curves was used to evaluate the glass-forming ability of this binary or ternary alloy. The calculated GFA results show good agreement with the experimental data in the compositional glass formation range of the investigated systems.
This chapter is intended to present systematically the methods and progress on the glass forming ability investigated by a thermodynamic approach or a combined thermodynamics and kinetics approach.
2. Calculation of GFA based on thermodynamics analysis
Usually, it is regarded the formation of metallic glasses is controlled by two factors, i.e., the cooling rate and the composition of the alloy. The critical cooling, which is the most effective gauge for GFA of the alloys, is hard to be measured experimentally. Hence, a great deal of efforts has devoted to the investigation on the correlation between the GFA and the composition of glass forming alloys. Inoue et al. (Inoue et al., 1998) and Johnson (Johnson, 1999) proposed the empirical rules to predict the element selection and compositional range of glass forming alloy. These rules have played an important role as a guideline for synthesis of BMGs for the last decade. However, recent experimental results have shown that the “confuse principal” and “deep eutectic rule” cannot be applicable to the Cu-Zr, Ni-Nb binary system (Xia et al., 2006) and Al-based ternary system (Guo et al., 2000). From a thermodynamic point of view, it is well known that crystallization is the only event that prevents the formation of an amorphous phase. During a melt–quenching process for metallic glass formation, the glass formation is exposed to crystallization competition of other crystalline phases from the undercooled melt between liquidus temperature Tl and glass transition Tg. The GFA of a melt is thus virtually determined by the stability of the undercooled melt and the competing crystalline phases, which can be analyzed by thermodynamic analysis. In this section, several GFA calculation based on thermodynamics analysis were introduced.
2.1. Calculation of the GFA of alloys based on T0 curve
2.1.1. Method
Generally, a glass can be formed during cooling when crystallization is avoided up to the occurrence of the glass transition. Thus, in order to predict the tendency to glass formation in a system and the composition regions where it is most probable, nucleation of crystals in the undercooled melt must be considered. The GFR will be the region of composition where nucleation of crystalline phases is less likely. Various models have been developed to analyze the GFA of alloys in the literature, as will be discussed in the following section, with different levels of approximation. T0 curve is one of the approaches used to estimate the GFA of the alloys.
A T0 curve is the locus of the compositions and temperatures where the free energies of two phases are equal. Thus, T0 curves can be calculated provided that their Gibbs free energy is known, i.e. an assessment of the system is available. The T0 curve between the liquid and a solid phase determines the minimum undercooling of the liquid for the partitionless formation of a crystalline solid with the same composition (Boettinger & Perepezko, 1993). Fig. 1 showed one example for a simple eutectic system. Alloys with T0 curves plunge steadily at low temperatures (dashed line in Fig. 1a), there will be no driving force for partitionless transformation in the composition region between them. If the equilibrium crystalline phases are not prone to nucleation, the glass can thus form. On the contrary, if T0 curves that are only slightly depressed below the stable liquidus curves are good candidate for partionless transformation of crystalline phases in the entire composition range (dashed line in Fig. 1b).
Figure 1.
Hypothetical T0 curves for a binary eutectic A–B system. (a) T0 curves drop to low temperature: glass formation is possible. (b) T0 curves intersect at low temperature: partitionless crystalline phase formation occurs (redrawn from Boettinger & Perepezko, 1993).
2.1.2. Application of T0 curve
Predictions of GFR based on T0 curves have been performed on several binary and some ternary systems. The construction of T0 curves for alloy glass needs a precise knowledge of thermodynamic properties of the supercooled liquid alloy and the introduction of the transition to the glassy state (Kim, et al., 1998). The quality of these predictions depends critically on the accuracy of the thermodynamic description and the introduction of the excess specific heat contribution is expected to improve the quality of results (Palumbo & Battezzati, 2008). However, as pointed out by Schwarz and co-workers (Schwarz et al., 1987), some discrepancies have been observed between the prediction and experimental results. For example, even when using the most recent thermodynamic assessment (Kumar, 1996) to calculate T0 curves in the Cu–Ti system, the results are not agreement with the reported experimental GFR. In fact, T0 curves for terminal solid solutions do not plunge at low temperatures as expected for glass forming systems (Kumar et al. 1996). Battezzati and co-workers (Battezzati, et al., 1990) have shown that in the Cu–Ti system the contribution of the excess specific heat is essential for describing the glass forming ability. An excess specific heat contribution has also been considered in the Al–Ti system (Cocco, et al., 1990) and the Fe–B system (Palumbo, et al., 2001).
2.2. Calculation of the GFA of alloys based on Miedema’s model
2.2.1. Method
Miedema’s model is an empirical theory for calculating heat of mixing in various binary systems both for the solid state (Miedema et al., 1975) and liquid (Boom et al., 1976). This model involves the calculations of the formation enthalpy of metallic glasses (amorphous phase) (ΔHamor), solid solutions (ΔHSS), and intermetallic compounds (ΔHinter) according to the following equations (Bakker 1988; Boer et al. 1988).
ΔHamor=ΔHchem(amor)+ΔHtopoE1
ΔHss=ΔHchem(SS)+ΔHelastic+ΔHstructureE2
and
ΔHinter=ΔHchem(inter)E3
where ΔHchem(amor) is the chemical mixing enthalpy of the amorphous state, ΔHtopo is the topology enthalpy of a glass, ΔHchem(SS) is the chemical mixing enthalpy of a solid solution, ΔHelastic is the elastic enthalpy of the solid solution calculated based on the continuous elastic model proposed by Friedel (Friedel, 1954) and Eshelby (Eshelby, 1954& 1956), ΔHstructure is the structure enthalpy induced by the structural changes, and ΔHchem(inter) is the chemical mixing enthalpy of an intermetallic compound. The formation enthalpy ΔHinter of a composition between two adjacent intermetallic compounds can be calculated using the level principle.
The chemical contribution of enthalpy of mixing of solid solution can be written as
ΔHchem=xAxB[xAΔHBinASS+xBΔHAinBSS]E4
where xA and xB represent the mole fraction of A and B atoms and ΔHSS is the enthalpy of solution of one element in another at infinite dilution. The data have been taken from Niessen et al. (Niessen, et al., 1983).
The elastic term in the enthalpy of formation originates from the atomic size mismatch, which can be expressed as
ΔHelastic=xAxB[xAΔHBinAelastic+xBΔHAinBelastic]E5
TheΔHiinjelastichas been obtained by using the formalism by Simozar and Alonso(Simozar & Alonso, 1984) as
ΔHiinjelastic=2μj(Vi−Vj)2Vj(3+4μjKi)E6
where μj is the shear modulus of the solvent, Vi and Vj are the molar volumes of the solute and the solvent, respectively and Ki is the compressibility of the solute.
The structural contribution of enthalpy for solid solution originates from the valence and the crystal structure of the solute and the solvent atom. It is found to have a very minor contribution and it is difficult to calculate. Hence, the structural contribution to enthalpy has been usually neglected (Basu, et al., 2008). In the case of the elastic and structural contributions are absent, thus the formation enthalpy of glasses can be calculated as
ΔHamor=ΔHchem(amor)+3.5∑1nxiTm,iE7
Where xi represents the mole fraction of component i atom, Tm,i is the melting temperature of the component i.
According to the Miedema’s model, an amorphous phase can be formed if the enthalpy of formation of the amorphous phase is less than that of the solid solution phase. The heat of formation in alloys generally arises from the interactions among the constituent atoms where the interfacial energy plays a major role. The interfacial energy mainly comes from the atomic size difference. It has also been postulated that the number of intermetallic phases appearing in an alloy system is a strong function of the heat of mixing. The number of intermetallic phase in an alloy system increases with the increase in the heat of mixing. This model can be directly used to determine the glass forming range in binary alloy systems and can be extended to ternary systems by neglecting the ternary interactions.
2.2.2. Calculation of the GFA for the binary alloy systems
Since the metallic glass formation process is controlled by thermodynamic factors, Miedema’s model was firstly used to predict the composition range of amorphous binary transition metal alloys (Kolk et al., 1988; Coehoorn et al., 1988; Murty, et al., 1992; Basu, et al., 2008). It is found that the predicted glass forming composition ranges are in good agreement with the experimental results. In the work of Takeuchi and Inoue (Takeuchi & Inoue, 2000), this approach has been used to calculate the mixing enthalpy and mismatch entropy of a number of bulk metallic glass alloy systems. It has been observed that the mixing enthalpy and normalised mismatch entropy for glass forming alloys vary within a certain range.
Figure 2.
Enthalpy–composition curves for binary Ti–Ni, Zr–Ni, Hf–Ni, Ti-Cu alloy systems (a–d). The curve with (•) and the curve with (Δ) represent amorphous and solid solution phase, respectively. The enthalpy values are in J/mol (From Basu, et al., 2008).
As shown in Fig. 2, in the work of Basu et al. (Basu, et al., 2008), glass forming range (GFR) has been determined for different binary (Ti–Ni, Zr–Ni, Hf–Ni, Ti–Cu, Zr–Cu, Hf–Cu) in (Zr, Ti, Hf)–(Cu, Ni) alloys based on the mixing enthalpy and mismatch entropy calculations. Though copper and nickel appear next to each other in the periodic table, the glass forming ability of the copper and nickel bearing alloys is different. Thermodynamic analysis reveals that the glass forming behaviour of Zr and Hf is similar, whereas it is different from that of Ti. The smaller atomic size of Ti and the difference in the heat of mixing of Ti, Zr, Hf with Cu and Ni leads to the observed changes in the glass forming behaviour. Enthalpy contour plots can be used to distinguish the glass forming compositions on the basis of the increasing negative enthalpy of the composition. This method reveals the high glass forming ability of binary Zr–Cu, Hf–Cu, Hf–Ni systems over a narrow composition.
In the recent work performed by Xia (Xia, et al., 2006), the GFA of an alloy is considered that the formation of the meta-stable amorphous state should include two aspects: (1) the driving force for the glass formation, i.e., −ΔHamor, and (2) the resistance of glass formation against crystallization, i.e. the difference between the driving force for glass phase and for the intermetallic compound formation ΔHamor−ΔHinter. When two glass forming alloys have the same −ΔHamor but different ΔHamor−ΔHinter, their GFA can then be dominated by ΔHamor−ΔHinter. The lower the value of ΔHamor−ΔHinter, the higher the GFA of the alloy. On the other hand, when two glass forming alloys have the same ΔHamor−ΔHinter but different ΔHamor, their GFA is dominated by −ΔHamor. The higher the value of −ΔHamor, the better the GFA. Since the contribution from entropies is much smaller as compared with that from the formation enthalpy of solid compounds (Delamare, et al., 1994), the GFA is expressed in terms of formation enthalpy alone. Based on this thermodynamic consideration, a new parameter γ* to evaluate GFA for glass formation was proposed by Xia et al. (Xia, et al., 2006) and expressed as
γ*=GFA∝ΔHamorΔHinter−ΔHamorE8
where ΔHamor and ΔHinter are the enthalpies for glass and intermetallic formation, respectively. Both ΔHamor and ΔHinter are calculated by Miedema’s macroscopic atom model.
This parameter has been successfully used to predict the GFR and the best GFA alloy compositions in Zr-Cu and Ni-Nb system by comparing the value of γ* of various alloy systems, respectively (Xia, et al., 2006).
Figure 3.
Calculated dependence of the parameter γ* on Zr and Ni concentration in Cu-Zr (a) and Ni-Nb (b) binary alloys, respectively (from Xia, et al., 2006).
Fig. 3 shows the calculated dependence of the parameter γ* on Zr and Ni concentration in Cu-Zr (a) and Ni-Nb (b) binary alloys, respectively, suggesting that the alloys Cu64Zr36 and Cu50Zr50 in Cu-Zr system, and Ni61.5Nb38.5 in Ni-Nb system are the best glass former, respectively. These predicted results are in good agreement with the experimentally reported Cu64.5Zr35.5 and Cu50Zr50, and Ni62Ni38 that could be made into bulk metallic glass rods with 2 mm in diameter, indicating that γ* is an effective parameter in identifying the best glass former in the Zr-Cu and Ni-Nb binary system.
Similarly, considering both the stability of liquid employing ΔHliq/ΔHinter, and the competition of glass and crystal using ΔHamor/ΔHinter, Ji et al. (Ji, et al., 2009) proposed a new parameter γ’ of GFA as
γ'=GFA∝ΔHliq•ΔHamor(ΔHinter)2E9
As Ji et al. described, this parameter γ’ is not only verified in five different binary bulk metallic glasses (Cu–Hf, Ni–Nb, Cu–Zr, Ca–Al, Pd–Si) but also showed wider application range comparing with the former model, but also have a better GFA estimation on the different composition than the parameter γ* because it is including ΔHliq in the evaluation expression. The predicated results are in good agreement with the experiments in all five different kinds of binary BMG systems and the biggest deviation of the peak of γ’ from the best current GFA composition is only about 6 at.% in Ca–Al alloy. Comparing with former GFA parameter γ*, γ’ takes account of liquid stability and shows more universal for evaluation GFA in different kinds of binary alloys (Ji, et al., 2009). Recently, Wang et al. also made a modification to Xia’s proposal and it works more convenient to describe the GFA of transition metal systems (Wang, et al., 2009).
2.2.3. Calculation of the GFA for the multicomponent alloy systems
Miedema’s approach has been extensively used by Nagarajan and Ranganathan (Nagarajan & Ranganathan, 1994), Takeuchi and Inoue (Takeuchi & Inoue, 2001&2004) and other researchers (Murty, et al., 1992; Rao, et al.; 2007; Basu, et al, 2008; Wang & Liu, 2009; Sun, et al., 2010) to determine the glass forming composition range (GFR) in a number of ternary and multicomponent systems. In the work performed by Takeuchi and Inoue (Takeuchi & Inoue, 2001), the amorphous-forming composition range (GFR) was calculated for 338 ternary amorphous alloy systems on the basis of the database given by Miedema's model in order to examine the applicability of the model, to analyze the stability of the amorphous phase, and to determine the dominant factors influencing the ability to form an amorphous phase. The mixing enthalpies of amorphous and solid solution phases were expressed as a function of alloy compositions on the basis of chemical enthalpy. The GFR was calculated for 335 systems except for the Al-Cu-Fe, Al-Mo-Si and Au-Ge-Si systems. The calculated results are in agreement with the experimental data for Cu-Ni- and Al-Ti-based systems. For typical amorphous alloy systems exemplified by the Zr-, La-, Fe- and Mg-based systems, it was recognized that the calculated GFR had been overestimated as a result of the model being simplified. It is found that the elastic enthalpy term arising in a solid solution phase stabilizes the amorphous phase, and the stabilization mechanism is particularly notable in Mg-based amorphous alloy systems. Short-range order plays an important role in the formation of Al-, Fe- and Pd-metalloid based systems (Takeuchi & Inoue, 2001).
Based on Miedema’s model and Alonso’s method, the glass forming ability/range (GFA/ GFR) of the Fe–Zr–Cu system was studied by thermodynamic calculation. It is found that when the atomic concentration of Zr is between 34% and 56%, no matter what the atomic concentrations of Fe and Cu are, amorphous phase could be obtained, thus the atomic mismatch playing a dominating role in influencing the GFA. While the atomic concentration of Zr is out of the above range, the GFA is highly influenced by the immiscibility between Fe and Cu (Wang & Liu, 2009).
Glass forming composition range for ternary Zr–Ti–Ni, Zr–Hf–Ni, Ti–Hf–Ni, Zr–Ti–Cu, Zr–Hf–Cu and Ti–Hf–Cu systems has been determined by extending the Miedema’s model to ternary alloy systems and by neglecting the ternary interaction parameter (Basu, et al., 2008). In their calculations, solid pure metals have been chosen to be the standard state and their enthalpy has been assigned to be zero. It is seen that the glass forming composition range for Ni bearing alloys is higher than that of the Cu bearing alloys, as heat of mixing of Ni is higher than that of Cu with Ti, Zr and Hf. In these ternary (Zr, Ti, Hf)–(Cu, Ni) alloys mixing enthalpy and mismatch entropy varies between (−13) and (−42) kJ/mol and 0.13 and 0.25, which is within the range predicted for glass formation (Basu, et al., 2008).
In the work of Oliveira et al., the γ* parameter proposed by Xia et al. was extended to the ternary Al–Ni–Y system. The calculated γ* isocontours in the ternary diagram are compared with experimental results of glass formation in that system. Despite some misfitting, the best glass formers are found quite close to the highest γ* values, leading to the conclusion that this thermodynamic approach can be extended to ternary systems, serving as a useful tool for the development of new glass-forming compositions (Oliveira et al., 2008).
Rao et al. (Rao et al. 2007) identified the composition with highest glass forming ability in Zr-Ti-Ni-Cu-Al quinary systems with the Gibbs-energy change between the amorphous and solid solution phases as the thermodynamic parameter by calculating the Gibbs-energy change with the help of Miedema, Miracle, mismatch entropy and configurational entropy models. ΔG shows the strong correlation with the reduced glass transition temperature (Tg/Tl) in Zr-based metallic glasses. Thus, ΔG can be used as a predictive GFA parameter to identify compositions with the highest GFA. The compositions with the highest GFA have been identified in a number of quinary systems by iso-free energy contour maps by representing quandary systems as quasiternary systems (Rao et al. 2007). The best glass forming composition has been identified by drawing iso-Gibbs-energy change contours by representing quinary systems as pseudo-ternary ones. Attempts have been made to correlate the Gibbs-energy change with different existing glass forming criteria and it is found that the present thermodynamic parameter has good correlation with the reduced glass transition temperature. Further, encouraging correlations have been obtained between the energy required for amorphization during mechanical alloying to the Gibbs-energy change between the amorphous and solid solutions.
2.3. Calculation of the GFA of alloys based on driving force criterion
2.3.1. Method
The basic underlying concept to predict the compositions of alloys having high GFA using the thermodynamic approach is that the compositions exhibiting the local melting minimum points favour amorphous phase formation. Thermodynamic approach of driving force criterion is based on a different concept. During a melt–quenching process for metallic glass formation, the glass formation is exposed to crystallization competition of other crystalline phases from the undercooled melt between liquidus temperature Tl and glass transition Tg. It is well known that the crystallization is the only event that prevents the formation of amorphous phase. Considering that crystallization is usually through the nucleation and growth process, the high GFA can be inversely predicted by searching a condition where the nucleation and growth of crystalline phases can be retarded. There are three dominating factors for the kinetics, (i) chemical driving force, (ii) interfacial energy, as an energy barrier, between the amorphous phase and the crystalline phases, (iii) the atomic mobility for rearrangement or transport of the partitioning atoms. According to the classical nucleation theory, the driving force of formation of the crystalline phases and interfacial energy, among other things, affects the nucleation rate of product phases. The interfacial energy between liquid and crystalline phases is known to be small compared to surface energy or grain boundary energy (Porter & Easterling, 1992), and therefore the role of interfacial energy in the nucleation kinetics of crystalline phases would be small. Then, the driving force of formation becomes the major factor that affect the nucleation kinetics of crystalline phases from amorphous alloy melts. It is believed that alloys with lower driving force for the formation of crystalline phases under the supercooled liquid state suggest higher GFA in the glass forming range. Therefore, Kim and co-workers proposed the minimum driving force criterion as a new thermodynamic calculation scheme to evaluate the composition dependence of the GFA (Kim et al., 2004). The driving force for the crystalline phases can be calculated using the critical assessed thermodynamic parameters by the CALPHAD method (Kaufman & Bernstein, 1970). In the CALPAHD method, the Gibbs energies of individual phases are described using thermodynamic models. Then, the model parameters are optimized considering relevant experimental information on phase equilibria or the other thermodynamic properties. The calculation of phase equilibrium is performed based on the minimum Gibbs energy criterion.
2.3.2. Application of the driving force criterion
Driving force criterion has been successfully used to explain the composition dependence of GFA in several glass forming alloys with unique GFA, such as Cu-Zr-Ti (Kim, et al., 2004), Mg-Cu-Y (Kim, et al., 2005), Al-Ce-Ni (Tang, et al. 2010), and Al-Cu-Zr (Bo, et al, 2010) systems, by calculating the driving force of formation of crystalline phases under metastable supercooled liquid states and by searching the local minima of the driving forces for crystallization. The calculated results are in good agreement with the experimental results. It has been indicated that the driving force criterion can be used as a new thermodynamic scheme to estimate the composition dependence of GFA in multicompoent alloy systems for the development of bulk amorphous alloys.
Figure 4.
Calculated driving forces of crystalline phases for Cu55Zr45–xTix alloys, versus Ti content at (a) 1073 K, (b) 973 K, and (c) 873 K (from Kim, et al., 2004).
For the Cu-Zr-Ti system, among a series of ternary alloys Cu60Zr40–xTix (x = 10, 20, 30), the alloy with the highest GFA should be the alloy Cu60Zr20Ti20 according to the maximum Trg criterion (Turnbull, 1969), while experiments (Inoue, et al., 2001) show it is Cu60Zr30Ti10. Although the other alloys, Cu55Ti35Zr10 (Lin & Johnson, 1995) and Cu47Ti33Zr11Ni8Si1 (Choi, et al., 1998) based on the Cu–Ti–Zr ternary system but at different region, have been published as alloys with high GFA, there is no empirical rule or factor that can explain why the high GFA is obtained at certain compositions in the Cu–Ti–Zr system (roughly Zr: Ti=3:1 and Zr: Ti=1:3). As already the thermodynamic parameters for all phases in this system obtained by the CALPHAD method, the GFA is estimated by calculating the driving forces of all crystalline phases under the undercooled liquid state. Fig. 4 shows the calculated driving forces of individual crystalline phases as a function of Ti content in a temperature range (600–800 oC) where the alloys correspond to supercooled liquids state. As shown in this figure, along the composition line Cu55Zr45–xTix with varying Ti content, the driving forces of crystalline phases show two local minimums, one at Zr-rich region (x = 7–10) and the other at Ti-rich region (x= 28–29). According to the driving force criterion, the two local minimum points in Fig. 4 are the compositions where the GFA is expected to be higher than other compositional region. In a sense that the Zr:Ti ratios in the two local minimum points are toughly close to 3:1 and 1:3, it can be said that the former is close to the Inoue’s composition and the latter is close to Johnson’s composition (Kim et al., 2004). This finding indicates that the composition dependency of the GFA in the Cu-Zr-Ti ternary alloy system can be explained by calculating the driving forces of formation of crystalline phases under metastable supercooled liquid states and by searching the local minima of the driving forces for crystallization (Kim, et al. 2004).
Similarly, Al–based amorphous, which was discovered in 1988 (He, et al., 1988& Inoue, et al., 1988), is also of particular interest because of its low density, good bending ductility and high tensile strength. It was found that, however, most of above mentioned parameters and rules capable of searching metallic glasses with high GFA are not applicable to Al–based amorphous (Guo, et al., 2000, Hackenberg, et al., 2002, Gao, et al., 2003, Zhu, et al. 2004).Al–Ce–Ni system is an unique Al–based system, which can be synthesized into a strong, flexible metallic glass with the widest GFR covering 2–15 at.% Ce and 1–30 at.% Ni (Inoue, 1998, Kawazoe et al., 1997). The alloys with high GFA are situated away from the eutectic point. Experimental results of the Al–Ce–Ni bulk amorphous alloys prepared with copper mold casting indicate that the amorphous sheets with 5 mm width and 0.2 mm thickness are obtained in Al86Ce4Ni10 and Al88Ce6Ni6 alloys without appreciable glass transition. On contrary, alloys Al82Ce8Ni10 and Al80Ce6Ni14 with △Tx values of 20 and 21K consist mainly of crystalline phases (Inoue, 1998). After a thermodynamic assessment of the Al-Ce-Ni system in the Al-rich corner was performed, a set of consistent thermodynamic parameters were obtained, and the thermodynamic properties of the Al-Ce-Ni amorphous alloys were calculated. The calculated results indicated that the alloys with high GFA in the Al–Ce–Ni system are far from the eutectic point, and the heats of mixing are from –15 to –49 kJ/mol of atom for the observed amorphous alloys (Tang et al., 2010).
As shown in Fig. 5, the relatively smaller nucleation driving forces for the formation of crystalline phases for the Al–10Ce based alloys (Fig. 5a) are generally indicative of their higher GFA with a reportedly wider GFR (1–30at.% Ni) (Kawazoe et al., 1997). In contrast, the relatively larger driving forces in the Al–10Ni based alloys (Fig. 7b) are associated with their poorer GFA and narrower GFR (2–10 at.% Ce) (Kawazoe et al., 1997). This finding is further confirmed by the melt spinning (Tang, et al., 2010) and the copper mold casting experimental results (Inoue, 1998).
Based on the experimental enthalpies of mixing of ternary liquid and undercooled liquid alloys as well as the evaluated isothermal sections, the Al–Cu–Zr ternary system has been assessed using the CALPHAD method. Most of the calculated results show good agreement with the experimental thermodynamic data and the reported phase diagrams. By employing the driving force criterion with the present thermodynamic description, the observed glass-forming ability in the Al–Cu–Zr system can be accounted for satisfactorily (Bo, et al., 2010).
Figure 5.
Calculated normalized nucleation driving force (per mole of atoms) for crystalline phases from undercooled (a) Al–10Ce–Ni, (b) Al–Ce–10Ni and (c) Al-6Ce-Ni metastable liquid at 800 oC (from Tang, et al., 2010).
2.4. Calculation of the GFA of alloys based on the other thermodynamic approach
By treating the glass transition as a second-order phase transformation from liquid phase (Palumbo, et al., 2001, Shao et al., 2005), which can give good predictions of all important GFA indicators such as the reduced glass transition temperature and the thermodynamic stability of the amorphous phase, Shao et al. established a full thermodynamic database for glass forming ability (GFA). The resultant thermodynamic database can be used to produce all major temperature-related GFA indicators such as Tg/Tl, Tg/Tm and Tx/(Tg+Tl). It is indicated that together with phase diagram prediction, such an extensive CALPHAD approach is a powerful tool for designing alloys with large GFA (Shao et al., 2005).
By using the computational thermodynamic approach exhibiting low-lying liquidus surfaces coupled with the reduced glass transition temperature criterion of Turnbull, regions of alloy composition suitable for experimental tests for glass formation of Zr–Ti–Ni–Cu–Al system were identify rapidly by Cao et al. (Cao, et al., 2006). The glass forming ability of the alloys we studied can be understood in terms of the relative liquidus temperature in a thermodynamically calculated temperature vs. composition section through a multicomponent phase diagram. It does not follow several other proposed thermodynamic or topological criteria.
A thermodynamic parameter (ΔHchem ×Sσ/kB) in the configuration entropy (Sconfig/R) range of 0.8-1.0 has been developed to identify excellent BMG composition using enthalpy of chemical mixing (ΔHchem), the mismatch entropy normalized by Boltzmann’s constant (Sσ/kB) and the configurational entropy (Sconfig/R) by Bhatt et al. and it has been demonstrated for the Zr-Cu-Al based ternary system. It is found that this approach can be used to predict the best BMG composition more closely than the earlier models (Bhatt, et al., 2007).
Based on the undercooling theory resulting from the existence of multicomponent chemical short-range order (CSRO) domains, the glass forming range (GFR) in Zr-Ni-Ti alloy system was predicted by thermodynamic calculation. The GFR predicted by the thermodynamic calculation is consistent with the experiment results (Liu, et al. 2008).
One of the ways to predict the possible bulk glass formation composition is the phase diagram calculation with suppression of the formation of intermetallic phases. The formation of stoichiometric intermetallic compounds which have the ordered structure of atoms into specific lattice sites can take a time for the rearrangement of atoms from liquid state. Thus, the formation of intermetallic compounds can be suppressed during the fast solidification process normally applied to the bulk glass production. Combining the obtained the thermodynamic database and the above concept, the amorphous formation diagram of the Cu–Zr–Ag system with the suppression of all binary and ternary intermetallic phases has been proposed by Kang and Jung (Kang & Jung, 2010)
3. Calculation of the GFA of alloys based on a combined thermodynamics and kinetics approach
As discussed above, the thermodynamic approach is useful since the thermodynamic parameters can be used to calculate the GFR in binary alloys and can also be used to predict the GFR in ternary systems based on the constituent binaries. One of the limitations of a purely thermodynamic approach is that it does not give the critical cooling rates for the glass formation. A combined thermodynamic and kinetic treatment, based on time-temperature-transformation curves (TTT) in the manner of Uhlmann and Davies has been presented (Saunders & Miodownik, 1986& 1988). This combined approach takes the thermodynamic parameters obtained from the phase diagram calculations and derives values for the free energy barrier for nucleation, free energy driving forces, and melting points used in kinetic equations. The combined approach has been successfully used to calculate the glass forming ability (GFA) of a wide range of binary and ternary alloy systems (Saunders & Miodownik, 1988). The calculated glass forming ranges for a wide number of binary and ternary alloy systems are in good agreement with experiment. A significant advantage of the combined approach is that data from binary alloy systems, often with little or no ternary modification, can be used to calculate the necessary thermodynamic input for the kinetic equations in higher order systems. This section outlines briefly the combined thermodynamic and kinetics method used for the calculation of the GFA for alloy systems.
3.1. Method
Critical cooling rates for glass formation can be obtained by Johnson-Mehl-Avrami isothermal transformation kinetics using the equation
X=1−exp[(−π/3)IvUc3t4]E10
where X is the volume fraction of material transformed, Iv is the nucleation frequency, Uc is the crystal growth rate, and t is the time taken to transform X. In the early stages of transformation the value of X approximates to
X≅πIvUc3t4E11
For homogeneous nucleation without pre-existing nuclei, the nucleation frequencyIvh is given by
Ivh=DnNva02exp(−ΔG*/kT)E12
where Dn is the diffusion coefficient necessary for crystallisation, Nv is the number of atoms per unit volume. a0 is an atomic diameter. k is Boltzmann’s constant. T is the transformation temperature, and ΔG*is the free energy barrier for nucleation of a spherical nucleus given by the expression
ΔG*=16π3(σ3/Gv2)E13
where σ is the liquid/crystal interfacial energy and Gv is the change in free energy per unit volume on solidification. An equation for Uc can be written as
Uc=fDga0[1−exp(−ΔGm/RT)]E14
where Dg is the diffusion coefficient for the atomic motion necessary for liquid to crystal growth, ΔGmis the molar free energy driving force for liquid to crystal growth, and R is the universal gas constant, and f is a structural constant denoted the fraction of sites on the interface where atoms may preferentially be added or removed, and is given by the following expression (Uhlmann, 1972)
f=0.2(T−mT)/TmE15
where Tm is the liquidus temperature. By assuming that Dn=Dg=the bulk liquid diffusivity and invoking the Stokes-Einstein relationship between diffusivity and viscosity η, equation (12) and (14) can be derived to give the time t needed to form a volume fraction X of transformed crystalline phase in an undercooled liquid as following
where t is the time taken to transformation volume fraction X of crystalline solid. η is the viscosity of liquid, a0 is an atomic diameter, f is a structural constant, Nv is the number of atoms per unit volume, ΔG*is the Gibbs energy barrier to nucleation and ΔGm is the Gibbs energy driving force for the liquid-crystal transformation. The constants have been typically taken as X=10-6, a0=0.28×10-9m, f=0.1 and Nv=5×1028 atoms/m3. In order to apply this equation to a real alloy system, it is necessary to derive or estimate the parameters η,ΔG*, and ΔGm (Saunders & Miodownik, 1988).
3.2. Estimation of η,ΔG*, and ΔGm
Since it is very difficult to measure experimentally the viscosity of supercooled liquid, there have been few measurements of it. In this case, the viscosity can be generally described as being between the liquidus temperature Tm and the glass transition Tg using a Doolittle-type expression involving the relative free volume fT (Ramachandrarao, et al., 1977) as
η=Aexp(B/fT)E17
where
fT=Cexp(−EH/RT)E18
EH is the hole formation energy and A, B, and C are constants. Because of the lack of experimental data and the EH value was estimated by means of a direct relationship from Tg (Ramachandrarao, et al., 1977). Assuming B is unity and fT and η are 0.03 and 1012Ns/m2, respectively at Tg, A and C have been approximated at 3.33×10-3 and 10.1, respectively. If Tg values are not available, crystallisation temperatures Tx are used as a first approximation.
Figure 6.
The construction used in calculating the driving force, ΔGm, for the crystallization of compound AB2 from a liquid of composition x1 in the A-B system.
For the crystallization of compounds AB2 from a liquid of composition x1 in the A-B system (Fig. 6), the molar free energy driving force for liquid to crystal growth, ΔGm, represents the Gibbs energy required to form one mole of crystalline phase from the liquid of composition x1, which can be obtained from thermodynamic phase diagram calculations to give explicitly molar heats of fusion Hmf and driving force that can be used toΔG*. Values of Hmfand ΔGm are calculated from partial molar Gibbs energies of elements A and B and free energy values. Therefore, ΔGmis expressed by the following equation
ΔGm=xAG¯AL+xBG¯BL−GcrystE19
where xA and xB are the mole fractions of elements A and B in the precipitating crystalline phase, respectively. G¯AL and G¯BL are the partial molar Gibbs energies of elements A and B in the liquid phase, respectively, and Gcryst is the integral free energy of the precipitating crystalline phase. The Gibbs energy functions in Eq. (19) can be obtained from the thermodynamic model parameters evaluated in the literature. In an A-B alloy system, a liquid composition x1 becomes unstable with respect to the compounds AB2 at the liquidus temperature Tm. At a given temperature T1, there is a driving force for the precipitation of the compound AB2 given by G1 (Fig. 6), where G1 is defined as the driving force to form one mole of compound AB2 in a liquid of composition x1. In all cases here ΔGm is equal to ΔG1. By using heats of formation in place of free energy values, Hmfcan be similarly evaluated.
The Gibbs energy barrier to nucleation of a spherical nucleus ΔG* can be described as
ΔG*=16π3N(σm3/ΔGm2)E20
where N is Avogadro’s number and the σm the molar liquid/crystal interfacial energy. σm is directly related to the molar enthalpy of fusion Hmf and expressed as
σm=αHmfE21
where α is a proportional constant. Hmfcan be obtained in a similar way to evaluate ΔGm based on bond energy values across the interface (Turnbull, 1950). Saunders and Miodownik empirically evaluated the constant α to be 0.41 (Saunders & Miodownik, 1988).
3.3. Calculation of critical cooling rates below T0 of disordered solid phases
The expression for t in equation (16) is derived assuming that the kinetics of the liquid to crystal transformation are limited by the bulk diffusivity, which is appropriate when the crystal composition differs from that of the liquid, or at compound compositions where substantial diffusion is necessary before the correct spatial relationships that define the ordered structure of the compound are achieved. However, at the temperatures below the T0 temperature of a disordered solid solution phase, the liquid becomes unstable with respect to a molecularly simple phase of the same composition. Consequently, no long range diffusion is necessary for the liquid to crystal and the kinetics are governed by atom motions of less than one atom in diameter. Then, the transformation is considered to be extremely difficult to suppress and this forms the T0 criteria for GFA. In such cases, it has been suggested that the rate limiting step for crystal growth is proportional to the rate at which atoms collide at the liquid/crystal interface, and an expression for the crystal growth rate is then given (Boettinger et al., 1984) by
Uc=fV0[1−exp(−ΔGm/RT)]E22
where V0 is the velocity of sound in the liquid metal. This is the same form as equation (14), but with V0 replacing the Dg/a0. Replacing the Dn/a0 in equation (12) with V0 and rearranging equation (12) and (14), an expression for t is derived as
t≈1V0{Xa0πf3Nvexp(ΔG*/kT)[1−exp(−ΔGm/RT)]3}1/4E23
The value for V0 has been taken as 1000 m/s by Saunders and Miodownik (Saunders & Miodownik, 1988), close to the a value used by Boettinger et al. (Boettinger et al., 1984) and no transformation is considered to occur below Tg.
From equations (10) to (23), the time-temperature-transformation (TTT) curve can be obtained. The critical cooling rate Rc necessary for amorphous phase formation with a melt quenching method can be evaluated from TTT curve calculated and approximated as follows
Rc=Tm−Tn5tnE24
where Τm and tn are the temperature and time at the nose of the TTT curve, respectively, since the cooling rate calculated directly from the isothermal transformation curve is somewhat overestimated compared with that from the CCT (continuous cooling transformation) curve, the right side of Equation (24) has divided by a factor of 5 to emulate continuous cooling. In the composition range with Rc˂1×10-7 K/s, which has been generally known to be a maximum available cooling rate for melting quenching, the amorphous phase formation may be possible.
3.4. Evaluation of glass forming ranges in alloy systems
The combined thermodynamic and kinetic approach has been undertaken to evaluate the GFA of a wide number of binary and ternary alloy systems since the pioneering work performed by Sanders and Miodownik (Saunders & Miodownik, 1988; Shim et al., 1999; Clavaguera-Mora, 1995; Tokunaga, et al., 2004; Abe, et al., 2006; Ge, et al., 2008; Palumbo, & Battezzati, 2008; Mishra & Dubey, 2009). They calculated the free energy driving forces from the thermodynamic databases, free energy barrier for nucleation and melting points, and employed this data to kinetic calculation. There is in good agreement between the predicted glass forming ranges and those experimentally observed. It is indicated that the approach has the potential to predict glass forming ability in multicomponent alloys using mainly binary input data.
The first attempts to couple kinetic models with reliable thermodynamic data using the CALPHAD methodology was performed by Saunders and Miodownik (Saunders & Miodownik, 1988). In his work, the combined thermodynamics and kinetics approach was presented in detail and undertaken to evaluate the GFA of a wide range of binary (Au-Si, Pd-Si, Ti-Be, Zr-Be, Hf-Be, Cu-Ti, Co-Zr, Ni-Zr, Cu-Zr, Ni-P, Pd-P) and ternary (Ni-Pd-P, Cu-Pd-P, Co-Ti-Zr, Zr-Be-Hf, Ti-Be-Hf) alloy systems (Saunders & Miodownik, 1988). The TTT curves and the critical cooling rate for glass formation Rc were estimated. There is excellent agreement between the predicted and observed GFRs of binary systems, apart from the discrepancies in the Ti-Be and Cu-Ti systems. The approach was then extended to give predications for critical cooling rates in ternary and multicomponent alloys using mainly binary information. The results would appear to indicate that the combined approach takes into account a number of the major effects that govern glass formation and has the potential to predict GFA in multicomponent systems (Saunders & Miodownik, 1988).
In the work performed by Ge et al., (Ge, et al., 2008), the glass forming ability (GFA) of nine compositions of Cu-Zr and thirteen of Cu-Zr-Ti alloys in terms of critical cooling rate and fragility were evaluated by combining CALPHAD technique with kinetic approach. The driving forces for crystallization from the undercooled liquid alloys were calculated by using Turnbull and Thompson-Spaepen (TS) Gibbs free energy approximate equations, respectively. As shown in Fig. 7, time-temperature-transformation (TTT) curves of these alloys were obtained with Davies-Uhlmann kinetic equations based on classical nucleation theory. With Turnbull and TS equations, the critical cooling rates are calculated to be in the range of 9.78 ×103-8.23×105 K/s and 4.32 ×102-3.63×104 K/s, respectively, for Cu-Zr alloys, and 1.38×102-7.34×105 K/s and 0.64-1.36×104 K/s, respectively, for Cu-Zr-Ti alloys (Ge, et al., 2008).
Figure 7.
Calculated TTT curve of Cu-Zr (left) and Cu-Zr-Ti (right) alloys by (a) Turnbull model and (b) TS model (from Ref. Ge, et al., 2008).
Based on topological, kinetic and thermodynamic considerations, Yang et al. (Yang, et al., 2010) have discussed the existence of the multiple maxima in GFA in a single eutectic system in Al–Zr–Ni system. It is apparent that, when taken alone, none of the factors seemed to be able to fully explain this phenomenon we have observed. It is suggested that glass formation is an intricate balance of kinetic, thermodynamic and also topological factors. Perhaps in good glass formers, all factors could come to a consensus at one composition or one compositional zone, where the best glass former(s) are located. However, for marginal glass formers like Al-based alloys, each of these factors could point to a different alloy composition, where conditions are best suited for glass formation.
Recently, Considering chemical short-range ordering and metastability of undercooled melts, Zhu and co-workers have applied a simplified quasi-kinetic approach in order to predict the GFR in binary Al-rare earth (Zhu, et al., 2004) and Al-based Al-Gd-Ni(Fe) ternary system (Zhu, et al., 2004), using CALPHAD databases. They derived the derived an expression for the reduced time t’ = t/tmin for the formation of a minimal quantity of crystalline solid, where t represents the composition-dependent time needed for the transformation and tmin the minimum transformation time at a certain optimum composition:
t'≈(exp(ΔG*/kT)[1−exp(−ΔGm/RT)]3)1/4E25
This formula is in fact equivalent to Eqs. (10) to (23), except that the effect of parameters related to atomic transport is neglected. The calculated reduced times for various solid crystalline phases are then used to predict the GFR, i.e. the region of the composition space where these times are higher. A qualitative satisfactory agreement can be observed (Zhu, et al., 2004). The ability to predict the GFR of candidate metallic glass systems indicates a simple but effective approach for reducing reliance on extensive experimental trial and error in the search for new metallic glass systems (Zhu, et al., 2004).
Search for new bulk metallic glasses (BMGs) system or composition by predicting the GFA of an alloy system is of interesting and theoretical and practical significance. In this chapter, the progress on the calculation or predication the glass forming ability by thermodynamics approach or a combined thermodynamics and kinetics approach have been reviewed. It is found that a good agreement between the predicated glass forming ability and those experimentally observed has been obtained. It is indicated that the thermodynamic approach developed in the literature has proved useful to predict the glass forming ability of a number of alloys system. It has revealed that the combined thermodynamics and kinetics approach has the advantage to predict the glass forming ability of the multicomponent alloys using the reliable database of binary system assessed by CALPHAD method. It has been accepted that the thermodynamic approach and/or the combined thermodynamic and kinetic approach are effective ways for the prediction of the GFA of metallic glass alloys.
References
1.AbeT.ShimonoM.OdeM.OnoderaH.2006Thermodynamic modeling of the undercooled liquid in the Cu-Zr system. Acta Materialia, 544 (February 2006), 909915 , 1359-6454
2.BakkerH.1988Enthalpies in Alloys: Miedema’s Semi-Empirical Model. Trans Tech, 978-0-87849-783-6 Uetikon-Zurich, Switzerland
3.BasuJ. MurtyB. S.RanganathanS.2008Glass forming ability: Miedema approach to (Zr, Ti, Hf)-(Cu, Ni) binary and ternary alloys. Journal of Alloys and Compounds. 4651/2 , (October 2008), 163172 , 0925-8388
4.BattezzatiL.BariccoM.RiontinoG.SolettaI.1990Thermodynamic evaluation of the Cu-Ti system in view of solid state amorphization reactions. Journal de Physique Colloques. 51C4 (1990), 7985 , 0449-1947
5.BhattJ.JiangW.XiaJ.WangQ.ChuangD.MurtyB. S.2007Optimization of bulk metallic glass forming compositions in Zr-Cu-Al system by thermodynamic modeling. Intermetallics, 155/6, (May/June 2007), 716721 , 0966-9795
6.BoH.WangJ.JinS.QiH. Y.YuanX. L.LiuL. B.JinZ. P.2010Thermodynamic analysis of the Al-Cu-Zr bulk metallic glass system. Intermetallics, 1812 (December 2010), 23222327 , 0966-9795
7.BoomR. De BoerM. R.MiedemaA. R.1976On the heat of mixing of liquid alloys.Journal of Less-Common Metetals, 452 (March 1976), 237245 , 0925-8388
8.BoettingerW. J.PerepezkoJ. H.1993 in: Rapidly Solidified Alloys, Liebermann, H. H., (Ed.), Dekker M., 17978-0-82478-951-2 New York,
9.Clavaguera-MoraM. T.1995 Glass materials: thermodynamic and kinetic quantities. Journal of Alloys and Compounds, 2201/2, (April 1995), 197205 , 0925-8388
10.CoccoG.SolettaI.BattezzatiL.BariccoM.EnzoS.1990Mechanical alloying of the Al-Ti system. Philosophical Magazine B, 614 (April 1990), 473479 , 1364-2812
11.CoehoornR.Van der KolkG. J. Van denBroek. J. J. MinemuraT.MiedemaA. R.1988Thermodynamics of the stability of amorphous alloys of two transition metals.Journal of Less-Common Metetals, 1401/2, (June 1988), 307316 , 0925-8388
12.ChenQ. J.ShenJ.FanH. B. SunJ. F. HuangY. J.MccartneyD. G.2005Glass-forming ability of an Iron-based alloy enhanced by Co addition and evaluated by a new criterion. Chinese Physics Letters, 227 (July 2005), 17361738 , 0025-6307X
13.ChoiY. H. BuschR.JohnsonW. L.1998The effect of silicon on the glass forming ability of the Cu47Ti34Zr11Ni8 bulk metallic glass forming alloy during processing of composites. Journal of Applied Physics, 8312 (June 1998), 79937997 . 0021-8979
14.DelamareJ. LemarchandD.VigierP.1994, Structural investigation of the metastable compound A1 in an as-cast Fe---Nd eutectic alloy. Journal of Alloys and Compouds, 2162 (January 1994), 273280 , 0925-8388
15.De BoerF. R. BoomR. MattensW. C. M. MiedemaA. R.NiessenA. K.1988 Cohesion in Metals: Transition Metals Alloys. North-Holland, 978-0-44487-098-8 Amsterdam, The Netherlands
16.EshelbyJ. D.1954 Distortion of a Crystal by Point Imperfections. Journal of Applied Physics, 252 (February 1954), 255261 . 0021-8979
17.EshelbyJ. D.1956 The continuum theory of lattice defects. Solid State Physics : Advances in Research and Applications. Frederick Seitz and David Turnbull, (Ed.), 3 (1956), 79144 , Elsevier, 978-0-12374-292-6 Amsterdam, the Netherlands
18.FriedelJ.1954Electronic structure of primary solid solutions in metals. Advances in Physics, 312 (October 1954), 446507 , 0001-8732
19.FanG. J.ChooH.LiawP. K.2007 A new criterion for the glass-forming ability of liquids. Journal of Non-Crystalline Solids, 3531 (January 2007), 102107, 0022-3093
20.GaoM. C.HackenbergR. E.ShifletG. J.2003Journal of Alloys and Compounds, 3531/2 , (April 2003), 114123 , 0925-8388
21.GeL.HuiX.WangE. R.ChenG. L.ArroyaveR.LiuZ. K.2008Prediction of the glass forming ability in Cu-Zr binary and Cu-Zr-Ti ternary alloys. Intermetallics, 161 (January 2008), 2733 , 0966-9795
22.GuoF. Q.PoonS. J.ShifletG. J.2000Investigation of glass formability in Al-based multinary alloys. Scripta Materialia, 4312 (November 2000), 10891095 , 1359-6462
30.InoueA.OhteraK.TsaiA. P.MasumotoT.1988Aluminum-based amorphous alloys with tensile strength above 980 MPa (100 kg/mm2). Japanese Journal of Applied Physics, 274 (April 1988), L479L482 , 0021-4922
31.InoueA.1998 Amorphous, nanoquasicrystalline and nanocrystalline alloys in Al-based systems. Progress in Materials Science, 435 (December 1998), 365520 , 0079-6425
32.JiX. L.PanN. Y.2009A thermodynamic approach to assess glass-forming ability of bulk metallic glasses. Transaction of Nonferrous Metals Society of China, 195 October 2009), 12711279 , 1003-6326
33.JiX. L.PanY.NiF.2009A thermodynamic criterion for predicting glass-forming ability in binary metallic glasses. Materials and Design, 303 (March 2009), 842845 , 0261-3069
35.KaufmanL.BernsteinH.1970 Computer calculation of phase diagrams. Academic Press; 978-1-84569-215-5 NY, USA
36.KangD. H.JungI. H.2010 Critical thermodynamic evaluation and optimization of the Ag-Zr, Cu-Zr and Ag-Cu-Zr systems and its applications to amorphous Cu-Zr-Ag alloys. Intermetallics, 185 (May 2010), 815833 , 0966-9795
37.KawazoeY.MasumotoT.SuzukiK.InoueA.TsaiA. P.YuJ. Z.AiharaJ. T.NakanomyoT.1997Nonequilibrium phase diagrams of ternary amorphous alloys. In: Landolt- Bönstein, New Series III: Condensed Matter, Y. Kawazoe, J. Z. Yu, A. P. Tsai & T. Masumoto, (Ed.), 37A112978-3-54065-664-7
38.KimD.LeeB. J.KimN. J.2004Thermodynamic approach for predicting the glass forming ability of amorphous alloys. Intermetallics. 1210-11 , (October-November 2004), 11031107 , 0966-9795
39.KimD. LeeB. J.KimN. J.2005 Prediction of composition dependency of glass forming ability of Mg-Cu-Y alloys by thermodynamic approach. Scripta Materialia, 52 (2005), 969972 , 1359-6462
40.KimY. K.SohJ. R.KimH. S.LeeH. M.1998Thermodynamic prediction of glass forming range in Al-Mg-REM ternary system. CALPHAD. 222 (June 1998). 2212301 , 0364-5916
41.KumarK. C. H.AnsaraI.WollantsDelaeyP.L.1996Thermodynamic optimisation of the Cu-Ti system. Zeitschrift für Metallkunde. 878 (August 1996), 666672 , 0004-4093
43.LinX. H.JohnsonW. L.1995Formation of Ti-Zr-Cu-Ni bulk metallic glasses.Journal of Applied Physics, 7811 (December 1995), 65146519 . 0021-8979
44.LuZ. P.LiuC. T.2002A new glass-forming ability criterion for bulk metallic glasses. Acta Materialia, 5013 (August 2002), 35013512 , 1359-6454
45.LiuX. J.HuiX. D.ChenG. L.2008Thermodynamic calculation and experimental investigation of glass formation in Zr-Ni-Ti alloy system. Intermetallics, 162 (February 2008), 262266 , 0966-9795
46.MiedemaA. R. BoomR.De BoerM. R.1975On the heat of formation of solid alloys. Journal of Less-Common Metals, 412 July 1975), 283298 , 0925-8388
47.MishraR. K. DubeyS.2009Glass forming ability of materials: A thermodynamic approach. Journal of Non-Crystalline Solids. 35545-47 (November 2009), 21992204 , 0022-3093
48.MurtyB. S.RanganathanS.RaoM. M.1992Solid state amorphization in binary Ti-Ni, Ti-Cu and ternary Ti---Ni---Cu system by mechanical alloying. Materials Science and Engeering A, 1492 (January 1992), 231240 , 0921-5093
49.NagarajanR.RanganathanS.1994A study of the glass-forming range in the ternary Ti-Ni-Al system by mechanical alloying. Materials Science and Engeering A, 179/180, 1 (May 1994), 168172 , 0921-5093
50.NiessenA. K. BoerF. R.de .BoomR. ChatelP. F.de .MatternW. C. M.MiedemaA. R.Model predictions for the enthalpy of formation of transition metal alloys II. CALPHAD, 771 (Januay-March 1983), 5170 , 0364-5916
51.OliveiraM. F.AliagaL. C. BolfariniR. C. BottaW. J.KiminamiC. S.2008Thermodynamic and topological instability approaches for forecasting glass-forming ability in the ternary Al-Ni-Y system. Journal of Alloys and Compounds, 4641/2, (September 2008), 118121 , 0925-8388
52.PorterD. A.EasterlingK. E.1992Phase transformation in metals and alloys. Chapman & Hall, 978-0-74875-741-1 London, UK
53.PalumboM.CacciamaniG.BoscoE.BariccoM.2001Thermodynamic analysis of glass formation in Fe-B system, CALPHAD, 254 (December, 2001), 625637 . 0364-5916
54.PalumboM.BattezzatiL.2008Thermodynamics and kinetics of metallic amorphous phases in the framework of the CALPHAD approach. CALPHA, 322 (June 2008), 295314 , 0364-5916
55.RaoB. S.BhattJ.MurtyB. S.2007Identification of compositions with highest glass forming ability in multicomponent systems by thermodynamic and topological approaches. Materials Science and Engineering A, 449-451 , 1 (March 2007), 211214 , 0921-5093
56.RamachandraraoP.CantorB.CahnR. W.1977Free volume theories of the glass transition and the special case of metallic glasses. Journal of Materials Science, 1212 (December 1977), 24882502 , 0022-2461
57.SaundersN.MiodownikA. P.1988 Evaluation of glass forming ability in binary and ternary metallic alloy systems- an application of thermodynamic phase diagram calculation. Materials Science and Technology, 49 (September 1988), 768777 , 0267-0836
58.SaundersN.MiodownikA. P.1986Thermodynamic aspects of amorphous phase formation. Journal of Materials Research, 11 (January 1986), 3846 . 0884-2914
59.SaundersN.MiodownikA. P.1985The use of free energy vs composition curves in the prediction of phase formation in codeposited alloy thin films. CALPHAD, 93 (September 1985), 283290 , 0364-5916
60.SchwarzR. B.NashP.TurnbullD.1987The use of thermodynamic models in the prediction of the glass-forming range of binary alloys. Journal of Materials Research. 24 (July 1987), 456460 , 0884-2914
62.ShimJ. H.ChungH.J.LeeD. N.1999Calculation of phase equilibria and evaluation of glass-forming ability of Ni-P alloys.Journal of Alloys and Compounds, 2821/2, (January 1999), 145181 , 0925-8388
63.SimozarS.AlonsoJ. A. On the Factors Controlling Glass Forming Ability of Metallic Alloys Formed by Fast Liquid Quenching. Physica Status Solidi (a), 8111 (January 1984), 5561 , 1862-6319
64.SunS. P.YiD. Q.LiuH. Q.ZangB. JiangY.2010Calculation of glass forming ranges in Al-Ni-RE (Ce, La, Y) ternary alloys and their sub-binaries based on Miedema’s model. Journal of Alloys and Compounds, 5061 (September, 2010) 377387 , 0925-8388
65.TakeuchiA.InoueA.2000Calculations of Mixing Enthalpy and Mismatch Entropy for Ternary Amorphous Alloys. Materials Transactions JIM, 4111 (November 2000), 13721378 , 0916-1821
66.TakeuchiA.InoueA.2001 Calculations of Amorphous-Forming Composition Range for Ternary Alloy Systems and Analyses of Stabilization of Amorphous Phase and Amorphous-Forming Ability. Materials Transactions JIM , 427 July 2001) 1435 EOF1444, 0916-1821
67.TakeuchiA.InoueA.2001Thermodinamic and kinetic analyses of formation of amorphous and nanocrystalline alloys with the aid of computer and database.Journal of Optoelectronics and Advanced Materials, 62 (June 2004), 533540 , 1454-4164
68.TangC. Y.DuY.WangJ.ZhouH. Y.ZhangL. J.ZhengF.LeeJ.YaoQ.2010Correlation between the thermodynamics and the glass forming ability in the Al-Ce-Ni system, Intermetallics, 185 (May 2010), 900906 , 0966-9795
69.TokunagaT.OhtaniH.HasebeM.2004Thermodynamic evaluation of the phase equilibria and glass-forming ability of the Fe-Si-B system. CALPHAD,284 (December 2004), 354362 . 0364-5916
70.TurnbullD.1950 Formation of crystal nuclei in liquid metals. Journal of Applied Physics, 2111 (November 1950), 10221028 , 0021-8979
71.TurnbullD.1969 Under what conditions can a glass be formed? Contemporary Physics, 105 (October 1969), 473488 , 0010-7514
72.UhlmannD. R.1972A kinetic treatment of glass formation. Journal of Non-Crystalline Solids, 74 (May 1972), 337348 , 0022-3093
73.Van der KolkG. J. MiedemaA. R.NiessenA. K.1988 On the composition range of amorphous binary transition metal alloys. Journal of Less-Common Metetals, 1451/2, (December 1988), 117 , 0925-8388
74.WangT. L.WangW. C.LiJ. H.LiuB. X.2010Glass formation of the Fe-Hf system studied by thermodynamic calculation and ion beam mixing. Journal of Alloys and Compounds, 4931/2 . (March 2010), 154157 , 0925-8388
75.WangT. L.LiJ. H.LiuB. X.2009Proposed thermodynamic method to predict the glass formation of the ternary transition metal systems. Physical Chemistry Chemical Physics, 1114 (February 2009), 23712373 , 1463-9076
76.WangT. L.LiuB. X.2009Glass forming ability of the Fe-Zr-Cu system studied by thermodynamic calculation and ion beam mixing. Journal of Alloys and Compounds, 4811/2, (July 2009), 156160 , 0925-8388
77.XiaL.FangS. S.WangQ.DongY. D.LiuC. T.2006Thermodynamic modeling of glass formation in metallic glasses. Applied Physic Letters, 8817 (April 2006), 171905171905 -3, 0003-6951
79.YangH.LimK. Y.LiY.2010Multiple maxima in glass-forming ability in Al-Zr-Ni system. Journal of Alloys and Compounds, 4891 (January 2010), 183187 , 0925-8388
80.YuanZ. Z.BaoS. L. LuY.ZhangD. P.YaoL.2008A new criterion for evaluating the glass-forming ability of bulk glass forming alloys. Journal of Alloys and Compouds, 4591/2, ( July 2008), 251260 , 0925-8388
81.ZhuA.ShifletG. J.MiracleD. B.2004Glass forming ranges of Al-rare earth metal alloys: thermodynamic and kinetic analysis. Scripta Materialia, 507 (April 2004), 987991 , 1359-6462
82.ZhuA.PoonS. J.ShifletG. J.2004On glass formability of Al-Gd-Ni (Fe). Scripta Materialia, 5012 (June 2004), 14511455 , 1359-6462
Written By
Chengying Tang and Huaiying Zhou
Submitted: 10 November 2010Published: 15 September 2011