《中文版Exploring-Chemistry-with-Electronic-Structure-Methos-.docx》由会员分享,可在线阅读,更多相关《中文版Exploring-Chemistry-with-Electronic-Structure-Methos-.docx(66页珍藏版)》请在第一文库网上搜索。
1、Exp1oringChemistrywithE1ectronicStructureMethodSecondEdithionJamesB.ForesmanAe1eenFrischGaussian,IncPittsburgh,PA2002年9月25日特别声明本文转自南开大学BBS网站,在此对译者表示衷心感谢!!用GaiISSiai1研究化学问题说明接触GaUSSian已经很久了,但真正用GaUSSian做东西还是临近博士毕业时的事情。当时做计算的时候,就特别希望有一本具体怎么使用从头算的书,可惜一直没有找到。来到这里后,在新买的GaUSSian98包中发现了这本书,感觉如获至宝,也希望能够提供给想
2、用Gaussian做东西的朋友。我不是专门做量化的,很多术语不清楚怎么翻译,手头又没有中文的资料,错误的地方,只能希望内行来指点了。其实这本书里面介绍的东西,不止限于GaUSSian程序的。对于从事从头算研究的都有帮助。内容中有很多计算实例,都是在GaUSSian94,98程序中提供的。节译自Exp1oringChemistrywithE1ectronicStructureMethos,SecondEdition,作者JameSBoForesman,E1eenFrisch出版社GaUSSian,Inc,USA,1996特别声明1用Ga1ISSian研究化学问题1说明1前言1运行Gaussian
3、2UniX/1inux平台2WindoWS平台2输出文件2第一章计算模型41.1 计算化学概述4分子力学理论4电子结构理论4密度泛函(DenSityFunctiona1Methods)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 耦合簇(COUPIedC1USter)和二次结构相关(QUadratiCCOnfigUration1nteraCtiOn)方法.27密度泛函方法271.4 资源的使用29第七章高精度能量模型317.1 预测热化学31原子化能31电子亲和势31离子化能31质子亲和能317.2 理论模型的评价327.3 G2分子基(MO1eCU1eSet)以及缺陷及对缺陷的解释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理论背景;附录BGaUSSiaI1输入方法简介前言GaUSSian可以做很多事情,具体包括: 分子能量和结构研究 过渡态的能量和结构研究 化学键以及反应的能量 分子轨道 偶极矩和多极矩 原子电荷和电势 振动频率 红外和拉曼光谱 核磁 极化率和超极化率 热力学性质 反应途径计算可以模拟在气相和溶液中的体系,模拟
8、基态和激发态。Gaussian是研究诸如取代效应,反应机理,势能面和激发态能量的有力工具。运行GaussianUnix/1inux平台运行Gaussian前要设置好运行参数,比如在CShe11中,需要加这两句setenvg94rootdirectory/directory指程序的上级目录名source$g94root/g94/bsd/g94.1ogin然后运行就可以了。比如有输入文件,采用CShe11时的运行格式是:g94h2o.1og在UniX系统中,输入文件是.com为扩展名的,输出文件为.1og;Windows平台图形界面就不用多说了输入输出文件介绍在WindOWS系统中,输入文件是.g
9、jf为扩展名,输出文件为.out。下面是一个输入文件#TRHF/6-31G(d)TestMyfirstGaussianjob:watersing1epointenergy010-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-18Jou1esWc=2625.500kJmo1=27.2116eV;再后面是布居分析,有分子轨道情况,各个轨道的本征值(能量),各个原子的电荷,偶极距。然后是整个计算结果的一个总结,各小节之间用分开,所要的东西基本在里面了。然后是一句格言,随机有GaIISSian程序从它的格言库里选出的(在19999.exe中,想看的可以用文本格式打开这个文件,自己去找,学英语的好机会)。1然后是CPU时间,注意这不是真正的运行时间,是CPU运行的时间,真正的时间要长一些。如果几个工作一起做的话(WindOW下
12、好像不可能,UniX/1inix下可以同时做多个工作),实际计算时间就长很多了。最后一句话,Norma1terminationofGaussian94很关键,如果没有这句话,说明工作是失败的,肯定在什么地方出错误了。这是这里应该有出错信息。根据输入文件的设置,输出文件还要多一些内容,上面的是基本的东西。第一章计算模型计算化学的方法主要有分子力学理论(MO1eCU1arMeChaniCS)和电子结构理论(E1ectronicStructureTheory)。两者的共同点是: 计算分子的能量,分子的性质可以根据能量按照一定的方法得到。 进行几何优化,在起始结构的附近寻找具有最低的能量的结构。几何优
13、化是根据能量的一阶导数进行的。 计算分子内运动的频率。计算依据是能量的二阶导数。1.1 计算化学概述分子力学理论分子理论采用经典物理对分子进行处理,可以在MM3,HyperChem,Quanta,Syby1,AIChemy等软件中看到。根据所采用的力场的不同,分子理论又分为很多种。分子力学理论方法很便宜(做量化的经常用贵和便宜来描述计算,实际上就是计算时间的长短,因为对于要花钱上机的而言,时间就是金钱;对于自己有机器的,要想算的快,也要多在机器上花钱),可以计算多达几千个原子的体系。其缺点是每一系列参数都是针对特定原子得出的。没有对于原子各个状态的统一参数。计算中忽略了电子,只考虑键和原子,自
14、然就不能处理有很强电子效应的体系,比如不能描述键的断裂。电子结构理论这一理论基于薛定鄂方程,采用量子化学方法对分子进行处理。主要有两类:1 .半经验方法包括AMI,MIND03,PM3,常见的软件包有MOPAC,AMPAC,HyperChem,以及Gaussian0半经验方法采用了一些实验得来的参数,来帮助对薛定鄂方程的求解。2 .从头算从头算,在解薛定鄂方程的过程中,只采用了几个物理常数,包括光速,电子和核的质量,普朗克常数,在求解薛定鄂方程的过程中采用一系列的数学近似,不同的近似也就导致了不同的方法。最经典的是Hartree-FOCk方法,缩写为HF。从头算能够在很广泛的领域提供比较精确的信息,当然计算量要比前面讲的方法大的多,就是贵得多了。密度泛函(De1ISityFunctiona1Methods)密度泛函是最近几年兴起的第三类电子结构理论方法。它采用泛函(以函数为变量的函数)对薛定鄂方程进行求解,由于密度泛函包涵了电子相关,它的计算结果要比HF方法好,计算速度也快。1.2化学模型(MOde1Chemistries)Gaussian认为所谓理论是,一个理论模型,必须适用于任何种类和大小体系,它的应用限制只应该来自于计算条件的限制。这里包括两点: 一个理