Ab initio molecular dynamics simulations of low energy recoil events in wurtzite AlN have been performed to determine threshold displacement energies, defect production and evolution mechanisms, role of partial charge transfer during the process, and the influence of irradiation-induced defects on the properties of AlN. The results show that the threshold displacement energies, Ed, along the direction parallel to the basal planes are smaller than those perpendicular to the basal planes. The minimum Ed values are determined to be 19 eV and 55 eV for N and Al atom, respectively, which occur along the [ ] direction. In general, the threshold displacement energies for N are smaller than those for Al atom, indicating the N defects would be dominant under irradiation. The defect production mechanisms have been analyzed. It is found that charge transfer and redistribution for both the primary knock-on atom and the subsequent recoil atoms play a significant role in defect production and evolution. Similar to the trend in oxide materials, there is a nearly linear relationship between Ed and the total amount of charge transfer at the potential energy peak in AlN, which provides guidance on the development of charge-transfer interatomic potentials for classic molecular dynamics simulations. Finally, the response behavior of AlN to low energy irradiation is qualitatively investigated. The existence of irradiation-induced defects significantly modifies the electronic structure, and thus affects the magnetic, electronic and optical properties of AlN. These findings further enrich the understanding of defects in the wide bandgap semiconductor of AlN.
Due to lattice mismatch between epitaxial films and substrates, in-plane strain fields are produced in the thin films, with accompanying structural distortions, and ion implantation can be used to controllably engineer the strain throughout the film. Because of the strain profile, local defect energetics are changed. In this study, the effects of in-plane strain fields on the formation and migration of oxygen vacancies in KTaO3 are investigated using first-principles calculations. In particular, the doubly positive charged oxygen vacancy (V2+O) is studied, which is considered to be the main charge state of the oxygen vacancy in KTaO3. We find that the formation energies for oxygen vacancies are sensitive to in-plane strain and oxygen position. The local atomic configuration is identified, and strong relaxation of local defect structure is mainly responsible for the formation characteristics of these oxygen vacancies. Based on the computational results, formation-dependent site preferences for oxygen vacancies are expected to occur under epitaxial strain, which can result in orders of magnitude differences in equilibrium vacancy concentrations on different oxygen sites. In addition, all possible migration pathways, including intra- and inter-plane diffusions, are considered. In contrast to the strain-enhanced intra-plane diffusion, the diffusion in the direction normal to the strained plane is impeded under the epitaxial strain field. These anisotropic diffusion processes can further enhance site preferences.