REVIEW 3 major objections 5 minor 1 cited by
Nonlinear valley Hall effect in a bilayer transition metal dichalcogenide
T0 review · 3 major / 5 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read A centrosymmetric bilayer can host a valley Hall effect in the nonlinear regime.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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).
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Model calculations of a MX2 bilayer, Eq. (3)] 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.
- [Figs. 3(a) and 4(b)] 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.
- [Eq. (5) and model calculations] 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.
minor comments (5)
- [Near Eq. (3)] Typo: 'Possion ratio' should be 'Poisson ratio'.
- [Eq. (5)] 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.
- [Fig. 2(a)] 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.
- [Fig. 1(c) caption] 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'.
- [Introduction and first-principles section] 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.
Circularity Check
No significant circularity: the NVH conductivity is computed from independently fitted Hamiltonian parameters and cross-checked by first-principles DFT; no fitted input is renamed as a prediction.
full rationale
The derivation chain is self-contained: (i) symmetry analysis shows which nonlinear conductivity tensor components survive under P, T, C3z, C2x, and Mx; (ii) a k.p Hamiltonian (Eqs. 2-4) is constructed with monolayer and interlayer parameters taken from Refs. [22-24], which are external DFT-based band-structure parameters; (iii) the NVH conductivity is evaluated via the Berry connection polarizability dipole, Eqs. (5)-(7); and (iv) first-principles DFT calculations independently reproduce the qualitative and semiquantitative behavior (e.g., chi_xyy ~ 0.28 vs 0.33 from the model). No Hamiltonian parameter (v_F, Delta, lambda_c, lambda_v, t_cc, t_vv, gamma_1-3) is fitted to the target observable chi_NVH, and no equation defines chi_NVH in terms of itself. The valley-contrasting sign follows from the tau-dependent strain coupling and is verified by momentum-resolved BCP patterns rather than imposed by construction. The zero-Poisson-ratio approximation in the strain Hamiltonian is an acknowledged modeling limitation ('the strain terms are obtained by assuming a zero Possion ratio', Eq. 3), which may shift quantitative values under realistic uniaxial stress, but it is an accuracy caveat, not a circular input. Similarly, the nonmonotonic strain dependence above 4% in DFT versus the model is attributed to neglected higher-order q terms, an approximation caveat rather than a logical circularity. There is no load-bearing self-citation chain: Refs. [22-24] are independent prior works, and the first-principles calculation provides an external cross-check. Therefore the central claim is not forced by definition or by fitted inputs; it is an independent model-plus-DFT prediction.
Assumptions & free parameters
free parameters (8)
- vF (Fermi velocity)
- Δ (crystal-field band gap)
- δ (energy shift)
- λc (conduction-band SOC) =
-1.5 meV (MoS2 value λc,0)
- λv (valence-band SOC)
- tcc (interlayer conduction hopping) =
70.6 meV·Å (tcc,0)
- tvv (interlayer valence hopping)
- γ1, γ2, γ3 (strain coupling parameters)
assumptions (5)
- domain assumption The k·p Hamiltonian (Eqs. 2-4) with the given basis and parameters accurately describes the low-energy electronic structure of strained bilayer MoS2 near the K± valleys.
- domain assumption The nonlinear Hall conductivity is given by the BCP dipole formula in Eqs. 5-7, which assumes a relaxation-time approximation and the clean limit as in Ref. [19].
- ad hoc to paper Uniaxial strain is captured by linear terms in ε with a zero Poisson ratio (Eq. 3).
- domain assumption Interlayer coupling is limited to tcc and tvv hoppings in Eq. 4.
- standard math Fermi-Dirac distribution at zero temperature (Eq. 5).
Cite this review
Pith. "Pith review of Nonlinear valley Hall effect in a bilayer transition metal dichalcogenide." pith.science (2026). https://pith.science/paper/BRRBHI6L
@misc{pith2026241219502,
author = {Pith},
title = {Pith review of: Nonlinear valley Hall effect in a bilayer transition metal dichalcogenide},
year = {2026},
howpublished = {\url{https://pith.science/paper/BRRBHI6L}},
note = {Machine review of arXiv:2412.19502}
}
abstract
Valley-contrasting Hall transport conventionally relies on the inversion symmetry breaking in two-dimensional systems, which greatly limits the selection range of valley materials. In particular, while monolayer transition metal dichalcogenides have been widely utilized as a well-known class of valley materials in valleytronics, the centrosymmetric nature hinders the realization of valley-contrasting properties in the bilayer counterparts. Here, taking MoS$_{2}$ as an example, we discover valley-contrasting transport in bilayer transition metal dichalcogenides by exploring nonlinear transport regime. Using effective models and first-principles calculations, our work demonstrates that nonvanishing nonlinear valley Hall conductivities emerge in a uniaxially strained MoS$_{2}$ bilayer, owing to strain-induced band tilts of Dirac fermions. With the aid of small spin-orbit-coupling induced band splittings, the conduction bands generate much remarkable nonlinear valley Hall conductivity. Moreover, the nonlinear conductivities are highly tunable through modulating the strength and the direction of the strain, chemical potential, and interlayer gap. Our findings not only expands material choices for valleytronic applications, but also provides opportunities for designing advanced electronic devices that leverage nonlinear valley transports.
Figures
Forward citations
Cited by 1 Pith paper
-
Sliding-tuned Quantum Geometry in Moir\'e Systems: Nonlinear Hall Effect and Quantum Metric Control
Interlayer sliding in multi-twist moiré systems acts as a tuning knob for Berry curvature and quantum metric, enabling a sliding-driven nonlinear Hall effect and quantum-metric control for testing fractional Chern ins...
Reference graph
Works this paper leans on
-
[19]
K. Das, K. Ghorai, D. Culcer, and A. Agarwal, Phys. Rev. Lett. 132, 096302 (2024)
2024
-
[1]
D. Xiao, W. Yao, and Q. Niu, Phys. Rev. Lett. 99, 236809 (2007)
2007
-
[2]
J. R. Schaibley, H. Yu, G. Clark, P. Rivera, J. S. Ross, K. L. Seyler, W. Yao, and X. Xu, Nature Reviews Ma- terials 1, 16055 (2016)
work page 2016
-
[3]
W. Yao, D. Xiao, and Q. Niu, Phys. Rev. B 77, 235406 (2008)
2008
-
[4]
Xiao, G.-B
D. Xiao, G.-B. Liu, W. Feng, X. Xu, and W. Yao, Phys. Rev. Lett. 108, 196802 (2012)
2012
-
[5]
T. Cao, G. Wang, W. Han, H. Ye, C. Zhu, J. Shi, Q. Niu, P. Tan, E. Wang, B. Liu, and J. Feng, Nature Commu- nications 3, 887 (2012)
work page 2012
-
[6]
H. Zeng, J. Dai, W. Yao, D. Xiao, and X. Cui, Nature Nanotechnology 7, 490 (2012)
work page 2012
-
[7]
S. Wu, J. S. Ross, G.-B. Liu, G. Aivazian, A. Jones, Z. Fei, W. Zhu, D. Xiao, W. Yao, D. Cobden, and X. Xu, Nature Physics 9, 149 (2013)
work page 2013
Show all 30 references
-
[8]
M. Sui, G. Chen, L. Ma, W.-Y. Shan, D. Tian, K. Watan- abe, T. Taniguchi, X. Jin, W. Yao, D. Xiao, and Y. Zhang, Nature Physics 11, 1027 (2015)
2015
-
[9]
Shimazaki, M
Y. Shimazaki, M. Yamamoto, I. V. Borzenets, K. Watan- abe, T. Taniguchi, and S. Tarucha, Nature Physics 11, 1032 (2015)
2015
-
[10]
Ideue and Y
T. Ideue and Y. Iwasa, Annual Review of Condensed 6 Matter Physics 12, 201 (2021)
2021
-
[11]
Z. Z. Du, H.-Z. Lu, and X. C. Xie, Nature Reviews Physics 3, 744 (2021)
2021
-
[12]
Bandyopadhyay, N
A. Bandyopadhyay, N. B. Joseph, and A. Narayan, Ma- terials Today Electronics 8, 100101 (2024)
2024
-
[13]
Sodemann and L
I. Sodemann and L. Fu, Phys. Rev. Lett. 115, 216806 (2015)
2015
-
[14]
Z. Z. Du, C. M. Wang, H.-Z. Lu, and X. C. Xie, Phys. Rev. Lett. 121, 266601 (2018)
2018
-
[15]
Y. Gao, S. A. Yang, and Q. Niu, Phys. Rev. Lett. 112, 166601 (2014)
2014
-
[16]
C. Wang, Y. Gao, and D. Xiao, Phys. Rev. Lett. 127, 277201 (2021)
2021
-
[17]
H. Liu, J. Zhao, Y.-X. Huang, W. Wu, X.-L. Sheng, C. Xiao, and S. A. Yang, Phys. Rev. Lett. 127, 277202 (2021)
2021
-
[18]
J. Wang, H. Zeng, W. Duan, and H. Huang, Phys. Rev. Lett. 131, 056401 (2023)
2023
-
[20]
Strachan, A
J. Strachan, A. F. Masters, and T. Maschmeyer, ACS Appl. Energy Mater. 4, 7405 (2021), publisher: American Chemical Society
2021
-
[21]
Xiao, M.-C
D. Xiao, M.-C. Chang, and Q. Niu, Rev. Mod. Phys. 82, 1959 (2010)
2010
-
[22]
Liu, W.-Y
G.-B. Liu, W.-Y. Shan, Y. Yao, W. Yao, and D. Xiao, Phys. Rev. B 88, 085433 (2013)
2013
-
[23]
S. Fang, S. Carr, M. A. Cazalilla, and E. Kaxiras, Phys. Rev. B 98, 075106 (2018)
2018
-
[24]
Korm´ anyos, V
A. Korm´ anyos, V. Z´ olyomi, V. I. Fal’ko, and G. Burkard, Phys. Rev. B 98, 035408 (2018)
2018
-
[25]
K. Kang, T. Li, E. Sohn, J. Shan, and K. F. Mak, Nature Materials 18, 324 (2019)
2019
-
[26]
Wang and X
H. Wang and X. Qian, npj Computational Materials 5, 119 (2019)
2019
-
[27]
Q. Chen, Y. Wei, P.-B. Zhai, and Y.-J. Gong, Rare Met- als 43, 6125 (2024)
2024
-
[28]
Manzeli, D
S. Manzeli, D. Ovchinnikov, D. Pasquier, O. V. Yazyev, and A. Kis, Nature Reviews Materials 2, 17033 (2017)
2017
-
[29]
W. Choi, N. Choudhary, G. H. Han, J. Park, D. Akin- wande, and Y. H. Lee, Materials Today 20, 116 (2017)
2017
-
[30]
M. Wu, Y. Xiao, Y. Zeng, Y. Zhou, X. Zeng, L. Zhang, and W. Liao, InfoMat 3, 362 (2021)
2021
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.