multiwfn计算拉普拉斯键级时，键级数值会受到基组的影响吗？

机智猴

各位老师好，学生有一个关于m062xd3+不同基组会导致两个化学键的拉普拉斯键级差值的符号（正负号）不同的问题，如下：1.体系是含氮氧体系，20原子左右，存在弱相互作用（氢键）。
2.计算方法是明尼苏达泛函M062X+D3校正，由于之前是打算用来计算mayer键级，因此没用pople基组来优化结构（不加弥散不合适，加了弥散又不适合mayer键级的计算），因此分别使用了def2TZVP/def-TZVP/cc-pvtz/cc-pvdz（好基组搭配密度泛函确实会不太明智，耗时更多，但是初衷是想计算更好的结构），结果最后发现相同方法搭配不同基组计算出的某两个键的拉普拉斯键级的差值竟然出现了不同符号（化学键1的拉普拉斯键级-化学键2的拉普拉斯键级）。


3.def2-TZVP和cc-pvdz计算得出差值的符号相同，def-TZVP和cc-PVTZ计算得出的差值符号相同。
4.求各位老师帮忙分析可能出现的问题。

机智猴

忘了说细节，补充一下。
1.体系中只有C、N、O、H三种元素，不存在第四周期的重原子。
2.这两个键一个键是C-N键，一个键是N-H键，N-H键和一个氧原子形成了氢键。我看sob老师的文章中，尽管是用不同的方法\不同的基组，得到的键级的变化趋势也是相同的，所以不太懂这是为什么。

sobereva

含氢键的体系用诸如6-311G**优化完全没问题，你对基组的选用有所误解

2 DFT一般用def2-TZVP就行了。不用考虑相关一致性基组。如果本身两个键的LBO键级就差不多，不同基组下符号有异不奇怪

3 暗示cc-pVDZ级别可能太低。>=def2-TZVP档次的基组得到的结论应当是基本一样的

机智猴

sobereva 发表于 2019-7-9 03:49
含氢键的体系用诸如6-311G**优化完全没问题，你对基组的选用有所误解

2 DFT一般用def2-TZVP就行了。不用 ...

感谢sob老师的回答。确实是对基组的理解不到位，我再认真学习一下基组。