求助：怎样将GAMESS的MO数据转成ORCA轨道初猜文件呢？

13799725819

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

zjxitcc

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

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

复制代码

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

1. formchk h2o.chk h2o.fch
   
2. fch2inp h2o.fch

复制代码

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

1. dat2fch h2o.dat h2o.fch

复制代码

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

1. fch2mkl h2o.fch
   
2. 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，操作如下（选择一种即可）

1. orca_2mkl h2o -mkl
   
2. 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计算最好使用全电子基组。如果实在对赝势有需求，可以再提问。

13799725819

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

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