Gaussian 讲解67-Gaussian 09 G1-G4 方法
G1,G2,G2MP2,G3,G3MP2,G3B3,G3MP2B3,G4,G4MP2
说明
这些方法关键词请求用 Gaussian-1(也就是常说的 G1)[Pople89, Curtiss90],Gaussian-2(G2) [Curtiss91],Gaussian-3(G3)[Curtiss98]和 Gaussian-4(G4)[Curtiss07]方法计算非常精确的能 量。G2MP2 进行修改版的 G2,也就是 G2(MP2)的计算,它对基组展开校正使用 MP2 而不 是 MP4 [Curtiss93],不仅具有接近完全 G2 方法的精度,而且节省了计算时间。G3MP3 进行做 了类似修改的 G3(MP2) [Curtiss99]。可以用 G3B3 和 G3MP2B3 关键词调用 G3 的变体,其中使用 B3LYP 的结构和频率[Baboul99]。G4 和 G4MP2 使用第四代 Gaussian 方法[Curtiss07, Curtiss07a]。
所有这些方法都是复合能量计算,其中对指定的分子体系作多个预定义的计算。在指定 了其中一个关键词后,自动执行所有这些计算步骤,并在输出中显示最终计算的能量值。这 些关键词不需要定义基组关键词。 在选择这些方法之前,用户应当考虑其它的高精度方法。CBS-QB3 同样精确但明显更 快,而 W1U 更准确(但更慢)。 Opt=Maxcyc=n,QCISD=Maxcyc=n 或 CCSD=Maxcyc=n 关键词可以与这些关键词结 合使用,分别用于定义优化的最大步数和 QCISD 或 CCSD 循环的最大步数。
选项
ReadIsotopes
该选项允许指定不同于默认值(分别为 298.15 K,1 大气压,无换算,和最 丰同位素)的温度,压强,频率换算因子,和/或同位素。该选项可用于重新 进行分析,但使用不同于检查点文件数据的参数。
但是需要注意,所有这些参数也可以在执行路径部分(关键词 Temperature, Pressure,和 Scale)和分子说明部分(Iso=参数)指定,如这个例子:
#T Method/6-31G(d) JobType Temperature=300.0 ...
...
0 1
C(Iso=13)
...
ReadIsotopes 输入具有下面的格式:
temp pressure [scale] 数值必须为实数。
原子 1 的同位素质量
原子 2 的同位素质量
...
原子 n 的同位素质量
其中的 temp,pressure,和 scale 分别为所需的温度,压强,和当用于热化 学分析时可选的频率数据换算因子(默认不换算)。其余的行是分子中各原子的同位素质量,按照分子说明部分出现的相同顺序排列。如果原子质量指 定为整数,程序将自动使用相应的真实同位素质量(例如,对 18O 指定为 18, 那么 Gaussian 使用的值为 17.99916)。
Restart 从检查点文件恢复部分完成的计算。当与 ReadIso 选项同时使用时,这个选 项允许用不同的热化学参数和/或同位素选择,快速计算能量。
例子
计算的输出总结。在任务的所有组成部分输出之后,Gaussian 为这些方法打印结果的表格。以下是 G2 计算的输出:
开始是温度和压强,之后是用于计算 G2 能量的各种成分。输出的最后是在 0 K 和在指 定温度时的 G2 能量(后者包括完整的热校正,而不仅仅是零点能校正),以及(在输出的 最后一行)G2 理论预测的焓和吉布斯自由能(都用热校正的 G2 能量计算)。(注意,这个 总结部分还打印了在 G1 级别预测的同样的量)。
能量标记有如下含义(作为例子,使用 G2):
G2 (0 K) 零点校正电子能量:E0 = Eelec + ZPE
G2 Energy 热校正能量:E = E0 + Etrans + Erot + Evib
G2 Enthalpy 用 G2 预测的能量计算的焓:H = E + RT
G2 Free Energy 用 G2 预测的能量计算的 Gibbs 自由能:G = H - TS
在不同温度下重新进行计算。下面的两步任务演示了在不同温度下(快速)运行第二个 G2 计算的方法。这个任务先在 298.15 K 计算 G2 能量,接下来在 300 K:
推荐
-
-
QQ空间
-
新浪微博
-
人人网
-
豆瓣