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

bangbangbang

各位老师好，想要分析S1态下，分子的芳香性，使用NICS方法。但不知道CASSCF如何计算。
根据论坛中Dalton的相关学习，将激发态的高斯计算文件gif,在Multiwfn中生.dal文件，编写了关于CASSCF的.mol文件。
想请教和麻烦老师，看一下这两个文件是可以进行CASSCF计算嘛，我提交后出现了报错，不知道是哪里出了问题？

zjxitcc

（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. 
   
5. mokit{NMR,root=1}
   
6. 
   
7. 0 1
   
8. C      0.00000000    0.00000000    0.00000000
   
9. C      0.00000000    0.00000000    1.34900000
   
10. C      1.56200000    0.00000000    1.34894100
    
11. C      1.56200000   -0.00023900   -0.00005900
    
12. H     -0.76022400   -0.00000200   -0.76288200
    
13. H     -0.76022400    0.00013500    2.11188200
    
14. H      2.32222400    0.00000000    2.11182300
    
15. 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

bangbangbang

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

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)。

sobereva

利用Multiwfn令Dalton计算时使用其它程序产生的轨道作为初猜
http://sobereva.com/740（http://bbs.keinsci.com/thread-52653-1-1.html）

此文里就有个CASSCF例子