To carry out an MD simulation, the reported Vashishta-Rahman form of interatomic potential is employed. This method is becoming popular recently, as it describes structural and dynamical properties of a series of honeycomb layered oxides, namely Na 2 M 2 TeO 6 (M = Ni, Zn, Co and Mg), as well as NASICONs. A Na-Na = 0 eV, A Na-Li = 80000 eV, A Na-Fe = 27000 eV, A Na-Te = 7000 eV, n Na-Na = n Na-Li = n Na-Fe = n Na-Te = 11 a The value used from the reported literature .