第10章 量子化学计算简介.ppt

2020/6/3,第十章计算量子化学简介,2020/6/3,2020/6/3,计算机的发展,经典力学,量子力学,分子和物质结构问题,分子力学,分子动力学,量子化学,计算化学,1687Newton,Heisenberg1926SchrdingerBorn,1946WestheimerMayer,1927HeitlerLondon,统计力学,Before1900Boltzmann,1957AlderWainwright,10.1引论,2020/6/3,基础理论,应用,量子化学,多体理论,化学键理论密度矩阵理论传播子理论,分子轨道理论价键理论配位场理论,计算方法,从头计算法abinitio密度泛函法DFT超自洽场方法半经验法,量子有机量子无机量子生物量子药物分子设计材料中的应用,10.1引论,2020/6/3,1.类氢原子轨道的具体形式,10.1.1原子轨道,10.1.1原子轨道,2020/6/3,10.1.1原子轨道,2020/6/3,进行量子化学从头计算首先应选定一组合适的原子轨道基函数,最常见的是STO-GTO系基函数.Slater函数的形式为：,该函数适合于描述电子云的分布，在半经验法计算中常采用。轨道指数，决定原子轨道的衰减速度。,2.Slater型原子轨道(STO),径向函数,轨道指数,10.1.1原子轨道,2020/6/3,特点：原子轨道计算中常用，有很好的渐近性质，与真实的轨道形状接近；在近核部分表现很好；STO相互不正交；对于3中心和4中心积分难以解析求算。,10.1.1原子轨道,2020/6/3,3.Gauss型原子轨道(GTO),10.1.1原子轨道,2020/6/3,的意义：角量子数,10.1.1原子轨道,2020/6/3,角度部分,10.1.1原子轨道,2020/6/3,在离核距离远处衰减太快近核的地方不够尖锐能很方便地解析求算3中心和4中心积分,10.1.1原子轨道,2020/6/3,10.1.1原子轨道,2020/6/3,10.1.2基组及其选择,实际计算中采用STO-GTO基组，以STO作为自洽场的基函数,而每个STO用若干个GTO来拟合,解出的分子轨道用STO的线性组合来表示.计算中常采用的基组：(1)STO-NG型极小基组:每个占据的原子轨道只对应1个STO,而每个STO则用N个GTO来拟合,N一般为16.,拟合时采用最小二乘法,并结合归一化条件来确定N个GTO组合成1个STO的组合系数和轨道指数。例如一个标准STO-3G的结果为：,10.1.2基组及其选择,2020/6/3,STO-NG的特点是，从1s轨道到3d、4f轨道，每个轨道均用相同几个GTO拟合，其中STO-3G应用最广，它是量化计算中极小基的首选基组，常用于几何结构优化、寻找反应过渡的试探函数、较大体系的量化计算等。例如，C的原子轨道为1s2s2px2py2px,先用5个STO拟合5个原子轨道，再用15个GTO拟合5个STO。,10.1.2基组及其选择,2020/6/3,(2)分裂价基N-31G和3-21G:每个原子内层轨道用N个GTO拟合单的STO；而价层分为内轨(I)和外轨(O),对应指数为z，和z的两个STO基函数,分别用N，个GTO和N”个GTO来拟合。这种基组称为分裂价基.6-31G：内层轨道用6个GTO拟合1个STO，1个STO拟合1个与原子轨道；价层分为内外轨，即双z，I层用3个GTO拟合1个STO(I)，O层用1个GTO(O)拟合1个STO(O)，而1个STO(I)和1个STO(O)拟合1个价层原子轨道。分裂价基基组有、与等可选。它比STO-NG仅价轨道增加一倍，但计算精度提高不少。,10.1.2基组及其选择,2020/6/3,(3)双扩展基组:例如,6-31G*或6-31(d,p)表示在6-31G基组的基础上,再增加极化函数来描述原子轨道,在此是用p型函数来描述s型轨道的极化,用d型函数来描述p型轨道的极化.这样原子轨道包括内存轨道、价层轨道，再加上角量子数更高的原子轨道。N-31G*表示体系中非氢元素的原子都加上极化函数(即价轨道为s、p、d，则分别加上p、d、f极化函数)。N-31G*表示体系所有元素(包括氢)的原子都加上极化函数。现以乙烯分子为例，说明采用不同基组时原子轨道数目。,10.1.2基组及其选择,2020/6/3,STO-3GC:1s2s2px2py2px，H:1s,共14个原子轨道、42个初始GTO。3-21GC:2s、2p为双层，H:1s为双层,共2(1+24)+24=26个STO、42个初始GTO。6-31G*C:2s、2p层外加d极化函数，共2(1+24+45)+24=66个STO、100个初始GTO.6-31G*在6-31G*基础上增加H1s所外加得p极化函数，共78个STO、112个初始GTO.随着基组增大，原子轨道数逐步增加，而计算机容量与时间则与基组数的四次方成正比。,10.1.2基组及其选择,2020/6/3,上述的基组还可以对价层轨道进行三指数基组分裂，如6-311G，其内层轨道由6个GTO组成，价层轨道被分成3个部分，分别由3个、1个和1个GTO组成。此外，除了极化函数外，还可增加弥散函数来增加计算精度。通常是在“G”前添加“+”和“+”来表示。“+”表示弥散到s轨道，“+”表示弥散到s轨道和s轨道.,其它基组：有效核势ECP和赝势基组(将内层电子和核的势能用函数代替)，如Lanl2dz，SDD等;相关一致基组(考虑了相关能，基函数非常大，计算结果接近HF极限，本身包含有极化函数，还可以加上弥散函数)，如cc-pVDZ,cc-pVTZcc-pV6Z，aug-cc-pVTZ等.,10.1.2基组及其选择,2020/6/3,基组总结,MOLCAOAO(STO)STO-NGTO极小基：STO-NGTO双基：使用两个收缩基来近似一条轨道分裂价基：内层轨道使用收缩基，价层使用多基极化基：比轨道更高角动量的基函数弥散基：轨道指数较小，更弥散的基函数,10.1.2基组及其选择,2020/6/3,2020/6/3,1.基本近似基本常数：Plank常数、电子质量、电子电量、原子序数依据的：（1）非相对论近似:认为电子质量等于其静止质量。（2）Born-Oppenherimer近似(绝热近似):认即将核的运动和电子运动分离开来进行处理。由于原子核质量一般约比电子的质量大103105，分子中核的运动要比电子的运动慢近千倍。因此在考察电子运动时，可把核近似看成不动。（3）轨道近似(又称单电子近似):原子的单电子波函数称为原子轨道，分子的单电子波函数称为分子轨道。蛋电子近似是把N个电子体系的总波函数写成N个单电子波函数的乘积：,10.2.1从头计算法,10.2计算方法,2020/6/3,从头计算法(abinitio)就是基于上述3个基本近似后,不在做任何近似来求解体系的Schrodinger方程.这3个基本近似即为从头计算法的“头”.而引入3个基本近似后,Schrodinger方程的具体表达形式为Hartree-Fock-Roothaan方程-LCAO法.,HF方法中的矩阵元求算，采取Fock算符的本征函数为基Roothaan提出使用原子轨道作为基函数，从而将微分方程求解问题简化为代数方程求解问题,(1,2,N)=1(1)2(2)N(N),10.2.1从头计算法,2020/6/3,2.从头计算的误差来源及其校正三个基本近似给自洽场从头算带来了一定的误差.(1)相对论误差对于内蹭电子,尤其是重原子及其所组成的体系的内层电子,相对论效应尤为突出.(2)轨道近似误差(即电子相关能的误差)这是由单组态(即独立电子模型)带来的,特别是把自旋相反的交换积分视为零,相关能就考虑不够.(3)B-O近似误差相对较小.相对论误差的校正可采用微扰法或相对论的H-F理论.相关能的校正有组态相互作用法(CI)和自旋非限制Roothaan方程法.,10.2.1从头计算法,2020/6/3,二.半经验计算法,1.概论在3个层次上对Roothaan方程进行简化:波函数、积分、Hamilton算符。由此形成3类近似方法：(1)单电子近似完全不考虑双电子的作用而挑选等效的Hamilton量，如EHMO法。(2)用统计平均模型计算交换位能的X方法。(3)以零微分重叠(ZDO)近似为基础的计算方法.如CNDO,PM3等,10.2.2半经验计算法,2020/6/3,ZDO(零微分重叠近似),根据上述近似可以得到ZDO下的Fock矩阵元对角元：非对角元：,基本近似：,(1).不同轨道之间的重叠积分等于零:(2).对于双电子积分,10.2.2半经验计算法,2020/6/3,ZDO存在的问题1）.基变化下不能保证双电子积分仍等于零.例如：,2）.对于电子与其他核实的相互作用考虑得不够。,10.2.2半经验计算法,2020/6/3,Dewar-理论,参数化一般原则：单中心电子排斥积分：拟合原子能级的实验数据双中心电子排斥积分：结合单中心排斥积分和核间距计算得到其他积分：从参数化近似式计算得到，依据近似级别不同采用不同的拟合方式,10.2.2半经验计算法,2020/6/3,取可精确得到的实验数据做参照分子生成热,参数化策略：,分子几何构型,偶极矩,第一电离电势.,10.2.2半经验计算法,2020/6/3,2.EHMO(ExtendedHuckelMolecularOrbitalMethod)是对HMO的扩展,由Hoffmann在1963年所提出.改进点:.可用于非共轭体系方程中要考虑全部价电子,在LCAO中应包括所有的价轨道，认为内层电子仍保持原来的原子轨道状态；忽略所有的双电子项。.对所有重叠积分Sij都要进行计算计算时原子轨道一般用Slater函数.,10.2.2半经验计算法,2020/6/3,(3).考虑了所有原子之间的相互作用：每个价层轨道使用一个STO表示,所有重叠积分直接计算，交换积分用下式计算:,常数K取1.752.00之间.,(4).矩阵元Hmm：实由验测得的原子轨道的电离势*不适用于优化结构计算，z目前仅用于固体和有机聚合物的能带结构计算,10.2.2半经验计算法,2020/6/3,3.CNDO(CompleteNeglectofDifferentialOverlapMethod),采用了忽略微分重叠近似,对所有三中心和四中心积分都加以忽略,只剩下单中心和双中心积分.相邻原子间的核实积分当做参数,用半经验方法加以确定.只考虑价电子而把内层电子和原子核一起视为不变的原子实.,对于双电子积分作如下处理：m和l属于同一原子Am和l属于不同原子A和B,10.2.2半经验计算法,2020/6/3,CNDO近似下的Fock矩阵元：,不再简单地忽略不同原子上轨道之间的共振积分：,A原子中m轨道上的电子在A核实作用下的能量,A原子中m轨道上的电子在B核实作用下的能量,10.2.2半经验计算法,2020/6/3,CNDO下的FOCK矩阵元,10.2.2半经验计算法,2020/6/3,CNDO计算所需参数:重叠积分S,实Hamilton量元素U、VAB，电子排斥积分AB和成键参量AB的数据。CNDO按参量方案的不同可分为：CNDO/1、CNDO/2、CNDO/SW、CNDO/BW、CNDO/S。,10.2.2半经验计算法,2020/6/3,在CNDO/2的基础上不再忽略单中心双电子积分,Slater-condon参量F0，仍使用STO直接计算G，F2通过与光谱得到的原子能级拟合得到,m,l来自同一原子,4.INDO方法(IntermediateNeglectofDifferentialOverlapMethod),10.2.2半经验计算法,2020/6/3,INDO下的FOCK矩阵元,10.2.2半经验计算法,2020/6/3,同一原子，同一原子轨道：单中心单电子积分同一原子，不同原子轨道：单中心双电子积分不同原子，不同原子轨道：双中心双电子积分,10.2.2半经验计算法,2020/6/3,5.MINDO(ModifiedIntermediateNeglectofDifferentialOverlapMethod),一些基本假定：保留全部单中心排斥积分；双中心电子排斥积分选成=gAB.单中心核矩阵元为:,双原子核实积分取为:,6.MNDO(ModifiedNeglectofDiatomicOverlapMethod),与INDO的主要区别在于对双中心排斥积分和双中心核实-电子吸引积分的处理.属NDDO法.,10.2.2半经验计算法,2020/6/3,7.AM1(AustinModel1),只对MNDO中的排斥函数(CRF)进行了修正:,8.PM3(Parametric3),是一种基于MNDO模型的新参量化方法.对每一个非氢原子有由18个参数组成的一个参数集,氢原子只有11个参数.,10.2.2半经验计算法,2020/6/3,计算方法总结与评价,10.2.2半经验计算法,2020/6/3,10.2.2半经验计算法,2020/6/3,10.2.3密度泛函理论（DFT）Hohenberg和Kohn第一定理：Hohenberg和Kohn第二定理能量表示式Kohn-Sham方程交换函数与相关函数,10.2.3DFT法,2020/6/3,2020/6/3,目前常用的量子化学程序,Gaussian(98,03，09)HyperchemGAMESSMolproCrystalMOPAC,10.3计算程序简介,2020/6/3,10.3.1.Hyperchem简介,Hyperchem是一款以质量高，灵活易操作而闻名的分子模拟软件。通过3D对量子化学计算，分子力学及动力学进行模拟。Hyperchem可提供比其它Windows软件更多的模拟工具：图形界面，半经验计算方法（AM1，PM3等），UHF，RHF和CI，7.0版新增加的密度泛函。可进行单点能，几何优化，分子轨道分析，蒙特卡罗和分子力学计算，预测可见-紫外光谱。目前，已至8.010版本。,1.概述,10.3.1Hyperchem,2020/6/3,2.程序的主要功能：1)、建立、显示分子骨架模型该程序包含有实验数据构成的数据库（键长、键角、二面角等），可以迅速构成分子体系骨架，建立起3D、2D分子模型。2)、分子化学反应性质可以给出分子的静态性质，包括势能、势能的梯度、静电势、分子轨道能量、基态、激发态分子轨道系数、振动分析、振动频率、振动模式、红外吸收强度、自旋密度-ESR光谱的耦合常数、预测化学反应的位置、说明化学反应的途径和机制、解释分子的动力学行为等等。,10.3.1Hyperchem,2020/6/3,3)、计算结果图形化分子轨道、静电势的图形、总电荷密度的分布、总自旋密度的图形等等。,4)、分子构型优化SteepestDescentConjugateGradientFletcher-ReevesPolak-RibiereBlock-Digonal、Newton-RaphsonBloackDiggonalEigenvectotFollowing,10.3.1Hyperchem,2020/6/3,5)、计算过渡态理论上预测化学反应的过渡结构，活化能可由实验测定，用于说明化学键的形成和断裂过程。过渡态是沿反应途径势能面的最高点，相应的结构有一个负本征值、一个虚频率，二阶能量导数相对于几何参数的矩阵（Hessian）在一个方向极大其它方向极小仅有一个负本征值。,6)、计算方法量子化学计算，可以包括100个原子的分子体系，与计算机的内存、硬盘和计算时间有关。从头算类方法从头算已经引入了优选的基组和算法，在LCAO-MO的基础上，没有引入其它近似。,10.3.1Hyperchem,2020/6/3,HyperChem使用S-G型轨道作为基函数STO-1G(H-He)，STO-2G(H-Xe)，STO-3G(H-Xe)，STO-4G(H-Xe)，STO-5G(H-Xe)，STO-6G(H-Xe)，3-21G、3-21G*、3-21G*(H-Ar)，4-21G、4-21G*、4-21G*(H-Ne)，6-21G、6-21G*、6-21G*(H-Ar)，4-31G、4-31G*、4-31G*(H-Ne)，5-31G、5-31G*、5-31G*(H-F)，6-31G、6-31G*、6-31G*(H-Ar)，6-311G、6-311G*、6-311G*(H-Ar),重原子加f轨道，氢原子加D轨道极化。可以自行扩展基组，当然基组参数的调整应当是科学的。可以对不同的原子使用不同的基函数.,10.3.1Hyperchem,2020/6/3,半经验类方法半经验方法适合于有机化合物，可以得到定性的、有意义的大分子量化计算结果。但对于无机化合物、主族元素甚至最好的半经验方法也难以得到好的结果。EHMO、CNDO、INDO、MINDO/3、MNDO、AM1、PM3、ZINDO/1、ZINDO/2EHMO是最粗糙的半经验近似方法。可以提供分子轨道的形状和能量形成近似的电子密度图。ZINDO也是修正的INDO版本，可以包括过渡金属，给出参数化的紫外光谱。NDDO是MNDO、AM1、PM3方法的基础。CNDO、INDO的参数化方案是由Pople发展起来的，MINDO/3、MNDO、AM1是由Dewar等人发展起来的，比CNDO、INDO更接近实验值。,10.3.1Hyperchem,2020/6/3,PM3基于忽略双原子近似，基于AM1方法。Z-INDO适合于第一周期、第二周期过渡元素。Moller-Plesset微扰理论是说明相关能的另一种方法，MP2.MP3MP4可以使用MP2方法。计算MP2相关能-仅对单点适用。半经验、从头算方法均提供：能量及能量梯度、偶极矩、总电子密度、总电子自旋密度（闭壳层体系为零）、静电势、生成热、轨道能级、简正振动模式和频率、红外光谱强度、紫外吸收频率和强度、（log、hin文件）、,10.3.1Hyperchem,2020/6/3,7)、使用程序步骤首先准备初始构型：初始结构可由HyperChem数据库提供选择计算方法：从头算、半经验从头算基组的选择电荷和自旋多重度：闭壳层体系的多重度应该为1具有未成对电子的体系多重度应大于1。收敛限、迭代次数的选择。UHF、RHF、ROHF的选择：RHF适合于闭壳层体系，UHF适合于所有体系。CI计算：闭壳层单组态基态、半经验、从头算、双重态、三重态开壳层基态。MINDO/3MNDOAM1PM3的参数中已经包括了CI的影响。LOG文件作为计算结果的输出。,10.3.1Hyperchem,2020/6/3,2.程序基本使用方法以Hyperchem7.0为例，说明该程序的基本使用方法。(1).启动程序,chem,在屏幕上，双击绿色烧杯可得右边的工作界面：,10.3.1Hyperchem,2020/6/3,窗口各部分功能简介标题名称：最大、最小化、退出按钮菜单条：,工具条,工作区状态行,FilePreferences:对Hyperchem程序进行初始化设置,10.3.1Hyperchem,2020/6/3,(2)打开已存在的数据文件FileOpen文件夹文件名(Hyperchem7/Samples/Aromatic/C60为例)（或点击工具条中的打开工具图标）,选择分子图形的显示方式Display-Labels可以选择原子、化学键等标记方式；Dispay-Rendering选择分子结构显示方式,(3).建立计算分子的数据文件以对硝基苯酚为例：选择Build-DefaultElement可以显示指定元素的基本性质：,10.3.1Hyperchem,2020/6/3,选择绘图工具后，得到碳碳骨架。左键、右键增减键序,把2D骨架转换为3D图像,选择其它原子进行相应的碳原子替换。,BuildH2/17=6,18=5,40=1/2;3/5=1,6=6,11=9,16=1,25=1,30=1/1,2,3;4/1;5/5=2,32=1,38=5/2;6/7=2,8=2,9=2,10=2,28=1/1;99/5=1,9=1/99;LeaveLink1atFriApr0801:21:372005,MaxMem=0cpu:1.0,重新打印作业输入和标准途径,10.3.2Gaussian程序,2020/6/3,作业执行过程,1/38=1/1;2/17=6,18=5,40=1/2;3/5=1,6=6,11=9,16=1,25=1,30=1/1,2,3;4/1;5/5=2,32=1,38=5/2;6/7=2,8=2,9=2,10=2,28=1/1;99/5=1,9=1/99;,L101L202L301L302L303L401L502L601L9999,L101读取标题和分子输入部分；L202坐标重新定向，计算对称性，检查变量；L301产生基组信息；L302计算重叠积分，动能积分和势能积分；L303计算多极积分；L401形成分子轨道初猜；L502自洽场迭代求解；L601布居数和相关分析；L9999完成计算和输出,10.3.2Gaussian程序,2020/6/3,(EnterC:G03Wl101.exe)-h2osp-SymbolicZ-matrix:Charge=0Multiplicity=1O0.0.0.H0.6620.0.468H-0.6620.0.468IsotopesandNuclearProperties:Atom123IAtWgt=1611AtmWgt=15.99491461.00782501.0078250IAtSpn=011AtZEff=0.00000000.00000000.0000000AtQMom=0.00000000.00000000.0000000AtGFac=0.00000002.79284602.7928460LeaveLink101atFriApr0801:21:382005,MaxMem=6291456cpu:0.0,10.3.2Gaussian程序,2020/6/3,(EnterC:G03Wl202.exe)Inputorientation:-CenterAtomicAtomicCoordinates(Angstroms)NumberNumberTypeXYZ-1800.0000000.0000000.0000002100.6620000.0000000.468000310-0.6620000.0000000.468000-Distancematrix(angstroms):1231O0.0000002H0.8107210.0000003H0.8107211.3240000.000000StoichiometryH2OFrameworkgroupC2VC2(O),SGV(H2)Deg.offreedom2FullpointgroupC2VNOp4LargestAbeliansubgroupC2VNOp4LargestconciseAbeliansubgroupC2NOp2Standardorientation:-CenterAtomicAtomicCoordinates(Angstroms)NumberNumberTypeXYZ-1800.0000000.0000000.0936002100.0000000.662000-0.3744003100.000000-0.662000-0.374400-Rotationalconstants(GHZ):1289.0084222572.1186153396.2468433LeaveLink202atFriApr0801:21:402005,MaxMem=6291456cpu:1.0,10.3.2Gaussian程序,2020/6/3,(EnterC:G03Wl301.exe)Standardbasis:6-31G(6D,7F)Thereare7symmetryadaptedbasisfunctionsofA1symmetry.Thereare0symmetryadaptedbasisfunctionsofA2symmetry.Thereare2symmetryadaptedbasisfunctionsofB1symmetry.Thereare4symmetryadaptedbasisfunctionsofB2symmetry.Integralbufferswillbe262144wordslong.Raffenetti1integralformat.Two-electronintegralsymmetryisturnedon.13basisfunctions,30primitivegaussians,13cartesianbasisfunctions5alphaelectrons5betaelectronsnuclearrepulsionenergy10.8432720049Hartrees.IExCor=0DFT=FEx=HFCorr=NoneExCW=0ScaHFX=1.000000ScaDFX=1.0000001.0000001.0000001.000000IRadAn=0IRanWt=-1IRanGd=0ICorTp=0NAtoms=3NActive=3NUniq=2SFac=2.76D+00NAtFMM=60Big=FLeaveLink301atFriApr0801:57:332005,MaxMem=6291456cpu:1.0,10.3.2Gaussian程序,2020/6/3,(EnterC:G03Wl302.exe)NPDir=0NMtPBC=1NCelOv=1NCel=1NClECP=1NCelD=1NCelK=1NCelE2=1NClLst=1CellRange=0.0.One-electronintegralscomputedusingPRISM.One-electronintegralsymmetryusedinSTVIntNBasis=13RedAO=TNBF=7024NBsUse=131.00D-06NBFU=7024LeaveLink302atFriApr0801:57:352005,MaxMem=6291456cpu:1.0(EnterC:G03Wl303.exe)DipDrv:MaxL=1.LeaveLink303atFriApr0801:57:362005,MaxMem=6291456cpu:0.0,10.3.2Gaussian程序,2020/6/3,(EnterC:G03Wl401.exe)HarrisfunctionalwithIExCor=205diagonalizedforinitialguess.ExpMin=1.61D-01ExpMax=5.48D+03ExpMxC=8.25D+02IAcc=1IRadAn=1AccDes=1.00D-06HarFok:IExCor=205AccDes=1.00D-06IRadAn=1IDoV=1ScaDFX=1.0000001.0000001.0000001.000000HarrisEn=-76.0428456360852Initialguessorbitalsymmetries:Occupied(A1)(A1)(B2)(A1)(B1)Virtual(A1)(B2)(B1)(A1)(B2)(B2)(A1)(A1)Theelectronicstateoftheinitialguessis1-A1.LeaveLink401atFriApr0801:57:372005,MaxMem=6291456cpu:0.0,10.3.2Gaussian程序,2020/6/3,(EnterC:G03Wl502.exe)Warning!Cutoffsforsingle-pointcalculationsused.ClosedshellSCF:RequestedconvergenceonRMSdensitymatrix=1.00D-04within128cycles.RequestedconvergenceonMAXdensitymatrix=1.00D-02.Requestedconvergenceonenergy=5.00D-05.Nospecialactionsifenergyrises.UsingDIISextrapolation,IDIIS=1040.Integralsymmetryusagewillbedecideddynamically.KeepR1integralsinmemoryincanonicalform,NReq=424251.IEnd=18776IEndB=18776NGot=6291456MDV=6285583LenX=6285583SymmetrynotusedinFoFDir.MinBra=0MaxBra=1Meth=1.IRaf=0NMat=1IRICut=1DoRegI=TDoRafI=FISym2E=0JSym2E=0.,10.3.2Gaussian程序,2020/6/3,Cycle5Pass1IDiag1:E=-75.9273519590932Delta-E=-0.000003762089Rises=FDamp=FDIIS:error=5.91D-05atcycle5NSaved=5.NSaved=5IEnMin=5EnMin=-75.9273519590932IErMin=5ErrMin=5.91D-05ErrMax=5.91D-05EMaxC=1.00D-01BMatC=4.02D-08BMatP=1.89D-06IDIUse=1WtCom=1.00D+00WtEn=0.00D+00Coeff-Com:0.711D-02-0.177D-010.103D+00-0.460D+000.137D+01Coeff:0.711D-02-0.177D-010.103D+00-0.460D+000.137D+01Gap=0.746Goal=NoneShift=0.000RMSDP=4.58D-05MaxDP=1.68D-04DE=-3.76D-06OVMax=2.54D-04SCFDone:E(RHF)=-75.9273519591A.U.after5cyclesConvg=0.4576D-04-V/T=1.9895S*2=0.0000KE=7.673560164965D+01PE=-2.025358608090D+02EE=3.902963519532D+01LeaveLink502atFriApr0801:57:412005,MaxMem=6291456cpu:2.0,10.3.2Gaussian程序,2020/6/3,(EnterC:G03Wl601.exe)CopyingSCFdensitiestogeneralizeddensityrwf,ISCF=0IROHF=0.*PopulationanalysisusingtheSCFdensity.*Orbitalsymmetries:Occupied(A1)(A1)(B2)(A1)(B1)Virtual(A1)(B2)(B2)(B1)(A1)(A1)(B2)(A1)Theelectronicstateis1-A1.Alphaocc.eigenvalues-20.52887-1.44410-0.81350-0.57281-0.51130Alphavirt.eigenvalues-0.234900.331091.135361.163731.17135Alphavirt.eigenvalues-1.436361.500541.68694Condensedtoatoms(allelectrons):1231O8.0907140.3123390.3123392H0.3123390.374540-0.0445753H0.312339-0.0445750.374540Mullikenatomiccharges:11O-0.7153922H0.3576963H0.357696SumofMullikencharges=0.00000Atomicchargeswithhydrogenssummedintoheavyatoms:11O0.0000002H0.0000003H0.000000SumofMullikencharges=0.00000,10.3.2Gaussian程序,2020/6/3,Electronicspatialextent(au):=17.2272Charge=0.0000electronsDipolemoment(field-independentbasis,Debye):X=0.0000Y=0.0000Z=-2.4602Tot=2.4602Quadrupolemoment(field-independentbasis,Debye-Ang):XX=-6.8696YY=-4.2696ZZ=-6.1388XY=0.0000XZ=0.0000YZ=0.0000TracelessQuadrupolemoment(field-independentbasis,Debye-Ang):XX=-1.1103YY=1.4898ZZ=-0.3795XY=0.0000XZ=0.0000YZ=0.0000Octapolemoment(field-independentbasis,Debye-Ang*2):XXX=0.0000YYY=0.0000ZZZ=-1.6189XYY=0.0000XXY=0.0000XXZ=-0.5144XZZ=0.0000YZZ=0.0000YYZ=-0.9389XYZ=0.0000Hexadecapolemoment(field-independentbasis,Debye-Ang*3):XXXX=-4.8330YYYY=-4.5460ZZZZ=-5.2984XXXY=0.0000XXXZ=0.0000YYYX=0.0000YYYZ=0.0000ZZZX=0.0000ZZZY=0.0000XXYY=-1.6915XXZZ=-1.7089YYZZ=-1.5145XXYZ=0.0000YYXZ=0.0000ZZXY=0.0000N-N=1.084327200490D+01E-N=-2.025355816034D+02KE=7.673560164965D+01SymmetryA1KE=6.825101673560D+01SymmetryA2KE=0.000000000000D+00SymmetryB1KE=4.561871027519D+00SymmetryB2KE=3.922713886524D+00NoNMRshieldingtensorssonospin-rotationconstants.LeaveLink601atFriApr0801:57:442005,MaxMem=6291456cpu:2.0,2020/6/3,(EnterC:G03Wl9999.exe)1|1|UNPC-UNK|SP|RHF|6-31G|H2O1|PCUSER|08-Apr-2005|0|#PHF/6-31G|h2osp|0,1|O,0,0.,0.,0.|H,0,0.662,0.,0.468|H,0,-0.662,0.,0.468|Version=x86-Win32-G03RevB.02|State=1-A1|HF=-75.927352|RMSD=4.576e-005|Dipole=0.,0.,0.9679067|PG=C02VC2(O1),SGV(H2)|ERRORISAHARDYPLANT.ITFLOURISHETHINEVERYSOIL.MARTINFARQUHARTUPPERJobcputime:0days0hours0minutes18.0seconds.Filelengths(MBytes):RWF=11Int=0D2E=0Chk=7Scr=1NormalterminationofGaussian03atFriApr0801:57:462005.,10.3.2Gaussian程序,2020/6/3,Gaussian中的分子结构优化过程与输出文件示例乙烯分子优化,10.3.2Gaussian程序,2020/6/3,#PRHF/6-31G(d)OptTestEthyleneGeometryOptimization01CC1CCH1CH2HCCH1CH2HCC3180.H2CH1HCC3180.H2CH1HCC4180.Variables:CC=1.31CH=1.07HCC=121.5,输入文件,10.3.2Gaussian程序,2020/6/3,优化作业执行过程,1/18=20,38=1/1,3;2/9=110,17=6,18=5,40=1/2;3/5=1,6=6,7=1,11=1,16=1,25=1,30=1/1,2,3;4/7=1/1;5/5=2,38=5/2;6/7=2,8=2,9=2,10=2,28=1/1;7/1,2,3,16;1/18=20/3(1);99/99;2/9=110/2;3/5=1,6=6,7=1,11=1,16=1,25=1,30=1/1,2,3;4/5=5,7=1,16=3/1;5/5=2,38=5/2;7/1,2,3,16;1/18=20/3;2/9=110/2;6/7=2,8=2,9=2,10=2,19=2,28=1/1;99/9=1/99;,L101L103L202L301L302L303L401L502L601L701L702L703L716L9999,10.3.2Gaussian程序,2020/6/3,优化开始的标志,第一次做L103初猜Hessian矩阵,优化结束的标志,10.3.2Gaussian程序,2020/6/3,10.3.2Gaussian程序,2020/6/3,10.3.2Gaussian程序,2020/6/3,10.3.2Gaussian程序,2020/6/3,Hessian矩阵,10.3.2Gaussian程序,2020/6/3,10.3.2Gaussian程序,2020/6/3,10.3.2Gaussian程序,2020/6/3,优化时的坐标选项：OPT=RedundantOPT=AddRedundantOPT=Z-matrixOPT=Cartesian部分冻结优化：#HF/6-31G*Op