Thermodynamic Properties of Li-Sb Liquid Solution by QAM

The optimization of the thermodynamic properties of a liquid Li-Sb system is performed based on all available literature data. The description of the excess Gibbs energy as a function of temperature and composition is presented using a qualitative associate model (QAM), described in our earlier work on the Li-Pb system. The modeled molar mixing enthalpy changes, excess Gibbs energy and the molar excess entropy values are characterized by strong negative deviations from the ideal solution in the entire concentration range. It was found that the best correlation between experimental data and those calculated using QAM is observed if the Li3Sb and LiSb associates are assumed to be in the liquid Li-Sb solution. Calculations of the microscopic function, i.e., the concentration-concentration in the long-wavelength limit, (Scc(0)) excess stability function (ES) and short-range-order parameter (SRO(α)) show the extreme (maximum or minimum) values almost exactly for the concentrations of Li3Sb associate (and intermetallic phase in solid alloys) and LiSb. The quotient of real Scc(0) to this for the ideal solution (Scc(0)id) is characterized by a similar property to the mentioned structure factors.


I. INTRODUCTION
IT is currently necessary to find an electrode material with a high reversible capacity and low overall cost of production. The Li and Sb-containing alloys are of main interest because of their good energy storage properties and ability to be used in the future generation of Li-Ion or high-temperature liquid batteries, with high operational parameters. [1][2][3][4][5] Preparation of such electrode materials requires information about their thermodynamic properties and phase equilibria. Thus, reliable knowledge about the thermodynamic properties of the Li-Sb system is very important.
The thermodynamic investigation of the liquid phase existing in the Li-Sb system was published in 2015 by Kane et al. [6] In their paper, [6] the activities of Li in Li-Sb liquid and liquid-solid solutions were obtained by electromotive force measurements (EMF) in the temperature range from 698 K to 1048 K and concentration range of x Li between 0.008 and 0.75. Based on the obtained EMF results, Kane et al. [6] calculated the partial enthalpies and entropies of Li in liquid Li-Sb alloys. Then, Terlicka et al. [7] and later Beutl et al. [8] and Li et al. [9] presented the mixing enthalpy of liquid Li-Sb alloys. Also, Trybula et al. [10] investigated the thermodynamic properties of Li-Sb liquid solutions using EMF measurements. One should also notice that the differential thermal analysis measurements conducted by Beutl et al. [8] showed the different shape of the liquidus temperature compared to that proposed by Sangster and Pelton. [11] Recently, the phase diagram of the Li-Sb system was calculated in 2017 by Zhang et al. [12] and also by Li et al., [9] both using the association model.
Liquid solutions of Li-Sb are characterized by a significant negative deviation from the ideal one. In this case, the change in the enthalpy of mixing has a minimum near the concentration corresponding to one of the intermetallic phases (compounds) in the solid state. Its shape on both sides of this concentration is similar to the linear function. A very similar situation occurs inter alia for systems such as Li-Bi, [13] Li-Pb [14][15][16] and Li-Sn. [17] What is more, the partial mixing enthalpies are characterized by the S-shape (the inflection point occurs at the concentration close to that for the phase) with a very rapid increase (decrease) of the value in the vicinity of the intermetallic phase (compound). As is shown by, e.g., Saboungi et al. [18] and Gąsior et al., [19] this shape of the partial and molar function is very difficult to describe with one of the different often used types of polynomials (models), e.g., the power series or Redlich-Kister. In such cases, as shown in Saboungi et al. [18] and Gąsior et al., [19] it is necessary to use two polynomials (or the graphical interpretation of experimental data): one describing a given thermodynamic function in the range of solutions with Li concentration greater than the intermetallic phase and the other for the remaining concentration range. However, in such a case it is necessary to keep the condition that derivatives calculated using both polynomials are equal at the point of inflection of a partial function. Another difficulty is that, by increasing the polynomial degree to better approximate the experimental data, the approximation function begins to wave, i.e., it has many inflection points or local minima and maxima. A completely different behavior is seen in the model (formalism) proposed in References 14 and 20 called QAM. It can describe very rapid changes of the thermodynamic isotherm, in the vicinity of the concentration characteristic of one of the intermetallic phases observed in the solid state, by using only one or two parameters calculated based on experimental data. Comparing the results of the elaboration of thermodynamics of the liquid phase by the association model, the authors of References 9 and 12 used 11 and 12 parameters for the description of thermodynamic properties of the Li 3 Sb associate and the excess Gibbs energy of three binary solutions, Li-Li 3 Sb, Li-Sb and Li 3 Sb-Sb, by the Redlich-Kister model. [21] Moreover, the reaction equilibrium constant and the additional equation have to be applied to determine the mole number of Li 3 Sb, Sb and Li. As can be seen, the use of the association model for the description of the composition dependence of excess Gibbs energy of binary solutions needs a large number of coefficients (11 (12) + equilibrium constant), and additionally a procedure for the determination of a number of moles based on the equation of the equilibrium constant, unless the equation for calculation of the number of moles is given (additional parameters).
Compared to the cited earlier publications, [9,12] where the association model was used, we would like to show that the interpretation of thermodynamic properties for systems with a large negative mixing enthalpy change and excess Gibbs energy is possible with a much lower number of coefficients and with acceptable standard deviation using the QAM model presented in References 14 and 20.
Therefore, the main aim of this work was to optimize the thermodynamic properties of a liquid Li-Sb system based on all literature reports on the mixing enthalpy change and the partial excess Gibbs energy of lithium using a qualitative associate model (QAM). The correlation between experimental and calculated data for the assumed type of associates was done. Moreover, the calculation of the structure factors based on the QAM was also made and compared with the literature data.

II. OPTIMIZATION PROCEDURE
The qualitative associate model can be applied to determine the molar and partial mixing excess Gibbs energy for binary systems, which have strong negative deviations from the ideal solutions. In such a case, the curves describing thermodynamic functions, on both sides of the concentration characteristic for the assumed associate, have very high curvature (the radius of curvature is so large that the graph resembles almost straight lines). If the considered binary system has significant negative deviations in thermodynamic properties, the model used to develop such properties must take into account the above-mentioned linearity. The association model presented by Sommer is one such model. [22] However, as was shown in References 9 and 12, in this case the number of optimized thermodynamic parameters must be large and additional calculations of a number of moles of solution components have to be conducted. Thus, only the use of specialized programs will allow the determination of thermodynamic data. It is completely different for other methods where knowledge of optimized thermodynamic parameters allows calculation of the thermodynamic properties of solutions (e.g., Redlich-Kister, [21] power series or orthonormal polynomials models). Accepting this concept, a new approximation method was worked out [20] and applied in our earlier paper. [14] Because it applies only the type of associate without quantitative relations, it was named the qualitative associate model (QAM). This model can be also used to elaborate the excess Gibbs energy based on the data from the association model.
The starting point for working out this model was the observation that the description of the data discussed above is possible if the equation for the excess Gibbs energy is a modification of the following one: The plot of the concentration function (f(X)) in Eq.
[1] is drawn in Figure 1 for different values of N and M (1 À Y, Y) parameters. Notably, the shape of curves is very similar to the plots of the molar mixing thermodynamic functions of the systems (e.g., the enthalpy and Gibbs free energies of mixing) with large negative deviation from the ideal behavior.
If M equals 2, 3 and 4 and N equals 1, the maximum of the f(X) has the same value as for N = 1 and M = 2, 3 and 4, however, for X higher than 0.5 (2/3, 3/4, 4/5). An equation suitable for the development of experimental data cannot be obtained only by the change of the exponent in Eq. [1] from 0.5 to another one. Additionally, the concentration function f(X) had to be modified. The simplest form of the equation for the description of thermodynamic properties of binary solutions with strong negative deviations and with one associate deduced by the author is as follows: ) on the m parameter and X was also described in our previous paper [14] (the m parameter is simply an adjustable parameter that is dependent on the experimental data by the meaning of the best fit). The partial excess Gibbs energy of components can be calculated by a combination of the Gibbs-Du¨hem relation and Eq. [2]. The integral and partial functions may be written as follows: :

½5
The meaning of all symbols in Eqs. [3] through [5] has been explained above in this chapter. If, instead of N and M, the 1 À Y and Y coefficient of the A 1ÀY B Y associate are used, then Y + (1 À Y) = 1, and Eqs. [3] through [5] are transformed as follows: If two associates are necessary to describe the thermodynamic properties of the solution, Eqs. [2] through [8] should be a sum of two expressions representing the excess Gibbs energy of each associate.
Notably, if m = 0.5, the function of composition f(X) in Eqs. [1] and [6] has the same shape as the that in Figure 1, depending on the value of N and M. When has the value 0.5 and is equal to the maximal mole number of A 1 B 1 associate, which can be formed in the binary A-B solution (the equilibrium reaction constant of associate formation K as equals infinity). For the other values of N and M or 1 À Y and Y, the maximal values of f(X) can be calculated from the equation: and the maximal mole number of A 1ÀY B Y associates equals Y. As a result, it changes linearly from Y to 0 or from Y to 1 when the composition of the solution changes from 0 to Y or from Y to 1 (X = Y for the A 1ÀY B Y associate type).
In the case of one associate, the calculations of the a, b (c and d if necessary) and m parameter are very easy using the minimization procedures described in the next chapter. The problem is more complicated when 2 (or 3) associates occur, because a minimized function may possess many minima, and the obtained a i , b i , c i , d i and m i (i = 1, 2, …) parameters could be improper; this means, for example, that the calculated excess molar entropy may have very high positive or very low negative values.
To avoid this problem, additional observations were made. It was noticed that the smaller the negative value of the excess Gibbs energy (or the enthalpy of mixing), the more its shape is similar to that in Figure 1. This means that the lower the negative value of the thermodynamic function, the closer to 0.5 the value of m should be. This fact is easy to observe when analyzing, for example, the change of the enthalpy of mixing for Li-Sb and Li-Pb [14] systems (for Li-Sb, m @ 0.5, D mix H @ À 71 kJ mol À1 and for Li-Pb, m @ 0.59, D mix H @ À 30 kJ mol À1 ).
The second observation concerns the difference between m and 0.5. It was noticed that Dm = m 1 À 0.5 is close to the following: where: where a, b, c and d are the coefficients of temperature function f(T) in Eqs. [1] through [8], T av is the average temperature calculated using all experimental data, and e is the base of the natural logarithm (e = 2,71828…). T av is calculated just after loading all experimental data into the computer and before the optimization process (minimization of the residual sum of squares). Parameters m 2 or m 3 , characteristic for the second or third associate (if necessary), are calculated as follows: where G 2ð3Þ is the value of the temperature function . Using Eqs. [12] and [13] for the m 1(2,3) calculations, and in the case of linear dependence of f(T) function, the a i , b i parameters have only one value for the given set of experimental data for which the residual sum of squares has the minimum.

III. RESULTS AND DISCUSSION
Using the available thermodynamic data such as the mixing enthalpy change from the solution calorimetric method [7,9] and the partial excess Gibbs energy of Li in Li-Sb liquid alloys measured by the EMF study, [6,10] the adjustable parameters a, b of the f(T) function in Eqs. [1] through [5] were calculated by the linear regression method, and the residual sum of squares (RSS), used as an optimal criterion, is calculated according to the following equation: where G cal i or H cal j and G exp i or ðH exp j Þ denote the excess Gibbs energy of Li ðG ex Li Þ or mixing enthalpy change (D mix H), calculated and experimental, respectively, and k and l are the number of experimental data of the cited thermodynamic functions.
The a i and b i parameters for which the residual sum of squares (RSS) has the minimum was calculated using the SIMPLEX method (computer program).
The calculations were conducted in three variants. In the first one (I), it was assumed that only the Li 3 Sb associate occurs in the liquid Li-Sb solutions. In the second (II) and third (III) variants, the Li 3 Sb and Li 2 Sb and the Li 3 Sb and LiSb associates were assumed to exist in the Li-Sb solution.
The results of calculations are presented in Table I, where the computed standard deviations for the enthalpy of mixing (SD H ), partial excess Gibbs energy of Li (SD GE ), global one (SD G ) and coefficients a, b and m are shown.
For the first variant, with only one associate (Li 3 Sb), the standard deviation was the highest. In the case of the second and third variant (Li 3 Sb and Li 2 Sb, Li 3 Sb and LiSb associates, respectively), all calculated standard deviations are comparable. However, because the Li 2 Sb intermetallic phase decomposes at about 754 K, it was assumed that the Li 3 Sb and LiSb associates exist in the liquid phase. Thus, the results of the calculation in this study are presented for variant (III) (bold).
The integral molar mixing enthalpies and entropies, as well as excess Gibbs energy of mixing for the liquid phase calculated using the qualitative associate model for 973 K, are collected in Table II  All presented molar thermodynamic functions have an almost triangle shape (which means that the radius of curvature is so large that the graph resembles almost straight lines and the equilibrium reaction constant of associate formation K as equals infinity). This means that the partial ones must be of S-shape with the inflection point for x Li = 0.75 (Li 3 Sb), which is confirmed in Figures 2 through 4. Presented in Figure 2, the calculated DG E Li data and experimental values [6,10] are generally in good agreement, and the observed deviation are < 6 pct.
The DG E Li and DG E Sb at 973 are negative in the entire concentration range, and their minimal values are À 79947 and À 222,595 J mol À1 , respectively, in contrast to the partial mixing enthalpy change and partial excess entropy of Li and Sb, which are characterized by both negative and positive values. Such behavior of the mentioned partial functions is the effect of the slightly convex nature of the change of enthalpy of mixing D mix H for the alloys of compositions < x Li = 0.55. The limiting values of partial enthalpy and the partial excess entropy of Li at 973 K are -86,599 J mol À1 and À 6.84 J (mol K) À1 , respectively. The same functions for antimony reach À 318657 J mol À1 and À 98.73 J (mol K) À1 , respectively. Both partial enthalpies and excess entropies have minima and maxima. The Li-Sb alloys have large negative deviations of thermodynamic properties from the ideal solution, which has been already described in References 7 through 10. Thus, it is probable that in the liquid phase the associates are forming as the result of the short-range ordering. This problem was discussed, inter alia, by Ruppersberg et al. [23,24] and Bhatia and Thornton [25] based on the results of the intensity I(k) of the scattered radiation for the liquid Li-Pb solutions. Its correlation with the fluctuation theory is shown below.
The intensity of the scattered radiation of the x-ray or neutrons is connected with the static structure factor by the relation: where I(k) is the intensity of the scattered radiation, S(k) is the static structure factor, N is the number of scattering points (for N = 1 mole it is Avogadro's number), f is the atomic scattering factor, a is half of the scattering angle, and k is the incident   [7] at T = 1024 K and [9] at T = 922 K to 979 K. wavelength.In turn, the static structure factor S(k) is the function of the pair distribution function gðrÞ (and vice versa), which measures the probability of finding an atom at a distance r from another atom at r = 0. For two kinds of atoms (binary solution), the probability of finding three pair of atoms, A-A, B-B and A-B, at distance r must be considered. Therefore, in this case there are three partial structure factors: S AA (k), S BB (k) and S AB (k) [also denoted as a AA (a 11 ), a BB (a 22 ) and a AB (a 12 )] and three pair distribution functions: g AA ðrÞ; g AB ðrÞ and g BB ðrÞ. As Faber and Ziman [26] or Enderby et al. [27] and Bhatia and Thornton [25] showed, the intensity of the scattered radiation I(k) can be expressed in the form of a linear equation of the structure factors depending on the radial distribution functions [26] or the structure factors depending on the fluctuation of concentration S CC ðkÞ, the number of atoms S NN ðkÞ [25] and the partial factor S NC ðkÞ being the correlation function of S CC ðkÞ and S NN ðkÞ. So, I(k) can be expressed by equations:

the mole fraction of A-B alloy components and f A and f B are atomic scattering factors of A and B atoms.
As Bhatia and Thornton [25] show, for the long wavelength limit, this means that for k = 0, S cc (0), S NN (0) and S NC (0) can be calculated directly with the thermodynamic and physical properties of the alloy according to the equations: where j T is the isothermal compressibility, V 0 is the molar volume of alloys, V A ; V B are the partial molar volumes of the A and B components, respectively, R is the universal gas constant, T is temperature, N is Avogadro's number, ðDXÞ 2 is the mean-square fluctuation in concentration, ðDNÞ 2 is the mean-square fluctuation in a number of particles, DNDX is the correlation between the two fluctuations, and D mix G is the molar Gibbs free energy change.
As the above equations show, S CC (0) can only be calculated based on the thermodynamics of a solution.  Points at the temperature range T = 873 K to 973 K: d, [6] and . [10]  The discussion of S CC (0) for the different solutions was presented, inter alia, by Ruppersberg et al., [23,24] Bhatia and Thornton [25] and Shing et al. [28] The general conclusion is that for the alloys the negative deviations from the ideal, once S CC (0) shows lower values than for the ideal solutions, and the minimum value of S CC (0) is observed for the concentration, one of the intermetallic compounds occurred in a given system. However, for many systems characterized by the negative deviation from ideality, S CC (0) calculated based on the available equations describing thermodynamics (excess Gibbs free energy) of liquid solutions shows an inflection, not a minimum for the composition of the intermetallic phase observed in the solid solutions. The reason for this is the use of an improper model(s) to elaborate the experimental thermodynamic data. Continuing the discussion on the structure of alloys with high negative deviations of D mix G E from the ideal solution, note that S CC (0) can be calculated not only using the molar Gibbs energy change but also the activity and activity coefficient, applying the following relations: where D mix G is the mixing Gibbs energy (D mix G ¼ D mix G id þ D mix G E ); a A(B) and X A(B) are the activities and the mole fractions of components A and B (Li and Sb), respectively; T is the temperature in Kelvin, R is the universal gas constant, and S CC (0) i is the limit value of S CC (k) in the case of ideal solutions. The S CC (0) for the real and ideal liquid solution for the Li-Sb system is presented in Figure 5, which also shows the difference DS CC (0) [22] between the ideal and real concentration factors. The latter was suggested to be more proper for the identification of the type of associate in the case of the lack of evident minimum at the S CC (0) curve. As could be expected, the S CC (0) dependence on the concentration of Sb at 1500 K, visible in Figure 5, presents a clear minimum located almost exactly for x Li = 0.75 (Li 3 Sb), which in the literature is thought to be the associate's creation in the liquid phase. The presence of this associate in the liquid Li-Sb solutions was expected because of the Li 3 Sb intermetallic phase in solid solutions and deep minima of thermodynamic functions analyzed above in the vicinity of composition x Li = 0.75. Consequently, the assumption of this type of associate was assumed in the elaboration of thermodynamic properties by Eqs. [3] and [5], and, in this case, the minimum at the S CC (0) curve for x Li = 0.75 must have appeared. Figure 5 shows that the DS CC (0) maximum is for the alloy of x Li = 0.73 and not for Li 3 Sb. The second weak minimum is for x Li = 0.5, which corresponds to LiSb associate. However, the minimum value for Li 3 Sb is much lower than that for LiSb.
The influence of temperature on the S CC (0) and DS CC (0) is shown in Figure 6. As is clearly seen, the increase of temperature changes the shape of both functions and the dominant influence of Li 3 Sb at 1000 K is lowered at 1500 K and 2000 K; a greater influence of LiSb is also observed. As a result, the lower value of S CC (0) of the LiSb associate and higher value of  S CC (0) of the Li 3 Sb associate is observed at 2000 K compared to that at 1000 K and 1500 K. This means that the amount of LiSb associates increases and of Li 3 Sb associates decreases with increasing temperature. A similar conclusion can be made by analyzing the DS CC (0).
The second function, which suggests the existence of associates in the liquid phase, is the excess stability function (ES) proposed by Darken. [29] This function is directly related to the partial structure factor S CC (0) and is defined as the second derivative of excess Gibbs energy (D mix G E ) of the concentration of the component (x i ): Substituting the D mix G E as a sum of partial excess Gibbs energy of components (D G E i ) multiplied by its mole fractions (x i ) and performing appropriate mathematical transformations, the following equation is reached: where c i is the activity coefficient of component i.
The ES function usually presents the maximum near or at the concentration of the high-temperature intermetallic phase. When the above-mentioned maximum exists in the studied system, the existence of associates in the liquid phase is highly probable. The calculated ES function for the Li-Sb system at three different temperatures (1000 K, 1500 K and 2000 K) is presented in Figure 7. A clear maximum exists for x Li = 0.75, which suggests the occurrence of associates of the chemical formula Li 3 Sb. A slight maximum is also noticeable for x Li = 0.50 (LiSb). Moreover, the influence of temperature on the structure of Li-Sb liquid alloys is obvious, and the same conclusions can be drawn as in the case of S CC (0). The next parameters indicating the occurrence of the associate(s) in the liquid phase is the Warren-Cowley short-range-order parameter (SRO(a)). [28,[30][31][32][33] Its values vary from 1 to À 1. If it is negative, this means that in the liquid solutions the interaction between different atoms is much stronger than between the same atoms and the associates can be formed. If the value of SRO(a) is positive, the interaction is stronger between the same type of atoms, so the formation of associates with different atoms is impossible, and such solutions will tend to segregate.
The calculation of the SRO(a) was done using the following formula: where: where x 1 , x 2 is the molar fraction of the binary system components (Li, Sb), and z is the coordination number of liquid solutions. Based on the calculation presented in, [33] the z equals 10 for the liquids. The values of the chemical short-range-order parameter SRO(a) for liquid Li-Sb solutions are presented in Figures 8(a) and (b). They are negative in the entire concentration range, with the strong, clear minimum at x Li = 0.75. This decrease near the composition of the high-temperature intermetallic phase confirms that, in the liquid phase of Li-Sb system, the Li 3 Sb associate exists. For comparison, Figure 8(a) also shows the values obtained by Trybula et al. [10] from the elaboration of experimental data by the Sommer's association model. [22] The great differences observed are the consequence of the interpretation of the experimental data. The development of SRO(a) with the temperature increase is presented in Figure 8(b). Considering that stronger interaction of atoms should be manifested by the lower values of SRO(a), we can state that the higher the temperature, the worse the interaction between components of Li 3 Sb is and the more LiSb associates exist in the liquid solution. It is because with the increase of temperature the solutions lead to the ideal ones.
Finishing the discussion on the structure of liquid Li-Sb solutions, the authors would like to draw attention to the property of the S parameter in the equation of SRO(a) (Eq. [28]), used also by Reference 34 to model the diffusion coefficients (the inverse of the Q parameter is used in the calculation of viscosity and diffusion coefficients by Singh and Sommer [34] ), which is the quotient of the real S CC (0) to S CC (0) id of ideal solutions. Figure 9 shows the same three temperatures (1000 K, 1500 K, 2000 K) as presented above for the structure factors.
Q is also exactly characterized by the deep minimum for the proper composition of Li 3 Sb. Its values change from 0 to 1 for the systems with a negative deviation from the ideal, and those with positive ones have values > 1. Hence, it also can be a kind of indicator of the existence of associates in liquid solutions. One can also notice that Q is independent of the coordination number z as a, and it is calculated directly only from the thermodynamic properties of liquid solutions similar to S CC (0) and ES. The conclusions that can be reached from the analysis of Q at different temperatures are the same as in the case of S CC (0), ES and SRO(a).
The changes in the enthalpy of the mixing of liquid Li-Sb alloys calculated using the QAM model at 1000 K are illustrated in Figure 10(a) and compared with the experimental [7,9] and modeled data presented by Li et al. [9] via an associate-based model at 1000 K and Trybula et al. [10] by associate-based and sub-regular models, both computed at 1000 K. There is a very good correlation between the values of the enthalpy of mixing obtained experimentally and those obtained by the proposed QAM model for Li-Sb alloys. This is contrary to the models used by Li et al. [9] and Trybula et al., [10] where the experimental values differ more significantly from those obtained through the experimental methods. It is especially noticeable for alloys rich in Li (0.9<x Li < 0.8). Moreover, the predicted D mix H by References 9 and 10 is less negative compared to that presented in this work model and in Reference 12. The shape of D mix H of all visible thermodynamic models is generally asymmetrical. However, the minimum of this function differs from about x Li = 0.68 for both models presented by Trybula et al. [10] (blue and red lines) to x Li = 0.70 for Li et al. [9] (green line) up to x Li = 0.75 for the QAM model and [12] (black line).
The comparison of the partial enthalpy change of Li and Sb is shown in Figure 10(b), where the experimental results of References 6, 7 and 9, are drawn along with those calculated using the QAM model presented in this work and the association model by Reference 12 used in the calculations of the Li-Sb phase diagram. The experimental values [9] of partial enthalpy of Li and Sb agree very well with those obtained by the QAM model. Comparing the modeled data with the experimental values, as was done by Zhang et al. [12] shows better compatibility of DH Li values and somewhat worse agreement of DH Sb . The difference between the measured and calculated values of partial enthalpy of Sb is because literature data were probably measured in the solid-liquid range. [7,9] In Figures 11 and 12 the partial Gibbs energy change of Li and Sb is compared with those from Reference 10 at three temperatures. In Figure 13, the same functions are shown only at 1000 K and compared with the experimental data. [6,7] The observed differences between experimental and modeled values at 1000 K ( Figure 13) are the lowest for the QAM model presented in this work (maximum 2.5 kJ mol À1 ), and the highest is for the sub-regular model (about 50 kJ mol À1 ). Maximum deviations of the values calculated by the association  model [10] and from the experimental results are about 10 kJ mol À1 for the solutions from the concentration range between x Sb @ 0.6 to 0.8. Data obtained by Reference 9 for x Sb @ 0.4 to 1 and 0 to 0.15 are close to those gained from QAM model and experiment. [6,7] Between x Sb @ 0.15 to 0.4 the courses of lines are completely different. Analyzing the shape of functions of the partial molar Gibbs energy of Li and Sb calculated applying the QAM and association model [10] shows that, first, their inflection points are located in different concentrations. In the case of the QAM model, the inflection is almost exact for x Sb = 0.25, independent of temperature, and for Reference 10 the inflection is for x Sb @ 0.85. The shapes of the partial molar Gibbs energy of Li and Sb are also different.

IV. CONCLUSIONS
Based on the newest available thermodynamic property data on the Li-Sb system with strong negative deviations from ideal solutions, a description of the integral molar and partial quantities such as mixing enthalpies, molar excess Gibbs energy and molar excess entropy was proposed by means of a qualitative associate model.
The calculated integral molar mixing enthalpies agree very well with the available experimental data. Moreover, the elaborated partial excess Gibbs energies are in  [9] ; sub-regular model of Trybula et al. (red line) [10] ; association model of Trybula et al. (blue line). [10] (b) Partial enthalpy of Li and Sb calculated at 973 K by a new QAM model (DH Li black solid line, DH Sb blue solid line) together with the literature data: association model presented by Zhang et al. [12] (DH Li : black dashed line, DH Sb : blue dashed line) and experimental data: d [6] ; [7] at T = 922 K, m [7] at T = 820 K to 1024 K, j [9] at T = 929 K to 1080 K (Color figure online).  The two microscopic functions S CC (0), and SRO(a), together with the excess stability function (ES), obtained in this work, present a clear maximum and/or minimum located almost exactly at the composition proper for the intermetallic phase of higher melting temperature from the Li-Sb system. These facts suggest the prevalence of Li 3 Sb associates in the liquid phase, at least near the melting points of alloys. Moreover, the above-mentioned functions also present extrema (minimum or maximum) for x Li = 0.50 at higher temperatures, which suggests the existence of a second LiSb associate in a liquid Li-Sb solution. Additionally, the analysis of microscopic functions for different temperatures shows that the higher the temperature is, the higher the values of S CC (0), SRO(a) and Q and the lower the ES for the concentration near Li 3 Sb. This suggests that the stability of Li 3 Sb associates decreases with increasing temperature. The reverse trend is observed for the LiSb associate.

CONFLICT OF INTEREST
There are no conflicts to declare.

OPEN ACCESS
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.