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

zjxitcc

我试图使用G09的GVB(1)计算氟化氢分子的单点能，在调换轨道那一步出错了（No more virtuals for natural orbitals! Error termination），不是很理解如何构造GVB的输入文件，感觉可能是对称性没指定，但不知道如何指定。请版上各位大大指导一下，gjf和log文件见附件。PS：我alter了占据分子轨道3和5的顺序，是因为3是F的2pz与H的1s混合而成的成键轨道，6是对应的反键轨道，想着应该配成一对，才是GVB(1)；下行写的2是看了高斯的手册，GVB需要写2，表示“GVB对中的每一个都有2个自然轨道”，不过个人觉得这行很鸡肋，都叫GVB了，当然是2，为什么还要写2？无奈GVB计算例子太少，没学习出正确的输入文件格式。

zjxitcc

已解决。关于GVB计算网上没有多少解疑和实例，最后自己再看了几遍高斯说明书GVB计算部分，慢慢摸索出来了。先说结果：轨道我没调换错，不指定对称性也行（即GVB选项里不需要lowsym关键词，至少我算过的几个简单算例都不需要）。trick在于GVB-PP（完美配对）计算的话，gjf文件末尾需要一行“ 2”输入；GVB(n)有n对，就要写n个2，且不能顶行写，需空1格。附件我上传了H-F和CH4两个格式正确的gjf文件。

另外，还发现GAMESS和Gaussian09的波函数有些时候会收敛到不同的局域最小点上，导致能量有明显差别，这一点与CASSCF计算很像。打算以后有空写一个帖子讨论GVB的计算格式和收敛。还学会了以G09的结果为初始提供给GAMESS使用，这样两个软件就能收敛到同一个较低能量的波函数上了。

liusy

请问为什么是轨道3和5调换？

zjxitcc

liusy 发表于 2019-1-21 12:15
请问为什么是轨道3和5调换？

我更新了2#中附的两个分子的输入文件，更为简洁和严谨。

在H-F分子这个例子里，3号轨道是H-F成键，5号是F的2p孤电子对，6号像H-F反键（还不确定，可以通过优化得到更好的结果）。我们做GVB要求初始轨道尽量是一个成键配一个反键（这像极了做CAS(2,2)，如果你熟悉CAS计算的话），所以我们要交换3与5，使得3与6处在配对的位置。
注意：高斯里活性轨道是从HOMO的标号往下数，LUMO的标号往上数，所以要把活性轨道挪往HOMO-LUMO线附近。

在CH4这个例子里，我故意把一根C-H键缩短，目的是使其轨道更容易被辨认出。同理，这根C-H键的成键轨道在3号，反键轨道在9号，要把它们挪往HOMO-LUMO线附近，所以alter了3与5，6与9.

这是早期我刚开始学习GVB计算时的做法，这一套做法现在来看已经过时了，只能对简单小分子玩一玩。一般地，还要结合局域化方法、键解离时考虑RHF不稳定等等，如果对GVB感兴趣可以看我们前段时间发的文章DOI: 10.1021/acs.jctc.8b00854
《Automatic Construction of the Initial Orbitals for Efficient Generalized Valence Bond Calculations of Large Systems》

liusy

好的，非常感谢，我稍后看下文章，请问你是否现在用的gamess算gvb？

zjxitcc

liusy 发表于 2019-1-21 20:02
好的，非常感谢，我稍后看下文章，请问你是否现在用的gamess算gvb？

是的，GAMESS的SOSCF收敛算法，结合我们产生的良好的初始轨道，收敛不是问题。高斯做过渡金属的GVB计算简直是灾难。另外，GAMESS最大只能做12对，想做更大的需要改源代码重新编译。

liusy

好的，谢谢指点，目测论坛有人谈过将GVB 12对变成30对，我现在做的主要是过渡金属的了，应该要用到的超过12对

zjxitcc

liusy 发表于 2019-1-21 20:22
好的，谢谢指点，目测论坛有人谈过将GVB 12对变成30对，我现在做的主要是过渡金属的了，应该要用到的超过12 ...

那是我实验室里的前辈，当年他最后并没有改成功，有时候不能work。后来我改成了，一个shell脚本简便易用

liusy

哇哦，太厉害了，你们一般用哪个版本的gamess

zjxitcc

liusy 发表于 2019-1-21 21:39
哇哦，太厉害了，你们一般用哪个版本的gamess

一般用2017版本的，因为我开始做的时候出了2017。2016和2018版我想应该也没啥问题。更早期的话，有些变量名不一样，因此修改之处会略微有些不同。

liusy

安装GAMMES的编译环境是需要编译器和数学库：ifort 12.1.0+MKL,请问你那边还有安装包没？如果有的话，可以发我一份吗，谢谢！

zjxitcc

liusy 发表于 2019-1-23 19:04
安装GAMMES的编译环境是需要编译器和数学库：ifort 12.1.0+MKL,请问你那边还有安装包没？如果有的话，可以 ...

Intel编译器（含MKL）学术版免费，官网可下。此论坛亦有，可自行搜索解决。

liusy

got，感谢

liusy

zjxitcc 发表于 2019-1-23 19:41
Intel编译器（含MKL）学术版免费，官网可下。此论坛亦有，可自行搜索解决。

请问你可以发一个用GAMESS计算gvb的input文件给我参考下吗？

zjxitcc

liusy 发表于 2019-2-24 20:00
请问你可以发一个用GAMESS计算gvb的input文件给我参考下吗？

GAMESS安装包下本来就有例子，在gamess/tests/standard目录下用grep 'SCFTYP=GVB' -in *.inp命令搜索，就可以看到哪些文件是GVB任务。