在orca中如何在已有分子激发态能量的情况下切换泛函计算soc？

adreamy

比如用B3LYP/def2svp计算了水分子的激发态能量之后，如何在B3LYP/def2svp得到的激发态能量的基础上切换泛函（如M062x）去计算soc？求教。希望大佬提供一个满足要求输入文件。谢谢

wzkchem5

首先，一般不建议TDDFT激发能和SOC用不同泛函。因为适合算SOC的泛函一般没有理由不适合激发能，反之亦然，如果发现一个算得准另外一个就算得不准，那一般是因为其他方面考虑得不准，例如溶剂效应、H-T效应、非谐性效应等，而不是说明激发能和SOC需要用不同泛函。
其次，如果一定要用不同泛函算SOC和激发能的话，那么直接按正常的输入文件写法，做一个M06-2X的SOC计算就行了，这个计算也会输出激发能，但你不报这个激发能而报B3LYP的激发能就行了。没有任何办法根本不得到M06-2X的激发能就得到M06-2X的SOC，只能是先得到M06-2X的激发能再弃而不用。但是这样做有两个前提条件：（1）B3LYP和M06-2X的激发态成分必须定性符合，不能不检查激发态成分就假设B3LYP的S1态对应M06-2X的S1态；（2）正如我所说，用不同泛函算激发能和SOC几乎肯定是不合理的做法，如果没有十足的理由，即使能这么算，也不要这么算。

adreamy

wzkchem5 发表于 2023-5-22 21:06
首先，一般不建议TDDFT激发能和SOC用不同泛函。因为适合算SOC的泛函一般没有理由不适合激发能，反之亦然， ...

好的，谢谢老师。

zjxitcc

如果你觉得你的建模、计算步骤无误，也恰当考虑了环境的影响，计算结果仍然不满足期望/与实验结果相差过多，那你应该用更严谨的理论方法（如EOM-CCSD或CASSCF-NEVPT2）算激发态和SOC，而非换泛函

adreamy

zjxitcc 发表于 2023-5-23 14:53
如果你觉得你的建模、计算步骤无误，也恰当考虑了环境的影响，计算结果仍然不满足期望/与实验结果相差过多 ...

好的，谢谢老师。