Pith. sign in

REVIEW 1 major objections 2 minor 44 references

Arsenic substitution for molybdenum in MoS2 monolayers shifts the Fermi level to produce p-type behavior while interstitial arsenic produces n-type behavior.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.3

2026-07-01 06:40 UTC pith:FQCNZAHB

load-bearing objection Standard DFT defect run on As in MoS2 that infers p-type and n-type from neutral-supercell Fermi shifts without the formation-energy curves needed to confirm equilibrium doping. the 1 major comments →

arxiv 2606.22119 v2 pith:FQCNZAHB submitted 2026-06-20 cond-mat.mtrl-sci

First-principles study of the impact of As doping on the structural and electronic properties of MoS₂ monolayer

classification cond-mat.mtrl-sci
keywords MoS2 monolayerarsenic dopingDFT calculationsdefect statesp-type dopingn-type dopingFermi level shiftelectronic properties
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper performs density functional theory calculations on a molybdenum disulfide monolayer to determine how arsenic atoms and vacancies change the atomic arrangement and the positions of energy bands. Defect states appear in the middle of the original band gap for all modified structures. When arsenic replaces a molybdenum or sulfur atom the Fermi level moves down into the valence band region, which the calculations identify as p-type doping. When arsenic sits in an interstitial position the Fermi level moves up into the conduction band region, identified as n-type doping. These shifts are presented as routes to improved performance in photocatalysis, photovoltaics, and field-effect transistors.

Core claim

Introduction of arsenic defects in MoS2 monolayers generates midgap states; substitution at the Mo site or S site moves the Fermi level toward the valence band (p-type), while interstitial placement moves it toward the conduction band (n-type), with the substitution cases proposed for photocatalysis and photovoltaics and the interstitial case proposed for field-effect transistors.

What carries the argument

Density functional theory calculations of formation energies, relaxed geometries, and electronic density of states for vacancy and arsenic-doped configurations in the MoS2 monolayer.

Load-bearing premise

The electronic band positions and Fermi level locations computed for ideal defect structures in an isolated monolayer will match the carrier type and device behavior that appear in real samples after growth and contact formation.

What would settle it

Hall-effect or Seebeck-coefficient measurements on experimentally fabricated As-doped MoS2 monolayers that show the opposite carrier sign from the DFT-predicted p-type or n-type behavior.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Substitution of As for Mo creates p-type doping that can supply holes for photocatalytic reactions.
  • The same substitution is claimed to support high-efficiency photovoltaic devices through modified band alignment.
  • Interstitial As creates n-type doping that can increase electron density in field-effect transistor channels.
  • Midgap defect states are expected to influence recombination rates and therefore carrier lifetimes in all doped cases.

Where Pith is reading between the lines

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

  • The same defect-induced carrier-type switch might be testable in other monolayer transition-metal dichalcogenides by replacing the chalcogen or metal site with a group-V element.
  • Substrate interactions or finite-temperature effects omitted from the monolayer model could alter the predicted Fermi-level positions in actual devices.
  • Direct comparison of calculated defect formation energies with measured incorporation rates during growth would test whether the lowest-energy configurations are actually realized.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

1 major / 2 minor

Summary. The manuscript reports DFT calculations on the structural and electronic properties of MoS₂ monolayers containing S or Mo vacancies and As dopants (substitutional at Mo or S sites, or interstitial). Defect states appear in the mid-gap for all cases. The Fermi level shifts downward into the valence band for V_S, V_Mo, As_Mo and As_S (interpreted as p-type) and upward into the conduction band for As interstitial (interpreted as n-type). These shifts are used to propose applications in photocatalysis, high-efficiency photovoltaics, and enhanced FET performance.

Significance. If the equilibrium carrier types were correctly identified, the results would add to the literature on defect engineering of 2D TMDs for optoelectronics and electronics. The computational framework is conventional, but the doping-type assignments rest on an incomplete analysis that does not establish the thermodynamic equilibrium Fermi level.

major comments (1)
  1. [Abstract / Electronic properties results] The classification of p-type (V_S, V_Mo, As_Mo, As_S) and n-type (As interstitial) behavior is based solely on the Fermi-level position in the neutral-supercell DOS. Standard defect physics requires formation-energy curves E_f(q, E_F) for multiple charge states q, solution of the charge-neutrality condition, and 2D image-charge corrections to locate the equilibrium E_F. The abstract and reported results give no indication that these steps were performed; the neutral-DOS shift alone does not determine carrier type when defect levels are deep or when formation energies cross inside the gap. This directly undermines the application claims in the final paragraph.
minor comments (2)
  1. [Abstract] Computational details (XC functional, plane-wave cutoff, k-mesh, supercell size, convergence criteria) are not summarized even at the level of the abstract; these must be stated explicitly for reproducibility.
  2. [Results] No error bars, convergence tests with respect to supercell size, or comparison to available experimental defect levels or doping data are mentioned.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their thorough review and valuable comments on our manuscript. We address the major comment regarding the classification of carrier types below.

read point-by-point responses
  1. Referee: [Abstract / Electronic properties results] The classification of p-type (V_S, V_Mo, As_Mo, As_S) and n-type (As interstitial) behavior is based solely on the Fermi-level position in the neutral-supercell DOS. Standard defect physics requires formation-energy curves E_f(q, E_F) for multiple charge states q, solution of the charge-neutrality condition, and 2D image-charge corrections to locate the equilibrium E_F. The abstract and reported results give no indication that these steps were performed; the neutral-DOS shift alone does not determine carrier type when defect levels are deep or when formation energies cross inside the gap. This directly undermines the application claims in the final paragraph.

    Authors: We agree with the referee that determining the equilibrium carrier type rigorously requires computing the formation energies as a function of charge state and Fermi level, solving the charge neutrality condition, and applying appropriate corrections for the 2D system. Our study relied on the Fermi level position within the neutral supercell density of states, which provides an initial indication of the doping behavior commonly used in the literature for similar systems. However, we recognize that this approach has limitations, especially for deep midgap states. We will revise the abstract and the discussion section to explicitly state that the p-type and n-type assignments are based on the neutral supercell Fermi level shift and to moderate the application claims, indicating that a more comprehensive charged-defect analysis would be required to confirm the thermodynamic equilibrium carrier concentrations. revision: partial

Circularity Check

0 steps flagged

No circularity; standard forward DFT simulation of defect DOS

full rationale

The paper reports direct outputs from DFT calculations on MoS2 supercells with vacancies and As substitutions/interstitials: optimized geometries, total energies, and neutral-supercell DOS plots from which Fermi-level positions are read off. No equations, fitted parameters, or self-citations are used to derive the reported p-type/n-type assignments; the shifts are literal computational results. The study contains no derivation chain that reduces any claim to its own inputs by construction, satisfying the self-contained benchmark criterion.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only the abstract is available; typical first-principles studies rest on standard DFT approximations whose validity cannot be checked here.

axioms (1)
  • domain assumption Density functional theory with the chosen functional and pseudopotentials yields reliable defect formation energies and electronic level positions for the modeled supercells.
    This is the implicit foundation of any DFT defect study; the abstract provides no details on functional choice or convergence.

pith-pipeline@v0.9.1-grok · 5718 in / 1241 out tokens · 45448 ms · 2026-07-01T06:40:50.880778+00:00 · methodology

0 comments
read the original abstract

This study is aimed at exploring the structural and electronic properties of doped MoS$_2$ monolayers, including Mo and S vacancies and As doped systems, employing DFT calculations. The electronic properties were analyzed to understand how these modifications affect the behavior of the material. Introduction of defects generates new defect states in the midgap. In the S-vacancy (V$_\text{S}$), Mo-vacancy (V$_{\text{Mo}}$), As-Mo (As substituting Mo), and As-S (As substituting S) doped systems, the downward shift of the Fermi level to the valence band indicates a $p$-type behavior. In the As interstitial system the Fermi level shifts to the conduction band, suggesting an $n$-type semiconductor. The results highlight that doping MoS$_2$ with As, particularly at the Mo site, can be used in photocatalysis and high-efficiency photovoltaics. Additionally, the As interstitial system demonstrates an enhanced performance in field-effect transistors (FETs).

Figures

Figures reproduced from arXiv: 2606.22119 by A. Daouadi, M. L. Benkhedir.

Figure 1
Figure 1. Figure 1: (Colour online) Top and side views of geometric structures of 3×3×1 supercell with 27 atoms: (a) undoped MoS2, (b) VMo system, (c) VS system, (d) As-Mo doped system, (e) As-S doped system and (f) As-interstitial system. blue, yellow and red are Mo, S, and As atoms. 23702-4 [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (Colour online) Band structure (a), total and partial DOS of undoped monolayer MoS2 (b). Fermi level is set at 0 eV. The difference between the band gap values from the band structure and DOS mainly results from an insufficiently dense 𝑘-mesh [36] and a large smearing parameter that blurs DOS features and reduces the apparent gap [37], which tends to blur the sharp features of the DOS by introducing a tail… view at source ↗
Figure 3
Figure 3. Figure 3: (Colour online) Band structure (a), total and partial DOS of Mo-vacancy in monolayer MoS2 (b). Fermi level is set at 0 eV. Γ M K Γ -2 -1 0 1 2 Energy (eV) 0 10 20 30 40 DOS (states/eV) -2 -1 0 1 2 TDOS Mo-4d S-3p (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (Colour online) Band structure (a), total and partial DOS of S-vacancy in monolayer MoS2 (b). Fermi level is set at 0 eV [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (Colour online) Band structure (a), total and partial DOS of monolayer MoS2 with As doping at the Mo site (b). Fermi level is set at 0 eV. explanation of the origin of defects, the contribution from the As atom was not observed due to the PDOS contained contributions from all the S and Mo atoms, owing to the presence of only one As atom. To gain a better understanding, we plotted PDOS for the S and Mo atom… view at source ↗
Figure 6
Figure 6. Figure 6: (Colour online) Geometric structures and PDOS of atoms S and Mo neighboring As: As-Mo doped system. Fermi level is set at 0 eV [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Band structure (a), total and partial DOS of monolayer MoS2 with As doping at the S site (b). Fermi level is set at 0 eV. -4 -3 -2 -1 0 1 2 3 4 Energy (eV) 0 1 2 3 4 5 6 7 8 9 DOS (States/eV) 0 Mo-d Mo-p Mo-s As-p As-s S-p S-s -4 -3 -2 -1 0 1 2 3 4 Energy (eV) 0 1 2 3 4 5 6 7 8 9 DOS (States/eV) 0 Mo-d Mo-p Mo-s As-p As-s S-p S-s [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: (Colour online) Geometric structures and PDOS of atoms S and Mo neighboring As: As-S doped system. Fermi level is set at 0 eV. Γ M K Γ -3 -2 -1 0 1 2 Energy (eV) 0 10 20 30 40 DOS (states/eV) -3 -2 -1 0 1 2 TDOS Mo-4d S-3p As-4p As-4s (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: (Colour online) Band structure (a), total and partial DOS of monolayer MoS2 with As interstitial doping (b). Fermi level is set at 0 eV. as detailed in table 1. No clear and direct contribution from the As orbitals is observed in the PDOS. However, more significant changes in electronic density become more pronounced when analyzing the 23702-8 [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 9
Figure 9. Figure 9: (Colour online) Band structure (a), total and partial DOS of monolayer MoS2 with As interstitial doping (b). Fermi level is set at 0 eV. -3 -2 -1 0 1 2 3 Energy (eV) 0 1 2 3 4 5 6 7 8 9 10 DOS (States/eV) 0 Mo-d Mo-p Mo-s S-p S-s As-p As-s [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: (Colour online) Geometric structures and PDOS of atoms S and Mo neighboring As: As interstitial doping system. Fermi level is set at 0 eV. PDOS of the sulfur and molybdenum atoms neighboring the As atom (see figure 8), where a shoulder that appears attached to the valence band edge located at 0 eV, arises from the contribution of the Mo-4𝑑 and As-4𝑝 and As-4𝑠 orbitals. Moreover, the shift of the Fermi lev… view at source ↗

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

44 extracted references · 43 canonical work pages · 1 internal anchor

  1. [1]

    S., Geim A

    Novoselov K. S., Geim A. K., Morozov S. V., Jiang, Zhang Y., Dubonos S. V., Grigorieva I. V., Firsov A. A., Science, 2004,306, No. 5696, 666–669, doi:10.1126/science.1102896

  2. [2]

    Y., Özyilmaz B., Zhang Y., Kim P., Phys

    Han M. Y., Özyilmaz B., Zhang Y., Kim P., Phys. Rev. Lett., 2007,98, No. 20, 206805, doi:10.1103/PhysRevLett.98.206805

  3. [3]

    W., Hone J., Science, 2008,321, No

    Lee C., Wei X., Kysar J. W., Hone J., Science, 2008,321, No. 5887, 385–388, doi:10.1126/science.1157996

  4. [4]

    CastroNetoA.H.,GuineaF.,PeresN.M.R.,NovoselovK.S.,GeimA.K.,Rev.Mod.Phys.,2009,81,109–162, doi:10.1103/RevModPhys.81.109

  5. [5]

    A., Yoffe A., Adv

    Wilson J. A., Yoffe A., Adv. Phys., 1969,18, No. 73, 193–335, doi:10.1080/00018736900101307

  6. [6]

    H., Kalantar-Zadeh K., Kis A., Coleman J

    Wang Q. H., Kalantar-Zadeh K., Kis A., Coleman J. N., Strano M. S., Nat. Nanotechnol., 2012,7, No. 11, 699–712, doi:10.1038/nnano.2012.193

  7. [7]

    DasS.,ZhangW.,DemarteauM.,HoffmannA.,DubeyM.,RoelofsA.,NanoLett.,2014,14,No.10,5733–5739, doi:10.1021/nl5025535

  8. [8]

    Kuc A., Heine T., Chem. Soc. Rev., 2015,44, 2603–2614, doi:10.1039/C4CS00276H

  9. [9]

    11, 2254–2260, doi:10.1016/j.physb.2011.03.044

    Ding Y., Wang Y., Ni J., Shi L., Shi S., Tang W., Physica B, 2011,406, No. 11, 2254–2260, doi:10.1016/j.physb.2011.03.044

  10. [10]

    Yang J., Kawai H., Wong C. P. Y., Goh K. E. J., J. Phys. Chem. C, 2019,123, No. 5, 2933–2939, doi:10.1021/acs.jpcc.8b10496

  11. [11]

    Feng L.-p., Sun H.-q., Li A., Su J., Zhang Y., Liu Z.-t., Mater. Chem. Phys., 2018,209, 146–151, doi:10.1016/j.matchemphys.2018.01.015

  12. [12]

    Mater., 2014,26, No

    FengQ.,ZhuY.,HongJ.,ZhangM.,DuanW.,MaoN.,WuJ.,XuH.,DongF.,LinF.,JinC.,WangC.,ZhangJ., Xie L., Adv. Mater., 2014,26, No. 17, 2648–2653, doi:10.1002/adma.201306095

  13. [13]

    Noh J.-Y., Kim H., Park M., Kim Y.-S., Phys. Rev. B, 2015,92, No. 11, 115431, doi:10.1103/PhysRevB.92.115431

  14. [14]

    Cheng Y., Zhu Z., Mi W., Guo Z., Schwingenschlögl U., Phys. Rev. B, 2013,87, No. 10, 100401, doi:10.1103/PhysRevB.87.100401

  15. [15]

    Yue Q., Chang S., Qin S., Li J., Phys. Lett. A, 2013,377, No. 19-20, 1362–1367, doi:10.1016/j.physleta.2013.03.034

  16. [16]

    Electron Devices, 2018,65, No

    Tsai Y.-C., Li Y., IEEE Trans. Electron Devices, 2018,65, No. 2, 733–738, doi:10.1109/TED.2017.2782667

  17. [17]

    M., Mathew S., ACS Omega, 2023,8, No

    Balan B., Xavier M. M., Mathew S., ACS Omega, 2023,8, No. 29, 25649–25673, doi:10.1021/acsomega.3c02084

  18. [18]

    TongX.,AshalleyE.,LinF.,LiH.,WangZ.M.,NanoMicroLett.,2015,7,203–218,doi:10.1007/s40820-015- 0034-8

  19. [19]

    Nanotechnol., 2013,8, No

    Lopez-Sanchez O., Lembke D., Kayci M., Radenovic A., Kis A., Nat. Nanotechnol., 2013,8, No. 7, 497–501, doi:10.1038/nnano.2013.100

  20. [20]

    K., Friedman A

    Perkins F. K., Friedman A. L., Cobas E., Campbell P., Jernigan G., Jonker B. T., Nano Lett., 2013,13, No. 2, 668–673, doi:10.1021/nl3043079

  21. [22]

    L., Cococcioni M., Dabo I., et al., J

    Giannozzi P., Baroni S., Bonini N., Calandra M., Car R., Cavazzoni C., Ceresoli D., Chiarotti G. L., Cococcioni M., Dabo I., et al., J. Phys.: Condens. Matter, 2009,21, No. 39, 395502, doi:10.1088/0953- 8984/21/39/395502

  22. [23]

    Kaltenbacher, W

    Giannozzi P., Andreussi O., Brumme T., Bunau O., Nardelli M. B., Calandra M., Car R., Cavazzoni C., Ceresoli D., Cococcioni M., et al., J. Phys.: Condens. Matter, 2017,29, No. 46, 465901, doi:10.1088/1361- 648X/aa8f79

  23. [24]

    Baerends E., Theor. Chem. Acc., 2001, 265–269, doi:10.1007/s002140050031

  24. [25]

    Generalized gradient approximation made simple,

    PerdewJ.P.,BurkeK.,ErnzerhofM.,Phys.Rev.Lett.,1996,77,No.18,3865,doi:10.1103/PhysRevLett.77.3865. 23702-10 First-principles study of the impact of as doping

  25. [26]

    P., Ruzsinszky A., Csonka G

    Perdew J. P., Ruzsinszky A., Csonka G. I., Vydrov O. A., Scuseria G. E., Constantin L. A., Zhou X., Burke K., Phys. Rev. Lett., 2008,100, No. 13, 136406, doi:10.1103/PhysRevLett.100.136406

  26. [27]

    Catal., 2000,189, No

    Raybaud P., Hafner J., Kresse G., Kasztelan S., Toulhoat H., J. Catal., 2000,189, No. 1, 129–146, doi:10.1006/jcat.1999.2698

  27. [28]

    WangY.,TsengL.-T.,MurmuP.P.,BaoN.,KennedyJ.,IonescM.,DingJ.,SuzukiK.,LiS.,YiJ.,Mater.Des., 2017,121, 77–84, doi:10.1016/j.matdes.2017.02.037

  28. [29]

    W., Nørskov J

    Bollinger M., Jacobsen K. W., Nørskov J. K., Phys. Rev. B, 2003,67, No. 8, 085410, doi:10.1103/PhysRevB.67.085410

  29. [30]

    B., ACS Nano, 2011,5, No

    Johari P., Shenoy V. B., ACS Nano, 2011,5, No. 7, 5903–5908, doi:10.1021/nn201698t

  30. [31]

    V., Phys

    Komsa H.-P., Krasheninnikov A. V., Phys. Rev. B, 2015,91, No. 12, 125304

  31. [32]

    M., Yakobson B

    Zhou W., Zou X., Najmaei S., Liu Z., Shi Y., Kong J., Lou J., Ajayan P. M., Yakobson B. I., Idrobo J.-C., Nano Lett., 2013,13, No. 6, 2615–2622, doi:10.1021/nl4007479

  32. [33]

    D., Sanvito S., Phys

    Dolui K., Rungger I., Pemmaraju C. D., Sanvito S., Phys. Rev. B, 2013,88, 075420, doi:10.1103/PhysRevB.88.075420

  33. [34]

    Appl., 2019,3, No

    Zhong M., Shen C., Huang L., Deng H.-X., Shen G., Zheng H., Wei Z., Li J., npj 2D Mater. Appl., 2019,3, No. 1, 1, doi:10.1038/s41699-018-0083-1

  34. [35]

    Lebegue S., Eriksson O., Phys. Rev. B, 2009,79, No. 11, 115409, doi:10.1103/PhysRevB.79.115409

  35. [36]

    Setyawan W., Curtarolo S., Comput. Mater. Sci, 2010,49, No. 2, 299–312, doi:10.1016/j.commatsci.2010.05.010

  36. [37]

    Toriyama M., Ganose A., Dylla M., Anand S., Park J., Brod M., Munro J., Persson K., Jain A., Snyder G., Preprint arXiv:2103.03469, 2021

  37. [38]

    Sci., 2016,651, 215–221, doi:10.1016/j.susc.2016.05.003

    Salehi S., Saffarzadeh A., Surf. Sci., 2016,651, 215–221, doi:10.1016/j.susc.2016.05.003

  38. [39]

    C., Addou R., Wallace R

    Santosh K., Longo R. C., Addou R., Wallace R. M., Cho K., Nanotechnology, 2014,25, No. 37, 375703, doi:10.1088/0957-4484/25/37/375703

  39. [40]

    G., Ullah S., Phys

    Menezes M. G., Ullah S., Phys. Rev. B, 2021,104, No. 12, 125438, doi:10.1103/PhysRevB.104.125438

  40. [41]

    J., Chen Y., Ko C., Jang C., Sun Y., Sinclair R., Chang J., Tongay S., Wu J., Nano Lett., 2014,14, No

    Suh J., Park T.-E., Lin D.-Y., Fu D., Park J., Jung H. J., Chen Y., Ko C., Jang C., Sun Y., Sinclair R., Chang J., Tongay S., Wu J., Nano Lett., 2014,14, No. 12, 6976–6982, doi:10.1021/nl503251h

  41. [42]

    Xie W.-J., Li X., Zhang F.-J., Chem. Phys. Lett., 2020,746, 137276, doi:10.1016/j.cplett.2020.137276

  42. [43]

    BeshirB.T.,ObodoK.O.,AsresG.A.,RSCAdv.,2022,12,No.22,13749–13755,doi:10.1039/D2RA00775D

  43. [44]

    LiuQ.,YouJ.,XiongY.,LiuW.,SongM.,RenJ.,XueQ.,TianJ.,ZhangH.,WangX.,J.ColloidInterfaceSci., 2024,675, 772–782, doi:10.1016/j.jcis.2024.07.044

  44. [45]

    Mano P., Jitwatanasirikul T., Roongcharoen T., Takahashi K., Namuangruk S., Appl. Surf. Sci., 2025,688, 162339, doi:10.1016/j.apsusc.2025.162339. 23702-11 A. Daouadi, M. L. Benkhedir Дослiдження впливу легування арсеном на структурнi та електроннi властивостi моношару MoS2 з перших принципiв А. Даудi, М. Л. Бенхедiр Лабораторiя теоретичної та прикладної ф...