计算化学公社

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

[GAMESS-US] 请教老师GKS-EDA分析出现NaN错误怎么解决?

[复制链接 Copy URL]

31

帖子

0

威望

585

eV
积分
616

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
GAMESS的输入文件是frag_guess_wfn生成的,在gaussian下的计算都没有任何问题。但是在GAMESS算到MONMER 2 WITH ALL BASIS时就会出现NaN的诡异问题:                                                         ITER EX      TOTAL ENERGY        E CHANGE  DENSITY CHANGE    DIIS ERROR      VIR. SHIFT       DAMPING        INTEGRALS    SKIPPED          * * *   INITIATING DIIS PROCEDURE   * * *
   1  0                 NaN              NaN   0.000000000           NaN     0.000000000     1.000000000    2636094216012012497669
   2  1                 NaN              NaN   0.000000000           NaN     0.000000000     1.000000000              013334506471

而MONMER 2 WITH OWN BASIS可以正常收敛:
ITER EX      TOTAL ENERGY        E CHANGE  DENSITY CHANGE    DIIS ERROR      VIR. SHIFT       DAMPING        INTEGRALS    SKIPPED
          * * *   INITIATING DIIS PROCEDURE   * * *
   1  0    -1522.2388327210 -1522.2388327210   0.000285946   0.000033622     0.000000000     1.000000000     1175542717  116847328
   2  1    -1522.2388310343     0.0000016866   0.000060958   0.000011627     0.000000000     1.000000000      774311246  141396303
NEXT ENERGY RISES: SWITCHING DIIS OFF FOR  5 ITERS, NEW ETHRSH= 1.05E-05
   3  2    -1522.2388310343     0.0000000000   0.000060958   0.000010297     0.000000000     0.062500000      653380427  147293203
   4  3    -1522.2388310343    -0.0000000000   0.000068830   0.000011627     0.000000000     0.000000000      653380427  147293203
   5  4    -1522.2388310299     0.0000000044   0.000115906   0.000016205     0.000000000     0.000000000      660769987  146942035
   6  5    -1522.2388309995     0.0000000304   0.000205264   0.000029385     0.000000000     0.000001068      682458636  145894549
   7  0    -1522.2388309117     0.0000000878   0.000133993   0.000052813     0.000000000     0.000035740      714458274  144348315
          * * *   INITIATING DIIS PROCEDURE   * * *
   8  1    -1522.2388310431    -0.0000001314   0.000008957   0.000001699     0.000000000     1.000000000      687475045  145664756
   9  2    -1522.2388310426     0.0000000005   0.000005931   0.000002881     0.000000000     1.000000000      546431377  152356585
  10  3    -1522.2388310423     0.0000000003   0.000000597   0.000000310     0.000000000     1.000000000      516969751  153706404


GAMESS的版本是2021-R2,其他相似体系都没问题,可以正常收敛,输出文件太大(150mb)就不附上了


3K-M3-OPT-TPSSH-FromPaper-TZVP.zip

9.66 MB, 阅读权限: 20, 下载次数 Times of downloads: 6

3809

帖子

4

威望

8000

eV
积分
11889

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2022-11-23 22:05:24 | 只看该作者 Only view this author
能把各个元素用的基组说一下么?试试看去掉不必要的弥散函数。
自动做多参考态计算的程序MOKIT

31

帖子

0

威望

585

eV
积分
616

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-11-23 23:17:19 | 只看该作者 Only view this author
zjxitcc 发表于 2022-11-23 22:05
能把各个元素用的基组说一下么?试试看去掉不必要的弥散函数。

Fe和K是def2tzvp,C,H,O,N用的是def2svp

3809

帖子

4

威望

8000

eV
积分
11889

Level 6 (一方通行)

MOKIT开发者

4#
发表于 Post on 2022-11-23 23:29:55 | 只看该作者 Only view this author
liyaohx 发表于 2022-11-23 23:17
Fe和K是def2tzvp,C,H,O,N用的是def2svp

感觉是GAMESS的bug。能否把相关文件(gjf,out,5个fch文件,inp, gms)打包压缩发到我邮箱njumath[at]sina.cn,我看看?
自动做多参考态计算的程序MOKIT

31

帖子

0

威望

585

eV
积分
616

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2022-11-24 00:26:35 | 只看该作者 Only view this author
zjxitcc 发表于 2022-11-23 23:29
感觉是GAMESS的bug。能否把相关文件(gjf,out,5个fch文件,inp, gms)打包压缩发到我邮箱njumath[at]si ...

您好,已发送

49

帖子

0

威望

1136

eV
积分
1185

Level 4 (黑子)

6#
发表于 Post on 2023-9-28 15:14:04 | 只看该作者 Only view this author

你好,请问你的问题解决了吗?我也遇到了相同的问题,想问一下你最后是怎么解决的

31

帖子

0

威望

585

eV
积分
616

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2023-12-26 12:22:43 | 只看该作者 Only view this author
zh7729 发表于 2023-9-28 15:14
你好,请问你的问题解决了吗?我也遇到了相同的问题,想问一下你最后是怎么解决的

目前还没有解决,感觉是BUG

本版积分规则 Credits rule

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

GMT+8, 2024-11-25 17:59 , Processed in 0.186820 second(s), 23 queries , Gzip On.

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