韩正

（北京化工大学，北京100029）

研究与开发

非水电解质溶液的状态方程

韩正

（北京化工大学，北京100029）

以甲醇为溶剂的非水电解质溶液作为研究对象，构造出了能够准确描述非水电解质溶液热力学性质的状态方程。状态方程中包含2个可调节参数，在不改变可调节参数的情况下，该状态方程对于不同温度、浓度的非水电解质溶液的热力学性质保持着良好的预测能力；其中的亥姆霍斯自由能由微扰理论计算得到，并阐述了非水电解质溶液的蒸汽压、溶剂活度的求解过程。

非水电解质溶液；状态方程；热力学

根据电解质溶液中溶剂的不同，人们把以水作为溶剂的电解质溶液称作“水电解质溶液”，而把以其它物质（如甲醇、乙醇等）作为溶剂的电解质溶液称作“非水电解质溶液”。用来描述电解质溶液热力学性质的模型主要有2种，他们分别是活度系数模型和状态方程模型，而本次研究所开发的模型正是针对以甲醇为溶剂的非水电解质溶液所开发的状态方程模型。

状态方程（Equation of State，EOS）是表征流体压强、密度和温度3个热力学参量的函数关系式，在确定了流体的状态方程之后，就可以确定流体在不同温度下的压力和密度。电解质溶液的状态方程模型法是一种首先利用热力学理论确定电解质溶液体系的亥姆霍兹自由能，然后再将该亥姆霍兹自由能与电解质溶液的其他热力学性质进行关联——如活度系数、渗透系数、蒸汽压和密度等，从而在数值上对电解质溶液热力学性质进行准确描述的模型化处理方法。

电解质溶液在化学工业生产过程之中十分常见，深入研究电解质溶液的热力学性质对于化工生产来说具有十分重要的意义。对于电解质溶液的蒸汽压预测工作在许多领域中扮演着非常重要的角色，这些领域包括热泵吸附、精馏、有机溶剂循环利用和盐湖中贵金属的萃取等等。在过去的几十年里，人们对于电解质溶液热力学性质的测量和模型化处理已经在许多方面取得了很多阶段性成果，但是其中绝大多数的研究成果都是局限于水电解质溶液范围之内。到目前为止，人们对于非水电解质溶液的研究成果还是屈指可数。尽管对于水电解质溶液来说已经有很多实验数据存在于文献之中以供查阅，但是与之相较而言，人们可以查阅到的关于非水电解质溶液的实验数据还是少之又少。因此，开发出简单合理的预测性模型来描述非水电解质溶液的热力学性质对于化工生产来说是非常有意义的。

1 理论基础

1.1 状态方程简介

本研究对象是以碱金属卤化物为溶质、以甲醇为溶剂的非水电解质溶液。因为盐溶质完全解离，所以溶液中存在3种粒子——（盐）阳离子、（盐）阴离子和甲醇分子。盐离子被看成是带电的Lennard-Jones（LJ）球体，它们身上具有可与甲醇分子形成氢键的缔合点，1个甲醇分子被看成是带有1个点偶极和3个缔合点的LJ球体。

在温度T和体积V条件下，电解质溶液系统由N个粒子组成，第i种物质的粒子数为Ni。运用微扰理论可以将整个电解质溶液系统总的亥姆霍兹自由能A微扰展开成如下形式[1]：

式中，k是玻尔兹曼常数，上角标hs、LJ、elect、assoc和chain分别代表来自硬球排斥、Lennard-Jones、静电、缔合和成链作用的贡献项，其中静电作用项（elect）包括离子-离子（cc）、离子-偶极（cd）和偶极-偶极（dd）3项。利用电介质溶液的相关热力学理论和计算机程序可以将上述各个亥姆霍斯自由能贡献项的数值计算出来，然后再将各个贡献项相加，最终就可以确定电解质溶液系统总的亥姆霍兹自由能[2-4]。

从压缩因子与密度和亥姆霍兹自由能的关系中可以得到电解质溶液系统的状态方程：

式中，ρ是粒子数密度，组分k的化学位由下式得到：

将上述状态方程运用于非水电解质溶液，就可以计算出非水电解质溶液的蒸汽压、溶剂活度和密度。

1.2 状态方程的计算功能

1.2.1溶剂活度的计算

由活度的定义出发可以将溶液中的溶剂活度定义为：

式中，as为溶剂活度，fs为溶剂的活度系数，xs为溶剂的摩尔分数；下标s代表溶剂。而溶剂的活度系数fs与剩余化学位μres的关系可以由下式确定[5]：

式中，下标0表示纯态物质。

将上式带入式（5）即可得到电解质溶液溶剂活度的表达式：

1.2.2蒸汽压的计算

将P=ρkT带入上式，上式可以变为：

由于溶剂是电解质溶液的主要组成部分，因此可以将溶剂的蒸汽压看做是电解质溶液的蒸汽压，由此可以得到电解质溶液蒸汽压的表达式：

2 状态方程的确立

状态方程中包含了一些参数，其中有些是可调节参数，有些是不可调节参数，不可调节参数主要包括甲醇溶剂的参数和盐溶质的阴阳离子参数。

2.1 不可调节参数的确定

2.1.1甲醇溶剂的参数

研究针对的非水电解质溶液的溶剂是甲醇。甲醇分子是一种带有1个点偶极和3个缔合点的LJ球体。通过本状态方程可以确定甲醇的偶极距为8.30×10-30C·m，甲醇的介电常数ε为32.49，这与298.15 K下的实验值32.7非常之接近。

对于甲醇来说，除了偶极距以外，在本状态方程中还有4个与温度无关的参数需要确定，这其中包括2个LJ参数（软球直径σ和能量参数ε/k）和2个缔合参数（缔合体积参数κassoc和缔合能量参数εassoc/k）。通过本状态方程可以确定这4个参数，结果见表1。

在确定了甲醇溶剂的4个参数的同时，本状态方程可以计算出甲醇的饱和液体在不同浓度下的密度和蒸汽压，结果见图1。

图1 不同温度下甲醇的蒸汽压和密度Fig 1 Vapor pressure and density of methanol at different temperature

表1 甲醇的状态方程参数和性质数据（298.15 K、0.1 MPa）Tabl 1 Methanol’s equation of state parameters and property data (298.15 K and 0.1 MPa)

由图1可知，该计算值（曲线）与实验值（点）吻合程度很高，其温度是283.15～333.15 K，密度和压力的平均绝对偏差分别为0.043%和0.087%。计算所得的甲醇在298.15 K时的密度为785.7 kg/m，这实验值786.2 kg/m接近。

2.1.2盐溶质的阴阳离子参数

每种盐溶质所包含的阴阳离子有2个LJ参数：离子直径σion和离子能量参数εion/k，其中离子直径σion可以从Pauling离子半径得到，离子能量参数εion/ k可以通过Mavroyannis和Stephen离散理论计算得到[6-7]。

这里αion是从文献中得到的离子极化率[8]；ne是离子中的电子数目。所有这些阴阳离子的参数都列在表2中。粒子间所有的LJ相互作用都被考虑了进去，一共包括甲醇-甲醇、甲醇-离子和离子-离子3种作用。

表2 阴阳离子参数Tab 2 Parameters of cations and anions

可以通过离子大小来预测每个离子上的缔合点数S。在通常情况下，离子越大缔合点数就越多。在本工作中，对于阳离子Li+、Na+、K+的缔合点数都取7，Rb+取8，Cs+取10；对于阴离子来说，Cl-、Br-和I-的缔合点数分别取9、10和12。每个离子的缔合体积参数为κii=0.001。

2.2 可调节参数的确定

本状态方程包括2个可调节参数：第1个参数是计算式(2)中Acd时所用到的有效平均离子直径σi，它对于所有种类的盐溶质来说都是可调节的；另1个参数就是对于每1种离子来说都不可忽略的离子-甲醇缔合能量参数εassoc/k，这个参数也可以通过本状态方程计算得到。本次工作忽略了阴离子和甲醇的缔合作用，所以在本状态方程中不存在阴离子-甲醇缔合项，因此在本状态方程中只存在2个参数。

本状态方程通过了9种碱金属卤化物的甲醇电解质溶液的测试，溶液的共同状态是0.1 MPa、298.15 K。所得溶液蒸汽压的平均绝对偏差（AAD）的计算：

式中，NP是试验点数，f是计算对象（p和a），上标cal和exp表明分别计算值和实验值。

AAD和可调节参数列在表3中，用以判断本状态方程的准确度。

表3 状态方程的可调节参数以及计算所得的蒸汽压p和溶剂活度a的平均绝对偏差（AADs）Tab 3 Adjustable parameters of this equation of state(EOS) and average absolute deviations(AADs)for vapor pressure(p)and solvent activity(a)calculated

从表3可看出，本状态方程计算出来的非水电解质溶液的蒸汽压和溶剂活度的数值与实验值的吻合度很高，二者总的平均绝对偏差分别是1.112%和0.106%。因此可以判断，本状态方程能够准确地描述非水电解质溶液的热力学性质。

除了可以计算活度和蒸汽压之外，本状态方程还可以计算非水电解质溶液的密度，与活度系数模型相比这也是本状态方程的另一个优势功能所在。但是由于历史上对于非水电解质溶液的研究工作进行得太少，因此暂时还无法在历史文献中查询到相关的密度实验值作为本状态方程所计算出的密度计算值的参照物，所以在此处就不再展示计算结果。

2.3 可调节参数的通用性

本状态方程所包含的2个可调节参数具有很高的通用性，在不改变可调节参数的前提下，针对同一种非水电解质溶液可以计算出它在不同温度和浓度下的蒸汽压和溶剂活度。举例来说，图2和图3展示了本状态方程在使用表3中可调节参数的情况下计算出的相同非水电解质溶液在不同温度下的蒸汽压和溶剂活度（可调节参数298.15 K时），计算出来的数值是更高温度状态下的数值。

图2 甲醇电解质溶液的预测蒸汽压Fig 2 Predicted vapor pressure of methanol electrolyte solutions

由图2和图3可知，此时本状态方程依然能够准确地描述非水电解质溶液的热力学性质，其平均绝对偏差列于表4中。

值得注意的是，本状态方程的可调节参数的通用性还不仅仅局限于上述情况。在不改变可调节参数的情况下，把某种非水电解质溶液的浓度提高、甚至是将浓度和溶液温度同时提高，对于其溶液压力和溶剂活度的计算值相对于实验值而言依然保持着很高的吻合度。图4展示了本状态方程使用表3中NaI-甲醇电解质溶液的可调节参数所计算出来的其在更高的浓度或者更高的温度和浓度下的溶剂活度和溶液蒸汽压，平均绝对偏差见表4所示（可调节参数用的是质量浓度0.024～0.755 mol/ kg下的参数，计算结果是0～4.338 mol/kg下）。

图3 甲醇电解质溶液的预测溶剂活度Fig 3 Predicted solvent activities of methanol electrolyte solutions

表4不改变可调节参数计算出的蒸汽压p和溶剂活度a的平均绝对偏差(AADs）（0.1 MPa）Tab 4 Average absolute deviations(AADs)about vapor pressure(p)and solvent activity(a)calculated by adjustable parameters not change(0.1 MPa)

3 结论

目前关于电解质溶液的研究成果有很多，但是大多数研究针对的对象都是以水为溶剂的水电解质溶液，而本状态方程所研究的对象是以甲醇为溶剂的非水电解质溶液。因此，对于电解质溶液领域来说具有填补空白的意义。

状态方程包含了2个可调节参数，它能够准确地计算以甲醇为溶剂、以碱金属卤化物为溶质的非水电解质溶液的蒸汽压、密度和溶剂活度等热力学性质数据。所包含的2个可调节参数具备很强的通用性——在不改变可调节参数的情况下，对于不同温度、浓度的非水电解质溶液的热力学性质依然具备十分准确的预测性。

图4 高含量下的NaI-甲醇电解质溶液的预测溶剂活度和蒸汽压Fig 4 Predicted solvent activity and vapor pressure of NaIMethanol electrolyte solution at higher molality

表5 不改变可调节参数计算出的蒸汽压p和溶剂活度a的平均绝对偏差（0.1 MPa）Tab 5 Average absolute deviations(AADs)about vapor pressure (p)and solvent activity(a)calculated by adjustable parameters not change(0.1 MPa)

[1]Boyle E J,Mcquarrie D A.Perturbation-theory of electrolyte -solutions[J].J Chem Phys,1985,83(6):3077-3082.

[2]Liu Z P,Wang W C,Li Y G.An equation of state for electrolyte solutions by a combination of low-density expansion of non-primitive mean spherical approximation and statistical associating fluid theory[J].Fluid Phase Equilibria,2005,227 (2):147-156.

[3]Liu W B,Li Y G,Lu J F,A new equation of state for real aqueous ionic fluids based on electrolyte perturbation theory, mean spherical approximation and statistical associating fluid theory[J].Fluid Phase Equilibria,1999,158-160:595-606.

[4]Liu W B,Liu Z P,Li Y G,Comparison of Perturbation Theory and Mean Spherical Approximation for Polar Fluids and Ion Dipole Mixtures Based on Molecular Simulation Data[J].Fluid Phase Equilibria,2001,178(1/2):45-71

[5]韩正.非水电解质溶液系统的状态方程研究[D].北京:北京化工大学,2008:79-81.

[6]Pauling L.The NatureofChemicalBond[M].3rd ed.Ithaca,NY: Cornell University Press,1960.

[7]Mavroyannis C,Stephen M J.Dispersion forces[J].Mol Phys, 1962,5(3):629-638.

[8]CokerH.Empiricalfree-ion polarizabilitiesofthe alkalimetal, alkaline earth metal,and halide ions[J].The journalofphysical chemistry,1976,80(19):2078-2084.

TQ013.1

A DOI 10.3969/j.issn.1006-6829.2016.04.006

2016-05-27