zjxitcc 发表于 2021-4-1 15:22 老师打扰了,刚才我没有说清楚,我就是想基于RHF/RDFT 检验RHF->UHF稳定性,从而得到开壳层单重态 |
paramecium86 发表于 2021-4-1 15:19 谢谢老师,感恩 |
498746012 发表于 2021-4-1 15:07 有,但是功能有限,只能基于RHF/RDFT 检验RHF->UHF稳定性,不能基于UHF/UDFT 检验UHF的内部不稳定性(而高斯都支持)。看GAMESS手册里的$SCF章节,UHFCHK关键词。 你这一开始算的就是UB3LYP,在GAMESS里不支持检验。 |
498746012 发表于 2021-4-1 14:41 不是。 如果想先做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 好用。基本上这是我自己拿来做计算时候的经验之谈,自己琢磨着试出来的一些经验。 |
zjxitcc 发表于 2021-4-1 14:01 老师您好,学生想请教一下gamess可以做stable=opt的计算吗,有没有相关的关键词呀 |
paramecium86 发表于 2021-4-1 13:35 学生是第一次接触gamess,老师您说的先DIIS再SOSCF, 是算完DIIS,读check文件接着算SOSCF,还是如下一个任务就搞定呀?$SCF DIIS=.t. SOSCF=.t. $END对吗 |
zjxitcc 发表于 2021-4-1 14:01 谢谢老师,感恩 |
498746012 发表于 2021-4-1 12:36 (1)你这个例子 与 帖子一开始说的例子 不是同一个体系,应该讲清楚。乍一眼看上去咋有机体系还能不收敛。。 (2)贴输入文件时最好直接上传文件,或用论坛 插入代码 功能贴出,这样 空格 之类的格式不会乱。你这样直接复制粘贴,有些格式消失了,也不知道你到底写对了没有。 (3)假设你写的输入文件是对的,就按#7楼给你的建议。 |
paramecium86 发表于 2021-4-1 13:35 谢谢老师,太感谢了 |
498746012 发表于 2021-4-1 12:36 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各种组合都试试 |
zjxitcc 发表于 2021-4-1 10:21 谢谢老师,太感谢了 |
498746012 发表于 2021-4-1 10:10 那很简单,GAMESS有提供与高斯guess=mix等价的关键词 $GUESS MIX=.TRUE. $END 用不着IORDER,JORDER。 另外,这只能保证你这个单点能算对,至于EDA里效果如何,还得具体再看。 |
zjxitcc 发表于 2021-4-1 09:59 老师您好,我想算开壳层单重态的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.