目录

内反应坐标IRC

前提

用户已经完成过渡态搜索,得到过渡态结构。IRC会分别得到从过渡态结构出发,沿着反应物、产物方向的变化路径,共同的起点是过渡态结构。本例以$CH_4-F_2$取代反应的过渡态为例,计算IRC。

参数设置

泛函、基组、精度设置应该与过渡态搜索一致,便于比较在相同参数下的能量变化。过渡态搜索、IRC计算,参数标准可以采用结构优化相同精度的参数。

IRC计算,需要设置的参数非常少,分子结构为过渡态结构,Main面板选择IRC以及正确的基组、泛函等参数,IRC计算并不需要很高精度,只是为了校验过渡态对应的反应物、产物物种,所以选择了较小的基组DZP,并且设置了Frozen core为large。具体参考:ADF参数设置详解

IRC的详细参数,需要设置的也不多,Direction默认选择了Both,也就是包含从过渡态-反应物、过渡态-产物。

但注意: Forward path、Backward path,与反应物、产物并没有一一对应的关系,因为对程序来说,反应物与产物只是相对而言的。 如果计算过频率,这会对IRC计算有帮助,因此可以在Initial Hessian这里读取之前计算频率生成的*.results/adf.rkf文件:

计算过程中,可以使用SCM > Movie 监控结构的变化。IRC实际上相当于一系列的限制性结构优化。这些优化里面,没有收敛的结构,一般我们并不关心,所以可以只显示收敛的点:Movie窗口 > View > Accepted Geometries Only,如此显示较为光滑的IRC曲线。计算完毕后,程序会重组一下数据,让反应物、过渡态、产物,连接成连续的变化的过程:

当然这时候,Movie需要读取生成的*.results/ams.rkf文件。

AMS软件提供免费试用,试用申请方式参见AMS免费试用