
3
n, [G]ab(k) = 2eRe{Pn′̸=nAnn′
a(k)An′n
b(k)/[ϵn(k)−
ϵn′(k)]}, with Ann′
b=⟨un(k)|i∂kbun′(k)⟩. Here, a, b de-
notes Cartesian coordinates and G(k) is even under PT.
For clarity, we concentrate on the intraband limit and
focus on scalar impurities that preserve PT symmetry
and conserve spin. The distribution function in Eq. (4)
can be directly computed via a spatially uniform kinetic
equation, ∂tfs(k, t)−eE(t)·∂kfs(k, t)/¯h=I{fs(k, t)}
where [27,32,33]
I{fs(k, t)}=X
k′Ws
k′→kfs(k′, t)−Ws
k→k′fs(k, t)(6)
describes the spin-dependent collision integral.
The distribution function can be solved in the standard
perturbative fashion: in powers of Eand for weak skew-
scattering by using the relaxation time approximation,
see SI for a detailed derivation [28]. As such, we expand
fs(k, t) as
fs(k, t) = f0(k) + X
ℓ,m
f(m)
ℓ,s (k, t),(7)
where the second term captures the deviation of the dis-
tribution function from the equilibrium distribution func-
tion, f0(k). Here subscript ℓ= 1,2,· · · denote order in E
and the superscript m= 0,1,· · · denote its dependence
on skew-scattering rate; m= 0 captures the purely sym-
metric part independent of skew-scattering.
Note f0(k) is the same for both ↑and ↓due to PT
symmetry. Similarly, even as m= 0 contributions
depend on the (transport) relaxation time: (τs)−1=
⟨Pk′wS
s,k′,k(1 −cos θvv′)⟩,PT symmetry in Eq. (2) en-
sure f(0)
ℓ,s (k, t) are the same for ↑and ↓since τ↑=τ↓=τ.
Here θvv′is the angle between v(k) and v(k′), and ⟨· · · ⟩
indicates an average over an energy contour. In contrast,
skew-scattering (m= 1) contributions to the distribution
function f(1)
ℓ,s (k, t) have opposite signs for opposite spins,
see Fig. 1b: a property key to ASN.
Writing E(t) = Eeiωt +c.c. and substituting the dis-
tribution functions into Eq. (4) enables to directly dis-
cern the nonlinear Hall responses. Amongst the possi-
ble second-order nonlinear Hall responses obtained (see
Table I), two are PT-even; the rest are odd. The first
PT-even response is the intrinsic nonlinear Hall (INH)
effect [16,17,31] obtained by combining the second
term in Eq. (5) with f0(k). This yields an INH current
[jINH]a(t) = Re(j0
a+j2ω
aei2ωt) with j0
a=χINH
abc [Eb]∗Ec
and j2ω
a=χINH
abc EbEc, where χINH
abc [16,17,31] depends
only on band geometric quantities. χINH
abc is independent
of τand insensitive to ωin the semiclassical limit.
The second PT-even nonlinear Hall response, ASN, is
the main result of our work. This nonlinear Hall effect
arises from combining Ωs(k)×E(t) with the skew dis-
tribution function f(1)
1,s (k, t). This produces a nonlin-
ear Hall response: [jASN]a(t) = Re(j0
a+j2ω
aei2ωt) with
j0
a=χASN
abc [Eb]∗Ecand j2ω
a=χASN
abc EbEcwith
χASN
abc = 2e3εadb
¯h2X
k,k′,s
Ωs
d(k)˜τ2
ωwA
s,k,k′∂f0(k′)
∂k′c
,(8)
where Ωs
d(k) denotes the dcomponent of Ωs(k), ˜τω=
τ/(1 + iωτ) and εadb is the Levi-Civita symbol. Since
Ωs
dand wA
sare both odd under PT , their product is
even producing a finite extrinsic nonlinear Hall effect.
Importantly, χASN
abc scales as τ2wAfor ωτ ≪1. As a
result, χASN
abc is expected to dominate the nonlinear Hall
response in clean systems. At finite ω,χASN
abc displays a
characteristic ωdependence varying rapidly on the scale
1/τ (see below); this ωdependence distinguishes it from
both the ωinsensitive χINH as well as interband effects
that have characteristic ωdependence on the scale of
interband transition energy ϵn−ϵm.
Symmetry, scattering, and chiral photocurrents: ASN
has several striking attributes. Due to its Berry curva-
ture roots, χASN
abc is antisymmetric in its first two indices
yielding a nonlinear Hall effect [34] always transverse to
the applied electric field. This antisymmetric nature
imposes additional point-group symmetry constraints
as compared to conventional skew-scattering nonlineari-
ties [7,11,13]. For example, in 2D, antisymmetric non-
linear χabc requires broken rotational symmetry [3,17].
ASN’s antisymmetric behavior contrasts with that of
another PT-even nonlinear response that arises from
combining v(k) with f(0)
2,s (k, t) [16,18] to produce a clas-
sical nonlinearity, χDrude
abc . Importantly, χDrude
abc has a sus-
ceptibilty that is completely symmetric when its indices
are permuted yielding a response that need not always be
transverse as required of Hall type responses [34]. Exper-
imentally, this fully symmetric nonlinear Drude response
can be weeded out via interchanging the directions of
driving field and response: symmetric χabc is even under
exchange, whereas nonlinear Hall responses are odd.
Perhaps most striking is how ASN produces a
helicity-dependent chiral photocurrent: [j⟲]a=
i
2Im[χabc](E∗
bEc−EbE∗
c).ASN chiral photocurrent arises
from its part quantum geometric and part skew-
scattering origins. First, since ASN depends on skew
scattering χASN
abc possesses both real and imaginary com-
ponents arising from the complex valued ˜τ2
ωin Eq. (8).
Second, because ASN proceeds from the anomalous ve-
locity Ω×E, its susceptibility is asymmetric allowing for
a non-zero j⟲after both band cindices are summed.
Importantly, ASN’s combination of geometric nature
and scattering processes is essential. For instance, even
as symmetric scattering alone enables a nonlinear Drude
conductivity χDrude
abc [16,18] that has an imaginary com-
ponent, it nevertheless is completely symmetric under
any interchange of indices yielding a zero j⟲. Similarly,
while χINH
abc is also asymmetric, it nevertheless is purely
real producing helicity blind photocurrents. As a result,
to our knowledge, χASN
abc is the only intraband nonlinearity