计算化学公社

标题: GAMESS-US计算过渡金属不收敛 [打印本页]

作者
Author: jiayuhan 时间: 2020-3-14 18:34
标题: GAMESS-US计算过渡金属不收敛
各位好，请问在使用GAMESS-US计算能量分解时遇到不收敛的问题怎么处理。
我的体系是Rh团簇与H2O的配合物，MP2， SDD for Rh, 6-311G** for O,H 方法。
看总能量变化感觉是能量震荡导致的不收敛，请教各位大神们这种情况该如何修正使其收敛
谢谢

作者
Author: sobereva 时间: 2020-3-14 18:44
尝试DIIS=.TRUE. ETHRSH=2.0 SWDIIS=0.005
写$DFT SWITCH=1.0 $END没意义，当前又不是DFT计算
也可尝试用DAMP设置阻尼、SHIFT设置能级移动

也可以用Gaussian等其它程序计算得到fch/molden文件，载入Multiwfn，用主功能100的子功能2导出GAMESS-US输入文件时选择把当前波函数也写入文件作为初猜，并且挪到你的输入文件里读取当初猜用

另：过渡金属团簇用MP2描述非常垃圾，耗时还高，强烈建议改用DFT。

作者
Author: jiayuhan 时间: 2020-3-14 18:55

sobereva 发表于 2020-3-14 18:44
尝试DIIS=.TRUE. ETHRSH=2.0 SWDIIS=0.005
写$DFT SWITCH=1.0 $END没意义，当前又不是DFT计算
也可尝试用 ...

感谢sob老师
请问用DFT做GAMESS的输入文件该怎么写

作者
Author: sobereva 时间: 2020-3-14 19:37

jiayuhan 发表于 2020-3-14 18:55
感谢sob老师
请问用DFT做GAMESS的输入文件该怎么写

论坛里以往GAMESS分类里的帖子就有现成的例子

作者
Author: jiayuhan 时间: 2020-3-14 23:00

sobereva 发表于 2020-3-14 19:37
论坛里以往GAMESS分类里的帖子就有现成的例子

感谢Sob老师
老师还有一个问题，就是GAMESS-US可以续算，或是降低收敛标准吗
如果体系在200步的时候马上就要收敛但没收敛该如何解决呢

作者
Author: zjxitcc 时间: 2020-3-14 23:14

jiayuhan 发表于 2020-3-14 23:00
感谢Sob老师
老师还有一个问题，就是GAMESS-US可以续算，或是降低收敛标准吗
如果体系在200步的时候马 ...

可以续算。第200步的结果应该存在.dat文件里，你可以把它改成.inp文件（注意修改关键词）再算。

作者
Author: jiayuhan 时间: 2020-3-14 23:25

zjxitcc 发表于 2020-3-14 23:14
可以续算。第200步的结果应该存在.dat文件里，你可以把它改成.inp文件（注意修改关键词）再算。

非常感谢！！！我试试去

作者
Author: jiayuhan 时间: 2020-3-16 12:05

sobereva 发表于 2020-3-14 19:37
论坛里以往GAMESS分类里的帖子就有现成的例子

Sob老师你好，在LMOEDA的结果中并没有Eorb 一项。得到的Epolar 可以等效于其他EDA程序得到的Eorb吗

作者
Author: sobereva 时间: 2020-3-17 13:18

jiayuhan 发表于 2020-3-16 12:05
Sob老师你好，在LMOEDA的结果中并没有Eorb 一项。得到的Epolar 可以等效于其他EDA程序得到的Eorb吗

我不知道你说的其它的是哪个
字面上有一定对应关系，但也有的程序还单独拆分出CT项

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