计算化学公社

标题: gamess自定义基组出错求助 [打印本页]

作者
Author: idoit-y 时间: 2023-1-15 19:50
标题: gamess自定义基组出错求助
在gamess中，我用m062x/def2tzvp进行能量分解分析，基组设置参考的是Basis Set Exchange (BSE) 和[原创]Gamess自定义基组与赝势的输入格式——basnam - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com) 和求gamess能量分解输入文件的格式 - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com) ，已经检索了类似的贴子，实在不知道问题出在哪儿了，想求一个大佬解惑。（也祝各位大佬小年快乐！）

作者
Author: paramecium86 时间: 2023-1-16 18:18
输入文件格式问题。最后那个$END 之前没有空格。会提示楼里附件里的报错。而且设置了PP=READ 但是没有给出$ECP也是会报错的

作者
Author: 123qwertybobo 时间: 2023-1-16 21:28
诚如paramecium86所说，写了pp=read，后面要去找金属的小核赝势（MWB28）。还需注意一行不能太长，好像是小于64字符，具体我忘记了。自定义基组我忘了要不要按照第一次出现顺序。大概率还会不收敛，需要加大格点。反正大概是这样，你再试试。

作者
Author: zjxitcc 时间: 2023-1-16 22:40
本帖最后由 zjxitcc 于 2023-1-16 22:43 编辑

正好又是一个展示frag_guess_wfn小程序便捷性的例子。写一个高斯gjf文件，划分好片段，提交（假设叫11.gjf）

frag_guess_wfn 11.gjf >& 11.out &

复制代码

泛函为M062X，H, C, N使用6-311G(d,p)基组，Cr使用def2TZVP基组，O用6-311+G(d)基组，48核耗时4.5 min。算完后产生GKS-EDA输入文件，交给xeda做计算

xeda 11.inp 00 48 >& 11.gms &

复制代码

GAMESS每一步均极速收敛，共耗时3 min。整个过程不用考虑各种格式问题、SCF大概率不收敛、赝势书写等一大堆缺点，不用反复调试、试错，一次成功，计算时间还快得多得多。相关阅读《GKS-EDA计算简介》
(, 下载次数 Times of downloads: 8)

作者
Author: idoit-y 时间: 2023-1-21 10:28

zjxitcc 发表于 2023-1-16 22:40
正好又是一个展示frag_guess_wfn小程序便捷性的例子。写一个高斯gjf文件，划分好片段，提交（假设叫11.gjf ...

好的谢谢老师，我去研究一下。

作者
Author: idoit-y 时间: 2023-1-26 17:59

123qwertybobo 发表于 2023-1-16 21:28
诚如paramecium86所说，写了pp=read，后面要去找金属的小核赝势（MWB28）。还需注意一行不能太长，好像是小 ...

谢谢您，我删掉了pp=read，然后检查了$end前面的所有空格，都存在，也确认了一行不超过64字符，还是报错了，可以麻烦您再看看吗

作者
Author: idoit-y 时间: 2023-1-26 18:00
修改后的输入文件

作者
Author: idoit-y 时间: 2023-1-26 18:02

paramecium86 发表于 2023-1-16 18:18
输入文件格式问题。最后那个$END 之前没有空格。会提示楼里附件里的报错。而且设置了PP=READ 但是没有给 ...

您好，我删除了pp=read，也再次确认了end前面的空格，依旧报错，可以麻烦您再帮我看看吗？

作者
Author: beefly 时间: 2023-1-26 21:02

idoit-y 发表于 2023-1-26 18:02
您好，我删除了pp=read，也再次确认了end前面的空格，依旧报错，可以麻烦您再帮我看看吗？

这个错误信息很明确，H原子的基组错了，某个S收缩函数没有写NGAUSS：
*** ERROR, ILLEGAL BASIS FUNCTION TYPE=1.000000 NGAUSS= 0
ERROR OCCURED READING ATOM NO. 2, INPUT ATOM NO. 2, NAME=H
CHARGE= 1.0 X,Y,Z= -6.82971952 0.67727201 0.34468891

作者
Author: idoit-y 时间: 2023-2-1 19:47

beefly 发表于 2023-1-26 21:02
这个错误信息很明确，H原子的基组错了，某个S收缩函数没有写NGAUSS：
*** ERROR, ILLEGAL BASIS FUNCTI ...

谢谢老师，按照老师的提示我修改后进行了计算，没有再进行报错，但得到的结果没有all basis set和own basis set，只有final energy,感觉是没有进行能量分解，能麻烦老师再帮我看一下吗。

作者
Author: idoit-y 时间: 2023-2-1 19:49
计算得到的输出文件及其输入文件

作者
Author: idoit-y 时间: 2023-2-7 09:59
问题已经解决，非常感谢各位老师的回复。

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