GVB(1)计算H-F分子输入文件求教

liusy

好的，非常感谢

liusy

zjxitcc 发表于 2019-1-21 21:11
那是我实验室里的前辈，当年他最后并没有改成功，有时候不能work。后来我改成了，一个shell脚本简便易用{ ...

要将gvb的12对变成30对，请问你用的什么脚本？能否告知一下，谢谢！

pyscf

请问GVB相比于普通的DFT计算有啥优势？

zjxitcc

pyscf 发表于 2019-3-23 07:01
请问GVB相比于普通的DFT计算有啥优势？

在能量和性质方面，GVB波函数本身没有动态相关，相对能量等准确度被杂化泛函（甚至是纯泛函）吊打，没有优势（事实上CASSCF大多数时候也是这样）。所以做GVB其实只是万里长征第一步，后面发展类似于CASPT2、NEVPT2等带动态相关的高精度方法才是真正目的。

在波函数方面，由于GVB轨道大多是局域在化学键上的（Boys和PM局域轨道很多时候也这样），定性图像是对的（例如对双自由基）而且比较清晰，讨论分析问题方便，这是屈指可数的用处之一（有别的替代方法，并非特长）。另一个用处是可以作为CASSCF的初始轨道，使之快速收敛（有别的替代方法，并非特长）。

在计算量上比CASSCF小很多，也是仅存的优点之一，便于发展带动态相关的高精度方法。

使用难点只有一个：对过渡金属体系超难收敛（比CAS还难），且操作较为繁琐。4#贴出的paper里我们解决了这个问题。

zjxitcc

liusy 发表于 2019-3-21 15:46
要将gvb的12对变成30对，请问你用的什么脚本？能否告知一下，谢谢！

自己写的shell脚本，暂时不开源，因为说明文档还没写，估计可能今年会开源（加上4楼paper的实现代码）。

看你前面做基本的GVB计算还不是很熟，先熟练做三四对的计算吧···如果对开发方法感兴趣可以先看看W. A. Goddard和Martin Head Gordon等人有关GVB的文章。如果是实用性考虑，那其实没多大用处，很多时候还不如做个CASSCF或DMRG。

如果有过渡金属体系想合作可以私信，不在此一一介绍计算细节，因为后续还有很多问题，时间原因无法逐一回复

pyscf

zjxitcc 发表于 2019-3-23 09:57
在能量和性质方面，GVB波函数本身没有动态相关，相对能量等准确度被杂化泛函（甚至是纯泛函）吊打，没有 ...

谢谢你的回复和解答