In this work, we demonstrate that β-Ga2O3 shows orientation-dependent polarized photoluminescence (PL) emission and give a comprehensive insight into gallium oxide's PL spectral properties. We characterized the polarization and spectral dependencies of both the incident and emitted light for (−201) unintentionally doped (UID) as well as (−201) and (010) Sn-doped and Fe-doped crystals. We observed for UID and Sn-doped samples that the electron to self-trapped hole and native defect-related emission bands are linearly polarized with polarized emission intensities ordered as E || c (and c*) > E || a (and a*) > E || b. Furthermore, the spectral shape of emission does not change between the UID and Sn-doped samples; instead, the Sn-doping quenches the total PL spectral intensity. For Fe-doped samples, polarized red emission caused by unintentional Cr3+ doping generates emission intensities ordered E || b > E || c (and c*) > E || a (and a*). It is also observed that in some circumstances, for some doped crystals, the PL spectra can show variations not only in intensity but also in spectral shape along different polarization directions. As an example, the PL emission band for emission along c is blueshifted relative to that along a in Sn-doped β-Ga2O3.

Gallium oxide (β-Ga2O3) is a material of interest in both electrical and optical applications due to its unique properties. It is a wide bandgap material with an indirect, fundamental bandgap near 276 nm (4.5 eV), though optical transition energies range from 282 to 253 nm (4.5–4.9 eV), depending on the crystallographic orientation with E||c < E||a < E||b.1–3 Its beta phase, which has a monoclinic structure as seen in Fig. 1, is thermally and chemically stable and can be synthesized through inexpensive melt-growth methods, which allow for large, high-quality bulk substrates. Furthermore, β-Ga2O3 has a large Baliga's figure of merit estimated to be >3000 times that of Si.2 That said, to control the material's electrical and optical properties, understanding and characterizing its defects is a priority.

FIG. 1.

Crystal structure of β-Ga2O3 showing the (−201) plane and the [102] direction along (a) a orientation, (b) b orientation, and (c) c orientation.

FIG. 1.

Crystal structure of β-Ga2O3 showing the (−201) plane and the [102] direction along (a) a orientation, (b) b orientation, and (c) c orientation.

Close modal

One method to do this is through photoluminescence (PL). For β-Ga2O3, the PL spectrum is complex due to the multiple potential defects that could constitute recombination paths as well as vibrational broadening, electron-phonon coupling, etc. This causes deconvolving the spectrum and defect identification to be quite challenging.4 In this context, the literature generally describes the PL spectrum by separating it into regions. The highest energy region is the UV region (centered around 380–390 nm or 3.26–3.18 eV), which is generally described as being caused by the relaxation from the conduction band to self-trapped holes (STHs).5–11 The next region is called the blue region (centered around 410–480 nm or 3.02–2.58 eV) and is believed to be caused by the recombination of donor-acceptor pairs or donor to STH.5,9,11–14 In general, gallium and oxygen vacancies or interstitials are typically identified as the cause of the blue emission, though it has also been found that extended defects can cause blue emission in poor-quality samples.4,15 The next region is the green region (around ∼500–600 nm or 2.48–2.07 eV), which is not as well understood as the UV and blue regions but is typically attributed to defects, recombination of donor-acceptor pairs caused by gallium and oxygen vacancies or interstitials using OII or OI to GaII, though dopants may also cause this.9,11,16,17 Last, there is also a red region (within 650–850 nm or 1.91–1.46 eV) that is typically attributed to Cr3+, which is an unintentional dopant and comes from the iridium crucible during boule growth.18–20 It is theorized that Cr3+ causes a red emission only when the sample is insulating.21 

However, previous literature for β-Ga2O3 typically only discusses results for unpolarized and excitation-polarized PL. In this context, there have been a handful of prior works analyzing polarization properties of the PL emission. From our understanding, Yamaga et al. first measured polarized emission using unintentionally doped (UID) and Si-doped β-Ga2O3, whose results showed weak emission from b orientation.6 Later, Cho et al. characterized the polarized emission in (−201) bulk UID β-Ga2O3 and showed dominant emission from [102] as well as weak emission from b. However, this measurement was a transmission measurement, meaning the emitted PL spectra would transmit through the bulk of the sample, potentially affecting it.10 Here, we verify their results and explore the polarization of emission on a broader set of bulk crystals containing a wide set of dopants and orientations, demonstrating that the PL emission in β-Ga2O3 is polarization-dependent due to anisotropy of the crystal structure. As such, our results give a comprehensive insight into gallium oxide's PL spectral properties.

Several samples of bulk β-Ga2O3 substrates purchased from Novel Crystal Technology were used in this study. The samples include (−201) UID β-Ga2O3, (−201) Sn-doped β-Ga2O3, (010) Sn-doped β-Ga2O3, (−201) Fe-doped β-Ga2O3, and (010) Fe-doped β-Ga2O3. Details on the substrate properties can be found on the Novel Crystal Technology (NCT) website.22 Additionally, elemental analysis for Ga2O3 single crystal growths can be found in previous literature19,23–28 and PL characterization of (−201) UID β-Ga2O3 was reported in our previous work.4 

Photoluminescence was performed using ultrafast (fs) pulses from a wavelength-tunable (690–1040 nm or 1.8–1.2 eV) Ti:sapphire (Coherent Chameleon Vision Ultra) laser, which passed through a third-harmonic generator (Coherent Harmonics). The laser was then polarized using a linear polarizer (Glan-Laser alpha-BBO polarizer prism, 210–450 nm or 5.9–2.76 eV) followed by a zero-order half-wave plate to control the excitation polarization angle of the laser.4 The polarized laser then reflected off a beam splitter (Semrock 266 nm beam splitter) and through a lens followed by a pinhole hitting the sample and exciting it. This was followed by emission from the sample, which was focused back through the pinhole and lens and transmitted through the beam splitter into a polarizer (Thorlabs). The emission then went through another lens, focusing the emission into an optical fiber. The PL spectra were collected at room temperature using an optical fiber connected to a broadband spectrometer in the range of 300–800 nm or 4.13–1.55 eV (Avasoft AvaSpec dual-channel spectrometer).

Samples were measured using an excitation of 240 nm (5.16 eV) and 267 nm (4.64 eV). Excitation polarization was altered from 0° to 180° with a 15° step size. The emission polarization was measured from 0° to 180° with a 15° step size as well (or left out of the setup for unpolarized emission measurements). The collected data were corrected to remove the response caused by the spectrometer and eliminate the grating and detector response to extract the response of the sample itself. The data collected were corrected for the spectrometer spectra using the Ocean Insight HL-3 plus visible-near infrared (VIS-NIR) light source. The spectrum for this light source was calibrated by the Shanghai Calibration Laboratory. The calibrated blackbody radiant energy spectrum from the laboratory was divided by the spectrum collected by our spectrometer to get a correction factor. This correction factor was then applied to the spectra collected for all samples. All the corrected spectra were then normalized for the integration time.

Photoluminescence (PL) was first done on (−201) UID bulk β-Ga2O3. Initially, polarized excitation along b and [102] was done with total, unpolarized emission detected [Fig. 2(a)]. It was found that the total, unpolarized emission depends on the excitation's polarization. Furthermore, the maximum intensity between orientations changes depending on the excitation wavelength. This is because of the anisotropy of the absorption edge for β-Ga2O3, with E || b > E || [102] (energy-wise).3,29 As such, the PL spectral intensity will depend on how the excitation is polarized to β-Ga2O3 orientation and how close this is to the absorption edge for that orientation. This is seen in Fig. 2(b), which was made by integrating the PL spectra of E || b and subtracting it from the integrated PL spectra of E || [102]. The difference is then divided by the sum of the two integrated orientations. It is seen that the PL spectrum has maximum intensity along E || b below 260 nm (4.77 eV) excitation. The maximum intensity flips to E || [102] when the excitation is above 260 nm (4.77 eV) excitation.

FIG. 2.

PL spectra of (−201) UID β-Ga2O3. (a) PL spectra at different polarization excitations while emission is not polarized. Excitation is at 240 nm (5.16 eV). The inset is the same plot normalized. (b) PLE of the integrated PL spectra for Ein || b – Ein || [102] and then divided by the sum of the two. The points circled indicate the excitation wavelengths used for PL spectral measurements. (c) PL spectra at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || b. The inset is the same plot normalized. (d) PL spectra at different polarization emissions excited at 267 nm (4.64 eV) with excitation Ein || [102]. Missing data in the graph are from the 2nd harmonic of the excitation laser. The inset is the same plot normalized.

FIG. 2.

PL spectra of (−201) UID β-Ga2O3. (a) PL spectra at different polarization excitations while emission is not polarized. Excitation is at 240 nm (5.16 eV). The inset is the same plot normalized. (b) PLE of the integrated PL spectra for Ein || b – Ein || [102] and then divided by the sum of the two. The points circled indicate the excitation wavelengths used for PL spectral measurements. (c) PL spectra at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || b. The inset is the same plot normalized. (d) PL spectra at different polarization emissions excited at 267 nm (4.64 eV) with excitation Ein || [102]. Missing data in the graph are from the 2nd harmonic of the excitation laser. The inset is the same plot normalized.

Close modal

That said, when a polarizer was placed in front of the spectrometer to determine the orientation of the light emission with respect to the sample, emission polarization was observed to be always most intense along [102], no matter what excitation is used or excitation polarization is used. This is in agreement with the results seen in previous literature.6,10 In fact, there are very little, if not any, PL spectra coming from b as seen in Figs. 2(c) and 2(d). Another interesting trait is that the overall PL spectral shape does not change, no matter what polarization (excitation or emission) is used. This points out likely to the same radiative-recombination paths occurring along different orientations.

While we are uncertain of the reason for this polarized emission dependence, we theorize that it is a phenomenon likely arising from the intrinsic crystal structure. The conduction band is made up of Gallium orbitals, while the valence band is made up of oxygen 2p orbitals. Self-trapped holes (STHs) are also connected to the oxygen atom in the crystal structure.30–32 The radiative recombination of an electron from the conduction band to STH could be an electron moving from a gallium site (atom or vacancy) to an oxygen site (atom or vacancy).10,33 The path the electron takes during recombination lies along [102], which could be along the a, c, or a combination of both. So, the lack of PL emission along b means the electron's path is not along that vector. Previous literature on STHs discus that the threefold O(I) and O(II) sites, which run along a and c, respectively, are more favorable for STHs. While the fourfold O(III) site along b is less favorable energy-wise for STHs.30–32 This means that, energy-wise, the relaxation path along the O(I) and O(II) are more favorable than O(III),30–32 or more favorable along a and c while less favorable along b.

(−201) Sn-doped bulk β-Ga2O3 has similar results to the (−201) UID bulk β-Ga2O3 sample. Figure 3(a) shows that the excitation polarization dependence is weaker for Sn-doped β-Ga2O3 than with the UID β-Ga2O3. Doping will cause a shift in the absorption edge energy along different orientations.34,35 This is likely why, when excited at 240 nm (5.16 eV), the excitation polarization dependence is weaker due to the absorption edge for b-orientation redshifting. That is also why the excitation polarization dependence is stronger when excited at 267 nm (4.64 eV) in Fig. 3(c) Also, similar to the UID β-Ga2O3 sample, the excitation polarization dependence flips depending on excitation energy, with E || b more intense at lower wavelengths (higher energies) and E || [102] more intense at higher wavelengths (lower energies). This is again due to the location of the absorption edge for different orientations; it is just redshifted due to the Sn-doping.

FIG. 3.

PL spectra of (−201) Sn-doped β-Ga2O3. (a) PL spectra at different polarization excitations and emission is not polarized, excited at 240 nm (5.16 eV). (b) PL spectra at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || b. (c) PL spectra at different polarization excitations and emission is not polarized, excited at 267 nm (4.64 eV). (d) PL spectra at different polarization emissions excited at 267 nm (4.64 eV) with excitation Ein || [102]. Missing data were purposefully emitted because of the 2nd harmonic of the excitation laser. All insets are normalized versions of the plot within which they are embedded.

FIG. 3.

PL spectra of (−201) Sn-doped β-Ga2O3. (a) PL spectra at different polarization excitations and emission is not polarized, excited at 240 nm (5.16 eV). (b) PL spectra at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || b. (c) PL spectra at different polarization excitations and emission is not polarized, excited at 267 nm (4.64 eV). (d) PL spectra at different polarization emissions excited at 267 nm (4.64 eV) with excitation Ein || [102]. Missing data were purposefully emitted because of the 2nd harmonic of the excitation laser. All insets are normalized versions of the plot within which they are embedded.

Close modal

Further similarities between UID β-Ga2O3 and Sn-doped β-Ga2O3 can be seen by looking at where the emission originates from. This is observed in Figs. 3(b) and 3(d) where strong emission along [102] occurs and there is a very weak emission along b. This arises at any excitation energy and any excitation polarization used, just like with UID β-Ga2O3.

Taking a closer look by plotting both the (−201) UID bulk β-Ga2O3 sample and (−201) Sn-doped bulk β-Ga2O3 together in Fig. 4 shows that the Sn-doped sample has a weaker emission compared to the UID sample. However, the actual shape of the overall PL emission does not change between the two samples. This is in agreement with the PL signature for pure and Si-doped Ga2O3 samples measured by Yamaga et al.6 This means that the doping did not generate any new radiative-recombination paths nor cause a change in the probability of the radiative-recombination paths. Instead, Sn appears to have quenched the overall radiative-recombination paths in the sample. Specifically, the ratio of emission along [102] between exciting along the E || [102] and E || b for UID is much larger than for the Sn-doped sample by a factor of around 2.

FIG. 4.

PL spectra comparison of (−201) UID and Sn-doped Ga2O3. (a) Comparison of samples when 240 nm (5.16 eV) excitation polarization is changed, and emission is measured along [102]. The inset is the same plot normalized. (b) Comparison of samples when emission along different orientations is changed, and 240 nm (5.16 eV) excitation polarization is kept E || b. The inset is the same plot normalized.

FIG. 4.

PL spectra comparison of (−201) UID and Sn-doped Ga2O3. (a) Comparison of samples when 240 nm (5.16 eV) excitation polarization is changed, and emission is measured along [102]. The inset is the same plot normalized. (b) Comparison of samples when emission along different orientations is changed, and 240 nm (5.16 eV) excitation polarization is kept E || b. The inset is the same plot normalized.

Close modal

Looking at a different crystal orientation, Fig. 5 shows the PL emission results of a (010) Sn-doped β-Ga2O3 sample. Figure 5(a) shows the results for the total emission when the excitation at 240 nm (5.16 eV) is polarized. As seen, the PL spectrum does not change between excitation polarizations. This is again due to the Sn-doping, which causes a red shift in the absorption edge, allowing the emission along a to match c. To determine the orientation of the sample, transmission was done to find the bandgap and, therefore, sample orientation as seen in Fig. S1.36 

FIG. 5.

PL spectra of Sn-doped β-Ga2O3. (a) PL spectra of (010) Sn-doped β-Ga2O3 at different polarization excitations, and emission is not polarized, excited at 240 nm (5.16 eV). (b) PL spectra of (010) Sn-doped β-Ga2O3 at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || b. The inset is the same plot normalized. (c) Comparison of (010) and (−201) Sn-doped sample orientations. The inset is the same plot normalized.

FIG. 5.

PL spectra of Sn-doped β-Ga2O3. (a) PL spectra of (010) Sn-doped β-Ga2O3 at different polarization excitations, and emission is not polarized, excited at 240 nm (5.16 eV). (b) PL spectra of (010) Sn-doped β-Ga2O3 at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || b. The inset is the same plot normalized. (c) Comparison of (010) and (−201) Sn-doped sample orientations. The inset is the same plot normalized.

Close modal

Interestingly, there is emission polarization dependence between the sample orientations, as seen in Fig. 5(b). c emits more strongly than a. Furthermore, in this sample, contrary to the (−201), we observe differences in PL spectra for different emission directions. As seen in the Fig. 5(b) inset, there is a blueshift when comparing the emission from a to the emission from c. This suggests that the radiative recombination that occurs along a generates a more significant fraction of higher energy photons. Furthermore, the low energy tail along a is less pronounced. That said, the overall PL shape remains very similar in shape with dominant emission in the UV, which is characteristic of samples with low density of structural defects.4 So, when comparing light being emitted along different polarizations, these differences in emission spectra point to radiative-recombination paths being different for emission between a and c.

Figure 5(c) compares a (−201) Sn-doped β-Ga2O3 sample with a (010) Sn-doped β-Ga2O3 sample. When comparing the overall intensity, we see that the emission along c > [102] > a > b, which is in the same order as the wavelength absorption edge for each orientation. These results match with previous literature on Si-doped β-Ga2O3 samples.6 

Whenever a β-Ga2O3 sample is insulating, it causes a quenching of the UV, blue, and green PL emission. Instead, unintentionally doped Cr3+ which originates from the crucible during growth causes the presence of a red emission within the Ga2O3 sample. The Fe-doping pins the Fermi level causing Cr3+, which has orbital transitions that generate red emission when excited.18,19,21 This is seen in Fig. 6 where the PL emission for Fe-doped β-Ga2O3 is plotted and shows a red emission.

FIG. 6.

PL spectra of (−201) Fe-doped β-Ga2O3. (a) PL spectra at different polarization excitations and emission is not polarized, excited at 240 nm (5.16 eV). (b) PL spectra at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || [102]. (c) PL spectra at different polarization excitations and emission is not polarized, excited at 267 nm (4.64 eV). (d) PL spectra at different polarization emissions excited at 267 nm (4.64 eV) with excitation Ein || b. All insets are of the same plot normalized.

FIG. 6.

PL spectra of (−201) Fe-doped β-Ga2O3. (a) PL spectra at different polarization excitations and emission is not polarized, excited at 240 nm (5.16 eV). (b) PL spectra at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || [102]. (c) PL spectra at different polarization excitations and emission is not polarized, excited at 267 nm (4.64 eV). (d) PL spectra at different polarization emissions excited at 267 nm (4.64 eV) with excitation Ein || b. All insets are of the same plot normalized.

Close modal

Figures 6(a) and 6(c) are the plots of the polarized excitation and total, un-polarized emission at 240 nm (5.16 eV) and 267 nm (4.64 eV), respectively. They both show polarization dependence, which is dependent on excitation. Furthermore, the intensity of the emission generated by the polarized excitation is also dependent on the excitation energy. Polarization transmission measurements were performed to determine the sample orientation as seen in Fig. S2.36 Exciting along [102] will generate a stronger emission at 240 nm (5.16 eV) excitation, while exciting along b will generate a stronger emission at 267 nm (4.64 eV) excitation. This is the opposite of what was observed for UID and Sn-doped β-Ga2O3, though the recombination process differs for Fe-doped β-Ga2O3. There might be a correlation between the emission from Cr3+ and its position within the crystal structure versus the emission from band transitions. For example, the Cr3+ is more likely to occur along orientations where emission from band transitions is less likely.

In Figs. 6(b) and 6(d) are the plots of the polarized emission at 240 nm (5.16 eV) and 267 nm (4.64 eV), respectively. As seen, there is weak polarization dependence from the emission with E || b generating a slightly stronger intensity. Again, this is the opposite of what was observed for UID and Sn-doped β-Ga2O3. Furthermore, there was no change within the PL shape when the emission plots were normalized, suggesting no change in radiative-recombination paths between orientations.

(010) Fe-doped β-Ga2O3 shows very different results. As seen in Figs. 7(a) and 7(c), which are the plots of the polarized excitation and total, un-polarized emission at 240 nm (5.16 eV) and 267 nm (4.64 eV), respectively, there is no or very weak polarization dependence due to excitation polarization. But, as seen in Figs. 7(b) and 7(d), which are the plots of the polarized emission at 240 nm (5.16 eV) and 267 nm (4.64 eV), respectively, there is polarization dependence depending on the emission's orientation from the sample. Again, polarization transmission measurements were performed to determine the sample orientation as seen in Fig. S3.36 The emission from c has a dominant emission around 690 nm (1.8 eV). But that peak is quenched when looking at the emission from a. That said, it is worth mentioning that the loss of the 690 nm (1.8 eV) peak might not always occur. When a lamp and monochromator were used for excitation rather than a pulsed laser, the polarized emission between Eout || a and Eout || c both showed a 690 nm (1.8 eV) peak. The causes behind this observation are currently unknown but are being investigated.

FIG. 7.

PL spectra of (010) Fe-doped β-Ga2O3. (a) PL spectra at different polarization excitations and emission is not polarized, excited at 240 nm (5.16 eV). (b) PL spectra at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || a. (c) PL spectra at different polarization excitations and emission is not polarized, excited at 267 nm (4.64 eV). (d) PL spectra at different polarization emissions excited at 267 nm (4.64 eV) with excitation Ein || c.

FIG. 7.

PL spectra of (010) Fe-doped β-Ga2O3. (a) PL spectra at different polarization excitations and emission is not polarized, excited at 240 nm (5.16 eV). (b) PL spectra at different polarization emissions excited at 240 nm (5.16 eV) with excitation Ein || a. (c) PL spectra at different polarization excitations and emission is not polarized, excited at 267 nm (4.64 eV). (d) PL spectra at different polarization emissions excited at 267 nm (4.64 eV) with excitation Ein || c.

Close modal

We characterized bulk, single crystal samples of (−201) and (010) β-Ga2O3 of various dopants using polarized excitation and looking at the polarization of the emission. UID and Sn-doped β-Ga2O3 both have a strong emission along c, which is slightly weaker along a and very weak along b. There is no change in the PL shape between orientations for (−201) UID and Sn-doped β-Ga2O3 samples. That said, (010) Sn-doped β-Ga2O3 shows a blueshift for the emission along a compared to c, suggesting a difference in radiative-recombination paths.

Fe-doped β-Ga2O3 exhibits a red emission instead of a UV, blue, or green emission. This red emission is caused by Cr3+, which is unintentionally doped during the boule growth. For the (−201) Fe-doped β-Ga2O3 sample, there is weak polarization dependence for the emission, with b generating a slightly stronger emission intensity. For (010) Fe-doped β-Ga2O3, under our experimental conditions, we observe a strong polarization dependence with a loss of the peak around 690 nm for the emission along a.

This material is based upon work supported by the Air Force Office of Scientific Research under Award No. FA9550-21-1-0507 (program manager: Dr. Ali Sayir). Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the United States Air Force.

The authors have no conflicts to disclose

Jacqueline Cooke: Data curation (equal); Formal analysis (equal); Investigation (equal); Writing – original draft (equal); Writing – review & editing (equal). Minhan Lou: Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal). Michael A. Scarpulla: Conceptualization (equal); Funding acquisition (equal); Methodology (equal); Project administration (equal); Supervision (equal); Visualization (equal); Writing – review & editing (equal). Berardi Sensale-Rodriguez: Conceptualization editing (equal); Funding acquisition editing (equal); Project administration editing (equal); Resources editing (equal); Supervision editing (equal); Visualization editing (equal); Writing – review & editing editing (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
J. A.
Spencer
,
A. L.
Mock
,
A. G.
Jacobs
,
M.
Schubert
,
Y.
Zhang
, and
M. J.
Tadjer
,
Appl. Phys. Rev.
9
,
011315
(
2022
).
2.
S. J.
Pearton
,
J.
Yang
,
P. H.
Cary
,
F.
Ren
,
J.
Kim
,
M. J.
Tadjer
, and
M. A.
Mastro
,
Appl. Phys. Rev.
5
,
011301
(
2018
).
3.
F.
Ricci
,
F.
Boschi
,
A.
Baraldi
,
A.
Filippetti
,
M.
Higashiwaki
,
A.
Kuramata
,
V.
Fiorentini
, and
R.
Fornari
,
J. Phys.: Condens. Matter
28
,
224005
(
2016
).
4.
J.
Cooke
,
P.
Ranga
,
J.
Jesenovec
,
J. S.
McCloy
,
S.
Krishnamoorthy
,
M. A.
Scarpulla
, and
B.
Sensale-Rodriguez
,
Sci. Rep.
12
,
3243
(
2022
).
5.
K.
Shimamura
,
E. G.
Víllora
,
T.
Ujiie
, and
K.
Aoki
,
Appl. Phys. Lett.
92
,
201914
(
2008
).
6.
M.
Yamaga
,
T.
Ishikawa
,
M.
Yoshida
,
T.
Hasegawa
,
E. G.
Villora
, and
K.
Shimamura
,
Phys. Status Solidi C
8
,
2621
(
2011
).
7.
Q. D.
Ho
,
T.
Frauenheim
, and
P.
Deák
,
Phys. Rev. B
97
,
115163
(
2018
).
9.
Y.
Berencén
,
Y.
Xie
,
M.
Wang
,
S.
Prucnal
,
L.
Rebohle
, and
S.
Zhou
,
Semicond. Sci. Technol.
34
,
035001
(
2019
).
10.
J.
Bin Cho
,
G.
Jung
,
K.
Kim
,
J.
Kim
,
S. K.
Hong
,
J. H.
Song
, and
J. I.
Jang
,
J. Phys. Chem. C
125
,
1432
(
2021
).
11.
N.
Zhang
,
I. M.
Kislyakov
,
C.
Xia
,
H.
Qi
,
J.
Wang
, and
H. F.
Mohamed
,
Opt. Express
29
,
18587
(
2021
).
12.
H.
Wakai
,
Y.
Sinya
, and
A.
Yamanaka
,
Phys. Status Solidi C
8
,
537
(
2011
).
13.
T.
Onuma
,
S.
Fujioka
,
T.
Yamaguchi
,
M.
Higashiwaki
,
K.
Sasaki
,
T.
Masui
, and
T.
Honda
,
Appl. Phys. Lett.
103
,
041910
(
2013
).
14.
O. M.
Bordun
,
B. O.
Bordun
,
I. Y.
Kukharskyy
, and
I. I.
Medvid
,
J. Appl. Spectrosc.
84
,
46
(
2017
).
15.
J.
Cooke
et al,
Opt. Mater. Express
12
,
4341
(
2022
).
16.
H.
Oda
,
N.
Kimura
,
D.
Yasukawa
,
H.
Wakai
, and
A.
Yamanaka
,
Phys. Status Solidi A
214
,
1600670
(
2017
).
17.
C.
Liu
et al,
Semicond. Sci. Technol.
33
,
095022
(
2018
).
18.
C.
Remple
,
L. M.
Barmore
,
J.
Jesenovec
,
J. S.
McCloy
, and
M. D.
McCluskey
,
J. Vac. Sci. Technol. A
41
,
022702
(
2023
).
19.
R.
Sun
,
Y. K.
Ooi
,
P. T.
Dickens
,
K. G.
Lynn
, and
M. A.
Scarpulla
,
Appl. Phys. Lett.
117
,
052101
(
2020
).
20.
Y.
Tokida
and
S.
Adachi
,
J. Appl. Phys.
112
,
063522
(
2012
).
21.
V.
Vasyltsiv
,
A.
Luchechko
,
Y.
Zhydachevskyy
,
L.
Kostyk
,
R.
Lys
,
D.
Slobodzyan
,
R.
Jakieła
,
B.
Pavlyk
, and
A.
Suchocki
,
J. Vac. Sci. Technol. A
39
,
033201
(
2021
).
22.
Novel Crystal Technologies Inc.
, Datasheets: standard specifications for β-Ga2O3 substrates., see https://www.novelcrystal.co.jp/eng/wp-content/uploads/2022/04/5cfefe5fd718f6d5ef0f66c39ee4ebe3.pdf (last accessed 11 December 2023)
23.
Z.
Galazka
et al,
J. Cryst. Growth
404
,
184
(
2014
).
24.
Z.
Galazka
et al,
ECS J. Solid State Sci. Technol.
6
,
Q3007
(
2017
).
25.
J. R.
Ritter
,
K. G.
Lynn
, and
M. D.
McCluskey
,
J. Appl. Phys.
126
,
225705
(
2019
).
26.
Z.
Galazka
,
J. Appl. Phys.
131
,
031103
(
2022
).
27.
J. S.
McCloy
,
J.
Jesenovec
,
B. L.
Dutton
,
C.
Pansegrau
,
C.
Remple
,
M. H.
Weber
,
S.
Swain
,
M.
McCluskey
, and
M.
Scarpulla
,
Proc. SPIE
12002
,
1200205
(
2022
).
28.
A.
Kuramata
,
K.
Koshi
,
S.
Watanabe
,
Y.
Yamaoka
,
T.
Masui
, and
S.
Yamakoshi
,
Jpn. J. Appl. Phys.
55
,
1202A2
(
2016
).
29.
J.
Cooke
,
P.
Ranga
,
A.
Bhattacharyya
,
X.
Cheng
,
Y.
Wang
,
S.
Krishnamoorthy
,
M. A.
Scarpulla
, and
B.
Sensale-Rodriguez
,
J. Vac. Sci. Technol. A
41
,
013406
(
2023
).
30.
Y. K.
Frodason
,
K. M.
Johansen
,
L.
Vines
, and
J. B.
Varley
,
J. Appl. Phys.
127
,
075701
(
2020
).
31.
B. E.
Kananen
,
N. C.
Giles
,
L. E.
Halliburton
,
G. K.
Foundos
,
K. B.
Chang
, and
K. T.
Stevens
,
J. Appl. Phys.
122
,
215703
(
2017
).
32.
J. B.
Varley
,
A.
Janotti
,
C.
Franchini
, and
C. G.
Van De Walle
,
Phys. Rev. B
85
,
081109(R)
(
2012
).
33.
P.
Deák
,
Q.
Duy Ho
,
F.
Seemann
,
B.
Aradi
,
M.
Lorke
, and
T.
Frauenheim
,
Phys. Rev. B
95
,
075208
(
2017
).
34.
X.
Ma
,
Y.
Zhang
,
L.
Dong
, and
R.
Jia
,
Res. Phys.
7
,
1582
(
2017
).
35.
S.
Gao
,
W.
Li
,
J.
Dai
,
Q.
Wang
, and
Z.
Suo
,
Mater. Res. Express
8
,
025904
(
2021
).
36.
See supplementary material online for supporting figure content.

Supplementary Material