计算化学公社

标题: 求助：怎样将GAMESS的MO数据转成ORCA轨道初猜文件呢？ [打印本页]

作者
Author: 13799725819 时间: 2024-5-13 23:36
标题: 求助：怎样将GAMESS的MO数据转成ORCA轨道初猜文件呢？
本帖最后由 13799725819 于 2024-5-14 17:33 编辑

我想使用ORCA5.0.3进行CASSCF方面的计算，之前已经有了GAMESS的CASSCF计算输出的MO数据（*.dat）；
那么，我怎样能在使用ORCA做CASSCF计算时，不需要再次进行UNO计算，而直接利用GAMESS的CASSCF计算输出的MO数据呢？

作者
Author: zjxitcc 时间: 2024-5-13 23:49
本帖最后由 zjxitcc 于 2024-5-14 00:21 编辑

GAMESS传给ORCA：
首先要明确，GAMESS支持球谐型和Cartesian型基函数两种，而ORCA只支持球谐型基函数，因此我们要在GAMESS中也用前者进行计算。首先我们写一个gjf文件，例如

%chk=h2o.chk
%mem=4GB
%nprocshared=2
#p RHF/6-31G(d,p) 5D 7F nosymm int=nobasistransform
title
0 1
O 0.00000000 0.00000000 0.06200700
H 0.00000000 -0.78397600 -0.49205200
H 0.00000000 0.78397600 -0.49205200

复制代码

这里特地写了5D 7F表示使用球谐型基函数。6-31G(d,p)只是一个示例，你可以用其他基组如cc-pVTZ，也可以用赝势/自定义/混合基组。做完这个简单的Gaussian计算后有h2o.chk文件，运行

formchk h2o.chk h2o.fch
fch2inp h2o.fch

复制代码

产生GAMESS输入文件h2o.inp，内含坐标，基组信息和HF关键词，修改成你要的CASSCF关键词，或者你的多步CASSCF计算步骤，提交GAMESS计算，最终获得h2o.dat文件，运行

dat2fch h2o.dat h2o.fch

复制代码

这个h2o.fch前面我们就有了，这里拿来覆盖里面的轨道，这样h2o.fch文件中就是GAMESS收敛的CASSCF轨道，运行

fch2mkl h2o.fch
orca_2mkl h2o_o -gbw

复制代码

产生h2o_o.inp, h2o_o.mkl, h2o_o.gbw三个文件，在h2o_o.inp中添加你想要的关键词（坐标和基组都写好了，不用改，加关键词就行）。经过传轨道，三个量化软件的HF和CASSCF均可在1圈左右收敛，避免重复计算。这些步骤的精髓是经过一个fch文件进行中转；注意不能随便搞一个h2o.fch文件过来欺骗dat2fch小程序，那样基函数不对应，传轨道是错误的。

ORCA传给GAMESS：
我们知道ORCA的orca_2mkl小程序可以从gbw文件产生mkl文件或molden文件，比如有一文件h2o.gbw，操作如下（选择一种即可）

orca_2mkl h2o -mkl
orca_2mkl h2o -molden

复制代码

方法一：使用mkl2inp小程序，例如
mkl2inp h2o.inp

方法二：使用mkl2fch和fch2inp小程序，例如
mkl2fch h2o.mkl
fch2inp h2o.fch

方法三：使用molden2fch和fch2inp小程序，例如
molden2fch h2o.molden -orca
fch2inp h2o.fch

产生的GAMESS inp文件关键词是RHF/ROHF的，基组信息写好了，不需要自己添加基组。自己打开inp文件添加一下CASSCF的关键词，交上去算就行了。如果每一步都合理操作，GAMESS的CASSCF应该在3圈内收敛。mkl2inp, molden2fch, fch2inp这些都是开源程序MOKIT中的小程序，MOKIT还有十几个小程序用于各种量化软件传轨道，可以将常见量化软件联用。

注意，mkl和molden文件中均不含赝势信息，因此若你的ORCA计算使用赝势，那转化产生的inp文件也不含赝势（难为无米之炊，mkl和molden都不含，无法凭空变出）。为方便起见，你的ORCA计算最好使用全电子基组。如果实在对赝势有需求，可以再提问。

作者
Author: 13799725819 时间: 2024-5-14 11:08

zjxitcc 发表于 2024-5-13 23:49
GAMESS传给ORCA：
首先要明确，GAMESS支持球谐型和Cartesian型基函数两种，而ORCA只支持球谐型基函数，因 ...

谢谢！真是太全面了，我去试一试。

作者
Author: 413 时间: 2024-12-17 15:19
请问您用的哪个版本的gamess？

作者
Author: 13799725819 时间: 2024-12-24 19:59

413 发表于 2024-12-17 15:19
请问您用的哪个版本的gamess？

Gamess (2021-09-30 R2 Patch 1)

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