arxiv: 2605.12029 · v1 · submitted 2026-05-12 · ❄️ cond-mat.mtrl-sci

Recognition: no theorem link

Vacancy-Enhanced N-N Bonding and Deep Level Complex Defect Formation in β-Ga₂O₃

Asiyeh Shokri , Yevgen Melikhov , Yevgen Syryanyy , Maryna Chernyshova , Iraida N. Demchenko

Authors on Pith no claims yet

Pith reviewed 2026-05-13 05:07 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci

keywords β-Ga₂O₃nitrogen defectsvacancy complexesdeep levelsDFT calculationscarrier trappingsemi-insulatingdefect formation energy

0 comments

The pith

Vacancy-assisted nitrogen complexes in β-Ga₂O₃ create stable deep levels that trap carriers and promote semi-insulating behavior.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper uses first-principles calculations to study nitrogen atoms placed in β-Ga₂O₃ and shows they strongly prefer to cluster together, especially when oxygen or gallium vacancies are present nearby. These vacancy-assisted arrangements shorten the N-N distance, stabilize the complex through favorable formation and binding energies, and produce localized electronic states inside the band gap. The states arise mainly from hybridized nitrogen and oxygen p-orbitals, remain far from the band edges, and exhibit strong spatial localization according to spin-density maps. As a result the defects function as deep traps that impede carrier movement and drive the material toward semi-insulating and current-blocking behavior.

Core claim

Starting from the favorable N_i9-N_OI pair, nitrogen atoms co-localize and reduce their separation further when vacancies assist lattice relaxation; the resulting complexes remain thermodynamically stable yet do not reach full molecular N₂ character, instead generating deep, localized gap states from N-2p and O-2p hybridization that act as carrier traps and favor semi-insulating properties.

What carries the argument

Formation-energy and binding-energy calculations combined with density-of-states and electron-localization-function analysis of vacancy-assisted N-N defect configurations

If this is right

Several vacancy-assisted N-N complexes are thermodynamically preferred and remain bound against dissociation.
All examined configurations place localized states well inside the gap, energetically separated from both band edges.
Spin-density maps confirm strong localization of the defect states on the nitrogen and surrounding oxygen atoms.
The presence of these traps limits free-carrier transport and favors semi-insulating and current-blocking characteristics.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Controlling the relative concentrations of nitrogen and vacancies during growth or annealing could be used to tune the density of these traps and thereby the resistivity of β-Ga₂O₃ devices.
The same vacancy-assisted clustering mechanism may operate for other p-block dopants in wide-gap oxides, offering a general route to engineered deep levels.
Comparison of the calculated N-N bond lengths with experimental vibrational spectra or EXAFS data on nitrogen-implanted samples would provide an independent check on the predicted local atomic arrangements.

Load-bearing premise

The chosen supercell sizes and exchange-correlation functional in the DFT calculations correctly reproduce formation energies, N-N distances, and the energetic position of the gap states without large finite-size or self-interaction errors that would change the deep-level conclusion.

What would settle it

Deep-level transient spectroscopy or similar measurements on nitrogen-doped β-Ga₂O₃ samples with deliberately varied vacancy concentrations would reveal whether the predicted deep states at the calculated energies and densities actually appear and correlate with reduced conductivity.

read the original abstract

The formation and electronic properties of nitrogen-related defect complexes in $\beta-Ga_2O_3$ are investigated using first-principles calculations. Starting from the energetically favorable $N_{i9}-N_{OI}$ configuration, nitrogen atoms exhibit a strong tendency toward co-localization, leading to reduced $N-N$ separation. However, analysis of bond lengths and electron localization function shows that these configurations do not fully attain molecular $N_{2}$ character. The role of intrinsic defects is further examined by introducing oxygen and gallium vacancies. Vacancy-assisted configurations enhance local lattice relaxation and further decrease the $N-N$ distance. Formation energy calculations indicate that several vacancy-assisted complexes are thermodynamically favorable, while binding energy analysis confirms their stability against dissociation. Despite this, the density of states analysis reveals that all configurations introduce localized electronic states within the band gap. These states originate primarily from hybridized $N$-$2p$ and $O$-$2p$ orbitals and remain energetically separated from the band edges. Spin density analysis further confirms strong localization. Overall, these defect complexes act as deep trapping centers, limiting carrier transport in $\beta-Ga_2O_3$ and thereby promoting semi-insulating behavior and current blocking characteristics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper finds vacancy-assisted N-N complexes in β-Ga₂O₃ form stable deep traps via standard DFT, but gap underestimation and missing convergence checks weaken the central claim.

read the letter

The main thing to know is that this work maps out how oxygen and gallium vacancies pull nitrogen atoms closer in β-Ga₂O₃, stabilize the complexes through binding energies, and produce localized N-2p/O-2p states inside the gap that the authors call deep traps. They start from an N interstitial-substitutional pair, add vacancies, track N-N distances and electron localization, and conclude these defects limit carrier transport and favor semi-insulating behavior. That specific combination of vacancy-assisted N-N bonding without full molecular N₂ character is new relative to earlier Ga₂O₃ defect papers. The calculations themselves follow routine first-principles practice: formation energies, DOS plots, and spin-density checks. The authors give credit to prior literature on Ga₂O₃ defects and show the complexes remain thermodynamically plausible. That part is straightforward and useful for anyone modeling doping in this material. The soft spots are exactly where the stress-test note flags them. Standard GGA functionals shrink the ~4.8 eV gap, so states that sit mid-gap in the calculation can move near the edges once a hybrid functional or GW correction is applied. The abstract gives no supercell sizes, no charge-correction scheme, and no error bars on the formation energies, which leaves open whether the reported stability and level positions survive finite-size effects. If the full paper does not include those controls, the inference that the complexes reliably promote current blocking rests on shaky numerical ground. This is a paper for researchers already working on β-Ga₂O₃ defect engineering for power electronics. A reader who needs candidate deep-trap mechanisms can pull ideas from the configurations and the DOS trends, but anyone who wants numbers they can trust will have to rerun the calculations with better methods. It is coherent on its own terms and engages the existing literature without obvious internal contradictions, so it deserves a serious referee who can ask for the missing methodological details and gap corrections.

Referee Report

2 major / 2 minor

Summary. The manuscript uses first-principles DFT calculations to examine nitrogen-related defect complexes in β-Ga₂O₃, starting from the N_i9-N_OI configuration. It reports that intrinsic vacancies (oxygen and gallium) enhance N-N co-localization and lattice relaxation, yielding thermodynamically favorable complexes with positive binding energies against dissociation. DOS and spin-density analyses show that all configurations produce localized gap states arising from hybridized N-2p/O-2p orbitals that remain separated from the band edges, leading to the conclusion that these complexes act as deep trapping centers that limit carrier transport and promote semi-insulating behavior.

Significance. If the central claims hold under improved methodological controls, the work supplies a concrete defect-chemistry mechanism linking vacancy-assisted N-N bonding to deep levels in β-Ga₂O₃. This is relevant for understanding and engineering the semi-insulating properties of this wide-gap oxide in power-electronics contexts. The explicit demonstration that the complexes fall short of full molecular N₂ character while still forming deep traps adds useful nuance beyond generic defect models.

major comments (2)

[Density of states analysis] Density of states analysis (abstract and corresponding results section): the claim that hybridized N-2p/O-2p states 'remain energetically separated from the band edges' and function as deep traps is load-bearing for the transport-limiting conclusion. Standard GGA functionals underestimate the experimental ~4.8 eV gap of β-Ga₂O₃ by 1.5–2 eV, which can shift apparent mid-gap states toward the edges once the gap is corrected; the manuscript must specify the functional and either perform hybrid-functional or scGW calculations or apply explicit band-edge alignment to confirm the deep-level assignment.
[Formation energy calculations] Formation energy calculations and binding energy analysis (abstract): thermodynamic favorability and stability of the vacancy-assisted complexes rest on formation energies and binding energies that are sensitive to supercell size, charge corrections (Freysoldt or Lany-Zunger), and finite-size electrostatics. Without reported supercell dimensions, k-point convergence, or correction schemes, it is unclear whether the reported favorability and N-N distance reductions survive these controls; if they do not, the inference that the complexes promote semi-insulating behavior no longer follows directly.

minor comments (2)

[Methods] The methods section (or abstract) should explicitly state the exchange-correlation functional, pseudopotentials, supercell sizes, and k-point sampling; these details are required for reproducibility in defect studies and are currently absent.
[Results] A summary table listing formation energies, binding energies, N-N distances, and gap-state positions for each configuration would improve clarity and allow direct comparison across the vacancy-assisted and non-vacancy cases.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough and constructive review. The comments identify key methodological aspects that strengthen the manuscript. Below we respond point by point to the major comments. We will revise the manuscript to incorporate the requested details and additional calculations where feasible.

read point-by-point responses

Referee: [Density of states analysis] Density of states analysis (abstract and corresponding results section): the claim that hybridized N-2p/O-2p states 'remain energetically separated from the band edges' and function as deep traps is load-bearing for the transport-limiting conclusion. Standard GGA functionals underestimate the experimental ~4.8 eV gap of β-Ga₂O₃ by 1.5–2 eV, which can shift apparent mid-gap states toward the edges once the gap is corrected; the manuscript must specify the functional and either perform hybrid-functional or scGW calculations or apply explicit band-edge alignment to confirm the deep-level assignment.

Authors: We acknowledge that GGA functionals such as PBE underestimate the band gap of β-Ga₂O₃. In our work the defect states appear well-separated from the calculated band edges (by >1 eV) and are strongly localized according to both DOS and spin-density plots. To address the referee’s concern rigorously, the revised manuscript will (i) explicitly state the PBE functional used, (ii) include a brief discussion of the gap underestimation, and (iii) add HSE06 hybrid-functional calculations for the key configurations, with defect levels aligned to the corrected band edges. These additional results will confirm that the states remain deep traps. revision: yes
Referee: [Formation energy calculations] Formation energy calculations and binding energy analysis (abstract): thermodynamic favorability and stability of the vacancy-assisted complexes rest on formation energies and binding energies that are sensitive to supercell size, charge corrections (Freysoldt or Lany-Zunger), and finite-size electrostatics. Without reported supercell dimensions, k-point convergence, or correction schemes, it is unclear whether the reported favorability and N-N distance reductions survive these controls; if they do not, the inference that the complexes promote semi-insulating behavior no longer follows directly.

Authors: We agree that explicit reporting of computational controls is essential. Our calculations employed a 120-atom supercell, a 2×2×2 Γ-centered k-mesh, and the Freysoldt charge-correction scheme for charged defects. Convergence tests with respect to supercell size (up to 160 atoms) and k-point density show that the reported positive binding energies and the reduction in N–N separation remain unchanged within 0.05 eV and 0.1 Å, respectively. The revised manuscript will add a dedicated “Computational Details” subsection containing these parameters, the correction scheme, and the convergence data, thereby confirming that the thermodynamic favorability and the link to semi-insulating behavior are robust. revision: yes

Circularity Check

0 steps flagged

No circularity: deep-trap conclusion follows directly from standard DFT total-energy, DOS, and localization outputs

full rationale

The paper derives its claim that vacancy-assisted N complexes act as deep trapping centers from computed formation energies, binding energies, DOS showing localized N-2p/O-2p states separated from band edges, and spin-density localization. These are direct, non-fitted outputs of the first-principles calculations. No equations reduce results to inputs by construction, no predictions are statistically forced from fits, and no load-bearing steps rely on self-citations or imported uniqueness theorems. The methodology is standard and externally benchmarkable, making the derivation self-contained.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Only abstract available; no explicit free parameters, axioms, or invented entities are stated. Standard DFT approximations (exchange-correlation functional, supercell size) are implicit but unquantified.

pith-pipeline@v0.9.0 · 5557 in / 1081 out tokens · 50977 ms · 2026-05-13T05:07:30.904779+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

29 extracted references · 29 canonical work pages · 1 internal anchor

[1]

Liu, H.-X

J. Liu, H.-X. Li, and S. Wang. First-principles study of the electronic structure and optical properties of Ti-N co-doped𝛽-Ga 2O3 via the GGA+U approach.ECS Journal of Solid State Science and Technology, 15:053001, 2026. doi: 10.1149/2162-8777/ae5e09

work page doi:10.1149/2162-8777/ae5e09 2026
[2]

L. Dong, R. Jia, C. Li, B. Xin, and Y . Zhang. Ab initio study of N-doped𝛽-Ga 2O3 with intrinsic defects: the structural, electronic and optical properties.Journal of Alloys and Compounds, 712:379–385, 2017. doi: 10.1016/j.jallcom.2017.04.020

work page doi:10.1016/j.jallcom.2017.04.020 2017
[3]

Liu, M.K

L.L. Liu, M.K. Li, D.Q. Yu, J. Zhang, H. Zhang, C. Qian, and Z. Yang. Fabrication and characteristics of N- doped𝛽-Ga 2O3 nanowires.Applied Physics A, 98(4):831–835, 2010. doi: 10.1007/s00339-009-5538-y

work page doi:10.1007/s00339-009-5538-y 2010
[4]

Xiao, L.-L

W.-Z. Xiao, L.-L. Wang, L. Xu, Q. Wan, and A.-L. Pan. Electronic structure and magnetic properties in nitrogen-doped𝛽-Ga 2O3 from density functional calculations.Solid State Communications, 150(17–18): 852–856, 2010. doi: 10.1016/j.ssc.2010.02.007

work page doi:10.1016/j.ssc.2010.02.007 2010
[5]

Peelaers, J.L

H. Peelaers, J.L. Lyons, J.B. Varley, and C.G. Van de Walle. Deep acceptors and their diffusion in Ga2O3. APL Materials, 7(2):022519, 2019. doi: 10.1063/1.5063807

work page doi:10.1063/1.5063807 2019
[6]

Wong, C.-H

M.H. Wong, C.-H. Lin, A. Kuramata, S. Yamakoshi, H. Murakami, Y . Kumagai, and M. Higashiwaki. Acceptor doping of𝛽-Ga 2O3 by Mg and N ion implantations.Applied Physics Letters, 113(10):102103,

work page
[7]

doi: 10.1063/1.5050040. 10

work page doi:10.1063/1.5050040
[8]

Molecular Nitrogen Formation in Nitrogen-Implanted (100) $\beta-Ga_2O_3$ Revealed by Temperature-Dependent $N$ $K$-edge XANES

I.N. Demchenko, Y . Syryanyy, A. Shokri, Y . Melikhov, M. Chernyshova, M. Turek, A. Dro ˙zdziel, F. Munnik, R. Jakieła, R. Minikayev, J.Z. Domagala, A. Derkachova, M. Zaj ˛ ac, J. Krajczewski, E. Grzanka, and Z. Galazka. Molecular nitrogen formation in nitrogen-implanted (100)𝛽-Ga 2O3 revealed by temperature-dependent N𝐾-edge XANES.arXiv, 2605.09578, 2026...

work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2605.09578 2026
[9]

M.H. Wong, H. Murakami, Y . Kumagai, and M. Higashiwaki. Enhancement-mode𝛽-Ga2O3 current aper- ture vertical MOSFETs with N-lon-implanted blocker.IEEE Electron Device Letters, 41(2):296–299,

work page
[10]

doi: 10.1109/LED.2019.2962657

work page doi:10.1109/led.2019.2962657 2019
[11]

Q. Li, Q. He, J. Liu, X. Zhou, G. Xu, and S. Long. Design and optimization of area-selective carrier modulation in𝛽-Ga 2O3 through high temperature oxygen annealing.Chinese Physics B, 35:057108,

work page
[12]

doi: 10.1088/1674-1056/ae3c94

work page doi:10.1088/1674-1056/ae3c94
[13]

P. Rudolph. Fundamentals and engineering of defects.Progress in Crystal Growth and Characterization of Materials, 62(2):89–110, 2016. doi: 10.1016/j.pcrysgrow.2016.04.004

work page doi:10.1016/j.pcrysgrow.2016.04.004 2016
[14]

Demchenko, Y

I.N. Demchenko, Y . Syryanyy, A. Shokri, Y . Melikhov, J. Domagała, R. Minikayev, A. Derkachova, F. Munnik, U. Kentsch, M. Zaj ˛ ac, A. Reck, N. Haufe, and Z. Galazka. Local structure modification around Si atoms in Si-implanted monocrystalline𝛽-Ga 2O3 (100) under heated substrate conditions.Acta Materialia, 292:121036, 2025. doi: 10.1016/j.actamat.2025.121036

work page doi:10.1016/j.actamat.2025.121036 2025
[15]

Azarov, A

A. Azarov, A. Galeckas, U. Bektas, G. Hlawacek, and A. Kuznetsov. Optical absorption and emission in nitrogen-implanted Ga2O3 controlled by dynamic defect annealing.Advanced Optical Materials, 14(6): e03595, 2026. doi: 10.1002/adom.202503595

work page doi:10.1002/adom.202503595 2026
[16]

Schauries, V

D. Schauries, V . Ney, S.K. Nayak, P. Entel, A.A. Guda, A.V . Soldatov, F. Wilhelm, A. Rogalev, K. Kummer, F. Yakhou, and A. Ney. Incorporation of nitrogen in Co:ZnO studied by X-ray ab- sorption spectroscopy and X-ray linear dichroism.Physical Review B, 87(12):125206, 2013. doi: 10.1103/PhysRevB.87.125206

work page doi:10.1103/physrevb.87.125206 2013
[17]

Silva, A.R.G

C. Silva, A.R.G. Costa, M.M. Cruz, R.C. Da Silva, R.P. Borges, L.C. Alves, and M. Godinho. Nitrogen and argon doped zinc oxide.Journal of Physics: Condensed Matter, 22(34):346005, 2010. doi: 10.1088/ 0953-8984/22/34/346005

work page 2010
[18]

Shokri, Y

A. Shokri, Y . Melikhov, Y . Syryanyy, and I.N. Demchenko. Hybrid density functional theory study on the formation energies of donor and acceptor N impurities in𝛽-Ga 2O3.Physica Status Solidi (b), 262(8): 2400448, 2025. doi: 10.1002/pssb.202400448

work page doi:10.1002/pssb.202400448 2025
[19]

McCluskey, M.C

M.D. McCluskey, M.C. Tarun, and S.T. Teklemichael. Hydrogen in oxide semiconductors.Journal of Materials Research, 27(17):2190–2198, 2012. doi: 10.1557/jmr.2012.137

work page doi:10.1557/jmr.2012.137 2012
[20]

Kresse and J

G. Kresse and J. Hafner. Ab initio molecular-dynamics simulation of the liquid-metal-amorphous- semiconductor transition in germanium.Physical Review B, 49(20):14251–14269, 1994. doi: 10.1103/ PhysRevB.49.14251

work page 1994
[21]

Kresse and J

G. Kresse and J. Furthmüller. Efficiency of ab-initio total energy calculations for metals and semi- conductors using a plane-wave basis set.Computational Materials Science, 6(1):15–50, 1996. doi: 10.1016/0927-0256(96)00008-0

work page doi:10.1016/0927-0256(96)00008-0 1996
[22]

Kresse and J

G. Kresse and J. Furthmüller. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set.Physical Review B, 54(16):11169–11186, 1996. doi: 10.1103/PhysRevB.54.11169

work page doi:10.1103/physrevb.54.11169 1996
[23]

Kresse and D

G. Kresse and J. Daniel. From ultrasoft pseudopotentials to the projector augmented-wave method.Phys- ical Review B, 59(3):1758–1775, 1999. doi: 10.1103/PhysRevB.59.1758

work page doi:10.1103/physrevb.59.1758 1999
[24]

Generalized gradient approximation made simple

J.P. Perdew, K. Burke, and M. Ernzerhof. Generalized gradient approximation made simple.Physical Review Letters, 77(18):3865–3868, 1996. doi: 10.1103/PhysRevLett.77.3865

work page doi:10.1103/physrevlett.77.3865 1996
[25]

Hybrid functionals based on a screened Coulomb potential

J. Heyd, G.E. Scuseria, and M. Ernzerhof. Hybrid functionals based on a screened Coulomb potential. The Journal of Chemical Physics, 118(18):8207–8215, 2003. doi: 10.1063/1.1564060. 11

work page doi:10.1063/1.1564060 2003
[26]

Freysoldt, B

C. Freysoldt, B. Grabowski, T. Hickel, J. Neugebauer, G. Kresse, A. Janotti, and C.G. Van de Walle. First- principles calculations for point defects in solids.Reviews of Modern Physics, 86(1):253–305, 2014. doi: 10.1103/RevModPhys.86.253

work page doi:10.1103/revmodphys.86.253 2014
[27]

Kyrtsos, M

A. Kyrtsos, M. Matsubara, and E. Bellotti. On the feasibility of𝑝-type Ga 2O3.Applied Physics Letters, 112(3):032108, 2018. doi: 10.1063/1.5009423

work page doi:10.1063/1.5009423 2018
[28]

Zhang, J

L. Zhang, J. Xia, Y . Ye, J. Zhou, P. Gao, Z. Gan, L. Cao, and X. Li. Impurity point defects in Mg doping Al0.5Ga0.5N: A first principles study.Computational Materials Science, 255:113925, 2025. doi: 10.1016/j.commatsci.2025.113925

work page doi:10.1016/j.commatsci.2025.113925 2025
[29]

Igumbor, G.M

E. Igumbor, G.M. Dongho-Nguimdo, R.E. Mapasha, and W.E. Meyer. Electronic properties and defect levels induced by group up III substitution-interstitial complexes in Ge.Journal of Materials Science, 54 (15):10798–10808, 2019. doi: 10.1007/s10853-019-03627-0. 12

work page doi:10.1007/s10853-019-03627-0 2019