{"id":"8bf94555-71b2-4fdd-be75-30e984971351","arxiv_id":"2506.01464","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"DFT calculations predict that Mo1-xWxS2 monolayer alloys retain a direct bandgap at the K point that tunes from 1.696 eV (MoS2) to 1.858 eV (WS2), with optical absorption shifting to higher energy as tungsten content increases.","lead":"This paper uses computer simulations to map how mixing molybdenum and tungsten atoms in a single layer of MoS2 changes the material's electronic and optical properties across all mixture ratios. It finds that the alloy keeps a direct bandgap that shifts from about 1.70 to 1.86 eV as tungsten content rises, which matters for making color-tunable sensors and light emitters.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Thermodynamic stability claim rests on a formation energy relative to isolated atoms; the composition-property map also lacks configurational averaging.","rationale":"The reader correctly identified the single-configuration-per-composition issue as the weakest assumption for the bandgap map, and I agree that the missing configurational averaging undermines the precision of the composition-property trend. However, I see the stability claim as the more load-bearing weak point, because it is stated as an absolute conclusion ('thermodynamically stable') in both the abstract and the conclusion, yet Eq. (8) measures cohesion relative to isolated atoms rather than alloy stability relative to the parent compounds. A negative Eb is a necessary but not sufficient condition for thermodynamic stability of the alloy against decomposition, and the paper does not provide the relevant decomposition energy. This is a concrete, checkable gap rather than a mere disagreement with consensus. I would keep the verdict CONDITIONAL: the corrections needed are well-defined and affordable, and the central tunable direct-gap trend is likely to survive, but as written the paper overstates thermodynamic stability and the statistical basis of its composition-property map. My proposed test directly addresses the most load-bearing concern and would either confirm or correct the headline claims.","tokens_in":22327,"tokens_out":1363,"duration_ms":13477,"concrete_test":"Recompute the alloy energetics using the standard mixing enthalpy, ΔH_mix(x) = E(alloy, x) - [(1-x)E(MoS2) + xE(WS2)], for at least 3-5 distinct 4x4 supercell configurations per composition (or SQS structures) and check whether the minimum and average ΔH_mix are positive or negative. Additionally, recompute the direct K-point gap and the fitted bowing parameter for the same configurations to quantify the configuration-to-configuration spread.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's headline claims that Mo1-xWxS2 monolayer alloys are 'thermodynamically stable' across the full composition range. The only quantitative support is the formation energy Eb in Eq. (8), referenced to isolated atoms, with values from about -6.91 to -7.37 eV/atom. These large negative values mostly reflect the strong Mo-S and W-S bonds and do not probe alloy stability, which is the energy of the alloy relative to the competing parent phases (MoS2 and WS2) or to phase separation. The paper itself states in Section III A that 'an exhaustive search for the absolute ground-state configuration was not performed,' so the relevant mixing enthalpy is never computed. A positive mixing enthalpy would not contradict the negative Eb values, but it would contradict the stability claim as stated. The same single-configuration limitation affects the bandgap map: the direct-gap, near-monotonic trend, and bowing parameter b=0.1564 are derived from one hand-picked arrangement per composition, with no SQS construction or ensemble averaging, even though the paper acknowledges in Section II that properties 'might exhibit some dependence on the specific configuration chosen.' Both issues point to overstatement of what the calculations establish, not to an obvious internal contradiction in the DFT results themselves.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports DFT-PBE with DFT-D3 dispersion corrections (Quantum ESPRESSO) for monolayer Mo1-xWxS2 alloys represented by 4x4x1 (48-atom) supercells at ten compositions (x = 0, 0.0625, 0.125, 0.25, 0.375, 0.5, 0.625, 0.75, 0.875, 1), with one hand-picked configuration per composition. Calculated quantities include lattice parameters and bond lengths, formation energies Eb relative to isolated atoms (Eq. 8), K-point bandgaps, projected densities of states, and RPA optical spectra (epsilon1, epsilon2, alpha, n, k, R, sigma) up to 8 eV. The central claims are that the alloys are thermodynamically stable across the full composition range, that the direct K-point gap is preserved for all compositions, that the gap increases near-monotonically from 1.696 eV (MoS2) to 1.858 eV (WS2) with a downward bowing parameter b = 0.1564, and that the main optical features blueshift with increasing W content.","tokens_in":22512,"tokens_out":28530,"duration_ms":246193,"significance":"If the quantitative claims were properly supported, the systematic full-composition map would be a useful reference for bandgap and optical engineering in 2D TMDC alloys; the retention of the direct gap and the small bowing are consistent with experimental PL studies reporting near-linear composition tuning in CVD-grown Mo1-xWxS2 [23]. The paper's strengths are its systematic coverage of ten compositions (versus the two or three studied in most prior work), a consistent supercell methodology, explicit acknowledgment of key limitations (neglect of excitonic effects and configurational averaging), and falsifiable predictions for the composition dependence of the bandgap and optical spectra that can be benchmarked against experiment. The endpoint Kohn-Sham gaps (1.696 and 1.858 eV) and the general bandgap trend are plausible and consistent with earlier PBE calculations for the pure monolayers, so the electronic-structure core is likely sound. However, the thermodynamic-stability claim is not established by the computed formation energies, and the optical-response results contain internal inconsistencies that must be resolved.","major_comments":[{"comment":"The stability claim in the abstract and conclusion that the alloys are 'thermodynamically stable throughout the full compositional range' is not supported by the quantity actually computed. The formation energy Eb in Eq. (8) references isolated atoms; it shows only that the compounds are bound relative to free atoms. The relevant criterion is the mixing enthalpy relative to the parent compounds, Delta_H_mix/atom = Eb(x) - [(1-x)Eb(0) + xEb(1)], where the Eb values from Table I correspond to the same 48-atom supercells so the atomic-reference terms cancel. Using the paper's own Table I values, this difference is positive for x >= 0.375 (e.g., about +0.010 eV/atom at x = 0.5 and +0.033 eV/atom at x = 0.875), meaning the alloy configurations are higher in energy than the phase-separated mixture of MoS2 and WS2, the opposite of the stated conclusion. Because the mixing entropy at typical growth temperatures is of the same order as these energy differences, a defensible stability statement requires the proper enthalpy reference and at least a finite-temperature T*Delta_S_mix discussion. Please compute and report Delta_H_mix and revise the stability claim accordingly.","section":"Section III A, Eq. (8), Table I"},{"comment":"The composition-property map is built from a single hand-picked configuration per composition, and the manuscript itself acknowledges in Section II that the properties 'might exhibit some dependence on the specific configuration chosen' and in Section III A that 'an exhaustive search for the absolute ground-state configuration was not performed.' Consequently, the reported gap values, the 'near-monotonic' trend, and the bowing parameter are not averaged over configurations or computed for ground-state configurations, and no uncertainty is attached to them. The magnitude of configurational fluctuations is not bounded by the calculations: the non-monotonic step in the reported gaps (1.741 eV at x = 0.375 versus 1.737 eV at x = 0.5) is only 4 meV, while the residuals of the quadratic fit reach about 21 meV (next comment), both well above the stated 0.01 eV/atom energy-convergence threshold. Please quantify the configuration sensitivity, for example by computing several distinct configurations or special quasirandom structures at each composition and reporting the spread, and adjust the strength of the quantitative claims accordingly.","section":"Section II; Section III A; Section III B"},{"comment":"The quadratic fit with b = 0.1564 shown in Figure 3 does not describe the reported data well. Evaluating Eg(x) = (1-x)*1.696 + x*1.858 - 0.1564*x*(1-x) at the compositions in Table I gives residuals of about +16 meV at x = 0.25 (fit 1.707 eV versus data 1.723 eV) and about +21 meV at x = 0.375 (fit 1.720 eV versus data 1.741 eV). The effective bowing required to pass through individual data points ranges from about 0.07 at x = 0.375 to about 0.18 at x = 0.875, so the single-bowing-parameter description is not a faithful summary of the data. Please report the fit residuals; if the scatter is configurational in origin it supports the concern in the preceding comment, and if it is intrinsic, the composition dependence of the gap is not simply quadratic.","section":"Section III B, Figure 3"},{"comment":"The optical results contain internal inconsistencies that affect the claimed systematic blueshift. Table III lists Epeak (the energy of the extinction-coefficient peak 'near the A exciton region') as 2.14, 2.36, 2.36, 2.08, 2.15, 2.04, 2.04, 2.40, 2.12 and 2.26 eV for the ten compositions; these values are non-monotonic and do not track the bandgap trend. The figure captions (for example Figs. 5, 7, 9 and 11) report epsilon2 A-feature positions that blueshift systematically (about 2.14 to 2.25 eV for H1 through H14), but these disagree with Table III for the same compositions (e.g., H4: 2.17 versus 2.08 eV; H8: 2.20 versus 2.04 eV; H12: 2.23 versus 2.40 eV). In addition, Section III C states that the Epeak shift corresponds to a wavelength tuning of about 30 nm, while the Conclusion claims 'approximately 100 nm (~550-650 nm based on Epeak shifts)'; the endpoint Epeak shift from 2.14 to 2.26 eV corresponds to 579 to 549 nm, a 30 nm span, not 100 nm. Please reconcile the two sets of peak positions and report one consistent value for the wavelength tuning range.","section":"Section III C, Table III; Conclusion"}],"minor_comments":[{"comment":"The text refers to a 'Projector Augmented Wave (PAW) pseudopotential,' but PAW is an augmentation method rather than a pseudopotential in the norm-conserving or ultrasoft sense; the corresponding objects in Quantum ESPRESSO are PAW data sets. Please correct the terminology and specify which PSLibrary data-set flavors were used for the geometry/electronic and optical runs.","section":"Section II"},{"comment":"The sentence stating that Eb reaches 'a minimum value of -7.280 eV/atom ... near the W-rich end (x = 0.875)' is confusing because the quoted range of Eb is from -6.910 to -7.370 eV/atom and the most negative value in Table I is that of WS2. Please reword, for example to 'the most negative value among the alloy configurations.'","section":"Section III A"},{"comment":"The parenthetical in Section III A reading '(as per Table I value for WS2, previously -7.24, then -7.399)' appears to be leftover editing notes; please remove it and ensure the text agrees with Table I.","section":"Section III A"},{"comment":"The axis labels in the optical figures contain stray characters (for example 'xx x 10^8' in the absorption panels), and the numeric annotations ('1', '2', '3', ...) attached to the curves are not explained in the captions. Please repair the labels and define or remove the annotations.","section":"Figures 5-12"},{"comment":"The note stating that missing kpeak data are denoted by '-' is never used in the table; either populate the missing entries or remove the note.","section":"Table III"},{"comment":"For the RPA optical spectra, please state the k-point grid and Gaussian broadening used in the epsilon.x calculation and report a brief convergence check for the peak positions in epsilon2, since the optical trends are the basis of the claims in Section III C.","section":"Section II"},{"comment":"The statement that the computed 0.162 eV tuning range 'aligns remarkably well' with the experimental PL range of about 0.18 eV (Ref. [23]) should be framed as a trend-level comparison, because PBE yields Kohn-Sham gaps while PL probes excitonic transitions and the agreement of the differences could be partly fortuitous.","section":"Section III B"}],"recommendation":"major_revision","confidential_remarks":"I believe the paper is within the journal's scope and the bandgap trend is likely correct, but the manuscript currently overstates the thermodynamic stability, and the optical tables and figures need reconciliation. The main risk to acceptance is the stability claim, which is directly contradicted by the paper's own formation energies when converted to a mixing enthalpy; this is fixable with additional total-energy comparisons and a softened wording. I would also ask the authors to check the novelty claim ('first comprehensive DFT study'), since prior computational work on Mo1-xWxS2 alloys exists beyond Ref. [17]; a more measured phrasing would be safer. For reproducibility, I recommend that the Quantum ESPRESSO input files be deposited rather than made available only upon request."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the useful part: this paper gives a consistent PBE-D3 map of monolayer Mo1-xWxS2 over the entire composition range, with lattice parameters, band gaps, and RPA optical spectra at eight intermediate compositions. That's an extension over earlier work that mostly hit x=0, 0.5, 1. The band gap tuning range they compute (1.696 to 1.858 eV, about 0.16 eV) matches the experimental PL tuning from the CVD alloy paper (Ref. 23) well, and the structural parameters look sensible. The paper is also fairly transparent about its own limits: Section II states that properties 'might exhibit some dependence on the specific configuration chosen,' and Section III A admits that no exhaustive search for the ground-state configuration was performed. That honesty is to their credit.\n\nThe soft spots are real, though. The abstract and conclusion say the alloys are 'thermodynamically stable,' but the only quantitative support is the formation energy in Eq. (8), referenced to isolated atoms. That quantity mostly reflects the strength of the metal–chalcogen bonds; it says nothing about whether the alloy is stable against phase separation into MoS2 and WS2. To make the stability claim stick, they need a mixing enthalpy relative to the parent compounds, which they have not computed. That is a headline claim, not a side remark.\n\nSecond, the entire band gap/composition curve and the bowing parameter b=0.1564 eV rest on one hand-picked configuration per composition. Without an ensemble average or SQS, the small dip at x=0.5 could easily be a configurational artifact. The trend is probably robust, but the precision implied by a fitted bowing parameter is not supported.\n\nThird, the optical section labels peaks as A and B excitons, but the calculation is scalar-relativistic and has no spin-orbit coupling. The physical A/B splitting in these monolayers is spin-orbit in origin. The peaks they see are just interband transitions; the nomenclature is misleading, even though they do note the RPA lacks excitonic effects.\n\nWho is this for? A reader who wants a quick engineering reference for how lattice parameters, band gap, and optical constants vary with W content in this alloy family. It's not a methodological advance, and it doesn't establish new physics. It deserves peer review rather than a desk reject, because the underlying data are reproducible and the trend is useful, but the authors need to fix the stability reference and either add configurational averaging or soften the claims accordingly. I'd send it to a referee with a note that those two issues need to be addressed.","headline":"Useful full-range DFT map of Mo1-xWxS2, but the stability claim uses the wrong reference and the per-composition configurational sampling is too thin to support the fitted bowing.","tokens_in":23079,"tokens_out":4493,"would_cite":false,"duration_ms":47270,"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":"Direct bandgap survives every MoS2-WS2 alloy mix, 1.696 to 1.858 eV","keywords":["Mo1-xWxS2","monolayer transition metal dichalcogenides","direct bandgap","bandgap engineering","density functional theory","optical absorption","bandgap bowing"],"falsifier":"Compute the electronic structure for many random configurations or special quasirandom structures at $x=0.5$: if any of them puts the valence band maximum or conduction band minimum away from K, or shifts the gap by more than the roughly 0.1 eV scale of the claimed trend, the central claim fails. Alternatively, photoluminescence measurements on an alloy series could test whether the emission stays direct and continuously tunable at every intermediate composition.","tokens_in":22108,"feed_emoji":"💡","tokens_out":7420,"duration_ms":71227,"temperature":0.7,"pith_summary":"The paper argues that mixing MoS2 and WS2 within a single monolayer produces a thermodynamically stable alloy whose direct bandgap survives across the entire composition range. The gap rises near-monotonically from 1.696 eV in MoS2 to 1.858 eV in WS2, with a small downward bowing of 0.1564 eV at mid-composition. Because the band edges stay at the K point, alloying does not sacrifice the direct-gap character that makes monolayers useful for light absorption and emission. The calculated optical spectra shift to higher energy as tungsten content increases, so the same alloy series offers continuous control over both the electronic gap and the optical response.","feed_headline":"Direct bandgap survives every MoS2-WS2 alloy mix, 1.696 to 1.858 eV","feed_subtitle":"Composition sets the bandgap and absorption edge, keeping tunable light emission in one monolayer.","key_machinery":"The central object is the series of $4\\times4\\times1$ supercell alloy models in which 1, 2, 4, 6, 8, 10, 12, or 14 of the 16 metal sites are filled with W atoms, giving configurations H1 through H14. The mechanism carrying the argument is the isoelectronic substitution of Mo by W: both are group-6 metals, so no charge carriers are added, the 2H lattice is preserved with bond lengths near 3.18 Å, and the band edges are controlled by metal $d$ states. As W $5d$ states replace Mo $4d$ states at the valence and conduction edges, the K-point gap widens, and the bowing parameter $b = 0.1564$ eV quantifies the small deviation from linear interpolation.","core_discovery":"Using density functional theory on $4\\times4\\times1$ supercells with progressively more W atoms on the Mo sublattice, the paper finds that every composition from pure MoS2 to pure WS2 keeps a direct gap at the K point. The Kohn-Sham gap increases from 1.696 to 1.858 eV in a near-monotonic way, fitted by $E_g(x) = (1-x)\\,1.696 + x\\,1.858 - 0.1564\\,x(1-x)$ eV. Formation energies are negative for every alloy, and bond lengths and lattice parameters barely change, so the 2H structure remains intact. The projected density of states shows W $5d$ states progressively replacing Mo $4d$ states at both band edges, which the paper identifies as the origin of the gap increase. RPA optical spectra up to 8 eV show the A/B absorption features blueshifting with W content, with a total absorption-edge tuning range of roughly 100 nm.","pith_inferences":["If configurational disorder is averaged over, the bowing parameter may differ from 0.1564 eV, and the near-monotonic gap might develop a spread comparable to the total tuning range; the paper's single-configuration data should be read as a baseline rather than a final distribution.","The reported enhancement of optical absorption at 6.25 and 12.5% W suggests a local Mo-W interaction that could be tested by computing the joint density of states for multiple low-W configurations.","Since PBE underestimates gaps and RPA omits excitons, the absolute peak energies are shifted; a GW-BSE calculation at a few compositions would show whether the roughly 100 nm tuning range survives many-body corrections."],"forward_implications":["A single monolayer can act as a wavelength-tunable emitter or detector, with the absorption onset moving across roughly 100 nm as $x$ goes from 0 to 1.","Because the gap remains direct at every composition, alloying should not suppress radiative recombination the way indirect-gap materials would, supporting light-emitting applications.","The refractive index, reflectance, and optical conductivity evolve predictably with composition, allowing optical component design without changing material family.","The small bowing of 0.1564 eV means simple linear interpolation is a good first approximation for device design at intermediate compositions.","Thermodynamic stability across the full range suggests these alloys can be synthesized as homogeneous monolayers rather than phase-separated domains."],"supporting_citations":[{"why":"CVD-grown monolayer alloy study that supplies the experimental photoluminescence tuning range used for validation.","marker":"[23]"},{"why":"Experimental feasibility and synthesis context for atomically substitutional TMDC alloying, plus high-accuracy DFT comparison for 2D materials.","marker":"[22]"},{"why":"Provides GW-BSE reference values and exciton-binding estimates used to interpret the optical spectra.","marker":"[20]"},{"why":"Experimental optical gap of monolayer MoS2 used to benchmark the calculated PBE gap.","marker":"[39]"},{"why":"Experimental exciton binding and optical gap data for monolayer WS2 used to benchmark the tungsten-rich end.","marker":"[33]"},{"why":"Effective band structure theory of random alloys that grounds the interpretation of bandgap bowing.","marker":"[38]"},{"why":"Many-body optical spectrum of MoS2 used to motivate why excitonic effects must be kept in mind when reading the RPA peaks.","marker":"[34]"}],"fun_headline_variants":["MoS2-WS2 alloy monolayers: direct gap tunes from 1.696 to 1.858 eV","All MoS2-WS2 monolayer alloys retain direct bandgap, tunable 1.696-1.858 eV","First-principles study: Mo1-xWxS2 monolayers show tunable direct gap","MoS2-WS2 alloy monolayers: bandgap continuously tunable, optics blue-shift"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Each composition is represented by one hand-picked arrangement of W atoms, and the paper assumes that arrangement behaves like the real random alloy well enough that configurational disorder would not change the direct-gap character, the gap trend, or the bowing.","fun_headline_variants_meta":{"raw":{"variants":["MoS2-WS2 alloy monolayers: direct gap tunes from 1.696 to 1.858 eV","All MoS2-WS2 monolayer alloys retain direct bandgap, tunable 1.696-1.858 eV","First-principles study: Mo1-xWxS2 monolayers show tunable direct gap","MoS2-WS2 alloy monolayers: bandgap continuously tunable, optics blue-shift"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000533,"raw_usage":{"total_tokens":2623,"prompt_tokens":1060,"completion_tokens":1563,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":676,"completion_tokens_details":{"reasoning_tokens":1452}},"tokens_in":676,"tokens_out":1563,"duration_ms":12591,"temperature":1.0,"reasoning_tokens":1452,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:40:23.513727+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the electronic structure for many random configurations or special quasirandom structures at $x=0.5$: if any of them puts the valence band maximum or conduction band minimum away from K, or shifts the gap by more than the roughly 0.1 eV scale of the claimed trend, the central claim fails. Alternatively, photoluminescence measurements on an alloy series could test whether the emission stays direct and continuously tunable at every intermediate composition.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Effective band structure theory of random alloys that grounds the interpretation of bandgap bowing."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"CVD-grown monolayer alloy study that supplies the experimental photoluminescence tuning range used for validation."},{"cited_title":"Ganatra and Q","cited_arxiv_id":null,"evidence_quote":"Experimental feasibility and synthesis context for atomically substitutional TMDC alloying, plus high-accuracy DFT comparison for 2D materials."},{"cited_title":"Komsa and A","cited_arxiv_id":null,"evidence_quote":"Provides GW-BSE reference values and exciton-binding estimates used to interpret the optical spectra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental optical gap of monolayer MoS2 used to benchmark the calculated PBE gap."},{"cited_title":"Grimme, Semiempirical gga-type density functional constructed with a long-range dispersion correction, Journal of Computational Chemistry 27, 1787 (2006)","cited_arxiv_id":null,"evidence_quote":"Experimental exciton binding and optical gap data for monolayer WS2 used to benchmark the tungsten-rich end."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Many-body optical spectrum of MoS2 used to motivate why excitonic effects must be kept in mind when reading the RPA peaks."}],"review_version":1}