计算化学公社

标题: Gamess 程序中对于cc-pvtz-pp基组的赝势部分书写报错 [打印本页]

作者
Author: xyzabc666 时间: 2020-11-30 19:31
标题: Gamess 程序中对于cc-pvtz-pp基组的赝势部分书写报错

Gamess 程序中对于cc-pvtz-pp基组的赝势部分书写报错（1.inp/1.log），按照论坛中http://bbs.keinsci.com/thread-671-1-1.html以及赝势的函数形式以及在量子化学程序中定义的方式http://sobereva.com/188进行了修正（2.inp/2.log），还是出现错误。赝势部分具体怎么改才能运行正确，恳求各位老师帮帮忙。谢谢！

作者
Author: zjxitcc 时间: 2020-11-30 19:50
本帖最后由 zjxitcc 于 2020-11-30 19:52 编辑

我看BSE基组库上并没有Th元素的cc-pVTZ-PP基组或赝势数据，你是从哪里复制的数据？

作者
Author: xyzabc666 时间: 2020-11-30 19:59

zjxitcc 发表于 2020-11-30 19:50
我看BSE基组库上并没有Th元素的cc-pVTZ-PP基组或赝势数据，你是从哪里复制的数据？

sob老师博文：谈谈赝势基组的选用 http://sobereva.com/373的3.3部分推荐的http://tyr0.chem.wsu.edu/~kipeters/basis.html

作者
Author: zjxitcc 时间: 2020-11-30 20:31

xyzabc666 发表于 2020-11-30 19:59
sob老师博文：谈谈赝势基组的选用 http://sobereva.com/373的3.3部分推荐的http://tyr0.chem.wsu.edu/~ki ...

你这不太对啊，你用了cc-pVTZ-PP的基组数据，赝势部分官网上提醒说
these -PP basis sets should be used ONLY with the ECP60MDF Stuttgart/Koeln pseudopotentials
按照这个提示，赝势你得用ECP60MDF，而你却选了ECP60MWB。

作者
Author: xyzabc666 时间: 2020-11-30 20:51
本帖最后由 xyzabc666 于 2020-11-30 21:00 编辑

zjxitcc 发表于 2020-11-30 20:31
你这不太对啊，你用了cc-pVTZ-PP的基组数据，赝势部分官网上提醒说
these -PP basis sets should be use ...

多谢提醒。找到了它的基组然后就基组库下载一个赝势为60的版本，没注意。按照http://www.tc.uni-koeln.de/PP/clickpse.en.html修正为ECP60MDF了。另外我想问一下，对于Th元素高斯软件用赝势SDD关键词，是这两种的哪一种类型。谢谢。

作者
Author: zjxitcc 时间: 2020-11-30 21:00

xyzabc666 发表于 2020-11-30 20:51
多谢提醒。找到了它的基组然后就基组库下载一个赝势为60的版本，没注意。按照http://www.tc.uni-koeln.de ...

高斯官网（http://gaussian.com/pseudo/）写了对于Th的SDD默认是MWB60。你也可以用gfinput等关键词查看SDD在log文件中的基组/赝势数据。

作者
Author: xyzabc666 时间: 2020-11-30 21:05
好的，谢谢您

作者
Author: zjxitcc 时间: 2020-11-30 21:31

xyzabc666 发表于 2020-11-30 21:05
好的，谢谢您

可能你后面还会上来问，换成ECP60MDF咋还是一样的报错&怎么解决。你还得把g-h改成g-ul，s-h改成s-g，p/d/f类似，因为此时最高角动量变了。然后还得把基组部分的H角动量两行删掉。

作者
Author: xyzabc666 时间: 2020-11-30 22:20
Gamess的ECP60MDF数据版本我还没找到，能推荐个网址吗？还有g-h改成g-ul，s-h改成s-g，p/d/f类似具体是怎么改呀，恳求大佬发个附件。

作者
Author: zjxitcc 时间: 2020-11-30 22:32
本帖最后由 zjxitcc 于 2020-11-30 22:35 编辑

xyzabc666 发表于 2020-11-30 22:20

Gamess的ECP60MDF数据版本我还没找到，能推荐个网址吗？还有g-h改成g-ul，s-h改成s-g，p/d/f类似具体是怎么 ...

就在赝势网站上啊，很显眼
(, 下载次数 Times of downloads: 43)

另外，附上gjf和inp输入文件，如果你想自己手写inp文件，可以参考这个。但其实这个inp文件并不是手写的，而是由高斯fchk文件产生的，用我写的fch2inp小程序（https://gitlab.com/jxzou/mokit）即可。这样做的好处是只需要编辑gjf文件，高斯算完后，用fch2inp小程序从fchk文件产生inp文件，连分子轨道都写在inp文件里了，GAMESS读进去1圈收敛，免了做SCF的时间（GAMESS的SCF远不如高斯稳健，容易不收敛或收敛到鞍点波函数上）。
(, 下载次数 Times of downloads: 21)

作者
Author: xyzabc666 时间: 2020-12-1 09:22
好的，感谢您的帮助！

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