An electrostatic spectral neighbor analysis potential for lithium nitride

Machine-learned interatomic potentials based on local environment descriptors represent a transformative leap over traditional potentials based on rigid functional forms in terms of prediction accuracy. However, a challenge in their application to ionic systems is the treatment of long-ranged electrostatics. Here, we present a highly accurate electrostatic Spectral Neighbor Analysis Potential (eSNAP) for ionic α-Li3N, a prototypical lithium superionic conductor of interest as a solid electrolyte or coating for rechargeable lithium-ion batteries. We show that the optimized eSNAP model substantially outperforms traditional Coulomb–Buckingham potential in the prediction of energies and forces, as well as various properties, such as lattice constants, elastic constants, and phonon dispersion curves. We also demonstrate the application of eSNAP in long-time, large-scale Li diffusion studies in Li3N, providing atomistic insights into measures of concerted ionic motion (e.g., the Haven ratio) and grain boundary diffusion. This work aims at providing an approach to developing quantum-accurate force fields for multi-component ionic systems under the SNAP formalism, enabling large-scale atomistic simulations for such systems.

[1]  U. Holland,et al.  Local ionic motion in the superionic conductor Li3N , 1978 .

[2]  R. Buckingham,et al.  The Classical Equation of State of Gaseous Helium, Neon and Argon , 1938 .

[3]  B. Boukamp,et al.  Fast ionic conductivity in lithium nitride , 1978 .

[4]  Gábor Csányi,et al.  Accuracy and transferability of Gaussian approximation potential models for tungsten , 2014 .

[5]  K. Differt,et al.  NMR study of diffusion in Li3N , 1981 .

[6]  A Modified Embedded-Atom Method Interatomic Potential for Ionic Systems: 2NN MEAM-Qeq , 2016 .

[7]  Aidan P Thompson,et al.  Extending the accuracy of the SNAP interatomic potential form. , 2017, The Journal of chemical physics.

[8]  I. Tanaka,et al.  First principles phonon calculations in materials science , 2015, 1506.08498.

[9]  J. Ziegler The stopping and range of ions in solids vol 1 : The stopping and ranges of ions in matter , 2013 .

[10]  R. Kondor,et al.  Gaussian approximation potentials: the accuracy of quantum mechanics, without the electrons. , 2009, Physical review letters.

[11]  Cormac Toher,et al.  Charting the complete elastic properties of inorganic crystalline compounds , 2015, Scientific Data.

[12]  Michele Parrinello,et al.  Generalized neural-network representation of high-dimensional potential-energy surfaces. , 2007, Physical review letters.

[13]  Christian Trott,et al.  Spectral neighbor analysis method for automated generation of quantum-accurate interatomic potentials , 2014, J. Comput. Phys..

[14]  A. Rabenau,et al.  Ionic conductivity in Li3N single crystals , 1977 .

[15]  Shyue Ping Ong,et al.  Accurate Force Field for Molybdenum by Machine Learning Large Materials Data , 2017, 1706.09122.

[16]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[17]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[18]  Allen Pei,et al.  Robust Pinhole-free Li3N Solid Electrolyte Grown from Molten Lithium , 2017, ACS central science.

[19]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[20]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[21]  P. Hohenberg,et al.  Inhomogeneous Electron Gas , 1964 .

[22]  P. P. Ewald Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[23]  R. Bratschitsch,et al.  Assignment of the NV 0 575-nm zero-phonon line in diamond to a 2 E- 2 A 2 transition , 2013, 1301.3542.

[24]  Nongnuch Artrith,et al.  High-dimensional neural-network potentials for multicomponent systems: Applications to zinc oxide , 2011 .

[25]  P. Rinke,et al.  Multiscale model for disordered hybrid perovskites: The concept of organic cation pair modes , 2017, Physical Review B.

[26]  J. Behler Atom-centered symmetry functions for constructing high-dimensional neural network potentials. , 2011, The Journal of chemical physics.

[27]  G. Henkelman,et al.  A climbing image nudged elastic band method for finding saddle points and minimum energy paths , 2000 .

[28]  Blöchl,et al.  Projector augmented-wave method. , 1994, Physical review. B, Condensed matter.

[29]  Anubhav Jain,et al.  Python Materials Genomics (pymatgen): A robust, open-source python library for materials analysis , 2012 .

[30]  A. Rabenau,et al.  Re-evaluation of the lithium nitride structure , 1976 .

[31]  R. Kondor,et al.  On representing chemical environments , 2012, 1209.3140.

[32]  M. Baskes,et al.  Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals , 1984 .

[33]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[34]  C. Catlow,et al.  Potential models for ionic oxides , 1985 .

[35]  Jian Luo,et al.  Quantum-accurate spectral neighbor analysis potential models for Ni-Mo binary alloys and fcc metals , 2018, Physical Review B.

[36]  Alexander V. Shapeev,et al.  Moment Tensor Potentials: A Class of Systematically Improvable Interatomic Potentials , 2015, Multiscale Model. Simul..

[37]  Yizhou Zhu,et al.  Strategies Based on Nitride Materials Chemistry to Stabilize Li Metal Anode , 2017, Advanced science.