---
title: Nonlinear Valley Nernst Effect (NVNE)
url: https://www.emergentmind.com/topics/nonlinear-valley-nernst-effect-nvne
type: topic
---

# Nonlinear Valley Nernst Effect (NVNE)

The nonlinear valley Nernst effect (NVNE) is a transverse, second-order thermoelectric response in multivalley crystals, especially two-dimensional systems with inequivalent \(K\) and \(K'\) valleys, in which a temperature gradient generates a valley current or, equivalently in some formulations, a valley-contrasting transverse magnetization current that scales quadratically with the thermal drive. Current work distinguishes an intrinsic NVNE, independent of the relaxation time, from Berry-curvature-driven and Drude-like nonlinear thermoelectric terms. In the 2025 literature, the intrinsic effect is formulated either through the quantum metric and temperature-gradient-induced orbital magnetization, or through the Berry connection polarizability dipole; in both cases it survives in settings where linear valley Nernst or valley Hall responses are symmetry-forbidden [2501.15835] [2509.01961] [2508.19586].

## 1. Definition and response structure

A standard valley-current definition is
\[
\mathbf{J}^v \equiv \mathbf{J}^{\mathrm K}-\mathbf{J}^{\mathrm K'},
\]
and the second-order valley thermoelectric response can be written as
\[
j_i^{v}=\alpha^{v(2)}_{ijk}\,(\nabla_jT)(\nabla_kT),
\]
or, in the thermal-field notation \(E^T=-\nabla T/T\),
\[
j_\beta^{(2)}=\alpha_{\beta;\gamma\delta}^T\,E^T_\gamma E^T_\delta.
\]
The defining feature is that the transverse response is valley-odd: the valley-resolved component changes sign between \(\tau=\pm1\), corresponding to \(K\) and \(K'\) [2509.01961] [2501.15835].

NVNE is distinct from the linear valley Hall effect (LVHE) and linear valley Nernst effect (LVNE). LVHE is a first-order transverse valley current driven by an electric field and controlled by Berry curvature \(\Omega_n(\mathbf{k})\). LVNE is the thermal analogue, likewise first order and Berry-curvature-based. In PT-symmetric systems, \(\Omega(\mathbf{k})=0\) identically, so those linear responses vanish, whereas NVNE can remain finite because its intrinsic part is built from PT-even band-geometric quantities [2501.15835].

| Response | Driving field | Geometric control |
|---|---|---|
| LVHE | Electric field | Berry curvature \(\Omega_n(\mathbf{k})\) |
| LVNE | \(\nabla T\) | Berry curvature \(\Omega_n(\mathbf{k})\) or orbital angular momentum |
| NVNE | \((\nabla T)^2\) or \(E^T E^T\) | Quantum metric, thermoelectric orbital magnetization, or Berry connection polarizability dipole |

In strained bilayer graphene, the response is explicitly Nernst-type: with strain along \(x\), the surviving element is
\[
\alpha^{v(2)}_{xyy}\neq0,\qquad \alpha^{v(2)}_{yxx}=0,
\]
so a temperature gradient along \(y\) induces a transverse valley current along \(x\) [2509.01961].

## 2. Intrinsic band geometry and microscopic mechanisms

One 2025 formulation identifies the intrinsic NVNE with the quantum metric \(g_{ij}\), the symmetric part of the quantum geometric tensor. For band \(n\),
\[
g_{ij}^{(n,\mathbf{k})}
=\operatorname{Re}\langle \partial_i u_n|[1-|u_n\rangle\langle u_n|]|\partial_j u_n\rangle
=\operatorname{Re}\sum_{m\neq n}\mathcal{R}_{nm,i}\mathcal{R}_{mn,j},
\]
with \(\mathcal{R}_{mn,i}=i\langle u_m|\partial_i u_n\rangle\). In the quantum-kinetic treatment, the second-order current takes the form
\[
\mathbf{j}=-e\,\mathrm{Tr}[\hat v\,\rho]+\sum_q \mathbf{E}^T\times \mathbf{M}^\Omega_q,
\]
and solving the Liouville equation in the relaxation-time approximation yields an intrinsic, \(\tau_s\)-independent conductivity tensor \(\alpha_{\beta;\gamma\delta}^T\) expressible through band-resolved quantum-metric objects \(\widetilde{\mathcal G}\) and \(\mathcal G\) [2501.15835].

A closely related semiclassical formulation writes the intrinsic nonlinear thermoelectric tensor as
\[
\alpha^{\mathrm{in}}_{abc}
=\frac{e}{\hbar}\sum_n\int[d\mathbf{k}]\,
\frac{(\varepsilon^n_{\mathbf{k}}-E_f)^2}{T^2}
\frac{\partial f_0}{\partial \varepsilon^n_{\mathbf{k}}}
\,\chi^n_{abc}(\mathbf{k}),
\]
with
\[
\chi^n_{abc}(\mathbf{k})
=\sum_{m\neq n}\hbar
\left[
\frac{2v_a^n\mathcal G^{mn}_{bc}-v_b^n\mathcal G^{mn}_{ac}-v_c^n\mathcal G^{mn}_{ab}}
{\varepsilon^n_{\mathbf{k}}-\varepsilon^m_{\mathbf{k}}}
\right],
\]
and
\[
\mathcal G^{mn}_{ab}(\mathbf{k})
=\operatorname{Re}\!\left[\mathcal A_a^{nm}(\mathbf{k})\mathcal A_b^{mn}(\mathbf{k})\right],\qquad
g^n_{ab}(\mathbf{k})=\sum_{m\neq n}\mathcal G^{mn}_{ab}(\mathbf{k}).
\]
This makes the quantum-metric origin explicit and separates the intrinsic contribution from Berry-curvature and Drude terms, which scale as \(\tau\) and \(\tau^2\), respectively [2509.01961].

A second 2025 formulation expresses intrinsic NVNE through the Berry connection polarizability dipole (BCPD). The gauge-invariant tensor
\[
G_{ab}^n(\mathbf{k})
=
2\,\operatorname{Re}\sum_{m\neq n}
\frac{\mathcal A_a^{nm}\mathcal A_b^{mn}}
{\varepsilon_n-\varepsilon_m}
\]
enters the intrinsic valley tensor as
\[
\beta^v_{abc}
=
2\int_K[d\mathbf{k}]\,
\big(v_a^nG_{bc}^n-v_b^nG_{ac}^n\big)\,
\frac{(\varepsilon_n-\mu)^2}{T^2}\,
f_0'(\varepsilon_n).
\]
This formulation emphasizes that the intrinsic response is a Fermi-surface property of band quantum geometry and is independent of scattering [2508.19586].

The thermoelectric orbital-magnetization channel supplies the valley discriminator in PT-symmetric systems. The temperature-gradient-induced orbital magnetic moment \(m^{\nabla T}_{n,a}\) in strained bilayer graphene and the orbital magnetic moment \(m_{q,\beta}^{T;o}\) in the PT-symmetric quantum-kinetic treatment both acquire opposite signs in the two valleys. In that sense, the thermoelectric correction to orbital magnetization plays a role analogous to orbital angular momentum in the linear valley Hall picture [2509.01961] [2501.15835].

At low temperature, nonlinear Mott-type relations connect thermal and electric second-order tensors. One relation is
\[
\alpha^{\mathrm{in}}_{abc}(E_f,T)\simeq
\frac{\pi^2 k_B^2T^2}{3e}\,
\frac{\partial}{\partial E_f}\chi^{\mathrm{in}}_{abc}(E_f),
\]
while another is
\[
\beta^v_{abc}=\mathcal L\,\chi^v_{abc},\qquad
\mathcal L=\frac{\pi^2k_B^2}{3e^2}.
\]
These results show that the intrinsic thermal response is systematically linked to nonlinear valley Hall transport, although the two papers emphasize different geometric kernels [2509.01961] [2508.19586].

## 3. Symmetry conditions and tensor selection rules

A central symmetry result is that NVNE is allowed under conditions that suppress linear valley responses. In PT-symmetric systems,
\[
T[\Omega(\mathbf{k})]=-\Omega(-\mathbf{k}),\qquad
P[\Omega(\mathbf{k})]=-\Omega(-\mathbf{k}),
\]
so PT enforces \(\Omega(\mathbf{k})=0\), killing LVHE and LVNE. By contrast, the quantum metric satisfies
\[
T[\mathcal G(\mathbf{k})]=\mathcal G(-\mathbf{k}),\qquad
P[\mathcal G(\mathbf{k})]=\mathcal G(-\mathbf{k}),
\]
hence remains PT-even and can support a finite intrinsic second-order response [2501.15835].

In globally inversion- and time-reversal-symmetric systems, the Berry-curvature contribution vanishes pointwise,
\[
\Omega^n_{ab}(\mathbf{k})=0\quad\Rightarrow\quad \alpha^{\mathrm B}_{abc}=0,
\]
but the intrinsic valley difference can remain finite because time reversal maps \(K\) to \(K'\) without forcing the single-valley tensor to vanish. Thus \(\alpha^{\mathrm{in},K}_{abc}\) and \(\alpha^{\mathrm{in},K'}_{abc}\) may be equal in magnitude and opposite in sign, so the total charge current cancels while a pure valley current survives [2509.01961].

Local rather than global symmetry controls the allowed tensor components. In two dimensions, the largest local symmetry compatible with a nonzero intrinsic NVNE is a single mirror. With strain along \(x\) in bilayer graphene, the surviving local symmetry is \(M_y\), and the only allowed intrinsic valley component is
\[
\alpha^{\mathrm{in,valley}}_{xyy}\neq0.
\]
Additional rotational symmetries such as \(C_3\) or \(C_n\) with \(n>2\) forbid the response [2509.01961].

A complementary selection-rule formulation uses the pseudotensor
\[
\mathcal X^v_{cd}\equiv \frac12\epsilon_{abc}\beta^v_{abd},
\]
which transforms under a point-group operation \(\mathcal O\) as
\[
\mathcal X^v\to \eta_v\,\det(\mathcal O)\,\mathcal O\,\mathcal X^v\,\mathcal O^{-1},
\]
with \(\eta_v=+1\) for valley-preserving and \(\eta_v=-1\) for valley-switching operations. In point group \(C_s\), this allows \(\beta^v_{xyy}=-\beta^v_{yxy}\), whereas three- and fourfold rotational symmetries typically forbid the same components [2508.19586].

Anisotropy or tilt is therefore not a minor perturbation but a symmetry-enabling ingredient. In the tilted Dirac model of the PT-symmetric study, \(\alpha_t\to0\) kills NVNE. In strained bilayer graphene, the temperature gradient must be applied perpendicular to the strain direction. These results establish NVNE as a symmetry-selected anisotropic transport coefficient rather than a generic consequence of valley degeneracy [2501.15835] [2509.01961].

## 4. Canonical models and material realizations

The prototypical PT-symmetric model is an anisotropic tilted Dirac semimetal with
\[
H_\tau(\mathbf q)=v_F(\tau\sigma_x q_x+\sigma_y q_y)+\tau\alpha_t q_y\sigma_0+\Delta\sigma_z.
\]
Its Berry curvature vanishes under PT symmetry, but the band-resolved quantum metric remains finite. For the valence band,
\[
\begin{bmatrix}
\widetilde{\mathcal G}^{vc}_{xx} & \widetilde{\mathcal G}^{vc}_{xy}\\
\widetilde{\mathcal G}^{vc}_{yx} & \widetilde{\mathcal G}^{vc}_{yy}
\end{bmatrix}
=
-\frac{v_F^2}{8\varepsilon^5}
\begin{bmatrix}
v_F^2q_y^2+\Delta^2 & -v_F^2q_xq_y\\
-v_F^2q_xq_y & v_F^2q_x^2+\Delta^2
\end{bmatrix}.
\]
The intrinsic NVNE tensor component is obtained analytically as
\[
\alpha_{y;xx}^{T,(2)}(\zeta)
=
\zeta\frac{e\tau\alpha_t}{64\pi\hbar T^2}
\left[
17+29\frac{\Delta^4}{\mu_c^4}-46\frac{\Delta^2}{\mu_c^2}
\right],
\]
with \(\zeta=\pm1\) for conduction and valence bands. In the gapless limit,
\[
\alpha_{y;xx,\max}^{T,(2)}(\zeta)
=
\zeta\frac{17e\alpha_t}{64\pi\hbar T^2},
\]
independent of \(\mu_c\). The valley-dependent thermoelectric orbital magnetic moment is
\[
\frac{m^{T,o}_{v,z}}{\nabla_xT}
=
-\frac{e v_F^4 q_y}{8\hbar \varepsilon_0^4}
\left[
1+2\tau\alpha_t
\frac{v_F^2q_y^2+\Delta^2}{\varepsilon_0 v_F^2 q_y}
\right],
\]
and the integrated valley orbital magnetization is
\[
M_v^T=
\frac{e\tau\alpha_t}{16\pi\hbar}
\left[
\frac1{\varepsilon_F}+\frac{3\Delta^2}{\varepsilon_F^3}
\right]\nabla_xT.
\]
These expressions make the valley sign, tilt dependence, and band dependence explicit [2501.15835].

The same paper extends the analysis to PT-symmetric bilayer graphene with trigonal warping and to the organic conductor \(\alpha(\mathrm{BEDT\!-\!TTF})_2\mathrm I_3\). In bilayer graphene, integrating the thermoelectric orbital magnetization yields
\[
M_z^T(K)\approx -444.5\times (e\nabla_yT/\hbar),\qquad
M_z^T(K')\approx +444.5\times (e\nabla_yT/\hbar),
\]
while in \(\alpha(\mathrm{BEDT\!-\!TTF})_2\mathrm I_3\),
\[
M_v^T(K)\approx +28.14\times (e\nabla_yT/\hbar),\qquad
M_v^T(K')\approx -28.14\times (e\nabla_yT/\hbar).
\]
In both materials the sign reverses between valleys, and in the organic compound the sign of the NVNE also reverses between valence and conduction bands [2501.15835].

A distinct but related realization is uniaxially strained AB-stacked bilayer graphene, described by
\[
\hat H_{\tau_v}(\mathbf k)=N_1(\mathbf k)\sigma_x+N_2(\mathbf k)\sigma_y,
\]
with
\[
N_1=-\frac{\hbar^2}{2m^*}(k_x^2-k_y^2)+\tau_v v_3\hbar k_x+w,\qquad
N_2=-\tau_v\frac{\hbar^2}{m^*}k_xk_y-v_3\hbar k_y.
\]
For strain along \(x\), only \(\alpha^{\mathrm{in,valley}}_{xyy}\) survives, so \(\nabla T\parallel \hat y\) is required to generate \(j_x^v\). The intrinsic coefficients satisfy
\[
\alpha^{\mathrm{in},K}_{xyy}=-\alpha^{\mathrm{in},K'}_{xyy},
\]
leaving zero charge current and a pure valley current. Changing from compressive (\(w<0\)) to tensile (\(w>0\)) strain reverses the sign of \(\alpha^{v(2)}_{xyy}\), with peaks near Lifshitz transitions around \(w\simeq -1\,\mathrm{meV}\) and \(w\simeq +3\,\mathrm{meV}\) [2509.01961].

First-principles work on bilayer \(\mathrm{WTe}_2\) provides an inversion-asymmetric realization. With point group \(C_s\), two well-separated valleys on the \(k_x\) axis, and \(\beta^v_{xyy}=-\beta^v_{yxy}\) allowed by symmetry, the intrinsic NVNE component \(\beta^v_{xyy}\) peaks near \(\mu\approx0.09\,\mathrm{eV}\) with magnitude on the order of \(\mathrm{nA\cdot nm/K^2}\) at \(T=10\,\mathrm K\) [2508.19586].

## 5. Scaling laws, tunability, and experimental signatures

Parameter dependence is unusually transparent in the tilted Dirac and strained bilayer graphene models. In the PT-symmetric tilted Dirac case, the response grows with the tilt or anisotropy parameter \(\alpha_t\), vanishes for \(\alpha_t\to0\), decreases as the gap \(\Delta\) increases, and changes sign between conduction and valence bands through the factor \(\zeta=\pm1\). For finite \(\Delta\), the analytical expression depends on \(\Delta^2/\mu_c^2\) and \(\Delta^4/\mu_c^4\); in the gapless limit it becomes independent of \(\mu_c\) [2501.15835].

Temperature scaling requires care because different papers normalize the response differently. In the \(E^T=-\nabla T/T\) convention, the intrinsic coefficient in the PT-symmetric Dirac model scales as \(1/T^2\). In the \((\nabla T)^2\) convention used for strained bilayer graphene and for the BCPD theory, the low-temperature intrinsic coefficients obey Sommerfeld-type relations and scale as \(T^2\). This difference reflects tensor conventions rather than a contradiction in the existence of the effect [2501.15835] [2509.01961] [2508.19586].

Several experimental observables have been proposed. Direct transverse thermovoltage measurements apply \(\nabla T\) longitudinally and read out a transverse signal; a nonzero response in globally PT-symmetric samples, where linear Nernst transport is absent, is a primary signature. Nonlocal valley transport uses valley filters or injectors and tests sign reversal under valley polarization switching. Scanning SQUID or NV-center magnetometry can image edge magnetization currents, exploiting \(j_{\mathrm{mag}}=\nabla\times M\). Kerr rotation, circular dichroism, and nonreciprocal directional dichroism have been proposed as valley-sensitive optical probes of edge accumulation or magnetic-quadrupole structure [2501.15835] [2509.01961] [2508.19586].

The nonlocal theory for intrinsic NVNE predicts a second-harmonic signal with characteristic \(\rho^2\) scaling,
\[
\frac{\delta V_{2\omega}}{(\Delta T)^2}=\rho^2(c_1+c_2T^2),
\]
and a ratio
\[
\frac{c_2T^2}{c_1}=\frac{Q^2}{\mathcal L},
\]
where \(Q\) is the thermopower and \(\mathcal L\) is the Lorenz number. In bilayer \(\mathrm{WTe}_2\), using \(\rho\approx10^4\,\Omega\), \(\sigma^v\approx6\,\mu\mathrm S\), \(x\sim \ell_v\sim1\,\mu\mathrm m\), and \(40\,\mathrm{K}/\mu\mathrm m\), the estimated second-harmonic nonlocal voltage is
\[
\delta V_{2\omega}\sim0.01\,\mu\mathrm V,
\]
and at \(\mu\approx0.09\,\mathrm{eV}\) the ratio \(c_2T^2/c_1\) is \(\sim0.001\), so the signal is dominated by the direct NVNE pathway [2508.19586].

The PT-symmetric study gives order-of-magnitude estimates in a different language. With \(\alpha_t\) of order \(0.1\)–\(0.3\,\mathrm{eV\cdot \AA}\) and \(T\approx50\)–\(150\,\mathrm K\), the gapless NVNE coefficient can generate measurable transverse voltages under strong thermal gradients. For \(\nabla T\approx1\)–\(10\,\mathrm{K/mm}\) and device widths of \(1\)–\(10\,\mu\mathrm m\), the resulting magnetization current density is stated to correspond to microampere-per-meter scale transverse currents and microvolt-to-millivolt thermovoltages, although exact SI conversion depends on band parameters and geometry [2501.15835].

## 6. Relation to adjacent nonlinear Nernst effects and current limitations

A recurring misconception is that all nonlinear Nernst effects are equivalent. Earlier work on strained transition-metal dichalcogenides analyzed a nonlinear anomalous Nernst current in time-reversal-symmetric but inversion-broken systems, controlled by Berry curvature near the Fermi surface and by the pseudovector
\[
\Lambda_d^T=-\frac1{T^2}\int[d\mathbf k]\,
\frac{\partial f_0}{\partial k_d}
(E_{\mathbf k}-\mu_e)^2\Omega_z(\mathbf k).
\]
There the valley mechanism is explicit, but the two valley contributions add to a net transverse charge current rather than a pure valley current [1903.00814].

Likewise, the 2019 study of bilayer \(\mathrm{WTe}_2\) treated the nonlinear anomalous Nernst effect in an inversion-broken, time-reversal-symmetric material and found a Berry-curvature-controlled second-order charge current proportional to a Fermi-surface quantity \(\Lambda^T\). Valley-resolved contributions can be defined, but the reported observable is the net charge response, not a symmetry-protected pure valley current [1905.09814].

The 2020 quantum-kinetic treatment introduced a thermoelectric Berry-curvature dipole \(D_i^T\) and showed that, in inversion-broken crystals, second-order thermoelectric transport may contain two intrinsic Berry-curvature-linked terms. That analysis also stated that charge and pure valley responses depend on whether the relevant anisotropy or tilt is valley-odd or valley-even: valley-odd tilt tends to make charge currents add, whereas valley-even deformations can leave \(j_x=0\) but \(j_x^v\neq0\) [2010.15340]. This suggests that the modern NVNE literature sits at the intersection of two threads: Berry-curvature-based nonlinear anomalous thermoelectric transport in inversion-broken systems, and quantum-metric-based or BCPD-based pure valley transport that survives global inversion and time-reversal symmetry.

The present theories also share limitations. The PT-symmetric quantum-kinetic and semiclassical treatments focus on the intrinsic, clean-limit contribution, neglecting disorder-dependent skew-scattering and side-jump corrections in the main result. Strong intervalley scattering can reduce the net NVNE by relaxing valley imbalance. The strained bilayer-graphene analysis assumes weak, slowly varying temperature gradients, a low-energy \(k\cdot p\) regime, and no intervalley scattering. The first-principles nonlocal theory similarly assumes long intervalley diffusion length and uses transport-scaling arguments to separate intrinsic and extrinsic terms [2501.15835] [2509.01961] [2508.19586].

Taken together, these works define NVNE as a genuine nonlinear valley caloritronic effect rather than a simple extension of linear Nernst physics. Its most distinctive forms are the intrinsic, relaxation-time-independent responses permitted by band quantum geometry, selected by local mirror symmetry, amplified by anisotropy or trigonal warping, and observable even when Berry curvature and all linear valley Hall or Nernst signals vanish globally.

Source: https://www.emergentmind.com/topics/nonlinear-valley-nernst-effect-nvne