MODIFIED HELGESON, KIRKHAM AND FLOWERS EQUATION OF STATE FOR STANDARD PARTIAL MOLAL PROPERTIES OF IONS
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.