Concentration Dependence of Thermodynamic Properties of NaPb Liquid Alloy ()
1. Introduction
Most of the binary alloys are grown from their corresponding liquid phase. In order to improve the quality of growth and to understand the energetic of the formation of alloys, it is important to have as good an understanding as possible of the properties of the liquid alloys. The thermodynamic properties of liquid Na-Pb alloy in molten state exhibits anomalous behaviour as a function of concentration [1]. In present paper, we have used regular associated solution model to understand the thermodynamic properties of Na-Pb alloys in the liquid state at 700 K. Regular associated solution model has been proved to study the thermodynamic and structural properties of weakly, moderately and highly interacting liquid alloys [2-6] by assuming the formation of complex. Such assumptions have been used in different models [7-13]. We have assumed Na and Pb atoms are energetically favoured to form NaPb complex in the molten state.
Theoretical formalism is given in Section 2; Section 3 deals with the numerical result and discussion. Conclusions are provided in Section 4.
2. Theory
Let one mole of binary solution of NaPb alloy comprising of x1 mole of A (=Na) atoms and x2 moles of B (=Pb) atoms. The presence of AB (=NaPb) type complex in the solution results in a depletion of concentration of free atoms of the components of A and B. The liquid solution is thus composed of three species namely free atoms A and B and the complex AB. As a result of associations, the thermodynamic behaviour of the components A and B is governed by the true mole fractions xA and xB rather than the gross mole fraction x1 and x2. Thus it is convenient to operate with two frames of references, one referring to gross mole fractions x1 and x2 and other referring to actual mole fractions of each species (xA, xB and xAB). Further, it is assumed that there are n1 moles of species A, n2 moles of species B and n3 moles of species AB per mole of the binary solution. From the conservation of mass, the two frames of reference can be interrelated as follows:
, ![](https://www.scirp.org/html/6-4800040\ec490437-dee9-49cc-82b8-9e8585732d15.jpg)
and
(1)
, ![](https://www.scirp.org/html/6-4800040\2d6f7339-9892-4673-8b00-a73ae5d5f076.jpg)
and
(2)
Here,
(3)
Now using Equation (3) in Equation (2), we have,
, ![](https://www.scirp.org/html/6-4800040\8de654dd-9280-4522-b827-ddbad8d7ca83.jpg)
(4)
For the sake of convenience one or more of these frames of reference may be used.
Now xA, xB and xAB can be inter-related with each other as follows
![](https://www.scirp.org/html/6-4800040\d442e419-d29b-435e-bed4-e0f5e6bd9da8.jpg)
Using
and
, we get
![](https://www.scirp.org/html/6-4800040\3b27da71-ba30-4840-9cf6-213824d881fd.jpg)
After performing some algebraic operations and rearranging the terms, we obtain
![](https://www.scirp.org/html/6-4800040\47f0365a-b040-4644-adc9-01769a46eabd.jpg)
and
(5a)
Similarly we can obtain
(5b)
In regular associated solution, the gross chemical potentials of components 1 and 2 are equal to the chemical potentials of the monomeric species A and B [2-7]. The activity coefficients
,
and
of monomers and complex can be expressed in terms of pairwise interaction energies through [3-5]
(6a)
(6b)
(6c)
where
,
and
are interaction energies for the species A, B; A, AB and B, AB respectively, T the temperature and R stands for the universal gas constant. The equilibrium constant in a regular associated can be obtained [3,4] as
(7)
Now using the equations listed above the free energy G is given by
(8)
Once the expressions for G are obtained, other thermodynamic and microscopic functions follow readily. Enthalpy of mixing (H), entropy of mixing (SM) are related to G through standard thermodynamic relations
(9)
(10)
In a regular associated solution
and
, where
and
are respective gross activity coefficients of components 1 and 2. Thus
(11a)
and
(11b)
The pairwise interaction energies, the equilibrium constants and the activity coefficients at infinite dilution can be written as [5]
(12a)
(12b)
where
and
are the respective activity coefficients of component A and that of B at zero concentrations.
Solving Equations (6a) and (6b) we obtain
(13a)
(13b)
where a1 and a2 are activities of the components Na and Pb Using Equations (7), (12) and (13) we can derive
(14a)
(14b)
At equiatomic composition, where
, we have
(15)
3. Results and Discussion
We have calculated the concentration of complex, equilibrium constant and pairwise interaction energies following the method employed by Lele and Ramchancrarao [2] and D. Adhikari et al [3,5,6] using Equations (11)-(15) by iterative procedure. The equilibrium constant and interaction energies for the alloy Na-Pb in liquid state at 700 K are found to be K = 0.0118,
= –27.747 kJ·mol–1,
= –21.930 kJ·mol–1 and
= –7.562 kJ·mol–1.The calculated and observed value of integral excess free energy of mixing (
) is in good agreement (Figure 1). The calculated integral excess free energy of mixing is minimum (–12.22 kJ·mol–1) at xNa = 0.6 which almost matches with the observed value [1]. The observed asymmetry in integral excess free energy of mixing is well explained by our theoretical model.
We have observed that if the interaction energies are supposed to be independent of temperature, i.e.
, then H and
so obtained are in very poor agreement with experimental data. This simply suggests importance of the dependence of interaction energies on temperature. On using Equation (10) and observed values of H [11] we have chosen the following values for the given parameters as the best fit values for the heat of formation of NaPb complex
Jmol–1·K–1,
Jmol–1·K–1,
Jmol–1·K–1
and
43530
1200 J·mol–1.The dependence of energy parameters on temperature can be observed from the study of H and SM. It is found that the pairwise interaction energies to be considerably dependent on temperature. It is found from the analysis that the heat of mixing is negative at all concentration. Our theoretical calculation shows that the minimum value of the heat of mixing is –18.02 kJ·mol–1 at
= 0.6. The observed minimum value is also at
= 0.6 [1]. The calculated values are in reasonable agreement with the observed values (Figure 2). The concentration dependence of asymmetry in H is well explained.
We have calculated entropy of mixing of NaPb alloy in liquid state using Equation (10). The calculated values always match in sign with observed values. The calculated values and experimental values are in reasonable agreement at all concentration range (Figure 3).
![](https://www.scirp.org/html/6-4800040\caf7edeb-a379-45d3-9902-2fc55a0963cb.jpg)
Figure 1. Integral excess free energy of mixing (Gxs) versus xNa (concentration of Na) in the liquid NaPb solution (700 K); (––––) theory, (ooo) experiment [1].
![](https://www.scirp.org/html/6-4800040\598bc1db-e1da-4796-91d7-eb72140e358a.jpg)
Figure 2. Heat of mixing (H) versus xNa of liquid NaPb solution (700 K); (––––) theory, (ooo) experiment [1]
![](https://www.scirp.org/html/6-4800040\36c1ed05-f4b5-40ec-ae42-a9041eee9a31.jpg)
Figure 3. Entropy of mixing (SM) versus xNa of liquid NaPb (700 K ) (––––) theory, (ooo) experiment [1]
4. Conclusions
The asymmetry in thermodynamic properties like, integral excess free energy of mixing, heat of mixing and entropy of mixing of NaPb liquid alloys at 700 K is well explained by the present theoretical model. Present theoretical analysis shows that there exist complexes NaPb in NaPb alloy in molten state and the pairwise interaction energies are temperature dependent.