计算化学公社

标题: 求教结果文件出错Operation on file out of range. [打印本页]

作者
Author: 小学徒 时间: 2021-4-8 18:20
标题: 求教结果文件出错Operation on file out of range.
本帖最后由小学徒于 2021-4-8 18:20 编辑

优化分子的s1态，取用的基态优化好的结构且基态结构无虚频。s1态第一遍优化的输入文件#p opt td(nstates=5) b3lyp/6-311g(d,p)
显示出现震荡，优化不收敛

取能量最低结构，修改命令行
#p opt(gdiis,maxstep=5,notrust) td=(nstates=5) b3lyp/6-311g(d,p)
结果“Convergence failure -- run terminated.”显示SCF不收敛
看了一下老师的这篇文章http://sobereva.com/61，重新修改了命令行，
#p Stable=Opt SCF=QC td=(nstates=5) b3lyp/6-311g(d,p)
结果文件出错，显示“Operation on file out of range.”

想要请教一下，遇到这种情况要如何修改？（分别用g09和g16都计算过，都出现了最后的问题）
附件中是这几次计算的结果文件

作者
Author: sobereva 时间: 2021-4-10 06:30
stable和TD根本就没法一起用

尝试#p opt(gdiis,maxstep=5,notrust) td=(nstates=5) b3lyp/6-311g(d,p)是得当的，SCF不收敛该怎么解决就怎么解决，该尝试的都需要尝试，仔细看博文开头的红字。
文中都已经明确说了不要随便用SCF=QC干嘛还用

上传相关文件时说清楚每个文件都是对应什么的，没人愿意一个个去看然后去猜

作者
Author: 小学徒 时间: 2021-4-13 17:42
好的谢谢老师，针对SCF不收敛我换了降低SCF收敛限的方法进行计算，命令行如下
#p opt SCF=conver=6 td=(nstates=5) b3lyp/6-311g(d,p)
计算结果出现震荡，位移低于收敛限，但是受力高于收敛限，我想问一下老师这种结果是说明我修改优化不收敛活着SCF不收敛的方法不对吗？还是分子的结构有问题？
我后续要如何改进？相关文件已重新上传，文件为每一步的结果文件，文件名字已经改为每一步的目的和操作

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