2009 Portland GSA Annual Meeting (18-21 October 2009)

Paper No. 13
Presentation Time: 11:25 AM

MODIFIED HELGESON, KIRKHAM AND FLOWERS EQUATION OF STATE FOR STANDARD PARTIAL MOLAL PROPERTIES OF IONS


WALTHER, John V., Huffington Department of Earth Sciences, Southern Methodist University, Dallas, TX 75225, walther@smu.edu

An equation of state for standard thermodynamic properties of aqueous species at P and T above the critical point of H2O was developed by Helgeson and Kirkham (1976) and Helgeson, Kirkham and Flowers (1981) and later revised (revised HKF). The electrostatics in the revised HKF equations are based on Born charging from a vacuum to a dielectric medium. This produces a linear dependence of standard Gibbs energy on the inverse of the ion’s radius, 1/ri. In order to fit the data available ri of cations and anions are a function of P and T. Also ri for cations and anions of the same crystallographic size are different which contravenes Born charging theory. Because of the quadrupole properties of H2O it is proposed that in addition to a Born charging term the Gibbs energy of ion-water-dipole and ion-water-quadrupole interactions Goion-q be explicitly considered:

Goion-q = - niNo|Zi|eoμH2O/(ri + rH2O)2 + niNoZieoqH2O/2(ri + rH2O)3

In the equation ni is the number of H2O molecules in the primary hydration sphere of ion i, No and eo are Avogadro's number and the electron charge, Zi gives the ion's charge, μH2O stands for the dipole moment of a water molecules, rH2Odenotes the radius of a water molecule and qH2O is its quadrupole moment. As the equation indicates negatively charged ions will have a greater stability than positive ions of the same charge and size as experimental heat of hydration experiments indicate. The electrostatic contribution to standard Gibbs energy is then

Goi = (No(Zieo)2/2(ri + 2rH2O)(1/ε-1) - niNo|Zi|eoμH2O/(ri + rH2O)2 + niNoZieoqH2O/2(ri + rH2O)3

where ε is the dielectric constant of water. A term for induced dipoles which vary as 1/(ri+rH2O)4 could also be added which gives a "fit" equation in terms of 1/(ri+rH2O)nj , nj = (1, 2...). Provisions for transition metal ions with unsymmetrical 3d orbitals might also be needed as they have larger heats of solvation than s-orbital ions of similar charge and size.