- Email: [email protected]
Entropy of mixing of liquid metal alloys
174
Physica l14B (1982) 174-180 North-Holland Publishing Company
ENTROPY OF MIXING O F LIQUID M E T A L ALLOYS K.N. K H A N N A and P. S I N G H Department of Physics, V.S.S.D. College, Kanpur, India Received 1 March 1982
An analysis of the dependence of entropy of mixing on the concentration of the constituents in an alloy has been presented, based on the hard-sphere reference system. The parameters of the hard-sphere system, as calculated by fitting the excess entropy at equiatomic concentration, provide results comparable with those obtained by the variational method. The alloying properties are found to be sensitive to the atomic composition of the alloy.
1. Introduction The hard-sphere model has been successfully used to calculate several thermodynamical properties of liquid metals and their alloys [1-7]. Recently, the variational method based on the Gibbs-Bogolyubov inequality has been applied by several workers [7-10] to study the entropy of mixing of equiatomic binary alloys. Hoshino and Young [11] have developed a theory for the entropy of mixing of a pseudo binary alloy in which the formation of molecules has been assumed. The basic quantities are the packing fraction, hard-sphere diameters and effective volume of the binary liquid mixtures. These quantities, and hence the thermodynamical properties, change significantly when the atomic volumes of alloying elements are far from equality. In the literature the hard-sphere parameters such as packing fraction and excess volume are determined by minimising the free energy with respect to the hard-sphere diameters as (OF/tgo.i)a,v = 0. The present work examines an alternative procedure of fitting these parameters to match the experimental value of the excess entropy at equiatomic concentration and using an empirical relation to get the parameters for general concentrations. Following the success of our earlier calculations [12] for monovalent
liquid alloy systems, we now intend to examine the applicability of the new fitting procedure to a wider class of homovalent and heterovalent liquid binary mixtures and to elucidate the concentration dependence of various contributions to the entropy of mixing.
2. Hard-sphere description A liquid binary mixture with components of atomic concentrations Cl and c2 comprises ClN hard spheres with diameter o-~ and c2N with diameter o-2. Following Umar et al. [2] the expression for the entropy in the hard-sphere reference system is given by Shs = ~gas -I'- ~7 q'- ~c "q'-~ a ,
(1)
where Sc is the ideal entropy of mixing and $~ corresponds to the mismatch between the hard spheres with different diameters. With O as the atomic volume of the system (V/N), the packing fraction ~7 is given by 77"
'1~ = ~
3
(CIO"1 Q- C20r23) .
The expressions for various S contributions are
0378-4363/82/0000-0000/$02.75 O 1982 North-Holland
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
then given by [6, 7]
NkBS~=2.5+ In [~(2(-m ~ B s0
NkB
T ) 3/2] ,
(2)
- ( c l In c1 + c2 In c2)
/VkBS~= _(so _ 1)(se + 3),
(3) 1
with ~: = - --' 17/
(4)
s~
Nka = clc2(trl - tr2)2{(Xl + x2)[se(~:- 1) - In ~] + 3x,(~ - 1)},
(5)
with x, = (~1 + ~ 2 ) / ( c ~ + c 2 ~ ) ,
175
Thus knowing the value of excess entropy at any concentration, the entropy of mixing can be calculated straightforwardly by including the term Sc- strictly concentration-dependent term. In the binary liquid mixture, formed by mixing constituents of different atomic sizes, the change in the packing fraction and volume on alloying plays a significant role in the determination of thermodynamical properties, such as excess entropy. In the present computation we have considered that the actual volume of an alloy is concentration dependent and is different from the ideal volume of an alloy. The effective volume has been considered as /2 = J-~ideal + ClC2 AJ-~o,
(11)
X2 = 0r10"2(C10"21+ C20"2)/(C10"31+ C20"3)2 .
where The hard-sphere formula for the entropy of mixing is given by AShs = Shs-- C l S l - c2S2,
(6)
where
Thus the entropy of mixing becomes [6] A S m = ASgas + A S , / + So- + S e ,
(7)
where
(12)
The same value of the effective volume has been used in the determination of packing fraction as that used to determine AS~ and S= for the binary mixture, viz. 7/" 77 = ~ (ClOt31+ C2O'3)
ASgas = S g a s - ClSgas,1 -- C2Sgas,2
O1
02
= C I ' O I ' ~ - + C2'02"~ ,
12
(13)
using the packing fraction for the separate liquids as ~/i = 7r~3/6/2i.
and
AS~ = S, - clS,,- c2S~.
(9)
In eq. (7) Sc is the ideal entropy of mixing. Thus, the excess entropy can be written as
NkB
and the fractional change in the volume on mixing is defined as m / 2 = (/2 -- /2ide, al)/~'~ .
Si = Sgas,i + S~,i.
AS e
f~ideal = Cl~'~I + C2~'~2
AS~ + AS, + So-.
(10)
3. Results and discussion
In the present work the atomic volume /2 of the liquid mixture at equiatomic concentration is so fitted as to satisfy eq. (10) with the experi-
176
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
mental excess entropy [13]. That leads to an evaluation of AO0 from eq. (11) and the same value of A/20 is then used to calculate 12 and 77 at different concentrations, and thereby to evaluate the three contributions on the right-hand side of eq. (10) at different concentrations. These are plotted in fig. l(a to f) for six binary liquid systems. While S~ is always positive (see eq. (5)), the present computation predicts that the contributions AS~ and A S ~ to excess entropy can be
"o2
either positive or negative. The term ASh which depends upon the relative magnitude of the packing fractions of the constituent elements, takes the positive or negative values as following: AS~ > 0 for ~7 < ~7c and AS~ < 0 for ~7 > ~1c, where the critical packing fraction is determined by S,c = c]S m + c2S~. Thus AS~ becomes negative in PbSn, AIMg and MgZn and positive in InZn, CdTI and Cdln. It is interesting to note that the plot of AS~ with respect to concentration is
' o6o
Sr
'ol.
?° ~
.o~$
~
.
.
s
..6 ~
,.o
"0~0
- "o|
' o65
l
d~ <3 - "05 I
. .o 4
'oAo o'~
.'~,,,.~+'o~0
-- "05"
I 'ogo -'@7
.oIS 'Olo
- .OI~
'00~
- '10
-
o
.1l
"1
'2
'5
'4 '~ '~, --C I It
"7
'8
'9
Po
Fig. 1. (b)
Fig. 1. (a)
Ca.'t't.
",4
t+
"lg "1o
'lO
1.,8
z
"o8
+'1
<3
I
"o6
'*4 'o2
"e2
.i
'R
"3
'A
"s
CI----~.
Fig. 1. (c)
.&
"7
"8
"q
i.o
o
•+.
.2
"5
'A
.5 Cf--,"
Fig. 1. (d)
"G
"T
"8
.,I
I.o
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
177
Pb S~
.oo'5 "2
A
~.0
t~
L-.o¢ -2
Fig. 1. (e)
Fig. 1. (f)
Fig. 1. (a)-(f)Concentration dependences of AS~, A S ~ and S,, which contribute to the entropy of mixing for six liquid binary mixtures. Contribution of So in CdTI is very small and is not shown in the figure.
asymmetric about the mean concentration value (c1 = c2 = 0.5). This p h e n o m e n o n is closely associated with the volume of the constituent elements. With suffixes chosen such that /21//22 > 1 the maximum occurs at cl < 0 . 5 and the shift becomes more dominant with the increase of volume ratio. This is particularly so in the case of NaCs where the ratio is as large as three [12]. Table II exhibits the behaviour of the shifting in the AS~ curve. The term A S ~ depends upon the relative magnitudes of the atomic volumes of the constituent elements. In most of the cases ASg, remains positive even for heterovalent alloys, contradicting the views of Hafner [3], but becomes negative for large volume contraction such as in A1Mg and MgZn. Here, too, the shifting in the maximum of the ASg~ curve is observed as described earlier but is less affected in comparison to AS,. Examples of numerical results are 12% for AIMg and 4% for CdIn. The mismatch term S~ increases with increasing Io-2O-ll and remains positive at all concentrations for
all alloys. In this case also the shifting depends upon the ratio of the two volumes but is again small as in ASs,. As regards the comparison of the three contributions, we find AS~ > ASg-*> S~ in CdTl and CdIn alloys and AS~ > S~ > ASs,, in InZn, PbSn, A1Mg and MgZn. Thus, as remarked by other workers [3, 10], AS~ dominates in the determination of excess entropy, contrary to pure metals where ASs-* dominates. Our packing fractions [eq. (13)] at different concentrations, as shown in fig. 2, are in reasonable agreement with those derived by Hoshino [6] only for the alloys Table I Parameters used in the calculations Alloy
,(21(a.u.)
/22(a.u.)
7/1
172
AS*/kB
ALMg CdIn InZn CdTl MgZn PbSn
127.99 158.8 188.86 158.8 171.49 222.10
171.49 188.86 111.59 205.21 111.59 194.6
0.527 0.451 0.446 0.451 0.483 0.421
0.483 0.446 0.489 0.375 0,489 0.431
-0.099 0.066 0.144 0.122 -0.3 0.0
178
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
Table II Effective volume, packing fraction and volume of mixing for six binary alloys at equiatomic concentration
Alloy
Effective volume /2(a.u.) Ours Ref. [10]
Packing fraction ~7 Ours Ref. [10]
Excess volume A~(%) Ours Ref. [10]
Shifting in ASh Shifting (%) (~Q1/122)
A(trA+ o'a) (trn + o'a) (%)
AIMg
147.196
146.78
0.5104
0.504
-1.7
-2.7
+20
0.746
Cdln InZn CdTI MgZn PbSn
174.826 149.26 181.67 137.22 208.07
175.78 151.02 183.82 134.82 -
0.4457 0.4649 0.4089 0.5006 0.4262
0.441 0.466 0.412 0.486 -
0.6 -0.6 -0.18 -3.15 -0.13
0.7 0.5 1.0 -6.3 -
+6 -14 +10 -12 -4
0.840 1.69 0.774 1.54 1.14
-0.50 (-0.22)") 0.2 -0.09 0.07 -0.9 -0.036
") From ref. [3].
h a v i n g small excess v o l u m e , b u t d e v i a t e considerably for the alloys possessing large volume c o n t r a c t i o n such as A I M g a n d M g Z n . T h i s is d u e to t h e s e r i o u s a s s u m p t i o n O = Oid,~l c o n s i d e r e d
.~!
~
in t h e d e r i v a t i o n of A*/ in t h e following e q u a t i o n [61: * / - (C1./1 + c2./2) = clc2(*/2- */1)(02C101 "F C202
01)
If w e c o n s i d e r t h e a c t u a l v o l u m e (/2), w e find
~s~
* / - (Cl*/1 + c2./2) = [CLC2(*/2- */1)(O2 - O1) - (c1./1 + c2./2) A O ' ] / O . '4g
•~,s
In view of o u r e m p i r i c a l r e l a t i o n v o l u m e (eq. (11)), eq. (15) b e c o m e s
ea.I~
. . . . . .
(14)
-
(15)
of excess
* / - (c1./1 + c2./2) = c,c2[(*/2 - rh)(O2 - 01) - */iae~lAOo]/g2.
(16)
./~
'42 .41 ./,.
~C4 rFig. 2. Concentration dependence of the packing fraction (~/) of binary liquid alloys. Full and broken curves refer to the values calculated by using eq. (13) (our result) and eq. (14) [6], respectively.
N o w it s h o u l d b e n o t e d t h a t 7/ can p o s s e s s a v a l u e > o r <*/ideal d e p e n d i n g u p o n t h e sign a n d m a g n i t u d e of ADo, as is o b v i o u s f r o m eq. (16), a n d t h u s in o u r analysis t h e p a c k i n g d e n s i t y d e c r e a s e s in I n Z n , CdT1 a n d C d l n a n d i n c r e a s e s in A I M g a n d M g Z n . Thus, t h e c o n d i t i o n 7 / < Cl*/1 + C2./2 for */1 > */2, ~~1 < ~-~2a n d */1 < *12, J~l > O2 s t a t e d b y H o s h i n o [6] is n o t f o l l o w e d in A1Mg a n d M g Z n . T h e v a r i a t i o n b e t w e e n */ a n d Cl is o b s e r v e d t o b e m a x i m u m for C d T I a n d m i n i m u m for C d I n , a n d is p r o p o r t i o n a l t o ]*/1 - */2]. A m o r e i m p o r t a n t f e a t u r e of t h e p a c k i n g f r a c t i o n is t h a t t h e c o n c e n t r a t i o n d e p e n d e n c e of A*/ e x h i b i t s a n o n l i n e a r b e h a v i o u r , as d i s p l a y e d in fig. 3. H e r e
K.N. Khanna and P. Singh / Entropy of mixing o/ liquid metal alloys
179
'olk
6.0
'or2 .o|o -o~j
"OOG.
I
°
'004 .I
~S,,t Cl---~
/
"2
'~
"/~
'5
,¢,,
.7
"8
"*1
I.o
Fig. 4. Concentration dependence of hard spheres of liquid metal alloys; CdTI (open circles), CdIn (black circles), A1Mg (crosses); A, B, C, D, E represent Cd, In, T1, AI, Mg, respectively.
cd. c-&TL .o0~ Fig. 3. C o n c e n t r a t i o n d e p e n d e n c e of A17 of six binary liquid alloys.
again the shift in the maximum value from mean concentration depends upon the ratio of the two volumes. W e now turn to the hard-sphere diameters. In a homovalent alloy like MgZn, the diameter of the lighter element is increased whereas the diameter of the heavier element is reduced, as was also reported by other workers [9, 10]. For the heterovalent alloys our analysis reveals that the hard-sphere diameters of a more electronegative element expands and that of a more electropositive element shrinks except in AIMg, which is in agreement with other theoretical results [3, 10]. T h e change in hard-sphere diameters is observed to be minimum at about 50-50% composition that are reported in table II and increases on either direction. Fig. 4 explains the change in hard-sphere diameters with concentration; only three alloys have been shown in order to avoid the complexity of the figure. The excess volume, an important feature in the determination of thermodynamical properties [14, 15], is plotted at different concentrations in fig. 5. Alloys AIMg and MgZn show a large
l
Ca.I~.
"o61
~ -'°°4 _,oo 8
-'ol~. --'ol~
--'0~.0 --"0~.1
Fig. 5. Concentration dependence of volume of mixing (AD) of six liquid metal alloys. volume contraction, in agreement with other results [10, 16]. This behaviour of AIMg in comparison to alloys CdTI and CdIn possessing approximately same volumes ratio may indicate that hard-sphere mixture of AIMg is more densely packed than the other two alloys. This argument is also supported by the result ~TA~g> cam or ~Tc,m. Again the maximum value of the A/2 curve shifts from the mean concentration in the manner described earlier. In table II some of the quantities at 50-50%
18o
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
c o m p o s i t i o n a r e c o m p a r e d with t h e r e c e n t calc u l a t i o n [10] using a v a r i a t i o n a l t e c h n i q u e . T h e comparison shows reasonable agreement.
4. Conclusion T h e p r e s e n t s t u d y e x p l a i n s t h e c h a n g e in excess e n t r o p y with c o n c e n t r a t i o n , f r o m alloy t o alloy; a d e c r e a s e in excess e n t r o p y in t h e alloy w h e r e t h e c o n t r i b u t i o n s in A S , a n d S~ t e n d to c a n c e l e a c h o t h e r a n d vice versa. T h e r a t i o of t h e v o l u m e s o f t h e c o n s t i t u e n t e l e m e n t s is a significant f e a t u r e t h a t c h a n g e s t h e t h e r m o d y n a m i c a l p r o p e r t i e s of an alloy s h a r p l y . T h e v a r i a t i o n of A/2 a n d A,/ with c o n c e n t r a t i o n a n d t h e shift of t h e m a x i m u m v a l u e f r o m m e a n concentration show the associative tendency bet w e e n u n l i k e a t o m s is s t r o n g at t h e p a r t i c u l a r concentration.
Acknowledgement T h e a u t h o r s a r e g r a t e f u l to P r o f e s s o r D . P . K h a n d e l w a l f o r s t i m u l a t i n g discussions a n d helpful comments.
References [1] I.H. Umar and W.H. Young, J. Phys. F4 (1974) 525. [2] I.H. Umar, I. Yokoyama and W.H. Young, Phil. Mag. 34 (1976) 535. [3] J. Hafner, Phys. Rev. A16 (1977) 351. [4] W.H. Young, Proc. Int. Conf. on Liquid Metals (Bristol, 1976), Inst. Phys. Conf. Ser. 30, p. 1. [5] E.G. Visser, W. van der Lugt and J. Th. M. de Hosson, F10 (1980) 1681. [6] K. Hoshino, J. Phys. F10 (1980) 2157. [7] I. Yokoyama, A. Meyer, M.J. Stott and W.H. Young, Phil. Mag. 35 (1977) 1021. [8] N.W. Ashcroft and D. Stroud, Solid State Phys. 33 (1978) 1. [9] R.N. Singh, J. Phys. F10 (1980) 1411. [10] R.N. Singh and R.B. Choudhary, J. Phys. F l l (1981) 1577. [11] Kozo Hoshino and W.H. Young, J. Phys. F10 (1980) 1365. [12] K.N. Khanna and D.P. Khandelwal, Phys. Stat. Sol. (b) 106 (1981) 715. [13] R. Hultgren, P.D. Desai, D.T. Hawkins, M. Gleiser and K.K. Kelley, Selected Values of the Thermodynamic Properties of Binary Alloys (1973), (Metals Park, Ohio: American Society of Metals). [14] A.B. Bhatia and R.N. Singh, Phys. letters 78 (1980) 460. [15] M. Shimoji, Liquid Metals (Academic Press, New York, 1977). [16] J. Harrier, Inst. Phys. Conf. Ser. No. 30 (Bristol, 1976) p. 107.
Physica l14B (1982) 174-180 North-Holland Publishing Company
ENTROPY OF MIXING O F LIQUID M E T A L ALLOYS K.N. K H A N N A and P. S I N G H Department of Physics, V.S.S.D. College, Kanpur, India Received 1 March 1982
An analysis of the dependence of entropy of mixing on the concentration of the constituents in an alloy has been presented, based on the hard-sphere reference system. The parameters of the hard-sphere system, as calculated by fitting the excess entropy at equiatomic concentration, provide results comparable with those obtained by the variational method. The alloying properties are found to be sensitive to the atomic composition of the alloy.
1. Introduction The hard-sphere model has been successfully used to calculate several thermodynamical properties of liquid metals and their alloys [1-7]. Recently, the variational method based on the Gibbs-Bogolyubov inequality has been applied by several workers [7-10] to study the entropy of mixing of equiatomic binary alloys. Hoshino and Young [11] have developed a theory for the entropy of mixing of a pseudo binary alloy in which the formation of molecules has been assumed. The basic quantities are the packing fraction, hard-sphere diameters and effective volume of the binary liquid mixtures. These quantities, and hence the thermodynamical properties, change significantly when the atomic volumes of alloying elements are far from equality. In the literature the hard-sphere parameters such as packing fraction and excess volume are determined by minimising the free energy with respect to the hard-sphere diameters as (OF/tgo.i)a,v = 0. The present work examines an alternative procedure of fitting these parameters to match the experimental value of the excess entropy at equiatomic concentration and using an empirical relation to get the parameters for general concentrations. Following the success of our earlier calculations [12] for monovalent
liquid alloy systems, we now intend to examine the applicability of the new fitting procedure to a wider class of homovalent and heterovalent liquid binary mixtures and to elucidate the concentration dependence of various contributions to the entropy of mixing.
2. Hard-sphere description A liquid binary mixture with components of atomic concentrations Cl and c2 comprises ClN hard spheres with diameter o-~ and c2N with diameter o-2. Following Umar et al. [2] the expression for the entropy in the hard-sphere reference system is given by Shs = ~gas -I'- ~7 q'- ~c "q'-~ a ,
(1)
where Sc is the ideal entropy of mixing and $~ corresponds to the mismatch between the hard spheres with different diameters. With O as the atomic volume of the system (V/N), the packing fraction ~7 is given by 77"
'1~ = ~
3
(CIO"1 Q- C20r23) .
The expressions for various S contributions are
0378-4363/82/0000-0000/$02.75 O 1982 North-Holland
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
then given by [6, 7]
NkBS~=2.5+ In [~(2(-m ~ B s0
NkB
T ) 3/2] ,
(2)
- ( c l In c1 + c2 In c2)
/VkBS~= _(so _ 1)(se + 3),
(3) 1
with ~: = - --' 17/
(4)
s~
Nka = clc2(trl - tr2)2{(Xl + x2)[se(~:- 1) - In ~] + 3x,(~ - 1)},
(5)
with x, = (~1 + ~ 2 ) / ( c ~ + c 2 ~ ) ,
175
Thus knowing the value of excess entropy at any concentration, the entropy of mixing can be calculated straightforwardly by including the term Sc- strictly concentration-dependent term. In the binary liquid mixture, formed by mixing constituents of different atomic sizes, the change in the packing fraction and volume on alloying plays a significant role in the determination of thermodynamical properties, such as excess entropy. In the present computation we have considered that the actual volume of an alloy is concentration dependent and is different from the ideal volume of an alloy. The effective volume has been considered as /2 = J-~ideal + ClC2 AJ-~o,
(11)
X2 = 0r10"2(C10"21+ C20"2)/(C10"31+ C20"3)2 .
where The hard-sphere formula for the entropy of mixing is given by AShs = Shs-- C l S l - c2S2,
(6)
where
Thus the entropy of mixing becomes [6] A S m = ASgas + A S , / + So- + S e ,
(7)
where
(12)
The same value of the effective volume has been used in the determination of packing fraction as that used to determine AS~ and S= for the binary mixture, viz. 7/" 77 = ~ (ClOt31+ C2O'3)
ASgas = S g a s - ClSgas,1 -- C2Sgas,2
O1
02
= C I ' O I ' ~ - + C2'02"~ ,
12
(13)
using the packing fraction for the separate liquids as ~/i = 7r~3/6/2i.
and
AS~ = S, - clS,,- c2S~.
(9)
In eq. (7) Sc is the ideal entropy of mixing. Thus, the excess entropy can be written as
NkB
and the fractional change in the volume on mixing is defined as m / 2 = (/2 -- /2ide, al)/~'~ .
Si = Sgas,i + S~,i.
AS e
f~ideal = Cl~'~I + C2~'~2
AS~ + AS, + So-.
(10)
3. Results and discussion
In the present work the atomic volume /2 of the liquid mixture at equiatomic concentration is so fitted as to satisfy eq. (10) with the experi-
176
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
mental excess entropy [13]. That leads to an evaluation of AO0 from eq. (11) and the same value of A/20 is then used to calculate 12 and 77 at different concentrations, and thereby to evaluate the three contributions on the right-hand side of eq. (10) at different concentrations. These are plotted in fig. l(a to f) for six binary liquid systems. While S~ is always positive (see eq. (5)), the present computation predicts that the contributions AS~ and A S ~ to excess entropy can be
"o2
either positive or negative. The term ASh which depends upon the relative magnitude of the packing fractions of the constituent elements, takes the positive or negative values as following: AS~ > 0 for ~7 < ~7c and AS~ < 0 for ~7 > ~1c, where the critical packing fraction is determined by S,c = c]S m + c2S~. Thus AS~ becomes negative in PbSn, AIMg and MgZn and positive in InZn, CdTI and Cdln. It is interesting to note that the plot of AS~ with respect to concentration is
' o6o
Sr
'ol.
?° ~
.o~$
~
.
.
s
..6 ~
,.o
"0~0
- "o|
' o65
l
d~ <3 - "05 I
. .o 4
'oAo o'~
.'~,,,.~+'o~0
-- "05"
I 'ogo -'@7
.oIS 'Olo
- .OI~
'00~
- '10
-
o
.1l
"1
'2
'5
'4 '~ '~, --C I It
"7
'8
'9
Po
Fig. 1. (b)
Fig. 1. (a)
Ca.'t't.
",4
t+
"lg "1o
'lO
1.,8
z
"o8
+'1
<3
I
"o6
'*4 'o2
"e2
.i
'R
"3
'A
"s
CI----~.
Fig. 1. (c)
.&
"7
"8
"q
i.o
o
•+.
.2
"5
'A
.5 Cf--,"
Fig. 1. (d)
"G
"T
"8
.,I
I.o
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
177
Pb S~
.oo'5 "2
A
~.0
t~
L-.o¢ -2
Fig. 1. (e)
Fig. 1. (f)
Fig. 1. (a)-(f)Concentration dependences of AS~, A S ~ and S,, which contribute to the entropy of mixing for six liquid binary mixtures. Contribution of So in CdTI is very small and is not shown in the figure.
asymmetric about the mean concentration value (c1 = c2 = 0.5). This p h e n o m e n o n is closely associated with the volume of the constituent elements. With suffixes chosen such that /21//22 > 1 the maximum occurs at cl < 0 . 5 and the shift becomes more dominant with the increase of volume ratio. This is particularly so in the case of NaCs where the ratio is as large as three [12]. Table II exhibits the behaviour of the shifting in the AS~ curve. The term A S ~ depends upon the relative magnitudes of the atomic volumes of the constituent elements. In most of the cases ASg, remains positive even for heterovalent alloys, contradicting the views of Hafner [3], but becomes negative for large volume contraction such as in A1Mg and MgZn. Here, too, the shifting in the maximum of the ASg~ curve is observed as described earlier but is less affected in comparison to AS,. Examples of numerical results are 12% for AIMg and 4% for CdIn. The mismatch term S~ increases with increasing Io-2O-ll and remains positive at all concentrations for
all alloys. In this case also the shifting depends upon the ratio of the two volumes but is again small as in ASs,. As regards the comparison of the three contributions, we find AS~ > ASg-*> S~ in CdTl and CdIn alloys and AS~ > S~ > ASs,, in InZn, PbSn, A1Mg and MgZn. Thus, as remarked by other workers [3, 10], AS~ dominates in the determination of excess entropy, contrary to pure metals where ASs-* dominates. Our packing fractions [eq. (13)] at different concentrations, as shown in fig. 2, are in reasonable agreement with those derived by Hoshino [6] only for the alloys Table I Parameters used in the calculations Alloy
,(21(a.u.)
/22(a.u.)
7/1
172
AS*/kB
ALMg CdIn InZn CdTl MgZn PbSn
127.99 158.8 188.86 158.8 171.49 222.10
171.49 188.86 111.59 205.21 111.59 194.6
0.527 0.451 0.446 0.451 0.483 0.421
0.483 0.446 0.489 0.375 0,489 0.431
-0.099 0.066 0.144 0.122 -0.3 0.0
178
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
Table II Effective volume, packing fraction and volume of mixing for six binary alloys at equiatomic concentration
Alloy
Effective volume /2(a.u.) Ours Ref. [10]
Packing fraction ~7 Ours Ref. [10]
Excess volume A~(%) Ours Ref. [10]
Shifting in ASh Shifting (%) (~Q1/122)
A(trA+ o'a) (trn + o'a) (%)
AIMg
147.196
146.78
0.5104
0.504
-1.7
-2.7
+20
0.746
Cdln InZn CdTI MgZn PbSn
174.826 149.26 181.67 137.22 208.07
175.78 151.02 183.82 134.82 -
0.4457 0.4649 0.4089 0.5006 0.4262
0.441 0.466 0.412 0.486 -
0.6 -0.6 -0.18 -3.15 -0.13
0.7 0.5 1.0 -6.3 -
+6 -14 +10 -12 -4
0.840 1.69 0.774 1.54 1.14
-0.50 (-0.22)") 0.2 -0.09 0.07 -0.9 -0.036
") From ref. [3].
h a v i n g small excess v o l u m e , b u t d e v i a t e considerably for the alloys possessing large volume c o n t r a c t i o n such as A I M g a n d M g Z n . T h i s is d u e to t h e s e r i o u s a s s u m p t i o n O = Oid,~l c o n s i d e r e d
.~!
~
in t h e d e r i v a t i o n of A*/ in t h e following e q u a t i o n [61: * / - (C1./1 + c2./2) = clc2(*/2- */1)(02C101 "F C202
01)
If w e c o n s i d e r t h e a c t u a l v o l u m e (/2), w e find
~s~
* / - (Cl*/1 + c2./2) = [CLC2(*/2- */1)(O2 - O1) - (c1./1 + c2./2) A O ' ] / O . '4g
•~,s
In view of o u r e m p i r i c a l r e l a t i o n v o l u m e (eq. (11)), eq. (15) b e c o m e s
ea.I~
. . . . . .
(14)
-
(15)
of excess
* / - (c1./1 + c2./2) = c,c2[(*/2 - rh)(O2 - 01) - */iae~lAOo]/g2.
(16)
./~
'42 .41 ./,.
~C4 rFig. 2. Concentration dependence of the packing fraction (~/) of binary liquid alloys. Full and broken curves refer to the values calculated by using eq. (13) (our result) and eq. (14) [6], respectively.
N o w it s h o u l d b e n o t e d t h a t 7/ can p o s s e s s a v a l u e > o r <*/ideal d e p e n d i n g u p o n t h e sign a n d m a g n i t u d e of ADo, as is o b v i o u s f r o m eq. (16), a n d t h u s in o u r analysis t h e p a c k i n g d e n s i t y d e c r e a s e s in I n Z n , CdT1 a n d C d l n a n d i n c r e a s e s in A I M g a n d M g Z n . Thus, t h e c o n d i t i o n 7 / < Cl*/1 + C2./2 for */1 > */2, ~~1 < ~-~2a n d */1 < *12, J~l > O2 s t a t e d b y H o s h i n o [6] is n o t f o l l o w e d in A1Mg a n d M g Z n . T h e v a r i a t i o n b e t w e e n */ a n d Cl is o b s e r v e d t o b e m a x i m u m for C d T I a n d m i n i m u m for C d I n , a n d is p r o p o r t i o n a l t o ]*/1 - */2]. A m o r e i m p o r t a n t f e a t u r e of t h e p a c k i n g f r a c t i o n is t h a t t h e c o n c e n t r a t i o n d e p e n d e n c e of A*/ e x h i b i t s a n o n l i n e a r b e h a v i o u r , as d i s p l a y e d in fig. 3. H e r e
K.N. Khanna and P. Singh / Entropy of mixing o/ liquid metal alloys
179
'olk
6.0
'or2 .o|o -o~j
"OOG.
I
°
'004 .I
~S,,t Cl---~
/
"2
'~
"/~
'5
,¢,,
.7
"8
"*1
I.o
Fig. 4. Concentration dependence of hard spheres of liquid metal alloys; CdTI (open circles), CdIn (black circles), A1Mg (crosses); A, B, C, D, E represent Cd, In, T1, AI, Mg, respectively.
cd. c-&TL .o0~ Fig. 3. C o n c e n t r a t i o n d e p e n d e n c e of A17 of six binary liquid alloys.
again the shift in the maximum value from mean concentration depends upon the ratio of the two volumes. W e now turn to the hard-sphere diameters. In a homovalent alloy like MgZn, the diameter of the lighter element is increased whereas the diameter of the heavier element is reduced, as was also reported by other workers [9, 10]. For the heterovalent alloys our analysis reveals that the hard-sphere diameters of a more electronegative element expands and that of a more electropositive element shrinks except in AIMg, which is in agreement with other theoretical results [3, 10]. T h e change in hard-sphere diameters is observed to be minimum at about 50-50% composition that are reported in table II and increases on either direction. Fig. 4 explains the change in hard-sphere diameters with concentration; only three alloys have been shown in order to avoid the complexity of the figure. The excess volume, an important feature in the determination of thermodynamical properties [14, 15], is plotted at different concentrations in fig. 5. Alloys AIMg and MgZn show a large
l
Ca.I~.
"o61
~ -'°°4 _,oo 8
-'ol~. --'ol~
--'0~.0 --"0~.1
Fig. 5. Concentration dependence of volume of mixing (AD) of six liquid metal alloys. volume contraction, in agreement with other results [10, 16]. This behaviour of AIMg in comparison to alloys CdTI and CdIn possessing approximately same volumes ratio may indicate that hard-sphere mixture of AIMg is more densely packed than the other two alloys. This argument is also supported by the result ~TA~g> cam or ~Tc,m. Again the maximum value of the A/2 curve shifts from the mean concentration in the manner described earlier. In table II some of the quantities at 50-50%
18o
K.N. Khanna and P. Singh / Entropy of mixing of liquid metal alloys
c o m p o s i t i o n a r e c o m p a r e d with t h e r e c e n t calc u l a t i o n [10] using a v a r i a t i o n a l t e c h n i q u e . T h e comparison shows reasonable agreement.
4. Conclusion T h e p r e s e n t s t u d y e x p l a i n s t h e c h a n g e in excess e n t r o p y with c o n c e n t r a t i o n , f r o m alloy t o alloy; a d e c r e a s e in excess e n t r o p y in t h e alloy w h e r e t h e c o n t r i b u t i o n s in A S , a n d S~ t e n d to c a n c e l e a c h o t h e r a n d vice versa. T h e r a t i o of t h e v o l u m e s o f t h e c o n s t i t u e n t e l e m e n t s is a significant f e a t u r e t h a t c h a n g e s t h e t h e r m o d y n a m i c a l p r o p e r t i e s of an alloy s h a r p l y . T h e v a r i a t i o n of A/2 a n d A,/ with c o n c e n t r a t i o n a n d t h e shift of t h e m a x i m u m v a l u e f r o m m e a n concentration show the associative tendency bet w e e n u n l i k e a t o m s is s t r o n g at t h e p a r t i c u l a r concentration.
Acknowledgement T h e a u t h o r s a r e g r a t e f u l to P r o f e s s o r D . P . K h a n d e l w a l f o r s t i m u l a t i n g discussions a n d helpful comments.
References [1] I.H. Umar and W.H. Young, J. Phys. F4 (1974) 525. [2] I.H. Umar, I. Yokoyama and W.H. Young, Phil. Mag. 34 (1976) 535. [3] J. Hafner, Phys. Rev. A16 (1977) 351. [4] W.H. Young, Proc. Int. Conf. on Liquid Metals (Bristol, 1976), Inst. Phys. Conf. Ser. 30, p. 1. [5] E.G. Visser, W. van der Lugt and J. Th. M. de Hosson, F10 (1980) 1681. [6] K. Hoshino, J. Phys. F10 (1980) 2157. [7] I. Yokoyama, A. Meyer, M.J. Stott and W.H. Young, Phil. Mag. 35 (1977) 1021. [8] N.W. Ashcroft and D. Stroud, Solid State Phys. 33 (1978) 1. [9] R.N. Singh, J. Phys. F10 (1980) 1411. [10] R.N. Singh and R.B. Choudhary, J. Phys. F l l (1981) 1577. [11] Kozo Hoshino and W.H. Young, J. Phys. F10 (1980) 1365. [12] K.N. Khanna and D.P. Khandelwal, Phys. Stat. Sol. (b) 106 (1981) 715. [13] R. Hultgren, P.D. Desai, D.T. Hawkins, M. Gleiser and K.K. Kelley, Selected Values of the Thermodynamic Properties of Binary Alloys (1973), (Metals Park, Ohio: American Society of Metals). [14] A.B. Bhatia and R.N. Singh, Phys. letters 78 (1980) 460. [15] M. Shimoji, Liquid Metals (Academic Press, New York, 1977). [16] J. Harrier, Inst. Phys. Conf. Ser. No. 30 (Bristol, 1976) p. 107.