We report the simulation of alpha recoil cascades in zirconolite (CaZrTi(2)O(7)) using the molecular dynamics method under a variety of simulation conditions. Interatomic potentials have been developed through fitting to crystalline and melted structures and their energy surfaces obtained from ab initio calculations. Single and double cascades have been modelled, showing that damage recovery in damaged crystalline material is greatly impaired when there are multiple events. The calcium ions are significantly easier to displace than other ions, which is consistent with a recent ab initio study of defect energies.