In this work, we report on the growth of high-mobility β-Ga2O3 homoepitaxial thin films grown at a temperature much lower than the conventional growth temperature window for metalorganic vapor phase epitaxy. Low-temperature β-Ga2O3 thin films grown at 600 °C on Fe-doped (010) bulk substrates exhibit remarkable crystalline quality, which is evident from the measured room temperature Hall mobility of 186 cm2/V s for the unintentionally doped films. N-type doping is achieved by using Si as a dopant, and a controllable doping in the range of 2 × 1016–2 × 1019 cm−3 is studied. Si incorporation and activation is studied by comparing the silicon concentration from secondary ion mass spectroscopy and the electron concentration from temperature-dependent Hall measurements. The films exhibit high purity (low C and H concentrations) with a very low concentration of compensating acceptors (2 × 1015 cm−3) even at this growth temperature. Additionally, an abrupt doping profile with a forward decay of ∼ 5 nm/dec (10 times improvement compared to what is observed for thin films grown at 810 °C) is demonstrated by growing at a lower temperature.

β-Ga2O3 is the emerging material of choice for power device applications owing to its material properties that translate to the predicted Baliga figure of merit (BFOM) several times larger than that of SiC and GaN.1 To reach this intrinsic potential, a material with low background doping and high mobility is necessary. Carrier transport in β-Ga2O3 suffers from severe polar optical phonon scattering that limits its theoretical room temperature mobility to 200 cm2/V s.2 Metalorganic vapor phase epitaxy (MOVPE) has emerged as the leading growth technique that allows for the growth of epitaxial films with mobility values closer to this theoretical limit.3 So far, the homoepitaxy of β-Ga2O3 using MOVPE has been demonstrated with different reactor geometries, different precursors and dopants, and different substrate orientations. All the growth and doping studies done on these homoepitaxial β-Ga2O3 thin films are typically performed at high temperatures (>800 °C).4–12 The lower temperature regime < 800 °C still remains unexplored and could help in furthering the understanding of the growth window for β-Ga2O3 epitaxy.

The low temperature epitaxial growth (LT) of β-Ga2O3 is interesting to study, especially for the development of power and high-frequency (RF) devices. Dopant incorporation, desorption, and segregation in semiconductor growth are strongly correlated with the growth temperature and could be better controlled if films are grown at lower temperatures. The effect of these phenomena could be very important in the development of 2D channel RF devices such as delta-doped and modulation-doped FETs where extremely sharp doping profiles become critical. Another key application of MOVPE LT epitaxy could be the Ohmic contact regrowth process to buried channel structures that could expand the capability of the MOVPE technique in terms of device processing.13,14 The regrowth process if performed at high temperatures could lead to dopant diffusion, creation, and migration of defects in the active region that could result in channels with reduced mobility. Here, we explore the low temperature regime of the MOVPE growth window by tuning the growth temperature down to 600 °C. We show that β-Ga2O3 thin films with electronic properties similar to thin films grown at higher temperatures are achievable. This allows us to expand the growth window of MOVPE by almost more than 200 °C. Furthermore, we demonstrate a sharp silicon doping profile enabled by low temperature epitaxy.

β-Ga2O3 thin films were grown on Fe-doped (010) semi-insulating bulk substrates (Novel Crystal Technology, Japan) in an MOVPE reactor (Agnitron Agilis) using triethylgallium (TEGa), O2, and diluted silane as precursors and argon as the carrier gas. Growths were carried out at substrate temperatures (Tg) of 550–810 °C. Prior to loading the substrates into the growth chamber, the substrates were cleaned with acetone, methanol, and DI water followed by a 30 min dip in a diluted HF solution.3 Films with thicknesses of 400–1000 nm were grown. The uniformity of the doping profile was characterized by capacitance–voltage (CV) measurements on lateral Schottky diode structures with Ni/Au Schottky pads and Ti/Au Ohmic contacts.

The films were grown at growth temperatures of Tg = 550 °C, 600 °C, 700 °C, and 810 °C. The gas flow values were set to 800  sccm and 1100  sccm for O2 and Ar, respectively, and the chamber pressure was maintained at 60 Torr. The TEGa bubbler was maintained at 16 °C, and the flow rate was set to 15.74 μmol/min. By keeping all other growth conditions identical, Tg was the only growth parameter varied. Figure 1 shows the atomic force microscopy (AFM) scans of the surface morphology of representative films grown at (a) 700 °C, (b) 600 °C, and (c) 550 °C. The films grown at Tg = 600 °C and 700 °C show extremely smooth surfaces with RMS roughness values of 0.38 nm and 0.53 nm, respectively. The film grown at Tg = 550 °C shows a rougher surface morphology (RMS roughness = 2.2 nm), indicating the onset of enhanced nucleation because of the reduced Ga adatom mobility.

FIG. 1.

5 μm × 5μm AFM (Bruker dimension icon) scans of the β-Ga2O3 thin films grown at substrate temperatures of (a) 700 °C (film thickness = 430 nm), (b) 600 °C (film thickness = 380 nm), and (c) 550 °C (film thickness = 220 nm). (d) Growth rate as a function of the growth temperature.

FIG. 1.

5 μm × 5μm AFM (Bruker dimension icon) scans of the β-Ga2O3 thin films grown at substrate temperatures of (a) 700 °C (film thickness = 430 nm), (b) 600 °C (film thickness = 380 nm), and (c) 550 °C (film thickness = 220 nm). (d) Growth rate as a function of the growth temperature.

Close modal

The step-flow growth condition could be achieved because the diffusion coefficient of Ga adatoms is much higher on a Ga2O3 surface as compared to that on a GaN surface which, along with ammonia cracking constraints, limits GaN growth to higher temperatures. Early modeling and experimental work reveal that Ga adatom surface diffusivity on a Ga2O3 substrate could be six orders higher than in GaN.8 The growth of homoepitaxial Ga2O3 films at lower temperatures is further supported by the fact that the cracking temperature of the precursor (TEGa) is lower than the growth temperatures used in this work.15 

Figure 1(d) shows the growth rate of the films with different Tg. The growth rate was measured and corroborated using secondary-ion mass spectroscopy (SIMS) and cross-sectional scanning electron microscopy (SEM) imaging of co-loaded c-plane sapphire substrates. The films grown at Tg = 700 °C and 810 °C showed a similar growth rate of ∼7 nm/min, indicating that, in this regime, the growth is mass transport limited. The films grown at Tg = 550 °C and 600 °C show reduced growth rates of ∼3 nm/min and ∼6 nm/min, respectively. This indicates that the growth rate is now surface kinetics limited. Therefore, further low temperature growth experiments were performed on films grown at 600 °C.

N-type doping is achieved in these films by flowing diluted silane into the chamber during the growth, and doping in the range of 2 × 1016–2 × 1019 cm−3 was studied. To understand the incorporation of Si and other unintentional impurities (C, H) in the films grown at 600 °C, secondary-ion mass spectroscopy (SIMS) measurements were performed on a stacked structure with different silane flows (see the supplementary material). Figure 2 shows the measured Si atom density (left Y axis) as a function of the silane flow. The Si atom density scales linearly with the silane flow, suggesting complete incorporation. Carbon (C) is always found to be near the SIMS detection limits (5 × 1016 cm−3) throughout the growth temperature range (supplementary material, Fig. SP2). We hypothesize that the low C incorporation in the films is likely due to the growth under an oxidizing atmosphere, wherein CO2 can easily desorb from the surface during the growth. At lower silane flows, hydrogen (H) is found to be near the SIMS detection limits (2 × 1016 cm−3) for films grown at 600 °C (supplementary material, Fig. SP2). However, the H concentration is orders of magnitude lower than the Si concentration even at higher silane flows and the active carrier concentration corresponds well with the Si concentration. Thus, high purity β-Ga2O3 thin films with controllable doping are achieved and further verified by the excellent transport properties as discussed later.

FIG. 2.

Comparison of Si density from SIMS and carrier density extracted from Hall measurements in the 600 °C grown β-Ga2O3 thin films as a function of silane flow.

FIG. 2.

Comparison of Si density from SIMS and carrier density extracted from Hall measurements in the 600 °C grown β-Ga2O3 thin films as a function of silane flow.

Close modal

Hall effect measurements (Ecopia HMS 7000) were performed to study the transport properties of these films. Ohmic contacts were formed by depositing Ti (50 nm)/Au (100 nm) on the four corners of the as-grown unintentionally doped (UID) and Si-doped β-Ga2O3 thin films using DC sputtering. The contacts were further annealed in N2 at 450 ∘C for 1.5 minutes so that the contacts retained their Ohmic nature even at low temperatures. Figure 2 (right Y axis) shows the measured carrier densities as a function of the silane flow for the films grown at 600 °C. The carrier density matches closely with the measured Si density from SIMS, which indicates complete activation of the Si atoms. The doping profiles were further extracted by CV profiling (see the supplementary material) and are found to be very uniform throughout the film thicknesses. Room temperature (RT) Hall mobilities as a function of carrier density are shown in Fig. 3. The UID film grown at 600 °C exhibits an impressive room temperature mobility of 186 cm2/V s (for ND+ − NA ∼ 2 × 1016 cm−3), which is close to the highest reported mobility values of high-temperature MOVPE-grown UID films with similar background doping. Our UID films grown at 810 °C exhibit a room-temperature Hall mobility of 184 cm2/V s for similar background doping (2 × 1016 cm−3). This clearly indicates that our LT-grown film is of high crystalline quality (low defects/compensation) equivalent to those grown at high temperatures.

FIG. 3.

Room temperature Hall mobility vs carrier density in the 600 °C grown β-Ga2O3 thin films.

FIG. 3.

Room temperature Hall mobility vs carrier density in the 600 °C grown β-Ga2O3 thin films.

Close modal

To further quantify the compensation and the Si donor behavior for films grown at 600 °C, temperature-dependent Hall measurements (80 K–340 K) were performed and were compared to films grown at 810 °C. Lightly doped films were grown at 600 °C and 810 °C by keeping the silane flows and all other growth conditions the same. Figure 4(a) shows the temperature-dependent mobility of the two films. Both films exhibit high RT mobilities of 140 cm2/V s (810 °C) and 133 cm2/V s (600 °C) for carrier densities of 1.75 × 1017 cm−3 and 1.2 × 1017 cm−3, respectively. The lower mobility for the LT grown film could be due to higher Si contaminants at the epilayer/substrate interface, which forms a low mobility parasitic channel.3 For low-doped films, the contribution from this parasitic channel could be non-trivial.3 The effect of this is also evident from the lower LT peak mobility of the LT grown film.

FIG. 4.

(a) Temperature-dependent Hall mobility of a lightly doped β-Ga2O3 thin film grown at 600 °C and compared with a film grown at 810 °C having nominally the same growth conditions (the inset shows the mobility values on the logarithmic scale). (b) Hall carrier density vs 1000/T of these two films. (Colored diamonds represent the measured values and the black dashed lines represent the fitting to the 2 donor model.)

FIG. 4.

(a) Temperature-dependent Hall mobility of a lightly doped β-Ga2O3 thin film grown at 600 °C and compared with a film grown at 810 °C having nominally the same growth conditions (the inset shows the mobility values on the logarithmic scale). (b) Hall carrier density vs 1000/T of these two films. (Colored diamonds represent the measured values and the black dashed lines represent the fitting to the 2 donor model.)

Close modal

A two donor and one acceptor charge model16,17 was used to analyze the temperature-dependent carrier densities in the two films as shown in Fig. 4(b), and the extracted parameters are summarized in Table I. Both films show a shallow donor state with an activation energy of 17–19 meV that could be attributed to Si. Two different deep donor states are also observed in these films. These donor energy levels have also been observed in films grown by different growth techniques, but their origin still remains to be understood. Deep donor states in β-Ga2O3 could originate from extrinsic traps/defects such as interstitials, antisites, and the incorporation of unintentional impurities.17 Also, the presence of residual Si impurities at the epilayer/substrate interface complicates the analysis of charge and transport characteristics. To date, no standard processing technique has been developed that could remove the accumulated residual Si impurities, which makes it an uncontrolled parameter. Nevertheless, the extracted concentrations of compensating acceptors in both the films are very low, of the order ∼2 × 1015 cm−3. This shows that β-Ga2O3 thin films with high crystalline quality (low defects) can now be grown over a wider growth temperature window of >200 °C using the MOVPE growth technique. The low temperature epitaxy approach using MOVPE now allows for the selective area regrowth of Ohmic contacts on β-Ga2O3 epitaxial films, which was so far limited to MBE, PLD, and ion-implantation techniques.13,14,18

TABLE I.

Summary of free parameters extracted by fitting the two donor one acceptor charge model to the experimental values. Nd1 and Nd2 are the concentrations of the two donors, Na the concentration of compensating acceptors, and Ed1 and Ed2 the donor energies. An effective mass of m* = 0.28mo was used to estimate Nc analytically.

TgNd1Ed1Nd2Ed2Naμ (RT)
(°C)(cm−3)(meV)(cm−3)(meV)(cm−3)(cm2/Vs)
600 1.1 × 1017 19 5.5 × 1016 100 2.5 × 1015 133 
810 1.3 × 1017 17 7.7 × 1016 65 2.0 × 1015 140 
TgNd1Ed1Nd2Ed2Naμ (RT)
(°C)(cm−3)(meV)(cm−3)(meV)(cm−3)(cm2/Vs)
600 1.1 × 1017 19 5.5 × 1016 100 2.5 × 1015 133 
810 1.3 × 1017 17 7.7 × 1016 65 2.0 × 1015 140 

Next, we studied the effect of the growth temperature on abruptness of the doping profiles in these Si-doped β-Ga2O3 thin films. Stacks with intentionally-doped layers sandwiched between UID layers (see the supplementary material) were grown at different temperatures, and only the silane flow was switched without any growth interruption between the layers. Both doping and depth profiles were analyzed by SIMS measurements. Figure 5 shows the doping vs depth profiles of the representative structures with similar silane flows grown at three different temperatures of Tg = 810 °C, 700 °C, and 600 °C. It can be seen that the forward Si decay (toward the growth direction) is strongly dependent on the growth temperature and becomes steeper when the growth temperature is lowered. By reducing the growth temperature from 810 °C to 600 °C, the forward Si decay changed from ∼50 nm/dec to ∼5 nm/dec. Note that the decay of the forward Si tail improved monotonically as the growth temperature was lowered. The corresponding forward Si decays for three different growth temperatures are found to be consistent over a wide range of Si densities in the doped layer—3 × 1017 cm−3 to 2 × 1019 cm−3 (Fig. 6).

FIG. 5.

Profile of Si, C, and H in the films grown at (a) 810 °C, (b) 700 °C, and (c) 600 °C analyzed by SIMS measurements. Black dashed lines indicate the silane switching times and the expected ideal doping profile.

FIG. 5.

Profile of Si, C, and H in the films grown at (a) 810 °C, (b) 700 °C, and (c) 600 °C analyzed by SIMS measurements. Black dashed lines indicate the silane switching times and the expected ideal doping profile.

Close modal
FIG. 6.

Forward Si decay rates for various doping densities in films grown at different temperatures.

FIG. 6.

Forward Si decay rates for various doping densities in films grown at different temperatures.

Close modal

It was observed that the reverse Si decay (toward the substrate) was nearly independent of the growth temperature as well as the Si density in the doped layer. This suggests that the forward Si tail that extends into the undoped layers could not be explained due to bulk diffusion of Si alone19,20 because Si in β-Ga2O3 does not diffuse (bulk diffusion) until heated to very high temperatures of 1000–1100 °C.18,21 A possible explanation could be that at higher growth temperatures, a significant amount of Si segregates and accumulates at the growth surface with lower incorporation in the preceding monolayers and continues to “ride” the surface even after the silane flow is switched off.19,22 This understanding is further strengthened by the observation of a significant increment of spatial delay in reaching the steady-state doping with higher growth temperatures after the silane flow is switched on as shown in Fig. 5. The characteristic delay length (ΔR) is found to decrease as the growth temperature is reduced. This kind of doping profile smearing has been observed in other semiconductor growths as well and has been attributed to variation in dopant segregation, desorption, and incorporation rates with growth temperature.23–29 In our case, the effect of growth temperature on Si segregation at the growth surface could be more pronounced than Si desorption. In general, the incorporation rate during semiconductor growth depends not only on the dopant desorption rates, dopant surface segregation (surface coverages), and the incoming dopant flux but also on the growth rate. Some additional strategies such as growth interruption coupled with purge cycles could be useful to achieve even sharper doping profiles. To quantitatively explain the doping profile smearing observed at elevated growth temperatures, modeling and experimental verification of dopant desorption, incorporation, segregation, and bulk diffusion for MOVPE β-Ga2O3 growth are necessary.

In conclusion, we demonstrate the growth of UID and Si-doped β-Ga2O3 thin films on Fe-doped (010) bulk substrates using the MOVPE technique at a growth temperature of 600 °C, which is significantly lower than the conventional MOVPE growth temperatures reported. High purity films with excellent transport properties are achieved, evident from the high RT Hall mobility of 186 cm2/V s measured on UID films as well as high mobilities in the Si-doped film. Extremely smooth films with atomically flat surfaces are achieved at this growth temperature that was attributed to higher Ga adatom mobility on a Ga2O3 surface. Si density from SIMS and carrier density from Hall measurements match closely and scale linearly with the silane flow, indicating complete incorporation and activation of Si. Compensating acceptors and other impurities such as C and H are found to be low and comparable to films grown at higher temperatures. Abrupt doping profiles are achieved with forward Si decay as sharp as 5 nm/dec. This was achieved due to the suppression of Si “surface riding” by growing at a lower temperature. This demonstration widens the MOVPE growth temperature window for β-Ga2O3 homoepitaxy by more than 200 °C. High-quality β-Ga2O3 thin films with sharp doping profiles could be beneficial for the development of β-Ga2O3-based RF devices.

See the supplementary material for CV-extracted doping profiles and profiles of Si, C, and H from SIMS scans on 600 °C and 810 °C grown stacks.

This work was supported by the Air Force Office of Scientific Research under Award No. FA9550-18-1-0507 (Program Manager: Dr. Ali Sayir). We also acknowledge the II–VI foundation Block Gift Program for financial support. This work was performed in part at the Utah Nanofab sponsored by the College of Engineering and the Office of the Vice President for Research.

The data that support the findings of this study are available within the article and its supplementary material.

1.
M.
Higashiwaki
and
G. H.
Jessen
, “
Guest editorial: The dawn of gallium oxide microelectronics
,”
Appl. Phys. Lett.
112
,
060401
(
2018
).
2.
N.
Ma
,
N.
Tanen
,
A.
Verma
,
Z.
Guo
,
T.
Luo
,
H. G.
Xing
, and
D.
Jena
, “
Intrinsic electron mobility limits in β-Ga2O3
,”
Appl. Phys. Lett.
109
,
212101
(
2016
).
3.
Z.
Feng
,
A. F. M. A. U.
Bhuiyan
,
Z.
Xia
,
W.
Moore
,
Z.
Chen
,
J. F.
McGlone
,
D. R.
Daughton
,
A. R.
Arehart
,
S. A.
Ringel
,
S.
Rajan
, and
H.
Zhao
, “
Probing charge transport and background doping in MOCVD grown (010) β-Ga2O3
,”
Phys. Status. Solidi. RRL
14
,
2000145
(
2020
).
4.
Z.
Feng
,
A. F. M.
Anhar Uddin Bhuiyan
,
M. R.
Karim
, and
H.
Zhao
, “
MOCVD homoepitaxy of Si-doped (010) β-Ga2O3 thin films with superior transport properties
,”
Appl. Phys. Lett.
114
,
250601
(
2019
).
5.
F.
Alema
,
B.
Hertog
,
A.
Osinsky
,
P.
Mukhopadhyay
,
M.
Toporkov
, and
W. V.
Schoenfeld
, “
Fast growth rate of epitaxial β-Ga2O3 by close coupled showerhead MOCVD
,”
J. Cryst. Growth
475
,
77
82
(
2017
).
6.
Y.
Zhang
,
F.
Alema
,
A.
Mauze
,
O. S.
Koksaldi
,
R.
Miller
,
A.
Osinsky
, and
J. S.
Speck
, “
MOCVD grown epitaxial β-Ga2O3 thin film with an electron mobility of 176 cm2/V s at room temperature
,”
APL Mater.
7
,
022506
(
2019
).
7.
S.
Bin Anooz
,
R.
Grüneberg
,
C.
Wouters
,
R.
Schewski
,
M.
Albrecht
,
A.
Fiedler
,
K.
Irmscher
,
Z.
Galazka
,
W.
Miller
,
G.
Wagner
,
J.
Schwarzkopf
, and
A.
Popp
, “
Step flow growth of β-Ga2O3 thin films on vicinal (100) β-Ga2O3 substrates grown by MOVPE
,”
Appl. Phys. Lett.
116
,
182106
(
2020
).
8.
R.
Schewski
,
M.
Baldini
,
K.
Irmscher
,
A.
Fiedler
,
T.
Markurt
,
B.
Neuschulz
,
T.
Remmele
,
T.
Schulz
,
G.
Wagner
,
Z.
Galazka
, and
M.
Albrecht
, “
Evolution of planar defects during homoepitaxial growth of β-Ga2O3 layers on (100) substrates-A quantitative model
,”
J. Appl. Phys.
120
,
225308
(
2016
).
9.
M.
Baldini
,
M.
Albrecht
,
A.
Fiedler
,
K.
Irmscher
,
D.
Klimm
,
R.
Schewski
, and
G.
Wagner
, “
Semiconducting Sn-doped β-Ga2O3 homoepitaxial layers grown by metal organic vapour-phase epitaxy
,”
J. Mater. Sci.
51
,
3650
3656
(
2016
).
10.
P.
Ranga
,
A.
Bhattacharyya
,
A.
Rishinaramangalam
,
Y. K.
Ooi
,
M. A.
Scarpulla
,
D.
Feezell
, and
S.
Krishnamoorthy
, “
Delta-doped β-Ga2O3 thin films and β-(Al0.26Ga0.74)2O3/β-Ga2O3 heterostructures grown by metalorganic vapor phase epitaxy
,”
Appl. Phys. Express
13
,
045501
(
2020
).
11.
G.
Wagner
,
M.
Baldini
,
D.
Gogova
,
M.
Schmidbauer
,
R.
Schewski
,
M.
Albrecht
,
Z.
Galazka
,
D.
Klimm
, and
R.
Fornari
, “
Homoepitaxial growth of β-Ga2O3 layers by metal-organic vapor phase epitaxy
,”
Phys. Status Solidi A
211
,
27
33
(
2014
).
12.
P.
Ranga
,
A.
Rishinaramangalam
,
J.
Varley
,
A.
Bhattacharyya
,
D.
Feezell
, and
S.
Krishnamoorthy
, “
Si-doped β-(Al0.26Ga0.74)2O3 thin films and heterostructures grown by metalorganic vapor-phase epitaxy
,”
Appl. Phys. Express
12
,
111004
(
2019
).
13.
K. D.
Chabak
,
K. D.
Leedy
,
A. J.
Green
,
S.
Mou
,
A. T.
Neal
,
T.
Asel
,
E. R.
Heller
,
N. S.
Hendricks
,
K.
Liddy
,
A.
Crespo
 et al, “
Lateral β-Ga2O3 field effect transistors
,”
Semicond. Sci. Technol.
35
,
013002
(
2020
).
14.
Z.
Xia
,
C.
Joishi
,
S.
Krishnamoorthy
,
S.
Bajaj
,
Y.
Zhang
,
M.
Brenner
,
S.
Lodha
, and
S.
Rajan
, “
Delta doped β-Ga2O3 field effect transistors with regrown Ohmic contacts
,”
IEEE Electron Device Lett.
39
,
568
571
(
2018
).
15.
R.
Miller
,
F.
Alema
, and
A.
Osinsky
, “
Epitaxial β-Ga2O3 and β-(AlxGa1−x)2O3/β-Ga2O3 heterostructures growth for power electronics
,”
IEEE Trans. Semicond. Manuf.
31
,
467
474
(
2018
).
16.
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
 et al, “
Donors and deep acceptors in β-Ga2O3
,”
Appl. Phys. Lett.
113
,
062101
(
2018
).
17.
A. T.
Neal
,
S.
Mou
,
R.
Lopez
,
J. V.
Li
,
D. B.
Thomson
,
K. D.
Chabak
, and
G. H.
Jessen
, “
Incomplete ionization of a 110 meV unintentional donor in β-Ga2O3 and its effect on power devices
,”
Sci. Rep.
7
,
13218
(
2017
).
18.
K.
Sasaki
,
M.
Higashiwaki
,
A.
Kuramata
,
T.
Masui
, and
S.
Yamakoshi
, “
Si-Ion implantation doping in β-Ga2O3 and its application to fabrication of low-resistance Ohmic contacts
,”
Appl. Phys. Express
6
,
086502
(
2013
).
19.
A.
Mauze
,
Y.
Zhang
,
T.
Mates
,
F.
Wu
, and
J. S.
Speck
, “
Investigation of unintentional Fe incorporation in (010) β-Ga2O3 films grown by plasma-assisted molecular beam epitaxy
,”
Appl. Phys. Lett.
115
,
052102
(
2019
).
20.
S.
Iyer
,
R.
Metzger
, and
F.
Allen
, “
Sharp profiles with high and low doping levels in silicon grown by molecular beam epitaxy
,”
J. Appl. Phys.
52
,
5608
5613
(
1981
).
21.
R.
Sharma
,
E.
Mark
,
C.
Fares
,
M.
Tadjer
,
F.
Ren
,
A.
Kuramata
, and
S. J.
Pearton
, “
The role of annealing ambient on diffusion of implanted Si in β-Ga2O3
,”
AIP Adv.
9
,
085111
(
2019
).
22.
P.
Ranga
,
A.
Bhattacharyya
,
A.
Chmielewski
,
S.
Roy
,
N.
Alem
, and
S.
Krishnamoorthy
, “
Delta-doped
β-Ga2O3 films with low FWHM charge profile grown by metalorganic vapor-phase epitaxy
,” arXiv:2008.12934 (
2020
).
23.
J.-E.
Sundgren
,
J.
Knall
,
W.-X.
Ni
,
M.-A.
Hasan
,
L.
Markert
, and
J.
Greene
, “
Dopant incorporation kinetics and abrupt profiles during silicon molecular beam epitaxy
,”
Thin Solid Films
183
,
281
297
(
1989
).
24.
J.
Greene
,
S.
Barnett
,
A.
Rockett
, and
G.
Bajor
, “
Modeling of dopant incorporation, segregation, and ion/surface interaction effects during semiconductor film growth by molecular beam epitaxy and plasma-based techniques
,”
Appl. Surf. Sci.
22–23
,
520
544
(
1985
).
25.
A.
Rockett
,
J.
Klem
,
S.
Barnett
,
J.
Greene
, and
H.
Morkoc
, “
Si incorporation probabilities and depth distributions in Ga1xAlxAs films grown by molecular-beam epitaxy
,”
J. Appl. Phys.
59
,
2777
2783
(
1986
).
26.
S.
Barnett
and
J.
Greene
, “
Si molecular beam epitaxy: A model for temperature dependent incorporation probabilities and depth distributions of dopants exhibiting strong surface segregation
,”
Surf. Sci.
151
,
67
90
(
1985
).
27.
K.
Köhler
,
R.
Gutt
,
J.
Wiegert
, and
L.
Kirste
, “
Diffusion of Mg dopant in metal-organic vapor-phase epitaxy grown GaN and AlxGa1−xN
,”
J. Appl. Phys.
113
,
073514
(
2013
).
28.
K.
Tomita
,
K.
Itoh
,
O.
Ishiguro
,
T.
Kachi
, and
N.
Sawaki
, “
Reduction of Mg segregation in a metalorganic vapor phase epitaxial grown GaN layer by a low-temperature AlN interlayer
,”
J. Appl. Phys.
104
,
014906
(
2008
).
29.
O.
Elleuch
,
K.
Lekhal
,
Y.
Guan
, and
T. F.
Kuech
, “
Highly tin doped GaAs at low growth temperatures using tetraethyl tin by metal organic vapor phase epitaxy
,”
J. Cryst. Growth
507
,
255
259
(
2019
).

Supplementary Material