计算化学公社

标题: Dipole 计算是不是必须将分子质心置于原点？ [打印本页]

作者
Author: renpj 时间: 2025-1-26 15:45
标题: Dipole 计算是不是必须将分子质心置于原点？
我在测试使用 Multiwfn 计算时发现dipole的数值跟分子的位置相关，似乎需要将分子的质心放置于原点才能得到比较正确的值。比如 N2 的偶极应该是 0 ，如果手动将 mwfn 文件中的坐标改成质心不在0 点则最后计算的 dipole 不再是 0。
我研究了一下 mwfn 文件的格式，（T. Lu, Q. Chen, 2022, DOI 10.26434/chemrxiv-2021-lt04f-v6）没有找到必须将分子质心置于原点的说明。

作者
Author: sobereva 时间: 2025-1-27 14:11
要说清楚在Multiwfn里怎么操作的，Multiwfn多个功能都能输出偶极矩，而且算法不同。并且给出截图，体现出你读的是什么数据，别把比如核电荷对偶极矩的贡献当成了偶极矩（因为还有电子对偶极矩的贡献）。

这是个中心不在原点的N2的mwfn文件（Gaussian做HF/3-21G单点计算后Multiwfn转换出的mwfn文件）
(, 下载次数 Times of downloads: 1)

Multiwfn载入之，输入
300
5
可见给出的偶极矩就是为0：
(, 下载次数 Times of downloads: 34)

作者
Author: renpj 时间: 2025-1-30 17:01

sobereva 发表于 2025-1-27 14:11
要说清楚在Multiwfn里怎么操作的，Multiwfn多个功能都能输出偶极矩，而且算法不同。并且给出截图，体现出你 ...

谢谢卢老师！Multiwfn 的计算过程跟你是一样的，我现在怀疑是mwfn 文件的生成有问题。我自己再好好检查一下

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