计算化学公社

标题: irc两端结构优化收敛问题 [打印本页]

作者
Author: jing-quan63 时间: 2019-4-15 21:10
标题: irc两端结构优化收敛问题
各位老师好：
附件是我进行IRC计算得到附件中的曲线，其中的两端结构以及过渡态结构都符合预期，计算所用命令为：#p IRC=(CALCFC,STEPSIZE=5,MAXPOINT=25) M062X/6-31g(d,p)，采用相同方法和基组，取右端反应物结构进行几何优化时，出现SCF不收敛的情况，我的问题如下： 1：对irc曲线端点结构进行几何优化时，采用的方法与基组是否必须与IRC计算时相同？
2：我按照帖子中的方法，尝试过加上int=ultrafine,或者更改为B3LYP泛函进行优化，能够解决收敛问题，请问这时的结构是否有意义，势能面是否不同了呢？

作者
Author: sobereva 时间: 2019-4-16 00:15
1 原理上不需要。但如果用的不同，道理上说不通
2 此时相当于和IRC不在一个势能面上，如果要把IRC的数据和极小点的数据拿到一起说事，说不通

作者
Author: jing-quan63 时间: 2019-4-16 08:36

sobereva 发表于 2019-4-16 00:15
1 原理上不需要。但如果用的不同，道理上说不通
2 此时相当于和IRC不在一个势能面上，如果要把IRC的数据和 ...

那请问老师，我可以不用IRC两端结构来优化吗，比如我重新再gview里构建反应物初始结构，再选用与IRC计算相同的理论方法进行优化，这样可以吗？但我感觉这样得到的势垒就不在一个反应路径上了。

作者
Author: 哈哈哈 时间: 2019-4-16 08:39

sobereva 发表于 2019-4-16 00:15
1 原理上不需要。但如果用的不同，道理上说不通
2 此时相当于和IRC不在一个势能面上，如果要把IRC的数据和 ...

那遇到2这种情况如何处理呢？（类似于楼主的情况，优化的右端反应物有虚频，加tight int=ultrafine 优化后解决。但和IRC就不在一个势能面上了）
a 先int=ultrafine 优化rc 再拿该结构去int=ultrafine优化，再和IRC数据做比较。（这种有时去int=ultrafine 能得到正常的结构，但有时仍有虚频）
b ts 更精确所有的都在int=ultrafine 下做优化
我想了这两种，不确定能不能这样做。（但除了当下优化的这部分还有其他正常优化的数据，一起比较感觉还是有些不合理）

作者
Author: sobereva 时间: 2019-4-16 11:12

jing-quan63 发表于 2019-4-16 08:36
那请问老师，我可以不用IRC两端结构来优化吗，比如我重新再gview里构建反应物初始结构，再选用与IRC计算 ...

IRC两端的结构肯定比你自己随便摆的更接近于极小点，这么做没实际意义

作者
Author: sobereva 时间: 2019-4-16 11:13

哈哈哈发表于 2019-4-16 08:39
那遇到2这种情况如何处理呢？（类似于楼主的情况，优化的右端反应物有虚频，加tight int=ultrafine 优化 ...

找过渡态、走IRC也用ultrafine
如果需要和IRC上的点对比能量，至少优化完了重新用ultrafine算个单点

作者
Author: jing-quan63 时间: 2019-4-16 12:51

sobereva 发表于 2019-4-16 11:13
找过渡态、走IRC也用ultrafine
如果需要和IRC上的点对比能量，至少优化完了重新用ultrafine算个单点

好的，谢谢sober老师，非常感谢您解惑！

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