计算化学公社

标题: 加入弥散函数会对RESP电荷造成不利影响吗？ [打印本页]

作者
Author: h840473807 时间: 2019-7-3 15:46
标题: 加入弥散函数会对RESP电荷造成不利影响吗？
出于跑分子动力学的目的，计算了如图所示的分子的RESP电荷，在B3LYP-D3/def2svp优化后分别用图示的三个等级算单点，然后用Multiwfn计算RESP电荷，发现6-311G(d,p)加上弥散后，醚氧的负电荷增加很多，都比酯基氧的还大了，def2-TZVP也有这样的趋势，但不如6-311+G(d,p)明显。这就让我很疑惑，按我的理解，对这种不带负电荷的体系计算RESP电荷，加不加弥散应该结果差不多才对，但结果却差距很大。请问加弥散会对RESP电荷计算带来什么不利影响吗，在这种情况下，我应不应该带上弥散来计算RESP电荷呢？

作者
Author: sobereva 时间: 2019-7-3 22:00
由于6-311G**还不是很大，对价电子分布描述还不算特别充分，所以弥散函数对范德华表面外静电势分布还是会有可查觉的影响的，进而影响拟合静电势电荷。

如果你对比ma-TZVP（带最低程度弥散的def2-TZVPP）和def2-TZVP，你应当会发现结果差异不会那么大。而如果你去对比比如ma-QZVP和def2-QZVP，那差异就微乎其微了。

带弥散函数只会改进拟合静电势电荷计算结果，而不会带来任何不利影响。只不过如果本身基组已经质量挺高了，那么再加弥散函数就没必要了。

作者
Author: h840473807 时间: 2019-7-5 09:11
谢谢社长，尝试了ma-TZVP，结果确实如此。但我还有一点疑问，从multiwfn计算出的静电势极小值点的数值看（如下图），三个等级计算出值非常接近，而且羰基氧上的数值明显负于醚基氧，而计算出的这两种O的RESP电荷却因基组不同而相对大小都变化了。因为我是用MD研究聚合物电解质中羰基和醚基对Li离子的配位情况，两者相对配位能力的大小对计算结果非常重要。请问6-311+G(d,p)算出来醚基氧原子电荷反比羰基氧更负的情况合理吗？
@sobereva
(, 下载次数 Times of downloads: 45)

作者
Author: wuzhiyi 时间: 2019-9-9 18:26

sobereva 发表于 2019-7-3 22:00
由于6-311G**还不是很大，对价电子分布描述还不算特别充分，所以弥散函数对范德华表面外静电势分布还是会有 ...

请问sob老师，我最近也在研究弥散对RESP的影响，感觉对我的系统我的结果也和您说的类似

6-311+G**和6-311G**有一定差别

但ma-TZVP和def2-TZVP差别不大

但想问一下您这个结论是如何得到的呢？我自己查发现做这类benchmark的RESP文献里没看看，在[原子电荷计算方法的对比]里也没有对于弥散对于RESP的影响，想问一下sob老师有啥推荐的文献讨论了弥撒对RESP的影响。

作者
Author: sobereva 时间: 2019-9-9 18:41

wuzhiyi 发表于 2019-9-9 18:26
请问sob老师，我最近也在研究弥散对RESP的影响，感觉对我的系统我的结果也和您说的类似

6-311+G**和6- ...

之前我做过专门的基组对静电势影响的测试，测试了一批基组，没有发表而已

作者
Author: wuzhiyi 时间: 2019-9-9 19:01

sobereva 发表于 2019-9-9 18:41
之前我做过专门的基组对静电势影响的测试，测试了一批基组，没有发表而已

谢谢sob老师，不知道老师有计划公开这些数据吗？比如在预印本之类的地方？

作者
Author: sobereva 时间: 2019-9-10 06:41

wuzhiyi 发表于 2019-9-9 19:01
谢谢sob老师，不知道老师有计划公开这些数据吗？比如在预印本之类的地方？

没，就是用几个简单体系做的测试，还不足以构成一篇文章
以后是有打算写成文章，不仅讨论对拟合静电势电荷的影响，还讨论对Multiwfn做的静电势的定量分子表面分析的那些量的影响，筛出几个最划算的算静电势相关问题的基组

作者
Author: wuzhiyi 时间: 2019-9-10 20:35

sobereva 发表于 2019-9-10 06:41
没，就是用几个简单体系做的测试，还不足以构成一篇文章
以后是有打算写成文章，不仅讨论对拟合静电势电 ...

謝謝sob老師，希望可以看到文章儘早面世。

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