计算化学公社

标题: 利用Multiwfn将Gaussian生成的波函数读入到Gamess中计算，Gamess报波函数读取错误。 [打印本页]

作者
Author: 叫个啥名字呢 时间: 2020-7-2 16:22
标题: 利用Multiwfn将Gaussian生成的波函数读入到Gamess中计算，Gamess报波函数读取错误。
本帖最后由叫个啥名字呢于 2020-7-2 16:30 编辑

我在用gamess计算一个体系一直不收敛，读帖子看到用Gaussian产生的波函数作为初猜，于是用Multiwfn将Gaussian产生的波函数读入到Gamess做计算，运行出错，报错信息以及输入输出信息在附件中。我在Gamess以及Gaussian里面用了赝势，是否两个软件的赝势不匹配。其中两个软件中赝势均用的Stuttgart，其他基组泛函均相同。

作者
Author: zjxitcc 时间: 2020-7-2 16:40
GAMESS读入的时候只支持Cartesian型基函数，你得在gjf里加6D 10F，你现在的gjf默认是5D 7F。另外，最好还要加上int=nobasistransform。

PS：体系很奇怪，还对N用赝势

作者
Author: 叫个啥名字呢 时间: 2020-7-2 16:47

zjxitcc 发表于 2020-7-2 16:40
GAMESS读入的时候只支持Cartesian型基函数，你得在gjf里加6D 10F，你现在的gjf默认是5D 7F。另外，最好还要 ...

多谢大佬！
用Molcas和Gamess比较多，很少用gaussian，是把6D 10F直接写在关键词中就行吗？
PS：特殊用途。

作者
Author: zjxitcc 时间: 2020-7-2 16:48

叫个啥名字呢发表于 2020-7-2 16:47
多谢大佬！
用Molcas和Gamess比较多，很少用gaussian，是把6D 10F直接写在关键词中就行吗？
PS：特殊用 ...

是的。另外，MOKIT （https://gitlab.com/jxzou/mokit）中提供了高斯fchk文件生成OpenMolcas和GAMESS输入文件和轨道的小程序，分别是fch2inporb和fch2inp。

作者
Author: sobereva 时间: 2020-7-2 22:50
都用stuttgart赝势就不会赝势不匹配。但注意Gaussian内置的Stuttgart赝势不够全面，版本也老，因此如果两个程序间互传波函数，stuttgart赝势来源应当统一，比如要么同样来自BSE，要么同样来自Stuttgart赝势官网。

除了Gaussian里用6d 10f外，建议你把GAMESS输入文件里的ispher=1先去掉
我建议你先留意一下Gaussian显示的基函数数目和GAMESS-US显示的是否一致，这个先确保一致了再说收敛的事

另外注意GAMESS-US里与Gausssian对应的B3LYP应当写B3LYPV1R关键词

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