计算化学公社

标题: Laplacian键级计算结果不正常 [打印本页]

作者
Author: lyhnh 时间: 2017-5-15 23:16
标题: Laplacian键级计算结果不正常
本帖最后由 lyhnh 于 2017-5-15 23:17 编辑

Laplacian键级与键能有相当不错的相关性，想用Laplacian键级来预测键能。算了几个小分子，并用multiwfn计算Laplacian键级，却发现结果不太正常。
b3lyp,6-31(d,p)下优化甲醇，乙醛肟，氨基甲醛肟分子，并计算Laplacian键级。
甲醇的结果
# 1: 1(C ) 2(H ):  0.900845
# 2: 1(C ) 3(H ):  0.924333
# 3: 1(C ) 4(H ):  0.900902
# 4: 1(C ) 5(O ):  0.408599
# 5: 5(O ) 6(H ):  0.648910
乙醛肟的结果
# 1: 1(C ) 2(H ):  0.915290
# 2: 1(C ) 3(C ):  1.135302
# 3: 1(C ) 7(N ):  1.812945
# 4: 3(C ) 4(H ):  0.874389
# 5: 3(C ) 5(H ):  0.894044
# 6: 3(C ) 6(H ):  0.874348
# 7: 7(N ) 8(O ):  0.281471
# 8: 8(O ) 9(H ):  0.658438
氨基甲醛肟的结果
# 1: 1(C ) 2(H ):  0.925252
# 2: 1(C ) 3(N ):  1.785166
# 3: 1(C ) 6(N ):  1.118353
# 4: 3(N ) 4(O ):  0.229705
# 5: 4(O ) 5(H ):  0.660558
# 6: 6(N ) 7(H ):  0.821961
# 7: 6(N ) 8(H ):  0.831011
氧氢键键能分别为437.7，410.9，371.5 kJ/mol，与Laplacian键级的变化正好相反。而且两个肟的氮氧键键级也不正常，只有0.25左右
为什么会这样

作者
Author: sobereva 时间: 2017-5-15 23:48
“氮氧键键级也不正常”这是误解了拉普拉斯键级的本质。一方面本身这是极性键，又有一定电荷转移键特征，本身在它们之间电子密度聚集程度就低，拉普拉斯键级数值很小是正常的，没法和Mayer键级的数量级（一般接近形式键级）去对比，因为本质完全不同。

至于O-H键，也计算Mayer键级，以及考察键长进一步对比。

作者
Author: lyhnh 时间: 2017-5-15 23:57

sobereva 发表于 2017-5-15 23:48
“氮氧键键级也不正常”这是误解了拉普拉斯键级的本质。一方面本身这是极性键，又有一定电荷转移键特征，本 ...

谢谢sob老师

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