wjcmpWorld Journal of Condensed Matter Physics2160-69272160-6919Scientific Research Publishing10.4236/wjcmp.2026.163004wjcmp-153582ArticlePhysicsMathematicsTheoretical Predictions of Thermodynamic Behavior for Cd-Tl Melts at Different TemperaturesRoyNabin Kumar1SinghAnadi Prasad2ChaudharyRajendra Prasad3MandalJagdhar4JhaIndu Shekhar3 University Department of Chemistry, T.M. Bhagalpur University, Bhagalpur, India Marwari College, T.M. Bhagalpur University, Bhagalpur, India Department of Physics, M.M.A.M. Campus, Biratnagar, Nepal University Department of Physics, T.M. Bhagalpur University, Bhagalpur, India
The authors declare no conflicts of interest regarding the publication of this paper.
The molecular interaction volume model is utilized to analyze the composition dependence of thermodynamic properties like Gibbs excess free energy of mixing, free energy of mixing, activity coefficients, activity, microscopic properties such as concentration fluctuations, short-range order parameter and diffusivity ratio of Cd-Tl melts at 750 K. The required interaction parameters are evaluated by Newton-Raphson method using the observed data of activity coefficients of constituent atoms in the infinite dilute solution region. The computed data for excess free energy of mixing, Gibbs free energy of mixing, activity, activity coefficients and concentration fluctuations exhibit excellent harmony with the analogous experimental data for Cd-Tl melts at 750 K. The observed asymmetries in excess free energy of mixing, free energy of mixing and concentration fluctuations for Cd-Tl melts at 750 K are successfully explained. The short-range order parameter and diffusivity relative to the composition for Cd-Tl melts at 750 K are theoretically predicted. The results illustrate the segregating character of Cd-Tl melts at 750 K. Further, the thermodynamic, microscopic and transport properties of Cd-Tl melts at 650 K, 750 K, 850 K, 950 K and 1050 K are predicted theoretically. The segregating behavior of Cd-Tl melts decreases with the elevation in temperature from 650 - 1050 K.
Excess Free Energy of MixingActivity CoefficientsConcentration FluctuationsShort-Range Order ParameterMIVM Model1. Introduction
The Cd-Tl and Cd-Tl-based multicomponent alloys like Sn-Cd-Tl, Bi-Cd-Tl, Pb-Cd-Bi-Tl etc., have potential applications in industries, as semiconductors, as ceramic compounds, as microsensor in technological process, etc. [1][2]. The phase diagram [3] shows that the composition dependence of thermodynamic functions like excess free energy of mixing, G M E , free energy of mixing, G M and enthalpy of mixing, H M are asymmetric while the entropy of mixing, S M is symmetric for Cd-Tl melts at 750 K at x = 0.5. In the entire composition region x = 0.1 to 0.9, G M E and H M have positive values at each composition. The activity of both the components i.e. a C d and a T l illustrates the positive departures from the ideal mixing behaviour. Therefore, Cd-Tl system at 750 K in molten state has a segregating character [4].
The thermodynamic properties and wide range of applications indicate the importance of the Cd-Tl melts for the researchers. Nevertheless, there is scarcity of data related to the study of thermophysical properties of Cd-Tl in the literature. Earlier, some attempts have been made to explore the thermodynamic properties like H M and activity of the components Cd and Tl, and phase diagram by the optimization process on using CALPHAD approach [5] and G M , concentration fluctuations, S c c ( 0 ) , diffusivity, D and viscosity of Cd-Tl melts at 750 K by quasi-lattice theory [6]. Despite these, Cd-Tl melts at 750 K still requires further investigation.
In present work, MIVM model is utilized for the analysis of G M E , G M , activity coefficients, γ υ , activity, S c c ( 0 ) , short-range order parameter, α 1 (SRO), excess stability function, E X S and transport property like diffusivity ratio, D M D i d of Cd-Tl melts at 650 K, 750 K, 850 K, 950 K and 1050 K. The MIVM model [7] is a fluid-based model based on the statistical thermodynamics. It is necessary to point out that MIVM model is successfully employed for the analysis of the thermodynamic properties of several binary alloys (Zn-Bi, Al-Ga, Ni-Pd, Al-In, Zn-Cd, etc.) and multi-component alloys (Au-Sn-Bi, Au-Sn-Zn, Zn-In-Sn, Al-Sn-Zn, Zn-Cu-Sn-In, etc.) [8]-[15].
MIVM model is established by Tao [7] which is a fluid-based model and deduced from statistical thermodynamics, the fundamental notion of non-random interchange of liquid molecules and fluid-phase equilibria. For a binary melt i − j , the molecular excess free energy of mixing is expressed as [7]
where, the compositions of the element i and j represented by x i and x j , A i j and A j i refer to pair-potential energy interaction parameters, respectively and the first coordination numbers of the constituents i and j are respectively Z i and Z j . The pair potential energy interaction parameters may be expressed as [7]
Aji=exp[−εji−εiikT]
where, ε i i , ε j j and ε j i stand for potential energies for the pairs i − i , j − j and i − j respectively, temperature is represented by T and k represents Boltzmann constant. It is noted that ε i j = ε j i .
The partial molar energy and excess molar free energy are related as
where, γ υ ( υ = i , j ) represents the activity coefficient of binary melts.
It is noted that Equation (3) follows the condition that for i ≠ N , δ = 1 ; and i = N , δ = 0 . Again, the variables x j and x N determined from x [ i , N ] with x N = 1 − ∑ j = 1 N − 1 x j . Hence, the activity coefficients of components i and j of a binary melt are respectively, represented by
where, T m i represents the melting temperature in kelvin, Δ H m i , the enthalpy at melting temperature, ρ i = N i V i = 0.6022 V m i represents the molecular number density, r m i and r o i are initial value and first peak value for radial distribution function of metal i in molten state near T m i , R stands for molar gas constant and Z C represents the co-ordination number for closed packed structure which is usually taken to be 12. The radial distances may be represented by
roi=0.918dcov,i
where, d cov represents atomic covalent diameter and σ i the atomic diameter.
In infinite dilute solution region, such as x i or x j → 0 , the activity coefficients of the constituents i and j may be represented by
The S c c ( 0 ) and α 1 are found to be very important microscopic parameters to explore the interatomic interactions in a liquid system. If S c c ( 0 ) < S c c i d ( 0 ) then the system is ordered i.e. dissimilar ( i − j ) atoms or molecules coupled together as the closest neighbors. Similarly, S c c ( 0 ) > S c c i d ( 0 ) describes the segregating character of the binary melt and like (i.e. i − i or j − j ) atoms or molecules link together as the closest neighbors [4].
The G M and S c c ( 0 ) are related as [4][16]
Scc(0)=RT(∂2GM∂xi2)T,P,N−1=RT(∂2GM∂xj2)T,P,N−1
where,
GM=GME+GMid⇒GM=GME+RT[xilnxi+xjlnxj]
Using Equations (1) and (13) in (11), one can obtain
The excess stability function, E X S initiated by Darken [21][22] may be expressed as
EXSRT=Scc(0))−1−(Sccid(0))−1
Obviously, E X S = positive for ordering alloy, E X S = negative for segregating alloy and E X S = 0 for ideal alloy depending upon the conditions that S c c ( 0 ) < S c c i d ( 0 ) , S c c ( 0 ) > S c c i d ( 0 ) , and S c c ( 0 ) = S c c i d ( 0 ) respectively.
The transport property, such as diffusivity is very useful to describe the mixing behavior of a binary solution at the microscopic level. For a binary liquid system, the diffusivity ratio is expressed as [4][23]
DMDid=Sccid(0)Scc(0)
where, D M and D i d represent the mutual diffusion coefficient of the system and intrinsic diffusion coefficient of ideal system respectively. Now, D i d is represented by
Did=xiDj+xjDi
where, D i and D j are self-diffusion coefficients of pure elements i and j respectively.
Equation (20) illustrates that D M / D i d > 1 for ordered system, D M / D i d < 1 for segregating system and D M / D i d = 1 for ideal system.
The analytical equations for G M E , G M , γ μ and a μ ( μ = i , j ) are employed to assess the thermodynamic properties of Cd-Tl melts at 750 K. The necessary parameters of pure Cd and Tl are presented in Table 1 [24]. The coordination numbers Z i and Z j of the constituents of the alloy are determined from Equation (6) which are presented in Table 2. The parameters A j i and A i j are evaluated on solving Equations (8) and (9) simultaneously by Newton-Rapson method using the values of infinite dilute activity coefficients i.e. γ i ∞ (1.932) and γ j ∞ (3.826) [3]. The evaluated values of A j i and A i j are 0.7107 and 0.8815 respectively. The approximated values of A j i and A i j are slightly altered in order to achieve a good harmony between the theory and experiment [3] for G M E of Cd-Tl melts at 750 K. The reasonable values of A j i and A i j i.e. 0.7120 and 1.0990 respectively, are also incorporated in Table 2.
Table 1. Input parameters for the pure metals [24].
Metal,
i
ΔHmi
(KJ/mole)
σi
(×10
−8
cm)
roi
(×10
−8
cm)
Vmi
(cm
3
/mole)
Cd
6.40
3.00
2.54
14.0 [1 + 1.6 × 10
−4
(
T
-594)]
Tl
4.31
3.22
2.74
18.0 [1 + 1.15 × 10
−4
(
T
-575.65)]
Table 2. Values of A i j , A j i , Z i and Z j for the constituents of Cd-Tl melts at required temperature.
i−j
T
(K)
Aij
Aji
Zi
Zj
650
0.6757
1.1151
9.86
9.53
750
0.7120
1.0990
9.71
9.42
Cd-Tl
850
0.7410
1.0869
9.56
9.32
950
0.7648
1.0774
9.41
9.22
1050
0.7846
1.0698
9.27
9.12
The values of G M E / R T calculated from Equation (1) as a function of Cd for Cd-Tl alloys at 750 K show excellent agreement with the corresponding experimental values [3] as presented in Figure 1. In composition region, x C d = 0.1 to 0.9, the value of G M E is positive at each concentration. The computed and experimental [3] values of G M E / R T are in excellent harmony with maximum values at x C d = 0.6 i.e. 0.2478 (Theory) and 0.2504 (Experiment). The RMS error is found to be 0.00533% on employing Equation (22). Hence, the present theoretical investigation successfully reveals the symmetry in G M E / R T and segregating nature of Cd-Tl melts at 750 K.
RMS Error=∑i=1n(xi−yi)2n
where,
x i = theoretical (calculated) values;
y i = experimental values;
n = no. of data points.
Figure 1. G M E / R T vs x C d for Cd-Tl melts at 750 K.
Equation (13) is used to compute the values of G M / R T for Cd-Tl melts at 750 K which are illustrated in Figure 2 along with experimental values [3]. A well agreement is seen. In the whole composition region, x C d = 0.1 to 0.9, G M / R T is negative with minimum values at x C d = 0.4 i.e. G M / R T ( min m ) = − 0.4481 (Theory) and −0.4518 (Experiment). The RMS error computed from Equation (22) is found to be 0.00542%. Thus, the symmetry in G M / R T is successfully described by MIVM model.
Figure 2. G M / R T vs x C d for Cd-Tl melts at 750 K.
The activity coefficients γ C d and γ T l are evaluated from Equations (4) and (5) respectively in terms of the concentration of Cd. The theoretical and experimental [3] data of γ C d and γ T l are in excellent agreement as shown in Figure 3. The RMS errors for γ C d and γ T l computed from Equation (22) are found to be 0.029441% and 0.022264% respectively.
Figure 3. ( γ C d & γ T l ) vs x C d for Cd-Tl melts at 750 K.
The theoretical data of γ C d and γ T l are used to calculate the concentration dependence of a C d and a T l respectively from Equations (10a) and (10b) for Cd-Tl melts at 750 K, which are presented in Figure 4.
Figure 4. ( a C d & a T l ) vs x C d for Cd-Tl melts at 750 K.
A well concord is noticed between theory and experiment [3] with RMS errors for a C d and a T l computed from Equation (22) are found to be 0.008040% and 0.008330% respectively. The positive departures of a C d and a T l from ideal behavior refer to the segregating character of Cd-Tl melts at 750 K. It should be pointed out that the activity is a fortunate thermodynamic function which can be directly measured from the experiment. Therefore, the well harmony between theoretical and experimental data for a C d and a T l validates the reasonable values of A j i and A i j .
Equations (14) and (15) are employed to compute the concentration dependence of S c c ( 0 ) for Cd-Tl melts at 750 K. The theoretical data of S c c ( 0 ) are incorporated in Figure 5 along with S c c i d ( 0 ) and experimental values [3].
Figure 5. S c c ( 0 ) vs x C d for Cd-Tl melts at 750 K.
A well uniformity is noticed between the theory and experiment having the RMS error of 0.0118705% which is computed from Equation (22). It is observed that S c c ( 0 ) > S c c i d ( 0 ) in the composition region, 0.1 ≤ x C d < 0.95 and S c c ( 0 ) exhibits ideal behavior of Cd-Tl melts at 750 K in the regions 0 ≤ x C d ≤ 0.1 and 0.95 ≤ x C d < 1.0 . The maximum values of S c c ( 0 ) are at x C d = 0.60 i.e. S c c ( 0 ) ( max m ) = 0.520 (Theory) and 0.518 (Experiment). Hence, the present investigation reveals that Cd-Tl melt at 750 K possesses segregating character in the region 0.1 ≤ x C d < 0.95 .
Equation (18) is employed to achieve the concentration dependent values of α 1 for Cd-Tl melts at 750 K on using Z = 10. The theoretical data of α 1 in terms of the concentration of Cd are presented in Figure 6. In the full range of composition i.e. x C d = 0.1 to 0.9, α 1 is positive at each composition and having maximum value of 0.0576 at x C d = 0.7 , which reveals the segregating behavior of Cd-Tl melts at 750 K [25].
The concentration dependence of E X S / R T for Cd-Tl alloys at 750 K computed from Equation (19) is illustrated in Figure 7. Obviously, E X S / R T is negative at each composition of Cd in the concentration region x C d = 0.10 to 0.90 for Cd-Tl melts at 750 K. This exhibits the segregating nature of Cd-Tl melts at 750 K [22][26].
Figure 6. α 1 vs x C d for Cd-Tl melts at 750 K.
Figure 7. E X S / R T vs x C d for Cd-Tl melts at 750 K.
The diffusivity ratio, D M / D i d as a function of the composition of Cd for Cd-Tl melts at 750 K is calculated from Equation (20) and the theoretical data of D M / D i d are depicted in Figure 8. It is noticed that the diffusivity ratio, D M / D i d has minimum value of 0.4554 at 750 K at x C d = 0.7 . This confirms the segregation in Cd-Tl melts at 750 K [4].
Figure 8. D M / D i d vs x C d for Cd-Tl melts at 750 K.
For the theoretical analysis of the mixing behavior of Cd-Tl melts at different temperatures, the needed input parameters are presented in Table 1 and Table 2. The values of A j i and A i j at desired temperature, T K are determined via the following equations [7].
The computed data of G M E / R T and G M / R T for Cd-Tl melts assessed via Equations (1) and (13) at 650 K, 750 K, 850 K, 950 K and 1050 K are presented in Figure 9 and Figure 10 respectively. Figure 9 illustrates that the segregating nature of Cd-Tl melts degrades when temperature is upgraded from 650 K to 1050 K. The maximum values of G M E / R T are 0.2767 (at 650 K), 0.2478 (at 750 K), 0.2234 (at 850 K) all x C d = 0.6 while G M E / R T ( max m ) = 0.2033 (at 950 K) and 0.1863 (at 1050 K) both at x C d = 0.5 . Thus, Cd-Tl melts exhibit asymmetric behavior at 650 K, 750 K and 850 K, and symmetric character at 950 K and 1050 K relative to the composition of Cd. Again, it is noticeable from Figure 10 that the negative values of G M / R T for Cd-Tl melts increase due to rise in temperature from 650 K to 1050 K. The values of G M / R T are minimum at x C d = 0.4 for temperatures 650 K and 750 K i.e. −0.4256 and −0.4481, respectively while G M / R T shows minima at x C d = 0.5 for temperatures 850 K, 950 K and 1050 K i.e. −0.4703, −0.4899 and −0.5068 respectively in terms of x C d . Thus, Cd-Tl melts indicate asymmetry at 650 K and 750 K and symmetry at 850 K, 950 K and 1050 K with respect to x C d .
Figure 9. G M E / R T vs x C d for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K & 1050 K.
Figure 10. G M / R T vs x C d for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K & 1050 K.
The thermodynamic activities of Cd and Tl i.e. a C d and a T l are computed from Equations (10a) and (10b) respectively on applying the theoretical data of γ C d and γ T l determined from Equations (4) and (5) for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K and 1050 K relative to x C d . The theoretical results for a C d and a T l at the temperatures ranging from 650 K to 1050 K are shown in Figure 11. The results reveal the decrease in segregation for rise in temperature from 650 K to 1050 K for Cd-Tl melts. However, the decrease in the segregation is weak in the temperature region 650 - 1050 K.
Figure 11. ( a C d & a T l ) vs x C d for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K & 1050 K.
The temperature dependence of S c c ( 0 ) for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K and 1050 K relative to composition Cd is presented in Figure 12. For the calculation of S c c ( 0 ) , Equations (14) and (15) are applied. The S c c ( 0 ) - x C d isotherms indicate the degradation in segregating character of Cd-Tl melts for the rise in temperature from 650 K to 1050 K. The peak value of S c c ( 0 ) is noticed at x C d = 0.6 i.e. 0.618, 0.520, 0.466 and 0.427 respectively at 650 K, 750 K, 850 K, 950 K while a minor shift in peak value of S c c ( 0 ) is observed at 1050 K i.e. S c c ( 0 ) = 0.405 at x C d = 0.53 .
Figure 12. S c c ( 0 ) vs x C d for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K & 1050 K.
The α 1 - x C d isotherms for Cd-Tl melts are depicted in Figure 13. For this, the theoretical data of α 1 are computed with the help of Equation (18) at temperatures 650 - 1050 K. Obviously, the segregation in Cd-Tl melts goes on decreasing for the rise in temperature from 650 K to 1050 K. Again, the magnitude of α 1 is small which illustrates that Cd-Tl melts is a weakly segregating system [4][26] in the temperature range 650 - 1050 K.
Figure 13. α 1 vs x C d for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K & 1050 K.
The excess stability function of several binary alloys at a certain temperature has been discussed earlier [27][28]. Nonetheless, the limited data is available regarding the theoretical exploration of E X S at different temperatures [22][26], [29] on employing the optimization procedure. In present study, MIVM model is applied to explore the stability of Cd-Tl melts at various temperatures.
The variation of E X S / R T with respect to composition of Cd at 650 K, 750 K, 850 K, 950 K and 1050 K for Cd-Tl melts determined via Equation (19) is presented in Figure 14. Obviously, E X S / R T exhibits negative value at each composition in the region x C d = 0.10 to 0.90 at each temperature under consideration. Again, the negative values of E X S / R T declines with the upgradation in temperature from 650 K to 1050 K. This is a sign of decrease in the segregation of Cd-Tl melts [22][26] from 650 K to 1050 K.
The D M / D i d relative to the concentration of Cd for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K and 1050 K are computed using Equation (20) and the theoretical data are depicted in Figure 15. It is noticed that D M / D i d < 1 at the temperatures under consideration which confirms the segregating character of Cd-Tl melts. Again, D M / D i d has minimum values at x C d = 0.7 i.e. 0.3679 and 0.4554 respectively for 650 K and 750 K and again, has minimum values at x C d = 0.6 i.e. 0.5145, 0.5621 and 0.6022 for 850 K, 950 K and 1050 K respectively. This indicates lowering of segregating character of Cd-Tl melts at the undertaken temperatures.
Figure 14. E X S / R T vs x C d for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K & 1050 K.
Figure 15. D M / D i d vs x C d for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K & 1050 K.
4. Conclusions
In present work, the concentration dependent properties such as G M E , G M activity coefficient, activity, S c c ( 0 ) , α 1 , E X S and D M / D i d for Cd-Tl melts at 750 K have been successfully explored. The results exhibit the segregating nature of Cd-Tl melts at 750 K. The asymmetries in G M E , G M and S c c ( 0 ) for Cd-Tl melts at 750 K are successfully explained. Again, MIVM model is employed for the theoretical predictions of G M E , G M activity, S c c ( 0 ) , α 1 , E X S and D M / D i d for Cd-Tl melts at 650 K, 750 K, 850 K, 950 K and 1050 K. The results indicate that the segregating character of Cd-Tl melts degrades when temperature rises from 650 K to 1050 K. Further, Cd-Tl melt is a feebly segregating alloy in the temperature range 650 - 1050 K.
Therefore, MIVM model is a relevant and reliable model to explore the thermodynamic behavior of segregating alloys at different temperatures.
Author Contributions
Conceptualization, I. S. Jha and J. Mandal; methodology, N. K. Roy; validation, I. S. Jha, J. Mandal and R. P. Chaudhary; formal analysis; R. P. Chaudhary; investigation, A. P. Singh; writing original draft preparation, A. P. Singh; writing-review and editing, N. K. Roy; visualization, I. S. Jha; supervision, J. Mandal. All authors have read and agreed to published version of the manuscript.
ReferencesDinsdale, A.T. (1991) SGTE Data for Pure Elements. Calphad, 15, 317-425. https://doi.org/10.1016/0364-5916(91)90030-n 10.1016/0364-5916(91)90030-nhttps://doi.org/10.1016/0364-5916(91)90030-nDinsdale, A.T.1991SGTE Data for Pure ElementsCalphad59169110.1016/0364-5916(91)90030-nPelton, A.D. (1988) The Cd-Na (Cadmium-Sodium) System. Bulletin of Alloy Phase Diagrams, 9, 41-46. https://doi.org/10.1007/bf02877458 10.1007/bf02877458https://doi.org/10.1007/bf02877458Pelton, A.D.1988The Cd-Na (Cadmium-Sodium) SystemBulletin of Alloy Phase Diagrams910.1007/bf02877458Hultgren, R., Desai, P.D., Hawkins, D.T., Gleiser, M. and Kelley, K.K. (1973) Selected Values of the Thermodynamic Properties of Binary Alloys. Metal Park (Ohio), ASM.Hultgren, R.Desai, P.D.Hawkins, D.T.Gleiser, M.Kelley, K.K.1973Selected Values of the Thermodynamic Properties of Binary AlloysMetal Park (Ohio)Singh, R.N. and Sommer, F. (1997) Segregation and Immiscibility in Liquid Binary Alloys. Reports on Progress in Physics, 60, 57-150. https://doi.org/10.1088/0034-4885/60/1/003 10.1088/0034-4885/60/1/003https://doi.org/10.1088/0034-4885/60/1/003Singh, R.N.Sommer, F.1997Segregation and Immiscibility in Liquid Binary AlloysReports on Progress in Physics6010.1088/0034-4885/60/1/003Liu, Y.J., Yu, D. and Ge, Y. (2009) Thermodynamic Optimization of the Cd-Na and Cd-Tl Binary Systems. Journal of Alloys and Compounds, 473, 60-64. https://doi.org/10.1016/j.jallcom.2008.05.106 10.1016/j.jallcom.2008.05.106https://doi.org/10.1016/j.jallcom.2008.05.106Liu, Y.J.Yu, D.Ge, Y.2009Thermodynamic Optimization of the Cd-Na and Cd-Tl Binary SystemsJournal of Alloys and Compounds47310.1016/j.jallcom.2008.05.106Awe, O.E. (2018) Thermodynamic Investigation of Thermophysical Properties of Thallium-Based Liquid Alloys. Physics and Chemistry of Liquids, 57, 296-310. https://doi.org/10.1080/00319104.2018.1443453 10.1080/00319104.2018.1443453https://doi.org/10.1080/00319104.2018.1443453Awe, O.E.2018Thermodynamic Investigation of Thermophysical Properties of Thallium-Based Liquid AlloysPhysics and Chemistry of Liquids5710.1080/00319104.2018.1443453Tao, D.P. (2000) A New Model of Thermodynamics of Liquid Mixtures and Its Application to Liquid Alloys. ThermochimicaActa, 363, 105-113. https://doi.org/10.1016/s0040-6031(00)00603-1 10.1016/s0040-6031(00)00603-1https://doi.org/10.1016/s0040-6031(00)00603-1Tao, D.P.2000A New Model of Thermodynamics of Liquid Mixtures and Its Application to Liquid AlloysThermochimica Acta60310010.1016/s0040-6031(00)00603-1Ding, G.H., Liang, L.M., Meng, P., Wang, L., Wang, Y.X., Wang, Y.J., et al. (2019) Application of the Molecular Interaction Volume Model to the Assessment of Thermodynamic Properties of Liquid Zn-Bi Alloys. MetallicMaterials, 57, 55-60. https://doi.org/10.4149/km_2019_1_55 10.4149/km_2019_1_55https://doi.org/10.4149/km_2019_1_55Ding, G.H.Liang, L.M.Meng, P.Wang, L.Wang, Y.X.Wang, Y.J.2019Application of the Molecular Interaction Volume Model to the Assessment of Thermodynamic Properties of Liquid Zn-Bi AlloysMetallic Materials5710.4149/km_2019_1_55Yang, H.W., Tao, D.P. and Zhou, Z.H. (2008) Prediction of the Mixing Enthalpies of Binary Liquid Alloys by Molecular Interaction Volume Model. ActaMetallurgicaSinica ( EnglishLetters), 21, 336-340. https://doi.org/10.1016/s1006-7191(08)60056-3. 10.1016/s1006-7191(08)60056-3https://doi.org/10.1016/s1006-7191(08)60056-3Yang, H.W.Tao, D.P.Zhou, Z.H.2008Prediction of the Mixing Enthalpies of Binary Liquid Alloys by Molecular Interaction Volume ModelActa Metallurgica Sinica (English Letters)71910810.1016/s1006-7191(08)60056-3Liang, L.M., Ding, G.H., Wang, Y.X. and Liu, Y. (2021) Thermodynamic Properties of Liquid Zn-Cd Alloys Investigated by the Molecular Interaction Volume Model. PhysicsandChemistryofLiquids, 59, 706-715. https://doi.org/10.1080/00319104.2020.1808655 10.1080/00319104.2020.1808655https://doi.org/10.1080/00319104.2020.1808655Liang, L.M.Ding, G.H.Wang, Y.X.Liu, Y.2021Thermodynamic Properties of Liquid Zn-Cd Alloys Investigated by the Molecular Interaction Volume ModelPhysics and Chemistry of Liquids5910.1080/00319104.2020.1808655Awe, O.E. and Oshakuade, O.M. (2017) Theoretical Prediction of Thermodynamic Activities of Liquid Au-Sn-X (X=Bi, Sb, Zn) Solder Systems. Physica B: Condensed Matter, 507, 84-94. https://doi.org/10.1016/j.physb.2016.11.035 10.1016/j.physb.2016.11.035https://doi.org/10.1016/j.physb.2016.11.035Awe, O.E.Oshakuade, O.M.Bi, S2017Theoretical Prediction of Thermodynamic Activities of Liquid Au-Sn-X (X=Bi, Sb, Zn) Solder SystemsPhysica B: Condensed Matter50710.1016/j.physb.2016.11.035Roy, N.K., Vrind, N., Sahu, R.K., Mandal, J. and Jha, I.S. (2025) Prediction of the Alloying Behaviour of Ni-Pd Melts by Molecular Interaction Volume Model. IOP Conference Series: Materials Science and Engineering, 1334, Article 012005. https://doi.org/10.1088/1757-899x/1334/1/012005 10.1088/1757-899x/1334/1/012005https://doi.org/10.1088/1757-899x/1334/1/012005Roy, N.K.Vrind, N.Sahu, R.K.Mandal, J.Jha, I.S.2025Prediction of the Alloying Behaviour of Ni-Pd Melts by Molecular Interaction Volume ModelIOP Conference Series: Materials Science and Engineering133401200510.1088/1757-899x/1334/1/012005Sah, S.K., Jha, I.S. and Koirala, I. (2023) Theoretical Examination of Some Thermodynamic Properties in In-Bi-Sn Liquid Alloy and Its Sub-Binary Systems. Journal of Electronic Materials, 52, 6350-6359. https://doi.org/10.1007/s11664-023-10571-y 10.1007/s11664-023-10571-yhttps://doi.org/10.1007/s11664-023-10571-ySah, S.K.Jha, I.S.Koirala, I.2023Theoretical Examination of Some Thermodynamic Properties in In-Bi-Sn Liquid Alloy and Its Sub-Binary SystemsJournal of Electronic Materials5210.1007/s11664-023-10571-yOdusote, Y.A., Popoola, A.I., Adedayo, K.D. and Ogunjo, S.T. (2017) Thermodynamic Properties of Al in Ternary Lead-Free Solder Al-Sn-Zn Alloys. Materials Science- Poland, 35, 583-593. https://doi.org/10.1515/msp-2017-0080 10.1515/msp-2017-0080https://doi.org/10.1515/msp-2017-0080Odusote, Y.A.Popoola, A.I.Adedayo, K.D.Ogunjo, S.T.2017Thermodynamic Properties of Al in Ternary Lead-Free Solder Al-Sn-Zn AlloysMaterials Science-Poland3510.1515/msp-2017-0080Sah, S.K., Jha, I.S. and Koirala, I. (2024) Thermodynamic Activity in Zn-Cu-Sn-In Liquid Solder Alloys: A Comprehensive Analysis Using the Molecular Interaction Volume Model. Welding International, 38, 265-276. https://doi.org/10.1080/09507116.2024.2325151 10.1080/09507116.2024.2325151https://doi.org/10.1080/09507116.2024.2325151Sah, S.K.Jha, I.S.Koirala, I.2024Thermodynamic Activity in Zn-Cu-Sn-In Liquid Solder Alloys: A Comprehensive Analysis Using the Molecular Interaction Volume ModelWelding International3810.1080/09507116.2024.2325151Bhatia, A.B. and Thornton, D.E. (1970) Structural Aspects of the Electrical Resistivity of Binary Alloys. Physical Review B, 2, 3004-3012. https://doi.org/10.1103/physrevb.2.3004 10.1103/physrevb.2.3004https://doi.org/10.1103/physrevb.2.3004Bhatia, A.B.Thornton, D.E.1970Structural Aspects of the Electrical Resistivity of Binary AlloysPhysical Review B210.1103/physrevb.2.3004Bhatia, A.B. and Hargrove, W.H. (1974) Concentration Fluctuations and Thermodynamic Properties of Some Compound Forming Binary Molten Systems. Physical Review B, 10, 3186-3196. https://doi.org/10.1103/physrevb.10.3186 10.1103/physrevb.10.3186https://doi.org/10.1103/physrevb.10.3186Bhatia, A.B.Hargrove, W.H.1974Concentration Fluctuations and Thermodynamic Properties of Some Compound Forming Binary Molten SystemsPhysical Review B1010.1103/physrevb.10.3186Warren, B.E. (1969) X-Ray Diffraction. Addison-Wesley Publishing Company.Warren, B.E.1969X-Ray DiffractionCowley, J.M. (1950) An Approximate Theory of Order in Alloys. Physical Review, 77, 669-675. https://doi.org/10.1103/physrev.77.669 10.1103/physrev.77.669https://doi.org/10.1103/physrev.77.669Cowley, J.M.1950An Approximate Theory of Order in AlloysPhysical Review7710.1103/physrev.77.669Singh, R.N. (1987) Short-Range Order and Concentration Fluctuations in Binary Molten Alloys. Canadian Journal of Physics, 65, 309-325. https://doi.org/10.1139/p87-038 10.1139/p87-038https://doi.org/10.1139/p87-038Singh, R.N.1987Short-Range Order and Concentration Fluctuations in Binary Molten AlloysCanadian Journal of Physics6510.1139/p87-038Darken, L.S. (1967) Thermodynamics of Binary Metallic Solutions. Transactions of the Metallurgical Society of AIME, 239, 80-89.Darken, L.S.1967Thermodynamics of Binary Metallic SolutionsTransactions of the Metallurgical Society of AIME239Singh, A.P., Chaudhary, R.P., Kumari, R., Jha, I.S. and Mandal, J. (2022) Stability of NaK Liquid Alloys at Different Temperatures. Materials Today: Proceedings, 66, 2251-2258. https://doi.org/10.1016/j.matpr.2022.06.097 10.1016/j.matpr.2022.06.097https://doi.org/10.1016/j.matpr.2022.06.097Singh, A.P.Chaudhary, R.P.Kumari, R.Jha, I.S.Mandal, J.2022Stability of NaK Liquid Alloys at Different TemperaturesMaterials Today: Proceedings6610.1016/j.matpr.2022.06.097Darken, L.S. (1948) Diffusion, Mobility and Their Interrelation through Free Energy in Binary Metallic Systems. Transactions of AIME, 175, 184-201.Darken, L.S.Diffusion, M1948Diffusion, Mobility and Their Interrelation through Free Energy in Binary Metallic SystemsTransactions of AIME175Iida, T. and Guthrie, R.I.L. (1988) The Physical Properties of Liquid Metals. Oxford University Press.Iida, T.Guthrie, R.I.L.1988The Physical Properties of Liquid MetalsChaudhary, R.P., Mandal, J. and Jha, I.S. (2022) Concentration Dependence of the Bulk and Surface Properties of Liquid AlZn Alloys at 1000 K. MaterialsToday: Proceedings, 59, 397-404. https://doi.org/10.1016/j.matpr.2021.11.358 10.1016/j.matpr.2021.11.358https://doi.org/10.1016/j.matpr.2021.11.358Chaudhary, R.P.Mandal, J.Jha, I.S.2022Concentration Dependence of the Bulk and Surface Properties of Liquid AlZn Alloys at 1000 KMaterials Today: Proceedings5910.1016/j.matpr.2021.11.358Singh, A.P., Chaudhary, R.P., Hussain, M.M., Mandal, J. and Jha, I.S. (2025) Theoretical Study of the Stability of CdSn Liquid Alloy at Different Temperatures. Diffusion Foundations and Materials Applications, 37, 37-51. https://doi.org/10.4028/p-efgqf5 10.4028/p-efgqf5https://doi.org/10.4028/p-efgqf5Singh, A.P.Chaudhary, R.P.Hussain, M.M.Mandal, J.Jha, I.S.2025Theoretical Study of the Stability of CdSn Liquid Alloy at Different TemperaturesDiffusion Foundations and Materials Applications3710.4028/p-efgqf5Jha, I.S., Singh, R.N., Srivastava, P.L. and Mitra, N.R. (1990) Stability of HgNa and HgK Liquid Alloys. Philosophical Magazine B, 61, 15-24. https://doi.org/10.1080/13642819008208649 10.1080/13642819008208649https://doi.org/10.1080/13642819008208649Jha, I.S.Singh, R.N.Srivastava, P.L.Mitra, N.R.1990Stability of HgNa and HgK Liquid AlloysPhilosophical Magazine B6110.1080/13642819008208649Abodunrin, O.W. and Ajayi, A.A. (2022) Calculation of Excess Stability Functions of Four Binary Alloys. Jordan Journal of Physics, 14, 193-200.Abodunrin, O.W.Ajayi, A.A.2022Calculation of Excess Stability Functions of Four Binary AlloysJordan Journal of Physics14Shrestha, G.K., Singh, B.K., Jha, I.S., Singh, B.P. and Adhikari, D. (2017) Optimization Method for the Study of the Properties of Al-Sn Binary Liquid Alloys. PhysicaB: CondensedMatter, 514, 1-7. https://doi.org/10.1016/j.physb.2017.03.005 10.1016/j.physb.2017.03.005https://doi.org/10.1016/j.physb.2017.03.005Shrestha, G.K.Singh, B.K.Jha, I.S.Singh, B.P.Adhikari, D.2017Optimization Method for the Study of the Properties of Al-Sn Binary Liquid AlloysPhysica B: Condensed Matter51410.1016/j.physb.2017.03.005