Pith. sign in

REVIEW 5 major objections 4 minor 32 references

Winding Berry dipole on uniaxially strained graphene/hBN/hBN moir\'e trilayers

T0 review · 5 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read In a graphene/hBN/hBN moiré trilayer, the local Berry dipole winds through the moiré and points along only six symmetry-allowed directions, opening a route to longitudinal nonlinear currents driven by quantum geometry alone, even in…

desk verdict A plausible but overclaimed proposal: the winding Berry dipole map is new and worth discussing, but the six-direction quantization is not established and the paper needs more validation before the central claim can be accepted. read the letter →

arxiv 2412.10584 v1 pith:IPGYEGU6 submitted 2024-12-13 cond-mat.mes-hall

Angiolo Huaman , Salvador Barraza-Lopez This is my paper · ORCID
classification cond-mat.mes-hall
keywords Berrycurvaturedipolemoirégraphene/hBNnonlinearHalleffectsemiclassicalapproximationtight-bindingwindingvectorfieldsecond-ordercurrent
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that the moiré pattern formed by placing graphene on a bilayer of hexagonal boron nitride (hBN) turns the local Berry curvature dipole D into a winding vector field. Ab initio calculations of the electrostatic potential from the substrate are folded into a tight-binding model of graphene's π electrons, and a semiclassical prescription assigns a local D(r) to each registry between the two lattices. The central result is that D(r) points along only six symmetry-allowed directions and winds around the high-symmetry ABA and ABC stackings. If correct, this means the direction of second-order nonlinear currents can be steered by the moiré texture alone, including longitudinal currents that normally require impurity scattering, even in ultraclean samples.

What carries the argument

The central object is the local Berry dipole D(r), the Fermi-weighted integral of k-derivatives of the Berry curvature over the Brillouin zone, evaluated from a semiclassically constructed tight-binding Hamiltonian of graphene under the substrate potential. The moiré is treated as a parameter space of local registries: for each registry, the potential ΔV(r) renormalizes on-site energies and hopping integrals, giving a local unit-cell Hamiltonian whose Berry curvature yields D. The winding of D(r) follows from how the relative mass (εA−εB) and the pattern of largest hoppings vary across the moiré, with mirror symmetry enforcing Dx=0 for the high-symmetry stackings.

What would settle it

Compute the Berry curvature dipole directly from Bloch states of the full (54,55) moiré supercell—without the local-registry decomposition—and compare the resulting direction and magnitude with the winding map in Fig. 6(d). If the full-supercell D does not point along the six predicted directions or does not wind around ABA/ABC, the semiclassical claim fails. A complementary experimental falsifier would be a second-harmonic transport measurement that finds current directions outside the six allowed angles.

Watch

Extended reading notes

Core claim

The authors set out to show that the moiré between graphene and an AB-stacked hBN bilayer provides a mechanism for generating second-order nonlinear currents that is purely quantum-geometric. Starting from the ab initio potential difference ΔV(r) produced by the substrate, they compute renormalized on-site energies εA, εB and first-nearest-neighbor hopping integrals tj for each local registry of the moiré. The resulting local Hamiltonian yields a Berry dipole D whose direction is constrained by the trigonal symmetry of the lattice to six angles (±90°, ±30°, ±150° relative to the x-axis), and whose map across the (54,55) moiré supercell winds around the ABA and ABC stacking points. They conclude that moiré textures are a knob for homogenizing and scattering nonlinear currents, opening the possibility of longitudinal second-order currents based solely on the Berry dipole.

Load-bearing premise

The entire construction treats the moiré as if each local patch has its own periodic Hamiltonian, assuming the local registry persists over a few lattice constants; if that semiclassical picture breaks down, the winding map of D would not represent the actual nonlinear response.

Editorial extensions

If this is right

  • The local direction of the Berry dipole—and hence the direction of the second-order nonlinear Hall current—is set by the moiré registry and can be tuned by stacking configuration and strain.
  • Longitudinal second-order currents, normally attributed to skew scattering and side jumps, can in principle arise from the Berry dipole alone when its direction tilts within the moiré, even in ultraclean samples.
  • The winding of D(r) around ABA and ABC stacking points makes the moiré a topological texture for nonlinear transport, analogous to winding polar fields in other moiré systems.
  • Only six directions of D are symmetry-allowed, so the current response is quantized in direction rather than continuously variable.
  • The semiclassical local-registry Hamiltonian provides a parameter-space method to compute D for any large moiré supercell without diagonalizing the full superlattice.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A direct test would be a full moiré-supercell calculation of the Berry curvature dipole without the semiclassical local-registry approximation; if the winding pattern survives, the semiclassical map is validated.
  • The six-direction constraint likely generalizes to any moiré with trigonal lattice symmetry, so other graphene/substrate systems may show the same quantized winding rather than continuous rotation.
  • If the winding is robust, it could be probed by second-harmonic generation or photocurrent imaging: the local current direction should rotate by 60° steps across the moiré, which scanning photocurrent microscopy might resolve.
  • The authors' sampling of 82 registries is sparse; a denser map (or machine-learning interpolation over the parameter space) would reveal whether the winding number around ABA/ABC is +1 or −1 and whether domain walls carry the winding.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

5 major / 4 minor

Summary. The paper studies the Berry curvature dipole D in a moiré formed by graphene on an AB-stacked hBN homobilayer under uniaxial strain. The authors compute the moiré-induced local potential ΔV(r) from DFT, use it to renormalize the on-site energies and nearest-neighbor hoppings of a graphene tight-binding model via Eq. (5), and then evaluate D from Eq. (3) for different local registries. For the highly symmetric ABA and ABC stackings, D is purely transverse. The paper's central claim is that at generic moiré positions, D tilts and takes only six symmetry-allowed directions, leading to a winding vector field D(r) across the moiré. This is presented in Figs. 5 and 6(d), and the authors argue that such winding could enable longitudinal as well as transverse second-order currents solely from the Berry dipole.

Significance. If the six-direction winding claim is correct, the work introduces a concrete, materials-realistic mechanism by which moiré texture can control the local direction of the nonlinear Hall response, potentially enabling longitudinal second-order currents without extrinsic scattering. The DFT-to-tight-binding pipeline is standard, the assumption of a local unit-cell Hamiltonian is stated explicitly, and the prediction is falsifiable by direct computation or experiment. These are genuine strengths. However, the load-bearing evidence for the quantization and winding is currently incomplete, and the semiclassical local approximation is not validated. The significance of the paper depends crucially on whether the winding survives for generic hopping asymmetries and in a full moiré treatment.

major comments (5)
  1. [Fig. 4 and the text following it (six-direction quantization)] The six-direction quantization of D is inferred from panels in Fig. 4 in which only one hopping differs from the other two (e.g., t1 ≠ t2 = t3). At a generic moiré registry all three hoppings Δt1, Δt2, Δt3 are distinct, no mirror symmetry pins D onto a bond axis, and D is a continuous function of the hopping parameters except at gap closures. The statement that 'slight asymmetries barely distort this picture' is not quantified, and no calculation is shown for a generic (Δt1, Δt2, Δt3) triple. To support the central claim, the authors should compute D for several generic registries inside the moiré, show that the direction indeed stays on the six axes (or identify the conditions under which it does not), and specify how the discrete directions are selected by the pattern of hoppings.
  2. [Introduction and Fig. 6(d) (semiclassical local-unit-cell approximation)] The paper explicitly states that the main assumption is constructing D as a local vector field via locally commensurate unit cells. This assumption is not validated against any full moiré supercell calculation. The (54,55) supercell contains 17,714 atoms, which makes a full DFT calculation impractical, but a model-level check (for example, computing the Berry curvature of a full moiré tight-binding Hamiltonian at a few representative positions, or comparing with a continuum-model average) would provide a direct test of whether the local-unit-cell D(r) represents the actual response. Without such a check, the winding map in Fig. 6(d) remains an extrapolation from parametrized local cells.
  3. [Equation (3) and the numerical methods paragraph] No convergence tests are reported for the k-space integration in Eq. (3), and no error estimates are given for the D values. The text only says that Gauss-Legendre quadrature is used in the region where the Berry curvature is largest. Because the reported D magnitudes are small (around 1 Å) and the chemical-potential features are below 1 meV (temperature 10^-6 K), the results may be sensitive to the integration grid, the size of the integration region, and the smearing. The authors should provide a convergence study, for example a plot of D versus the number of quadrature points for one representative registry, and state the numerical uncertainty on the plotted vectors in Fig. 6(d).
  4. [Eq. (5) and the parameter-space mapping (82 registries)] The procedure by which Δtj and ΔE are obtained for the 82 moiré positions is not described. It is unclear whether DFT calculations were performed for every local registry or whether the results were interpolated from a small set of configurations (ABA, ABC, and the two asymmetric cells in Fig. 5). This matters because the six-direction quantization and the winding pattern depend on the functional form of Δtj as a function of local displacement. The authors should state the protocol explicitly: which registries were computed ab initio, how the remaining ones were generated, and what error is introduced by the interpolation.
  5. [Conclusion and Eq. (4) (longitudinal-current claim)] The conclusion states that the winding 'opens the possibility of longitudinal (as well as transversal) second order currents based solely on the Berry dipole.' However, Eq. (4) relates D to transverse conductivities (σyxx and σxyy) in a homogeneous system. For a moiré with spatially varying D, one must specify how the local currents are averaged or scattered to produce a measurable longitudinal component. The authors acknowledge that 'a subsequent scattering that was not calculated here' is needed, but this leaves the headline application without a concrete derivation. The winding claim itself is independent of this, but the framing should be moderated or the current formula for a spatially modulated D should be derived.
minor comments (4)
  1. [Abstract and Introduction] The phrase 'homobilayer of hexagonal boron nitride' is used, but the system is graphene on an hBN bilayer, which is a trilayer; consider a clearer wording such as 'graphene on an AB-stacked hBN bilayer.'
  2. [Fig. 4 caption] The caption says 'the (0,0) point in all subplots of Fig. 4 to correspond with the location of shifted the K−point'; this should read 'the shifted K−point'.
  3. [Fig. 6(d)] Unlike panels (a)–(c), panel (d) has no visible color scale or legend for the vector magnitudes. Adding a reference arrow and a scale (or stating that all vectors are normalized) would improve readability.
  4. [References] Reference [14] (Sethna, Statistical Mechanics) is a general textbook; citing it for the notion of 'parameter space' of local registries is somewhat loose. A more specific reference to the moiré registry parameterization, or a short explanation in the text, would help the reader.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: D(r) is computed from ab initio ΔV and tight-binding parameters, not fitted or reduced to self-cited inputs.

full rationale

The derivation of D(r) is self-contained. ΔV(r) is obtained from DFT, Eq. (5) converts it into renormalized hoppings and on-site energies with t = -2.83 eV from the literature, Eq. (6) defines the local Hamiltonian, and Eq. (3) yields the Berry dipole without reference to any target D value. The winding map in Fig. 6(d) is a post-processing of these computed D vectors over 82 registries, and the six-direction quantization is argued from the trigonal symmetry of the hopping model (Fig. 4) with external Ref. [13], not from a self-cited uniqueness theorem. The only self-citation, Ref. [15] to a review on strained graphene by one of the present authors, is used only as an analogy for the semiclassical strain-engineering approximation; the same assumption is supported by external Ref. [7]. The claim about longitudinal currents is explicitly flagged as a possibility ('a subsequent scattering that was not calculated here'), not presented as a derived consequence. No fitted parameter is relabeled as a prediction, and no load-bearing step reduces by construction to an input. The unquantified 'slight asymmetries barely distort' statement and the sparse 82-point sampling affect the robustness of the six-direction claim, but they are correctness concerns, not circularity. Hence no circularity.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The calculation relies on standard tight-binding and DFT inputs. The only hand-chosen parameters are the 0.5% strain and the artificial 0.5 eV gap used for illustration in Fig. 4. No new physical entities are postulated.

free parameters (2)
  • Uniaxial strain epsilon = 0.5%
    Graphene is elongated along the x direction by 0.5% to break threefold symmetry; this is a chosen model parameter, not fitted to the target result.
  • Gap offset DeltaE in Fig. 4 = -0.5 eV
    Used to make Berry curvature maps more visible; not the physical 10 meV gap and not used in the main D(r) maps.
assumptions (5)
  • standard math Berry dipole formula from Sodemann and Fu and the semiclassical transport relation for nonlinear currents
    Equations (3) and (4) are taken from the literature and define the computed quantity and its physical meaning.
  • domain assumption Tight-binding pi-electron model with only first-nearest-neighbor hoppings
    Equations (5) and (6) neglect higher-order hoppings and other orbitals; this is a standard but restrictive model for graphene.
  • domain assumption Semiclassical local field approximation
    The paper states that constructing D as a local vector field is its main assumption; this is central to interpreting the maps in Fig. 6.
  • domain assumption PBE-GGA density functional theory with van der Waals corrections gives a reliable effective potential
    The DFT potential DeltaV is the input to the tight-binding model; functional approximations affect the quantitative values of on-site and hopping corrections.
  • domain assumption Moiré supercell is commensurate and twist-free
    The (54,55) cell with homogeneous strain is chosen to represent the graphene/hBN interface, as in prior work on this system.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Winding Berry dipole on uniaxially strained graphene/hBN/hBN moir\'e trilayers." pith.science (2026). https://pith.science/paper/IPGYEGU6

@misc{pith2026241210584,
  author       = {Pith},
  title        = {Pith review of: Winding Berry dipole on uniaxially strained graphene/hBN/hBN moir\'e trilayers},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IPGYEGU6}},
  note         = {Machine review of arXiv:2412.10584}
}
abstract

Nonlinear Hall-like currents can be generated by a time-periodic alternating bias on two-dimensional (2D) materials lacking inversion symmetry. To hint that the moir\'e between graphene and its supporting substrate contributes to the homogeneity of nonlinear currents, the change in the local potential $\Delta V(r)$ around horizontally strained graphene due to a homobilayer of hexagonal boron nitride (hBN) was obtained from ab initio calculations, and corrections to on-site energies and hopping matrix elements on graphene's tight-binding electronic dispersion of $\pi-$electrons were calculated. Relying on a semiclassical approximation, Berry dipoles $D$ are seen to change orientation and wind throughout the moir\'e lattice.

Figures

Figures reproduced from arXiv: 2412.10584 by the authors.

Figure 1
Figure 1. FIG. 1. (a) Radial part (times [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Berry dipole [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. FIG. 4. Berry curvature for the valence band around the [PITH_FULL_IMAGE:figures/full_fig_p003_4.png] view at source ↗
Figures from the paper (2 more)
Figure 5
Figure 5. Figure 5: FIG. 5. (a) Moir [PITH_FULL_IMAGE:figures/full_fig_p004_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. (a) Relative displacement between the graphene unit cells and the closest hBN ones, for over 82 local registries within the moir [PITH_FULL_IMAGE:figures/full_fig_p005_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 23 canonical work pages

  1. [1]

    Sodemann and L

    I. Sodemann and L. Fu, Quantum nonlinear Hall effect induced by Berry curvature dipole in time-reversal invariant materials, Phys. Rev. Lett. 115, 216806 (2015)

    2015

  2. [2]

    Battilomo, N

    R. Battilomo, N. Scopigno, and C. Ortix, Berry curvature dipole in strained graphene: a Fermi surface warping effect, Phys. Rev. Lett. 123, 196403 (2019)

    work page 2019

  3. [3]

    Xiao, M.-C

    D. Xiao, M.-C. Chang, and Q. Niu, Berry phase effects on electronic properties, Rev. Mod. Phys. 82, 1959 (2010)

    2010

  4. [4]

    Zhang, J

    C.-P. Zhang, J. Xiao, B. T. Zhou, J.-X. Hu, Y .-M. Xie, B. Yan, and K. T. Law, Giant nonlinear Hall effect in strained twisted bilayer graphene, Phys. Rev. B 106, L041111 (2022)

    work page 2022

  5. [5]

    P. He, G. K. W. Koon, H. Isobe, J. Y . Tan, J. Hu, A. H. C. Neto, L. Fu, and H. Yang, Graphene moir ´e superlattices with giant quantum nonlinearity of chiral Bloch electrons, Nat. Nanotechnol. 17, 378 (2022)

    work page 2022

  6. [6]

    Z. Z. Du, H.-Z. Lu, and X. C. Xie, Nonlinear Hall effects, Nat. Rev. Phys. 3, 744 (2021)

    2021

  7. [7]

    J. Jung, E. Laksono, A. M. DaSilva, A. H. MacDonald, M. Mucha-Kruczy ´nski, and S. Adam, Moir ´e band model and band gaps of graphene on hexagonal boron nitride, Phys. Rev. 6 B 96, 085442 (2017)

    work page 2017

  8. [8]

    Cazeaux, D

    P. Cazeaux, D. Clark, R. Engelke, P. Kim, and M. Luskin, Relaxation and domain wall structure of bilayer moir´e systems, J. Elast. 154, 443 (2023)

    work page 2023

Show all 32 references
  1. [9]

    Bennett, G

    D. Bennett, G. Chaudhary, R. J. Slager, E. Bousquet, and P. Ghosez, Polar meron-antimeron networks in strained and twisted bilayers, Nat. Commun. 14, 1629 (2023)

    2023

  2. [10]

    S. Li, Z. Sun, N. J. McLaughlin, A. Sharmin, N. Agarwal, M. Huang, S. H. Sung, H. Lu, S. Yan, H. Lei, R. Hovden, H. Wang, H. Chen, L. Zhao, and C. R. Du, Observation of stacking engineered magnetic phase transitions within moir ´e supercells of twisted van der Waals magnets, N...

    2024

  3. [11]

    C. Hu, M. H. Naik, Y .-H. Chan, J. Ruan, and S. G. Louie, Light-induced shift current vortex crystals in moir ´e heterobilayers, Proc. Natl. Acad. Sci. U.S.A.120, e2314775120 (2023)

    2023

  4. [12]

    A. O. Fumega and J. L. Lado, Ferroelectric valley valves with graphene/mote2 van der waals heterostructures, Nanoscale 15, 2181 (2023)

    2023

  5. [13]

    Engelke, H

    R. Engelke, H. Yoo, S. Carr, K. Xu, P. Cazeaux, R. Allen, A. M. Valdivia, M. Luskin, E. Kaxiras, M. Kim, J. H. Han, and P. Kim, Topological nature of dislocation networks in two-dimensional moir´e materials, Phys. Rev. B 107, 125413 (2023)

    2023

  6. [14]

    J. P. Sethna, Statistical Mechanics: Entropy, Order Parameters, and Complexity, 1st ed. (Oxford U. Press, Oxford, U.K., 2006)

    2006

  7. [15]

    G. G. Naumis, S. Barraza-Lopez, M. Oliva-Leyva, and H. Terrones, Electronic and optical properties of strained graphene and other strained 2D materials: a review, Rep. Prog. Phys. 80, 096501 (2017)

    2017

  8. [16]

    Kresse and J

    G. Kresse and J. Furthm ¨uller, Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set, Phys. Rev. B 54, 11169 (1996)

    1996

  9. [17]

    Kresse and J

    G. Kresse and J. Furthm¨uller, Efficiency ofab initio total energy calculations for metals and semiconductors using a plane-wave basis set, Comput. Mater. Sci. 6, 15 (1996)

    1996

  10. [18]

    Kresse and D

    G. Kresse and D. Joubert, From ultrasoft pseudopotentials to the projector augmented-wave method, Phys. Rev. B 59, 1758 (1999)

    1999

  11. [19]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

    1996

  12. [20]

    J. c. v. Klime ˇs, D. R. Bowler, and A. Michaelides, Van der Waals density functionals applied to solids, Phys. Rev. B 83, 195131 (2011)

    2011

  13. [21]

    Klime ˇs, D

    J. Klime ˇs, D. R. Bowler, and A. Michaelides, Chemical accuracy for the van der Waals density functional, J. Phys.: Condens. Matter 22, 022201 (2009)

    2009

  14. [22]

    Sinha, P

    S. Sinha, P. C. Adak, A. Chakraborty, K. Das, K. Debnath, L. D. V . Sangani, K. Watanabe, T. Taniguchi, U. V . Waghmare, A. Agarwal, and M. M. Deshmukh, Berry curvature dipole senses topological transition in a moir´e superlattice, Nat. Phys. 18, 765 (2022)

    2022

  15. [23]

    A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, The electronic properties of graphene, Rev. Mod. Phys. 81, 109 (2009)

    2009

  16. [24]

    R. W. Jansen and O. F. Sankey, Ab initio linear combination of pseudo-atomic-orbital scheme for the electronic properties of semiconductors: results for ten materials, Phys. Rev. B 36, 6520 (1987)

    1987

  17. [25]

    O. F. Sankey and D. J. Niklewski, Ab initio multicenter tight-binding model for molecular-dynamics simulations and other applications in covalent systems, Phys. Rev. B 40, 3979 (1989)

    1989

  18. [26]

    J. D. Bernal, The structure of graphite, Proc. R. Soc. Lond. A 106, 749 (1924)

    1924

  19. [27]

    Ma, S.-Y

    Q. Ma, S.-Y . Xu, H. Shen, D. MacNeill, V . Fatemi, T.-R. Chang, A. M. Mier Valdivia, S. Wu, Z. Du, C.-H. Hsu, S. Fang, Q. D. Gibson, K. Watanabe, T. Taniguchi, R. J. Cava, E. Kaxiras, H.-Z. Lu, H. Lin, L. Fu, N. Gedik, and P. Jarillo-Herrero, Observation of the nonlinear Hall...

    2019

  20. [28]

    S.-Y . Xu, Q. Ma, H. Shen, V . Fatemi, S. Wu, T.-R. Chang, G. Chang, A. M. M. Valdivia, C.-K. Chan, Q. D. Gibson, J. Zhou, Z. Liu, K. Watanabe, T. Taniguchi, H. Lin, R. J. Cava, L. Fu, N. Gedik, and P. Jarillo-Herrero, Electrically switchable Berry curvature dipole in the mono...

    2018

  21. [29]

    K. Kang, T. Li, E. Sohn, J. Shan, and K. F. Mak, Nonlinear anomalous Hall effect in few-layer WTe 2, Nat. Mater. 18, 324 (2019)

    2019

  22. [30]

    D. K. Efimkin and A. H. MacDonald, Helical network model for twisted bilayer graphene, Phys. Rev. B 98, 035404 (2018)

    2018

  23. [31]

    Zhang, A

    F. Zhang, A. H. MacDonald, and E. J. Mele, Valley chern numbers and boundary modes in gapped bilayer graphene, Proc. Natl. Acad. Sci. U.S.A. 110, 10546 (2013)

    2013

  24. [32]

    J. D. Verbakel, Q. Yao, K. Sotthewes, and H. J. W. Zandvliet, Valley-protected one-dimensional states in small-angle twisted bilayer graphene, Phys. Rev. B 103, 165134 (2021)

    2021

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.