计算化学公社

标题: 高斯输入文件信息咨询 [打印本页]

作者
Author: 北极天99 时间: 2015-6-3 21:01
标题: 高斯输入文件信息咨询

本人初学者，是一枚菜鸟，问老师们几个基础的问题：

1.我在用qst3计算过渡态，开始的时候命令是这样写的# opt=(calcall,qst3) ub3lyp/6-311+g(d,p) geom=connectivity但是途中L502错误，我就看sob老师给出的方法加上了int=ultrafine即输入文件这样写# opt=(calcall,qst3) ub3lyp/6-311+g(d,p) int=ultrafine geom=connectivity最后是计算正常结束了，计算出来的东西没有虚频。是不是证明我的输入命令哪里出问题了？
2.我计算的时候经常出现SCF不收敛的情况，我看到过sob老师的帖子，所以我基本上都是加上int=ultrafine，还有或者是加上SCF=vshift=300.我就是不太清楚，像这样的命令应该加在哪个地方？是不是需要加在opt的括号内？
3.对于开壳层和闭壳层问题的计算，有的时候体系可能是开壳层，比如用mp2方法计算，我命令直接写了mp2，忘记加上u了，或者有的我不确定到底是开壳层或者闭壳层，一律用的是mp2，而不是ump2，请问老师，这样是不是计算出来的结果就不准了?或者说，我如果不加u，体系是开壳层的情况下，高斯程序会不会自动判断出来是开壳层，计算的时候自动按照开壳层计算？

多谢多谢，感激不尽！

作者
Author: 北极天99 时间: 2015-6-3 22:00
4.我还有一个计算同一个过渡态的，命令这样写的# opt=qst3 ub3lyp/6-31+g(d,p) geom=connectivity，最后出现了错误，请教老师，我应该怎么解决啊？非常感谢！
(, 下载次数 Times of downloads: 22)

作者
Author: sobereva 时间: 2015-6-4 07:56
1 不建议用QST3，比较麻烦，比起opt=TS也并无明显优势
2 比如# M062X/6-311G** int=ultrafine scf=vshift=300
3 对于自旋多重度>1的体系，写MP2、HF、B3LYP等等，自动就会默认为UMP2、UHF、UB3LYP。
自旋多重度=1的体系（单重态），高斯默认用闭壳层方式计算。但这样的体系也可能有开壳层的情况（诸如双自由基），要想合理计算就必须明确写U了（还得同时结合guess=mix之类关键词破坏初猜轨道对称性才起作用）。
4 到最后已经发生震荡难以收敛了，最终超过了步数上限。此时建议拿出最后的结构，用opt(calcall,TS)继续做，还不行再把gdiis,maxstep=4,notrust加进opt里。

PS: 略大的输出文件（>3MB）建议压缩一下再作为附件发，这样可以节省论坛储存空间，下载也更快，我帮你更新了。

作者
Author: 北极天99 时间: 2015-6-4 09:19
好的，非常感谢老师，我再计算一下试试。

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