Please use this identifier to cite or link to this item: https://doi.org/10.1063/1.4960790
Title: An efficient algorithm for finding the minimum energy path for cation migration in ionic materials
Authors: Rong, Z
Kitchaev, D
Canepa, P 
Huang, W
Ceder, G
Keywords: Computation theory
Density functional theory
Diffusion barriers
Energy barriers
Intercalation compounds
Interpolation
Ions
Lithium-ion batteries
Positive ions
Reaction kinetics
Computational expense
General applications
High throughput studies
Intermediate image
Intermediate structures
Linear Interpolation
Minimum energy paths
Nudged elastic band
Computational efficiency
Issue Date: 2016
Publisher: American Institute of Physics Inc.
Citation: Rong, Z, Kitchaev, D, Canepa, P, Huang, W, Ceder, G (2016). An efficient algorithm for finding the minimum energy path for cation migration in ionic materials. Journal of Chemical Physics 145 (7) : 74112. ScholarBank@NUS Repository. https://doi.org/10.1063/1.4960790
Rights: Attribution 4.0 International
Abstract: 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. © 2016 Author(s).
Source Title: Journal of Chemical Physics
URI: https://scholarbank.nus.edu.sg/handle/10635/179608
ISSN: 0021-9606
DOI: 10.1063/1.4960790
Rights: Attribution 4.0 International
Appears in Collections:Elements
Staff Publications

Show full item record
Files in This Item:
File Description SizeFormatAccess SettingsVersion 
10_1063_1_4960790.pdf4.11 MBAdobe PDF

OPEN

NoneView/Download

Google ScholarTM

Check

Altmetric


This item is licensed under a Creative Commons License Creative Commons