计算化学公社

标题: 请教老师GKS-EDA分析出现NaN错误怎么解决？ [打印本页]

作者
Author: liyaohx 时间: 2022-11-23 11:10
标题: 请教老师GKS-EDA分析出现NaN错误怎么解决？
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）就不附上了

作者
Author: zjxitcc 时间: 2022-11-23 22:05
能把各个元素用的基组说一下么？试试看去掉不必要的弥散函数。

作者
Author: liyaohx 时间: 2022-11-23 23:17

zjxitcc 发表于 2022-11-23 22:05
能把各个元素用的基组说一下么？试试看去掉不必要的弥散函数。

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

作者
Author: zjxitcc 时间: 2022-11-23 23:29

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，我看看？

作者
Author: liyaohx 时间: 2022-11-24 00:26

zjxitcc 发表于 2022-11-23 23:29
感觉是GAMESS的bug。能否把相关文件（gjf，out，5个fch文件，inp, gms）打包压缩发到我邮箱njumath[at]si ...

您好，已发送

作者
Author: zh7729 时间: 2023-9-28 15:14

liyaohx 发表于 2022-11-24 00:26
您好，已发送

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

作者
Author: liyaohx 时间: 2023-12-26 12:22

zh7729 发表于 2023-9-28 15:14
你好，请问你的问题解决了吗？我也遇到了相同的问题，想问一下你最后是怎么解决的

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

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