电子管理与计算化学

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

2、孔 hf方法基本上反映了fermi孔(由于波函数的反对称性), 完全没有反映columb孔!正六角形h6分子的相关函数(a)r2在氢核上(b)r2在氢核中间hf的相关能是精确的(exact)! 由于hf方法主要未考虑columb相关, 忽略了电子之间的排斥作用(排斥使得能量升高), 因此, 相关能是负值, 将导致体系计算能量更加负 lowdin的定义: 指定的一个hamilton量的某个本征态的电子相关能, 是指该hamilton量在该态的精确本征值与它的限制的hartree-fock极限期望值之差。相关能的定义无限(完备)基组变分计算有限基组变分计算近似hf极限值hf极限值相关能理论值总能量

3、实验值精确非相对论极限值相对论效应校正基组误差近似相关能实验值基组误差相对论效应校正(近似)近似相关能计算值近似非相对论极限值各种能量值的关系相关能在化学问题中的重要性 理论和实验能量差值的比较 实验误差: 1kcal/mol0.002hartree 体系总能量 1015hartree 体系相关能 10-11hartree 总能量差值受相关能影响很大。 弱相互作用（比如分子间作用） 氢键键能310kcal/mol一般处理方法 把 r12 包含在波函数中 适用于非常小的体系 有很多难计算的积分 he的hylleras 波函数 用一组多电子波函数来表示体系的波函数 hartree-fock 行列式

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

5、的一种 因为把电子都限制在了hf的占据轨道上, 所以使得不同自旋的电子过分接近, 导致排斥能升高 为了改进这种硬性限制, 允许电子在所有可能的轨道上自由填充, 再确定出每种填充方式的几率。组态相互作用00111111111|(,)aaabababcabciiijijijkijkiaijabijkabcnaiiainabijiaijbjniajbttt 参考行列式(hartree-fock波函数)单激发态行列式(把占据轨道i用未占据的轨道a代替)双激发态行列式等等用变分原理来确定ci系数cis 包括所有单激发态 用于激发态计算, 不是用于基态相关能cisd 包括所有单激发和双激发态 对计算基态的

6、相关能非常有用 o2v2 个行列式 (o= 占据轨道的数目, v=未占据轨道的数目)cisdt 单重, 双重和三重激发态 仅限于小分子, 大概 o3v3 个行列式完全 ci 所有可能的激发态 (o+v)!/o!v!)2 个行列式 在基组给定的情况下, 得到精确的相关能 大概限于14个轨道上分布14个电子组态相互作用tddetttabcijkijkabcabcijkabijijababijaiiaai respect towith /minimize*0h组态相互作用 这是很大的求本征值的问题, 可以用迭代法求解 ci是行列式的线性组合 得到的是精确能量的上限(变分性) 可应用于激发态 比起非变

7、分方法, 梯度可以简单计算出来tthe 对距离很远的两个体系, 有e(a+b) =e(a)+e(b) 对有2个基函数的he原子, cid为 对两个he原子分别计算得到 e=e(a)+e(b), = ab 两个he原子距离很远时, cid为(没有包括四重激发态)大小一致性aaaaat22110bbaabbaabbbbaaaabbbbbaaaaabababatttttt222211112211221102211022110)(bbbbaaaabatt221122110多组态自洽场方法( mcscf ) ci中采用的轨道是自洽场的分子轨道, (占据轨道以及未占据的轨道), 只对c0, c1等进行了变

8、分, 并没有保证使用的分子轨道是最优的 多组态自洽场方法, 不仅对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): 把所有自洽场分子轨道分为内层(永远双占据), 外

9、层(永远不占据)和活性轨道, 只对活性轨道进行mcscf计算, 这样, 同时减小了k和n。能量0cas-scf内层活性空间只有价轨道只有homo+lumo 把hamilton量分为精确可解的部分和微扰部分 按l展开hamilton量, 能量和波函数 把l幂次相同的项合并在一起微扰方法vhhhh00000iiiee221022100llllleeeevhhdeeeedeeee10202112012000101100100000vvhvvhh 选择 h0,其本征函数就是分子轨道波函数 按照hartree-fock行列式, 单激发, 双激发和多重激发行列式的方式展开微扰波函数 对能量的微扰校正为ml

10、ler-plesset 微扰理论ifh0bajijiaaabijhfhfhfmphfdedeeeeddeee,2010220000010vvvhabcijkijkabcabcijkabijijababijaiiaaiaaa1hf方法是mbpt方法的一级近似 hf方程的hamilton量就是mbpt的零级表示 hf方程中的交换作用, 就是mbpt的一级表示occupied(0)0011;niiiihfe(1)011111(0)(1)hf00111();( |)( | )2nnnnnjjij ijij iijijvjkeii jjij ijrreee 在每个阶次都是大小一致的 mp2 二阶计算相对

11、容易(只需要双激发态) 在hartree-fock是一个相当好的起点时, 二阶将可以得到很大一部分相关能 实际使用中可以到四阶(单, 双, 三和四重激发态) mp4 能够得到大部分的相关能 整个系列有点振荡(偶数阶更低一些) 如果出现了严重的自旋污染或者hartree-fock不是一个好的起点, 相关能的收敛很困难mller-plesset 微扰理论hfmp2mp3mp4mp5mp6mbpt方法的说明 mbpt可以有mp2, mp3, mp4等各级修正, 主要应用的是mp2, mp4, 而mp3给出的修正往往很差 mp2是从头算最经常使用的可承受包含相关能的方法,最常用的组合就是mp2/6-3

12、1g*, 但是在基函数超过500后就很慢了 mbpt方法有可能给出不体系能量低的值, 因为它没有变分原理保证其高于体系真实能量 cisd 可以写作 t1 和 t2 生成了相应的系数组合的所有可能的单激发态和双激发态 耦合簇理论的波函数耦合簇理论021)1 (ttcisdabijbajiabijt ,02t021)1exp(ttccsd ccsd 能量和波幅可以通过解下式得到耦合簇理论0)(0)(0)(0dededeccccabijccccaicccchhh cc 方程在波幅上是二次的 必须用迭代法求解 每次迭代的计算成本类似于 mp4sdq cc 是大小一致的 ccsd 单和双激发态校正, 直

13、到无穷大阶次 ccsd(t) 把三重激发用微扰的方式加上 qcisd 组态相互作用, 且添加足够的二次项以使得其大小一致(等价于ccsd的截断形式)耦合簇理论总结 ci是最直接的方法 mp是最常用的方法 耦合簇方法(cc)是最有效的计算相关能的方法分子模拟与计算化学分子模拟与计算化学模型化学model chemistries模型图(pople图)minimal sto-3g hf mp2 mp3 mp4 qcisd(t) . . . full ci split-valence 3-21g polarized 6-31g* 6-311g* diffuse 6-311+g* high ang. m

14、om. 6-311+g(3df,p) n2 n5 n6 n7 schrdinger equation basis模型化学 模型化学理论是一个完整的算法, 用来计算任何分子体系的能量 它在应用过程中没有任何主观因素 它必须是大小一致的, 这样每个分子的能量才能被唯一确定 简单的模型化学使用单一的理论方法和基组 复合的模型化学把几个理论模型和基组联合起来, 在低的计算成本下得到更精确的结果 模型化学的作用在于用我们容许的计算量对某些种类的分子得到最精确的结果确定方法的可信度 实验数据为确定某个理论方法的可信度提供了条件 g2和g3方法已经可以达到1kcal/mol的精度(化学精度) g2包括了原子

15、化能, 电离能, 电子亲和能和质子亲和能等热化学数据, 体系中包含1-2个非氢原子 这些数据集合可以用来建立衡量理论水平怎样把我们的方法用于未知体系 如果我们研究的体系没有实验数据, 就需要把理论方法逐级改进, 直到达到可信的结果 从模型图的左上角开始, 逐级向右下角改进, 进行几种方法的计算 当计算结果没有什么改进了, 我们的结果就可信了 注意: 这种方法不能用于dft理论, 虽然它在基组方面可以改进确定模型的理论方法 几何结构很容易改进到可信范围, 而能量比较困难(因为它取决于相关能的改进) 因此, 一般在低理论水平优化结构而在高理论水平计算能量值 采用的方法是 能量计算方法/基组/优化方

16、法/基组 例如: 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.2 22.052.27.06.8 108个主族分子在6-31g(d,p)下计算得到的原子化能与实验值的平均偏差模型化学比较(电子亲和能) f的电子亲和能计算值(ev)hfmp2b3lypsto-3g-10.16 -10.16-9.013-21g-1.

17、98-1.22-0.866-31g(d)-0.39+1.07+1.056-311+g(2df,p)+1.20+3.44+3.466-311+g(3df,2p)+1.19+3.54+3.46实验值+3.48分子模拟与计算化学分子模拟与计算化学密度泛函理论density functional theory密度泛函理论 从头算理论从schrodinger方程的波函数出发, 对原子和分子体系进行了近似描写 但是, 波函数并不是可观测量, 只有电子密度是可观测的, 用电子密度来描写体系当然是可能的 随着电子数目的增加, 波函数增大, hf以及mpn, ccsd(t), ci等的计算量也急剧增加, 但是电

18、子密度r(x,y,z)并不增加, 仅仅是空间范围有所扩大动机 波函数本身在本质上是无法解释的 减小体系的尺寸: n个电子的体系的波函数包含了4n个坐标 对于大体系, 基于波函数的方法很快就无能为力了, 虽然计算能力在持续提高, 但是电子之间的运动是相关的 希望使用某些物理可观测量, 而不是几率振幅电子密度与波函数的关系 电子密度可以用波函数描写, 电子密度是三维的, 波函数是3n维的, 因此对多电子体系, 用电子密度比波函数要简便 电子密度必然失去了波函数的一些信息21( , , )*( , , )( , , )|niiix y zx y zx y znr 电子能量的组成 总电子能量可以划分为

19、e = et + ene +ej + ex +ecet= 电子动能ene= 电子与核的库仑吸引能ej= 电子之间的库仑排斥能ex= 交换能, 对电子自排斥能的校正ec= 不同自旋的电子运动的相关能 et, ene, & ej 包含了 e 的大部分 ex ec 能量是电子密度r的函数 核-电子吸引和电子-电子排斥 动能的thomas-fermi 近似 交换能的slater近似3/223/5)3(103rrcdctexr cr4 /3dc 9a831/3a 1enerzaarradjrr(1)r(2)r12d1d2thomas-fermi-dirac (tfd) 模型xa 模型 tfd 不

20、能预测成键, 总能量的误差在15-50 %. 如果slater的ex中的a被看作参数, 就能迭代更好的结果 在xa 模型中, a = 3/4. 称为hartree-fock-slater方法 虽然xa 模型被现代泛函所取代, 但是它在无机体系和初步计算中仍然被使用 hohenberg 和 kohn (1964) 能量是密度的泛函 er 泛函是普遍适用的, 于体系无关 精确的密度对应于 er的最小值 只能应用于基态理论基础泛函和函数 函数: 把一个数变成一个数 泛函: 把一个函数变成一个数, 例如, 积分是f(x)的泛函: ff(x)= f(x)dx3320( )( )24( )302811(

21、)44f xxf xxf x dxxf xx kohn-sham方法的基本思想 kohn 和 sham (1965) 局域泛函的变分方程 r 可以写作单slater行列式, 但是其轨道与hartree-fock不同, 其轨道仅仅相对于hartree的轨道, 没有包含电子的交换作用! exc 包括了电子相关和交换作用 近似逼近的思想: 把分子能量表示成一些项的和, 在这些项中, 只有一项, 此项是未知的却是相对小的 自洽场的思想: 初始猜测一个电子密度, 从它来逼近真实的电子密度, 从最后得到的电子密度再计算体系的性质kohn-sham能量 电子在势场中的总能量(在假想的无相互作用的参考态rr(

22、x,y,z)下): 到现在为止没有引入任何近似, 但是也无法实际用上式进行计算, 因为动能和电子排斥能的形式未知00002201 nuclei a10000( )( ) ( )neeennaneiiiiaeeetvvzvv rretr v r drvrrrrrrr kohn-sham能量的分析 对未知形式的动能和电子排斥能, 把主要部分用参考态表示出来, 只剩下小部分是未知形式的!000022102101020120120102000121200()( )( )12( )( )1( ) ( )2xcrnriieeeereeettttdrrrvdrdrvrrrer v r drtdrdrrtvr

23、rrrrrrrrrrrrrr er tsr ener jr tr excr vnn交换相关泛函exc的意义1.exc包括了动能相关能和势能相关能和交换能, 2.每一项都是负值, 导致体系能量更低3.每一项的准确函数都是未知的, 只能用近似方法来逼近4.其中的交换能是大的(hf的一级近似), 相关能是小的(hf的二级及以上近似), ex ec!5.这里的交换相关能与hf方法中的不完全一样kohn-sham(ks)方程 同样要把多电子方程转换为单电子方程组, 类似于hf方法, 得到对应于hf方程的kohn-sham方程222122201( )1(1)(1)(1)2(1) ( )( )( )ksks

24、ksaixciiinucleiaianksriixcxczrdrvrrervrrr rrrr kohn-sham方程的解法 ks方程与hf方程完全类似, 因此其解法也与hf基本相同, 即给定一组基组, 猜测基组系数, 构成分子轨道, 计算电子密度, 得到电子势场, 再由此得到新的电子密度和电子能量; 由新电子密度再得到新的势场, 进行自洽场计算 不同之处: ks的fock算符比较复杂, 其矩阵元(特别是, 无法用解析方法求解, 因此一般用格点法或拟合函数法近似逼近其矩阵元的值交换相关泛函 交换相关泛函exc是dft理论中主要问题 在ks方程中出现的是交换相关势场vxc , 它是交换相关能泛函的

25、微商 只要知道交换相关能泛函exc, 就知道交换相关势场泛函vxc, 它们都是密度的函数; 在分子中, 电子密度在空间中各个点不同, 它们的值在空间各个点也不同 交换相关能有相当明确的物理意义和图像, 因此比较好进行近似估计 xcxcevrrr构造密度泛函 精确的形式无法知道, hohenberg-kohn 定理仅仅是证明其存在 密度泛函的一般形式为: 对 lsda: a = b = c = 0 对纯泛函: a = 0 有可能系统地改进泛函, 但是因为对其限制条件和性质尚未搞清楚, 需要更进一步的研究localhflocalxxxcc (1)(1)ea eb eaec ec er 局域自旋密度

26、泛函 (lsda) 只与密度有关 从电子气中推导出来 广义梯度近似 (gga) 与 |r|/r4/3有关 meta-gga 与 有关 杂化泛函 混进一些hartree-fock交换作用(精确的交换作用)密度泛函的类型(r) 12i(r)i2精确性与计算成本化学精度增加化学精度增加计算成本降低计算成本降低交换泛函, exlsdaggameta-ggaxa1951dirac1930g96b86b88pw91pbe1996rpbe1999revpbe1998xpbe2004pw86mpwtpss2003br89pkzb1999cs1975相关泛函, eclsdaggameta-ggaw38xpbe2

27、004pw86pbe1996pw91lyp1988b95tpss2003pkzb1999b88vwn1980pz81pw92ca data1980体系电子密度的分布及其描写内层电子gga也可能无法完全描写价层电子必须使用gga金属固体价层lda就可以了 交换-相关泛函必须用数值方法来积分 没有用公式积分那样强大 能量和梯度花费的时间是hartree-fock方法的1-3倍 频率计算是hf的2-4倍 有些计算成本可以弥补回来, 因为对纯密度泛函, 进行库仑相互作用计算时可以采用密度拟合方法近似计算exc 项的计算dft 的性能 要得到精确的能量, 需要一些精确的交换泛函, 虽然系统改进的泛函(比

28、如tpss)已经胜过杂化泛函了 在能量方面, gga 泛函已经对lsda做出了很大的改进, 从gga到meta-gga的改进没有从lsda到gga的幅度大 开壳层dft计算出现自旋污染的情况不多 多组态问题和范德华作用还很难用dft研究dft 术语 交换-相关孔密度, nxc(r,r) 描写的是在点r处出现有关电子时, 它减小其它电子的总密度在r处的分布 线性响应函数, (r,r;) 描写的是在r处出现有关微扰势, 在 r出的总密度随频率的变化 能量密度, x 或 c 每个粒子的能量密度dft 策略 象通常一样, 在用很大的弥散和极化基组前先用最小或dz基组计算 做初步计算或确定起始结构时, 用lsda而不要用hf, 这对包含金属的体系特别重要 用纯泛函优化结构可以使用密度拟合来加快速度 要注意 的值, 必要时要进行稳定性检验 如果需要用更精确的积分格子, 那么先用默认的格子优化, 再用更精确的格子进行进一步优化gaussian中的dft 交换泛函和相关泛函都必须给出. tpsstpss 不是所有的交换泛函和相关泛函都可以组合. lsda 泛函: lsda, svwn5, xalp