计算化学公社

标题: 求助gamess-us的CASSCF参数设置 [打印本页]
作者
Author: tjeremiah    时间: 2023-9-19 06:55
标题: 求助gamess-us的CASSCF参数设置
看到文献里写用gamess可以计算两个分子的CASSCF数据，以两个乙烯分子为例计算，但对CASSCF数据不太会设置。
(, 下载次数 Times of downloads: 23)
比如我想用乙烯分子做CAS(4,4)应该如何设置参数？
我没找到应该设置（4，4）的部分，个人理解是应该是设置$DRT部分，但只知道NMCC=8 NDOC=8 NVAL=12，其他的就不会设置了，看手册也有点模糊。
(, 下载次数 Times of downloads: 25)
输入文件见附件，请问应该如何修改？
作者
Author: zjxitcc    时间: 2023-9-19 11:28
本帖最后由 zjxitcc 于 2023-9-19 11:40 编辑

这未免有点麻烦。只需一个gjf文件，使用符合直觉的简单写法，MOKIT便可自动调用GAMESS做CASSCF计算，活性空间也不用指定，自动算出来就是CASSCF(4,4)。gjf文件如下
  1. %mem=8GB
  2. %nprocshared=4
  3. #p CASSCF/cc-pVDZ

  5. mokit{CASSCF_prog=GAMESS}

  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打开可视化。计算文件见附件 (, 下载次数 Times of downloads: 12)

输出文件里有各种能量和未成对电子数
  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
作者
Author: tjeremiah    时间: 2023-9-19 16:04
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对应的部分找到关键词，不知道这之间应该怎么转换？
作者
Author: beefly    时间: 2023-9-19 17:15
gamess输入太繁琐了，它处理cc-pvnz类型基组的积分效率也很低，casscf这部分的算法太陈旧，并且手册里很多陈年旧账干扰初学者。建议改用免费的opemmolcas做casscf（或者还有最新版orca？好几年没用orca做casscf计算了，不便评论），它的casscf结合密度拟合以及Cholesky分解，能算很大的分子体系。

个人不推荐gamess的另一个原因，是它的一个作者mike schmidt擅长玩政治正确，五毒俱全。也许某天看中国用户不顺眼，把你骂一顿，不光是带点敏感内容啥的，那居高临下的姿态也让人反感。不过，Mark Gordon是个好大爷。
作者
Author: chrinide    时间: 2023-9-20 06:38
本帖最后由 chrinide 于 2023-10-20 11:47 编辑

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



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3