计算化学公社

标题: GAMESS-US计算含过渡金属体系单点读入Gaussian轨道信息计算异常 [打印本页]
作者
Author: xrpeng    时间: 2023-9-17 20:25
标题: GAMESS-US计算含过渡金属体系单点读入Gaussian轨道信息计算异常
原本计算的体系较大，难以收敛，使用小体系（FeO）测试（未进行优化和频率计算）
使用Gaussian16在UDFT水平下计算单点，使用Multiwfn将*.fchk转化为GAMESS-US输入文件，只取轨道部分并只保留了占据轨道，读入GANESS-US的*.inp文件中，计算发现第一步能量与Gaussian收敛后的能量差距较大，后续迭代过程也很奇怪
考虑不同程序对赝势的处理可能不同，使用O2对比分析，没有上述情况，所以问题很有可能在赝势这里，但两个程序Fe原子的赝势基组均为读入，保持一致，所以不知道是哪里出了问题。测试文件在压缩包中。
作者
Author: zjxitcc    时间: 2023-9-17 21:24
本帖最后由 zjxitcc 于 2023-9-17 21:34 编辑

建议使用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

  8.           -----------------
  9.           DENSITY CONVERGED
  10.           -----------------
复制代码

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

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用球谐型基函数做计算。
作者
Author: xrpeng    时间: 2023-9-18 09:41
zjxitcc 发表于 2023-9-17 21:24
建议使用MOKIT的fch2inp小程序将fch文件转化为inp文件，无需你写基组名、手动添加赝势数据、添加PP=READ等 ...

谢谢！我这就试试！
作者
Author: xrpeng    时间: 2023-9-18 15:32
zjxitcc 发表于 2023-9-17 21:24
建议使用MOKIT的fch2inp小程序将fch文件转化为inp文件，无需你写基组名、手动添加赝势数据、添加PP=READ等 ...

老师，我还是想知道我的输入文件哪里出了问题，然后我发现是您的输入文件中Fe的基组和我写的的不一样，请问这是因为Gaussian和GAMESS处理过渡金属时不一样吗，为什么我直接从基组网站中获取的基组不能用呢
作者
Author: zjxitcc    时间: 2023-9-18 15:45
本帖最后由 zjxitcc 于 2023-9-18 15:52 编辑
xrpeng 发表于 2023-9-18 15:32
老师，我还是想知道我的输入文件哪里出了问题，然后我发现是您的输入文件中Fe的基组和我写的的不一样，请 ...
我不知道你的文件哪里有问题，因为我一直使用fch2inp，不手写gamess输入文件，从不会碰到这种问题。Gaussian和GAMESS处理过渡金属一样。

关于基组和赝势：我在2L上传的文件里的基组和赝势不是我手写的，是fch2inp小程序生成的。传轨道时本来就应该连同基组和赝势数据一起传过去，不能（至少不建议）使用其他程序内置基组/赝势，也不需要手写。例如叫cc-pVTZ名称的基组，在不同量化程序里具体基组数据经常不一样。只有同时传基组和赝势数据，才能保证两个程序的计算严格对应（当然，还有一些其他因素也要对应）。
作者
Author: xrpeng    时间: 2023-9-18 16:01
zjxitcc 发表于 2023-9-18 15:45
我不知道你的文件哪里有问题，因为我一直使用fch2inp，不手写gamess输入文件，从不会碰到这种问题。Gaussia ...

好的，谢谢！



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3