The trajectory evolutions simulated by the enhanced sampling DeePMD are shown in ). We start with a model with one Co ion already present in a Li layer leaving a Co vacancy in a TM layer, we then employ the CI-NEB method to calculate the kinetic barriers of the intralayer Co migration order (the right panel in Figure ) and the interlayer Co migration order (the left panel in Figure ). The former case exhibits a higher barrier (1.57 eV) than that of the latter case (1.08 eV), meaning that the interlayer Co migration dynamic path indeed is likely to occur ahead of the intralayer dynamic path, which is consistent with our DeePMD trajectory's prediction.