Prediction and characterisation of radiation damage in fluorapatite

Molecular dynamics simulations, used in conjunction with a set of classical pair potentials, have been employed to examine simulated radiation damage cascades in the fluorapatite structure. Regions of damage have subsequently been assessed for their ability to recover and the effect that damage has on the important structural units defining the crystal structure, namely phosphate tetrahedra and calcium meta-prisms. Damage was considered by identifying how the phosphorous coordination environment changed during a collision cascade. This showed that PO 4 units are substantially retained, with only a very small number of under or over coordinated phosphate units being observed, even at peak radiation damage. By comparison the damaged region of the material showed a marked change in the topology of the phosphate polyhedra, which polymerised to form chains up to seven units in length. Significantly, the fluorine channels characteristic of the fluorapatite structure and defined by the structure's calcium meta-prisms stayed almost entirely intact throughout. This meant that the damaged region could be characterised as amorphous phosphate chains interlaced with regular features of the original undamaged apatite structure. Displacement cascade simulations in fluorapatite show selective amorphisation, with the phosphate sub-structure forming amorphous chains interlaced with much more damage-resistant calcium meta-prisms.