计算化学公社

标题: 稀有气体二聚体对积分格点较敏感 [打印本页]

作者
Author: sobereva 时间: 2017-2-14 17:40
标题: 稀有气体二聚体对积分格点较敏感
IJQC刚出了篇文章，测试不同带色散校正的泛函和CCSD(T)计算稀有气体二聚体势能面
(, 下载次数 Times of downloads: 37)
结论是B3LYP-D3(zero)相对来说最理想，表现差强人意。而通常认为普通泛函中算弱相互作用最好之一的M06-2X-D3对此类体系则不甚理想。

最值得一说的是测试发现对于这类体系，int=ultrafine对于M06-2X-D3也不够用，即图中绿线有波折，而且对于越轻的稀有气体波折越轻，fine时则是烂得一塌糊涂，根本没法用。图中都是带-D3的，F、UF、SF代表Fine、Ultrafine、Superfine。
(, 下载次数 Times of downloads: 31)

一般认为B3LYP这样非meta泛函用Fine足够了，不过对于Ne2，fine的时候波动也很明显，对于Ar、Kr的情况，fine的势能曲线和ultrafine有显著偏差。所以对于算稀有气体二聚体时候要注意积分格点的选取。

此例也反映出，靠色散校正，是没法明显削弱积分格点精度不足情况下PES的皱纹特征的。

作者
Author: helpme 时间: 2017-2-14 19:40
谢谢sob推荐有意思的好文。
考察的R都很大，作用力又很弱，是不是就会这样？比如考虑 CH4 。。。CH4，是不是也会这样表现？

作者
Author: 我本是个娃娃 时间: 2017-2-14 20:34
测试的是零阻尼情况下的计算结果，不知道BJ阻尼结果如何。另外，都是气体分子，没有其他分子，如H2-H2这样的体系的测试

作者
Author: sobereva 时间: 2017-2-14 21:49

我本是个娃娃发表于 2017-2-14 20:34
测试的是零阻尼情况下的计算结果，不知道BJ阻尼结果如何。另外，都是气体分子，没有其他分子，如H2-H2这样 ...

一样

作者
Author: sobereva 时间: 2017-2-14 21:53

helpme 发表于 2017-2-14 19:40
谢谢sob推荐有意思的好文。
考察的R都很大，作用力又很弱，是不是就会这样？比如考虑 CH4 。。。CH4，是不 ...

也不排除是因为原子数太少导致此现象比较明显，因为积分格点是从各个原子中心发散分布的，原子数越多格点数也会越多，因此甲烷二聚体估计比其等电子体Ne2的这个问题要弱化

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