中性分子添加背景电荷后如何正确计算偶极矩呢

pikachuupup · 发表于 Post on 2024-3-30 10:04:59

各位老师好，我使用ORCA5.0.3进行偶极矩计算时，使用的方法为B3LYP/def2-TZVPD，计算了单个分子比如氢气分子在背景电荷（比如电量为1的正电荷）下存在的总偶极矩，但是我发现此时偶极矩的数值依赖于坐标的选择，当把体系与背景电荷同时简单的平移时，体系的偶极矩就会发生明显的变化。那么如何才能正确计算中性分子添加背景电荷后的偶极矩呢？或者说利用ORCA或Gaussian可以通过一些关键词的限定来实现吗？谢谢各位老师！

pikachuupup · 发表于 Post on 2024-3-30 10:25:02

另外学生换用Gaussian16计算时发现，关键词中无论是否添加nosymm，无论如何改变体系的坐标都不影响中性分子在背景电荷下总偶极矩的数值，为何会与ORCA有这方面的差异呢？还请诸位老师赐教

sobereva · 发表于 Post on 2024-3-30 16:40:26

如果背景电荷是固定不动的，中性分子平移时由于被背景电荷极化情况的不同，分子的偶极矩肯定会变化
如果背景电荷和中性分子一起平移，原理上中性分子的偶极矩不可能变化。如果你在Gaussian里就是这么操作的，那么Gaussian给你的结果没问题。如果ORCA里反常，可能是其计算偶极矩的规则有什么特殊性，此时你可以用Multiwfn载入ORCA的molden文件算偶极矩，计算方式和标准定义一样，进主功能300的子功能5就能得到（建议先进入主功能0看当前体系里的原子信息，确认背景电荷并没被当做一个原子记录进去）。

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

pikachuupup · 发表于 Post on 2024-3-30 17:02:01

sobereva 发表于 2024-3-30 16:40
如果背景电荷是固定不动的，中性分子平移时由于被背景电荷极化情况的不同，分子的偶极矩肯定会变化
如果背 ...

谢谢sob老师，这次疏忽了，一定注意~

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

[ORCA] 中性分子添加背景电荷后如何正确计算偶极矩呢