计算化学公社

 找回密码 Forget password
 注册 Register
Views: 8808|回复 Reply: 6
打印 Print 上一主题 Last thread 下一主题 Next thread

[波函数分析求助] 使用Multiwfn算Mayer键级为负数

[复制链接 Copy URL]

63

帖子

0

威望

620

eV
积分
683

Level 4 (黑子)

这是我算高斯的输入文件
%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键级,出来的结果如图所示

这是Total valences and free valences defined by Mayer的结果

请问:这里负值是否有物理意义?如有,代表什么?
谢谢各位大神!

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

2#
发表于 Post on 2021-5-7 03:16:49 | 只看该作者 Only view this author
去掉savenbos
另外,最好上传gjf文件,要么确保你发帖的时候输入文件部分中的空行和实际gjf里一致,否则很容易造成误会
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

63

帖子

0

威望

620

eV
积分
683

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-5-7 10:48:02 | 只看该作者 Only view this author
本帖最后由 shiguang 于 2021-5-7 10:50 编辑
sobereva 发表于 2021-5-7 03:16
去掉savenbos
另外,最好上传gjf文件,要么确保你发帖的时候输入文件部分中的空行和实际gjf里一致,否则很 ...

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

1043

帖子

0

威望

4188

eV
积分
5231

Level 6 (一方通行)

4#
发表于 Post on 2021-5-7 13:32:39 来自手机 | 只看该作者 Only view this author
Mayer键级是基于Mulliken布居的,正如Mulliken布居可能为负(不合理),Mayer键级也无法保证是正定的。建议用Wiberg键级,因为基于正交基的键级一定是正定的

35

帖子

0

威望

633

eV
积分
668

Level 4 (黑子)

5#
发表于 Post on 2022-3-22 15:11:21 | 只看该作者 Only view this author
granvia 发表于 2021-5-7 13:32
Mayer键级是基于Mulliken布居的,正如Mulliken布居可能为负(不合理),Mayer键级也无法保证是正定的。建议 ...

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

1043

帖子

0

威望

4188

eV
积分
5231

Level 6 (一方通行)

6#
发表于 Post on 2022-3-26 20:49:43 来自手机 | 只看该作者 Only view this author
邓苏微 发表于 2022-3-22 15:11
我想问一下,那Mayer键级如果用Multiwfn计算出来确实为负数,是否合理呢?我今天也做一个小分子的计算, ...

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

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

7#
发表于 Post on 2022-3-26 22:22:02 | 只看该作者 Only view this author
邓苏微 发表于 2022-3-22 15:11
我想问一下,那Mayer键级如果用Multiwfn计算出来确实为负数,是否合理呢?我今天也做一个小分子的计算, ...

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

另外,6#提到的基于Lowdin正交化轨道的Wiberg键级在Multiwfn的键级分析功能里直接有相应的选项。但仍然不适合带弥散函数的情况。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 23:42 , Processed in 0.233861 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list