IRC问题求助

alan

root 发表于 2015-12-29 21:44
从NO+O’出发，随着N与O’距离靠近，势能一路下降就能到达NO2的势阱里，
也就是说，反应路径上不会出现 ...

你好，我最近做了relaxed scan然后求过渡态，虽然求出了一个虚频不大的过渡态，但是做IRC没有得到理想的结果，我怀疑我scan找的鞍点不对（by the way,我的柔性扫描没有正常结束，但看到一个弱弱的鞍点，就拿来求过渡态了）。可以麻烦你看看我的scan为什么错误结束吗？以及，图中倒数第二点是否可以认为是鞍点？

root

alan 发表于 2015-12-29 22:47
你好，我最近做了relaxed scan然后求过渡态，虽然求出了一个虚频不大的过渡态，但是做IRC没有得到理想的 ...

能量几乎是一直升高，到倒数第二个构型的时候，S与C的距离都已经超过4 $\AA$了，十有八九这个解离过程也是没有过渡态的。

最后一个点与倒数第二个构型相比在空间上的取向更稳定一些，你用倒数第二个构型做初始猜想找出来的过渡态，得到的也只会是连接了分子在空间相对取向上的两个局域极小值。

4 Angstrom已经可以算是解离产物了，再远算的意义不大。至于最后错误的结束，是因为最后一步的优化没结束，本质原因是两个碎片距离太远了，之间的相对作用力太小

alan

root 发表于 2015-12-29 23:15
能量几乎是一直升高，到倒数第二个构型的时候，S与C的距离都已经超过4 $\AA$了，十有八九这个解离过程也 ...

所以可以在文章里说明，这个断键通过relaxed scan后，不考虑过渡态，即不考虑此断键反应路径了，是吧？谢谢~~~（这个过渡态算了好久。。）

话说，顺便提一个奇怪的现象，我求的过渡态做IRC（两个方向同时做），得到的轨迹如下图。根本不是身为过渡态的一阶鞍点图形，这个IRC路径是不是也可以从侧面反应出此过渡态的不存在？谢谢！

root

alan 发表于 2015-12-29 23:30
所以可以在文章里说明，这个断键通过relaxed scan后，不考虑过渡态，即不考虑此断键反应路径了，是吧？谢 ...

是的，我觉得是没有过渡态的。

对于IRC，先不看这个曲线的形状，从两端的构型上看，就是空间取向改变了一下，能量略有下降。也就是说这个“过渡态”与解离反应无关。
再说形状，从IRC=0到附近的极大值点，能量变化其实很小，只有2个波数，是因为优化过渡态的时候收敛标准不严格，以及这个构型本来就比较软，在一般的优化步长下难以准确走到驻点。所以这个IRC也不算奇怪。

alan

root 发表于 2015-12-30 01:07
是的，我觉得是没有过渡态的。

对于IRC，先不看这个曲线的形状，从两端的构型上看，就是空间取向改变 ...

太谢谢你了~