计算化学公社

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

[GAMESS-US] gamess不收敛

[复制链接 Copy URL]

442

帖子

0

威望

1500

eV
积分
1942

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 413 于 2022-11-1 10:34 编辑

请问下面这个错误咋解决呢?我搜了gamess forum上,那个有个建议,但是没解决呢
我还是了guess=hcore,这个错误是避开了,但是200圈的时候scf没收敛

“gamess forum的建议http://classic.chem.msu.su/cgi-b ... HW-8870-1309+00.htm
As a workaround, try to use: $SCF DIRSCF=.TRUE. FDIFF=.FALSE. NCONV=6  EXTRAP=.TRUE. DAMP=.TRUE. SHIFT=.TRUE. RSTRCT=.FALSE.   DIIS=.TRUE. SOSCF=.FALSE. $END

                    * * * * * WARNING * * * * *
POSSIBLE LINEAR DEPENDENCE PROBLEMS DETECTED, INPUT QMTTOL= 1.0E-06,
WITH     7 OVERLAP EIGENVALUES LYING BETWEEN 1.0E-05 AND QMTTOL.
IF YOU EXPERIENCE SCF CONVERGENCE DIFFICULTIES WITH THE DENSITY ALMOST BUT
NOT QUITE CONVERGED, YOU MAY NEED TO LOOSEN -QMTTOL- BY A FACTOR 5 TO 10.
**** ERROR, ORBITAL   1 IS MISSING FROM REORDERING INSTRUCTIONS
    0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
    0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
    0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
    0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0
SOMETIMES THE ABOVE ERROR IS CASED BY AN INABILITY TO
PROJECT THE HUCKEL GUESS ONTO THE CURRENT BASIS SET,
E.G. WHEN THE 1S HUCKEL ORBITAL OF A H ATOM CANNOT
BE PROJECTED ONTO A PURELY P BASIS SET BECAUSE YOU
ARE TRYING TO DO EXCITED STATE H ATOMS.  IN SUCH A
CASE, PLEASE USE GUESS=HCORE OR GUESS=MOREAD INSTEAD.
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT 10:23:14 01-NOV-2022

4289

帖子

4

威望

9541

eV
积分
13910

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2022-11-1 11:22:17 | 只看该作者 Only view this author
本帖最后由 zjxitcc 于 2022-11-1 11:24 编辑

如果你是做常规HF/DFT计算,推荐使用MOKIT中的fch2inp小程序从Gaussian向GAMESS传轨道;
如果你是做EDA计算,推荐使用GKS-EDA方法,并且使用frag_guess_wfn小程序产生带片段和整体轨道信息的.inp文件。建议阅读《GKS-EDA计算简介
不建议尝试所列举的英文建议。GAMESS的SCF收敛性远不如Gaussian,所以不建议在GAMESS中换关键词试来试去(关键词加多了 自己也没明白是啥意思),而且即使好不容易收敛了 也没法检验波函数稳定性。
使用上述两种技巧,论坛上几乎所有GAMESS单点SCF不收敛、xxx-EDA过程中SCF不收敛的帖子都可得到解决,彻底终结此类问题。

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

442

帖子

0

威望

1500

eV
积分
1942

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2022-11-1 12:09:24 | 只看该作者 Only view this author
zjxitcc 发表于 2022-11-1 11:22
如果你是做常规HF/DFT计算,推荐使用MOKIT中的fch2inp小程序从Gaussian向GAMESS传轨道;
如果你是做EDA计 ...

谢谢提醒哈
其实我就是用的fch2inp传的轨道,后来自己手动改了命令行,把MOREAD抹去了,结果搞的不收敛
现在已经重新提交上去了,谢谢谢谢。有了结果我会及时反馈

442

帖子

0

威望

1500

eV
积分
1942

Level 5 (御坂)

4#
 楼主 Author| 发表于 Post on 2022-11-1 21:21:26 | 只看该作者 Only view this author
zjxitcc 发表于 2022-11-1 11:22
如果你是做常规HF/DFT计算,推荐使用MOKIT中的fch2inp小程序从Gaussian向GAMESS传轨道;
如果你是做EDA计 ...

收敛了,加上了moread这个关键词

本版积分规则 Credits rule

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

GMT+8, 2026-2-21 12:46 , Processed in 0.171353 second(s), 20 queries , Gzip On.

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