计算化学公社

标题: 用QC方法解决SCF不收敛，stable关键词的应用问题？ [打印本页]
作者
Author: 一个平凡之人    时间: 2020-11-12 09:10
标题: 用QC方法解决SCF不收敛，stable关键词的应用问题？
SCF不收敛的输入文件已经上传，请各位老师帮忙分析一下输入文件中有什么问题，导致不收敛，谢谢各位老师。关于sob老师提到的IOpt的选项是在高斯09用户参考手册的第224页。

作者
Author: wzkchem5    时间: 2020-11-12 10:33
需要贴输出文件
很多时候光看输入文件看不出来，除非重新跑一遍（因为SCF不收敛并不必然意味着输入文件写错了，没有人能只看输入文件就预测这个任务能不能收敛）。而既然你已经跑过了，最好连输出文件一起贴
作者
Author: 一个平凡之人    时间: 2020-11-12 11:52
wzkchem5 发表于 2020-11-12 10:33
需要贴输出文件
很多时候光看输入文件看不出来，除非重新跑一遍（因为SCF不收敛并不必然意味着输入文件写 ...

现在正在运行加了QC的文件，等它运行完有问题我再贴出来。谢谢您的建议。
作者
Author: zjxitcc    时间: 2020-11-12 12:41
本帖最后由 zjxitcc 于 2020-11-12 12:52 编辑

你没（认真）看上次回复时推荐的两篇博文。即使是预优化，B3LYP/STO-3G也是过于粗糙 以至于没有意义。加QC的任务也不用算了，建议直接kill掉。

你这个问题的合理输入文件示例（y, z值请根据自己机器的实际配置修改）
  1. %chk=xxx.chk
  2. %mem=yGB
  3. %nprocshared=z
  4. #p UB3LYP/6-31G(d) em=GD3BJ nosymm stable=opt

  6. title

  8. -3 3
  9. [坐标]

  11. --Link1--
  12. %chk=xxx.chk
  13. %mem=yGB
  14. %nprocshared=z
  15. #p opt=calcfc freq UB3LYP/6-31G(d) em=GD3BJ nosymm guess=read geom=allcheck
复制代码

另外，建议到linux做计算，不要在Windows下算。不要怕接触linux，安装VMware虚拟机 或 Windows的Ubuntu子系统 十分简单，在里面跑linux版高斯，节约时间，时间=生命。安装教程很多，如《在VMware 15中安装CentOS 7.6的完整过程视频演示》http://sobereva.com/454
再另外，如果不支持em=GD3BJ，说明你的高斯版本过于老旧。
作者
Author: 一个平凡之人    时间: 2020-11-12 13:05
zjxitcc 发表于 2020-11-12 12:41
你没（认真）看上次回复时推荐的两篇博文。即使是预优化，B3LYP/STO-3G也是过于粗糙 以至于没有意义。加QC ...

好的，谢谢您，我会认真考虑您的建议。
作者
Author: sobereva    时间: 2020-11-13 03:22
本来STO-3G算中性都没法忍受，算带-3这么多负电荷体系还用STO-3G简直是不可思议。
仔细把此文看了
谈谈量子化学中基组的选择
http://sobereva.com/336（http://bbs.keinsci.com/thread-3545-1-1.html）

也不清楚你这个体系到底是算的什么问题，我都怀疑这个古怪的模型本身都没法体现实际要研究的问题。普遍是越古怪的结构SCF越难收敛。初学者遇到SCF不收敛问题时，在使用我博文里提到的那些解决SCF不收敛关键词之前，先好好考虑一下当前的结构、电荷、自旋多重度设定等等是否确实有意义，要不然就算收敛了也是白算。

手册里有IOp，也有opt，但绝对没有IOpt。描述的时候一定要严格准确。

作者
Author: zjxitcc    时间: 2020-11-13 09:38
看了sob老师的回答，也想提醒你一句，检查你的电荷设置是否有问题，感觉大概率有误。
作者
Author: 一个平凡之人    时间: 2020-11-13 14:36
sobereva 发表于 2020-11-13 03:22
本来STO-3G算中性都没法忍受，算带-3这么多负电荷体系还用STO-3G简直是不可思议。
仔细把此文看了
谈谈量 ...

嗯嗯，谢谢您的解答，我再认真考虑一下我的结构。
作者
Author: 一个平凡之人    时间: 2020-11-13 14:36
zjxitcc 发表于 2020-11-13 09:38
看了sob老师的回答，也想提醒你一句，检查你的电荷设置是否有问题，感觉大概率有误。

嗯嗯，好的，谢谢您的提醒，我再认真考虑一下。



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3