第四章_Gaussian软件的使用

1、 分子的能量和结构分子的能量和结构 过渡态的能量和结构过渡态的能量和结构 振动频率振动频率 红外和拉曼光谱（包括红外和拉曼光谱（包括预共振拉曼）预共振拉曼） 热化学性质热化学性质 成键和化学反应能量成键和化学反应能量 化学反应路径化学反应路径 分子轨道分子轨道 原子电荷原子电荷 电多极矩电多极矩 NMR NMR 屏蔽和磁化系数屏蔽和磁化系数 自旋自旋- -自旋耦合常数自旋耦合常数 振动圆二色性强度振动圆二色性强度 电子圆二色性强度电子圆二色性强度 g g 张量和超精细光谱的张量和超精细光谱的其它张量其它张量 旋光性旋光性 振动振动- -转动耦合转动耦合 非谐性的振动分析和振非谐性的振动分析和振

2、动动- -转动耦合转动耦合 电子亲和能和电离势电子亲和能和电离势 极化和超极化率（静态极化和超极化率（静态的和含频的）的和含频的） 各向异性超精细耦合常各向异性超精细耦合常数数 静电势和电子密度静电势和电子密度 Gjf 计算任务文计算任务文件件 Chk 临时文件临时文件 Rwf 读写文件读写文件 Int 双电子积分双电子积分 D2e 双电子积分双电子积分导数文件导数文件 PDB， ENT MOL， HIN CDX， C3D Schrodinger方程方程 Born-Oppenheimer近似近似 Hellmann-Feyman定理定理 维里定理维里定理基组的影响 基组是体系内轨道的数学描述。大

3、的基组由于对电子在空间上有小的限制而具有更大的精确性。 用于电子结构计算的标准的基组使用线性的高斯函数来模拟轨道。最小基组 最小基组包含了描述轨道的最少的函数数量。 H：1s C：1s，2s，2px，2py，2pz STO-3G 是最小基组分裂基组 增大基组的第一个方法就是增加每个原子基函数的数量。分裂基组，比如3-21G 和6-31G，对于价键轨道都用两个函数来进行描述，比如 H：1s，1s C：1s，2s，2s，2px，2py，2pz，2px，2py，2pz双zeta 基组，如Dunning-Huzinaga 基组(D95)，采用每个原子的两种不同大小的函数的线性组合来描述分子轨道。同样的

4、，三重分裂基组，如6-311G，采用三个不同大小的收缩函数来描述轨道。极化基组 分裂基组允许轨道改变其大小，但不能改变形状。极化基组则取消了这样的限制，增加了角动量。 比如在碳原子上增加d 轨道的成分，在过渡金属上增加f 轨道成分。有些在氢原子上增加p 轨道成分 6-31G(d)、6-31G(d，p)弥散函数(Diffuse Functions) 弥散函数是s 和p 轨道函数的大号的版本。他们允许轨道占据更大的空间。对于电子相对离原子核比较远的体系，如含有孤对电子的体系，负离子，以及其他带有明显负电荷的体系，激发态的体系，含有低的离子化能的体系，以及纯酸的体系等，弥散函数都有重要的应用。 6-

5、31+G(d)基组表示的是6-31G(d)基组在重原子上加上弥散基组，6-31G+(d)基组表示对于氢原子也加上弥散函数。高角动量基组 6-31G(2d)就是在6-31G 基础上增加两个d 轨道的函数，而6-311+G(3df，3pd)则增加了更多的极化函数，包括三个分裂的价键基组，在重原子和氢原子上加的弥散函数，在重原子上加的三个d 函数和一个f 函数，在氢原子上加的三个p 函数和一个d 函数。这样的基组在电子相关方法对于描述电子之间的作用有很重要意义。第4周期以后的原子的基组 第4周期以上的原子的基组很难处理。由于存在非常大的核，原子核附近的电子通过有效核电势方法(ECP)进行了近似，这一

6、处理同时也包含了相对论效应相对论效应。这其中，LANL2DZ 是最有名的基组。GaussianOverlay1Overlay9,10,11,99L101L102L122Overlay0L0L001 Gaussian 程序的结构程序的结构link0: 初始化程序，控制初始化程序，控制overlaylink1: 读入并处理读入并处理Route Section，建立要执行的建立要执行的link列表列表link9999: 终止计算终止计算Overlay99L9999-M- 2MW-#- MaxDisk=400MBGaussian程序使用的内存单位W是双精度字，相当于8字节 2MW=16MB设置方法：将

7、Default.r1文件改成default.rouPreferences:对对Gaussian程序进行初程序进行初始化设置始化设置自定义外部文字编辑器，用来打开自定义外部文字编辑器，用来打开.out文文件件link.exe所在的文件夹所在的文件夹临时文件存放文件夹临时文件存放文件夹开始作业开始作业 暂停当前作业暂停当前作业 当前当前link后暂停后暂停 终止当前作业和批处理终止当前作业和批处理 恢复当前作业恢复当前作业 在当前作业完成后终止批处理在当前作业完成后终止批处理 终止当前作业终止当前作业 编辑或建立批处理编辑或建立批处理 批处理专用批处理专用打开外部编辑器打开外部编辑器 编辑输出文件

8、编辑输出文件 将.chk文件转换为.fch文件，这种 文件可以使用图形软件打开 将.fch文件还原为.chk文件 从指定.chk文件中显示作业的 route section和title 将.chk文件转换成文本格式 将以前版本的Gaussian产生的 .chk文件转换为G03的.chk文件 编辑批处理作业文件 转换不同格式的分子结构文件 读取.fch文件中的数据并生成 三维空间网格图 利用.chk文件中的分子轨道， 生成电子密度和静电势的空间 分布网格图 从.chk文件中打印出频率 和热化学数据NewZmat工具界面工具界面H2O分子单点计算后的分子单点计算后的.chk文件文件转换成转换成.f

9、ch文件后，在文件后，在Chem3D中所作分子轨道中所作分子轨道(HOMO)，电荷电荷密度和静电势空间分布图密度和静电势空间分布图% Section 设定作业运行的环境变量 Route Section 设定作业的控制项Title 作业题目Charge & Multipl 电荷与自旋多重度Molecule Specification 分子说明%chk=water.chk % Section： 行首以行首以%开始，段开始，段后无空行后无空行%rwf=water.rwf#p hf/6-31g scfcyc=250 scfcon=8 Route Section：行首以行首以# #开始，段后加空行开始，

10、段后加空行 Water ennergy title：作业的简要描述，段后加空行作业的简要描述，段后加空行 0 1 Molecular Specification：O 分子说明部分，段后通常加空行分子说明部分，段后通常加空行H 1 R1H 1 R1 2 a1R1=1.04a1=104.0 Gaussian作业的格式作业的格式 定义计算过程中的临时文件%chk=name.chk.chk文件在计算中记录分子几何构型，分子轨道，力常数矩阵等信息%rwf=name.rwf.rwf文件主要在作业重起时使用，当计算量比较大时，.rwf文件通常会非常大，此时需要将之分割保存%int=, %d2

11、e=name.d2e .int文件在计算过程中存储双电子积分, .d2e文件在计算过程中存储双电子积分的二阶导数 内存使用控制%mem=n 控制运行过程中使用内存的大小,可以以W或者MB，GB为单位default：6000000W48MB 综合考虑到计算的需要和硬件水平，内存并非给得越多越好，最有效率的方法是根据作业类型估算所需要内存的大小M + 2NB2 M：不同类型作业需要的最小内存， NB ：计算所使用基函数的数目 作业类型作业类型 f g h i j SCF能量能量 4 MW 4 MW 9 MW 23 MW 60 MW SCF梯度梯度 4 MW 5 MW 16 MW 38 MWSCF振

12、动分析振动分析 4 MW 9 MW 27 MW MP2能量能量 4 MW 5 MW 10 MW 28 MW 70 MW MP2梯度梯度 4 MW 6 MW 16 MW 38 MWMP2振动分析振动分析 6 MW 10 MW 28 MW Route Section以# 开始，# 控制作业的输出#N 正常输出；默认 (没有计算时间的信息)#P 输出更多信息。包括每一执行模块在开始和结束时与计算机系统有关的各种 信息 (包括执行时间数据，以及SCF计算的收敛信息)#T 精简输出：只打印重要的信息和结果。 Route Section主要由方法，基组，任务类型三部分组成方法与基组后续课程专门介绍Gaus

13、sian程序能完成的任务类型：Sp单点能量计算（默认任务类型）ADMP&BOMP动力学计算动力学计算Opt分子几何结构优化Force计算核的受力计算核的受力Freq振动分析Stable波函数稳定性测试波函数稳定性测试Irc反应途径计算Volume计算分子体积计算分子体积IrcMax在指定反应途径上找能量最大值在指定反应途径上找能量最大值Density=CheckPoint仅计算布居分析仅计算布居分析Scan势能面扫描势能面扫描Guess=Only仅作分子轨道初猜仅作分子轨道初猜Polar极化率和超极化率计算极化率和超极化率计算ReArchive从从.chk文件中提取存档文件中提取存档Gauss

14、ian 程序中的部分关键词程序中的部分关键词 Guess Guess=read Guess=mix Guess=alter SCFCYC=n default=64 SCFCON=n SCFDM SCF=DM SCFQC SCF=QC direct (default) 优化平衡态 opt=z-matrix 优化过渡态 opt=(ts,z-matrix,noeigentest) optcyc=n Iop(1/8=6) Freq=noraman Freq=ReadIsotopes Route Section采用自由格式，大小写不敏感 同一行不同项之间可以使用空格，逗号和“/”连接； 例：#p hf/

15、6-31g scfcyc=230 scfcon=8 #p,hf/6-31g,scfcyc=230,scfcon=8 关键词可以通过 keyword =option,keyword(option), keyword=(option1, option2, .)，keyword( option1, option2, .)指定 例：#p HF/STO-3G opt #p HF/STO-3G opt=(TS,readfc) Title部分必须输入，但是程序并不执行，起标识和说明作用 输入分子的电荷和自旋多重度 例： 电荷 多重度2s+1 H2O 0 1 H3O+ 1 1 NO 0 2 O2 0 2 电荷

16、多重度部分通常也算作分子说明 分子说明部分主要用来定义分子核相对位置分子说明部分主要用来定义分子核相对位置分子核相对位置可以用笛卡尔坐标，内坐标分子核相对位置可以用笛卡尔坐标，内坐标(Z-matrix)，或者是二者混合，或者是二者混合 表示表示笛卡尔坐标是内坐标的一种特殊形式笛卡尔坐标是内坐标的一种特殊形式分子坐标的格式为分子坐标的格式为: 元素符号, x, y, z 元素符号 (n) 原子1 键长 原子2 键角 原子3 二面角 格式代码 (0,180) 以原子在分子中的序数表示 通过右手规则确定例例1：使用使用HF方法，优化方法，优化 H2O2分子分子%chk=h2o2.chk%rwf=h2

17、o2.rwf#p hf/6-31g optH2O2 energy calculation0 1 HO 1 0.9O 2 1.4 1 105.0H 3 0.9 2 105.0 1 120.0%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1H 0.000 0.000 0.000 O 0.000 0.900 0.000 O 1.350 1.262 0.000 H 1.464 1.742 -0.752 内坐标表示内坐标表示 笛卡尔直角坐标表示笛卡尔直角坐标表示上述两个作业将对上述两个作业将对H2O2分子的结构进行完全

18、优化，包括所有的键长键角和二面角。分子的结构进行完全优化，包括所有的键长键角和二面角。通过将内坐标定义成变量，可以对分子结构进行部分优化通过将内坐标定义成变量，可以对分子结构进行部分优化%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1 HO 1 r1O 2 r2 1 a1H 3 r1 2 a1 1 d1r1=0.9 r2=1.4a1=105.0d1=120.0 通过对两个键长和键角使用同通过对两个键长和键角使用同一变量定义可以控制分子的对一变量定义可以控制分子的对称性称性%chk=h2o2.chk%rwf=h2o2.rwf#p hf/6-31g optH2O2 energy calculation0 1 HO 1 r1O 2 r2 1 a1H 3 r1 2 a1 1 d1r1=0.9a1=105.0 变量部分变量部分 d1=120.0r2=1.4 常量部分常量部分通过设定常量，可以使程序通过设定常量，可以使程序只对分子结构进行只对分子结构进行部分部分优化优化例例2： C2H4分子分子 平面形分子平