计算化学公社

标题: 做IRC时过渡态不在能量最高点 [打印本页]

作者
Author: wch990122 时间: 2021-10-15 16:15
标题: 做IRC时过渡态不在能量最高点
各位老师好，我在寻找图一两个构型间的过渡态时遇到了问题。我的步骤如下：
1.首先对键角进行柔性扫描。
柔性扫描命令如下：
#p opt=modredundant ub3lyp/6-31g(d) nosymm geom=check guess=read
扫描结果如图2所示。

2.选取柔性扫描能量最高点为过渡态初猜结构，寻找过渡态。
TS命令如下：
#p opt=(TS,calcfc,noeigen,tight) freq ub3lyp/6-31g(d) em=gd3 int=ultrafine

程序正常结束，虚频为-102.21，振动方向合理。

3.IRC计算
命令如下：
#p irc=(calcfc,maxpoint=150,lqa,stepsize=5) ub3lyp/6-31g(d) int=ultrafine geom=check guess=read
程序正常结束。IRC见图3

问题如下：IRC正常结束后看图3时，发现过渡态能量不在最高点，并且在一个方向上走了一步就停止了。
请问各位老师，出现该问题的原因是过渡态搜寻不够精确吗？可否给一些建议，谢谢各位。
关键输入输出文件已上传。

作者
Author: wzkchem5 时间: 2021-10-15 16:57
参见http://sobereva.com/400

作者
Author: sobereva 时间: 2021-10-16 02:53
你怎么能找TS时候加D3，跑IRC不加D3？这都不是同一个计算级别了
要用B3LYP-D3就对所有过程都用，要么都不用

作者
Author: wch990122 时间: 2021-10-20 17:06

sobereva 发表于 2021-10-16 02:53
你怎么能找TS时候加D3，跑IRC不加D3？这都不是同一个计算级别了
要用B3LYP-D3就对所有过程都用，要么都不 ...

谢谢老师，已经发现问题了。还有两个问题。
1.我对比了IRC加上D3和IRC没有加D3的两种情况，发现两种情况IRC生成物端的结构相同，对生成物端进行优化后的结构完全一致，那么下一步的过渡态还需要再重新算一遍吗（原来已经算过一遍了，算的时候加了D3）。
2.进行IRC计算时，用不用geom=check读取TSchk文件的区别大不大，直接输入过渡态结构进行计算会不会出现计算错误。

作者
Author: zjxitcc 时间: 2021-10-20 17:12
本帖最后由 zjxitcc 于 2021-10-20 17:13 编辑

wch990122 发表于 2021-10-20 17:06
谢谢老师，已经发现问题了。还有两个问题。
1.我对比了IRC加上D3和IRC没有加D3的两种情况，发现两种情况 ...

IRC时，自己复制进过渡态结构与从chk文件里读取结构，这两种做法是等价的，所以二者不会有区别。

但你这类反应，要时刻记得读取上一步波函数（guess=read）进行计算，若你不读取，高斯重新开始SCF迭代，很容易收敛到别的SCF解上去。假设你找过渡态后检验了波函数稳定性、确保波函数稳定，但你做IRC时忘记读取TS的稳定波函数进来，大概率会算出不稳定的、错误的态，导致后续全是错的。

作者
Author: sobereva 时间: 2021-10-21 01:01

wch990122 发表于 2021-10-20 17:06
谢谢老师，已经发现问题了。还有两个问题。
1.我对比了IRC加上D3和IRC没有加D3的两种情况，发现两种情况 ...

1 不管加D3的影响是明显还是可以忽略不计，横向对比的数据必须保证影响电子能量的所有因素完全统一

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