请问，如果想要算不同CT态下的电荷，是需要先得到每个CT态的fchk文件，再计算拟合静电势电荷吗？因为我的体系比较大，如果想得到每个CT态的fchk文件，需要对每个CT态进行结构优化，这样很耗时，有其他简单点的方法吗

对D+和A-分别用Multiwfn计算拟合静电势电荷（参考Multiwfn手册4.7.1节计算CHELPG电荷的例子。文献里用的Hirshfeld电荷非常不适合此目的，对静电势重现性很差，http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml里我专门说了），并要求导出记录原子电荷和原子坐标的chg文件，然后通过文本编辑器将两个chg文件合并，再将此chg文件（对应D+A-状态）载入Multiwfn，进入主功能7，选-2 Calculate interaction energy between fragments based on atomic charges，定义两个片段分别为D+和A-部分，即得到二者间的相互作用能。