计算化学公社

标题: 求助能量分解的问题，计算出的数很大 [打印本页]
作者
Author: 小翟w    时间: 2024-4-2 09:46
标题: 求助能量分解的问题，计算出的数很大
老师们好，我用xeda提交inp文件计算能量分解，计算出的数很大，下面是inp文件，out文件，想请老师帮忙看看，这种问题要怎么解决，谢谢老师！
作者
Author: zjxitcc    时间: 2024-4-2 14:07
本帖最后由 zjxitcc 于 2024-4-2 14:20 编辑

很简单，至少有一处SCF不收敛。GKS-EDA对于3个片段的复合物需进行7次不同的SCF计算，任何一次SCF不收敛都会导致这个结果。这个问题是GAMESS的老毛病，甚至GKS-EDA开发者也深受其扰，在这个论坛上的GAMESS-US板块你可以找到无数这种帖子http://bbs.keinsci.com/forum.php ... ter=typeid&typeid=9

（几乎）彻底解决这个问题的办法，也是基本上解决上面所有这类帖子问题的办法：用frag_guess_wfn小程序生成GKS-EDA输入文件
我们写一个gjf文件，例如MO24.gjf，内容如下
  1. %mem=128GB
  2. %nprocshared=64
  3. #p wB97XD/def2TZVP guess(fragment=3)

  5. {gks}

  7. 0 1 0 1 0 1 0 1
  8. C(fragment=1)   -0.00101100    0.00124700   -0.43694500
  9. C(fragment=1)   -0.00136600    0.01059200    0.89704900
  10. H(fragment=1)   -0.01620900    0.92064900   -1.00842900
  11. H(fragment=1)    0.01491600   -0.92604200   -0.99551100
  12. H(fragment=1)    0.00802700   -0.91192000    1.46254900
  13. H(fragment=1)   -0.01130000    0.94073400    1.44981300
  14. Mo(fragment=2)    3.08135100   -0.00290900   -0.08727100
  15. F(fragment=2)    2.67718100    1.59537400    0.75412500
  16. F(fragment=2)    2.63018600   -1.75611000    0.29835400
  17. O(fragment=2)    4.47609500    0.09568300   -0.95496500
  18. Mo(fragment=3)   -3.08087600    0.00152100   -0.08735400
  19. F(fragment=3)   -2.63724400    1.75809500    0.29129700
  20. F(fragment=3)   -2.66987900   -1.59235800    0.75916300
  21. O(fragment=3)   -4.47650800   -0.10583000   -0.95269300

复制代码
核数，内存和基组/赝势可以根据自己的需求修改。这个gjf可以用GaussView打开可视化，点击Tools->Atom Groups查看分片情况，一目了然。提交计算，即运行
  1. frag_guess_wfn MO24.gjf >MO24.out 2>&1
复制代码
frag_guess_wfn会调用Gaussian完成7次SCF计算，耗时仅2min 30s，生成GKS-EDA输入文件（含所有必须关键词，坐标，基组和赝势信息），并将7套轨道写入MO24.inp文件。提交GKS-EDA任务
  1. $HOME/software/gamess_xeda/rungms MO24.inp 00 64 >MO24.gms 2>&1
复制代码
耗时不到1 min完成计算，7个SCF依次极速收敛。我这里是在本地机器使用XEDA，用户还可以选择把这个inp文件上传到XACS云平台上做计算。整个过程不需要费心基组/赝势的书写，不需要关心SCF收敛问题。

推荐阅读：
关于frag_guess_wfn小程序的更多EDA例子https://jeanwsr.gitlab.io/mokit-doc-mdbook/chap5-3.html
仅使用frag_guess_wfn小程序的解压即用方式https://jeanwsr.gitlab.io/mokit- ... want-frag_guess_wfn
这个页面上有使用frag_guess_wfn小程序自动完成GKS-EDA计算后发表的文章https://jeanwsr.gitlab.io/mokit-doc-mdbook/citing.html

PS：
（1）当前只是一种分片方式，要注意左右两侧的过渡金属分子是否为开壳层的可能，从而形成反铁磁耦合的对称破缺单重态。这种计算在GAMESS里极其痛苦，还无法检验波函数稳定性，而frag_guess_wfn则可以调用Gaussian进行波函数稳定性检验，确保波函数稳定。
（2）当然，会有人推荐你尝试其他办法，但目前来说在GKS-EDA框架下没有一种办法比这层楼提供的办法更简单省时、有效和智能化。
作者
Author: 小翟w    时间: 2024-4-7 11:37
本帖最后由 小翟w 于 2024-4-7 11:40 编辑

老师，我拿frag_guess_wfn小程序算输入文件，SCF不收敛了，老师能帮忙看看吗 (, 下载次数 Times of downloads: 0)
作者
Author: zjxitcc    时间: 2024-4-7 14:54
小翟w 发表于 2024-4-7 11:37
老师，我拿frag_guess_wfn小程序算输入文件，SCF不收敛了，老师能帮忙看看吗

换用Linux g16、多给点内存和核数即可解决，不要用32位Win版G09拖拉机。
作者
Author: 小翟w    时间: 2024-4-7 16:17
zjxitcc 发表于 2024-4-7 14:54
换用Linux g16、多给点内存和核数即可解决，不要用32位Win版G09拖拉机。

好的，谢谢老师，我试试
作者
Author: 小翟w    时间: 2024-5-7 09:58
zjxitcc 发表于 2024-4-7 14:54
换用Linux g16、多给点内存和核数即可解决，不要用32位Win版G09拖拉机。

老师，我换了g16算了，结果数字还是很大，老师有空帮我看看



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