Chen Changguo，Zhang Hongmin，Huang Zongqing
(Department of Applied Chemistry， Chongqing
University，Chongqing 400044)
Key words Quantum Electrochemistry, Electrode/Solution Interface
摘要 从电极表面的簇模型、电极/溶液界面上的吸附模型和吸附体系、双电层电容、电极电势的费米球概念及理论计算、电催化氧化还原等方面综述了近年来量子化学理论和方法在电化学研究中所取得的进展。
关键词 量子电化学 电极/溶液界面
电极/溶液界面的量子化学研究进展
陈昌国 张洪敏 黄宗卿
(重庆大学应用化学系 重庆
1 电极表面的簇模型
对于金属电极，量子化学研究中广泛采用金属簇模型，即用一个或几个按特定几何构型排布的基质原子组成的簇来类比电极表面[2]。如图1所示的簇模型在不同程度上模拟了铂金属表面[3]。合理的簇模型应能体现不同晶面的表面结构和吸附位结构，应能反映大块金属特殊的电子性质[4]。簇模型的核心在于将固体的电子结构作定域化描述。较好地定域化在于能否从解有限簇模型的Hartree-Fock方程中得出无限的固体所具有的一些特性，这也是簇模型的优点，它使固体的电子结构和局部对称性简化为确定分子的对称性和电子结构这类量子化学易于处理的问题[5]。
图1铂金属表面的簇模型
用簇模型进行量子化学研究时，存在两种误差：一种来自簇模型对金属表面的近似；一种来自量子化学计算方法的近似。金属表面原子的电子云重叠非常大，不能忽略簇原子以外周围原子的作用。原则上，原子簇越大越能体现金属性，但加大原子簇需付出更多的计算时间，甚至为了得到计算结果不得不采用精确度降低的近似计算方法。小簇虽能采用精确度较高的计算方法，但其与表面之间的共性值得怀疑。一种解决办法是采用徐昕等人[2]提出的“金属态”原子概念，对自由原子的Slater指数xa进行修正，导出金属中原子的Slater指数xm。如Cu原子的xa、xm值可选择为表1所示。对CO在铜上的吸附，当用一个Cu原子来模拟铜表面，其计算结果也能较好地解释UPS、XPS等电子能谱数据。解决簇模型本体近似问题的另一种办法是将原子簇分为活性位附近的内部区及其周围的外部区。内部区和外部区的原子采用不同的基集和近似处理来兼顾本体性质与计算精度。
| 表1 Cu原子的xa、xm[6] | ||||||||||||||||||||
|
由于电极表面常有氧化膜存在，量子化学研究中也采用均匀分布的点电荷来模拟晶体中Madelung势场的影响[7]。
总之，簇模型理论本身还在不断地发展，它为量子化学处理具有无限结构的固体提供了一条有利途径。
2 电极/溶液界面上的吸附模型
离子、分子等在电极表面的吸附位一般有顶位(on Top)、桥位(on
Bridge)和穴位(on Hole)等三种基本形式(图2)。在不同的晶面上，
相同吸附位的环境(对称性)也不一样， 如顶位吸附就有图3的几种情况。
图2 吸附位的三种基本形式
图3 在不同晶面上顶位吸附的对称性
电极表面吸附的量子化学处理与一般固体表面相比还须考虑两个因素：界面强电场和溶剂极性的影响。前者一般通过模拟电极电势的变化来实现。田昭武等[4]根据关系式：
提出了“给簇加上过剩电荷dQ来移动其费米能级ef，按-def =e.dE定量模拟电极电位的变化dE”的方法，并计算了CO在Pt上的吸附。现场红外光谱测得CO的谱峰频率与电位的关系定量地符合量子化学计算结果。
电解质和溶剂影响的处理，或以共吸附、偶极取向、势场模型引入计算，或者在某些场合作为一级近似忽略之。近几年来，溶剂效应的量子化学计算方法也逐步发展起来，主要有两种模型[8]：自洽场理论模型SCF和量子力学-分子力学相结合的QM-MM模型。前者反映的是溶剂影响的整体效应，后者可将溶质-溶剂相互作用分解为静电作用、范德华作用和氢键作用等。后面将介绍的水分子和离子在电极表面上的吸附研究也有助于溶剂作用的计算机模拟。
3 电极/溶液界面吸附的研究
3.1 水分子的吸附态
由于水溶液电解质居多，水分子在电极表面的吸附研究就显得尤为重要。Muller[9]采用Khan-Sham方法发现H2O通过O在Al簇顶位吸附时，H-O-H平面与Al簇法线夹角为55°，且O距Al簇为3.9a.u.时，能量最低。Seller[10]采用了相对论核芯势(RECP)ab initio和微扰理论校正，得出了一个水分子在汞表面的顶位、桥位和穴位吸附(图4)的键能分别为13.1，12.2和11.6
kcal/mol，而O与Hg表面平衡距离分别为5.23，4.89和4.86a.u.。
图4 水分子在汞表面上的吸附
q 为水分子偶极矩(m )与汞表面法线(n)之间的夹角。
Renat[11]进一步用Gaussian 90对Hg7-H2O吸附体系进行了计算，得到吸附水分子的偶极矩为2.18D，接近于实验值1.82D，并证明了水分子在汞表面吸附时偶极矩垂直向外时能量最低。
3.2离子的吸附
量子化学对离子吸附的研究近来也较活跃。Blance[12]用ab initio和CNDO/2相结合对比研究了卤素原子X和离子X- 在汞表面的吸附，对内部区和外部区的汞原子分别采用了不同的基集。结果表明卤原子和离子的吸附非常相似，只是离子吸附具有较大的键长和较小的键能(F- 除外)。
德国学者[13]用ab initio计算获得了I-、Li+与Pt表面作用的解析势函数：
VPt-I- =6.30x106exp(-4.59r)- 4.66x106exp(-0.80r)
VPt-Li+ =1.87x103exp(-1.33r)-
1.23x103exp(-1.02r)+1.06x104/ r8.350
其中r为I-、Li+与Pt表面的距离。同时还发现I-在Pt(100)上的吸附键能比Li+大，两者的优先吸附位都是穴位。
林文锋[14]等用电荷自洽离散变分Xa方法(SCC-DV-Xa)研究了CN-在Ag4(1,3)上的吸附，得到了电子集居数和电荷分布(见表2)。每个CN-离子反键轨道向Ag电极传递了0.16个电子电荷而形成C-Ag吸附键，同时增强了C-N间的键合。电位负移则抑制电荷传递而削弱了吸附作用。
| 表2 CN- 的电子集居数和电荷分布 | ||||||||||||||||||||||||||||||||||||||||||||
|
3.3 缓蚀剂的研究
除硫脲之外，目前对缓蚀剂的量子化学研究主要集中在芳环及杂环中含氮、氧、硫等杂原子的衍生物以及脂肪醇、脂肪胺等化合物上。
3.4氢等小分子气体的吸附
电极过程常常涉及气体的吸附和析出，尤其是H2、O2、CO等小分子气体在燃料电池以及电催化氧化还原反应中举足轻重，人们因此也越来越重视从微观结构上阐明它们的吸附机理。文献[17]采用EHMO方法得到了H2在(111)、(110)和(100)等不同镍晶面上离解吸附的势函数曲线都呈Morse函数形式，根据结合能还确定了(110)和(111)为优先吸附面，而顶位、桥位则是优先吸附位。
为了搞清在银表面起催化作用的是分子态氧还是原子态氧， 王文宁等[18]分别用原子簇模型(CM)和浸入吸附原子簇模型(DAM)的从头计算法研究了Ag6O2吸附体系。结果表明两个最低能态是1A1和3A2，它们分别对应在[110]槽位吸附的过氧分子O22-和在[001]方向上吸附的超氧分子O2-，理论优化的吸附构型和实验测量结果一致。
CO在金属表面上的吸附由于与有机物的催化氧化还原密切相关，而且也是碳一化学的重点研究对象，报道较多，这里不再叙述。
3.5表面共吸附
由于电化学体系涉及的微粒物种众多，共吸附在所难免。量子化学研究共吸附能够详细地描绘各个物种之间的相互作用。以CO为探针分子，可以研究S和Bi在Pt表面催化反应中起的作用。研究表明，S原子从Pt电极上得到电子，降低了Pt对CO的电子给予能力，使CO的氧化变得容易进行。另外，S与CO之间有排斥作用，这两者都使CO吸附减弱，nco上升，与现场红外光谱的实验结果一致。但Bi的情况则恰恰相反[19]。
Bonzel等[20]用ab initio对K+与H2O在Pt表面共吸附模型Pt10-K-H2O进行了计算。结果表明，在与K+共吸附时，H2O分子平面倾斜160° 时最稳定。这种不寻常的构型是由于H2O-K的近程作用和静电作用引起的。
4 双电层电容
van den Eeden等[21]根据簇模型用ab initio计算方法研究了金属引起的电容CM。双电层电荷密度的分布是场强E的函数。在电场存在下，体系的哈密顿算符变为：
其中：Zm 是电场方向上第m个电子的坐标。算符EZm加在单电子算符上，由此计算出的波函数反映了体系受外界的所有影响。CM的计算式为：
其中：r是在电场E作用下r处的电荷密度。 研究结果表明线性Li簇的CM值随原子数n增大而逐渐收敛，而三维簇模型用不同的近似所得到的CM大致相同。 他们进一步研究了Li(111)20(7-3-3-7)，Li(100)22(9-4-9)，Li(110)22(9-4-9)三种簇， 发现金属晶面对表面电荷分布和CM有较大影响。
5 电极电势的费米球概念以及理论计算
电化学体系的特点就是当金属与不同电解质接触时，其费米能级ef发生不同程度的变化。史美伦[22]从量子统计处理固体电子状态的费米面和费米球概念出发，对电极电势的本质进行了讨论。所谓零电荷电势(PZC)相应于费米球外电子完全被剥离的状态，即固相的裸费米球。电极的热力学电势与平衡费米球的半径相联系，动力学性质与费米球半径的扩大和缩小以及球外电子溅入费米海洋或费米海洋表面电子逸出费米球的数量相联系。
Jaleli-Heravi等[23]用ab initio和AM1二种方法计算了苯醌和萘醌类七种物质的氧化还原电势， 结果表明它们的电极电势7%的变化来自于分子相互之间的空穴作用，15%的变化来自于扩散作用。分子轨道对电极电势影响电化学反应的解释是[24]:随着电势改变，电极价带移动，当电极表面功函数的变化等于电势变化时就达到了平衡。该理论可以解释一些吸附体系的性质随电极电势变化的关系。
6 电催化氧化还原过程
过渡金属配合物中的CoTPP和FeTPP对氧在水溶液中的电化学还原催化活性较好，但电化学方法难以区分电子是转移到中心离子还是卟啉环上。文献[25]用限制性CNDO/2研究了不同还原阶段体系的HOMO和LUMO的组成。结果表明Fe(III)TPPCl的HOMO轨道中dxz和dyz占95%，故第一个电子被Fe的d轨道接受，而FeTPPCl-的LUMO主要由卟啉环上N、C的轨道组成，故第二个电子传递到卟啉环上。同样的分析表明，CoTPP还原时二个电子均与Co的d轨道发生作用。
用量子化学方法研究电化学界面能够从电子水平上探索电化学过程。近十余年来，我国在量子电化学领域的理论研究已有一些报道， 亦取得了一些高水平的成果。但对于电极/溶液界面来说，很多复杂的问题一时还难以用量子化学方法来解决，或者是在研究中引入了过多的假设或简化。因此，如何提高量子化学研究的准确性以及推广量子化学方法在电化学中的应用值得进一步研究。
7 参考文献
[1] Bockris J O’M, Khan S U M. Surface Electrochemistry:A Molecular Level approach. Plenum，New York and London，1993.
[2] 徐昕，王南钦，张乾二. 金属上化学吸附的簇模型量子化学研究. 化学通报，1994，(4)24-21.
[3] Lin W F， Sun S G，Tian Z Q et al. Quantum chemistry and in situ FT-IR spectroscopic studies on potential-dependent properties of coadsorbed on Pt electrode. Elctrochim Acta，1993，38(8):1107-1109.
[4] 林文锋，孙世钢，田中群 等. 电化学界面吸附微观模型和量子化学处理.科学通报，1993，38(24):2252-2254.
[5] Jaroslav K, Piercarlo F. Theoretical aspects of metal atom cluster.Chem Rev，1986，86:539-545.
[6] Xu X，Wang N Q, Zhang Q E. Studies on surface analogy:ab initio calculations for the CO/Cu chemisorption system. Surf Sci，1992，274:386-390.
[7] 吕鑫，徐昕，王南钦等. NiO(100)表面小分子(Co，NO，O2)吸附DV-Xa嵌入簇模型研究。高等学校化学学报，1995，16(3):424-427.
[8] 蒋华良，陈凯先，嵇汝远. 溶剂效应量子化学研究进展. 化学通报，1995，(4):1-5.
[9] Muller J E, Harris. Cluster study of the interaction of a water molecule with aluminum surface. J. Phys Rev Lett， 1984，53(26):2493-2496.
[10] Sellers H, Sudhakar P V. The interaction between water and the liquid-mercury surface. J. Chem Phys, 1992, 97(9):6644-6648.
[11] Renat P，Nazmutdinov S, Probst M et al. Quantum chemical study of an H2O molecule on an uncharged mercury surface. J. Electroanal Chem，1994，369:227-231.
[12] Blanco M, Rubio J, Lllas F. An initio cluster model approach to the chemisorption on mercury:halogen and halides. J Electroanal Chem，1989，261:39-46.
[13] Seitz-beywl J, Poxleitner M, Probst M et al. On the interaction of ions with a platinum metal surface. J Quantum Chem，1992，42:1141-1147.
[14] Lin W F，Tian Z Q, Sun S G et al. Quantum chemistry studies on electronic properties of CN-. Electrochim Acta， 1992，37(2):211-216.
[15] 董俊华，宋光铃，曹楚南等. 铁电极上硫脲及衍生物的缓蚀作用研究. 物理化学学报，1996，12(3):252-257.
[16] Romanowski S, Pietrza T M. On the role of hydrogen bonds in thiourea condensates-the correlation of quantum-chemical approach with the experimental results. J Electroanal Chem，1993，345:83-90.
[17] 王文宁，范康年，邓景发. 分子氧在Ag(110)面上的吸附构型、吸附态和吸附能的CM和DAM研究. 化学学报，1995，53:1000-1006.
[18] 苏文煅，陈子云，杨旭辉等. 氢在镍晶面上的吸附态研究。厦门大学学报(自然科学版)，1987，26(3):317-321.
[19] Lin W F, Sun S G, Tian Z W. Investigations of coadsorption of carbon monoxide with S or Bi adatoms at a platinum electrode by in-situ FT-IR spectroscopy and quantum chemistry analysis. J Electroanal Chem, 1994, 364:1-5.
[20] Bonzel H P, Pirug G, Muller J E. Reversible H2O adsorption on Pt(111)-K:work-function changes and molecular orientation. Phys Rev Lett, 1987, 58(20):2138-214.
[21] Van den Eeden A L G , Sluyters J H and van Lenthe J H.Quantum chemical calculations on the electric double layer at the mental/vacum interface.Part one:Hartree-Fock cluster description. J Electroanal Chem, 1986, 208: 243-248.
[22] 史美伦. 费米面与电化学(一).同济大学学报，1981，(3):16-21.
[23] Jalali-Heravi M, Namazian M, Peacock T E . Theoretical studies of electrode potentials in aqueous solution:Investigation of individual contributions from electrostatic，cavity and dispersion interactions to redoxpotential. J Electroanal Chem，1995，385:1-7.
[24] Anderson A B. The influence of electrochemical potential on chemistry at electrode surface modeeled by MO theory. J Electroanal Chem， 1990，280:37-41.
[25] 郭孝纯，黄楚宝，郭富春 等. FeTPP和CoTPP电化学氧化还原的量子化学研究. 高等学校化学学报，1991，12(9):1242-1246.
陈昌国 男，
40岁，教授。主要从事电化学和结构化学研究。