2012年10月24日.ppt

上传人:ga****84 文档编号:455707 上传时间:2018-10-09 格式:PPT 页数:6 大小:260.50KB
下载 相关 举报
2012年10月24日.ppt_第1页
第1页 / 共6页
2012年10月24日.ppt_第2页
第2页 / 共6页
2012年10月24日.ppt_第3页
第3页 / 共6页
2012年10月24日.ppt_第4页
第4页 / 共6页
2012年10月24日.ppt_第5页
第5页 / 共6页
点击查看更多>>
资源描述

1、电子结构方法复习,方法目的:给定一组原子构型,计算电子轨道波函数,从而得到系统能量、原子间 作用力及其它物理、化学性质。 对应量子方程:非相对论不含时情况下,为多体薛定谔方程。其中包含所有原子核 及电子自由度。,数值近似求解:因为多体薛定谔方程无解析精确解,需要在很多近似的基础上进行 数值求解。 Born-Oppenheimer 近似:对原子核及电子自由度分离变量,固定原子核的位 置,求解电子波函数。,Slater 行列式:多体电子波函数的薛定谔方程依旧无法数值求解。为了能够求解, 进一步分离电子轨道波函数的相关性,忽略电子波函数中因为库仑 排斥力引起的关联,只考虑电子作为费米子所具有的交换反

2、对称和 不可区分性,写出 Slater 行列式作为多体电子轨道波函数的解的函 数形式,是对真实多体波函数的近似。,变分原理:任意假设的波函数所对应的能量都大于等于真实的波函数解所对应的能 量,并且等号仅当波函数为真实解时成立。利用变分原理,可以先任意 选取一个电子轨道波函数的初始解,然后迭代调整这个解,使系统能量 逐渐变小,直到找到能量最优解。注意这个能量最小化的过程是非物理 的,只是一个数值求解的过程。 Fock 算符:利用拉格朗日未定乘子法,得到 Fock 算符,因为Fock 算符与轨道波函数有关,需要迭代求解:,一旦解出了 Fock 算符及对应的轨道能量,就可以得到系统总能量,基组:因为

3、每个分子轨道的具体函数形式 还是未知的,所以还要用一组给定 函数形式的基函数 (basis functions) 组成的基组 (basis set) 来展开分子轨道。,理论上任意形式的基组(基函数要满足正交归一化条件)都可以。实际求解时需要 权衡精度和计算量的要求适当选取。常用的基组: STO-3G:计算量小,精度一般。 3-21G,6-31G,6-311G, 6-311+G (3df,3pd):计算量适中,精度还可以。 cc-pVDZ, cc-pVTZ, cc-pVQZ: 精度高,计算量大。,赝势:用经典的解析函数代替内部电子自由度,只有价电子用量子方法计算。可以 节省大量计算时间,并可以包

4、含相对论效应。,电子关联方法:HF 引入 Slater 行列式,忽略了电子间的库仑排斥力对电子波函数 的贡献。电子关联方法试图基于 HF 计算结果的基础上,对电子轨 道波函数进行修正,以加入库仑排斥力的贡献。,激发态的 Slater 行列式:为了基于 HF 计算结果的基础上,对电子轨道波函数进行 修正,需要在 RHF 的 Slater 行列式的基础上,加入 Slater 行列式的激发态作为对电子波函数的修正项。如果 有一个电子调整到更高能量的轨道上,称为 S-type 激发 态,以此类推,有 D-type,T-type,Q-type,5-type, 6-type, 注意 Slater 行列式的

5、激发态只是一种数学上的 形式,没有直接对应的物理意义。,Configuration Interaction (CI) 方法: 在电子轨道波函数的表达式中加入激发态的 Slater 行列式,再用变分原理和拉格朗 日未定乘子法求出对应的系数。,Mller-Plesset (MP) 方法: 基于量子力学的微扰论。MP1 就是 HF。常用的是 MP2 和 MP4。,Coupled Cluster (CC) 方法: 用激发算符的指数算符使对电子波函数的修正一次包括无穷多级。,计算精度大致为:HF MP2 CISD MP4(SDQ) CCSD MP4 CCSD(T),关键在于根据计算问题的要求,权衡计算精度和计算量,适当选取计算方法和基组。,

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 学术论文资料库 > 毕业论文

Copyright © 2018-2021 Wenke99.com All rights reserved

工信部备案号浙ICP备20026746号-2  

公安局备案号:浙公网安备33038302330469号

本站为C2C交文档易平台,即用户上传的文档直接卖给下载用户,本站只是网络服务中间平台,所有原创文档下载所得归上传人所有,若您发现上传作品侵犯了您的权利,请立刻联系网站客服并提供证据,平台将在3个工作日内予以改正。