计算化学公社

 找回密码 Forget password
 注册 Register

求助gamess计算开壳层单重态总是收敛到闭壳层

查看数: 16257 | 评论数: 15 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-4-1 09:45

正文摘要:

各位老师打了,用gamess计算开壳层单重态总是收敛到闭壳层,是不是写的关键词有问题呀?下边是输入文件   $CONTRL SCFTYP=UHF RUNTYP=ENERGY ICHARG=0 MULT=1      DFTTYP=B3LYPV1R ...

回复 Reply

498746012 发表于 Post on 2021-4-1 16:14:55
zjxitcc 发表于 2021-4-1 15:22
有,但是功能有限,只能基于RHF/RDFT 检验RHF->UHF稳定性,不能基于UHF/UDFT 检验UHF的内部不稳定性(而 ...

老师打扰了,刚才我没有说清楚,我就是想基于RHF/RDFT 检验RHF->UHF稳定性,从而得到开壳层单重态
498746012 发表于 Post on 2021-4-1 16:07:30
paramecium86 发表于 2021-4-1 15:19
不是。
如果想先做DIIS 在做 SOSCF大概写法如下
$scf    dirscf=.true. fdiff=.f. soscf=.f.

谢谢老师,感恩
zjxitcc 发表于 Post on 2021-4-1 15:22:29
498746012 发表于 2021-4-1 15:07
老师您好,学生想请教一下gamess可以做stable=opt的计算吗,有没有相关的关键词呀

有,但是功能有限,只能基于RHF/RDFT 检验RHF->UHF稳定性,不能基于UHF/UDFT 检验UHF的内部不稳定性(而高斯都支持)。看GAMESS手册里的$SCF章节,UHFCHK关键词。

你这一开始算的就是UB3LYP,在GAMESS里不支持检验。

paramecium86 发表于 Post on 2021-4-1 15:19:23
498746012 发表于 2021-4-1 14:41
学生是第一次接触gamess,老师您说的先DIIS再SOSCF, 是算完DIIS,读check文件接着算SOSCF,还是如下一个任 ...

不是。
如果想先做DIIS 在做 SOSCF大概写法如下
$scf    dirscf=.true. fdiff=.f. soscf=.f.
         diis=.t. ethrsh=2.0 shift=.t.
         swdiis=0.001 damp=.t. $end  

其中重要的是 swdiis 这个阈值。 达到这个阈值之后就停止DIIS 开始soscf 而在设置上 还是diis=.t.  soscf=.f.

gamess的scf里面里面选项很多。你最好每个都看看什么意思。不收敛的时候可能哪个就有用。不过我并不知道为什么有时候diis好用 有时候 soscf好用。有时候就是 先diis后soscf 好用。基本上这是我自己拿来做计算时候的经验之谈,自己琢磨着试出来的一些经验。
498746012 发表于 Post on 2021-4-1 15:07:30
zjxitcc 发表于 2021-4-1 14:01
(1)你这个例子 与 帖子一开始说的例子 不是同一个体系,应该讲清楚。乍一眼看上去咋有机体系还能不收敛 ...

老师您好,学生想请教一下gamess可以做stable=opt的计算吗,有没有相关的关键词呀
498746012 发表于 Post on 2021-4-1 14:41:09
paramecium86 发表于 2021-4-1 13:35
SCF不容易收敛 GAMESS的老毛病了。
如果DIIS不行 就换成SOSCF 如果soscf不行 ,就先DIIS再SOSCF。或者干 ...

学生是第一次接触gamess,老师您说的先DIIS再SOSCF, 是算完DIIS,读check文件接着算SOSCF,还是如下一个任务就搞定呀?$SCF DIIS=.t. SOSCF=.t. $END对吗
498746012 发表于 Post on 2021-4-1 14:11:47
zjxitcc 发表于 2021-4-1 14:01
(1)你这个例子 与 帖子一开始说的例子 不是同一个体系,应该讲清楚。乍一眼看上去咋有机体系还能不收敛 ...

谢谢老师,感恩
zjxitcc 发表于 Post on 2021-4-1 14:01:51
498746012 发表于 2021-4-1 12:36
老师我的开壳层单重态不收敛,刚开始能量一直下降到-2060.792,后来能量升高的很厉害,有没有关键词可以 ...

(1)你这个例子 与 帖子一开始说的例子 不是同一个体系,应该讲清楚。乍一眼看上去咋有机体系还能不收敛。。
(2)贴输入文件时最好直接上传文件,或用论坛 插入代码 功能贴出,这样 空格 之类的格式不会乱。你这样直接复制粘贴,有些格式消失了,也不知道你到底写对了没有。
(3)假设你写的输入文件是对的,就按#7楼给你的建议。
498746012 发表于 Post on 2021-4-1 13:54:49
paramecium86 发表于 2021-4-1 13:35
SCF不容易收敛 GAMESS的老毛病了。
如果DIIS不行 就换成SOSCF 如果soscf不行 ,就先DIIS再SOSCF。或者干 ...

谢谢老师,太感谢了
paramecium86 发表于 Post on 2021-4-1 13:35:21
498746012 发表于 2021-4-1 12:36
老师我的开壳层单重态不收敛,刚开始能量一直下降到-2060.792,后来能量升高的很厉害,有没有关键词可以 ...

SCF不容易收敛 GAMESS的老毛病了。
如果DIIS不行 就换成SOSCF 如果soscf不行 ,就先DIIS再SOSCF。或者干脆二者都别用,先用noconv算一个。然后再读
比如可以先试试soscf
$scf    dirscf=.true. fdiff=.f. noconv=.f. soscf=.t. diis=.f. ethrsh=0.002
         shift=.t. damp=.f. $end   

遇到不好收敛的 noconv soscf diis damp 这几个关键词true false各种组合都试试
498746012 发表于 Post on 2021-4-1 10:41:26
zjxitcc 发表于 2021-4-1 10:21
那很简单,GAMESS有提供与高斯guess=mix等价的关键词
$GUESS MIX=.TRUE. $END
用不着IORDER,JORDER。
...

谢谢老师,太感谢了
zjxitcc 发表于 Post on 2021-4-1 10:21:24
498746012 发表于 2021-4-1 10:10
老师您好,我想算开壳层单重态的EDA,无法读取分子轨道,我在调关键词看看能不能算出开壳层单重态

那很简单,GAMESS有提供与高斯guess=mix等价的关键词
$GUESS MIX=.TRUE. $END
用不着IORDER,JORDER。

另外,这只能保证你这个单点能算对,至于EDA里效果如何,还得具体再看。
498746012 发表于 Post on 2021-4-1 10:10:24
zjxitcc 发表于 2021-4-1 09:59
我问个问题哈。。。你都用MOKIT的fch2inp小程序产生GAMESS输入文件了,那算出开壳层单重态就是几秒钟的事, ...

老师您好,我想算开壳层单重态的EDA,无法读取分子轨道,我在调关键词看看能不能算出开壳层单重态

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

GMT+8, 2026-2-25 17:25 , Processed in 0.177261 second(s), 25 queries , Gzip On.

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