In this study, an ab initio molecular dynamics method is employed to study low energy recoil events in zirconate pyrochlores (A 2 Zr 2 O 7 , A = La, Nd and Sm). It shows that both cations and anions in Nd 2 Zr 2 O 7 and Sm 2 Zr 2 O 7 are generally more likely to be displaced than those in La 2 Zr 2 O 7 . The damage end states mainly consist of Frenkel pair defects, and the Frenkel pair formation energies in Nd 2 Zr 2 O 7 and Sm 2 Zr 2 O 7 are lower than those in La 2 Zr 2 O 7 . These results suggest that the order-disorder structural transition more easily occurs in Nd 2 Zr 2 O 7 and Sm 2 Zr 2 O 7 resulting in a defect-fluorite structure, which agrees well with experimental observations. Our calculations indicate that oxygen migration from 48f and 8b to 8a sites is dominant under low energy irradiation. A number of new defects, including four types of cation Frenkel pairs and six types of anion Frenkel pairs, are revealed by ab initio molecular dynamics simulations. The present results may help to advance the fundamental understanding of the irradiation response behavior of zirconate pyrochlores.