{"id":"acc1fa1d-8790-4602-bbca-62bfb0d39495","arxiv_id":"2607.29084","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Intrinsic tensile strain from the asymmetric Janus stack selects a 2×2 Nb-trimerized charge-density wave in NbSSiAs2; 1% compression makes one configuration Z2-topological without destroying ~6 K superconductivity.","lead":"First-principles calculations predict that the Janus monolayer NbSSiAs2—NbS2 with one sulfur layer replaced by an As–Si–As stack—develops a 2×2 charge-density-wave pattern selected by an internal 7.35% tensile strain, and that 1% external compression turns one pattern into a topological, still-superconducting state. Why read it: it proposes a built-in structural knob, rather than external strain, for choosing among competing electronic orders in 2D materials.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'primarily strain' attribution is underdetermined: the two reference models do not separate strain from chemistry in the actual Janus substitution, and the NbS2 baseline is protocol-dependent.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the two clamped reference models assume that in-plane lattice constant, varied independently, captures the strain channel, and that the NbS2 baseline is robust. My analysis adds specificity about why the decomposition is underdetermined: in the actual Janus substitution, strain and chemistry are generated by the same structural change, so the two reference calculations are projections of a single coupled perturbation rather than independent physical knobs. The fixed-lattice Janus result shows chemistry alone shifts q_min from 0.77 to ~0.67 ΓM, demonstrating that the channels are not cleanly separable. The acknowledged protocol dependence of the NbS2 off-M position further weakens the 'redirection' framing. These points support the reader's CONDITIONAL verdict rather than ACCEPT. They do not refute the paper's main computational findings—the 2×2 CDW, the bistability, the topological transition, or the superconductivity estimates—because those are less dependent on the causal attribution. The proposed test—a lattice-constant sweep for the Janus composition plus a PBEsol cross-check of pristine NbS2—would directly test whether strain is necessary and sufficient for the M-point selection. No change to the reader's verdict is required; it should remain CONDITIONAL pending such a test.","tokens_in":12871,"tokens_out":4309,"duration_ms":50857,"concrete_test":"Compute the soft-mode minimum q-vector q_min for the Janus NbSSiAs2 composition with in-plane lattice constants swept from 3.346 Å to 3.592 Å in steps of ~0.05 Å, fully relaxing internal coordinates at each a, using the same QE/PBE protocol. Additionally recompute the pristine NbS2 phonon dispersion with PBEsol. If q_min(a) for Janus reaches M only near the relaxed lattice constant, the intrinsic-strain attribution is supported; if q_min stays near 0.67ΓM until chemistry is changed, or if PBEsol puts pristine NbS2's minimum at M, then the claim that strain is the primary selector does not hold.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central causal claim—that Janus-induced intrinsic tensile strain redirects the parent soft-mode instability to M and thereby selects the 2×2 CDW—rests on the two constrained reference structures in Figs. 1d–1e. In the strain-only model, NbS2 is clamped to the Janus lattice constant (a = 3.592 Å) and the soft mode moves to M; in the fixed-lattice Janus model, the As-Si-As substitution is made at the NbS2 lattice constant and the soft mode sits at ~0.67 ΓM. From these two points the paper concludes that strain, not chemistry/polarity, is the primary selector. This inference requires that the actual Janus substitution decomposes additively into an independent lattice-constant channel and a chemistry channel, but no such decomposition is demonstrated. The 7.35% lattice expansion is not an external control parameter in the real material; it is a consequence of the same substitution whose 'chemistry' is supposedly separated out. The fixed-lattice Janus result itself shows that chemistry alone moves the instability from ~0.77 to ~0.67 ΓM, so the two channels are not orthogonal. Moreover, the paper concedes (text near Fig. 1c) that the pristine NbS2 off-M position 'can vary with computational details' [66,67]. If another same-protocol functional places the NbS2 minimum at M, the redirection framing loses its baseline; the claimed shift from off-M to M is then an artifact of the chosen baseline, not a robust physical effect. The cDFPT result in Fig. 2b shows that Nb-d screening is necessary for the M instability, but it does not attribute that screening to the strain channel; it is equally compatible with chemistry-driven changes in the low-energy electronic structure. Thus the load-bearing assertion—that intrinsic strain, rather than the polar/chemical asymmetry, is what selects the 2×2 order—is underdetermined by the presented two-point decomposition.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports a first-principles study of monolayer Janus NbSSiAs2, derived from NbS2 by replacing one S layer with an As-Si-As trilayer. The authors find that the 7.35% in-plane lattice expansion relative to NbS2 displaces the parent soft-phonon instability from an off-M position (~0.77 Gamma-M) to the M point, selecting a 2x2 CDW reconstruction. Using constrained DFT perturbation theory they trace the M-point softening to momentum-selective electron-phonon coupling of Nb-d states, and they identify two nearly degenerate 1+3 Nb-trimerized 2x2 configurations separated by 0.56 meV/f.u. External biaxial strain is shown to reverse their relative stability; at 1% compression the 1+3-hollow configuration becomes a dynamically stable Z2-nontrivial state with gapless edge states, while Allen-Dynes estimates give Tc ~6.2 K. Both unstrained CDW phases are predicted to superconduct at ~6.6-7.1 K. The paper includes phonon stability, 300 K AIMD, three exchange-correlation functionals for the energy landscape, and Wilson-loop/topological edge-state checks.","tokens_in":13132,"tokens_out":8135,"duration_ms":81691,"significance":"If correct, the paper would establish a clean internal mechanism - Janus-induced intrinsic strain as a structural selector of CDW wavevector - and a concrete material platform in which CDW order, topology, and superconductivity coexist. The technical execution is generally careful: same-protocol comparisons, a cDFPT control that removes the Nb-d screening channel and eliminates the M-point instability, a NEB double-well profile, phonon and AIMD stability checks, and Wannier-based topology diagnostics. These controlled tests make the microscopic origin of the instability and the metastability of the two CDW configurations credible. The main caveat is that the central strain-versus-chemistry attribution is not as directly demonstrated as the rest of the paper, and the baseline on which it rests is acknowledged to be protocol-sensitive. With additional decomposition calculations or moderated claims, the paper would be a valuable contribution.","major_comments":[{"comment":"The central causal attribution is underdetermined. The two constrained reference models (Figs. 1d and 1e) are read as separating an 'intrinsic strain' channel from a 'chemical/polarity' channel, but the 7.35% lattice expansion in the real Janus structure is itself a consequence of the As-Si-As substitution. Clamping NbS2 to the Janus lattice constant therefore includes the chemical effect on the lattice, while the fixed-lattice Janus model at the NbS2 lattice constant still contains the substitution and its internal relaxations. Moreover, Fig. 1e shows that the substitution alone moves the soft-mode minimum from ~0.77 to ~0.67 Gamma-M, so the two channels are not orthogonal. The conclusion 'intrinsic strain primarily controls' is an inference, not a decomposition. A quantitative test is needed: e.g., Janus phonons at several intermediate in-plane lattice constants, or a force-constant de","section":"Intrinsic-Strain-Driven CDW Instability (Figs. 1b-1e)"},{"comment":"The redirection narrative depends on a protocol-specific baseline. The text near Fig. 1c concedes that the pristine NbS2 off-M position 'can vary with computational details' [66,67]. Since only one functional (PBE) is used for the phonon wavevector selection, the claimed shift from off-M to M may reflect the chosen baseline rather than a physical redirection. The authors should verify the NbS2 and Janus soft-mode positions with at least one additional same-protocol functional or show convergence with respect to pseudopotentials and q-grid; otherwise the conclusion should be framed as a commensuration/stabilization effect relative to this PBE baseline. This is load-bearing because the M-point selection is the starting point for the entire CDW phase diagram.","section":"Intrinsic-Strain-Driven CDW Instability (text near Fig. 1c)"}],"minor_comments":[{"comment":"The second data row is labeled '1+3-filled' but from the text and Fig. 5 it should be the unstrained '1+3-hollow' configuration (Tc = 7.09 K). Please correct the label.","section":"Table 1"},{"comment":"The notation q_CDW = Gamma-M is confusing: Gamma-M denotes a line, while the instability is at the M point. Use q_CDW = M or state explicitly 'along Gamma-M'.","section":"Figs. 1 and 2"},{"comment":"The functional labeled 'PZ' in Figs. 3c and 3d is not defined. Please name the exchange-correlation functional and cite it.","section":"Fig. 3 caption"},{"comment":"In the Methods paragraph, the phrase 'see also refs. [60,61] therein' is unclear. Please specify what the Supporting Information contains.","section":"Method"}],"recommendation":"major_revision","confidential_remarks":"The paper has strong technical controls and a plausible central claim, but the strain-versus-chemistry attribution needs additional evidence before publication. The Table 1 row label is a simple fix. I do not see grounds for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new thing here is NbSSiAs₂, a Janus monolayer with a computed 2×2 CDW selected at the M point, two nearly degenerate Nb-trimerized configurations, strain-tunable switching between them, and a claim that 1% compression produces a Z₂-nontrivial state that still superconducts near 6–7 K in Allen–Dynes estimates. If the material can be made, that is a nice platform for coexisting order. The authors do the work properly: same-protocol comparison against pristine NbS₂, a strain-only reference model, a fixed-lattice Janus model, cDFPT showing the M-point instability dies when Nb-d screening is removed, NEB double-well profile, three functionals, phonon stability and 300 K AIMD for both CDW phases, and Wilson-loop plus edge-state checks for topology. That is credible, careful first-principles work.\n\nThe soft spot is the causal attribution. The paper concludes that intrinsic strain, not chemistry or polarity, is the primary selector of the M-point CDW. That relies on the two constrained reference models. The strain-only NbS₂ model does reproduce the M point, but the fixed-lattice Janus model also moves the instability from ~0.77 to ~0.67 ΓM, so the chemical substitution is not silently sitting out. Because the lattice expansion is itself a consequence of the substitution, the decomposition is a constructed control rather than a direct measurement in the real material. Worse, the pristine NbS₂ off-M baseline is protocol-dependent, as the authors admit; if that baseline moved to M under another functional, the 'redirection' framing would lose its anchor. This does not sink the paper—the biaxial-strain scan on the actual material is direct evidence for strain tunability—but it does mean the 'intrinsic strain as design principle' claim should be softened.\n\nSmaller issues: the 0.56 meV/f.u. energy difference and strain crossover sit at DFT precision limits; the superconductivity numbers are estimates for a hypothetical material with μ*=0.10; and Table 1 mislabels a row (the 7.09 K entry is assigned to 1+3-filled instead of 1+3-hollow) and cites inconsistent Tc values between text and table. These are fixable.\n\nAll in, the paper deserves a serious referee. The central prediction of a bistable 2×2 CDW in this Janus monolayer is well supported; the mechanism attribution needs more careful treatment of the strain-chemistry separation, ideally with an additional functional for the phonon baseline or a discussion of why the constrained models are valid. I'd send it to review and ask for those clarifications, not a rejection. I'd also bring it to a reading group and would cite it if I were working on CDWs in 2D materials.","headline":"Solid computational prediction of a strain-selected 2x2 CDW in a new Janus monolayer, with a plausible but not fully resolved strain-vs-chemistry story; worth refereeing.","tokens_in":13862,"tokens_out":3934,"would_cite":true,"duration_ms":35068,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A monolayer's built-in strain selects a 2×2 charge-density-wave order, and 1% compression turns it topological while superconductivity survives.","keywords":["two-dimensional Janus materials","charge density wave","intrinsic strain","NbSSiAs2","NbS2","electron-phonon coupling","Z2 topology","superconductivity"],"falsifier":"A temperature-dependent electron diffraction or STM experiment on free-standing or substrate-supported monolayer NbSSiAs2 that finds no 2×2 CDW superlattice, or finds an ordering wave vector away from M, would falsify the selection claim. A computational check would be to recompute the strained-NbS2 reference phonons with several functionals: if the M-point minimum disappears with a different functional while the pristine NbS2 off-M instability remains, the redirection conclusion would not be robust.","tokens_in":12620,"feed_emoji":"⚛️","tokens_out":5712,"duration_ms":54194,"temperature":0.7,"pith_summary":"The paper argues that Janus asymmetry in monolayer NbSSiAs2 acts as an internal, built-in tensile strain of about 7.35%, and that this intrinsic strain — not the chemical polarity of the asymmetric layers — is the main control that shifts the parent NbS2 soft phonon from an off-M, incommensurate wave vector to the M point. That wave-vector selection makes a commensurate 2×2 CDW reconstruction the ground state, 31 meV/f.u. below the competing 1×2 stripe. Within the 2×2 phase, two nearly degenerate Nb-trimerized configurations form a double-well energy landscape, and external biaxial strain reverses which one is lower. Under 1% compressive strain, the 1+3-hollow configuration becomes Z2-nontrivial with gapless edge states while retaining phonon-mediated superconductivity with Tc ≈ 6.2 K. If correct, this gives a design rule: internal Janus strain can select CDW order, and strain can push the same lattice into a coexisting topological superconducting CDW state.","feed_headline":"A monolayer's built-in strain picks a 2×2 charge-density-wave order","feed_subtitle":"NbSSiAs2 inherits NbS2's soft phonon, but its internal tension pins it at the M point — and 1% compression makes the CDW topological.","key_machinery":"The load-bearing mechanism is the Janus-induced intrinsic tensile strain: replacing one S layer with an As–Si–As trilayer expands the in-plane lattice constant from 3.346 Å to 3.592 Å (7.35%). Two constrained reference calculations isolate the strain channel: strained NbS2 at the Janus lattice constant (which reproduces the M-point instability) and fixed-lattice Janus NbSSiAs2 at the NbS2 lattice constant (which only broadens the branch). The electronic origin is identified by constrained density functional perturbation theory and momentum-resolved fluctuation diagnostics of the phonon self-energy, which resolve the k-resolved contributions of the electron-phonon vertex g² and bare susceptib","core_discovery":"The central discovery is that intrinsic strain from Janus structural asymmetry redirects the CDW instability. In pristine NbS2, the lowest soft mode sits at q ≈ 0.77 ΓM, pointing to an incommensurate near-3×3 distortion. In NbSSiAs2, the same soft mode appears at the M point (q = ΓM), and constrained calculations show that clamping NbS2 to the Janus lattice constant reproduces this M-point selection, while keeping the Janus chemistry but restoring the NbS2 lattice constant shifts the minimum only to ≈0.67 ΓM. The authors conclude that the 7.35% in-plane expansion is the dominant structural degree of freedom selecting a commensurate 2×2 order. cDFPT then shows that the M-point softening disap","pith_inferences":["If lattice constant is the dominant selector, then other Janus or intercalation variants of NbS2 with different in-plane expansion can be screened by computing only the strain response of the parent soft mode, without full phonon calculations of each chemistry.","A direct test would be to grow monolayer NbSSiAs2 on substrates with different lattice constants: a substrate that compresses the lattice by ~1% should flip the CDW configuration from 1+3-filled to 1+3-hollow and, if the topological prediction is right, switch on edge conductance.","Uniaxial strain or biaxial strain beyond the ±3% window might stabilize the competing 1×2 stripe phase or an incommensurate state, extending the phase diagram the paper sketches for biaxial strain.","The momentum-selective electron-phonon coupling mechanism suggests that selective doping of the Nb-d bands, rather than strain, might also move the ordering vector, giving a complementary electronic route to CDW phase selection."],"forward_implications":["If the strain-redirection claim is right, the CDW ordering vector of Janus NbSSiAs2 should be exactly at M, giving a commensurate 2×2 superlattice observable by STM or electron diffraction.","The near degeneracy of the two 1+3 configurations means small external biaxial strain can switch between filled and hollow trimers, effectively acting as a switch on the CDW texture.","Compressive strain drives the 1+3-hollow configuration through a band inversion, yielding a Z2=1 phase with gapless edge states; because Tc remains ≈6.2 K, the same material is a platform where CDW, topology, and superconductivity coexist.","Since both CDW configurations retain Tc on the 6–7 K scale, CDW order here does not destroy superconductivity, unlike cases where CDW gaps the whole Fermi surface.","The functional-dependence checks (PBE, PBEsol, PZ) all show strain reversal of ΔE, so the bistability and strain response are not artifacts of one exchange-correlation functional."],"fun_headline_variants":["Janus strain locks in 2×2 CDW order in NbSSiAs2","Built-in strain selects CDW phase, compression adds topology","Strain picks CDW in NbSSiAs2; compression adds topology","Compression makes CDW topological, preserving superconductivity","Janus strain redirects soft mode, selecting 2×2 CDW"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The central claim rests on the assumption that the two constrained reference structures cleanly separate intrinsic strain from chemical/polarity effects, and that the parent NbS2 soft-mode position is computed reliably enough to define the 'off-M' baseline.","fun_headline_variants_meta":{"raw":{"variants":["Janus strain locks in 2×2 CDW order in NbSSiAs2","Built-in strain selects CDW phase, compression adds topology","Strain picks CDW in NbSSiAs2; compression adds topology","Compression makes CDW topological, preserving superconductivity","Janus strain redirects soft mode, selecting 2×2 CDW"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001239,"raw_usage":{"total_tokens":4942,"prompt_tokens":783,"completion_tokens":4159,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":527,"completion_tokens_details":{"reasoning_tokens":4064}},"tokens_in":527,"tokens_out":4159,"duration_ms":28899,"temperature":1.0,"reasoning_tokens":4064,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T13:56:06.574885+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A temperature-dependent electron diffraction or STM experiment on free-standing or substrate-supported monolayer NbSSiAs2 that finds no 2×2 CDW superlattice, or finds an ordering wave vector away from M, would falsify the selection claim. A computational check would be to recompute the strained-NbS2 reference phonons with several functionals: if the M-point minimum disappears with a different functional while the pristine NbS2 off-M instability remains, the redirection conclusion would not be robust.","supporting_citations":[],"review_version":1}