The luminescent and recombination properties of V-pit defects in p-GaN(Mg) grown by metalorganic chemical vapor deposition (MOCVD) were studied by scanning electron microscopy (SEM) in the secondary electron, cathodoluminescence (CL), and electron beam induced current (EBIC) modes, combined with CL spectra measurements and EBIC collection efficiency measurements. Similar studies were performed on low-dislocation-density freestanding n-GaN crystals. For MOCVD p-GaN films, the SEM investigations were supplemented by capacitance–voltage, current–voltage, deep level transient spectroscopy analysis with Ni Schottky diode, and Ohmic contacts. These experiments show that V-pits in p-GaN increase the leakage current of Schottky diodes, as in n-GaN films and crystals. EBIC imaging and EBIC collection efficiency results suggest that in the region of V-pits, a parasitic p–n junction is formed. We also observe that, in V-pits, the CL spectra the contribution of the 3.2 eV defect band is strongly enhanced compared to the 3 eV blue CL band that dominates the spectra.

GaN and its solid solutions are currently widely used for fabrication of light emitting devices (LEDs) operating in the visible and UV spectral range and in high-power devices, such as rectifiers and field effect transistors.1,2 However, in spite of rapid progress in the development of GaN growth technology with low dislocation densities, many optical and electrical properties of extended defects in this material are not fully understood.3 One of the most common defects in GaN-based structures is a so-called V-pit, a hexagonal pit formed by threading dislocations when they cross the surface of the growing layer.4 The facets of these V-pits are formed by six {10−11} planes inclined in respect to the (0001) growth surface at 61° angle.5–8 Often, but not always, these V-pits can have a flat bottom.6–9 These defects are believed to be partly responsible for the GaN-based LEDs demonstrating good quantum efficiency, despite a very high density of threading dislocations exceeding ∼108 cm−2 when grown on alien substrates (sapphire, SiC, Si).4−10 The explanation commonly offered is that the quantum wells (QWs) in LEDs are inclined on the facets of the V-pit and are thinner than on the flat surface and thus have a higher bandgap. As a result, recombination mainly occurs in the narrow-bandgap QWs portions on the flat surface and the impact of nonradiative recombination via threading dislocations is strongly diminished.4−10 

For that reason, it is often considered advantageous to have as many threading dislocations as possible to be terminated with V-pits and the dimensions of the pits to be large.4−11 The situation is not as straightforward as that: In concentration (and hence the bandgap of the GaN/InGaN QW) on the facets of the V-pits can be measurably different from the flat portion of the structure, electrically active impurities can be segregated on the facets of the V-pits, the hole injection efficiency into the QWs can be different in the vicinity of V-pits, and the density of nonradiative recombination defects can vary near the V-pits.10,12−15 Detailed understanding of the properties of V-pits is necessary to optimize the performance of GaN-based LEDs.

The properties of V-pits have been studied in detail for n-GaN, where it has been shown that the pit facets are easily decorated with oxygen and thus show a higher shallow donor density than in the flat portions of the films. The result is a strong increase of the leakage current of the Schottky diodes.13 It was also observed that enhanced shallow donor density on the facets of the V-pits leads to a red shift of the band edge photoluminescence (PL) peak, more pronounced for the facets than for the flat bottom.6,7 For GaN/InGaN QW structures, it has been noted that the In concentration on the facets of the V-pits can be higher than in the flat (0001) oriented part of the film, which can overpower in some cases the increase of the bandgap due to the thinning of the QW on the facets and result in either blue shift or red shift of the QW PL peak near the V-pit and to possible negative effects of the V-pits on QW luminescence efficiency.10−12 Among the growth factors, it has been noted that the addition of In surfactant can influence the density and dimensions of the V-pits in QW structures,16 while increasing the temperature of the barrier growth can increase the dimensions of the V-pits.17 However, virtually nothing is known about the optical and electrical properties of the V-pits in p-type GaN. In this paper, we present the results of V-pit studies performed by cathodoluminescence (CL) and electron beam induced current (EBIC) measurements. For comparison, the results of similar measurements are presented for V-pits in freestanding n-GaN crystals.

The samples used were (a) a p-GaN film grown by metalorganic chemical vapor deposition (MOCVD) and (b) an unintentionally doped (UID) n-GaN freestanding crystal grown by hydride vapor phase epitaxy (HVPE). The p-GaN film was grown on basal plane sapphire and consisted of 1.6 μm-thick unintentionally doped (UID) n-GaN buffer, and 2 μm-thick layer doped with Mg to the Mg concentration of (3–4) × 1019 cm−3 determined by secondary ion mass spectrometry (SIMS) calibration. The film was grown in the same way as the p-GaN emitter layer in our standard blue LED process fabrication. The electrical properties of the UID n-GaN buffer were measured after a calibration run with no p-GaN deposition. These measurements showed that the net shallow donor concentration at the surface, as determined by capacitance–voltage (C–V) profiling in UID n-GaN buffer, was below 1016 cm−3.

The HVPE sample was acquired from Kyma Technologies Inc. (USA) and had thickness of 650 μm, the net donor concentration of 5 × 1016 cm−3, and the threading dislocation density below 107 cm−2. Electrical properties and deep electron and hole traps spectra were studied previously and are described in detail elsewhere.18,19

The CL studies were carried out in a scanning electron microscope JSM-6490 (Jeol, Japan) with the MonoCL-3 system. A Hamamatsu photomultiplier was used as a detector. Measurements were carried out with beam energy of 10 keV and beam current of 1–10 nA in the temperature range from 80 to 300 K. EBIC measurements were carried out at room temperature in the scanning electron microscope JSM-840 (Jeol, Japan) with beam current about 0.1 nA and beam energy in the range from 15 to 39 keV. The EBIC collection efficiency was taken as the EBIC current Ic normalized by the product of the probing electrons energy Eb and beam current Ib, Ic/(Eb × Ib). Experimental setups are described in detail elsewhere.3,20

For the MOCVD p-GaN sample, we also performed electrical and deep trap spectra studies including (a) current–voltage (I–V) measurements in the dark and under illumination with a high-power (optical output power density 250 mW/cm2) GaN-based LED with peak wavelength of 365 nm and (b) capacitance vs frequency (C–f) measurements (frequency range 20 Hz–1 MHz) in the dark and under illumination with high-power LEDs with peak wavelengths from 950 to 365 nm and optical power density 250 mW/cm2, capacitance vs voltage (C–V) measurements in the frequency range 20 Hz–1 MHz in the dark and under illumination, deep level transient spectroscopy with electrical (DLTS) or optical (ODLTS) injection. For these measurements; for the EBIC characterization, the contacts were prepared by e-beam evaporation through a shadow mask of Ni/Au (50/200 nm) and rapid thermal annealing at 600°C for 30 s in nitrogen, followed by room temperature deposition of 30-nm-thick, 1 mm in diameter Ni dots. The Ni/Au provides the Ohmic contact, while Ni forms a Schottky diode to p-GaN. Experimental setups are described elsewhere.21−23 

I–V characteristics of the MOCVD p-GaN measured at 295, 130, and 400 K are shown in Fig. 1. The figure also shows the I–V characteristics measured at room temperature and 130 K with 365 nm LED excitation (photon energy 3.4 eV) (at 400 K the photocurrent was strongly masked by the increased dark current). The plus sign in the voltage axis in this figure corresponds to positive voltage applied to the Ni Schottky contact. There are three key points from the data: (a) the current for two voltage directions is similar, (b) the temperature dependence of the dark current is quite weak (the activation energy is close to 0.1 eV, the actual dependence not shown for the lack of space), and (c) the open-circuit voltage is positive, as would be the case if the Ni Schottky diodes were shorted and the signals were determined by the two p-GaN/n-GaN junctions turned face-to-face, with the contribution of the junction with the thinner Ni electrode dominant because of the lower area and better contact transparency (in the case of I–V measurements under illumination). As shown in Sec. III B, the V-pit density in the p-GaN film is quite high and a credible explanation of the lack of rectification at the Ni/p-GaN Schottky diode seems to be the strong leakage via the V-pits, as in the case of n-type GaN.13 

FIG. 1.

Dark I–V characteristics measured at 130 K (blue line), 297 K (black line), 400 K (red line); also shown are I–V characteristics measured with 365 nm LED excitation at 130 K (magenta line), 297 K (orange line), and 400 K (red line).

FIG. 1.

Dark I–V characteristics measured at 130 K (blue line), 297 K (black line), 400 K (red line); also shown are I–V characteristics measured with 365 nm LED excitation at 130 K (magenta line), 297 K (orange line), and 400 K (red line).

Close modal

Capacitance–voltage results also confirm the lack of rectification by the Ni Schottky diode and indicate the charged center concentration to be 9.3 × 1017 cm−3, with the built-in voltage close to 3 eV, as is often observed in GaN p–n junctions.23 Illumination with high-power 365 nm wavelength LED increased the photo-concentration to 1018 cm−3 (Fig. 2). The measured concentration is too high to correspond to the space charge boundary movement with applied voltage in the n-GaN buffer of the film, so that the measured concentration is most likely due to the p-GaN side of the junction. The reason is that the n-side is almost depleted down to the ∼(0.5–1) μm adjacent to the sapphire substrate where the donor concentration is very high23 so that the space charge region boundary movement with changing voltage mostly occurs in the p-type portion of the p–n junction.

FIG. 2.

1/C2 vs voltage plots for the p-GaN MOCVD sample measured at 10 kHz in the dark (black line) and under 365 nm LED illumination (violet line).

FIG. 2.

1/C2 vs voltage plots for the p-GaN MOCVD sample measured at 10 kHz in the dark (black line) and under 365 nm LED illumination (violet line).

Close modal

The DLTS spectrum of the sample measured at 10 kHz to alleviate problems due to the high series resistance21 is shown in Fig. 3. The dominant feature of the spectrum is the peak near 300 K due to a prominent majority carriers trap (the capacitance decreases with time after the application of the forward bias pulse).24 The activation energy is 0.85 eV and capture cross section is 8.9 × 10−14 cm2 if the signal comes from the p-GaN region, and the trapped charges are holes, as suggested by I–V and C–V measurements. The parameters of these deep acceptors are close to the characteristics of gallium-vacancy-related deep acceptor complexes with shallow donors (VGa-D)2− or to C acceptors on N site, with CN expected to be major deep hole traps in n-GaN.23 The VGa-related defects are not expected to be present in high concentrations in p-GaN because of their high formation energy when the Fermi level is close to the valence band edge,25,26 but the CN acceptor traps have been detected in p-GaN by DLTS.27 Another possibility28 is a deep donor with a transition level near Ev + 0.85 eV, which has been predicted for a complex of nitrogen vacancy donors VN3+ with Mg acceptors, (VN-Mg)2+ and can be easily formed in p-type material. In principle, these two scenarios can be distinguished by measuring CL spectra: theory predicts a broad luminescence line in the yellow–green region of 2.3–2.5 eV for the CN,25,26 whereas for (VN-Mg)2+ deep donor, the respective luminescence band should be lying in the red spectral region and centered near 1.8 eV.28 However, detailed CL spectra measurements allowing to decide which of the defects is observed have not yet been done for the red–green spectral range.

FIG. 3.

DLTS (blue line) and ODLTS (red line) spectra measured with bias −1V, time windows 1.75/17.5 s, and excitation pulse length 5 s at 10 kHz; in DLTS, the forward bias pulse was +1 V, and in ODLTS the excitation pulse was provided by illumination with 365 nm LED.

FIG. 3.

DLTS (blue line) and ODLTS (red line) spectra measured with bias −1V, time windows 1.75/17.5 s, and excitation pulse length 5 s at 10 kHz; in DLTS, the forward bias pulse was +1 V, and in ODLTS the excitation pulse was provided by illumination with 365 nm LED.

Close modal

The same hole traps could also be clearly seen in ODLTS spectra measured with an excitation pulse provided by illumination with high-power LED with photon energy 3.4 eV (Fig. 3). In addition, an electron trap with a peak near 210 K and the level near Ec − 0.45 eV and electron capture cross section of 1.6 × 10−16 cm2 could be detected. Such traps have been previously reported in some n-GaN films.23 

Secondary electron (SE) imaging of the MOCVD p-GaN sample and the freestanding HVPE n-GaN sample are shown in Figs. 4 and 5. For the p-GaN, there is a high density (106–107 cm−2) of hexagonal V-pits with flat bottom and in-plane dimensions of about 1 μm, which allows us to calculate their depth as being below 1 μm (Fig. 4). In freestanding HVPE n-GaN, the V-pit density does not exceed 102 cm−2; however, their size is larger than 100 μm (Fig. 5).

FIG. 4.

SE image of V-pits in p-GaN.

FIG. 4.

SE image of V-pits in p-GaN.

Close modal
FIG. 5.

SE image of V-pit in HVPE n-GaN.

FIG. 5.

SE image of V-pit in HVPE n-GaN.

Close modal

The room temperature CL spectra measured in HVPE n-GaN are shown in Fig. 6. On the facets of the V-pits, the CL intensity increases and the CL peak position shifts to lower energies, which correlates with the results published in Ref. 6. The spectrum measured far from a defect can be described with one Gaussian peak centered at 3.4 eV with width of 50 meV. The spectrum in the V-pit center can be fitted with two Gaussian peaks exp[−(E − Em)2/2σ2] with energies Em of 3.375 and 3.32 eV and σ equal to 45 and 80 meV, respectively. The spectrum in the V-pit side facet can be fitted also with two Gaussian peaks with Em of 3.25 and 3.05 eV and σ equal to 95 and 100 meV, respectively. The Schottky diodes made on V-pits of the HVPE sample demonstrate a practically Ohmic behavior, i.e., they are shorted by V-pits. This can indicate an increased donor concentration in V-pits and correlates with oxygen donor segregation on the surface framing the V-pit.6 The origin of the emission bands observed is unknown; however, different bands are dominant for the V-pit-free “bulk” area and for the facets and bottom regions of V-pits. Thus, it can be assumed that different defects are accumulated in the V-pit facet and bottom regions. The emission band at 3.05 eV is close to that produced by dislocations gliding at room temperature.15,20

FIG. 6.

Room temperature CL spectra measured in HVPE n-GaN in the center of V-pit (red line) and on its side (olive line). The spectrum measured far from the V-pit is shown with the blue line. The CL intensity on the side is two and five times higher than in the center and far from defect, respectively. The results of spectra deconvolution are shown with dashed lines with the same color coding as the measured spectra.

FIG. 6.

Room temperature CL spectra measured in HVPE n-GaN in the center of V-pit (red line) and on its side (olive line). The spectrum measured far from the V-pit is shown with the blue line. The CL intensity on the side is two and five times higher than in the center and far from defect, respectively. The results of spectra deconvolution are shown with dashed lines with the same color coding as the measured spectra.

Close modal

Due to a smaller size of V-pits in the MOCVD p-GaN sample, it was difficult to separate the spectra from their edges and centers; therefore, the spectra measured present some average ones. Typical CL spectra obtained in p-GaN on the V-defect and far from it, obtained at 80 and 300 K, are shown in Fig. 7. In the CL spectra of bulk p-GaN, there is no clearly defined luminescence line associated with the near-band-edge transition. This is in good agreement with the results of Refs. 29 and 30 at high magnesium concentrations. A small peak close to the near-band-edge transition is present in some V-pit spectra; however, it can come from the n-GaN buffer because in pits, the e-beam can penetrate deeper into the sample. The spectra measured on V-pits are shifted to higher energies when compared to those on the defect-free region. This trend is opposite to the case of n-GaN. However, fitting the spectra with Gaussian peaks shows that in the range from 2.5 to 3.6 eV, they can be described by superposition of three peaks with Em of 3.28 and 3.2 eV for both V-pits and defect-free region, and the most intense peak at 3.03 eV for V-pits and 2.98–3 eV for the defect-free region. Em for all peaks is practically independent of temperature in the range 80–300 K. The relative intensities of Gaussian peaks normalized by 3–3.03 eV peak intensity for V-pit and two defect-free regions are shown in Fig. 8. While the positions of Gaussian peaks in defect-free regions are the same across the sample, their intensities varied in different regions. Therefore, the data for defect-free regions are shown in Fig. 8 for two typical cases.

FIG. 7.

CL spectra measured on the p-GaN MOCVD sample at 80 (a) and 300 K (b) on the V-pit (magenta) and far from it (blue).

FIG. 7.

CL spectra measured on the p-GaN MOCVD sample at 80 (a) and 300 K (b) on the V-pit (magenta) and far from it (blue).

Close modal
FIG. 8.

The intensity of Gaussian peaks normalized by the intensity of 3–3.03 V peak at 80 (a) and 300 K (b). The results are shown for the V-pit region (magenta) and for two representative places in the defect-free region (blue).

FIG. 8.

The intensity of Gaussian peaks normalized by the intensity of 3–3.03 V peak at 80 (a) and 300 K (b). The results are shown for the V-pit region (magenta) and for two representative places in the defect-free region (blue).

Close modal

Thus, it can be assumed that the emission in both defect-free and V-pit regions is determined by the same defects. However, in V-pits, the relative intensities of high energy peaks are higher, which explains the overall spectral shift to higher energies. The emission peaks with similar Em have been already observed in GaN with high Mg concentration.29 Nevertheless, some differences remain. As reported in Ref. 29, the intensity of the UV band at 3.2 eV practically disappears at room temperature, and at 300 K, the blue band at 3 eV dominates in the spectrum. As seen in Fig. 8, the opposite behavior is observed in our work, i.e., relative intensity of UV bands increases at room temperature. Contrary to observations in Ref. 29, the position of the blue band is practically independent of the temperature. The origin of the blue and UV luminescence bands in p-GaN (Mg) is still under debate. Theoretical analysis25,26 suggests that the blue line near 3 eV is due to a strongly lattice-coupled quasi-deep isolated Mg acceptor center, whereas the UV line near 3.2 eV could be related to a transition involving a donor state of the Mg–H complex with the transition level near Ev + 0.13 eV. One can easily envisage the change in relative concentrations of such centers in the V-pit-free and V-pit regions of the p-GaN film.

In the EBIC mode, V-pits in p-GaN produce a bright contrast, i.e., the collected current increases at them (Fig. 9). The current direction was the same as would be observed for the p-GaN/n-GaN junction with the Ni Schottky diode effectively shorted, most likely by the V-pits, as already inferred in Sec. III A from I–V and C–V measurements. The induced current is collected by the p–n junction between the n-type buffer layer and the p-type epilayer. As follows from Fig. 9, V-pits essentially enhance the excess currier collection. Pits usually increase the collected current due to the suppression of backscattering electron emission, which in turn increases the absorbed energy. However, in this case, the current increase cannot exceed 15%–20%. In the samples studied, the contrast exceeds 50% at 35 keV and the increase is even larger at lower beam energies; therefore, this purely “geometric” mechanism can be excluded.

FIG. 9.

EBIC image of V-pits in p-GaN. Eb = 35 keV; the two bright spots are due to V-pits.

FIG. 9.

EBIC image of V-pits in p-GaN. Eb = 35 keV; the two bright spots are due to V-pits.

Close modal

The dependence of collected current Ic normalized by the product of beam energy (Eb) and current (Ib) on electron range R is shown in Fig. 10. The electron range is calculated as R = 0.47 × RB, where RB = 13.2 × Eb1.75 nm is the electron Bethe range.31,32 Far from V-pits, the collected current starts to increase for R > 1.5 μm. That means at smaller R, the electrons cannot reach the depletion region of p–n junction, which is in agreement with the p-layer thickness. With the increase of R, the number of excess carriers generated near the depletion region increases; therefore, the collected current also increases. To check if the collected current increase is determined by the layer thickness decrease in the pits, the curve in Fig. 10 is shifted by 1000 nm (the maximum pit depth). The collected current values in the pit are still larger than far from the pit, and the normalized EBIC current dependence on depth is different from the one observed in defect-free regions. Thus, the EBIC results allow us to assume that in V-pits, conductivity changes to n-type and the p–n junction formed around V-pits increases the collection efficiency. The outer radius of this n-type cylinder can be quite small; therefore, when the excess carriers are generated far from the buffer layer, the probability to reach the buffer layer decreases, which determined the current decrease at small R. The p–n junction formation around V-pit defects well correlates with the increase of oxygen concentration near them because oxygen is known to be a shallow donor in GaN.

FIG. 10.

The dependences of normalized collected current on the electron range for V-pit (red solid squares) and far from it (blue open squares). The dependence for V-pit shifted by 1000 nm (approximate depth of the V-pit) along the R axis is shown with olive stars.

FIG. 10.

The dependences of normalized collected current on the electron range for V-pit (red solid squares) and far from it (blue open squares). The dependence for V-pit shifted by 1000 nm (approximate depth of the V-pit) along the R axis is shown with olive stars.

Close modal

Our results strongly suggest that V-pits in p-GaN can short Schottky diodes, similarly to what is observed in n-GaN. The reason is likely related to oxygen donor segregation on the side facets and the bottom of the V-pit funnel. In p-GaN, judging by the EBIC collection efficiency results, this is accompanied by the formation of a parasitic p–n junction around the V-pit. In GaN/InGaN multi quantum well (MQW) LEDs, the thickness of the p-GaN emitter layer is typically below 0.2 μm1 and thus lower than the possible depth of the V-pit. Therefore, the high density of V-pits in such layers can lead to excessive leakage of the p–n junction and considerable loss of injection efficiency.

We also observe that, in the vicinity of V-pits, the relative intensity of Mg-related centers giving rise to the blue band near 3 eV and the UV band near 3.2 eV is seriously redistributed in favor of the latter, suggesting, based on theoretical predictions, that the V-pit region is enriched by the Mg–H complexes that have survived the Mg activation annealing.

The work at NUST MISiS was financially supported by the Ministry of Science and Higher Education of the Russian Federation in the framework of Increase Competitiveness Program of NUST “MISiS” (No. К2-2020-040). The work at IMT RAS was supported in part by the State Task No 075-00355-21-00. The work at Korea University was supported by the Technology Innovation Program (No. 20004479, Micro LED chip and module fabrication technology with integrated field effect transistor for IoT application) funded by the Ministry of Trade, Industry & Energy (MOTIE, Korea). The work at UF was sponsored by NSF DMR (No. 1856662 (James Edgar)].

The data that support the findings of this study are available within the article.

1.
S.
Nakamura
,
Angew. Chem. Int. Ed.
54
,
7770
7788
(
2015
).
2.
Power GaN Devices: Materials, Applications and Reliability
, edited by
M.
Meneghini
,
G.
Meneghesso
, and
E.
Zanoni
(
Springer International Publishing
,
2017
), p.
379.
3.
E. B.
Yakimov
,
A. Y.
Polyakov
,
I.-H.
Lee
, and
S. J.
Pearton
, J. Appl. Phys.
123
,
161543
(
2018
).
4.
A.
Hangleiter
,
F.
Hitzel
,
C.
Netzel
,
D.
Fuhrmann
,
U.
Rossow
,
G.
Ade
, and
P.
Hinze
,
Phys. Rev. Lett.
95
,
127402
(
2005
).
5.
J.
Kim
,
Y.-H.
Cho
,
D.-S.
Ko
,
X.-S.
Li
,
J.-Y.
Won
,
E.
Lee
,
S.-H.
Park
,
J.-Y.
Kim
, and
S.
Kim
,
Opt. Express
22
,
A857
A866
(
2014
).
6.
T.
Sekiguchi
,
W.
Lee
,
X.
Luo
,
T.
Kimura
, and
Y.
Cho
,
Phys. Status Solidi C
14
,
1700054
(
2017
).
7.
M.
Zhang
,
D.
Cai
,
Y.
Zhang
,
X.
Su
,
T.
Zhou
,
M.
Cui
,
C.
Li
,
J.
Wang
, and
K.
Xu
,
Mater. Lett.
198
,
12
(
2017
).
8.
J.-D.
Gao
,
J.-L.
Zhang
,
X.
Zhu
,
X.-M.
Wu
,
C.-L.
Mo
,
S.
Pan
,
J.-L.
Liu
, and
F.-Y.
Jiang
,
J. Appl. Cryst.
52
,
637
(
2019
).
9.
M.
Kim
,
S.
Choi
,
J.-H.
Lee
,
C.
Park
,
T.-H.
Chung
,
J. H.
Baek
, and
Y.-H.
Cho
,
Sci. Rep.
7
,
42221
(
2017
).
10.
F. C.-P.
Massabuau
,
P.
Chen
,
M. K.
Horton
,
S. L.
Rhode
,
C. X.
Ren
,
T. J.
O'Hanlon
,
A.
Kovács
,
M. J.
Kappers
,
C. J.
Humphreys
,
R. E.
Dunin-Borkowski
, and
R. A.
Oliver
,
J. Appl. Phys.
121
,
013104
(
2017
).
11.
U.
Jahn
,
O.
Brandt
,
E.
Luna
,
X.
Sun
,
H.
Wang
,
D. S.
Jiang
,
L. F.
Bian
, and
H.
Yang
,
Phys. Rev. B
81
(
12
),
125314
(
2010
).
12.
X.
Wu
,
J.
Liu
, and
F.
Jiang
,
J. Appl. Phys.
118
,
164504
(
2015
).
13.
S.
Besendorfer
,
E.
Meissner
,
A.
Tajalli
,
M.
Meneghini
,
J. A.
Freitas
, Jr.
,
J.
Derluyn
,
F.
Medjdoub
,
G.
Meneghesso
,
J.
Friedrich
, and
T.
Erlbacher
,
J. Appl. Phys.
127
,
015701
(
2020
).
14.
S.
Kurai
,
S.
Higaki
,
N.
Imura
,
K.
Okawa
,
R.
Makio
,
N.
Okada
,
K.
Tadatomo
, and
Y.
Yamada
,
Phys. Status Solidi B
255
,
1700358
(
2018
).
15.
O.
Medvedev
,
O.
Vyvenko
,
E.
Ubyivovk
,
S.
Shapenkov
,
A.
Bondarenko
,
P.
Saring
, and
M.
Seibt
,
J. Appl. Phys.
123
,
161427
(
2018
).
16.
D.
Won
,
X.
Weng
, and
J. M.
Redwing
,
J. Appl. Phys.
108
,
093511
(
2010
).
17.
A. Y.
Polyakov
,
N. B.
Smirnov
,
I. V.
Shchemerov
,
E. B.
Yakimov
,
E. E.
Yakimov
,
K. C.
Kim
, and
I.-H.
Lee
,
ECS J. Solid State Sci. Technol.
7
,
Q80
(
2018
).
18.
I.-H.
Lee
,
A. Y.
Polyakov
,
N. B.
Smirnov
,
E. B.
Yakimov
,
S. A.
Tarelkin
,
A. V.
Turutin
,
I. V.
Shemerov
, and
S. J.
Pearton
,
J. Appl. Phys.
119
,
205109
(
2016
).
19.
I.-H.
Lee
,
A. Y.
Polyakov
,
E. B.
Yakimov
,
N. B.
Smirnov
,
I. V.
Shchemerov
,
S. A.
Tarelkin
,
S. I.
Didenko
,
K. I.
Tapero
,
R. A.
Zinovyev
, and
S. J.
Pearton
,
Appl. Phys. Lett.
110
,
112102
(
2017
).
20.
P. S.
Vergeles
,
V. I.
Orlov
,
A. Y.
Polyakov
,
E. B.
Yakimov
,
T.
Kim
, and
I.-H.
Lee
,
J. Alloys Comp.
776
,
181
(
2019
).
21.
A. Y.
Polyakov
,
N. B.
Smirnov
,
I.-H.
Lee
, and
S. J.
Pearton
,
J. Vac. Sci. Technol. B
33
,
061203
(
2015
).
22.
A. Y.
Polyakov
,
C.
Haller
,
N. B.
Smirnov
,
A. S.
Shiko
,
I. V.
Shchemerov
,
S. V.
Chernykh
,
L. A.
Alexanyan
,
P. B.
Lagov
,
Y. S.
Pavlov
,
J.-F.
Carlin
,
M.
Mosca
,
R.
Butté
,
N.
Grandjean
, and
S. J.
Pearton
,
J. Appl. Phys.
126
,
125708
(
2019
).
23.
A. Y.
Polyakov
and
I.-H.
Lee
,
Mat. Sci. Eng.
94
,
1
(
2015
).
24.
Capacitance Spectroscopy of Semiconductors
, edited by
J. V.
Li
and
G.
Ferrari
(
Pan Stanford Publishing Pte Ltd
,
Singapore
,
2018
), p.
437
.
25.
J. L.
Lyons
,
A.
Janotti
, and
C. G.
Van de Walle
,
Phys. Rev. Lett.
108
,
156403
(
2012
).
26.
J. L.
Lyons
,
A.
Janotti
, and
C. G.
Van de Walle
,
Phys. Rev. B
89
,
035204
(
2014
).
27.
T.
Narita
,
Y.
Tokuda
,
T.
Kogiso
,
K.
Tomita
, and
T.
Kachi
,
J. Appl. Phys.
123
,
161405
(
2018
).
28.
Q.
Yan
,
A.
Janotti
,
M.
Scheffler
, and
C. G.
Van de Walle
,
Appl. Phys. Lett.
100
,
142110
(
2012
).
29.
M. A.
Reshchikov
and
H.
Morkoç
,
J. Appl. Phys.
97
,
061301
(
2005
).
30.
S. F.
Chichibu
,
K.
Shima
,
K.
Kojima
,
S.
Takashima
,
M.
Edo
,
K.
Ueno
,
S.
Ishibashi
, and
A.
Uedono
,
Appl. Phys. Lett.
112
,
211901
(
2018
).
31.
E. B.
Yakimov
,
S. S.
Borisov
, and
S. I.
Zaitsev
,
Semiconductors
41
,
411
(
2007
).
32.
E. B.
Yakimov
,
J. Alloys Comp.
627
,
344
(
2015
).