|
Gaussian 讲解71-Gaussian 09 IRCMax说明 使用 Petersson 及合作者[Eyring35, Truhlar70, Truhlar71, Garrett80, Skodje82, Malick98, Petersson98, Schwartz98, Petersson98a]的方法进行 IRCMax 计算。这一计算类型用过渡结构做为输入，用 GS2 算法沿着 指定的反应路径寻找最大的能量[Gonzalez89, Gonzalez90]（细节参见 IRC=GS2）。 必须用标准方法为 IRCMax 任务指定可选的同位素。 必需的输入 IRCMax 需要两个模型化学作为选项，中间用冒号隔开：IRCMax(模型 2:模型 1)。这 里是一个执行路径的例子： # IRCMax(B3LYP/6-31G(d,p):HF/6-31G(d,p)) 这个计算在 HF/6-31G(d,p)反应路径上寻找点，该点使 B3LYP/6-31G(d,p)能量位于最大 值。 选项 Zero 产生用于零曲率变分过渡态理论（ZC-VTST）[Eyring35, Truhlar70, Truhlar71, Garrett80, Skodje82, Petersson98]的数据。考虑下面的执行路径： # IRCMax(MP2/6-311G(d):HF/6-31G(d),Zero,Stepsize=15,CalcAll) 这个任务将从 HF/6-31G(d)过渡态开始，沿着 HF/6-31G(d)的 IRC 以步长 0.15 amu1/2 Bohr 进 行寻找，直到找到 MP2/6-311G(d)能量的最大值（含 HF/6-31G(d)零点能）为止。对 MP2/6-311G(d)过渡态沿着 HF/6-31G(d)的 IRC 进行优化。用 ZC-VTST 版本绝对速率理论进 行反应速率计算所需的全部量都包含输出中：过渡态的转动惯量，所有的实振动频率 （HF/6-31G(d)），隧穿的虚频（拟合到 MP2/6-311G(d)，加上 ZPE），以及过渡态的 MP2/6-311G(d)加上 ZPE 的总能量。注意，如果未使用 CalcAll，那么所有这些量（ZPE，频 率，等）只在 HF/6-31G(d)级别进行计算，并且这些量用于 IRCMax 路径上所有的点。 ZC-VTST 选项 Zero 在 IRCMax 计算中包含零点能。 选则路径选项 Forward 只沿着反应路径的正向进行跟踪。 Reverse 只沿着反应路径的反向进行跟踪。 ReadVector 读入矢量进行跟踪。格式为 Z-矩阵（FFF(I), I=1,NVAR），以（8F10.6）读入。 MaxPoints=N 沿反应路径进行检查的点数（沿着正反两个方向，如果计算中都考虑的话）。 默认是 6。 StepSize=N 反应路径上的步长，单位是 0.01 amu1/2-Bohr。默认是 10。 MaxCyc=N 设定每个几何优化中的最大步数。默认是 20。 选择坐标系的选项 MassWeighted 跟踪质量加权内坐标（Z-矩阵）的路径（和跟踪质量加权笛卡尔坐标的路径 是等价的）。MW 是 MassWeighted 的同义词。这是默认选项。 Internal 跟踪内坐标（Z-矩阵）路径，不使用质量加权。 Cartesian 跟踪笛卡尔坐标路径，不使用质量加权。 收敛有关的选项 VeryTight 在反应路径每一点的优化中，使用严格的收敛标准。该选项对反应路径上很 小的步长是必要的。 产生初始力常数选项 CalcFC 指定在第一个点计算力常数。 CalcAll 指定在每一个点计算力常数。在路径上的每一优化点，对运动计算投影的振 动频率[Baboul97]。注意，这些投影的振动频率仅对质量加权内坐标的反应路径 计算有效。 FCCards 指定从分子说明后面的输入串读入笛卡尔坐标的力和力常数。该选项使用内 部命令 FCList，可以从 Quantum Chemistry Archive 读取力常数。输入的格式 为： 能量 （格式 D24.16） 笛卡尔力 （格式为 6F12.8 的几行） 力常数 （格式为 6F12.8 的几行） 力常数为倒三角形式：((F(J,I),J=1,I),I=1,3Natoms)，其中 3Natoms 是笛卡尔坐标 数量。如果同时指定 FCCards 和 ReadIsotopes，则在能量，笛卡尔梯度和/ 或笛卡尔力常数之前输入原子量。 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 重新开始没有完成的 IRCMax 计算，或对已完成的 IRCMax 计算在反应路径 上加入新的点继续计算。 适用性 对计算的 IRC 部分（上面的模型 1）需要解析梯度。对模型 2 可以使用任何非复合的能 量方法和基组。 有关的关键词 IRC，Opt，Freq |
