计算化学公社

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

[波函数分析求助] Mayer键级与模糊键级计算结果不同

[复制链接 Copy URL]

168

帖子

0

威望

644

eV
积分
812

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
老师,我想问一下,我做的是过渡金属催化有机物反应,阴离子体系,优化用的是b3lyp/6-31+G(d),过渡金属是SDD。在计算键级时,发现Mayer键级与模糊键级计算出来的结果正好相反。
Mayer:M-C1为0.294,M-C2为0.334,M-C1<M-C2;模糊键级:M-C1为1.05,M-C2为0.946,M-C1>M-C2。


由于Mayer键级不能使用弥散函数,又在M06-2X/def2-SVP/SDD水平上优化了一次,得到的结果依然不一样。
Mayer:M-C1为0.607,M-C2为0.75,M-C1<M-C2;模糊键级:M-C1为1.121,M-C2为0.997,M-C1>M-C2。
想问一下老师,我该取哪个键级更正确一点?

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

2#
发表于 Post on 2022-3-10 22:06:33 | 只看该作者 Only view this author
完全没必要重新优化。而且M06-2X非常不适合算过渡金属配合物这是基本常识,而且def2-SVP才2-zeta还不带弥散函数明显不适合表现带显著负电荷体系的电子结构。用PBE0/def2-TZVP算单点多好
好好看
简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html
谈谈量子化学中基组的选择
http://sobereva.com/336http://bbs.keinsci.com/thread-3545-1-1.html

在这个级别重算后再说
北京科音自然科学研究中心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

168

帖子

0

威望

644

eV
积分
812

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2022-3-11 10:06:54 | 只看该作者 Only view this author
好的,谢谢老师

168

帖子

0

威望

644

eV
积分
812

Level 4 (黑子)

4#
 楼主 Author| 发表于 Post on 2022-3-13 14:15:26 | 只看该作者 Only view this author
sobereva 发表于 2022-3-10 22:06
完全没必要重新优化。而且M06-2X非常不适合算过渡金属配合物这是基本常识,而且def2-SVP才2-zeta还不带弥散 ...


老师,我又用PBE0/def2-TZVP/SDD计算了一下单点,得到的结果是
Mayer:M-C1为0.34,M-C2为0.60,M-C1<M-C2;模糊键级:M-C1为1.05,M-C2为0.95,M-C1>M-C2
想问一下老师,取哪个键级更好呢

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

5#
发表于 Post on 2022-3-14 02:10:50 | 只看该作者 Only view this author
化小白 发表于 2022-3-13 14:15
老师,我又用PBE0/def2-TZVP/SDD计算了一下单点,得到的结果是
Mayer:M-C1为0.34,M-C2为0.60,M-C1M ...

算一下BCP的电子密度,如果M-C2比M-C1大,拿Mayer键级说事
北京科音自然科学研究中心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

168

帖子

0

威望

644

eV
积分
812

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2022-3-14 11:37:22 | 只看该作者 Only view this author
老师,我算出来的结果是
M-C1为Density of all electrons:  0.1043503553E+00,M-C2为Density of all electrons:  0.1056431862E+00,M-C1<M-C2
老师,确实是M-C2大一点,那我就用Mayer键级分析,谢谢老师

168

帖子

0

威望

644

eV
积分
812

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2022-3-14 11:44:56 | 只看该作者 Only view this author
sobereva 发表于 2022-3-14 02:10
算一下BCP的电子密度,如果M-C2比M-C1大,拿Mayer键级说事

老师,我想问一下,如果M-C2键级>M-C1键级,可以说明在断键的时候,M-C1比M-C2先断键吗?并且算出来的断M-C2键的能垒确实比断M-C1键的能垒高

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

8#
发表于 Post on 2022-3-15 02:22:07 | 只看该作者 Only view this author
化小白 发表于 2022-3-14 11:44
老师,我想问一下,如果M-C2键级>M-C1键级,可以说明在断键的时候,M-C1比M-C2先断键吗?并且算出来的断M ...

可以这么解释
北京科音自然科学研究中心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

168

帖子

0

威望

644

eV
积分
812

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2022-3-15 16:18:58 | 只看该作者 Only view this author

好的,谢谢老师

168

帖子

0

威望

644

eV
积分
812

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2022-3-15 16:23:20 | 只看该作者 Only view this author

老师,计算Mayer键级所用的泛函和基组可以和计算单点时所用的泛函和基组不一样吗?我计算单点用的是M06/ma-def2-tzpv /SDD,而且已经全部算完了

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

11#
发表于 Post on 2022-3-16 01:44:37 | 只看该作者 Only view this author
化小白 发表于 2022-3-15 16:23
老师,计算Mayer键级所用的泛函和基组可以和计算单点时所用的泛函和基组不一样吗?我计算单点用的是M06/m ...

可以
那叫ma-def2-TZVP
北京科音自然科学研究中心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

168

帖子

0

威望

644

eV
积分
812

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2022-3-17 17:40:12 | 只看该作者 Only view this author
sobereva 发表于 2022-3-16 01:44
可以
那叫ma-def2-TZVP

好的,谢谢老师提醒,老师,我想问一下,过渡金属M-C和M-O之间的键级也是可通过Mayer键级来比较吧?不是说只有M-C1和M-C2之间才可以相互比较的吧?
老师,可以用RESP拟合静电势电荷计算得到过渡金属M为0.597,C1为-0.594,C2为-0.513,来说明M-C1更难断吗

6万

帖子

99

威望

6万

eV
积分
125145

管理员

公社社长

13#
发表于 Post on 2022-3-18 06:55:18 | 只看该作者 Only view this author
化小白 发表于 2022-3-17 17:40
好的,谢谢老师提醒,老师,我想问一下,过渡金属M-C和M-O之间的键级也是可通过Mayer键级来比较吧?不是 ...

可以

原子电荷和键的强弱没直接关系
北京科音自然科学研究中心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

168

帖子

0

威望

644

eV
积分
812

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2022-3-18 11:12:07 | 只看该作者 Only view this author
sobereva 发表于 2022-3-18 06:55
可以

原子电荷和键的强弱没直接关系

好的,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-22 10:02 , Processed in 0.188849 second(s), 19 queries , Gzip On.

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