Molecular dynamics simulation of energetic uranium recoil damage in zircon

Defect production and amorphisation due to energetic uranium recoils in zircon (ZrSiO₄), which is a promising ceramic nuclear waste form, is studied using molecular dynamics (MD) simulations with a partial charge model. An algorithm that distinguishes between undamaged crystal, crystalline defects and amorphous regions is used to develop a fundamental understanding of the primary damage state. The amorphous cascade core is separated from the surrounding crystal by a defect-rich region. Small, chemically inhomogeneous amorphous clusters are also produced around the core. The amorphous regions consist of under-coordinated Zr and polymerised Si leading to amorphisation and phase separation on a nanometer scale into Zr- and Si-rich regions. This separation could play an important role in the experimentally observed formation of nanoscale ZrO₂ in ZrSiO₄ irradiated at elevated temperatures.