计算化学公社

标题: 使用Multiwfn算Mayer键级为负数 [打印本页]

作者
Author: shiguang 时间: 2021-5-6 18:01
标题: 使用Multiwfn算Mayer键级为负数
这是我算高斯的输入文件
%rwf=2.rwf
%NoSave
%chk=2.chk
%nprocshared=64
%mem=5GB
#p opt freq genecp em=gd3bj pbe1pbe pop=(nboread,savenbos)

2

1 1
[原子坐标]

Pd 0
SDD
****
H C O N 0
def2TZVP
****

Pd 0
SDD

$nbo bndidx $end
我用Multiwfn想计算N原子与Pd之间的Mayer键级，出来的结果如图所示
(, 下载次数 Times of downloads: 33)
这是Total valences and free valences defined by Mayer的结果
(, 下载次数 Times of downloads: 35)
请问：这里负值是否有物理意义？如有，代表什么？
谢谢各位大神！

作者
Author: sobereva 时间: 2021-5-7 03:16
去掉savenbos
另外，最好上传gjf文件，要么确保你发帖的时候输入文件部分中的空行和实际gjf里一致，否则很容易造成误会

作者
Author: shiguang 时间: 2021-5-7 10:48
本帖最后由 shiguang 于 2021-5-7 10:50 编辑

sobereva 发表于 2021-5-7 03:16
去掉savenbos
另外，最好上传gjf文件，要么确保你发帖的时候输入文件部分中的空行和实际gjf里一致，否则很 ...

谢谢sob老师！这个我是直接从gjf文件里粘贴的，没想到在贴子里会有两个空行

作者
Author: granvia 时间: 2021-5-7 13:32
Mayer键级是基于Mulliken布居的，正如Mulliken布居可能为负（不合理），Mayer键级也无法保证是正定的。建议用Wiberg键级，因为基于正交基的键级一定是正定的

作者
Author: 邓苏微 时间: 2022-3-22 15:11

granvia 发表于 2021-5-7 13:32
Mayer键级是基于Mulliken布居的，正如Mulliken布居可能为负（不合理），Mayer键级也无法保证是正定的。建议 ...

我想问一下，那Mayer键级如果用Multiwfn计算出来确实为负数，是否合理呢？我今天也做一个小分子的计算，确实为负数，只不过很小，-0.03

作者
Author: granvia 时间: 2022-3-26 20:49

邓苏微发表于 2022-3-22 15:11
我想问一下，那Mayer键级如果用Multiwfn计算出来确实为负数，是否合理呢？我今天也做一个小分子的计算， ...

这是Mayer键级定义本身所致，与使用什么软件计算无关（你也可以用Gaussian的IOP选项来计算Mayer键级，结果应该与Multiwfn一致）。建议使用Wiberg键级，用NBO程序可以算（也可以基于Lowdin正交化原子轨道计算，但得自己编个小脚本）

作者
Author: sobereva 时间: 2022-3-26 22:22

邓苏微发表于 2022-3-22 15:11
我想问一下，那Mayer键级如果用Multiwfn计算出来确实为负数，是否合理呢？我今天也做一个小分子的计算， ...

轻微负值直接当成0
明显负值时注意看你的基组带没带弥散函数，带的话Mayer键级没法用，可以用模糊键级（通常键的极性越弱和Mayer键级相符程度越高），或者去掉弥散函数再计算。

另外，6#提到的基于Lowdin正交化轨道的Wiberg键级在Multiwfn的键级分析功能里直接有相应的选项。但仍然不适合带弥散函数的情况。

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