电解质溶液理论模型及活度计算方法（学校XX XX省 XX 邮编XX）摘要:近年来，电解质溶液活度的研究日益活跃，本文回顾了电解质溶液活度理论研 究的发展进程，列举了现目前电解质溶液一些理论模型及活度预测及计算的方法。关键词:电解质溶液 热力学模型 活度系数 计算方法Electrolyte Solution Theory Model And Calculation Of MethodLiu Bangchao(Panzhihua University, Sichuan, Panzhihua, 617000)Abstract： In recent years, study of electrolyte solution activity has become increasingly active ,this paper reviews the theories of electrolyte solution and the process of development,list the current theoretical models of electrolyte solution and a number of activity predictions and calculation of method .Keywords: Electrolyte solution Thermodynamic model Activity coefficient Calculation弓丨言:玄冶金过程中，实际的体系绝大部分为非理想体系，因此在热力学性质的计算中必 须考虑以活度代替浓度，以便对体系的热力学行为进行准确的分析，因此建立一个活度 数据库就显得十分重要。活度的概念是1908年由Lewis提出的，早期的研究工作主要集中在实验测定方面。 由于体系的数量众多，依靠实验已远远不能满足冶金、化工、材料等工艺过程的要求。 唯一的解决办法就是在有限的实验数据的基础上进行理论上的估算，即热力学模型的方 法。本文在回顾许多资料的基础之上，总结了近几年以来的电解质溶液理论模型及活度 计算方法。一、电解质溶液活度计算方法的回顾11、德拜一休格尔理论1923年,荷兰化学家德拜（Debye）和他的助手休格尔（H iicke1 ）提出了微观电解质溶液 理论，既当电解质酸、碱或盐溶于水中时，分子就离解为带正电荷的阳离子和带负电荷 的阴离子。他们提出的物理模型是：一个阳离子（中心离子）最近的周围有较多的阴离子 形成一种阴离子氛。如果考虑离子的直径，公式ln Yiz 2 De 2kiD2 DkT 1 + kaK 2 =5 D工zDkT 107C DZ2i i其中：f为离子的活度系数,Z为离子的电价，三为质子电荷；D为介质的介电常数,aii为正、负离子的有效半径之和,k为玻尔兹曼常数,N为阿弗加德罗常数。如果把离子看成点电荷：公式为：（极限公式）lg 十-AZ |JT(B：JT I)条件：适用于任何溶剂中的极稀的强电解质溶液。lgV = lg7 lg(A 0. OOVmM)对于实用的活度系数：+1电解质溶液的活度系数是随着离子强度的增大而斌小的，且离子的价数越高及溶剂 的介电常数越小，则离子平均活度系数的减小就越显著。12、戴维斯公式(1 0.20 I)1 p适用于25C电解质水溶液活度系数的计算。1961年，戴维斯又修改了上式为：二 0.50 Z Z ( 0.30 I)lg 丫=0.50 |Z Z+ +lg 丫+不但修改德拜一休格尔公式，还有许多化学家设法改善这个理论，使它能够用到高 浓度的电解质溶液中去。如1970年Scatched建立了Satched理论；1973年Pitzer建立了Pitzer 理论。1.3、以Ptzer电解质溶液理论为基础的二个改进型方程131、Pitzer从电解质水溶液的径向分布函数出发，提出了溶液的总过量自由能表达式， 再导出了渗透系数与活度系数的计算公式。并在1973年提出了计算电解质溶液渗透系数 和溶液活度系数的方程皿2。渗透系数的表达式如下：cA011/22V V 小 口仆(2V V )3/20 1 = Z Z + mm_丄(p (0) + p (1) e-ai1/2) + m 2m_xC 0M X 1 + bi 1/2VMX MXVMX活度系数的表达式如下：对1：1价型电解质有：lnJMX = fr+mPMX + m2CMX，上式计算出的活度 系数是以无限稀释为参考态的。132、PitzcrLi方程.对于高浓度电解质溶液，固其质量摩尔浓度m可能很大，特别是时于纯熔盐，更有 m T 0，因此前面所说的Pitzer方程误差较大，为此Pitzcr和李以圭提出了以摩尔分数为 浓度单位，从过量Gibbs自由能出发得到的溶剂和溶质的活度系数计算公式：d V G ex / RT 丿 2 A 13/ 2宀、t ln Y = = x x + 3 X 2对溶剂：1d n1 + p 11/22X%G/rT =a 2n(l+P1f) + (叟-2仔)对溶质：lnn2二叫X p 1 +p富-213/2 )21 + pI 1/21X此方程已成功地用于550C及108Pa时的计算。可以看到，计算的温度、压力范围 已有较大幅度扩大，明显优于经典的计算公式。当把这个方程用于混合电解质水溶液的 活度系数的计算时，其符合度在可接受范围。14、Clegg-Pitzer 方程Clegg和Pitzer提出了一个计算电解质活度系数的方程，其主要构想为：把体系 中所有的组分都作为作用粒子，并用摩尔分数X表示浓度。对单一对称性电解质水溶 液，Clegg-Pitzer方程的表示式如下：2A I3/2lny =xx -XXB exp(oIi/2)+X2 w+(X -X)V+XX2 (2-3X )V1 1+pIl/2M X MXXMX 1MXMX 1 11 MXMX 1MXX2I1/2(1-2I ) Xlny y =-2A _ln(+p/1/2)+JZ)+ mxB g(aI1/2)+(1-X )exp/1/2)m 2X pX1+pI1/22 MXXMXXX+2X2W + 2X V + X (2-3X )V -2W1 1MXMX 1MXMXMX 1MX1MXClegg-Pitzer方程已被成功地用于计算0 % -100%HNO3水溶液体系的热力学性质计 算以及HNO3-H O体系的活度系数的计算7。以上应用结果表明，Clegg-Pitzer方程的适 用范围及计算精度均优于原Pitzer方程，不足之处是由于其是一个四参数的方程，使得 计算过程较为复杂。二、热力学模型在电解质溶液理论中的研究现状 21、高阶亚正规模型在前人提出的一些算法的基础上9，王之昌等10提出了高阶亚正规模型的一种 新的解析算法，将已知的组元活度表达为组元的多项式函数，以子二元系或多元系的边 界条件为基础，计算其它组元全浓度范围内组元的活度系数。具体的解法如下：对一个由N组元组成的体系,可用N-1个独立的变量y来表示，iy与各组元摩尔分i数关系为：ix - niy (ii-y)ixjiy, x=y (1-y), x=yy对ij k三元系则有：ij jik/Y = 1x , Zxk取浓度变换：ix+ xjk如果已知三元系中组元i的活度，根据亚正规溶液模型的数学解析法，ln y)可以 表示为多项式函数：In 丫iA Yjk(1)2Ajkl ny= 1j2lnykj-12 0 2+j - 1dlLAk +工呂1 j - k Y (1- j)1 j -k+7(1-T)A YjZkjk0YjZkjk0.(3)(2)式中Y、Z已知，A 由(1)式中In y拟合定值。将拟合出的A 代入式(2)、 JkiJk(3)中就可求出其它两组分的活度系数。此方程既利用了二元系边界条件又利用了三元系的边界条件，使得它的三元、二元 和纯组元的数据热力学严格一致，适用于实验点无规律分布或某些边界条件未知的体 系，很多学者运用此模型1112 13进行三元及多元系活度的计算，均取得了很好的效果。 此方法应用时必须已知其中一组元的热力学性质，这对热力学数据缺乏的体系来说，应 用将受到很大的限制。预测功能弱，需要二元及三元系甚至采用多元实测数据才能预测 多元系的热力学数据。22、积分方程理论积分方程理论是以Ornstein-Zemike积分方程为基础，从求解体系中各种粒子的径 向分布函数g(r)着手，建立统计力学理论模型14。OZ方程将直接相关函数c(r)和总相 关函数h(r)=g(r)-1用积分方程表示如下：h (r ) = c (r ) + J p (r ) c (r ) h (r ) dr “、1212313233，(1)式中，p为分子数密度。因为函数h(r)和c(r)均属未知，所以求解上述积分方程必 须引入近似方法。工程上常用的近似方法为PercusYevick(PY)近似和平均球近似法(MSA)。采用PY近似，己求解出具有下列位能函数(r)的硬球流体:(r d)(2)的解析解。phs 1 + g+ g 2 g 3 pkT =(1 -g )3Ahs AidNkTghs (d)=4g- 3g 271-g)(1 -g /2) (1 -g)3 9.(3)采用平均球近似：f h(r) = 1(r d)式中，d为硬球直径，p、A分别为压力和亥氏自由能，g为堆积因子(g二p血3/6)。 式(3)为纯硬球流体的状态方程(即著名的Camahan Starting方程)，式(4)为纯硬球流体 的亥氏自由能公式，式(5)为在硬球直径d处的径向分布函数。由于硬球流体现已普遍 用做实际流体的参考基准，故以上公式在工程上己获得广泛应用。对于具有下列位能函数的电解质离子水溶液:(r d)也已获得解析解15:A Ahs式中,u (r)= d)p z2兀 Qp 2(7) 2NkTpDkTi r为屏蔽系数，需从下列六次方程中解出： 4r 2 =a 2 Z p zi_(兀 /学比2 Pn 2,i，1 r 31+rd +古+环i1+rdii式中，d为离子i的水化直径，为电子电荷，D为介电常数。上式己用于任何价态 的单一和混合电解质水溶液热力学计算，其预测精度高于Pitzer方程。采用平均球近似，对于具有下列位能函数的离子偶极混合物：ua (r) = z z e2 / r,i jucd(r) = z |LX e(r p )/r2,i d