|
|
本帖最后由 啊不错的飞过海 于 2023-4-1 13:32 编辑
ORCA做一般的UHF/UKS计算时，可以通过!UNO关键字简便地让程序在SCF结束后自动生成其自然轨道，还可以通过%print指定使其输出自然轨道占据数等信息，非常方便。进行使用%scf Brokensym选项做BS-UHF或BS-UKS计算时，ORCA还默认在对称破缺态下生成其UCO(UHF Corresponding Orbital，不会译)并输出UCO重叠情况。
然而，结合使用%scf Brokensym选项与!UCO关键字时，ORCA只会输出高自旋态的自然轨道，此时在%print段指定输出的自然轨道占据数等信息也是基于高自旋轨道的；想要得到对称破缺态的自然轨道信息，就只能新开一个非限制单重态任务（不使用%scf Brokensym选项），然后用对称破缺的轨道给它做初猜，重做一遍SCF，比较麻烦。ORCA手册中似乎也没有写怎么让使用%scf Brokensym选项的任务输出其自旋破缺单重态下自然轨道的写法。
有没有什么办法让ORCA直接在一个任务中做自旋破缺计算并输出自旋破缺单重态自然轨道呢？
===========================================
2023.4.1更新：
方法找到了，开启波函数稳定性分析时会默认给BS轨道算UNO/UCO。缺陷是BS-波函数稳定性分析联用时高自旋轨道的稳定性也会被检查，用时上比较亏。
ORCA这样做默认设定的意图似乎是在一个任务中直接算出来BS任务的高低自旋能隙，高自旋的做一次，BS的再做一次，保证输出的能隙是准的。我的任务倒是确实需要算这个能隙，对于BS-HF/BS-DFT来讲做一下也还算方便，对于用双杂化泛函的任务就有些画蛇添足了。
===========================================
附-只输出高自旋态下UNO的ORCA输入文件头样例：
! wB97M-V UNO UCO def2-TZVP def2/J def2-TZVP/C RIJCOSX tightSCF miniprint nopop
%maxcore 1000
%pal nprocs 8 end
%scf
BrokenSym 1,1
end
%output
printlevel mini
print[P_natpop] 1
print[P_UNO_OccNum] 1
print[P_OrbEn] 0
end
|
|
发表于 Post on 2023-3-23 13:57:29
收藏 Add to favorites
楼主 Author