|
Gaussian 讲解64-Gaussian 09 Freq(2)简正模式有关的选项 HPModes 在频率输出中除了正常的三个数字输出外,还包括高精度格式(五个数字) 的振动频率本征矢。 InternalModes 把模式打印成冗余内坐标的位移。IntModes 是这个选项的同义词。 SaveNormalModes 把所有模式保存到检查点文件。SaveNM 是该选项的同义词。默认为 NoSaveNormalModes 或 NoSaveNM。 ReadNormalModes 从检查点文件读取保存的模式。ReadNM 是该选项的同义词。默认为 NoReadNormalModes 或 NoReadNM。 SelectNormalModes 读取输入,选择特定的模式用于显示。SelectNM 是该选项的同义词。默认为 NoSelectNormalModes或NoSelectNM。AllModes表示在输出中包含所有模式。 输入部分格式的说明见上。注意,这个选项不影响 SaveNormalModes 的运行, 它总是把所有模式保存到检查点文件。 SortModes 在输出中,按照 ONIOM 层对模式进行排序。 ModelModes 在 ONIOM 计算中,只显示涉及最小模型体系的模式。 MiddleModes 在三层 ONIOM 计算中,只显示涉及两个模型体系的模式。 PrintFrozenAtoms 在模式的输出中,默认不打印冻结原子的零位移。该选项要求列出所有的原 子。 修改分子说明的选项 ModRedundant 读入对冗余内坐标的修改(即,和InternalModes合用)。注意,在Opt Freq 计算中,对优化和模式分析使用相同的坐标,这与Opt=ModRedundant相同。 关于输入数据格式的细节,参见Opt关键词的说明。 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)。 关于算法和流程的选项 Analytic 指定解析计算能量二阶导。该选项仅能用于 RHF,UHF,CIS,CASSCF,MP2 和全部 DFT 方法,它对这些方法是默认的。 Numerical 请求对能量的一阶解析导进行数值导数计算,得到能量的二阶数值导数。它 可以用于任何具有梯度的方法,并且对那些只能计算解析梯度而不能计算解 析二阶导数的方法是默认的。Freq=Numer 可以和 Polar=Numer 合用于同一个 任务步骤中。 DoubleNumer 请求对能量进行二次数值插分,产生力常数。对于没有解析梯度的方法,这 是默认且唯一的选择。EnOnly 是 DoubleNumer 的同义词。 Cubic 进行解析二阶导数的数值差分,产生三阶导数。仅用于有解析频率但没有解 析三阶导的方法。 Step=N 定义数值差分的步长为 0.0001*N (长度单位:埃,除非定义了 Units=Bohr)。 如果同时使用了 Freq=Numer 和 Polar=Numer,N 也用来定义电场的步长。对 于 Hartree-Fock 和关联方法的 Freq=Numer,默认是 0.001 Å,对于 GVB 和 CASSCF 的 Freq=Numer,默认是 0.005 Å,对 Freq=EnOnly 默认是 0.01 Å。 对 Freq=Anharmonic 或 Freq=VibRot 默认是 0.025 Å。 Restart 该选项用于在最后一个完成的结构计算后,重新开始数值频率计算。失败的 数值频率任务可以从检查点文件重新开始,只需简单地重复原始任务的执行 路径,并对 Freq=Numer 关键词添加 Restart 选项。不需要其它输入。 如果命名并保存了读写文件,也可以用 Restart 关键词从失败的任务重新开 始解析频率计算。 DiagFull 对角化整个(3Natoms)2 力常数矩阵——包括平移和转动自由度——并报告最低 的振动频率,进行频率计算数值稳定性的测试。这个计算在简正频率分析之 前进行,因为后者将投影掉这些模式。它的输出报告最低的九个模式,其中 较高的三个对应于常规频率分析的三个最低模式。理想情况下,这个分析得 到的最低六个模式具有非常小的数值。如果它们显著非零,则可能表示平移 和转动是该体系的重要模式,需要更好的积分网格,结构未收敛,等。应对 该体系作进一步研究,以获得更准确的频率。参见下面的例子部分关于该选 项的输出。 默认为 DiagFull;NoDiagFull 表示跳过这个分析。 ReadFC 从检查点文件读取前一频率计算中获得的力常数,并使用不同温度、压强或 同位素,用最少的计算花费重新进行模式分析和热化学分析。注意,因为基 组从检查点文件读取,因此不用输入基组。如果前一任务中还指定了 Raman 选项,如果还想使用这个选项也无需再指定它。 TwoPoint 在计算数值导数时,在每个方向做两次位移。这是默认的。FourPoint 将做 四次位移,但只能用于链接 106(Freq=Numer)。对 Freq=DoubleNumer 无效。 NFreq=N 请求用 Davidson 对角化求解最低的 N 个频率。目前,该选项仅能用于 ONIOM(QM:MM)模型化学。 |