CalculationofMolecularIntegrals分子积分的计算课件.ppt

上传人(卖家):三亚风情 文档编号:3376308 上传时间:2022-08-25 格式:PPT 页数:16 大小:360.01KB
下载 相关 举报
CalculationofMolecularIntegrals分子积分的计算课件.ppt_第1页
第1页 / 共16页
CalculationofMolecularIntegrals分子积分的计算课件.ppt_第2页
第2页 / 共16页
CalculationofMolecularIntegrals分子积分的计算课件.ppt_第3页
第3页 / 共16页
CalculationofMolecularIntegrals分子积分的计算课件.ppt_第4页
第4页 / 共16页
CalculationofMolecularIntegrals分子积分的计算课件.ppt_第5页
第5页 / 共16页
点击查看更多>>
资源描述

1、3.3 SCF Procedure(自洽场计算)Roothaan方程FC=SCE是广义的矩阵本征值方程lklkijijkjilklijDhF,)|()|(2klkllkDCCDocc,在实基函数空间 F=F(D)=F(C)欲得到正确的C、E,必须借助自洽叠代处理。构建分子选定基组分子对称性分析计算分子积分MO 初始猜测自洽叠代求解Roothaan方程 Fc=Sc电子布居分析计算偶极矩分子性质计算绘制 MO、电子密度和自旋密度分布等图形电子相关校正(CI 或 MP 微扰处理)Hatree-Fock calculationPost-Hartree-Fock procedure从从头头算算法法基基本

2、本框框图图ab initioSCFprocedure3.3.1 Logical Structure of SCF calculations(自洽场计算的逻辑结构)选取一组合适的MO试探波矢(trial vectors)Roothaan方程借助自洽叠代求解。基本步骤:)0()0(2)0(1(0),mcccT用以形成试探矩阵)0()0()0()0()0()0(GhFTTD解 Roothaan 方程 F(0)c(1)=(1)Sc(1)得)1()1()1(2)1(1)1()1()1(2)1(1,nmnmccccoccupied 取 作为试探波矢)1()1(2)1(1(1),mcccT用以形成试探矩阵)

3、1()1()1()1()1()1(GhFTTD解 Roothaan 方程 F(1)c(2)=(2)Sc(2)得)2()2()2(2)2(1)2()2()2(2)2(1,nmnmccccoccupied 重复步骤 直至 E、T 或 D 达到收敛为止。计 算双电子积分(i j|k l)存于存于硬盘硬盘计 算单电子积分形成密度矩阵D=TT+计算电子总能量E=Tr D(h+F)hij,Sij 存于存于硬盘硬盘(i j|k l)形成Fock矩阵 F密度矩阵初始猜测初始猜测hij解Roothaan方程 Fc=SCESij 自洽场计算的 逻辑结构 1(非直接法从头算)(非直接法从头算)能量或密度矩阵是否达到

4、预定收敛精度?SCF完完 成成是否开始下一次叠代Indirect ab initio形成密度矩阵D=TT+计算电子总能量E=Tr D(h+F)计 算双电子积分(i j|k l)形成Fock矩阵 F密度矩阵初始猜测初始猜测计 算单电子积分hij解Roothaan方程Fc=SCESij 自洽场计算的 逻辑结构 2(直接法从头算)(直接法从头算)能量或密度矩阵是否达到预定收敛精度?SCF完完 成成是否开始下一次叠代Direct Ab InitioSixty-Four-Dollar Problems in SCF(自洽场计算中的几个关键问题)双电子积分是从头算的“bottleneck”。SCF叠代可占

5、总计算量 1/31/2。提高效率亦至关重要The follwing problems are vital:How to choose the MO trial vectors?How to recognize the convergence of iterations?How to accelerate the iteration process?3.3.2 MO Initial Guess(MO初始猜测)数学上要求方程的解是存在和唯一的,故自洽解应与试探函数的选择无关。但实践表明,如试探MO选择不当则SCF收敛慢甚至失败试探函数c(0)的选择称为“MO的初始猜测”对c(0)的基本要求:对称性

6、(不可约表示构成)与精确对称性(不可约表示构成)与精确MO解一致解一致 在在此前提下令前提下令 c(0)尽可能接近于精确解尽可能接近于精确解当使用当使用core hamiltonian收敛遇困难时,可尝试改收敛遇困难时,可尝试改用用第第类猜测类猜测。取单电子哈密顿的本征矢,适于大多数体系且简单易行。(i.e.core hamiltonian)常用的MO初始猜测)0(c 取半经验计算方法所得MO向从头算的基集投影作为价MO试探本征矢,加内层自洽AO作为初始猜测。在 Gaussian 与 HyperChem 中提供:Projected Hckel(i.e.EHMO);Projected CNDO;

7、Projected INDO)0()0()0(Schc 将分子碎片(molecular fragments)的MO拼装成大分子的试探函数。e.g.联苯分子(PhPh)的初始猜测可用苯的分子轨道拼装。但通用软件未提供此技术,需自行输入HyperChem中提供的4种初始猜测选择CycleTrialvector3.3.3 Recognition of SCF Convergence(自洽场收敛的判别)自然叠代相邻两次循环的关系自洽:第第q次与第次与第q+1次叠代结果之差达到了要求精度次叠代结果之差达到了要求精度q-th:(q+1)-th:)1(qD)1(qF)1(qc)1(qE)()(qq,c)(q

8、D)(qF)(qc)(qE)1()1(qq,cScFc ScFc(公式:E=TrD(h+F))常用以下三类判据(criteria)Total energy criterion(总能量判据)主要用于半经验法(CNDO、AM1、PM3 etc.),在一些简易从头算软件中亦有时采用通用标准:10-810-9 a.u.10-510-6 kcal/mol 1 kcal/mol 4.184 kJ/mol 1.594 10-3 a.u.HyperChem 的缺省选值取 10-3 kcal/mol 不妥应取 10-5 kcal/mol )1()(qqEEHyperChem的对话框 MO waefunction

9、 criterion(波函数判据)比较对应MO本征矢每个分量的最大误差,比总能量判据要求严格通用标准:10-6实践表明:E 收敛到10-8 a.u.时,c 收敛到10-510-6(两者精度约差两者精度约差 2 个量级个量级)1()(Max qqcc Density matrix criterion(密度矩阵判据)通用标准:10-6与波函数判据精度等级相同,D TT。但因D是实对称阵,独立矩阵元数为n(n-1)/2+1。故在从头算和DFT的自洽叠代中多采用D判据)1()(Max qqDD为节约机时,G98的几何优化作业中在E10-5的远离平衡几何点取 10-4,优化到E10-5后改取 10-6在

10、G98单点计算(single point calculation,关键词SP)作业中,缺省值取 10-4。如欲得到更精确的能量和波函数,应在命令行添加关键词选项 SCF=Tight3.3.4 Acceleration of Convergence(自洽叠代的收敛和加速)按自然叠代处理(即不加任何附加措施)一般收敛较慢,有些分子甚至不能收敛快收敛 q 20慢收敛 q=2064不收敛 振荡Iteration qE三类典型的叠代过程vibrationalslowfast消除振荡和加速收敛的措施假想分子 HCNO最难收敛。70年代初后曾长期用作攻克收敛技术的测试样本Pulay 于1980年提出“叠代空间直接反演”处理技术(DIIS,即 Direct Inversion in Iterative Space),妥善解决了收敛问题,使多数分子可在叠代20次以内达到高精度自洽H C:N:O:.DIIS 是迄今提出过的最有效的保证和加速收敛技术。已被各种通用软件包所采用八十年代末,DIIS 法被推广应用于共轭梯度法的分子几何优化,因后者也是通过反复叠代来完成432

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 办公、行业 > 各类PPT课件(模板)
版权提示 | 免责声明

1,本文(CalculationofMolecularIntegrals分子积分的计算课件.ppt)为本站会员(三亚风情)主动上传,163文库仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。
2,用户下载本文档,所消耗的文币(积分)将全额增加到上传者的账号。
3, 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(发送邮件至3464097650@qq.com或直接QQ联系客服),我们立即给予删除!


侵权处理QQ:3464097650--上传资料QQ:3464097650

【声明】本站为“文档C2C交易模式”,即用户上传的文档直接卖给(下载)用户,本站只是网络空间服务平台,本站所有原创文档下载所得归上传人所有,如您发现上传作品侵犯了您的版权,请立刻联系我们并提供证据,我们将在3个工作日内予以改正。


163文库-Www.163Wenku.Com |网站地图|