求助，SCF并行化方案

hebrewsnabla

zha23 发表于 2022-3-31 19:24
pyscf中有一个参数我不明白，其积分中有一个为cache缓存的参数，但我不知道其作用是什么，以及怎么用
(*in ...

你可以去issue区提问

zjxitcc

zha23 发表于 2022-3-31 17:22
谢谢各位，我所在的课题组是做EDA的，老师觉得GAMESS的效率拖累了EDA的效率，所以想重新写一个自洽场。

GAMESS慢 倒还不是主要问题，你看PSI4也没快到哪里去，但是能免费做SAPT这个特点就是很香；GAMESS的SCF收敛性差，收敛关键词太少，收敛功能太弱、不支持检验波函数稳定性等等缺点限制了GKS-EDA的用途。这些我跟苏老师交流过。

论坛上有很多人做EDA过程中SCF不收敛，为了临时地解决这个问题，我去年在MOKIT里写了frag_guess_wfn小程序(https://gitlab.com/jxzou/mokit)，此小程序可以自动调用Gaussian做HF/DFT计算，然后生成GKS-EDA输入文件（内含坐标、基组和片段及整体的轨道，及GKS-EDA关键词），这样在GAMESS里的GKS-EDA时间就基本可以忽略不计了，一般算3圈就很可能收敛了。

我倒是建议你直接基于PySCF或ORCA，看能否调用他们的接口，加上你自己GKS-EDA做分解部分的代码，得到GKS-EDA分解结果，这样不仅计算速度快，开发程序也快。如果你还能整合进pyscf主分支，或添加进ORCA开源脚本区https://github.com/ORCAQuantumChemistry/CompoundScripts，那就造福大量人了。当然，如果你的抱负是开发一整个量化程序，那从最基本的开始写是可以的。

zha23

zjxitcc 发表于 2022-4-2 00:54
GAMESS慢 倒还不是主要问题，你看PSI4也没快到哪里去，但是能免费做SAPT这个特点就是很香；GAMESS的SCF收 ...

多谢您的建议，我这一部分代码是寒假为了复习之前知识写的，主要我对这方面有一定兴趣，以及除了EDA，XMVB计算的SCF效率也较低，老师也想自己做一套独立的，正好也可以当作明年研究生开题的课题。

喵星大佬

zha23 发表于 2022-4-2 19:23
多谢您的建议，我这一部分代码是寒假为了复习之前知识写的，主要我对这方面有一定兴趣，以及除了EDA，XMV ...

那毕业估计比较危险了

hebrewsnabla

你的上一个帖子不是已经解决了效率问题吗，为什么在issue里面还说efficiency much lower than pyscf

wzkchem5

zha23 发表于 2022-4-2 12:23
多谢您的建议，我这一部分代码是寒假为了复习之前知识写的，主要我对这方面有一定兴趣，以及除了EDA，XMV ...

那我倒觉得不如把XMVB的源码要来，然后在上面改，改善XMVB的效率，这样受益的不止你一个，还有所有XMVB的用户。相比之下，你自己写SCF代码哪怕确实比XMVB快了，也肯定很难有比XMVB更多的用户。
量化软件也是国外卡中国脖子的方向之一，如果每个有志向写量化程序的人都从SCF从头做起，那中国永远也做不出能匹敌国外的软件，等于有100个地基但是没有一座高楼。