关于IRC问题

Hao

1. 我在计算过渡态的IRC时，输入文件设置了maxpoint=20（irc=(rcfc,forward,maxpoint=20,stepsize=10)） ，但是计算结果仅跑了二个点，而且正常结束，并且提示PES minimum detected on this side of the pathway.。可是第二个点的能量和过渡态的相近，但是结构不是我预计的稳定体。请教：怎么判断这个结果正确与否呢？

SUMMARY OF CORRECTOR INTEGRATION:

   Predictor End Point Energy = -4459.517180

   Old End Point Energy       = -4459.517180

   Corrected End Point Energy = -4459.517180

   Predictor End-Start Dist.  =     0.317673

   Old End-Start Dist.        =     0.316790

   New End-Start Dist.        =     0.316823

   New End-Old End Dist.      =     0.002495

CORRECTOR INTEGRATION CONVERGENCE:

   Recorrection delta-x convergence threshold:    0.010000

   Delta-x Convergence Met

Point Number:   1          Path Number:   1

   CHANGE IN THE REACTION COORDINATE =    0.31682

   NET REACTION COORDINATE UP TO THIS POINT =    0.31682

  # OF POINTS ALONG THE PATH =   1

  # OF STEPS =   2

PES minimum detected on this side of the pathway.

    Magnitude of the gradient =      0.0000673

Calculation of FORWARD path complete.

Reaction path calculation complete.

Energies reported relative to the TS energy of       -4459.517165

--------------------------------------------------------------------------

    Summary of reaction path following

--------------------------------------------------------------------------

                        Energy    RxCoord

   1                    0.00000   0.00000

   2                   -0.00002   0.31682

--------------------------------------------------------------------------

2. 如果我想进行几何优化这个第二个点，请问：我可不可以直接读取这个IRCf.chk（不解的地方是： 直接geom=check 是表示读取的是最后那点的几何坐标吗？）还是需要我把最后一个点的几何坐标从IRCf结果文件中取出来再放入输入文件中进行优化？

sobereva

1 这样的结果不对。通常是因为过渡态结构定位得不够精确引起的。用int=ultrafine结合tight几何优化收敛限重新定位过渡态，然后再试。还不行的话calcall也用上。

2 此情况不建议用geom=check，免得搞乱。

comchem

这种情况下一般是过渡态有问题，可先检查一下过渡态的虚频对应的振动情况是否合适。