
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39266571
52314
10.1038/s41467-024-52314-5
Article
Intrinsic dipole Hall effect in twisted MoTe2: magnetoelectricity and contact-free signatures of topological transitions
Fan Feng-Ren 12
http://orcid.org/0000-0002-4843-070X
Xiao Cong 23
http://orcid.org/0000-0003-2883-4528
Yao Wang wangyao@hku.hk

12
1 https://ror.org/02zhqgq86 grid.194645.b 0000 0001 2174 2757 New Cornerstone Science Laboratory, Department of Physics, University of Hong Kong, Hong Kong, China
2 grid.194645.b 0000000121742757 HKU-UCAS Joint Institute of Theoretical and Computational Physics at Hong Kong, Hong Kong, China
3 grid.437123.0 0000 0004 1794 8068 Institute of Applied Physics and Materials Engineering, University of Macau, Taipa, Macau China
12 9 2024
12 9 2024
2024
15 799713 2 2024
30 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
We discover an intrinsic dipole Hall effect in a variety of magnetic insulating states at integer fillings of twisted MoTe2 moiré superlattice, including topologically trivial and nontrivial ferro-, antiferro-, and ferri-magnetic configurations. The dipole Hall current, in linear response to in-plane electric field, generates an in-plane orbital magnetization M∥ along the field, through which an AC field can drive magnetization oscillation up to THz range. Upon the continuous topological phase transitions from trivial to quantum anomalous Hall states in both ferromagnetic and antiferromagnetic configurations, the dipole Hall current and M∥ have an abrupt sign change, enabling contact-free detection of the transitions through the magnetic stray field. In configurations where the linear response is forbidden by symmetry, the dipole Hall current and M∥ appear as a crossed nonlinear response to both in-plane and out-of-plane electric fields. These magnetoelectric phenomena showcase fascinating functionalities of insulators from the interplay between magnetism, topology, and electrical polarization.

The interplay between magnetism, topology, and electrical polarization is a fascinating research direction in modern condensed matter physics. Here, the authors discover an intrinsic dipole Hall effect in various magnetic insulating states of twisted MoTe2—enabling contact-free detection of topological phase transitions.

Subject terms

Topological matter
Spintronics
Research Grant Council of Hong Kong SAR China (AoE/P-701/20, HKU SRFS2122-7S05, A-HKU705/21) and New Cornerstone Science Foundation.Research Grant Council of Hong Kong SAR China (AoE/P-701/20, HKU SRFS2122-7S05, A-HKU705/21) and New Cornerstone Science Foundation.Research Grant Council of Hong Kong SAR China (AoE/P-701/20, HKU SRFS2122-7S05, A-HKU705/21) and New Cornerstone Science Foundation. UM Start-up Grant (SRG2023-00033-IAPME).issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

In magnetic insulators, in-depth exploration of magnon transport has unveiled promising opportunities for low-power-consumption information technologies1,2. As the electrical counterpart of magnon, a charge-neutral elementary excitation carrying electric dipole can also transport energy, momentum, as well as electrical polarization in an insulator3. Such a concept has caught limited attention until very recently in the context of ferroelectrics where the dipole excitation corresponds to certain phonon3–9. Notably, flow of electric dipole in directions perpendicular to its moment is detectable via the accompanied magnetic stray field9.

Rhombohedral (R) homobilayers of transition metal dichalcogenides (TMDs) are a versatile platform that hosts ferroelectricity10,11, intrinsic magnetism12, and nontrivial topology13–16. With inversion symmetry broken, the electron affinity difference between atoms leads to an interlayer electrical polarization having opposite signs at the MX and XM stacking configurations that are related by an interlayer sliding (Fig. 1a)10,11,17–20. Incommensurate and marginally twisted homobilayers, switching of the polarization accompanied by sliding can be achieved by a sizable out-of-plane electric field10,11,20–23, forming the basis of ferroelectric functionalities.Fig. 1 Schematic of the dipole Hall effect in twisted MoTe2.

a Schematic of twisted R-stack TMDs bilayer, featuring opposite out-of-plane electrical polarizations at MX and XM stacking regions. b FM and AFM insulating states at filling ν = − 1 and −2. Dashed circles denote the moiré orbitals of doped carriers, polarized in opposite layers at MX and XM sites. The colored area signifies both the filling fraction and spin. c, d Dipole Hall response to in-plane electric field E∥ for (c) QAH and (d) trivial magnetic insulating state. The red dashed arrows denote the interlayer tunneling current due to annihilation of accumulated dipole on edges, and the red solid arrows denote the layer counter flows corresponding to the bulk dipole Hall current. Together, they form a current loop that generates the in-plane magnetization in bulk. Upon the topological phase transition controlled by perpendicular field E⊥, the dipole Hall current and associated orbital magnetization M∥ (blue arrows) have an abrupt sign change. e Equilibrium (black) and total (red) magnetization in an AC field E∥(ω).

Twisting the homobilayers by a modest angle results in a moiré pattern where stacking registries alternate between MX and XM with few nm periodicities (Fig. 1a). The stacking-dependent electrical polarization then becomes an antiferroelectric background pinned by the moiré, in which doped carriers experience a hexagonal superlattice with two degenerate moiré orbitals at MX and XM regions polarized in opposite layers respectively24,25. Berry phase from such layer texture of carrier manifests as an emergent magnetic field of quantized flux per moiré cell25,26, which underlies nontrivial topology of low-energy minibands24. With the intrinsic ferromagnetism from Coulomb exchange between moiré orbitals12, this system has become an exciting platform for exploring quantum anomalous Hall (QAH) effects, where both integer and fractional QAH effects were observed in twisted MoTe2 (tMoTe2)13–16. Through the antiferroelectric nature of the carrier wave function, both the magnetism and topology can be manipulated by a modest perpendicular electric field (E⊥)12,25. At filling factor ν = − 1, i.e., one hole per moiré cell, experiment has reported a continuous topological phase transition from ferromagnetic (FM) QAH to a trivial FM, with the increase of E⊥14. Hartree-Fock calculations further suggest the existence of an antiferromagnetic (AFM) state at ν = − 2, featuring a topological transition from trivial to an AFM QAH upon the increase of E⊥27–29.

Here we discover an intrinsic dipole Hall effect generally present in a variety of magnetic insulating states at integer filling factors in twisted R stack TMDs homobilayers. In the insulating bulk, a pure flow of interlayer dipole excitation of the doped carrier is generated due to the dipole Berry curvature in superlattice minibands. Such quantum geometric origin allows linear dipole current response to an in-plane electric field E∥, rather than the usual field gradient. On top of the equilibrium magnetization out-of-plane, the dipole Hall current corresponds to an in-plane orbital magnetization M∥ along E∥. Through this magnetoelectric response, an AC electrical field can thus drive magnetization oscillations up to the terahertz range. Remarkably, upon the continuous topological quantum phase transitions tuned by E⊥ in both the ν = −1 FM and ν = −2 AFM states, the dipole Hall conductivity and the associated M∥ have an abrupt change, enabling contact-free detection of the transitions through the magnetic stray field. In the ν = −1 ferro- and a ν = −3 ferri-magnetic configurations where this linear response is forbidden at E⊥ = 0 by the C2yT symmetry, the dipole Hall current and M∥ appear as a nonlinear response to both E⊥ and E∥.

Results

Quantum geometric origin

For electrons in a coupled bilayer, the dipole current operator reads j^a=e2{v^a,σ^z}. Here σ^z is the Pauli matrix in the layer index subspace, and represents the interlayer charge dipole p^=ed0σ^zz^, where d0 is the interlayer distance. By the semiclassical theory, the dipole current contributed by an electron is given by (see Supplementary Note 1):1 Jan(k)=jan(k)+ϒabn(k)E∥,b,

where n and k are the band index and the wave vector, respectively, and summation over repeated Cartesian indices a, b is implied. The first term on the right side is the expectation value of j^a. The second term is the anomalous dipole current induced by field, where2 ϒabn(k)=2eℏIm∑n′≠njann′(k)vbn′n(k)[εn(k)−εn′(k)]2

can be termed, in the same spirit of the spin Berry curvature30, as dipole Berry curvature. Here εn(k) is the band energy, and the numerator involves the interband matrix elements of interlayer-dipole current and velocity operators.

Summing over k in the filled bands of the insulator yields the total dipole current density. We find an intrinsic dipole current in linear response to E∥: jD,a = σabE∥,b, where σab=∑n∫[dk]f0(k)ϒabn(k), f0 the equilibrium Fermi distribution, and [dk] is shorthand for dk/(2π)2. For the off-diagonal component of dipole conductivity σab, the threefold rotation symmetry C3z forbids its symmetric part, thus σxy = −σyx, and the transverse transport is described by a Hall conductivity σD=σxy−σyx/2. The dipole Hall current reads3 jD=σDE∥×z^,

where the dipole Hall conductivity is given by4 σD=∑n∫[dk]f0(k)ϒDn(k),

with ϒDn(k)=ϵabϒabn(k)/2 being the antisymmetric part of the dipole Berry curvature. Here ϵ is the Levi-–Civita symbol. As a time-reversal odd pseudoscalar (zz component of a pseudotensor), σD is allowed by rotation and primed improper rotations (which are combinations of time reversal with inversion, mirror, or roto-reflection).

Longitudinal orbital magnetoelectric response

The interlayer-dipole Hall current is also manifested as an in-plane orbital magnetization parallel to E∥, as shown by the schematics in Fig. 1c and d. This is physically intuitive as a traveling charge dipole can induce an orbital magnetic moment31. To see this connection formally, one notices that the in-plane orbital magnetic moment operator m^=14(p^×v^−v^×p^) can be recast into m^=d02z^×j^, where j^ is the aforementioned dipole current operator (details in Supplementary Note 2). This form of the in-plane orbital magnetic moment operator can also be obtained by rigorous quantum mechanical treatment of magnetic field effect in the continuum model of coupled twisted bilayers32,33, as shown in Supplementary Note 2.

The equilibrium in-plane magnetization is prohibited by the C3z symmetry. Thus, the in-plane orbital magnetization appears from the first order of electric field M∥,a = χabE∥,b. For the examples to be discussed, C3z symmetry renders χxx = χyy, whereas MxT symmetry forbids χxy, χyx, leaving the magnetoelectric response a longitudinal form,5 M∥=d02σDE∥.

Namely, this longitudinal magnetoelectric response is equivalent to the intrinsic dipole Hall effect (details in Supplementary Note 2). Therefore, measuring the in-plane magnetization, by magneto-optical means or via the magnetic stray field, allows a contact-free detection of the dipole Hall effect. Moreover, an AC E∥ field will drive an oscillating M∥ that adds on top of the equilibrium magnetization (Fig. 1e), such that the net magnetization can process with a frequency upper bound by the charge gap, reaching THz range for the examples below.

Before proceeding to specific behaviors of the proposed effects, some comments are in order.

First, the orbital magnetization in linear response to E field in 2D insulators is a boundary-independent bulk thermodynamic quantity. Based on the in-plane orbital moment operator obtained by the quantum mechanical treatment and perturbation calculations on its response to in-plane electric field, we determine the in-plane intrinsic orbital magnetoelectric response unambiguously (see Supplementary Note 2). Consistently, one can also obtain this response coefficient by using the Maxwell relation, which states that it also quantifies the linear response of in-plane electric polarization to in-plane magnetic field. We have performed direct perturbation calculation of the latter response and get the expected same result (see Supplementary Note 2), further corroborating our theory for the in-plane orbital magnetoelectric response.

Second, the measurable longitudinal orbital magnetoelectric susceptibility is found to be quantitatively equivalent to the bulk dipole Hall conductivity in terms of the conventional dipole current definition employed here, although there can be alternative definition of dipole current operator in the presence of interlayer coupling (similar to the case of defining spin current in the presence of spin-orbit coupling34). On the one hand, it is reminded that our formulated magnetoelectric response is unambiguous, independent of the definition of dipole Hall current. On the other hand, the quantitative equivalence means that the in-plane magnetoelectric response renders a way to measure the bulk dipole Hall conductivity of the conventional dipole current.

Third, the interlayer tunneling plays a critical role in making relevant an in-plane orbital magnetoelectricity in the context of a double-layer 2D system. It is the interplay of the interlayer tunneling and intralayer moiré potentials that underlies this magnetoelectric response while setting either to zero will lead to a null response (Supplementary Fig. 3).

Dipole Hall effect in the FM insulator at ν = −1

We first demonstrate the dipole Hall conductivity at ν = −1 for 3.9° tMoTe2. Figure 2a shows the variation of Hall conductivities with interlayer bias, from a Hartree-Fock mean-field calculation. The black-circle (orange-triangle) curve represents the charge (dipole) Hall conductivity. At small interlayer bias E⊥, the FM is in a QAH state with Chern number 1, while the dipole Hall conductivity vanishes at E⊥ = 0 due to the C2yT symmetry. A finite E⊥ breaks the C2yT and switches on the dipole Hall effect. The FM transits from QAH to a topologically trivial state at a critical E⊥, upon which the dipole Hall conductivity undergoes an abrupt change.Fig. 2 Dipole Hall effect at ν = − 1 in tMoTe2.

a Charge (black circles) and dipole (orange triangles) Hall conductivity as a function of E⊥. Twist angle θ = 3. 9∘. b–d Berry curvatures (upper panels) and quasiparticle band dispersions (lower panels) along k path γ1 to γ2 (c.f. inset for moiré Brillouin zone). The band dispersion is color-coded with the layer polarization: PL = 1 (−1) for Bloch states fully polarized on layer T(B). e Dipole Hall conductivity as a function of twist angle and E⊥, where the abrupt sign changes mark the topological phase boundaries.

The simultaneous abrupt change in the charge and dipole Hall conductivities is not a coincidence. The band geometric properties arising from layer pseudospin textures underly both quantities, which become noncontinuous at the band inversion induced by the interlayer bias. To see this, we introduce layer-projected Berry curvature ΩL (L = T, B) which satisfies ΩT + ΩB = Ω and ΩT −ΩB = ΥD (see Supplementary Note 1), where Ω is the usual k-space Berry curvature. Figure 2b and c show the Berry curvatures (upper panels) and band dispersion (lower panels) before and after the topological phase transition, respectively. Notably, at these positive interlayer biases, the filled band predominantly occupies layer B, and ΩB dominates over ΩT in magnitude. When E⊥ crosses the transition point, band inversion occurs in the vicinity of  −κ point, where the Berry curvatures all get reversed. As a result, an abrupt change of σD accompanies the change of Chern number. The reversal of the associated in-plane orbital magnetization (Fig. 1c and d) enables contact-free detection of this topological transition.

The phase diagram for the dipole Hall conductivity as a function of twist angle and bias is presented as Fig. 2e. One observes that as the twist angle decreases, the critical bias for dipole-Hall jump decreases. This is because as the twist angle decreases, the energy band also narrows, and the critical bias decreases accordingly.

Dipole Hall effect in the AFM insulator at ν = −2

Compared to the case of ν = −1, the ν = −2 AFMz state exhibits a different symmetry, where C2yT is replaced by C2y (Fig. 1b). Therefore, a pronounced dipole Hall effect is present at E⊥ = 0 (Fig. 3c). As the interlayer bias is increased, the dipole Hall conductivity undergoes a sudden jump (Fig. 3c) upon the topological transition to the AFMz QAH state, reminiscent of the finding in the ν =−1 FM configuration. We separately examine the contributions from the spin-down and spin-up channels to charge and dipole Hall conductivities, as shown in Fig. 3a and b, respectively. As the bias increases in the positive (negative) direction, the transition from trivial to QAH state occurs in the spin-down (spin-up) channel, whose contribution to the dipole Hall conductivity has a sudden sign change. The underlying picture is similar to the case of ν = −1, where the band inversion in a spin channel leads to an abrupt change in both σD and Chern number (Supplementary Fig. 2).Fig. 3 Dipole Hall effect in the AFMz insulating state at ν = − 2 in tMoTe2.

a–c Dipole Hall (orange triangles) and charge Hall conductivities (black circles) as functions of interlayer bias for (a) spin-down, (b) spin-up, and (c) total contributions, respectively. The twist angle θ = 3. 9∘. d Dipole Hall conductivity as a function of twist angle (θ) and interlayer bias.

The dependence of dipole Hall conductivity on twist angle and interlayer bias is shown in Fig. 3d. The topological phase diagram is complementary to that of the ν = −1 case, where critical E⊥ decreases with twisting angle27. With contributions from both spin up and down carriers, σD has larger magnitude here compared to the ν = −1 case. Near zero interlayer bias, we find a modest E∥ = 107 V/m can generate a sizable in-plane orbital magnetization M∥ ∼ 0.01 of μB/nm2.

Crossed nonlinear dipole Hall effect

The study of nonlinear Hall effect is another recent focus of condensed matter physics35–37. With σD symmetry forbidden in the ν =  −1 state at E⊥ = 0, its E⊥ dependence in Fig. 2a implies a type of nonlinearity—intrinsic nonlinear dipole Hall effect (Fig. 4a). We can define ΔDn(k)=∂ϒDn(k)/∂E⊥ as the dipole Berry curvature polarizability with respect to E⊥, which is also a band geometric quantity:6 ΔDn=eℏ2Im∑m≠n2jnm×vmnpn−pmεn−εm3−∑ℓ≠njℓm×vmn+vℓm×jmnpnℓεn−εℓεn−εm2+∑ℓ≠mjℓn×vnm+vℓn×jnmpmℓεm−εℓεn−εm2

and its flux through filled bands,7 α(E⊥)=∑n∫[dk]f0(k)ΔDn(k).

Substituting them into Eq. (3), we get a crossed nonlinear dipole current response,8 jD=α(0)E∥×E⊥.

The concomitant nonlinear orbital magnetoelectric response is given by9 M∥=d02α(0)E⊥E∥.

As shown in Fig. 4a, when the two electric fields are AC, M∥ has sum-frequency and difference-frequency components.Fig. 4 Crossed nonlinear dipole Hall effect at ν = − 1 and  −3 in tMoTe2.

a Schematic of the nonlinear orbital magnetoelectric response to two AC electric fields in-plane (E∥) and out-of-plane (E⊥), in which d0 is the interlayer spacing, α the flux of dipole Berry curvature polarizability. b The flux of dipole Berry curvature polarizability as a function of E⊥ at ν = − 1. c Distinct distributions of dipole Berry curvature polarizability in the moiré Brillouin zone at three different E⊥ in the QAH phase. d Quasiparticle band dispersion of ν = − 3 ferrimagnetic state, for which the twist angle θ = 3.5∘. The blue-solid (red-dashed) lines represent spin-up (spin-down) bands. The black-dashed line denotes the Fermi level (EF). e Layer-resolved carrier distribution of this ferrimagnetic state in a moiré supercell. Two spin-up holes occupy the B and C moiré orbitals respectively, and one spin-down hole occupies the layer-hybridized A orbital. f Dipole Hall conductivity as a function of interlayer bias. The black-dashed line is generated using σD = α(0)E⊥.

Surprisingly, Fig. 2a shows the simple scaling of the crossed nonlinear response given in Eqs. (8) and (9) are protected over the entire QAH phase, until a sudden jump occurs as the signature of the topological transition. Over such broad range of E⊥, the variation of the minibands (from Fig. 2d to b) does lead to significant changes in the distribution of ΔDn(k) in the moiré Brillouin zone, as shown in Fig. 4c. Nevertheless, their flux over the filled Chern band remains a constant in the QAH phase, in stark contrast to that in the topological trivial phases (Fig. 4b). This protected scaling throughout the QAH phase is observed at other twisting angles, and also for the AFMz configuration, while the actual flux can vary modestly with angle (details in Supplementary Fig. 4).

The nonlinear dipole Hall response is not limited to QAH state. We showcase another example at filling factor ν = −3, which hosts a topologically trivial ferrimagnetic insulating state, with two spin-up holes in the layer-polarized B and C orbitals and one spin-down hole in the layer-hybridized A orbital (Fig. 4d and e). Details of the Hartree-Fock calculation are given in Supplementary Note 4. Possessing C2yT symmetry as well (Fig. 4e), the linear dipole Hall effect is forbidden at zero interlayer bias, whereas the crossed nonlinear dipole Hall response is anticipated. This is confirmed by the calculated dipole Hall conductivity as a function of E⊥ given in Fig. 4f, which shows that the response is well captured by Eq. (8) for E⊥ range up to 40 mV/nm.

Discussion

We can also provide an alternative picture for the magnetoelectricity here. The bulk dipole Hall current tends to accumulate interlayer dipole on the boundary, which will get annihilated through interlayer tunneling. This tunneling current is inversely proportional to the dipole relaxation time and proportional to the dipole density. In steady state, the interlayer tunneling current on the edges will balance with the bulk dipole Hall current that corresponds to layer counter flows. Together they form a current loop that generates the in-plane magnetization in bulk.

Compared to the known magnetoelectric phenomena38,39, the magnetoelectricity from the intrinsic dipole Hall effect here has several features. The in-plane magnetoelectricity in 2D insulators here is solely of orbital origin, enabled by the spin-conserved interlayer tunneling across the twisted interface40,41. As an intrinsic effect in an insulator, the magnetoelectric response here is protected by a charge gap. This allows the phenomena to be explored in the AC regime, where the generated magnetization can adiabatically follow an AC electric field, generating magnetization oscillations at high frequencies upper bound by the charge gap. This could potentially provide a mechanism to induce magnetic oscillations of ferromagnets at faster timescales. The form of crossed nonlinear magnetoelectric susceptibility further points to novel device functionalities, where the nonlinearity underlies logic operations, second harmonic generation, rectification, etc.

In comparison with the magnon-based functionalities also based on magnetic insulators, the intrinsic dipole Hall effect exhibits several complementary features. The magnon transport is typically driven by temperature gradient and magnetic field gradient42, while the dipole Hall current is a response directly to an in-plane electric field through band geometric effect. The response can also be dramatically tuned by the out-of-plane electric field, including the magnitude of the susceptibility and the abrupt sign change when tuning across the topological phase transition point. Exploring the transport of electric dipole and that of magnon in magnetic insulators represent distinct and complementary approaches to low-power-consumption information technologies, which could potentially be combined to exploit the advantages of both.

Methods

Hartree–Fock mean-field method

To determine the ground magnetic states of tMoTe2 at different fillings, we employed a self-consistent k-space Hartree-Fock mean-field approach. Starting from the continuum model for twisted R-stacking MoTe2, Coulomb interactions were incorporated at the mean-field level27. Self-consistent solution of the mean-field Hamiltonian yielded the ground state for each filling factor. A detailed exposition of the methodology is provided in Supplementary Note 4.

Supplementary information

Supplementary Information “Intrinsic dipole Hall effect in twisted MoTe2: magnetoelectricity and contact-free signatures of topological transitions”

Peer Review File

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52314-5.

Acknowledgements

This work is supported by the Research Grant Council of Hong Kong SAR China (AoE/P-701/20, HKU SRFS2122-7S05, A-HKU705/21), the National Key R&D Program of China (2020YFA0309600), and New Cornerstone Science Foundation. C.X. acknowledges support by the UM Start-up Grant (SRG2023-00033-IAPME).

Author contributions

W.Y. conceived and initiated the research. F.-R.F. performed the Hartree–Fock calculations and analysis, based on the formula derived by C.X. All authors discussed the results, contributed to the interpretations, and wrote the manuscript.

Peer review

Peer review information

Nature Communications thanks the anonymous, reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

The numerical data generated by the custom codes for this study that support the findings are available from the corresponding author on request.

Code availability

The custom codes prepared for this study that support the findings are available on Code Ocean (10.24433/CO.0067402.v1).

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Feng-Ren Fan, Cong Xiao.
==== Refs
References

1. Chumak AV Vasyuchka VI Serga AA Hillebrands B Magnon spintronics Nat. Phys. 2015 11 453 10.1038/nphys3347
Chumak, A. V., Vasyuchka, V. I., Serga, A. A. & Hillebrands, B. Magnon spintronics. Nat. Phys. 11, 453 (2015).10.1038/nphys3347
2. Yu, H., Xiao, J. & Pirro, P. Magnon spintronics. Perspectives on Magnon Spintronics. J. Magn. Magn. 450, 1 (2018). https://www.sciencedirect.com/science/article/pii/S0304885317338118.
3. Bauer G Polarization transport in ferroelectrics Phys. Rev. Appl. 2023 20 050501 10.1103/PhysRevApplied.20.050501
Bauer, G. et al. Polarization transport in ferroelectrics. Phys. Rev. Appl. 20, 050501 (2023).10.1103/PhysRevApplied.20.050501
4. Shen K Magnon-Ferron coupling mediated by dynamical Dzyaloshinskii-Moriya interaction in a two-dimensional multiferroic model Phys. Rev. B 2022 106 104411 10.1103/PhysRevB.106.104411
Shen, K. Magnon-Ferron coupling mediated by dynamical Dzyaloshinskii-Moriya interaction in a two-dimensional multiferroic model. Phys. Rev. B 106, 104411 (2022).10.1103/PhysRevB.106.104411
5. Chotorlishvili L Khomeriki R Sukhov A Ruffo S Berakdar J Dynamics of localized modes in a composite multiferroic chain Phys. Rev. Lett. 2013 111 117202 10.1103/PhysRevLett.111.117202 24074117
Chotorlishvili, L., Khomeriki, R., Sukhov, A., Ruffo, S. & Berakdar, J. Dynamics of localized modes in a composite multiferroic chain. Phys. Rev. Lett. 111, 117202 (2013).24074117 10.1103/PhysRevLett.111.117202
6. Bauer GEW Iguchi R Uchida K-i Theory of transport in ferroelectric capacitors Phys. Rev. Lett. 2021 126 187603 10.1103/PhysRevLett.126.187603 34018785
Bauer, G. E. W., Iguchi, R. & Uchida, K.-i Theory of transport in ferroelectric capacitors. Phys. Rev. Lett. 126, 187603 (2021).34018785 10.1103/PhysRevLett.126.187603
7. Chen G Lan J Min T Xiao J Narrow waveguide based on ferroelectric domain Wall Chin. Phys. Lett. 2021 38 087701 10.1088/0256-307X/38/8/087701
Chen, G., Lan, J., Min, T. & Xiao, J. Narrow waveguide based on ferroelectric domain Wall. Chin. Phys. Lett. 38, 087701 (2021).10.1088/0256-307X/38/8/087701
8. Li Q Subterahertz collective dynamics of polar vortices Nature 2021 592 376 10.1038/s41586-021-03342-4 33854251
Li, Q. et al. Subterahertz collective dynamics of polar vortices. Nature 592, 376 (2021).33854251 10.1038/s41586-021-03342-4
9. Tang P Iguchi R Uchida K-i Bauer GEW Thermoelectric polarization transport in ferroelectric ballistic point contacts Phys. Rev. Lett. 2022 128 047601 10.1103/PhysRevLett.128.047601 35148138
Tang, P., Iguchi, R., Uchida, K.-i & Bauer, G. E. W. Thermoelectric polarization transport in ferroelectric ballistic point contacts. Phys. Rev. Lett. 128, 047601 (2022).35148138 10.1103/PhysRevLett.128.047601
10. Weston A Interfacial ferroelectricity in marginally twisted 2D semiconductors Nat. Nanotechnol. 2022 17 390 10.1038/s41565-022-01072-w 35210566
Weston, A. et al. Interfacial ferroelectricity in marginally twisted 2D semiconductors. Nat. Nanotechnol. 17, 390 (2022).35210566 10.1038/s41565-022-01072-w
11. Wang X Interfacial ferroelectricity in rhombohedral-stacked bilayer transition metal dichalcogenides Nat. Nanotechnol. 2022 17 367 10.1038/s41565-021-01059-z 35039684
Wang, X. et al. Interfacial ferroelectricity in rhombohedral-stacked bilayer transition metal dichalcogenides. Nat. Nanotechnol. 17, 367 (2022).35039684 10.1038/s41565-021-01059-z
12. Anderson E Programming correlated magnetic states with gate-controlled moiré geometry Science 2023 381 325 10.1126/science.adg4268 37347950
Anderson, E. et al. Programming correlated magnetic states with gate-controlled moiré geometry. Science 381, 325 (2023).37347950 10.1126/science.adg4268
13. Cai J Signatures of fractional quantum anomalous hall states in twisted MoTe2 Nature 2023 622 63 10.1038/s41586-023-06289-w 37315640
Cai, J. et al. Signatures of fractional quantum anomalous hall states in twisted MoTe2. Nature 622, 63 (2023).37315640 10.1038/s41586-023-06289-w
14. Park H Observation of fractionally quantized anomalous Hall effect Nature 2023 622 74 10.1038/s41586-023-06536-0 37591304
Park, H. et al. Observation of fractionally quantized anomalous Hall effect. Nature 622, 74 (2023).37591304 10.1038/s41586-023-06536-0
15. Xu F Observation of integer and fractional quantum anomalous hall effects in twisted bilayer MoTe2 Phys. Rev. X 2023 13 031037
Xu, F. et al. Observation of integer and fractional quantum anomalous hall effects in twisted bilayer MoTe2. Phys. Rev. X 13, 031037 (2023).
16. Zeng Y Thermodynamic Evidence of Fractional Chern Insulator in Moiré MoTe2 Nature 2023 622 69 10.1038/s41586-023-06452-3 37494955
Zeng, Y. et al. Thermodynamic Evidence of Fractional Chern Insulator in Moiré MoTe2. Nature 622, 69 (2023).37494955 10.1038/s41586-023-06452-3
17. Li L Wu M Binary compound bilayer and multilayer with vertical polarizations: two-dimensional ferroelectrics, multiferroics, and nanogenerators ACS Nano 2017 11 6382 10.1021/acsnano.7b02756 28602074
Li, L. & Wu, M. Binary compound bilayer and multilayer with vertical polarizations: two-dimensional ferroelectrics, multiferroics, and nanogenerators. ACS Nano 11, 6382 (2017).28602074 10.1021/acsnano.7b02756
18. Zhao P Xiao C Yao W Universal superlattice potential for 2D materials from twisted interface inside h-BN substrate NPJ 2D Mater. Appl. 2021 5 1 10.1038/s41699-021-00221-4
Zhao, P., Xiao, C. & Yao, W. Universal superlattice potential for 2D materials from twisted interface inside h-BN substrate. NPJ 2D Mater. Appl. 5, 1 (2021).10.1038/s41699-021-00221-4
19. Liang J Optically probing the asymmetric interlayer coupling in rhombohedral-stacked MoS2 bilayer Phys. Rev. X 2022 12 041005
Liang, J. et al. Optically probing the asymmetric interlayer coupling in rhombohedral-stacked MoS2 bilayer. Phys. Rev. X 12, 041005 (2022).
20. Rogée L Ferroelectricity in untwisted heterobilayers of transition metal dichalcogenides Science 2022 376 973 10.1126/science.abm5734 35617404
Rogée, L. et al. Ferroelectricity in untwisted heterobilayers of transition metal dichalcogenides. Science 376, 973 (2022).35617404 10.1126/science.abm5734
21. Vizner Stern M Interfacial ferroelectricity by van der Waals sliding Science 2021 372 1462 10.1126/science.abe8177
Vizner Stern, M. et al. Interfacial ferroelectricity by van der Waals sliding. Science 372, 1462 (2021).10.1126/science.abe8177
22. Woods CR Charge-polarized interfacial superlattices in marginally twisted hexagonal boron nitride Nat. Commun. 2021 12 347 10.1038/s41467-020-20667-2 33436620
Woods, C. R. et al. Charge-polarized interfacial superlattices in marginally twisted hexagonal boron nitride. Nat. Commun. 12, 347 (2021).33436620 10.1038/s41467-020-20667-2
23. Yasuda K Wang X Watanabe K Taniguchi T Jarillo-Herrero P Stacking-engineered ferroelectricity in bilayer boron nitride Science 2021 372 1458 10.1126/science.abd3230
Yasuda, K., Wang, X., Watanabe, K., Taniguchi, T. & Jarillo-Herrero, P. Stacking-engineered ferroelectricity in bilayer boron nitride. Science 372, 1458 (2021).10.1126/science.abd3230
24. Wu F Lovorn T Tutuc E Martin I MacDonald AH Topological insulators in twisted transition metal dichalcogenide homobilayers Phys. Rev. Lett. 2019 122 086402 10.1103/PhysRevLett.122.086402 30932597
Wu, F., Lovorn, T., Tutuc, E., Martin, I. & MacDonald, A. H. Topological insulators in twisted transition metal dichalcogenide homobilayers. Phys. Rev. Lett. 122, 086402 (2019).30932597 10.1103/PhysRevLett.122.086402
25. Yu H Chen M Yao W Giant magnetic field from moiré induced Berry phase in homobilayer semiconductors Natl Sci. Rev. 2020 7 12 10.1093/nsr/nwz117 35296065
Yu, H., Chen, M. & Yao, W. Giant magnetic field from moiré induced Berry phase in homobilayer semiconductors. Natl Sci. Rev. 7, 12 (2020).35296065 10.1093/nsr/nwz117
26. Zhai D Yao W Theory of tunable flux lattices in the homobilayer moiré of twisted and uniformly strained transition metal dichalcogenides Phys. Rev. Mater. 2020 4 094002 10.1103/PhysRevMaterials.4.094002
Zhai, D. & Yao, W. Theory of tunable flux lattices in the homobilayer moiré of twisted and uniformly strained transition metal dichalcogenides. Phys. Rev. Mater. 4, 094002 (2020).10.1103/PhysRevMaterials.4.094002
27. Fan F-R Xiao C Yao W Orbital chern insulator at ν = −2 in twisted MoTe2 Phys. Rev. B 2024 109 L041403 10.1103/PhysRevB.109.L041403
Fan, F.-R., Xiao, C. & Yao, W. Orbital chern insulator at ν = −2 in twisted MoTe2. Phys. Rev. B 109, L041403 (2024).10.1103/PhysRevB.109.L041403
28. Jiang K Zhou S Dai X Wang Z Antiferromagnetic Chern insulators in noncentrosymmetric systems Phys. Rev. Lett. 2018 120 157205 10.1103/PhysRevLett.120.157205 29756861
Jiang, K., Zhou, S., Dai, X. & Wang, Z. Antiferromagnetic Chern insulators in noncentrosymmetric systems. Phys. Rev. Lett. 120, 157205 (2018).29756861 10.1103/PhysRevLett.120.157205
29. Liu X Gate-tunable antiferromagnetic Chern insulator in twisted bilayer transition metal dichalcogenides Phys. Rev. Lett. 2024 132 146401 10.1103/PhysRevLett.132.146401 38640385
Liu, X. et al. Gate-tunable antiferromagnetic Chern insulator in twisted bilayer transition metal dichalcogenides. Phys. Rev. Lett. 132, 146401 (2024).38640385 10.1103/PhysRevLett.132.146401
30. Sun Y Zhang Y Felser C Yan B Strong intrinsic spin Hall effect in the TaAs family of weyl semimetals Phys. Rev. Lett. 2016 117 146403 10.1103/PhysRevLett.117.146403 27740815
Sun, Y., Zhang, Y., Felser, C. & Yan, B. Strong intrinsic spin Hall effect in the TaAs family of weyl semimetals. Phys. Rev. Lett. 117, 146403 (2016).27740815 10.1103/PhysRevLett.117.146403
31. Xiao C Niu Q Conserved current of nonconserved quantities Phys. Rev. B 2021 104 L241411 10.1103/PhysRevB.104.L241411
Xiao, C. & Niu, Q. Conserved current of nonconserved quantities. Phys. Rev. B 104, L241411 (2021).10.1103/PhysRevB.104.L241411
32. Lee JY Theory of correlated insulating behaviour and spin-triplet superconductivity in twisted double bilayer graphene Nat. Commun. 2019 10 5333 10.1038/s41467-019-12981-1 31767862
Lee, J. Y. et al. Theory of correlated insulating behaviour and spin-triplet superconductivity in twisted double bilayer graphene. Nat. Commun. 10, 5333 (2019).31767862 10.1038/s41467-019-12981-1
33. Qin W MacDonald AH In-plane critical magnetic fields in magic-angle twisted trilayer graphene Phys. Rev. Lett. 2021 127 097001 10.1103/PhysRevLett.127.097001 34506202
Qin, W. & MacDonald, A. H. In-plane critical magnetic fields in magic-angle twisted trilayer graphene. Phys. Rev. Lett. 127, 097001 (2021).34506202 10.1103/PhysRevLett.127.097001
34. Culcer D Semiclassical spin transport in spin-orbit-coupled bands Phys. Rev. Lett. 2004 93 046602 10.1103/PhysRevLett.93.046602 15323781
Culcer, D. et al. Semiclassical spin transport in spin-orbit-coupled bands. Phys. Rev. Lett. 93, 046602 (2004).15323781 10.1103/PhysRevLett.93.046602
35. Ma Q Observation of the nonlinear hall effect under time-reversal-symmetric conditions Nature 2019 565 337 10.1038/s41586-018-0807-6 30559379
Ma, Q. et al. Observation of the nonlinear hall effect under time-reversal-symmetric conditions. Nature 565, 337 (2019).30559379 10.1038/s41586-018-0807-6
36. Kang K Li T Sohn E Shan J Mak KF Nonlinear anomalous hall effect in few-layer WTe2 Nat. Mater. 2019 18 324 10.1038/s41563-019-0294-7 30804510
Kang, K., Li, T., Sohn, E., Shan, J. & Mak, K. F. Nonlinear anomalous hall effect in few-layer WTe2. Nat. Mater. 18, 324 (2019).30804510 10.1038/s41563-019-0294-7
37. Du ZZ Lu H-Z Xie X Nonlinear hall effects Nat. Rev. Phys. 2021 3 744 10.1038/s42254-021-00359-6
Du, Z. Z., Lu, H.-Z. & Xie, X. Nonlinear hall effects. Nat. Rev. Phys. 3, 744 (2021).10.1038/s42254-021-00359-6
38. Dong S Liu J-M Cheong S-W Ren Z Multiferroic materials and magnetoelectric physics: symmetry, entanglement, excitation, and topology Adv. Phys. 2015 64 519 10.1080/00018732.2015.1114338
Dong, S., Liu, J.-M., Cheong, S.-W. & Ren, Z. Multiferroic materials and magnetoelectric physics: symmetry, entanglement, excitation, and topology. Adv. Phys. 64, 519 (2015).10.1080/00018732.2015.1114338
39. Spaldin NA Ramesh R Advances in magnetoelectric multiferroics Nat. Mater. 2019 18 203 10.1038/s41563-018-0275-2 30783227
Spaldin, N. A. & Ramesh, R. Advances in magnetoelectric multiferroics. Nat. Mater. 18, 203 (2019).30783227 10.1038/s41563-018-0275-2
40. Stauber T Low T Gómez-Santos G Chiral response of twisted bilayer graphene Phys. Rev. Lett. 2018 120 046801 10.1103/PhysRevLett.120.046801 29437442
Stauber, T., Low, T. & Gómez-Santos, G. Chiral response of twisted bilayer graphene. Phys. Rev. Lett. 120, 046801 (2018).29437442 10.1103/PhysRevLett.120.046801
41. Zhai D Chen C Xiao C Yao W Time-reversal even charge hall effect from twisted interface coupling Nat. Commun. 2023 14 1961 10.1038/s41467-023-37644-0 37029110
Zhai, D., Chen, C., Xiao, C. & Yao, W. Time-reversal even charge hall effect from twisted interface coupling. Nat. Commun. 14, 1961 (2023).37029110 10.1038/s41467-023-37644-0
42. Bauer GE Tang P Iguchi R Uchida K-i Magnonics vs. Ferronics J. Magn. Magn. Mater. 2022 541 168468 10.1016/j.jmmm.2021.168468
Bauer, G. E., Tang, P., Iguchi, R. & Uchida, K.-i Magnonics vs. Ferronics. J. Magn. Magn. Mater. 541, 168468 (2022).10.1016/j.jmmm.2021.168468
