{"id":"ba5ca2d5-c527-4300-9287-0cc490c9d559","arxiv_id":"2607.18895","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"An ensemble Monte Carlo model of monolayer MoS2 reproduces measured transistor and gas-sensing data, but its absolute mobility and sensor curves rely on fitted scattering and adsorption parameters.","lead":"This paper presents a Monte Carlo simulation framework for atomically thin MoS2 that connects first-principles material parameters to transistor and gas-sensor measurements. A generalist might read it to see how far physics-based simulation can go in designing 2D semiconductor devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The single calibrated acoustic deformation potential (Eq. 5) is the load-bearing free parameter; Sec. 5.6 explicitly concedes it is not first-principles-validated, so the downstream 'predictions' are partly consistency checks.","rationale":"The reader's weakest-assumption analysis points to the same load-bearing issue: D_ac is calibrated, not derived from first principles, and the paper explicitly concedes this in Sec. 5.6. My read of the full text confirms that this parameter is not merely cosmetic. It enters Eq. (5) quadratically and propagates into the mobility benchmark, the velocity-field curve, the device I–V, and the ambient/sensor baseline. Because the conclusion's strongest claim emphasizes a physically grounded and reusable basis, an unverified calibration at the root of the transport chain is a genuine soft spot. I considered whether the compact-model substitution (Sec. 4.1) is more load-bearing, but that substitution affects only the device-level claim and is self-identified; D_ac affects all transport results, including the intrinsic mobility and substrate ordering that are central to the paper. The paper deserves credit for transparency: Table 3 and Sec. 5.6 clearly separate predicted from fitted quantities, and the DFPT verification of seven electronic-structure and polar-phonon parameters is a real, parameter-free check. The sensor section also uses per-device adsorption constants while leaving the transduction shape as an EMC output. These features support a conditional rather than hostile reading. My recommendation is therefore unchanged from the reader's conditional verdict, with the concrete test being an independent first-principles determination of the acoustic deformation potential and a no-refit rerun of the pipeline.","tokens_in":12427,"tokens_out":3947,"duration_ms":37827,"concrete_test":"Obtain the effective acoustic deformation potential for monolayer MoS2 from a published or newly computed EPW/Wannier electron-phonon calculation (or directly compute the acoustic-limited mobility from first principles). Insert that value into the EMC kernel in place of the calibrated D_ac, refit nothing, and recompute the temperature-dependent mobility in Fig. 2 and the device transfer characteristic in Fig. 5(b). If the room-temperature mobility shifts by more than ~20% or the power-law exponent γ changes by more than ~0.1, the single global calibration is not a benign scale factor and the central predictive claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Equation (5) sets λ_ac ∝ D_ac², and D_ac enters every transport result: μ(T) in Figs. 2–3, the velocity-field characteristic feeding the device model in Fig. 5, and the mobility baseline of the ambient/sensor extension in Eq. (13). The paper's own Sec. 5.6 flags that the frozen-strain ≈13 eV deformation potential is 'a different quantity from the effective scattering deformation potential that enters the mobility,' and that the DFPT comparison 'confirms the electronic-structure and polar-phonon inputs but not the acoustic scattering coupling.' Thus the central claim of a physically grounded, non-empirical framework rests on one calibrated constant whose value is not reported and whose transferability across densities, temperatures, and devices is assumed rather than demonstrated. Because D_ac changes the relative weight of acoustic versus polar-optical and remote-phonon scattering, a calibration chosen to match room-temperature mobility could absorb errors in other rates and still produce the right low-field value while giving incorrect temperature dependence or substrate ordering. The device and sensor reproductions then inherit this calibration, so they do not independently validate the first-principles basis. The relative trends remain plausible, but the strongest claim of full predictive grounding is weaker than stated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a self-consistent ensemble Monte Carlo (EMC) framework for monolayer MoS2, implemented in ViennaEMC, with a full scattering stack (acoustic and intervalley deformation-potential phonons, polar-optical Fröhlich and piezoelectric coupling, remote substrate phonons, screened charged impurities, surface roughness, and adsorbates), degenerate statistics, and self-consistent Poisson coupling. The intrinsic electronic-structure and polar-phonon inputs are compared to DFPT calculations, and seven quantities are claimed to match the reference parametrization. The EMC velocity-field characteristic feeds a one-dimensional velocity-saturation channel model, used to reproduce the transfer characteristic of a CVD monolayer device in air and vacuum. An ambient/adsorbate extension is then used to reproduce oxygen partial-pressure conductivity, NO2, and NH3 gas-sensing data from several devices. The manuscript explicitly separates predicted from fitted quantities and concludes that the framework is a physically grounded, reusable basis for 2D-material FET and gas-sensor modelling.","tokens_in":12756,"tokens_out":7064,"duration_ms":68368,"significance":"The paper makes a useful open-source contribution: it ships an EMC transport kernel plus DFPT input decks and simulation data, and it provides a genuine DFPT check of seven electronic-structure and polar-phonon parameters. Several outputs are genuine EMC predictions rather than fits: the temperature exponents and substrate ordering in Figs. 2 and 4, the non-monotonic density dependence with a peak near 10^13 cm^-2, the shape of the velocity-field characteristic v(E), and the shape S(theta) of the adsorbate transduction curve. These are real strengths. However, the load-bearing acoustic deformation potential D_ac in Eq. (5) is a single calibrated constant that the paper itself concedes is not DFPT-verified, and the device and sensor validations depend on additional fitted parameters (V_th, SS, R_c, DeltaN, K). The claim that the framework is 'physically grounded rather than empirical' is therefore only partially supported: the relative trends are plausible, but the absolute mobility and the quantitative sensor response are not first-principles predictions in the strong sense stated in the abstract.","major_comments":[{"comment":"The calibrated D_ac is load-bearing and undermines the 'essentially parameter-free' claim. Eq. (5) makes the acoustic scattering rate proportional to D_ac^2, and this rate enters every transport result: mu(T) in Figs. 2-3, the velocity-field characteristic feeding Fig. 5, and the mobility baseline in Eq. (13). Section 5.6 explicitly concedes that the DFPT comparison 'confirms the electronic-structure and polar-phonon inputs but not the acoustic scattering coupling,' and that the frozen-strain ~13 eV value is a different quantity. Yet Table 3 lists mu(T) as resting on a 'single global calibration' while the text describes the backbone as 'essentially parameter-free.' Calibrating D_ac at room temperature can absorb errors in other scattering rates and still give the correct low-field mobility while producing incorrect temperature dependence or substrate ordering. The authors should report","section":"Sec. 5.6, Eq. (5), Table 3"},{"comment":"The device-level validation is substantially weaker than the abstract implies. The manuscript states that direct EMC terminal-current extraction 'does not yield a clean Id-Vg' and was abandoned; the transfer characteristic is instead computed from a 1D gradual-channel compact model with fitted V_th, SS, and R_c, using the bulk EMC v(E) as input. The agreement in Fig. 5(b) therefore tests the compact model plus fitted device parameters more than the scattering stack. To support the multiscale claim, the authors should quantify the sensitivity of the reproduced Id-Vg to V_th/SS/R_c, and ideally show that the EMC v(E) differs from a generic Caughey-Thomas form in a way that measurably affects the transferred characteristic. Without that, the device panel does not independently corroborate the transport kernel.","section":"Sec. 4.1, Eqs. (10)-(12)"},{"comment":"The gas-sensing 'validation' is partly a curve-fitting exercise. The EMC output is the transduction shape S(theta), but the mapping theta(C)=KC/(1+KC) uses K fitted per device (and per dataset), DeltaN is fitted per gas/device, and N_max_c is introduced without being listed in Table 3. Consequently, the NO2 and NH3 concentration curves in Figs. 7-8 are not predictions from the framework; they are fits to each dataset. A meaningful predictive test would be to fix DeltaN and N_max_c from one device and K from an independent adsorption measurement or from one pressure point, and then predict another device's concentration response. As written, the sensor section demonstrates interpolation within a flexible model, not predictive power.","section":"Sec. 4.2, Eq. (13), Table 3"},{"comment":"The language of the abstract and conclusion overstates the evidential basis. The conclusion says the framework 'establishes a physically grounded and reusable basis,' while Sec. 5.6 correctly acknowledges the main limitations: the acoustic scattering coupling is unverified and the gate-tunable few-layer sensing mechanism is outside scope. Since the only parameter controlling the absolute mobility is calibrated, and direct device simulation failed, the strong claim that the model is 'physically grounded rather than empirical' should be softened to something like 'microscopic scattering-based with transparent parameter accounting.' The paper would be more credible if the calibrated and fitted elements were described as such throughout, rather than only in Sec. 5.6.","section":"Abstract, Conclusion, Sec. 5.6"}],"minor_comments":[{"comment":"The manuscript has several encoding/typesetting problems: 'I d-Vg', 'Fr¨ ohlich', and other LaTeX artifacts appear in the extracted text. Please ensure the final PDF uses consistent math typesetting.","section":"Throughout"},{"comment":"The annotation 'EMC intrinsic (this work) 194' beside the experimental bars is unclear. State explicitly that this is the room-temperature intrinsic mobility at the same density used for the Kaasbjerg limit, and how it relates to the 135-195 cm2/Vs range quoted in Sec. 5.1.","section":"Fig. 3"},{"comment":"The 'in-plane epsilon_inf (slab)' entry says it 'needs 2D unfold.' Since this quantity enters screening in Eq. (7), spell out how the slab value was used, or whether a literature value was substituted.","section":"Table 2"},{"comment":"Define n0, DeltaN, and N_max_c in physical units and clarify the donor/acceptor sign convention in the text before Eq. (13). The upper/lower sign description is terse.","section":"Eq. (13)"},{"comment":"The multilayer and few-layer results (Late et al., Mamun et al.) are described in Sec. 5.6 as corroboration only. Consider labeling them as such in the main text and figure captions as well, or moving them to supplementary material.","section":"Sec. 5.5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest in Sec. 5.6 about its limitations, which is commendable, but the abstract and conclusion systematically overclaim first-principles grounding. The central issue is the single calibrated D_ac: it controls all absolute transport results and is explicitly not DFPT-verified. This can be addressed either by computing D_ac with a full electron-phonon (EPW) treatment or by reframing the paper as a semi-empirical multiscale model and adding sensitivity analyses. I recommend major revision rather than rejection because the open framework, the DFPT parameter checks, and the genuinely predicted relative trends (temperature exponents, substrate ordering, density dependence, S(theta) shape) are valuable contributions. Please ask for a sensitivity analysis around D_ac and the device/sensor fitted parameters, and for a consistent reframing of the 'physically grounded' claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: this paper builds a full ensemble Monte Carlo transport kernel for monolayer MoS2 with a proper scattering stack, DFPT checks on the structural and polar-phonon inputs, and multiscale links to FET and gas-sensor models. It is transparent about what is fitted and what is predicted — Table 3 is exactly the sort of accounting you wish more papers had. The catch, which the paper itself concedes, is that the acoustic deformation potential D_ac in Eq. (5) — the parameter that scales every acoustic-scattering rate in every result — is set by a single global calibration, not computed, and its value is not reported. The device I-V also comes from a 1D gradual-channel model with fitted Vth, SS, and Rc after direct EMC current extraction failed (Sec 4.1). So the phrase \"physically grounded rather than empirical\" is doing more work than the evidence supports.\n\nWhat is genuinely new: no one has put this full 2D scattering stack — acoustic, intervalley, Fröhlich, piezoelectric, remote phonons, charged impurities, roughness — into an open EMC solver and coupled it to device and sensor response in one workflow. The framework is reusable and well documented, and the relative physical trends are real EMC outputs. The substrate ordering (suspended>hBN≈SiO2>Al2O3≈CaF2>HfO2), the non-monotonic density peak around 1e13 cm-2, and the rising temperature exponent on softer dielectrics are all consistent with the physics and would survive changes in the acoustic calibration. The DFPT agreement for seven quantities (lattice constant, mass, gap, LO energy, Born charges, sound velocities) is real and worth having.\n\nThe soft spots are as follows. D_ac is load-bearing: it sets the relative weight of acoustic versus polar-optical scattering, so a room-temperature fit can absorb errors in other rates and still reproduce the correct low-field mobility while distorting temperature dependence or substrate ordering. The sensor section fits K and DeltaN per dataset; the shape S(theta) is an EMC output, but the \"predicted\" concentration curves are not independent. The device reproduction inherits the compact-model fits, so it is not a first-principles validation of transport. None of these are fatal — the paper openly labels them — but they should be acknowledged in the framing. A sensitivity analysis over D_ac and reporting the fitted values would materially strengthen what is already a useful paper.\n\nWho it's for: anyone modeling MoS2 FETs or 2D-material gas sensors, and people standing up EMC workflows for 2D channels. It deserves a serious referee, not a desk reject. I'd send it out with a request for the sensitivity analysis and a tightened claim about predictive status.\n\nCandidly,\n[You]","headline":"Honest, useful EMC framework for MoS2, but the 'physically grounded' claim rests on one calibrated acoustic deformation potential and a compact-model device stand-in.","tokens_in":13225,"tokens_out":2716,"would_cite":true,"duration_ms":49285,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A single ensemble Monte Carlo transport kernel, calibrated once, reproduces measured transport, transistor, and gas-sensing behaviour in monolayer MoS2.","keywords":["ensemble Monte Carlo","transition-metal dichalcogenide","MoS2","electron transport","gas sensing","multiscale modelling","density-functional perturbation theory"],"falsifier":"A first-principles electron-phonon (EPW) calculation of the acoustic deformation potential in monolayer MoS2 that differs materially from the calibrated value, or a measured suspended-device mobility that deviates from the predicted T^(-1.49) power law, would falsify the transport backbone. A simpler check: if the model's v(E) fails to reproduce the transfer characteristic of a device with a different channel length or contact design, the bulk-to-channel coupling would be invalid.","tokens_in":12287,"feed_emoji":"","tokens_out":8118,"duration_ms":60182,"temperature":0.7,"pith_summary":"This paper aims to establish that one semiclassical ensemble Monte Carlo transport framework can span the full range from electron-phonon coupling to terminal current and chemical sensing in monolayer MoS2. The claim matters because it would give device designers a physically grounded, reusable simulation basis for 2D-material transistors and gas sensors, replacing purely empirical mobility models and lumped compact models. The author shows that with a single global calibration—the acoustic deformation potential—the same Boltzmann-transport kernel reproduces substrate-specific mobility power laws, the velocity-field characteristic, the measured transfer curves of a CVD device in air and vacuum, and the concentration-dependent response to O2, NO2, and NH3 across several independent devices. If true, this would unify transport and sensing physics under one computational umbrella.","feed_headline":"One Monte Carlo model reproduces MoS2 devices and gas sensors","feed_subtitle":"One kernel spans electron-phonon coupling, transistor current, and gas sensing in monolayer MoS2.","key_machinery":"The central object is the ensemble Monte Carlo Boltzmann-transport kernel, which stochastically propagates carrier ensembles through field-driven flights and self-scattering events, coupled self-consistently to Poisson's equation with a degenerate 2D closure. Its output velocity-field characteristic v(E), fitted to the Caughey-Thomas form, feeds a one-dimensional gradual-channel model for the transistor drain current, while an adsorbate model couples coverage-dependent charge transfer and Coulomb scattering into the same conductivity. The key identity is the transduction S(θ)=σ0/σ(θ)−1, computed directly by the EMC, with only the per-device adsorption constant K fitted.","core_discovery":"The central discovery is that a self-consistent ensemble Monte Carlo kernel with a full scattering stack—acoustic deformation potential, intervalley phonons, polar-optical Fröhlich and piezoelectric coupling, remote substrate phonons, screened charged impurities, and surface roughness—verified against density-functional perturbation theory for seven parameters, reproduces a broad set of measurements without per-result fitting. The temperature-dependent mobility follows a substrate-specific power law whose exponent rises from γ=1.49 for suspended films to γ=1.79 on HfO2; the multiscale coupling to a one-dimensional velocity-saturation channel model reproduces the measured transfer characteris","pith_inferences":["Because the single calibrated acoustic deformation potential is not first-principles-verified, a future EPW calculation that disagrees would require re-testing the model's transferability across densities and temperatures.","The coverage-dependent transduction curve could be probed with time-resolved or noise measurements to separate the charge-transfer and scattering contributions directly.","The same framework should extend to other monolayer TMDs (WS2, WSe2) with only material parameters changed, offering a strong test of the mechanism.","The bulk v(E) replacing a full gated-channel EMC may break down for very short channels or significantly different contact resistances; a full device simulation with proper terminal extraction would be needed to confirm the multiscale coupling."],"forward_implications":["Substrate choice is predictive: each dielectric yields a distinct mobility power-law exponent, and stiff-mode insulators like CaF2 retain high mobility despite a large dielectric discontinuity.","The intrinsic phonon-limited mobility provides an upper bound that quantifies the extrinsic (contact, defect, substrate) limitation of any measured device.","The same velocity-field characteristic explains the linear-to-saturation transition in a real gated channel and the air→vacuum mobility recovery ratio.","The sensor response shape is an output of the transport kernel, so transistor and sensor behaviour share a common microscopic basis.","The transport backbone requires only one global calibration; remaining fits are confined to device parameters (Vth, SS, Rc) and per-gas adsorption constants."],"fun_headline_variants":["One Monte Carlo model reproduces MoS2 devices and gas sensors","Self-consistent Monte Carlo nails MoS2 transport and sensing","Multiscale MC matches MoS2 transistor and sensor experiments","Ensemble Monte Carlo unifies MoS2 electronics and chemosensing","Single simulation stack reproduces MoS2 transistors and gas response"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the effective acoustic deformation potential, fixed by a single global calibration rather than computed from first principles, remains valid across all carrier densities, temperatures, and devices—a premise the paper itself flags as unverified for the acoustic scattering coupling.","fun_headline_variants_meta":{"raw":{"variants":["One Monte Carlo model reproduces MoS2 devices and gas sensors","Self-consistent Monte Carlo nails MoS2 transport and sensing","Multiscale MC matches MoS2 transistor and sensor experiments","Ensemble Monte Carlo unifies MoS2 electronics and chemosensing","Single simulation stack reproduces MoS2 transistors and gas response"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000235,"raw_usage":{"total_tokens":1377,"prompt_tokens":824,"completion_tokens":553,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":568,"completion_tokens_details":{"reasoning_tokens":464}},"tokens_in":568,"tokens_out":553,"duration_ms":4655,"temperature":1.0,"reasoning_tokens":464,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T13:59:29.376223+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A first-principles electron-phonon (EPW) calculation of the acoustic deformation potential in monolayer MoS2 that differs materially from the calibrated value, or a measured suspended-device mobility that deviates from the predicted T^(-1.49) power law, would falsify the transport backbone. A simpler check: if the model's v(E) fails to reproduce the transfer characteristic of a device with a different channel length or contact design, the bulk-to-channel coupling would be invalid.","supporting_citations":[],"review_version":1}