Journal of Thermodynamics

Journal of Thermodynamics / 2011 / Article

Research Article | Open Access

Volume 2011 |Article ID 874979 |

W. A. Oates, H. Wang, D. S. Xu, R. Yang, "Evaluation of the Component Chemical Potentials in Analytical Models for Ordered Alloy Phases", Journal of Thermodynamics, vol. 2011, Article ID 874979, 4 pages, 2011.

Evaluation of the Component Chemical Potentials in Analytical Models for Ordered Alloy Phases

Academic Editor: Brian J. Edwards
Received17 Nov 2010
Accepted27 Jan 2011
Published29 Mar 2011


The component chemical potentials in models of solution phases with a fixed number of sites can be evaluated easily when the Helmholtz energy is known as an analytical function of composition. In the case of ordered phases, however, the situation is less straightforward, because the Helmholtz energy is a functional involving internal order parameters. Because of this, the chemical potentials are usually obtained numerically from the calculated integral Helmholtz energy. In this paper, we show how the component chemical potentials can be obtained analytically in ordered phases via the use of virtual cluster chemical potentials. Some examples are given which illustrate the simplicity of the method.

1. Introduction

Chemical potentials for the components in alloy phases are often required: they are useful, for example, in phase diagram calculations. In the case of a binary substitutional alloy model, (A,B), which uses a fixed number of sites, 𝑁𝑆, containing 𝑁𝐴 and 𝑁𝐡 atoms of element A and B, respectively, functions like (πœ•πΉ/πœ•π‘π΄)𝑁𝑆 are not chemical potentials. They are equal to the difference between the component chemical potentials, for example, (πœ•πΉ/πœ•π‘π΅)𝑁𝑆=πœ‡π΅βˆ’πœ‡π΄. This difference is usually referred to as the diffusion potential [1].

This problem with definition does not mean that the individual component chemical potentials are unobtainable in solution phases with a fixed number of sites. They can be obtained from the partial derivative of the calculated Helmholtz energy, F, for example, πœ‡π΄=(πœ•πΉ/πœ•π‘π΄)𝑉,𝑇,𝑁𝐡 in circumstances where F can be expressed as a function of 𝑁𝐴 and 𝑁𝐡. In a completely disordered solution the Helmholtz energy and the component chemical potentials are well defined Ξ”mixπΉπ‘šξ“=𝑅𝑇𝑖π‘₯𝑖lnπ‘₯𝑖,Ξ”mixπœ‡π‘–=𝑅𝑇lnπ‘₯𝑖.(1) Analytical expressions for the chemical potentials can also be derived for models of single lattice phases which take into account deviations from random mixing. Results for the pair quasichemical (Q-C) approximation [2, 3] and for a four-point cluster in the same approximation [4, 5] have been reported.

It is the evaluation of the component chemical potentials in ordered (or antiferromagnet) phases which poses a problem, because there is no longer an explicit relation between F and 𝑁𝐴, 𝑁𝐡. Instead, F is a functional involving internal order parameters. It is because of this that the usual way of obtaining the component chemical potentials has been to numerically differentiate the calculated integral Helmholtz energy.

In this paper, we show how component chemical potentials can be easily obtained in any cluster approximation in either ordered or the single lattice state via the use of virtual chemical potentials (VCPs). VCPs are defined in Section 2. Previously, only point VCPs appear to have been used, but cluster VCPs are also definable and, as we show, are equally useful. The use of cluster VCPs in calculating the equilibrium distribution of clusters and species in partially ordered phases is discussed in Section 3. In Section 4, we show how the component chemical potentials are simply related to the VCPs, and in Section 5, we present the results for some example model calculations.

2. Virtual Chemical Potentials

In their original treatment for calculating the equilibrium distribution of lattice defects in solids, Wagner and Schottky used the law of mass action [6], but later, Schottky gave a more formal treatment of this approach in terms of point VCPs [7]. These point VCPs were used extensively by KrΓΆger [8] in discussing defect equilibria in ionic and semiconductor compounds.

Schottky distinguished between two types of constituent of a solution phase, building units and structural elements. The building units can be regarded as the normal components, while the structural elements are the majority and defect species occurring on the sublattice sites. When a structural element is created, the number of complementary structural elements cannot be kept constant due to the requirement of a definite site ratio; it is, therefore, not possible to assign a true chemical potential to a structural element, nor can they be accessed experimentally. It is possible, however, to define a point VCP for a species A on a sublattice i as πœ‡π΄π‘–=ξƒ©πœ•πΉπœ•π‘π΄(𝑖)ξƒͺ𝑉,𝑇,𝑁1(𝑗),(2) where F is the Helmholtz energy, 𝑁𝐴(𝑖) the number of species or constituents of type A on sublattice i, and 𝑁1(𝑗) to all other point species on all sublattices. We have used the notation πœ‡π΄π‘– here rather than πœ‡π΄(𝑖) since the latter is often used to denote the chemical potential of a component A in a phase i.

The concept of VCPs is readily extended to consider larger clusters than the point. For example, the following can be defined for pair and four-point clusters: πœ‡π΄1𝐡2=ξƒ©πœ•πΉπœ•π‘(12)𝐴𝐡ξƒͺ𝑉,𝑇,𝑁2(12),πœ‡π‘ƒ1𝑄2𝑅3𝑆4=ξƒ©πœ•πΉπœ•π‘(1234)𝑃𝑄𝑅𝑆ξƒͺ𝑉,𝑇,𝑁4(1234),(3) where 𝑁2(12) and 𝑁4(1234) refer to all other pairs and four-point clusters, respectively.

3. Equilibrium Distribution of Species in Ordered Phases

Schottky showed that the value of point VCPs lies in their computational convenience in a modeling context. This can be illustrated by considering a model for an ordered phase comprising two elements, A and B, distributed between two sublattices, 1 and 2. If the sublattices are assumed to be of equal size, then this ordered phase can be represented as (A, B) : (A, B).  

We will consider this phase in the nearest neighbor pair Q-C approximation. If we consider a closed system, Lagrangian multipliers can be assigned to the mass balances πœ†π΄π‘§π‘π΄=2𝑁(12)𝐴𝐴+𝑁(12)𝐴𝐡+𝑁(12)𝐡𝐴,πœ†π΅π‘§π‘π΅=2𝑁(12)𝐡𝐡+𝑁(12)𝐡𝐴+𝑁(12)𝐴𝐡.(4) Minimization of the Lagrangian 𝐿=𝐴+πœ†π΄ξ‚€π‘§π‘π΄ξ‚€π‘βˆ’2(12)𝐴𝐴+𝑁(12)𝐴𝐡+𝑁(12)𝐡𝐴+πœ†π΅ξ‚€π‘§π‘π΅ξ‚€π‘βˆ’2(12)𝐡𝐡+𝑁(12)𝐡𝐴+𝑁(12)𝐴𝐡,(5) followed by the elimination of the Lagrangian multipliers gives the following equilibrium relations between the pair VCPs πœ‡π΄1𝐴2+πœ‡π΅1𝐡2=πœ‡π΄1𝐡2+πœ‡π΅1𝐴2,πœ‡π΄1𝐡2=πœ‡π΅1𝐴2.(6) The solution of these equations, subject to normalization and mass balance constraints, leads to the equilibrium values for the pair probabilities.

In the Q-C approximation the relation between the VCPs and the pair probabilities can be obtained from Ξ”mixπ‘ˆπ‘š=𝑧2𝑖≠𝑗𝑝(12)π‘–π‘—π‘Š(12)𝑖𝑗,Ξ”(7)mixπ‘†π‘š=𝑧2𝑆2βˆ’(π‘§βˆ’1)𝑆1,(8) where π‘Šπ‘–π‘—=πœ€π‘–π‘—βˆ’0.5(πœ€π‘–π‘–+πœ€π‘—π‘—) is the pair exchange energy (the πœ€(12)𝑖𝑗 are the bond energies), z is the coordination number, and 𝑝(12)𝑖𝑗 is the mean pair probability.

The dimensionless pair and point entropies in (8), 𝑆2, and 𝑆1 are given by 𝑆2=βˆ’π‘–π‘—π‘(12)𝑖𝑗ln𝑝(12)𝑖𝑗,𝑆1=βˆ’π΄ξ“π‘–12𝑝𝐴(𝑖)ln𝑝𝐴(𝑖),(9) where 𝑝𝐴(𝑖) is the mean probability or sublattice mole fraction of the species A on sublattice i.

The pair VCPs may then be obtained from the Helmholtz energy minimization (dimensionless) πΉπ‘š=π‘ˆπ‘šβˆ’π‘†π‘š. For example, πœ‡π΅1𝐡2=πœ€(12)𝐡𝐡+ln𝑝(12)π΅π΅βˆ’ξ‚€π‘§βˆ’1𝑧ln𝑝𝐡(1)𝑝𝐡(2).(10) Substitution of such expressions for the VCPs into (6) then leads to the solution for the equilibrium pair distribution.

It should be noted that this use of VCPs is not the only, nor necessarily the most convenient, method to calculate equilibrium cluster distributions. Many using the CVM, for example, use the natural iteration method [9] to calculate these distributions.

4. Component Chemical Potentials

A principal advantage of VCPs lies in their relation to the component chemical potentials. We will first consider the same example as was used in the previous section and then present analogous relations for other examples.

Consider a system which is open to the component B. We lose the mass balance constraint for B and must now consider a Lagrangian based on the grand potential, Ξ©=π΄βˆ’πœ‡π΅π‘π΅: 𝐿=π΄βˆ’πœ†π΄ξ‚€π‘§π‘π΄βˆ’ξ‚€2𝑁(12)𝐴𝐴+𝑁(12)𝐴𝐡+𝑁(12)π΅π΄ξ‚ξ‚βˆ’πœ‡π΅ξ‚€π‘π΅(1)+𝑁𝐡(2)=π΄βˆ’πœ†π΄ξ‚€π‘§π‘π΄βˆ’ξ‚€2𝑁(12)𝐴𝐴+𝑁(12)𝐴𝐡+𝑁(12)π΅π΄ξ‚ξ‚βˆ’πœ‡π΅ξ‚€1𝑧2𝑁(12)𝐡𝐡+𝑁(12)𝐴𝐡+𝑁(12)𝐡𝐴,(11) from which we can obtain, πœ•πΏπœ•π‘(12)𝐡𝐡=πœ‡π΅1𝐡2βˆ’2π‘§πœ‡π΅=0,(12) so that in this case, the component chemical potential is related to just the one pair VCP πœ‡π΅=𝑧2πœ‡π΅1𝐡2.(13) Similar simple expressions are readily obtained for other cluster models. The following lists some examples (n.n. refers to nearest neighbor interactions and n.n.n. to next nearest neighbor interactions).

Four-sublattices, Bragg-Williams (B-W) approxn. Ξ”mixπœ‡π΅=πœ‡1(𝐡),(14) bcc, n.n., Q-C approxn. Ξ”mixπœ‡π΅=4πœ‡2(𝐡)βˆ’7πœ‡1(𝐡),(15) bcc, n.n. & n.n.n, Q-C approxn. Ξ”mixπœ‡π΅=4πœ‡ξ…ž2(𝐡)+3πœ‡2ξ…žξ…ž(𝐡)βˆ’13πœ‡1(𝐡),(16) fcc CVM-T approxn. Ξ”mixπœ‡π΅=2πœ‡4(𝐡)+6πœ‡2(𝐡)βˆ’5πœ‡1(𝐡),(17) bcc CVM-T approxn. Ξ”mixπœ‡π΅=6πœ‡4(𝐡)βˆ’12πœ‡3(𝐡)+4πœ‡ξ…ž2(𝐡)+3πœ‡2ξ…žξ…ž(𝐡)βˆ’πœ‡1(𝐡),(18) where πœ‡4(𝐡)=π‘›ξ“π‘–π‘—π‘˜π‘™1𝑛ln𝑝(π‘–π‘—π‘˜π‘™)𝐡𝐡𝐡𝐡,πœ‡3(𝐡)=π‘›ξ“π‘–π‘—π‘˜1𝑛ln𝑝(π‘–π‘—π‘˜)𝐡𝐡𝐡,πœ‡2(𝐡)=𝑛𝑖𝑗1𝑛ln𝑝(𝑖𝑗)𝐡𝐡,πœ‡1(𝐡)=𝑛𝑖1𝑛ln𝑝𝐡(𝑖).(19) Here, n is the number of different types of cluster or subcluster; for example, n = 4 and n = 2, respectively, for the number of types of n.n. and n.n.n. clusters in the bcc n.n. and n.n.n. Q-C approximation.

It should be noted that there is no relation similar to those given for chemical potentials which permit the analytical calculation of partial molar energies or entropies.

5. Example Calculations

In the examples shown in Figures 1 and 2 for a solution phase A-B, the molar integral Helmholtz mixing energy, Ξ”mixπΉπ‘š, has been calculated from the integral mixing energy and integral mixing entropy. For example, in the Q-C n.n. two-sublattice approximation, the following equations have been used: Ξ”mixπ‘ˆπ‘š=𝑧2ξ“π‘–β‰ π‘—π‘π‘–π‘—π‘Šπ‘–π‘—,Ξ”mixπ‘†π‘šξ“=βˆ’π‘–π‘—π‘π‘–π‘—ln𝑝𝑖𝑗+(π‘§βˆ’1)2𝑖12𝑦𝐴(𝑖)ln𝑦𝐴(𝑖).(20) The chemical potentials shown in the figures have been obtained from (14) and (18). Explicitly, the chemical potentials of A shown in Figure 2 have been obtained from πœ‡π΄=6ln𝑝(1234)π΄π΄π΄π΄βˆ’124ξ‚€ln𝑝(123)𝐴𝐴𝐴+ln𝑝(124)𝐴𝐴𝐴+ln𝑝(134)𝐴𝐴𝐴+ln𝑝(234)𝐴𝐴𝐴+44ξ‚€ln𝑝(13)𝐴𝐴+ln𝑝(14)𝐴𝐴+ln𝑝(23)𝐴𝐴+ln𝑝(24)𝐴𝐴+32ξ‚€ln𝑝(12)𝐴𝐴+ln𝑝(34)π΄π΄ξ‚βˆ’14ξ‚€ln𝑝𝐴(1)+ln𝑝𝐴(2)+ln𝑝𝐴(3)+ln𝑝𝐴(4).(21) Here, the n.n.n. have been taken to involve the sublattices 1-2 and 3-4.

The chemical potentials calculated from the VCPs agree well with those obtained numerically from the independently calculated integral quantity from CVM. Slight difference is due to the n.n. approximation in the VCP calculation, which can obviously be overcome by a straightforward employment of the n.n.n. approximation in the present method.

Besides the simplicity in definition, the use of VCPs also reduces the number of independent variables in the calculation of chemical potentials. In the CVM calculations for n-component alloy, there are necessarily 2𝑛 independent variables, whereas it is significantly decreased to 2n through the definition of VCPs.

6. Conclusion

Component chemical potentials are easily obtained in analytical forms by virtue of cluster VCPs in ordered alloy phases, instead of the usual numerical calculations from the integral Helmholtz energy. The example calculation based on pair quasichemical approximation is compared with the CVM calculation with a four-point cluster in the same approximation, illustrating the simplicity of the method.

Furthermore, the use of VCPs benefits direct comparisons with simulation results, in which systems are always restricted to a fixed number of total sites.


The support of the Ministry of Science and Technology of China under Grant no. 2006CB605104 and the Natural Science Foundation of China under Grant no. 50631030 is gratefully acknowledged. W. A. Oates wishes to acknowledge support for a short visit to the Institute of Metal Research, CAS, Shenyang, China.


  1. F. C. LarchΓ© and J. W. Cahn, β€œThe interactions of composition and stress in crystalline solids,” Acta Metallurgica, vol. 33, no. 3, pp. 331–357, 1985. View at: Google Scholar
  2. G. S. Rushbrooke, Introduction to Statistical Mechanics, Clarendon Press, Oxford, UK, 1949.
  3. E. A. Guggenheim, Mixtures, Clarendon Press, Oxford, UK, 1952.
  4. E. A. Guggenheim and M. L. McGlashan, β€œApproximations relating to strictly regular mixtures,” Molecular Physics, vol. 5, no. 13, pp. 433–445, 1962. View at: Google Scholar
  5. E. A. Guggenheim and M. L. McGlashan, β€œGeneralization of quasiehemical formulas,” The Journal of Chemical Physics, vol. 42, no. 7, pp. 2544–2547, 1965. View at: Google Scholar
  6. C. Wagner and S. Schottky, β€œTheorie der geordneten Mischphasen,” Zeitschrift fΓΌrPhysikalische Chemie B, vol. 11, p. 163, 1930. View at: Google Scholar
  7. W. Schottky, in Halbleiter Probleme, W. Schottky, Ed., vol. 4, pp. 235–268, Fr. Vieweg und Sohn, Braunschweig, Germany, 1958.
  8. F. A. KrΓΆger, The Chemistry of Imperfect Crystals, North-Holland, Amsterdam, The Netherlands, 1973.
  9. T. Kikuchi, β€œSuperposition approximation and natural iteration calculation in cluster-variation method,” Journal of Chemical Physics, vol. 60, p. 1071, 1974. View at: Google Scholar

Copyright Β© 2011 W. A. Oates et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.