Next Article in Journal
Broadband Continuous Transverse Stub (CTS) Array Antenna for High-Power Applications
Previous Article in Journal
Controllable Carrier Doping in Two-Dimensional Materials Using Electron-Beam Irradiation and Scalable Oxide Dielectrics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Temperature and Incidence Angle on the Irradiation Cascade Effect of 6H-SiC: Molecular Dynamics Simulations

Key Laboratory for Wide Band Gap Semiconductor Materials and Devices of Education, School of Microelectronics, Xidian University, Xi’an 710071, China
*
Author to whom correspondence should be addressed.
Micromachines 2023, 14(11), 2126; https://doi.org/10.3390/mi14112126
Submission received: 25 October 2023 / Revised: 16 November 2023 / Accepted: 17 November 2023 / Published: 19 November 2023
(This article belongs to the Section D1: Semiconductor Devices)

Abstract

:
SiC devices have been typically subjected to extreme environments and complex stresses during operation, such as intense radiation and large diurnal amplitude differences on the lunar surface. Radiation displacement damage may lead to degradation or failure of the performance of semiconductor devices. In this paper, the effects of temperature and incidence angle on the irradiation cascade effect of 6H-SiC were investigated separately using the principles of molecular dynamics. Temperatures were set to 100 K, 150 K, 200 K, 250 K, 300 K, 350 K, 400 K and 450 K. The incidence direction was parallel to the specified crystal plane, with angles of 8°, 15°, 30°, 45°, 60° and 75° to the negative direction of the Z-axis. In this paper, the six types of defects were counted, and the microscopic distribution images and trajectories of each type of defect were extracted. The results show a linear relationship between the peak of the Frenkel pair and temperature. The recombination rate of Frenkel pairs depends on the local temperature and degree of aggregation at the center of the cascade collision. Increasing the angle of incidence first inhibits and then promotes the production of total defects and Frenkel pairs. The lowest number of total defects, Frenkel pairs and antisite defects are produced at a 45° incident angle. At an incidence angle of 75°, larger size hollow clusters and anti-clusters are more likely to appear in the 6H-SiC.

1. Introduction

SiC material is a typical third-generation semiconductor material. It has been used in space exploration, nuclear reactions and electrical equipment at high altitudes on plateaus. 3C-SiC is mainly used in photovoltaic applications [1]. 4H-SiC materials are mainly used to fabricate high-frequency and high-power devices [2]. 6H-SiC materials are mainly used for the manufacture of photoconductive detectors [3]. Researchers typically use Rutherford Backscattering Spectroscopy (RBS) and Transmission Electron Microscope (TEM) to probe the irradiation damage of SiC materials [4]. Wenhui Dong et al. studied different doses of Si, H ion injection experiments using Scanning Electron Microscopy (SEM) and Transmission Electron Microscopy. The results show that the irradiation damage effect of 6H-SiC is highly dependent on the type of irradiated ions, the amount of ion injection and the irradiation sequence [5,6]. W. J. Weber investigated defect clusters and amorphization induced by high-dose neutron irradiation in SiC materials [7]. The generation and recrystallization behavior of SiC amorphous crystals has been achieved by low-energy He and Fe ion injection [8,9]. C. Hemmingsson et al. conducted an investigation into the deep-level defects induced by electron irradiation in 4H-SiC using deep-level transient spectroscopy. They observed multiple electron traps and a single hole trap [10]. Concurrently, A. Castaldini et al. identified that among the five defect levels introduced by irradiation, three exhibited similar electrical characteristics [11]. Katoh et al. investigated the composite materials of SiC fibers and pyrocarbon under neutron irradiation conditions. They found that the density of defects significantly influences the post-neutron irradiation electrical properties of SiC semiconductors [12]. Wang Pengfei et al. found that the type and structure of irradiation-induced defects in 6H-SiC materials affect the extent of neutron irradiation damage. They observed that annealing temperature affects the recombination of defects and lattice recovery, ultimately influencing the resistivity of 6H-SiC materials [13]. Other researchers also found that the radiation damage caused by P ion injection partially recovered or even disappeared at different annealing temperatures and annealing times [14]. In addition, several researchers have investigated the displacement threshold energies and the laws of defect aggregation in SiC lattice atoms [15,16].
Silicon vacancies in silicon carbide (SiC) have been proposed as attractive candidates for quantum technology applications, such as quantum sensing and quantum repeaters [17]. As the use of 6H-SiC becomes more widespread, the complex environment poses new challenges for the application of 6H-SiC materials. Further understanding of the defect formation mechanism in 6H-SiC will help broaden the application fields of 6H-SiC. For example, in SiC device manufacturing technology, proton irradiation can be used to create localized regions of high resistance in a semiconductor [18]. Based on the principles of molecular dynamics (MD), this paper simulated the cascade effect of 6H-SiC at different temperatures and incidence directions. The evolution of point defects and clusters in 6H-SiC materials was studied under specific stress conditions.

2. Methods

In this paper, irradiation experiments on 6H-SiC targets are modelled and simulated using the Large-scale Atomic\Molecular Massive Parallel Simulator (LAMMPS) [19]. 6H-SiC is a hexagonal crystal. The lattice constants are a = b = 3.08101 Å and c = 15.1248 Å. Table 1 shows the atomic coordinates of the unit cell of 6H-SiC [20]. Figure 1a,b show the unit cell and the supercell of the model for 6H-SiC, respectively. The supercell is obtained by expanding the unit cell.
When low-energy particles irradiate materials, the main mode of energy loss is elastic collisions between atoms. Displacement cascades from cascade collisions are the main form of displacement damage. Figure 2 shows an irradiation model with fully periodic boundary conditions. Periodic boundary conditions ensure the integrity of the forces on the simulated system boundary atoms. The blue area in the model is a 1.5 nm thick constant temperature area. The temperature of the thermostatic area is controlled by the velocity rescaling (temp/rescale) method. The method is used to simulate the dissipation of energy. Collision area, the red area of the model, can be used to simulate the cascade collision effect of 6H-SiC. The Si-PKA is set in the inner part of the collision area, and the place is 2 nm from the thermostatic area on the upper surface. The aim is to avoid displaced atoms from cascade collisions entering the thermostat area. The incident direction of the PKA is parallel to the x-z plane and at an angle of 8° along the negative direction of the Z-axis.
The potential function describes the potential energy of the atoms in a material as a function of the positions of the atoms, i.e., it describes the forces between the atoms. The accuracy of the potential function parameters determines the accuracy of the simulation results. The potential function used in this paper is the Tersoff/ZBL potential function, whose parameters were confirmed in the experiments of Erhart et al. [21,22]. According to molecular dynamics, the atoms in the system obey the classical mechanics. The velocity of motion of each atom can be derived from the definition of velocity, i.e., the atoms in the system satisfy the set of Equations (1). As shown in Table 2, the simulation system was first relaxed by constant-pressure, constant-temperature ensemble (NPT) for 30 ps. And the micro-canonical ensemble (NVE) was used during the irradiation process. Table 3 shows the experimental design. Experiment 1 parameters are kept constant except temperature. Experiment 2 parameters are kept constant except for the angle of incidence.
F i ( t ) = m i ν i ( t ) t ν i = r i ( t ) t .
F i ( t ) is the force on particle i at time t . m i , ν i and r i are the mass, velocity and coordinates of particle i, respectively.
Figure 3 is a schematic representation of the Wiger–Seitz method. The Open Visualisation Tool (OVITO) is used to identify and count the types and clusters of defects [23]. Frenkel pairs, silicon vacancy (VSi), carbon vacancy (VC), silicon interstitial (ISi), carbon interstitial (IC), silicon antisite (SiC), carbon antisite (CSi) and cluster can be shown in the OVITO. Among them, clusters have a search radius of 0.22 nm [16,24].

3. Results and Discussion

3.1. Validation of the Model

The density, radial distribution function, atomic angular distribution function and coordination number of the model are calculated in this section. In addition, to verify the reliability of the model, the experimental data from other researchers are compared in the following content.
Table 4 shows that the lattice constants after NPT optimization are not significantly different from the experimental values. Figure 4a shows that the bond length of C-Si is 1.8645 Å in the simulated system. Figure 4b shows that the coordination number of the 6H-SiC crystal is 12, consistent with the coordination number of the actual hexagonal crystal structure. Figure 4c shows angles of 109.5° and 109.2° for C-Si-C and Si-C-Si, respectively.
Table 5 shows the error statistics information. Compared with the experimental value [20], the density error of the model in this paper is about 3.302%. The error of the Si-C bond length is about −1.3388%. The atomic angle errors of Si-C-Si and C-Si-C are about −0.2439% and 0.0759%, respectively.

3.2. Effect of Temperature on 6H-SiC Cascade Collisions

In the field of space exploration, extreme environments pose a severe impact on equipment. For example, the temperature on the lunar surface can reach −183 °C or 127 °C. This subsection focuses on the variation in the 6H-SiC irradiation cascade effect at ambient temperatures of 100 K, 150 K, 200 K, 250 K, 300 K, 350 K, 400 K and 450 K.
Figure 5a shows that different temperatures have no effect on the overall trend of defect evolution in 6H-SiC. The peak formation times of the Frenkel pairs are basically the same. The duration of cascade collisions becomes longer with increasing temperature. Figure 5b shows that the fitted curve between the peak of the Frenkel pair and the temperature satisfies a linear relationship, which can be expressed by Equation (2). The reason is that the transition from low to high temperatures enhances the thermal motion of the atoms. The energy of the lattice atoms is more likely to meet the displacement threshold energy. The permanent Frenkel pairs and their recombination rates have no obvious regularity. The paper further explains this phenomenon through the temperature nephogram.
N m a x = a + k * T ,
where N m a x is the peak of the Frenkel pair, a is the intercept ( a = 111.94048 ± 2.99957 ). k is the slope ( k = 0.00976 ± 0.01007 ), and T is the temperature in Kelvin.
Figure 6 shows the nephogram at different temperatures. A peak represents a cascade collision centre, corresponding to a recombination centre. As the temperature increases from 100 to 200 K, the number of recombination centres increases and the dispersion decreases. The result is an increase in the recombination rate of Frenkel pairs. As the temperature rises from 200 K to 350 K, the dispersion of the recombination centre increases, leading to a decrease in the recombination rate of Frenkel pairs. A larger cascade of atomic collisions occurs when the temperature is raised from 350 K to 450 K, corresponding to a local temperature of over 1000 K. Further, high temperature annealing promotes the recombination of defects [14].
The barrier potential for various types of point defects in SiC is given in Table 6 [25,26,27]. Figure 7 shows that the irradiation defect types are mainly in the form of point defects at different temperatures. These clusters are mainly vacancy clusters and antisite clusters. The reason is that IC and ISi have relatively low migration potential. It is not easy to form interstitial clusters.

3.3. Effect of Incidence Angle on 6H-SiC Cascade Collisions

The incident angle of irradiated ions affects the particle deposition trajectory and hence the transient current of the device, a phenomenon that has been reported [28,29]. The 6H-SiC is an anisotropic crystal. Irradiated particles with different incidence angles may form different cascade collisions in 6H-SiC, leading to different levels of irradiation damage. Therefore, this section focuses on the effect of different incidence angles on 6H-SiC cascade collisions.
Figure 8a shows that the number of permanent defects produced by cascade collisions tends to decrease and then increase as the angle of incidence increases. When the incidence angles are 8° and 15°, the total number of defects is roughly the same. The number of total defects, SiC and CSi, are the least at the 45° incident angle. Among them, the total defect consists of VC, VSi, IC, ISi, SiC and CSi. Figure 8b shows a similar regularity for the number of permanent Frenkel pairs generated by cascade collisions. Figure 8c indicates that the number of permanent CSi tends to decrease and then increase as the angle of incidence increases. When the incidence angle is 15° and 30°, the difference in CSi quantity is very small. Figure 8d shows that there is no obvious regularity in the number of SiC. To further explain the differences in the number of defects in each category, we provide reactions between the various types of point defects studied and their reaction barrier (ERe) and reaction distances (r) [27,30]. In Table 7, the reaction barriers for CC and SiSi are smaller than those of the other forms. Therefore, during the defective recombination phase, a large number of recombination occurs between VC and IC and between VSi and ISi. This leads to a large reduction in Frenkel pairs as shown in Figure 8b. However, by comparison with Figure 8b,c, it can be observed that the amount of CSi does not change significantly during the defect recombination phase. The recombination of Si atoms in SiC and VSi is the reason for the decrease in the amount of SiC in Figure 8d. This also suggests that the SiC produced at the beginning of irradiation is not stable SiC, but a combination of VC and ISi. As shown in Table 8 [26,27], the formation energy of SiC is lower than that of CSi, so the peak defects in Figure 8c are lower than those in Figure 8d. The results in Figure 8e,f show that the variation in the incidence angle has no significant effect on the final percentage of vacancy defects and antisite defects. This shows that the effect of the variation in the PKA incidence angle on the electrical properties of the electronic devices is more in the regional distribution of defect clusters of vacancy defects and interstitial atomic defects generated within the material.
Figure 9 shows the microscopic distribution of various defect types. The trajectories in the figure are produced by atoms with displacements greater than 4 Å. Defects are distributed in the shape of grape bunches near the PKA trajectory. As the angle of incidence increases, local cascade collisions become more intense, and the aggregation of defects becomes more pronounced at the end of the PKA trajectory. Consistent with the above results, the local amorphization is more obvious within the 6H-SiC material. A distinct region characterized by the aggregation of vacancy defects and interstitial atomic defects emerges irrespective of alterations in the incident angle. Furthermore, this region undergoes modifications as the incident angle varies.
Figure 10 shows that vacancy clusters and antisite clusters with an atomic number of two are the most common. The largest size of vacancy clusters and antisite clusters appear at a 75° incident angle. At an incident angle of 60°, there are no antisite clusters and no vacancy clusters greater than two. However, it can be concluded from Table 6 that the migration potentials of vacancy and antisite defects are much higher than those of interstitial defects, and therefore the clusters formed are mainly vacancy and antisite clusters.

4. Conclusions

During irradiation of 6H-SiC crystals, there is a monotonically increasing linear relationship between the peak of the Frenkel pair and temperature. The peak of the Frenkel pair increases with increasing temperature, but the amount of change is small. The recombination rate of Frenkel pairs depends on the local temperature and the degree of defective recombination center aggregation. As the temperature increases from 300 K to 450 K, the cascade collision distribution tends to concentrate, and local high temperatures are more likely to occur. As the angle of incidence increases, the number of both total defects and Frenkel pairs exhibits a trend of first decreases and then increases. Among them, the number of defects is minimal at a 45° incident angle. At a 75° incident angle, larger vacancy clusters and antisite clusters are produced. As the angle of incidence increases, secondary trajectories from irradiation cascade effect converge towards the end of the PKA trajectory.

Author Contributions

Conceptualization, Y.C.; Data curation, Y.C., C.Y. and H.W.; Formal analysis, C.Y. and H.W.; Funding acquisition, H.L.; Investigation, C.Y. and H.W.; Methodology, Y.C.; Project administration, H.L.; Resources, H.L.; Software, Y.C.; Supervision, H.L.; Validation, Y.C.; Visualization, H.L.; Writing—original draft, Y.C.; Writing—review and editing, C.Y. and H.W. All authors have read and agreed to the published version of the manuscript.

Funding

National Natural Science Foundation of China (Grant No. U2241221).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sameera, J.N.; Islam, M.A.; Islam, S.; Hossain, T.; Sobayel, M.; Akhtaruzzaman; Amin, N.; Rashid, M.J. Cubic silicon carbide (3C–SiC) as a buffer layer for high efficiency and highly stable CdTe solar cell. Opt. Mater. 2022, 123, 111911. [Google Scholar] [CrossRef]
  2. Yang, F.; Tang, X.; Wei, X.; Sang, L.; Liu, R.; Bai, S.; Peng, T.; Zhao, G.; Yang, P.; Yang, T.; et al. Demonstrations of high voltage SiC materials, devices and applications in the solid state transformer. J. Cryst. Growth 2023, 604, 127059. [Google Scholar] [CrossRef]
  3. Liu, Q.; Zhou, D.; Cai, X.; Qi, M.; Xu, W.; Chen, D.; Ren, F.; Zhang, R.; Zheng, Y.; Lu, H. Effect of very high-fluence proton radiation on 6H-SiC photoconductive proton detectors. IEEE Electron. Device Lett. 2019, 40, 1929–1932. [Google Scholar] [CrossRef]
  4. Gloginjić, M.; Erich, M.; Kokkoris, M.; Liarokapis, E.; Fazinić, S.; Karlušić, M.; Luketić, K.T.; Petrović, S. The quantitative 6H-SiC crystal damage depth profiling. J. Nucl. Mater. 2021, 555, 153143. [Google Scholar] [CrossRef]
  5. Wang, X.; Ji, Y.; Zhang, M.; Zhao, Y.; Chen, Y.; Zhao, Z.; Pan, S.; Wang, H. Damage effects in 6H-SiC single crystals by Si&H dual ion irradiation: A combined Raman and XRD study. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 2020, 485, 20–25. [Google Scholar]
  6. Dong, W.; Shen, Q.; Wei, M.; Lei, P.; Song, L.; Chang, Q.; Ye, C. Research on the surface damage of Si+ and H+ co-implanted 6H-SiC before and after annealing. Nucl. Instrum. Methods Phys. Res. Sect. B Beam Interact. Mater. At. 2023, 538, 81–86. [Google Scholar] [CrossRef]
  7. Weber, W.J.; Gao, F. Irradiation-induced defect clustering and amorphization in silicon carbide. J. Mater. Res. 2010, 25, 2349–2353. [Google Scholar] [CrossRef]
  8. Li, B.; Liu, H.; Shen, T.; Xu, L.; Wang, J.; Zhao, F.; Peng, D.; Li, J.; Sheng, Y.; Xiong, A. Irradiation-induced microstructure damage in He-irradiated 3C-SiC at 1000 °C. J. Eur. Ceram. Soc. 2020, 40, 1014–1022. [Google Scholar] [CrossRef]
  9. Daghbouj, N.; Li, B.; Callisti, M.; Sen, H.; Karlik, M.; Polcar, T. Microstructural evolution of helium-irradiated 6H–SiC subjected to different irradiation conditions and annealing temperatures: A multiple characterization study. Acta Mater. 2019, 181, 160–172. [Google Scholar] [CrossRef]
  10. Hemmingsson, C.; Son, N.T.; Kordina, O.; Bergman, J.P.; Janzén, E.; Lindström, J.L.; Savage, S.; Nordell, N. Deep level defects in electron-irradiated 4H SiC epitaxial layers. J. Appl. Phys. 1997, 81, 6155–6159. [Google Scholar] [CrossRef]
  11. Castaldini, A.; Cavallini, A.; Rigutti, L.; Nava, F.; Ferrero, S.; Giorgis, F. Deep levels by proton and electron irradiation in 4H–SiC. J. Appl. Phys. 2005, 98, 053706. [Google Scholar] [CrossRef]
  12. Konopka, A.; Greulich-Weber, S.; Dierolf, V.; Jiang, H.; Gerstmann, U.; Rauls, E.; Sanna, S.; Schmidt, W. Microscopic structure and energy transfer of vacancy-related defect pairs with Erbium in wide-gap semiconductors. Opt. Mater. 2011, 33, 1041–1044. [Google Scholar] [CrossRef]
  13. Wang, P.F.; Ruan, Y.F.; Huang, L.; Zhu, W. Nitrogen-promoted formation of graphite-like aggregations in SiC during neutron irradiation. J. Appl. Phys. 2012, 111, 063517. [Google Scholar] [CrossRef]
  14. Gu, J.J.; Zhao, J.H.; Bu, M.Y.; Wang, S.M.; Fan, L.; Huang, Q.; Li, S.; Yue, Q.-Y.; Wang, X.-L.; Wei, Z.-X.; et al. Study on the damage evolution of 6H-SiC under different phosphorus ion implantation conditions and annealing temperatures. Results Phys. 2022, 43, 106127. [Google Scholar] [CrossRef]
  15. Devanathan, R.; De La Rubia, T.D.; Weber, W.J. Displacement threshold energies in β-SiC. J. Nucl. Mater. 1998, 253, 47–52. [Google Scholar] [CrossRef]
  16. 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]
  17. Singh, H.; Anisimov, A.N.; Baranov, P.G.; Suter, D. Identification of different silicon vacancy centers in 6H-SiC. arXiv 2022, arXiv:2212.10256. [Google Scholar]
  18. Lebedev, A.A.; Veinger, A.I.; Davydov, D.V.; Kozlovski, V.V.; Savkina, N.S.; Strel’chuk, A.M. Doping of n-type 6H–SiC and 4H–SiC with defects created with a proton beam. J. Appl. Phys. 2000, 88, 6265–6271. [Google Scholar] [CrossRef]
  19. Thompson, A.P.; Aktulga, H.M.; Berger, R.; Bolintineanu, D.S.; Brown, W.M.; Crozier, P.S.; in ‘tVeld, P.J.; Kohlmeyer, A.; Moore, S.G.; Nguyen, T.D.; et al. LAMMPS—A flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comput. Phys. Commun. 2022, 271, 108171. [Google Scholar]
  20. Capitani, G.C.; Di Pierro, S.; Tempesta, G. The 6 H-SiC structure model: Further refinement from SCXRD data from a terrestrial moissanite. Am. Mineral. 2007, 92, 403–407. [Google Scholar] [CrossRef]
  21. Erhart, P.; Albe, K. Analytical potential for atomistic simulations of silicon, carbon, and silicon carbide. Phys. Rev. B 2005, 71, 035211. [Google Scholar]
  22. Tersoff, J. Modeling solid-state chemistry: Interatomic potentials for multicomponent systems. Phys. Rev. B 1989, 39, 5566. [Google Scholar]
  23. Stukowski, A. Visualization and analysis of atomistic simulation data with OVITO—The Open Visualization Tool. Model. Simul. Mater. Sci. Eng. 2009, 18, 015012. [Google Scholar] [CrossRef]
  24. Jortner, J. Cluster size effects. Z. Für Phys. D At. Mol. Clust. 1992, 24, 247–275. [Google Scholar]
  25. Shrader, D.; Khalil, S.M.; Gerczak, T.; Allen, T.R.; Heim, A.J.; Szlufarska, I.; Morgan, D. Ag diffusion in cubic silicon carbide. J. Nucl. Mater. 2011, 408, 257–271. [Google Scholar]
  26. Kondo, S.; Katoh, Y.; Snead, L.L. Analysis of grain boundary sinks and interstitial diffusion in neutron-irradiated SiC. Phys. Rev. B 2011, 83, 075202. [Google Scholar] [CrossRef]
  27. Zheng, M.J.; Swaminathan, N.; Morgan, D.; Szlufarska, I. Energy barriers for point-defect reactions in 3 C-SiC. Phys. Rev. B 2013, 88, 054105. [Google Scholar] [CrossRef]
  28. Roche, P.; Autran, J.L.; Gasiot, G. Technology downscaling worsening radiation effects in bulk: SOI to the rescue. In Proceedings of the 2013 IEEE International Electron Devices Meeting, Washington, DC, USA, 9–11 December 2013. [Google Scholar]
  29. Bartra, W.C.; Vladimirescu, A.; Reis, R. Bulk and FDSOI Sub-micron CMOS transistors resilience to single-event transients. In Proceedings of the 2015 IEEE International Conference on Electronics, Circuits, and Systems (ICECS), Cairo, Egypt, 6–9 December 2015. [Google Scholar]
  30. Roma, G.; Crocombette, J.-P. Evidence for a kinetic bias towards antisite formation in SiC nano-decomposition. J. Nucl. Mater. 2010, 403, 32–41. [Google Scholar]
Figure 1. (a) Unit cell model for 6H-SiC; (b) Supercell model for 6H-SiC.
Figure 1. (a) Unit cell model for 6H-SiC; (b) Supercell model for 6H-SiC.
Micromachines 14 02126 g001
Figure 2. Schematic diagram of the irradiation model.
Figure 2. Schematic diagram of the irradiation model.
Micromachines 14 02126 g002
Figure 3. Schematic diagram of the Wiger–Seitz method.
Figure 3. Schematic diagram of the Wiger–Seitz method.
Micromachines 14 02126 g003
Figure 4. Structural data for the models in this paper. (a) Radial distribution function; (b) coordination number; (c) atomic angle distribution function.
Figure 4. Structural data for the models in this paper. (a) Radial distribution function; (b) coordination number; (c) atomic angle distribution function.
Micromachines 14 02126 g004
Figure 5. Production of a Frenkel pair at different temperatures. (a) Number of Frenkel pairs versus time; (b) Peak number of Frenkel pairs and permanent defects at different temperatures.
Figure 5. Production of a Frenkel pair at different temperatures. (a) Number of Frenkel pairs versus time; (b) Peak number of Frenkel pairs and permanent defects at different temperatures.
Micromachines 14 02126 g005
Figure 6. The nephogram at different temperatures. (a) 100 K; (b) 150 K; (c) 200 K; (d) 250 K; (e) 300 K; (f) 350 K; (g) 400 K; (h) 450 K.
Figure 6. The nephogram at different temperatures. (a) 100 K; (b) 150 K; (c) 200 K; (d) 250 K; (e) 300 K; (f) 350 K; (g) 400 K; (h) 450 K.
Micromachines 14 02126 g006
Figure 7. The number of clusters by size. (a) Vacancy clusters; (b) Interstice clusters; (c) Antisite clusters.
Figure 7. The number of clusters by size. (a) Vacancy clusters; (b) Interstice clusters; (c) Antisite clusters.
Micromachines 14 02126 g007
Figure 8. (a) Number of total defects versus time curves; (b) Number of Frenkel pairs versus time curves; (c) Number of CSi versus time curves; (d) Number of SiC versus time curves; (e) Percentage of vacancy; (f) Percentage of antisite.
Figure 8. (a) Number of total defects versus time curves; (b) Number of Frenkel pairs versus time curves; (c) Number of CSi versus time curves; (d) Number of SiC versus time curves; (e) Percentage of vacancy; (f) Percentage of antisite.
Micromachines 14 02126 g008
Figure 9. Microscopic distribution of point defects. (a) 8°; (b) 15°; (c) 30°; (d) 45°; (e) 60°; (f) 75°.
Figure 9. Microscopic distribution of point defects. (a) 8°; (b) 15°; (c) 30°; (d) 45°; (e) 60°; (f) 75°.
Micromachines 14 02126 g009
Figure 10. The number of clusters by size. (a) Vacancy clusters; (b) Antisite clusters.
Figure 10. The number of clusters by size. (a) Vacancy clusters; (b) Antisite clusters.
Micromachines 14 02126 g010
Table 1. Atomic coordinates for the 6H-SiC moissanite [15].
Table 1. Atomic coordinates for the 6H-SiC moissanite [15].
Atomxyz
Si 1000
Si 20.33330.66670.1664
Si 30.66670.33330.3329
Si 40.33330.66670.0412
Si 50.66670.33330.2080
Si 6000.3746
Table 2. Parameters related to the time step.
Table 2. Parameters related to the time step.
StageTime Step (fs)Number of MD StepsTime (ps)
Relaxation phase (NPT)130,00030
Cascade collision phase (NVE)0.0150,0000.5
Defect compounding stage (NVE)0.120,0002
Steady-state phase (NVE)120,00020
Table 3. The design of the experiment.
Table 3. The design of the experiment.
Experimental Group 1Temperature
(K)
Angle of Incidence
(°)
Velocity Component of the y-Axis 2 (Å/ps)Velocity Component of the z-Axis (Å/ps)
Experiment 1: Different ambient temperatures100, 150, 200, 250,
300, 350, 400, 450
257.8014−1834.3527
Experiment 2: Different angles of incidence300257.8014−1834.3527
30015°479.4312−1789.2616
30030°926.18997−1604.2081
30045°1309.8304−1309.8304
30060°1604.2081−926.1900
30075°1789.2616−479.4312
1 The dimensions of the model are 65a × 38b × 13c. Total atomic number is 770,640. 2 The velocity component of the x-axis is 0.
Table 4. Lattice constants.
Table 4. Lattice constants.
Lattice Constants 1
(Å)
This Work
(before NPT)
This Work
(after NPT)
Reference [15]
a3.08103.0809993.0810
b5.33645.3364005.3364
c15.124815.12479915.1248
1 The hexagonal lattice constants in the table were all converted to orthogonal lattice constants.
Table 5. Error statistics table.
Table 5. Error statistics table.
This Work 1Reference [15]Average Error 2 (%)
Density (g/mol)3.3163.213.302
Bond lengths (Å)1.86451.8898−1.3388
Angles (°)109.2 (Si-C-Si)109.467 (Si-C-Si)−0.2439
109.5 (C-Si-C)109.417 (C-Si-C)0.0759
1 The temperature is 300 K. 2 ± indicates that the value in this paper is greater or less than the experimental value.
Table 6. Potential barriers for various types of point defects.
Table 6. Potential barriers for various types of point defects.
VCVSiICISiSiCCSi
Migration potential/eV3.663.200.671.4811.6011.70
Table 7. The recombination pattern of point defects, reaction barriers (ERe) and reaction distances (r) [27,30].
Table 7. The recombination pattern of point defects, reaction barriers (ERe) and reaction distances (r) [27,30].
Recombination PatternReaction Barriers (ERe)/eVReaction Distances (r)/Å
IC + VC → CC0.433.08
ISi + VSi → SiSi0.175.34
IC + VSi → CSi1.253.30
ISi + VC → SiC1.113.70
IC + SiC → CC + ISi1.344.36
ISi + CSi → SiSi + IC0.644.36
Table 8. The formation energy of defects [26,27].
Table 8. The formation energy of defects [26,27].
VCVSiICISiSiCCSi
Formation energy/eV4.194.976.958.753.564.03
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

Chen, Y.; Liu, H.; Yan, C.; Wei, H. Influence of Temperature and Incidence Angle on the Irradiation Cascade Effect of 6H-SiC: Molecular Dynamics Simulations. Micromachines 2023, 14, 2126. https://doi.org/10.3390/mi14112126

AMA Style

Chen Y, Liu H, Yan C, Wei H. Influence of Temperature and Incidence Angle on the Irradiation Cascade Effect of 6H-SiC: Molecular Dynamics Simulations. Micromachines. 2023; 14(11):2126. https://doi.org/10.3390/mi14112126

Chicago/Turabian Style

Chen, Yaolin, Hongxia Liu, Cong Yan, and Hao Wei. 2023. "Influence of Temperature and Incidence Angle on the Irradiation Cascade Effect of 6H-SiC: Molecular Dynamics Simulations" Micromachines 14, no. 11: 2126. https://doi.org/10.3390/mi14112126

APA Style

Chen, Y., Liu, H., Yan, C., & Wei, H. (2023). Influence of Temperature and Incidence Angle on the Irradiation Cascade Effect of 6H-SiC: Molecular Dynamics Simulations. Micromachines, 14(11), 2126. https://doi.org/10.3390/mi14112126

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