计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2047|回复 Reply: 11
打印 Print 上一主题 Last thread 下一主题 Next thread

[GAMESS-US] gamess自定义基组出错求助

[复制链接 Copy URL]

12

帖子

0

威望

521

eV
积分
533

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
在gamess中,我用m062x/def2tzvp进行能量分解分析,基组设置参考的是Basis Set Exchange (BSE)[原创]Gamess自定义基组与赝势的输入格式——basnam - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com)求gamess能量分解输入文件的格式 - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com) ,已经检索了类似的贴子,实在不知道问题出在哪儿了,想求一个大佬解惑。(也祝各位大佬小年快乐!)

11.inp

9.62 KB, 下载次数 Times of downloads: 16

11.gms

21.47 KB, 下载次数 Times of downloads: 7

12

帖子

0

威望

521

eV
积分
533

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2023-2-7 09:59:10 | 只看该作者 Only view this author
问题已经解决,非常感谢各位老师的回复。

12

帖子

0

威望

521

eV
积分
533

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2023-2-1 19:49:01 | 只看该作者 Only view this author
计算得到的输出文件及其输入文件

11.inp

9.62 KB, 下载次数 Times of downloads: 8

11.gms

106.34 KB, 下载次数 Times of downloads: 1

12

帖子

0

威望

521

eV
积分
533

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2023-2-1 19:47:25 | 只看该作者 Only view this author
beefly 发表于 2023-1-26 21:02
这个错误信息很明确,H原子的基组错了,某个S收缩函数没有写NGAUSS:
*** ERROR, ILLEGAL BASIS FUNCTI ...

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

744

帖子

21

威望

5351

eV
积分
6515

Level 6 (一方通行)

9#
发表于 Post on 2023-1-26 21:02:14 | 只看该作者 Only view this author
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

12

帖子

0

威望

521

eV
积分
533

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2023-1-26 18:02:52 | 只看该作者 Only view this author
paramecium86 发表于 2023-1-16 18:18
输入文件格式问题。最后那个$END 之前没有空格。会提示楼里附件里的报错 。 而且设置了PP=READ 但是没有给 ...

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

12

帖子

0

威望

521

eV
积分
533

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2023-1-26 18:00:30 | 只看该作者 Only view this author
修改后的输入文件

11.gms

26.22 KB, 下载次数 Times of downloads: 4

11.inp

9.62 KB, 下载次数 Times of downloads: 5

12

帖子

0

威望

521

eV
积分
533

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2023-1-26 17:59:33 | 只看该作者 Only view this author
123qwertybobo 发表于 2023-1-16 21:28
诚如paramecium86所说,写了pp=read,后面要去找金属的小核赝势(MWB28)。还需注意一行不能太长,好像是小 ...

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

12

帖子

0

威望

521

eV
积分
533

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2023-1-21 10:28:29 | 只看该作者 Only view this author
zjxitcc 发表于 2023-1-16 22:40
正好又是一个展示frag_guess_wfn小程序便捷性的例子。写一个高斯gjf文件,划分好片段,提交(假设叫11.gjf ...

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

4289

帖子

4

威望

9541

eV
积分
13910

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2023-1-16 22:40:03 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2023-1-16 22:43 编辑

正好又是一个展示frag_guess_wfn小程序便捷性的例子。写一个高斯gjf文件,划分好片段,提交(假设叫11.gjf)
  1. 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做计算
  1. xeda 11.inp 00 48 >& 11.gms &
复制代码

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

评分 Rate

参与人数
Participants 2
eV +8 收起 理由
Reason
hebrewsnabla + 3 GJ!
123qwertybobo + 5 其实简单的可以写个小脚本,不过mokit确实g.

查看全部评分 View all ratings

自动做多参考态计算的程序MOKIT

26

帖子

0

威望

1103

eV
积分
1129

Level 4 (黑子)

3#
发表于 Post on 2023-1-16 21:28:18 | 只看该作者 Only view this author
诚如paramecium86所说,写了pp=read,后面要去找金属的小核赝势(MWB28)。还需注意一行不能太长,好像是小于64字符,具体我忘记了。自定义基组我忘了要不要按照第一次出现顺序。大概率还会不收敛,需要加大格点。反正大概是这样,你再试试。

334

帖子

1

威望

6153

eV
积分
6507

Level 6 (一方通行)

2#
发表于 Post on 2023-1-16 18:18:51 | 只看该作者 Only view this author
输入文件格式问题。最后那个$END 之前没有空格。会提示楼里附件里的报错 。 而且设置了PP=READ 但是没有给出$ECP也是会报错的

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 02:02 , Processed in 0.188399 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list