1. Introduction
The In-Tl system represents a shape-memory alloy [1] possessing low melting point i.e. 578 K [2] which is less than that of pure In (i.e. 702.76 K) and Tl (i.e. 800 K). Again, In-Tl system has potential applications in metallurgy, thermostats, electric circuits, hydraulic lines, etc. [3]. The mixing properties of In-Tl molten alloys exhibit interesting behavior, i.e. the excess free energy of mixing,
, entropy of mixing,
and free energy of mixing,
are symmetric while heat of mixing,
is asymmetric about C = 0.5 for InTl melts at 723 K [2]. The activity,
shows positive departures from linearity confirming that InTl melts at 723 K is a demixing alloy i.e. segregating alloy [4]. Thus, the structure and thermophysical properties of InTl melts have drawn attention of several investigators [1] [3] [5]. In the past, the thermophysical properties like diffusivity ratio,
,viscosity,
,
and concentration fluctuations,
have been investigated [3] by quasi-lattice theory [4]. In quasi-lattice theory, the size effect i.e. size ratio of the constituent atoms, is coupled with the entropic and enthalpic effects to explain the thermodynamics of molten binary alloys and hence two model parameters namely size ratio and interaction energy are utilized to investigate the thermodynamic properties of the alloys. The simple statistical model [6] has been applied to study the segregating nature of InTl melts at 723 K [5] in terms of thermodynamic properties like
,
,
and
; microscopic properties like
as well as short-range order parameter,
; surface properties like surface concentration,
as well as surface tension,
and transport properties like
and
for InTl melts at 723 K. In simple statistical model, the entropic and enthalpic effects are utilized to explain the thermodynamics of binary molten alloys on assuming the constituent atoms are comparable in size (volume) and hence there is only one model parameter, known as interaction energy. However, the study of the stability function and surface concentration fluctuations
of InTl melts at 723 K is not available in the literature. It should be pointed out that for InTl system, the size factor (
i.e.
) and electronegativity difference (i.e.
) are small to account for the mixing properties of the alloy. The potential applications as shape memory alloys in medical devices, robotics etc., as well as the interesting properties of mixing ensure that InTl system is very crucial for the further theoretical investigation. Thus in present investigation, the theoretical explanation of the thermophysical properties of InTl melts is taken into account in the framework of self-association model. In this model, the size (i.e. volume) of the constituent atoms of the binary alloy is supposed to be comparable. Again, it is assumed that the atoms of the constituent elements form like-atoms clusters or self-associates. Therefore, the present model has two important parameters like ratio of number of two atoms in the self-associates and the interaction energy [4].
The self association model i.e. SAM [4] is employed to describe the thermodynamic functions like
,
,
and
,
and
; and microscopic functions like
well as
and stability i.e. stability function,
for InTl melts at 723 K. Again, Butler model [7] in the framework of SAM is employed to analyze surface properties like
,
and
for InTl melts at 723 K. Further, the Moelwyn-Hughes formula [8] is used in the light of SAM model to explain the transport properties like
and
for InTl melts at 723 K. In the past, SAM model is successfully used to analyze the thermophysical properties of several alloys like GaSn, AlSn, FeSn, CdGa, AlTl, ZnCd, BiSn, GaZn, etc. [9]-[14].
2. Formalism
In the substructure of SAM [4], the required equations are formulated for the bulk properties namely
,
,
,
,
,
as well as
for a binary melt at a specified temperature. Again, the analytical expressions for the surface properties like
,
and
for a binary solution are developed by Butler model [7] in conjunction with SAM model. Further, the expressions for transport properties namely diffusivity,
as well as viscosity,
are deduced using Moelwyn-Hughes model [8] in the framework of SAM model for a binary solution.
2.1. Bulk Properties: Thermodynamic and Microscopic Properties of Binary Liquid Alloys
A binary system, A-B is supposed to be consisted of
as well as
atoms of metal A and metal B respectively where
and
denote the concentrations of constituents A and B, respectively. In SAM model, it is considered that the atoms are situated on equivalent lattice sites with Z nearest neighbors. The effective interactions are short-ranged and exist only between nearest neighbors. Further, consider that sizes of A and B are comparable and there exist like-atom clusters or self-associates like
and
so that
(1)
where i and j = number of atoms in self-associates types A and type B matrices respectively.
In SAM model [4] assuming
[15], the expression for
is given by
(2)
where,
and
(3)
with
and
= universal gas constant and absolute temperature respectively,
, denotes interchange energy while
,
, ratio of number of two atoms in self-associates. It is pointed out that
and
are independent of composition and dependent on the pressure and temperature.
The equation of activity [4] for a binary melt is represented by
(4a)
(4b)
Equations (2) and (4) provide
(5)
(6)
Again,
for a binary melt is represented by [4]
(7)
Equations (2) and (7) provide
(8)
where,
= temperature derivative of
.
The
for a binary melt is expressed as
Equations (2), (8) and (9) provide
(9)
The structural functions like
[16] and
[17], [18] are very useful microscopic functions for acquiring the knowledge of the nature as well as degree of local arrangement of atoms in the binary melts. The
and
are related as [16]
(10)
Expressions (2) and (10) yield
(11)
where,
(12)
For
, expression (12) reduces to
(13)
which represents the ideal value of
.
The activity and
are related as [16]
(14)
The Equation (14) provides the experimental data of
for the binary melts if the experimental data for
and
are utilized.
Again,
in the nearest-neighbor shell for a binary melt [4], [17], [18] and
are related as
with
(15)
where,
denotes the coordination number, which is taken to be 10.
The
is an important parameter used to describe the stability of a binary melt. In terms of
, excess stability function,
[19] is represented by [4], [20], [21]
(16)
2.2. Surface Properties of Binary Liquid Alloys
The properties like
,
as well as surface concentration fluctuations,
describe the characteristics like energy of phase alternation, mechanical behavior, thin film, etc. Hence, the present investigation is concerned with the study of the surface properties of binary melts by Butler model [7].
According to Butler model [7], the surface tension,
of a binary melt is expressed as
(17a)
(17b)
where,
= surface tension of ith component in molten state,
= partial excess free energy at the surface for component i,
= partial excess free energy in the bulk for component i,
= molar surface area of
component, and
and
= surface and bulk compositions respectively. Now,
may be expressed as [9]
(18)
where,
= atomic volume of component i,
= Avogadro’s number.
Now, the expression for
is given by [4]
(19)
where,
= activity coefficient of ith component.
Again,
and
are related as [22]
(20)
with,
[23] where
= surface coordination number and
= bulk coordination number. It is necessary to proclaim that the value of
lies between 0.50 to 0.84 [24]. In present study, the value of
is taken to be 0.818 [25].
Now,
which deals with the arrangement of atoms at the surface of binary solution. The expression for
is given by [26]
(21)
with,
(22)
where,
(23)
and
(24)
The ideal
is represented by [26]
(25)
2.3. Transport Properties of Binary Liquid Alloys
The diffusivity ratio and
is related as [4], [14]
(26)
where,
denotes the mutual diffusion coefficient while
represents the intrinsic diffusion coefficient of an ideal mixture.
Again,
is represented by
(27)
where,
and
= self-diffusion coefficients of pure elements A and B respectively.
The viscosity is helpful to explain the interatomic interactions between the elements of a binary melt. From Moelwyn-Hughes approach [8],
can be expressed as
(28)
where,
= viscosity of the ideal melt which is represented as
(29)
where,
and
= viscosities of pure elements A and B respectively.
Now,
may be represented by [4]
(30)
Equations (28) and (30) provide
(31)
Equation (31) indicates that
depends upon the magnitude of
.
3. Results and Discussions
The analytical expressions are employed to compute the thermophysical properties of InTl melts at 723 K.
3.1. Bulk Properties:
,
,
,
,
,
,
and
for InTl Melts at 723 K
The model parameters like
and
are required for the theoretical investigation of the bulk properties of InTl melts at 723 K. The values of
and
are estimated by the method of successive approximation approach [11], [13], [14] on taking the experimental data of
[2] in Equation (2) for InTl melts at 723 K. The suitable values are
The theoretical data for
of InTl melts at 723 K are presented in Figure 1 along with the correlated experimental data [2] relative to the composition of In which displays well agreement in the whole composition region of In,
to 0.9. The values of
are negative in the whole composition range with minima at
i.e.
(Theory) and −0.5669 (Experiment). Thus, the symmetry in
for InTl melts at 723 K is well explained theoretically.
Figure 1.
and
versus
for InTl melts at 723 K.
The
for InTl melts at 723 K is computed by the relation
(32)
A well agreement is observed between theory and experiment [2] for
of InTl melts at 723 K as depicted in Figure 1. The positive values of
are found in whole region,
to 0.9 having maxima i.e.
(Theory) and 0.1267 (Experiment) both at
. Hence, the symmetry in
for InTl melts at 723 K is explained satisfactorily. Again, the magnitude as well as sign of
reveal that InTl melts at 723 K is a feebly segregating alloy [4].
Now, the activities
and
of In and Tl respectively, computed via Equations (5) and (6) are displayed in Figure 2 along with the corresponding experimental results [2] for InTl melts at 723 K. An excellent harmony is noticed. Again, the small positive departures from Raoult’s law are noticed for
and
in the complete composition range i.e.
to 0.9, which signifies the feebly segregating nature [4] of InTl melts at 723 K. Thus, the well agreement between theory and experiment for
as well as
provides the authenticity of the evaluated values of n and W.
Figure 2.
versus
for InTl melts at 723 K.
For the computations of
and
for InTl melts 723 K via Equations (8) and (9), the values of
and
are required. The evaluated values of
and
by least square method [13], [14] on employing the experimental values of
[2] in expression (8) are −1 × 10−5 K−1 and 0.333 R, respectively.
Figure 3 shows the comparison between the theoretical result of
determined from Equation (8), and the correlated experimental data [2] for InTl melts at 723 K and exhibits excellent agreement in the region,
to 0.9. The peak values are observed at
i.e.
(Theory) and 0.6596 (Experiment). This successfully reproduces the symmetry in
for InTl melts at 723 K.
Again, the theoretical values of
in terms of the composition of In for InTl melts at 723 K, computed from Equation (9) are in reasonable agreement with the corelated experimental data [2] as mentioned in Figure 3. However, a minor discrepancy is noticed i.e.
(Theory) at
and 0.0947 (Experiment) at
, which shows the asymmetric behaviour of InTl melts at 723 K. Further, the values of
are positive in the complete composition region of In, which indicates the demixing nature of InTl melts at 723 K [4].
Figure 3.
and
versus
for InTl melts at 723 K.
The theoretical data of
for InTl melts at 723 K in terms of composition of In, computed via Equations (11) and (12) are mentioned in Figure 4. The experimental data of
[2] and ideal
are also displayed in Figure 4. A well agreement is noticed between theory and experiment with some departures in the region
so that maximum error is about 5.38% at
. In the complete range of composition,
,
having the peak values of 0.315 (Theory) and 0.318 (Experiment) both at
. This advocates the segregating nature of InTl melts at 723 K [4]. Further, the difference between the peak values of
and
i.e.
for InTl melts at 723 K is found to be 0.085 which is very close to 0.0846 as available in the literature [3].
Equation (15) is utilized to compute the values of
relative to the composition of In for InTl melts at 723 K. Figure 4 also displays the
isotherm for InTl melts at 723 K. It should be pointed out that no experimental data of
for InTl melts at 723 K are available in the literature for comparison. It is obvious from Figure 4 that
is positive at each composition throughout the composition range of In form 0.1 to 0.9. The sign and magnitude of
indicate that InTl melt is a feebly segregating alloy.
Figure 4.
and
versus
for InTl melts at 723 K.
Figure 5.
versus
for InTl melts at 723 K.
The excess stability function is useful to study the nature of atomic interactions in a binary melt. The negative value of
signifies the segregating nature of binary melts because in this case,
. Similarly,
becomes positive for ordered alloys due to
[4]. The theoretical values of
are utilized in Equation (16) to explain the concentration dependence of
for InTl melts at 723 K. Figure 5 displays the
isotherm for InTl
melts at 723 K. Obviously, InTl melt at 723 K represents the segregating behaviour because its value is negative in the complete composition region of In i.e.
to 0.9. Again, the magnitude of the negative values of
advocates the weakly segregating nature of InTl melts at 723 K.
3.2. Surface Properties:
,
and
for InTl Melts at
723 K
The simultaneous solution of Equations (17a) and (17b) provides the values of
and
relative to the bulk concentration of In i.e.
(
). The necessary parameters are evaluated from Equations (18), (19) and (20). The values of
for pure In and pure Tl at 723 K are achieved from the expression given below
(33)
where,
.
= melting temperature;
= surface tension of molten component i in pure state at T K;
= surface tension of pure component i at
K;
= derivative of surface tension for ith component with respect to temperature.
The values of input parameters are incorporated in Table 1 [27].
Table 1. Surface tension and temperature coefficient of surface tension for Tl and In [27].
Metal |
Melting point
(K) |
Surface tension (
in mN·m−1) |
Tempr. coefficient of surface tension (
in mN·m−1·K−1) |
In |
429.55 |
556 |
−0.09 |
Tl |
575.65 |
464 |
−0.08 |
The values of
and
at 723 K are found to be 529.590 mN·m−1 and 452.212 mN·m−1, respectively.
The values of
for pure components i.e. In and Tl at 723 K, are determined from the equations given below
(34a)
with
(34b)
where,
= density of the ith component at 723 K;
= density of the ith component at melting temperature;
= atomic mass of the ith component;
= temperature coefficient of the density for the ith component.
The required physical quantities are presented in Table 2 [27].
Table 2. Data for the physical properties of In and Tl [27].
Metal |
M.P. (
in K) |
Atomic mass (
in Kg) |
Density (
in Kg·m−3) |
Tempr. Coefft. of Density (
in Kg·m−3·K−1) |
In |
429.55 |
0.11482 |
7030 |
−0.68 |
Tl |
575.65 |
0.20437 |
11350 |
−1.30 |
The theoretical values of
relative to the bulk concentration of In i.e.
are depicted in Figure 6. It is clear that
rises when the bulk composition of In increases. Again, in the complete composition range of In,
to 0.9,
denotes the negative deviation from linearity which advocates that atoms of Tl segregate on the surface for InTl melts at 723 K.
Figure 7 displays the theoretical values of
in terms of the bulk composition of In for InTl melts at 723 K. Obviously, the value of
increases when concentration of In increases from 0.0 to 1.0. This arises due to the fact that pure In has larger value of surface tension than pure Tl.
Now, equations (21) and (22) are utilized to study the composition dependency of
relative to the bulk composition of In for InTl melts at 723 K. The theoretical result is illustrated in Figure 8 along with the ideal values of
. Obviously,
exhibits positive deviation from
in the entire bulk composition of In. The
has maximum value of 0.3301 at
which indicates the maximum possibility of phase separation in InTl melts at 723 K close to
[13].
Figure 6.
versus
for InTl melts at 723 K.
Figure 7.
versus
for InTl melts at 723 K.
Figure 8.
versus
for InTl melts at 723 K.
3.3. Transport Properties:
and
for InTl Melts at 723 K
Figure 9 displays the theoretical values of
relative to the composition of In for InTl melts at 723 K. For computations, the theoretical values of
for InTl melts at 723 K are used in Equation (26). It is evident from Figure 9 that
is less than 1 in the complete composition region of In with the minimum value of 0.7468 at
, which is very close to 0. 7465 as available in the literature [5]. This indicates that
is greater than
in the region,
to 0.9. The minimum value of
shows that the InTl melt at 723 K is a weakly segregating alloy [28].
The computed data of
are employed to compute the values if
relative to the composition of In for InTl melts at 723 K. The values of
for pure In and pure Tl melts at 723 K are evaluated from the data available in the literature [29], which are found to be 0.9872 (mPa s) and 1.8678 (mPa s) respectively. Figure 10 illustrates the composition dependency of
for InTl melts at 723 K in the region
. Obviously, the viscosity of InTl melts at 723 K declines with the elevation of bulk composition of In.
Figure 9.
versus
for InTl melts at 723 K.
Figure 10.
versus
for InTl melts at 723 K.
4. Conclusions
The thermophysical properties of InTl melts at 723 K are successfully analyzed by SAM model. The experimentally observed symmetries in
,
and
as well as asymmetry in
are successfully analyzed. The composition dependency of
,
,
,
,
,
,
,
,
suggests that InTl melt at 723 K is a feebly segregating alloy. Again, the concentration dependence of
advocates that the atoms of Tl isolate at the surface for InTl melts at 723 K. For InTl melts at 723 K, the surface tension increases when bulk concentration of In increases.
Therefore, SAM model is a reliable and an appropriate model for the theoretical assessment of the demixing (or segregating) liquid alloys at a specified temperature. Further, the model may be employed to study the concentration dependence of the thermophysical properties of segregating molten alloys at various temperatures in the framework of optimization procedure. However, the present work has some limitations, such as scanty of experimental data and the application of self-association model in phase diagram calculations.
Author Contributions
Conceptualization: I.S. Jha and J. Mandal; methodology: R.P. Chaudhary; validation: I.S. Jha; J. Mandal and N.K. Roy; formal analysis: N.K. Roy; investigation: M.M. Hussain; writing-original draft preparation: M.M. Hussain; writing-review and editing: R.P. Chaudhary; visualization: I.S. Jha; supervision: J. Mandal.