计算化学公社

标题: 二聚体设定体系电荷为0，结构优化后一个分子电荷为0.2另一分子电荷为-0.2是否正常 [打印本页]

作者
Author: Yaohaohao 时间: 2023-7-13 17:26
标题: 二聚体设定体系电荷为0，结构优化后一个分子电荷为0.2另一分子电荷为-0.2是否正常
本帖最后由 Yaohaohao 于 2023-7-13 17:26 编辑

各位老师好，我用高斯将一个分子进行结构优化（文件名为1），优化结束后在gview里复制分子建立了二聚体体系（文件名为dimer)继续结构优化，结束后用Multiwfn计算了Hirshfeld charge，单分子和二聚体体系电荷和自旋多重度都设置为0 3。但是我发现将dimer体系中的两个分子的电荷分别加和，一个分子是0.1985，另一个分子是-0.1985，一个分子带正电另一个分子带负电体系总和为0。请问各位老师这种情况是否正常？我设想的是dimer体系中的2个分子应该都为中性分子，电荷近似为0。请各位老师不吝赐教，感谢。 (, 下载次数 Times of downloads: 4)

作者
Author: wzkchem5 时间: 2023-7-13 18:01
正常的，两个彼此靠近的、不太小的分子，有0.2e量级的电荷转移还是比较常见的

作者
Author: Yaohaohao 时间: 2023-7-14 08:38

wzkchem5 发表于 2023-7-13 18:01
正常的，两个彼此靠近的、不太小的分子，有0.2e量级的电荷转移还是比较常见的

请问老师有没有什么计算办法，能让两个分子间不发生电荷转移呢，就是在这个二聚体体系里确保两个分子各自的电荷都为0。谢谢

作者
Author: zjxitcc 时间: 2023-7-14 10:54

Yaohaohao 发表于 2023-7-14 08:38
请问老师有没有什么计算办法，能让两个分子间不发生电荷转移呢，就是在这个二聚体体系里确保两个分子各自 ...

有，用constrained DFT方法。NWChem，Q-Chem和PSI4都支持。不过要注意这算出的是假想状态，不是真实情况。
《谈谈约束性DFT(CDFT)》http://sobereva.com/271

作者
Author: Yaohaohao 时间: 2023-7-14 11:00

zjxitcc 发表于 2023-7-14 10:54
有，用constrained DFT方法。NWChem，Q-Chem和PSI4都支持。不过要注意这算出的是假想状态，不是真实情况 ...

好的，谢谢老师，我学习一下

欢迎光临计算化学公社 (http://bbs.keinsci.com/)