计算化学公社

标题: 高斯16用alter关键词交换轨道不成功 [打印本页]

作者
Author: BinWang 时间: 2022-11-8 17:03
标题: 高斯16用alter关键词交换轨道不成功
本帖最后由 BinWang 于 2022-11-8 20:22 编辑

各位老师：大家好！我用高斯做CASSCF时用alter交换轨道不成功怎么办，输入文件如下：
%nproc=32
%mem=50GB
%chk=TS.chk
# hf/sto-3g nosym

formodehide

0 1
C 6.278396925488155 0.3296048074971284 0.1612059285632868
C 7.788403787818916 0.5894418134350702 -0.1311618416620294
C 7.762495217934132 2.139420268042496 -0.3533564496557764
C 5.867968047740744 1.788195467503076 0.06338917681079459
C 6.196403650152662 2.476439572059121 -1.102198903619217
C 5.467184660592673 2.565136345128419 1.26644842529743
H 5.802371083583121 -0.2989234239358279 -0.587199286107845
H 6.091460354252916 -0.08753628353557222 1.14870194114183
H 8.117514641046018 0.07483849552960756 -1.030672031911053
H 8.447328911457912 0.2974198555700084 0.6865224178579373
H 8.372922069894065 2.489716671980996 -1.17789521060085
H 7.991582628918603 2.728250960785596 0.5295288821588113
H 5.604255692400966 3.633070420302092 1.094183981563241
H 4.399192830935121 2.391604210781828 1.437408801334604
H 5.995501971207139 2.242717735017672 2.167543925466545
F 6.180259124906911 1.624486184194161 -2.207691519403577

--Link1--
%oldchk=TS.chk
%chk=TS_Boys.chk
%mem=50GB
%nprocshared=32
#p hf/sto-3g guess(read,local,only,save) geom=allcheck nosym

0 1

接着有：
%oldchk=TS_Boys.chk
%chk=TS_cas_avdz_change_obi.chk
%mem=50GB
%nprocshared=32
#p cas(12,9)/sto-3g guess(read,alter) geom=Allcheck nosym

0 1

30,44
9,23
8,22

作者
Author: zjxitcc 时间: 2022-11-8 19:25
本帖最后由 zjxitcc 于 2022-11-8 19:29 编辑

cas(12,9)/sto-3g guess(read,alter)表示交换轨道后进行CASSCF(12,9)计算。你的轨道确实被交换了，只不过它要继续进行CASSCF轨道优化，至于轨道优化结果如何，并非一定如你所愿，跟你的初始轨道、轨道优化算法有关。如果你写guess(read,alter,only,save)，它在交换轨道后会退出程序、不做计算，可以看到你的轨道真的被交换了，进而说明高斯程序没有问题。另外顺便提醒，CASSCF活性轨道互相之间是没有顺序的，注意不要比较活性轨道谁序号小、谁序号大，这是CASSCF方法的酉不变性带来的。

这个从RHF/STO-3G开始，局域化、手动看轨道，调换轨道，期望小基组收敛，最后给大基组读取，步骤繁琐，人为因素参与多。不如直接用MOKIT算一下，输入文件如下

%mem=48GB
%nprocshared=48
#p CASSCF/cc-pVDZ
mokit{}
0 1
坐标

复制代码

计算过程是自动的，不用中途停下来看、挑，直接看结果：程序自动确定当前结构下活性空间为CASSCF(4,4)，同时获得自然轨道文件*_CASSCF_NO.fch，4个自然轨道为
(, 下载次数 Times of downloads: 8)

看起来是C-F成键、反键轨道，以及C-C键和C=C键的混合，轨道下方数值为自然轨道占据数。这个体系看起来(4,4)或(6,6)就足够了，(12,9)未免过于耗时。想算(6,6)就把上述输入文件关键词改成CASSCF(6,6)，程序会自动给出结果。

计算相关文件压缩包 (, 下载次数 Times of downloads: 2)

作者
Author: BinWang 时间: 2022-11-8 20:18

zjxitcc 发表于 2022-11-8 19:25
cas(12,9)/sto-3g guess(read,alter)表示交换轨道后进行CASSCF(12,9)计算。你的轨道确实被交换了，只不过它 ...

非常感谢老师的回复

欢迎光临计算化学公社 (http://bbs.keinsci.com/)