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 →
First-principles study of the impact of As doping on the structural and electronic properties of MoS₂ monolayer
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
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.
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
- 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.
Referee Report
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)
- [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)
- [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.
- [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
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
-
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
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
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.
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
Reference graph
Works this paper leans on
-
[1]
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
work page internal anchor Pith review doi:10.1126/science.1102896 2004
-
[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]
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]
CastroNetoA.H.,GuineaF.,PeresN.M.R.,NovoselovK.S.,GeimA.K.,Rev.Mod.Phys.,2009,81,109–162, doi:10.1103/RevModPhys.81.109
-
[5]
Wilson J. A., Yoffe A., Adv. Phys., 1969,18, No. 73, 193–335, doi:10.1080/00018736900101307
-
[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]
DasS.,ZhangW.,DemarteauM.,HoffmannA.,DubeyM.,RoelofsA.,NanoLett.,2014,14,No.10,5733–5739, doi:10.1021/nl5025535
-
[8]
Kuc A., Heine T., Chem. Soc. Rev., 2015,44, 2603–2614, doi:10.1039/C4CS00276H
-
[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]
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]
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]
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]
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]
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]
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]
Tsai Y.-C., Li Y., IEEE Trans. Electron Devices, 2018,65, No. 2, 733–738, doi:10.1109/TED.2017.2782667
-
[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]
TongX.,AshalleyE.,LinF.,LiH.,WangZ.M.,NanoMicroLett.,2015,7,203–218,doi:10.1007/s40820-015- 0034-8
-
[19]
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]
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
-
[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
-
[23]
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
-
[24]
Baerends E., Theor. Chem. Acc., 2001, 265–269, doi:10.1007/s002140050031
-
[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
-
[26]
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
-
[27]
Raybaud P., Hafner J., Kresse G., Kasztelan S., Toulhoat H., J. Catal., 2000,189, No. 1, 129–146, doi:10.1006/jcat.1999.2698
-
[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
-
[29]
Bollinger M., Jacobsen K. W., Nørskov J. K., Phys. Rev. B, 2003,67, No. 8, 085410, doi:10.1103/PhysRevB.67.085410
-
[30]
Johari P., Shenoy V. B., ACS Nano, 2011,5, No. 7, 5903–5908, doi:10.1021/nn201698t
-
[31]
V., Phys
Komsa H.-P., Krasheninnikov A. V., Phys. Rev. B, 2015,91, No. 12, 125304
2015
-
[32]
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
-
[33]
Dolui K., Rungger I., Pemmaraju C. D., Sanvito S., Phys. Rev. B, 2013,88, 075420, doi:10.1103/PhysRevB.88.075420
-
[34]
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
-
[35]
Lebegue S., Eriksson O., Phys. Rev. B, 2009,79, No. 11, 115409, doi:10.1103/PhysRevB.79.115409
-
[36]
Setyawan W., Curtarolo S., Comput. Mater. Sci, 2010,49, No. 2, 299–312, doi:10.1016/j.commatsci.2010.05.010
- [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
-
[39]
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
-
[40]
Menezes M. G., Ullah S., Phys. Rev. B, 2021,104, No. 12, 125438, doi:10.1103/PhysRevB.104.125438
-
[41]
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
-
[42]
Xie W.-J., Li X., Zhang F.-J., Chem. Phys. Lett., 2020,746, 137276, doi:10.1016/j.cplett.2020.137276
-
[43]
BeshirB.T.,ObodoK.O.,AsresG.A.,RSCAdv.,2022,12,No.22,13749–13755,doi:10.1039/D2RA00775D
-
[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
-
[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я теоретичної та прикладної ф...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.