{"id":"f25b9b94-d77b-4962-9551-42118c85f14a","arxiv_id":"2412.19502","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Uniaxial strain makes a centrosymmetric MoS2 bilayer exhibit a sizable nonlinear valley Hall conductivity, as shown by effective-model and first-principles calculations.","lead":"Bilayer molybdenum disulfide, usually too symmetric for valley Hall effects, develops a valley-dependent transverse current when squeezed along one direction. The work predicts a measurable nonlinear Hall signal from these centrosymmetric bilayers, expanding the materials available for valleytronic devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Zero-Poisson-ratio strain protocol in Eq. (3) supports the quantitative NVH claim; realistic uniaxial stress (ν≈0.22 in MoS2) could shift the quoted 0.33 peak and tuning curves, even though the symmetry-based existence is not at risk.","rationale":"The paper's central claim is the existence and size of a nonlinear valley Hall effect in strained bilayer MoS2, and the symmetry analysis provides a robust existence argument: uniaxial strain breaks C3z, P and T allow the BCP-dipole valley contrast, and C2x/Mx select χxyy while forbidding χyxx. The effective model and first-principles calculation agree at 2% strain (0.33 vs 0.28 in the same units), which is genuine independent support. The weakest load-bearing point is the quantitative reliability of the strain Hamiltonian. I agree with the reader that the zero-Poisson-ratio assumption is the key vulnerability: it directly controls the strain-induced tilts and therefore the size of the BCP dipole, and the DFT check does not fully remove the issue if the DFT calculation uses the same fixed-transverse-cell protocol. The high-order q truncation and the unavailable Supplementary Material are secondary; they affect high-strain behavior and derivability, but not the small-strain peak that anchors the headline number. A protocol-controlled model/DFT recomputation can settle the issue. Thus the reader's CONDITIONAL verdict is appropriate; I do not see a reason to move the verdict to ACCEPT or REJECT. If the proposed test shows the peak shift is small, the quantitative claim would be robust to the Poisson-ratio ambiguity and the manuscript could be strengthened by stating the strain protocol explicitly.","tokens_in":9633,"tokens_out":24760,"duration_ms":267279,"concrete_test":"Recompute the model curves in Fig. 2(a) and Fig. 3(a,b) with the realistic strain tensor for uniaxial stress: u_xx=ε(cos²θ−ν sin²θ), u_yy=ε(sin²θ−ν cos²θ), u_xy=ε(1+ν) sinθ cosθ, using ν=0.22 for MoS2, re-expressing the Eq. (3) strain terms as functions of u_ij (so γ3 is effectively multiplied by 1+ν and the trace terms by 1−ν). If the χxyy peak at 2% strain stays within 10–20% of 0.33 and remains near μ≈3 meV, the concern is minor. Independently, run a DFT calculation under uniaxial stress (relax the transverse cell to zero stress rather than fixing it) and compare the χxyy(μ) peak; a shift beyond about 20% would confirm that the quoted values are specific to the zero-Poisson-ratio protocol.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The single most load-bearing assumption is the zero-Poisson-ratio treatment of strain in Eq. (3). For a uniaxial strain along direction θ, the symmetry-allowed strain Hamiltonian is built from the pseudogauge combination u_xx−u_yy ∝ ε(1+ν)cos2θ and u_xy ∝ ε(1+ν)sin2θ, plus the trace u_xx+u_yy ∝ ε(1−ν). Setting ν=0, as stated in the text, fixes a specific strain protocol in which the transverse lattice constant is held fixed. Real suspended or substrate-strained MoS2 under uniaxial stress has ν≈0.22, so the effective γ3 coupling is roughly 22% larger and the γ1/γ2 terms change in the opposite direction through the trace. Because the reported NVH conductivity arises from small strain-induced band tilts whose BCP contribution is amplified by the tiny conduction-band SOC splitting (near-degenerate denominator in Eq. (7)), even a 20% change in the strain coupling can shift the peak magnitude and its chemical-potential location, quoted as 0.33 (e3/ℏ) Å meV−1 at μ=3 meV for 2% strain. The first-principles verification in Fig. 4 appears to use the same u_yy=0 strain protocol, so it does not independently validate the uniaxial-stress case. This concern is quantitative rather than existential: the symmetry analysis correctly guarantees a nonvanishing χNVH once C3z is broken, but the advertised numbers and the strain- and gap-tuning curves are protocol-dependent until the Poisson-ratio issue is resolved.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper predicts a nonlinear valley Hall (NVH) effect in centrosymmetric 2H-MoS2 bilayers under uniaxial strain. The authors first show by symmetry analysis that the simultaneous presence of inversion and time-reversal symmetry forbids the Berry curvature dipole but allows the Berry connection polarizability (BCP) to produce a nonlinear Hall current, and that C3z symmetry must be broken for the response to be nonvanishing. They then construct a k.p model with spin-orbit coupling and interlayer hopping, calculate the NVH conductivity, and find a large chi_xyy peak of 0.33 (e^3/hbar) Angstrom meV^-1 at 2% strain near the conduction-band edge, about 8.03 nm*mA/V^2. The paper demonstrates tunability of the response with strain strength, strain direction, chemical potential, SOC strength, and interlayer hopping, and corroborates the main results with first-principles calculations, which give a peak of 0.28 in the same units at the same strain.","tokens_in":1476,"tokens_out":1482,"duration_ms":78531,"significance":"The central symmetry argument is sound and the paper addresses an important gap: extending valley Hall physics to centrosymmetric bilayer systems. The combination of an effective model and first-principles calculations with reasonable agreement at small strain is a strength, as is the explicit identification of the BCP mechanism and the proposal of the conduction-band SOC splitting as a tunable knob and as a detectable signal. If the quantitative predictions hold for realistic experimental strain protocols, the work will be of considerable interest to the valleytronics community and could motivate experiments on strained bilayer TMDs using nonlocal resistance measurements.","major_comments":[{"comment":"The strain Hamiltonian in Eq. (3) is derived assuming a zero Poisson ratio, explicitly stated in the text. For real uniaxial stress in MoS2 with nu approximately 0.22, the pseudogauge strain combinations u_xx-u_yy and u_xy are multiplied by (1+nu) while the trace term u_xx+u_yy is multiplied by (1-nu). Because the reported NVH conductivity arises from small strain-induced band tilts and is amplified by the near-degenerate denominator in Eq. (7) associated with the small conduction-band SOC splitting, a roughly 20% change in the effective strain couplings can shift the advertised peak magnitude and its chemical-potential location, and can modify the tuning curves in Fig. 3(a). The first-principles calculation in Fig. 4 appears to use the same fixed-transverse-lattice protocol, so it does not independently validate the uniaxial-stress case. I request that the authors either incorporate a realistic Poisson ratio into the strain parameterization or clearly reframe the predictions as valid for a specific fixed-transverse-lattice-constant protocol, and discuss how the peak values and tuning curves would change under a typical experimental uniaxial-stress setup.","section":"Model calculations of a MX2 bilayer, Eq. (3)"},{"comment":"The model calculation in Fig. 3(a) predicts that chi_xyy increases monotonically with strain up to 8%, while the first-principles result in Fig. 4(b) shows a decrease for strain beyond 4%. The paper attributes this discrepancy to neglecting q-related high-order terms in the Hamiltonian. Because the abstract and text emphasize that the effect is highly tunable through modulating the strength of strain, this discrepancy directly limits the predictive power of the model for the strain-strength-dependence claim. The authors should explicitly restrict their strain-tunability claim to strains of about 4% or below, or extend the effective model to include the relevant higher-order q terms.","section":"Figs. 3(a) and 4(b)"},{"comment":"The integral in Eq. (5) is written over the full Brillouin zone, but the k.p Hamiltonian in Eq. (2) is valid only in a neighborhood of the K+/- valleys. The model calculation must be restricted to a finite integration window around each valley; otherwise the numerical value depends on an uncontrolled cutoff. Please specify the integration cutoff (or matching procedure) used to produce Figs. 2 and 3, and demonstrate that the reported peak values and the strain-direction dependence are insensitive to the chosen window size.","section":"Eq. (5) and model calculations"}],"minor_comments":[{"comment":"Typo: 'Possion ratio' should be 'Poisson ratio'.","section":"Near Eq. (3)"},{"comment":"The notation (-1)^tau inside the BZ integral is ambiguous because tau is a valley index, not a function of k. It would be clearer to write the conductivity as a sum over valleys, with the valley label attached to the Hamiltonian and the integral over each valley patch.","section":"Eq. (5)"},{"comment":"The axis label appears to list both chi_xyy and chi_yxx without subscripts; please format the labels cleanly so that the two curves are unambiguously identified.","section":"Fig. 2(a)"},{"comment":"The phrase 'q is given in units of the reciprocal of the lattice constant, a' should read 'q is given in units of 1/a'.","section":"Fig. 1(c) caption"},{"comment":"When comparing with the nonlinear Hall conductivity of few-layer WTe2, give the specific value used for the comparison and the reference from which it is taken.","section":"Introduction and first-principles section"}],"recommendation":"major_revision","confidential_remarks":"The paper makes a solid symmetry-based case for the NVH effect in centrosymmetric bilayers and provides a two-pronged model/DFT study. The main concern is that all quantitative predictions are tied to a zero-Poisson-ratio strain protocol, which is a special case rather than the generic experimental condition. This is fixable, but it requires either new calculations with nu approximately 0.22 or a careful reframing of the claims as protocol-specific. The paper fits the journal scope and the central existence claim appears correct."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a good paper and the central claim holds up. The authors show that uniaxial strain breaks C3z in centrosymmetric bilayer MoS2 and generates a nonlinear valley Hall effect from Berry connection polarizability, with opposite signs for the two valleys. The symmetry analysis is convincing and the effective model matches DFT near small strain (0.33 vs 0.28 in the same units). That is real, independent cross-checking.\n\nWhat is genuinely new: the BCP dipole mechanism is extended from strained graphene and organic monolayers (Ref. [19]) to bilayer TMDs. The extra ingredients — interlayer coupling, tunable van der Waals gap, and the small conduction-band SOC splitting — are not just decoration; the paper shows the conductivity is highly sensitive to the tiny conduction-band splitting, which is a useful way to detect small spin-orbit splittings, and that the van der Waals gap is a practical tuning knob. The strain-direction dependence (cos/sin 2θ) is a clean prediction.\n\nSoft spots, in order of importance. First, the quantitative predictions rest on the strain Hamiltonian in Eq. (3), which the text explicitly says is built assuming zero Poisson ratio. That is a real limitation. Real MoS2 under uniaxial stress has ν≈0.22, which changes the pseudogauge couplings by roughly 20% and shifts the quoted peak value and its chemical-potential location. The DFT verification in Fig. 4 uses the same u_yy=0 protocol, so it does not independently validate the uniaxial-stress case. This is a quantitative concern, not an existential one: the symmetry-based existence of the NVH effect does not depend on the Poisson ratio. But the advertised peak of 0.33 (e3/ℏ) Å meV−1 is protocol-specific until this is addressed. Second, the model neglects higher-order q terms, which the paper itself acknowledges is why the model and DFT diverge for strains above 4% (Fig. 4b). That is an honest, stated limitation. Third, the analytic derivation is in the Supplementary Material, which was not available for this review, and no code or data are released. Those are addressable, not fatal.\n\nThe paper is worth a serious referee. It is a solid theoretical contribution with a clear symmetry argument, model-DFT agreement at small strain, and experimentally falsifiable predictions. I would send it to review and ask the authors to discuss the Poisson-ratio dependence and release the SM and, ideally, the calculation code.\n\nWho is this for? Researchers in valleytronics and nonlinear Hall physics. It does not reshape the field, but it removes a symmetry constraint that was assumed necessary for valley-contrasting transport and opens a specific material class for experiments.","headline":"Solid symmetry-based prediction of strain-induced nonlinear valley Hall effect in bilayer MoS2, with quantitative numbers that are protocol-dependent because the strain model assumes zero Poisson ratio.","tokens_in":10582,"tokens_out":1421,"would_cite":true,"duration_ms":16649,"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 centrosymmetric bilayer can host a valley Hall effect in the nonlinear regime.","keywords":["nonlinear valley Hall effect","Berry connection polarizability","bilayer MoS2","uniaxial strain","valleytronics","transition metal dichalcogenides","second harmonic transport"],"falsifier":"A second-harmonic transport measurement on a Hall bar of bilayer MoS2 under 2% uniaxial strain along x, with the chemical potential tuned about 3 meV above the conduction band minimum, should show a transverse current quadratic in the applied field with magnitude near $0.3\\,(e^3/\\hbar)\\,\\text{Å}\\,\\text{meV}^{-1}$ and with opposite signs at $K_+$ and $K_-$, while $\\chi_{yxx}$ stays zero; a null transverse second-harmonic signal at that doping would contradict the central claim.","tokens_in":9408,"feed_emoji":"⚡","tokens_out":8455,"duration_ms":77124,"temperature":0.7,"pith_summary":"Bilayer molybdenum disulfide is centrosymmetric, so it shows no linear valley Hall effect, and valley physics in it has been considered blocked. This paper shows that the block disappears in the nonlinear transport regime: under uniaxial strain, the Dirac bands at the $K_+$ and $K_-$ valleys tilt, and the Berry connection polarizability generates opposite second-order transverse currents from the two valleys. The predicted nonlinear valley Hall conductivity reaches about $0.33\\,(e^3/\\hbar)\\,\\text{Å}\\,\\text{meV}^{-1}$ near the conduction band edge, which is much larger than in strained graphene and comparable to the organic conductor that has been studied for the same effect. The response is tunable by strain strength and direction, chemical potential, and the interlayer van der Waals gap, and first-principles calculations confirm the model's main features. If right, the paper makes uniaxially strained centrosymmetric TMD bilayers viable valleytronic materials despite having inversion symmetry.","feed_headline":"Bilayer MoS2 gains a nonlinear valley Hall effect under strain","feed_subtitle":"Uniaxial strain tilts the Dirac bands, yielding a tunable transverse E² current in a centrosymmetric bilayer.","key_machinery":"The central object is the Berry connection polarizability (BCP), defined per band as $$$G^{{\\alpha\\beta}}$_n(\\mathbf{k})=2\\,\\mathrm{Re}\\sum_{m\\neq n}\\frac{$A^{{nm}}$_\\$\\alpha$(\\mathbf{k})$A^{{mn}}$_\\$\\beta$(\\mathbf{k})}{\\varepsilon_{n\\mathbf{k}}-\\varepsilon_{m\\mathbf{k}}},$$ with $A^{nm}_\\alpha=\\langle u_{n\\mathbf{k}}|i\\partial_\\alpha|u_{m\\mathbf{k}}\\rangle$, combined with the BCP dipole $\\Lambda_{\\alpha\\beta\\gamma}=\\sum_n \\lambda^{\\alpha\\beta\\gamma}_n\\,\\partial f/\\partial \\varepsilon$, where $\\lambda$ involves the band velocity times $G$. The BCP is the geometric quantity that survives when both inversion and time-reversal symmetry are present, unlike the Berry curvature and Berry curvature dipole. In the model, the strain term $\\epsilon[\\gamma_1\\sigma_0+\\gamma_2\\sigma_z+\\gamma_3(\\cos 2\\theta\\,\\sigma_x-\\tau\\eta\\sin2\\theta\\,\\sigma_y)]$ tilts the low-energy Dirac bands, and the small conduction-band spin-orbit splitting $\\lambda_c$ makes the energy denominator in $G$ small, strongly amplifying the response. The machinery converts an applied field squared into a valley-antisymmetric transverse current.","core_discovery":"On its own terms, the paper's discovery is that a $\\mathcal P$- and $\\mathcal T$-symmetric bilayer can still have a valley Hall effect when the current is quadratic in the applied field. Symmetry analysis shows that $\\mathcal P$ and $\\mathcal T$ together forbid the Berry curvature dipole but not the Berry connection polarizability dipole, so the only requirement is breaking the threefold rotation $C_{3z}$. A uniaxial strain does this and tilts the Dirac cones; the tilt makes the BCP dipole $\\Lambda_{xyy}$ mirror-symmetric at each valley and opposite between valleys, giving a nonzero $\\chi^{\\mathrm{NVH}}_{xyy}$ while $\\chi^{\\mathrm{NVH}}_{yxx}$ vanishes for strain along $x$. Because the conduction-band spin-orbit splitting is only a few meV, the BCP around the conduction band edge is large, and the conductivity peaks at $0.33\\,(e^3/\\hbar)\\,\\text{Å}\\,\\text{meV}^{-1}$ at $\\mu=3$ meV. First-principles DFT gives the same qualitative behavior with a peak value of $0.28\\,(e^3/\\hbar)\\,\\text{Å}\\,\\text{meV}^{-1}$ at 2% strain, supporting the claim.","pith_inferences":["Inference: because the symmetry argument depends only on $\\mathcal P$, $\\mathcal T$, and broken $C_{3z}$, the same mechanism should appear in other centrosymmetric 2D bilayers with gapped Dirac-like bands, not just MoS2, whenever uniaxial strain creates band tilts.","Inference: the $\\cos 2\\theta$/ $\\sin 2\\theta$ strain-direction dependence suggests a simple device control—rotating the strain axis by 45 degrees transfers the transverse response from $\\chi_{xyy}$ to $\\chi_{yxx}$, effectively rotating the Hall current direction without reversing the field.","Inference: since the BCP grows as the inverse cube of the interband splitting, materials with even smaller conduction-band spin-orbit splitting than MoS2 should show substantially larger nonlinear valley Hall conductivities; the paper does not calculate this, but it follows from Eq. (7)."],"forward_implications":["A measurable second-harmonic transverse voltage should appear in uniaxially strained bilayer MoS2, detectable through nonlocal resistance measurements.","The predicted peak value is about 0.33 (e3/ℏ) Å meV−1 at a chemical potential of 3 meV, far exceeding the nonlinear valley Hall response of strained graphene and comparable to organic-conductor values.","Strain strength, strain direction, carrier doping, and the van der Waals gap each provide an independent tuning knob for the magnitude and sign of the nonlinear valley Hall current.","The effect is naturally largest near small conduction-band spin-orbit splittings, so the nonlinear valley Hall response can serve as a detector of those splittings.","First-principles calculations reproduce the model's sign structure and peak value at 2% strain, indicating the effect is not an artifact of the low-energy model."],"supporting_citations":[{"why":"Introduces the Berry curvature dipole mechanism and the nonlinear Hall effect formalism that this paper extends to the valley-contrasting BCP case.","marker":"[13]"},{"why":"Shows that nonlinear valley Hall responses can appear in centrosymmetric monolayers, the direct precedent for the bilayer result here.","marker":"[19]"},{"why":"One of the two sources for the monolayer k·p Hamiltonian and spin-orbit parameters used to build the bilayer model in Eq. (3).","marker":"[22]"},{"why":"The other source for the monolayer k·p Hamiltonian and strain/band parameters entering Eq. (3).","marker":"[23]"},{"why":"Supplies the interlayer coupling Hamiltonian in Eq. (4) that mixes the two monolayers.","marker":"[24]"},{"why":"Experimental observation of the nonlinear Hall effect in few-layer WTe2, used as the benchmark that the predicted conductivity surpasses.","marker":"[25]"},{"why":"Computational study of WTe2 nonlinear Hall conductivity, used as another benchmark comparison for the predicted magnitude.","marker":"[26]"}],"fun_headline_variants":["Strain unlocks nonlinear valley Hall effect in bilayer MoS2","Nonlinear valley Hall current appears in strained bilayer MoS2","Strained bilayer MoS2 shows quadratic valley Hall response","Bilayer MoS2 gets nonlinear valley Hall under uniaxial strain","Nonlinear Hall effect in bilayer MoS2 via strain-induced tilt"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The size of the predicted conductivity rests on a strain Hamiltonian derived by setting the Poisson ratio to zero, while the real material contracts in the perpendicular direction under uniaxial strain; that approximation could shift the tilts and the numerical value of the conductivity even though the symmetry-based existence of the effect would remain.","fun_headline_variants_meta":{"raw":{"variants":["Strain unlocks nonlinear valley Hall effect in bilayer MoS2","Nonlinear valley Hall current appears in strained bilayer MoS2","Strained bilayer MoS2 shows quadratic valley Hall response","Bilayer MoS2 gets nonlinear valley Hall under uniaxial strain","Nonlinear Hall effect in bilayer MoS2 via strain-induced tilt"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000505,"raw_usage":{"total_tokens":2503,"prompt_tokens":1022,"completion_tokens":1481,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":1393}},"tokens_in":638,"tokens_out":1481,"duration_ms":10292,"temperature":1.0,"reasoning_tokens":1393,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T00:16:12.015003+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A second-harmonic transport measurement on a Hall bar of bilayer MoS2 under 2% uniaxial strain along x, with the chemical potential tuned about 3 meV above the conduction band minimum, should show a transverse current quadratic in the applied field with magnitude near $0.3\\,(e^3/\\hbar)\\,\\text{Å}\\,\\text{meV}^{-1}$ and with opposite signs at $K_+$ and $K_-$, while $\\chi_{yxx}$ stays zero; a null transverse second-harmonic signal at that doping would contradict the central claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The other source for the monolayer k·p Hamiltonian and strain/band parameters entering Eq. (3)."},{"cited_title":"Korm´ anyos, V","cited_arxiv_id":null,"evidence_quote":"Supplies the interlayer coupling Hamiltonian in Eq. (4) that mixes the two monolayers."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental observation of the nonlinear Hall effect in few-layer WTe2, used as the benchmark that the predicted conductivity surpasses."},{"cited_title":"Wang and X","cited_arxiv_id":null,"evidence_quote":"Computational study of WTe2 nonlinear Hall conductivity, used as another benchmark comparison for the predicted magnitude."}],"review_version":1}