A method for predicting the variation with temperature of the solubilities of gases in non-polar liquids

Abstract
An expression for the chemical potential due to Kirkwood & Boggs is adapted to give rigorous expressions for Henry’s coefficient (H) for the solubility of a gas in a liquid and for the temperature dependence of this coefficient, in terms of radial distribution functions (g) and a molecular coupling parameter. If the solute-solvent and solvent-solvent molecular interactions are similar in strength the expression forTdlnH/dTreduces toTdlnH/ dT=L/RT+ (1+αT) InP° /H(i) whereL, aandp° are the molar latent heat, the coefficient of thermal expansion and the vapour pressure of the pure solvent. Equation (i) is closely obeyed by the simple systems Ar-CH4, Ar-O2and Ar—N2, though it becomes markedly less accurate when applied to the solubilities of common gases in liquids. This is to be expected since the solute-solvent and solvent—solvent intermolecular force fields are then very different. By assuming these force fields to be of the Lennard-Jones type and making simplifying assumptions relatinggfor the solute in the solvent togfor the pure solvent, the equationTdlnH/ dT=L/RT+ (1 +αT) InP° /H-Q0(1 -ε°αβσ3αβ/ε°ββσ3ββ(ii) is then obtained in whichQ0=L/RT- 1 +αT(1 +αT) InP°Vβ/RT, whereVβis the molar volume of the solvent,ε°ββ,σββ,ε°αβandσαβare the Lennard-Jones force constants for the solvent-solvent and solute—solvent interactions respectively. This equation is found to predictTdlnH/ dTfor gases dissolved in common liquids with sufficient accuracy to be of practical value The equationTdlnH/ dT= 2-αT+ (1 +αT) InRT/VβH, valid at solvent reduced temperatures between about 0.5 and 0.65, is found in practice to provide a useful approximation to (ii) both for simple systems and for the permanent gases dissolved in common solvents. Expression (i) is shown to be related to an expression previously developed by Longuet-Higgins.