计算化学公社

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

[GAMESS-US] gamess-us 能量分解的SCF能量从负值突然变到正值

[复制链接 Copy URL]

581

帖子

0

威望

1700

eV
积分
2281

Level 5 (御坂)

各位老师好,想请教一下为什么gamess的能量分解SUPERMOLECULE SCF能量从负值变到了正值,最终导致SCF不收敛,如何解决这种不收敛问题呀?输入文件如下
$CONTRL SCFTYP=UHF RUNTYP=EDA ICHARG=0 MULT=1 NOSYM=1 ICUT=11 PP=READ
  DFTTYP=B3LYPV1R MAXIT=200 ISPHER=1 $END
$SYSTEM MWORDS=1500 MEMDDI=9000  $END
$DFT NRAD=99 NLEB=590 $END
$SCF  dirscf=.t. diis=.t. damp=.t. $END
$GUESS MIX=.TRUE. $END
$lmoeda matom(1)=53,19 mcharg(1)=0,0 mmult(1)=1,1 $end
$DATA
GAMESS inp format file produced by MOKIT, nbf=767


问题11.JPG (572.5 KB, 下载次数 Times of downloads: 25)

问题11.JPG

334

帖子

1

威望

6155

eV
积分
6509

Level 6 (一方通行)

2#
发表于 Post on 2021-4-4 12:29:30 | 只看该作者 Only view this author
本帖最后由 paramecium86 于 2021-4-4 15:05 编辑

GAMESS的SCF不稳 抽风什么的属于常规现象 。DIIS不行就换SOSCF试一试比如写成

$scf    dirscf=.true. fdiff=.f. noconv=.f. soscf=.t. diis=.f. ethrsh=0.002
         shift=.t. damp=.f. $end   

或者先DIIS再SOSCF
$scf    dirscf=.true. fdiff=.f. soscf=.f.
         diis=.t. ethrsh=2.0 shift=.t.
         swdiis=0.001 damp=.t. $end  

各种控制SCF收敛的方法都试一试先。
嗯。。。我才发现前两天的帖子回复的就是你呀。。。 那看来你这个以上方法是解决不了你这个问题

581

帖子

0

威望

1700

eV
积分
2281

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2021-4-5 10:36:12 | 只看该作者 Only view this author
paramecium86 发表于 2021-4-4 12:29
GAMESS的SCF不稳 抽风什么的属于常规现象 。DIIS不行就换SOSCF试一试比如写成

$scf    dirscf=.true. fd ...

谢谢老师,感恩

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 02:23 , Processed in 0.325532 second(s), 23 queries , Gzip On.

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