It is therefore reasonable to apply the bias on Co ions to drive their migrations and observe the associated instability of the oxygen sublattice in our enhanced sampling MD simulations. Here we note that the simulation time in the biased MD trajectories does not correspond to a realistic time scale due to the employment of the enhanced sampling approach, however, the dynamic correlations between the TM migrations and the oxygen dimer formations still can be reflected by our MD results, which are qualitatively equivalent to accelerated trajectories induced by the applied bias on the collective variable (CV).