大家在orca中使用双杂化泛函wB2PLYP的时候，w值的调控怎么操作

南京吴彦祖

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


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


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

谢谢大家。

sobereva

这种事应该发邮件问作者

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

wzkchem5

贴文献原文和你的输入输出文件上来，不然没人能解答为什么数据不一样。
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

南京吴彦祖

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

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

风起~

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

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

1. --------------------------------------------------------------------------------
   
2.     ROOT    Mult      Excitation energy[Eh]           [cm-1]          [eV]
   
3. --------------------------------------------------------------------------------
   
4.      0       3            0.00000000                      0.00          0.000
   
5.      1       3            0.02912859                   6392.99          0.793
   
6.      2       1            0.02989529                   6561.26          0.813
   
7.      3       3            0.04034258                   8854.17          1.098
   
8.      4       3            0.04390579                   9636.21          1.195
   
9.      5       1            0.05093341                  11178.59          1.386
   
10.      6       3            0.05329215                  11696.28          1.450
    
11.      7       1            0.05696040                  12501.36          1.550
    
12.      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

wzkchem5

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

这 ...

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

风起~

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

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

wzkchem5

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

什么叫能量最低点？

风起~

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

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

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

wzkchem5

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

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

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

风起~

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

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

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

复制代码

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

wzkchem5

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

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

风起~

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

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

1. ! PBE0 ROKS def2-svp def2-svp/c cpcm(toluene) tightSCF miniprint
   
2. %pal nprocs 36
   
3. end
   
4. %maxcore 8000
   
5. %rocis
   
6.   DoLowerMult true
   
7.   nroots 5
   
8.   DoDFTCIS true
   
9.   doRI true
   
10.   soc true
    
11. end
    
12. * xyz 0 3

复制代码

wzkchem5

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

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

风起~

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

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