|
还有另一种途径是不学习Molpro各种语法和功能,直接采用MOKIT程序黑箱、自动式的多参考计算,这里以两根O-H键均被拉长的水分子为例,写一个高斯gjf文件
%mem=4GB
%nprocshared=2
#p CASPT2/cc-pVDZ
mokit{CASPT2_prog=Molpro}
0 1
O -0.23497692 0.90193619 -0.068688
H 1.26502308 0.90193619 -0.068688
H -0.73568721 2.31589843 -0.068688
输入文件使用最常用的gjf格式,想要什么理论方法直接写(哪怕高斯不支持该方法也行,因为并不是仅调用高斯算的),符合书写直觉。中间CASPT2_prog=Molpro表示调用Molpro做CASPT2计算,一样符合直觉。用automr提交,命令如下
automr h2o.gjf >h2o.out 2>&1
仅20秒结束,不用肉眼看轨道、挑选轨道、做轨道局域化、大小基组投影、担心收敛性、调试关键词等各种曲折操作,直接算出来就是CASSCF(4,4)和CASPT2(4,4),程序会自动确定活性轨道是两根O-H的成键、反键轨道。h2o.out输出文件里有各个方法的电子能量、双自由基特征、未成对电子数等大量信息,直观清楚
E(CASCI) = -75.90802109 a.u.
E(CASSCF) = -75.90823063 a.u.
----------------------- Radical index -----------------------
biradical character (2c^2) y0= 0.165
tetraradical character(2c^2) y1= 0.139
Yamaguchi's unpaired electrons (sum_n n(2-n) ): 1.121
Head-Gordon's unpaired electrons(sum_n min(n,(2-n))): 0.607
Head-Gordon's unpaired electrons(sum_n (n(2-n))^2 ): 0.315
-------------------------------------------------------------
Leave subroutine do_cas at Sat Aug 12 14:46:05 2023
Enter subroutine do_mrpt2...
CASPT2 based on CASSCF orbitals.
Frozen_core = F, CASPT2 using program molpro
$molpro -t 1 -n 2 -m 250m h2o_uhf_gvb4_CASSCF_CASPT2.com
E(ref) = -75.90823031 a.u.
IP-EA shift = 0.25 (default)
E(corr) = -0.13847648 a.u.
E(CASPT2) = -76.04670678 a.u.
Leave subroutine do_mrpt2 at Sat Aug 12 14:46:16 2023
另有fch文件可用GaussView和Multiwfn直接打开(不管你用什么量化程序算,都会给出fch文件,十分通用)
CASSCF初始轨道 h2o_uhf_uno_asrot2gvb4_s.fch
CASSCF收敛轨道,也是自然轨道 h2o_uhf_gvb4_CASSCF_NO.fch
MOKIT程序免费、开源,见https://gitlab.com/jxzou/mokit
中文简介见https://gitlab.com/jxzou/mokit/-/blob/master/README_zh.md |
|