请教如何通过分析波函数计算bond dipole moment？

Ninilol

sob老师您好。我把图片中的体系opt之后，通过Multiwfn分析fchk，计算分子片段的偶极矩（C12-N15两个原子，分析的结果在附件的out.txt）。我想计算的是C12-N15的键偶极矩，结果显示Molecular dipole moment (a.u.):      -1.402155      0.208681      0.039040。展示在vmd里面，这个dipole moment的方向和C12-N15这根键是有一定角度的。所以想请教在C12-N15这根键上的偶极矩是否就是这个分子片段的dipole moment在键上的投影？

sobereva

如果你要表征键偶极矩，这样算不合理。因为C-N这个片段是带净电荷的，因此此片段的偶极矩是依赖于坐标原点选取的，而当前这俩原子又明显远离坐标原点。
相对合理的计算键的偶极矩的方法看下文的相关部分
Multiwfn支持的分析化学键的方法一览
http://sobereva.com/471

Ninilol

谢谢社长的回复！我根据http://sobereva.com/471去尝试计算键偶极矩。我遇到了几个问题：

（1）以下图的这个体系，图片是过渡态结构（C1-C13成键），把ts.fchk放到Multiwfn分析，设置如下
-10 Set threshold for determining one- and two-center LMO, current: 100   0 %
-9 Set the method for calculating orbital composition, current: Hirshfeld
-8 If calculating center position and dipole moment of LMOs, current: Yes
-6 Set localization method, current: Pipek-Mezey with Mulliken population
-5 If also localizing core orbitals, current: Yes
-4 If calculating and printing orbital energies, current: No
-3 If reloading newly generated .fch file, current: Yes
-2 Set criterion of convergence, current:  0.00010000
-1 Set maximum cycles, current:  80
比较奇怪的是，输出的结果里面Alpha LMOs的Two-center bond dipole moments有13C  -   1C的dipole，而Beta LMOs里面则没有这一对的dipole。这是为什么呢？

（2）有些原子对在Alpha LMOs和Beta LMOs中都有Two-center bond dipole moments，比如C13 - C5。那这一根键的dipole可以理解成Alpha和Beta的bond dipole的加和吗？
C:\Users\35247\Desktop\Snipaste_2024-11-14_23-35-04.png