请问文中这种不同CT态的Ecoul值是怎么算的

zhangxx · 发表于 Post on 2023-11-16 09:20:51

请问文中这种不同CT态的Ecoul值是怎么算的。这是这篇文章的DOI号，DOI: 10.1021/acs.jpcc.6b09927

sobereva · 发表于 Post on 2023-11-17 02:46:06

对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-部分，即得到二者间的相互作用能。

zhangxx · 发表于 Post on 2023-11-18 09:45:04

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[波函数分析求助] 请问文中这种不同CT态的Ecoul值是怎么算的