求助gamess-us的CASSCF参数设置

tjeremiah

看到文献里写用gamess可以计算两个分子的CASSCF数据，以两个乙烯分子为例计算，但对CASSCF数据不太会设置。

比如我想用乙烯分子做CAS(4,4)应该如何设置参数？
我没找到应该设置（4，4）的部分，个人理解是应该是设置$DRT部分，但只知道NMCC=8 NDOC=8 NVAL=12，其他的就不会设置了，看手册也有点模糊。

输入文件见附件，请问应该如何修改？

zjxitcc

这未免有点麻烦。只需一个gjf文件，使用符合直觉的简单写法，MOKIT便可自动调用GAMESS做CASSCF计算，活性空间也不用指定，自动算出来就是CASSCF(4,4)。gjf文件如下

1. %mem=8GB
   
2. %nprocshared=4
   
3. #p CASSCF/cc-pVDZ
   
4. 
   
5. mokit{CASSCF_prog=GAMESS}
   
6. 
   
7. 0 1
   
8. C   -2.3280701600     -0.1699583400     -1.8514691000
   
9. H   -1.7949064100     -1.0976632600     -1.8514691000
   
10. H   -3.3980701600     -0.1699583400     -1.8514691000
    
11. C   -2.4237513400      0.6322979400      0.6618613600
    
12. H   -1.8905876000     -0.2954069800      0.6618613600
    
13. H   -3.4937513400      0.6322979400      0.6618613600
    
14. C   -1.6527958500      1.0050189600     -1.8514691000
    
15. H   -2.1859596000      1.9327238800     -1.8514691000
    
16. H   -0.5827958500      1.0050189600     -1.8514691000
    
17. C   -1.7484770400      1.8072752300      0.6618613600
    
18. H   -2.2816407800      2.7349801500      0.6618613600
    
19. H   -0.6784770400      1.8072752300      0.6618613600
    

复制代码

提交

1. automr c2h4_dimer.gjf >c2h4_dimer.out 2>&1

复制代码

中间无需停下来看轨道、调换轨道等操作。4核并行，2min算完，当前目录下就有fch文件c2h4_dimer_uhf_gvb12_CASSCF_NO.fch，内含CASSCF自然轨道及其占据数，可直接用GaussView或Multiwfn打开可视化。计算文件见附件 c2h4_dimer_GMS.tar.gz (1.34 MB, 下载次数 Times of downloads: 11)

2023-9-19 11:28 上传 Uploaded

点击下载
Click to download

输出文件里有各种能量和未成对电子数

1. E(CASCI)  =      -156.09059822 a.u.
   
2. E(CASSCF) =      -156.09202343 a.u.
   
3. ----------------------- Radical index -----------------------
   
4. biradical character   (2c^2) y0=  0.111
   
5. tetraradical character(2c^2) y1=  0.095
   
6. Yamaguchi's unpaired electrons  (sum_n n(2-n)      ):  0.781
   
7. Head-Gordon's unpaired electrons(sum_n min(n,(2-n))):  0.412
   
8. Head-Gordon's unpaired electrons(sum_n (n(2-n))^2  ):  0.155
   
9. -------------------------------------------------------------

复制代码

若不写CASSCF_prog=GAMESS，文件更简单，默认自动调用PySCF做CASSCF计算，更快（结果还是严格一样的）。激发态CASSCF计算等更多例子见https://jeanwsr.gitlab.io/mokit-doc-mdbook/chap5-1.html

tjeremiah

zjxitcc 发表于 2023-9-19 11:28
这未免有点麻烦。只需一个gjf文件，使用符合直觉的简单写法，MOKIT便可自动调用GAMESS做CASSCF计算，活性空 ...

谢谢老师的回复。
我主要是想用这篇文章后面的4-fold 方法得到透热态另外输出各态之间的耦合项。但gamess这部分只有CISTEP=GUGA可以采用，这个的关键词是$DRT;但我看cas(4,4)是用的$det和CISTEP=ALDET，其中的NELS=4 NACT=4没在GUGA对应的部分找到关键词，不知道这之间应该怎么转换？

beefly

gamess输入太繁琐了，它处理cc-pvnz类型基组的积分效率也很低，casscf这部分的算法太陈旧，并且手册里很多陈年旧账干扰初学者。建议改用免费的opemmolcas做casscf（或者还有最新版orca？好几年没用orca做casscf计算了，不便评论），它的casscf结合密度拟合以及Cholesky分解，能算很大的分子体系。

个人不推荐gamess的另一个原因，是它的一个作者mike schmidt擅长玩政治正确，五毒俱全。也许某天看中国用户不顺眼，把你骂一顿，不光是带点敏感内容啥的，那居高临下的姿态也让人反感。不过，Mark Gordon是个好大爷。

chrinide

Mike已经被清除出GA队伍了 @beefly