二聚体设定体系电荷为0，结构优化后一个分子电荷为0.2另一分子电荷为-0.2是否正常

Yaohaohao · 发表于 Post on 2023-7-13 17:26:50

本帖最后由 Yaohaohao 于 2023-7-13 17:26 编辑

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

1.zip (5.25 MB, 下载次数 Times of downloads: 4)

wzkchem5 · 发表于 Post on 2023-7-13 18:01:32

正常的，两个彼此靠近的、不太小的分子，有0.2e量级的电荷转移还是比较常见的

Yaohaohao · 发表于 Post on 2023-7-14 08:38:48

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

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

zjxitcc · 发表于 Post on 2023-7-14 10:54:21

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

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

Yaohaohao · 发表于 Post on 2023-7-14 11:00:39

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

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

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

[新手求助] 二聚体设定体系电荷为0，结构优化后一个分子电荷为0.2另一分子电荷为-0.2是否正常