Gamma/Beta-Ga2O3 Polymorph Heterostructures
- Gamma/Beta-Ga2O3 polymorph heterostructures are homo-interfaces that join metastable cubic and stable monoclinic Ga2O3 phases without changing stoichiometry.
- Formation methods such as ion irradiation and MOCVD allow precise control over layer thickness, interface abruptness, and defect profiles to tailor functional properties.
- Atomic-scale and thermal studies reveal strain-sensitive band offsets with type-II alignment, offering design principles for high-power electronic and optoelectronic devices.
Searching arXiv for papers on γ/β-Ga2O3 polymorph heterostructures to ground the article in the cited literature. Gamma/Beta-GaO polymorph heterostructures are homo-interfaces formed between two polymorphs of the same compound, namely metastable cubic defective-spinel -GaO and stable monoclinic -GaO. Unlike conventional semiconductor heterostructures, where interfacial functionality is produced by changing chemical composition, these structures modulate properties across an interface while retaining identical chemistry on both sides. In the recent literature, such stacks are generated either by disorder-induced ordering in -GaO0 under ion irradiation or, in a distinct epitaxial context, by exploiting a narrow MOCVD window for 1-Ga2O3 on MgAl4O5, with 6 inclusions or conversion emerging under thermal treatment (Abdullaev et al., 2024, Azarov et al., 2023, Tang et al., 2023).
1. Definition, structural motif, and distinction from conventional heterostructures
The defining feature of a 7-Ga8O9 polymorph heterostructure is that it joins two crystal structures of Ga0O1 across a sharply defined interface without changing stoichiometry. In the ion-irradiation-derived form, a top layer of 2-Ga3O4 transforms into 5-Ga6O7 once a disorder threshold is reached, leaving a crystalline 8 film atop the remaining 9 substrate and thereby creating a homo-interface between polymorphs of the same oxide (Abdullaev et al., 2024, Azarov et al., 2023). This configuration differs from conventional heterostructures because the interfacial contrast arises from symmetry, strain, and defect topology rather than from a compositional discontinuity.
The two constituent phases are structurally distinct. 0-Ga1O2 is monoclinic and anisotropic, whereas 3-Ga4O5 is a cubic defective spinel with a well-ordered oxygen FCC or cubic close-packed sublattice and a partially occupied Ga sublattice (Tang et al., 2023, Abdullaev et al., 9 Sep 2025). In diffraction, the 6 phase is associated with a broad peak near 7 assigned to the 8 reflection of the cubic spinel structure in irradiation-derived stacks, and with an out-of-plane 9 peak at 0 in MOCVD-grown films on MgAl1O2(100) (Abdullaev et al., 2024, Tang et al., 2023).
A central point in the literature is interface abruptness. Rutherford backscattering spectroscopy in channeling mode shows a “box-like” signal characteristic of the disorder-induced 3 transition, while prior TEM work and modeling support a sharply defined 4 interface; in thermal-transport modeling, an interface thickness of 5 nm was adopted conservatively (Abdullaev et al., 2024). High-resolution transmission electron microscopy likewise confirms abruptness and crystallographic registry in the band-offset study, where the interface is treated as atomically abrupt and formed with preservation of the oxygen sublattice and Ga migration (Liu et al., 2024).
This structural paradigm has broader implications. Because the interface is chemically homogeneous yet crystallographically discontinuous, it provides a platform for tuning thermal transport, optical signatures, radiation response, and band alignment through orientation, strain state, thickness, and disorder history rather than by alloy design. A plausible implication is that 6-Ga7O8 occupies an intermediate conceptual position between classical heteroepitaxy and phase-boundary engineering in functional oxides.
2. Formation routes and crystallographic realization
The most established route is disorder-induced ordering in 9-Ga0O1. In this approach, ion irradiation at room temperature converts the near-surface region of 2-Ga3O4 into 5-Ga6O7 once sufficient disorder has accumulated (Azarov et al., 2023). In thermal-transport studies, Ga8 irradiation at a fluence of 9 cm0 and energies of 0.5, 1.0, and 1.7 MeV yielded 1-layer thicknesses of approximately 350, 650, and 1000 nm, respectively, on 2, 3, and 4 5-Ga6O7 substrates (Abdullaev et al., 2024). In optical studies, 8 layers of approximately 9m were generated from 0-Ga1O2 films on sapphire using 1.5 MeV 3Ni4 at a dose of 5 cm6, and 300 nm auxiliary 7 layers were formed on 8 and 9 0-Ga1O2 wafers to produce double-polymorph bilayers (Galeckas et al., 2024). A related annealing study used 400 keV Ni3 at 4 cm5, with 7° off-axis incidence, to fabricate a 6 nm 7 film on a 8 9-Ga0O1 substrate (Azarov et al., 2024).
These irradiation-generated structures exhibit orientation relationships and low-mismatch registry. High-resolution ADF-STEM in the radiation-tolerance work showed 2 alignment across a sharp interface (Azarov et al., 2023). In the band-offset study, experimentally observed or modeled interface orientations included 3, 4, and 5, with detailed modeling focused on 6 and 7 (Liu et al., 2024).
A distinct realization is epitaxial 8-Ga9O00 growth on MgAl01O02(100) by MOCVD. There, nominally phase-pure 03-Ga04O05 is obtained only within a narrow temperature window centered near 06C; at 07C mixed 08 occurs, and at 09C nominally phase-pure 10-Ga11O12 is obtained (Tang et al., 2023). The epitaxy is cube-on-cube, with 13 and 14 out of plane. In that system, 15-Ga16O17 inclusions can form inside the 18 matrix, and annealing at 19C promotes 20 recrystallization while preserving the ordered oxygen sublattice of the 21 matrix (Tang et al., 2023).
The two formation routes are complementary rather than interchangeable. Ion irradiation directly generates 22 homo-interfaces in bulk or wafer substrates, whereas MOCVD on MgAl23O24 exposes atomic-scale defect physics, epitaxial relationships, and thermal conversion pathways in an epitaxial 25-dominated environment. This suggests that the field uses “heterostructure” in two closely related senses: a deliberately formed 26-on-27 stack and a 28 matrix containing oriented 29 domains or inclusions.
3. Atomic-scale structure, defects, and interface topology
At the atomic scale, 30-Ga31O32 is treated as a defective spinel with a fully ordered cubic close-packed oxygen sublattice and a partially occupied cation sublattice (Tang et al., 2023). In STEM along 33, image contrast scales with Ga column density, and the ideal 34 motif appears as hexagonal arrays with six tetrahedrally coordinated Ga columns at the edges and five octahedral Ga columns, including a central octahedral column with approximately double the Ga column density and therefore higher HAADF intensity (Tang et al., 2023). This ideal motif is especially visible in APB-free regions near the interface.
A high density of antiphase boundaries is a defining feature of relaxed 35-Ga36O37 grown on MgAl38O39(100). These APBs arise from lattice shifts in the Ga sublattice while preserving oxygen close packing. A 40 nm fully strained transition layer at the film/substrate interface is APB-free, after which APBs nucleate as the film relaxes (Tang et al., 2023). The planar-defect model extended from 41-Al42O43 explains these APBs as glide on oxygen-preserving 44 planes with cation sublattice shifts that reproduce the diagonal high-intensity stripes, tri-column low-intensity patterns, and alternating octahedral/tetrahedral intensities seen in HAADF-STEM (Tang et al., 2023).
Within 45, 46-Ga47O48 inclusions exhibit specific epitaxial and rotational-domain relationships. XRD 49-scans after annealing reveal four in-plane rotational domains of 50, with the relationship 51, and 180° rotational 52 domains form via twinning on the 53 plane (Tang et al., 2023). The prevalence of 180° domain pairs is rationalized by the comparatively low twin-boundary formation energy on 54, while 90° 55 rotations are discussed as likely facilitated by the higher symmetry of 56 (Tang et al., 2023).
In the irradiation-derived 57 stacks, the interface is also structurally distinctive. Rutherford backscattering/channeling never reaches the random level, indicating that the transformation does not proceed through full amorphization during 58-layer formation (Abdullaev et al., 2024). In the band-offset work, nanoscale strain maps recorded along 59 and 60 yielded 61 in-plane and 62 out-of-plane, with 63 strained relative to 64 by 65 and 66 at the interface (Liu et al., 2024).
A common misconception is that disorder-induced formation of 67-Ga68O69 implies an amorphous or poorly crystallized layer. The combined diffraction, channeling, and microscopy evidence does not support that interpretation. Rather, the 70 layer is crystalline, exhibits a defective-spinel topology, and can support either abrupt homo-interfaces or coherent inclusion relationships, even though its cation sublattice is intrinsically disordered (Abdullaev et al., 2024, Azarov et al., 2023, Tang et al., 2023).
4. Thermal transport and anisotropic heat-flow engineering
Thermal transport is one of the clearest experimentally quantified functionalities of 71-Ga72O73 polymorph heterostructures. The cross-plane thermal conductivity was measured by femtosecond laser-based time-domain thermoreflectance with modulation frequencies from 0.73 to 10 MHz, using the thermal penetration depth
74
which enabled depth profiling across 75 stacks (Abdullaev et al., 2024). The spot size was approximately 76m, larger than the thermal diffusion length, so the measurements were restricted to cross-plane conductivity.
The room-temperature conductivity contrast between the two phases is strong. For 77-Ga78O79, cross-plane thermal conductivity was measured as 80–81 W m82 K83, independent of the initial 84-substrate orientation (Abdullaev et al., 2024). For pristine 85-Ga86O87, the measured cross-plane thermal conductivity was anisotropic: 88 W m89 K90 for 91, 92 W m93 K94 for 95, and 96 W m97 K98 for 99 (Abdullaev et al., 2024). The 00 values exceed the amorphous Ga01O02 minimum, reported as the Cahill limit of approximately 03 W m04 K05, which confirms the crystalline character of the 06 phase (Abdullaev et al., 2024).
For 07 layers formed by 0.5, 1.0, and 1.7 MeV Ga08 irradiation, three-layer fits yielded 09, 10, and 11 W m12 K13, respectively, with corresponding 14 values of 15, 16, and 17 MW m18 K19 (Abdullaev et al., 2024). In the model, the 20 thermal boundary conductance was fixed at 250 MW m21 K22, consistent with the assumption of an abrupt interface and limited sensitivity to this parameter within the measurement window (Abdullaev et al., 2024).
The most striking thermal result is the order-of-magnitude conductivity step across the interface when 23 is formed on 24-25. There, 26–27 W m28 K29 contrasts with 30–20 W m31 K32, producing a step-like change in cross-interface heat flow (Abdullaev et al., 2024). Depth profiling showed that for the 33 nm 34 layer, the effective conductivity began to increase near 35 MHz, where 36 nm, matching the 37-layer thickness extracted from RBS-C. For thicker 38 layers of 39 and 40 nm, no frequency dependence was observed because the thermal penetration depth at the lowest modulation frequency did not exceed 41 nm, so the measurement probed only the 42 region (Abdullaev et al., 2024).
Molecular-dynamics simulations provide a partially convergent but not identical picture. For 43-Ga44O45, equilibrium MD gave 46 W m47 K48 with the Born–Mayer–Huggins potential and 49 W m50 K51 with tabGAP, while reverse NEMD extrapolation yielded bulk 52 W m53 K54 and a phonon mean free path 55 nm (Abdullaev et al., 2024). The approximately twofold discrepancy relative to experiment was attributed to residual biaxial strain after 56 transformation and/or ion-induced defects; applying strain 57, 58 in MD reduced 59 to values comparable to experiment (Abdullaev et al., 2024).
These results establish a design principle already emphasized in the literature: select 60-61 to maximize the conductivity contrast and tune the 62-layer thickness via irradiation energy to control the depth at which heat encounters the high-63 64 substrate (Abdullaev et al., 2024). This suggests an avenue for thermal energy conversion and phononic functionality within a chemically uniform oxide platform.
5. Optical signatures, phase identification, and annealing-induced transformation
Optical characterization has been used both to build a comparative library of Ga65O66 polymorph signatures and to monitor the evolution of 67 bilayers. In the optical library study, diffuse-reflectance spectroscopy with Kubelka–Munk analysis, transmittance spectroscopy, low-temperature photoluminescence, and nano-FTIR were applied under uniform analysis conditions across polymorphs and heterostructures (Galeckas et al., 2024). For band-gap extraction, the work used
68
with 69 for direct-allowed and 70 for indirect-allowed transitions, and in practice followed the conventional direct-gap analysis across polymorphs (Galeckas et al., 2024). The diffuse-reflectance form was
71
while the Urbach tail was analyzed using
72
For standalone 73-Ga74O75 thin films, the reported room-temperature values were 76 eV and 77 meV from DRS, and 78 eV and 79 meV from transmittance (Galeckas et al., 2024). For 80-Ga81O82 thin film 83, the values were 84 eV and 85 meV by DRS, and 86 eV and 87 meV by transmittance (Galeckas et al., 2024). For the 88 double-polymorph bilayer, DRS yielded 89 eV and 90 meV (Galeckas et al., 2024). The study explicitly noted that the bilayer value is lower than those of standalone 91 and 92 films, reflecting near-surface probing depth, possible spectral mixing across the interface, and the tendency of DRS to underestimate 93 relative to transmittance in thick or stacked structures (Galeckas et al., 2024).
Near-field nano-FTIR provided spatially resolved identification of the heterointerface. A cross-sectional line scan over 1500 nm with 50 nm step size showed distinct near-field amplitude spectra in 94 and 95 domains, enabling phase identification with 96 nm spatial resolution (Galeckas et al., 2024). For 97-Ga98O99, near-field modes were observed around 666, 670, 696, and 728 cm00, and the 01 phase exhibited a distinct near-field spectral envelope sufficient for domain discrimination (Galeckas et al., 2024). This confirms that the 02 interface is optically recognizable on the nanoscale even where detailed atomic-scale interface characterization is not the focus.
Low-temperature photoluminescence is dominated by self-trapped-hole emission. Across crystalline polymorphs, including 03 and 04, intrinsic UV emission in the 2.9–3.3 eV range is attributed to recombination of free electrons with self-trapped holes (Galeckas et al., 2024). Orientation-dependent 05-Ga06O07 spectra show a main UV peak at 3.28 eV with FWHM 0.54 eV for 08, and 3.22 eV with FWHM 09–0.64 eV for 10 (Galeckas et al., 2024). Optical anisotropy in 11-Ga12O13 is negligible for 14, with 15 eV, but strong for 16, where polarization along 17 produces a 18 eV blue shift of the absorption edge (Galeckas et al., 2024).
Annealing studies connect these optical signatures to structural transformation. In 19 bilayers with 20 nm 21 on 22 23-Ga24O25, isochronal air anneals for 30 min from 300 to 26C revealed a two-stage kinetic behavior (Azarov et al., 2024). Stage I, below 27C, consists of epitaxial 28 regrowth at the buried interface; the 29 layer shrinks, for example from 30 nm to 31 nm after 32C, while its XRD 33 peak shifts closer to the tabulated 34, indicating improved crystallinity (Azarov et al., 2024). Stage II, beginning at 35–36C, is a non-planar transformation driven by misoriented 37 nucleation near the surface; the 38 peak vanishes by 39C and is replaced by 40 and 41 reflections at approximately 42 and 43, respectively (Azarov et al., 2024).
The optical changes mirror the structural transition. DRS Tauc analysis of the bilayers gave absorption edges around 4.38–4.43 eV at low-to-moderate annealing temperatures, with a threshold temperature range beyond which the bilayer becomes optically single-phase 44-like (Azarov et al., 2024). PL at 10 K showed the self-trapped-hole band near 3.2 eV, implantation-induced green luminescence near 2.5 eV, red luminescence near 2.0 eV, and an extrinsic Cr45 R-line doublet near 1.78 eV (Azarov et al., 2024). The red band intensified up to 46C and then shut off rapidly at higher temperatures, tracking the disappearance of 47 and indicating that the RL band is associated with defective spinel structure and its defects (Azarov et al., 2024).
These observations make phase identification in 48 heterostructures explicitly multimodal. XRD and RBS-C establish structure and thickness, DRS gives near-surface optical band-edge behavior, PL tracks intrinsic and defect-related recombination, and nano-FTIR resolves polymorph domains at the 49 nm scale (Galeckas et al., 2024, Azarov et al., 2024).
6. Radiation tolerance, disorder recovery, and polymorphic transformation under irradiation
Radiation response is a defining theme of this material system. The irradiation-derived 50 heterostructure was introduced as a crystalline product of disorder accumulation in 51-Ga52O53, rather than an endpoint of amorphization (Azarov et al., 2023). In room-temperature implantation experiments, 54 double-polymorph Ga55O56 structures tolerated disorder equivalent to hundreds of displacements per atom without severe degradation of crystallinity, with the 57-Ga58O59 channeling fingerprint persisting up to 265 dpa for 400 keV 60Ni61 implantation (Azarov et al., 2023). Classical MD correlated the 62 transformation with approximately 200 Ga Frenkel pairs per simulation cell, corresponding to 63 dpa (Azarov et al., 2023).
The mechanistic explanation is sublattice-specific. In 64-Ga65O66, the Ga sublattice is intrinsically defective and therefore relatively insensitive to added point defects, while the oxygen sublattice exhibits a strong tendency to recover FCC stacking after collision cascades (Azarov et al., 2023). Pearson-correlation analysis of partial radial distribution functions and O–O bond-angle distributions showed that the O sublattice retains fcc-characteristic peaks under Ni or Au excess but loses them for Ga excess, which drives local amorphization through charge transfer and disruption of ionic cohesion (Azarov et al., 2023). This establishes an important limitation: the exceptional radiation tolerance is not chemically unconditional.
Under swift heavy ion irradiation, the divergence between 67 and 68 becomes even more striking. In 69 heterostructures with 70 thicknesses of 350, 650, and 1000 nm on 71 72-Ga73O74, irradiation with Kr75 ions at 50, 100, and 147 MeV and fluences of 76 and 77 cm78 produced crystalline 79-phase tracks inside 80-Ga81O82, but no observable tracks in 83-Ga84O85 itself under identical conditions (Abdullaev et al., 9 Sep 2025). For 147 MeV Kr, the electronic stopping powers were essentially the same in the two phases, 86 eV/\AA\ in 87 and 2150 eV/\AA\ in 88, so the difference is not attributable to markedly different deposited electronic energy density (Abdullaev et al., 9 Sep 2025).
In 89-Ga90O91, tracks appear as crystalline 92-phase cylinders of approximately 3–4 nm diameter by STEM; MD showed initial amorphous cores of 93 nm at 200 ps that recrystallized into final tracks of 94 nm after long-term evolution (Abdullaev et al., 9 Sep 2025). In 95-Ga96O97, MD showed initial amorphous regions of 98 nm at 200 ps, but these were erased upon long-term evolution, leaving only residual point defects below STEM detection; even larger initial amorphous cores of 99 nm recovered (Abdullaev et al., 9 Sep 2025). Arrhenius extrapolation of oxygen-sublattice recovery yielded room-temperature half-lives of 00 s for 01 and 02 s for 03, with similar activation barriers of 04 and 05 eV, respectively (Abdullaev et al., 9 Sep 2025).
The thermal-transport signature is consistent with the microscopy. In irradiated bulk 06-Ga07O08, normalized thermal conductivity 09 dropped markedly and by more than half at the highest energy, whereas in 10-Ga11O12 13 remained nearly constant, indicating maintained homogeneity and the absence of extended defects (Abdullaev et al., 9 Sep 2025). This is counter to the usual correlation between low thermal conductivity and track formation, since 14-Ga15O16 has the lower thermal conductivity yet shows no observable tracks (Abdullaev et al., 9 Sep 2025).
A frequent oversimplification is to attribute radiation hardness here purely to reduced heat transport or to a generic “disordered structure.” The published interpretation is more specific: 17-Ga18O19 rapidly recovers disorder because its defective spinel cation sublattice offers multiple energetically accessible configurations, while the oxygen FCC backbone supplies a robust template for reordering (Abdullaev et al., 9 Sep 2025, Azarov et al., 2023). This suggests that radiation tolerance in these heterostructures is inseparable from polymorphic topology and sublattice dynamics.
7. Band alignment, strain sensitivity, and device-relevant considerations
Electronic band alignment at the 20-Ga21O22 interface has been addressed by density functional theory informed by experimental atomistic interface analysis (Liu et al., 2024). Because both sides consist of the same compound, core-level XPS alignment is not applicable, so the study used potential-lineup and vacuum-alignment methods. The valence-band offset in the lineup method is
23
and the conduction-band offset is obtained from the band gaps: 24 In vacuum alignment, absolute band edges are referenced to vacuum through the slab average potential (Liu et al., 2024).
The principal conclusion is that the offsets are small, type-II, and highly sensitive to strain and orientation. For Interface I, 25, the lineup was 26 eV; for Interface II, 27, 28 eV (Liu et al., 2024). In terms of conduction-band offsets 29, Interface I gave approximately 30 meV for strained 31 and 32 meV for unstrained 33, while Interface II gave approximately 34 meV for strained 35 and 36 meV for unstrained 37 (Liu et al., 2024). The corresponding valence-band offsets were also orientation- and strain-dependent, but in all computed cases electrons and holes preferred opposite sides of the interface, yielding type-II alignment (Liu et al., 2024).
Strain is decisive because the 38-phase conduction-band minimum has a strong deformation potential. The study reported 39 eV for the 40-oriented slab and cited 41 eV for 42 (Liu et al., 2024). When strain is applied to 43 only,
44
so compressive volumetric strain in 45 can move the CBO by several tenths of an eV (Liu et al., 2024). For the experimentally GPA-informed interfacial state, 46, and the extrapolated CBO shift was approximately 47 eV relative to the unstrained case, implying a net 48 eV for Interface I under that strain state (Liu et al., 2024). The authors cautioned that this value is indicative because strain relaxation and defects may modify both the average electrostatic potential and the band edges.
These results motivate discussion of a possible two-dimensional electron gas at the interface. The work states that 2DEG formation is theoretically possible, especially when 49 is sufficiently positive so that electrons localize on the 50 side, but that gradual strain relaxation with distance from the interface may shift the location of any electron gas or reduce its overall probability of formation (Liu et al., 2024). The study does not solve a Poisson–Schrödinger problem, but provides the standard triangular-well expression for subband energies as a context for confinement (Liu et al., 2024).
From a device standpoint, several implications recur across the literature. Type-II alignment and small offsets in the few-hundred-meV range suggest potential use in HEMTs, modulation-doped structures, tunnel barriers, or resonant transport elements within a single-chemistry platform, provided orientation and strain can be controlled (Liu et al., 2024). At the same time, the thermal and structural studies impose constraints. Annealing above 51–52C drives non-planar 53 conversion and texture evolution in bilayers (Azarov et al., 2024), while the radiation studies indicate that 54 is best placed on the exposed surface if track-free behavior under electronic-energy-dominant irradiation is desired (Abdullaev et al., 9 Sep 2025). The optical study further indicates that DRS with focal-plane control and nano-FTIR provide practical metrology for phase identification in such stacks (Galeckas et al., 2024).
A plausible synthesis of these findings is that 55-Ga56O57 heterostructures are most technically distinctive when their polymorphic contrast is used deliberately rather than incidentally: 58 can provide low thermal conductivity, nanoscale optical identifiability, and exceptional radiation tolerance, whereas 59 supplies anisotropic heat transport, structural stability, and the conventional electronic backbone. The outstanding challenge, as the band-offset study makes explicit, is that interfacial electronic functionality is not fixed by chemistry and therefore depends critically on local strain, orientation, and defect topology (Liu et al., 2024).