Thermodynamic Optimization of the Ternary Ga-Sn-Te System Using Modiﬁed Quasichemical Model

: Thermoelectric (TE) materials are of great interest to many researchers because they directly convert electric and thermal energy in a solid state. Various materials such as chalcogenides, clathrates, skutterudites, eutectic alloys, and intermetallic alloys have been explored for TE applications. The Ga-Sn-Te system exhibits promising potential as an alternative to the lead telluride (PbTe) based alloys, which are harmful to environments because of Pb toxicity. Therefore, in this study, thermodynamic optimization and critical evaluation of binary Ga-Sn, binary Sn-Te, and ternary Ga-Sn-Te systems have been carried out over the whole composition range from room temperature to above liquidus temperature using the CALPHAD method. It is observed that Sn-Te and Ga-Te liquids show the strong negative deviation from the ideal solution behavior. In contrast, the Ga-Sn liquid solution has a positive mixing enthalpy. These different thermodynamic properties of liquid solution were explicitly described using Modiﬁed Quasichemical Model (MQM) in the pair approximation. The asymmetry of ternary liquid solution in the Ga-Sn-Te system was considered by adopting the toop-like interpolation method based on the intrinsic property of each binary. The solid phase of SnTe was optimized using Compound Energy Formalism (CEF) to explain the high temperature homogeneity range, whereas solid solution, Body-Centered Tetragonal (BCT) was optimized using a regular solution model. Thermodynamic properties and phase diagram in the Ga-Sn-Te and its sub-systems were reproduced successfully by the optimized model parameters. Using the developed database, we also suggested several ternary eutectic compositions for designing TE alloy with improved properties.


Introduction
Thermoelectric (TE) materials are being continuously developed for application in waste heat recovery and electronic cooling because of their ability to interchange between thermal energy and electrical energy in a solid-state [1,2].As these materials can harvest waste energy and directly convert it to electricity, they can be applied as green and sustainable energy sources.However, the low conversion efficiency of the TE devices is one of the major limiting factors for commercial applications.This conversion efficiency is described by a dimensionless thermoelectric figure of merit (ZT) represented as S 2 σT/(k), where S, σ, k, and T are Seebeck coefficient, electrical conductivity, thermal conductivity, and temperature, respectively.The Seebeck coefficient of a material is represented as the ratio of an induced thermoelectric voltage to the temperature difference across that material [3].Generally, the Seebeck coefficient is inversely proportional to the number of charge carriers, whereas electrical conductivity is directly proportional to the number of charge carriers [4,5].Therefore, it is challenging to maximize the collective term, S 2 σ.Another approach is to reduce thermal conductivity by introducing phonon scattering sites such as solid solutions, multi-phase microstructure, precipitations, or nanostructuring.entropy at 298.15 K.As the thermodynamic principle, the Gibbs free energies of pure elements can be calculated as follows: where H • 298.15K , S • 298.15K , and C p are standard enthalpy of formation from the stable pure elements at 298.15 K, entropy at 298.15 K, and heat capacity, respectively.The heat capacity expression of each compound can be determined by fitting heat capacity data.The Neumann-Kopp rule [33] can be used to predict the C p expression in case of no available reliable data on C p .

Liquid Solution
The Modified Quasichemical Model (MQM) [27,28], which accounts for the shortrange ordering of the nearest neighbor atoms, was used to describe the liquid solutions.Compared to the conventional Bragg-Williams Random Mixing Model, the MQM provides a more realistic description of the entropy of the solution.In the MQM, Gibbs energy of a pair formation can be expressed in a polynomial of a pair fraction instead of a component fraction.Besides, the coordination number is allowed to vary with the composition.These added features of the MQM provide greater flexibility in reproducing the binary experimental data and combining optimized binary liquid parameters into a larger database for a multicomponent system.
In the MQM, the reaction of pair exchanging in a binary A-B liquid solution can be expressed by the distribution of A and B atoms over the sites of the quasi lattice as follows: (A-A) + (B-B)= 2(A-B); ∆g AB (4) where (i-j) represents the First-Nearest Neighbor (FNN) pair between components i and j, and ∆g AB is the Gibbs energy change of forming two moles of (A-B) pairs.The Gibbs energy of the liquid solution can be estimated by the following equation: where n A , n B , and n AB are the moles of A, B atoms, and A-B pair, respectively.g • A and g • B are the molar Gibbs energies of the pure A and B components.∆S config is the configurational entropy of randomly mixing (A-A), (B-B), and (A-B) pairs.
where n ij is the number of moles of (i-j) pairs.Z A and Z B are the coordination numbers of A and B, respectively.Pair fraction X ii , mole fraction X i , and coordination equivalent fraction Y i are defined as follows: Equation ( 10) can be obtained by mass balance of (i-j) pairs.The ∆g AB is the model parameter and can be expanded as a polynomial in terms of the pair fraction as follows:

∆g
• AB , ∆g i0 AB , and ∆g 0j AB are the model parameters, and that can be a function of temperature.The composition dependent coordination numbers are given as follows: where Z A AA and Z A AB are the value Z A when all nearest-neighbor atoms of atom A are A and Bs, respectively.Z B BB and Z B BA are defined analogously.The composition of maximum SRO in each binary subsystem is estimated by ratio of the coordination numbers Z B /Z A .In this study, all the coordination numbers are set to 6, except Z Te GaTe = 4 to have a maximum SRO composition, X Te = 0.6 in the Ga-Te binary system [29].
The Gibbs energy of the ternary liquid solution can be predicted using the interpolation technique based on the model parameter of the sub-binary systems.One of the advantages of the MQM is the flexibility to choose the interpolation method depending on the nature of each binary system.Toop-like interpolation with Te as an asymmetric component was used to describe the asymmetry of the ternary liquid with the positive enthalpy of mixing in the Ga-Sn system and the strong negative enthalpy of mixing in the Ga-Te and Sn-Te systems.

Solid Solution
In the Ga-Sn-Te ternary system, there are only two existing solid solution phases, SnTe in Sn-Te binary and Body-Centered Tetragonal (BCT) in the Ga-Sn binary system.The Compound Energy Formalism (CEF) [34] is used to describe the Gibbs energy of SnTe to explain the high temperature homogeneity range of the compound.The homogeneity range of the SnTe phase is described by a large intrinsic vacancy at high temperature; thus, it was formulated as (Sn, Va) (Te).The Gibbs energy of the solid solution can be expressed as: G sol. s = y Sn y Te G Sn:Te + y Va y Te G Va:Te + RT y Sn lny Sn + y Va lny Va + ∑ i,j,k y i y j y k L i,j:k + ∑ i,j,k y k y i y j L k:i,j where y i and y i are the site fraction of species, i, in the corresponding sublattice.G i:j is the Gibbs free energy of an end-member, and the L i,j:k and L k:i,j are adjustable interaction parameters between the species in related sublattices.BCT solid solution was optimized using one lattice regular solution model.This solution model is one lattice version of the CEF.The Gibbs energy expression of the solution using this model is as follows: where L i Sn,Ga is the adjustable model parameter.The self-consistent thermodynamic parameter from the present optimization is given in Table 1.

Ga-Sn Binary System
The Ga-Sn system has been thermodynamically assessed by Anderson and Ansara [35].They proposed the Ga-Sn phase diagram based on the extensive review of the thermodynamic and phase diagram data.Recently, thermodynamic optimization of this system was done by Kulawik et al. [36].They used a regular solution model with exponential model parameters for the liquid solution.Their proposed thermodynamic description can satisfactorily represent available thermodynamic and phase diagram data.However, to keep the consistency with other binaries, Ga-Sn liquid solution was optimized using MQM in this study.According to the previous optimization, the Ga-Sn is a simple binary eutectic system that includes three phases: Othorhombic_Ga, BCT_A5, and liquid.The liquidus line was determined by Delcroix et al. [37], Predel [38], and Puschin et al. [39] using thermal analysis and by Shurai et al. [40], Trebukhov et al. [41], and Zivkovic et al. [42] using differential thermal analysis (DTA).The solidus line in Sn rich side was determined by Predel [38] using X-ray and thermal analysis.By measuring electrical resistance, Audit et al. [43] and Bosio et al. [44] have reported extrapolated Sn solubility in Ga at the eutectic temperature 0.016 at.% Sn and 0.027 at.% Sn, respectively.This reported solubility of Sn in Ga is very small.Therefore, it has not been included in the present optimization.The BCT phase was described using one lattice regular solution model.No thermodynamic property was reported for the BCT solid solution.Therefore, model parameters were obtained to fix the solidus data reported by Predel [38].The calculated phase diagram of the Ga-Sn system is shown in Figure 1 along with experimental data [37][38][39][40][41][42].Using MQM, only one parameter was required to describe the thermodynamic properties of the liquid solution.
This reported solubility of Sn in Ga is very small.Therefore, it has not been included in the present optimization.The BCT phase was described using one lattice regular solution model.No thermodynamic property was reported for the BCT solid solution.Therefore, model parameters were obtained to fix the solidus data reported by Predel [38].The calculated phase diagram of the Ga-Sn system is shown in Figure 1 along with experimental data [37][38][39][40][41][42].Using MQM, only one parameter was required to describe the thermodynamic properties of the liquid solution.Danilin and Yatsenko [45] and Katayama et al. [46] determined activity in the liquid solution using the Electromotive Force (EMF) method at 1063 K and 1073 K, respectively.Zivkovic et al. [42] determined the activity in the liquid solution using DTA.Mixing enthalpy of liquid Ga-Sn alloy was measured by Bros et al. [47] at 742 K, Bros and Laffitte [48] at 513 K, Gambino et al. [49] at 723 K, and Predel et al. [50] at 623 K using calorimetry.Skoropanov et al. [51] and Zivkovic et al. [42] determined the mixing enthalpy of liquid using DTA at 743 K and 1000-1200 K, respectively.The calculated activity and mixing enthalpy of liquid solution are shown in Figure 2a,b, respectively.These calculated thermodynamic properties with a single positive parameter are in good agreement with the respective experimental data [42,[45][46][47][48][49][50][51].Danilin and Yatsenko [45] and Katayama et al. [46] determined activity in the liquid solution using the Electromotive Force (EMF) method at 1063 K and 1073 K, respectively.Zivkovic et al. [42] determined the activity in the liquid solution using DTA.Mixing enthalpy of liquid Ga-Sn alloy was measured by Bros et al. [47] at 742 K, Bros and Laffitte [48] at 513 K, Gambino et al. [49] at 723 K, and Predel et al. [50] at 623 K using calorimetry.Skoropanov et al. [51] and Zivkovic et al. [42] determined the mixing enthalpy of liquid using DTA at 743 K and 1000-1200 K, respectively.The calculated activity and mixing enthalpy of liquid solution are shown in Figure 2a,b, respectively.These calculated thermodynamic properties with a single positive parameter are in good agreement with the respective experimental data [42,[45][46][47][48][49][50][51].
This reported solubility of Sn in Ga is very small.Therefore, it has not been included in the present optimization.The BCT phase was described using one lattice regular solution model.No thermodynamic property was reported for the BCT solid solution.Therefore, model parameters were obtained to fix the solidus data reported by Predel [38].The calculated phase diagram of the Ga-Sn system is shown in Figure 1 along with experimental data [37][38][39][40][41][42].Using MQM, only one parameter was required to describe the thermodynamic properties of the liquid solution.Danilin and Yatsenko [45] and Katayama et al. [46] determined activity in the liquid solution using the Electromotive Force (EMF) method at 1063 K and 1073 K, respectively.Zivkovic et al. [42] determined the activity in the liquid solution using DTA.Mixing enthalpy of liquid Ga-Sn alloy was measured by Bros et al. [47] at 742 K, Bros and Laffitte [48] at 513 K, Gambino et al. [49] at 723 K, and Predel et al. [50] at 623 K using calorimetry.Skoropanov et al. [51] and Zivkovic et al. [42] determined the mixing enthalpy of liquid using DTA at 743 K and 1000-1200 K, respectively.The calculated activity and mixing enthalpy of liquid solution are shown in Figure 2a,b, respectively.These calculated thermodynamic properties with a single positive parameter are in good agreement with the respective experimental data [42,[45][46][47][48][49][50][51].

Sn-Te Binary System
The Sn-Te phase diagram consists of two eutectic reactions and one intermediate congruently melting phase, SnTe.This binary system has been thermodynamically assessed by various authors [52- the model parameters for this binary system.Biltz and Mecklenburg [58], Kobayashi [59], Umeda et al. [60], and Harries et al. [61] used thermal analysis to measure the liquidus temperature.Brebrick and Strauss [62] measured the partial pressure of Te 2 over SnTe and determined the liquidus temperature near the SnTe region.Rakotomavo et al. [63] and Le Bouteiller et al. [64] determined the liquidus temperature using EMF measurements and DTA analysis.Rakotomavo et al. [63] reported Eutectic temperature 507±1 K on Sn rich side and 673 ± 1 K on Te rich side based upon DTA measurements.Le Bouteiller et al. [64] reported eutectic temperature 675 ± 2 K in Te rich side.Kurosawa et al. [65] measured the electrical conductivity of the liquid Sn-Te system to determine the liquidus.Brebrick [66,67] determined the homogeneity range of the intermetallic compound SnTe by measuring the partial pressure of Te 2 over SnTe and lattice parameter measurements.Shelimova and Abrikosov [68], Sha and Brebrick [69], and Rogacheva et al. [70,71] determined the homogeneity range by hardness, partial pressure measurements, and lattice parameter as well as hardness measurements, respectively.The calculated phase diagram of the Sn-Te system is shown in Figure 3a along with experimental data [58][59][60][61][62][63][64][65] and the most recent optimization [56].The calculated enlarged homogeneity range of the SnTe compound is shown in Figure 3b along with experimental data [66][67][68][69][70][71].As seen in Figure 3b, compared to the previous calculation by Liu et al. [56], the calculated homogeneity range of the SnTe is in good agreement with the experimental data.
All proposed thermodynamic descriptions can satisfactorily reproduce reported experimental thermodynamic and phase diagram data in this system.However, Hsieh et al. [52] and Kattner et al. [53] used old thermodynamic descriptions of the pure elements, and Kuznetsov [54] and Sharma and Chang [57] obtained model parameters that were not consistent with other CALPHAD type modeling.They have extensively reviewed the available thermodynamic data and those data were employed in the present study to determine the model parameters for this binary system.Biltz and Mecklenburg [58], Kobayashi [59], Umeda et al. [60], and Harries et al. [61] used thermal analysis to measure the liquidus temperature.Brebrick and Strauss [62] measured the partial pressure of Te2 over SnTe and determined the liquidus temperature near the SnTe region.Rakotomavo et al. [63] and Le Bouteiller et al. [64] determined the liquidus temperature using EMF measurements and DTA analysis.Rakotomavo et al. [63] reported Eutectic temperature 507±1 K on Sn rich side and 673 ± 1 K on Te rich side based upon DTA measurements.Le Bouteiller et al. [64] reported eutectic temperature 675 ± 2 K in Te rich side.Kurosawa et al. [65] measured the electrical conductivity of the liquid Sn-Te system to determine the liquidus.Brebrick [66,67] determined the homogeneity range of the intermetallic compound SnTe by measuring the partial pressure of Te2 over SnTe and lattice parameter measurements.Shelimova and Abrikosov [68], Sha and Brebrick [69], and Rogacheva et al. [70,71] determined the homogeneity range by hardness, partial pressure measurements, and lattice parameter as well as hardness measurements, respectively.The calculated phase diagram of the Sn-Te system is shown in Figure 3a along with experimental data [58][59][60][61][62][63][64][65] and the most recent optimization [56].The calculated enlarged homogeneity range of the SnTe compound is shown in Figure 3b along with experimental data [66][67][68][69][70][71].As seen in Figure 3b, compared to the previous calculation by Liu et al. [56], the calculated homogeneity range of the SnTe is in good agreement with the experimental data.  of SnTe are listed in Table 2. Blachnik et al. [81], Medzhidov and Rasulov [82], and Yamaguchi et al. [83] measured heat content (HT-H298.15K)using calorimetry.The calculated heat content (HT-H298.15K) is shown in Figure 5 along with experimental data [81][82][83].Recently, Pashinkin et al. [84] measured the heat capacity of SnTe using calorimetry in the temperature range 350-600 K.They reported entropy of SnTe at 298.15 K (S298.15K)as 102 ± 4.2 J/mol-K, which is in good agreement with the present optimized value 100.37 J/mol-K.Formation enthalpy of SnTe at 298.15 K (∆H f 298.15K ) was measured by Mcateer and Harry [76] using potentiometric analysis, by Robinson and Bever [77], Vecher et al. [78] and Misra [79] using calorimetry, and by Colin and Drowart [80] using Mass spectrometry.The ∆H f 298.15K of SnTe are listed in Table 2. Blachnik et al. [81], Medzhidov and Rasulov [82], and Yamaguchi et al. [83] measured heat content (H T -H 298.15K ) using calorimetry.The calculated heat content (H T -H 298.15K ) is shown in Figure 5 along with experimental data [81][82][83].Recently, Pashinkin et al. [84] measured the heat capacity of SnTe using calorimetry in the temperature range 350-600 K.They reported entropy of SnTe at 298.15 K (S 298.15K ) as 102 ± 4.2 J/mol-K, which is in good agreement with the present optimized value 100.37 J/mol-K.

Ga-Te Binary System
The optimized result of the Ga-Te binary system by the authors' recent study [29] was directly adopted in the present study.The optimized phase diagram and the integral enthalpy of mixing of liquid solution were shown in Figure 6a,b, respectively.The liquid immiscibility at the Ga-rich side and the asymmetric enthalpy curve were described well by using MQM with the modification of the coordination number ratio between Ga and Te as Z Ga GaTe /Z Te GaTe = 6/4.

Ga-Te Binary System
The optimized result of the Ga-Te binary system by the authors' recent study [29] was directly adopted in the present study.The optimized phase diagram and the integral enthalpy of mixing of liquid solution were shown in Figure 6a,b, respectively.The liquid immiscibility at the Ga-rich side and the asymmetric enthalpy curve were described well by using MQM with the modification of the coordination number ratio between Ga and Te as / = 6/4.

Ga-Te Binary System
The optimized result of the Ga-Te binary system by the authors' recent study [29] was directly adopted in the present study.The optimized phase diagram and the integral enthalpy of mixing of liquid solution were shown in Figure 6a,b, respectively.The liquid immiscibility at the Ga-rich side and the asymmetric enthalpy curve were described well by using MQM with the modification of the coordination number ratio between Ga and Te as / = 6/4.

Ga-Sn-Te Ternary System
The description of the ternary Ga-Sn-Te is not known over the whole composition and temperature range.There are some investigations on vertical sections Ga 2 Te 3 -SnTe [85,86] and GaTe-SnTe [86,87].The three ternary compounds, Ga 2 SnTe 3 , GaSnTe 2 [87,88], and Ga 6 SnTe 10 [85,86,89,90], have been reported in this system.Rustamov et al. [87] investigated the GaTe-SnTe system using the ultrasonic technique.They reported two ternary compounds, Ga 2 SnTe 3 and GaSnTe 2 , corresponding to the peak in velocity-composition isotherm obtained from the ultrasonic technique.However, they confirmed only one compound (GaSnTe 2 ) using thermographic and metallographic techniques, the detection of other peaks could be due to structural defects such as Sn vacancy.Guittard et al. [89] reported the ternary phase Ga 6 SnTe 10 using XRD and thermal analysis for 15 days at 600 • C.They reported crystal structure of the ternary phase Ga 6 SnTe 10 as rhombohedral (space group R32) and lattice parameters a = 10.207Å, α = 87.74• .Dedegkaev et al. [85] studied the Ga 2 Te 3 -SnTe system using XRD and reported the ternary phase Ga 6 SnTe 10 with a cubic primitive unit cell of lattice parameter 10.2 Å. Gather et al. [86] measured excess enthalpy using a heat flow calorimeter for the ternary liquid Ga-Sn-Te.Moreover, 2021, 11, 1363 10 of 15 they have determined the Ga 2 Te 3 -SnTe and GaTe-SnTe quasi binary phase diagram and confirmed the existence of ternary phase Ga 6 SnTe 10 using thermal analysis.The existence of ternary compounds Ga 2 SnTe 3 and GaSnTe 2 are still controversial [86].Therefore, ternary compound Ga 6 SnTe 10 was only considered in the present calculation.The Gibbs energy of the compound was estimated using the Neumann-Kopp rule with the slight modification of excess enthalpy term for best fitting the phase diagram in the quasi-binary Ga 2 Te 3 -SnTe.The liquidus line is measured by [87,89] for GaTe-SnTe and by [86,89] for Ga 2 Te 3 -SnTe.
As a first approximation of Gibbs free energy of the ternary liquid solution in the Ga-Sn-Te system, "Toop-like" asymmetric approximation [91] with Te as the asymmetric component was used.One adjustable ternary model parameter was added to avoid unexpected expansion of liquid miscibility gap toward the Sn-rich side.It was solved by increasing the stability of the ternary liquid solution with small negative parameter, g 001 GaTe(Sn) = −8900.The calculated vertical section GaTe-SnTe and Ga 2 Te 3 -SnTe are shown in Figure 7a,b along with experimental data [86][87][88][89], respectively.The predicted liquidus is in good agreement with the respective experimental data [86][87][88][89].The calculated mixing enthalpy of liquid Ga x Sn y -Te is shown in Figure 8 along with experimental data [86].The enthalpy of mixing decreases with increasing Ga/Sn ratio due to the higher contribution of Ga-Te, which is more negative than liquid Sn-Te.The same tendency was observed in the experimental data measured at 1203 K with Ga/Sn ratio.The calculated liquidus projection over the entire composition range of the ternary system Ga-Sn-Te is shown in Figure 9.The black solid lines in the liquidus projection represent the univariant equilibria.The colored lines are isothermal lines calculated at the interval of 50 K.The different types of invariant reactions in this system are given in Table 3.
detection of other peaks could be due to structural defects such as Sn vacancy.Guittard et al. [89] reported the ternary phase Ga6SnTe10 using XRD and thermal analysis for 15 days at 600 °C.They reported crystal structure of the ternary phase Ga6SnTe10 as rhombohedral (space group R32) and lattice parameters a = 10.207Å, α = 87.74°.Dedegkaev et al. [85] studied the Ga2Te3-SnTe system using XRD and reported the ternary phase Ga6SnTe10 with a cubic primitive unit cell of lattice parameter 10.2 Å. Gather et al. [86] measured excess enthalpy using a heat flow calorimeter for the ternary liquid Ga-Sn-Te.Moreover, they have determined the Ga2Te3-SnTe and GaTe-SnTe quasi binary phase diagram and confirmed the existence of ternary phase Ga6SnTe10 using thermal analysis.The existence of ternary compounds Ga2SnTe3 and GaSnTe2 are still controversial [86].Therefore, ternary compound Ga6SnTe10 was only considered in the present calculation.The Gibbs energy of the compound was estimated using the Neumann-Kopp rule with the slight modification of excess enthalpy term for best fitting the phase diagram in the quasi-binary Ga2Te3-SnTe.The liquidus line is measured by [87,89] for GaTe-SnTe and by [86,89] for Ga2Te3-SnTe.
As a first approximation of Gibbs free energy of the ternary liquid solution in the Ga-Sn-Te system, "Toop-like" asymmetric approximation [91] with Te as the asymmetric component was used.One adjustable ternary model parameter was added to avoid unexpected expansion of liquid miscibility gap toward the Sn-rich side.It was solved by increasing the stability of the ternary liquid solution with small negative parameter, g ( ) = −8900.The calculated vertical section GaTe-SnTe and Ga2Te3-SnTe are shown in Figure 7a,b along with experimental data [86][87][88][89], respectively.The predicted liquidus is in good agreement with the respective experimental data [86][87][88][89].The calculated mixing enthalpy of liquid GaxSny-Te is shown in Figure 8 along with experimental data [86].The enthalpy of mixing decreases with increasing Ga/Sn ratio due to the higher contribution of Ga-Te, which is more negative than liquid Sn-Te.The same tendency was observed in the experimental data measured at 1203 K with Ga/Sn ratio.The calculated liquidus projection over the entire composition range of the ternary system Ga-Sn-Te is shown in Figure 9.The black solid lines in the liquidus projection represent the univariant equilibria.The colored lines are isothermal lines calculated at the interval of 50 K.The different types of invariant reactions in this system are given in Table 3.

Application of the Database
Three vertical sections calculated using the present database in the Ga-Sn-Te system that have great potential as TE materials are shown in Figure 10a-c.The solid red circle in the figure indicates the prospective compositions of the potential eutectic alloys that may be investigated for TE application.

Application of the Database
Three vertical sections calculated using the present database in the Ga-Sn-Te system that have great potential as TE materials are shown in Figure 10a-c.The solid red circle in the figure indicates the prospective compositions of the potential eutectic alloys that may be investigated for TE application.

Application of the Database
Three vertical sections calculated using the present database in the Ga-Sn-Te system that have great potential as TE materials are shown in Figure 10a-c.The solid red circle in the figure indicates the prospective compositions of the potential eutectic alloys that may be investigated for TE application.

Conclusions
In this study, binary Sn-Te and Ga-Sn and ternary Ga-Sn-Te systems are thermodynamically assessed based on available experimental information using the CALPHAD method.The liquid solution in binary Sn-Te was reassessed by using the Modified Quasichemical Model (MQM) to account for the short-range ordering (SRO) exhibited by the liquid solution.To keep the consistency with other sub-binaries, the liquid Ga-Sn solution was optimized using MQM.The Gibbs free energy of thermoelectric compound, SnTe was optimized using Compound Energy Formalism to explain the hightemperature homogeneity range of the phase.BCT in the Ga-Sn binary system was optimized using a regular solution model.A set of self-consistent thermodynamic descriptions was obtained for all the phases in the binary system Ga-Sn, Sn-Te, and that of Ga-Te binary system was taken from our previous optimization.The available thermodynamic and phase diagram data were reproduced within the experimental error limits.Combining the optimized sub-binary systems, thermodynamic properties of the ternary Ga-Sn-Te system were predicted well with only one adjustable ternary model parameter.Several eutectic compositions were predicted, and these might be used for thermoelectric applications.Furthermore, the optimized database will be extended to a larger telluride-based system by adding elements such as In and Bi to explore the multiphase eutectic systems for thermoelectric (TE) applications.

Figure 2 .
Figure 2. (a) Calculated activity of Ga (l) and Sn (l) in liquid at 1073 K along with experimental data from [42,45,46] and (b) mixing enthalpy of liquid at 723 K along with experimental data from [42,47-51].
57].Hsieh et al.[52], Kattner et al.[53],Sharma and Chang [57],Kuznetsov [54], and Liu et al.[56] used an associate solution model with SnTe as an associate component, whereas Merzhanov et al.[55] used a formal mathematical model.All proposed thermodynamic descriptions can satisfactorily reproduce reported experimental thermodynamic and phase diagram data in this system.However, Hsieh et al.[52] and Kattner et al.[53] used old thermodynamic descriptions of the pure elements, andKuznetsov [54]  and Sharma and Chang [57] obtained model parameters that were not consistent with other CALPHAD type modeling.They have extensively reviewed the available thermodynamic data and those data were employed in the present study to determine Metals 2021, 11, 1363 7 of

Figure 3 .Figure 3 .
Figure 3. (a) Calculated phase diagram of Sn-Te system along with experimental data and (b) enlarged homogeneity range of SnTe compound along with experimental data from [58-71].Le Bouteiller et al. [64], Rakotomavo et al. [63], and Nakamura et al. [72] have performed EMF measurements on galvanic cells to determine thermodynamic propertiesFigure 3. (a) Calculated phase diagram of Sn-Te system along with experimental data and (b) enlarged homogeneity range of SnTe compound along with experimental data from [58-71].Le Bouteiller et al. [64], Rakotomavo et al. [63], and Nakamura et al. [72] have performed EMF measurements on galvanic cells to determine thermodynamic properties such as activity, partial enthalpy in liquid phase at 673-913 K, 673-1123 K, and at 1100 K, respectively.Predel et al. [73] determined the activity of Te (l) by partial pressure measurements.Partial enthalpy of Sn (l) and Te (l) at high dilution in liquid Sn-Te was measured by Amzil et al. [74] using reaction calorimetry.The mixing enthalpy and entropy of the liquid solution were measured using calorimetry by Blachnik and Gather [75] and Rakotomavo et al. [63], and EMF by Nakamura et al. [72].The calculated activity of Sn (l) and Te (l), mixing enthalpy, mixing entropy, and partial enthalpy of Sn (l) and Te (l) in the liquid solution at 1100 K are shown in Figure 4a-d along with respective experimental data [63,64,72-75].The calculated thermodynamic properties are in good agreement with the experimental data.

Figure 4 .
Figure 4. Calculated thermodynamic properties at 1100 K (a) activity of Sn (l) and Te (l), (b) mixing enthalpy of liquid, (c) mixing entropy of liquid, and (d) partial enthalpy of Sn (l) and Te (l) in liquid solution along with respective experimental data from [63,64,72-75].

Figure 6 .
Figure 6.Calculated (a) Ga-Te phase diagram and (b) integral enthalpy of mixing at 1123 K using thermodynamic description of Kumar et al. [29].

Figure 5 .
Figure 5. Calculated heat content (H T -H 298.15K ) of SnTe as function of temperature along with experimental data from [81-83].

Figure 6 .
Figure 6.Calculated (a) Ga-Te phase diagram and (b) integral enthalpy of mixing at 1123 K using thermodynamic description of Kumar et al. [29].

Figure 6 .
Figure 6.Calculated (a) Ga-Te phase diagram and (b) integral enthalpy of mixing at 1123 K using thermodynamic description of Kumar et al. [29].

Figure 8 .
Figure 8. Calculated enthalpy of mixing in liquid GaxSny-Te at 1203 K along with experimental data from [86].

Figure 9 .
Figure 9. Calculated liquidus projection in the Ga-Sn-Te system.

Figure 8 . 16 Figure 8 .
Figure 8. Calculated enthalpy of mixing in liquid Ga x Sn y -Te at 1203 K along with experimental data from [86].

Figure 9 .
Figure 9. Calculated liquidus projection in the Ga-Sn-Te system.

Figure 9 .
Figure 9. Calculated liquidus projection in the Ga-Sn-Te system.

Table 1 .
Optimized model parameter of the Ga-Sn-Te system (All values are in J/mol).

Table 2 .
Formation enthalpy of SnTe at 298.15 K.

Table 2 .
Formation enthalpy of SnTe at 298.15 K.

Table 3 .
Invariant reactions in Ga-Sn-Te system.

Table 3 .
Invariant reactions in Ga-Sn-Te system.