计算化学公社

标题: ORCA进行自旋破缺SCF计算时，有无直接输出自旋破缺单重态自然轨道的方法？ [打印本页]

作者
Author: 啊不错的飞过海 时间: 2023-3-23 13:57
标题: ORCA进行自旋破缺SCF计算时，有无直接输出自旋破缺单重态自然轨道的方法？
本帖最后由啊不错的飞过海于 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

作者
Author: wzkchem5 时间: 2023-3-23 17:44
可能是没有的，但可以考虑用compound job功能把broken symmetry计算和UNO/UCO计算写到一个输入文件里

作者
Author: hebrewsnabla 时间: 2023-3-23 17:57
或许你可以直接提feature request让他们把此种情况的UNO改成基于目标态而不是高自旋态的。这样比较符合直觉。

作者
Author: 啊不错的飞过海 时间: 2023-3-23 21:35

wzkchem5 发表于 2023-3-23 17:44
可能是没有的，但可以考虑用compound job功能把broken symmetry计算和UNO/UCO计算写到一个输入文件里

如果想写进一个文件的话倒是$new_job也可以，至于compound job，感觉不如python......灵便（逃）

作者
Author: 啊不错的飞过海 时间: 2023-3-23 21:51

hebrewsnabla 发表于 2023-3-23 17:57
或许你可以直接提feature request让他们把此种情况的UNO改成基于目标态而不是高自旋态的。这样比较符合直觉 ...

有道理，这就去他们论坛发帖

作者
Author: wzkchem5 时间: 2023-3-23 22:24

啊不错的飞过海发表于 2023-3-23 14:35
如果想写进一个文件的话倒是$new_job也可以，至于compound job，感觉不如python......灵便（逃）

对，对于这种简单的任务$new_job也可以，但是$new_job将来有可能会不支持，因为它默认把第一个任务的大部分选项carry over到第二个任务里，容易导致一些意外问题。compound job不会做这个carry over的步骤。

作者
Author: 啊不错的飞过海 时间: 2023-4-1 13:30

wzkchem5 发表于 2023-3-23 17:44
可能是没有的，但可以考虑用compound job功能把broken symmetry计算和UNO/UCO计算写到一个输入文件里

现在看来是有的，做波函数稳定性分析时BS波函数的自然轨道会被算一遍，具体情况更新在主楼了

作者
Author: 13799725819 时间: 2023-10-11 16:43
请问，为什么说“对于BS-HF/BS-DFT来讲做一下也还算方便，对于用双杂化泛函的任务就有些画蛇添足了”？麻烦给解释一下，谢谢~

作者
Author: zjxitcc 时间: 2023-10-11 17:00
本帖最后由 zjxitcc 于 2023-10-11 17:01 编辑

13799725819 发表于 2023-10-11 16:43
请问，为什么说“对于BS-HF/BS-DFT来讲做一下也还算方便，对于用双杂化泛函的任务就有些画蛇添足了”？麻烦 ...

双杂化任务没法直接检验波函数稳定性，需要单独提交一个任务，检验双杂化的自洽部分的波函数稳定性（无论Gaussian还是ORCA都是这样）。
https://jeanwsr.gitlab.io/mokit- ... 5.html#4524-fch2mkl
这里有一个Gaussian中检验PWPB95自洽部分的波函数稳定性的关键词示例。

作者
Author: 13799725819 时间: 2023-10-12 00:45

zjxitcc 发表于 2023-10-11 17:00
双杂化任务没法直接检验波函数稳定性，需要单独提交一个任务，检验双杂化的自洽部分的波函数稳定性（无 ...

谢谢！我看一看

作者
Author: 啊不错的飞过海 时间: 2023-10-17 14:33

13799725819 发表于 2023-10-11 16:43
请问，为什么说“对于BS-HF/BS-DFT来讲做一下也还算方便，对于用双杂化泛函的任务就有些画蛇添足了”？麻烦 ...

ORCA会在BS任务做完自洽场计算后给出其SCF高低自旋态能隙。
BS-HF/BS-DFT的能量就是SCF的能量，这个能隙是能用的；双杂化的能量还要加一个微扰——但做完SCF时微扰还没算，其SCF能量也没啥意义，此时输出的这个能隙就没什么用，只是浪费算力。

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