Thermodynamic, Structural, Surface and Transport Properties of In-Tl Liquid Alloy at Different Temperatures ()
1. Introduction
Various properties, such as hardness, corrosion resistance, wear resistance, mechanical strength, and fatigue strength, melting and boiling temperatures etc. of an alloy are different from its individual components. The properties of alloys in solid state are usually determined by thermodynamic, surface, structural, transport and electrical properties of alloys in liquid state. But there is complexity in determining the properties of liquid alloys due to lack of long range atomic order. Therefore several models [1] - [9] have long been proposed to understand the properties of alloys in liquid state.
The binary In-Tl alloy has widely been studied in the past because of its wide range of practical applications [10] [11] [12] . It is a shape memory alloy [13] having low melting temperature (fusible alloy). The indium based fusible alloy has the ability to wet and adhere to ceramic and glass surfaces. This characteristic can be useful in sealing and bonding of sensitive electronic devices, glass, ceramic and temperature-sensitive materials.
In this work we have used regular solution model [14] to explain thermodynamic and structural properties of In-Tl at 723 K. In regular solution model, the interaction energy is considered as input parameter and is determined by fitting observed free energy of mixing at different concentrations. Thermodynamic behavior of the alloy is studied by computing free energy of mixing (GM), activity (a), entropy of mixing (SM) and heat of mixing (HM). Structural properties is interpreted by the concentration fluctuation in the long wavelength limit (Scc(0)) and chemical short range order parameter (α1). The surface properties are studied by using Butler’s model [15] and viscosity is studied with the help of Moelwyn-Hughes equation [16] . In order to compute all these properties at higher temperatures optimization procedure is used. Basic theoretical formulation of regular solution model has been presented in Section (2), result and discussion in Section (3) and conclusion are dealt in Section (4).
2. Theoretical Formalism
2.1. Thermodynamic Properties
Regular solution model [14] deals with the A-B alloy in liquid state in which the constituent atoms A and B are sufficiently similar in size and shape so that they are interchangeable on the lattice or quasi-lattice, and the configuration is no longer independent of the mutual disposition of the two or more kinds of molecules. We suppose a liquid binary mixture A-B having cA (
c) mole of A and cB {
} mole of B respectively, where cA and cB are the mole fractions of A (
In) and B (
Tl) in the binary liquid alloy A-B.
The expression for free energy of mixing (GM) can be derived on the basis of regular soultion model. It is given as [14]
(1)
where R is the universal gas constant, T the temperature, and ω the interaction energy parameter to be determined.
The term
represents the heat of mixing in the frame work of regular solution model, i.e.
(2)
The activity aA of the element A in the binary alloys A-B is given by the standard relation
(3)
Using Equations (1) and (3) we get the activity of element A as
(4)
Similarly, the activity of the element B is given by
(5)
The temperature derivative of GM provides an expression for integral entropy of mixing (SM) which is given by
(6)
From Equations (1) and (6), we get
(7)
The relation between free energy of mixing (GM), entropy of mixing (SM), and heat of mixing (HM) is given by a standard thermodynamic relation,
(8)
Using Equations (1), (7) and (8), we get
(9)
2.2. Structural Properties
The concentration structure factor or concentration fluctuation in the long wavelength limit (Scc(0)) which has been used to study the nature of atomic order in binary liquid alloy [17] can also be deduce using regular solution model and is given by using relation,
(10)
From Equations (1) and (10), we get
(11)
The observed value of concentration fluctuation in the long wavelength limit Scc(0) is obtained from observed data of the activities of the constituent species of the binary liquid alloys from the relation
(12)
where, aA and aB are the activities of the component of A and B respectively.
The degree of local order in the liquid mixture is quantified by the Warren-Cowley [18] [19] short range order parameter (α1). It provides insight into the local arrangement of the atoms in the molten alloys. The theoretical values of (α1) can be evaluated as
(13)
where,
,
and Z is the coordination number, which is taken 10 (= Z) for our purpose.
2.3. Surface Property
Butler derived an expression for the surface tension of liquid binary mixture on the basis of the assumption that there exists a mono atomic layer at the surface of a liquid solution, called surface monolayer, as a separate phase that is in thermodynamic equilibrium with the bulk phase [15] [20] . The expression for surface tension, Γ on the basis of Butler equation for binary A-B solution at temperature T is given by
(14)
where, R is universal gas constant.
and
are mole fraction of component “
” in the surface and bulk respectively. Γ1 and Γ2 are the surface tension of the pure component 1 and 2 respectively.
and
(
= 1.2) are partial excess free energy of component “
” in the surface and the bulk respectively. The molar surface area of the component “
” can be computed by using the relation
(15)
where, K (=1.091) is geometrical factor for the liquid alloy, NA is Avogadro’s number, Vi (=
) is the molar volume of the component i.
where,
and
= atomic volume of In at melting point
= atomic volume of Tl at melting point
= melting temperature of In
= melting temperature of Tl
= volume coefficient of In at constant temperature
= volume coefficient of Tl at constant temperature
2.4. Transport Properties
The diffusion coefficients also provide the insight into the mixing behavior of an alloy in microscopic level. The relation between diffusion coefficient and concentration fluctuation is given as [21]
(16)
With
(17)
where, DM is the mutual diffusion coefficient and Did is the intrinsic diffusion coefficient for an ideal mixture; D1 and D2 are the self-diffusivities of pure components A and B respectively.
In term of energy order parameter ω, the diffusion coefficient can be expressed as [22]
(17) If
, then it indicates compound formation and if
then there is tendency of phase separation. For ideal mixing,
approaches 1.
The mixing behavior of liquid alloys at microscopic level can also be understood in terms of viscosity. We have employed the Moelwyn-Hughes equation [16] in order to observe the atomic transport behavior of In-Tl alloy. The Moelwyn-Hughes equation for viscosity of liquid mixture is given as
(18)
where,
is the viscosity of pure component K and for most liquid metals, it can be calculated from Arrhenius type equation [21] at temperature T as
(19)
where,
is constant (in unit of viscosity) and En is the energy of activation of viscous flow for pure metal (in unit of energy per mole).
2.5. Optimization of Free Energy of Mixing (GM), Activity (a), Heat of Mixing (HM), Concentration Fluctuations (Scc(0)), Surface Tension and Viscosity
The model parameters which can be used to optimize thermodynamic, structural and transport properties of the liquid mixture can be obtained by thermodynamic description based on statistical thermodynamics or polynomial expressions. The adjustable coefficients, used in the process, are estimated by least square method which gives an idea to extrapolate into temperature and concentration region in which the direct observed determination is unavailable.
Redlich-Kister polynomial equation gives the composition dependence of excess free energy of mixing and is given by
(20)
with
(21)
The coefficients Kl are the function of the temperature.
The values for the free energy of mixing (GM) of In-Tl liquid alloy at different temperatures can be calculated from Equation (1) by knowing the values of ordering energy parameter (ω) at different temperatures from the relation [23]
(22)
where, ω(T) is the order energy parameter at the temperature 723 K, ω(TK) is order
energy parameter at required temperature TK, and
represents temperature
derivative of order energy parameter which has already been estimated from observed data of entropy of mixing (SM) of In-Tl liquid alloys at 723 K from Hultgren et al. 1973 [24] .
3. Results and Discussion
3.1. Theoretical Investigation of Different Properties of In-Tl Liquid Alloy at 723 K
The value of main input parameters i.e. interchange energy (ω) and
used for the calculation of thermodynamic properties of In-Tl liquid alloy at 723 K were computed from Equations (1) and (9) by using observed values [24] of GM and HM respectively. The best fit parameters were found to be
The positive value of interaction energy indicates that In-Tl is segregating in nature. Using above interchange energy (ω), we computed free energy of mixing (GM), activity (a) entropy of mixing (SM), heat of mixing (HM), concentration fluctuation (Scc(0)), ratio of diffusion coefficient (DM/Did), short-range order parameter (α1), surface tension (Γ1) and viscosity (η) at 723 K using regular solution model.
The plot of observed and computed values of free energy of mixing with respect to the concentration of Tl is shown in Figure 1. The computed and observed values [24] of free energy of mixing are in good agreement in the entire concentration range. Both the computed and observed value of free energy of mixing are minimum at
, which indicates that In-Tl alloys in liquid state at 723 K is symmetric at equiatomic composition.
We computed the activities of In and Tl in the liquid In-Tl alloys using the same energy parameter (ω) as used in Equation (1), and then compared them with observed values as depicted in Figure 2. The plot shows that there is a good agreement between the sets of computed and observed results. Activity coefficient represents the measure of tendency of a component to leave the solution and it provides correlation of the behavior of the systems.
The heat of mixing and entropy of mixing of In-Tl binary liquid alloy at 723K were computed from Equations (9) and (7) respectively. To determine heat of mixing (HM) and entropy of mixing (SM) we need temperature derivatives of energy parameters. The observed values of heat of mixing [24] were utilized to obtain the temperature derivatives by the successive approximation. The best fit
parameters was found to be
. Both computed and observed values
of heat of mixing and entropy of mixing are shown in Figure 3 and Figure 4. Both the computed and observed value of heat of mixing and entropy of mixing are maximum at
. The positive value of heat of mixing indicates the segregating nature of the alloy which is in agreement with the sign of the interaction energy ω. The value of entropy of mixing is positive in whole range of concentration as depicted in Figure 4.
We computed concentration fluctuation in the long wavelength limit Scc(0) using Equation (11) with the help of same energy parameter ω for the consistency as used for the computation of free energy of mixing, activity and heat of mixing. The observed Scc(0) were computed from Equation (12). It was found that the computed and observed values of
at all compositions. This indicates that segregation is favored in the In-Tl alloy at 723 K. The computed values of Scc(0) are in good agreement with the observed values in entire concentration range. The theoretical value is maximum at
(i.e.
) as depicted in Figure 5.
Figure 1. Compositional dependence of (Gm/RT)exp and (Gm/RT)Th of In-Tl liquid alloy at 723 K.
Figure 2. Compositional dependence of lna(i)exp and lna(i)Th of In-Tl liquid alloy at 723 K.
Figure 3. Compositional dependence of (Hm/RT)exp and (Hm/RT)Th of In-Tl liquid alloy at 723 K.
Figure 4. Compositional dependence of (Sm/R)exp and (Sm/R) of In-Tl liquid alloy at 723 K.
Figure 5. Compositional dependence of (Scc(0))exp and (Scc(0))Th of In-Tl liquid alloy at 723 K.
3.2. Theoretical Investigation of In-Tl at 623 K, 723 K, 800 K and 923 K by Optimization Method
By using the best fit value of
and ω(T) at the given temperature T = 723 K
in Equation (22), the value of ω(TK) at temperature TK = 623 K, 800 K, 923 K are estimated and listed in Table 1.
The values of free energy of mixing (GM) of In-Tl liquid alloy at different temperatures (i.e. 623,723,800 and 923 K) computed by using the corresponding values of ω(T) in Equation (22) over the entire range of concentration and then they are used to calculate the corresponding excess free energy of mixing (
) of the alloy at temperature of study by using Equation (20) conjugation with Equations (21) and (22).
The least-square method was used to calculate the parameters involved in Equation (21) and then the optimized coefficients for the alloy are computed which are listed in the Table 2.
We used the parameters i.e. k0, k1, k2 and k3 to obtain partial excess free energy. The partial excess free energy of mixing (
) of the components A (=In) and B (=Tl) were computed and tabulated below using equations [22]
(23)
and
(24)
The partial excess free energy of mixing of both the components In and Tl involved in In-Tl liquid alloy at different temperatures (i.e. = 623 K, 723 K, 800 K, 923 K) were calculated separately over the entire concentration range by Equations (23) and (24) with the help of optimized coefficients. With the help of this optimized partial excess free energy of mixing of both the components in the alloy have been used to calculate the corresponding excess free energy of the alloy at different temperatures over the entire range of concentration from the relation
(25)
The optimized values of excess free energy of mixing for the alloys at different temperatures (i.e. = 623 K, 723 K, 800 K, 923 K) over the entire concentration range are shown in Figure 6.
Now, the activity coefficients (Υi), (i = In or Tl) at different temperature over the entire range of concentration for In or Tl components have been computed
from the relation
, with
, where, ai and ci be the activity
and concentration of the component respectively of In-Tl liquid alloy at corresponding temperature.
The calculated optimized values partial excess free energy of mixing, the corresponding activity coefficients and corresponding activity of both the components involved in In-Tl liquid alloy in the entire concentration range at the temperature T = 623 K, 723 K, 800 K, 923 K are shown in Table 3, Table 4, Table 5 and Table 6.
The concentration fluctuations in long wavelength limit
of the alloy at different temperatures in entire concentration range have been calculated using Equation (14) with the help of the optimized values of the activity of both the components which is shown in Figure 7. The values of concentration fluctuation decreases as the temperature of the alloy increases correspond to the concentration of Tl component. Also,
at entire concentration range of Tl at all temperature of investigation which indicates segregating nature of the alloy.
Figure 6. Compositional dependence of
of In-Tl liquid at different temperatures.
Figure 7. Graph for Scc(0) at different temperatures versus CTl of In-Tl liquid alloy (the ideal values in the graph are taken at 723 K).
Table 1. Estimated values of order energy parameter at different temperatures in In-Tl liquid alloy.
Table 2. Calculated values of optimized coefficients Al, Bl, Cl and Dl (l = 0 to 3) in liquid alloy In-Tl.
Table 3. Optimized values of partial excess free energy of mixing, activity coefficients and activity of both the components involved in In-Tl liquid alloys at 723 K.
Table 4. Optimized values of partial excess free energy of mixing, activity coefficients and activity of both the components involved in In-Tl liquid alloys at 623 K.
Table 5. Optimized values of partial excess free energy of mixing, activity coefficients and activity of both the components involved in In-Tl liquid alloys at 800 K.
Table 6. Optimized values of partial excess free energy of mixing, activity coefficients and activity of both the components involved in In-Tl liquid alloys at 923 K.
The natural logarithms of optimized values of activity of both the components in the alloy at different temperatures in the entire concentration range are shown in Figure 8.
The Warren-Cowley chemical short-range order parameter α1 is used to explain the degree of local arrangement of atoms in the mixture. The chemical short-range order parameter was computed from Equation (13) using coordination number, Z = 10 for liquid the alloy at 723 K. The value of short range order parameter has been found positive in all concentration range which indicates that the alloy is segregating at all compositions. The value of short range order parameter has been found maximum at CTl = 0.5 at 723 K as shown in Figure 9. We have computed chemical short-range order parameter at temperatures 723 K, 623 K, 800 K and 923 K. The value of short range order parameter has been found maximum at CTl = 0.5 at all temperature of study as shown in Figure 9. The plot depicts that α1 is positive in the entire range of concentration showing that α1 in In-Tl is segregating.
The surface tension of the liquid alloys can be computed using Equation (14) conjugation with Equations (15) and (16) using Buttler’s model [15] [20] [25] . The ratio of partial excess Gibbs energy in the bulk and that in the surface can be
expressed as
, β has been taken as 0.83 [26] [27] . We have taken the
surface tension of In-Tl and temperature coefficients for pure In and Tl components from the reference [28] . The surface tension of the pure component at temperature of investigation has been calculated using the relation.
, where,
(= −0.09 mNm−1K−1 for In, −0.08 mN
m−1K−1 for Tl) is temperature coefficient of surface tension, Tm (= 430 K for In and 577 K for Tl) is melting temperature and T = 723 K.
The partial excess free energy of mixing of the pure components was taken from Hultgren et al. [24] . It was found from the analysis that the computed surface tension for In-Tl system at 723 K is less than ideal value (= X1Γ1 + X2Γ2) at all the concentration of In i.e. there is negative departure of surface tension from ideality at 723 K. It was found that surface tension increases with increase in the concentration of component Tl on the alloy as shown in the Figure 10.
To find the surface tension (Γ) of the alloy at different temperature we have used the values of optimized partial excess free energy i.e.
from the Tables 3-6, and using Buttler’s Equation (16) conjugation with Equation (17). The calculated values of surface tension at different temperatures are shown in Figure 10. The surface tension of the alloy in the concentration range of Tl component showed that as the temperature increases, the surface tension decreases.
Viscosity of pure component In and Tl at 723 K are calculated using Equation (19) which is required to compute viscosity of In-Tl alloy at 723 K with the help of the value of the constants E and ηok for the metals [28] . To compute viscosity using Moelwyn-Hughes equation [16] , heat of mixing HM is required which is taken from regular solution model calculations. The result of viscosity at 723 K of the alloy with the ideal values (
) is shown in Figure 11. The graph shows that with the increase in concentration of Tl component in the alloy, the viscosity of the alloy decreases (Figure 11).
The optimized viscosity of the alloys is calculated by the Moelwyn-Hughes Equation (22) conjugation with Equation (23). The values of optimized HM/RT is used in whole concentration range obtained using Equation (11) with the help of
optimized values of ω(T) and
at temperature of study. The values of
parameters on Equation (23) is taken from [28] . The plot exhibits that the viscosity show small negative deviation from ideality at all compositions at 723 K. The computed values of viscosity at different temperatures is as shown in Figure 11.
The diffusion coefficients
of In-Tl liquid alloy at temperatures of study
in the entire concentration range can also be calculated using the optimized data from the relation (19) which are shown in the Figure 12. The
which indicates that there is tendency of phase separation.
Figure 8. Graph for lna at different temperatures versus CTl of In-Tl liquid.
Figure 9. Graph for short range order parameter i.e. α1 at different temperatures versus CTl of In-Tl liquid alloy (the ideal values in the graph are taken at 723 K).
Figure 10. Graph for surface tension i.e. Γ at different temperatures versus CTl of In-Tl liquid alloy (the ideal values in the graph are taken at 723 K).
Figure 11. Graph for viscosity i.e. η at different temperatures versus CTl of In-Tl liquid alloy (the ideal values in the graph are taken at 723 K).
Figure 12. Graph for diffusion coefficient i.e.
at different temperatures versus CTl of In-Tl liquid alloy.
4. Conclusion
From the theoretical investigation following conclusions can be drawn:
1) Symmetry is observed in free energy of mixing, heat of mixing and entropy of mixing at all temperatures of investigations.
2) Activity decrease with the increase in concentration of Tl component in all temperatures of study i.e. 623 K, 723 K, 800 K and 923 K. And, with compare to temperature of study, it decreases as temperatures of investigation increases.
3) At all temperatures of investigation, the surface tension decrease with the increase of bulk concentration of Tl in Tl-In liquid alloy. In context of different temperatures, as temperature of study increases, it decreases.
4) Viscosity increases as the concentration of Tl component increases in the alloy. And, as temperature of study increases viscosity also increases.
5) The diffusion coefficients i.e.
at all compositions which indicates that there is tendency of phase separation
Acknowledgements
Authors are grateful to University Grants Commission, Nepal for providing financial support to pursue this work under Faculty Research Grant Program-2017.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.