计算化学公社

标题: 求助：关于不同体系的NBO轨道作用的疑问 [打印本页]

作者
Author: suosuosky 时间: 2019-10-8 02:23
标题: 求助：关于不同体系的NBO轨道作用的疑问
请问sob老师，我在计算两个体系 HCN和Ph-CN分别与另一个分子中的P作用，两个体系的作用能分别为 17和 27kcal/mol。

在进行NBO分析（NBO6）时发现，两者的轨道作用能差别很大，其中两个体系中较大的轨道作用能如下：

HCN:
  18. LP ( 1) N  5          30. BD*( 2) P  1- O  3    20.85 0.91 0.123
18. LP ( 1) N  5          32. BD*( 2) P  1- O  4    20.82 0.91 0.123

Ph-CN:
19. LP ( 1) N 12          53. LV ( 1) P 16       169.20 0.49 0.258
19. LP ( 1) N 12          74. BD*( 1) P 16- O 18    3.98 0.93 0.054
19. LP ( 1) N 12          75. BD*( 1) P 16- O 19    3.98 0.93 0.054

请问老师：
1-为什么HCN这个体系中NBO输出文件直接没有找到LV(P)这种轨道？
2-两个体系NBO差别这么大，可能是什么原因造成的？是NBO分析结果不准确吗？
非常感谢老师。

附件是NBO6进行NBO分析的输出文件。

作者
Author: sobereva 时间: 2019-10-8 07:19
看不到你的结构没法判断
对于弱相互作用范畴，靠NBO判断基本没有意义

作者
Author: suosuosky 时间: 2019-10-8 15:57

sobereva 发表于 2019-10-8 07:19
看不到你的结构没法判断
对于弱相互作用范畴，靠NBO判断基本没有意义

非常感谢老师的回复。
我在计算两个体系的时候，发现NBO分析时这两个体系差别较大，不清楚是什么原因造成的。所以想咨询一下老师，大概是什么原因，为什么HCN体系没有LV(P)的轨道作用。

附件是两个体系的构型，因为优化后的构型log文件比较大上传不了，上传的是优化后计算单点能的输出文件。麻烦老师看一下。谢谢！

作者
Author: sobereva 时间: 2019-10-10 00:46
我在你的另一个帖子里已经说了你这样用B2PLYP波函数是错的，先把这个问题修正了再说

作者
Author: suosuosky 时间: 2019-10-16 00:20

sobereva 发表于 2019-10-10 00:46
我在你的另一个帖子里已经说了你这样用B2PLYP波函数是错的，先把这个问题修正了再说

谢谢老师的回复。

您在我之前的帖子里回复我，如果使用B3LYP-D3/6-311+g**水平上获得的波函数进行ELF分析时可以的。我对体系进行NBO分析也是在此水平上获得的波函数进行分析的，B2PLYP-D3方法只是用来获得结合能和作用能信息用的。

请问老师，两个体系的具有较大NBO轨道作用的两个轨道之间相差较大可能是什么原因造成的。我应该在什么角度分析比较合适？谢谢老师！

作者
Author: sobereva 时间: 2019-10-16 00:27

suosuosky 发表于 2019-10-16 00:20
谢谢老师的回复。

您在我之前的帖子里回复我，如果使用B3LYP-D3/6-311+g**水平上获得的波函数进行ELF ...

NBO轨道搜索方式存在极大的任意性，其算法很不优雅，有时几何或电子结构特征稍微变化一点就导致判断出的NBO轨道显著不同，所以不要太把NBO的E2当回事，特别是模棱两可的时候
到底输出的有没有意义，用Multiwfn看一下轨道图形进行判断。不要NBO输出什么，就想办法对其予以解释，很可能输出的数据本来就是缺乏化学意义的。
对于弱相互范畴的作用，我是极力不建议拿NBO分析说事的

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