Studies of Defect Structure in Epitaxial AlN/GaN Films Grown on (111) 3C-SiC

Several aspects such as the growth relation between the layers of the GaN/AlN/SiC heterostructure, the consistency of the interfaces, and elemental diffusion are achieved by High Resolution Transmission Electron Microscopy (HR-TEM). In addition, the dislocation densities together with the defect correlation lengths are investigated via High-Resolution X-ray Diffraction (HR-XRD) and the characteristic positron diffusion length is achieved by Doppler Broadening Spectroscopy (DBS). Moreover, a comparative analysis with our previous work (i.e., GaN/AlN/Si and GaN/AlN/Al2O3) has been carried out. Within the epitaxial GaN layer defined by the relationship F4¯3m (111) 3C-SiC || P63mc (0002) AlN || P63mc (0002) GaN, the total dislocation density has been assessed as being 1.47 × 1010 cm−2. Compared with previously investigated heterostructures (on Si and Al2O3 substrates), the obtained dislocation correlation lengths (Le = 171 nm and Ls =288 nm) and the mean distance between two dislocations (rd = 82 nm) are higher. This reveals an improved crystal quality of the GaN with SiC as a growth template. In addition, the DBS measurements upheld the aforementioned results with a higher effective positron diffusion length LeffGaN2 = 75 ± 20 nm for the GaN layer.


Introduction
Gallium nitride (GaN) is an important wide band-gap semiconductor used in the fabrication of nitride based heterostructure devices for photodetectors [1], electronics [2,3], and light emitting diodes [4,5] and is considered as a perspective material for positron moderation [6,7]. In this sense, understanding GaN's defect structure becomes vital since any impurities, vacancy defects, or strains can influence its unique physical and optical characteristics [3,5,8].
Several growth techniques have recently been studied for the preparation of large-area GaN films but a major obstacle is still represented by the lack of a suitable material that could be used as a substrate for the nitrides while being compatible from both thermal and structural points of view [7,9]. The ideal solution would be the use of single-crystal GaN substrates but as of yet there is no viable method for producing such substrates with a low cost and on a large area, on account of the reduced nitrogen solubility and diffusion in liquid gallium [10]. Consequently, sapphire or silicon have been used as substrates for GaN epitaxy despite their respective limitations. For instance, although Si is used as a promising high-quality and low-cost substrate, it still presents a large lattice mismatch

Materials and Methods
The GaN/SiC sample used in this work with dimensions of 10 × 10 × 0.35 mm 3 was produced at the NTT Advanced Technology Corporation (Kanagawa, Japan) for highelectron-mobility transistors. As presented in our previous work [13,14], the wafers' exact growth conditions were not made available by the producer.
Three complementary techniques were applied for the characterization of the material in order to study the microstructure, the interfaces' characteristics, and the defects distribution. The samples were investigated via HR-TEM, HR-XRD, and DBS. The HR-TEM measurements were carried out under 200 kV with a Cs-corrected Titan Themis 200 microscope (Thermo Fisher Scientific, Waltham, MA, USA) designed with a high-brightness electron source. Elemental line profiling was performed using a Scanning Transmission Electron Microscopy (STEM) detector and a Super-X Energy Dispersive Spectroscopy (EDS) detector for Selected Area Electron Diffraction (SAED). EDS data processing was performed using ImageJ software (National Institutes of Health and the Laboratory for Optical and Computational Instrumentation, Madison, WI, USA) [23], while SingleCrystal ® (Oxford, UK) and CrystalMaker ® (Oxford, UK) were used for simulating the SAED data and crystal structures [24].
The HR-XRD data were collected using a high-resolution Rigaku SmartLab X-ray diffractometer (Neu-Isenburg, Germany) featuring a 9 kW rotating anode (CuK α = 1.5418 Å) and a high-resolution hybrid detector (HyPix-3000). The ω scans, taken at a fixed 2θ Bragg angle corresponding to selected (h k l) planes, were recorded in a double-axis configuration using Cross Beam Optics (CBO) and Ge (220) × 4 monochromator, resulting in an axial divergence of 0.003 • in the vertical diffraction plane of the goniometer. In order to prevent the samples' curvature from influencing the measurements, the incidence beam was limited by an incidence slit of 1 mm, while on the detector side 4 mm and 38.5 mm receiving slits were used (open detector configuration). The samples were aligned with the (111) plane of the SiC substrate prior to the measurement of the selected GaN planes in order to remove any asymmetry due to optics or sample misalignment. DBS measurements were performed at room temperature using a beam of incident positrons with controllable energy at the slow positron beam line of the High Energy Physics Institute in Beijing, China. The intensity of slow positrons obtained by solid neon moderation was~2 × 10 6 s −1 , forming a beam with a Ø 5.5 mm cross-section. A high purity Ge detector with an efficiency of 27% (model GEM20P4, ORTEC, Zoetermeer, The Netherlands) was located 20 cm away from the samples in a perpendicular orientation with respect to the axis of the positron beam. The approximate resolution of this detector was FWHM = 0.97 keV for the 511 keV line. Each spectrum was acquired over an 8 min cycle using fixed incident positron energy (E + ) in the range of 0.5 to 25 keV. The resulting statistics of individual spectra is~5 × 10 5 counts for the 511 keV peak region.

Results and Discussion
The evaluation of the TEM cross-section images ( Figure 1) reveals a buffer layer of aluminum nitride (AlN); therefore, the GaN/SiC sample can be regarded as a GaN/AlN/SiC heterostructure. Usually, such a layer aims at reducing the number of defects in the GaN film on account of the low lattice mismatch between AlN and GaN of 2.4% [25].
positrons with controllable energy at the slow positron beam line of the High Energy Physics Institute in Beijing, China. The intensity of slow positrons obtained by solid neon moderation was ~2 × 10 6 s −1 , forming a beam with a Ø 5.5 mm cross-section. A high purity Ge detector with an efficiency of 27% (model GEM20P4, ORTEC, Zoetermeer, the Netherlands) was located 20 cm away from the samples in a perpendicular orientation with respect to the axis of the positron beam. The approximate resolution of this detector was FWHM = 0.97 keV for the 511 keV line. Each spectrum was acquired over an 8 min cycle using fixed incident positron energy (E+) in the range of 0.5 to 25 keV. The resulting statistics of individual spectra is ~5 × 10 5 counts for the 511 keV peak region.

Results and Discussion
The evaluation of the TEM cross-section images ( Figure 1) reveals a buffer layer of aluminum nitride (AlN); therefore, the GaN/SiC sample can be regarded as a GaN/AlN/SiC heterostructure. Usually, such a layer aims at reducing the number of defects in the GaN film on account of the low lattice mismatch between AlN and GaN of 2.4% [25]. Figure 1a depicts HR-TEM micrographs and SAED patterns of the studied samples, revealing the arrangement of the SiC and AlN planes with respect to their interface. Below the interface, the planes of atoms parallel to the interface featured an interlayer distance of 2.52 Å , corresponding to the (111) lattice planes of 43 cubic SiC (ICDD 00-029-1129). Additionally, the (200) lattice planes of cubic SiC were highlighted, showing an interlayer distance of 2.18 Å . Above the interface, the planes of atoms parallel to it featured an interlayer distance of 2.49 Å , which corresponds to the (0002) lattice planes of 63 hexagonal AlN (ICDD 00-025-1133). Furthermore, interlayer distances of 2.37 Å were also measured, corresponding to the (2022) lattice planes of hexagonal AlN. The TEM analysis confirms epitaxial growth with the relationship 43 (111) SiC || 63 (0002) AlN, as can be deduced using the overlays of simulated crystals in Figure 1b.  In Figure 2a, HR-TEM and SAED reveal the structure of AlN and GaN below and above their interface. Below the interface, the planes of atoms parallel to it show an interlayer distance corresponding to the (0002) lattice planes of P63mc hexagonal AlN. In addition, the (2022) lattice planes of hexagonal AlN are also depicted with an interlayer distance of 2.37 Å. Above the interface, the planes of atoms parallel to it featured an interlayer distance of 2.59 Å, which corresponds to the (0002) lattice planes of P63mc hexagonal GaN (ICDD 00-050-0792). Finally, interlayer distances of 1.89 Å were also measured, corresponding to the (1012) lattice planes of hexagonal GaN. It could thus be deduced that through a coherent interface, (0002) hexagonal GaN had grown over the (0002) hexagonal AlN with the relationship P63mc AlN (0002) || P63mc (0002) GaN, as simulated in Figure 2b. In Figure 2a, HR-TEM and SAED reveal the structure of AlN and GaN below and above their interface. Below the interface, the planes of atoms parallel to it show an interlayer distance corresponding to the (0002) lattice planes of 63 hexagonal AlN. In addition, the (2022) lattice planes of hexagonal AlN are also depicted with an interlayer distance of 2.37 Å . Above the interface, the planes of atoms parallel to it featured an interlayer distance of 2.59 Å , which corresponds to the (0002) lattice planes of 63 hexagonal GaN (ICDD 00-050-0792). Finally, interlayer distances of 1.89 Å were also measured, corresponding to the (1012) lattice planes of hexagonal GaN. It could thus be deduced that through a coherent interface, (0002) hexagonal GaN had grown over the (0002)    For the assessment of the threading dislocation density, ρ d and the correlation length, L, several ω scans (rocking curves) were performed for selected (h k l) planes (see Figure 4). In order to determine the edge dislocation density, ρ e d and edge correlation length, L e , the (1015) GaN plane was used, while for the screw values ρ s d and L s , the (0004) plane was measured. For the assessment of the threading dislocation density, d and the correlation length, L, several ω scans (rocking curves) were performed for selected (h k l) planes (see Figure 4). In order to determine the edge dislocation density, d e and edge correlation length, e , the (1015) GaN plane was used, while for the screw values d s and s , the (0004) plane was measured.
The resulted rocking curve scans were processed applying the equation used by Kaganer et al. [26]: where Ii represents the integrated peak intensity, Ibackgr represents the background intensity, A and B are parameters that describe the dislocation density and dislocation correlation range, respectively, and x is the arbitrary direction along which the correlations are measured. The dislocation density and correlation length are associated with the A and B parameters. Their values have been extracted from a nonlinear least square fit, minimizing the difference between the observed intensity and the calculated one. They are defined as where b is the Burgers vector with the chosen values of the edge b e = 0.32 nm, while the screw b s = 0.52 nm lengths were given by the lattice parameters a and c of nominally strainfree hexagonal GaN. The dimensionless parameters f and g depend on the skew geometry of the diffraction setup. They are defined as e = 0.7 cos 2 ( ) cos 2 ( ) 4 cos 2 ( B ) ; s = 0.5 sin 2 ( ) cos 2 ( ) 4 cos 2 ( B ) ; e = 2 cos( B ) cos( ) cos( ) ; s = 2 cos( B ) cos( ) sin( ) where represents the angle between the scattering vector and the sample surface, represents the angle between the sample surface and the incident vector, and B represents the Bragg angle.  The obtained threading dislocation densities corresponding to the GaN layer grown on SiC are d e = 1.37 × 10 10 cm −2 and d s = 1.07 × 10 9 cm −2 . Consequently, the total threading dislocation density calculated as the sum of those two component densities is d t = 1.47 × 10 10 cm −2 . With the latter value, the mean distance between two dislocations was determined as being 82 nm according to . The values of the dislocation correlation lengths, i.e., e = 171 nm and s = 288 nm, extracted from the peak profiles, can  [13,14]. The resulted rocking curve scans were processed applying the equation used by Kaganer et al. [26]: where I i represents the integrated peak intensity, I backgr represents the background intensity, A and B are parameters that describe the dislocation density and dislocation correlation range, respectively, and x is the arbitrary direction along which the correlations are measured. The dislocation density and correlation length are associated with the A and B parameters. Their values have been extracted from a nonlinear least square fit, minimizing the difference between the observed intensity and the calculated one. They are defined as where b is the Burgers vector with the chosen values of the edge b e = 0.32 nm, while the screw b s = 0.52 nm lengths were given by the lattice parameters a and c of nominally strain-free hexagonal GaN. The dimensionless parameters f and g depend on the skew geometry of the diffraction setup. They are defined as where ψ represents the angle between the scattering vector and the sample surface, ϕ represents the angle between the sample surface and the incident vector, and θ B represents the Bragg angle. The obtained threading dislocation densities corresponding to the GaN layer grown on SiC are ρ e d = 1.37 × 10 10 cm −2 and ρ s d = 1.07 × 10 9 cm −2 . Consequently, the total threading dislocation density calculated as the sum of those two component densities is ρ t d = 1.47 × 10 10 cm −2 . With the latter value, the mean distance between two dislocations was determined as being 82 nm according to r d = 1/ ρ t d 1/2 [27]. The values of the dislocation correlation lengths, i.e., L e = 171 nm and L s = 288 nm, extracted from the peak profiles, can be ascribed to the low lattice mismatch between the substrate and the film. The uncertainty of the presented values is within the least significant digit. DBS is one of the most sensitive techniques used for the detection of open-volume defects in materials [28]. Whenever positrons are implanted into condensed matter, they annihilate with electrons in less than 10 −9 s. While in a defect-free material the positrons are in a Bloch state, i.e., delocalized, in the presence of defects, both negatively charged and neutral, the positrons are trapped into them before their annihilation. Thus, the annihilation gamma-rays carry information of the electrons associated with the defects. [29]. This phenomenon will then cause narrowing in the DBS Spectrum. The shape of the broadened Doppler spectrum is characterized by the S and W parameters, which will ultimately enable the assessment of the aforementioned defects. The S and W parameters reflect the changes caused by the annihilation of positrons with electrons with low and high momentum, respectively. While higher S values imply greater positron entrapment within defects, the W parameter is influenced by annihilation of positrons with core electrons and is used to detect the presence of impurities [21,30].
The DBS data were analyzed using the VEPFIT software and characterized through the sharpness parameter, S, calculated as the ratio between the counts in the annihilation peak's central area (|E γ −511 keV| < 0.78 keV) and the total peak counts (500-522 keV). Considering the available TEM information on the wafer's structure, a positron implantation profile with different layer densities was used [31]: with z = z ρ 0 ρ(ζ)/ρ 0 dζ, where z represents the depth at which the positron is located and ρ 0 stands for the density of the substrate. For performing the data analysis, the densities of GaN, AlN, and SiC were chosen as 6.15 g cm −3 , 3.26 g cm −3 , and 3.21 g cm −3 , respectively. When solving the positron transport problem, the implantation, diffusion, trapping, and annihilation of free positrons were considered. Positron surface trapping, epithermal positrons, thermal positrons diffusing back toward the surface, and positronium (Ps) emission were all integrated using the VEPFIT model: where S e , S s , and The relative estimated triple state of Ps emitted from the surface is given by F Ps (E + ), which is the ratio between the valley area counts of the energy spectrum (400-500 keV) and the total peak counts. F Ps (E + ) was fitted concomitantly on the same VEPFIT model. The effective positron diffusion length is one of the parameters derived from data analysis and is described as 1/2 , where D + represents the positron diffusion coefficient, λ b corresponds to the bulk annihilation rate, k t stands for the positron trapping rate, and n t is the defect density. In order to not influence the uncertainty of the above parameters, the diffusion length of the substrate was set to 184 nm [32]. The thicknesses of the layers were also fixed to the values determined via a TEM analysis. Figure 5 shows the depth profiles S(E + ) and F Ps (E + ) obtained for the GaN/AlN/SiC heterostructure along with the layer boundary depths calculated using the mean penetration depth z m = (36/ρ)E 1.62 + . The best fit results are shown in Table 1. Our attempts to fit the experimental data using a three-layer model resulted in physically incorrect data, so a four-layer model was proposed instead, thereby splitting the GaN film into two sublayers (i.e., GaN1 and GaN2). A detailed explanation why the three-layer model results are incorrect is given in references [13,14].  Figure 5, a strong decrease in S was observed for E + The DBS data were analyzed using the VEPFIT software and characterized through the sharpness parameter, S, calculated as the ratio between the counts in the annihilation peak's central area (|Eγ−511 keV| < 0.78 keV) and the total peak counts (500-522 keV). Considering the available TEM information on the wafer's structure, a positron implantation profile with different layer densities was used [31]: with = ( )/ d , where z represents the depth at which the positron is located and stands for the density of the substrate. For performing the data analysis, the densities of GaN, AlN, and SiC were chosen as 6.15 g cm −3 , 3.26 g cm −3 , and 3.21 g cm −3 , respectively. When solving the positron transport problem, the implantation, diffusion, trapping, and annihilation of free positrons were considered. Positron surface trapping, epithermal positrons, thermal positrons diffusing back toward the surface, and positronium (Ps) emission were all integrated using the VEPFIT model: where , which is the ratio between the valley area counts of the energy spectrum (400-500 keV) and the total peak counts. FPs (E+) was fitted concomitantly on the same VEPFIT model. The effective positron diffusion length is one of the parameters derived from data analysis and is described as where D + represents the positron diffusion coefficient, λb corresponds to the bulk annihilation rate, kt stands for the positron trapping rate, and nt is the defect density. In order to not influence the uncertainty of the above parameters, the diffusion length of the substrate was set to 184 nm [32]. The thicknesses of the layers were also fixed to the values determined via a TEM analysis. Figure 5 shows the depth profiles S(E+) and FPs(E+) obtained for the GaN/AlN/SiC heterostructure along with the layer boundary depths calculated using the mean penetration depth zm = (36/ρ) . . The best fit results are shown in Table 1. Our attempts to fit the experimental data using a three-layer model resulted in physically incorrect data, so a four-layer model was proposed instead, thereby splitting the GaN film into two sublayers (i.e., GaN1 and GaN2). A detailed explanation why the three-layer model results are incorrect is given in references [13,14].  Figure 5, a strong decrease in S was observed for E+ ≲ 1 keV due to the reduced ability of thermalized positrons to diffuse back toward the surface. This phenomenon is induced by the local electric field resulting from the band bending near the surface. Consequently, less Ps emission takes place at the surface and thus the positron diffusion length is reduced. At E+ ≳ 1 keV, S increases slowly with E+ in the GaN film area, then 1 keV due to the reduced ability of thermalized positrons to diffuse back toward the surface. This phenomenon is induced by the local electric field resulting from the band bending near the surface. Consequently, less Ps emission takes place at the surface and thus the positron diffusion length is reduced. At E + The DBS data were analyzed using the VEPFIT software and characterized throu the sharpness parameter, S, calculated as the ratio between the counts in the annihilati peak's central area (|Eγ−511 keV| < 0.78 keV) and the total peak counts (500-522 keV Considering the available TEM information on the wafer's structure, a positron implan tion profile with different layer densities was used [31]: with = ( )/ d , where z represents the depth at which the positron is located an stands for the density of the substrate. For performing the data analysis, the densiti of GaN, AlN, and SiC were chosen as 6.15 g cm −3 , 3.26 g cm −3 , and 3.21 g cm −3 , respective When solving the positron transport problem, the implantation, diffusion, trapping, a annihilation of free positrons were considered. Positron surface trapping, epithermal po itrons, thermal positrons diffusing back toward the surface, and positronium (Ps) em sion were all integrated using the VEPFIT model: , which the ratio between the valley area counts of the energy spectrum (400-500 keV) and t total peak counts. FPs (E+) was fitted concomitantly on the same VEPFIT model. The effective positron diffusion length is one of the parameters derived from da analysis and is described as where D + represents the positr diffusion coefficient, λb corresponds to the bulk annihilation rate, kt stands for the positr trapping rate, and nt is the defect density. In order to not influence the uncertainty of t above parameters, the diffusion length of the substrate was set to 184 nm [32]. The thic nesses of the layers were also fixed to the values determined via a TEM analysis. Figure 5 shows the depth profiles S(E+) and FPs(E+) obtained for the GaN/AlN/S heterostructure along with the layer boundary depths calculated using the mean penetr tion depth zm = (36/ρ) . . The best fit results are shown in Table 1. Our attempts to fit t experimental data using a three-layer model resulted in physically incorrect data, so four-layer model was proposed instead, thereby splitting the GaN film into two sublaye (i.e., GaN1 and GaN2). A detailed explanation why the three-layer model results are i correct is given in references [13,14]. Table 1. Best fit parameters obtained via VEPFIT software using the S(E+) and FPs(E+) depth profiles. Values without any error margins represent fixed parameters.

Sample
GaN/SiC S  Figure 5, a strong decrease in S was observed for E+ ≲ 1 keV due to t reduced ability of thermalized positrons to diffuse back toward the surface. This phenom enon is induced by the local electric field resulting from the band bending near the surfa Consequently, less Ps emission takes place at the surface and thus the positron diffusi length is reduced. At E+ ≳ 1 keV, S increases slowly with E+ in the GaN film area, th 1 keV, S increases slowly with E + in the GaN film area, then seemingly reaches a plateau in the AlN layer (for E + The DBS data were analyzed using the VEPFIT softw the sharpness parameter, S, calculated as the ratio between peak's central area (|Eγ−511 keV| < 0.78 keV) and the to Considering the available TEM information on the wafer's tion profile with different layer densities was used [31]: , = / ( , with = ( )/ d , where z represents the depth at w stands for the density of the substrate. For performing of GaN, AlN, and SiC were chosen as 6.15 g cm −3 , 3.26 g cm When solving the positron transport problem, the implan annihilation of free positrons were considered. Positron su itrons, thermal positrons diffusing back toward the surfa sion were all integrated using the VEPFIT model: where Se, Ss, and Si represent the characteristic parameter annihilation of epithermal positrons (Se) and thermalized annihilation of thermalized positrons within the bulk of t (Si). Fe(E+) stands for the fraction of epithermal positrons a corresponds to the fraction of surface-annihilated thermali sents the fraction of epithermal positrons annihilated with relative estimated triple state of Ps emitted from the surfa the ratio between the valley area counts of the energy sp total peak counts. FPs (E+) was fitted concomitantly on the The effective positron diffusion length is one of the analysis and is described as corresponds to the bulk annihilatio trapping rate, and nt is the defect density. In order to not i above parameters, the diffusion length of the substrate wa nesses of the layers were also fixed to the values determin Figure 5 shows the depth profiles S(E+) and FPs(E+) heterostructure along with the layer boundary depths calc tion depth zm = (36/ρ) . . The best fit results are shown in experimental data using a three-layer model resulted in four-layer model was proposed instead, thereby splitting t (i.e., GaN1 and GaN2). A detailed explanation why the th correct is given in references [13,14]. According to Figure 5, a strong decrease in S was obs reduced ability of thermalized positrons to diffuse back tow enon is induced by the local electric field resulting from the Consequently, less Ps emission takes place at the surface a length is reduced. At E+ ≳ 1 keV, S increases slowly with 17 keV), before finally achieving full saturation at E + Nanomaterials 2021, 11, x FOR PEER REVIEW 7 of 11 The DBS data were analyzed using the VEPFIT software and characterized through the sharpness parameter, S, calculated as the ratio between the counts in the annihilation peak's central area (|Eγ−511 keV| < 0.78 keV) and the total peak counts (500-522 keV). Considering the available TEM information on the wafer's structure, a positron implantation profile with different layer densities was used [31]: with = ( )/ d , where z represents the depth at which the positron is located and stands for the density of the substrate. For performing the data analysis, the densities of GaN, AlN, and SiC were chosen as 6.15 g cm −3 , 3.26 g cm −3 , and 3.21 g cm −3 , respectively. When solving the positron transport problem, the implantation, diffusion, trapping, and annihilation of free positrons were considered. Positron surface trapping, epithermal positrons, thermal positrons diffusing back toward the surface, and positronium (Ps) emission were all integrated using the VEPFIT model: where represents the fraction of epithermal positrons annihilated within the bulk of the i-th layer. The relative estimated triple state of Ps emitted from the surface is given by FPs (E+), which is the ratio between the valley area counts of the energy spectrum (400-500 keV) and the total peak counts. FPs (E+) was fitted concomitantly on the same VEPFIT model.
The effective positron diffusion length is one of the parameters derived from data analysis and is described as represents the positron diffusion coefficient, λb corresponds to the bulk annihilation rate, kt stands for the positron trapping rate, and nt is the defect density. In order to not influence the uncertainty of the above parameters, the diffusion length of the substrate was set to 184 nm [32]. The thicknesses of the layers were also fixed to the values determined via a TEM analysis. Figure 5 shows the depth profiles S(E+) and FPs(E+) obtained for the GaN/AlN/SiC heterostructure along with the layer boundary depths calculated using the mean penetration depth zm = (36/ρ) . . The best fit results are shown in Table 1. Our attempts to fit the experimental data using a three-layer model resulted in physically incorrect data, so a four-layer model was proposed instead, thereby splitting the GaN film into two sublayers (i.e., GaN1 and GaN2). A detailed explanation why the three-layer model results are incorrect is given in references [13,14].  Figure 5, a strong decrease in S was observed for E+ ≲ 1 keV due to the reduced ability of thermalized positrons to diffuse back toward the surface. This phenomenon is induced by the local electric field resulting from the band bending near the surface. Consequently, less Ps emission takes place at the surface and thus the positron diffusion length is reduced. At E+ ≳ 1 keV, S increases slowly with E+ in the GaN film area, then 21 keV within the SiC substrate layer. seemingly reaches a plateau in the AlN layer (for E+ ≳ 17 keV), before finally achieving full saturation at E+ ≳ 21 keV within the SiC substrate layer. Since Ps is not formed within the bulk GaN layers, the inequality SGaN1 = 0.4576 ± 0.0004 < SGaN2 = 0.4615 ± 0.0004 indicates either a smaller (in size) or lesser number of defects as well as a longer Leff within the GaN1 sublayer [33]. Given the resulting disagreement in eff GaN1 = 13 ± 0.4 nm being smaller than eff GaN2 = 75 ± 20 nm, the underlying reason could be related to the presence of a local electric field directed toward the surface. The results which summarize the characteristics of the GaN layers listed in this study along with previously reported results [13,14] are shown in Table 2. Additionally, Figure  4c,d show the omega scans for (0004) and (1015) planes for GaN films grown on (111) 3C-SiC, (0001) Al2O3, and (111) Si substrates [13,14]. Comparing the results, one can observe that GaN exhibits an improved crystal quality in the GaN/SiC sample, with the lowest structural defect density of d t = 1.47 × 10 10 cm −2 and the highest value of eff GaN2 = 78 ± 20 nm. All of the investigated samples (i.e., GaN/Al2O3 [13], GaN300/Si, GaN700/Si [14], and GaN/SiC) demonstrated shorter effective positron diffusion lengths than that of defectfree GaN (i.e., eff DF =135 nm), with the closest value being obtained for GaN/SiC. It is important to note that the GaN/Al2O3 mean distance between two dislocations was calculated within this paper and the positron data was revised by considering the relative estimated triple state of Ps emitted from the surface. Since Ps is not formed within the bulk GaN layers, the inequality S GaN1 = 0.4576 ± 0.0004 < S GaN2 = 0.4615 ± 0.0004 indicates either a smaller (in size) or lesser number of defects as well as a longer L eff within the GaN1 sublayer [33]. Given the resulting disagreement in L GaN1 eff = 13 ± 0.4 nm being smaller than L GaN2 eff = 75 ± 20 nm, the underlying reason could be related to the presence of a local electric field directed toward the surface.
The results which summarize the characteristics of the GaN layers listed in this study along with previously reported results [13,14] are shown in Table 2. Additionally, Figure 4c,d show the omega scans for (0004) and (1015) planes for GaN films grown on (111) 3C-SiC, (0001) Al 2 O 3 , and (111) Si substrates [13,14]. Comparing the results, one can observe that GaN exhibits an improved crystal quality in the GaN/SiC sample, with the lowest structural defect density of ρ t d = 1.47 × 10 10 cm −2 and the highest value of L GaN2 eff = 78 ± 20 nm. All of the investigated samples (i.e., GaN/Al 2 O 3 [13], GaN300/Si, GaN700/Si [14], and GaN/SiC) demonstrated shorter effective positron diffusion lengths than that of defect-free GaN (i.e., L DF eff = 135 nm), with the closest value being obtained for GaN/SiC. It is important to note that the GaN/Al 2 O 3 mean distance between two dislocations was calculated within this paper and the positron data was revised by considering the relative estimated triple state of Ps emitted from the surface.
The smallest mean distance between two dislocations, r d = 15 nm, which is a consequence of the high value of the GaN300/Si wafer dislocation density (ρ t d = 4.37 × 10 11 cm −2 ), implies a lower quality of the GaN layer as compared to the GaN/Al 2 O 3 , GaN700/Si, and GaN/SiC wafers. The conclusion of the lower crystal quality is also sustained by the linear dislocations and Al diffusion [14] near the AlN/GaN interface, as well as by the highest FWHM in omega scans for this GaN300/Si layer, as shown in Figure 4c,d. Even though GaN700/Si displays more point defects and linear dislocations at the interface, the decreased lengths of elemental diffusion relative to the GaN layer thickness induces a better quality of the top film, lowering the amount of threading dislocation densities [14]. Table 2. Total defect densities, defect correlation lengths, and effective positron diffusion lengths in GaN, where d is the layer thickness, ρ t d is the total threading dislocation density, r d is the mean distance between two dislocations, L is the defect correlation length (edge and screw), and L eff is the effective positron diffusion length.

Conclusions
Three methods were used for evaluating the relationship between the substrate, buffer, and film heterostructure and the consistency of its interfaces, layer thicknesses, dislocation densities, and positron-characteristic diffusion lengths in the GaN/AlN/SiC heterostructure. Edge and screw dislocations were assessed within the epitaxially grown layers defined by the relationship P63mc (0002) GaN || P63mc (0002) AlN || F43m (111) 3C-SiC. The total dislocation density has been assessed as being 1.47 × 10 10 cm −2 . Compared with previously investigated heterostructures (on Si and Al 2 O 3 substrates), the obtained dislocation correlation lengths (L e = 171 nm and L s = 288 nm) and the mean distance between two dislocations (r d = 82 nm) are higher. This reveals an improved crystal quality of the GaN with SiC as a growth template. In addition, the DBS measurements upheld the aforementioned results, with higher effective positron diffusion length L GaN2 eff = 75 ± 20 nm for the GaN layer.