arxiv: v1 [cond-mat.other] 15 Jan 2013

Transcription

1 axiv: v1 [cond-mat.othe] 15 Jan 2013 Abstact Modified Associate Fomalism without Entopy Paadox: Pat I. Model Desciption Dmity N. Saulov,a, Igo G. Vladimiov a, A. Y. Klimenko a a School of Engineeing, The Univesity of Queensland, AUSTRALIA A Modified Associate Fomalism is poposed fo themodynamic modelling of solution phases. The appoach is fee fom the entopy paadox descibed by Lück et al. Z. Metallkd pp The model is consideed in its geneal fom fo an abitay numbe of solution components and an abitay size of associates. Asymptotic behaviou of chemical activities of solution components in binay dilute solutions is also investigated. Key wods: 1. Intoduction Themodynamic modeling D, Entopy C The associate model in its vaious modifications has been successfully used fo modelling solution phases of metallugical and chemical engineeing inteest [1 7]. Besmann and Spea [3], fo example, utilised the modified associated species model fo glasses used in nuclea waste disposal. Recently, Yazhenskikh et al. [5 7] have successfully applied the associate species model to model melting behaviou of coal ashes which is an impotant poblem in coal gasification technologies. Good ageement between model pedictions and available expeimental data was epoted. Accoding to the classical associate model descibed by Pigogine and Defay [8], the stong inteactions esult in the fomation of stable configuations of mixing paticles, the socalled association complexes o, biefly, associates. Those paticles that ae not involved in the fomation of associates ae called fee paticles o, intechangeably, monopaticles. The associated solution is then consideed to be an ideal solution of monopaticles and diffeent associates. Fo example, a binay associated solution of components A and B, in which only AB-associates ae fomed, is consideed to be a tenay ideal solution of the A-monopaticles, B-monopaticles and ABassociates. Coesponding autho: School of Engineeing, The Univesity of Queensland, St. Lucia, QLD 4072, Austalia. Phone: Fax: addess: d.saulov@uq.edu.au Dmity N. Saulov. The Gibbs fee enegy G of the associated solution of n A moles of the solution component A and n B moles of B is then given by G=n A1 g A1 + n B1 g B1 + n AB g AB TS conf, 1 whee the configuational entopy of mixing S conf is expessed as S conf = Rn A1 lnx A1 + n B1 lnx B1 + n AB lnx AB. 2 Hee, n A1, n B1, n AB ae the mole numbes and g A1, g B1, g AB ae the mola Gibbs fee enegies of the A-monopaticles, B- monopaticles and AB-associates, espectively; T is the absolute tempeatue and R is the univesal gas constant. The mola factions x A1, x B1 and x AB ae defined in a usual way. Fo example, x A1 = n A1 /n A1 + n B1 + n AB. The othe mola factions ae defined similaly. The equilibium values of the mole numbes n A1, n B1 and n AB ae detemined by minimising the Gibbs fee enegy G given by Eq. 1, subject to the mass balance constaints n A = n A1 + n AB, n B = n B1 + n AB. The adjustable paamete of the model is the mola Gibbs fee enegy g AB of the eaction 3 A+B AB. 4 Besmann and Spea [3] consideed an ideal mixtue of associate species instead of monopaticles and associates. The stoichiomety of the associate species was specified so that all the species contain two non-oxygen atoms pe fomula unit. In

2 S conf [J mol 1 K 1 ] Associated solution Ideal solution x A B B Fig. 1. The configuational entopy of mixing fo the associated solution, in which only the AB-associates ae fomed, in the high tempeatue limit compaed with the ideal entopy of mixing. this appoach, contibutions of diffeent species to the configuational entopy of mixing ae equally weighted. Lück et al. [9] descibed a emakable featue of the configuational entopy of mixing given by Eq. 2. This featue was efeed to as an entopy paadox and is biefly descibed in the next section. 2. Entopy paadox Lück et al. consideed the high tempeatue limit fo the associated solution. The tempeatue T is assumed to be so high that the entopy tem plays a dominant ole in the Gibbs fee enegy of the solution and the enthalpy changes on foming diffeent associates can be neglected. The configuational entopy of mixing given by Eq. 2 in the high tempeatue limit is highe than that used in the egula solution model the entopy of ideal mixing of the solution components. Fo example, the configuational entopy of mixing of the associated solution in which only the AB-associates ae fomed is compaed with the ideal entopy of mixing in Fig. 1. At the same time, the fomation of associates models shot ange odeing in the solution. As pointed out by Lück et al. [9], it is a paadoxical esult that the configuational entopy, which is a measue of disode, appeas to be highe in a solution with odeing than in a completely disodeed solution. Anothe intepetation of this entopy paadox was descibed by Pelton et al. [10]. Conside the binay associated solution whee only AB-associates ae fomed and assume that thee ae no Gibbs fee enegy changes on foming AB-associates fom monopaticles, so that g AB = 0. In this case, the configuational entopy of mixing of the solution should be equal to that of an ideal solution, since no inteactions between mixing paticles ae assumed. Eq. 2 howeve, leads to highe values fo the configuational entopy of mixing that educes to the ideal configuational entopy of mixing only when g AB =+ ; see also Fig. 1. The oveestimation of the configuational entopy can esult in eithe undeestimation of the non-configuational entopy o oveestimation of the enthalpy of mixing o both. This, in ten, can undemine the pedictive capabilities of the model. As pointed out by Lück et al. [9] and late by Pelton et al. [10], the expession fo the configuational entopy of mixing used in the quasichemical model does educe to the ideal configuational entopy of mixing when g AB = 0. In the next section, we popose the Modified Associate Fomalism which is fee of the entopy paadox. The paadox is esolved by distinguishing between all possible spatial aangements of paticles in an associate that have not been taken into account in pevious associate models. 3. Model assumptions The model poposed in this pape is based on the following assumptions. 1 Similaly to the classical associate model [8], we assume that inteactions between mixing paticles esult in the fomation of associates which ae in a stable dynamic equilibium with each othe. The associates of the model ae undestood as a tool fo modelling shot-ange inteactions between mixing paticles. The associates of the model, howeve, may epesent eal associated complexes pesent in solution phases. 2 We assume that the associates do not inteact with each othe and ae unifomly distibuted ideally mixed ove a lattice. Equivalently, the occupancies of the sites of the associate lattice ae stochastically independent and have identical pobability distibutions. 3 In contast to the classical model, all pue solution components and the chemical solution of these components ae teated in a unified way. Moe pecisely, we assume that the solution and all its components consist of noninteacting associates of the same size, so that the associates ae composed of the same numbe of paticles. In this appoach, the contibutions of diffeent associates to the configuational entopy of mixing ae equally weighted. 4 Following the convention see, fo example, Ref. [11] fo moe details, we also assume that paticles of the same type ae indistinguishable, while paticles of diffeent types and paticle sites within an associate ae distinguishable. 4. Model desciption Fo simplicity of exposition, we exemplify the model by consideing a binay solution A B whee the associates ae composed of thee paticles. Let n A and n B be the mole numbes of A and B paticles in the solution. Since, by assumption 4 above, the paticle sites within an associate ae distinguishable, we also assume that they ae numbeed. If the 1st and 2nd sites in an associate ae both occupied by A-paticles, while the 3d one is occupied by a B-paticle, such an associate is said to be of type [AAB]. Othe types of associates ae defined similaly.

3 Thus, in the solution consideed, thee ae 2 3 = 8 diffeent types of associates,[aaa],[baa],[aba],[aab],[bba],[bab],[abb], [BBB]. It is impotant to note that all these types should be taken into account in calculating the configuational entopy. Let n [i jk] be the mole numbe of [i jk]-associates, whee the tiplet of symbolic indices i, j,k = A,B specify the associate type. Then the mass balance constaints take the fom n A = 3n [AAA] + 2n [BAA] + n [ABA] + n [AAB] + n [BBA] + n [BAB] + n [ABB], n B = 3n [BBB] + 2n [ABB] + n [BAB] + n [BBA] + n [AAB] + n [ABA] + n [BAA]. By a standad combinatoial agument, the numbe of available micostates Ω is Ω= N i, j,k=a,b i, j,k=a,b n [i jk]! 5 N n [i jk]!, 6 whee N is the Avogado numbe. Unde assumption 2 of ideal mixing of associates, the configuational entopy of mixing S conf is S conf = k B lnω, 7 whee k B is Boltzmann s constant. Applying the Stiling fomula to Eqs. 6 and 7, whee S conf = R x [i jk] = i, j,k=a,b n [i jk] lnx [i jk], 8 n [i jk] n [i i, j,k j k ] =A,B is the mola faction of the[i jk]-associates. Theefoe, the Gibbs fee enegy of the solution is given by G= i, j,k=a,b n [i jk] g[i jk] + lnx [i jk], 9 whee g [i jk] is the mola Gibbs fee enegy of [i jk]-associates. In one mole of pue solution component A thee ae 1/3 moles of[aaa]-associates. Hence, g [AAA] = 3g A, whee g A is the mola Gibbs fee enegy of the pue solution component A. Similaly, g [BBB] = 3g B, whee g B is the mola Gibbs fee enegy of the pue solution component B. The associates, which consist of both A and B paticles, will be efeed to as mixed associates. Fo the binay solution consideed, thee ae 8 2 = 6 types of mixed associates. Thei mola Gibbs fee enegies g [i jk], with [i jk] [AAA],[BBB], ae adjustable paametes of the model. These paametes can depend on the tempeatue T as g [i jk] = h [i jk] Ts [i jk], 10 whee h [i jk] and s [i jk] ae the mola enthalpy and entopy of[i jk]- associates, espectively. Altenatively, as adjustable paametes of the model, one can employ the Gibbs fee enegies of the eactions of foming mixed associates fom pue [AAA] and [BBB]-associates. Fo example, instead of g [AAB], the Gibbs fee enegy g [AAB] of the following eaction can be used, so that 2 3 [AAA]+ 1 [BBB] [AAB], 11 3 g [AAB] = g [AAB] 2 3 g [AAA] 1 3 g [BBB]. 12 The othe Gibbs fee enegies g [i jk] ae defined similaly. Note that Eqs. 5 define n [AAA] and n [BBB] as linea functions of the mole numbes of the solution components n A and n B and of the mole numbes n [i jk] of the mixed associate types with [i jk] [AAA],[BBB]. The mole numbes of the latte associates, which consist of both A and B paticles, ae the intenal vaiables of the model. These ae detemined by minimising the Gibbs fee enegy of the solution at constant n A and n B, subject to the mass balance constains of Eqs. 5. The equilibium values of n [i jk], with [i jk] [AAA],[BBB], ae found fom G = 0, 13 n [i jk] n A,n B,n [i j k ] whee the deivative in n [i jk] is calculated fo fixed n A, n B and fixed five vaiables n [i j k ], with[i j k ] [i jk],[aaa],[bbb]. By the chain ule, Eq. 13 eads G n [i jk] n [BBB] n [BBB] n [i jk] n [i jk] = Theefoe, combining the last equation with Eqs. 5 and 9 gives n [i jk] n [AAA] α [i jk] /3 n [BBB] β = exp g [AAB]. 15 [i jk] /3 Hee, α [i jk] and β [i jk] stand fo the numbes of A and B paticles in the[i jk]-associate, espectively. Fo example, α [AAB] = 2 and β [AAB] = 1. Recall that Eqs. 5 define n [AAA] as a linea function of n A and n [i jk] s, and n [BBB] as a linea function of n B and n [i jk] s, whee [i jk] [AAA],[BBB]. The chemical potential µ A of the solution component A is calculated as follows G µ A = n A n B G n [i jk] = n [i jk] n A [i jk] n [BBB] n[aaa] n A n[bbb] n A = G n A n [BBB] n [BBB] n [i jk] n [i jk] n [i jk] + n [i jk] [i jk] n A n [BBB] n [i jk] + n [i jk] [i jk] n A G + n [i jk] [AAA] n [i jk] n[i jk]. n A 16

4 Hee, the sums ae taken ove i, j,k = A,B such that [i jk] [AAA],[BBB]. Using Eqs. 9 and 14, one veifies that µ A = G n A = g A lnx [AAA]. 17 The chemical potential µ B of the solution component B is calculated in a simila way. Note that [AAB], [ABA] and [BAA] associates diffe in the spatial aangement of the constituent paticles, though they have the same chemical composition A 2 B. The same distinction holds fo [BBA], [BAB] and [ABB] associates of common composition AB 2. Let us assume, fo a moment, that only one spatial aangement of associates is allowed fo each composition, while the othe aangements ae pohibited enegetically, fo example, g [BAA] = g [ABA] =+ and g [BAB] = g [ABB] =+. In this case, the poposed fomalism educes to the associate species model [3], in which associate species ae composed of thee paticles. Unlike pevious modifications of the associate model, we take into account all possible spatial aangements of paticles in an associate. If two o moe associates of diffeent types ae spatially symmetic to each othe, then thei mola Gibbs fee enegies ae consideed equal. Fo example, 2-paticle associates [AB] and [BA] ae symmetic and theefoe ae endowed with equal Gibbs enegies g [AB] = g [BA]. In moe complex cases, howeve, Gibbs fee enegy levels ae ascibed to associates depending on the spatial aangement of paticles in them, so that both enegy splits and multiple levels may occu fo associates of common chemical composition. Fo example, if paticle sites in 3-paticle associates ae aanged linealy as shown in Fig. 2a, then g [AAB] = g [BAA] while g [ABA] can be diffeent. Altenatively, if the paticle sites ae aanged as in Fig. 2b, the thee associate types ae all symmetic to each othe, and hence, g [AAB] = g [ABA] = g [BAA]. Note that no assumptions on spatial aangements of paticle sites within an associate have been made so fa in the famewok of the poposed model. In geneal, it is impossible to detemine in advance the numbe of enegy levels and thei multiplicities fo the associates of a paticula composition. Futhemoe, such associates may have the same Gibbs fee enegy of fomation, even if they ae not spatially symmetic. As a easonable initial appoximation, all associates of a given composition can be endowed with the same Gibbs enegy of fomation. This assumption can be efined subsequently in the pocess of themodynamic model optimisation fo eal chemical systems, if use of seveal enegy levels appeas to povide a bette fit to expeimental data. Fo the est of this section and also in Sections 5 and 6, we assume, fo simplicity, that the mola Gibbs fee enegy of an associate type is completely specified by its chemical composition, g [AAB] = g [ABA] = g [BAA] g 2,1, g [BBA] = g [BAB] = g [ABB] g 1,2. 18 Hee, the subscipt 2,1 signifies that the associate consists of two A-paticles and one B-paticle, with 1,2 and simila a Paticle sites ae aanged along a hoizontal line. The associates [BAA] and [AAB] ae mio eflections of each othe about axis a; they ae not symmetic, howeve, to [ABA]. b Paticle sites ae aanged at the vetices of an equilateal tiangle. All associates ae symmetic to each othe by otation about axis b, pependicula to the plane of the tiangle. Fig. 2. Schematic views of diffeent aangements of paticle sites. indices undestood appopiately. Moe geneal case is consideed in Section 7. Using Eqs. 15 and 18, n [AAA] n 3,0, n [AAB] = n [ABA] = n [BAA] n 2,1 /3, n [BBA] = n [BAB] = n [ABB] n 1,2 /3, n [BBB] n 0,3. In tems of mola factions, Eqs. 19 ead x [AAA] x 3,0, x [AAB] = x [ABA] = x [BAA] x 2,1 /3, x [BBA] = x [BAB] = x [ABB] x 1,2 /3, x [BBB] x 0,3. Substitution of Eqs. 19 and 20 into Eq. 9 gives 19 20

5 G = n 3,0 g 3,0 + n 0,3 g 0,3 + n 2,1 g 2,1 + n 1,2 g 1,2 + n 3,0 lnx 3,0 + n 0,3 lnx 0,3 x2,1 + n 2,1 ln + n 1,2 ln 3 x1,2 3. The mass balance constaints of Eqs. 5 educe to n A = 3n 3,0 + 2n 2,1 + n 1,2, n B = 3n 0,3 + 2n 1,2 + n 2, The equilibium values of the associate mole numbes ae detemined by = 3exp g 2,1, n 2,1 n 2/3 3,0 n1/3 0,3 n 1,2 n 1/3 3,0 n2/3 0,3 5. Case g 2,1 = 0 and g 1,2 = 0 = 3exp g 1,2. 23 Recalling Eq. 18, conside the situation whee thee ae no Gibbs fee enegy changes of foming diffeent associates, that is, g 2,1 = 0 and g 1,2 = 0. In this case, the A and B paticle species mix ideally. The composition of a andomly selected tiplet of paticles follows the binomial distibution, well-known in pobability theoy; see, fo example, Ref. [12]. Moe pecisely, the pobabilities of choosing an associate of paticula compositions, o, equivalently, the mola factions of appopiate associates, ae x 3,0 = x 3 A, x 2,1 = 3x 2 A x B, x 1,2 = 3x A x 2 B, x 0,3 = x 3 B. 24 Hee, x A and x B ae the mola factions of A and B paticles. Since the total mole numbe of associates n tot is equal to n A + n B /3, then the mole numbes of associates ae given by n 3 A n 3,0 = 3n A + n B 2, n 2 A n 2,1 = n B n A + n B 2, n A n 2 B n 1,2 = n A + n B 2, n 3 B n 0,3 = 3n A + n B By diect inspection, the mole numbes given by Eqs. 25 satisfy Eqs. 22 and Eqs. 23 with g 2,1 = 0 and g 1,2 = 0. Since the Gibbs fee enegy G given by Eq. 21 is a stictly convex function of the mole numbes of associates, then Eq. 23 has no othe solution satisfying the mass balance. Futhemoe, if g 2,1 = 0 and g 1,2 = 0, then the mola Gibbs fee enegies of associates ae given by g 3,0 = 3g A, g 2,1 = 2g A + g B, g 1,2 = g A + 2g B, g 0,3 = 3g B. 26 Finally, substitution of Eqs. 25 and 26 into Eq. 21 yields G=n A g A + n B g B + n A lnx A + n B lnx B. 27 Thus, the poposed model coectly educes to the ideal solution model in the case whee thee ae no inteactions between mixing paticles. 6. Dilute solutions Pelton et al. [10] pointed out anothe inteesting featue of the associate model that occus in dilute solutions. They consideed the associated solution of monopaticles A and B and associates A 2 B. The highly odeed solution ich in component B consists pimaily of B-monopaticles and A 2 B-associates. Accoding to the authos, the chemical activity a B of the component B behaves asymptotically as1 x A /2 athe than1 x A fo small x A. Thus, lim da B/dx A = 1/2. Note, howeve, x A 0 that this behaviou of a B is obseved only in the limiting case g 2,1 =. Fo any finite value, lim da B/dx A = 1. We x A 0 pesent the poof of this esult fo the model poposed in this pape. Using the pesented technique, simila esult can be established fo the example consideed by Pelton et al. [10]. As obtained in Section 4, a B = x 1/3 0,3. The equied deivative is then calculated as follows. Fist, da B = 1 dx 0,3 dx A 3 x 2/3 0,3 dx A = 1 3 x 2/3 0,3 1 dxa. 28 dx 0,3 Secondly, assuming g 2,1 and g 1,2 finite and denoting the ight hand sides of Eqs.?? by ε 2,1 3exp g 2,1, ε 1,2 3exp g 1,2, 29 we obtain x 2,1 = ε 2,1 x 2/3 3,0 x1/3 0,3, x 1,2 = ε 1,2 x 1/3 3,0 x2/3 0,3. 30 Substitution of Eqs. 30 into the mass balance constaints of Eqs. 22, with the latte witten in tems of mola factions, gives 1=x 3,0 + x 0,3 + ε 2,1 x 2/3 3,0 x1/3 0,3 + ε 1,2x 1/3 3,0 x2/3 0,3, 31 1 x A = x 0, ε 2,1x 2/3 3,0 x1/3 0, ε 1,2x 1/3 3,0 x2/3 0,3. 32 Now, Eq. 31 implicitly defines x 0,3 as a function of x 3,0, while Eq. 32 defines x A as a function of x 3,0 and x 0,3. Diffeentiating Eq. 31 gives dx 3,0 = ε 1,2t+ 3 1ε 2,1t 2 dx 0, ε 1,2t ε 2,1t 2, 33

6 whee t x3,0 x 0,3 1/3. 34 Diffeentiation of Eq. 32 with espect to x 0,3 and substitution of Eq. 34 into the esultant expession yields dx A = 1 4 dx 0,3 9 ε 1,2t 1 9 ε 2,1t ε1,2 t 2 + ε 2,1 t ε 1,2t+ 3 1ε 2,1t ε 1,2 t ε 2,1t 1 = ε 1,2t 1 9 ε 2,1t ε 1,2 + ε 2,1 t ε 1,2t+ 1 3 ε 2,1t 2 9 t ε 1, ε 2,1t. 35 Note that x 0,3 1, x 3,0 0 and t 0 as x A 0. Theefoe, substituting Eq. 35 into Eq. 28 and taking the limit as x A 0, one veifies that lim da B/dx A = 1 fo any finite values x A 0 of g 2,1 and g 1,2. Now conside the limiting case whee g 2,1, while g 1,2 emains finite. Recalling Eqs. 29 and taking the limit in Eq. 35 give lim ε 2,1 + dx A = 2 dx 0, ε 1,2t 1 6 t3. 36 Finally, substituting Eq. 36 into Eq. 28 and taking the limit as x A 0, we obtain lim da B/dx A = 1/2 as g 2,1. x A 0 7. Model equations fo abitay numbe of solution components and abitay size of associates Conside an -component solution A 1... A and assume that the components and the solution itself all consist of m- paticle associates. Thee ae m distinguishable types of associates. Some of these types have the same compositions. Now conside an m-paticle associate that consists of paticles of types A 1,...,A, espectively. Thus, the composition of the associate is specified by the -tuple of nonnegative integes satisfying k k = m. Omitting the dependence on, the set of such tuplets, which epesent all possible compositions of m-paticle associates in the -component solution, is denoted by S m. Its cadinality, that is, the numbe N comp of diffeent compositions is computed as N comp = m+ 1! m! 1!. 37 Conside a pue solution component A i. Its m-paticle associates all have the composition σ i 0,...,0,m,0,...,0. 38 }{{}}{{} i 1 i The compositions σ 1,...,σ and coesponding associate types ae efeed to as pue. The complementay set S m \{σ 1,...,σ } of mixed compositions, containing two o moe diffeent paticle species, is witten biefly as S m. Thus, S m is constituted by those -tuples fom S m with at least two nonzeo enties. In one mole of the solution component A i, thee ae 1/m moles of m-paticle associates of composition σ i. Thus, thei mola Gibbs fee enegy g σi is g σi = mg i, 39 whee g i is the mola Gibbs fee enegy of the solution component A i. In geneal, the total numbe N k1,...,k of distinguishable types of associates which have the same composition is descibed by the multinomial coefficient N k1,...,k = m!. 40 We assume that all distinguishable types of the associates of composition ae endowed with J k1,...,k diffeent values of the mola Gibbs fee enegy. Let d [ j] be the numbe of the associate types at the jth enegy level g [ j], that is, the multiplicity of the level, so that k i! J k1,...,k d [ j] = N k1,...,k. 41 All the d [ j] associate types have the same mola Gibbs fee enegy of fomation g [ j] accoding to the eaction k 1 m A 1 m k m A m A 1 k1...a k [ j]. 42 Moe pecisely, the Gibbs fee enegy of the eaction is defined by g [ j] = g [ j] k 1 m g σ 1... k m g σ. 43 The mola Gibbs fee enegies g [ j] of mixed associates of composition S m, o altenatively, the coesponding fomation enegies g [ j] fom Eq. 43 ae adjustable paametes of the model. Now let n [ j] denote the mole numbe of the associates of composition at the jth enegy level. The mass balance constaints then ead J k1,...,k n i = k i n [ j],,...,, 44 S m whee n i is the mole numbe of A i paticles in the solution. One veifies that n i = m S m J k1,...,k n [ j] mn tot. 45 Note that the total numbe of associates n tot is independent of composition of the solution since they ae assumed to be of equal size. The Gibbs fee enegy of the solution is computed as G = S m J k1,...,k n [ j] g [ j] + J k1,...,k n [ j] S m ln j] x[ d [ j], 46

7 whee x [ j] = n[ j] n tot 47 is the mola faction of an appopiate associate type. Eqs. 44 define n σi as a linea function of the mole numbe n i of the solution component A i and of the mole numbes of the mixed associates n [ j], whee S m. The equilibium values of n [ j], which ae intenal vaiables of the model, ae detemined by minimising the Gibbs fee enegy of the solution at constant n 1,...,n, subject to the mass balance constaints of Eqs. 44. The minimum is found by setting G n1,...,n = 0 48 fo all S m and fo all [ j]=1,...,j k1,...,k. Recalling that n σi ae functions of n [ j], whee S m, and using the chain ule, G + k 1,...,k G = d [ j] n σi k i/m exp = Substitution of Eqs. 44 and 46 into Eq. 49 gives n [ j] j] k 1,...,k g[. 50 Using Eqs. 44, one veifies that n σi is a function of n i and of n [ j], whee S m. The chemical potential µ 1 of the solution component A 1 can now be calculated as follows: G µ 1 = n 1 n 2,...,n J k1,...,k G k = 1,...,k S m n 1 + G n σ i + n 1 = G n σ1 + n σ1 n 1 + G S m S m J k1,...,k j] n[ n 1 J k1,...,k G n 1 = G n σ1 n σ1 n 1 51 Hee, Eq. 49 have been used. Using Eqs. 39 and 44, one veifies that µ 1 = g m lnx σ The chemical potentials of the othe components ae calculated in a simila way. Now assume that all the associate types of composition have the same Gibbs fee enegy of fomation fo any S m, so that J k1,...,k = 1 and d 1 = m!/ k i!. This assumption, which can be subsequently efined in themodynamic model optimisation of eal chemical systems, allows substantially educe the numbe of adjustable paametes of the model. Since the Gibbs fee enegies of pue associates ae fixed, the numbe of adjustable paametes is equal to N comp, whee N comp is given by Eq. 37. Then, omitting the supescipt [ j], Eq. 46 educes to G = g x k1,...,k k i! + ln m! S m n k1,...,k 53 Similaly to Section 5, conside the case when all Gibbs fee enegies g k1,...,k of fomation of associates ae equal to zeo, so that paticles of the solution components ae mixed ideally. In this case, the composition of m andomly selected paticles follows the multinomial distibution; see Ref. [12] fo moe details. The mola faction of the associates with composition S m is expessed as x k1,...,k = m! and hence, thei mole numbe is n k1,...,k = m 1! m 1 n i x k i i k i! 54 n k i i k i!. 55 Hee, n i and x i ae the mole numbe and the mola faction of the solution component A i. One veifies that the mole numbes given by Eq. 55 descibe the solution of Eq. 50 subject to the mass balance constaints of Eqs. 44 in the case g k1,...,k = 0. The uniqueness of the solution is ensued by the stict convexity of the Gibbs fee enegy G given by Eq. 53 in the vaiables n k1,...,k. Fom Eqs. 39 and 43, the mola Gibbs fee enegy of the associates of the composition is given by g k1,...,k = k i g i, 56 whee, g i is the mola Gibbs fee enegy of the solution component A i. Substitution of Eqs into Eq. 53 and a staightfowad, though lengthy, veification shows that the poposed model coectly educes to the -component ideal solution model in this case. That is, G= n i g i + lnx i. 57 Fom this eduction, it is immediately follows that the model with associates of size can be epoduced by the models with associates of size 2, 3 and so on. Fo example, conside the model with associates of size 2. Any 2-associate is a combination of two -associate. Now assume that thee is no Gibbs fee enegy change on foming any 2-associate fom coesponding two -associates. As demonstated above, the model with 2-associate educes to the ideal mixtue of -associates.

8 Thus, the model with 2-associates is moe geneal and include the model with -associates as a paticula case. 8. Effective adjustable paametes of the model The adjustable paametes of the model elated to associates of the composition ae g [ j], whee [ j]= 1,...,J k1,...,k. One should also take into account the multiplicity of the enegy levels d [ j]. In geneal, the numbe of adjustable paametes of the model inceases exponentially with the incease in size of associates m. Howeve, the numbe of adjustable paametes can be substantially educed without loss of geneality of the model as descibed below. Conside associates with the composition. In geneal, J k1,...,k enegy levels ae possible fo these associates. The mole numbe n k1,...,k of all associates with the composition is given by J k1,...,k n k1,...,k = n [ j]. 58 The mass balance constaints of Eqs. 44 take the fom n i = k i n k1,...,k,,...,. 59 S m Summing up ove [ j] in Eq. 50, one veifies that the equilibium value of n k1,...,k is calculated as whee n k1,...,k J k1,...,k Z k1,...,k = Z k1 n σi k,...,k 60 i/m d [ j] exp j] g[. 61 The function Z k1,...,k defined by Eq. 61 plays a ole of the patition function which descibes the distibution of associates of the composition ove enegy levels. That is, n [ j] n k1,...,k = x [ j] x k1,...,k = d [ j] Z k1,...,k exp j] g[ 62 In fact, Z k1,...,k is a single effective adjustable paamete elated to associates of the composition. The othe themodynamic paametes of the model can be expessed in tems of Z k1,...,k, whee S m, and thei deivatives. Indeed, using Eqs. 39, 43, 59 and 62, one veifies that the Gibbs fee enegy of the solution given by Eq.46 takes the fom xk1,...,k G= n i g i + n k1,...,k ln. 63 Z S k1 m,...,k Accoding to Eq. 61, Z k1,...,k vaies fom zeo to infinity. When the associates of the composition ae pohibited enegetically, that is g [ j] + fo all [ j], Z k1,...,k appoaches zeo. If the associates of the composition is highly pefeable, Z k1,...,k appoaches infinity. Simila to the Gibbs fee enegies of associates, the optimal values of the effective adjustable paametes Z k1,...,k can be detemined by the tial-and-eo pocedue which is conventionally used in themodynamic model optimisation of eal chemical systems. 9. Excess Gibbs enegy tems Similaly to the modified associate species model [4], egula o, in geneal, polynomial, excess Gibbs fee enegy tems can be included into the poposed fomalism. These tems take into account inteactions between associates. A pobabilistic intepetation of the polynomial excess Gibbs fee enegy tems is pesented in Ref. [13]. In fact, this intepetation povides a theoetical justification fo such tems. Teating the associates as paticles, the esults of Ref. [13] ae applicable to the model pesented in this pape. Note that it is desiable to use the excess tems only fo fine tuning of the model, while the main adjustable paametes ae the mola Gibbs fee enegies of associates. The following two conditions on the absolute value of the inteaction paametes should be satisfied; see, fo example, the monogaph by Pigogine and Defay [8] fo moe details. i The absolute values of the inteaction paametes should be small compaed with the mola Gibbs fee enegies of the associates. As pointed out by Pigogine and Defay [8], if the inteaction between associates, say C1 and C2, is sufficiently stong to alte the vibational and otational states of the associates, then the associate C1C2 is included into the set of associates by the definition of the classical associate model. In the famewok of the poposed fomalism, one should conside the associates of lage size. ii The values should also be small in compaison with. Othewise, the assumption of ideal mixing of associates is less justified. Again, lage associates should be taken into account. Note, howeve, that the second condition is sometimes elaxed in ode to fit the expeimental data available fo eal solutions. This is the case, fo example, fo the solutions with immiscibility, which is the esult of elatively weak, compaed with the Gibbs fee enegies of associates, epulsive inteactions between the associates. 10. Discussion on applicability of the model The suggested modified associate fomalism belongs to associate-type models. As a esult, the ange of applicability of the fomalism is at least the same as that of the classical associate model o the associate species model. Thee ae, howeve, some distinctions in using the suggested fomalism and the pevious modifications of the associate model. These distinctions ae discussed below. Since the pesent study is intended as a theoetical intoduction to the modified associate fomalism, the esults on its application to eal chemical systems will be epoted elsewhee.

9 In contast to the pevious modifications, whee associates of abitay compositions can be included into the set of associates, the compositions of associates ae defined by thei size in the famewok of the modified associate fomalism. Theefoe, a modelle should pay special attention to selection of the size of associates. If expeimental infomation about compositions of associates that pesent in the solution phase is available, this infomation should clealy be taken into account. The size of associates should be lage enough to incopoate those compositions. When the suggested fomalism is applied fo themodynamic desciption of multicomponent solution phases, the size of associate should be lage enough to incopoate the associates of the compositions, which coincide with those of maximum odeing in all binding binay systems. Fo example, conside the tenay system A B C and assume that in the binay system B C the composition of maximum odeing is that of the associate B 2 C, while in the systems A B and A C maximum odeing occus at the compositions of the associates AB and AC, espectively. To incopoate the equied composition, the size of associates should be divisible by 2 and by 3. Theefoe, the size of associates fo the tenay system A B C should be at least 6. It is desiable to use 6-paticle associates fo themodynamic model optimisation fo all the binay systems. Thee is, howeve, no need to use 6-paticle associates fom the beginning. Let us assume that the systems A-B and A-C ae initially optimised with 2-paticle associate, while the system B-C is optimised with 3-paticle associate. As discussed in the end of Section 7, the model with associates of size can be epoduced by the models with associates of size 2, 3 and so on. Using this popety of the poposed fomalism, the desciptions of the binay systems with 2- and 3-paticle associates can be eplaced by the equivalent desciptions with 6-paticle associates in a staightfowad way. Then, the binay systems can be combined into the tenay one. As a demonstational example, conside the system A B. When the system is optimised with 2-paticle associates the Gibbs fee enegies of the associates g [AA], g [BB], g [AB] and g [BA] ae known. Due to spatial symmety, g [AB] = g [BA]. In the equivalent desciption of the system A B, 6-paticle associates ae fomed fom thee 2-paticle associates with no Gibbs fee enegy changes on such fomations. The associate of the composition A 6 is fomed fom thee [AA]-associates with g A6 = 3g [AA]. The associate of the composition A 5 B is fomed ethe fom two [AA]-associates and[ab]-associates o fom two[aa]-associates and [BA]-associates. Then, ecalling that g [AB] = g [BA], A 5 B- associates have only one enegy level g A5 B= 2g [AA] +g [AB] with multiplicity 6. The A 4 B 2 -associate can be fomed fom two [AA]-associates and one [BB]-associate with the enegy level g [1] A 4 B 2 = 2g [AA] + g [BB] and the multiplicity d [1] A 4 B 2 = 3. Altenatively, A 4 B 2 -associate can be fom fom one[aa]-associate and eithe two [AB]-associates o two [BA]-associates o one [AB]- associate and one [BA]-associate. In this case, the enegy level is g [2] A 4 B 2 = g [AA] + 2g [AB], and its multiplicity is d [2] A 4 B 2 = 12. In a simila way, one veifies that A 3 B 3 -associates have two enegy levels: g [1] A 3 B 3 = g [AA] +g [BB] +g [AB] with the multiplicity d [1] A 3 B 3 = 12 and g [2] A 3 B 3 = 3g [AB] with the multiplicity d [2] A 3 B 3 = 8. A 2 B 4 - associate also have two enegy levels: g [1] A 2 B 4 = g [AA] + 2g [BB] with the multiplicity d [1] A 2 B 4 = 3 and g [2] A 2 B 4 = g [BB] + 2g [AB] with the multiplicity d [2] A 2 B 4 = 12, while AB 5 -associates have only one enegy level g AB5 = 2g [BB] + g [AB] with the multiplicity equals to 6. B 6 -associates also have only one enegy level g B6 = 3g [BB] with the multiplicity equals to 1. The binay systems A C and B C ae teated similaly. Pevious modifications of the associate model allow moe flexibility in fitting expeimental data compaed with the suggested appoach, since stoichiomety of associates can be abitaily alteed to impove the fit. This additional flexibility can be helpful fo modelling a binay system with high degee of odeing, whee the descibed paadox is of elatively low pactical impotance. Howeve, when such a binay system is combined with othe binaies to fom multicomponent model, abitay selection of associate stoichiomety can esult in the negative consequences of the entopy paadox. Simila to the pevious modifications, a miscibility gap can not be epoduced without excess Gibbs fee enegy tems in the poposed fomalism. Howeve, as demonstated by Besmann et el. [4] fo the associate species model, immiscibility, which is the esult of epulsive inteactions between associate o associate species, can be accuately epesented by the associatetype models with the polynomial excess Gibbs fee enegy. Teating associates o associate species as paticles and using the esults of Ref. [13], one veifies that the coefficients of the polynomial excess Gibbs fee enegy explicitly elate to enegies of inteactions between associates. It is impotant to note that, in the case of the modified quasichemical model [10], immiscibility is descibed by empiical polynomial expansions of the Gibbs fee enegy of the quasichemical eaction. In contast to associate-type models, physical meaning of the coefficients of such expansions is not clea. Anothe distinction of the poposed fomalism fom the pevious modifications is elated to multicomponent associates, that is, the associates which consist of paticles of thee o moe types. In the pevious modifications, no multicomponent associates ae initially consideed. This implies that all multicomponent associates ae assumed to have infinite positive Gibbs fee enegy. Multicomponent associates ae usually included into the model, when such associates ae equied to fit available expeimental data in multicomponent systems. In contast, the compositions of multicomponent associates ae pescibed by the selected size of associate in the poposed fomalism. It would be beneficial to develop a method fo estimating the Gibbs fee enegies of the multicomponent associates fom those of binay associates. The development of such a method, howeve, is the topic fo sepaate study. As an initial appoximation, the Gibbs fee enegies of fomation of the multicomponent associates fom pue associates can be set to zeo. If expeimental data ae available in multicomponent systems, the Gibbs fee enegies of the multicomponent associates can be adjusted in a conventional tial-and-eo pocedue.

10 A possibility to select compositions of multicomponent associates abitaily, which exists in the pevious modifications, could povide moe flexibility in data fitting compaed with the fixed set of compositions. It could be the case, that abitaily selected additional multicomponent associates povide bette fitting to paticula expeimental data. One should ecognise, howeve, that the pice fo the possibility to select the compositions of associates abitaily is the entopy paadox, which undemines the fundamental pedictive capabilities of the model. In ou opinion, the necessity of additional multicomponent associates fo fitting the available expeimental data indicates that the size of associates should be inceased to incopoate the equied compositions. 11. Conclusion In this pape, the Modified Associate Fomalism has been poposed fo themodynamic modelling of solutions. The pesented appoach is fee fom the entopy paadox and coectly educes to the ideal solution model, whee thee ae no Gibbs fee enegy changes on foming associates of diffeent types. Asymptotic behaviou of chemical activities of solution components in binay dilute solutions has been investigated. It has been demonstated that the deivative of the chemical activity of a solution component in its mola faction at teminal composition has the expected value fo any finite value of the Gibbs fee enegies of fomation of associates. Acknowledgment [6] E. Yazhenskikh, K. Hack, M. Mulle, Citical themodynamic evaluation of oxide systems elevant to fuel ashes and slags. Pat 2: Alkali oxide - alumina systems, Calphad - Compute Coupling of Phase Diagams and Themochemisty [7] E. Yazhenskikh, K. Hack, M. Mulle, Citical themodynamic evaluation of oxide systems elevant to fuel ashes and slags. Pat 3: Silicaalumina system, Calphad - Compute Coupling of Phase Diagams and Themochemisty In pess. [8] I. Pigogine, R. Defay, Chemical Themodynamics, Longmans, Geen and Co., [9] R. Lück, U. Geling, B. Pedel, An entopy paadox of the association model, Zeitschift fu Metallkunde [10] A. D. Pelton, S. A. Degteov, G. Eiksson, C. Robelin, Y. Dessueault, The modified quasichemical model I - Binay solution, Metallugical and Mateials Tansactions B: Pocess Metallugy and Mateials Pocessing Science 31B [11] L. D. Landau, E. M. Lifshitz, Satistical Physics Pat 1, Vol. 5 of Couse of Theoetical Physics, Pegamon Pess, [12] W. Felle, An Intoduction to Pobability Theoy and Its Applications, Vol. 1, John Wiley & Sons, [13] D. Saulov, On the multicomponent polynomial solution models, Calphad - Compute Coupling of Phase Diagams and Themochemisty The pesent wok has been suppoted by the Austalian Reseach Council. Refeences [1] F. Somme, Association model fo the desciption of the themodynamic fuctions of liquid alloys II Numeical teatment and esults, Zeitschift fu Metallkunde [2] J. W. Hastie, D. W. Bonnell, A pedictive phase equilibium model fo muticomponent oxide mixtues Pat II. Oxides of Na-K-Ca-Mg-Al-Si, High Tempeatue Science [3] T. M. Besmann, K. E. Spea, Themochemical modeling of oxide glasses, Jounal of the Ameican Ceamic Society [4] T. M. Besmann, N. S. Kulkani, K. E. Spea, Themodynamical analysis and modeling of the Al 2 O 3 -C 2 O 3, C 2 O 3 -SiO 2 and Al 2 O 3 -C 2 O 3 -SiO 2 systems elevant to efactoies, Jounal of the Ameican Ceamic Society [5] E. Yazhenskikh, K. Hack, M. Mulle, Citical themodynamic evaluation of oxide systems elevant to fuel ashes and slags. Pat 1: Alkali oxide - silica systems, Calphad - Compute Coupling of Phase Diagams and Themochemisty