计算化学公社

标题: 使用PSI4和sobEDA做能量分解，尝试了许多方法SCF都不收敛 [打印本页]

作者
Author: NGC626 时间: 2024-1-18 12:22
标题: 使用PSI4和sobEDA做能量分解，尝试了许多方法SCF都不收敛
本帖最后由 NGC626 于 2024-1-18 12:32 编辑

各位老师好，我建立了一个碳原子在几层铱原子上吸附的模型，并用CP2K在PBE/DZVP-MOLOPT-SR-GTH下做了结构优化，之后我尝试用了PSI4和sobEDA对它做能量分解，但是一直提示不收敛，我查阅并尝试了sob老师的博文(http://sobereva.com/61)中的内容，以及帖子(http://bbs.keinsci.com/thread-15720-1-1.html、http://bbs.keinsci.com/thread-6995-1-1.html、http://bbs.keinsci.com/thread-4570-1-1.html)中的方法，折腾了快一个月都没有解决，以下是我用CP2K对该结构做opt和计算能量的输入输出文件，以及用sobEDA和PSI4做能量分解的输入输出文件，恳请各位老师指导。

作者
Author: zjxitcc 时间: 2024-1-18 14:33
（1）这不是任何一种EDA的问题，是建模问题，浪费了快一个月时间。一开始压根就不应该提交这种计算，即使侥幸收敛了也没有意义。底下的Ir原子数太少了，仔细阅读《使用量子化学程序基于簇模型计算金属表面吸附问题》http://sobereva.com/540。
（2）碳原子基态是三重态，你现在把它当单重态算，完全错了。退一步讲，就算“碳原子吸附在Ir原子表面”这个复合物体系基态是单重态，你也得给出理由来。

作者
Author: NGC626 时间: 2024-1-19 13:00

zjxitcc 发表于 2024-1-18 14:33
（1）这不是任何一种EDA的问题，是建模问题，浪费了快一个月时间。一开始压根就不应该提交这种计算，即使侥 ...

老师您好，我尝试了重新建模，但是又发生了报错，请问可能是什么原因导致的呢？输入输出文件如下

作者
Author: zjxitcc 时间: 2024-1-19 15:34

NGC626 发表于 2024-1-19 13:00
老师您好，我尝试了重新建模，但是又发生了报错，请问可能是什么原因导致的呢？输入输出文件如下

如果计算本身不要求，但你自己强行加了scf(conver=6)和IOp(7/127=-99)，表示您已误入歧途。其他问题：
（1）int=superfine没有必要
（2）def2SVP基组做能量分解没有太大的意义，你这体系全用def2TZVP也能算的动
（3）没有加色散校正，急需阅读《谈谈“计算时是否需要加DFT-D3色散校正？”》http://sobereva.com/413
（4）最好用一个杂化泛函，当前是纯泛函
（5）SCF不收敛按此文解决《解决SCF不收敛问题的方法》http://sobereva.com/61
在（1）~（4）都解决之后，如果你遇到（5）的问题解决不了，可再来提问。待你成功算出结果后，需用对称破缺DFT再算一遍，看结果是否有所变化。

作者
Author: NGC626 时间: 2024-1-21 13:05
本帖最后由 NGC626 于 2024-1-22 11:00 编辑

zjxitcc 发表于 2024-1-19 15:34
如果计算本身不要求，但你自己强行加了scf(conver=6)和IOp(7/127=-99)，表示您已误入歧途。其他问题：
...

老师您好，我尝试了更换参数，但是SCF仍然不收敛，输入输出文件如下

作者
Author: wzkchem5 时间: 2024-1-21 23:29

NGC626 发表于 2024-1-21 06:05
老师您好，我尝试了更换参数，但是计算发生了报错，请问可能是什么原因导致的呢？输入输出文件如下

以后不要用“计算发生了报错”这种措辞，得先自己分析是什么类别的报错，并且在问问题的时候指明这个类别。这样显得你思考了，问问题有诚意。
我相信你不会不知道以下两行是说明这个计算是因为什么原因报错了
>>>>>>>>>> Convergence criterion not met.
SCF Done: E(RPBE1PBE) = -5011.01716865 A.U. after 129 cycles
我也相信你不会不知道看sob老师的哪篇博文解决这个问题（或者说在sob老师博客的搜索框里搜索什么关键词）

作者
Author: NGC626 时间: 2024-1-22 09:45
本帖最后由 NGC626 于 2024-1-22 10:59 编辑

wzkchem5 发表于 2024-1-21 23:29
以后不要用“计算发生了报错”这种措辞，得先自己分析是什么类别的报错，并且在问问题的时候指明这个类别 ...

谢谢老师提醒，和此前的错误还是一样，SCF不收敛，已经尝试了许多论坛里的方法和sob老师帖子里的建议，折腾了很久

作者
Author: NGC626 时间: 2024-1-22 11:47
本帖最后由 NGC626 于 2024-1-22 12:02 编辑

wzkchem5 发表于 2024-1-21 23:29
以后不要用“计算发生了报错”这种措辞，得先自己分析是什么类别的报错，并且在问问题的时候指明这个类别 ...

老师您好，我想到一个方法，就是先用高斯或者ORCA算一次这个体系的单点，收敛之后，再做能量分解，请问这样的方法是否具有合理性呢？

作者
Author: wzkchem5 时间: 2024-1-22 16:36

NGC626 发表于 2024-1-22 04:47
老师您好，我想到一个方法，就是先用高斯或者ORCA算一次这个体系的单点，收敛之后，再做能量分解，请问这 ...

你本来也要算这个体系的单点啊，不算单点怎么做能量分解，只不过是你自己跑一个单点和sobEDA替你跑一个单点的区别

作者
Author: NGC626 时间: 2024-1-22 18:31
本帖最后由 NGC626 于 2024-1-22 18:34 编辑

wzkchem5 发表于 2024-1-22 16:36
你本来也要算这个体系的单点啊，不算单点怎么做能量分解，只不过是你自己跑一个单点和sobEDA替你跑一个单 ...

老师好，我用ORCA在PBE0/def2SV(P)下算了一下这个体系的单点，发现是可以收敛的，但是用高斯在同样的级别在下算却不能收敛，在高斯里面换成def2TZVP之后仍然不能收敛，以下是这两次计算的输入输出文件，请问可能是recalc_hess的原因吗？如果是的话，如何才能在高斯里算单点的时候用类似的功能呢？如果不是的话，在ORCA里能收敛的根本原因是什么呢？

作者
Author: wzkchem5 时间: 2024-1-22 20:04

NGC626 发表于 2024-1-22 11:31
老师好，我用ORCA在PBE0/def2SV(P)下算了一下这个体系的单点，发现是可以收敛的，但是用高斯在同样的级别 ...

如果经检查高斯的输入文件没有问题（尤其注意赝势等），尝试用Multiwfn或MOKIT把ORCA的波函数传给高斯作为初猜。
一个软件能不能让一个体系能收敛，有一定的偶然性，不是所有情况都能说出为什么

作者
Author: NGC626 时间: 2024-1-23 11:07

wzkchem5 发表于 2024-1-22 20:04
如果经检查高斯的输入文件没有问题（尤其注意赝势等），尝试用Multiwfn或MOKIT把ORCA的波函数传给高斯作 ...

谢谢老师的建议，我在论坛搜索之后，找到了把高斯收敛的波函数作为ORCA的初猜的方法，请问如何才能把ORCA的波函数作为高斯的初猜呢？

作者
Author: wzkchem5 时间: 2024-1-23 16:45

NGC626 发表于 2024-1-23 04:07
谢谢老师的建议，我在论坛搜索之后，找到了把高斯收敛的波函数作为ORCA的初猜的方法，请问如何才能把ORCA ...

仔细读http://sobereva.com/379
把orca波函数用orca_2mkl转成molden格式，读进Multiwfn（或MOKIT），导出成fch，用高斯的unfchk转换成chk格式，用guess=read读取

作者
Author: NGC626 时间: 2024-1-23 18:03

wzkchem5 发表于 2024-1-23 16:45
仔细读http://sobereva.com/379
把orca波函数用orca_2mkl转成molden格式，读进Multiwfn（或MOKIT），导 ...

谢谢老师

作者
Author: gaojiahao 时间: 2024-11-25 16:08
老师，您好。我再用sobEDA进行计算时，因为牵扯到Am,Cm等元素。所以我没有采用D3校正。但在计算Pu元素时，我用b3lyp方法和混合基组进行计算时，加入D3校正可以计算出结果，不加D3却不收敛。请问这是什么原因呢？

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