Multiple cascade radiation damage simulations of pyrochlore

Abstract

ABSTRACT We report molecular dynamics simulations of multiple radiation damage cascades in the pyrochlores Gd2Ti2O7 and Gd2Zr2O7 and the solid solution Gd2(Zr x Ti1–x )2O7 (x = 0.25, 0.50, 0.75). Using a simulation cell of 360,448 atoms, for each compound 2200 decay events are simulated over a total time of 10 ns, with each recoiling uranium atom (primary-knock-on atom) assigned initial kinetic energy of 5 keV. The structures generated are analysed using Steinhardt local order parameters. There is a large increase in volume for the Ti pyrochlore associated with a transition to an amorphous structure which resembles the melt while preserving the local environment of the Ti. The calculated dose for amorphisation is 20 eV atom−1 which compares well with experiment overlap of cascade and damage accumulation drives the amorphisation suppressessing the healing mechanisms. The behaviour of the zirconate is different – the substantial anion disorder produced by each recoil event is followed by healing and reversion to the parent pyrochlore. In the solid solution the onset of amorphisation is delayed to later times on increasing the Zr concentration and overall swelling reduced. Our simulations highlight the importance of ion mobility, associated with the weaker Zr–O bonds, in healing.