pith. sign in

arxiv: 2604.15175 · v1 · submitted 2026-04-16 · ❄️ cond-mat.mtrl-sci

Optimal spin-qubit hallmarks of sulfur-vacancy defects in 4H-SiC: Design from first principles

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

classification ❄️ cond-mat.mtrl-sci
keywords 4H-SiCspin qubitsulfur vacancydefect engineeringfirst principlesGW calculationsBethe-Salpeter equationoptical control
0
0 comments X

The pith

A sulfur-vacancy defect in 4H-SiC combines a stable triplet ground state with isolated gap levels and near-infrared optical transitions to enable optically controlled spin qubits.

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

The paper applies first-principles calculations to a defect formed by a silicon vacancy paired with a sulfur atom on a carbon site, called VSiSC, inside 4H-silicon carbide. It reports that every configuration of this defect maintains a spin-triplet ground state that lies below higher-energy singlet states, remains dynamically and thermodynamically stable, and produces sharp isolated electronic peaks inside the host band gap according to GW results. Bethe-Salpeter equation solutions further show intense optical excitations in the near-infrared range, while the constituent elements possess high-abundance zero-nuclear-spin isotopes that support long coherence times. A sympathetic reader would care because these properties together address the combined needs for optical addressability and minimal decoherence in solid-state spin qubits.

Core claim

The authors establish that the VSiSC defect in 4H-SiC has a spin-triplet ground state lower in energy than singlet states for all considered configurations, forms isolated electronic states within the band gap as shown by GW calculations, and exhibits intense near-infrared optical excitations according to Bethe-Salpeter equation results, making it an excellent candidate for an optically controlled spin qubit with high spin-coherence time ensured by zero-nuclear-spin isotopes of the constituent elements.

What carries the argument

The VSiSC defect (silicon vacancy paired with sulfur substituting carbon) whose electronic states and optical transitions are computed using GW and Bethe-Salpeter equation methods to confirm triplet stability and qubit functionality.

If this is right

  • The triplet ground state lower than singlets enables the spin polarization cycle required for qubit initialization and readout.
  • Sharp isolated defect peaks inside the band gap reduce unwanted mixing with host electronic states.
  • Intense near-infrared optical excitations provide efficient optical control of the spin state.
  • High-abundance zero-nuclear-spin isotopes of silicon, carbon, and sulfur support long spin coherence without isotopic enrichment.

Where Pith is reading between the lines

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

  • The same vacancy-dopant pairing strategy could be applied to screen other wide-bandgap hosts for additional qubit candidates.
  • Natural isotopic abundances of the constituent atoms may allow fabrication without special purification steps.
  • The near-infrared optical window aligns with common photonic infrastructure, suggesting easier integration into existing optical networks.

Load-bearing premise

The GW approximation and Bethe-Salpeter equation accurately describe the defect's ground-state spin, excitation energies, and stability in real material without substantial errors from supercell size or functional approximations.

What would settle it

Experimental measurement of the ground-state spin multiplicity and the energies of optical absorption lines in a sample containing the VSiSC defect; absence of the predicted triplet ground state or mismatch with the calculated near-infrared transitions would falsify the qubit suitability claim.

Figures

Figures reproduced from arXiv: 2604.15175 by Marisol Alc\'antara Ortigoza, Sergey Stolbov.

Figure 1
Figure 1. Figure 1: in an A-B-C-B fashion (Fig.1). Each bilayer is made of one flat monolayer of pure C and another flat monolayer of pure Si. Each bilayer is internally such that one species (say, Si) sits directly above the other species (say, C), with respect to the c￾axis, as shown in [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: Spin density iso-surfaces calculated for the triplet and singlet states of the VSiSC defect in the hk and kh configurations. Black, green, and yellow balls mark the Si, C, and S atom, respectively. Red and blue lobes correspond to the spin-up and spin-down density distribution [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
read the original abstract

By applying our methodology, we propose a defect in 4H-SiC which combines a Si vacancy and a C atom substituted with S (VSiSC) to have a spin-triplet ground state with the spin qubit functionality. Our calculations confirm that all configurations of the defect have a dynamically and thermodynamically stable triplet ground state and higher energy singlet states, essential for the spin-qubit polarization cycle. From GW calculations, we found that the electronic states associated with the defect form sharp and isolated peaks within the band gap for both triplet and singlet states. Further Bethe-Salpeter-equation calculations show that all considered configurations have intense optical excitations in the near infrared spectrum range. Analysis of the excitation energies and rates indicate that the VSiSC defect can be an excellent optically controlled spin qubit. Crucially, the host elements and the dopant have high-abundance isotopes with zero nuclear spin ensuring high spin-coherence time of the qubit.

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 / 3 minor

Summary. The manuscript proposes the VSiSC defect (Si vacancy paired with S substituting C) in 4H-SiC as a candidate for an optically controlled spin qubit. First-principles calculations establish that all configurations exhibit a stable triplet ground state with higher-lying singlets, isolated defect levels inside the gap (from GW), and intense near-infrared optical transitions (from Bethe-Salpeter equation) that enable spin initialization and readout. Thermodynamic and dynamic stability are confirmed, and the zero-nuclear-spin isotopes of the constituent elements are highlighted as ensuring long coherence times.

Significance. If the computed defect properties hold, the work identifies a promising qubit platform in a technologically mature wide-bandgap semiconductor. The combination of a triplet ground state, isolated gap levels, strong NIR transitions, and nuclear-spin-free host/dopant isotopes addresses several key requirements for scalable spin qubits. The use of GW and BSE methods for defect electronic structure and optical spectra provides a more rigorous treatment than standard DFT, strengthening the design claim.

major comments (2)
  1. [§3.2] §3.2 and Fig. 4: The claim that the optical transitions are 'intense' and suitable for qubit control rests on BSE oscillator strengths, but no quantitative comparison is provided to the known VSi or other SiC defects (e.g., divacancy) whose experimental transition rates are available. Without this benchmark, it is difficult to judge whether the computed rates are competitive.
  2. [§4.1] §4.1, Table II: Supercell-size convergence for the GW quasiparticle levels is shown only for the 4×4×1 cell; the shift between 3×3×1 and 4×4×1 is ~0.2 eV for the triplet defect state. This residual finite-size error could affect the predicted isolation of the levels inside the gap and should be extrapolated or corrected before asserting 'sharp and isolated peaks'.
minor comments (3)
  1. [Abstract] The abstract states that 'all configurations' are stable, but the main text discusses only three; a brief enumeration of the full set of configurations considered (and why others were discarded) would improve clarity.
  2. [§2] Notation for the defect (VSiSC) is introduced without an explicit structural diagram in the main text; adding a ball-and-stick figure of the relaxed triplet geometry would aid readers unfamiliar with SiC defect literature.
  3. [§5] The zero-nuclear-spin argument is stated qualitatively; a short table listing natural abundances and nuclear spins of 28Si, 12C, and 32S would make the coherence-time advantage quantitative.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of our work and the recommendation for minor revision. We address each major comment below and will update the manuscript accordingly.

read point-by-point responses
  1. Referee: [§3.2] §3.2 and Fig. 4: The claim that the optical transitions are 'intense' and suitable for qubit control rests on BSE oscillator strengths, but no quantitative comparison is provided to the known VSi or other SiC defects (e.g., divacancy) whose experimental transition rates are available. Without this benchmark, it is difficult to judge whether the computed rates are competitive.

    Authors: We agree that a direct quantitative benchmark would strengthen the assessment of competitiveness. In the revised manuscript we will add a comparison of our computed BSE oscillator strengths and implied transition rates to available experimental and theoretical values for the VSi defect and the divacancy in 4H-SiC, placed in §3.2 near the discussion of Fig. 4. revision: yes

  2. Referee: [§4.1] §4.1, Table II: Supercell-size convergence for the GW quasiparticle levels is shown only for the 4×4×1 cell; the shift between 3×3×1 and 4×4×1 is ~0.2 eV for the triplet defect state. This residual finite-size error could affect the predicted isolation of the levels inside the gap and should be extrapolated or corrected before asserting 'sharp and isolated peaks'.

    Authors: We thank the referee for highlighting this point. Although the shift between the two supercells is modest relative to the gap width, we will perform an extrapolation of the GW levels to the infinite-supercell limit (or compute one additional larger cell) and include the corrected values in the revised Table II and §4.1 to rigorously confirm that the defect states remain sharply isolated. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The paper's central claims rest on standard first-principles workflows (DFT ground-state relaxations, GW quasiparticle corrections, and Bethe-Salpeter equation optical spectra) applied to candidate defect configurations in 4H-SiC. These methods are benchmarked against established SiC band gaps and known defect levels rather than fitted to the target qubit metrics (triplet ground state, NIR transitions, coherence). No equation reduces the predicted spin-qubit hallmarks to a parameter fit or self-citation chain; the stability, level positions, and transition intensities are computed outputs, not inputs. The derivation is therefore self-contained against external benchmarks and receives the default non-circularity finding.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work relies on standard first-principles approximations for defect calculations in semiconductors; no additional free parameters or new physical entities are introduced beyond the defect model itself.

axioms (1)
  • domain assumption Standard DFT, GW, and Bethe-Salpeter approximations accurately describe defect electronic structure and optical excitations in 4H-SiC
    Invoked to predict ground-state spin, stability, and excitations from the abstract methodology description.

pith-pipeline@v0.9.0 · 5474 in / 1174 out tokens · 46843 ms · 2026-05-10T10:22:23.990975+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

37 extracted references · 37 canonical work pages

  1. [1]

    Gruber, A

    A. Gruber, A. Drabenstedt, C. Tietz, L. Fleury, J. Wrachtrup, C. von Borczyskowski, Scanning Confocal Optical Microscopy and Magnetic Resonance on Single Defect Centers, Science 276 2012 (1997)

  2. [2]

    M. W. Doherty, N. B. Manson, P. Delaney, F. Jelezko, J. Wrachtrup, and L. C.L. Hollenberg, The nitrogen-vacancy colour centre in diamond. Physics Reports 528, 1 (2013) and references therein

  3. [3]

    Gali, Ab initio theory of the nitrogen-vacancy center in diamond, Nanophotonics; 8, 1907–1943 (2019)

    A. Gali, Ab initio theory of the nitrogen-vacancy center in diamond, Nanophotonics; 8, 1907–1943 (2019)

  4. [4]

    Koehl, B

    W. Koehl, B. Buckley, F. Heremans, G. Calusine, and D. D. Awschalom, Room temperature coherent control of defect spin qubits in silicon carbide. Nature 479, 84–87 (2011)

  5. [5]

    G. Bian, G. Thiering, and A. Gali, Theory of optical spin-polarization of axial divacancy and nitrogen-vacancy defects in 4H-SiC. Phys. Rev. Research 7, 013320 (2025)

  6. [6]

    Gottscholl, M

    A. Gottscholl, M. Kianinia, V. Soltamov, S. Orlinskii, G. Mamin, C. Bradac, C. Kasper, K. Krambrock, A. Sperlich, M. Toth, I. Aharonovich, and V. Dyakonov, Initialization 16 and read-out of intrinsic spin defects in a van der Waals crystal at room temperature. Nat. Mat. 19, 540 (2020)

  7. [7]

    J. R. Reimers, J. Shen, M. Kianinia, C. Bradac, I. Aharonovich, M. J. Ford, and P. Piecuch, Photoluminescence, photophysics, and photochemistry of the 𝑉𝐵 − − defect in hexagonal boron nitride. Phys. Rev. B 102, 144105 (2020)

  8. [8]

    Mathur, A

    N. Mathur, A. Mukherjee, X. Gao, J. Luo, B. A. McCullian, T. Li, A. N. Vamivakas, and G. D. Fuchs, Excited-state spin-resonance spectroscopy of 𝑉𝐵 − defect centers in hexagonal boron nitride. Nat. Comm. 13, 3233 (2022)

  9. [9]

    Stolbov and M

    S. Stolbov and M. Alcántara Ortigoza, Aluminum Vacancy/Sulfur Complex in Wurtzite AlN as an Optically Controllable Spin Qubit. Phys. Rev. B 109, L241108 (2024)

  10. [10]

    Stolbov and M

    S. Stolbov and M. Alcántara Ortigoza, Spin qubit properties of the boron-vacancy/carbon defect in the two-dimensional hexagonal boron nitride. J. Phys.: Condens. Matter 37 385503 (2025)

  11. [11]

    Balasubramanian, P

    G. Balasubramanian, P. Neumann, D. Twitchen, M. Markham, R. Kolesov, N. Mizuochi, J. Isoya, J. Achard, J. Beck, J. Tissler, V. Jacques, P. R. Hemmer, F. Jelezko, and J. Wrachtrup, Ultralong spin coherence time in isotopically engineered diamond. Nat. Mat. 8, 383 (2009)

  12. [12]

    Bar-Gill, L.M

    N. Bar-Gill, L.M. Pham, A. Jarmola, D. Budker, and R.L. Walsworth, Solid-state electronic spin coherence time approaching one second. Nat. Communs. 4, 1743 (2013)

  13. [13]

    H. Seo, A. L. Falk, P. V. Klimov, K. C. Miao, G. Galli, and D. D. Awschalom. Quantum decoherence dynamics of divacancy spins in silicon carbide. Nat. Communs. 7, 12935 (2016)

  14. [14]

    and Ohshima, T., 2022

    Anderson, C.P., Glen, E.O., Zeledon, C., Bourassa, A., Jin, Y., Zhu, Y., Vorwerk, C., Crook, A.L., Abe, H., Ul-Hassan, J. and Ohshima, T., 2022. Five-second coherence of a single spin with single-shot readout in silicon carbide. Science advances, 8(5), 5912 (2022)

  15. [15]

    V. H. Rodgers, L. B. Hughes, M. Xie, P. C. Maurer, S. Kolkowitz, A. C. Bleszynski Jayich, and N. P. de Leon. Materials challenges for quantum technologies based on color centers in diamond. MRS Bulletin 46, 1 (2021)

  16. [16]

    N. T. Son, C. P. Anderson, A. Bourassa, K. C. Miao, C. Babin C, M. Widmann, M. Niethammer, J. U. Hassan, N. Morioka, I. G. Ivanov, F. Kaiser, J. Wrachtrup, and D. D. Awschalom, Developing silicon carbide for quantum spintronics. Appl. Phys. Lett. 116, 190501 (2020)

  17. [17]

    J. Heyd, G. E. Scuseria, and M. Ernzerhof, Hybrid functionals based on ascreened Coulomb potential, J. Chem. Phys. 118, 8207 (2003)

  18. [18]

    Qian, X.-G

    T. Qian, X.-G. Gong, and J.-H. Yang, First-principles studies of hydrogen irradiation effects on the photoluminescence properties of nitrogen-vacancy centers in 4H-SiC. Phys. Rev. Mat. 9, 116201 (2025)

  19. [19]

    Iwamoto, T

    S. Iwamoto, T. Shimura, H. Watanabe, and T. Kobayashi, Oxygen-related defects in 4H- SiC from first principles. Appl. Phys. Express 17, 051008 (2024). 17

  20. [20]

    S. G. Carter, I. Tathfif, C. Burgess, B. VanMil, S. Debata, and P. Dev, Influence of nitrogen doping and annealing on the silicon vacancy in 4H-SiC. Phys. Rev. B 112, 085209 (2025)

  21. [21]

    Qian and J.-H

    T. Qian and J.-H. Yanga, First-principles studies of nitrogen doping in 4H-SiC. J. Appl. Phys. 137, 215701 (2025)

  22. [22]

    L. Sun, W. Zhang, B. Peng, and W. Zhang, Investigation of PCVSi- defect in 4H–SiC as a candidate for a qubit. J. Phys. D: Appl. Phys. 58 025306, (2025)

  23. [23]

    J. R. Maze, A. Gali, E. Togan, Y. Chu, A. Trifonov, E. Kaxiras. and M D Lukin, Properties of nitrogen-vacancy centers in diamond: the group theoretic approach. New J. Phys. 13 025025 (2011)

  24. [24]

    M. W. Doherty, N. B. Manson, P. Delaney, and L. C. L. Hollenberg, The negatively charged nitrogen-vacancy centre in diamond: the electronic solution. New Journal of Physics 13 025019 (2011)

  25. [25]

    Shishkin, M

    M. Shishkin, M. Marsman, and G. Kresse, Accurate Quasiparticle Spectra from Self- Consistent GW Calculations with Vertex Corrections. Phys. Rev. Lett. 99, 246403 (2007)

  26. [26]

    Onida, L

    G. Onida, L. Reining, and A. Rubio, Electronic excitations: density-functional versus many-body Green’s-function approaches. Rev. Mod. Phys. 74, 601 (2002)

  27. [27]

    Kresse, G

    G. Kresse, G. J. Furthmüller, Efficient Iterative Schemes for Ab Initio Total-Energy Calculations Using a Plane-Wave Basis Set. Comput. Mat. Sci., 6, 15 (1996)

  28. [28]

    Kresse, J

    G. Kresse, J. Joubert, From Ultrasoft Pseudopotentials to the Projector Augmented-Wave Method. Phys. Rev. B, 59, 1758 (1999)

  29. [29]

    J. P. Perdew, S. Burke, M. Ernzerhof, Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 77, 3865 (1996)

  30. [30]

    Kokalj, XCrySDen—a new program for displaying crystalline structures and electron densities, J

    A. Kokalj, XCrySDen—a new program for displaying crystalline structures and electron densities, J. Mol. Graphics Modelling, 17, 176 (1999)

  31. [31]

    B. R. McClain S. Stolbov, and M. Alcántara Ortigoza, Ab initio evaluation of the electronic and optical properties of VBCB defect in wurtzite boron nitride as promising single-photon emitter. Phys. Letts. A 564, 131085 (2025)

  32. [32]

    Davidsson, V

    J. Davidsson, V. Ivády, R. Armiento, N. T. Son, A. Gali, and I. A. Abrikosov, First principles predictions of magneto-optical data for semiconductor point defect identification: the case of divacancy defects in 4H–SiC. New J. Phys. 20, 023035 (2018)

  33. [33]

    Ivády, I

    V. Ivády, I. A. Abrikosov, and A. Gali, First principles calculation of spin-related quantities for point defect qubit research. npj Computational Materials 4:76 (2018)

  34. [34]

    J. R. Weber, W. F. Koehl, J. B. Varley, A. Janotti, B. B. Buckley, C. G. Van de Walle, and D. D. Awschalom, Quantum computing with defects. PNAS, 107 8513 (2010)

  35. [35]

    Y.-Q. Su, L. Zhang, Y . Wang, J.-X. Liu, V . Muravev, K. Alexopoulos, I. A. W. Filot, D. G. Vlachos, and E. J. M. Hensen, Stability of heterogeneous single-atom catalysts: a scaling law mapping thermodynamics to kinetics. npj Comput. Mater. 6, 144 (2020). 18

  36. [36]

    M. K. Mainali, P. Dulal, B. Shrestha, E. Amonette, A. Shan, and N. J. Podraza, Optical properties of 4H-SiC and 6H-SiC from infrared to vacuum ultraviolet spectral range ellipsometry (0.05–8.5eV). Surf. Sci. Spectra 31, 026003 (2024)

  37. [37]

    Bockstedte, F

    M. Bockstedte, F. Schütz, T. Garratt, V . Ivády, and A. Gali, Ab initio description of highly correlated states in defects for realizing quantum bits. npj Quantum Materials 3, 31 (2018)