计算化学公社

标题: 关于Mayer键级判断原子是否成键的问题 [打印本页]

作者
Author: Are 时间: 2024-6-4 21:02
标题: 关于Mayer键级判断原子是否成键的问题
本人想通过Mayer键级判断图（已隐藏C上的H）中K与O（选中的四个原子为K与3个O）是否成键，利用http://sobereva.com/414中的方法，发现并未输出K（82号原子）与O的键级，猜测键级应该小于0.05 (, 下载次数 Times of downloads: 15) (, 下载次数 Times of downloads: 13)
目前实验上已经证明有这种结构的存在 (, 下载次数 Times of downloads: 16) ，我在PBE0泛函下计算，基组选择6-31G（d,p）和LANL2DZ，优化后的键长和实验差0.03埃
本来的想法是K-O键为离子键，不能用Mayer键级判断，但是看到老师分析NaCl是可以拿到0.8的键级，所以对一开始的想法产生怀疑

请问老师我这里可以判定K-O未成键吗？还是有别的方法可以说明这里是否成键

作者
Author: sobereva 时间: 2024-6-5 02:31
算一下K+的结合能，看数量级就知道算不算
让Multiwfn导出bndmat.txt键级矩阵，K相关的键级是多少一看便知，不用猜

作者
Author: Are 时间: 2024-6-5 09:09

sobereva 发表于 2024-6-5 02:31
算一下K+的结合能，看数量级就知道算不算
让Multiwfn导出bndmat.txt键级矩阵，K相关的键级是多少一看便知 ...

感谢老师回复。老师，您说的算K+的结合能指的是把K+和团簇分成两个片段，然后用gaussian里的counterpoise关键词计算吗？

作者
Author: sobereva 时间: 2024-6-6 00:40

Are 发表于 2024-6-5 09:09
感谢老师回复。老师，您说的算K+的结合能指的是把K+和团簇分成两个片段，然后用gaussian里的counterpoise ...

是
counterpoise不是必须的，自己手动求差也可以

作者
Author: Are 时间: 2024-6-11 16:46
本帖最后由 Are 于 2024-6-11 16:47 编辑

sobereva 发表于 2024-6-6 00:40
是
counterpoise不是必须的，自己手动求差也可以

非常感谢老师的回复
老师，我计算出Li+与[Ti2V4O5(OCH3)14]-的相互作用能为-169.81 kcal/mol， K+与[Ti2V4O5(OCH3)14]-的相互作用能为-112.57 kcal/mole，NaCl中Na+与Cl-的相互作用能为-132.32 kcal/mol，在同一数量级。所以，请问老师：
1：这种情况我们是否可以认为Li和K都与多酸中O成键
2：计算得到的与多酸的相互作用能Li+大于K+的，但实验中本体电解表示K[Ti2V4O5(OCH3)14]更稳定，且用软硬酸碱理论解释稳定性，认为K+与多酸阴离子相互作用更强，所以想请问老师造成计算的相互作用能大小顺序和实验上不一致的原因是什么?当前情况下有无更好的方法去表征金属阳离子和多酸阴离子之间的相互作用能呢？或者有什么方法去解释本体电解中KPOVs比LiPOVs更稳定的原因吗？实验文章：DOI: 10.1039/d3ta06432h
(, 下载次数 Times of downloads: 18) (, 下载次数 Times of downloads: 17)

作者
Author: sobereva 时间: 2024-6-12 02:50

Are 发表于 2024-6-11 16:46
非常感谢老师的回复
老师，我计算出Li+与[Ti2V4O5(OCH3)14]-的相互作用能为-169.81 kcal/mol， K+与[Ti2 ...

1 可以
2 Li+的正电荷密度更高，极化作用更强，有相互作用更强的道理。可以做个sobEDA能量分解对差异做更深入讨论
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

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