An efficient algorithm for finding the minimum energy path for cation migration in ionic materials
The Nudged Elastic Band (NEB) is an established method for finding minimumenergy 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 nondiffusing species, causing instabilities in the ensuing calculation. Here, 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 ofmore »
 Authors:

^{[1]};
^{[1]};
^{[2]}
;
^{[1]};
^{[3]}
 Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States). Dept. of Materials Science and Engineering
 Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States). Dept. of Materials Science and Engineering; Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Materials Sciences Division
 Massachusetts Inst. of Technology (MIT), Cambridge, MA (United States). Dept. of Materials Science and Engineering; Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Materials Sciences Division; Univ. of California, Berkeley, CA (United States). Dept. of Materials Science and Engineering
 Publication Date:
 Grant/Contract Number:
 AC0205CH11231; OCI1147503; 3F31144; EDCBEE
 Type:
 Published Article
 Journal Name:
 Journal of Chemical Physics
 Additional Journal Information:
 Journal Volume: 145; Journal Issue: 7; Related Information: © 2016 Author(s).; Journal ID: ISSN 00219606
 Publisher:
 American Institute of Physics (AIP)
 Research Org:
 Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
 Sponsoring Org:
 USDOE Office of Science (SC), Basic Energy Sciences (BES) (SC22). Joint Center for Energy Storage Research (JCESR); National Science Foundation (NSF)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 25 ENERGY STORAGE; 36 MATERIALS SCIENCE; 97 MATHEMATICS AND COMPUTING; 37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY
 OSTI Identifier:
 1298334
 Alternate Identifier(s):
 OSTI ID: 1416455; OSTI ID: 1474902
Rong, Ziqin, Kitchaev, Daniil, Canepa, Pieremanuele, Huang, Wenxuan, and Ceder, Gerbrand. An efficient algorithm for finding the minimum energy path for cation migration in ionic materials. United States: N. p.,
Web. doi:10.1063/1.4960790.
Rong, Ziqin, Kitchaev, Daniil, Canepa, Pieremanuele, Huang, Wenxuan, & Ceder, Gerbrand. An efficient algorithm for finding the minimum energy path for cation migration in ionic materials. United States. doi:10.1063/1.4960790.
Rong, Ziqin, Kitchaev, Daniil, Canepa, Pieremanuele, Huang, Wenxuan, and Ceder, Gerbrand. 2016.
"An efficient algorithm for finding the minimum energy path for cation migration in ionic materials". United States.
doi:10.1063/1.4960790.
@article{osti_1298334,
title = {An efficient algorithm for finding the minimum energy path for cation migration in ionic materials},
author = {Rong, Ziqin and Kitchaev, Daniil and Canepa, Pieremanuele and Huang, Wenxuan and Ceder, Gerbrand},
abstractNote = {The Nudged Elastic Band (NEB) is an established method for finding minimumenergy 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 nondiffusing species, causing instabilities in the ensuing calculation. Here, 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 highthroughput 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 longstanding obstacle to the computational screening of intercalation compounds for Liion and multivalent batteries.},
doi = {10.1063/1.4960790},
journal = {Journal of Chemical Physics},
number = 7,
volume = 145,
place = {United States},
year = {2016},
month = {8}
}
Works referenced in this record:
A climbing image nudged elastic band method for finding saddle points and minimum energy paths
journal, December 2000
journal, December 2000
 Henkelman, Graeme; Uberuaga, Blas P.; Jónsson, Hannes
 The Journal of Chemical Physics, Vol. 113, Issue 22, p. 99019904
SelfConsistent Equations Including Exchange and Correlation Effects
journal, November 1965
journal, November 1965
 Kohn, W.; Sham, L. J.
 Physical Review, Vol. 140, Issue 4A, p. A1133A1138
Efficient iterative schemes for ab initio totalenergy calculations using a planewave basis set
journal, October 1996
journal, October 1996
 Kresse, G.; Furthmüller, J.
 Physical Review B, Vol. 54, Issue 16, p. 1116911186
Alloys of platinum and early transition metals as oxygen reduction electrocatalysts
journal, September 2009
journal, September 2009
 Greeley, J.; Stephens, I. E. L.; Bondarenko, A. S.
 Nature Chemistry, Vol. 1, Issue 7, p. 552556