REVIEW 1 major objections 6 minor 7 references
Why Do Thick MOCVD-Grown beta-Ga2O3 Epilayers on (001) Substrates Crack: Crystallographic Origin
T0 review · 1 major / 6 minor · reviewed 2026-07-14 · grok-4.5
Pith's one-line read Thick MOCVD beta-Ga2O3 films on (001) substrates grow as (-401) layers and crack from ~1.8 um because of a ~4.1% tensile mismatch along [104].
desk verdict Solid crystallographic explanation for why thick MOCVD Ga2O3 on (001) cracks: the film is actually (−401) with ~4.1% tensile mismatch along one in-plane axis. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The crystallographic orientation relationship between the (-401) epilayer and (001) substrate: [010]epi || [010]sub and [104]epi || [100]sub. It converts an intended homoepitaxial growth into a lattice-mismatched heteroepitaxy whose anisotropic strain drives cracking.
What would settle it
Cross-sectional TEM or STEM of the earliest stages of growth (or of the annealed-only control) that either identifies low-index facets matching the (-401)/(001) relation or shows continuous (001) registry would confirm or refute the proposed selection mechanism.
Extended reading notes
Core claim
Despite growth on nominally (001) substrates under conventional oxygen-rich MOCVD conditions, the epilayers adopt a predominantly (-401)-oriented structure from the earliest stages. The resulting epitaxial alignment [010]epi || [010]sub and [104]epi || [100]sub produces zero lattice mismatch along [010] but a theoretical tensile in-plane strain of approximately +4.1% along [104], which relaxes by forming channeling cracks parallel to [010]epi once thickness exceeds ~1.8 um.
Load-bearing premise
The claim that oxygen-annealed faceting of the (001) surface is what selects the (-401) orientation rests on AFM topography and surface-energy arguments; the paper does not show the actual nucleation interface or facet indices by cross-sectional imaging.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript reports a systematic thickness series (0.3–3.5 µm) of MOCVD β-Ga₂O₃ epilayers grown at ~3.5 µm h⁻¹ on (001) native substrates. AFM shows striated morphology aligned with [010]sub and increasing RMS roughness beyond coalescence; optical microscopy shows channeling cracks parallel to [010]epi above ~1.8 µm. HRXRD ω–2θ scans establish that the epilayers are predominantly (−401)-oriented from the earliest stages (with −401/−802/−1203 intensity growing and no broad 002 epilayer peak), while φ-scans of the −201 reflection establish the single-domain in-plane relation [010]epi ∥ [010]sub and [104]epi ∥ [100]sub. Rocking-curve FWHMs indicate higher twist than tilt, both decreasing with thickness. The authors calculate a geometric tensile mismatch ε[104] ≈ +4.1% and argue that this anisotropic strain drives the observed cracking; they further hypothesize that oxygen-rich annealing produces a faceted (001) surface that templates (−401) growth.
Significance. If the orientation–mismatch–cracking chain holds, the work supplies a concrete crystallographic explanation for a known industrial bottleneck: thick MOCVD layers on the largest, lowest-cost (001) β-Ga₂O₃ wafers crack and roughen. That explanation is actionable—it identifies the epitaxial relationship itself, rather than generic thermal stress or growth-rate effects, as the limiting factor—and therefore guides nucleation strategies (temperature/ambient ramps, low-T nucleation layers) that could enable device-relevant thick (001) films without AlGaO buffers that impede vertical transport. Strengths include a well-controlled thickness series, mutually consistent AFM/optical/HRXRD datasets, and a transparent geometric strain calculation using standard lattice periodicities. The reconstruction hypothesis is weaker and is already labeled as such, but it is not required for the central claim.
major comments (1)
- The central claim (orientation + anisotropic mismatch → cracking) is well supported by the thickness series, ω–2θ, φ-scans, rocking curves, AFM striation orientation, and optical crack direction/critical thickness. No load-bearing inconsistency is present. The only material gap is the absence of cross-sectional TEM/STEM of the interface (Discussion / Fig. 7), which would confirm facet indices and the nucleation pathway; the authors already present the reconstruction as a hypothesis, so this is not required to accept the cracking-origin claim, but a brief statement of what TEM would (and would not) change would strengthen the Discussion.
minor comments (6)
- Abstract and Results: the phrase “in theory, a maximum tensile in-plane strain of approximately +4.1%” is slightly ambiguous; state explicitly that this is the geometric (unrelaxed) mismatch and that the actual elastic strain is partially relaxed by defects and cracks.
- Fig. 4(a) and associated text: the progressive −401 peak shift to lower 2θ is used as evidence of tensile relaxation; a short quantitative estimate of Δd/d (or residual strain) from the measured 2θ values would make the argument more precise.
- Fig. 5: the anomalous smaller FWHM of the −401 reflection measured coplanar with [104]epi for the 0.3 µm sample is noted but only briefly discussed; a sentence linking this to incomplete defect-mediated relaxation before cracking would help the reader.
- Supporting Information Fig. S2 is useful; consider moving a condensed version of the “no broad 002 epilayer peak” argument into the main text so that the single-domain (−401) conclusion is self-contained.
- Materials and Methods: growth rates are taken from Si-doped calibration runs; a short note that undoped layers under the same conditions give the same rate (or that SIMS marker runs were co-loaded) would remove any residual ambiguity.
- References: a few recent (001) MOCVD morphology papers are cited; ensuring the most recent HVPE vs MOCVD comparison reviews are included would help place the work for non-specialists.
Circularity Check
No significant circularity: orientation, epitaxial relation, and mismatch strain are measured or computed from independent crystallographic inputs, not forced by construction or self-citation.
full rationale
The paper's load-bearing chain is experimental: HRXRD ω–2θ scans (Figs. 3–4, S1–S2) establish the dominant (−401) epilayer orientation from the thinnest films onward; skew-symmetric φ-scans of the −201 reflection fix the single-domain in-plane relation [010]epi ∥ [010]sub and [104]epi ∥ [100]sub; optical microscopy and AFM then show cracks and striations oriented consistently with that relation; the theoretical tensile strain ε[104] ≈ +4.1 % is obtained by arithmetic on standard monoclinic lattice periodicities (L[104] = 23.464 Å, L[100] = 12.210 Å) that are not fitted to the cracking data. Progressive −401 peak shifts and rocking-curve narrowing are reported as observations, not as predictions renormalized to the same dataset. The faceted-reconstruction hypothesis (Fig. 7) is explicitly labeled as such and is not required for the mismatch–crack claim. Minor citations to the authors’ prior GaN buffer work appear only as background analogies for thermal-strain arguments that the paper itself discards; they do not underwrite the crystallographic origin. No self-definitional loop, fitted-input-as-prediction, uniqueness import, or ansatz smuggling is present. The derivation is therefore self-contained against external crystallographic benchmarks.
Assumptions & free parameters
free parameters (2)
- Lattice periodicities L[104] and L[100] used for geometric mismatch =
L[104]=23.464 Å, L[100]=12.210 Å → ε≈+4.1%
- Critical thickness for cracking (~1.8 µm) =
~1.8 µm
assumptions (4)
- domain assumption β-Ga2O3 is monoclinic with the standard lattice and the listed plane/direction relationships between (001) and (−401).
- domain assumption Griffith fracture criterion: cracking becomes favorable when elastic-energy release exceeds new-surface energy.
- domain assumption Surface energies of β-Ga2O3 are highly anisotropic, with (001) among the highest (~1.95 J m−2), so the surface reconstructs/facets under oxygen-rich conditions.
- ad hoc to paper Thermal-mismatch strain is not the dominant crack driver under the present conditions.
invented entities (1)
-
Oxygen-driven faceted (001) reconstruction that templates (−401) epilayer growth
Cite this review
Pith. "Pith review of Why Do Thick MOCVD-Grown beta-Ga2O3 Epilayers on (001) Substrates Crack: Crystallographic Origin." pith.science (2026). https://pith.science/paper/OJBFZSSP
@misc{pith2026260710378,
author = {Pith},
title = {Pith review of: Why Do Thick MOCVD-Grown beta-Ga2O3 Epilayers on (001) Substrates Crack: Crystallographic Origin},
year = {2026},
howpublished = {\url{https://pith.science/paper/OJBFZSSP}},
note = {Machine review of arXiv:2607.10378}
}
read the original abstract
Thick, defect free epitaxial layers grown using industry standard techniques are a fundamental requirement for the widespread adoption of fully vertical power devices based on ultra wide bandgap gallium oxide (Ga2O3). However, metal-organic chemical vapour deposition (MOCVD) of such layers on native beta-Ga2O3 substrates with the largest diameter (001) orientation remains relatively unexplored, and the origins of the reported surface roughening and cracking with increasing thickness are not yet fully understood. To address this, we report a systematic study of MOCVD grown beta-Ga2O3 epilayers deposited at growth rates of ~3.5 um/h, with thicknesses from 0.3 to 3.5 um. The epilayers exhibit a relatively smooth but striated surface morphology, with progressively increasing nanometre-scale roughness beyond coalescence and crack formation observed from ~1.8 um thickness. High resolution X ray diffraction reveals that, despite growth on (001) substrates, the epilayers adopt a predominantly (-401)-oriented structure from the earliest stages of growth. Rocking curve analysis further indicates a higher degree of in-plane twist than tilt, both decreasing with increasing epilayer thickness. While the epilayer and substrate are lattice-matched along the [010] in plane direction, the epitaxial alignment in the orthogonal epilayer [104] in plane direction imposes, in theory, a maximum tensile in-plane strain of approximately +4.1% arising from the underlying lattice mismatch, thereby driving crack formation perpendicular to this direction. Our results suggest that this epitaxial relationship is likely associated with faceted reconstruction of the (001) substrate surface during annealing, driven by the minimisation of surface energy under oxygen-rich MOCVD growth conditions.
Reference graph
Works this paper leans on
-
[1]
β‐Gallium Oxide Devices: Progress and Outlook,
1 M. Higashiwaki, “ β‐Gallium Oxide Devices: Progress and Outlook,” Physica Status Solidi (RRL) – Rapid Research Letters 15(11), (2021). 14 2 M.J. Tadjer, “Toward gallium oxide power electronics,” Science (1979). 378(6621), 724–725 (2022). 3 Y. Wang, M. Zhu, and Y. Liu, “Growth process, defects, and dopants of bulk β-Ga2O3 semiconductor single crystals,” ...
2021
-
[2]
Toward emerging gallium oxide semiconductors: A roadmap,
Zhao, Z. Jia, J. Ye, J. Zhang, S. Long, X. Tao, R. Zhang, and Y. Hao, “Toward emerging gallium oxide semiconductors: A roadmap,” Fundamental Research 1(6), 697–716 (2021). 6 M. Higashiwaki, K. Sasaki, K. Goto, K. Nomura, Q.T. Thieu, R. Togashi, H. Murakami, Y
2021
-
[3]
Ga2O3 Schottky barrier diodes with n -Ga2O3 drift layers grown by HVPE,
Kumagai, B. Monemar, A. Koukitu, A. Kuramata, and S. Yamakoshi, “Ga2O3 Schottky barrier diodes with n -Ga2O3 drift layers grown by HVPE,” in 2015 73rd Annual Device Research Conference (DRC), (IEEE, 2015), pp. 29–30. 7 Z. Hu, K. Nomoto, W. Li, N. Tanen, K. Sasaki, A. Kuramata, T. Nakamura, D. Jena, and H.G. Xing, “Enhancement-Mode Ga2O3 Vertical Transisto...
2015
-
[4]
Ultra-wide bandgap semiconductor Ga2O3 power diodes,
Kong, H. Zhou, and Y . Hao, “Ultra-wide bandgap semiconductor Ga2O3 power diodes,” Nat. Commun. 13(1), 3900 (2022). 11 W. Xu, Y. Chen, X. Zhou, X. Xu, X. Dong, H. Long, Z. Fang, Z. Yang, D. Chen, H. Qi, and K.H.L. Zhang, “Recent advancement and perspective of epitaxial growth and doping of β- Ga2O3 thin films for power electronics,” APL Mater. 13(10), (20...
2022
-
[5]
Homoepitaxial growth of β-Ga2O3 layers by halide vapor phase epitaxy,
Kumagai, M. Higashiwaki, A. Kuramata, S. Yamakoshi, B. Monemar, and A. Koukitu, “Homoepitaxial growth of β-Ga2O3 layers by halide vapor phase epitaxy,” Applied Physics Express 8(1), 015503 (2015). 15 15 J.H. Leach, K. Udwary, J. Rumsey, G. Dodson, H. Splawn, and K.R. Evans, “Halide vapor phase epitaxial growth of β-Ga2O3 and α-Ga2O3 films,” APL Mater. 7(2...
2015
-
[6]
High -quality (001) β-Ga2O3 homoepitaxial growth by metalorganic chemical vapor deposition enabled by in situ indium surfactant,
Mudiyanselage, H. Fu, and B. Zhang, “High -quality (001) β-Ga2O3 homoepitaxial growth by metalorganic chemical vapor deposition enabled by in situ indium surfactant,” Appl. Phys. Lett. 120(21), (2022). 20 Y. He, Y . Liu, X. Lu, Z. Wang, X. Song, Y . Zhou, X. Zheng, X. Ma, and Y . Hao, “(001) β- Ga2O3 epitaxial layer grown with in-situ pulsed Al atom assis...
2022
-
[7]
Evolution of planar defects during homoepitaxial growth of β-Ga2O3 layers on (100) substrates—A quantitative model,
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(22), (2016). 26 A.F.M. Anhar Uddin Bhuiyan, Z. Feng, J.M. Johnson, H.-L. Huang, J. Hwang, and H. Zhao, “MOCVD Epitaxy of Ultrawide Bandgap β-(AlxGa1–x)2O3 with High-Al Comp...
2016
Reviewed July 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.