Rong Ziqin, Kitchaev Daniil, Canepa Pieremanuele, Huang Wenxuan, Ceder Gerbrand
Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA.
J Chem Phys. 2016 Aug 21;145(7):074112. doi: 10.1063/1.4960790.
The Nudged Elastic Band (NEB) is an established method for finding minimum-energy paths and energy barriers of ion migration in materials, but has been hampered in its general application by its significant computational expense when coupled with density functional theory (DFT) calculations. Typically, an NEB calculation is initialized from a linear interpolation of successive intermediate structures (also known as images) between known initial and final states. However, the linear interpolation introduces two problems: (1) slow convergence of the calculation, particularly in cases where the final path exhibits notable curvature; (2) divergence of the NEB calculations if any intermediate image comes too close to a non-diffusing species, causing instabilities in the ensuing calculation. In this work, we propose a new scheme to accelerate NEB calculations through an improved path initialization and associated energy estimation workflow. We demonstrate that for cation migration in an ionic framework, initializing the diffusion path as the minimum energy path through a static potential built upon the DFT charge density reproduces the true NEB path within a 0.2 Å deviation and yields up to a 25% improvement in typical NEB runtimes. Furthermore, we find that the locally relaxed energy barrier derived from this initialization yields a good approximation of the NEB barrier, with errors within 20 meV of the true NEB value, while reducing computational expense by up to a factor of 5. Finally, and of critical importance for the automation of migration path calculations in high-throughput studies, we find that the new approach significantly enhances the stability of the calculation by avoiding unphysical image initialization. Our algorithm promises to enable efficient calculations of diffusion pathways, resolving a long-standing obstacle to the computational screening of intercalation compounds for Li-ion and multivalent batteries.
推挤弹性带(NEB)是一种用于寻找材料中离子迁移的最小能量路径和能垒的既定方法,但当与密度泛函理论(DFT)计算结合时,其巨大的计算成本阻碍了它的广泛应用。通常,NEB计算是从已知初始态和终态之间连续中间结构(也称为影像)的线性插值开始的。然而,线性插值带来了两个问题:(1)计算收敛缓慢,特别是在最终路径呈现明显曲率的情况下;(2)如果任何中间影像过于靠近非扩散物种,NEB计算就会发散,导致后续计算不稳定。在这项工作中,我们提出了一种新方案,通过改进路径初始化和相关的能量估计工作流程来加速NEB计算。我们证明,对于离子框架中的阳离子迁移,将扩散路径初始化为基于DFT电荷密度构建的静电势中的最小能量路径,能在0.2 Å的偏差内重现真实的NEB路径,并使典型NEB运行时间最多缩短25%。此外,我们发现从这种初始化得到的局部弛豫能垒能很好地近似NEB能垒,与真实NEB值的误差在20 meV以内,同时将计算成本降低多达5倍。最后,对于高通量研究中迁移路径计算的自动化至关重要的是,我们发现新方法通过避免非物理影像初始化显著提高了计算的稳定性。我们的算法有望实现扩散路径的高效计算,解决了锂离子和多价电池插层化合物计算筛选中一个长期存在的障碍。