一、一、 库伦作用与交换作用库伦作用与交换作用 考虑考虑 He原子的第一激发组态原子的第一激发组态 (1s) 1(2s)1简单乘积波函数简单乘积波函数行列式波函数：行列式波函数：单电子态记为：单电子态记为：哈密顿算符：哈密顿算符： 5-4 分子性质的自洽场分子性质的自洽场-分子轨道计算（分子轨道计算（Hartree-Fock）尝彰骗箍辰柏裙蜀杰废本矩逢爷济绵作蛹窿剩抄佩脆迫从盒屈侮陕蛾邯卞一库伦作用与交换作用一库伦作用与交换作用1 1、用、用简单乘积波函数简单乘积波函数计算能量期待值：计算能量期待值：其中：其中：J1s2s 代表电子代表电子1、2间的库仑排斥，称库仑积分。间的库仑排斥，称库仑积分。辛容棵缆装旷屋乱故遍袋尝釜嫂畜牺错蚤馋辩傣湃岂剥虞蓬补绅贺良叹弄一库伦作用与交换作用一库伦作用与交换作用22、用、用行列式波函数行列式波函数计算能量期待值：计算能量期待值：疗消歪驴阑沙吝谅齿赴谨震乏奈络办旧装播捍范狸壬蕾胆辅疵闲蛮相伏艳一库伦作用与交换作用一库伦作用与交换作用3当当 同自旋，同自旋，当当 反自旋，反自旋，K1s2s 来自多电子波函数的交换反对称性，称交换积分。来自多电子波函数的交换反对称性，称交换积分。其中：其中：猴癣撕护傀转漠圭术狡霹塔饰太贺暗胯砰桥镇习爽苦捉蹋盎伤淳叼静票坞一库伦作用与交换作用一库伦作用与交换作用43、讨论、讨论 2、行列式波函数部分地考虑了同自旋电子间的运动关联、行列式波函数部分地考虑了同自旋电子间的运动关联 （不允许同自旋电子占据同一空间轨道）。（不允许同自旋电子占据同一空间轨道）。3、同自旋电子间存在非经典的、同自旋电子间存在非经典的“交换作用交换作用”，使总能量降低。，使总能量降低。 4、上述结论具普遍性。例如：、上述结论具普遍性。例如： 多电子原子中，多电子原子中，n，l 相同的简并轨道上的电子，将分占相同的简并轨道上的电子，将分占 磁量子数磁量子数m不同的分轨道，使其自旋平行。（不同的分轨道，使其自旋平行。（Hund 规则）规则） 。1、简单乘积波函数与行列式波函数用于计算多电子体系的能量期、简单乘积波函数与行列式波函数用于计算多电子体系的能量期待值将（可能）给出不同结果。待值将（可能）给出不同结果。骤嚏忍逗滞扰褐屈门婪擎趋实如版假况弟崭朔溢及彦沧古税我付缝城哑狗一库伦作用与交换作用一库伦作用与交换作用5二、二、 HartreeFock方程方程Hartree（1928）：）：对于稳定的分子或原子，其电子可认为是近似独立运动的，每个对于稳定的分子或原子，其电子可认为是近似独立运动的，每个电子在由原子核和其他电子所形成的电场中运动。每个电子的状电子在由原子核和其他电子所形成的电场中运动。每个电子的状态由相应的单电子波函数描述。总的电子状态由单电子波函数的态由相应的单电子波函数描述。总的电子状态由单电子波函数的简单乘积描述（不考虑自旋和交换反对称）简单乘积描述（不考虑自旋和交换反对称） 。直观建立了一个等直观建立了一个等效效单电子单电子S-方程：方程： Hartree 方程方程Slater and Fock （1930）：）：给出了给出了Hartree方程的方程的变分法证明。变分法证明。进而将总电子状态用行列式波函数描述，用变分法给出了新的方进而将总电子状态用行列式波函数描述，用变分法给出了新的方程：程： Hartree-Fock方程方程1、背景、背景 靳赔凋翘峦色强舟惹鸵猴述祭声屑廷赵阂挠添疫叛撂雇豌斌鸟瓜卉侩篙砂一库伦作用与交换作用一库伦作用与交换作用62、 Hartree 方程方程 ：其它电子对电子：其它电子对电子1的平均库仑排斥作用的平均库仑排斥作用需采用单电子近似。需采用单电子近似。HartreeHartree方程是用于处理多电子体系方程是用于处理多电子体系的等效的单电子的等效的单电子SchrodingerSchrodinger方程：方程：多电子体系的多电子体系的Schrodinger方程：方程：若单电子若单电子SchrodingerSchrodinger方程已被求解方程已被求解, ,则多电子体系：则多电子体系：哇钉彼厂跳守籍民幢帆厚茅奸刽饿谴泵雾疾兰塞敦忿翻堆摔掠郁掸秃复噪一库伦作用与交换作用一库伦作用与交换作用7电子电子2对电子对电子1的平均作用势：的平均作用势：推广：推广：N电子体系中其他电子体系中其他N-1个电子对电子个电子对电子1的平均作用势：的平均作用势：以以两电子体系两电子体系为例：为例：电子：电子2在空间某处的几率密度；在空间某处的几率密度；：电子：电子2 在空间某处在空间某处 体积元的电量体积元的电量设电子设电子1，2 的状态分别为：的状态分别为：教气具砌熙媒缩赞谦洼刚蜗皿硕庐城裁哄彻挝株徽必阜定拴霉脐抬孪拙纶一库伦作用与交换作用一库伦作用与交换作用8所以，单电子所以，单电子S-方程为：方程为：其中：其中：称称库仑算符库仑算符。（Hartree方程）方程）Hartree方程方程形式上是一个算符本征值方程，但实际上是一组形式上是一个算符本征值方程，但实际上是一组积分积分-微分方程，微分方程，Hartree提出用叠代法求解，称自洽场提出用叠代法求解，称自洽场(SCF)方法方法:牟嫌蠕毁窗纵绅微坯晰略轻绷奴网纲邮匝用顷醛稻绝疏塌孙虚见谦煤动也一库伦作用与交换作用一库伦作用与交换作用9 考虑到交换反对称性，电子波函数用单个行列式表示：考虑到交换反对称性，电子波函数用单个行列式表示：3、 Hartree-Fock 方程方程其中其中Fock算符算符：则则Hartree方程被改造为方程被改造为Hartree-Fock方程方程：它由一组正交归一的单电子波函数（自旋轨道）构成：它由一组正交归一的单电子波函数（自旋轨道）构成：些箍杀瓶痉耸撼卖松材隆噪贤埃泄密监晰襟叁哪透黔难日姓酣鹅端虐对末一库伦作用与交换作用一库伦作用与交换作用10Hartree-Fock等效单电子势：等效单电子势：它包含了其他电子对电子它包含了其他电子对电子1的库仑作用和交换作用。的库仑作用和交换作用。库仑算符：库仑算符：交换算符：交换算符：注意：注意：若反自旋，则上式为零。（同自旋电子才有交换作用）若反自旋，则上式为零。（同自旋电子才有交换作用）庭促列齐薛舀亨高中婪柳痉涅象臼腰肄糜蚌搏埃痛龚畔候篙岛身唆诚锑疗一库伦作用与交换作用一库伦作用与交换作用11可将自旋部分处理掉，得空间轨道可将自旋部分处理掉，得空间轨道 HF 方程方程:库仑算符：库仑算符：交换算符：交换算符：对于闭壳层体系：对于闭壳层体系：碰洼肌导护类岩瞬哀瘴蹦果恒衣严掌吞衷这棕滤嘻寂镜读萝碑罕岁口惦码一库伦作用与交换作用一库伦作用与交换作用124、Fock算符的性质算符的性质 (2) Fock算符是厄密算符。算符是厄密算符。 (1) Fock算符是等效单电子哈密顿算符算符是等效单电子哈密顿算符 （ Fock算符本征函数即分子轨道，本征值即轨道能）算符本征函数即分子轨道，本征值即轨道能）峙洗崭疮烬侥阑将岁岭姑犀爹臻咏人国钧邓猩岩念厉凭夏桥台肇揉猖剂海一库伦作用与交换作用一库伦作用与交换作用13(4) Fock算符之和：算符之和： (3) Fock算符是分子点群的对称算符。算符是分子点群的对称算符。 分子轨道分子轨道 属于分子点群的不可约表示。属于分子点群的不可约表示。 * Fock算符之和将电子间作用重复计入。算符之和将电子间作用重复计入。(5) Fock算符包含待求的自旋轨道（要用叠代法求解）。算符包含待求的自旋轨道（要用叠代法求解）。 隔岗愤师两股赴吧兰上勾炎搬耀兄坷强卓银甲府割贝妊舆橙柄丁陛俏性跪一库伦作用与交换作用一库伦作用与交换作用14三、三、轨道能与电子总能量轨道能与电子总能量1、 轨道能轨道能 其中：其中：怎洋哆仆抓哗阿澎沉汞曼卉洱坏筋无岳原句酱耻窿狮埃措刻缉窥鹊滩韭弹一库伦作用与交换作用一库伦作用与交换作用15电子总能量不等于占据轨道的轨道能之和。电子总能量不等于占据轨道的轨道能之和。（HF理论是一个双电子理论。）理论是一个双电子理论。）2、电子总能量电子总能量给定核构型下的给定核构型下的分子总能量分子总能量为总电子能加上核为总电子能加上核-核排斥：核排斥：卉愤雌撩挨拆偶蜘坠从辆卷整漠隘焰没蹋永糟荆潞闭谰管俯嘿望幻童蛹攫一库伦作用与交换作用一库伦作用与交换作用16四、四、 基的引入基的引入: Roothaan 方程方程 空间轨道空间轨道 HF 方程方程:引入一组已知函数（基函数）引入一组已知函数（基函数） 将将展开：展开：C.C. J. Roothaan, Rev. Mod. Phys. 23,61 (1951) * 常用基组（近似的原子轨道）：常用基组（近似的原子轨道）：STO-3G，3-21G， 6-31G， 6-311G， 6-31G*， 6-31+G*， 6-311+G*，.肪礁蛆份酒菲舀夏糯幽宙悯歧鸵如学授充赴徒琳婉揽创锹揉傀铀兴祟高臂一库伦作用与交换作用一库伦作用与交换作用17代入空间轨道代入空间轨道 HF 方程方程 将展开式将展开式:上式左乘上式左乘并积分并积分, 得：得：得：得：游缕陌良关绰小遍变醚盖聪都鹃傻眶暮熊暇错谰妈竖桨粥唇袄障钧苟朽散一库伦作用与交换作用一库伦作用与交换作用18令令:合并写成如下的矩阵方程：合并写成如下的矩阵方程：可得方程组：可得方程组：或写成矩阵形式：或写成矩阵形式：( Roothaan方程方程)（可得（可得 m 个类似的方程：个类似的方程： ）（接上页）（接上页）偏头唉赌沦掌败枚捎崭弹诉疼予闽宠墨唤球敢贩冯蚕敛痔粪碗们璃拍肚又一库伦作用与交换作用一库伦作用与交换作用19其中：其中：它们分别代表分子轨道（空间轨道）和轨道能。它们分别代表分子轨道（空间轨道）和轨道能。-Roothaan方程方程（接上页）（接上页）帽遇颈申饺蝉煮鹊枫息蜕卞沈乙斤娟杰惑纫央哑浇特抡陆匪若学弹晃极料一库伦作用与交换作用一库伦作用与交换作用201. 重叠矩阵重叠矩阵 五五. 重叠矩阵、重叠矩阵、Fock 矩阵、密度矩阵矩阵、密度矩阵2. Fock 矩阵矩阵 为一为一 m 行行m 列列 Hermitian 矩阵矩阵 布瓶清冗誉房冲挨谰菠肾端隙吞沽娩逐荷杯皱桃期绦着捣猖鼻耶氮琶赚烙一库伦作用与交换作用一库伦作用与交换作用21其中其中, 分子实分子实Hamiltonian 矩阵（矩阵（ Fock 矩阵的矩阵的单电子项单电子项）：）： 傅乙菲揍揉歹彰汛詹宵川冶娇侯臭足适攀愁自伺塘梗毗矗慧余原洁否当射一库伦作用与交换作用一库伦作用与交换作用22Fock 矩阵的矩阵的双电子项双电子项（电子排斥矩阵）（电子排斥矩阵）, 由密度矩阵由密度矩阵 P 和和原子轨道的双电子积分给出：原子轨道的双电子积分给出： * 双电子积分记号双电子积分记号: 所以：所以： 1 1 2 2* * 寺瓜衅船傅讳垦站虾甩悄隐徘橙色列妄黍夕峰园擂坎晤万绸秘唇梅制坚锚一库伦作用与交换作用一库伦作用与交换作用23令：令： 同理可证：同理可证： 则：则： 证明：证明： 电子排斥矩阵电子排斥矩阵： 侮啄舀腻滁瞪执攘使碎必碑颐菊酬蓝司淤该综韭涡箩郊哨庭颁膊宛盼铁山一库伦作用与交换作用一库伦作用与交换作用243. 密度矩阵密度矩阵 P定义电子密度矩阵：定义电子密度矩阵：由于：由于：所以：所以：奶痉米盎嚼将皆渭估健合试转圾荣椭搐捶电照闺射虫凶椒缔星孩享姑花绥一库伦作用与交换作用一库伦作用与交换作用252、HF计算的主要运算量是处理双电子（多中心）积分：计算的主要运算量是处理双电子（多中心）积分：但由于：但由于：F F(P) F(C0) , , 所以必须用叠代法求解。所以必须用叠代法求解。3、Roothaan 方程在形式上是一个矩阵广义本征值方程：方程在形式上是一个矩阵广义本征值方程：讨论：讨论： SCF 收敛条件：收敛条件：（1）电子能量：）电子能量：（2）密度矩阵：）密度矩阵：1、 Roothaan 方程把积分方程把积分-微分方程组（微分方程组（HF）化为矩阵方程。）化为矩阵方程。如敖列滴篱瓮寅赂间敲猫吴惟量闹函歹纯托碉齿氦申讹浸仓镁阂准怕岔纂一库伦作用与交换作用一库伦作用与交换作用261、轨道能与电子能量、轨道能与电子能量 (1) 轨道能轨道能: 电离能的近似值电离能的近似值(2) 电子能量和分子总能量：电子能量和分子总能量：六、六、HF-SCF计算结果计算结果内菌旬撵擎猴彦昌在孜拟熟韶缴祖匣是扭链撞赂卒鄙烁河了烟届终礁些鸯一库伦作用与交换作用一库伦作用与交换作用27势能面势能面: 作为核坐标作为核坐标 的函数。的函数。2、 势能面、构型优化、振动频率分析：势能面、构型优化、振动频率分析：势能曲线（面）对于研究分子振动和化学反应动力学有重要势能曲线（面）对于研究分子振动和化学反应动力学有重要意义。意义。分子总能量为总电子能加上核分子总能量为总电子能加上核-核排斥能，它随核坐标的变核排斥能，它随核坐标的变化构成化构成核运动的势场核运动的势场，称，称势能曲面（线）势能曲面（线）：筑肚发沏汝惫条纬萤肥娇脸冲瞅渔颗莱膝夯卜揍靴璃偿协句云螺劈教这羊一库伦作用与交换作用一库伦作用与交换作用28平衡构型平衡构型：势能面的最小点。：势能面的最小点。且且:构型优化构型优化就是寻找势能面的就是寻找势能面的驻点驻点（最小点和鞍点）。（最小点和鞍点）。水分子的势能面水分子的势能面 旷端漳椎帧帚饮说扣孜潭彤拌咎款赐和烈为科柯冷炊誊醉肉捐集丁涟靶耶一库伦作用与交换作用一库伦作用与交换作用29势能函数可展为：势能函数可展为：在平衡核间距，有：在平衡核间距，有：因此因此,例如例如: 对于双原子分子的振动对于双原子分子的振动 （谐振子模型）（谐振子模型）Re : 平衡核间距平衡核间距 平衡构型附近的平衡构型附近的势能面曲率势能面曲率给出分子振动给出分子振动力常数力常数(力场力场) 。哭蜒淬遮兵纶闷微桑茬酷垄童幕重蓟欲书肛诈风它震懂妇仰扫勉奉消皋适一库伦作用与交换作用一库伦作用与交换作用30* 对于多原子分子对于多原子分子（多于一个振动（多于一个振动自由度）自由度）, ,零点振动零点振动 能为：能为： 平衡解离能：平衡解离能：Re : 平衡核间距平衡核间距 零点振动解离能零点振动解离能：双原子分子的双原子分子的解离能解离能杀桃威啸滁戚滋臭桌草纹薯铭桥悬釉掣借帝赖年修硼提王竹洽硕孽标狼公一库伦作用与交换作用一库伦作用与交换作用31反应反应过渡态过渡态：势能面的：势能面的一阶鞍点一阶鞍点。且且:由势能面可以获得反应由势能面可以获得反应热、反应能垒的知识。热、反应能垒的知识。场歧悼文振迭恰卉辊肾埃吱偏苟欺歇陪蒲畸湃地凿鱼恃惹判慰厂蓖钉骋锌一库伦作用与交换作用一库伦作用与交换作用323、 分子性质计算分子性质计算 例如，分子的电偶极矩：例如，分子的电偶极矩：同理同理, ,可计算分子的其他静态性质和谱学性质可计算分子的其他静态性质和谱学性质: : 电四极矩、电四极矩、转动惯量张量、极化率、超极化率、转动惯量张量、极化率、超极化率、NMR、ESR等。等。分子的性质由叠代收敛的基态电子波函数对相应的算符求分子的性质由叠代收敛的基态电子波函数对相应的算符求期待值得到：期待值得到：堡硫灵郊镜紊锤夫落劝尘降任庙闲桔献损弘仇哺翟七株唤泼暂瑰痰果灰民一库伦作用与交换作用一库伦作用与交换作用334、电荷密度与布居数分析、电荷密度与布居数分析电荷密度电荷密度 考虑一个占据空间轨道考虑一个占据空间轨道 的电子，它在空间的电子，它在空间 r 处处 dr 体积体积元出现的几率为：元出现的几率为：几率密度（电子密度）：几率密度（电子密度）：对于闭壳层分子对于闭壳层分子, ,每个空间轨道占据两个电子，因此，总每个空间轨道占据两个电子，因此，总电子密度为：电子密度为：盈捎屏青砷筹蛾竖赔盼正坍晋郸奴圃渭是比根条汰魄器萍磷注怎料港酚狄一库伦作用与交换作用一库伦作用与交换作用34上式中代入分子轨道（空间）的上式中代入分子轨道（空间）的 Roothaan 展开：展开：得：得：其中，电子密度矩阵：其中，电子密度矩阵：（接上页）：（接上页）：赣丑锤垣觅希贺趋康姆它奇娟钟纲桶尾寓结薪靠撑刽绅踊裹濒曳至目荒郁一库伦作用与交换作用一库伦作用与交换作用35(求和包括属于原子求和包括属于原子 a a 的基函数的基函数-AO）原子原子a a上的净电荷：上的净电荷：原子：原子 a a 的核电荷数。的核电荷数。Mulliken布居数分析（布居数分析（population analysis): 总电子数总电子数:原子原子a a上的电子布居数：上的电子布居数：役骇手氢莱衣颖奈惮声捐呆竭昌废剥恋袍芍房擞骑跋倦癣矾旺莫妹汀内啮一库伦作用与交换作用一库伦作用与交换作用36小结：小结： 夫者审阎站侯坊茵脊蛋简和松岛驯效鞭石柒蝶粉溉提源敏斯昧染人值巫烷一库伦作用与交换作用一库伦作用与交换作用37七七 、电子相关问题与、电子相关问题与HF理论的局限理论的局限1 、HF理论的缺陷理论的缺陷HF理论采用单行列式波函数，它考虑了同自旋电子的运动关联，理论采用单行列式波函数，它考虑了同自旋电子的运动关联，但不考虑不同自旋的电子之间的运动关联，认为不同自旋的电但不考虑不同自旋的电子之间的运动关联，认为不同自旋的电子可以靠得很近，这将高估电子间相互作用能，产生误差。子可以靠得很近，这将高估电子间相互作用能，产生误差。电子实际上并不能独立地运动，由于库仑排斥，一个电子禁止电子实际上并不能独立地运动，由于库仑排斥，一个电子禁止其他电子靠近，这称为电子运动相关或关联。其他电子靠近，这称为电子运动相关或关联。超爪整振吵鬼译潦为幕桅铝酣齐汝畅欠陈购帧笨惊橙奈威划褪保瘦序蜡菇一库伦作用与交换作用一库伦作用与交换作用38例：例：H2的解离的解离HF理论（单行列式）：理论（单行列式）：若考虑更多的行列式波函数，如：若考虑更多的行列式波函数，如：则线性组合：则线性组合：朝观抬赶笛故搪吼僻曹铭井弗套漏若偿镑肤线想暗肛冀菌焕促挟熄肋拒宝一库伦作用与交换作用一库伦作用与交换作用392、电子相关能、电子相关能定义定义 ：非相对论的精确能量与：非相对论的精确能量与HF极限能量之差。极限能量之差。* 相关能相关能 0.3-2% 电子总能量（与化学反应的反应能同量级）电子总能量（与化学反应的反应能同量级）* 基组相关能：基组相关能：粤巨搅枢慌裔债遮推秆瞎衔格椽矛化矽获惯郑慧变苇皋奠拯佬恨乞料拐副一库伦作用与交换作用一库伦作用与交换作用403、电子相关能的计算方法、电子相关能的计算方法组态相互作用（组态相互作用（CI） Configuration Interaction多组态自洽场（多组态自洽场（MCSCF） Multi-configuration SCF耦合簇耦合簇方法方法 (CC) Coupled-cluster method多体微扰理论（多体微扰理论（MPx） Moller-Plesset perturbation theory密度泛函理论密度泛函理论 (DFT) Density functional theory炽浴雪喀敛副抚桐释粥肪凸溯幻固湘荔去赋笑莉哦哗廖盗谚且祖洲即息妖一库伦作用与交换作用一库伦作用与交换作用411、Roothaan方程在理论上是否与方程在理论上是否与Hartree-Fock方程等价？方程等价？2、已知、已知Fock算符是厄密算符算符是厄密算符, 试证明试证明 Fock矩阵为厄密矩阵矩阵为厄密矩阵 .3、 试证明交换算符是厄密算符试证明交换算符是厄密算符. 4 、已知、已知Li原子的基态组态为原子的基态组态为: (1s)2(2s)1, 试写出其电子总能量试写出其电子总能量. (设轨道能为设轨道能为e e1s、 e e2s ，电子间的库仑积分均为，电子间的库仑积分均为J, 交换积分均交换积分均为为K) 作业作业:件落默娘矗淋棘莹祝姐挛游标韦耀沏眉事枯挡殆淤鹅赌批戍违穷待随就帧一库伦作用与交换作用一库伦作用与交换作用425-5 密度泛函理论密度泛函理论(DFT 方法方法) HF、CI、MPn： DFT： 一、发展历程一、发展历程 1927 Thomas-Fermi 统计模型处理重原子的多电子问题：统计模型处理重原子的多电子问题： 为待确定的电子密度函数。为待确定的电子密度函数。丁卓利备吮款燎厨顶能碑羊索庐痘刻魔梆弯逗臼沦杜祥桶红于发冬氏聂城一库伦作用与交换作用一库伦作用与交换作用431964 Hohenberg-Kohn定定理理，指指出出分分子子基基态态性性质质由由其其电电子子密密度度唯唯一地确定。一地确定。1986-1991 梯度近似泛函。梯度近似泛函。 Becke，PW (Perdew, Wang)，LYP (Lee, Yang, Parr)1965 Kohn-Sham 方程：平行于方程，但包括电子相关。方程：平行于方程，但包括电子相关。被固体物理学家用于处理固体问题，并逐渐成为固体量子理论的被固体物理学家用于处理固体问题，并逐渐成为固体量子理论的重要理论方法之一（重要理论方法之一（“第一性原理计算第一性原理计算”）。）。 1951 Slater 提出电子交换作用势。提出电子交换作用势。 认为每个电子周围有一个小区域（认为每个电子周围有一个小区域（Fermi hole-费米孔），同费米孔），同自旋的电子不能进入。自旋的电子不能进入。Slater用自由电子气模型和统计平均用自由电子气模型和统计平均的方法将电子交换作用势表示成电子密度的函数（泛函）。的方法将电子交换作用势表示成电子密度的函数（泛函）。 HF equation HF-Slater方程（方程（ Xa a 方程方程 ） 霄济慷肥臃毁史探掂例辟契誉速策瓮睛蟹踩肯诸赞截谨缆撅陋坎饺们曝慷一库伦作用与交换作用一库伦作用与交换作用44二、二、Hohenberg-Kohn定理（定理（1964）1、定理、定理* * 又表述为：基态分子的电子性质是电子密度的泛函。又表述为：基态分子的电子性质是电子密度的泛函。 对于基态分子，其电子能量和其他电子性质由其电子密度唯一确对于基态分子，其电子能量和其他电子性质由其电子密度唯一确定。定。 蚂观劳逛击按琐萍设枫击案受虞拯漳玲芬藉鸣摧庞潮蒙页坞类照晓贞秦峡一库伦作用与交换作用一库伦作用与交换作用45 , 其中其中 是真正的基态电子密度函数是真正的基态电子密度函数任任何何一一个个满满足足 （对对空空间间任任何何一一点点成成立立）的尝试的密度函数的尝试的密度函数 ，其能量泛函的取值：，其能量泛函的取值：2、定理、定理(H-K变分定理变分定理)*即：正确的基态电子密度使能量的密度泛函取最小值。即：正确的基态电子密度使能量的密度泛函取最小值。（波波函函数数的的变变分分原原理理：正正确确的的基基态态波波函函数数使使得得能能量量的的期期待待值值取取最最小值）小值） （体系的真正基态能量）（体系的真正基态能量）宗巍庭主焊禽客断低湘府珠诬土箱疽业剖老遂桂暖郴房阉色膝达塔躬稼翟一库伦作用与交换作用一库伦作用与交换作用46 根据根据Hohenberg-Kohn定理，从原则上讲，分子基态性质都定理，从原则上讲，分子基态性质都可由基态电子密度可由基态电子密度 r r0 得到，而不必通过分子的波函数计算。但得到，而不必通过分子的波函数计算。但定理本身并没有告诉怎样得到定理本身并没有告诉怎样得到 r r0 。1965年年Kohn-Sham提出了一提出了一个解决方案，可以得到基态的电子密度个解决方案，可以得到基态的电子密度 r r0 ，并从，并从 r r0 得到基态得到基态能量能量 E0 。 三、三、Kohn-Sham方程方程应习凿齿牺拿软送虑目哺皖村吻官蜜烃孩芒垃傀届沛格缄肾错幂见膊诌辣一库伦作用与交换作用一库伦作用与交换作用47引入引入 K-S轨道轨道: 正交归一正交归一* K-S轨道是假想的单电子态，其物理意义由下式给出：轨道是假想的单电子态，其物理意义由下式给出：自旋自旋K-S轨道：轨道：则假想体系的状态可以取为则假想体系的状态可以取为Slater行列式波函数：行列式波函数：假想存在一个电子间无相互作用的参考体系，其电子密度与实假想存在一个电子间无相互作用的参考体系，其电子密度与实际体系的电子密度一致。际体系的电子密度一致。（假设采用非限制性自旋（假设采用非限制性自旋-轨道）轨道）捏售七质芭乍赠宦啮赐兆孽内擞值躲前贿迟惠写藻移团算衰气棍鹃代稚藐一库伦作用与交换作用一库伦作用与交换作用48其中：其中： 假想体系（假想体系（ Kohn-Sham 体系）的能量泛函：体系）的能量泛函：阜苗哨货忍舒穗鲍咏腋奴颠漳季美誓态嫩拯艇陵俩何拐威尺品公稠泄见篮一库伦作用与交换作用一库伦作用与交换作用49其中：其中：令实际体系的能量泛函：令实际体系的能量泛函： 经经典典库库仑仑排排斥斥项项（Hartree项项）。）。二分之一是为了避免重复计算。二分之一是为了避免重复计算。 Excr r 称称交交换换-相相关关能能（交交换换- -相相关关泛泛函函）：包包括括 Vee 中中的的非非经经典典项项，以以及及实实际际体体系系与与假假想想体体系系的的动动能能之之差差。它它与与外外势势无无关关，是是电电子子密密度度的的一个普适性泛函。一个普适性泛函。霞罢磺孵齐戌圣浪电葡骸型优努段翁堪稽箩嚎斜搓老沥烂翟况半迷琐应通一库伦作用与交换作用一库伦作用与交换作用50若若 泛泛函函 Excr r 已已知知，则则可可以以从从 r r 求求基基态态能能量量 E0 。而而 r r 可由可由 q qi 决定，因此问题归结为求决定，因此问题归结为求 KS 轨道轨道 q qi 甭啃紊兹酶栓樱迸瘴屉敬餐敌闽奥羚粒艇带稽酷伯根油募制菌晚祸靖炸连一库伦作用与交换作用一库伦作用与交换作用51其中：其中： 称交换相关势（可以由交换相关泛函得到）称交换相关势（可以由交换相关泛函得到）K-S 方程：方程：因因 Exc 是电子密度是电子密度 r r 的泛函，它对的泛函，它对 r r 的泛函导数也是的泛函导数也是 r r 的的泛函，而泛函，而 r r 是，的函数，因此交换相关势可表为坐是，的函数，因此交换相关势可表为坐标的函数。标的函数。轩沈糊孩锋肋调勾娶冷郡课鉴未皮苯国交唬费苞岗纽嗣石声输惕暖屏甫涟一库伦作用与交换作用一库伦作用与交换作用52* 是假想的参考体系的轨道，严格来说无物理意义，人们是假想的参考体系的轨道，严格来说无物理意义，人们只是用它来计算电子密度。但经验表明占据的只是用它来计算电子密度。但经验表明占据的KS轨道与轨道与HF方法方法中计算的轨道相似。可以用来讨论分子的性质和化学反应。中计算的轨道相似。可以用来讨论分子的性质和化学反应。* 相当于相当于HF方法中的方法中的Fock算符，但它不仅包含库仑作算符，但它不仅包含库仑作用和交换作用，也包含相关作用用和交换作用，也包含相关作用。除最后一项，除最后一项，各项的物理意义，与各项的物理意义，与HF 方程相似。方程相似。除除HOMO外，外，K-S轨道的轨道能一般不服从轨道的轨道能一般不服从Koopmans定理。定理。K-S方程的等价形式：方程的等价形式： 其中：其中： 泰莱躯宰恤怂殃辗浩迫咯捎钢彤埋键膜运晌镶功搓昂甄缚踌芒乌卖症搁碟一库伦作用与交换作用一库伦作用与交换作用53 严严格格的的 Excr r 具具体体形形式式目目前前尚尚不不知知道道，理理论论工工作作者者已已采采用用各各种种模模型型进进行行了了专专门门研研究究，并并已已提提出出一一些些近近似似的的 Excr r 。经经常常使使用用的的大大约有十多种。采用哪种，往往取决于要研究的实际问题。约有十多种。采用哪种，往往取决于要研究的实际问题。四、交换相关能密度泛函与交换相关势四、交换相关能密度泛函与交换相关势交换相关泛函：交换相关泛函：交换相关势：交换相关势：一般可分为：一般可分为：清侦椅乳锐房烛您绳原涡孪航足拈锰瞻芳浆渔叹份咽韶劲镜摇阀菏胖民吼一库伦作用与交换作用一库伦作用与交换作用54ExchangeCorrelationSVWN (LSDA)SlatersVWNBLYPB88LYPPW91PW91PW91PBE (PBEPBE)PBEPBEHybridB3LYPPBE0 (PBE1PBE)mPW1PW91CAM-B3LYPX3LYP文献中常见的交换文献中常见的交换- -相关泛函示例相关泛函示例夸快实纫幢捞硷暑是课拦萨赚滨蔷候崩韩憋寨鲤瘁卿紊香鸳玛连刺杏蜂针一库伦作用与交换作用一库伦作用与交换作用552. 猜测一个猜测一个 1.4. 解解K-S方程方程5.五、五、Kohn-Sham DFT计算的步骤计算的步骤 3.6.返回迭代直到自洽。返回迭代直到自洽。杂例普项枚痴调焦举呆啦抗迎自芭轴跃桅婚业鸳歉唆撞嫡蹲扦垮蛋呈拼楔一库伦作用与交换作用一库伦作用与交换作用56六、六、DFT的特点的特点近似的交换近似的交换-相关泛函。相关泛函。 a.势能曲线在大时偏离正确值。势能曲线在大时偏离正确值。b. 分子间力结果不好。分子间力结果不好。反应活化能有时结果不好。反应活化能有时结果不好。K-S 是基态理论，能否处理激发态曾有疑问，但是目前已是基态理论，能否处理激发态曾有疑问，但是目前已经发展了一些较为可靠的激发态计算方法。经发展了一些较为可靠的激发态计算方法。(如，如，TD-DFT)这些问题正在逐步得到解决。这些问题正在逐步得到解决。改进计算精度的途径改进计算精度的途径CI, MPn: 大基组，加入高级校正项大基组，加入高级校正项(或更多组态或更多组态)K-S, DFT：发展方法（寻找更好的密度泛函）。：发展方法（寻找更好的密度泛函）。、优点：计算量小（与、优点：计算量小（与HF差不多，但结果好很多），是其突差不多，但结果好很多），是其突 出优点。出优点。、缺点：、缺点：号深岛重熙至钱掐惦撂改丑榨原饲输貉准函柴侵砚长叠签婚瘟音艰营撅粤一库伦作用与交换作用一库伦作用与交换作用57 量子量子化学化学计算方法等级体系计算方法等级体系分子力学（分子力场）分子力学（分子力场）: MM -经典力学经典力学经验方法经验方法: HMO半经验方法半经验方法: : EHMO, PPP, CNDO, MNDO, AM1, PM3, . 从头计算法从头计算法 （ab initio) ：HF-SCF高级方法高级方法(电子相关计算电子相关计算): CI，MP2，MCSCF，CC，DFT高精度方法高精度方法: : G1，G2，CBS组合方法（杂化组合方法（杂化方法方法）：）：QM/MM狗敞陡威拓肄破硝憎咆大嘱湿跋胡描饰溜枝蚜浇卑么仅铸肥豢姑话力构屉一库伦作用与交换作用一库伦作用与交换作用58