求助ORCA基于TD-双杂化泛函计算激发态与基态单点能计算E0,0差值异常大

曹丽佳

各位老师，请问：对20个左右原子且无重金属的有机小分子，先用Gaussian基于TD-PBE0-D3(BJ)/def2-TZVP进行基态和激发态结构的优化（第一张图片），后用ORCA使用TD-双杂化TD-DSD-PBEP86-D3(BJ)/def2-TZVP计算激发态的单点能，减去双杂化泛函DSD-PBEP86-D3(BJ)/def2-TZVP计算基态的单点能，得到E_0,0绝热跃迁能。
①发现ORCA计算出来的基态单点能（第二张图片）与激发态的单点能（第三张图片）相差非常大，即（激发态：-531.945）-（基态，-533.122）这是合理的吗（直觉不是很合理），是计算思路哪里出现了问题？
②对比使用Gaussian基于DSD-PBEP86-D3(BJ)/def2-TZVP计算基态的单点能得出的结果是-533.1008383，与ORCA计算出的基态单点能（-533.122）是差不多，说明基态分子结构坐标没有问题，是激发态的分子坐标的问题？
非常感谢各位老师的指正!