|
|
各位老师，我最近对一个过渡态进行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图形为：
，左边端点是和预期一样的反应物的结构，右边端点也接近于产物的结构，但irc图中靠近过渡态的几个点转折太大了，所以我接下来调整步长试了一下。
②跑irc的关键词是：# M062X/6-311G(d,p) irc(calcall,maxpoint=20,stepsize=3) int=ultrafine scrf=smd guess=read，得到的IRC图形为：
,看irc走的趋势是对的，是氢原子从反应物向产物运动。但因为两端结构离反应物和产物的结构相差太远，我又把步长和步数调整了一下。
③跑irc的关键词是：# M062X/6-311G(d,p) irc(calcfc,maxpoint=50,stepsize=10,LQA) int=ultrafine scrf=smd guess=read，总共跑了101个点，步长调到了10，得到的IRC图形为：
，左边端点结构是想要的反应物结构，右边最开始几步是往产物的方向跑的，但是最后几步就开始向反应物方向跑了。
现在就很疑惑，为什么把步长或者步数变大，IRC结构就变得特别奇怪。这是不是因为我这个过渡态不准确导致的？我是不是要重新优化过渡态呢？望各位大佬不吝赐教。（注：在IRC任务中加guess=read关键词是因为不加的时候IRC任务第一次输出的SCF Done能量和找过渡态最后一步的SCF Done能量明显不同）
|
|
发表于 Post on 2019-9-12 17:47:40
|
只看该作者 Only view this author
|只看大图 View big picture
|倒序浏览 Reverse view
|阅读模式 Reading model
收藏 Add to favorites
发表于 Post on 2019-9-13 06:57:35
楼主 Author