计算化学公社

标题: 利用Multiwfn计算BCP性质（AIM）比较化学键出现问题怎么分析 [打印本页]

作者
Author: qdluck 时间: 2024-9-2 09:42
标题: 利用Multiwfn计算BCP性质（AIM）比较化学键出现问题怎么分析
本帖最后由 qdluck 于 2024-9-2 09:52 编辑

（1）目标对象：化学键A（C-O）和化学键B（C-O），注：以上A、B两键是分别属于两个类似分子相同位置的共价键；（2）问题描述：利用Multiwfn分别做电子密度拓扑分析，获得各自BCP位置及性质，希望判断化学键强弱，根据AIM分析，电子密度、势能密度、能量密度、 Sign(lambda2)*rho等均支持实验结果，即A键＞B键，但根据IGMH手册中Delta-g (under Hirshfeld partition)数值，却出现相反结果，即B键＞A键，前者为0.3016，后者为0.2953，结果相差很小，请问这种情况，出现这种情况的原因，以及是否可以忽略Delta-g的差异？

作者
Author: sobereva 时间: 2024-9-2 10:35
Delta-g不是严格的判据。另外，此处Delta-g指代不明，没有唯一计算方式，诸如可以算BCP位置的值、delta_g_inter极大点位置的值等，应当说清楚。

作者
Author: qdluck 时间: 2024-9-2 10:46

sobereva 发表于 2024-9-2 10:35
Delta-g不是严格的判据。另外，此处Delta-g指代不明，没有唯一计算方式，诸如可以算BCP位置的值、delta_g_i ...

好的，谢谢老师，我是根据《Tutorial of using Multiwfn to perform IGMH analysis》手册中“1.4 Evaluating δg at bond critical points.”计算的，本文中的Delta-g (under Hirshfeld partition)是用Multiwfn做电子密度拓扑分析获得的BCP位置及性质中得到的，其他Delta-g的计算方式我还不太了解，再次谢谢老师

作者
Author: qdluck 时间: 2024-9-2 15:44
本帖最后由 qdluck 于 2024-9-2 15:49 编辑

sobereva 发表于 2024-9-2 10:35
Delta-g不是严格的判据。另外，此处Delta-g指代不明，没有唯一计算方式，诸如可以算BCP位置的值、delta_g_i ...

补充测试：老师，分别利用delta_g（BCP位置）和delta_g_inter（max极大值）计算得到了delta_g，结果出现了相反的情况，后者基于极大值得到的delta_g（化学键A为0.3711，化学键B为 0.3781，即A＞B）更符合实验结果，那前者基于BCP得到的delta_g（化学键A为0.2953，化学键B为0.3016，即A＜B）是否可忽略？这样出现相反情况的原因是什么呢？

作者
Author: sobereva 时间: 2024-9-3 11:31

qdluck 发表于 2024-9-2 15:44
补充测试：老师，分别利用delta_g（BCP位置）和delta_g_inter（max极大值）计算得到了delta_g，结果出现 ...

后者本来物理意义就更强，既然已经算过了，就没必要再考虑BCP处的delta-g

作者
Author: qdluck 时间: 2024-9-3 12:13
好的，收到，感谢老师指点，谢谢啦

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