计算化学公社

标题: 请问文中这种不同CT态的Ecoul值是怎么算的 [打印本页]

作者
Author:
zhangxx    时间: 2023-11-16 09:20
标题: 请问文中这种不同CT态的Ecoul值是怎么算的
请问文中这种不同CT态的Ecoul值是怎么算的。 这是这篇文章的DOI号,DOI: 10.1021/acs.jpcc.6b09927   



作者
Author:
sobereva    时间: 2023-11-17 02:46
对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-部分,即得到二者间的相互作用能。
作者
Author:
zhangxx    时间: 2023-11-18 09:45
请问,如果想要算不同CT态下的电荷,是需要先得到每个CT态的fchk文件,再计算拟合静电势电荷吗?因为我的体系比较大,如果想得到每个CT态的fchk文件,需要对每个CT态进行结构优化,这样很耗时,有其他简单点的方法吗




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3