sophieding 发表于 2020-10-16 14:29
好的，感谢指教。IRC关键词写这个主要是过渡态优化也有这个，所以一起写上了，后来改为irc=(calcfc,recal ...

嗯，刚说错了，irc里有recalc，只是没有gdiis

wzkchem5 发表于 2020-10-16 09:49
IRC后面不应该写(recalc=5,gdiis)。recalc和gdiis都是opt的专属选项，IRC没有这两个选项，所以报输入文件 ...

好的，感谢指教。IRC关键词写这个主要是过渡态优化也有这个，所以一起写上了，后来改为irc=(calcfc,recalc=5)已可运算。反应物的结构已进行了优化~

sophieding 发表于 2020-10-16 09:07
谢谢老师的指导。1.之前已经采用PM6对C5-C5进行了柔性扫描，得出的总能量曲线一直往下降。问过相关的老师 ...

IRC后面不应该写(recalc=5,gdiis)。recalc和gdiis都是opt的专属选项，IRC没有这两个选项，所以报输入文件语法错误。
你优化了反应物结构吗？可能过渡态和反应物结构非常接近，或甚至即使在DFT水平下反应也是没有势垒的。

wzkchem5 发表于 2020-10-15 22:38
1. 把初始结构的C5-C6缩短一些，再试。一般而言，如果过渡态优化收敛到类似产物的结构，那么把初猜结构往反 ...

谢谢老师的指导。1.之前已经采用PM6对C5-C5进行了柔性扫描，得出的总能量曲线一直往下降。问过相关的老师，老师说可能是过渡态能垒比较低，导致能垒曲线一直往下降。查过相关文献，这个分子的类似物在断C5-C6的能垒时确实很低。
2.IRC输入文件已上传，望老师帮忙查看。
3.不好意思，之前写错了，应该是C3-O4。想做C3-O4均裂为两个自由基键解离能的能垒图，与C5-C6断键的做一个能垒图作为比较，探究哪个路径优先。

1. 把初始结构的C5-C6缩短一些，再试。一般而言，如果过渡态优化收敛到类似产物的结构，那么把初猜结构往反应物方向变一变，反之则往产物方向变。如果仍然不行，可以对C5-C6键做一个relaxed scan，取最高点作为初猜结构。
2. 上传IRC文件的gjf，我们看看。L1就报错一般是因为关键词语法没写对。
3. 你指的是C3-O4吗？你会做C5-C6，那C3-O4也是同理啊