计算化学公社

 找回密码 Forget password
 注册 Register
Views: 673|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[GAMESS-US] 求助gamess-us的CASSCF参数设置

[复制链接 Copy URL]

9

帖子

0

威望

401

eV
积分
410

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
看到文献里写用gamess可以计算两个分子的CASSCF数据,以两个乙烯分子为例计算,但对CASSCF数据不太会设置。

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

输入文件见附件,请问应该如何修改?

c2h4cas.inp

94.45 KB, 下载次数 Times of downloads: 2

339

帖子

0

威望

4999

eV
积分
5338

Level 6 (一方通行)

5#
发表于 Post on 2023-9-20 06:38:06 来自手机 | 只看该作者 Only view this author
本帖最后由 chrinide 于 2023-10-20 11:47 编辑

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

689

帖子

21

威望

5019

eV
积分
6128

Level 6 (一方通行)

4#
发表于 Post on 2023-9-19 17:15:52 | 只看该作者 Only view this author
gamess输入太繁琐了,它处理cc-pvnz类型基组的积分效率也很低,casscf这部分的算法太陈旧,并且手册里很多陈年旧账干扰初学者。建议改用免费的opemmolcas做casscf(或者还有最新版orca?好几年没用orca做casscf计算了,不便评论),它的casscf结合密度拟合以及Cholesky分解,能算很大的分子体系。

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

9

帖子

0

威望

401

eV
积分
410

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-9-19 16:04:37 | 只看该作者 Only view this author
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对应的部分找到关键词,不知道这之间应该怎么转换?

3813

帖子

4

威望

8002

eV
积分
11895

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2023-9-19 11:28:17 | 只看该作者 Only view this author
本帖最后由 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

  4. mokit{CASSCF_prog=GAMESS}

  5. 0 1
  6. C   -2.3280701600     -0.1699583400     -1.8514691000
  7. H   -1.7949064100     -1.0976632600     -1.8514691000
  8. H   -3.3980701600     -0.1699583400     -1.8514691000
  9. C   -2.4237513400      0.6322979400      0.6618613600
  10. H   -1.8905876000     -0.2954069800      0.6618613600
  11. H   -3.4937513400      0.6322979400      0.6618613600
  12. C   -1.6527958500      1.0050189600     -1.8514691000
  13. H   -2.1859596000      1.9327238800     -1.8514691000
  14. H   -0.5827958500      1.0050189600     -1.8514691000
  15. C   -1.7484770400      1.8072752300      0.6618613600
  16. H   -2.2816407800      2.7349801500      0.6618613600
  17. 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)

输出文件里有各种能量和未成对电子数
  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
自动做多参考态计算的程序MOKIT

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 16:22 , Processed in 0.459226 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list