计算化学公社

标题: 过渡态正常收敛，有一个虚频且振动模式正确，但IRC有问题 [打印本页]

作者
Author: wyy015020 时间: 2019-9-12 17:47
标题: 过渡态正常收敛，有一个虚频且振动模式正确，但IRC有问题
各位老师，我最近对一个过渡态进行irc计算，参考sob老师的博文，通过调整关键词，但发现IRC总是跑不对，请问是不是过渡态是有问题的，需要重新优化。我试了以下几种情况：
①我优化的过渡态对应的是一个氢转移反应，虚频-500左右，优化过渡态的关键词是：# M062X/6-311G(d,p)  opt(calcall,ts,noeigen,maxstep=3,notrust)  freq scrf=smd  int=ultrafine，跑irc的关键词是：# M062X/6-311G(d,p) irc(calcall,maxpoint=20,stepsize=5) int=ultrafine  scrf=smd  guess=read，得到IRC图形为： (, 下载次数 Times of downloads: 30) ，左边端点是和预期一样的反应物的结构，右边端点也接近于产物的结构，但irc图中靠近过渡态的几个点转折太大了，所以我接下来调整步长试了一下。
②跑irc的关键词是：# M062X/6-311G(d,p) irc(calcall,maxpoint=20,stepsize=3)  int=ultrafine  scrf=smd  guess=read，得到的IRC图形为： (, 下载次数 Times of downloads: 26) ,看irc走的趋势是对的，是氢原子从反应物向产物运动。但因为两端结构离反应物和产物的结构相差太远，我又把步长和步数调整了一下。
③跑irc的关键词是：# M062X/6-311G(d,p) irc(calcfc,maxpoint=50,stepsize=10,LQA) int=ultrafine  scrf=smd  guess=read，总共跑了101个点，步长调到了10，得到的IRC图形为： (, 下载次数 Times of downloads: 25) ，左边端点结构是想要的反应物结构，右边最开始几步是往产物的方向跑的，但是最后几步就开始向反应物方向跑了。
现在就很疑惑，为什么把步长或者步数变大，IRC结构就变得特别奇怪。这是不是因为我这个过渡态不准确导致的？我是不是要重新优化过渡态呢？望各位大佬不吝赐教。（注：在IRC任务中加guess=read关键词是因为不加的时候IRC任务第一次输出的SCF Done能量和找过渡态最后一步的SCF Done能量明显不同）

作者
Author: Jasminer 时间: 2019-9-12 22:31
你第二次的计算结果挺好，需要的话增加步数就可以了。
为什么要改成第三次那样改算法加步长呢？

作者
Author: sobereva 时间: 2019-9-13 06:57
2的结果就是对的，增加步数让两边跑长点，自然就像你要的反应物和产物了（除非漏了中间体）

1和3八成都是因为步长还不够小，导致TS旁边的点的波函数没收敛到对应TS的状态，引起势能面突变。要么你也可以尝试把forward和reverse分别跑每次都用guess=read

作者
Author: 最爱喵星人 时间: 2019-9-13 08:01
请问如果步长如此影响irc曲线，那怎么确定得到irc是正确的呢？我的irc跑完后反应物与过渡态能量很相近，请问是怎么回事呢？

作者
Author: wyy015020 时间: 2019-9-13 11:22

Jasminer 发表于 2019-9-12 22:31
你第二次的计算结果挺好，需要的话增加步数就可以了。
为什么要改成第三次那样改算法加步长呢？

好的，谢谢老师，我把步数增加到maxpoint=50试试

作者
Author: wyy015020 时间: 2019-9-13 11:25

sobereva 发表于 2019-9-13 06:57
2的结果就是对的，增加步数让两边跑长点，自然就像你要的反应物和产物了（除非漏了中间体）

1和3八成都 ...

好的，谢谢sob老师，我现在把关键词设为# M062X/6-311G(d,p) irc(calcfc,maxpoint=30,stepsize=5,forward) int=ultrafine scrf=smd guess=read 试试

作者
Author: wyy015020 时间: 2019-9-13 11:26

最爱喵星人发表于 2019-9-13 08:01
请问如果步长如此影响irc曲线，那怎么确定得到irc是正确的呢？我的irc跑完后反应物与过渡态能量很相近，请 ...

我也不知道呢，抱歉哈

作者
Author: wyy015020 时间: 2019-9-13 12:56

sobereva 发表于 2019-9-13 06:57
2的结果就是对的，增加步数让两边跑长点，自然就像你要的反应物和产物了（除非漏了中间体）

1和3八成都 ...

sob老师，我还想请教一个问题，就是我想对这个过渡态进行单点能计算，是不是也应该加上guess=read呢？

作者
Author: wyy015020 时间: 2019-9-13 16:54

wyy015020 发表于 2019-9-13 12:56
sob老师，我还想请教一个问题，就是我想对这个过渡态进行单点能计算，是不是也应该加上guess=read呢？

我试了一下直接把过渡态的log文件保存成单点能的输入文件，不加guess=read和加guess=read关键词，这个过渡态的单点能差了13kcal/mol

作者
Author: 最爱喵星人 时间: 2019-9-13 21:41

sobereva 发表于 2019-9-13 06:57
2的结果就是对的，增加步数让两边跑长点，自然就像你要的反应物和产物了（除非漏了中间体）

1和3八成都 ...

请问如果步长如此影响irc曲线，那怎么确定得到irc是正确的呢？我的irc跑完后反应物与过渡态能量很相近，请问是怎么回事呢？

作者
Author: sobereva 时间: 2019-9-15 11:50

最爱喵星人发表于 2019-9-13 21:41
请问如果步长如此影响irc曲线，那怎么确定得到irc是正确的呢？我的irc跑完后反应物与过渡态能量很相近， ...

相关问题在此文里解释得很清楚
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

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