计算化学公社

标题: 类似石墨烯结构还原NO的反应，柔性扫描能量一直下降 [打印本页]

作者
Author: wch990122 时间: 2021-9-4 11:09
标题: 类似石墨烯结构还原NO的反应，柔性扫描能量一直下降
各位老师好：
我再计算一个类似石墨烯结构还原NO的反应时，遇到了如下问题。
首先确定了上一步反应的过渡态，并通过IRC进行了验证，IRC正常结束，并且连接着正确的反应物，想用IRC中得到的产物进行下一步的计算。下一步应该是N2在结构上的解吸，根据文献推断图一中的C-N键会断裂，于是使用柔性扫描C-N键长，但是在柔性扫描的过程中，却发现随着键长增加，能量始终在下降。
我想请问一下：这种情况是说明这个过程没有过渡态吗？在柔性扫描中能量始终下降，为何在IRC中计算到图一的产物就自动结束了呢？能量不是应该继续减小吗？

作者
Author: wch990122 时间: 2021-9-4 11:12
抱歉，IRC中反应物和产物标反了，应该是这样。

作者
Author: sobereva 时间: 2021-9-4 13:00

wch990122 发表于 2021-9-4 11:12
抱歉，IRC中反应物和产物标反了，应该是这样。

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了

作者
Author: sobereva 时间: 2021-9-4 13:04
IRC两端根本不对应于精确的极小点结构，哪怕IRC跑得再长，下文说了。必须对两端进行优化
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）
优化完了如果确实得到N2没有掉下去的极小点到时候再说

作者
Author: wch990122 时间: 2021-9-4 15:37

sobereva 发表于 2021-9-4 13:00
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员 ...

抱歉老师，我下次注意。

作者
Author: wch990122 时间: 2021-9-4 15:39

sobereva 发表于 2021-9-4 13:04
IRC两端根本不对应于精确的极小点结构，哪怕IRC跑得再长，下文说了。必须对两端进行优化
在Gaussian中计算 ...

老师，上面提供的图1就是优化完成后的结构。
这是优化关键词
#p opt freq ub3lyp/6-31g(d) em=gd3

作者
Author: wzkchem5 时间: 2021-9-4 19:55

wch990122 发表于 2021-9-4 08:39
老师，上面提供的图1就是优化完成后的结构。
这是优化关键词
#p opt freq ub3lyp/6-31g(d) em=gd3

确认产物没有虚频了吗？可能沿着IRC以及之后你优化找到的不是产物的稳定结构，而是过渡态。对于分叉反应，走一遍IRC不一定走到正确的产物，可能得到另一个过渡态，从那个过渡态继续走IRC才能得到产物

作者
Author: wch990122 时间: 2021-9-4 20:54

wzkchem5 发表于 2021-9-4 19:55
确认产物没有虚频了吗？可能沿着IRC以及之后你优化找到的不是产物的稳定结构，而是过渡态。对于分叉反应 ...

确定过了，没有虚频。

作者
Author: sobereva 时间: 2021-9-5 01:25
柔性扫描的时候guess=read读取之前做opt freq无虚频任务的chk再试，如果延续了之前任务的波函数，理应不会能量直接下降，否则优化的时候直接就会优化到其它结构去

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