Systematic calculation of threshold displacement energies: Case study in rutile