电子管理与计算化学

电子管理与计算化学电子相关能•在Hartree-Fock近似下,每个电子感受到的是所有其它电子的平均密度•但是,在某个瞬间,两个电子不能出现在同一地点•电子之间必然两两之间互相避免,即它们的运动是相关的•给定基组的情况下,相关能就是精确能量与Hartree-Fock能量之差•每个电子对的相关能大概是20kcal/mol电子相关的表现•对于闭壳层组态,电子相关性分解为两个部分:自旋相同的电子之间和自旋相反的电子之间•f(r1,r2)=[fab(r1,r2)+faa(r1,r2)]/2•faa(r1,r2):Fermi孔;fab(r1,r2):Columb孔•HF方法基本上反映了Fermi孔(由于波函数的反对称性),完全没有反映Columb孔!正六角形H6分子的相关函数(a)r2在氢核上(b)r2在氢核中间HF的相关能是精确的(Exact)!!!•由于HF方法主要未考虑Columb相关,忽略了电子之间的排斥作用(排斥使得能量升高),因此,相关能是负值,将导致体系计算能量更加负•Lowdin的定义:指定的一个Hamilton量的某个本征态的电子相关能,是指该Hamilton量在该态的精确本征值与它的限制的Hartree-Fock极限期望值之差。相关能的定义无限(完备)基组变分计算有限基组变分计算近似HF极限值HF极限值相关能理论值总能量实验值精确非相对论极限值相对论效应校正基组误差近似相关能实验值基组误差相对论效应校正(近似)近似相关能计算值近似非相对论极限值各种能量值的关系相关能在化学问题中的重要性•理论和实验能量差值的比较–实验误差:1kcal/mol~0.002Hartree–体系总能量101~5Hartree–体系相关能10-1~1Hartree–总能量差值受相关能影响很大。•弱相互作用（比如分子间作用）–氢键键能～3－10kcal/mol一般处理方法•把r12包含在波函数中–适用于非常小的体系–有很多难计算的积分–He的Hylleras波函数•用一组多电子波函数来表示体系的波函数–Hartree-Fock行列式和激发态行列式–有非常多的激发态行列式,收敛很慢–组态相互作用(CI)相关方法的目标•可以明确地给出(welldefined)–不需要人为选择就可以应用于所有分子–可用于构建模型化学•效率–不能仅限于非常小的体系•变分性–是精确能量的上限•大小一致性(sizeextensive)–E(A+B)=E(A)+E(B)–需要对热化学进行正确描写•计算成本与精确性的等级性–这样计算才能被系统改进组态相互作用方法•HF方法把电子按从低到高填充,是合理的,但是,有一定的人为性•从统计力学的角度来看,任何填充方式都是可能的,只不过其在真实体系中的几率不同,HF的填充方式是几率最大的一种•因为把电子都限制在了HF的占据轨道上,所以使得不同自旋的电子过分接近,导致排斥能升高•为了改进这种硬性限制,允许电子在所有可能的轨道上自由填充,再确定出每种填充方式的几率。组态相互作用00111111111||||||(,)aaabababcabciiijijijkijkiaijabijkabcnaiiainabijiaijbjniajbttt参考行列式(Hartree-Fock波函数)单激发态行列式(把占据轨道i用未占据的轨道a代替)双激发态行列式等等•用变分原理来确定CI系数•CIS–包括所有单激发态–用于激发态计算,不是用于基态相关能•CISD–包括所有单激发和双激发态–对计算基态的相关能非常有用–O2V2个行列式(O=占据轨道的数目,V=未占据轨道的数目)•CISDT–单重,双重和三重激发态–仅限于小分子,大概O3V3个行列式•完全CI–所有可能的激发态–((O+V)!/O!V!)2个行列式–在基组给定的情况下,得到精确的相关能–大概限于14个轨道上分布14个电子组态相互作用tddEtttabcijkijkabcabcijkabijijababijaiiaairespecttowith/ˆminimize**0H组态相互作用•这是很大的求本征值的问题,可以用迭代法求解•CI是行列式的线性组合•得到的是精确能量的上限(变分性)•可应用于激发态•比起非变分方法,梯度可以简单计算出来ttHE•对距离很远的两个体系,有E(A+B)=E(A)+E(B)•对有2个基函数的He原子,CID为•对两个He原子分别计算得到–E=E(A)+E(B),=AB•两个He原子距离很远时,CID为(没有包括四重激发态)大小一致性AAAAAt22110BBAABBAABBBBAAAABBBBBAAAAABABABAtttttt222211112211221102211022110))((BBBBAAAABAtt221122110多组态自洽场方法(MCSCF)•CI中采用的轨道是自洽场的分子轨道,(占据轨道以及未占据的轨道),只对c0,c1等进行了变分,并没有保证使用的分子轨道是最优的•多组态自洽场方法,不仅对c0,c1等进行变分,还对分子轨道内部的基组系数c11,c12等进行变分,这样就能保证分子轨道最优的情况下,求最优的行列式系数•当然,它的计算量要比CI还要大!组态相互作用方法的简化•CI方法的实现就是降低行列式的个数(2K!)/[N!(2K-N)!],即减小K和N•只允许低等级激发,一级激发(CIS,N(2K-N)项),二级激发(CID,N(N-1)(2K-N)(2K-N-1)/[2!2!]项);一级+二级激发(CISD)•完全活性空间自洽场方法(CASSCF):把所有自洽场分子轨道分为内层(永远双占据),外层(永远不占据)和活性轨道,只对活性轨道进行MCSCF计算,这样,同时减小了K和N。能量0CAS-SCF内层活性空间只有价轨道只有HOMO+LUMO•把Hamilton量分为精确可解的部分和微扰部分•按l展开Hamilton量,能量和波函数•把l幂次相同的项合并在一起微扰方法VHHHHˆˆˆˆˆ00000iiiEE221022100ˆˆˆlllllEEEEVHHdEEEEdEEEE10202112012000101100100000ˆˆˆˆˆˆˆVVHVVHH•选择H0,其本征函数就是分子轨道波函数•按照Hartree-Fock行列式,单激发,双激发和多重激发行列式的方式展开微扰波函数•对能量的微扰校正为Møller-Plesset微扰理论iFHˆˆ0bajijiaaabijHFHFHFMPHFdEdEEEEddEEE,2010220000010]ˆ[ˆˆˆVVVHabcijkijkabcabcijkabijijababijaiiaaiaaa1HF方法是MBPT方法的一级近似•HF方程的Hamilton量就是MBPT的零级表示•HF方程中的交换作用,就是MBPT的一级表示Occupied(0)0011;NiiiiHFE(1)011111(0)(1)HF00111();[(|)(|)]2NNNNNjjijijijiijijVJKEiijjijijrrEEE•在每个阶次都是大小一致的•MP2–二阶计算相对容易(只需要双激发态)•在Hartree-Fock是一个相当好的起点时,二阶将可以得到很大一部分相关能•实际使用中可以到四阶(单,双,三和四重激发态)•MP4能够得到大部分的相关能•整个系列有点振荡(偶数阶更低一些)•如果出现了严重的自旋污染或者Hartree-Fock不是一个好的起点,相关能的收敛很困难Møller-Plesset微扰理论HFMP2MP3MP4MP5MP6MBPT方法的说明•MBPT可以有MP2,MP3,MP4等各级修正,主要应用的是MP2,MP4,而MP3给出的修正往往很差•MP2是从头算最经常使用的可承受包含相关能的方法,最常用的组合就是MP2/6-31G*,但是在基函数超过500后就很慢了•MBPT方法有可能给出不体系能量低的值,因为它没有变分原理保证其高于体系真实能量•CISD可以写作•T1和T2生成了相应的系数组合的所有可能的单激发态和双激发态•耦合簇理论的波函数耦合簇理论021)ˆˆ1(TTCISDabijbajiabijt,02ˆT021)ˆˆ1exp(TTCCSD•CCSD能量和波幅可以通过解下式得到耦合簇理论0)ˆ(0)ˆ(0)ˆ(0dEdEdECCCCabijCCCCaiCCCCHHH•CC方程在波幅上是二次的•必须用迭代法求解•每次迭代的计算成本类似于MP4SDQ•CC是大小一致的•CCSD–单和双激发态校正,直到无穷大阶次•CCSD(T)–把三重激发用微扰的方式加上•QCISD–组态相互作用,且添加足够的二次项以使得其大小一致(等价于CCSD的截断形式)耦合簇理论总结•CI是最直接的方法•MP是最常用的方法•耦合簇方法(CC)是最有效的计算相关能的方法分子模拟与计算化学模型化学ModelChemistries模型图(Pople图)MinimalSTO-3GHFMP2MP3MP4QCISD(T)...FullCISplit-Valence3-21GPolarized6-31G*6-311G*Diffuse6-311+G*HighAng.Mom.6-311+G(3df,p)…N2N5N6N7SchrödingerEquationBasis模型化学•模型化学理论是一个完整的算法,用来计算任何分子体系的能量•它在应用过程中没有任何主观因素•它必须是大小一致的,这样每个分子的能量才能被唯一确定•简单的模型化学使用单一的理论方法和基组•复合的模型化学把几个理论模型和基组联合起来,在低的计算成本下得到更精确的结果•模型化学的作用在于用我们容许的计算量对某些种类的分子得到最精确的结果确定方法的可信度•实验数据为确定某个理论方法的可信度提供了条件•G2和G3方法已经可以达到1kcal/mol的精度(化学精度)•G2包括了原子化能,电离能,电子亲和能和质子亲和能等热化学数据,体系中包含1-2个非氢原子•这些数据集合可以用来建立衡量理论水平怎样把我们的方法用于未知体系•如果我们研究的体系没有实验数据,就需要把理论方法逐级改进,直到达到可信的结果•从模型图的左上角开始,逐级向右下角改进,进行几种方法的计算•当计算结果没有什么改进了,我们的结果就可信了•注意:这种方法不能用于DFT理论,虽然它在基组方面可以改进确定模型的理论方法•几何结构很容易改进到可信范围,而能量比较困难(因为它取决于相关能的改进)•因此,一般在低理论水平优化结构而在高理论水平计算能量值•采用的方法是–能量计算方法/基组//优化方法/基组•例如:CCSD(T)/6-311G(2d,p)//HF/6-31G模型化学比较(键长)•108个主族分子在6-31G(d,p)下计算得到的键长与实验值的平均偏差键长HFMP2LDAGGA杂化泛函误差(A)0.0210.0140.0160.0170.011模型化学比较(原子化能)AEHFMP2LDAGGA杂化误差(Kcal/mol)119.222.052.27.06.8•108个主族分子在6-31G(d,p)下计算得到的原子化能与实验值的平均偏差模型化学比较(电子亲和能)•F的电子亲和能计算值(eV)HFMP2B3LYPSTO-3G-10.16-10.16-9.013-21G-1.98-1.22-0.866-31G(d)-0.39+1.07+1.056-311+G(2df,p)+1.20+3.44+3.46