computational_chemistry-4.ppt_第1页
computational_chemistry-4.ppt_第2页
computational_chemistry-4.ppt_第3页
computational_chemistry-4.ppt_第4页
computational_chemistry-4.ppt_第5页
已阅读5页,还剩16页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、计算化学理论和应用 第四讲,2005 Computational Chemistry laboratory Beijing Normal university,计算化学理论和应用 第四讲,2005 Computational Chemistry laboratory Beijing Normal university,基组总结,MOLCAOAO(STO)STO-NGTO 收缩基: STO-NGTO 极小基: 使用一个收缩基来近似一条轨道 双z基: 使用两个收缩基来近似一条轨道 分裂价基: 内层轨道使用收缩基,价层使用多z基 极化基:比轨道更高角动量的基函数 弥散基:轨道指数较小,更弥散的基函数

2、,HFR方程的自洽场法求解,Fock矩阵 重叠矩阵 ei 组成的对角阵 系数列矩阵,需要将广义本征方程转化为标准方程,以方便求解,关键:将重叠矩阵转化为单位矩阵,因为重叠矩阵为厄米矩阵,因此总可以找到一个酉阵将它对角化,在此基础上,构造矩阵,可以得到,因此可以通过如下方式对广义本征方程进行转换,SCF 计算实例: HeH+,几何构型: He-H键长:R=1.4632a.u.; 核电荷:He 2, H 1; 电子总数:2 基组:STO-3G,一、指定分子几何构型和基组,轨道指数和系数的确定最小二乘法,通过调整Ci和ai可以使3个GTO与STO达到最佳拟合 刻度因子: 简化GTO指数的优化过程,优

3、化得到的He的1S和H的1S轨道,二、初猜,计算Fock矩阵矩阵元,四指标积分,共N4个,实际要计算的少一些,对于本例,N2,q6,最简单的初猜:,三、对角化S矩阵,矩阵对角化的方法,可以解得,四、求解得到新的密度矩阵,五、构造新的Fock矩阵,迭代,迭代过程,自洽场计算收敛问题及解决,两种原因造成的不收敛:数学上不收敛 化学上不合理,1,改变初猜,使用GUESS=CORE,HUCKEL,MIX关键词; 2,对于开壳层体系,先对于其相应的闭壳层离子计算到 收敛,然后使用GUESS=READ关键词从chk文件中读取 分子轨道初猜; 3,先用小基组得到分子轨道,然后同样使用GUESS=RE AD关键词从chk文件中读取; 4,使用能级提升法,将HOMO和LUMO的能量间隔加大;,5,将SCF=SLEAZY收敛限放松,或者SCFCON关键词 改变收敛限; 6,对于使用了弥散函数的情形,需要提高积分精度, 加关键词SCF=NoVarAcc; 7,略微改变分子几何结构,得到收敛的波函数后,再使 用原来的几何构型读分子轨道计算; 8,考虑使用不同级别的计算方法和基组(小一些的); 9,使用强制收敛方法SCF=QC或者SCF=DM。,Koopmans 原理,BSSE(基组重叠误差),1961年 Ransil使用HF方法计算发现 He-He势能曲线上存在

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论