GAMESS-US计算含过渡金属体系单点读入Gaussian轨道信息计算异常

xrpeng

原本计算的体系较大，难以收敛，使用小体系（FeO）测试（未进行优化和频率计算）
使用Gaussian16在UDFT水平下计算单点，使用Multiwfn将*.fchk转化为GAMESS-US输入文件，只取轨道部分并只保留了占据轨道，读入GANESS-US的*.inp文件中，计算发现第一步能量与Gaussian收敛后的能量差距较大，后续迭代过程也很奇怪
考虑不同程序对赝势的处理可能不同，使用O2对比分析，没有上述情况，所以问题很有可能在赝势这里，但两个程序Fe原子的赝势基组均为读入，保持一致，所以不知道是哪里出了问题。测试文件在压缩包中。

zjxitcc

建议使用MOKIT的fch2inp小程序将fch文件转化为inp文件，无需你写基组名、手动添加赝势数据、添加PP=READ等关键词，且保证一切数据严格对应，GAMESS的SCF立即收敛

1.                                                                                 NONZERO     BLOCKS
   
2. ITER EX      TOTAL ENERGY        E CHANGE  DENSITY CHANGE    DIIS ERROR      INTEGRALS    SKIPPED
   
3.           * * *   INITIATING DIIS PROCEDURE   * * *
   
4.    1  0     -198.3408673384  -198.3408673384   0.000077581   0.000049649          76676          0
   
5.    2  1     -198.3408673380     0.0000000004   0.000006555   0.000002949          76668          4
   
6.    3  2     -198.3408673380     0.0000000000   0.000005553   0.000005048          76657          8
   
7. 
   
8.           -----------------
   
9.           DENSITY CONVERGED
   
10.           -----------------

复制代码

使用方式：运行fch2inp FeO-uGAU.fchk。
结果见附件 FeO-uGAU.inp (50.94 KB, 下载次数 Times of downloads: 4)

2023-9-17 21:22 上传 Uploaded

点击下载
Click to download

FeO-uGAU.gms (88.88 KB, 下载次数 Times of downloads: 1)

2023-9-17 21:22 上传 Uploaded

点击下载
Click to download

MOKIT可联网安装、下载Linux预编译版直接使用、下载Win预编译版直接使用、自己编译各种方式，见
https://gitlab.com/jxzou/mokit/-/blob/master/README_zh.md
https://jeanwsr.gitlab.io/mokit-doc-mdbook/chap2_install.html

相关中文教程《利用MOKIT从Gaussian向其他量化程序传轨道》。顺带一提，虽然你在gjf文件里写了6D 10F；即使你写5D 7F，fch2inp小程序也支持传轨道，让GAMESS用球谐型基函数做计算。

xrpeng

zjxitcc 发表于 2023-9-17 21:24
建议使用MOKIT的fch2inp小程序将fch文件转化为inp文件，无需你写基组名、手动添加赝势数据、添加PP=READ等 ...

谢谢！我这就试试！

xrpeng

zjxitcc 发表于 2023-9-17 21:24
建议使用MOKIT的fch2inp小程序将fch文件转化为inp文件，无需你写基组名、手动添加赝势数据、添加PP=READ等 ...

老师，我还是想知道我的输入文件哪里出了问题，然后我发现是您的输入文件中Fe的基组和我写的的不一样，请问这是因为Gaussian和GAMESS处理过渡金属时不一样吗，为什么我直接从基组网站中获取的基组不能用呢

zjxitcc

xrpeng 发表于 2023-9-18 15:32
老师，我还是想知道我的输入文件哪里出了问题，然后我发现是您的输入文件中Fe的基组和我写的的不一样，请 ...

我不知道你的文件哪里有问题，因为我一直使用fch2inp，不手写gamess输入文件，从不会碰到这种问题。Gaussian和GAMESS处理过渡金属一样。

关于基组和赝势：我在2L上传的文件里的基组和赝势不是我手写的，是fch2inp小程序生成的。传轨道时本来就应该连同基组和赝势数据一起传过去，不能（至少不建议）使用其他程序内置基组/赝势，也不需要手写。例如叫cc-pVTZ名称的基组，在不同量化程序里具体基组数据经常不一样。只有同时传基组和赝势数据，才能保证两个程序的计算严格对应（当然，还有一些其他因素也要对应）。

xrpeng

zjxitcc 发表于 2023-9-18 15:45
我不知道你的文件哪里有问题，因为我一直使用fch2inp，不手写gamess输入文件，从不会碰到这种问题。Gaussia ...

好的，谢谢！