《计算化学及其应用05MO计算中的基组.ppt》由会员分享,可在线阅读,更多相关《计算化学及其应用05MO计算中的基组.ppt(48页珍藏版)》请在金锄头文库上搜索。
1、计算化学及其应用MO计算中的基组Basis Sets for Molecular Orbital CalculationsLCAO 近似Hartree-Fock轨道的数值解法只能对原子和双原子分子使用双原子分子的轨道类似与原子轨道的线性组合比如, H2中的Sigma键 1sA + 1sB对多原子分子, 其轨道必须用原子轨道的线性组合(LCAO)来近似表示基函数称为基函数其中心通常选在原子上当然, 也可以使用比原子轨道更具有一般性, 更灵活的基函数当基函数足够多足够好时, 得到的分子轨道就能相当准确基函数的一般形式径向部分角度部分基组的角度部分5D/6D与7F/10F5D/6D: x2, y2,
2、 z2 (x2+y2+z2)s轨道xy, yz, zx7F/10F: x3, xy2, xz2 x(x2+y2+z2)y3, y2z, y2x y(x2+y2+z2)z3, z2x, z2y z(x2+y2+z2)xyzp轨道径向部分: Slater型函数(STO)对氢原子, 它是精确的通常用于原子计算有正确的渐进性质(r)在核处满足尖点条件(r0)三中心和四中心双电子积分没有可用的公式STO与类氢离子的径向函数1s, 2p, 3d, 4f的指前因子为单项式, 它们之间相互正交为了保证其它的s轨道(2s, 3s,.)与1s正交, 其它的s轨道的指前因子都是多项式, 3p, 4p, ., 4d,
3、 5d, ., 5f, 6f, . 也是如此多项式在积分时*d会分解成N2项, 计算比较麻烦, 于是Slater引入了Slater型原子轨道, 它相当于把类氢离子径向波函数正交化限制去掉, 而要求其指前因子为单项式! (为分子计算提供条件)类氢离子径向分布函数径向部分: Gaussian型函数(GTF)Gaussian型函数的多中心积分GTO乘积定理: 双中心函数乘积=单中心函数STO和GTF的比较在r比较大时, 衰减太快在核处没有尖点优点: 所有的双电子积分都可以用公式直接计算收缩Gaussian基函数用GTF的线性组合来构造一个更适合进行分子计算的基函数一个GTF称为一个素基函数(prim
4、itive basis function)基组基组就是给一些原子确定出的一组指数和收缩系数基组的类型最小基组双zeta, 三zeta基组等阶层分裂基组极化函数弥散函数最小基组只有中性原子的轨道壳层才给定基函数比如对C元素有1s, 2s, 2px, 2py, 2pzSTO-3G3 GFT拟合一个STO STO的指数从原子计算得到, 用一些典型的分子来调整这也称为单zeta基组 (zeta, , 是Slater型轨道的指数)双Zeta基组 (DZ)最小基组中的每个基函数都变成两个一组比较紧(离原子核近, 指数大), 另一组比较松(离原子核较元素, 指数比较小)因此在描写电子云方面, 径向上有一定的
5、伸缩性如果原子在分子中失去一部分电子, 其密度就会收缩如果原子得到一些电子, 其密度就会膨胀双Zeta基组使得在径向上有一定伸缩性Zeta1 Zeta2C1 + C2远核 近核C1=1, C2=0C1=0, C2=1C11, C21+价层分裂基组只有价层的基组加倍(基函数越少, 计算越快, 计算成本越小)内层轨道用最小基组, 因为它们原子化合成分子的过程中几乎不变3-21G (1s用3个GTF, 近核的2s, 2p用2个GTF, 远离核的2s, 2p用1个GTF)6-31G (1s用6个GTF, 近核的2s, 2p用3个GTF, 远离核的2s, 2p用1个GTF)极化函数角动量更大的函数添加到
6、基组中时, 轨道在角度方向就有了更大的可变性比如对氢原子加上p函数, 对C原子加上d函数如果不加极化函数, 即使使用大基组, Hartree-Fock计算得到的NH3的结构是平面型的没有极化函数时, 环丙烷的张力太大6-31G(d) (也写作 6-31G*) 对非氢原子添加d函数6-31G(d,p) (也写作6-31G*) 对氢原子添加p函数且对非氢原子添加d函数DZP 带极化函数的DZEffect of Polarization Functions弥散函数把指数很小的函数添加到基组中需要的情形: 阴离子, 电负性高的原子F, Cl, 计算电子亲和能和气相酸度, 弱作用(氢键, 范德华作用)6
7、-31+G 对非氢原子添加弥散s和p函数6-31+G 对氢原子添加弥散s函数且对非氢原子添加弥散s和p函数相关一致基函数一系列逐渐增大的基组,考虑了电子相关效应其径向和角度两方面的近似是相互一致的, 可用它外推到基组极限cc-pVDZ DZ 且对非氢原子添加d函数, 对氢原子添加p函数cc-pVTZ 价层三重分裂函数, 且对非氢原子添加2个d函数和1个f函数, 对氢原子添加2个p函数和1个d函数cc-pVQZ, cc-pV5Z, cc-pV6Z还可以再添加弥散函数 (aug-cc-pVXZ)相关一致基组(第二周期)(sp)组极化函数cc-pVDZ3s2p1dcc-pVTZ4s3p2d1fcc-
8、pVQZ5s4p3d2f1g相关一致基组的精确度单Zeta双Zeta三Zeta极化函数Zeta数目1d2d1f3d2f1g6-31G*(中等基组)半经验基组0级1级2级1级2级3级近似等级赝势, 有效势基组在化学作用中, 内层轨道没有大的变化价层轨道同时受到核和内层电子的作用可以构造一个赝势, 来代替来自核和内层电子对价层电子的静电势可以用更少的基函数来描写原子(但是却引入了其它近似)对于重元素, 赝势可以把难以处理的相对论效应给包括进去, 改善了计算结果计算化学及其应用单电子性质One Electron Properties分子轨道的概念波函数的整体性: 体系的电子是整体分布在整个空间的波函
9、数的可分性: 电子在整个能量区间或空间是有相对分布的按照能量划分, 得到不同的分子轨道按分子轨道能量划分: 按照空间划分, 得到定域的分子轨道不同能量的分子轨道(离域)分子轨道是离域的, 分布在整个分子上正交归一的可以在图形界面下画出三维的分子轨道图, 其中|i(r)| = ci(r) 有正负值正负号只有相对的意义, 没有绝对的意义定域分子轨道定域分子轨道形成一个正交归一函数集合, 占据与非占据之间也正交仅仅分布在1到2个原子上, 至多3个原子上可以与化学上的电子分布的直观概念联系起来: Lewis结构, 杂化, 成键(, p键)可以用于改进算法, 因为它具有相当的可加性计算量不再随N4增加,
10、 而是相当线性地增加离域分子轨道和定域分子轨道可以通过一个酉变换联系起来NBO中的自然定域分子轨道乙烷的CH和*CH原子电荷在化学家中, 电子的转移和共享是最直接的化学成键图像电子转移对应于原子电荷, 共享对应于于共价键的键级在分子轨道理论中, 没有原子的概念, 每个电子是分布于整个分子体系中因此确定“分子中的原子”就是把分子轨道理论与化学直观联系起来的桥梁Mulliken集居数分析把总的电子密度分割到各个原子或基函数上每个轨道都是归一的 对于闭壳层分子, 每个占据轨道上有2个电子 (Ne 是电子总数)Mulliken集居数分析每个轨道是基函数的线性组合把它代入电子总数的表达式并且重新组合各项
11、Mulliken集居数分析密度矩阵, Smn 重叠矩阵按照基函数所属的原子把 分配到不同原子上Mulliken集居数分析矩阵集中在各个原子上原子电荷A-B间电子数: Mulliken键级各个原子上电子数Mulliken的PS矩阵Wiberg和和Mayer键级键级l共价键的键级可以定量地反映键的强弱lMulliken键级与我们的经验键级没有定量的关系lWiberg键级及其推广Mayer键级(PS)(PS)2lWiberg键级与经验键级基本一致, lMayer键级与经验键级完全一致l例如:H-H (1.0), C-H (0.99), C=C (2.01)l对不同基组有一定的稳定性l应用:近似估计原
12、子之间共价键成份lGaussian03的NBO中可以输出Wiberg键级分子中的原子轨道(AOIM)从分子轨道理论得到原子电荷和键级的关键的确定“分子中的原子”!分子中的原子类似与自由原子, 因此其原子轨道应该是最小原子轨道, 而不是计算必须使用的双Zeta或更高的基组相对于自由原子, 分子中的原子的原子轨道要收缩(失去电子时)或膨胀(得到电子时), 这由分子波函数来确定, 并不是事先指定的这样得到的原子轨道称为分子中的原子轨道, 把分子波函数用这种轨道投影出来, 再进行Mulliken集居数分析AOIM对Mulliken分析的改进自然键轨道自然键轨道(NBO)分析分析http:/www.ch
13、em.wisc.edu/nbo5/l使用自然正交的轨道(一阶约化密度的本征轨道)进行轨道分析l自然轨道与原子轨道只是定性相似l因为正交, 所以轨道之间没有重叠, 与化学键性质不符l所得到的电荷比较稳定, 与基组的关系不大l电荷, 键型, 杂化方向, 共振成份等2s和2p自然原子轨道(乙烷, B3LYP/6-311+G*)AIM(分子中的原子)方法电子密度与波函数不同, 是一个可观测的真实量体系中电子密度在空间中的变化是有规律的电子密度函数在空间中的变化可以用其对空间坐标(x,y,z)的一阶和二阶微商来描写根据电子密度函数的微商把电子划分为不同的区域, 归属于相应区域的核电子密度(乙烯)电子密度
14、的梯度按照电子密度梯度的零线划分电子密度给每个原子来自静电势的原子电荷(CHelpG)CHELPG (= CHarges from ELectrostatic Potentials using a Grid based method) 原子电荷要拟合出分子周围一些点处的静电势, 首先计算分子的静电势3pm取点按照分子大小再扩大28pm抛弃在范德华半径内的点对其它所有点进行拟合优化偶极矩对于 Hartree-Fock 波函数, 偶极就是偶极经典表示的期望值它可以用密度矩阵和基函数的一组积分来表示电子密度在rC处, 一个单位电荷感受到的势能静电势电子总能量的表示Kopmanns定理和Brilliouin定理Kopmanns定理: 本征值ek近似对应于第k个轨道的电离能的负值Brilliouin定理: Y0是体系的Slater行列式, Y0a是用一个虚轨道代替一个占据轨道后的行列式, 那么: =0。用HF行列式函数Y0为零级近似基态波函数计算出的单电子算符的期望值精确到一级。(在Moller-Plesset展开下)电负性光谱电负性: Allen, 1989电负性是原子价层电子的平均能量对主族元素: 可以用理论方法计算, 对于分子体系, 可以使用AOIM方法进行计算