Electrochemical methods such as potentiostatic and galvanostatic intermittent titration techniques (PITT and GITT) have been frequently utilized to experimentally quantify lithium diffusivity in NMC. These techniques need the knowledge of electrode particle shape and dimension, and uncertainty about these parameters leads to substantial errors in predicting the diffusion coefficient. In addition, because these techniques consider the response of the whole electrochemical cell, it is hard to distinguish the effect of different structural factors on Li diffusivity in a single NMC active material. Therefore, an appropriate method still needs to be developed to capture the structural effects on lithium diffusivity in NMC. For this purpose, a multi-level modelling from Density Functional Theory (DFT) to kinetic Monte Carlo (KMC) should be implemented.
In this study, we will use DFT to find the ground state energy of NMC at different lithium concentrations and configurations. Also, the minimum energy path of lithium migration and the related activation barrier will be found by Climbing Image-Nudge Elastic Band (CI-NEB) method. Then by implementing the configurational dependent activation barrier into the KMC simulation, the lithium diffusivity will be studied. This atomistic simulation gives insight about the structural effects on lithium diffusivity in NMC to further develop this cathode material for high performance LIBs.