计算化学公社

标题: 用MOKIT程序生成的输入文件 [打印本页]

作者
Author: wangxue123 时间: 2024-5-14 22:08
标题: 用MOKIT程序生成的输入文件
本帖最后由 wangxue123 于 2024-5-15 16:49 编辑

老师您好，
我用molpro计算casscf，希望将molpro计算得到的轨道转移到高斯中。按照手册，需要提供一个.fch文件，所以我用高斯做了RHF计算得到了一个.fch文件，molpro计算得到了.xml文件，然后执行xml2fch a.xml a.fch，有如下报错。mokit是用conda安装的.

(, 下载次数 Times of downloads: 22)

作者
Author: zjxitcc 时间: 2024-5-14 22:15
本帖最后由 zjxitcc 于 2024-5-14 22:36 编辑

需先用fch2com小程序从fch文件生成com文件，用这个com文件做计算（后缀名你可以改成inp，这个文件里关闭了点群对称性，因为xml2fch不支持点群对称性），获得同名xml文件，然后使用xml2fch传轨道回fch文件。一开始的gjf文件里记得写上nosymm int=nobasistransform。

你没有用fch2com生成com文件，而是自己写了一个com文件，并且用了Molpro内置基组；后面Gaussian的fch文件还是独立的（与Molpro计算不相关），这目前不被xml2fch所支持。

作者
Author: wangxue123 时间: 2024-5-15 16:46
本帖最后由 wangxue123 于 2024-5-15 16:48 编辑

zjxitcc 发表于 2024-5-14 22:15
需先用fch2com小程序从fch文件生成com文件，用这个com文件做计算（后缀名你可以改成inp，这个文件里关闭了 ...

老师您好，我用fch2com生成了一个.com文件，然后执行计算，有如下报错。请问高斯提供的.fch文件需要和molpro执行的计算需要一样吗？现在我想将molpro做casscf计算得到的轨道传给高斯，那提供的.fch文件需要用的是高斯做casscf计算所提供的吗，还是说我用高斯做RHF计算得到的.fch文件也行？
(, 下载次数 Times of downloads: 24)

作者
Author: zjxitcc 时间: 2024-5-15 16:52

wangxue123 发表于 2024-5-15 16:46
老师您好，我用fch2com生成了一个.com文件，然后执行计算，有如下报错。请问高斯提供的.fch文件需要和mol ...

用高斯做RHF/ROHF计算获得h2o.fch文件，然后fch2com产生h2o.com文件，用Molpro做目标计算，用xml2fch传轨道回h2o.fch文件
xml2fch h2o.xml h2o.fch（你给错参数了，这条命令与h2o.com无关）

作者
Author: wangxue123 时间: 2024-5-31 17:47
本帖最后由 wangxue123 于 2024-5-31 17:48 编辑

zjxitcc 发表于 2024-5-15 16:52
用高斯做RHF/ROHF计算获得h2o.fch文件，然后fch2com产生h2o.com文件，用Molpro做目标计算，用xml2fch传轨 ...

感谢老师的回答，我用.com文件进行sa-casscf计算，显示
error
cannot find USERDEF-JFIT-FIT in basis catalogue
the problem occurs in binput

(, 下载次数 Times of downloads: 21) (, 下载次数 Times of downloads: 24)

作者
Author: zjxitcc 时间: 2024-5-31 18:27
本帖最后由 zjxitcc 于 2024-5-31 18:30 编辑

wangxue123 发表于 2024-5-31 17:47
感谢老师的回答，我用.com文件进行sa-casscf计算，显示
error
cannot find USERDEF-JFIT-FIT in basis ...

这跟MOKIT关系不大，你用了密度拟合加速的CASSCF（DF-CASSCF），没有指定辅助基/Molpro没能自动指定辅助基。要么不用密度拟合加速，要么自己指定辅助基。

建议使用MOKIT的automr自动做CASSCF计算，一个最简单的gjf文件、一个关键词mokit{CASSCF_prog=Molpro}即可自动调用Molpro，而且活性空间由automr自动构造，不需要Molpro的AVAS。例子见https://jeanwsr.gitlab.io/mokit-doc-mdbook/chap5-1.html

作者
Author: wangxue123 时间: 2024-6-3 21:25

zjxitcc 发表于 2024-5-31 18:27
这跟MOKIT关系不大，你用了密度拟合加速的CASSCF（DF-CASSCF），没有指定辅助基/Molpro没能自动指定辅助 ...

感谢老师的指导，我修改了.com文件,molpro已经能够正确计算了,我把输出的.xml文件用mokit执行xml2fch a.xml a.fch，发生如下报错。

(, 下载次数 Times of downloads: 25)
最开始用高斯做的hf是关闭了对称性的，然后mokit产生的.com拿去用molpro计算的时候我又加上了C2V对称性,不知道是不是这个原因。

作者
Author: zjxitcc 时间: 2024-6-3 21:30

wangxue123 发表于 2024-6-3 21:25
感谢老师的指导，我修改了.com文件,molpro已经能够正确计算了,我把输出的.xml文件用mokit执行xml2fch a.x ...

大概率是的。xml2fch目前不支持点群对称性。

作者
Author: wangxue123 时间: 2024-6-3 22:29
本帖最后由 wangxue123 于 2024-6-5 18:02 编辑

zjxitcc 发表于 2024-6-3 21:30
大概率是的。xml2fch目前不支持点群对称性。

好的，感谢老师的解答

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