Redlich-Kister Expansion

This equation can be used to fit activities in two component mixtures over the entire concentration range, and to calculate the activities of both components. It is used E-AIM for aqueous solutions of single organic compounds, and is an alternative to UNIFAC for those compounds for which the necessary activity data exist.

The expansion is described by Prausnitz et al. (1986), and is also covered in other chemical engineering textbooks. McGlashan (1963) gives examples of the use of the equation to illustrate the types of activity curves that are observed for different two component mixtures.

In E-AIM the Redlich-Kister expansion can be used with up to 10 fitted parameters. Examples of the use of the equation, to represent the properties of aqueous solutions of dicarboxylic acids at 298.15 K, are given by Clegg and Seinfeld (2006). These acids are present in the public database of compounds, on the Available Compounds selection page, and can be included in model calculations.

Expressions for solute (S) and water (W) mole fraction activity coefficients are given below, for terms including the first 5 fitted parameters (C1 to C5).

ln(fS) = ge/RT + [(1 - xS) × d(ge/RT)/d(xS)]

ln(fW) = ge/RT - [xS × d(ge/RT)/d(xS)]

where the excess Gibbs energy per mol of substance, ge, is given by:

ge/RT = xS(1 - xS)[C1  +  C2(1 - 2xS)  +  C3(1 - 2xS)2  +  C4(1 - 2xS)3  +  C5(1 - 2xS)4 ...]

and its differential with respect to xS, d(ge/RT)/d(xS), is given by:

d(ge/RT)/d(xS) = C1(1 - 2xS)   +  C2[-2xS(1 - xS) + (1 - 2xS)2]

+  C3[-4xS(1 - xS) + (1 - 2xS)2](1 - 2xS) +  C4[-6xS(1 - xS) + (1 - 2xS)2](1 - 2xS)2

+  C5[-8xS(1 - xS) + (1 - 2xS)2](1 - 2xS)3 ...

where Ci are the fitted parameters. Both activity coefficients are relative to a reference state of the pure liquid. The logarithm of the activity coefficient of the solute can be adjusted to a reference state of infinite dilution in water (fS*) by subtracting the value obtained with the equation above for xS = 0, yielding:

ln(fS*) = ln(fS) - (ge/RT + d(ge/RT)/d(xS))

This expression is valid for any concentration but, as stated above, the final term in ge/RT and its differential is calculated for xS = 0.


S. L. Clegg and J. H. Seinfeld (2006) Thermodynamic models of aqueous solutions containing inorganic electrolytes and dicarboxylic acids at 298.15 K. 1. The acids as nondissociating components. J. Phys. Chem. A 110, 5692-5717.

M. L. McGlashan (1963) Deviations from Raoult's law. J. Chem. Educ. 40, 516-518.

J. M Prausnitz, R. N. Lichtenthaler, and E. Gomes de Azevedo (1986) Molecular Thermodynamics of Fluid Phase Equilibria, 2nd. Edn., Prentiss-Hall.