We report the continuous Si doping in β-Ga2O3 epitaxial films grown by plasma-assisted molecular beam epitaxy through the use of a valved effusion cell for the Si source. Secondary ion mass spectroscopy results exhibit that the Si doping profiles in β-Ga2O3 are flat and have sharp turn-on/off depth profiles. The Si doping concentration was able to be controlled by either varying the cell temperatures or changing the aperture of the valve of the Si effusion cell. High crystal quality and smooth surface morphologies were confirmed on Si-doped β-Ga2O3 epitaxial films grown on (010) and (001) substrates. The electronic properties of Si-doped (001) β-Ga2O3 epitaxial film showed an electron mobility of 67 cm2/Vs at the Hall concentration of 3 × 1018 cm−3.

β-Ga2O3 has attracted enormous attention due to its ultra-wide bandgap (∼4.8 eV)1 and high breakdown field (∼8 MV/cm).2 These exceptional properties make this material a promising candidate for next-generation power electronics. Compared to other ultra-wide bandgap semiconductors such as AlN and diamond, β-Ga2O3 has shallow hydrogenic donors3 such as Ge,4–6 Sn,7–9 and Si.10–12 In addition, β-Ga2O3 is the only ultra-wide bandgap semiconductor that can be mass-produced as high-purity bulk substrates by the edge-defined fed growth (EFG) method13 and the Czochralski method.14 The availability of cost-effective bulk substrates allows the growth of β-Ga2O3 epitaxial films with superior crystal quality,15–17 which leads to an exceptionally low background carrier concentration.18,19 To measure the potential of a material for power applications, Baliga’s Figure of Merit (BFOM) is widely used.20 While BFOM is based on the assumption that the donors are hydrogenic and the carrier concentrations are not affected by background compensation, this FOM can be modified to reflect the real activation energies of donors and background carrier concentrations. As a result, β-Ga2O3 significantly outperforms other wide-bandgap semiconductor materials for power applications when the doping abilities and materials purity limitation are taken into consideration.21 The crystal orientation of β-Ga2O3 bulk substrates for epitaxial growth is also critical toward the application of β-Ga2O3-based power electronics. Vertical power electron devices are typically large and, therefore,22,23 the scalability of the β-Ga2O3 bulk substrates is the key to determining which crystal orientation is most efficient for the epitaxial growth. While the [010] zone is a scalable orientation in β-Ga2O3 crystal system, both the (100) and (001) orientations are in the [010] zone, which means they are suitable for large-scale bulk substrates. Moreover, even 100 mm diameter of (001) β-Ga2O3 oriented wafers are commercially available. Therefore, the (001) orientation is promising for future large-size power electronics applications.24 

Among all the shallow hydrogenic donors of β-Ga2O3, the importance of Si doping in β-Ga2O3 has been proven by its outstanding electrical properties and, therefore, its application to lateral and vertical electronic devices.2,25–27 High-performance enhancement-mode metal-oxide semiconductor field effect transistors (MOSFETs)28 and current aperture vertical electron transistors (CAVETs)29,30 have been achieved by Si-doped β-Ga2O3 films. While ion implantation is a well-established technique to obtain Si-doping at certain areas in β-Ga2O3 films, damage from ion radiation is inevitable during this process.31 The Si-doping for β-Ga2O3 epitaxial film has been widely studied in different growth methods such as metal-organic chemical vapor deposition (MOCVD),11,32 halide vapor phase epitaxy (HVPE),33and pulsed laser deposition (PLD).34 Moreover, for the application of β-(AlxGa1−x)2O3/β-Ga2O3 heterostructures, Si-doping in the β-(AlxGa1−x)2O3 has revealed its importance by showing superior conductivity35 and record-low contact resistance.36 On the other hand, the rapid oxidation of the Si source in molecular beam epitaxy (MBE) has limited the ability of Si doping to δ-doping.37 Effusion cells with an endplate38 and Si suboxide (SiO) source39 have been used to achieve controllable Si-doping in β-Ga2O3 epitaxial films grown by MBE. Due to the existence of a high background O2 partial pressure, typically >10−5 Torr, and the presence of active oxygen (presumably atomic oxygen O) inside the MBE growth environment, the Si source material in the effusion cell needs to be protected from oxidation. In this work, we report the continuous doping of Si in β-Ga2O3 epitaxial films on the (010) and (001) substrates by utilizing a valved effusion cell.

β-Ga2O3 epitaxial films were grown on Fe-doped (010) and (001) semi-insulating β-Ga2O3 substrates from Novel Crystal Technology by plasma-assisted MBE (PAMBE). The growth temperature was varied from 700 to 800 °C. The MBE system was equipped with a Veeco Unibulb radio-frequency plasma source running at 200 W with a foreline pressure of 60 Torr for the generation of O plasma to produce 1%–2% atomic O from pure O2 gas. Prior to the epitaxial growth, the substrate was Ga polished and O polished at a substrate temperature of 800 °C to remove the impurities on the surface of β-Ga2O3 substrates. Ga polishing is proven to be capable of effectively removing impurities and chemical-mechanical polish (CMP)-related damage from the surface and smoothing the surface of both the (010) and (001) substrates.4,40 O polishing by O plasma can further remove impurities to make sure the surface is ready for epitaxial growth. For Si doping, the temperature of the Si effusion cell was varied from 900 to 1200 °C. Figure 1 shows the schematic of a valved effusion cell used in this study. The cell uses a needle valve to control the flux into the growth chamber. The effusion cell of the Si source is designed to prevent the exposure of Si to the high oxygen background in the chamber while providing Si flux at the same time. The needle valve on the top of the crucible minimizes the exposure of the crucible to the environment in the main growth chamber while sustaining the Si flux from the source material. When the needle valve is fully open (100% open), the needle and the body seats are separated so that the Si flux can flow out of the cell through the openings. When the valve is fully closed (0% open), the needle and the body seats have no gap between each other, which results in no flux between the effusion cell and the growth. The opening of the needle valve is controlled by a programmable motor with a precision valve opening of ∼2.5 µm level to further control the Si doping concentration. While the Si flux can be controlled by changing the cell temperature, this requires time for cooling down/heating up and stabilizing the Si flux, whereas the needle valve can instantly change the Si flux to achieve controllable Si doping. During the growth, the Ga flux was varied from 1.0 × 10−7 to 3.0 × 10−7 Torr for conventional growth. We further performed metal-oxide catalyzed epitaxy (MOCATAXY) to improve the crystal quality of Si-doped (010) and (001) β-Ga2O3 epitaxial films.41–43 During the MOCATAXY growth, In flux was provided simultaneously with Ga flux to serve as a catalyst for the β-Ga2O3 epitaxial growth.44 Previously, the MOCATAXY growth showed that the growth rate of β-Ga2O3 epitaxial films was significantly enhanced by introducing an In flux into the MBE growth environment.42 This improvement in the growth rate of β-Ga2O3 epitaxial films is due to the catalytic mechanism of In during the growth of β-Ga2O3. In conventional MBE growth, Ga can only access atomic O as the active oxygen flux. On the other hand, in the MOCATAXY growth environment, In catalyst layer can access more oxygen, presumably molecular O2. Subsequently, a metal-exchange mechanism between Ga and In then allows Ga to preferentially incorporate in the β-Ga2O3 epitaxial film over In. As all the In is replaced by Ga during the growth, pure β-Ga2O3 epitaxial films with high growth rates can be grown at sufficiently high Ga fluxes. In addition, the maximum growth temperature can also be enhanced due to the suppression of Ga2O suboxide desorption by the presence of In in the MBE growth environment.45 Here, the Ga flux was fixed at 2.5 × 10−7 Torr with In flux at 4.0 × 10−7 Torr for catalyzed growth. Prior to the growth of the Si-doped β-Ga2O3 epitaxial films, an unintentionally doped (UID) buffer layer was grown to prevent the impurities on the surface of the substrate from incorporating into the Si-doped layer. Followed by the UID buffer layer, a 2-nm-thick β-(Al1−xGax)2O3 spacer layer with a target al composition of 2% to 3% was grown to enable the observation of thickness fringes in high resolution x-ray diffraction (HRXRD) ω-2θ scans, which helps determine the thickness of the Si-doped layer and characterize the crystal quality. Reflection high energy electron diffraction (RHEED) was used to observe the surface structure during the MBE growth. Surface morphologies of the epitaxial films were characterized by atomic force microscopy (AFM). For the determination of Si concentrations in the β-Ga2O3 epitaxial films, secondary ion mass spectrometry (SIMS) was used with a CAMECA IMS 7f tool. Hall measurements were performed with Ti/Au electrodes deposited on the corners of the samples in Van der Pauw geometry.

FIG. 1.

(a) The schematic of the valved effusion cell at the needle valve aperture of 0% and (b) the needle valve aperture of 100%. The needle valve aperture can be varied continuously from 0% to 100% by a programmable motor to control the Si doping concentration.

FIG. 1.

(a) The schematic of the valved effusion cell at the needle valve aperture of 0% and (b) the needle valve aperture of 100%. The needle valve aperture can be varied continuously from 0% to 100% by a programmable motor to control the Si doping concentration.

Close modal

During the growth, RHEED patterns were monitored to confirm the surface structure of Si-doed β-Ga2O3 epitaxial films. Figures 2(a) and 2(b) show the RHEED patterns of conventional MBE growth and (c), (d) the MOCATAXY growth of Si-doped β-Ga2O3 epitaxial films. The thicknesses of the UID layers and Si-doped layers grown by both growth methods were 100 and 200 nm, respectively. The RHEED images of conventional MBE growth of the UID buffer layer and Si-doped β-Ga2O3 epitaxial films did not show clear streaky patterns indicative of rough surface morphologies. The pattern of the Si-doped film was even more blurry than that of the UID buffer layer, which indicates surface roughening during Si doping of β-Ga2O3 by conventional MBE growth. This suggests that Si can be an anti-surfactant in the β-Ga2O3 epitaxial growth, which is also confirmed by the following HRXRD and AFM results. On the other hand, both the UID buffer layer and the Si-doped film grown by MOCATAXY growth showed clear streaky patterns during the whole growth.

FIG. 2.

RHHED images of conventional MBE growth of β-Ga2O3 (a) UID buffer layer (thickness: 100 nm) and (b) Si-doped film (thickness: 200 nm). The patterns are not clear on the UID films and got even more blurry on the Si-doped film, indicative of surface roughness during the Si-doping of β-Ga2O3. On the other hand, MOCATXY-grown β-Ga2O3 (c) UID buffer layer (thickness: 100 nm) and (d) Si-doped film (thickness: 200 nm) show similarly clear streaky patterns. This suggests that MOCATAXY-grown Si-doped films have smooth surface morphologies and no surface roughening during Si-doping.

FIG. 2.

RHHED images of conventional MBE growth of β-Ga2O3 (a) UID buffer layer (thickness: 100 nm) and (b) Si-doped film (thickness: 200 nm). The patterns are not clear on the UID films and got even more blurry on the Si-doped film, indicative of surface roughness during the Si-doping of β-Ga2O3. On the other hand, MOCATXY-grown β-Ga2O3 (c) UID buffer layer (thickness: 100 nm) and (d) Si-doped film (thickness: 200 nm) show similarly clear streaky patterns. This suggests that MOCATAXY-grown Si-doped films have smooth surface morphologies and no surface roughening during Si-doping.

Close modal

Figure 3(a) shows the SIMS result of a stack of Si-doped β-Ga2O3 layers separated by UID layers grown by MOCATAXY on a (010) β-Ga2O3 substrate. The spike in Si concentration at the film/substrate interface was due to the residual impurities on the surface of the β-Ga2O3 substrate though Ga polishing has reduced the Si concentration by more than one order of magnitude.46 Continuous Si doping in β-Ga2O3 epitaxial films with flat doping profiles and sharp turn-on was achieved by exploiting the valved effusion cell. The Si doping concentration was able to be either controlled by the aperture of the needle valve or the temperature of the effusion cell with an upper limit of ∼1.0 × 1020 cm−3 and a lower limit below the background impurity level of ∼5 × 1016 cm−3, which is the detection limit of SIMS measurement. The temperature of the Si effusion cell was varied from 900 to 1200 °C. At the effusion cell temperature of 1000 °C, the Si concentration was ∼4.0 × 1017 cm−3 with a needle valve opening of 20%. When the needle valve was fully open (100%), the Si concentration increased to ∼2.0 × 1018 cm−3, which is five times larger than that of the needle opening at 20%. This result demonstrates that the linear control of Si concentration in β-Ga2O3 epitaxial film was realized by utilizing the valved effusion cell. On the other hand, some surface segregation effect was observed in the turn-off of the doping profiles, which indicates possible surface riding of Si or a “memory” effect during the growth of β-Ga2O3 epitaxial films on (010) β-Ga2O3 substrates.

FIG. 3.

(a) SIMS profile of a β-Ga2O3 (010) film with different Si doping levels separated by UID layers grown by MOCATAXY growth at a substrate temperature of 800 °C. The Si concentrations were controlled by varying effusion cell temperatures from 900 to 1200 °C and needle valve apertures at 20% opening and 100% opening. The Fe profile is an indicator of the interface between the epitaxial film and the Fe-doped (010) semi-insulating substrate. (b) SIMS profile of a Si-doped β-Ga2O3 (010) SIMS-stack sample with 240 nm Si-doped layers grown at different Si cell temperatures and separated by UID spacer layers. (c) Arrhenius plot of the cell temperatures and Si doping concentrations, showing an apparent activation energy of ∼1.3 eV.

FIG. 3.

(a) SIMS profile of a β-Ga2O3 (010) film with different Si doping levels separated by UID layers grown by MOCATAXY growth at a substrate temperature of 800 °C. The Si concentrations were controlled by varying effusion cell temperatures from 900 to 1200 °C and needle valve apertures at 20% opening and 100% opening. The Fe profile is an indicator of the interface between the epitaxial film and the Fe-doped (010) semi-insulating substrate. (b) SIMS profile of a Si-doped β-Ga2O3 (010) SIMS-stack sample with 240 nm Si-doped layers grown at different Si cell temperatures and separated by UID spacer layers. (c) Arrhenius plot of the cell temperatures and Si doping concentrations, showing an apparent activation energy of ∼1.3 eV.

Close modal

An additional Si-doped (010) β-Ga2O3 SIMS-stack sample was prepared to investigate the Si doping behaviors at high temperatures (>1000 °C) range. Each Si-doped layer was grown for 240 nm with a growth rate of 4 nm/min to confirm the stability of Si doping during the growth. The SIMS profile of Si-doped (010) β-Ga2O3 film shown in Fig. 3(b) reveals the cell temperature dependence of Si doping concentrations. The Si-doped layers with Si effusion cell temperatures at 1200, 1100, and 1000 °C showed flat Si doping profiles and sharp turn-on/off. Moreover, two doping regions were observed from a single Si-doped layer. The first region (Region A) of the doping has a slightly larger Si concentration than the following region (Region B). Considering the needle valve is closed during the growth of the UID layer (prior to the Si-doped layers), excessive gas species, including Si and suboxide SiO, that could be active dopants after being incorporated into the β-Ga2O3 epitaxial films, were accumulated inside the Si effusion cell. Since this phenomenon is more obvious when the Si effusion cell is at a higher temperature, this suggests that the continuous Si doping by the valved cell is conducted by gas species that are able to accumulate in an effusion cell.47,48 Figure 3(c) shows the apparent activation energy extracted from the Arrhenius plot of the cell temperatures and Si doping concentrations. The activation energy showed a value of about 1.3 eV, which was lower than that of the element Si flux (∼5.0 eV)49 or SiO flux (∼3.6 eV).39 However, the identification of the exact gas species coming out of the Si effusion cell needs to be further investigated at actual MBE growth conditions and is beyond the scope of this study.

The crystal quality of the Si-doped (010) β-Ga2O3 epitaxial films was investigated by HRXRD, as shown in Fig. 4. While the film with a film thickness of 200 nm and Si concentration under 1 × 1017 cm−3 show clear thickness fringes, indicating a high quality film and an abrupt substrate–film interface. Thickness fringes were not observed from films with Si concentrations beyond 6 × 1017 cm−3, either by MOCATAXY (film thickness: 200 nm) or by conventional PAMBE (film thickness: 150 nm). In addition, a peak appears next to the (020) β-Ga2O3 peak on the higher angle side at Si concentrations larger than 6 × 1017 cm−3. This can be attributed to the strain introduction while Si is doped into (010) β-Ga2O3 films. The coherency strain calculated from the peak separation50 from the β-Ga2O3 (020) peak showed a value of −0.023%, which indicates Si doping might introduce tensile stress into β-Ga2O3 epitaxial films during the growth. Further work is ongoing to understand the peak separation. Figure 5 shows the AFM images of Si-doped β-Ga2O3 films grown by conventional MBE and MOCATAXY. The surface morphologies of Si-doped β-Ga2O3 films with Si concentrations above the mid-1017 cm−3 range have larger RMS roughness than films with Si concentrations under 1016 cm−3. Hall measurement was performed on Si-doped β-Ga2O3 epitaxial films grown on (010) substrates by conventional MBE growth and MOCATAXY growth. While MOCATAXY-grown Si-doped (010) films showed 1.8 × 1018 cm−3 Hall concentration and 72 cm2/Vs mobility, conventional MBE-grown Si-doped (010) films showed mobility under 10 or even unmeasurable results. This suggests that MOCATAXY growth can enhance the electronic properties by improving the crystal quality of β-Ga2O3 epitaxial films. Therefore, the HRXRD, AFM, and Hall results suggest that Si is potentially an anti-surfactant during the growth of β-Ga2O3 epitaxial films causing microstructure degradation and surface roughening along with the previous studies showing Si is known as an anti-surfactant in GaN growth.51,52 On the other hand, MOCATAXY growth showed its potential to act as a surfactant to counteract the anti-surfactant effect of Si in β-Ga2O3 epitaxial growth.

FIG. 4.

(a) HRXRD 2θ-ω scan of Si-doped β-Ga2O3 (010) films grown by conventional MBE and (b), (c) MOCATAXY. In (b) and (c), the broad peak next to the (020) is the characteristic of the strain introduced to the β-Ga2O3 (010) films, while Si is doped in the >1017 cm−3 range. On the other hand, clear thickness was observed at the Si concentration <1017 cm−3.

FIG. 4.

(a) HRXRD 2θ-ω scan of Si-doped β-Ga2O3 (010) films grown by conventional MBE and (b), (c) MOCATAXY. In (b) and (c), the broad peak next to the (020) is the characteristic of the strain introduced to the β-Ga2O3 (010) films, while Si is doped in the >1017 cm−3 range. On the other hand, clear thickness was observed at the Si concentration <1017 cm−3.

Close modal
FIG. 5.

(a) AFM images of Si-doped β-Ga2O3 (010) films grown by conventional MBE and (b), (c) MOCATAXY. Grooved features along the [001] direction can be observed in the films grown by either conventional MBE growth or MOCCATAXY growth. The value of RMS roughness on MOCATAXY-grown β-Ga2O3 epitaxial films is lower than that of conventional MBE-grown films, which corresponds to the streaky patterns on MOCATAXY-grown films from the RHEED measurement.

FIG. 5.

(a) AFM images of Si-doped β-Ga2O3 (010) films grown by conventional MBE and (b), (c) MOCATAXY. Grooved features along the [001] direction can be observed in the films grown by either conventional MBE growth or MOCCATAXY growth. The value of RMS roughness on MOCATAXY-grown β-Ga2O3 epitaxial films is lower than that of conventional MBE-grown films, which corresponds to the streaky patterns on MOCATAXY-grown films from the RHEED measurement.

Close modal

The Si doping availability of β-Ga2O3 was further investigated on the scalable (001) substrates. Since the (001) substrates were co-loaded with the (010) substrates in the MBE system, the growth conditions of the (001) Si-doped films were identical to those on the (010) substrates, which were shown in Sec. II. By growing with MOCATAXY, co-loaded (010) and (001) β-Ga2O3 films exhibited the best crystal qualities at the same growth conditions, which were confirmed by Mauze et al.42, Figure 6 shows the HRXRD and AFM results of continuous Si-doped β-Ga2O3 films grown on (001) substrates by MOCATAXY growth. Clear thickness fringes in the XRD patterns indicate the Si-doped (001) β-Ga2O3 epitaxial films have superior crystal quality and smooth surface morphology even at the high Si concentration of 1018 cm−3 range. The growth rates on (001) substrates were 5.0 nm/min, which corresponds to the results from previous MOCATAXY-grown UID films.42,43 The RHEED patterns during the MOCATAXY growth of Si-doped β-Ga2O3 epitaxial films showed clear streaky patterns, along with the AFM image showing the smooth surface morphology. The results of Hall measurements of the Si-doped films on (001) substrates are shown in the figures. For the 200-nm-thick Si-doped film with a Hall electron concentration of 3 × 1018 cm−3, the electron mobility was 67 cm2/Vs. The electron mobility of this Si-doped film is superior to previous studies on Sn-doped and Ge-doped β-Ga2O3 epitaxial films grown on (001) substrates.24,42,53

FIG. 6.

(a) HRXRD results of Si-doped β-Ga2O3 films grown on (001) substrates by MOCATAXY growth at a substrate temperature of 800 °C. Thickness fringes were confirmed at a Si concentration of up to 8.3 × 1018 cm−3. (b) Clear streaky RHEED patterns were observed during the growth of Si-doped β-Ga2O3 film on the (001) substrate. (c) AFM image of Si-doped 200-nm-thick β-Ga2O3 films grown on (001) substrates with an RMS roughness as low as 1.3 nm.

FIG. 6.

(a) HRXRD results of Si-doped β-Ga2O3 films grown on (001) substrates by MOCATAXY growth at a substrate temperature of 800 °C. Thickness fringes were confirmed at a Si concentration of up to 8.3 × 1018 cm−3. (b) Clear streaky RHEED patterns were observed during the growth of Si-doped β-Ga2O3 film on the (001) substrate. (c) AFM image of Si-doped 200-nm-thick β-Ga2O3 films grown on (001) substrates with an RMS roughness as low as 1.3 nm.

Close modal

Figure 7 shows the SIMS depth profile of this 100-nm-thick Si-doped β-Ga2O3 epitaxial film grown on the (001) substrate with the Hall measured electron concentration of 9 × 1018 cm−3 and electron mobility of 49 cm2/Vs. This result shows that a flat Si doping profile of 8.3 × 1018 cm−3 with sharp turn-on was achieved throughout the Si-doped (001) β-Ga2O3 epitaxial film. Additionally, this indicates that the doped Si has a sufficiently high ionization efficiency, which is also proven by previous studies suggesting that Si is a shallow donor in the β-Ga2O3 materials system.3 Therefore, the continuous Si-doping capability in β-Ga2O3 epitaxial films grown on (001) substrates was achieved in the range from below the background impurity level to as high as ∼1.0 × 1020 cm−3.

FIG. 7.

SIMS depth profile of a 100-nm-thick Si-doped β-Ga2O3 film grown on a Fe-doped (001) substrate by MOCATAXY growth. A flat Si doping profile was obtained by utilizing the valved effusion cell. The spike in Si concentration at the interface between the UID layer and the substrate is due to the residual impurities on the surface of the substrate.

FIG. 7.

SIMS depth profile of a 100-nm-thick Si-doped β-Ga2O3 film grown on a Fe-doped (001) substrate by MOCATAXY growth. A flat Si doping profile was obtained by utilizing the valved effusion cell. The spike in Si concentration at the interface between the UID layer and the substrate is due to the residual impurities on the surface of the substrate.

Close modal

In conclusion, continuous Si doping in β-Ga2O3 epitaxial films on (010) and (001) substrates grown by PAMBE was achieved by utilizing the valved effusion cell for the Si source material. By varying the effusion cell temperatures and valve aperture, the Si doping concentration can be precisely controlled on a linear scale corresponding to the opening of the needle valve. The clear thickness fringes confirmed from the HRXRD results prove that the Si-doped MOCATAXY-grown β-Ga2O3 epitaxial films exhibited high crystal quality. In addition, AFM images showed that the surface morphology of Si-doped films was smooth with a value of RMS roughness of around 1.0 nm. Controllable continuous doping is crucial to the β-Ga2O3 toward the application of power electronics, especially vertical devices that require a large uniform donor doping area.

The authors acknowledge funding from AFOSR through the GAME MURI Program, Award No. FA9550-18-1-0479, under project manager Dr. Ali Sayir.

The authors have no conflicts to disclose.

Takeki Itoh: Conceptualization (lead); Data curation (lead); Formal analysis (lead); Investigation (lead); Methodology (lead); Visualization (lead); Writing – original draft (lead); Writing – review & editing (lead). Akhil Mauze: Data curation (supporting); Formal analysis (supporting); Investigation (supporting); Methodology (supporting); Validation (supporting). Yuewei Zhang: Formal analysis (supporting); Investigation (supporting); Methodology (supporting); Supervision (supporting); Validation (supporting). James S. Speck: Funding acquisition (lead); Project administration (lead); Resources (lead); Supervision (lead); Writing – review & editing (lead).

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

1.
M.
Orita
,
H.
Ohta
,
M.
Hirano
, and
H.
Hosono
,
Appl. Phys. Lett.
77
,
4166
(
2000
).
2.
M.
Higashiwaki
,
K.
Sasaki
,
A.
Kuramata
,
T.
Masui
, and
S.
Yamakoshi
,
Appl. Phys. Lett.
100
,
013504
(
2012
).
3.
A. T.
Neal
,
S.
Mou
,
S.
Rafique
,
H.
Zhao
,
E.
Ahmadi
,
J. S.
Speck
,
K. T.
Stevens
,
J. D.
Blevins
,
D. B.
Thomson
,
N.
Moser
,
K. D.
Chabak
, and
G. H.
Jessen
,
Appl. Phys. Lett.
113
,
062101
(
2018
).
4.
E.
Ahmadi
,
O. S.
Koksaldi
,
S. W.
Kaun
,
Y.
Oshima
,
D. B.
Short
,
U. K.
Mishra
, and
J. S.
Speck
,
Appl. Phys. Express
10
,
041102
(
2017
).
5.
F.
Alema
,
G.
Seryogin
,
A.
Osinsky
, and
A.
Osinsky
,
APL Mater.
9
,
091102
(
2021
).
6.
P.
Ranga
,
A.
Bhattacharyya
,
L.
Whittaker-Brooks
,
M. A.
Scarpulla
, and
S.
Krishnamoorthy
,
J. Vac. Sci. Technol., A
39
,
030404
(
2021
).
7.
K.
Sasaki
,
A.
Kuramata
,
T.
Masui
,
E. G.
Víllora
,
K.
Shimamura
, and
S.
Yamakoshi
,
Appl. Phys. Express
5
,
035502
(
2012
).
8.
A.
Mauze
,
Y.
Zhang
,
T.
Itoh
,
E.
Ahmadi
, and
J. S.
Speck
,
Appl. Phys. Lett.
117
,
222102
(
2020
).
9.
X.
Du
,
Z.
Li
,
C.
Luan
,
W.
Wang
,
M.
Wang
,
X.
Feng
,
H.
Xiao
, and
J.
Ma
,
J. Mater. Sci.
50
,
3252
(
2015
).
10.
S.
Krishnamoorthy
,
Z.
Xia
,
S.
Bajaj
,
M.
Brenner
, and
S.
Rajan
,
Appl. Phys. Express
10
,
051102
(
2017
).
11.
Z.
Feng
,
A. F. M.
Anhar Uddin Bhuiyan
,
M. R.
Karim
, and
H.
Zhao
,
Appl. Phys. Lett.
114
,
250601
(
2019
).
12.
K.
Sasaki
,
M.
Higashiwaki
,
A.
Kuramata
,
T.
Masui
, and
S.
Yamakoshi
,
Appl. Phys. Express
6
,
086502
(
2013
).
13.
A.
Kuramata
,
K.
Koshi
,
S.
Watanabe
,
Y.
Yamaoka
,
T.
Masui
, and
S.
Yamakoshi
,
Jpn. J. Appl. Phys.
55
,
1202A2
(
2016
).
14.
Z.
Galazka
,
R.
Uecker
,
K.
Irmscher
,
M.
Albrecht
,
D.
Klimm
,
M.
Pietsch
,
M.
Brützam
,
R.
Bertram
,
S.
Ganschow
, and
R.
Fornari
,
Cryst. Res. Technol.
45
,
1229
(
2010
).
15.
F.
Alema
,
Y.
Zhang
,
A.
Osinsky
,
N.
Valente
,
A.
Mauze
,
T.
Itoh
, and
J. S.
Speck
,
APL Mater.
7
,
121110
(
2019
).
16.
Y.
Zhang
,
F.
Alema
,
A.
Mauze
,
O. S.
Koksaldi
,
R.
Miller
,
A.
Osinsky
, and
J. S.
Speck
,
APL Mater.
7
,
022506
(
2019
).
17.
G.
Seryogin
,
F.
Alema
,
N.
Valente
,
H.
Fu
,
E.
Steinbrunner
,
A. T.
Neal
,
S.
Mou
,
A.
Fine
, and
A.
Osinsky
,
Appl. Phys. Lett.
117
,
262101
(
2020
).
18.
F.
Alema
,
Y.
Zhang
,
A.
Osinsky
,
N.
Orishchin
,
N.
Valente
,
A.
Mauze
, and
J. S.
Speck
,
APL Mater.
8
,
021110
(
2020
).
19.
F.
Alema
,
Y.
Zhang
,
A.
Mauze
,
T.
Itoh
,
J. S.
Speck
,
B.
Hertog
, and
A.
Osinsky
,
AIP Adv.
10
,
085002
(
2020
).
20.
B. J.
Baliga
,
J. Appl. Phys.
53
,
1759
(
1982
).
21.
Y.
Zhang
and
J. S.
Speck
,
Semicond. Sci. Technol.
35
,
125018
(
2020
).
22.
M. H.
Wong
and
M.
Higashiwaki
,
IEEE Trans. Electron Devices
67
,
3925
(
2020
).
23.
A. J.
Green
,
J.
Speck
,
G.
Xing
,
P.
Moens
,
F.
Allerstam
,
K.
Gumaelius
,
T.
Neyer
,
A.
Arias-Purdue
,
V.
Mehrotra
,
A.
Kuramata
,
K.
Sasaki
,
S.
Watanabe
,
K.
Koshi
,
J.
Blevins
,
O.
Bierwagen
,
S.
Krishnamoorthy
,
K.
Leedy
,
A. R.
Arehart
,
A. T.
Neal
,
S.
Mou
,
S. A.
Ringel
,
A.
Kumar
,
A.
Sharma
,
K.
Ghosh
,
U.
Singisetti
,
W.
Li
,
K.
Chabak
,
K.
Liddy
,
A.
Islam
,
S.
Rajan
,
S.
Graham
,
S.
Choi
,
Z.
Cheng
, and
M.
Higashiwaki
,
APL Mater.
10
,
029201
(
2022
).
24.
S.-H.
Han
,
A.
Mauze
,
E.
Ahmadi
,
T.
Mates
,
Y.
Oshima
, and
J. S.
Speck
,
Semicond. Sci. Technol.
33
,
045001
(
2018
).
25.
Y.
Zhang
,
C.
Joishi
,
Z.
Xia
,
M.
Brenner
,
S.
Lodha
, and
S.
Rajan
,
Appl. Phys. Lett.
112
,
233503
(
2018
).
26.
E.
Farzana
,
F.
Alema
,
W. Y.
Ho
,
A.
Mauze
,
T.
Itoh
,
A.
Osinsky
, and
J. S.
Speck
,
Appl. Phys. Lett.
118
,
162109
(
2021
).
27.
N. K.
Kalarickal
,
Z.
Xia
,
J. F.
McGlone
,
Y.
Liu
,
W.
Moore
,
A. R.
Arehart
,
S. A.
Ringel
, and
S.
Rajan
,
J. Appl. Phys.
127
,
215706
(
2020
).
28.
M. H.
Wong
,
Y.
Nakata
,
A.
Kuramata
,
S.
Yamakoshi
, and
M.
Higashiwaki
,
Appl. Phys. Express
10
,
041101
(
2017
).
29.
M. H.
Wong
,
K.
Goto
,
H.
Murakami
,
Y.
Kumagai
, and
M.
Higashiwaki
,
IEEE Electron Device Lett.
40
,
431
(
2019
).
30.
M. H.
Wong
,
H.
Murakami
,
Y.
Kumagai
, and
M.
Higashiwaki
,
IEEE Electron Device Lett.
41
,
296
(
2020
).
31.
A.
Nikolskaya
,
E.
Okulich
,
D.
Korolev
,
A.
Stepanov
,
D.
Nikolichev
,
A.
Mikhaylov
,
D.
Tetelbaum
,
A.
Almaev
,
C. A.
Bolzan
,
A.
Buaczik
,
R.
Giulian
,
P. L.
Grande
,
A.
Kumar
,
M.
Kumar
, and
D.
Gogova
,
J. Vac. Sci. Technol., A
39
,
030802
(
2021
).
32.
A. F. M. A. U.
Bhuiyan
,
Z.
Feng
,
L.
Meng
,
A.
Fiedler
,
H.-L.
Huang
,
A. T.
Neal
,
E.
Steinbrunner
,
S.
Mou
,
J.
Hwang
,
S.
Rajan
, and
H.
Zhao
,
J. Appl. Phys.
131
,
145301
(
2022
).
33.
K.
Goto
,
K.
Konishi
,
H.
Murakami
,
Y.
Kumagai
,
B.
Monemar
,
M.
Higashiwaki
,
A.
Kuramata
, and
S.
Yamakoshi
,
Thin Solid Films
666
,
182
(
2018
).
34.
K. D.
Leedy
,
K. D.
Chabak
,
V.
Vasilyev
,
D. C.
Look
,
J. J.
Boeckl
,
J. L.
Brown
,
S. E.
Tetlak
,
A. J.
Green
,
N. A.
Moser
,
A.
Crespo
,
D. B.
Thomson
,
R. C.
Fitch
,
J. P.
McCandless
, and
G. H.
Jessen
,
Appl. Phys. Lett.
111
,
012103
(
2017
).
35.
F.
Alema
,
T.
Itoh
,
S.
Vogt
,
J. S.
Speck
, and
A.
Osinsky
,
Jpn. J. Appl. Phys.
61
,
100903
(
2022
).
36.
F.
Alema
,
C.
Peterson
,
A.
Bhattacharyya
,
S.
Roy
,
S.
Krishnamoorthy
, and
A.
Osinsky
,
IEEE Electron Device Lett.
43
,
1649
(
2022
).
37.
N. K.
Kalarickal
,
Z.
Xia
,
J.
McGlone
,
S.
Krishnamoorthy
,
W.
Moore
,
M.
Brenner
,
A. R.
Arehart
,
S. A.
Ringel
, and
S.
Rajan
,
Appl. Phys. Lett.
115
,
152106
(
2019
).
38.
J. P.
McCandless
,
V.
Protasenko
,
B. W.
Morell
,
E.
Steinbrunner
,
A. T.
Neal
,
N.
Tanen
,
Y.
Cho
,
T. J.
Asel
,
S.
Mou
,
P.
Vogt
,
H. G.
Xing
, and
D.
Jena
,
Appl. Phys. Lett.
121
,
072108
(
2022
).
39.
A.
Ardenghi
,
O.
Bierwagen
,
A.
Falkenstein
,
G.
Hoffmann
,
J.
Lähnemann
,
M.
Martin
, and
P.
Mazzolini
,
Appl. Phys. Lett.
121
,
042109
(
2022
).
40.
Y.
Oshima
,
E.
Ahmadi
,
S.
Kaun
,
F.
Wu
, and
J. S.
Speck
,
Semicond. Sci. Technol.
33
,
015013
(
2018
).
41.
P.
Vogt
,
A.
Mauze
,
F.
Wu
,
B.
Bonef
, and
J. S.
Speck
,
Appl. Phys. Express
11
,
115503
(
2018
).
42.
A.
Mauze
,
Y.
Zhang
,
T.
Itoh
,
F.
Wu
, and
J. S.
Speck
,
APL Mater.
8
,
021104
(
2020
).
43.
P.
Mazzolini
,
A.
Falkenstein
,
C.
Wouters
,
R.
Schewski
,
T.
Markurt
,
Z.
Galazka
,
M.
Martin
,
M.
Albrecht
, and
O.
Bierwagen
,
APL Mater.
8
,
011107
(
2020
).
44.
P.
Vogt
,
O.
Brandt
,
H.
Riechert
,
J.
Lähnemann
, and
O.
Bierwagen
,
Phys. Rev. Lett.
119
,
196001
(
2017
).
45.
P.
Vogt
and
O.
Bierwagen
,
Appl. Phys. Lett.
108
,
072101
(
2016
).
46.
T. J.
Asel
,
E.
Steinbrunner
,
J.
Hendricks
, and
A. T.
Neal
,
J. Vac. Sci. Technol., A
2
,
043403
(
2020
).
47.
D.
Starodub
,
E. P.
Gusev
,
E.
Garfunkel
, and
T.
Gustafsson
,
Surf. Rev. Lett.
06
,
45
(
1999
).
48.
C.
Gelain
,
A.
Cassuto
, and
P.
Le Goff
,
Oxid. Met.
3
,
139
(
1971
).
49.
F.
Afroz Faria
,
J.
Guo
,
P.
Zhao
,
G.
Li
,
P.
Kumar Kandaswamy
,
M.
Wistey
,
H.
(Grace) Xing
, and
D.
Jena
,
Appl. Phys. Lett.
101
,
032109
(
2012
).
50.
Y.
Oshima
,
E.
Ahmadi
,
S. C.
Badescu
,
F.
Wu
, and
J. S.
Speck
,
Appl. Phys. Express
9
,
061102
(
2016
).
51.
J.
Neugebauer
,
Phys. Status Solidi C
0
,
1651
(
2003
).
52.
A. L.
Rosa
,
J.
Neugebauer
,
J. E.
Northrup
,
C.-D.
Lee
, and
R. M.
Feenstra
,
Appl. Phys. Lett.
80
,
2008
(
2002
).
53.
K.
Hoshikawa
,
T.
Kobayashi
,
E.
Ohba
, and
T.
Kobayashi
,
J. Cryst. Growth
546
,
125778
(
2020
).