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
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.
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
- 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
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.
Referee Report
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)
- [§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.
- [§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)
- [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] 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.
- [§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
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
-
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
-
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
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
axioms (1)
- domain assumption Standard DFT, GW, and Bethe-Salpeter approximations accurately describe defect electronic structure and optical excitations in 4H-SiC
Reference graph
Works this paper leans on
- [1]
-
[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
work page 2013
-
[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)
work page 1907
- [4]
-
[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)
work page 2025
-
[6]
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)
work page 2020
-
[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)
work page 2020
- [8]
-
[9]
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)
work page 2024
-
[10]
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)
work page 2025
-
[11]
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)
work page 2009
-
[12]
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)
work page 2013
-
[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)
work page 2016
-
[14]
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)
work page 2022
-
[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)
work page 2021
-
[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)
work page 2020
-
[17]
J. Heyd, G. E. Scuseria, and M. Ernzerhof, Hybrid functionals based on ascreened Coulomb potential, J. Chem. Phys. 118, 8207 (2003)
work page 2003
-
[18]
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)
work page 2025
-
[19]
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
work page 2024
-
[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)
work page 2025
-
[21]
T. Qian and J.-H. Yanga, First-principles studies of nitrogen doping in 4H-SiC. J. Appl. Phys. 137, 215701 (2025)
work page 2025
-
[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)
work page 2025
-
[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)
work page 2011
-
[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)
work page 2011
-
[25]
M. Shishkin, M. Marsman, and G. Kresse, Accurate Quasiparticle Spectra from Self- Consistent GW Calculations with Vertex Corrections. Phys. Rev. Lett. 99, 246403 (2007)
work page 2007
- [26]
- [27]
- [28]
-
[29]
J. P. Perdew, S. Burke, M. Ernzerhof, Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 77, 3865 (1996)
work page 1996
-
[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)
work page 1999
-
[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)
work page 2025
-
[32]
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)
work page 2018
- [33]
-
[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)
work page 2010
-
[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
work page 2020
-
[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)
work page 2024
-
[37]
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)
work page 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.