计算化学公社

 找回密码 Forget password
 注册 Register

如何使用Dalton做分子S1态的CASSCF计算?

查看数: 397 | 评论数: 4 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-4-8 14:58

正文摘要:

本帖最后由 bangbangbang 于 2025-4-22 09:52 编辑 各位老师好,想要分析S1态下,分子的芳香性,使用NICS方法。但不知道CASSCF如何计算。 根据论坛中Dalton的相关学习,将激发态的高斯计算文件gif,在Multiwfn中 ...

回复 Reply

sobereva 发表于 Post on 2025-4-9 03:51:14
利用Multiwfn令Dalton计算时使用其它程序产生的轨道作为初猜
http://sobereva.com/740http://bbs.keinsci.com/thread-52653-1-1.html

此文里就有个CASSCF例子
zjxitcc 发表于 Post on 2025-4-8 17:11:27
bangbangbang 发表于 2025-4-8 16:57
好的老师,十分感谢您,我装MOKIT尝试算一下呢。
有几个问题,想请教一下您,麻烦您看一下
(1)我已经 ...

安装MOKIT, GAMESS和PySCF。正常编译即可,这些软件的编译教程早就写得较为成熟了,可以参考
GAMESS编译教程》,《离线安装PySCF-2.x》,《离线安装量子化学软件Dalton

Gaussian可以有,也可以没有。如果不想用,可以看《MOKIT自动多参考计算不调用Gaussian的技巧》。

活空间的数值“CASSCF(4,4)” 字词混乱 看不懂。MOKIT可以自动构造活性轨道,自动确定活性空间。如果您没做过CASSCF计算,最好是让MOKIT自动确定,而非强行指定。上述环丁二烯的例子中,不写(4,4)也可以,MOKIT经过计算会自动判断出CASSCF(4,4)。
bangbangbang 发表于 Post on 2025-4-8 16:57:29
zjxitcc 发表于 2025-4-8 15:15
(1)是gjf,不是gif。
(2)您压缩包中的KS-DFT计算对于CASSCF计算几乎没有帮助。仅能提供体系的坐标信息 ...

好的老师,十分感谢您,我装MOKIT尝试算一下呢。
有几个问题,想请教一下您,麻烦您看一下
(1)我已经在虚拟机里装了Dalton,“只需写好gjf文件便可自动构造CASSCF活性轨道,自动调用PySCF做CASSCF轨道优化,传轨道给Dalton进行NMR计算”,我需要安装MOKIT的同时,需要安装PySCF吧?
(2)它们之间需要怎么样的编译,才可以调用嘛?
(3)还想请教您,对于活空间的数值“CASSCF(4,4)”,我应如何获取呢?
zjxitcc 发表于 Post on 2025-4-8 15:15:11
本帖最后由 zjxitcc 于 2025-4-8 17:17 编辑

(1)是gjf,不是gif。
(2)您压缩包中的KS-DFT计算对于CASSCF计算几乎没有帮助。仅能提供体系的坐标信息。
(3)CASSCF不是黑箱方法,它需要提供一组良好的初始轨道,往往无法从RHF轨道开始计算,也不能“提供一组轨道,写个活性空间大小就开始CASSCF计算”。
(4)Dalton的CASSCF轨道优化比较慢,用它做CASSCF计算比较费时费事,建议只用它做NMR计算(就是说算NICS),而把CASSCF轨道优化 外包给更强的程序算。

建议使用MOKIT自动做CASSCF NICS计算,该程序免费开源,可联网在线安装(见安装方式)。用户只需写好gjf文件便可自动构造CASSCF活性轨道,自动调用PySCF做CASSCF轨道优化,传轨道给Dalton进行NMR计算,使Dalton的CASSCF极速收敛,从而迅速进入NMR计算步骤,获得NICS。例如,以下是一个环丁二烯的CASSCF S1态NMR计算示例
  1. %mem=32GB
  2. %nprocshared=16
  3. #p CASSCF/cc-pVDZ

  4. mokit{NMR,root=1}

  5. 0 1
  6. C      0.00000000    0.00000000    0.00000000
  7. C      0.00000000    0.00000000    1.34900000
  8. C      1.56200000    0.00000000    1.34894100
  9. C      1.56200000   -0.00023900   -0.00005900
  10. H     -0.76022400   -0.00000200   -0.76288200
  11. H     -0.76022400    0.00013500    2.11188200
  12. H      2.32222400    0.00000000    2.11182300
  13. H      2.32222400   -0.00037400   -0.76294100
复制代码
提交任务,运行
  1. automr C4H4.gjf >C4H4.out 2>&1 &
复制代码
MOKIT会自动调用Gaussian进行RHF/UHF计算,对这个例子会产生UNO轨道,变换得到局域配对UNO;调用GAMESS进行GVB计算,判断哪些轨道重要(将作为CASSCF活性轨道);调用PySCF进行CASSCF轨道优化,传轨道给Dalton进行CASSCF NMR计算。过程中无需自己手写GAMESS/PySCF/Dalton输入文件,无需自己手动传轨道。MOKIT还支持CASSCF ICSS计算,更多例子可以看https://jeanwsr.gitlab.io/mokit-doc-mdbook/chap5-4.html

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

GMT+8, 2025-8-15 20:30 , Processed in 0.158356 second(s), 26 queries , Gzip On.

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