Molecular dynamics simulations of the primary irradiation damage in Zirconium

Highlights: • The arc-dpa model is more accurate in comparison with NRT model when predicting N_s. • Compared with T, E_PKA is more likely to cause damages in Zr. • T has a significant influence on A, while has little effect on B in the power law. • PKA incident direction has little effect on N_peak, t_peak and N_s except for t_s. Atom-scale numerical calculations were performed to investigate the damage behavior in close-packed hexagonal zirconium (HCP-Zr) by primary irradiation through a molecular dynamics (MD) study. The influences of Primary Knock-on Atom (PKA) energy (E_PKA), the PKA incident direction and the ambient temperature (T) on the cascade collision were studied comprehensively. The results show that the athermal recombination corrected displacements per atom (arc-dpa) model is more accurate in comparison with the Norgett-Robinson-Torrens (NRT) model when predicting the steady vacancy (N_s). It was found that the vacancy peak (N_peak), the peak time (t_peak), and the steady time (t_s) increase as E_PKA and T increase. The steady vacancy also increases with the increase of E_PKA. It was also found that N_s increases and subsequently decreases by increasing T, suggesting that there is a suitable temperature to maximize N_s when E_PKA is constant. In addition, it was discovered that the PKA incident direction has little effect on vacancy dynamic history. It was proved that E_PKA can worsen the irradiation damage in crystal Zr, but T can relieve such damage. In addition, the incident direction was found to have an insignificant effect on the damage. This study highlights that the steady vacancy concentration (C) can characterize the material irradiation damage.