过渡态正常收敛，有一个虚频且振动模式正确，但IRC有问题

查看数: 12419 | 评论数: 10 | 收藏 Add to favorites 2

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

发布时间: 2019-9-12 17:47

各位老师，我最近对一个过渡态进行irc计算，参考sob老师的博文，通过调整关键词，但发现IRC总是跑不对，请问是不是过渡态是有问题的，需要重新优化。我试了以下几种情况： ①我优化的过渡态对应的是一个氢转移反应 ...

回复 Reply

sobereva 发表于 Post on 2019-9-15 11:50:00

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

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

最爱喵星人 发表于 Post on 2019-9-13 21:41:55

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

1和3八成都 ...

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

wyy015020 发表于 Post on 2019-9-13 16:54:57

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

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

wyy015020 发表于 Post on 2019-9-13 12:56:24

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

1和3八成都 ...

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

wyy015020 发表于 Post on 2019-9-13 11:26:29

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

我也不知道呢，抱歉哈

wyy015020 发表于 Post on 2019-9-13 11:25:39

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 试试

wyy015020 发表于 Post on 2019-9-13 11:22:49

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

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

最爱喵星人 发表于 Post on 2019-9-13 08:01:48

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

sobereva 发表于 Post on 2019-9-13 06:57:35

2的结果就是对的，增加步数让两边跑长点，自然就像你要的反应物和产物了（除非漏了中间体）

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

Jasminer 发表于 Post on 2019-9-12 22:31:17

你第二次的计算结果挺好，需要的话增加步数就可以了。
为什么要改成第三次那样改算法加步长呢？

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