Part XII – December 1969 – Papers - Series Representation of Thermodynamic Functions of Binary Solutions

Analytical representation of the thermodynamics of solutions is highly desirable from the standpoint of accuracy, compactness, and numerical manipulations. In particular, computer calculations are greatly implemented. Mathematical considerations show that previous expressions have one or more serious defects. This investigation shows a Fourier series to be satisfactory but that it is also possible to derive a new series which fits certain additional conditions. Included examples show the value of analytical expressions in giving a simple characterization of each system using some two to five parameters, the elimination of the Gibbs-Duhem integration, and the es timation of the error for the experimental function as well as derived functions. It is further shown that the present characterization provides easy comparison between systems. IN the past, thermodynamic calculations have depended to a considerable extent on tabular and graphical methods. As the volume and precision of such data increase such methods become less satisfactory. Specifically, the selection of the optimum representation and the estimation of errors require statistical methods which in turn require analytical representation. The utilization of such data require further manipulations which are best done analytically for maximum precision. For example, phase equilibria are determined by common tangents to free-energy curves: a graphical determination is normally of low accuracy. As computers are increasingly used analytical representations become almost mandatory. Insufficient mathematical consideration has been given previously to the selection of empirical expressions. Those expressions having some theoretical justification are generally too inflexible and mathematically unattractive. We consider the problem in some detail and show that a Fourier series can be effectively used. Also a new series is defined which has certain advantages. ANALYSIS We wish to consider the analytical representation of the heat of mixing, AH, the excess free energy, ?Gxs, and the excess entropy, ?sXS, as a function of composition, X, for binary solutions relative to the pure components in the same state. When a distinction is not required, we use W to denote any one of the above functions. One may use a Taylor expansion around X = 0 to generate a power series. As the derivatives are un- known we represent the series as W = A + BX + CX2 + DX3 + EX4 + ... [l] where the constants A , B, C , ..- are to be selected to provide some optimum fit. For the extremes of composition W is necessarily zero so it follows that A = 0 [2a] B +C + D + E +••• = 0 [2b] Nonelectrolytes, which we are considering, appear to satisfy the condition that d3W/dx3 = 0 [3] in the terminal regions. This is the basis of the a, ß, and Q functions used by Hultgren et al.' and others. While this condition does not have a strong theoretical basis it does appear desirable that any analytical relation should satisfy this condition. Darken2 and Turk-dogan and Darken3 have shown that many systems exhibit this behavior over an extended range from each terminal region, departure being restricted to a limited intermediate region. Since we have no a priori knowledge as to where this transition occurs we can require that this condition be satisfied only as a limit at the extreme compositions as a general condition. We will show later how more restricted conditions can be included in specific solutions. Darken2 has called this behavior the quadratic formalism; we call our application the limiting quadratic formalism, LQF. This condition applied to the above power series requires that D = 0 [4a] 4-3-2E +5-4-3_F + 6 • 5 . 4G + ••• =0 [4b] The form of the power series normally used, due to Margules,4 is W=X(1-X)(A + BX + CX2 + DX3 + EX4 + •••) [5] where A, B, C, --. are a new set of constants. (Guggenheim5 has given a variation of this expression in a more desirable form. Since, however, it is contained in the above expression it does not require separate consideration.) This form is precisely what results by incorporating the conditions in Eq. [2] into the power series and regrouping the constants. The LQF requires that B =C [6a] and 4.3.2(D-C) +5-4-3(E-D) + ••• =0 [6b] Thus, the correct form of the Margules expression with two adjustable parameters is w =X(1-X)[A + B +X2-2/3x3)] 171 and the EX4 term must be included before three adjustable parameters are permitted.