计算化学公社

标题: ORCA和Gaussian中TDDFT三重态计算结果相差大 [打印本页]

作者
Author: zhouoh 时间: 2021-1-24 23:03
标题: ORCA和Gaussian中TDDFT三重态计算结果相差大
各位老师们好，我最近在ORCA和Gaussian中分别用M062X泛函计算了蒽分子的三重激发态能级，但是两个软件的结果相差很大。但是单重态的结果却是相近的，所以比较疑惑，请各位看看是哪里出了问题。谢谢！

作者
Author: niobium 时间: 2021-1-25 09:15
我觉得这差别也不至于用“很大”来描述。真要比的话把ORCA的RIJCOSX去掉再比一比。这么小的体系没必要用RIJCOSX。

作者
Author: zhouoh 时间: 2021-1-25 09:51

niobium 发表于 2021-1-25 09:15
我觉得这差别也不至于用“很大”来描述。真要比的话把ORCA的RIJCOSX去掉再比一比。这么小的体系没必要用RIJ ...

谢谢回复！我去掉RIJCOSX之后重新算了一遍，发现相差还是有点大。主要是Gaussian计算的三重态能量普遍比ORCA的高出接近0.2eV，这我感觉已经不在误差范围内了。

作者
Author: sobereva 时间: 2021-1-25 22:33
计算更多的根试试，并且只让Gaussian计算三重态试试

作者
Author: zhouoh 时间: 2021-1-26 10:20

sobereva 发表于 2021-1-25 22:33
计算更多的根试试，并且只让Gaussian计算三重态试试

谢谢sob老师回复！我又计算了20个态，同时让Gaussian只计算三重态，结果还是没有改善。

作者
Author: sobereva 时间: 2021-1-26 19:22

zhouoh 发表于 2021-1-26 10:20
谢谢sob老师回复！我又计算了20个态，同时让Gaussian只计算三重态，结果还是没有改善。

原因不好说。换个泛函比如PBE0，或者都用TDA再比比试试。始终不开RI

作者
Author: huahua0621 时间: 2022-2-9 22:20
想问下最后题主问题解决了吗，遇到了相同的问题

作者
Author: wzkchem5 时间: 2022-2-10 00:09

huahua0621 发表于 2022-2-9 15:20
想问下最后题主问题解决了吗，遇到了相同的问题

有可能是ORCA 4.2.1的bug，我用ORCA 5.0.2测楼主的输入文件，结果只和高斯差0.003eV。
可以下载ORCA 5.0.2试试，或者等几天ORCA 5.0.3出来以后下载5.0.3

作者
Author: huahua0621 时间: 2022-2-10 22:16

wzkchem5 发表于 2022-2-10 00:09
有可能是ORCA 4.2.1的bug，我用ORCA 5.0.2测楼主的输入文件，结果只和高斯差0.003eV。
可以下载ORCA 5.0 ...

我用的是5.0.2，和singlet一样的输入只是改了MULT和坐标前的多重度，T1计算的结果甚至比S1的绝热能都高算的QD-NEVPT2

作者
Author: zjxitcc 时间: 2022-2-10 22:42

huahua0621 发表于 2022-2-10 22:16
我用的是5.0.2，和singlet一样的输入只是改了MULT和坐标前的多重度，T1计算的结果甚至比S1的绝热能都高 ...

你这问题与贴主所问完全不是同一个问题，你的问题复杂多了，而且高斯也没有NEVPT2给你对比，不能用“遇到了相同的问题”来描述。

QD-NEVPT2要讲出你的计算过程（计算流程，例如先算啥，再算啥，怎么选活性轨道，SA-CASSCF做完之后是否轨道符合自己的期望？SA-CASSCF用了几个根做平均？QD用了几个态？所用基组？），你看这细节是TDDFT的四五倍，你一个都没讲

作者
Author: wzkchem5 时间: 2022-2-11 02:21

huahua0621 发表于 2022-2-10 15:16
我用的是5.0.2，和singlet一样的输入只是改了MULT和坐标前的多重度，T1计算的结果甚至比S1的绝热能都高 ...

T1比S1能量高的分子是存在的（虽然极其罕见），算出这样的结果并不一定代表算错了。我记得听过Aspuru-Guzik讲过这种体系的计算，你可以搜一下他的文章

作者
Author: huahua0621 时间: 2022-2-16 20:47

zjxitcc 发表于 2022-2-10 22:42
你这问题与贴主所问完全不是同一个问题，你的问题复杂多了，而且高斯也没有NEVPT2给你对比，不能用“遇到 ...

DFT优化，QD-NEVPT2做单点，用DFT的轨道做的初猜，先调的CAS轨道，15个根做平均，然后算的QD-NEVPT2,吸收光谱基本和实验吻合，读了单重态的轨道去算的三重态，QD也是15个根，基组是def2-svp 对Ir用的SDD 算三重态也用的15个根，轨道也基本和之前一致，就是能量差距很大

作者
Author: zjxitcc 时间: 2022-2-16 21:06

huahua0621 发表于 2022-2-16 20:47
DFT优化，QD-NEVPT2做单点，用DFT的轨道做的初猜，先调的CAS轨道，15个根做平均，然后算的QD-NEVPT2,吸收 ...

你可以另开一贴，晒出你的分子，你的计算结果，你期望的结果（或实验值），以便详细分析、讨论

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