计算化学公社

标题: 键偶极矩的计算 [打印本页]

作者
Author: 小菜鸟求关照 时间: 2015-3-24 19:25
标题: 键偶极矩的计算
我用gauss做键偶极矩的计算，关键词是pop=nboread 为什么有的键偶极矩没有被打印出来呢？谢谢老师！

作者
Author: 小范范1989 时间: 2015-3-24 22:17
不知道有没有用，以前sob的回复，在小木虫
http://emuch.net/bbs/viewthread.php?tid=6357069

作者
Author: sobereva 时间: 2015-3-25 00:30
有多种方法计算键偶极矩
(1)Multiwfn里做ELF、LOL盆分析，计算对应于成键的盆的电多极矩（包括偶极矩）。参考手册4.17.2节或《使用Multiwfn做电子密度、ELF、静电势、密度差等函数的盆分析》（http://sobereva.com/179）
(2)进入Multiwfn主功能200，选2 Calculate atomic and bond dipole moments in Hilbert space，原理在手册3.200.2节有详细介绍
(3)用Multiwfn的模糊空间分析功能，选8 Perform integration in fuzzy overlap region for a real space functions，被积函数选择自定义函数。在计算之前，应当通过settings.ini的iuserfunc将自定义函数切换至21~23，详见手册2.8节的相应介绍。
(4)用NBO程序的dipole关键词，可以输出NBO和NLMO的键偶极矩，键偶极矩加和就是分子偶极矩。对于高斯，直接写pop=nboread并不会输出键偶极矩，必须在后面空一行再写上$NBO dipole $END才行。

每种方法结果和原理都有很大差异，毕竟键偶极矩本来就没有确切定义。(1)的原理最清楚，(2)最简单快速，但可信度逊于(1)，且不能用于有弥散函数的情况。(4)也可以。(3)不推荐。

(1)也可以用于多中心键，但多中心键必须能够以ELF/LOL盆描述出来，比如H3+的三中心键。如果是比如芳香环的pi多中心键，就不适合了。

多中心键往往是通过一个或多个分子轨道（或AdNDP轨道）描述的，这种情况可以计算轨道偶极矩，见
通过Multiwfn计算各个轨道的偶极矩
http://sobereva.com/251

作者
Author: 小菜鸟求关照 时间: 2015-3-25 07:37

小范范1989 发表于 2015-3-24 22:17
不知道有没有用，以前sob的回复，在小木虫
http://emuch.net/bbs/viewthread.php?tid=6357069

谢谢不太管用

作者
Author: 小菜鸟求关照 时间: 2015-3-25 07:41

sobereva 发表于 2015-3-25 00:30
有多种方法计算键偶极矩
(1)Multiwfn里做ELF、LOL盆分析，计算对应于成键的盆的电多极矩（包括偶极矩）。 ...

对于4，我写上了后面关键词有的键没有输出来为什么呢

作者
Author: sobereva 时间: 2015-3-25 07:46
输出文件压缩下发上来我看看

作者
Author: 小菜鸟求关照 时间: 2015-3-25 09:11

sobereva 发表于 2015-3-25 07:46
输出文件压缩下发上来我看看

老师您好我把文件发到您qq上了

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