计算化学公社

标题: 关于IRC问题 [打印本页]

作者
Author: Hao 时间: 2014-11-6 02:14
标题: 关于IRC问题
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结果文件中取出来再放入输入文件中进行优化？

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

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

作者
Author: comchem 时间: 2014-11-6 16:31
这种情况下一般是过渡态有问题，可先检查一下过渡态的虚频对应的振动情况是否合适。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)