A new equation for the temperature dependence of the excess Gibbs energy of solution phases

During CALPHAD optimization procedures, the excess Gibbs energy of solution phases is usually described by the Redlich–Kister (RK) polynomial. The RK polynomial includes interaction parameters (Li), usually described as linear functions of temperature (Li=ai−biT), with semi-empirical parameters ai and bi usually having the same signs. Consequently, the excess Gibbs energy will change its sign at a certain high temperature and will have an infinite value at infinitely high temperature. This is in conflict with thermodynamic constraints and can also lead to artificial phase stabilization in some of the calculated phase diagrams, such as the appearance of an artificial miscibility gap having no maximum critical temperature. Such artifacts can usually be avoided, if the following new semi-empirical equation is used for the interaction energy of solution phases: image. Parameter h0i (J/mol) is the enthalpy part of Li at T=0 K, while parameter τ0i (K) is the temperature at which Li would change its sign if it were extrapolated linearly from T=0 K. Parameter h0i can have any sign, but parameter τ0i must be always positive. The new equation can be extended to describe more complex behavior in certain systems as: image, with aij being adjustable semi-empirical parameters for the given Li.