Gaussian软件使用方法基组1

合集下载

第三章Gaussian软件的使用-new

%chk=h2o2.chk %rwf=h2o2.rwf #p hf/6-31g opt H2O2 energy calculation 01 H O 1 B1 O 2 B2 1 A1 H 3 B1 2 A1 1 D1 B1=0.9 A1=105.0 D1=120.0 B2=1.4
变量部分
常量部分
通过设定常量，可以使程序只对分子结构进行部分优化
%chk=h2o2.chk %rwf=h2o2.rwf #p hf/6-31g opt H2O2 energy calculation 01 H O 1 B1 O 2 B2 1 A1 H 3 B1 2 A1 1 D1 B1=0.9 B2=1.4 A1=105.0 D1=120.0
通过对两个键长和键角使用同一变量定义可以控制分子的对称性
Title Section
Title部分必须输入，但是程序并不执行，起标识和说明作用
Charge & Multipl.
输入分子的电荷和自旋多重度例：
H2O H3O+ · NO
电荷 0 1 0 多重度2s+1 1 1 2
电荷多重度部分通常也算作分子说明
Molecular Specification
局域密度局域密度近似近似
价电子从头算价电子从头算
模拟从头算模拟从头算
自洽场自洽场从头算从头算 SCFSCFab initio ab initio
分子碎片法分子碎片法
浮动球浮动球高斯法从头算法高斯法
从头算法 Ab Initio Ab Initio
梯度近似梯度近似
组态组态耦合相互作用耦合相互作用微扰处理电子对电子对微扰处理 CI 独立 CI 多组态独立多组态 Hartree-Fock 电子对自洽场 Hartree-Fock 电子对自洽场方程方程

Gaussian计算软件的使用

.
11
接CIS优化
• # CIS(restart)/3-21g* guess=read geom=check POPT TEST
• 或者 • # CIS/3-21g* POPT=restart TEST
.12Biblioteka TD光谱计算• #P B3LYP/6-31G* TD(singlets,Nstates=30, root=1) TEST
• P0 • 6-31+g(d) • **** • Au 0 • CEP-121G • **** • • Au 0 • CEP-121G
.
8
ZINDO光谱计算
• ZINDO • #P ZINDO=(singlets, Nstates=30, root=1)
GFINPUT IOP(6/7=3) TEST
.
15
G03中的Counterpoise
• #p b3lyp/gen pseudo=read Counterpoise=2 GFINPUT IOP(6/7=3) test
• • 2005-4-29 6:51 • • 020302 • bq 0.0 0.0 0.0 0 • bq 1 xxxx2 0 • c 2 cxx3 1 90.000 2 • c 2 cxx4 1 90.000 3 dih4 0 2 • c 2 cxx4 1 90.000 3 -dih4 0 2 • o 2 oxx6 1 90.000 3 dih6 0 2 • o 2 oxx6 1 90.000 3 -dih6 0 2 • ti 2 tixx8 1 tixxxx8 3 dih8 0 1
命令行：定义中间文件名和资源限制
关键词行：指定工作类型、模型化学和选项
＊空行
＊说明行：计算的简要描述

Gaussian 高斯使用指南

和 Parrinello 的经验，ADMP 传递电子自由度，而不是求解每个核结构的 SCF 方程。与
Car-Parrinello 不同之处在于，ADMP 传递密度矩阵而不是 MO。如果使用了原子中心基组，
5.模拟在反应和分子特性中溶剂的影响
在气相和在溶液之间，分子特性和化学反应经常变化很大。例如，低位构像在气相和在（不
同溶剂的）溶液中，具有完全不同的能量，构像的平衡结构也不同，化学反应具有不同的路
径。Gaussian 03 提供极化连续介质模型(PCM)，用于模拟溶液体系。这个方法把溶剂描述为
极化的连续介质，并把溶质放入溶剂间的空穴中。
2．高斯自带练习 ……………………………………………………………71
第八章选择合适的理论模型 …………………………………………………72
1 简要介绍 …………………………………………………………………72
2．高斯自带的练习 …………………………………………………………73
第九章高精度能量模型 ………………………………………………………76
Gaussian 03 的 PCM 功能包含了许多重大的改进，扩展了研究问题的范围：可以计算溶剂中
的激发能，以及激发态的有关特性；NMR 以及其它的磁性能；用能量的解析二级导数计算
振动频率，IR 和 Raman 光谱，以及其它特性；极化率和超极划率；执行性能上的改善。
G03W 的界面和 G98W 相比，没有什么变化，G98W 的用户不需要重新熟悉界面。
1.2.2 输出的解释 ………………………………………………………61
2.高斯中自带的练习 …………………………………………………………64
第六章：单点能计算 ……………………………………………………………64

Gaussian的使用(一)

可使用，程序能自动识别。可包含多个输入行，段后需加一空白行。可选的附加段空白行。下面以 H2O 分子在（从头算）RHF/STO-3G 水平上的单点能量计算为例说明各输入段的内容与格式。用 Notepad 编辑的作业输入文件 H2O-3G.gjf 的内容为 %Chk=H2O-3G 检查点文件名为 H2O-3G.chk 根据某些作业的特殊需要而输入。段后一般也应加一
# RHF/STO-3G 指定理论模型（计算方法）和基组
Water Energy 作业的题目
0 O H 0.
1
闭壳层中性分子，自旋多重度为 1 0. 0. 原子笛卡儿坐标的 x、y、z 分量
0.784 -0.554 0.
H -0.784 -0.554 0.
可在元素符号后加上原子的序列编号，也可用原子序数（核电荷数）代替元素符号。下面两种形式的原子标注在 G98W 亦适用： %Chk=H2O-3G # RHF/STO-3G
主要计算项目包括： • 分子的能量与几何结构 • 化学反应过渡态的能量与几何结构 • 振动频率分析 • 红外与拉曼光谱 • 分子的热化学性质 • 键能与化学反应能 • 化学反应途径 • 分子轨道的能量与性质 • 原子电荷分布（电子布居分析）和自旋密度分布
• 分子的多极矩（永久偶极矩和四极至十六极矩） • NMR 屏蔽常数、化学位移及分子的磁化率
Gaussian 软件的使用
第一部分几何结构的优化和能量的精确计算
Gaussian03 软件是目前国际上最为流行的量子化学计算软件，可执行各类不同精度和理论档次的 MO 计算，包括 Hartree-Fock 水平从头算（HF）、Post-HF 从头算（各级 CI 和 MP）、MC-SCF 法、密度泛函理论（DFT），以及多种半经验量子化学方法，通过计算可以预测分子和化学反应的许多性质，具有计算精度高、适用范围广及容易操作的特点。

Gaussian的使用(一)

主要计算项目包括： • 分子的能量与几何结构 • 化学反应过渡态的能量与几何结构 • 振动频率分析 • 红外与拉曼光谱 • 分子的热化学性质 • 键能与化学反应能 • 化学反应途径 • 分子轨道的能量与性质 • 原子电荷分布（电子布居分析）和自旋密度分布
• 分子的多极矩（永久偶极矩和四极至十六极矩） • NMR 屏蔽常数、化学位移及分子的磁化率
用默认的外部浏览器显示分子结构的 3D 图形
用外部文本编辑器打开输出文件，阅读、分析和编辑计算结果
2.在“作业编辑对话框”内进行在线编辑打开软件，进入 Gaussian03 主窗口。在主窗口中打开下拉菜单“File”, 点击“New”选项后将跳出“作业编辑对话框” （如图 1 所示），框内工具条上各按钮的功能为
图 1 中已标出每一段的主要作用和功能。
图 2 乙炔分子编辑对话框
其中 OPT 表示优化分子结构，HF 表示采用从头算 Hartree-Fock 方法，6 31G * 表示采用的基函数，0 和 1 分别表示乙烯分子的电荷和自旋多重度。分子说明部分采用的是分子的内坐标，它主要是通过分子中原子的空间构型进行设置，具体格式为：元素标记，参考原子，键长，参考原子，键角，参考原子，二面角
包含多个输入行，每行的行首均用#号开始。段后需加一空白行。题目段作业计算内容与目的的简要描述，以便于输出文件的阅读。
可包含多个输入行，段后需加一空白行。电荷与自旋行给定分子所带的净电荷数（正负整数）及自旋多重度
（正整数）。限一行，行后不加空白行。分子说明段给定分子中各原子的坐标。笛卡儿坐标或分子内坐标均
Water Energy
0 O1 H2 0.
1 0. 0.
0.ห้องสมุดไป่ตู้84 -0.554 0.

gaussian教程

gaussian教程节译自Exploring Chemistry with Electronic Structure Methos, Second Edition,作者James B. Foresman, Eleen Frisch 出版社Gaussian, Inc, USA, 1996前言Gaussian可以做很多事情,具体包括分子能量和结构研究过渡态的能量和结构研究化学键以及反应的能量分子轨道偶极矩和多极矩原子电荷和电势振动频率红外和拉曼光谱核磁极化率和超极化率热力学性质反应途径计算可以模拟在气相和溶液中的体系,模拟基态和激发态.Gaussian是研究诸如取代效应,反应机理,势能面和激发态能量的有力工具.全书结构序言运行Gaussian第一部分基本概念和技术第一章计算模型第二章单点能计算第三章几何优化第四章频率分析第二部分计算化学方法第五章基族的影响第六章理论方法的选择第七章高精度计算第三部分应用第八章研究反应和反应性第九章激发态第十章溶液中的反应附录A 理论背景附录B Gaussian输入方法简介运行GaussianUnix/Linux平台:运行gaussian前要设置好运行参数,比如在C Shell中,需要加这两句setenv g94root directory / directory指程序的上级目录名source $g94root/g94/bsd/g94.login然后运行就可以了.比如有输入文件,采用C Shell时的运行格式是g94 h2o.logWindows平台:图形界面就不用多说了输入输出文件介绍在Unix系统中,输入文件是.com为扩展名的,输出文件为.log;在Windows系统中,输入文件是.gjf为扩展名,输出文件为.out.下面是一个输入文件#T RHF/6-31G(d) TestMy first Gaussian job: water single point energy0 1O -0.464 0.177 0.0H -0.464 1.137 0.0H 0.441 -0.143 0.0第一行以#开头,是运行的说明行,#T表示指打印重要的输出部分,#P表示打印更多的信息.后面的RHF表示限制性Hartree-Fock方法,这里要输入计算所选用的理论方法6-31G(d)是计算所采用的基组,就是用什么样的函数组合来描述轨道Test是指不记入Gaussian工作档案,对于单机版没有用处.第三行是对于这个工作的描述,写什么都行,自己看懂就是了.第二行是空行,这个空行以及第四行的空行都是必须的.第五行的两个数字,分别是分子电荷和自旋多重度.第六行以后是对于分子几何构性的描述.这个例子采用的是迪卡尔坐标.分子结构输入完成后要有一个空行.对于Windows版本,程序的图形界面把这几部分分得很清楚.输入的时候就不要再添空行了.输出文件输出文件一般很长,对于上面的输入文件,其输出文件中,首先是版权说明,然后是作者,Pople的名字在最后一个.然后是Gaussian读入输入文件的说明,再将输入的分子坐标转换为标准内坐标,这些东西都不用去管.当然,验证自己的分子构性对不对就要看这个地方.关键的是有SCF Done的一行,后面的能量可是重要的,单位是原子单位,Hartree,1 Hartree= 4.3597482E-18 Joules或=2625.500 kJ/mol=27.2116 eV再后面是布居分析,有分子轨道情况,各个轨道的本征值(能量),各个原子的电荷,偶极距.然后是整个计算结果的一个总结,各小节之间用\分开,所要的东西基本在里面了.然后是一句格言,随机有Gaussian程序从它的格言库里选出的(在l9999.exe中,想看的可以用文本格式打开这个文件,自己去找,学英语的好机会).然后是CPU时间,注意这不是真正的运行时间,是CPU运行的时间,真正的时间要长一些.如果几个工作一起做的话(Window下好像不可能,Unix/Linix下可以同时做多个工作),实际计算时间就长很多了.最后一句话,"Normal termination of Gaussian 94"很关键,如果没有这句话,说明工作是失败的,肯定在什么地方出错误了.这是这里应该有出错信息.根据输入文件的设置,输出文件还要多一些内容,上面的是基本的东西.第一章计算模型1.1 计算化学的方法主要有分子理论(Molecular Mechanics)和电子结构理论(Electronic Structure Theory).两者的共同点是1. 计算分子的能量,分子的性质可以根据能量按照一定的方法得到.2. 进行几何优化,在起始结构的附近寻找具有最低的能量的结构.几何优化是根据能量的一阶导数进行的.3. 计算分子内运动的频率.计算依据是能量的二阶导数.1.2 分子理论分子理论采用经典物理对分子进行处理,可以在MM3,HyperChem, Quanta, Sybyl, Alchemy等软件中看到.根据所采用的力场的不同,分子理论又分为很多种.分子理论方法很便宜(做量化的经常用贵和便宜来描述计算,实际上就是计算时间的长短,因为对于要花钱上机的而言,时间就是金钱;对于自己有机器的,要想算的快,也要多在机器上花钱),可以计算多达几千个原子的体系.其缺点是1. 每一系列参数都是针对特定原子得出的.没有对于原子各个状态的统一参数.2. 计算中忽略了电子,只考虑键和原子,自然就不能处理有很强电子效应的体系, 比如不能描述键的断裂.1.3 电子结构理论这一理论基于薛定鄂方程,采用量子化学方法对分子进行处理.主要有两类:1. 半经验方法,包括AM1, MINDO/3, PM3,常见的软件包有MOPAC, AMPAC, HyperChem, 以及Gaussian.半经验方法采用了一些实验得来的参数,来帮助对薛定鄂方程的求解.2. 从头算.从头算,在解薛定鄂方程的过程中,只采用了几个物理常数,包括光速,电子和核的质量,普朗克常数,在求解薛定鄂方程的过程中采用一系列的数学近似,不同的近似也就导致了不同的方法.最经典的是Hartree-Fock方法,缩写为HF.从头算能够在很广泛的领域提供比较精确的信息,当然计算量要比前面讲的方法大的多,就是贵得多了.1.4 密度泛函(Density Functional Methods)密度泛函是最近几年兴起的第三类电子结构理论方法.它采用泛函(以函数为变量的函数)对薛定鄂方程进行求解,由于密度泛函包涵了电子相关,它的计算结果要比HF方法好,计算速度也快.1.5 化学模型(Model Chemistries)Gaussian认为所谓理论是,一个理论模型,必须适用于任何种类和大小体系,它的应用限制只应该来自于计算这里包括两点,1. 一个理论模型应该对于任何给定的核和电子有唯一的定义,就是说,对于解薛定鄂方程来讲,分子结构本身就可以提供充分的信息.2. 一个理论模型是没有偏见的,指不依靠于任何的化学结构和化学过程.这样的理论可以被认为是化学理论模型(theoretical-model chemistry),简称化学模型(model chemistry)(这个翻译我可拿不准,在国内没听说过).1.6 定义化学模型Gaussian包含多种化学模型,比如计算方法Gaussian关键词方法HF Hartree-Fock自恰场模型B3L YP Becke型3参数密度泛函模型,采用Lee-Yang-Parr泛函MP2 二级Moller-Plesset微扰理论MP4 四级Moller-Plesset微扰理论QCISD(T) 二次CI具体在第六章讨论基组基组是分子轨道的数学表达,具体见第五章开壳层,闭壳层指电子的自旋状态,对于闭壳层,采用限制性计算方法,在方法关键词前面加R对于开壳层,采用非限制性计算方法,在方法关键词前面加U.比如开壳层的HF就是UHF.对于不加的,程序默认为是闭壳层.一般采用开壳层的可能性是1. 存在奇数个电子,如自由基,一些离子2. 激发态3. 有多个单电子的体系4. 描述键的分裂过程模型的组合高精度的计算往往要几种模型进行组合,比如用中等算法进行结构优化,然后用高精度算法计算能量.第二章单点能计算2.1 单点能计算是指对给定几何构性的分子的能量以及性质进行计算,由于分子的几何构型是固定不变的,只是"一个点",所以叫单点能计算.单点能计算可以用于:计算分子的基本信息可以作为分子构型优化前对分子的检查在由较低等级计算得到的优化结果上进行高精度的计算在计算条件下,体系只能进行单点计算单点能的计算可以在不同理论等级,采用不同基组进行,本章的例子都采用HF方法2.2 计算设置计算设置中,要有如下信息:计算采用的理论等级和计算的种类计算的名称分子结构方法设置这里设置了计算要采用的理论方法,采用的基组,所要进行的计算的种类等信息.这一行,以#开头,默认的计算种类为单点能计算,关键词为SP,可以不写.这一部分需要出现的关键词有,计算的理论,如HF(默认关键词,可以不写),B3PW91;计算采用的基组,如6-31G, Lanl2DZ;布局分析方法,如Pop=Reg;波函数自恰方法,如SCF=Tight.Pop=Reg只在输出文件中打印出最高的5条HOMO轨道和最低的5条LOMU轨道,而采用Pop=Full则打印出全部的分子轨道.SCF设置是指波函数的收敛计算时的设定,一般不用写,SCF=Tight设置表示采用比一般方法较严格的收敛计算.计算的名称一般含有一行,如果是多行,中间不能有空行.在这里描述所进行的计算.分子结构首先是电荷和自旋多重度电荷就是分子体系的电荷了,没有就是0,自旋多重度就是2S+1,其中S是体系的总自旋量子数,其实用单电子数加1就是了.没有单电子,自旋多重度就是1.然后是分子几何构性,一般可以用迪卡尔坐标,也可以用Z-矩阵(Z-Matrix)多步计算Gaussian支持多步计算,就是在一个输入文件中进行多个计算步骤.2.3 输出文件中的信息例2.1 文件e2_01 甲醛的单点能标准几何坐标.找到输出文件中Standard Orientation一行,下面的坐标值就是输入分子的标准几何坐标.能量找到SCF Done: E(RHF)= -113.863697598 A. U. after 6 cycles这里的树脂就是能量,单位是hartree.在一些高等级计算中,往往有不止一个能量值,比如下一行E2=-0.3029540001D+00 EUMP2=-0.11416665769315D+03这里在EUMP2后面的数字是采用MP2计算后的能量.MP4计算的能量输出就更复杂了分子轨道和轨道能级对于按照计算设置所打印出的分子轨道,列出的内容包括,轨道对称性以及电子占据情况,O表示占据,V表示空轨道;分子轨道的本征值,也就是分子轨道的能量,分子轨道的顺序就是按照能量由低到高的顺序排列的;每一个原子轨道对分子轨道的贡献.这里要注意轨道系数,这些数字的相对大小(忽略正负号)表示了组成分子轨道的原子轨道在所组成的分子轨道中的贡献大小.寻找HOMO和LUMO轨道的方法就是看占据轨道和非占据轨道的交界处.电荷分布Gaussian采用的默认的电荷分布计算方法是Mullikin方法,在输出文件中寻找Total atomic charges可以找到分子中所有原子的电荷分布情况.偶极矩和多极矩Gassian提供偶极矩和多极矩的计算,寻找Dipole momemt (Debye),下面就是偶极矩的信息,再下两行是四极矩偶极矩的单位是德拜CPU时间和其他Job cpu time : 0days 0 hours 0 minuites 9.1 seconds.这里是计算的时间,注意是CPU时间.2.4 核磁计算例2.2 文件e2_02 甲烷的核磁计算核磁是单点能计算中另外一个可以提供的数据,在计算的工作设置部分,就是以#开头的一行里,加入NMR关键词就可以了,如#T RHF/6-31G(d) NMR Test在输出文件中,寻找如下信息GIAO Magnetic shielding tensor (ppm)1 C Isotropic = 199.0522 Anisotropy = 0.0000这是采用上面的设置计算的甲烷的核磁结果,所采用的甲烷构形是用B3L YP密度泛函方法优化得到的.一般的,核磁数据是以TMS为零点的,下面是用同样的方法计算的TMS(四甲基硅烷)的结果1 C Isotropic = 195.1196 Anisotropy = 17.5214这样,计算所得的甲烷的核磁共振数据就是-3.9ppm,与实验值-7.0ppm相比,还是很接近的.2.5 练习练习2.1 文件2_01 丙烷的单点能练习要点:寻找分子的标准坐标,寻找单点能,偶极矩的方向和大小,电荷分布练习2.2 文件2_02a (RR), 2_02b (SS), 2_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_06 C60的单点能练习要点:分析C60最高占据轨道注意在收敛方法选择的时候,要有SCF=Tight,否则有收敛问题.练习2.7 文件2_07 计算大小的CPU资源比较本练习比较不同基组函数数量,SCF方法对CPU时间,资源的占用情况.比较传统SCF方法(SCF=Convern),直接SCF方法(Gaussian默认方法)传统SCF 直接SCF基组函数数量int文件大小(MB) CPU时间CPU时间23 2 8.6 12.842 4 11.9 19.861 16 23.2 38.880 42 48.7 72.199 92 95.4 122.5118 174 163.4 186.8137 290 354.5 268.0156 437 526.5 375.0175 620 740.2 488.0194 832 1028.4 622.1很显然,函数数量对资源占用和CPU时间都有很大影响,函数越多,资源占用越大,CPU时间越长.理论上来讲,认为CPU时间和函数数量的四次方成正比,但实际上没有这么高, 在本例中,基本上和函数数量的2.5次方成正比.一般的讲,直接SCF方法的效率要比传统SCF方法要好,在本例中,当函数数量比较大时, 可以看到这一点.练习2.8 文件2_08a (O2), 2_08b (O3) SCF稳定性计算本例中采用SCF方法分析分子的稳定性.对于未知的体系,SCF稳定性是必须要做的.当分子本身不稳定的时候,所得到的SCF结果以及波函数等信息就没有化学意义.SCF稳定性分析是寻找是不是存在比当前状态能量更低的分子状态.关键词有Stable 检验分子的稳定性,放松对分子的限制,比如由闭壳层改为开壳层等.Stable=OPT 这一选项设定,当发现不稳定的时候,对新的状态进行优化.这种做法一般是不推荐的,因为所得到的新的状态的几何形太接近原来的几何构形.本例中首先计算闭壳层的单重态的氧分子.很显然,闭壳层单重态的氧分子不应该是稳定的.在输出文件中,我们可以找到这样的句子:The wavefuction has an RHF --> UHF instability.这表明存在一个UHF的状态,其能量要比当前状态低.这说明可能,能量最低的状态是单重态的,但不是闭壳层的;存在有更低能量的三重态;所计算的状态不是能量最低点,可能是过渡态.在三重态情况下重新计算,也进行稳定性验证,可以看到如下的句子The wavefunction is stable under the perturbations considered.臭氧是单重态的,但有不一般的电子结构.采用RHF Stable=Opt可以发现一个RHF-->UHF的不稳定性,在所得到的UHF状态下进行稳定性检验,采用UHF Stable=Opt,发现体系仍然不稳定.The wavefunction has an inernal instability再在此基础上进行的优化,体系又回到了RHF的状态.这时,就需要在进行SCF前的构性初始电子状态猜测上进行改动,使用Guess=Mix,在初始猜测中混合HOMO 和LUMO轨道,从而消除空间对称性,然后进行的UHF Guess=Mix Stable 表明得到了稳定的结构.确定电子状态还可以采用Guess=Alter详见Gaussian User's Reference第三章几何优化前面讨论了在特定几何构型下的能量的计算,可以看出,分子几何构型的变化对能量有很大的影响.由于分子几何构型而产生的能量的变化,被称为势能面.势能面是连接几何构型和能量的数学关系.对于双原子分子,能量的变化与两原子间的距离相关,这样得到势能曲线,对于大的体系,势能面是多维的,其维数取决与分子的自由度.3.1势能面势能面中,包括一些重要的点,包括全局最大值,局域极大值,全局最小值,局域极小值以及鞍点.极大值是一个区域内的能量最高点,向任何方向的几何变化都能够引起能量的减小.在所有的局域极大值中的最大值,就是全局最大值;极小值也同样,在所有极小之中最小的一个就是具有最稳定几何结构的一点.鞍点则是在一个方向上具有极大值,而在其他方向上具有极小值的点.一般的,鞍点代表连接着两个极小值的过渡态.寻找极小值几何优化做的工作就是寻找极小值,而这个极小值,就是分子的稳定的几何形态.对于所有的极小值和鞍点,其能量的一阶导数,也就是梯度,都是零,这样的点被称为稳定点.所有的成功的优化都在寻找稳定点,虽然找到的并不一定就是所预期的点. 几何优化有初始构型开始,计算能量和梯度,然后决定下一步的方向和步长,其方向总是向能量下降最快的方向进行.大多数的优化也计算能量的二阶导数,来修正力矩阵,从而表明在该点的曲度.3.2 收敛标准当一阶导数为零的时候优化结束,但实际计算上,当变化很小,小于某个量的时候,就可以认为得到优化结构.对于Gaussian,默认的条件是力的最大值必须小于0.00045,均方根小于0.0003为下一步所做的取代计算为小于0.0018,其均方根小于0.0012这四个条件必须同时满足,比如,对于非常松弛的体系,势能面很平缓,力的值已经小于域值,但优化过程仍然有很长的路要走.对于非常松弛的体系,当力的值已经低于域值两个数量级,尽管取代计算仍然高于域值,系统也认为找到了最优点.这条规则用于非常大,非常松弛的体系.3.3 几何优化的输入Opt关键字描述了几何优化例3.1 文件e3_01 乙烯的优化输入文件的设置行为#R RHF/6-31G(d) Opt Test表明采用RHF方法,6-31G(d)基组进行优化3.4 输出文件优化部分的计算包含在两行相同的GradGradGradGradGradGradGradGradGradGrad...........之间,这里有优化的次数,变量的变化,收敛的结果等等.注意这里面的长度单位是波尔.在得到每一个新的几何构型之后,都要计算单点能,然后再在此基础上继续进行优化,直到四个条件都得到满足.而最后一个几何构型就被认为是最优构型.注意,最终构型的能量是在最后一次优化计算之前得到的.在得到最优构型之后,在文件中寻找--Stationmay point found.其下面的表格中列出的就是最后的优化结果以及分子坐标.随后按照设置行的要求,列出分子有关性质例3.2 文件e3_02 氟代乙烯的优化3.5 寻找过渡态Gaissian使用STQN方法确定反应过渡态,关键词是Opt=QST2例3.3 文件e3_03 过渡态优化例中分析的是H3CO --> H2COH 的变化,输入文件格式#T UHF/6-31G(d) Opt=QST2 TestH3CO --> H2COH Reactants0,2structure for H3CO0,2structure for H2COHGaussian也提供QST3方法,可以优化反应物,产物和一个由用户定义的猜测的过渡态.3.6 难处例的优化有一些系统的优化很难进行,采用默认的方法得不到结果,其产生的原因往往是所计算出的力矩阵与实际的相差太远.当默认方法得不到结果时,就要采用其他的方法. Gaussian提供很多的选择,具体可以看User's Reference.下面列举一些.Opt=ReadFC 从频率分析(往往是采用低等级的计算得到的)所得到的checkpoint文件中读取初始力矩阵,这一选项需要在设置行之前加入%Chk= filename 一句,说明文件的名称.Opt=CalCFC 采用优化方法同样的基组来计算力矩阵的初始值.Opt=CalcAll 在优化的每一步都计算力矩阵.这是非常昂贵的计算方法,只在非常极端的条件下使用.有时候,优化往往只需要更多的次数就可以达到好的结果,这可以通过设置MaxCycle 来实现.如果在优化中保存了Checkpoint文件,那么使用Opt=Restart可以继续所进行的优化.当优化没有达到效果的时候,不要盲目的加大优化次数.这是注意观察每一步优化的区别,寻找没有得到优化结果的原因,判断体系是否收敛,如果体系能量有越来越小的趋势,那么增加优化次数是可能得到结果的,如果体系能量变化没有什么规律,或者,离最小点越来越远,那么就要改变优化的方法.也可以从输出文件的某一个中间构型开始新的优化,关键词Geom=(Check,Step=n)表示在取得在Checkpoint文件中第n步优化的几何构型3.7 练习练习3.1 文件3_01a (180), 3_01b (0) 丙烯的优化从两种丙烯的几何异构体进行优化,一个是甲基的一个氢原子与CCH形成180度二面角,另一个是0.优化结果表明,二者有0.003Hartree的差别,0度的要低.练习3.2 文件3_02a (0), 3_02b (180), 3_02c (acteald.) 乙烯醇的优化乙烯醇氧端的氢原子与OCC平面的二面角可以为0和180,优化得到的结果时,0度的能量比180度的低0.003Hartree,但同时做的乙醛的优化表明,乙醛的能量还要低,比0度异构体低0.027hartree.练习3.3 文件3_03 乙烯胺的优化运行所有原子都在同一平面上的乙烯胺的优化.比较本章的例子和练习,可以看到不同取代基对乙烯碳碳双键的影响.练习3.4 文件3_04 六羰基铬的优化本例采用STO-3G和3-21G基组,在设置行中加入SCF=NoVarAcc对收敛有帮助.3-21G基组的优化结果要优于STO-3G练习3.5 文件3_05a (C6H6), 3_05b (TMS) 苯的核磁共振采用6-31G(d)基组,B3L YP方法优化几何构性,采用HF方法,6-311+G(2d,p)基组在优化的几何构型基础上计算碳的化学位移.注意,核磁共振的可靠程度依赖准确的几何结构和大的基组.输入文件如下%Chk=NMR#T B3L YP/6-31G(d) Opt TestOptmolecule specification--Link1--%Chk=NMR%NoSave#T RHF/6-311+G(2d,p) NMR Geom=Check Guess=Read TestNMRcharg & spin同样,还需要采用同样方法计算TMS.下面是计算结果绝对位移相对位移实验值TMS Benzene188.7879 57.6198 131.2 130.9练习3.6 文件3_06a (PM3), 3_06b (STO-3G) 氧化碳60的优化C60中有两种碳碳键,一是连接两个六元环的6-6键,另一是连接六元环和无元环的5-6键. 氧化C60就有两种异构体.本例采用PM3和HF/STO-3G方法来判断那种异构体是稳定的,以及氧化后的C-C键的变化.采用Opt=AddRedundant关键词可以在输出文件中打印所要求的键长,键角,这一关键词需要在分子构型输入结束后在增加关于所要键长键角的信息,键长用两个原子的序列号表示,键角则用三个原子表示.计算结果显示,6-6键的氧化,碳碳键仍然存在,接近环氧化合物,而5-6键已经打开.采用不同的方法,得到的几何结构相差不多,但在能量上有很大差异.在采用MNDO,PM3,HF/3-21G方法得到的能量数据中,5-6键氧化的异构体的能量低,但采用HF/STO-3G得到的结果,确实6-6键氧化的能量低.Raghavachari在其进行的上述研究中阐述动力学因素同样是重要的;实验上还没有发现那个是能量最低的异构体;应该进行更精确的计算练习3.7 文件3_07 一个1,1消除反应的过渡态优化分析反应SiH4 --> SiH2 + H2, 可以采用Opt=(QST2, AddRedundant)关键词来进行过渡态优化,同时特别关注过渡态结构中的某个键长练习3.8 文件3_08 优化进程比较采用下述三种方法优化二环[2,2,2]直接采用默认方式冗余内坐标优化Opt;采用迪卡尔坐标优化Opt=Cartesian;采用内坐标优化Opt=Z-Matrix结果显示,冗余内坐标优化的优化次数最短,内坐标优化的次数最多.第四章频率分析频率分析可以用于多种目的,预测分子的红外和拉曼光谱(频率和强度)为几何优化计算力矩阵判断分子在势能面上的位置计算零点能和热力学数据如系统的熵和焓4.1 红外和拉曼光谱几何优化和单点能计算都将原子理想化了,实际上原子一直处于振动状态.在平衡态,这些振动是规则的和可以预测的.频率分析的计算要采用能量对原子位置的二阶导数.HF方法,密度泛函方法(如B3L YP), 二阶Moller-Plesset方法(MP2)和CASSCF方法(CASSCF)都可以提供解析二阶导数.对于其他方法,可以提供数值二阶导数.4.2 频率分析输入Freq关键词代表频率分析.频率分析只能在势能面的稳定点进行,这样,频率分析就必须在已经优化好的结构上进行.最直接的办法就是在设置行同时设置几何优化和频率分析.特别注意的是,频率分析计算是所采用的基组和理论方法,必须与得到该几何构型采用的方法完全相同!例4.1 文件e4_01 甲醛的频率分析例中采用的是已经优化好的几何构型,输入格式# RHF/6-31G(d) Freq Test4.3 频率和强度频率分析首先要计算输入结构的能量,然后计算频率.Gaussian提供每个振动模式的频率,强度,拉曼极化率.以下是例4.1的输出文件中的前四个频率1 2 3 4B1 B2 A1 A1。

Gaussian_03_使用(1)

因此，它和体系的总电子自旋S有关：
自旋多重度=2S+1
1个和2个单电子时自旋多重度的确定
容易知道，闭壳层体系的自旋多重度为1。具有一个单电子的体系自旋多重度为2。余下类推。成对电子对自旋多重度没有贡献。
几个自旋多重度判断的例子
1. 一、二周期原子的基态的自旋多重度判断：
原子：
H, He, Li， Be, B, C, N,
基本原则：根据研究体系，选择公认合理的基组 1. 配位化合物中心原子应该加d型极化函数 2. 相关能计算应该用大的基组，包括较多的极化函数 3. 对于阴离子等对电子束缚较弱的体系，应该添加弥
散函数 4. 对于弱相互作用体系，如氢键，弱分子间配合物，
范德华相互作用等，应当添加弥散函数 5. 计算激发态，应该使用大的基组，添加弥散函数
n
STO diGTOi
原子轨道分层 i1
1s 2s2 p 3s3 p 内层价层
内层和价层同样采用单zeta拟合方式，则为最小基组STONG
内层采用单zeta，价层采用多zeta，则为劈裂基组 3-21G，6-31G，6-311G等
3-21G为例
Gaussian03内含的劈裂基组
• 3-21G • 6-21G • 4-31G • 6-31G • 6-311G • D95 • D95V
均匀的，没有考虑瞬时的电子相关。因此准确度不高。
MP2理论
• 体系的哈密顿由两部分组成：
H H0 V
H0 是体系无微扰时可以精确求解的哈密顿；V 是加于H0 的，很小的微扰项；
由于微扰的存在，体系波函数也要加上一些微扰项：
(0) (1) 2(2) 3(3) .....
体系能量也要加上一些微扰项：

Gaussian软件使用方法技巧1

经验
这个方法是计算中最常用到的方法。它只需要你提供在反应坐标上位于TS两边的结构就可以了，程序会拟合两个结构来计算过渡态。一般人说最好是反应物和产物，这里我不认同。我认为初始的两个species的反应坐标为TS+n和TS-n，且n越小越好，这样拟合的结果更为接近TS的结构。一般先要确定两个species一定是在反应坐标上位于TS两侧的。这一点很简单，只需要自由优化一下，看结构式跑向R还是跑向P就可以了。如果都跑向R/P，那么就改变一下反应的主坐标参数，这个确认步骤很简单，只需要跑几个点。确认完两个结构现在反应坐标上位于TS两侧后，分别固定两个species的反应主坐标参数进行优化，目的是为了让除了主坐标参数以外的其他部分都达到一个较为合理的结构，便于qst2/qst3的进行。这一步优化不用优化到完全收敛，只需要在能量曲线平坦之后即可取出结构做qst2/qst3，节省时间，通常来说20步足以。
H
2.47445600 0.12780100 0.00402400
C
0.66866000 1.30146700 -0.16992700
H
1.26438900 2.14276400 -0.50621500
C
-0.66890700 1.30134300 -0.16992700
H
-1.26479300 2.14252900 -0.50621700
BCOT minus
01
15 19 1.8 F
15和19号原子将其距离固定为1.8A后优化
随后将优化好的结构算过渡态 opt（calcfc , ts ,noeigen）
如果还没有找到的话，可以对你猜测的构象进行频
率分析，然后查看有没有你所要的虚频，如果有的话用GV vibrations 最下面的manual displacement 调节震动，选中你需要的虚频向左调0.2-0.4或者你不要的虚频向右调0.2-0.4以微调构象，然后接着套用 opt(ts,calcfc,noeigen)继续找过渡态。

gaussian 基组

gaussian 基组Gaussian基组是一种重要的分子计算方法，它可以用来计算各种分子的电子结构和性质。

在本文中，我们将介绍Gaussian基组的基本概念和用法。

第一步，理解基组的概念基组是分子计算中用来描述电子结构和性质的一组数学基函数的称谓。

这些函数可以被用来近似地描述原子轨道和化学键的特征。

Gaussian基组是一种常用的基组，其中，高斯型函数被用来描述电子波函数。

高斯型函数是一种常用的数学函数，它具有一个峰值和一个标准差。

第二步，学习Gaussian基组的类型Gaussian基组可以分为不同的类型，其中最常见的类型是STO和GTO。

STO（Slater-type orbitals）是一个Slater波函数的结果，它是由一个序列的指数和一组系数定义的。

STO-3G和STO-6G是两个常见的STO基组，它们分别包含3个和6个原子轨道。

GTO（Gaussian-type orbitals）基组是由高斯型基函数组成的，这种基组比STO基组更易于处理和计算。

如6-31G和6-31G*是两个常见的GTO基组，它们更为完备。

第三步，应用Gaussian基组进行计算Gaussian基组可以被用来计算各种分子的电子结构和性质。

在使用Gaussian基组计算分子时，应选用适当的基组，确保计算结果的准确度和可靠性。

对于小分子的计算，通常可以使用STO-3G或STO-6G 基组；对于大分子的计算，最好选择GTO基组，例如6-31G或6-31G*，可能需要包括多个极化和扩展基组。

第四步，理解Gaussian基组的相关问题虽然Gaussian基组是一种非常有用的计算工具，但是它也存在一些局限。

首先，Gaussian基组只是基函数的一种选择，不同的基组可能会影响到计算结果的准确度和可靠性。

其次，Gaussian基组计算的速度不如其他分子计算方法，因此在处理大分子的计算时，需要更高效的计算方法。

最后，在使用Gaussian基组进行分子计算时，需要合理的计算策略和参数选择。

Gaussian计算软件的使用

geom=check TEST
TD相关关键词
• Singlets 只算单重态 • Triplets 只算三重态 • 50-50 单重和三重态各占一半 • Root=N 研究TD计算的第几个激发态ect
加快计算速度，减小硬盘使用空间 • Read 从Chk文件中读 TD （G03有，好
像G98 A11.2也有） • SOS Do sum-over states polarizabilities,
Gaussian 的启动
工作环境初始化设置
G98W 安装后的首次启动时需进行工作环境初始化设置，规定主程序、检查点文件、输入输出文件等默认子目录和路径。其内容记录在G98W 根目录下的初始化文件g98w.ini 上。用户可以在启动G98W 前用Notepad 对该文件直接编辑修改，也可点击“File”下拉菜单中的“Preference”选项在程序界面上进行设定。下面介绍后一种初始化方式。在主窗口上打开下拉菜单File
计算速度，减小硬盘使用空间 • Restart 从Chk文件中重新开始 CIS 计算 • Read 从Chk文件中读 CIS • 分析激发态电荷情况用 • Density=current 或者Density=CIS
接CIS优化
• # CIS(restart)/3-21g* guess=read geom=check POPT TEST
root=1) TEST • CIS 优化 • #P CIS(singlets,Nstates=30, root=1)/6-
31G* POPT FREQ TEST
CIS相关关键词
• Singlets 只算单重态 • Triplets 只算三重态 • 50-50 单重和三重态各占一半 • Root=N 研究CIS计算的第几个激发态Direct 加快

第四章-Gaussian软件的使用PPT课件

• Chk 临时文件 • Rwf 读写文件 • Int 双电子积分 • D2e 双电子积分
导数文件
• PDB， ENT • MOL， HIN • CDX， C3D
2021/6/4
3
运行界面
2021/6/4
4
准备知识
• Schrodinger方程 • Born-Oppenheimer近似 • Hellmann-Feyman定理 • 维里定理
信息 (包括执行时间数据，以及SCF计算的收敛信息) ✓ #T 精简输出：只打印重要的信息和结果。
➢ Route Section主要由方法，基组，任务类型三部分组成
✓ 方法与基组后续课程专门介绍 ✓ Gaussian程序能完成的任务类型：
Sp
单点能量计算（默认任务类型）
Opt
分子几何结构优化
Freq 振动分析
，采用每个原子的两种不同大小的函数的线
性组合来描述分子轨道。同样的，三重分裂
基组，如6-311G，采用三个不同大小的收缩
2021/6/4
函数来描述轨道。
9
极化基组
• 分裂基组允许轨道改变其大小，但不能改变形状。极化基组则取消了这样的限制，增加了角动量。
• 比如在碳原子上增加d 轨道的成分，在过渡金属上增加f 轨道成分。有些在氢原子上增加p 轨道成分
MP2梯度 4 MW 6 MW 16 MW 38 MW
MP2振动分析 6 MW 10 MW 28 MW
j ≈60 MW
≈70 MW
2021/6/4
28
Route Section
➢ Route Section以# 开始，# 控制作业的输出
✓ #N 正常输出；默认 (没有计算时间的信息) ✓ #P 输出更多信息。包括每一执行模块在开始和结束时与计算机系统有关的各种

第四章 Gaussian软件的使用

内存使用控制
%mem=n
控制运行过程中使用内存的大小,可以以W或者MB，GB为单位 default：6000000W＝48MB 综合考虑到计算的需要和硬件水平，内存并非给得越多越好，最有效率的方法是根据作业类型估算所需要内存的大小
不同作业使用内存的估算方法
M + 2NB2
M：不同类型作业需要的最小内存， NB ：计算所使用基函数的数目作业类型 SCF能量 SCF梯度 SCF振动分析 MP2能量 MP2梯度 MP2振动分析
Title Section
Title部分必须输入，但是程序并不执行，起标识和说明作用
Charge & Multipl.
输入分子的电荷和自旋多重度例：
H2O H3O+ · NO O2
电荷 0 1 0 0 多重度2s+1 1 1 2 2
电荷多重度部分通常也算作分子说明
Molecular Specification
f g 4 MW 4 MW 4 MW 5 MW 4 MW 9 MW 4 MW 5 MW 4 MW 6 MW 6 MW 10 MW
h 9 MW 16 MW 27 MW 10 MW 16 MW 28 MW
i 23 MW 38 MW 28 MW 38 MW
j ≈60 MW ≈70 MW
Route Section
.chk文件在计算中记录分子几何构型，分子轨道，力常数矩阵等信息
%rwf=name.rwf
.rwf文件主要在作业重起时使用，当计算量比较大时，.rwf文件通常会非常大，此时需要将之分割保存
%int=name.int, %d2e=name.d2e
.int文件在计算过程中存储双电子积分, .d2e文件在计算过程中存储双电子积分的二阶导数

第四章 Gaussian软件的使用

内存使用控制
%mem=n
控制运行过程中使用内存的大小,可以以W或者MB，GB为单位 default：6000000W＝48MB 综合考虑到计算的需要和硬件水平，内存并非给得越多越好，最有效率的方法是根据作业类型估算所需要内存的大小
不同作业使用内存的估算方法
M + 2NB2
M：不同类型作业需要的最小内存， NB ：计算所使用基函数的数目作业类型 SCF能量 SCF梯度 SCF振动分析 MP2能量 MP2梯度 MP2振动分析
高角动量基组
• 6-31G(2d)就是在6-31G 基础上增加两个 d 轨道的函数，而6-311++G(3df，3pd) 则增加了更多的极化函数，包括三个分裂的价键基组，在重原子和氢原子上加的弥散函数，在重原子上加的三个d 函数和一个f 函数，在氢原子上加的三个p 函数和一个d 函数。这样的基组在电子相关方法对于描述电子之间的作用有很重要意义。
第4周期以后的原子的基组
• 第4周期以上的原子的基组很难处理。由于存在非常大的核，原子核附近的电子通过有效核电势方法 (ECP)进行了近似，这一处理同时也包含了相对论效应。这其中， LANL2DZ 是最有名的基组。
基组的选择
L0
Gaussian
Overlay0
L001 L101 L102 L122
Water ennergy
01 O H 1 R1 H 1 R1 2 a1 R1=1.04 a1=104.0
title：作业的简要描述，段后加空行
Molecular Specification：
分子说明部分，段后通常加空行
% Section(link 0)
定义计算过程中的临时文件

Gaussian程序使用入门

内坐标表示
笛卡尔直角坐标表示
上述两个作业将对H2O2分子的结构进行完全优化，包括所有的键长键角和二面角。通过将内坐标定义成变量，可以对分子结构进行部分优化
%chk=h2o2.chk %rwf=h2o2.rwf #p hf/6-31g opt
H2O2 energy calculation
01 H O 1 r1 O 2 r2 1 a1 H 3 r1 2 a1 1 d1
% Section(link 0)
➢ 定义计算过程中的临时文件
%chk=name.chk
.chk文件在计算中记录分子几何构型，分子轨道，力常数矩阵等信息
%rwf=name.rwf
.rwf文件主要在作业重起时使用，当计算量比较大时，.rwf文件通常会非常大，此时需要将之分割保存
%int=name.int, %d2e=name.d2e
r1=0.9 r2=1.4 a1=105.0 d1=120.0
通过对两个键长和键角使用同一变量定义可以控制分子的对称性
%chk=h2o2.chk %rwf=h2o2.rwf #p hf/6-31g opt
H2O2 energy calculation
01 H O 1 r1 O 2 r2 1 a1 H 3 r1 2 a1 1 d1
.int文件在计算过程中存储双电子积分, .d2e文件在计算过程中存储双电子积分的二阶导数
➢ 内存使用控制
%mem=n
控制运行过程中使用内存的大小,可以以W或者MB，GB为单位 default：6000000W＝48MB 综合考虑到计算的需要和硬件水平，内存并非给得越多越好，最有效率的方法是根据作业类型估算所需要内存的大小
Preferences: 对Gaussian程序进行初始化设置

Gaussian 高斯使用指南

和 Parrinello 的经验，ADMP 传递电子自由度，而不是求解每个核结构的 SCF 方程。与
Car-Parrinello 不同之处在于，ADMP 传递密度矩阵而不是 MO。如果使用了原子中心基组，
据提供了分子中各原子之间的连接信息。自旋-自旋耦合常数可用来帮助识别分子的特定构
像，因为它们依赖于分子结构的扭转角。
除了以前版本提供的 NMR 屏蔽和化学位移以外，Gaussian 03 还能预测自旋-自旋耦合常数。
通过对不同构像计算这些常数，并对预测的和观测的光谱做比较，可以识别观测到的特定构
像。另外，归属观测的峰值到特定的原子也比较容易。
2.2 安装软件………………………………………………………………8
2.3 Gaussian的补钉和升级………………………………………………9
2.4 定制软件运行时占用内存和硬盘的最大空间………………………9
2.5 工作环境初始化设置…………………………………………………10
第三章高斯输入文件的创建和程序运行………………………………………11
5.模拟在反应和分子特性中溶剂的影响
在气相和在溶液之间，分子特性和化学反应经常变化很大。例如，低位构像在气相和在（不
同溶剂的）溶液中，具有完全不同的能量，构像的平衡结构也不同，化学反应具有不同的路
径。Gaussian 03 提供极化连续介质模型(PCM)，用于模拟溶液体系。这个方法把溶剂描述为
极化的连续介质，并把溶质放入溶剂间的空穴中。
提供对一般分子力场(MM)的支持，包括读入和修改参数。包含了独立的 MM 优化程序。
支持任何 ONIOM 模拟的外部程序。
(2) 修改和增强了溶剂模块
改善和增强了连续介质模型(PCM)：

gaussian基本概念和用法

1 计算流程上面的图应该是从L301说起的，即图中的“基组”，L301的作用是产生基组信息。

（在这之前还有L1：处理计算执行路径，创建执行链接的列表，并初始化scratch文件；L101：读取标题和分子说明部分；L103：berny优化到最小值；L202：重新定位坐标，计算对称性，检查变量）“基组”之后是“分子结构”即L401，形成初始的分子轨道初猜。

第一个框中有“分子轨道”和“能量”两项，我想这应该指的是L502即迭代求解SCF方程。

下面红线圈的框中则是迭代求解的具体过程。

每次迭代完毕都得到一个能量。

即（SCF Done: E(RHF) = -xx.xxxx A.U. after xx cycles）然后是L103判断力和位移是否满足收敛条件。

如果满足（即四个yes），则进行布居分析等，并完成计算。

如果不满足，则继续调整分子结构（L401），再次进行迭代求解，反复这个过程，直到满足力和位移收敛条件为止。

如果不是HF方法而是多体微扰，CI等后-scf方法的话，在每次迭代（L502）完成后还要多一步相关能的计算，以MP2为例，迭代完成后多出:L801(双电子积分变换的初始化),L906(半直接的MP2)和L1002(迭代求解CPHF方程)等计算。

对于楼主的图，我觉得力和位移判断为NO以后，箭头应该指向最头上的“基组”，因为实际过程中，每次构型调整都会重新定位坐标和产生基组信息。

常见问题分析1．检查是否有初始文件错误在命令行中加入 %kJob L301 or %kJob L302 如果通过则一般初始文件ok。

常见初级错误：a. 自旋多重度错误; b. 变量赋值为整数; c. 变量没有赋值或多重赋值; d. 键角小于等于0度，大于等于180度e. 分子描述后面没有空行; f. 二面角判断错误，造成两个原子距离过近; g. 分子描述一行内两次参考同一原子,或参考原子共线2． SCF(自洽场)不收敛则一般是L502错误,省却情况做64个cycle迭代(G03缺省128 cycles)a. 修改坐标，使之合理;b. 改变初始猜 Guess=Huckel 或其他的，看Guess关键词;c. 增加叠代次数 SCFCYC=N （对小分子作计算时最好不要增加，很可能结构不合理）; d. iop（5/13=1）这样忽略不收敛，继续往下做。

Gaussian软件使用方法基组1

关键词: GEN和PSEUDO=READ(对赝势基组)
Mo 0 Lanl2mb **** CO0 6-31G* **** Mo 0 Lanl2mb
对于比较轻的元素，比如CHON之类，b3lyp/6-311++g**就很优秀，有人甚至说，Cl以前的元素用6-31g就好，如果有孤对电子加个d就行了，如果带负电荷就加个弥散。
•对于使用较多的lanl系列基组，需注意：对于H~Ne范围的原子，lanl1mb和lan2mb基组等价于 sto-3g最小基组；
lanl1dz和lanl2dz基组等价于d95基组;
对于含金属—金属键体系一般使用lanl2系列，而不用 lanl1系列；
3) 混合基组的使用
例如, 对于Mo(CO)2+化合物, C和O原子可采用6-31+G基组, 但 Mo原子不在6-31G基组的使用范围之内.
极化(polarization)函数的使用
由于极化轨道的使用增强价轨道的空间柔软程度，因此其使用
场合主要在于环状化合物，例如有机环状化合物以及含有桥联
结构的金属簇合物等。
使用方法:在标准基组后加“*”，“**”，(d)，(d,p)， (3df,3pd)等例：6-31G*, 6-31G(ndf,mpd)
对于一些中等重量的元素，LanL2DZ或者DZVP之类的基组就很不错，没必要加更高的机组了。Fe的LanL2DZ有人换成 SDD试过，换了以后得到不同基态结果，所以千万小心使用 SDD。第一行过渡金属，如果只有一个，那么6-311+g*的计算量还是可以承受的；两个及两个以上，非赝势很难算。 lanl2dz和sdd都是比较好的选择。除非做单点能计算，否则非要上mp2算就没太大必要了。
到的能量要较HF来得精确，它是目前最常用的量子化学计算方法。

Gaussian软件应用——基组的影响

Gauss‎ian软件‎应用——基组的影响‎第五章基组的影响‎基组是体系‎内轨道的数‎学描述.大的基组由‎于对电子在‎空间上有小‎的限制而具‎有更大的精确‎性.用于电子结‎构计算的标‎准的基组使‎用线性的高‎斯函数来模‎拟轨道.Gauss‎ian提供‎大量的已经‎定义好的基‎组.5.1 最小基组最小基组包‎含了描述轨‎道的最少的‎函数数量.H: 1sC: 1s, 2s, 2px, 2py, 2pzSTO-3G是最小‎基组(虽然不是可‎能的最小基‎组),每一个基本‎函数中含有‎三个高斯函‎数,于是就有了‎3G的名称‎.S TO代表‎S late‎r形的轨道‎,这样,STO-3G就表示‎采用三个高斯函数‎来描述Sl‎a ter轨‎道.5.2 分裂基组增大基组的‎第一个方法‎就是增加每‎个原子基函‎数的数量.分裂基组,比如3-21G和6-31G,对于价键轨‎道都用两个‎函数来进行‎描述,比如H: 1s, 1s'C: 1s, 2s, 2s', 2px, spy, spz, spx', spy', spz'其中的主要‎轨道和非主‎要轨道在大‎小上不同.双zeta‎基组,如Dunn‎ing-Huzin‎aga基组‎(D95),采用每个原‎子的两种不‎同大小的函‎数的线性组合‎来描述分子‎轨道.同样的,三重分裂基‎组,如6-311G,采用三个不‎同大小的收‎缩函数来描‎述轨道.5.3 极化基组分裂基组允‎许轨道改变‎其大小,但不能改变‎形状.极化基组则‎取消了这样‎的限制,增加了角动‎量.比如在碳原‎子上增加d‎轨道的成分‎,在过渡金属‎上增加f轨‎道成分.有些在氢原‎子上增加p‎轨道成分.一般的,常用的极化‎基组是6-31G(d),这个基组来‎源与6-31G基组‎,并在其基础‎上,对于重原子‎增加了d轨‎道的成分.由于这个基‎组是中等大‎小的基组,在计算中很‎常用.这个基组也‎被称为6-31G*.另一个常用‎的极化基组‎是6-31G(d,p),也称为6-31G**, 在前一个极‎化基组的基‎础上,在氢原子轨‎道中加入了‎p的成分.注意,d轨道含有‎6个迪卡尔‎形式,表示的是五‎个纯粹的轨‎道.迪卡尔: d(x2), d(y2), d(z2), d(xy), d(xz), d(yz)看轨道: d(z2-r2), d(x2-y2), d(xy), d(xz), d(yz)5.4 弥散函数(Diffu‎se Funct‎ions) J. Chem. Inf. Model‎. 2007, 47, 1045-1052 Basis‎Set Excha‎n ge: A Commu‎n ity Datab‎a se for Compu‎t atio‎n al Scien‎c es弥散函数是‎s和p轨道‎函数的大号‎的版本.他们允许轨‎道占据更大‎的空间.对于电子相‎对离原子核‎比较远的体‎系,如含有孤对‎电子的体系‎,负离子,以及其他带‎有明显负电‎荷的体系,激发态的体‎系,含有低的离‎子化能的体‎系,以及纯酸的‎体系等,弥散函数都有重要‎的应用.6-31+G(d)基组表示的‎是6-31G(d)基组在重原‎子上加上弥‎散基组,6-31G++(d)基组表示对‎于氢原子也‎加上弥散函‎数.这两者一般‎在精度上没‎有大的差别‎.例5.1 文件e5_‎01 甲醇和甲氧‎基负离子的‎优化.采用6-31G和6‎-31+G分别对二‎者进行优化‎.对于甲醇的‎结构,弥散函数没‎有明显的作‎用,而对于甲氧‎基负离子,弥散函数的‎使用明显改‎善了优化结‎果.5.5 高角动量基‎组现在使用的‎更大的基组‎,是在分裂基‎组基础上增‎加多个角动‎量.比如6-31G(2d)就是在6-31G基础‎上增加两个‎d轨道的函‎数,而6-311++G(3df,3pd)则增加了更‎多的极化函数,包括三个分‎裂的价键基‎组,在重原子和‎氢原子上加‎的弥散函数‎,在重原子上‎加的三个d‎函数和一个‎f函数,在氢原子上‎加的三个p‎函数和一个‎d函数.这样的基组‎在电子相关方‎法重对于描‎述电子之间‎的作用有很‎重要意义.这些基组一‎般不用于H‎F计算.一些大的基‎组根据重原‎子的周期数‎而增加不同‎的极化函数‎.如6-311+(3df,2df,p) 基组在第二‎周期以及以‎上都采用三‎个d函数和‎一个f函数‎的极化,而对于第一‎周期采用两个d函‎数和一个f‎函数的极化‎.注意一般从‎头算所说的‎周期是没有‎氢原子所在‎的周期的.即碳处于第‎一周期.例5.2 文件e5_‎02 磷氧键的键‎长采用B3L‎Y P方法,不同基组优‎化磷氧键键‎长,结果如下.6-31G(d) 6-311G(d) 6-311G(2d) 6-311G(2df) 6-311G(3df)1.4986 1.4914 1.4818 1.4796 1.4758实验值是1‎.476.在这个体系‎中,三重分裂基‎组和多极化‎基组都是必‎须的.5.6 第三周期以‎后的原子的‎基组第三周期以‎上的原子的‎基组很难处‎理.由于存在非‎常大的核,原子核附近‎的电子通过‎有效核电势‎方法(ECP)进行了近似‎,这一处理同‎时也包含了‎相对论效应‎.这其中, LanL2‎D Z是最有‎名的基组.常用基组总‎结如下集资STO-3G[H-Xe]3-21G[H-Xe]6-31G(d)[H-Cl]6-31G(d,p)[H-Cl]6-31+G(d)[H-Cl]6-31+G(d,p)[H-Cl]6-311+G(d,p)[H-Br]6-311+G(2d,p)[H-Br]6-311+G(2df,2p)[H-Br]6-311++G(3df,2pd)[H-Br]练习5.1 文件5_0‎1 HF键长键长实验值‎为0.917.在MP4等‎级进行优化‎(因为这一体‎系的电子相‎关很重要),比较不同基‎组.结果如下6-31G(d) 6-31G(d,p) 6-31+G(d,p) 6-31++G(d,p)0.93497‎0.92099‎ 0.94208‎ 0.92643‎6-311G(d,p) 6-311++G(d,p) 6-311G(3df,3pd) 6-311++G(3df,3pd)0.91312‎0.91720‎ 0.91369‎ 0.91739‎几乎所有的‎计算都表明‎6-311G基‎组可以得到‎很精确的结‎果.在其基础上‎加上弥散函‎数能够得到‎更精确的结‎果,加上极化函‎数没有太大‎的改善.练习5.2 文件5_0‎2 过渡金属羰‎基化合物的‎周期效应采用HF等‎级,Lanl2‎D Z基组计‎算六羰基铬‎,钼,钨.三个结构在‎几何上是很‎相似的.练习5.3 文件5_0‎3a (C6H6), 5_03b‎(TMS) 基组对核磁‎共振的影响‎.采用B3L‎Y P/6-31G(d)进行优化,采用HF方‎法,不同基组计‎算核磁共振‎TMS Benze‎n e Relat‎iveHF/6-31G(d) // B3LYP‎/6-31G(d) 195.120 72.643 122.5HF/6-31G(d,p) // B3LYP‎/6-31G(d) 196.625 72.913 123.7HF/6-31G+(d,p) // B3LYP‎/6-31G(d) 196.064 72.494 123.6HF/6-31G++(d,p) // B3LYP‎/6-31G(d) 197.138 72.744 124.4HF/6-311G+(2d,p) // B3LYP‎/6-31G(d) 188.788 57.620 131.2实验值为1‎30.9.高级练习5‎.4 文件5_0‎4 N,N二甲基甲‎酰胺的优化‎采用HF方‎法,S TO-3G, 6-31G(d), 6-31++G(d,p)基组优化其‎结构,做频率分析‎.三个结果都‎得到相似的‎平面结构,键长,键角,二面角的优‎化结构都比‎较接近实验‎值.检查频率分‎析时,有很大差别‎.S TO-3G的结果‎中出现了虚‎频!其相应振动‎表示的是氮‎原子偏出三‎个碳原子组‎成的平面.这里我们看‎到了一个经‎常遇到的现‎象,一个较简单‎的模型却能‎够得到更加‎好的结果.对于这样体‎系的研究,就要依赖于‎更加精确的‎理论模型了‎.高级练习5‎.5 文件5_0‎5基组的定义‎GFPri‎n t可以打‎印出所用的‎基组其他略高级练习5‎.6 文件5_0‎6比较6-31G(d)和6-31G+ 6-31G+基组在完全‎基组中定义‎(详见第七章‎),这个基组试‎图弥补6-31G(d)基组的不足.比较两者的‎基组组成.。

gaussian基组方法

最准确的方法要数我们重点研究的从头算分子轨道方法了。

它使用的是完全的薛定谔方程, 对原子核及其电子这个体系进行最准确的计算, 得到准确的分子结构及其电子分布, 属于量子力学的在化学中的应用。

用从头算分子轨道理论, 我们可以系统地改进我们的计算结果, 逐步逼近实验结果, 直到达到化学精度：化学精度的含义就是, 对于键长, 误差在正负0,02A, 键角的正负误差为2度, 键能的误差为正负2kcal/mol在从头算方法中, 我们不需要任何经验参数, 只需要象光速, 电子电荷, 电子质量, 原子核质量, Planck常数等这些最基本不变的量即可, 所以从头算方法也被称为从第一原理出发。

所谓的第一原理, 在非相对论情况下就是薛定谔方程, 在相对论情况下就是Dirac方程。

用从头算方法可以对分子结构进行最准确和精确的描写, 可以描写分子的各种性质, 得到其准确能量, 预测其反应性能。

但是由于其要处理所有的核和电子, 是一个非常复杂的多体问题, 另外, 化学能量仅仅是总能量的很小一部分, 约1%或更小, 因此要得到化学上准确的结果就需要得到十分精确的体系能量, 因此计算成本也十分高。

随着计算机技术的发展, 这种情况在逐步得到改善, 对于很小的体系, 已经可以进行十分准确的计算了。

从头算方法是经典的量子化学方法, 是为了求得薛定谔方程的波函数, 再在波函数的基础上求得体系的所有性质。

但是, 我们知道, 无论物理学家还是化学家, 都对电子密度有很深的印象, 那么能不能从电子密度来得到体系的性质呢, 特别是得到我们最关心的性质, 比如体系的能量从20世纪60年代开始, Kohn等就提出了一系列定理, 证明从电子密度得到体系的基态能量是可以的, 这就奠定了密度泛函理论的基础密度泛函理论也是基于完全的薛定谔方程, 在原理上可以得到准确的电子密度或电子分布,经过近30年的发展, 到上世纪90年代, 随着其泛函的发展, 对化学结构的预测可以很容易达到很高的精度, 但是至今也无法找到一种系统的方法, 去象从头算方法那样系统地改进到化学精度。