Structure and thermodynamic stability of β-Ga2O3 surfaces
Abstract
We present a comprehensive first-principles investigation of all symmetrically inequivalent low-index surfaces of β-Ga2O3, examining their structural properties and thermodynamic stability across experimentally relevant growth conditions. Using density-functional theory with both the semi-local functional PBEsol and a modified PBE0 hybrid functional with 26% exact exchange, denoted PBE0(0.26), we calculate surface free energies for the (010), (100), (001), (2̅01), (110), (111), and (111̅) orientations, including the effects of harmonic vibrational contributions and varying oxygen chemical potentials. We demonstrate that the energetic ordering remains consistent across computational approaches and that the vibrational contributions remain below 0.2 J/m2 up to temperatures of 1000 K. A coordination-based model that correlates surface stability with the density of under-coordinated atoms reveals that under-coordinated oxygen atoms and tetrahedral Ga sites substantially destabilize surfaces, while exposed under-coordinated octahedral Ga atoms serve as indicators of surface stability. Our thermodynamic analysis shows that stoichiometric terminations dominate over nearly the entire range of chemical potentials relevant for β-Ga2O3 stability, while non-stoichiometric terminations emerge only under extreme reducing or oxidizing conditions. Notably, we predict the formation of stable Ga-rich terminations resembling Ga adlayers for the (100) and (2̅01) surfaces under highly reducing conditions.
Versions
Related
- Thermodynamic stability of the new black hole solution in f(R,G,T) gravity
- Thermodynamic Stability of Classical ℤ2 Lattice Gauge Models
- Kernel moments and thermodynamic stability in nonlocal Nambu--Jona-Lasinio models
- Thermodynamic stability of twisted domains in AgCrSe2 thin films grown on lattice-matched YSZ(111) substrate
- Charting the thermodynamic stability of hybrid perovskite alloys with machine learning
BibTeX
@misc{lion2026structure,
title = {{Structure and thermodynamic stability of $β$-Ga$_2$O$_3$ surfaces}},
author = {Konstantin Lion and Claudia Draxl},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2607.27424},
eprint = {2607.27424},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2607.27424},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Structure and thermodynamic stability of β-Ga_2O_3 surfaces
AU - Konstantin Lion
AU - Claudia Draxl
PY - 2026
DA - 2026/09/19
PB - arXiv
DO - 10.48550/arxiv.2607.27424
UR - https://arxiv.org/abs/2607.27424
AB - We present a comprehensive first-principles investigation of all symmetrically inequivalent low-index surfaces of β-Ga_2O_3, examining their structural properties and thermodynamic stability across experimentally relevant growth conditions. Using density-functional theory with both the semi-local functional PBEsol and a modified PBE0 hybrid functional with 26% exact exchange, denoted PBE0(0.26), we calculate surface free energies for the (010), (100), (001), (2̅01), (110), (111), and (111̅) orientations, including the effects of harmonic vibrational contributions and varying oxygen chemical potentials. We demonstrate that the energetic ordering remains consistent across computational approaches and that the vibrational contributions remain below 0.2 J/m^2 up to temperatures of 1000 K. A coordination-based model that correlates surface stability with the density of under-coordinated atoms reveals that under-coordinated oxygen atoms and tetrahedral Ga sites substantially destabilize surfaces, while exposed under-coordinated octahedral Ga atoms serve as indicators of surface stability. Our thermodynamic analysis shows that stoichiometric terminations dominate over nearly the entire range of chemical potentials relevant for β-Ga_2O_3 stability, while non-stoichiometric terminations emerge only under extreme reducing or oxidizing conditions. Notably, we predict the formation of stable Ga-rich terminations resembling Ga adlayers for the (100) and (2̅01) surfaces under highly reducing conditions.
ER -