做IRC时过渡态不在能量最高点

wch990122 · 发表于 Post on 2021-10-15 16:15:47

各位老师好，我在寻找图一两个构型间的过渡态时遇到了问题。我的步骤如下：
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时，发现过渡态能量不在最高点，并且在一个方向上走了一步就停止了。
请问各位老师，出现该问题的原因是过渡态搜寻不够精确吗？可否给一些建议，谢谢各位。
关键输入输出文件已上传。

wzkchem5 · 发表于 Post on 2021-10-15 16:57:54

参见http://sobereva.com/400

sobereva · 发表于 Post on 2021-10-16 02:53:16

你怎么能找TS时候加D3，跑IRC不加D3？这都不是同一个计算级别了
要用B3LYP-D3就对所有过程都用，要么都不用

wch990122 · 发表于 Post on 2021-10-20 17:06:34

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

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

zjxitcc · 发表于 Post on 2021-10-20 17:12:30

本帖最后由 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的稳定波函数进来，大概率会算出不稳定的、错误的态，导致后续全是错的。

sobereva · 发表于 Post on 2021-10-21 01:01:49

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 做IRC时过渡态不在能量最高点