《中文版Exploring-Chemistry-with-Electronic-Structure-Methos-.docx》由会员分享，可在线阅读，更多相关《中文版Exploring-Chemistry-with-Electronic-Structure-Methos-.docx（67页珍藏版）》请在课桌文档上搜索。
1、ExploringChemistrywithElectronicStructureMethodSecondEdithionJamesB.ForesmanAeleenFrischGaussian,IncPittsburgh,PA2002年9月25日特别声明本文转自南开大学BBS网站，在此对译者表示衷心感谢！!用GaiISSiail研究化学问题说明接触GaUSSian已经很久了，但真正用GaUSSian做东西还是临近博士毕业时的事情。当时做计算的时候，就特别希望有一本具体怎么使用从头算的书，可惜一直没有找到。来到这里后，在新买的GaUSSian98包中发现了这本书，感觉如获至宝，也希望能够提供给想
2、用Gaussian做东西的朋友。我不是专门做量化的，很多术语不清楚怎么翻译，手头又没有中文的资料，错误的地方，只能希望内行来指点了。其实这本书里面介绍的东西，不止限于GaUSSian程序的。对于从事从头算研究的都有帮助。内容中有很多计算实例，都是在GaUSSian94,98程序中提供的。节译自ExploringChemistrywithElectronicStructureMethos,SecondEdition,作者JameSBoForesman,EleenFrisch出版社GaUSSian,Inc,USA,1996特别声明1用GalISSian研究化学问题1说明1前言1运行Gaussian
3、2UniX/Linux平台2WindoWS平台2输出文件2第一章计算模型41.1 计算化学概述4分子力学理论4电子结构理论4密度泛函(DenSityFunctionalMethods)41.2 化学模型(MOdeIChemiStrieS)5定义化学模型5模型的组合5第二章单点能计算62.1 能量计算设置6路径6计算的名称6分子结构6多步计算72.3 输出文件中的信息7标准几何坐标。7能量7分子轨道和轨道能级7电荷分布7偶极矩和多极矩7CPU时间和其他72.4 核磁计算8第三章几何优化103.1 势能面103.2 寻找极小值10收敛标准10几何优化的输入H检查优化输出文件113.3 寻找过渡态1
4、13.4 难处理的优化11第四章频率分析144.1 预测红外和拉曼光谱14频率计算的输入14频率和强度14矫正因子和零点能。15简正模式15热力学16零点能(ZeroPointEnergy)和内能(ThermaIEnergy)17极化率和超极化率174.2 表征稳定点17下面列出了需要描述稳定点时必须考虑的问题18第五章基组的影响215.1 最小基组215.2 分裂基组215.3 极化基组215.4 弥散函数(DiffiISeFImCtionS)225.5 高角动量基组225.6 第三周期以后的原子的基组22第六章选择合适的理论模型256.1 使用半经验方法251.1 方法的局限性251.2
5、电子相关和后SCF方法26Hartree-Fock理论的限制26MPn方法261.3 耦合簇(COUPIedCLUSter)和二次结构相关(QUadratiCCOnfigUrationlnteraCtiOn)方法.27密度泛函方法271.4 资源的使用29第七章高精度能量模型317.1 预测热化学31原子化能31电子亲和势31离子化能31质子亲和能317.2 理论模型的评价327.3 G2分子基(MOleCUleSet)以及缺陷及对缺陷的解释327.4 理论模型的相对精确性327.5 组合方法33Gaussian-I和GaUSSian-2理论33完全基组方法(ComPIeteBasisSetM
6、otheds,CBS)34练习35第八章研究化学反应和反应性368.1 预测电子密度368.2 计算反应焰变368.3 研究势能面378.4 势能面扫描388.5 反应路径分析38运行IRC388.6 势能面研究实例39甲醛的解离39优化过渡态39频率分析39IRC计算40计算活化能401,2氢迁移反应40IRC的注意事项408.7 等构反应(ISodeSmiCReactions)41等构反应的局限41第九章激发态计算479.1 运行激发态计算479.2 激发态优化和频率分析48第十章溶液中的模型系统5510.1 反应场模型55Onsager模型的局限5510.2 运行SCRF计算55附录A5
7、9附录B60全书结构前言运行Gaussian第一部分基本概念和技术第一章计算模型；第二章单点能计算；第三章几何优化；第四章频率分析:第二部分计算化学方法第五章基族的影响；第六章理论方法的选择；第七章高精度计算；第三部分应用第八章研究反应和反应性；第九章激发态；第十章溶液中的反应；附录附录A理论背景；附录BGaUSSiaIl输入方法简介前言GaUSSian可以做很多事情，具体包括： 分子能量和结构研究 过渡态的能量和结构研究 化学键以及反应的能量 分子轨道 偶极矩和多极矩 原子电荷和电势 振动频率 红外和拉曼光谱 核磁 极化率和超极化率 热力学性质 反应途径计算可以模拟在气相和溶液中的体系，模拟
8、基态和激发态。Gaussian是研究诸如取代效应，反应机理，势能面和激发态能量的有力工具。运行GaussianUnix/Linux平台运行Gaussian前要设置好运行参数，比如在CShell中，需要加这两句setenvg94rootdirectory/directory指程序的上级目录名source$g94root/g94/bsd/g94.login然后运行就可以了。比如有输入文件,采用CShell时的运行格式是：g94h2o.log在UniX系统中，输入文件是.com为扩展名的，输出文件为.log；Windows平台图形界面就不用多说了输入输出文件介绍在WindOWS系统中，输入文件是.g
9、jf为扩展名，输出文件为.out。下面是一个输入文件#TRHF/6-31G(d)TestMyfirstGaussianjob：watersinglepointenergy010-0.4640.1770.0H-0.4641.1370.0H0.441-0.1430.0第一行以#开头，是运行的说明行，#1表示指打印重要的输出部分，#表示打印更多的信息。后面的RHF表示限制性Hartree-Fock方法，这里要输入计算所选用的理论方法。6-31G(d)是计算所采用的基组，就是用什么样的函数组合来描述轨道。Test是指不记入Gaussian工作档案，对于单机版没有用处。第三行是对于这个工作的描述，写什么
10、都行，自己看懂就是了。第二行是空行，这个空行以及第四行的空行都是必须的。第五行的两个数字，分别是分子电荷和自旋多重度。第六行以后是对于分子几何构性的描述。这个例子采用的是迪卡尔坐标。分子结构输入完成后要有一个空行。对于WindOWS版本，程序的图形界面把这几部分分得很清楚。输入时就不要添空行了。输出文件输出文件一般很长，对于上面的输入文件，其输出文件中，首先是版权说明，然后是作第2页共66页者，POPIe的名字在最后一个。然后是GaUSSian读入输入文件的说明，再将输入的分子坐标转换为标准内坐标，这些东西都不用去管。当然，验证自己的分子构性对不对就要看这个地方。关键的是有SCFDOIIe的一
11、行，后面的能量可是重要的，单位是原子单位，Hartree,1Hartree=4.3597482E-18JoulesWc=2625.500kJmol=27.2116eV；再后面是布居分析，有分子轨道情况，各个轨道的本征值（能量），各个原子的电荷，偶极距。然后是整个计算结果的一个总结，各小节之间用分开，所要的东西基本在里面了。然后是一句格言，随机有GaIISSian程序从它的格言库里选出的（在19999.exe中，想看的可以用文本格式打开这个文件，自己去找，学英语的好机会）。1然后是CPU时间，注意这不是真正的运行时间，是CPU运行的时间，真正的时间要长一些。如果几个工作一起做的话（WindOW下
12、好像不可能，UniX/Linix下可以同时做多个工作）,实际计算时间就长很多了。最后一句话，NormalterminationofGaussian94很关键，如果没有这句话，说明工作是失败的，肯定在什么地方出错误了。这是这里应该有出错信息。根据输入文件的设置，输出文件还要多一些内容，上面的是基本的东西。第一章计算模型计算化学的方法主要有分子力学理论(MOleCUlarMeChaniCS)和电子结构理论(ElectronicStructureTheory)。两者的共同点是： 计算分子的能量，分子的性质可以根据能量按照一定的方法得到。 进行几何优化，在起始结构的附近寻找具有最低的能量的结构。几何优
13、化是根据能量的一阶导数进行的。 计算分子内运动的频率。计算依据是能量的二阶导数。1.1 计算化学概述分子力学理论分子理论采用经典物理对分子进行处理，可以在MM3,HyperChem,Quanta,Sybyl,AIChemy等软件中看到。根据所采用的力场的不同，分子理论又分为很多种。分子力学理论方法很便宜(做量化的经常用贵和便宜来描述计算，实际上就是计算时间的长短，因为对于要花钱上机的而言，时间就是金钱；对于自己有机器的，要想算的快，也要多在机器上花钱)，可以计算多达几千个原子的体系。其缺点是每一系列参数都是针对特定原子得出的。没有对于原子各个状态的统一参数。计算中忽略了电子，只考虑键和原子，自
14、然就不能处理有很强电子效应的体系，比如不能描述键的断裂。电子结构理论这一理论基于薛定鄂方程，采用量子化学方法对分子进行处理。主要有两类：1 .半经验方法包括AMI,MIND03,PM3,常见的软件包有MOPAC,AMPAC,HyperChem,以及Gaussian0半经验方法采用了一些实验得来的参数，来帮助对薛定鄂方程的求解。2 .从头算从头算，在解薛定鄂方程的过程中，只采用了几个物理常数，包括光速，电子和核的质量，普朗克常数，在求解薛定鄂方程的过程中采用一系列的数学近似，不同的近似也就导致了不同的方法。最经典的是Hartree-FOCk方法，缩写为HF。从头算能够在很广泛的领域提供比较精确的
15、信息，当然计算量要比前面讲的方法大的多，就是贵得多了。密度泛函(DelISityFunctionalMethods)密度泛函是最近几年兴起的第三类电子结构理论方法。它采用泛函(以函数为变量的函数)对薛定鄂方程进行求解，由于密度泛函包涵了电子相关，它的计算结果要比HF方法好，计算速度也快。1.2化学模型(MOdelChemistries)Gaussian认为所谓理论是，一个理论模型，必须适用于任何种类和大小体系，它的应用限制只应该来自于计算条件的限制。这里包括两点： 一个理论模型应该对于任何给定的核和电子有唯一的定义，就是说，对于解薛定鄂方程来讲，分子结构本身就可以提供充分的信息。 一个理论模型
16、是没有偏见的，指不依靠于任何的化学结构和化学过程。这样的理论可以被认为是化学理论模型(theoretical-modelchemistry),简称化学模型(modelChemiStry)(这个翻译我可拿不准，在国内没听说过)。定义化学模型GaUSSian包含多种化学模型，比如计算方法Gaussian关键词方法HFHartree-Fock自恰场模型B3LYPBecke型3参数密度泛函模型，采用Lee-Yang-Parr泛函MP2二级Moller-Plesset微扰理论MP4四级Moller-Plesset微扰理论QCISD(T)二次Cl基组基组是分子轨道的数学表达，具体见第五章开壳层，闭壳层指电
17、子的自旋状态，对于闭壳层，采用限制性计算方法，在方法关键词前面加R；对于开壳层，采用非限制性计算方法，在方法关键词前面加U。比如开壳层的HF就是UHF。对于不加的，程序默认为是闭壳层。一般采用开壳层的可能性是 存在奇数个电子，如自由基，一些离子 激发态 有多个单电子的体系 描述键的分裂过程模型的组合高精度的计算往往要几种模型进行组合，比如用中等算法进行结构优化，然后用高精度算法计算能量。第二章单点能计算单点能计算是指对给定几何构性的分子的能量以及性质进行计算，由于分子的几何构型是固定不变的，只是“一个点”，所以叫单点能计算。单点能计算可以用于： 计算分子的基本信息 可以作为分子构型优化前对分子
18、的检查 在由较低等级计算得到的优化结果上进行高精度的计算 在计算条件下，体系只能进行单点计算 单点能的计算可以在不同理论等级，采用不同基组进行，本章的例子都采用HF方法和中等级基组。2.1 能量计算设置计算设置中，要有如下信息： 计算采用的理论等级和计算的种类 计算的名称 分子结构路径这里设置了计算要采用的理论方法，采用的基组，所要进行的计算的种类等信息。这一行，以#开头，默认的计算种类为单点能计算，关键词为SP,可以不写。这一部分需要出现的关键词有： 计算的理论，如HF（默认关键词，可以不写），B3PW91； 计算采用的基组，如631G,LANL2DZ； 布局分析方法，如Pop=Reg; 波
19、函数自恰方法，SCF=TightoPop=Reg只在输出文件中打印出最高的5条HOMO轨道和最低的5条LOMU轨道，而采用Pop=Full则打印出全部的分子轨道。SCF设置是指波函数的收敛计算时的设定，一般不用写，SCF=Tight设置表示采用比一般方法较严格的收敛计算。计算的名称一般含有一行，如果是多行，中间不能有空行。在这里描述所进行的计算。分子结构首先是电荷和自旋多重度电荷就是分子体系的电荷了，没有就是0,自旋多重度就是2S+1,其中S是体系的总自旋量子数，其实用单电子数加1就是了。没有单电子，自旋多重度就是1.然后是分子几何构型，一般可以用迪卡尔坐标，也可以用Z-矩阵(Z-MatriX
20、)多步计算GaUSSian支持多步计算，就是在一个输入文件中进行多个计算步骤。(见附录B)2.3 输出文件中的信息例2.1文件e2_01甲醛的单点能标准几何坐标。找到输出文件中StandardOrientation一行，下面的坐标值即输入分子的标准几何坐标。能量找到SCFDone：E(RHF)=-113.863697598A.U.after6cycles这里的数值就是能量，单位是hartree。在一些高等级计算中，往往有不止一个能量值，比如下一行：E2=-0.3029540001DOOEUMP2=-0.11416665769315D+03这里在EUMP2后面的数字是采用MP2计算的能量。MP4
21、计算的能量输出就更复杂了。分子轨道和轨道能级对于按照计算设置所打印出的分子轨道，列出的内容包括，轨道对称性以及电子占据情况，0表示占据，V表示空轨道；分子轨道的本征值，也就是分子轨道的能量，分子轨道的顺序就是按照能量由低到高的顺序排列的；每一个原子轨道对分子轨道的贡献。这里要注意轨道系数，这些数字的相对大小(忽略正负号)表示了组成分子轨道的原子轨道在所组成的分子轨道中的贡献大小。寻找HOMO和LUMO轨道的方法就是看占据轨道和非占据轨道的交界处。电荷分布GaUSSian采用的默认的电荷分布计算方法是MUnikin方法,在输出文件中寻找Totalatomiccharges,可以找到分子中所有原子
22、的电荷分布情况。偶极矩和多极矩GaSSian提供偶极矩和多极矩的计算，寻找DiPoIemomemt(Debye),下面就是偶极矩的信息，再下两行是四极矩。偶极矩的单位是德拜CPU时间和其他Jobcputime：Odays0hours0minuites9.1seconds0这里是计算的时间，注意是CPU时间。2.4 核磁计算例2.2文件e2_02甲烷的核磁计算核磁是单点能计算中另外一个可以提供的数据，在计算的工作设置部分，就是以#开头的一行里，加入NMR关键词就可以了，如：#TRHF6-31G(d)NMRTest在输出文件中，寻找如下信息GIAOMagneticshieldingtensor(p
23、pm)1CIsotropic=199.0522Anisotropy=0.0000这是采用上面的设置计算的甲烷的核磁结果，所采用的甲烷构形是用B3LYP密度泛函方法优化得到的。一般的，核磁数据是以TMS为零点的，下面是用同样的方法计算的TMS(四甲基硅烷)的结果：ICIsotropic=195.1196Anisotropy=17.5214这样，计算所得的甲烷的核磁共振数据就是3.9ppm,与实验值7.0PPm相比，还是很接近的。2.5练习练习2.1文件2_01丙烷的单点能练习要点：寻找分子的标准坐标，寻找单点能，偶极矩的方向和大小，电荷分布练习2.2文件2_02a(RR),2_02b(SS),2
24、_02c(RS)1,2.二氯1,2二氟乙烷的能量练习要点：比较该化合物三个旋光异构体的能量和偶极矩差异练习2.3文件2_03丙酮和甲醛的比较练习要点：比较甲基取代氢原子后带来的影响，说明能量比较必须在有同样的原子种类和数量的情况下进行练习2.4文件2_04乙烯和甲醛的分子轨道练习要点：寻找HOMo和LUMo能级，并分析能级的组成情况练习2.5文件2_05a,2_05b,2_05c烷，烯,快的核磁共振比较练习2.6文件2_06C60的单点能练习要点：分析C60最高占据轨道注意在收敛方法选择的时候，要有SCF=Tight,否则有收敛问题。练习2.7文件2_07计算大小的CPU资源比较本练习比较不同
25、基组函数数量，SCF方法对CPU时间，资源的占用情况。比较传统SCF方法(SCF=COnVem),直接SCF方法(GaUSSian默认方法)传统SCF直接SCF基组函数数量int文件大小(MB)CPU时间CPU时间2328.612.842411.919.8611623.238.8804248.772.1999295.4122.5118174163.4186.8137290354.5268.0156437526.5375.0175620740.2488.01948321028.4622.1很显然，基函数数量对资源占用和CPU时间都有很大影响，基函数越多，资源占用越大,CPU时间越长。理论上来讲，
26、认为CPU时间和基函数数量的四次方成正比，但实际上没有这么高，在本例中，基本上和基函数数量的2.5次方成正比。一般的讲，直接SCF方法的效率要比传统SCF方法要好，在本例中，当基函数数量比较大时，可以看到这一点。练习2.8文件2J)8a(O2),2_08b(03)SCF稳定性计算本例中采用SCF方法分析分子的稳定性。对于未知的体系，SCF稳定性是必须要做的。当分子本身不稳定的时候，所得到的SCF结果以及波函数等信息就没有化学意义。SCF稳定性分析是寻找是不是存在比当前状态能量更低的分子状态。关键词有：Stable检验分子的稳定性，放松对分子的限制，比如由闭壳层改为开壳层等。Stable=OPT
27、这一选项设定，当发现不稳定的时候，对新的状态进行优化。这种做法一般是不推荐的，因为所得到的新的状态的几何形太接近原来的几何构形。本例中首先计算闭壳层的单重态的氧分子。很显然，闭壳层单重态的氧分子不应该是稳定的。在输出文件中，我们可以找到这样的句子：TheWavefuctionhasanRHFUHFinstabilityUHF的不稳定性，在所得到的UHF状态下进行稳定性检验，采用UHFStable=Opt,发现体系仍然不稳定。Thewavefunctionhasaninernalinstability再在此基础上进行的优化，体系又回到了RHF的状态。这时，就需要在进行SCF前的构性初始电子状态猜
28、测上进行改动，使用Guess=Mix,在初始猜测中混合HOMO和LUMO轨道，从而消除空间对称性，然后进行的UHFGUeSS=MiXStabIe表明得到了稳定的结构。确定电子状态还可以采用Guess=Alter详见GaussianUsersReference第三章几何优化前面讨论了在特定几何构型下的能量的计算，可以看出，分子几何构型的变化对能量有很大的影响。由于分子几何构型而产生的能量的变化，被称为势能面。势能面是连接几何构型和能量的数学关系。对于双原子分子，能量的变化与两原子间的距离相关，这样得到势能曲线，对于大的体系，势能面是多维的，其维数取决与分子的自由度。3.1 势能面势能面中，包括一
29、些重要的点，包括全局最大值，局域极大值，全局最小值，局域极小值以及鞍点。极大值是一个区域内的能量最高点，向任何方向的几何变化都能够引起能量的减小。在所有的局域极大值中的最大值，就是全局最大值；极小值也同样，在所有极小之中最小的一个就是具有最稳定几何结构的一点。鞍点则是在一个方向上具有极大值，而在其他方向上具有极小值的点。一般的，鞍点代表连接着两个极小值的过渡态。3.2 寻找极小值几何优化做的工作就是寻找极小值，而这个极小值，就是分子的稳定的几何形态。对于所有的极小值和鞍点，其能量的一阶导数，也就是梯度，都是零，这样的点被称为稳定点。所有的成功的优化都在寻找稳定点，虽然找到的并不一定就是所预期的
30、点。几何优化由初始构型开始，计算能量和梯度，然后决定下一步的方向和步长，其方向总是向能量下降最快的方向进行。大多数的优化也计算能量的二阶导数，来修正力矩阵，从而表明在该点的曲度。收敛标准当一阶导数为零的时候优化结束，但实际计算上，当变化很小，小于某个量的时候，就可以认为得到优化结构。对于GaUSSian,默认的条件是： 力的最大值必须小于0.045 均方根小于0.0003 为下一步所做的取代计算为小于0.0018 其均方根小于0.0012这四个条件必须同时满足，比如，对于非常松弛的体系，势能面很平缓，力的值已经小于域值，但优化过程仍然有很长的路要走。对于非常松弛的体系，当力的值已经低于域值两个
31、数量级，尽管取代计算仍然高于域值，系统也认为找到了最优点。这条规则用于非常大，非常松弛的体系。几何优化的输入Opt关键字描述了几何优化例3.1文件e3_01乙烯的优化输入文件的设置行为 RRHF6-31G(d)OptTest表明采用RHF方法，631G(d)基组进行优化检查优化输出文件优化部分的计算包含在两行相同的：GradGradGradGradGradGradGradGradGradGrad之间，这里有优化的次数，变量的变化，收敛的结果等等。注意这里面的长度单位是波尔。在得到每一个新的几何构型之后，都要计算单点能，然后再在此基础上继续进行优化，直到四个条件都得到满足。而最后一个几何构型就被
32、认为是最优构型。注意，最终构型的能量是在最后一次优化计算之前得到的。在得到最优构型之后，在文件中寻找-Stationmaypointfoundo其下面的表格中列出的就是最后的优化结果以及分子坐标。随后按照设置行的要求，列出分子有关性质例3.2文件e3_02氟代乙烯的优化3.3 寻找过渡态Gaissian使用STQN方法确定反应过渡态，关键词是Opt=QST2例3.3文件e3_03过渡态优化例中分析的是H3CO-H2COH的变化，输入文件格式#TUHF6-31G(d)0pt=QST2TestH3C0-H2C0HReactants0,2structureforH3C00,2structurefor
33、H2COHGaUSSian也提供QST3方法，可以优化反应物，产物和一个由用户定义猜测的过渡态。3.4 难处理的优化有一些系统的优化很难进行，采用默认的方法得不到结果，其产生的原因往往是所计算出的力矩阵与实际的相差太远。当默认方法得不到结果时，就要采用其他的方法。Gaussian提供很多的选择，具体可以看UsersReferenceo下面列举一些。Opt=ReadFC从频率分析(往往是采用低等级的计算得到的)所得到的checkpoint文件中读取初始力矩阵，这一选项需要在设置行之前加入Chk=filename一句，说明文件的名称。Opt=CaICFC采用优化方法同样的基组来计算力矩阵的初始值。
34、Opt=CalcAlI在优化的每一步都计算力矩阵。这是非常昂贵的计算方法，只在非常极端的条件下使用。有时候，优化往往只需要更多的次数就可达到好的结果，这可以通过设置MaxCycle来实现。如果在优化中保存了Checkpoint文件，那么使用OPt=ReStart可以继续所进行的优化。当优化没有达到效果的时候，不要盲目的加大优化次数。这是注意观察每一步优化的区别，寻找没有得到优化结果的原因，判断体系是否收敛，如果体系能量有越来越小的趋势，那么增加优化次数是可能得到结果的，如果体系能量变化没有什么规律，或者，离最小点越来越远，那么就要改变优化的方法。也可以从输出文件的某一个中间构型开始新的优化，关
35、键词GeOm=(Check,SteP=II)表示在取得在CheCkPOint文件中第n步优化的几何构型3.5 练习练习3.1文件3_01a(180),3_01b(0)丙烯的优化从两种丙烯的几何异构体进行优化,一个是甲基的一个氢原子与CCH形成180度二面角，另一个是0.优化结果表明，二者有0.003Hartree的差别，0度的要低。练习3.2文件3_02a(0),3_02b(180),3_02c(actealdo)乙烯醇的优化乙烯醇氧端的氢原子与OCC平面的二面角可以为0和180,优化得到的结果时，0度的能量比180度的低0.003Hartree,但同时做的乙醛的优化表明，乙醛的能量还要低，比
36、0度异构体低0.027hartreeo练习3.3文件3_03乙烯胺的优化运行所有原子都在同一平面上的乙烯胺的优化。比较本章的例子和练习，可以看到不同取代基对乙烯碳碳双键的影响。练习3.4文件3_04六段基铭的优化本例采用STO-3G和3-21G基组，在设置行中加入SCF=NoVarAcc对收敛有帮助。3-2IG基组的优化结果要优于STO-3G练习3.5文件3_05a(C6H6),3_05b(TMS)苯的核磁共振采用6-31G(d)基组，B3LYP方法优化几何构性，采用HF方法，6-311+G(2d,P)基组在优化的几何构型基础上计算碳的化学位移。注意，核磁共振的可靠程度依赖准确的几何结构和大的
37、基组。输入文件如下%Chk=NMR#TB3LYP/6-31G(d)OptTestOptmoleculespecification-Linkl-%Chk=NMR%NoSave#TRHF/6-311+G(2d,p)NMRGeom=CheckGuess=ReadTestNMRcharg&spin同样，还需要采用同样方法计算TMS。下面是计算结果绝对位移相对位移实验值TMSBenzene188.787957.6198131.2130.9练习3.6文件3_06a(PM3),3_06b(STO-3G)氧化碳60的优化C60中有两种碳碳键，一是连接两个六元环的6-6键，另一是连接六元环和无元环的5-6键。氧
38、化C60就有两种异构体。本例采用PM3和HF/STO-3G方法来判断那种异构体是稳定的，以及氧化后的C-C键的变化。采用OPt=AddRedlmdant关键词可以在输出文件中打印所要求的键长，键角，这一关键词需要在分子构型输入结束后在增加关于所要键长键角的信息，键长用两个原子的序列号表示，键角则用三个原子表示。计算结果显示，6-6键的氧化，碳碳键仍然存在，接近环氧化合物，而5-6键已经打开。采用不同的方法，得到的几何结构相差不多，但在能量上有很大差异。在采用MNDO,PM3,HF/3-21G方法得到的能量数据中，5-6键氧化的异构体的能量低,但采用HF/STO-3G得到的结果，确实6-6键氧化
39、的能量低。Raghavachari在其进行的上述研究中阐述 动力学因素同样是重要的； 实验上还没有发现那个是能量最低的异构体； 应该进行更精确的计算练习3.7文件3_07一个1,1消除反应的过渡态优化分析反应SiH4-SiH2+H2,可以采用OPt=(QST2,AddRedImdalIt)关键词来进行过渡态优化，同时特别关注过渡态结构中的某个键长练习3.8文件3_08优化进程比较采用下述三种方法优化二环2,2,2 直接采用默认方式冗余内坐标优化OPt 采用迪卡尔坐标优化OPt=CarteSiaIl 采用内坐标优化C)Pt=Z-MatriX结果显示，冗余内坐标优化的优化次数最短，内坐标优化的次数
40、最多。第四章频率分析频率分析可以用于多种目的： 预测分子的红外和拉曼光谱(频率和强度) 为几何优化计算力矩阵 判断分子在势能面上的位置 计算零点能和热力学数据如系统的燃和焰4.1 预测红外和拉曼光谱几何优化和单点能计算都将原子理想化了，实际上原子一直处于振动状态。在平衡态，这些振动是规则的和可以预测的。频率分析的计算要采用能量对原子位置的二阶导数。HF方法，密度泛函方法(如B3LYP),二阶Moller-Plesset方法(MP2)和CASSCF方法(CASSCF)都可以提供解析二阶导数。对于其他方法，可以提供数值二阶导数。频率计算的输入Freq关键词代表频率分析。频率分析只能在势能面的稳定点
41、进行，这样，频率分析就必须在已经优化好的结构上进行。最直接的办法就是在设置行同时设置几何优化和频率分析。特别注意的是，频率分析计算是所采用的基组和理论方法，必须与得到该几何构型采用的方法完全相同。例4.1文件e4_01甲醛的频率分析例中采用的是已经优化好的几何构型，输入格式#RHF6-31G(d)FreqTest频率和强度频率分析首先要计算输入结构的能量，然后计算频率。Gaussian提供每个振动模式的频率，强度，拉曼极化率。以下是例4.1的输出文件中的前四个频率1234BlB2AlAlFrequencies1336.00411383.64491679.58432028.0971Red.mas
42、ses-1.36891.34421.10397.2497IRIntensities0.369423.15898.6240150.1861RamanActivities0.76574.517012.85948.1124Depolarizations0.75000.75000.59080.3281由于HF方法忽略的电子相关，所得到的频率结果有系统误差，大约在10%-12%一般的，对于HF方法，采用计算的频率乘以矫正因子0.8929,从而与实验值吻合。本文的计算由于采用的是中等基组，误差要稍大一些，大约15%o计算得到的相对强度数值是可靠的。矫正因子和零点能。其他方法得到的频率同样存在系统误差，下面
43、列出对于频率和零点能的矫正因子方法频率矫正因子零点能矫正因子HF/3-21G0.90850.9409HF6-31G(d)0.89290.9135MP2(Full)6-31G(d)0.94270.9646MP2(FC)6-31G(d)0.94340.9676SVWN6-31G(d)0.98331.0079BLYP6-31G(d)0.99401.0119B3LYP6-31G(d)0.96130.9804频率的矫正因子和用于计算热力学数据的零点能矫正因子之间有小的差异，但一般处理上，可以采用同样的因子，就是频率的矫正因子。基组也对矫正因子产生影响，相关的数据就要差相应的参考书了。简正模式输出文件同时
44、还提供振动的简正模式。比如对于如下的标准坐标（例4.1）Standardorientation：CenterAtomicCoordinates(Angstroms)NumberNumberXYZ160.0000000.000000-0.542500280.0000000.0000000.677500310.0000000.940000-1.082500410.000000-0.940000-1.082500其第一个振动模式是AtomANXYZ160.1728-0.0431-0.7041-0.70在标准坐标中，0.000.000.000.000.000.000.000.00四个原子的X坐标都是零，分子在YZ平面上。对于振动模式，其熟知的正负比大小要重要的多。在本例中，两个氢原子向X轴负方向有显著的振动，表明其偏离了X平面。这个振动产生的红外光谱在1189Cm-I（经过矫正）热力学频率分析也包括对体系的热力学分析。默认情况下，系统计算在298.15K和Iatm下的热力学数值。下面是输出的计算热力学的参数ThermochemiStryTempeature298.15KelvinoPressure1.00000
