计算化学公社

标题: 大家在orca中使用双杂化泛函wB2PLYP的时候，w值的调控怎么操作 [打印本页]

作者
Author: 南京吴彦祖 时间: 2023-2-5 22:56
标题: 大家在orca中使用双杂化泛函wB2PLYP的时候，w值的调控怎么操作
本帖最后由南京吴彦祖于 2023-2-6 14:17 编辑

因为最近看了一篇文章，文章中使用的wB2PLYP，但是不知道什么原因，我直接使用wB2PLYP,计算的结果该文章中的数据不一样。
所以想着是不是文章中调控了w值，但是看了si以及正文都没有发现作者使用的值是多少。

其次，那篇文章的作者还使用了SA-SF-PBE50，这个泛函在orca中应该怎么输入。

最后，大家认为在计算七嗪类的分子，双杂化泛函中较好的较精确的泛函是什么呢？

谢谢大家。

作者
Author: sobereva 时间: 2023-2-6 00:54
这种事应该发邮件问作者

估计没有人对这么具体的分子特意做过双杂化泛函的横测。心里没数的话可以用DLPNO-STEOM-CCSD算

作者
Author: wzkchem5 时间: 2023-2-6 05:00
贴文献原文和你的输入输出文件上来，不然没人能解答为什么数据不一样。
SA-SF-PBE50应该不是个泛函，而可能是用SA-SF-CIS方法结合PBE50泛函算的，就好比TD-B3LYP指的是用B3LYP做TDDFT，不是说TD-B3LYP是一个新的泛函的名字。SA-SF-CIS只能用Q-Chem做，不过orca里的DFT/ROCIS方法可以代替SA-SF-CIS，结果不严格等价，但应该相似。关于SA-SF-CIS和DFT/ROCIS的区别和联系，参见https://aip.scitation.org/doi/10.1063/5.0025428

作者
Author: 南京吴彦祖 时间: 2023-2-6 14:18

wzkchem5 发表于 2023-2-6 05:00
贴文献原文和你的输入输出文件上来，不然没人能解答为什么数据不一样。
SA-SF-PBE50应该不是个泛函，而可 ...

好的，谢谢指正，我已经放上了文献原文，以及自己计算的输入输出文件，请您看一下，谢谢。

作者
Author: 风起~ 时间: 2024-12-19 11:50

wzkchem5 发表于 2023-2-6 05:00
贴文献原文和你的输入输出文件上来，不然没人能解答为什么数据不一样。
SA-SF-PBE50应该不是个泛函，而可 ...

请问一下，我在用rocis计算一个基态为单重态的分子时，是否入sf计算一样要将自旋多重度设置为3呢？

--------------------------------------------------------------------------------
ROOT Mult Excitation energy[Eh] [cm-1] [eV]
--------------------------------------------------------------------------------
0 3 0.00000000 0.00 0.000
1 3 0.02912859 6392.99 0.793
2 1 0.02989529 6561.26 0.813
3 3 0.04034258 8854.17 1.098
4 3 0.04390579 9636.21 1.195
5 1 0.05093341 11178.59 1.386
6 3 0.05329215 11696.28 1.450
7 1 0.05696040 12501.36 1.550
8 1 0.05932491 13020.31 1.614

复制代码

这里的root2 的单重态激发能是相对于单重基态还是三重态呢？因为0.813 eV这个激发能明显太小。
而且这个分子的单三态能差是否就是0.813-0.793=0.02 eV呢
麻烦您指点一下

相关输出文件在云盘，麻烦您了
通过百度网盘分享的文件：dabna-1.out
链接：https://pan.baidu.com/s/15AVtwqSwo4592IA_IGrYRA?pwd=zcr1
提取码：zcr1

作者
Author: wzkchem5 时间: 2024-12-19 16:31

风起~ 发表于 2024-12-19 04:50
请问一下，我在用rocis计算一个基态为单重态的分子时，是否入sf计算一样要将自旋多重度设置为3呢？

这 ...

看Mult那一栏就知道谁是单重态谁是三重态了，Mult是Multiplicity的缩写
是否需要把参考态自旋多重度设为3，取决于S0和T1哪个容易用DFT算准

作者
Author: 风起~ 时间: 2024-12-19 16:47

wzkchem5 发表于 2024-12-19 16:31
看Mult那一栏就知道谁是单重态谁是三重态了，Mult是Multiplicity的缩写
是否需要把参考态自旋多重度设为 ...

‘是否需要把参考态自旋多重度设为3’
所以这个激发能是相对于能量最低点而言的么

作者
Author: wzkchem5 时间: 2024-12-19 17:25

风起~ 发表于 2024-12-19 09:47
‘是否需要把参考态自旋多重度设为3’
所以这个激发能是相对于能量最低点而言的么

什么叫能量最低点？

作者
Author: 风起~ 时间: 2024-12-19 19:26
本帖最后由风起~ 于 2024-12-19 19:37 编辑

wzkchem5 发表于 2024-12-19 17:25
什么叫能量最低点？

就是自旋多重度为 1或3 的基态的能量更低点

rocis计算基于roks的情况下，不能设置为闭壳层。如果基于rhf，就是普通的cis。
所以我基于roks，设置自旋多重度为3
但是这个激发能是相对于那个一点而言自旋多重度为3的情况么还是能量更低的自旋多重度为1的情况

作者
Author: wzkchem5 时间: 2024-12-19 20:10

风起~ 发表于 2024-12-19 12:26
就是自旋多重度为 1或3 的基态的能量更低点

rocis计算基于roks的情况下，不能设置为闭壳层。如果基于r ...

问题是你用ROCIS算出来能量最低的态是三重态啊。你是如何得到这个分子的基态是单重态的结论的？

作者
Author: 风起~ 时间: 2024-12-19 20:28

wzkchem5 发表于 2024-12-19 20:10
问题是你用ROCIS算出来能量最低的态是三重态啊。你是如何得到这个分子的基态是单重态的结论的？

原来这个会自动判断基态的自旋多重度

--------------------------------------------------------------------------------
ROOT Mult Excitation energy[Eh] [cm-1] [eV]
--------------------------------------------------------------------------------
0 3 0.00000000 0.00 0.000

复制代码

这个分子是MR-TADF最经典的一个分子，居然出现这种情况真的是震惊

作者
Author: wzkchem5 时间: 2024-12-19 21:55

风起~ 发表于 2024-12-19 13:28
原来这个会自动判断基态的自旋多重度

是不是因为你用的不是DFT/ROCIS？
单用ROCIS，会严重低估T1能量，甚至导致T1低于S0，这个和HF低估T1能量是一个道理

作者
Author: 风起~ 时间: 2024-12-20 10:02

wzkchem5 发表于 2024-12-19 21:55
是不是因为你用的不是DFT/ROCIS？
单用ROCIS，会严重低估T1能量，甚至导致T1低于S0，这个和HF低估T1能量 ...

我用的确实是DFTROCIS，好奇怪的结果

! PBE0 ROKS def2-svp def2-svp/c cpcm(toluene) tightSCF miniprint
%pal nprocs 36
end
%maxcore 8000
%rocis
DoLowerMult true
nroots 5
DoDFTCIS true
doRI true
soc true
end
* xyz 0 3

复制代码

作者
Author: wzkchem5 时间: 2024-12-20 15:20

风起~ 发表于 2024-12-20 03:02
我用的确实是DFTROCIS，好奇怪的结果

是不是基态波函数没有收敛到能量最低的波函数？

作者
Author: 风起~ 时间: 2024-12-20 15:59

wzkchem5 发表于 2024-12-20 15:20
是不是基态波函数没有收敛到能量最低的波函数？

可能是的
那么rocis也太不稳定了吧

作者
Author: wzkchem5 时间: 2024-12-20 18:45

风起~ 发表于 2024-12-20 08:59
可能是的
那么rocis也太不稳定了吧

如果真是这样的话，那么问题不是出在ROCIS，而是出在ROKS。ROKS本来就容易收敛到不稳定波函数
可以以UKS的QRO轨道做初猜来做ROKS试试，乃至于可以直接基于UKS的QRO轨道做ROCIS

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