计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › 做IRC时过渡态不在能量最高点

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

查看数: 6426 | 评论数: 5 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

wch990122

发布时间: 2021-10-15 16:15

正文摘要:

各位老师好，我在寻找图一两个构型间的过渡态时遇到了问题。我的步骤如下： 1.首先对键角进行柔性扫描。柔性扫描命令如下： #p opt=modredundant ub3lyp/6-31g(d) nosymm geom=check guess=read 扫描结果如图2 ...

回复 Reply

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

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

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

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

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文件的区别大不大，直接输入过渡态结构进行计算会不会出现计算错误。

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

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

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

参见http://sobereva.com/400

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