求推荐 用ORCA计算三线态能级的合适泛函/基组

waker

计算 50~60 个原子的TADF 分子的 △Est， 感觉TDDFT RI-PBE0/def2-TZVP 的精度不够高， EOM-CCSD/cc-pVTZ 耗时又太高了。求推荐一款计算精度和时间在这两者之间的计算级别，谢谢各位大佬了

wzkchem5

用双杂化泛函（例如DSD-PBEP86之类，不过建议调研一下benchmark文章看看哪个双杂化泛函算 △Est最好，算吸收光谱最好的泛函和算 △Est最好的泛函未必一样），或者用DLPNO-STEOM-CCSD

sobereva

w调控的泛函（参考JCTC, 11, 3851），或者DLPNO-STEOM-CCSD

适当的时候也可以用JCTC, 9, 3872里的OHF方法

waker

wzkchem5 发表于 2022-11-10 14:57
用双杂化泛函（例如DSD-PBEP86之类，不过建议调研一下benchmark文章看看哪个双杂化泛函算 △Est最好，算吸 ...

谢谢老师，但是DSD-PBEP86计算的结果有点奇怪:
TD-DFT EXCITED STATES (SINGLETS)
--------------------------------

the weight of the individual excitations are printed if larger than 0.01

STATE  1:  E=   0.053866 au      1.466 eV    11822.3 cm**-1 <S**2> =   0.000000
   115a -> 124a  :     0.011711
   120a -> 124a  :     0.039021
   120a -> 125a  :     0.010770
   120a -> 126a  :     0.013939
   121a -> 124a  :     0.077703
   122a -> 124a  :     0.487945
   122a -> 125a  :     0.038192
   122a -> 131a  :     0.012851
   123a -> 124a  :     0.214839
   123a -> 125a  :     0.013183

STATE  2:  E=   0.116604 au      3.173 eV    25591.6 cm**-1 <S**2> =   0.000000
STATE  3:  E=   0.112692 au      3.066 eV    24733.0 cm**-1 <S**2> =   0.000000
STATE  4:  E=   0.070449 au      1.917 eV    15461.8 cm**-1 <S**2> =   0.000000
STATE  5:  E=   0.127391 au      3.466 eV    27959.2 cm**-1 <S**2> =   0.000000
TD-DFT EXCITED STATES (TRIPLETS)
--------------------------------

the weight of the individual excitations are printed if larger than 0.01

STATE  1:  E=   0.210781 au      5.736 eV    46261.0 cm**-1 <S**2> =   2.000000
   119a -> 124a  :     0.023481
   119a -> 127a  :     0.213869
   120a -> 124a  :     0.280465
   120a -> 125a  :     0.018828
   120a -> 126a  :     0.021303
   120a -> 128a  :     0.048942
   120a -> 129a  :     0.034150
   120a -> 131a  :     0.013389
   121a -> 124a  :     0.011493
   122a -> 124a  :     0.096781
   123a -> 124a  :     0.057408
   123a -> 125a  :     0.018727

STATE  2:  E=   0.230293 au      6.267 eV    50543.5 cm**-1 <S**2> =   2.000000
STATE  3:  E=   0.080199 au      2.182 eV    17601.6 cm**-1 <S**2> =   2.000000
STATE  4:  E=   0.135533 au      3.688 eV    29746.2 cm**-1 <S**2> =   2.000000
STATE  5:  E=   0.105578 au      2.873 eV    23171.7 cm**-1 <S**2> =   2.000000

这个能级顺序比较奇怪，请问是我的输入文件出了问题吗？文件如附件

wzkchem5

waker 发表于 2022-11-10 08:33
谢谢老师，但是DSD-PBEP86计算的结果有点奇怪:
TD-DFT EXCITED STATES (SINGLETS)
------- ...

双杂化泛函TDDFT是先做一个普通的TDDFT计算，再对激发能做一个校正(doubles correction)。这个校正会改变激发能的相对大小，但是为了方便和校正前的结果比较（因为振子强度、跃迁成分等都没有校正），程序没有对激发能做重新排序。所以双杂化泛函TDDFT激发能不按大小排序是正常的

waker

wzkchem5 发表于 2022-11-10 16:14
双杂化泛函TDDFT是先做一个普通的TDDFT计算，再对激发能做一个校正(doubles correction)。这个校正会改变 ...

谢谢老师，另外我在做DLPNO-STEOM-CCSD 时，程序报错了，是出了什么问题吗？Final no active IP roots:    5  [file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 13]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.

[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 9]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.


[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 2]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.


[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 8]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.


[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 6]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.

[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 4]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.


[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 10]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.


[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 14]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.


[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 3]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.


[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 1]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.

[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 5]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.

[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 0]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.

[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 7]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.

[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 12]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.


[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 11]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.

[file orca_mdci/mdci_rhf_dlpno_eom.cpp, line 3718, Process 15]: Error (ORCA_MDCI): Solution of the IP-EOM-DLPNO-CCSD equations failed.


job aborted:
[ranks] message

[0-15] process exited without calling finalize

---- error analysis -----

[0-15] on DESKTOP-QF559HM
C:\ORCA\orca_mdci_mpi.exe ended prematurely and may have crashed. exit code 125

---- error analysis -----

ORCA finished by error termination in MDCI
Calling Command: mpiexec -np 16  C:\ORCA\orca_mdci_mpi.exe DLPNO.mdciinp.tmp DLPNO
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run

wzkchem5

waker 发表于 2022-11-11 14:05
谢谢老师，另外我在做DLPNO-STEOM-CCSD 时，程序报错了，是出了什么问题吗？Final no active IP roots:   ...

仔细看一下输出文件是不是IP-EOM-DLPNO-CCSD迭代不收敛，如果是的话，按手册MDCI相关章节的介绍，加level shift、调整maxdiis等辅助收敛

huahua0621

wzkchem5 发表于 2022-11-10 16:14
双杂化泛函TDDFT是先做一个普通的TDDFT计算，再对激发能做一个校正(doubles correction)。这个校正会改变 ...

老师 那我想麻烦问下您，我用双杂化泛函TD-DFT计算了单点，但是计算后第一个电子态和之前不用双杂化用普通泛函计算的是一致的，我算了10个state，但是第十个才是能量最低的态，那我的最低激发态应该是哪个呢

wzkchem5

huahua0621 发表于 2023-8-9 02:17
老师 那我想麻烦问下您，我用双杂化泛函TD-DFT计算了单点，但是计算后第一个电子态和之前不用双杂化用普 ...

取第10个。但是要留意第10个态是真的本该能量低，还是因为微扰理论发散、doubles correction虚高而造成的能量虚低。参见Phys. Chem. Chem. Phys., 2009, 11, 4611–4620，在里面搜break-down

jingzhe111

waker 发表于 2022-11-10 15:33
谢谢老师，但是DSD-PBEP86计算的结果有点奇怪:
TD-DFT EXCITED STATES (SINGLETS)
------- ...

请教一下老师，像您这样用双杂化泛函进行tddft计算，同时计算单重态和三重态，想要得到该分子s1态和t1态的电子能量即单点能，应该在orca.out的哪个地方读取呢？谢谢老师