计算化学公社

标题: 色散校正在高温高压化学反应的MD模拟中是否有必要添加 [打印本页]

作者
Author:
Zt2020666    时间: 2024-8-31 23:39
标题: 色散校正在高温高压化学反应的MD模拟中是否有必要添加
本帖最后由 Zt2020666 于 2024-8-31 23:54 编辑

最近使用基于BO近似的AIMD模拟一个500原子的有机体系固相分解,温度比较高,2000到3000K,压力也得好几个GPa。使用的PBE泛函。然后体系中不可避免会存在极性分子。在考虑加不加色散校正。主要是考虑色散校正,计算速度肯定会慢一些。所以不知道在这种情况下,加与不加D3BJ色散作用到底差异大不大。不是很想加,毕竟原子体系较大,本来计算就慢,而且我只关注原子的运动轨迹。此外,之前算别的体系加了之后,从反应趋势上,没看出来啥差异,但没做过系统的测试。大家讨论一下,加与不加色散校正在高温高压反应的模拟中差异会有多大?



作者
Author:
blueyangliu    时间: 2024-9-1 09:10
色散是很弱的作用,在高温高压MD中没什么影响
作者
Author:
Zt2020666    时间: 2024-9-1 11:39
blueyangliu 发表于 2024-9-1 09:10
色散是很弱的作用,在高温高压MD中没什么影响

共轭梯度优化晶胞和原子位置,发现加了D3也没啥改善,甚至不加D3的晶格参数更接近实验值。
作者
Author:
sobereva    时间: 2024-9-2 07:16
D3校正对耗时几乎没有任何影响。D4倒是对耗时有不可忽视的影响,特别是目前的CP2K里的D4,这部分只能靠OpenMP并行

对于PBEsol这样本身优化大部分体系晶胞参数就已经很好的泛函,加色散校正非必须(甚至可能恶化结果),哪怕色散作用很显著。如果某泛函本身已经表现很好,加色散校正也不是必须的,加了有时反倒矫枉过正。另外,高压下一些原子间距离比平时研究的平衡距离可能明显更小,此时D3可能不理想,特意考虑了近距离相互作用情况的D3M校正此时在原理上更适合(但用于动力学模拟也未必就有更好结果)。

作者
Author:
Zt2020666    时间: 2024-9-4 18:39
sobereva 发表于 2024-9-2 07:16
D3校正对耗时几乎没有任何影响。D4倒是对耗时有不可忽视的影响,特别是目前的CP2K里的D4,这部分只能靠Open ...

噢噢,了解了。但是有些审稿人总会提出这些问题,说白了有些工作不得不做去堵审稿人的嘴




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3