计算化学公社

标题: IRC得到的产物进行优化后结构变化很大 [打印本页]

作者
Author: 维维最爱小黄人 时间: 2016-12-13 20:45
标题: IRC得到的产物进行优化后结构变化很大
采用TS找到过渡态，初猜了很多种结构，最后得到的过渡态都一样，虚频振动也几乎一样，所以就觉得应该找对了过渡态，进行IRC计算，一端结构确实是反应物，然后把另一端的进行结构优化，但是优化后的结构变化很大，和虚频的振动方向也不一样，我以为是没有走到极小点，就将IRC的产物端设了60个点，结果45步就停止，得到的结构和刚开始用IRC计算得到的没多大差别，再进行结构优化，结构又变化很大，这可能是什么原因呢？如果取优化的结构作为产物，又和虚频振动的方向偏差很远。老师，我应该怎么做呢？图1是IRC得到的产物，图2是产物优化后的结构

作者
Author: sobereva 时间: 2016-12-14 00:16
如果IRC向产物一端你设了60个点但45个点就停了（提示已找到极小点），拿这个结构优化结构变化又很大，那么在优化的时候注意把步长上限设小，比如opt(calcfc,maxstep=5,notrust)，看看结果如何。
如果变化还很大，那说明IRC跑得不很精确，过早认为到了极小点，可以IRC的时候用calcall看看是否能跑得更长。

作者
Author: 维维最爱小黄人 时间: 2016-12-14 08:35
谢谢老师

，我试试

作者
Author: 维维最爱小黄人 时间: 2016-12-18 19:26

sobereva 发表于 2016-12-14 00:16
如果IRC向产物一端你设了60个点但45个点就停了（提示已找到极小点），拿这个结构优化结构变化又很大，那么 ...

老师，我将优化的步长变小，得到的结构还是变化很大，可是按照虚频振动的方向来看，产物应该是没有优化前的产物，即IRC得到的产物，是因为过渡态错了么？有没有可能IRC得到的产物和优化后的结构之间存在一个过渡态？之所以这样猜测，我是看到优化能量曲线先平稳一段，然后突然增加，然后又降低至平稳

作者
Author: sobereva 时间: 2016-12-18 19:51

维维最爱小黄人发表于 2016-12-18 19:26
老师，我将优化的步长变小，得到的结构还是变化很大，可是按照虚频振动的方向来看，产物应该是没有优化前 ...

没具体信息说不清楚
如果走IRC时候用了calcall而且IRC足够长，末尾的点不应当和极小点偏离多大

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