计算化学公社

标题: 请教如何通过分析波函数计算bond dipole moment？ [打印本页]

作者
Author: Ninilol 时间: 2024-11-5 23:10
标题: 请教如何通过分析波函数计算bond dipole moment？
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在键上的投影？

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

作者
Author: Ninilol 时间: 2024-11-14 23:49
本帖最后由 Ninilol 于 2024-11-14 23:50 编辑

谢谢社长的回复！我根据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

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