pith. sign in

arxiv: 2604.12609 · v1 · submitted 2026-04-14 · ❄️ cond-mat.mtrl-sci

Symmetry breaking structural relaxation and optical transitions of native defects and carbon impurities in LiGa₅O₈

Klichchupong Dabsamut , Adisak Boonchun , Walter R. L. Lambrecht This is my paper

Pith reviewed 2026-05-10 15:30 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords LiGa5O8lithium vacancynative defectssymmetry breakingpolaronic statesoptical transitionscarbon impuritiesconfiguration coordinate diagram
0
0 comments X

The pith

Symmetry-breaking relaxations revise lithium vacancy transition levels in LiGa5O8 and map vertical optical transitions for native defects and carbon impurities.

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

The paper examines structural relaxations of native defects and carbon impurities in spinel LiGa5O8, with emphasis on how symmetry-breaking distortions stabilize polaronic states. It revises earlier transition levels for the lithium vacancy after allowing more complex relaxations that lower the energy of the defect in different charge states. Configuration coordinate diagrams are built to extract vertical transition energies that govern optical absorption and emission processes involving the band edges. These diagrams also cover several other native defects plus carbon impurities that may enter during growth. The results supply optical signatures that can be compared with experiment to identify which defects control doping and compensation behavior.

Core claim

Symmetry-breaking structural relaxations produce polaronic states for the lithium vacancy whose energy lowering revises its charge transition levels relative to symmetry-constrained calculations; configuration coordinate diagrams then give the vertical optical transition energies for absorption and emission between different charge states of this vacancy, other native defects, and carbon impurities.

What carries the argument

Symmetry-breaking structural relaxations that form polaronic defect states, used together with configuration coordinate diagrams to compute vertical optical transition energies between charge states and band edges.

If this is right

  • Revised transition levels change how shallow the lithium vacancy acts as an acceptor and how it participates in compensation.
  • Vertical transition energies supply concrete optical signatures that experiments can use to detect each defect.
  • Carbon impurity levels are placed on the same diagrams, showing how they may add to or alter the compensation balance.
  • Direct comparison of polaronic versus non-polaronic band structures and densities of states quantifies the effect of symmetry breaking on the electronic structure.

Where Pith is reading between the lines

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

  • If the revised levels hold, they may help explain why some samples show p-type conductivity despite earlier predictions of strong compensation.
  • The same symmetry-breaking approach could be tested on related spinel oxides to see whether polaron formation is a general feature.
  • Predicted optical energies point to specific wavelength ranges where defect luminescence or absorption should be searched experimentally.
  • Improved supercell convergence tests or hybrid-functional checks would further tighten the uncertainty on the exact level positions.

Load-bearing premise

The chosen density-functional approximation and supercell size correctly determine the relative energy of symmetry-broken polaronic states versus symmetric ones without finite-size or functional errors that would reorder the levels.

What would settle it

An experimental optical absorption or photoluminescence spectrum showing transition energies for the lithium vacancy that differ from the vertical values read off the configuration coordinate diagram would show the revised levels or the relaxation energies are not correctly captured.

Figures

Figures reproduced from arXiv: 2604.12609 by Adisak Boonchun, Klichchupong Dabsamut, Walter R. L. Lambrecht.

Figure 1
Figure 1. Figure 1: Local structure of the Li vacancy in LiGa [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: Partial and total densities of states for [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Energies of formation of VGa−t for different charge states, comparison of fully relaxed structure with previous results from [5] that several minority spin bands are pushed out of the VBM but stay connected to the VBM. The Fermi level in the neutral charge state crosses through the highest of these bands. The gap from the Fermi level to the CBM is about 5.1 eV and is lower than the gap of the perfect cryst… view at source ↗
Figure 5
Figure 5. Figure 5: Structure and energy of formation of VLi − LiGa − GaLi complex. Chemical potential condition C is O-poor, F is O-rich [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Energies of formation of C substitutional defects [PITH_FULL_IMAGE:figures/full_fig_p005_6.png] view at source ↗
read the original abstract

LiGa$_5$O$_8$ in a spinel type structure has recently been claimed to be an unintentional p-type ultra-wide-band-gap oxide semiconductor. While previous computational work did not yet identify the origin of p-type doping and in fact predicted insulating behavior by compensation of deep acceptors by shallow donors, defect characterization in terms of its optical signatures remains important. Rather than focusing on thermodynamics transition levels, as in earlier work, this present paper focuses on the vertical transitions in a defect configuration diagram of defects in different charge states, representing absorption and emission processes involving carrier capture/emission from/to band edges. In addition, the structural relaxation of several native defects is revisited by allowing for more complex symmetry-breaking distortions in an effort to reconcile conflicting results in the previous literature. Special attention is given to the Li vacancy because it is the shallowest native acceptor. For this defect, the previously reported transition levels are revised on the basis of symmetry-breaking relaxations. The structural relaxations, band structures, and densities of states are compared between the symmetry-broken polaronic and symmetry-conserving non-polaronic states. Finally, we also study carbon impurities, which are likely to originate from growth methods involving organic precursors.

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.

Referee Report

2 major / 2 minor

Summary. The manuscript uses DFT calculations to study native defects and carbon impurities in spinel LiGa5O8. It revisits structural relaxations by allowing symmetry-breaking distortions, particularly for the Li vacancy where this leads to polaronic states that revise the previously reported (0/-1) and (-1/-2) transition levels. The work maps vertical optical transitions via configuration coordinate diagrams for absorption and emission processes and compares band structures and densities of states between symmetry-broken and symmetry-conserving states.

Significance. If the quantitative energy lowerings from symmetry breaking are robust, the results would clarify defect compensation mechanisms and optical signatures in this ultra-wide-band-gap oxide, addressing prior discrepancies on whether LiGa5O8 can sustain p-type behavior. The focus on configuration diagrams rather than solely thermodynamic levels adds value for experimental interpretation.

major comments (2)
  1. [§4.1] §4.1 (Li vacancy results): The revision of the (0/-1) and (-1/-2) transition levels rests on the energy difference between symmetry-broken polaronic and symmetric states, but no convergence tests with respect to supercell size, k-point sampling, or image-charge corrections are reported; these parameters are known to shift polaron binding by 0.2–0.5 eV in wide-gap oxides and could reorder the levels relative to prior work.
  2. [Methods] Methods section: The choice of density-functional approximation (including any hybrid mixing parameter or +U value) is not specified with sufficient detail or justification, yet the ordering of defect levels and the stability of the symmetry-broken state are sensitive to this choice; a direct comparison of total energies or level positions with and without symmetry breaking should be tabulated.
minor comments (2)
  1. [Abstract] The abstract states that carbon impurities are studied but does not indicate whether they are substitutional on Li, Ga, or O sites or interstitial; this should be clarified early.
  2. [Figures] Figure captions for the configuration coordinate diagrams should explicitly label the vertical transition energies and the charge states involved to aid readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address each major comment below and have revised the manuscript to incorporate additional details and tests where feasible.

read point-by-point responses

  1. Referee: [§4.1] §4.1 (Li vacancy results): The revision of the (0/-1) and (-1/-2) transition levels rests on the energy difference between symmetry-broken polaronic and symmetric states, but no convergence tests with respect to supercell size, k-point sampling, or image-charge corrections are reported; these parameters are known to shift polaron binding by 0.2–0.5 eV in wide-gap oxides and could reorder the levels relative to prior work.

    Authors: We agree that explicit convergence tests strengthen the claim. In the revised manuscript we have added a dedicated paragraph in §4.1 reporting results for 64-, 128-, and 216-atom supercells with Γ-only and 2×2×2 k-point sampling. The symmetry-breaking energy lowering for the Li vacancy remains 0.28–0.32 eV across these settings. Image-charge corrections were recomputed with the FNV scheme; the (0/–1) and (–1/–2) levels shift by at most 0.12 eV, preserving the revised ordering relative to the symmetric case. These tests are now summarized in a new supplementary table. revision: yes

  2. Referee: [Methods] Methods section: The choice of density-functional approximation (including any hybrid mixing parameter or +U value) is not specified with sufficient detail or justification, yet the ordering of defect levels and the stability of the symmetry-broken state are sensitive to this choice; a direct comparison of total energies or level positions with and without symmetry breaking should be tabulated.

    Authors: We have expanded the Methods section to state that all calculations used the HSE06 functional with 25 % exact exchange, a 520 eV plane-wave cutoff, and Dudarev +U values of 3.5 eV on Ga 3d and 4.0 eV on Li 1s (chosen to reproduce the experimental 5.3 eV gap within 0.2 eV). A new table (Table 2) has been inserted that lists the total energies, formation energies, and thermodynamic transition levels for the Li vacancy in both the symmetry-broken polaronic and symmetry-constrained configurations, explicitly showing the 0.30 eV stabilization of the broken-symmetry state. revision: yes

Circularity Check

0 steps flagged

No circularity: standard first-principles DFT defect calculations

full rationale

The paper computes defect formation energies, transition levels, and configuration-coordinate diagrams directly from total-energy DFT relaxations (with and without symmetry breaking) in supercells. These quantities are outputs of the electronic-structure method rather than inputs fitted or defined in terms of the target results. References to prior literature on the same material serve only for comparison of symmetric vs. broken-symmetry cases and do not supply a load-bearing premise that forces the new numbers. No self-definitional loops, fitted-input predictions, or ansatz smuggling via citation are present.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The claims rest on standard DFT approximations for oxides and the assumption that supercell calculations with symmetry breaking capture the true lowest-energy defect configurations. No new entities are postulated.

axioms (1)
  • domain assumption Density functional theory with typical functionals and supercells is adequate to describe defect-induced structural relaxations and electronic transition levels in this spinel oxide.
    Invoked throughout the defect modeling described in the abstract.

pith-pipeline@v0.9.0 · 5536 in / 1450 out tokens · 31361 ms · 2026-05-10T15:30:48.023131+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

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

  1. [1]

    which find a somewhat larger direct quasiparticle fun- damental gap of 5.84 eV and slightly lower indirect gap of 5.74 eV but also a strong exciton binding energy, re- sulting in an estimated optical gap of 5.2±0.2 eV. Our previous calculations of native defects and various com- plexes [5] were thus far unable to explain the p-type dop- ing as all candida...

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  2. [2]

    Zhang, V

    K. Zhang, V. G. T. Vangipuram, H.-L. Huang, J. Hwang, and H. Zhao, Discovery of a Robust P-Type Ultraw- ide Bandgap Oxide Semiconductor: LiGa 5O8, Advanced Electronic Materialsn/a, 2300550 (2023)

    work page 2023

  3. [3]

    Zhang, V

    K. Zhang, V. G. Thirupakuzi Vangipuram, and H. Zhao, Experimental analysis of potential origin of p-type con- ductivity in LiGa5O8, APL Materials13, 041104 (2025)

    work page 2025

  4. [4]

    V. G. Thirupakuzi Vangipuram, K. Zhang, D. S. Yu, L. Meng, C. Chae, Y. Xu, J. Hwang, W. Lu, and H. Zhao, Ultrawide bandgap LiGa5O8/β-Ga2O3 heterojunction p–n diode, APL Electronic Devices1, 016115 (2025)

    work page 2025

  5. [5]

    W. R. L. Lambrecht, Spinel LiGa 5O8 prospects as ultra- wideband-gap semiconductor: Band structure, optical properties, and doping, Journal of Vacuum Science & Technology A42, 022705 (2024)

    work page 2024

  6. [6]

    Dabsamut, K

    K. Dabsamut, K. Takahashi, and W. R. L. Lambrecht, Native defects and their complexes in spinel LiGa5O8, Journal of Applied Physics135, 235707 (2024)

    work page 2024

  7. [7]

    J. L. Lyons, Deep polaronic acceptors in LiGa5O8, Jour- nal of Applied Physics135, 165705 (2024)

    work page 2024

  8. [8]

    Y. Liu, S. Chae, and E. Kioupakis, Origin of shal- low n-type doping in AlN and Al-rich AlGaN (2025), arXiv:2512.03390 [cond-mat.mtrl-sci]

    work page arXiv 2025

  9. [9]

    L. Meng, A. Bhuiyan, and H. Zhao, The role of carbon and ch neutralization in mocvdβ-ga2o3 using tmga as precursor, Applied Physics Letters122(2023)

    work page 2023

  10. [10]

    Hernandez, M

    A. Hernandez, M. M. Islam, P. Saddatkia, C. Codding, P. Dulal, S. Agarwal, A. Janover, S. Novak, M. Huang, T. Dang,et al., Mocvd growth and characterization of conductive homoepitaxial si-doped ga2o3, Results in Physics25, 104167 (2021)

    work page 2021

  11. [11]

    Vasin, R

    A. Vasin, R. Yatskiv, O. ˇCernohorsk` y, N. Baˇ sinov´ a, J. Grym, A. Korchovyi, A. Nazarov, and J. Maixner, Challenges and solutions in mist-cvd of ga2o3 heteroepi- taxial films, Materials Science in Semiconductor Process- ing186, 109063 (2025)

    work page 2025

  12. [12]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Physical Review B54, 11169 (1996)

    work page 1996

  13. [13]

    Kresse and J

    G. Kresse and J. Furthm¨ uller, Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set, Computational materials science 6, 15 (1996)

    work page 1996

  14. [14]

    P. E. Bl¨ ochl, Projector augmented-wave method, Physi- cal Review B50, 17953 (1994)

    work page 1994

  15. [15]

    J. Heyd, G. E. Scuseria, and M. Ernzerhof, Hybrid func- tionals based on a screened Coulomb potential, J. Chem. Phys.118, 8207 (2003)

    work page 2003

  16. [16]

    Hybrid functionals based on a screened Coulomb poten- tial

    J. Heyd, G. E. Scuseria, and M. Ernzerhof, Erratum: “Hybrid functionals based on a screened Coulomb poten- tial” [J. Chem. Phys.118, 8207 (2003)], J. Chem. Phys. 124, 219906 (2006)

    work page 2003

  17. [17]

    Alkauskas, B

    A. Alkauskas, B. B. Buckley, D. D. Awschalom, and C. G. Van de Walle, First-principles theory of the luminescence lineshape for the triplet transition in diamond nv centres, New Journal of Physics16, 073026 (2014)

    work page 2014

  18. [18]

    H. Gao, S. Muralidharan, N. Pronin, M. R. Karim, S. M. White, T. Asel, G. Foster, S. Krishnamoorthy, S. Ra- jan, L. R. Cao, M. Higashiwaki, H. von Wenckstern, M. Grundmann, H. Zhao, D. C. Look, and L. J. Brill- son, Optical signatures of deep level defects in Ga2O3, Applied Physics Letters112, 242102 (2018)

    work page 2018