zhb 发表于 2021-10-28 20:07
“对于TDDFT算的激发态，只能算激发态的静电势，没法算激发态的ALIE”老师这句话的意思是我用TDDFT算的激 ...

搞清楚TDDFT激发态波函数形式和ALIE的表达式便知没法对TDDFT激发态波函数算ALIE。ALIE是对基态占据轨道来算的。

sobereva 发表于 2021-10-28 13:56
你的讨论一半正确，一半不严谨。

如果你是想考察当前分子处于激发态情况，其它分子的质子容易往这个分子 ...

“对于TDDFT算的激发态，只能算激发态的静电势，没法算激发态的ALIE”老师这句话的意思是我用TDDFT算的激发态结构的波函数不能算激发态的ALIE，为什么呢？但是我做的ALIE是用的和静电势一个含激发态波函数文件，那算出来的不应是激发态的ALIE吗？谢谢老师

你的讨论一半正确，一半不严谨。

如果你是想考察当前分子处于激发态情况，其它分子的质子容易往这个分子上哪里转移，原理上应当用激发态波函数来做波函数分析。对于TDDFT算的激发态，只能算激发态的静电势，没法算激发态的ALIE。

当前体系不大，把质子放在各个可能的质子化位点上，分别优化然后对比能量，是很容易实现的。建议你以算出来的能量作为最终判断哪里最容易接受质子的判据。如果结果和静电势、ALIE判断的相符，那么可以通过静电势、ALIE解释能量差异；如果不符，就别强行用静电势/ALIE尬讨论。