用MOKIT程序生成的输入文件

wangxue123 · 发表于 Post on 2024-5-14 22:08:10

本帖最后由 wangxue123 于 2024-5-15 16:49 编辑

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

zjxitcc · 发表于 Post on 2024-5-14 22:15:17

本帖最后由 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所支持。

wangxue123 · 发表于 Post on 2024-5-15 16:46:17

本帖最后由 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文件也行？

zjxitcc · 发表于 Post on 2024-5-15 16:52:42

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无关）

wangxue123 · 发表于 Post on 2024-5-31 17:47:08

本帖最后由 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

zjxitcc · 发表于 Post on 2024-5-31 18:27:07

本帖最后由 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

wangxue123 · 发表于 Post on 2024-6-3 21:25:53

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

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

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

zjxitcc · 发表于 Post on 2024-6-3 21:30:55

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

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

wangxue123 · 发表于 Post on 2024-6-3 22:29:56

本帖最后由 wangxue123 于 2024-6-5 18:02 编辑

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

好的，感谢老师的解答

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[辅助/分析程序] 用MOKIT程序生成的输入文件