Dynamical Chiral Nernst Effect in Twisted Van der Waals Few Layers

The Nernst effect is a fundamental thermoelectric conversion phenomenon that was deemed to be possible only in systems with magnetic field or magnetization. In this work, we propose a novel dynamical chiral Nernst effect that can appear in two-dimensional van der Waals materials with chiral structural symmetry in the absence of any magnetic degree of freedom. This unconventional effect is triggered by time variation of an out-of-plane electric field, and has an intrinsic quantum geometric origin linked to not only the intralayer center-of-mass motion but also the interlayer coherence of electronic states. We demonstrate the effect in twisted homobilayer and homotrilayer transition metal dichalcogenides, where the strong twisted interlayer coupling leads to sizable intrinsic Nernst conductivities well within the experimental capacity. This work suggests a new route for electric control of thermoelectric conversion.


I. INTRODUCTION
The Nernst effect, in which an electric Hall current is generated in response to an applied temperature gradient, is not only a central ingredient for energy-harvesting thermoelectric devices, but also a fundamental physical effect for probing quantum geometric properties of solids.While the ordinary Nernst effect is triggered by external magnetic field, the anomalous Nernst effect originates from the momentum-space Berry curvature [1], an intrinsic quantum geometric quantity related to the center-ofmass motion of Bloch electrons [2,3].Such an intrinsic Nernst response eliminates the necessity of applying a magnetic field, and has attracted extensive experimental studies in recent years [4][5][6][7].On the other hand, the requirement for magnetic ordering raises the issue of stability of thermoelectric materials against external magnetic perturbations, as well as excludes nonmagnetic materials from the playground of intrinsic Nernst response.
In this work we propose an alternative scheme for triggering intrinsic Nernst response by electrical means without invoking any magnetic degrees of freedom.Our proposal is based on the gate tunability of two-dimensional (2D) van der Waals (vdW) layered materials and novel quantum geometric properties of layered electrons, whose wave functions are distributed across different layers.Specifically, the time evolution of gate field generates dynamically a positional shift of layered electrons within the 2D plane, which results in an induced Berry curvature pointing to the out-of-plane direction and enabling the intrinsic Nernst response.We reveal that the underlying intrinsic quantum geometric quantity is linked to not only the intralayer center-of-mass motion but also the interlayer coherence of layered electrons.We show that the effect can appear in nonmagnetic 2D vdW materials with chiral structural symmetry to provide the needed chirality of Hall transport, thus is dubbed as the dynamical chiral Nernst effect.
We demonstrate this effect in vdW layered structures with twisted stacking [8][9][10][11][12][13][14].In this class of materials the interlayer chiral coupling rotates the layer pseudospin about in-plane axes that are of topologically nontrivial textures in the twisted landscapes [15][16][17][18].The quantum layer degree of freedom activated by interlayer coherence of electronic wave functions [19,20] enables coupling the out-of-plane electric field to in-plane charge transport.As examples, we show sizable intrinsic chiral Nernst conductivities in twisted bilayer and trilayer transition metal dichalcogenides (TMDs).

II. QUANTUM GEOMETRIC ORIGIN OF THE EFFECT
The intrinsic contribution to Nernst effect is most readily derived within the semiclassical theory [2,3,21,22], which has been extended to include corrections to the band quantities due to external fields [23][24][25].In particular, in vdW layered materials, the out-of-plane field E = E ẑ is coupled to the interlayer dipole moment p in the form of −E p, hence the temporal variation of the gate field Ė modifies the Bloch state [26,27], inducing a k -space Berry connection (ℏk is the 2D crystal momentum) [28] A For a band with index n, we have whose numerator involves the interband matrix elements of the velocity operator and the interlayer charge dipole arXiv:2401.13278v1[cond-mat.mes-hall]24 Jan 2024 operator p.The specific form of p depends on the layer number [29].In a bilayer system, for example, p = ed 0 σz with σz the Pauli matrix in the layer index subspace and d 0 the interlayer distance, and ε n is the unperturbed band energy.One can see that G n favors interlayer hybridized electronic states: According to Eq. ( 2), if the state |u n ⟩ is fully polarized in a specific layer around some k, then G n (k) is suppressed.In addition, as G n is gauge invariant, A Ė can be identified physically as an in-plane positional shift of the electron [30,31] induced by the time evolution of the out-of-plane field.This in-plane positional shift implies a dynamical generation of k -space Berry curvature in linear order of Ė: which points to the out-of-plane direction.The band index n has been suppressed for simplicity.Ω (k) will play a role in various Berry-curvature related phenomena.Particularly, as a k -space effective magnetic field, it corresponds to an orbital magnetization of topological nature that is given by the integration of gΩ in k space [2], with being the grand potential density contributed by a particular state and k B is the Boltzmann constant.
To account for the anomalous Nernst effect, it is well established that the circulating orbital magnetization current should be discounted from the local charge current density so that the transport current takes the form of [2,23] where [dk] is shorthand for n d 2 k/(2π) 2 .One observes that the second term arises from the topological orbital magnetization due to the k -space Berry curvature.On the other hand, the first term is the conventional expression for the charge current, with ṙ being the semiclassical velocity of an electron wave packet and f (ε) the equilibrium Fermi-Dirac distribution function.We do not consider the off-equilibrium distribution function ∼ ∇T in the first order of in-plane temperature gradient, because it does not contribute to the Nernst current in nonmagnetic systems -the present focus-as required by time reversal symmetry.Moreover, since the temperature gradient ∇T does not enter into ṙ, the first term in Eq. ( 5) does not give any current in the linear order of ∇T .As a result, the Nernst current driven by ∇T arises totally from the second term of Eq. ( 5), which reads where the response coefficient is intrinsic to the band structure and is a Fermi surface property.
One finds that ϑ is a time-reversal even pseudoscalar, thus it is invariant under rotation, but flips sign under space inversion, reflection and rotoreflection symmetries.As such, ϑ is allowed only if the system possesses a chiral crystal symmetry [28,32].This new type of Nernst response, dubbed as the dynamical chiral Nernst effect, is the first scenario of intrinsic Nernst response without magnetic degree of freedom.2D chiral lattices, such as Tellurene [33,34] and twisted layers of honeycomb lattice [28,32,35], are suitable platforms for this effect.In particular, in twisted structures, the chirality is locked to the twist direction: Since twisted configurations with opposite twist angles are mirror images of each other, whereas the mirror reflection flips the sign of ϑ, the direction of the Nernst current is reversed by reversing the twist direction.This character renders a unique knob to control the effect.
For a sine AC gate field E = E 0 sin ωt, the oscillating Nernst current is given by Here ω is required to be below the threshold for direct interband transition in order to validate the semiclassical treatment, and α H has the dimension of Nernst conductivity and quantifies the Nernst response with respect to the in-plane temperature gradient.

III. APPLICATION TO TWISTED HOMOBILAYER TMDS
We employ the four-band model of twisted bilayer MoTe 2 in Ref. [15].The Hamiltonian for spin-up carriers in the K valley is where the subscripts t and b refer to the top and bottom layer, δ = θẑ × r is the local displacement vector at position r between the two layers.T l with l = t, b is the layer-dependent kinetic energy (10) in which l = 1 for t and l = −1 for b in the exponential terms, and ∆ g = 1.1 eV is the energy gap between conduction and valence bands of monolayer MoTe 2 .∆ T is the interlayer tunneling potential w vv , G i is generated from the counterclockwise (i − 1)π/3 rotation of G 1 = (0, 1)4π/( √ 3a 0 ), a 0 = 3.47 Å is the lattice constant of monolayer MoTe 2 , and w cc = −2 meV, w vv = −8.5 meV, w cv = w vc = 15.3 meV are the tunneling parameters (subscripts c and v represent conduction and valence bands, respectively).∆ l is the layer-dependent electrostatic potential with where V c = −11.94meV and V v = −16 meV.Fig. 1(b) shows the valence band structure, and the Nernst conductivity α H measured in units of α 0 = ek B /h ≈ 3.34 nA/K [36][37][38], in the K valley of 1.4 • twisted bilayer MoTe 2 .The color on the band structures denotes σ z n = ⟨u n |σ z |u n ⟩, which characterizes the strength of interlayer hybridization [28,32].Specifically, σ z n ∼ 0 (|σ z n | ∼ 1) indicates strong (weak) interlayer hybridization.Eq. ( 2) dictates that a large G, thus α H , requires strong interlayer hybridization and small interband energy separation.In Fig. 1(b), the first two energy bands are close to each other, and interlayer hybridization is strong except around the Brillouin zone corners.Consequently, one observes large conductivity peaks at low temperatures within the energy window of these two bands [red and blue curves, Fig. 1(b)].In contrast, the third energy band contributes much smaller conductivity due to the large separation to other bands, although the states are also strongly interlayer hybridized.
The conductivity has a strong dependence on the temperature.At low temperatures, ∂ ε f /T in the integrand of Eq. ( 7) is sharply peaked with large magnitude around the Fermi level, while it becomes broadened with reduced magnitude when the temperature increases.To illustrate this temperature modulation, Fig. 1(d) presents the density of α H in the momentum space at 0.4 K and 2 K when the Fermi level is at −2.2 meV [black arrow in Fig. 1(b)].Clearly, the density becomes more widely distributed as the temperature increases, which result in broadened conductivity peaks [red vs blue curves, Fig. 1(b)].For temperatures O(10) K, the conductivity remains observable and exhibits strong thermal smearing features [yellow curve, Fig. 1(b)].
Fig. 1(c) shows the results of twisted bilayer MoTe 2 at another twist angle of 1.8 • .As the second energy band moves closer to the third, the most pronounced conductivity peaks now locate in the energy window that covers these two bands.All the other features remain qualitatively the same as those of 1.4 • .

IV. APPLICATION TO TWISTED HOMOTRILAYER TMDS
Next we consider twisted homotrilayer TMDs.Specifically, we assume that the top and bottom layers are parallel, while the middle layer is misorientated by θ.We also assume that the top layer can be translated with respect to the bottom layer by δ 0 from the fully aligned configuration [see Fig. 2(a)].Such a trilayer geometry shares the same moiré period with the bilayer case.
The Hamiltonian for the twisted trilayer can be built from that of the bilayer [29]: where the subscripts t, m and b refer to the top, middle and bottom layer, δ ′ = θẑ × r + δ 0 (δ = θẑ × r) is the local displacement vector at position r between the top-middle (middle-bottom) layer.T l with l = t, m, b is given in Eq. ( 10), in which l = +1 for t and b, and l = −1 for m in the exponential terms.∆ T is given by Eq. ( 11).The layer-dependent electrostatic potential Q l is built with Eq. ( 12) as [39]: In addition, to evaluate the Nernst conductivity, the interlayer charge dipole operator in Eq. ( 2) is replaced by p = ed 0 diag(1, 0, −1) in the trilayer geometry [29].
When the top and bottom layers are aligned with δ 0 = 0, the structure has mirror symmetry in the outof-plane direction, thus the chiral Nernst effect is forbidden.As the top layer is translated, the mirror symmetry is broken and the Nernst effect is turned on.Fig. 2(b) shows the modulation of the valence bands and layer hybridization by δ 0 in the K valley of 3.5 • twisted trilayer MoTe 2 .For example, when δ 0 is tuned from 0.1(a 1 + a 2 ) to 0.2(a 1 + a 2 ) [see Fig. 2(a) inset for a 1 and a 2 ], the second valence band moves closer to the first, and the energy separations between the second and third bands are enlarged in general [Fig.2(b1) vs (b2)].For larger δ 0 , the separation between the first two bands can be even smaller at certain locations of the momentum space [Figs.2(b3) and (b4)].When δ 0 is small, one notices that the low-energy states are strongly polarized to the middle layer, see e.g., the dominant blue circles for |ε n | < 5 meV in Fig. 2(b1).As the translation increases, the lowenergy states become strongly layer hybridized, as shown in Figs.2(b2)-(b4).Such modulations of energy and layer hybridization affect the Nernst conductivity.For example, the large energy separations and weak layer hybridization at low energies in Fig. 2(b1) dictates that the corresponding conductivity is negligible, which is confirmed by the cyan curve in Fig. 2(c).In contrast, remarkably strong conductivity peaks can be found with a Fermi level around 10 meV for the intermediate translation of δ 0 = 0.2(a 1 + a 2 ) [orange curve, Fig. 2(c)], and large conductivity peaks can be achieved with very low doping near the topmost band edges for large translations of δ 0 = (a 1 + a 2 )/3 and δ 0 = 0.4(a 1 + a 2 ) [green and yellow curves, Fig. 2(c)].

V. DISCUSSION
We have proposed a novel effect, the dynamical chiral Nernst effect, which is unique to 2D vdW materials with chiral lattice structures.This is also the first Nernst response of an intrinsic nature determined solely by the band structures in nonmagnetic materials.We revealed that its quantum geometric origin is rooted in the interlayer coherence of electronic states endowed by strong interlayer quantum tunneling.The effect is shown to be sizable in typical twisted bilayer and trilayer TMDs, and can be feasibly tuned by twist angle and interlayer translation.
Compared to the intrinsic anomalous Nernst effect in magnetic systems caused by k -space Berry curvature, the intrinsic chiral Nernst effect in our work does not require magnetic field or internal magnetization, and shows similar or even larger magnitudes that are well within experimental capacity.Our calculations showed that the chiral Nernst conductivity can reach ∼ 0.5α 0 in twisted bilayer TMDs [i.e., 2.4 A/(K• m) with the thickness of bilayer MoTe 2 ] and trilayer TMDs.In comparison, the intrinsic anomalous Nernst conductivity can reach ∼ 3 × 10 −4 α 0 in monolayer TMDs placed on a magnetic insulating substrate [40], and the experiment on magnetic topological semimetal Co 2 MnGa reported anomalous Nernst conductivities around 3 A/(K• m) [6].This work opens a new route towards in-plane thermoelectric conversion by out-of-plane dynamical control of layered vdW structures [41].The study can be gen-eralized to chiral thermal Hall effects and thermoelectric Hall transport of spin and valley degrees of freedom.

FIG. 1 .
FIG. 1.(a) Schematics of the experimental setup for the dynamical chiral Nernst effect.(b, c) Energy bands and Nernst conductivity αH at various temperatures from the K valley at twist angle of 1.4 • (b) and 1.8 • (c).α0 = ekB/h ≈ 3.34 nA/K.(d) The density of Nernst conductivity in the momentum space at two different temperatures when the Fermi level is −2.2 meV [black arrow in (b)] and θ = 1.4 • .