量子化学-HFR方程预案.pptVIP

  1. 1、本文档共38页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Hartree-Fock-Roothaan方程;Hartree方程 n电子原子体系的Hamiltonian算符为: 单电子算符: 双电子算符: 总电子Hamiltonian算符:;当体系只有一个电子时(即氢原子和类氢离子 ),Schr?dinger方程为: 该方程可以严格求解。 当体系有两个或两个以上电子时,Schr?dinger方程为: 由于双电子算符?(i,j)中含有两个电子的坐标,难以分离变量,故此Schr?dinger方程难以求解。 ;; 单电子近似的结果是将原来需要求解含有n个电子坐标的体系总波函数的问题拆分为求解n个单电子波函数的问题。后者由于自变量数目大幅度减少,所以在数学上的求解难度也相应地大为降低。 ;Hartree自洽场方法 在n电子原子中,电子i的势能算符为: ?j(j)是描述电子j的波函数,对电子j的所有空间可能的位置取平均可以得到电子i和电子j的平均排斥能: 由于上式右边是对电子j的坐标积分,所以 只与电子i坐标有关,这样其他 (n-1)个电子对电子i的排斥算符(即Coulomb算符)则可表示为:;单电子Hamiltonian算符可改写为: 电子i的单电子Schr?dinger方程为: 这就是著名的Hartree方程。 要求解?i(i),必须先知道?j(j),这相对于要解这个方程必须先知道方程的解。为解决这个困难,可采用迭代法进行求解。;先假定一组N个函数作为波函数的初值: 将这组初值代入Hartree方程,得到一组N个新的波函数作为一级近似波函数: 再用这组一级近似波函数代入Hartree方程,可得到二级近似波函数,重复这个过程,直到最后求解的单电子波函数、对应的轨道能量或总能量收敛为止。;缺陷;波函数的反对称性和Slater行列式 电子属于费米子(Fermions,自旋量子数为半整数的微观粒子)。由费米子构成全同粒子体系服从费米统计规律,体系总波函数必须满足两个必要条件:反对称性和全同性。 Hartree方程得到的波函数为: 该波函数不满足多电子波函数的反对称性条件,并且没有考虑电子自旋问题。; 波函数反对称性:任意交换两个粒子的全坐标后波函数的绝对值相同,符号相反。即: 全同性:任意交换两个粒子的全坐标后粒子的总几率分布不变。即:;在轨道近似条件下,多电子体系反对称波函数的最简单表达形式是由分子轨道波函数构成的单Slater行列式: 注:行列式的行号:标注轨道 列号:标注电子;用单Slater行列式表示的波函数能够满足费米子特性的要求: 反对称性:交换两个电子的全坐标,相当于交换行列式的两列,行列式的数值不变但正负符号相反; 全同性:每个电子在每个MO上出现的机会均等。 人们将采用Slater行列式(考虑了自旋波函数)后的Hartree方程称为Hartree-Fock (HF)方程。;解方程思路;Hartree-Fock算符和方程: 1.电子i的Hartree-Fock方程为: ;能量表达 1.先求解单电子Fock算符的期望值:;其中Kij为电子i和电子j之间的交换双电子积分。;记 Gi=;Hartree-Fock-Roothaan方程;(一)、Roothaan方程的简单推导; 该方程为广义本征方程。其中ε相当于算符F的本征值, C相当于算符F属于本征值ε的本征向量。;;定义密度矩阵元 ;总能获得酉阵A,使得I=A+SA;解HFR方程:;一个具体的例子—氨的自洽场计算:;输出文件;;;简单的单电子乘积; 轨道近似忽略了电子相关,所以Hartree方程两种相关都未考虑进去,而HF方程或HFR方程只考虑了自旋相关,并未考虑库仑相关。;Koopmans (库普曼)定理;对于定域共价键,键级=1/2(成键电子数-反键电子数)。例如,??2的键级为1,O2为2,N2为3。 对于离域π键,定义键级 ; 静电势 将单位正电荷从无穷远处移到分子周围空间某点所需做的机械功。;电荷密度;Hellmann-Feyman 定理

文档评论(0)

a336661148 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档