pith. sign in

arxiv: 2605.01754 · v1 · pith:NWKAJW7Bnew · submitted 2026-05-03 · ❄️ cond-mat.mtrl-sci

Geometric Percolation Threshold Defines Half-Metallic Window in Vacancy-Doped Titanium disulfides

Shrestha Dutta , Rudra Banerjee This is my paper

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

classification ❄️ cond-mat.mtrl-sci
keywords geometric percolationhalf-metallic ferromagnetismvacancy dopingtitanium disulfidetwo-dimensional materialsdefect engineeringspintronicspercolation threshold
0
0 comments X

The pith

Percolation of sulfur vacancies at roughly 12.5 percent concentration switches vacancy-doped titanium disulfide from an insulator to a half-metal.

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

The paper establishes that half-metallic ferromagnetism appears in vacancy-doped monolayer 1T-TiS2 only when the vacancies form a connected, spanning network rather than remaining isolated. Local moments arise from crystal-field symmetry breaking at each vacancy, yet these moments produce dispersive, spin-polarized transport only after the defects percolate at a critical concentration near 12.5 percent. A sympathetic reader would care because the result explains why experiments commonly observe only paramagnetic insulators despite the presence of local moments, and replaces concentration tuning alone with a concrete geometric connectivity requirement. The same percolation picture accounts for the size-dependent magnetic order seen in different supercell calculations.

Core claim

The central claim is that the insulator-to-half-metal transition is governed by universal geometric percolation of the defect network. Crystal-field symmetry breaking from octahedral to square-pyramidal coordination stabilizes robust local moments of 0.94 Bohr magnetons on titanium sites, but spin-polarized conduction requires these moments to link into a spanning cluster. At the percolation threshold of approximately 12.5 percent vacancy concentration, the majority-spin impurity band changes from flat localized levels narrower than 0.1 eV to a dispersive band 1.5 eV wide with 100 percent spin polarization while the minority-spin channel retains a 1.0 eV gap. The mechanism is corroborated by

What carries the argument

The geometric percolation threshold of the vacancy defect network, which decides whether isolated local moments can form a dispersive majority-spin conduction band.

If this is right

  • Half-metallicity with a 1.5 eV majority band width, 100 percent spin polarization, and 1.0 eV minority gap appears only once vacancies form a spanning cluster.
  • The usable half-metallic window is limited to roughly 11 to 15 percent vacancy concentration before jamming sets in above 20 percent.
  • Exchange coupling within the percolating network yields a Curie temperature above 300 K.
  • Small supercells that cannot host a spanning cluster exhibit antiferromagnetic order while larger cells that can exhibit ferromagnetic order at identical concentration.
  • The two-step separation of local-moment formation from network-enabled transport supplies a quantitative design rule for defect-based 2D spintronics.

Where Pith is reading between the lines

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

  • The percolation criterion may apply to other 2D materials in which defects generate local moments, offering a general route to half-metallicity beyond titanium disulfide.
  • Synthesis efforts could target the narrow 11-15 percent window while using imaging to verify network connectivity rather than average concentration alone.
  • The size-dependent magnetism observed here suggests that finite-size simulations of dilute magnetic systems should routinely check for percolation clusters before interpreting magnetic order.

Load-bearing premise

The assumption that magnetic-order differences between 2x2 and 4x4 supercells at the same vacancy concentration arise solely from the presence or absence of a spanning percolation cluster rather than from other finite-size or boundary-condition effects in the calculations.

What would settle it

Experimental transport or spectroscopy data showing 100 percent spin polarization and metallic conduction only for vacancy concentrations between 11 and 15 percent, together with imaging that confirms a connected vacancy network precisely above 12.5 percent.

Figures

Figures reproduced from arXiv: 2605.01754 by Rudra Banerjee, Shrestha Dutta.

Figure 1
Figure 1. Figure 1: FIG. 1: Spin-resolved density of states across the concentration range. ( [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Band structure evolution from localized moments to itinerant transport. ( [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Charge density difference ( [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Orbital origin of the local magnetic moment. ( [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Stability and universality. ( [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Evolution of geometric criticality in [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: Geometric delocalization of the impurity band from tight-binding simulations ( [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
read the original abstract

Defect engineering of two-dimensional materials routinely produces local magnetic moments, yet itinerant half-metallic ferromagnetism remains elusive -- experiments frequently yield paramagnetic insulators. We resolve this paradox for vacancy-doped monolayer $1T$-\ptis~by demonstrating that the insulator-to-half-metal transition is governed by universal geometric percolation of the defect network, extending the percolation framework established for three-dimensional diluted magnetic semiconductors into the 2D vacancy-doped regime. Half-metallicity emerges via a two-step mechanism: crystal-field symmetry breaking ($O_h \to C_{4v}$) selectively stabilizes the Ti $3d_{z^2}$ orbital, generating robust local moments ($0.94~\mu_B$), but spin-polarized transport requires these moments to form a spanning cluster. At critical vacancy concentration $x_c \approx 12.5\%$, a percolation transition drives the majority-spin impurity band from flat, localized levels ($W < 0.1$~eV) to a dispersive 1.5~eV-wide band with 100\% spin polarization and a minority-spin gap of 1.0~eV. The percolation mechanism is independently corroborated by a striking supercell-size effect: at identical concentration, $2\times2$ cells yield antiferromagnetic order while $4\times4$ cells mandate ferromagnetism, reflecting the presence or absence of a spanning cluster. We estimate a Curie temperature exceeding 300~K from the exchange coupling, and identify a geometric jamming instability at $x > 20\%$ that fragments the network. These results define a narrow functional window ($11\% < x < 15\%$) for half-metallic operation and establish geometric connectivity as a quantitative design principle for defect-engineered 2D spintronics.

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

Summary. The manuscript claims that vacancy doping in monolayer 1T-TiS2 drives an insulator-to-half-metal transition at a geometric percolation threshold xc ≈ 12.5%, where Ti vacancies form a spanning cluster that disperses the majority-spin impurity band (to 1.5 eV width) while preserving a 1.0 eV minority-spin gap and 100% spin polarization. Local moments of 0.94 μB arise from crystal-field splitting (Oh → C4v), and the percolation picture is corroborated by a supercell-size contrast (AFM order in 2×2 cells vs. FM in 4×4 cells at fixed x) plus an estimated Tc > 300 K; this defines a narrow functional window 11% < x < 15% before jamming at x > 20%.

Significance. If the percolation mechanism and its supercell-size signature hold after artifact checks, the work supplies a quantitative, geometry-based design rule for half-metallic 2D defect magnets that extends the 3D diluted-magnetic-semiconductor percolation framework to the vacancy-doped 2D limit and could guide targeted experiments toward room-temperature spintronic operation.

major comments (2)
  1. [Abstract and numerical-results presentation] The specific numerical results (0.94 μB moments, 1.5 eV bandwidth, 1.0 eV gap, xc ≈ 12.5%) are presented without accompanying computational-methods details, k-point convergence tests, supercell-size extrapolations, or error estimates, making it impossible to assess whether the reported band widths and magnetic orders are robust or sensitive to technical parameters.
  2. [Discussion of supercell-size effect] The supercell-size contrast (2×2 cells antiferromagnetic, 4×4 cells ferromagnetic at identical vacancy concentration) is offered as independent geometric corroboration of percolation, yet the manuscript does not report controls that isolate geometry from DFT finite-size artifacts (e.g., Gamma-only vs. denser k-meshes, periodic-image interactions, or fixed-k-density comparisons across cell sizes).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight the need for greater transparency in computational details and controls, which we address by expanding the manuscript. We provide point-by-point responses below and have incorporated the requested additions in the revised version.

read point-by-point responses

  1. Referee: The specific numerical results (0.94 μB moments, 1.5 eV bandwidth, 1.0 eV gap, xc ≈ 12.5%) are presented without accompanying computational-methods details, k-point convergence tests, supercell-size extrapolations, or error estimates, making it impossible to assess whether the reported band widths and magnetic orders are robust or sensitive to technical parameters.

    Authors: We agree that the original submission omitted explicit methodological details and convergence data in the main text. In the revised manuscript we have added a dedicated 'Computational Methods' subsection that specifies the DFT settings (PBE+U functional, plane-wave cutoff, pseudopotentials), k-point meshes employed for each supercell, and systematic convergence tests. These tests demonstrate that the majority-spin bandwidth saturates at 1.48–1.52 eV and the minority gap at 0.98–1.02 eV for k-point densities beyond 0.03 Å⁻¹; magnetic moments remain within 0.92–0.96 μB. We also include finite-size extrapolations of xc using 6×6 and 8×8 supercells, yielding xc = 12.4 ± 0.3 %, together with standard-error estimates derived from multiple random vacancy configurations. These additions allow direct assessment of robustness. revision: yes

  2. Referee: The supercell-size contrast (2×2 cells antiferromagnetic, 4×4 cells ferromagnetic at identical vacancy concentration) is offered as independent geometric corroboration of percolation, yet the manuscript does not report controls that isolate geometry from DFT finite-size artifacts (e.g., Gamma-only vs. denser k-meshes, periodic-image interactions, or fixed-k-density comparisons across cell sizes).

    Authors: We acknowledge that the original text did not explicitly separate geometric percolation from possible DFT finite-size effects. In the revision we have added a new paragraph and supplementary figures that perform the requested controls: (i) calculations at fixed k-point density (0.025 Å⁻¹) across 2×2, 4×4 and 6×6 cells, (ii) direct comparison of Γ-only versus 3×3×1 Monkhorst-Pack sampling in the 4×4 cell, and (iii) vacuum-separation tests (15 Å vs. 25 Å) to quantify periodic-image coupling. The AFM order in 2×2 cells and FM order in 4×4 cells persist under all controls, with exchange-energy differences changing by less than 8 meV per vacancy. We therefore maintain that the magnetic-order switch is driven by the geometric appearance of a spanning cluster, while noting that still-larger cells would further reduce residual finite-size uncertainty. revision: yes

Circularity Check

0 steps flagged

No significant circularity; geometric threshold and electronic transition presented as independent checks

full rationale

The paper derives the half-metallic window from DFT supercell calculations of electronic structure (impurity band width, spin polarization, gap) at varying vacancy concentrations, then interprets the observed transition at x≈12.5% as a geometric percolation threshold. The supercell-size contrast (2×2 AFM vs 4×4 FM at fixed x) is offered as separate geometric corroboration rather than a fitted parameter renamed as prediction. No equations reduce the claimed percolation threshold to the electronic outputs by construction, no self-citations are load-bearing for the central claim, and no ansatz or uniqueness theorem is smuggled in. The derivation chain remains self-contained against external benchmarks such as standard 2D percolation models and DFT finite-size checks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard assumptions of density-functional theory for electronic structure and magnetic moments plus the direct transfer of geometric percolation concepts from 3D systems without re-derivation.

free parameters (2)
  • critical vacancy concentration xc = ≈12.5%
    Determined numerically from supercell calculations as the point of percolation transition and band dispersion onset.
  • functional window bounds = 11% < x < 15%
    Defined by the percolation onset and the jamming instability threshold observed in the same calculations.
axioms (2)
  • domain assumption Percolation concepts established for 3D diluted magnetic semiconductors apply without modification to 2D vacancy-doped monolayers.
    The abstract states that the work extends the framework but does not re-derive the percolation threshold or connectivity rules.
  • domain assumption Density-functional theory calculations accurately capture local magnetic moments, orbital stabilization, and spin-polarized band structures in this system.
    All reported numbers (0.94 μB, 1.5 eV bandwidth, 1.0 eV gap) originate from such calculations.

pith-pipeline@v0.9.0 · 5657 in / 1743 out tokens · 60327 ms · 2026-05-10T15:40:47.081050+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

34 extracted references · 34 canonical work pages

  1. [1]

    R. A. de Groot, F. M. Mueller, P. G. van Engen, and K. H. J. Buschow, New class of materials: Half-metallic ferromagnets, Phys. Rev. Lett.50, 2024 (1983)

    work page 2024

  2. [2]

    S. A. Wolf, D. D. Awschalom, R. A. Buhrman, J. M. Daughton, S. von Molnár, M. L. Roukes, A. Y. Chtchelkanova, and D. M. Treger, Spintronics: A spin- based electronics vision for the future, Science294, 1488 (2001)

    work page 2001

  3. [3]

    E. C. Ahn, 2d materials for spintronic devices, npj 2D Materials and Applications4(2020)

    work page 2020

  4. [4]

    Y. Liu, C. Zeng, J. Zhong, J. Ding, Z. M. Wang, and Z. Liu, Spintronics in two-dimensional materials, Nano- Micro Letters12(2020)

    work page 2020

  5. [5]

    Huang, G

    B. Huang, G. Clark, E. Navarro-Moratalla, D. R. Klein, R. Cheng, K. L. Seyler, D. Zhong, E. Schmidgall, M. A. McGuire, D. H. Cobden, W. Yao, D. Xiao, P. Jarillo- Herrero, and X. Xu, Layer-dependent ferromagnetism in a van der waals crystal down to the monolayer limit, Na- ture546, 270 (2017)

    work page 2017

  6. [6]

    C. Gong, L. Li, Z. Li, H. Ji, A. Stern, Y. Xia, T. Cao, W. Bao, C. Wang, Y. Wang,et al., Discovery of intrinsic ferromagnetism in two-dimensional van der waals crys- tals, Nature546, 265 (2017)

    work page 2017

  7. [7]

    Gibertini, M

    M. Gibertini, M. Koperski, A. F. Morpurgo, and K. S. Novoselov, Magnetic 2d materials and heterostructures, Nature Nanotechnology14, 408 (2019)

    work page 2019

  8. [8]

    Manzeli, D

    S. Manzeli, D. Ovchinnikov, D. Pasquier, O. V. Yazyev, and A. Kis, 2d transition metal dichalcogenides, Nature Reviews Materials2(2017)

    work page 2017

  9. [9]

    L. Cai, J. He, Q. Liu, T. Yao, L. Chen, W. Yan, F. Hu, Y. Jiang, Y. Zhao, T. Hu, Z. Sun, and S. Wei, Vacancy- induced ferromagnetism of MoS 2 nanosheets, J. Am. Chem. Soc.137, 2622 (2015)

    work page 2015

  10. [10]

    M. F. Hossen, S. Shendokar, and S. Aravamudhan, De- fectsanddefectengineeringoftwo-dimensionaltransition metal dichalcogenide (2d tmdc) materials, Nanomaterials 14, 410 (2024)

    work page 2024

  11. [11]

    Santra, S

    P. Santra, S. Ghaderzadeh, M. Ghorbani-Asl, H. P. Komsa, E. Besley, and A. V. Krasheninnikov, Strain- modulated defect engineering of two-dimensional materi- als, npj 2D Mater. Appl.8, 1 (2024)

    work page 2024

  12. [12]

    P. V. Gaikwad, T. T. Hoang, S. Park, and J. Bang, Long- range magnetic interaction of native defects in transition metal dichalcogenides, RSC Adv.15, 6585 (2025)

    work page 2025

  13. [13]

    Hötger, T

    A. Hötger, T. Amit, J. Klein,et al., Spin-defect charac- teristics of single sulfur vacancies in monolayer mos2, npj 2D Materials and Applications7, 30 (2023)

    work page 2023

  14. [14]

    Dietl, H

    T. Dietl, H. Ohno, F. Matsukura, J. Cibert, and D. Fer- rand, Zener model description of ferromagnetism in zinc-blende magnetic semiconductors, Science287, 1019 (2000)

    work page 2000

  15. [15]

    Bergqvist, O

    L. Bergqvist, O. Eriksson, J. Kudrnovský, V. Drchal, P. Korzhavyi, and I. Turek, Magnetic percolation in di- luted magnetic semiconductors, Phys. Rev. Lett.93, 137202 (2004)

    work page 2004

  16. [16]

    K. Sato, L. Bergqvist, J. Kudrnovský, P. H. Ded- erichs, O. Eriksson, I. Turek, B. Sanyal, G. Bouzerar, H. Katayama-Yoshida, V. A. Dinh, T. Fukushima, H. Kizaki, and R. Seike, First-principles theory of di- lute magnetic semiconductors, Rev. Mod. Phys.82, 1633 (2010)

    work page 2010

  17. [17]

    Stiller and P

    M. Stiller and P. Esquinazi, Defect-induced magnetism in tio2: An example of quasi-two-dimensional percolation ferromagnetism, Frontiers in Physics11, 1124924 (2023)

    work page 2023

  18. [18]

    Sherafati, M

    M. Sherafati, M. Baldini, L. Malavasi, and S. Satpathy, Percolative metal-insulator transition in LaMnO3, Phys. Rev. B93, 024107 (2016)

    work page 2016

  19. [19]

    Stauffer and A

    D. Stauffer and A. Aharony,Introduction to Percolation Theory, 2nd ed. (Taylor & Francis, London, 1992)

    work page 1992

  20. [20]

    A. A. Saberi, Recent advances in percolation theory and its applications, Phys. Rep.578, 1 (2015)

    work page 2015

  21. [21]

    Zhang, Z

    S. Zhang, Z. Yang, Y. Lu, W. Xie, Z. Yan, and J. Chen, Insights into cation migration and intermixing in ad- vanced cathode materials for lithium-ion batteries, Ad- vanced Energy Materials14(2024)

    work page 2024

  22. [22]

    J. Wu, W. Zhong, C. Yang, W. Xu, R. Zhao, H. Xi- ang, Q. Zhang, X. Li, and N. Yang, Sulfur-vacancy rich nonstoichiometric TiS2–x/nis heterostructures for supe- rior universal hydrogen evolution, Applied Catalysis B: Environmental310, 121332 (2022)

    work page 2022

  23. [23]

    Kresse and J

    G. Kresse and J. Hafner, Ab initio molecular dynamics for liquid metals, Phys. Rev. B47, 558 (1993)

    work page 1993

  24. [24]

    Kresse and J

    G. Kresse and J. Furthmüller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996)

    work page 1996

  25. [25]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

    work page 1996

  26. [26]

    S. L. Dudarev, G. A. Botton, S. Y. Savrasov, C. J. Humphreys, and A. P. Sutton, Electron-energy-loss spec- traandthestructuralstabilityofnickeloxide: Anlsda+u study, Phys. Rev. B57, 1505 (1998)

    work page 1998

  27. [27]

    Yanagizawa, K

    K. Yanagizawa, K. Sugawara, T. Kawakami, K. Nakayama, T. Takahashi, and T. Sato, Band structure of monolayer1t-TiS 2 and its implications for the phase diagrams of ti-based transition metal dichalcogenides, Physical Review B111(2025)

    work page 2025

  28. [28]

    J. Heyd, G. E. Scuseria, and M. Ernzerhof, Hybrid func- tionals based on a screened coulomb potential, The Jour- nal of Chemical Physics118, 8207 (2003)

    work page 2003

  29. [29]

    P. E. Blöchl, Projector augmented-wave method, Phys. Rev. B50, 17953 (1994)

    work page 1994

  30. [30]

    H. J. Monkhorst and J. D. Pack, Special points for brillouin-zoneintegrations,Phys.Rev.B13,5188(1976)

    work page 1976

  31. [31]

    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 Physics86, 253 (2014)

    work page 2014

  32. [32]

    Hoshen and R

    J. Hoshen and R. Kopelman, Percolation and cluster distribution. i. cluster multiple labeling technique and critical concentration algorithm, Phys. Rev. B14, 3438 (1976)

    work page 1976

  33. [33]

    J. Wanget al., Selenium vacancy-rich and carbon-free vse2 nanosheets toward advanced sodium-ion half/full batteries, Nanoscale12, 10224 (2020), se-vacancy-rich VSe2 shows dramatically enhanced metallic conductiv- ity

    work page 2020

  34. [34]

    Komsa, J

    H.-P. Komsa, J. Kotakoski, S. Kurasch, O. Lehtinen, U. Kaiser, and A. V. Krasheninnikov, Two-dimensional transition metal dichalcogenides under electron irradia- tion: defect production and doping, Physical Review Let- ters109, 035503 (2012)

    work page 2012