Lec03电子相关效应

《Lec03电子相关效应》由会员分享，可在线阅读，更多相关《Lec03电子相关效应（67页珍藏版）》请在金锄头文库上搜索。

1、分子模拟与计算化学,电子相关效应 Electron Correlation Effects,电子相关能,在Hartree-Fock近似下, 每个电子感受到的是所有其它电子的平均密度 但是, 在某个瞬间, 两个电子不能出现在同一地点 电子之间必然两两之间互相避免, 即它们的运动是相关的 给定基组的情况下, 相关能就是精确能量与Hartree-Fock能量之差 每个电子对的相关能大概是 20 kcal/mol,电子相关的表现,对于闭壳层组态, 电子相关性分解为两个部分: 自旋相同的电子之间和自旋相反的电子之间 f(r1,r2) =fab(r1,r2)+ faa(r1,r2)/2 faa(r1,r2

2、): Fermi孔; fab(r1,r2): Columb孔 HF方法基本上反映了Fermi孔(由于波函数的反对称性), 完全没有反映Columb孔!,正六角形H6分子 的相关函数 (a)r2在氢核上 (b)r2在氢核中间,HF的相关能是精确的(Exact)!,由于HF方法主要未考虑Columb相关, 忽略了电子之间的排斥作用(排斥使得能量升高), 因此, 相关能是负值, 将导致体系计算能量更加负 Lowdin的定义: 指定的一个Hamilton量的某个本征态的电子相关能, 是指该Hamilton量在该态的精确本征值与它的限制的Hartree-Fock极限期望值之差。,相关能的定义,无限(完备

3、)基组变分计算,有限基组变分计算,相关能理论值,各种能量值的关系,相关能在化学问题中的重要性,理论和实验能量差值的比较 实验误差: 1kcal/mol0.002Hartree 体系总能量 1015Hartree 体系相关能 10-11Hartree 总能量差值受相关能影响很大。 弱相互作用（比如分子间作用）氢键键能310kcal/mol,一般处理方法,把 r12 包含在波函数中 适用于非常小的体系 有很多难计算的积分 He的Hylleras 波函数 用一组多电子波函数来表示体系的波函数 Hartree-Fock 行列式和激发态行列式 有非常多的激发态行列式, 收敛很慢 组态相互作用 (CI),

4、相关方法的目标,可以明确地给出(well defined) 不需要人为选择就可以应用于所有分子 可用于构建模型化学 效率 不能仅限于非常小的体系 变分性 是精确能量的上限 大小一致性(size extensive) E(A+B) = E(A) + E(B) 需要对热化学进行正确描写 计算成本与精确性的等级性 这样计算才能被系统改进,组态相互作用方法,HF方法把电子按从低到高填充, 是合理的, 但是, 有一定的人为性 从统计力学的角度来看, 任何填充方式都是可能的, 只不过其在真实体系中的几率不同, HF的填充方式是几率最大的一种 因为把电子都限制在了HF的占据轨道上, 所以使得不同自旋的电子过

5、分接近, 导致排斥能升高 为了改进这种硬性限制, 允许电子在所有可能的轨道上自由填充, 再确定出每种填充方式的几率。,组态相互作用,参考行列式(Hartree-Fock波函数),单激发态行列式(把占据轨道fi用未占据的轨道fa代替),双激发态行列式,等等,用变分原理来确定CI系数CIS 包括所有单激发态 用于激发态计算, 不是用于基态相关能 CISD 包括所有单激发和双激发态 对计算基态的相关能非常有用 O2V2 个行列式 (O= 占据轨道的数目, V=未占据轨道的数目) CISDT 单重, 双重和三重激发态 仅限于小分子, 大概 O3V3 个行列式 完全 CI 所有可能的激发态 (O+V)!

6、/O!V!)2 个行列式 在基组给定的情况下, 得到精确的相关能 大概限于14个轨道上分布14个电子,组态相互作用,组态相互作用,这是很大的求本征值的问题, 可以用迭代法求解 CI是行列式的线性组合 得到的是精确能量的上限(变分性) 可应用于激发态 比起非变分方法, 梯度可以简单计算出来,对距离很远的两个体系, 有E(A+B) =E(A)+E(B) 对有2个基函数的He原子, CID为对两个He原子分别计算得到 E=E(A)+E(B), = AB两个He原子距离很远时, CID为(没有包括四重激发态),大小一致性,多组态自洽场方法( MCSCF ),CI中采用的轨道是自洽场的分子轨道, (占据

7、轨道以及未占据的轨道), 只对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

8、):把所有自洽场分子轨道分为内层(永远双占据), 外层(永远不占据)和活性轨道, 只对活性轨道进行MCSCF计算, 这样, 同时减小了K和N。,CAS-SCF,只有 价轨道,只有 HOMO+LUMO,把Hamilton量分为精确可解的部分和微扰部分按l展开Hamilton量, 能量和波函数把l幂次相同的项合并在一起,微扰方法,选择 H0,其本征函数就是分子轨道波函数按照Hartree-Fock行列式, 单激发, 双激发和多重激发行列式的方式展开微扰波函数对能量的微扰校正为,Mller-Plesset 微扰理论,HF方法是MBPT方法的一级近似,HF方程的Hamilton量就是MBPT的零级表示

9、HF方程中的交换作用, 就是MBPT的一级表示,在每个阶次都是大小一致的 MP2 二阶计算相对容易(只需要双激发态) 在Hartree-Fock是一个相当好的起点时, 二阶将可以得到很大一部分相关能 实际使用中可以到四阶(单, 双, 三和四重激发态) MP4 能够得到大部分的相关能 整个系列有点振荡(偶数阶更低一些) 如果出现了严重的自旋污染或者Hartree-Fock不是一个好的起点, 相关能的收敛很困难,Mller-Plesset 微扰理论,HF,MP2,MP3,MP4,MP5,MP6,MBPT方法的说明,MBPT可以有MP2, MP3, MP4等各级修正, 主要应用的是MP2, MP4,

10、 而MP3给出的修正往往很差 MP2是从头算最经常使用的可承受包含相关能的方法,最常用的组合就是MP2/6-31G*, 但是在基函数超过500后就很慢了 MBPT方法有可能给出不体系能量低的值, 因为它没有变分原理保证其高于体系真实能量,CISD 可以写作T1 和 T2 生成了相应的系数组合的所有可能的单激发态和双激发态耦合簇理论的波函数,耦合簇理论,CCSD 能量和波幅可以通过解下式得到,耦合簇理论,CC 方程在波幅上是二次的 必须用迭代法求解 每次迭代的计算成本类似于 MP4SDQ CC 是大小一致的 CCSD 单和双激发态校正, 直到无穷大阶次 CCSD(T) 把三重激发用微扰的方式加上

11、 QCISD 组态相互作用, 且添加足够的二次项以使得其大小一致(等价于CCSD的截断形式),耦合簇理论,总结,CI是最直接的方法 MP是最常用的方法 耦合簇方法(CC)是最有效的计算相关能的方法,分子模拟与计算化学,模型化学 Model Chemistries,模型图(Pople图),Basis,模型化学,模型化学理论是一个完整的算法, 用来计算任何分子体系的能量 它在应用过程中没有任何主观因素 它必须是大小一致的, 这样每个分子的能量才能被唯一确定 简单的模型化学使用单一的理论方法和基组 复合的模型化学把几个理论模型和基组联合起来, 在低的计算成本下得到更精确的结果 模型化学的作用在于用我

12、们容许的计算量对某些种类的分子得到最精确的结果,确定方法的可信度,实验数据为确定某个理论方法的可信度提供了条件 G2和G3方法已经可以达到1kcal/mol的精度(化学精度) G2包括了原子化能, 电离能, 电子亲和能和质子亲和能等热化学数据, 体系中包含1-2个非氢原子 这些数据集合可以用来建立衡量理论水平,怎样把我们的方法用于未知体系,如果我们研究的体系没有实验数据, 就需要把理论方法逐级改进, 直到达到可信的结果 从模型图的左上角开始, 逐级向右下角改进, 进行几种方法的计算 当计算结果没有什么改进了, 我们的结果就可信了 注意: 这种方法不能用于DFT理论, 虽然它在基组方面可以改进,

13、确定模型的理论方法,几何结构很容易改进到可信范围, 而能量比较困难(因为它取决于相关能的改进) 因此, 一般在低理论水平优化结构而在高理论水平计算能量值 采用的方法是 能量计算方法/基组/优化方法/基组 例如: CCSD(T)/6-311G(2d,p)/HF/6-31G,模型化学比较(键长),108个主族分子在6-31G(d,p)下计算得到的键长与实验值的平均偏差,模型化学比较(原子化能),108个主族分子在6-31G(d,p)下计算得到的原子化能与实验值的平均偏差,模型化学比较(电子亲和能),F的电子亲和能计算值(eV),分子模拟与计算化学,密度泛函理论 Density Functional

14、 Theory,密度泛函理论,从头算理论从Schrodinger方程的波函数出发, 对原子和分子体系进行了近似描写 但是, 波函数并不是可观测量, 只有电子密度是可观测的, 用电子密度来描写体系当然是可能的 随着电子数目的增加, 波函数增大, HF以及MPn, CCSD(T), CI等的计算量也急剧增加, 但是电子密度r(x,y,z)并不增加, 仅仅是空间范围有所扩大,动机,波函数本身在本质上是无法解释的 减小体系的尺寸: N个电子的体系的波函数包含了4N个坐标 对于大体系, 基于波函数的方法很快就无能为力了, 虽然计算能力在持续提高, 但是电子之间的运动是相关的 希望使用某些物理可观测量,

15、而不是几率振幅,电子密度与波函数的关系,电子密度可以用波函数描写, 电子密度是三维的, 波函数是3n维的, 因此对多电子体系, 用电子密度比波函数要简便 电子密度必然失去了波函数的一些信息,电子能量的组成,总电子能量可以划分为E = ET + ENE +EJ + EX +ECET = 电子动能 ENE = 电子与核的库仑吸引能 EJ = 电子之间的库仑排斥能 EX = 交换能, 对电子自排斥能的校正 EC = 不同自旋的电子运动的相关能 ET, ENE, & EJ 包含了 E 的大部分 EX EC,能量是电子密度的函数 核-电子吸引和电子-电子排斥动能的Thomas-Fermi 近似交换能的S

16、later近似,Thomas-Fermi-Dirac (TFD) 模型,X 模型,TFD 不能预测成键, 总能量的误差在15-50 %. 如果Slater的Ex中的被看作参数, 就能迭代更好的结果 在X 模型中, = 3/4. 称为Hartree-Fock-Slater方法 虽然X 模型被现代泛函所取代, 但是它在无机体系和初步计算中仍然被使用,Hohenberg 和 Kohn (1964) 能量是密度的泛函 E 泛函是普遍适用的, 于体系无关 精确的密度对应于 E的最小值 只能应用于基态,理论基础,泛函和函数,函数: 把一个数变成一个数 泛函: 把一个函数变成一个数, 例如, 积分是f(x)的泛函: Ff(x)= f(x)dx,Kohn-Sham方法的基本思想,Kohn 和 Sham (1965) 局域泛函的变分方程 可以写作单Slater行列式, 但是其轨道与Hartree-Fock不同, 其轨道仅仅相对于Hartree的轨道, 没有包含电子的交换作用! EXC 包括了电子相关和交换作用 近似逼近的思想: 把分子能量表示成一些项的和, 在这些项中, 只有一项, 此项是未知的却是相对小的 自洽场的思想: 初始猜测一个电子密度, 从它来逼近真实的电子密度, 从最后得到的电子密度再计算体系的性质,