11计算化学-6电子相关方法共54页

电子相关的物理图像在Hartree-Fock方法中，采用单电子近似（轨道近似）—假设一个电子在原子核与其他电子形成的平均势场中独立运动。这种假设考虑了电子之间平均的相互作用，但没有考虑电子间的瞬时相关。由于Pauli不相容原理，自旋平行的两个电子不可能在空间同一点出现，基本上反映出每个电子周围有一个Fermi孔。但单电子近似导致独立运动的两个自旋反平行的电子有可能在某一瞬间在空间同一点出现。然而，由于电子间的排斥作用，这是不可能的。因此电子并不是独立运动，而是相互间有一定的制约作用，即它们的运动是相关的。每个电子周围有一个Fermi孔，应该还有一个Coulomb孔，禁止其他电子进入。电子间的这种相互关系称为电子运动的瞬时相关性或称动态相关性。单组态自恰场方法没有考虑这部分相关作用，引起相关能误差。用数学语言表达：双粒子算符不能精确地用某种等效单粒子算符之和来代替。电子相关能我们可以从两个电子同时出现的概率角度来考虑电子相关。设1(r1)是任一电子在小体积dr1附近出现的概率。2(r1,r2)是任意两电子同时分别在dr1和dr2出现的概率。则：2(r1,r2)/1(r1)就是已知一电子在r1时，在r2附近发现另一电子的概率。且r12=越小，Fr2(r1)值越负，表明环绕在r2电子周围的相关Coulomb孔存在（若两电子是独立运动，则F=0）电子间有Coulomb排斥作用，所以：2(r1,r2)/1(r1)1(r2)函数Fr2(r1)=[2(r1,r2)/1(r1)]-1(r2)0•由于HF单组态自恰场未考虑Columb相关,忽略了电子之间的排斥作用(排斥使得能量升高),因此,相关能是负值,将导致体系计算能量更加负•Lowdin的定义:指定的一个Hamilton量的某个本征态的电子相关能,是指该Hamilton量在该态的精确本征值与它的限制的Hartree-Fock极限期望值之差。•因为相关能等的存在，计算得到的体系总能量的理论值永远高于其真实值。相关能的定义无限(完备)基组变分计算有限基组变分计算近似HF极限值HF极限值相关能理论值总能量实验值精确非相对论极限值相对论效应校正基组误差近似相关能实验值基组误差相对论效应校正(近似)近似相关能计算值近似非相对论极限值各种能量值的关系相关能在化学问题中的重要性相关能在体系的总能量中所占的比例并不大，大约只占0.5%左右，甚至小于非相对论带来的误差。从总能的相对误差来说，HF方法是一个相当好的近似。遗憾的是，化学和物理过程涉及的是能量差值，非相对论误差在这一过程中可相互抵消，总能量差值受相关能影响很大！相关能的数值与一般化学反应的生成热、活化能等属于同一数量级，对于弱相互作用，有时甚至大一个数量级。因此，计算化学反应、构建精确的势能面、分子间弱相互作用或对分子性质做精确计算时，都需要考虑相关能。解决相关能问题成为理论化学领域中最活跃的方向！！一般处理方法•把r12包含在波函数中–适用于非常小的体系–有很多难计算的积分•后SCF方法（Post-SCFMethods）–组态相互作用(CI)及其近似方法–微扰理论（PT）–耦合簇（CC）理论组态相互作用方法（ConfigurationInteraction,CI）背景：•HF方法把电子按从低到高填充，是合理的，但是有一定的人为性•从统计力学的角度来看，任何填充方式都是可能的，只不过其在真实体系中的几率不同，HF的填充方式是几率最大的一种•因为把电子都限制在了HF的占据轨道上，所以使得不同自旋的电子过分接近，导致排斥能升高•为了改进这种硬性限制，允许电子在所有可能的轨道上自由填充，再确定出每种填充方式的几率。CI是一种严格的线性变分方法。对于N电子体系，在给定基组下，HF计算产生一组单电子轨道和电子在轨道上的一种排布。考虑电子相关效应也应允许其它的排布，即其它组态的存在。一般用HF态作为零级波函数，允许电子任意激发，因此完备的组态相互作用（FCI）波函数可以写为：组态相互作用上式是用二次量子化符号标记的，表示在轨道a上产生一个电子，表示从占据轨道i上湮灭掉一个电子，C是系数。上式也可写为：aaia5-1其中，是Hatree-Fock组态，其它组态可认为是电子从H-F态中的占据轨道激发到空轨道而产生。例如，被认为是一个电子从中的占据轨道i激发到空轨道a产生的单激发组态是两个电子从中的占据轨道i、j激发到空轨道a、b产生的双激发组态，等等。如果包括所有可能的激发，则叫做完全组态相互作用（FCI）。ai)(0000abij)(0FCI能够得到全部的相关能，在给定的基组下，是目前最为精确的理论方法。FCI能量满足变分上限性质，具有大小一致性大小一致性是指体系的总能量等于其无相互作用的子体系的能量之和。即对距离很远的两个体系,有E(A+B)=E(A)+E(B)完全组态相互作用（FCI）由于FCI展开式中的项数过多，当基集较大、电子数较多时，激发组态很多，FCI计算的工作量非常大，所以FCI计算只能对电子数不多、基组不大时才能实现。FCI的限制和部分CI部分CI最直接限制CI展开项的方法是限制体系的激发能级。一般来说，双激发态对CI波函数贡献最大。因此作为近似，（5-1）式通常取到第三项截止，就是仅考虑单、双激发的CI计算。研究表明，仅包括单、双电子激发的CI计算，计算结果可获得分子平衡态95%的相关能。仅包括单激发态的CI计算记为CIS这部分对波函数或能量的影响不大。一般来说至少要作到双激发，记为CIDΨCID=a0ψ0+aabijψabijijaboccuociaoccuocΨCIS=a0ψ0+aaiψai还可以计算包括单、双激发的CISDΨCISD=a0ψ0+aaiψai+aabijψabijijaboccuociaoccuoc虽然单激发态单独对能量、波函数影响不大，但单、双激发态之间能量矩阵元不为零，因此计算得到的相关能比双激发态有所改进。CISD对计算基态的相关能非常有用CISD共有O2V2个行列式(O=占据轨道的数目,V=未占据轨道的数目)CISDT–单重,双重和三重激发态仅限于小分子,大概O3V3个行列式完全CI–所有可能的激发态((O+V)!/O!V!)个行列式在基组给定的情况下,得到精确的相关能大概限于14个轨道上分布14个电子由于组态函数随激发程度增加而迅速增加，一般只考虑到双重激发，最多能考虑三重、四重激发，计算已经相当困难了。几种CI计算的相关能百分比（单位：%）分子CISDCISDTCISDTQFCIBH94.91（568）99.97（28698）（132686）HF95.41（552）96.49（6712）99.86（48963）（944348）H7+96.36（1271）96.87（24468）99.96（248149）（2923933）H2O94.48（1311）95.85（27026）99.85（332491）（94165610）NH394.44（2443）95.43（52595）99.84（619235）（48642057）括号内为组态波函数CSF数目N电子体系基函数随CI激发态等级而急剧上升冻芯近似冻芯近似是讨论电子相关方法中最常用的方法。由于化学反应主要涉及原子的价层电子，那些处于能量最低部分，占据原子内层的电子可以冻结起来。对于第二周期元素，冻芯即冻结1s电子轨道，第三周期元素，则冻结1s、2s、2px、2py、2pz轨道。由冻芯引起的误差，对含有相同原子的体系有基本相似的误差。为简化计算，可把分子轨道基组决定的轨道空间分成几个不同类型的子空间，计算时冻结内层空间与外层虚空间；即不考虑内层电子激发，也不考虑外层虚空间电子占据。电子激发仅在价层电子占据与价层虚空间之间进行。实际计算表明，若组态函数选得好，只要取其中较少一部分，就可获得大部分的相关能。所以选择CI组态函数是进行CI计算的重要技巧，既不能把有重要作用的组态函数漏掉，又要使选用的组态函数尽可能地少。部分CI计算中组态函数选择的三种情况只考虑一部分能级较高的占据轨道和一部分能级较低的空轨道两种极端情况：中间是只计算HOMO与LUMO轨道，右边的只冻结内层核心轨道，其余轨道全部计算。只可用于仅含第一、二周期元素的小体系多体微扰理论（Many-BodyPerturbationTheory,MBPT）基于Hartree-Fock参考态的Möller-Plesset二级微扰理论（MP2）是计算电子相关能最简单、应用最广泛的方法。该方法对平衡位置附近的闭壳层体系是可靠的。1.Brillouin-Wigner微扰展开多体微扰理论将体系的Hamilton算符分解成两部分：式中，Ĥ0是无微扰的Hamilton算符，为微扰项。Schrödinger方程可以表示为：波函数向一个正交归一化集合展开:波函数与能量可以表示为：E=E0+E(1)+E(2)+…能量校正到一级微扰:二级微扰:Mller-Plesset微扰能BW微扰展开从形式上看是一个封闭公式，但实际运算时由于E是未知，跌代要有一初始值，如果用近似值，效果不能保证。Mller-Plesset在此基础上做Rayleigh-Schrödinger展开处理，现在从头算程序里的微扰能即以这种方法编程。即微扰项定义为MP理论中，零级微扰哈密顿取单电子Fock算符的和。本征值Es对应占据轨道单电子能量i之和微扰波函数和能量按展开，可表示为：iFHˆˆ0MP一级微扰能对应HF能，即：HF方程中的交换作用,就是MPn的一级表示HF方法是MPn微扰法的一级近似二级微扰能MP2为：E(2)=-∑∑(a+b-i-j)-1|(ij||ab)|2ijoccabuoc这里求和是对所有双激发态求和。MP2只需要部分双电子积分变换为自旋轨道基，计算速度比较快。MP3中的三级微扰能为：E(3)=-∑∑(E0–Es)-1(E0–Et)-1V0s(Vst-V00st)V0tsDtD求和是对所有双激发态两次求和，再加上所有三重激发态之和。在MP4中单激发、双激发、三重激发、四重激发态都出现，计算相当复杂。特别是三重激发态计算相当困难，所以有些计算程序只计算单、双、四重激发，关键词为MP4SDQ。在密度泛函方法得到广泛应用之前，MP2方法是考虑电子相关的最便宜的方法，它可以成功应用于很多领域，一般都能得到很精确的结果，是理论化学中非常有力的工具。当然，MP2也有不能适用的体系，一般来说体系越特殊，所需要的计算理论水平就越高。在MP2不能得到好的结果的时候，就可以使用高等级的MP方法。在实践上，只有MP4得到比较广泛的应用，MP3对于MP2处理不好的体系一般也没有好的结果。基于与MP3同样的原因，MP5也很少被使用。虽然MP4比MP2昂贵很多，但其在MP2所不能解决的问题方面，的确能得到很精确的结果。下面是对一些MP方法的比较HCNCN-CN.MP2-91.82033-91.07143-91.11411MP3-91.82242-91.06862-91.12203MP4-91.82846-91.07603-91.13538MP5-91.83129-91.07539-91.14221MP6-91.83233-91.07694-91.14855MP7-91.83264-91.07678-91.15276MP8-91.83289-91.07699-91.15666FullCI-91.83317-91.07706-91.17006dE0.001atMP6MP6MP19FullCI-MP4-2.96-0.65-21.76(kcal/mol)计算都采用STO-3G基组。对于HCN和CN-，MP方法迅速收敛，在MP6等级，能量差就小于0.001；而对于自由基CN，要在MP19的等级上，等级之间的能量差才小于0.001.•在每个阶次都是大小一致的•MP2–二阶计算相对容易(只需要双激发态)•在Hartree-Fock是一个相当好的起点时,二阶将可以得到很大一部分相关能•实际使用中可以到四阶(单