计算化学公社

标题: 求推荐用ORCA计算三线态能级的合适泛函/基组 [打印本页]

作者
Author: waker 时间: 2022-11-10 14:16
标题: 求推荐用ORCA计算三线态能级的合适泛函/基组
计算 50~60 个原子的TADF 分子的 △Est，感觉TDDFT RI-PBE0/def2-TZVP 的精度不够高， EOM-CCSD/cc-pVTZ 耗时又太高了。求推荐一款计算精度和时间在这两者之间的计算级别，谢谢各位大佬了

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

作者
Author: sobereva 时间: 2022-11-10 15:01
w调控的泛函（参考JCTC, 11, 3851），或者DLPNO-STEOM-CCSD

适当的时候也可以用JCTC, 9, 3872里的OHF方法
(, 下载次数 Times of downloads: 30)

作者
Author: waker 时间: 2022-11-10 15:33

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

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

作者
Author: wzkchem5 时间: 2022-11-10 16:14

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

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

作者
Author: waker 时间: 2022-11-11 21:05
本帖最后由 waker 于 2022-11-11 21:07 编辑

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

作者
Author: wzkchem5 时间: 2022-11-11 21:10

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

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

作者
Author: huahua0621 时间: 2023-8-9 09:17

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

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

作者
Author: wzkchem5 时间: 2023-8-9 16:01

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

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

作者
Author: jingzhe111 时间: 2024-7-29 19:08

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

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

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