Next Article in Journal
Structural Features of Fatigue Crack Propagation of a Forging Die Made of Chromium–Molybdenum–Vanadium Tool Steel on Its Durability
Previous Article in Journal
Experimental Study on the Shear Strength and Failure Mechanism of Cemented Soil–Concrete Interface
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Materials
Volume 16
Issue 12
10.3390/ma16124224
materials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Dynamics Study on the Mechanism of Gallium Nitride Radiation Damage by Alpha Particles

by
Yang Liu
1,2,3,*,
Zhenpeng Xiong
1 and
Xiaoping Ouyang
1,3
1
School of Nuclear Science and Engineering, North China Electric Power University, Beijing 102206, China
2
State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang 330013, China
3
State Key Laboratory of Intense Pulsed Radiation Simulation and Effect, Xi’an 710024, China
*
Author to whom correspondence should be addressed.
Materials 2023, 16(12), 4224; https://doi.org/10.3390/ma16124224
Submission received: 23 April 2023 / Revised: 23 May 2023 / Accepted: 25 May 2023 / Published: 7 June 2023
Browse Figures
Versions Notes

Abstract

:
In special applications in nuclear reactors and deep space environments, gallium nitride detectors are subject to irradiation by α-particles. Therefore, this work aims to explore the mechanism of the property change of GaN material, which is closely related to the application of semiconductor materials in detectors. This study applied molecular dynamics methods to the displacement damage of GaN under α-particle irradiation. A single α-particle-induced cascade collision at two incident energies (0.1 and 0.5 MeV) and multiple α-particle injections (by five and ten incident α-particles with injection doses of 2 × 1012 and 4 × 1012 ions/cm2, respectively) at room temperature (300 K) were simulated by LAMMPS code. The results show that the recombination efficiency of the material is about 32% under 0.1 MeV, and most of the defect clusters are located within 125 Å, while the recombination efficiency of 0.5 MeV is about 26%, and most of the defect clusters are outside 125 Å. However, under multiple α-particle injections, the material structure changes, the amorphous regions become larger and more numerous, the proportion of amorphous area is about 27.3% to 31.9%, while the material’s self-repair ability is mostly exhausted.

1. Introduction

In the future, deep space exploration will require multi-functional electronic devices [1,2]. As one of the cornerstones of the modern electronics industry, semiconductor materials have undergone a long development process. With significant improvements in manufacturing technology and application areas, researchers have further studied and developed third-generation semiconductors based on second-generation semiconductors, such as silicon carbide (SiC) and gallium nitride (GaN). Third-generation semiconductors have higher performance, lower power consumption, and broader application fields than the first two generations of semiconductors. Therefore, third-generation semiconductors are considered the cutting edge of semiconductor technology. Gallium nitride semiconductor materials, as essential third-generation semiconductor materials, have a large bandgap (3.1 eV), a stable high breakdown field (3.5 × 106 V/m), high electron saturation velocity (2 MV/cm), and high thermal conductivity (1.3 W/cm·K). As a result, they are widely used in modern electronic devices, detectors, and other fields. In addition, gallium nitride also has excellent heat- and radiation-resistant characteristics, making it widely applicable in military and aerospace domains [2,3,4].
When the material is working in a harshly irradiated environment, radiation will result in severe degradation of material properties. Many studies have observed this phenomenon in different materials [5,6]. Radiation particles interact with semiconductor materials, mainly through nuclear scattering or nuclear reactions, which deposit a large amount of energy within the semiconductor material and result in the deposition of recoil atoms or reaction products inside the material [7]. Nuclear energy development has two main types: nuclear fusion energy and nuclear fission energy. Compared to nuclear fission reactions, nuclear fusion produces many α-particles. As a third-generation semiconductor detector, studying the radiation damage of GaN caused by α-particles is an essential topic in semiconductor development. Using α-particles to irradiate semiconductor materials directly is the most direct and effective method to study radiation damage effects, and a large number of semiconductor materials irradiated by α-particles have been investigated through various domestic and international experimental reactors, particle accelerators, or other radioactive sources, accumulating a wealth of experimental data. Although current experimental equipment can observe or identify specific defects in semiconductors, it has certain limitations in observing the entire irradiation process [8]. Therefore, relying solely on irradiation experiments cannot simulate every detail of the interaction between particles and matter, nor can it be used to analyze the causes of these defects. Moreover, the lengthy experimental cycles and high costs pose significant research barriers. In addition to experiments, the interaction between particles and matter can be simulated by computer. At present many computer codes have been developed to study the radiation damage. However, due to the limitation of computing power, the more detailed the calculation results, the smaller the calculation scale. In this work, the molecular dynamics method is used, which can simulate the damage effect of large-scale irradiation at the molecular level, and can clearly observe the damage-generation process.
In recent years, many scientists have used computer simulations to study the radiation damage to materials. He et al. [7] studied the displacement damage of GaN after 1, 10, 100, and 500 MeV proton irradiations. They reported that the NIEL deposited per unit thickness decreased with increasing proton energy. Li [9] found that 30 keV helium ions injected into monocrystalline silicon would produce dislocation, and the energy level and band structure would be affected. He et al. [8,10] simulated Ga PKA bombardment of the GaN matrix at 300 K to 900 K and found that the displacement threshold energy of GaN increases with the increase of temperature, and the defect trend is similar, but the higher the temperature, the higher the recombination efficiency. Xie et al. [11] used the Monte Carlo method to study neutron irradiation of GaN under different irradiation environments. The results showed that the PKA energy was lower, and the distribution range was wider under atmospheric neutron spectra. Ullah et al. [12,13,14] studied the defect production process in GaN under 50 eV/amu single atomic and molecular ion irradiation. The study showed that the damage varied nonlinearly with increasing molecular weight. Khanal et al. [15] conducted proton irradiation experiments on GaN at 100 KeV and found that the device performance decreased with increasing proton flux. Kucheyev et al. [16,17] studied the cascade density of GaN under keV cluster ion bombardment at room temperature and liquid nitrogen temperature. The study showed that the cascade density increased with increasing temperature, and the damage efficiency of molecular ion irradiation of GaN was higher than that of its single constituent ions. Karaseov et al. [18] used MD simulations to study the formation of damage caused by P, Ag, and PF4 ions on GaN. The results showed that the damage caused by Ag was more significant, and the large defect clusters increased nonlinearly. Chen et al. [19] studied the production and evolution of defects in GaN under 500 eV–40 KeV PKA irradiation. The study showed that many atoms would undergo displacement in a compact cascade volume during the collision stage under proton irradiation. Still, a significant combination of vacancies co-occurred, called pseudo-metallic behavior (PMB). Hosseini et al. [20] used MD simulations to study the defect and temperature effects of Ar atom bombardment on GaN. The results showed that as the temperature increased from 300 K to 350 K, the atomic loss in GaN increased, and the atomic mechanical stability decreased. Li [21] used molecular dynamics simulations to study the injection of α-particles into single-crystal silicon, and the results showed that a dislocation structure was produced during the injection process. Fan et al. [22] used Monte Carlo calculations to simulate the irradiation of aluminum materials by protons and α-particles. The results showed that the damage caused by α-particles was more significant than that caused by protons. Pan et al. [23] used computer simulations to study the radiation damage effects of fusion α-particle bombardment on amorphous targets (diamond and iron targets). It was concluded that the peak of elastic energy deposition was 5 eV/Å in both targets. The number of vacancies was proportional to the elastic energy deposition and inversely proportional to the displacement threshold energy. Li et al. [24] simulated the epitaxial growth of GaN films at different particle incident energies. The results show that increasing the incident energy of particles can improve the epitaxial surface mass.
Molecular dynamics simulation (MD) is a popular and reliable computer simulation method that can describe molecules’ motion at the microscopic level. Recently, MD simulation has been used to calculate the atomic/mechanical properties of materials [25,26,27,28,29]. In addition to the MD methods described above, ab initio molecular dynamics studies (AIMD) have also been widely used in recent years and have been successfully applied to group IIIA nitrides and group IIIA binaries [30,31]. All these approaches have become popular recently and will become important in the future. Therefore, this article uses molecular dynamics simulation to study the radiation damage mechanism of α-particles on GaN materials, which can better understand the evolution of damage inside the material during collision cascades and further understand the damage mechanism of α-particle irradiation on GaN. The innovation point of this work is that the GaN model is an imperfect matrix with vacancy defects and doped Fe atoms. Secondly, the α-particle is selected as the incident particle, and the single incident particle is compared with the cumulative incident particle. This work is of great significance for studying the irradiation effects of GaN-based devices in space radiation environments or nuclear fusion reactor environments.

2. Simulation Method

The MD simulation method can be used to simulate the process of material defect generation [19]. In this article, all simulations were performed using the open-source parallel MD code LAMMPS, which has been used to simulate the radiation damage of materials [32,33,34]. The results were analyzed using the visualization software OVITO [35,36,37]. At room temperature, GaN consistently exhibits the crystal structure of wurtzite, which is a hexagonal crystal system. The lattice constants of wurtzite GaN are a = b = 3.216   Å and c = 5.240   Å . In the simulation, we constructed a perfect GaN matrix with a size of 50 a × 50 a × 50 a containing 500,000 atoms. Then, based on the perfect GaN matrix, we constructed an imperfect GaN matrix with a vacancy defect rate of 10% and doped it with iron (Fe) atoms at a doping density of 5 × 10 18   c m - 3 . To simulate the radiation of α-particles on the GaN matrix, we added two positively charged nitrogen atoms as α-particles above the system in the z-direction and gave the α-particles an initial velocity to shoot at the GaN matrix as PKA. To avoid channeling effects, a slight angle was added to the PKA incident direction [8]. In the MD simulation calculation, the particle trajectory was obtained by solving Newton’s equation. The most important factor was finding an appropriate energy relationship formula to describe the interaction between particles, usually represented by a potential function. Therefore, we used the Tersoff potential to describe the interaction between GaN particles [19,38], then used the Lennard-Jones (LJ) potential to describe the interaction forces between GaN and Fe atoms as well as α-particles [39,40], and finally smoothed the Tersoff and LJ potential with the Ziegler–Biersack–Littmark (ZBL) potential so that high-energy short-range interactions within the material could be better described [19,41]. See Table 1 for specific parameter settings.
Periodic boundary conditions were applied along the x- and y-directions for the simulation system, while fixed boundary conditions were used along the z-direction [20,42,43]. The entire simulation process can be divided into three stages. The first stage is the relaxation stage: we first relaxed the entire simulation system by conjugate gradient minimization to release the stress in the system [44]. Then we allowed the system to relax under the NVT ensemble for 10 ps, during which the system temperature was 300 K, and the time step (∆t) was set to 0.0001 ps until the system reached a stable state [45]. In this stage, the temperature and potential energy of the system can reach equilibrium, making the structure reach a metastable state, as shown in Figure 1.
In Figure 1, it can be observed that the temperature and potential energy change trends during the relaxation stage are opposite, and both eventually approach stability, indicating that the system has reached a stable state and can proceed to the next simulation step. The second stage is the collision stage: two perpendicular constant-temperature walls were used [8,45,46], as shown in Figure 2a. Figure 2b shows the view from the z-direction under periodic boundary conditions, with black dots representing the position of the incident α-particle.
The thickness of the constant-temperature walls was set to 2a, and the Langevin thermostat was used to maintain a constant temperature of 300 K for the atoms inside the walls so that excess kinetic energy introduced by PKA could be dissipated by the walls [47]. A fixed layer with a thickness of about 3a was selected as the bottom layer of the system to prevent the entire simulation system from drifting during the collision cascade process [18,21]. The remaining atoms were in the active region, and all atoms in the active region were restricted to adiabatic motion using the NVE ensemble [44,45]. In addition, to ensure that the atomic displacement in each time step did not exceed 0.1 Å, we used a variable time step (∆t) [44,47]. The final stage is the annealing stage: in this stage, the time step (∆t) was increased to 0.1 ps, and the system was annealed for a long time until there was no significant change in the total number of defects [19,45].
The Wigner–Seitz cell recognition method was used to identify defects caused by irradiation in the simulation system, such as vacancies, interstitials, and anti-site atoms [35,36]. A vacancy is a defect with no atom on a lattice site, while an interstitial defect is when two or more atoms are on a lattice site. An anti-site atom defect occurs when an atom of the wrong type occupies a lattice site. Defects within a cutoff distance of 1.5a were considered defect clusters [13], and the size of a cluster was defined as the number of internal defects. Different sizes of clusters were classified as small clusters (containing from one to six defects), medium-sized clusters (containing 6–30 defects), large clusters (containing more than 30 defects), and local amorphization (when a cluster contains more than 100 defects).
In this study, MD simulations were performed on GaN irradiated by α-particles. First, single α-particle irradiation on GaN was simulated, corresponding to an injection dose of 4 × 1011 ions/cm2. Two different energies were selected as the incident energy of the α-particle: 0.1 and 0.5 MeV. Then, multiple α-particles were simulated to accumulate irradiation damage in GaN, with 5 and 10 α-particles used for cumulative irradiation, corresponding to injection doses of 2 × 1012 and 4 × 1012 ions/cm2, respectively.

3. Results and Analysis

3.1. Simulation of Radiation Damage in GaN by Single α-Particle Irradiation

The relationship between the number of defects caused by the irradiation of single α-particles with different energies and times was simulated by the MD method. Figure 3a shows the variation in the number of defects with time under 0.1 MeV α-particle irradiation. Figure 3b shows the variation in the number of defects with time under 0.5 MeV α-particle irradiation.
It can be seen from the figure that under irradiation with either energy level of α-particles, the number of defects increases first and then decreases to a steady state, with a peak value, namely, the displacement peak. This is because α-particles impact GaN, transfer energy to GaN through collisions, and cause displacement effects in the GaN system; the energy deposited in GaN is dissipated by the Langevin thermostat of the constant temperature wall, and GaN has a self-repairing ability, so the number of defects begins to decrease after reaching the peak value and eventually stabilizes. Compared with the displacement peak formed under 0.5 MeV irradiation, the time of the displacement peak formed under 0.1 MeV is slightly earlier. It can also be observed that the change trend of vacancy defects is the same as that of interstitial defects, and the number of vacancy defects slightly exceeds that of interstitial defects. This is because some interstitial defects with two or more atoms are formed during irradiation. Compared with the vacancy and interstitial defects, the figure shows little change in the peak and steady-state stages of anti-site defects, indicating that anti-site defects are more challenging to repair.
The defect distribution in the x-direction is shown in Figure 4. Figure 4a,b show the defect distribution under 0.1 MeV, where (a) is the distribution at the displacement peak stage and (b) is the distribution at the final steady state. Figure 4c,d show the defect distribution under 0.5 MeV at the displacement peak stage and the steady-state stage, respectively. It can be seen from the figure that after long-time annealing, the number of defects is significantly reduced, and the defects caused by 0.1 MeV are primarily concentrated in shallow positions. In comparison, those caused by 0.5 MeV are concentrated in deep positions. In addition, it can be observed that local amorphization may occur in areas with high defect density. By comparing the results under 0.1 MeV and 0.5 MeV, the results show that the higher the energy level of α-particles, the more severe the damage caused by them.
To compare the self-recovery ability of GaN under different irradiation energies, we recorded its recombination efficiency in Table 2. As shown in Table 2, in GaN with a fixed size, under the irradiation of α-particles with different energies, the higher the energy is, the more severe the displacement damage caused and the more difficult it is for the material to self-repair. Under α-particle irradiation, we used the Wigner–Seitz cell recognition method to identify and analyze the type and quantity of each defect, as shown in Table 2.
The Ga-N (Ga) interstitial defect in Table 2 represents a vacancy containing both Ga and N atoms at the original Ga site, where (Ga) is the original Ga site, and Ga-N represents the Ga and N atoms. In the case of low-energy α-particle irradiation, the material has a stronger ability to repair Ga vacancies than N vacancies. In contrast, under high-energy α-particle irradiation, N vacancies have a stronger repair ability than Ga vacancies. In addition to interstitial defects containing two atoms, interstitial defects containing three and four atoms can also be produced under α-particle irradiation. Among the interstitial defects, Ga-N (Ga) and Ga-N (N) interstitial defects are easier to form, followed by Ga-Ga (Ga) and N-N (N) interstitial defects. Interstitial defects containing three atoms are relatively difficult to form, while interstitial defects containing four atoms are the most difficult to form and are not easily repaired by the material’s self-recovery ability once formed. Compared with 0.1 MeV irradiation, under 0.5 MeV irradiation, although the repair abilities of some types of defects increased and some decreased, the overall recombination efficiency decreased. This again indicates that the higher the energy of the incident α-particles, the more difficult it is to repair the displacement damage formed.
As shown in Figure 4, the defects formed after irradiation can also form defect clusters due to the aggregation of defects. When the defect clusters are too dense, they can form local amorphization. The presence of defect clusters and local amorphization in the material can seriously affect semiconductor materials’ mechanical and electrical properties. Figure 5a shows that more defect clusters were produced under 0.5 MeV α-particle irradiation, and the clusters were located deeper. More defect clusters produced under α-particle irradiation were small clusters, with more clusters concentrated within 125 Å under 0.1 MeV. Most positions under 0.5 MeV were at depths of more than 125 Å. Figure 5b shows the number and location of local amorphization under irradiation with different energies. Under irradiation with higher energy, more local amorphous regions will be produced within the crystal, and their positions will be deeper.

3.2. Radiation Damage Simulation of Multiple α-Particle Cumulative Irradiation on GaN

After simulating the single α-particle irradiation of the GaN matrix, we also studied the cumulative damage effects caused by multiple α-particle irradiations on GaN. In this case, the energy of the incident α-particles was set to 0.1 MeV, and the injection doses were 2 × 1012 and 4 × 1012 ions/cm2, corresponding to five and ten incident α-particles, respectively. Figure 6 shows the radial distribution function of the GaN material at two injection rates [48].
The red curve represents the radial distribution function in the initial state, while the green curve represents the radial distribution function in the final state. In the absence of α-particle irradiation, the number of first-nearest neighboring atoms is nearly 13, and the number of atoms from the second-nearest neighbor to the nth-nearest neighbor decreases slowly. At this time, the radial distribution function shows a long- and short-range order and is relatively ordered. After accumulating five α-particles, the number of first-nearest neighboring atoms in the GaN material decreased to about nine, as shown in Figure 6a, and the number of atoms from the second-nearest neighbor to the nth-nearest neighbor also decreased. Figure 6b shows that after accumulating ten α-particles, the number of first-nearest neighboring atoms in the GaN material also decreased to the value slightly exceeding nine, and the change in the number of atoms from the second-nearest neighbor follows the same pattern as in Figure 6a. At this time, regardless of the injection rate, the long-range order of the radial distribution function became disordered, indicating that under the cumulative irradiation of multiple α-particles, the internal amorphization of GaN material became increasingly severe.
After accumulating irradiation, we further analyzed the damage inside the GaN material by examining the IDS damage morphology. Figure 7 shows the IDS surface topography in the z-direction, and Table 3 shows the percentage of each component in the IDS damage analysis.
As seen in Figure 7, a large hole appears on the surface of the GaN matrix, indicating that some surface atoms of the GaN matrix have been released from the confinement of the matrix lattice and sputtered out. The sputtering rates under different simulation conditions are listed in Table 4.
Figure 7a shows the IDS surface topography after accumulating five α-particle irradiations, and Figure 7b shows the IDS surface topography after accumulating ten α-particle irradiations. In both figures, a large hole appears on the surface of the GaN material due to the irradiation of the α-particles. Moreover, the entire hole has turned into an amorphous state on the IDS damage surface topography, indicating that the damage to the material has increased under the accumulated irradiation. Table 3 shows that there was already a 16.3% amorphous region in the initial state due to the imperfect GaN material with defects and doping used in this simulation. With the increased number of incident α-particles, the proportion of the amorphous region in the GaN matrix also increases. Under accumulated irradiation, the proportions of the hexagonal diamond, hexagonal diamond (1st neighbor), and hexagonal diamond (2nd neighbor) all decrease, with the hexagonal diamond decreasing the fastest. This suggests that the hexagonal diamond structure is the most susceptible to damage under the irradiation of multiple α-particles. Although the proportions of the hexagonal diamond (1st neighbor) and hexagonal diamond (2nd neighbor) also decrease under irradiation, they decrease slowly since the hexagonal diamond structure can transform into a hexagonal diamond (1st neighbor) or hexagonal diamond (2nd neighbor) after damage.
With increasing incident α-particle energy, the number of sputtered atoms and the sputtering rate also increases. Under cumulative irradiation, the sputtering rate increases significantly and increases with the number of injected particles. Therefore, large and deep damage holes will appear in the surface topography under cumulative injection.
Table 5 shows the types and quantities of defects under different cumulative irradiation conditions. It can be seen that with the increase in the number of α-particle injections, the number of each type of defect also increases. Compared with single α-particle injection, more types of damage are formed when multiple α-particles are injected, including gap defects containing more than five atoms. This indicates that the displacement damage of the GaN material is severe under cumulative irradiation and is not easily repaired by the material’s self-repairing ability. The table also shows that the number of N vacancies is always greater than that of Ga vacancies, indicating that the formation of N vacancies is more likely under α-particle irradiation. The number of vacancy defects is much greater than that of interstitial defects, which is consistent with the surface topography of IDS damage shown in Figure 7, where many atoms have been ejected from the matrix, leading to the formation of a large number of vacancy defects inside the matrix.
After comparing Table 2 and Table 5, it can be observed that under α-particle irradiation, the number of vacancy defects is higher than that of interstitial defects, and N vacancies are more likely to form. Those containing two atoms are more likely to form among the interstitial defects, while those containing more atoms are less likely to form. However, only under accumulated irradiation will interstitial defects containing more than five atoms form. Among the interstitial defects containing two atoms, Ga-N (Ga) and Ga-N (N) defects are more likely to form. To further understand the defect distribution and amorphization degree of GaN under accumulated irradiation.
Figure 8 shows the distribution of defect clusters. Figure 8a shows the distribution of defect clusters after accumulated irradiation of five α-particles, where a superlarge amorphous region containing more than 30,000 defects is formed. Figure 8b shows the distribution of defect clusters after accumulated irradiation of ten α-particles. A superlarge amorphous region containing more than 50,000 defects and several small amorphous regions containing 200–700 defects are formed. Under accumulated irradiation, most defect clusters are small or medium-sized, while the number of large clusters increases with the number of injected particles. At the same time, the number and size of amorphous regions also increase with the number of injected particles, as shown in Figure 9, where the amorphous region in Figure 9b is larger than that in Figure 9a. Furthermore, it can be observed that the amorphous region produced under the injection of five α-particles in Figure 9a is not a whole piece, but composed of several parts. The amorphous region produced by injection of ten α-particles in Figure 9b is a single block. This indicates that, at low dose injection, amorphous regions combine to produce large volume damage, while at high dose injection, large volume damage will be directly caused. At higher doses of cumulative injection, greater damage can appear further down.

4. Conclusions

This paper studied the displacement damage effect of α-particles on GaN material using the MD simulation method. The differences in damage between single α-particle irradiation at different energies and multiple α-particle injections were discussed. This paper mainly studied and compared the quantity, type, and distribution of defects, degree of amorphization, and repairability. The results obtained made it possible to draw the following conclusions. Under single α-particle irradiation, the displacement damage effect increases with increasing irradiation energy. N-vacancy defects are easier to form than Ga-vacancy defects, indicating that N atoms are more likely to break free from the lattice sites under irradiation. The self-repair ability of the GaN material decreases as irradiation energy increases. Defect recombination rate decreased from 32% to 26%. In addition, there is little difference in the repairability of vacancy and interstitial defects in GaN material. The defects caused by 0.1 MeV irradiation are more concentrated within 125 Å, while the defects caused by 0.5 MeV irradiation are mostly outside 125 Å. Under single α-particle irradiation, the size and depth of defect clusters increase with increasing energy, and at higher energy injection, deeper locations form larger defect clusters. In the simulation of multiple α-particle injections, the amorphous regions become larger and more numerous with the increase in injected atoms. The amorphous region increased from 16.3% to 27.3% and then to 31.9%. The sputtering rate also increases with the number of injected particles, which increased to 0.512% and 0.638%. For the study of the mechanism of heavy ion irradiation material, especially the study of latent track, a variety of models have been proposed, such as the Coulomb explosion model, thermal peak model and exciton model, which can be used to simulate and analyze the damage caused by heavy ion irradiation. In addition, experimental verification can be added to combine the microscopic parameters of molecular dynamics with the macroscopic characteristics of future devices.

Author Contributions

Conceptualization, Y.L. and X.O.; methodology, Y.L.; software, Z.X.; validation, Y.L.; formal analysis, Y.L.; investigation, Y.L. and Z.X.; resources, Y.L.; data curation, Y.L.; writing—original draft preparation, Z.X.; writing—review and editing, Y.L.; visualization, Z.X.; supervision, X.O.; project administration, Y.L.; funding acquisition, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Key R&D Program of China (No. 2022YFB1902700); National Natural Science Foundation of China (No. 11875129); Fund of State Key Laboratory of Intense Pulsed Radiation Simulation and Effect (No. SKLIPR1810); Fund of Innovation Center of Radiation Application (No. KFZC2020020402); Fund of State Key Laboratory of Nuclear Physics and Technology, Peking University (No. NPT2020KFY08); Joint Innovation Fund of China National Uranium Co., Ltd.-State Key Laboratory of Nuclear Resources and Environment, East China University of Technology (No. 2022NRE-LH-02).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sternberg, D.; Essmiller, J.; Colley, C.; Klesh, A.; Krajewski, J. Attitude Control System for the Mars Cube One Spacecraft. In Proceedings of the 2019 IEEE Aerospace Conference, Big Sky, MT, USA, 2–9 March 2019. [Google Scholar]
  2. Kodheli, O.; Lagunas, E.; Maturo, N.; Sharma, S.K.; Shankar, B.; Montoya, J.; Duncan, J.; Spano, D.; Chatzinotas, S.; Kisseleff, S. Satellite Communications in the New Space Era: A Survey and Future Challenges. IEEE Commun. Surv. Tutor. 2020, 23, 70–109. [Google Scholar] [CrossRef]
  3. Wada, Y.; Mizobata, H.; Nozaki, M.; Kobayashi, T.; Hosoi, T.; Kachi, T.; Shimura, T.; Watanabe, H. Insight into interface electrical properties of metal–oxide–semiconductor structures fabricated on Mg-implanted GaN activated by ultra-high-pressure annealing. Appl. Phys. Lett. 2022, 120, 082103. [Google Scholar] [CrossRef]
  4. Herriman, J.E.; Hellman, O.; Fultz, B. Phonon thermodynamics and elastic behavior of GaN at high temperatures and pressures. Phys. Rev. B 2018, 98, 214105. [Google Scholar] [CrossRef] [Green Version]
  5. Kato, N.; Sugiyama, M. Electron irradiation resistance of NiO/ZnO visible-light-transparent solar cells. Jpn. J. Appl. Phys. 2020, 59, 101004. [Google Scholar] [CrossRef]
  6. Geremew, A.K.; Kargar, F.; Zhang, E.X.; Zhao, S.E.; Aytan, E.; Bloodgood, M.A.; Salguero, T.T.; Rumyantsev, S.; Fedoseyev, A.; Fleetwood, D.M. Proton-irradiation-immune electronics implemented with two-dimensional charge-density-wave devices. Nanoscale 2019, 11, 8380–8386. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. He, B.W.; He, C.H.; Shen, S.S.; Chenyuan, M.L. Geant4 Simulation of Proton Displacement Damage in GaN. At. Energy Sci. Technol. 2017, 51, 543–548. [Google Scholar]
  8. He, H.; He, C.; Zhang, J.; Liao, W.; Zang, H.; Li, Y.; Liu, W. Primary damage of 10keVGaPKA in bulk GaN material under different temperatures. Nucl. Eng. Technol. 2020, 52, 1537–1544. [Google Scholar] [CrossRef]
  9. Li, R.; Zhu, R.; Chen, S.; He, C.; Xu, J. Study of damage generation induced by focused helium ion beam in silicon. J. Vac. Sci. Technol. B 2019, 37, 031804. [Google Scholar] [CrossRef]
  10. He, H.; Liao, W.L.; Wang, Y.Z.; Liu, W.B.; Zang, H.; He, C.H. Stability and interaction of cation Frenkel pair in wurtzite semiconductor materials. Comput. Mater. Sci. 2021, 196, 110554. [Google Scholar] [CrossRef]
  11. Xie, F.; Zhang, H.; Liu, F.; He, H.; Liao, W.L.; Huang, Y. Simulation of displacement damage of gallium nitride under different neutron irradiation environments. Acta Phys. Sin. 2020, 69, 9. [Google Scholar] [CrossRef]
  12. Ullah, M.W.; Kuronen, A.; Nordlund, K.; Djurabekova, F.; Karaseov, P.A.; Titov, A.I. Atomistic simulation of damage production by atomic and molecular ion irradiation in GaN. J. Appl. Phys. 2012, 112, 479. [Google Scholar] [CrossRef] [Green Version]
  13. Ullah, M.W.; Kuronen, A.; Nordlund, K.; Djurabekova, F. Effects of defect clustering on optical properties of GaN by single and molecular ion irradiation. J. Appl. Phys. 2013, 114, 1. [Google Scholar] [CrossRef]
  14. Byggmastar, J.; Granberg, F.; Nordlund, K. Effects of the short-range repulsive potential on cascade damage in iron. J. Nucl. Mater. Mater. Asp. Fission Fusion 2018, 508, 530–539. [Google Scholar] [CrossRef]
  15. Khanal, M.P.; Uprety, S.; Mirkhani, V.; Wang, S.; Yapabandara, K.; Hassani, E.; Isaacs-Smith, T.; Ahyi, A.C.; Bozack, M.J.; Oh, T.S. Impact of 100 keV proton irradiation on electronic and optical properties of AlGaN/GaN high electron mobility transistors (HEMTs). J. Appl. Phys. 2018, 124, 215702. [Google Scholar] [CrossRef]
  16. Kucheyev, S.O.; Williams, J.S.; Zou, J.; Li, G.; Jagadish, C.; Titov, A.I. Effect of ion species on implantation-produced disorder in GaN at liquid nitrogen temperature. Nucl. Instrum. Methods Phys. Res. 2002, 190, 782–786. [Google Scholar] [CrossRef]
  17. Kucheyev, S.O.; Williams, J.S.; Titov, A.I.; Li, G.; Jagadish, C. Effect of the density of collision cascades on implantation damage in GaN. Appl. Phys. Lett. 2001, 78, 2694–2696. [Google Scholar] [CrossRef]
  18. Karaseov, P.A.; Karabeshkin, K.V.; Mongo, E.E.; Titov, A.I.; Ullah, M.W.; Kuronen, A.; Djurabekova, F.; Nordlund, K. Experimental study and MD simulation of damage formation in GaN under atomic and molecular ion irradiation. Vacuum 2016, 129, 166–169. [Google Scholar] [CrossRef]
  19. Chen, N.; Rasch, E.; Huang, D.; Heller, E.R.; Gao, F. Atomic-Scale Simulation for Pseudo-Metallic Defect Generation Kinetics and Effective NIEL in GaN. IEEE Trans. Nucl. Sci. 2018, 65, 1108–1118. [Google Scholar] [CrossRef]
  20. Hosseini, A.; Nasrabadi, M.N. Investigation of vacancy defects and temperature effects on the GaN bombarding with argon atoms: Molecular dynamics simulation. Mater. Chem. Phys. 2021, 271, 124854. [Google Scholar] [CrossRef]
  21. Ji, L.; Liu, L.; Xu, Z.; Song, Y.; Fang, F. Molecular dynamics simulation on the effect of dislocation structures on the retention and distribution of helium ions implanted into silicon. Nanotechnol. Precis. Eng. 2020, 3, 7. [Google Scholar] [CrossRef]
  22. Fan, X.H.; Chen, B.; Pan, S.F. Monte Carlo simulation of proton and α-particle transport in aluminum materials. J. At. Mol. Phys. 2008, 25, 255–259. [Google Scholar]
  23. Pan, Z.Y.; Zhu, B.X. Computer simulation of α-particle irradiation damage. Comput. Phys. 1993, 10, 7. [Google Scholar]
  24. Li, R.; Wu, G.; Liang, K.; Xue, L.; Wang, S.; Liu, S. Influence of atomic incident kinetic energy on crystalline quality of epitaxial GaN thin films: A molecular dynamics study. Mater. Sci. Semicond. Process. 2022, 150, 106979. [Google Scholar] [CrossRef]
  25. Albooyeh, A.R.; Dadrasi, A.; Mashhadzadeh, A.H. Effect of point defects and low-density carbon-doped on mechanical properties of BNNTs: A molecular dynamics study. Mater. Chem. Phys. 2019, 239, 122107. [Google Scholar] [CrossRef]
  26. Hosseini, A.; Nasrabadi, M.N.; Esfandiarpour, A. Investigation of primary radiation damage near free surfaces in iron nanofoam with a model cylindrical nanovoids structure. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 2018, 439, 43–50. [Google Scholar] [CrossRef]
  27. Jolfaei, N.A.; Jolfaei, N.A.; Hekmatifar, M.; Piranfar, A.; Toghraie, D.; Sabetvand, R.; Rostami, S. Investigation of thermal properties of DNA structure with precise atomic arrangement via equilibrium and non-equilibrium molecular dynamics approaches. Comput. Methods Programs Biomed. 2020, 185, 105169. [Google Scholar] [CrossRef]
  28. Jafari, S.N.; Hakimi, Y.; Rouhi, S. Molecular dynamics investigation of the mechanical properties of two different graphyne allotropes: α-graphyne and α 2 -graphyne. Phys. E Low Dimens. Syst. Nanostruct. 2020, 119, 114022. [Google Scholar] [CrossRef]
  29. Sabetvand, R.; Ghazi, M.E.; Izadifard, M. Studying temperature effects on electronic and optical properties of cubic CH3NH3SnI3 perovskite. J. Comput. Electron. 2020, 19, 70–79. [Google Scholar] [CrossRef]
  30. Kakanakova-Georgieva, A.; Ivanov, I.G.; Suwannaharn, N.; Hsu, C.-W.; Cora, I.; Pécz, B.; Giannazzo, F.; Sangiovanni, D.G.; Gueorguiev, G.K. MOCVD of AlN on epitaxial graphene at extreme temperatures. CrystEngComm 2021, 23, 385–390. [Google Scholar] [CrossRef]
  31. Lundgren, C.; Kakanakova-Georgieva, A.; Gueorguiev, G.K. A perspective on thermal stability and mechanical properties of 2D Indium Bismide from ab initio molecular dynamics. Nanotechnology 2022, 33, 335706. [Google Scholar] [CrossRef]
  32. Liu, C.; Szlufarska, I. Distribution of defect clusters in the primary damage of ion irradiated 3C-SiC. J. Nucl. Mater. 2018, 509, 392–400. [Google Scholar] [CrossRef] [Green Version]
  33. Sahi, Q.; Kim, Y.S. Primary radiation damage characterization of α-iron under irradiation temperature for various PKA energies. Mater. Res. Express 2018, 5, 046518. [Google Scholar] [CrossRef]
  34. Mjr, A.; Mwdc, B.; Bs, C.; Jas, A. Primary radiation damage on displacement cascades in UO2, ThO2 and (U0.5Th0.5)O2—ScienceDirect. Comput. Mater. Sci. 2018, 154, 508–516. [Google Scholar]
  35. Ping, Z.A.; Xian, C.A.; Xz, A.; Yw, B. Effects of cutting parameters on the subsurface damage of single crystal copper during nanocutting process—ScienceDirect. Vacuum 2020, 187, 109420. [Google Scholar]
  36. Li, J.; Dong, L.; Dong, X.; Zhao, W.; Liu, J.; Xiong, J.; Xu, C. Study on wear behavior of FeNiCrCoCu high entropy alloy coating on Cu substrate based on molecular dynamics. Appl. Surf. Sci. 2021, 570, 151236. [Google Scholar] [CrossRef]
  37. Fu, X.; Xu, Z.; He, Z.; Hartmaier, A.; Fang, F. Molecular dynamics simulation of silicon ion implantation into diamond and subsequent annealing. Nucl. Instrum. Methods Phys. Res. 2018, 450, 51–55. [Google Scholar] [CrossRef]
  38. Nord, J.; Albe, K.; Erhart, P.; Nordlund, K. Modelling of compound semiconductors: Analytical bond-order potential for gallium, nitrogen and gallium nitride. J. Phys. Condens. Matter 2003, 15, 5649–5662. [Google Scholar] [CrossRef]
  39. Ting, L.; Ping, Z.; Peng, Y.; Siping, Z. In-plane thermal transport in black phosphorene/graphene layered heterostructures: A molecular dynamics study. Phys. Chem. Chem. Phys. 2018, 20, 21151–21162. [Google Scholar] [CrossRef]
  40. Liang, T.; Zhou, M.; Zhang, P.; Yuan, P.; Yang, D. Multilayer in-plane graphene/hexagonal boron nitride heterostructures: Insights into the interfacial thermal transport properties. Int. J. Heat Mass Transf. 2020, 151, 119395. [Google Scholar] [CrossRef]
  41. Ziegler, J.F.; Ziegler, M.D.; Biersack, J.P. SRIM—The stopping and range of ions in matter (2010). Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 2008, 268, 1818–1823. [Google Scholar] [CrossRef] [Green Version]
  42. Liang, K.; Sun, X.; Wu, G.; Zhang, L.; Gan, Z. The investigation of molecular beam epitaxy growth of GaN by molecular dynamics simulation. Comput. Mater. Sci. 2019, 173, 109426. [Google Scholar] [CrossRef]
  43. Libin, Z.; Han, Y.; Guo, Z.; Sheng, L.; Zhiyin, G. Molecular dynamics simulation of aluminum nitride deposition: Temperature and N:Al ratio effects. R. Soc. Open Ence 2018, 5, 180629. [Google Scholar]
  44. Jin, Y.; Huang, H.; Zhong, Y.; Yuan, X.; Peng, Q. Role of interface on irradiation damage of Cudiamond composites using classical molecular dynamics simulations. Ceram. Int. 2022, 48, 16813–16824. [Google Scholar] [CrossRef]
  45. Huang, H.; Tang, X.; Chen, F.; Yang, Y.; Liu, J.; Li, H.; Chen, D. Radiation damage resistance and interface stability of copper–graphene nanolayered composite. J. Nucl. Mater. 2015, 460, 16–22. [Google Scholar] [CrossRef]
  46. Tian, J.; Zhou, W.; Feng, Q.; Zheng, J. Molecular dynamics simulations with electronic stopping can reproduce experimental sputtering yields of metals impacted by large cluster ions. Appl. Surf. Sci. 2018, 435, 65–71. [Google Scholar] [CrossRef]
  47. Palko, J.W.; Srour, J.R. Amorphous Inclusions in Irradiated Silicon and Their Effects on Material and Device Properties. IEEE Trans. Nucl. Sci. 2009, 55, 2992–2999. [Google Scholar] [CrossRef]
  48. Tafrishi, H.; Sadeghzadeh, S.; Ahmadi, R.; Molaei, F.; Hassanloo, H. Investigation of tetracosane thermal transport in presence of graphene and carbon nanotube fillers—A molecular dynamics study. J. Energy Storage 2020, 29, 101321. [Google Scholar] [CrossRef]
Materials 16 04224 g001 550
Figure 1. Temperature and potential energy of the simulation system during the relaxation stage as time functions.
Figure 1. Temperature and potential energy of the simulation system during the relaxation stage as time functions.
Materials 16 04224 g001
Materials 16 04224 g002 550
Figure 2. A schematic diagram of the two perpendicular constant-temperature walls: (a) the yellow area represents the constant-temperature walls; (b) constant-temperature walls under periodic boundary conditions.
Figure 2. A schematic diagram of the two perpendicular constant-temperature walls: (a) the yellow area represents the constant-temperature walls; (b) constant-temperature walls under periodic boundary conditions.
Materials 16 04224 g002
Materials 16 04224 g003 550
Figure 3. The relationship between the number of defects and time under the irradiation of α-particles with different energies: (a) 0.1 MeV; (b) 0.5 MeV.
Figure 3. The relationship between the number of defects and time under the irradiation of α-particles with different energies: (a) 0.1 MeV; (b) 0.5 MeV.
Materials 16 04224 g003
Materials 16 04224 g004 550
Figure 4. Defect distribution maps (in the x-direction): (a) peak stage at 0.1 MeV; (b) final stage at 0.1 MeV; (c) peak stage at 0.5 MeV; (d) final stage at 0.5. The green ball is the vacancy defect and the red ball is the interstitial defect.
Figure 4. Defect distribution maps (in the x-direction): (a) peak stage at 0.1 MeV; (b) final stage at 0.1 MeV; (c) peak stage at 0.5 MeV; (d) final stage at 0.5. The green ball is the vacancy defect and the red ball is the interstitial defect.
Materials 16 04224 g004
Materials 16 04224 g005 550
Figure 5. The size and distribution of defect clusters: (a) Number and location of large, medium, and small defect clusters; (b) Number and location of amorphous regions.
Figure 5. The size and distribution of defect clusters: (a) Number and location of large, medium, and small defect clusters; (b) Number and location of amorphous regions.
Materials 16 04224 g005
Materials 16 04224 g006 550
Figure 6. Radial distribution functions of GaN material under different dose rates: (a) accumulation of five injected α-particles; (b) accumulation of ten injected α-particles.
Figure 6. Radial distribution functions of GaN material under different dose rates: (a) accumulation of five injected α-particles; (b) accumulation of ten injected α-particles.
Materials 16 04224 g006
Materials 16 04224 g007 550
Figure 7. IDS surface topography of GaN materials in the z-direction under different injection rates: (a) accumulated injection of five α-particles; (b) accumulated injection of ten α-particles.
Figure 7. IDS surface topography of GaN materials in the z-direction under different injection rates: (a) accumulated injection of five α-particles; (b) accumulated injection of ten α-particles.
Materials 16 04224 g007
Materials 16 04224 g008 550
Figure 8. Size and distribution of defect clusters under different injection amounts: (a) accumulated injection of five α-particles; (b) accumulated injection of ten α-particles.
Figure 8. Size and distribution of defect clusters under different injection amounts: (a) accumulated injection of five α-particles; (b) accumulated injection of ten α-particles.
Materials 16 04224 g008
Materials 16 04224 g009 550
Figure 9. Distribution of defects under different injection amounts: (a) accumulated injection of five α-particles; (b) accumulated injection of ten α-particles. The green ball is the vacancy defect and the red ball is the interstitial defect.
Figure 9. Distribution of defects under different injection amounts: (a) accumulated injection of five α-particles; (b) accumulated injection of ten α-particles. The green ball is the vacancy defect and the red ball is the interstitial defect.
Materials 16 04224 g009
Table
Table 1. Set specific parameters.
Table 1. Set specific parameters.
ParameterValueParameterValue
Lattice constanta = b = 3.216 Å
c = 5.240 Å		Atomic number500,000
Defect rate10%Doping rate5 × 1018 cm−3
Crystal size50a × 50a × 50aTemperature300 K
Incident energy0.1 MeV
0.5 MeV		Injected dose4 × 1011 ions/cm2
2 × 1012 ions/cm2
4 × 1012 ions/cm2
Table
Table 2. Classification and quantity statistics of defects.
Table 2. Classification and quantity statistics of defects.
0.1 MeV0.5 MeV
Defect ClassificationQuantityRepair Capability (%)Recombination Efficiency (%)QuantityRepair Capability (%)Recombination Efficiency (%)
VacancyGa vacancy16641.3431.31127622.9526.5
N vacancy21721.38129329.8
Interstitial defectTwo atomsGa-N (Ga)19435.2533.52111626.8626.92
Ga-N (N)1201075
Ga-Ga (Ga)2292
Ga-Ga (N)28
N-N (Ga)013
N-N (N)2090
Three atomsN-N-N (N)018.19128.57
Ga-Ga-N (Ga)224
Ga-Ga-N (N)19
Ga-N-N (Ga)322
Ga-N-N (N)314
Four atomsGa-Ga-N-N (N)0/1/
Ga-Ga-Ga-N (Ga)01
Table
Table 3. Proportions of different components in IDS damage analysis.
Table 3. Proportions of different components in IDS damage analysis.
StructureInitial State (%)Five α-Particles (%)Ten α-Particles (%)
Other16.327.331.9
Hexagonal diamond18.114.012.8
Hexagonal diamond (1st neighbor)29.624.723.1
Hexagonal diamond (2nd neighbor)36.134.032.2
Table
Table 4. Sputtering yield under different simulation conditions.
Table 4. Sputtering yield under different simulation conditions.
ModeNumber of Sputtered AtomsSputtering Yield (%)
Energy of 0.1 MeV260.0058
Energy of 0.5 MeV330.0073
Accumulation of five α-particles23020.512
Accumulation of ten α-particles28710.638
Table
Table 5. Defect classification and quantity statistics under different cumulative irradiations.
Table 5. Defect classification and quantity statistics under different cumulative irradiations.
Five α-ParticlesTen α-Particles
Defect TypeQuantityQuantity
VacancyGa vacancy566912,093
N vacancy11,47316,780
InterstitialTwo atomsGa-N (Ga)897412,803
Ga-N (N)33838144
Ga-Ga (Ga)323712
Ga-Ga (N)254470
N-N (Ga)234426
N-N (N)153399
Three atomsN-N-N (N)916
Ga-Ga-N (Ga)384623
Ga-Ga-N (N)94159
Ga-N-N (Ga)125291
Ga-N-N (N)101313
Four atoms116256
Five atoms or more232324
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, Y.; Xiong, Z.; Ouyang, X. Molecular Dynamics Study on the Mechanism of Gallium Nitride Radiation Damage by Alpha Particles. Materials 2023, 16, 4224. https://doi.org/10.3390/ma16124224

AMA Style

Liu Y, Xiong Z, Ouyang X. Molecular Dynamics Study on the Mechanism of Gallium Nitride Radiation Damage by Alpha Particles. Materials. 2023; 16(12):4224. https://doi.org/10.3390/ma16124224

Chicago/Turabian Style

Liu, Yang, Zhenpeng Xiong, and Xiaoping Ouyang. 2023. "Molecular Dynamics Study on the Mechanism of Gallium Nitride Radiation Damage by Alpha Particles" Materials 16, no. 12: 4224. https://doi.org/10.3390/ma16124224

APA Style

Liu, Y., Xiong, Z., & Ouyang, X. (2023). Molecular Dynamics Study on the Mechanism of Gallium Nitride Radiation Damage by Alpha Particles. Materials, 16(12), 4224. https://doi.org/10.3390/ma16124224

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop