第一性原理计算ppt课件.ppt

4 3闭核系统的能量 在分子模型中 关心分子的基态能量他们大部分具有闭壳系统N 2个轨道 N个电子的闭壳系统每一个空间轨道 i具有两个自旋轨道 i i 1 Hiicore 电子在裸露原子核作用场中运动的能量 对于在轨道 i运动的电子 如果一个轨道具有两个电子 那么对N 2个电子而言单个电子能量为2Hiicore 2 电子与电子的作用项 i和 j中存在4个电子 在一个轨道上的2个电子以库仑作用的方式与另一个轨道上的两个电子发生作用 记为4Jij 但是在这些组合中 包含了成对电子的两种作用方式 给出了总的交换能 Kij 最后 在同一轨道上的一对电子库仑作用也要考虑在内 但是这些电子具有成对的自旋因此不存在交换能 3 jii Kii 4 5hartree Fock方程 变分方法真实函数近似计算的能量总是高于真实函数得到的能量 波函数越好 则能量越低 等能量最低时 得到最好的波函数 在最小点时 能量的一次微分 E 0 5 这一条件作用于能量表达式 得到Hartree Fock方程正交归一化条件可以写为交叠矩阵Sij能量最小的限制条件可以通过Lagrange因子解决最小化的函数的微分加上限制条件乘以Lagrange乘法因子的微分 其和设置为0在每一个正交归一条件下的Lagrange乘法因子记为 ij 6 Lagrange乘法因子被写为 ij 多电子系统和单电子系统的最大区别在于存在电子之间的作用 包括库仑和交换积分找到多电子系统最好的波函数 保留轨道图像找到一个解 当一个电子的轨道变化时由于电子耦合作用会影响其他自旋轨道上的电子的运动 在上述情况下 这个解能保证所考虑的电子的运动 7 假设一个电子位于原子核作用场中一个自旋轨道 I中 其他电子位于轨道 j中 芯部哈密顿量 库仑作用算符 交换算符 8 Fock算符 多电子系统有效的哈密顿算符对于闭核系统 9 除非i等于j 否则Lagrange乘法因子为0 每一个电子都被假设在包括原子核和其他电子的固定区域中运动每一个电子方程得到的解都会影响系统中其他电子的解自洽场假设解的思路得到一个H F本征方程的试探解 用来计算库仑和交还能 解H F方 给出第二套解 一次类推 SCF方法逐渐得到对应越来越低能量的单电子的解 直到一点 在这一点所有的电子的方程不再改变 10 5 1原子的H F计算以及Slater规则 如果假设电子的分布实球形对称的近似分析函数 氢原子得到的径向函数不能直接用于多电子原子 因为内壳层电子对原子核电荷具有屏蔽作用 如果考虑屏蔽效应把轨道指数作适当地修正 仍可以采用氢原子的波函数形式 11 Slater 12 完整的轨道 我们要用角度部分乘以径向函数部分 13 Slater给出了一系列的轨道指数的确定原则 Z为原子数 为屏蔽常数 n 为有效的主量子数 n n 1 2 3时取与n相同的值 n 4 5 6s时 取3 7 4 0 4 2 屏蔽常数按下列规则确定 把轨道划分如下组 1s 2s 2p 3s 3p 3d 4s 4p 4f 5s 5p 5d 14 对给定的轨道 屏蔽常数有下列贡献之和确定如果到原子核的距离比到上述各组远 取0 同一组的其他电子贡献为0 35 但若其他轨道为1s 则贡献为0 3 电子的主量子数比现有轨道小2或更多 则贡献为1 如果每一个电子具有的主量子数比现有轨道小1 现有轨道为d或f 取1 0 现有轨道为s p取0 85 15 硅的价电子的屏蔽常数计算电子结构为 1s2 2s22p6 3s23p2 规则b得到3 0 35 根据规则 c 得2 0 根据规则 d 得8 0 85 我们得到原子核数为14时 Z 为4 15 16 5 2H F方程中的原子轨道线性组合的方法 LCAO 直接求解H F方程 在实际中是行不通的 每一个自旋轨道写为单个电子轨道的线性组合 一个电子的波函数 v称为基函数 常对应原子轨道 将用希腊符号 来标记 共有K个基函数 我们希望得到K个分子轨道 17 根据变分思想 我们得到对应最低能量的系数Cvi通过改变系数得到波函数对于给定的基函数以及波函数的形式 最好的系数应满足 能量最小原理 18 5 3闭壳系统和Roothann Hall方程 具有N 2轨道N个电子的闭壳系统系统的H F方程最早由Roothaan和Hall独立提出所谓的Roothaan或Roothann Hall方程矩阵形式 它可以用标准的方法求解 19 引入了轨道 i的原子轨道展开形式 把相应的 i轨道带入H F方程 20 每一侧都乘以 u 1 也是一基函数 积分 两个基函数之间的交叠积分 写为Suv两个基函数的交叠项并不一定为0 Fock矩阵 21 闭壳系统的Fock矩阵元素 Fock矩阵地元素可以写为心部 库仑作用 和交换能项的综合 Huvcore每一个元素可以通过对应单个电子哈密顿项的动能积分和势能积分得到 库仑以及交换作用项 22 库仑算符Jj 1 由自旋轨道 j的交互作用引起 利用轨道函数的基函数线性组合的方法 交换项可以写为 23 通过引入密度矩阵 计算某一点的电子密度 24 闭壳系统Fock矩阵中每一项元素为 25 5 4解Roothann Hall方程 Fock矩阵微K乘k阶矩阵 如果运用实函数 则是对称的Roothaan Hall方程可以写为矩阵形式 FC SCE E为由轨道能量组成的对角矩阵 26 标准本征方法需要矩阵形式FC CE R H方程中只有交叠矩阵S取单位矩阵才为标准形式 交叠矩阵中非对角线元素可能不为0 为了用标准方法解R H方程 必须对之进行变换 这相当于变换基函数 使之正交归一化 寻找一个矩阵X 使XTSX I 对称正交归一化条件下交叠矩阵是对角化的 27 D为包含S本征值 得对角矩阵 U则包含了S的本征矢量 X UD 1 2UT X写为S 1 2 S 1 2SS 1 2 I R H方程两边预乘以S 1 2 S 1 2FC S 1 2SCE S1 2CE S 1 2S1 2加入单位矩阵 S 1 2F S 1 2S1 2 C S1 2CE 28 S 1 2FS 1 2 S1 2C S1 2C E F S 1 2FS 1 2 C S1 2C 得到F C C E 用标准化的方法解决 方程有解的条件为 行列式展开为多项式求解 对角化F 系数矩阵 C 是F 的本征矢量 基函数系数c可以由C 得到 C S 1 2C 29 求解R H方程的一般的方法 1计算积分以形成Fock矩阵F 2计算交叠积分S3对角化S4形成S 1 25猜想 和计算初始密度集镇P 6利用积分和密度矩阵形成Fock矩阵7生成F S 1 2FS 1 28通过对角化F 解 给出能量本征值和本征矢量C 9计算分子轨道系数C C S 1 2C 10从C中计算新的密度矩阵11检查计算结果 如果复合要求 停止 否则重复步骤6 30 需要猜想密度矩阵P 最简单