The Berry Phase Rectification Tensor and The Solar Rectification Vector

^{1}

Max Planck Institute for the Physics of Complex Systems, Dresden 01187, Germany

^{2}

Max Planck Institute for Chemical Physics of Solids, Dresden 01187, Germany and

^{3}

Centre for Theoretical Studies, Indian Institute of Technology Kharagpur, Kharagpur 721302, India

We introduce an operational definition of the Berry Phase Rectification Tensor as the second order change of polarization of a material in response to an ideal short pulse of electric field. Under time reversal symmetry this tensor depends exclusively on the Berry phases of the Bloch bands and not on their energy dispersions, making it an intrinsic property to each material which contains contributions from both the inter-band shift currents and the intra-band Berry Curvature Dipole. We also introduce the Solar Rectification Vector as a technologically relevant figure of merit for bulk photo-current generation which counts the number of electrons contributing to the rectified current per incoming photon under ideal black-body radiation in analogy with the classic solar cell model of Shockley and Queisser. We perform first principle calculations of the Berry Phase Rectification Tensor and the Solar Rectification Vector for the Weyl semi-metal TaAs and the insulator LiAsSe which features large shift currents close to the peak of solar radiation intensity. We also generalize the formula for the Glass coefficient to include the spectral distribution of the incoming radiation, the directionality dependence of the conductivity of the material and the reflectivity at its surface.

PACS numbers:

72.15 . - v, 72.20 . M y, 73.43 . - f, 03.65 . V f

1. INTRODUCTION

Several remarkable relations between the quantum geometry and topology of electronic states in crystalline solids and their non-linear optical and transport properties have been unraveled over the last few decades. One of the second order effects intimately intertwined with the geometry of electronic wave-functions is the shift current [1-10], which is a rectified current arising in response to AC electric fields that drive optical transitions in crystals without inversion symmetry. The shift current arises from band-off-diagonal components of the velocity operator, in contrast to the injection current, which arises from the band-diagonal difference of group velocities between empty and occupied states

[2, 3, 5 - 10]

. In materials with time reversal symmetry, the shift current differs from the injection current in that only the shift current can lead to rectified currents in response to linearly polarized electric fields, making it attractive for photovoltaic applications

[6, 11 - 14]

since it can lead to a DC current for electric fields with no net degree of circular polarization such as those arriving from the sun.

Another second order effect with an intimate connection to the geometry of electronic wave functions is the non-linear Hall effect driven by the Berry Curvature Dipole (BCD)

[10, 15 - 17]

, that has been recently observed experimentally [18-23]. Like the shift current, the BCD non-linear Hall effect can also give rise to a rectified current in response to a linearly polarized electric field, but with the crucial difference that this effect occurs only in metals at low frequencies in a Drude-like fashion [10]. These two distinct effects, however, appear to have a connection, as suggested by the quantum-rectification sum rule (QRSR) derived in Ref.[10]. The QRSR states that the frequency integral of the total rectification conductivity for any time-reversal-invariant material is com-

FIG. 1: The Berry Phase Rectification Tensor measures the second order change of polarization (which equals the area in the right panel) in response to delta function pulse of electric field (left panel).

pletely independent of energy dispersions and entirely controlled by the quantum geometric Berry connections of the bands. The contributions to the QRSR arise from the BCD effect, which accounts for all the intra-band weight in a Drude-like peak, and the shift current, which accounts for all the inter-band contributions. In contrast, injection currents do not contribute to the QRSR under time reversal invariant conditions.

In this work, we will show that the QRSR defines a three-index Berry Phase Rectification Tensor (BPRT), which depends only on the Berry phase geometry of a time-reversal-invariant material. In d-dimensions the BPRT has units of (length)

^{d - 3}

, and therefore, it is dimensionless for 3D materials. We will introduce an operational definition of the BPRT by considering the change of polarization in response to ideal short pulses of electric fields (see Fig.1) and compute it from first principle calculations for the Weyl semi-metal TaAs [24-27] and demonstrate that, remarkably, the net contribution from the BCD [28] is comparable to the inter-band contribution from shift currents. We will also compute the BPRT for the insulator

{L i A s S e}_{2}

, which is a promising photo-voltaic

FIG. 2: Crystal structures and Solar Rectification Vectors (SRV) for (a)

{L i A s S e}_{2}

and (b) TaAs.

material [6]. Additionally we will introduce the notion of the Solar Rectification Vector (SRV), whose magnitude serves as a technologically useful figure of merit for bulk photo-current generation, and is defined from the average DC current produced under ideal black-body solar radiation, analogous to the model employed by Shockley and Quessier (SQ) in their seminal study of

p - n

junction solar cells

[29]

. We will demonstrate that

{L i A s S e}_{2}

could generate photo-currents as large as

5 m A / {c m}^{2}

, which is about

5 %

of the ideal maximum from the

S Q

model [29]. The direction of the SRV for these materials is shown as a yellow vector in Fig.[2].

As we will illustrate, one of the appealing features of our formulas for computing the BPRT and the SRV, is that they can be readily combined with realistic first principle calculations of material band structures. Therefore, these formulas have the potential to be used to implement a more comprehensive and systematic search of materials with technologically desirable bulk photovoltaic properties.

2. THE BERRY PHASE RECTIFICATION TENSOR

Consider an electronic system which is in equilibrium for

t < 0

, and at

t = 0

is subjected to a quench of the vector potential, leading to a delta-function pulse in the electric field of the form:

A (t) = A_{0} Θ (t), E (t) = - A_{0} δ (t),

where

Θ (t)

is the heavy-side function. After the pulse, the net change of electric polarization can be obtained from:

Δ P = \int_{0}^{\infty} j (t) d t

This protocol is illustrated in Fig.[1]. To second order in electric fields the current is related to the linear and second order conductivity as:

\begin{aligned} j^{γ} (t) = \int_{- \infty}^{\infty} d t_{1} σ_{(1)}^{γ α} (t, t_{1}) E^{α} (t_{1}) + \\ + \int_{- \infty}^{\infty} d t_{1} \int_{- \infty}^{\infty} d t_{2} σ_{(2)}^{γ α β} (t, t_{1}, t_{2}) E^{β} (t_{1}) E^{α} (t_{2}) \end{aligned}

where Greek letters denote space indices and here and throughout the paper sum over repeated indices is understood unless otherwise stated. Therefore, the expansion of the polarization in powers of

A_{0}

is:

Δ P^{γ} = - σ^{γ α} A_{0}^{α} + \frac{e^{3}}{ℏ^{2}} R^{γ α β} A_{0}^{β} A_{0}^{α} + \dots

where

σ^{γ α}

is the zero frequency linear conductivity, and

R^{γ α β}

defines the Berry Phase Rectification Tensor (BPRT) of the system, which is given by (see Appendix A):

\frac{e^{3}}{ℏ^{2}} R^{γ α β} = \int_{- \infty}^{\infty} d t σ_{(2)}^{γ α β} (t, 0, 0) = \int_{- \infty}^{\infty} \frac{d ω}{2 π} σ_{(2)}^{γ α β} (- ω, ω)

Eq.(4) expresses the notion that the BRPT measures the net displacement of the average position of an electronic system to second order in the strength of delta function pulse in electric field. As demonstrated in Ref.[10] for the ground state of any non-interacting system with a time-reversal-invariant band structure,

R^{γ α β}

only contains information about the Berry connections of the bands. Specifically, in the clean limit of a time reversal invariant band structure, it is given by:

\begin{aligned} R^{γ α β} = \frac{1}{4} {⟨ \partial^{α} Ω^{β γ} ⟩ & + ⟨ [A^{α}, i \partial^{β} A^{γ}] ⟩ + \\ + ⟨ [A^{α}, [A^{β}, {\bar{A}}^{γ}]] ⟩ + (α \leftrightarrow β)} \end{aligned}

where

\partial^{α} \equiv \partial / \partial k^{α}, Ω_{n}^{α β} = \partial^{α} A_{n n}^{β} - \partial^{β} A_{n n}^{α}

is the abelian Berry curvature of the

n

-th band, and

A_{n m}^{α} = A_{n m}^{α} (1 -

δ_{n m})

is the off-diagonal non-Abelian Berry connection and average is defined as:

⟨ \dots ⟩ = \sum_{n} \int \frac{d^{d} k}{(2 π)^{d}} f_{n} ⟨ n | (\dots) | n ⟩

where

f_{n}

is the Fermi-Dirac occupation function of the

n

th Bloch band. The first term in Eq. (6) arises from the intra-band BCD while the two other terms arise from the inter-band shift current. Therefore the BPRT connects these two effects and allows for a concrete measure of the shift current effect, as a collective displacement of average position of the electron system in response to an electric field pulse. Such an interpretation of the BPRT suggests that time domain spectroscopic techniques

[30, 31]

are promising tools to measure the BPRT in materials. By construction the BPRT is symmetric in its last two indices, and the

B C D

does not contribute to fully diagonal components of the form

R^{α α α}

. For a related, but more restricted, sum rule see Ref.[32].

3. ${L i A s S e}_{2}$

FIG. 3: The non-zero conductivity components for LiAsSe as

_{2}

a function of frequency in the clean limit

(Γ = 0)

. The dashed line is the solar Planck intensity distribution (in arbitrary units) showing its good overlap with the

σ_{(2)}^{z z z}

rectification conductivity of LiAsSe2.

Notice that the intra-band contribution to the BPRT, namely the Berry dipole contribution proportional to

⟨ \partial^{α} Ω^{β γ} ⟩

in Eq. (6), is gauge invariant as it only depends on the gauge invariant Berry curvature. Since the full BPRT tensor is also gauge invariant [10], it follows that the remainder of the tensor is separately gauge invariant. This remainder contains the contributions of interband transitions arising from shift-current effects [10]. Therefore, there is no gauge ambiguity in separating the intra- and inter-band contributions to the BPRT, and we will compute them separately for the semimetal TaAs in Sec.V.

4. THE SOLAR RECTIFICATION VECTOR

In this section we introduce a figure of merit for bulk photovoltaic materials that captures their amount of DC current generated from sun-light under ideal conditions, and therefore it is a more relevant quantity for technological applications. Consider an statistical ensemble electric fields that is time-translationally invariant, whose average and two-time autocorrelation functions are given by:

\bar{E^{α} (t_{1})} = 0, \bar{E^{β} (t_{1}) E^{α} (t_{2})} = I^{β α} (t_{1} - t_{2})

Now, from Eq.(3), it follows that the ensemble averaged current is given by:

\bar{j^{γ} (t)} = \iint_{- \infty}^{\infty} d t_{1} d t_{2} σ_{(2)}^{γ α β} (t, t_{1}, t_{2}) I^{β α} (t_{1} - t_{2})

Then the average current will be time independent and given by:

\bar{j^{γ}} = \int_{- \infty}^{\infty} d ω σ_{(2)}^{γ α β} (- ω, ω) I^{β α} (ω)

Where

I^{β α} (ω)

is the Fourier transform of the electric field auto-correlation function. In the next subsection, we will compute such electric field auto-correlation function for the ensemble of solar black-body radiation.

5. A. Autocorrelations for black-body radiation

The photon eigenmode representation of the electric field operator in Heisenberg representation and residing in a periodic box of volume

V

is:

\begin{array}{r} E (r, t) = i \sum_{μ, k} \sqrt{\frac{ℏ ω_{k}}{2 ε_{0} V}} [e_{μ} (k) a_{k}^{μ} e^{- i (ω_{k} t - k \cdot r)} \\ - e_{μ}^{*} (k) a_{k}^{† μ} e^{i (ω_{k} t - k \cdot r)}] \end{array}

where

ω_{k} = c | k |, a_{k}^{μ}, a_{k}^{† μ}

are photon destruction and creation operators satisfying standard commutation relations,

μ = x, y

labels the directions of light polarization and the unit vectors satisfy:

\begin{matrix} e_{x} \cdot k = e_{y} \cdot k = e_{x} \cdot e_{y} = 0 \\ \sum_{μ} e_{μ}^{α} (k) e_{μ}^{β *} (k) = δ_{α β} - {\hat{k}}_{α} {\hat{k}}_{β} \end{matrix}

The quantum ensemble average of the photon number is given the by the Planck distribution:

⟨ a_{k}^{μ †} a_{k^{'}}^{μ^{'}} ⟩ = δ_{{k k}^{'}} δ_{μ μ^{'}} n_{k}, n_{k} = \frac{1}{e^{β ℏ ω_{k}} - 1},

where

β = 1 / k_{B} T

is the inverse temperature of the black body. To compute the auto-correlation function of the electric field we use the following prescription for quantum-to-classical conversion of the averages:

\bar{E^{α} (r, t) E^{β} (r^{'}, t^{'})} \leftrightarrow ⟨ : E^{α} (r, t) E^{β} (r^{'}, t^{'}) : ⟩

Where in the left-hand side we have the correlations denoted by an over-bar, in which the electric field is viewed as a classical field, and in the right-hand side we have quantum mechanical averages in which the electric field is viewed as the operator appearing in Eq.(11). We have added a normal ordering of the quantum electric fields which removes the vacuum Casimir-type fluctuations which should be negligible in the thermodynamic limit

V \to \infty

. With this we obtain:

\begin{aligned} \bar{E^{α} (r, t) E^{β} (r, t^{'})} = \\ = \frac{ℏ}{ε_{0} V} \sum_{k} (δ_{α β} - {\hat{k}}_{α} {\hat{k}}_{β}) ω_{k} n_{k} \cos [ω_{k} (t - t^{'})] \end{aligned}

In the above formula all modes of light contribute to produce a fully isotropic correlator of electric fields at a given point

r

. This would be the situation if the point of interest was fully immersed inside the isotropic blackbody radiation ensemble. However, when the radiation that arrives at

r

is highly directional we can modify the formula above by a one with an extra weighing function inside the momentum integral

f (\hat{k})

\begin{aligned} \bar{E^{α} (r, t) E^{β} (r, t^{'})} = \\ = \frac{ℏ}{ε_{0} V} \sum_{k} (δ_{α β} - {\hat{k}}_{α} {\hat{k}}_{β}) ω_{k} n_{k} \cos [ω_{k} (t - t^{'})] f (\hat{k}) \end{aligned}

which phenomenologically takes into account the preferred directionality of the incoming light at a given point r. In the special case in which the function

f (\hat{k}) = 1

only for a very small region of solid angles,

Δ Ω

, around a specific direction

\hat{k} = {\hat{k}}_{0}

, we have that:

\begin{aligned} \bar{E^{α} (ω_{1}) E^{β} (ω_{2})} = \\ = \frac{Δ Ω}{4 π} \frac{ℏ (δ_{α β} - {\hat{k}}_{0}^{α} {\hat{k}}_{0}^{β})}{(2 π)^{2} ε_{0} (β ℏ c)^{3}} \frac{{(ℏ β | ω_{1} |)}^{3}}{e^{h β | ω_{1} | - 1}} δ (ω_{1} + ω_{2}) \end{aligned}

For the light emitted from the sun and arriving on earth,

Δ Ω = 4 π {(R_{E} / R_{S})}^{2}

and

{\hat{k}}_{0}

would correspond to the unit vector defining the direction from sun to earth,

R_{S}

is the radius of the sun, and

R_{E}

is the distance of the sun to the earth.

The formula above is expected to describe the ideal solar radiation before entering the earth atmosphere, but the solar irradiance spectrum on the surface of earth can still be reasonably well approximated by it (see e.g. [33]). Upon entering the atmosphere the spectral intensity gets reduced at specific frequencies corresponding to absorption of molecules in the atmosphere. There will also be a randomization of the transversality of the electric field due to elastic scattering. This could be accounted for by averaging the projector over some distribution of

{\hat{k}}_{0}

, but for simplicity we will simply replace the transverse projector in Eq. (18) by a delta function with all components of the electric field, namely we take the electric field correlator to be:

\bar{E^{α} (ω_{1}) E^{β} (ω_{2})} = \frac{ℏ Δ Ω}{2 (2 π)^{3}} \frac{δ_{α β} δ (ω_{1} + ω_{2})}{ε_{0} (β ℏ c)^{3}} \frac{{(ℏ β | ω_{1} |)}^{3}}{e^{ℏ β | ω_{1} |} - 1}

As we will see in the subsection III B, as along as the material is oriented so that its "solar rectification vector" (to be defined in next subsection) is orthogonal to

{\hat{k}}_{0}

, one obtains the same current from the Eq.(19) as the one one would obtain from the Eq.(18) that takes into account the transverse nature of incident light.

Therefore, in conclusion, the electric field autocorrelation function associated with solar radiation, and which is needed in order to compute the average DC cur- rent in Eq.(10), can be written as:

\begin{matrix} I^{β α} (ω) = \frac{ℏ δ^{β α}}{(2 π)^{2} ε_{0} {(β_{S} ℏ c)}^{3}} {(\frac{R_{S}}{R_{E}})}^{2} I (ℏ β_{S} | ω |) \\ I (x) = \frac{x^{3}}{e^{x} - 1} \end{matrix}

where

I (x)

is the dimensionless Planck spectral distribution function,

β_{S} = 1 / (k_{B} T_{S})

is the inverse temperature of the

sun T_{S} \approx 5778 K, R_{S}

is the radius of the sun, and

R_{E}

is the sun-earth distance. As we will see in the next subsection, the geometric factor

{(R_{S} / R_{E})}^{2}

will disappear from the normalized figure of merit that measures the electric current normalized to the incoming photon flux. They are needed however if one would like to compute the actual magnitude of the induced photo-current.

6. B. Definition of Solar Rectification Vector

Eqs. (10) and (20) provide the total current density induced by exposing a material to solar radiation. To define a figure of merit for current generation, it is convenient to normalize such total current density (number of electrons per unit area per unit time) by the photon flux

F_{s}

(number of photons per unit area per unit time

[F_{s}] = m^{- 2} s^{- 1}

), so as to count how many current-carrying electrons are induced per incoming photon. Let us therefore also compute the photon flux on earth. The photon density,

n

is given by:

\begin{matrix} n^{μ} (r) = a^{μ †} (r) a^{μ} (r), \\ a^{μ} (r, t) = \frac{1}{\sqrt{V}} \sum_{k} a_{k}^{μ} e^{- i ω_{k} t + i k \cdot r} \end{matrix}

The photon current can be obtained from the continuity equation for photon density:

\partial_{t} n^{μ} = - \nabla \cdot F^{μ} .

Which leads to the following expression for the photon current per unit of area per unit time:

F_{x} = \sum_{μ} F^{μ} \cdot \hat{x} = \frac{c}{V} \sum_{μ, k, k_{x} > 0} (\hat{k} \cdot x) a_{k}^{† μ} a_{k}^{μ}

The expression above would give us the number of photons moving on a single direction (e.g. towards the rigth) on a given surface of a black-body, which can be pictured to be the surface of the Sun. On the surface of the Earth the photon flux is given by:

\begin{aligned} F_{s} = {(\frac{R_{S}}{R_{E}})}^{2} \frac{c}{V} \sum_{μ, k, k_{x} > 0} & (\hat{k} \cdot x) ⟨ a_{k}^{† μ} a_{k}^{μ} ⟩ = \\ = {(\frac{R_{S}}{R_{E}})}^{2} \frac{2 ζ (3)}{(2 π)^{2}} \frac{c}{(β ℏ c)^{3}} \end{aligned}

here

ζ

is the Riemann zeta-function. Therefore, we define a three-component vectorial figure of merit called the Solar Rectification Vector (SRV), and denoted by

η^{γ} (γ = x, y, z)

, as the ratio of the electric current density and the photon flux, given by:

η^{γ} \equiv \frac{\bar{j^{γ}}}{e F_{s}}

Notice that the SRV defined above depends only on the temperature of the sun and the intrinsic rectification properties of the material. By combining Eqs. (18) and Eq.(26) with the definition from Eq.(27), one arrives at the following explicit form of the SRV:

η^{γ} = \frac{ℏ}{2 ζ (3) e c ε_{0}} \int_{- \infty}^{\infty} d ω σ_{(2)}^{γ α α} (- ω, ω) I (ℏ β_{S} | ω |)

Where

I (x)

is defined in Eq. (21). We see that

η^{γ}

is a dimensionless vector for 3D materials and its magnitude,

η = \sqrt{\sum_{i} η_{i}^{2}}

, counts the number of electrons contributing to the rectified current per incoming photon, therefore serving as a dimensionless figure of merit for the ideal bulk photo-current generated by a material under shortcircuit configurations. Because the Planck distribution has vanishing intensity at small frequency, in the clean limit, the BCD Drude-like peak does not contribute to the rectified current and all the contributions arise from shift currents in time-reversal invariant materials, while injection current contributions vanish (see Appendix B).

It is instructive to contrast this figure of merit with a comparable figure of merit for current generation for the ideal model of a p-n junction solar-cell proposed by Shockley-Quessier (SQ) [29]. In the most ideal version of the SQ model one assumes that every photon with energy above the band-gap is absorbed and generates an electron and a hole that are successfully collected at opposite electrodes contributing to the current generation of the material. Therefore in the SQ model the total current density is:

j_{S Q} = e F_{s} η_{S Q}

Here

η_{S Q}

is the analogue of the magnitude of the SRV, but defined for a conventional solar cell within the SQ model, and simply measures fraction of incident photons that have energy above the band-gap and thus can be absorbed by the material. The flux of such photons with energies above the band-gap,

Δ_{0}

, can be computed in an analogous fashion to the total photon flux but inserting an extra factor to enforce that their energy is larger than

Δ_{0}

, as follows:

F (Δ_{0}) = {(\frac{R_{S}}{R_{E}})}^{2} \frac{c}{V} \sum_{\begin{matrix} μ, k \\ k x > 0 \end{matrix}} Θ (c k - Δ_{0}) (\hat{k} \cdot x) ⟨ a_{k}^{† μ} a_{k}^{μ} ⟩

where

Θ

is the Heaviside function. The ratio of absorbed photons,

η_{S Q}

, in Eq. (29), can thus be obtained as

η_{S Q} = F (Δ_{0}) / F_{s}

. SQ were primarily interested in optimizing the total power delivered by the solar cell and

FIG. 4: Illustration of the geometry for incident (I), transmited (T) and reflected (R) light at the solar cell. For simplicity, we considered the light to be perpendicular to the surface. At the interface (solid black line) between vacuum (region 1) and the bulk photovoltaic material that makes the solar cell (region 2) the incoming light (I) splits onto reflected (R) and transmitted (T). The transmitted component generates the current J.

μ_{0}

is magnetic permitivity and

ϵ_{1}, ϵ_{2}^{x, y}

are permittivities of vacuum and the bulk photo-voltaic material respectively.

not the total current. It is however clear that the maximum current density is attained for

η_{S Q} = 1

which is achieved in the limit of vanishing gap. As we will see LiAsSe

_{2}

has a value of

η = 0.05

, and thus can produce a photocurrent that is

5 %

of the ideal

S Q

limit, in which every incoming photon produces an electron and a hole that do not recombine and are successfully captured at opposite electrodes contributing to the net current.

7. MATERIAL PENETRATION

Our discussion of photo-current generation so far has assumed that light propagates inside the material as it would do in vacuum. This is a reasonable approximation only if the incoming radiation is not substantially modified inside the material. This can be justified for example if the material is much thinner than the decay length for light absorbtion in the material,

L \sim λ_{i n} / Im [n]

, where

λ_{i n}

is the typical wavelength of the incoming radiation and

n

is the complex index of refraction of the material. However, a more realistic calculation of the photocurrent would include the frequency dependence of light absorbtion and propagation in the material, which is often accounted for by the Glass coefficient [11, 34].

In this section we take into account that the intensity of electric field decays inside the material, and derive the expression for the total current. For simplicity, we assume that incident radiation is normal to the solar cell surface so that that the electric field is parallel to the surface, namely along the XY plane in Fig.[4], and that the induced photo-current also flows along this XY plane. To estimate the electric field penetration, we use linear response theory. The geometry of the problem is shown in Fig.[4], and by the standard procedure of boundary condition matching from Maxwell equations (see e.g. Ref.[35]) leads to the following equations:

\begin{matrix} E_{I} (z, t) = E_{I} e^{i (k_{1} z - ω t)}, E_{R} (z, t) = E_{R} e^{- i (k_{1} z + ω t)} \\ H_{I} (z, t) = H_{I} e^{i (k_{1} z - ω t)}, H_{R} (z, t) = H_{R} e^{- i (k_{1} z + ω t)} \\ E_{I, R} = (E_{I, R}^{x}, E_{I, R}^{y}, 0), H_{I, R} = (H_{I, R}^{x}, H_{I, R}^{y}, 0) \\ E_{T} (z, t) = (E_{T_{x}} e^{i (k_{2 x} z - ω t)}, E_{T_{y}} e^{i (k_{2 y} z - ω t)}, 0) \\ \nabla \times H = - i J - ε_{0} ω E, \nabla \times E = ω μ_{0} H \\ J^{x} (ω) = σ_{(1)}^{x x} (ω) E^{x} (ω), J^{y} (ω) = σ_{(1)}^{y y} (ω) E^{y} (ω) \end{matrix}

For simplicity we assume the magnetic permitivities of the vacuum and the bulk photovoltaic material to be frequency independent constants. We also assume that the material has a point group so that its linear conductivity is diagonal in its set of principal axes, and we take one these axes to be the z-axis that coincides with the direction of incoming light. The frequency dependence of these diagonal components of the linear conductivity of the material can be obtained from standard linear response theory Refs.

[10, 32]

and are given by:

\begin{aligned} σ_{(1)}^{β β} (ω) & = - \frac{e^{2}}{ℏ} \sum_{n m} \int \frac{d^{3} k}{(2 π)^{3}} {δ_{n m} \frac{\partial ϵ_{n}}{\partial k^{β}} \frac{i \frac{\partial f_{n}}{\partial k^{β}}}{ω + i Γ} + \\ + i (f_{n} - f_{m}) (ϵ_{n} - ϵ_{m}) \frac{A_{n m}^{β} A_{m n}^{β}}{ω - ϵ_{n} + ϵ_{m} + i Γ}} \end{aligned}

where

Γ

is relaxation rate and

ϵ_{n}

is the n-th Bloch band energy dispersion. In the expression above

β

index is not summed over. From the above one can obtain the dispersion relations of light outside and inside the material to be:

\begin{aligned} k_{1}^{2} = μ_{0} ω^{2} ε_{0}, k_{2 x}^{2} = μ_{0} ω^{2} ε_{2}^{x} (ω), k_{2 y}^{2} = μ_{0} ω^{2} ε_{2}^{y} (ω), \\ ε_{2}^{x} (ω) = ε_{0} + i \frac{σ_{(1)}^{x x} (ω)}{ω}, ε_{2}^{y} (ω) = ε_{0} + i \frac{σ_{(1)}^{y y} (ω)}{ω} \end{aligned}

By taking into account standard electromagnetic boundary conditions [35], we can obtain the following relation describing the transmitted,

T

, component of the electric field inside the material for a given incident component, I:

\begin{aligned} E_{T} (ω, z) = \\ = \frac{2 n_{1}}{n_{1} + n_{2}^{x} (ω)} (E_{I}^{x} (ω), 0, 0) e^{i (Re [k_{2 x}] z - ω t)} e^{- Im [k_{2 x}] z} + \\ + \frac{2 n_{1}}{n_{1} + n_{2}^{y} (ω)} (0, E_{I}^{y} (ω), 0) e^{i (Re [k_{2 y}] z - ω t)} e^{- Im [k_{2 y}] z} \\ n_{1} = c \sqrt{μ_{0} ε_{0}}, n_{2}^{β} = c \sqrt{μ_{0} ε_{2}^{β} (ω)}, k_{2 β} = ω n_{2}^{β} (ω) / c \end{aligned}

Since the above relation between transmitted and incoming light is linear, we can compute ensemble averages of the transmitted electric field from the known ensemble averages of the incoming electric field. We take the ensemble average of the incoming electric field to be that of the ideal black-body radiation from Eq.(18), where we keep track explicitly of the transverse character of incoming light. From this it follows that the electric field auto-correlation function inside the bulk of the solar cell is:

\begin{aligned} \bar{E_{T}^{α} (ω, z) E_{T}^{β} (- ω, z)} = \\ = & δ_{α β} {| \frac{2 n_{1}}{n_{1} + n_{2}^{β} (ω)} |}^{2} e^{- 2 \frac{ω}{c} Im [n_{2}^{β} (ω)] z} \bar{E_{I}^{β} (ω) E_{I}^{β} (- ω)}, \end{aligned}

where in the above expression the repeated indices are not being summed over and indices are understood to take only

(x, y)

components. Thus the total current is:

\begin{aligned} \frac{\bar{I^{γ}}}{W} = \sum_{β = x, y} \int_{- \infty}^{\infty} d ω {| T^{β} (ω) |}^{2} \frac{σ_{(2)}^{γ β β} (- ω, ω)}{α^{β} (ω)} I^{β β} (ω) \\ T^{β} (ω) = \frac{2 n_{1}}{n_{1} + n_{2}^{β} (ω)}, α^{β} (ω) = 2 \frac{ω}{c} Im [n_{2}^{β} (ω)] \end{aligned}

where

W

is the width of the sample,

n_{2}^{β} (ω)

is the complex index of refraction from Eq. (41), and

I^{β β} (ω)

is given in Eq.

(20)

. The formula above generalizes that from Refs.

[11, 34]

by adding the frequency dependence of the reflection coefficient at the material surface, the spectral distribution of incident light as well as the directional dependence of the conductivity of the material. We expect it to provide accurate estimates for realistic bulkphotovoltaic solar cell devices.

8. BULK RECTIFICATION OF TaAs AND ${L i A s S e}_{2}$

In this section we combine our general considerations with first principle calculations of the non-linear optoelectronic properties of the 3D Weyl semi-metal TaAs and the insulator

{L i A s S e}_{2}

. Because it is semi-metallic, TaAs has BCD and shift current contributions whereas

{L i A s S e}_{2}

only has shift current contributions. Our primary interest in TaAs is to investigate the interplay between the intraband BCD and interband shift-currents to its BPRT, although we also compute its SRV, while in the case of

{L i A s S e}_{2}

our primary interest is the estimation of its SRV due to its promise as a photovoltaic material, but we also compute its BPRT. Throughout this section we will be using the ideal formalism that neglects the modification of the incoming radiation inside the material, namely we will compute the BPRT and SRV from Eqs. (6) and (27), respectively.

The band structures of

{L i A s S e}_{2}

and TaAs were calculated by the full-potential local-orbital (FPLO)

[36, 37]

minimum-basis code with the experimental crystal structure. Details of the calculations of band structure, Berry phases and non-linear rectification conductivity in the clean limit are presented in Appendix C. The frequency dependent non-linear conductivity of

{L i A s S e}_{2}

is shown in Fig.[3] (for the one of TaAs see Appendix C). Our results are in agreement with previous reports

[6, 38]

. In particular, the

σ^{z x x}

and

σ^{z z z}

components for TaAs can exceed

300 μ A / V^{2}

at infrared frequencies, while the

σ^{z z z}

component of

{L i A s S e}_{2}

can exceed

100 μ A / V^{2}

in a wide range frequency from

\sim 1.25

\sim 1.65 e V

, which as we will see, leads to a sizable SRV. The space symmetry I4

4_{1}

md of TaAs leads to the following non-zero components of the BPRT:

\begin{matrix} R^{x x z} = R^{y y z} = R^{x z x} = R^{y z y} = 0.13 \\ R^{z x x} = R^{z y y} = 0.05, R^{z z z} = - 0.04 . \end{matrix}

The values above contain the contribution from BCD and shift currents. The BCD part of the contribution to BPRT is:

\begin{matrix} R_{B C D}^{x x z} = R_{B C D}^{y y z} = R_{B C D}^{x z x} = R_{B C D}^{y z y} = 0.09 \\ R_{B C D}^{z x x} = R_{B C D}^{z y y} = - 0.18 \end{matrix}

Therefore, remarkably, we see that in spite of having all its spectral weight localized in a delta function Drude-like peak, in TaAs the BCD contributes to the net rectification tensor with a magnitude that is comparable with the net contribution from shift currents, and sometimes they have oppposite signs as it is the case for the

R^{z x x} = R^{z y y}

. The SRV for TaAs has a single non-trivial component

η_{T a A s} = (0, 0, - 0.001)

, which is oriented along its polar axes as depicted in Fig.2.

On the other hand, the lower Cc space group symmetry of

{L i A s S e}_{2}

allows for a larger set of independent components:

\begin{array}{ccc} R^{x x x} = 0.06, & R^{x y y} = - 0.05, & R^{x z z} = - 0.09, \\ R^{z x x} = 0.01, R^{z y y} = - 0.05, & R^{z z z} = 0.18 \\ R^{x x z} = R^{x z x} = - 0.002, & R^{y y x} = R^{y x y} = - 0.04 \\ R^{y y z} = R^{y z y} = - 0.02, & R^{z z x} = R^{z x z} = - 0.8 \end{array}

The SRV has two components

η_{{L i A s S e}_{2}} = (- 0.02, 0, 0.04)

, with the higher SRV magnitude

η = 0.05

, or in other words, it produces a photocurrent that is is

5 %

of the ideal limit for p-n junction solar cells.

9. DISCUSSION

We have introduced an operational definition of the Berry Phase Rectification Tensor as the second order polarization response of a material to an ideal delta function pulse of applied electric field. In the case of a time reversal invariant band structure this tensor is an intrinsic property of the material which only depends on the Berry connections of the Bloch bands and not on their energies. This tensor contains intra-band contributions which arise from the Berry curvature dipole and inter-band contributions which arise from the shift currents. This operational definition of the tensor is in principle applicable in the presence of electron-electron interactions, which is an interesting direction for future research. We have computed these tensors for TaAs, where we find that remarkably the Berry Dipole contributions are comparable to the shift current contributions.

We have also introduced a technologically relevant figure of merit for bulk photo-current generation termed the Solar Rectification Vector, which counts the number of electrons contributing to the rectified current per incoming photon, under ideal solar black-body radiation, in analogy to the model of solar cells in the classic work of Shockley and Queisser [29]. Using this, we have found that

{L i A s S e}_{2}

can produce a large bulk current of about

0.05

electrons per incident photon. We have also generalized the formula of photo-current production defining the Glass coefficient

[11, 34]

to include the black-body spectral distribution of intensity, the reflection coefficient at the surface of the material and the directional dependence of the linear conductivity of the material, which should provide accurate estimates of realistic photo-current production in solar cells based on bulk rectification mechanisms.

We would like to mention in passing that the formulae for rectification conductivities of Refs. [1, 2] are incomplete and are missing terms even for the inter-band shift currents. Details of the discrepancy and comparison with other various formulae employed in the literature

[1, 2, 4, 5, 9, 32, 39 - 43]

are presented in (see Appendix D). For the complete correct formulae for the general second order conductivity of any band structure see Ref. [10].

10. ACKNOWLEDGEMENTS

We thank Qian Niu for pointing to us the operational interpretation of the QRSR in terms of electric field pulses. We also thank Qian Niu, Cong Xiao, Jorge Facio, Jeroen van-den-Brink, Dennis Wawrzik, and Jhih-Shih You for valuable discussions.

11. APPENDIX A: SECOND ORDER RESPONSE TO ELECTRIC FIELDS IN TIME AND FREQUENCY DOMAINS

In time domain, the current that is second order in electric fields can be written as:

j_{2}^{γ} (t) = \iint_{- \infty}^{\infty} d t_{1} d t_{2} E^{β} (t_{2}) E^{α} (t_{1}) σ_{(2)}^{γ β α} (t, t_{1}, t_{2})

The second order conductivity appearing above was derived in Ref.[10], and is given by:

σ_{(2)}^{γ β α} (t, t_{1}, t_{2}) = - \frac{e^{3}}{ℏ^{2}} e^{Γ (t_{2} - t)} ⟨ [r^{β} (t_{2}), [r^{α} (t_{1}), v^{γ} (t)]] ⟩

for

t \geq t_{1} \geq t_{2}

and 0 otherwise. This correlation function has global time translation invariance, namely:

σ_{(2)}^{γ β α} (t, t_{1}, t_{2}) = σ_{(2)}^{γ β α} (t - τ, t_{1} - τ, t_{2} - τ),

which implies it has only two independent arguments when we shift the arguments by

t (τ = t)

. Thus using the following convention for Fourier transforms:

\begin{aligned} f (t) = \int_{- \infty}^{\infty} d ω f (ω) e^{- i ω t} \\ f (ω) = \int_{- \infty}^{\infty} \frac{d t}{2 π} f (t) e^{i ω t} \end{aligned}

we introduce the two argument conductivity in a frequency domain as:

\begin{aligned} σ_{(2)}^{γ β α} (0, τ_{1}, τ_{2}) = \\ = \iint \frac{d ω_{1} d ω_{2}}{(2 π)^{2}} e^{- i ω_{1} τ_{1} - i ω_{2} τ_{2}} σ_{(2)}^{γ β α} (ω_{1}, ω_{2}) . \end{aligned}

The Eq.(53) is the definition of the

σ_{(2)}^{γ β α} (ω_{1}, ω_{2})

. Now, the expression above leads to the following expression for the current:

\begin{aligned} j_{2}^{γ} (ω) = \int_{- \infty}^{\infty} d ω_{1} E^{β} (ω - ω_{1}) & E^{α} (ω_{1}) \times \\ \times σ_{(2)}^{γ β α} (- ω_{1}, ω_{1} - ω) \end{aligned}

Therefore the net rectified or DC current can be obtained from the above as the current at frequency

ω = 0

, and is given by:

j_{2, D C}^{γ} = \int_{- \infty}^{\infty} d ω E^{β} (- ω) E^{α} (ω) σ_{(2)}^{γ β α} (- ω, ω)

By using the above definitions and conventions one can arrive at the Eq.(4) from the main text, namely, we have that:

\begin{array}{r} \int_{- \infty}^{\infty} d t σ_{(2)}^{γ β α} (t, 0, 0) = \int_{- \infty}^{\infty} d t σ_{(2)}^{γ β α} (0, - t, - t) = \\ = ∭_{- \infty}^{\infty} d t \frac{d ω_{1} d ω_{2}}{(2 π)^{2}} σ_{(2)}^{γ β α} (ω_{1}, ω_{2}) e^{i (ω_{1} + ω_{2}) t} = \\ = \int_{- \infty}^{\infty} \frac{d ω}{2 π} σ_{(2)}^{γ α β} (- ω, ω) \end{array}

Now, the total current of the system is given by:

\begin{aligned} j^{γ} (t) = \int_{- \infty}^{\infty} d t_{1} σ_{(1)}^{γ α} (t, t_{1}) E^{α} (t_{1}) + \\ + \int_{- \infty}^{\infty} d t_{1} \int_{- \infty}^{\infty} d t_{2} σ_{(2)}^{γ α β} (t, t_{1}, t_{2}) E^{β} (t_{1}) E^{α} (t_{2}) \end{aligned}

after averaging over ensemble, namely using relations:

\bar{E^{α} (t_{1})} = 0, \bar{E^{β} (t_{1}) E^{α} (t_{2})} = I^{β α} (t_{1} - t_{2}),

we obtain the following expressions:

\begin{aligned} \bar{j^{γ} (t)} = \iint_{- \infty}^{\infty} d t_{1} d t_{2} σ_{(2)}^{γ α β} (0, t_{1} - t, t_{2} - t) I^{β α} (t_{1} - t_{2}) = \\ = \int_{- \infty}^{\infty} d ω σ_{(2)}^{γ α β} (- ω, ω) I^{β α} (ω) \end{aligned}

12. APPENDIX B: SUM RULE COMPUTATIONS

The second order rectification conductivity in general is given by:

\begin{aligned} σ_{(2)}^{γ β α} (- ω, ω) = σ_{J}^{γ β α} & (- ω, ω) + σ_{B C D}^{γ β α} (- ω, ω) + \\ + σ_{I}^{γ β α} (- ω, ω) + σ_{S}^{γ β α} (- ω, ω) \end{aligned}

where J stands for "Jerk", BCD for "Berry curvature dipole", I for "Injection" and

S

for "shift current" respectively. We consider band structure with small relaxations

(\forall n, m, n \neq m : Γ ≪ ϵ_{n} - ϵ_{m})

. Each contribution is given by:

σ_{J}^{γ β α} (- ω, ω) = \frac{e^{3}}{ℏ^{2}} \int \frac{d k}{(2 π)^{3}} \sum_{n m} \frac{\frac{\partial ϵ_{n}}{\partial k^{γ}} \frac{\partial^{2}}{\partial k^{α} \partial k^{β}} f_{n} δ_{n m}}{ω^{2} + Γ^{2}}

\begin{aligned} σ_{B C D}^{γ β α} (- ω, ω) = - \frac{1}{2} \frac{e^{3}}{ℏ^{2}} \frac{1}{ω + i Γ} \int \frac{d k}{(2 π)^{3}} \times \\ \times \sum_{n m} {\hat{A}}_{m n}^{γ} {\hat{A}}_{n m}^{α} \frac{\partial}{\partial k^{β}} (f_{m} - f_{n}) + (\begin{array}{c} α \leftrightarrow β \\ ω \leftrightarrow - ω \end{array}) \end{aligned}

σ_{I}^{γ β α} (- ω, ω) = \frac{e^{3}}{ℏ^{2}} \int \frac{d k}{(2 π)^{3}} \times

\begin{array}{r} \times \sum_{n m} \frac{(f_{m} - f_{n}) {\hat{A}}_{n m}^{β} {\hat{A}}_{m n}^{α} (\frac{\partial}{\partial k_{γ}^{γ}} ϵ_{n} - \frac{\partial}{\partial k^{γ}} ϵ_{m})}{{(ω - ϵ_{n} + ϵ_{m})}^{2} + Γ^{2}} + \\ + (\begin{array}{c} α \leftrightarrow β \\ ω \leftrightarrow - ω \end{array}) \end{array}

\begin{aligned} σ_{S}^{γ β α} (- ω, ω) = \frac{1}{2} \frac{e^{3}}{ℏ^{2}} \int \frac{d k}{(2 π)^{3}} \times \\ \times \sum_{n m} {{\hat{A}}_{m n}^{γ} \frac{\partial}{\partial k^{α}} \frac{(f_{n} - f_{m}) {\hat{A}}_{n m}^{β}}{ω - ϵ_{n} + ϵ_{m} + i Γ} + \\ + i \frac{(f_{n} - f_{m}) {\hat{A}}_{n m}^{β}}{ω - ϵ_{n} + ϵ_{m} + i Γ} \sum_{c} [{\hat{A}}_{m c}^{α} {\hat{\bar{A}}}_{c n}^{γ} - {\hat{\bar{A}}}_{m c}^{γ} {\hat{A}}_{c n}^{α}]} + \\ + (\begin{array}{c} α \leftrightarrow β \\ ω \leftrightarrow - ω \end{array}) . \end{aligned}

For purposes of ensemble averaging over incoming random light, we need to consider only diagonal components of the nonlinear rectification conductivity, namely

σ_{(2)}^{γ α α}

To simplify the notation for further analysis we write:

Q_{n m}^{γ β α} = {\hat{A}}_{n m}^{β} \frac{\partial}{\partial k^{α}} {\hat{A}}_{m n}^{γ} - i {\hat{A}}_{n m}^{β} \sum_{c} [{\hat{A}}_{m c}^{α} {\hat{A}}_{c n}^{γ} - {\hat{A}}_{m c}^{γ} {\hat{A}}_{c n}^{α}]

Thus, the conductivity which is the one that enters into the SRV, is:

σ_{J}^{γ α α} (- ω, ω) = \frac{e^{3}}{ℏ^{2}} \int \frac{d k}{(2 π)^{3}} \sum_{n m} \frac{\frac{\partial ϵ_{n}}{\partial k^{γ}} \frac{\partial^{2}}{\partial k^{α} \partial k^{α}} f_{n} δ_{n m}}{ω^{2} + Γ^{2}}

σ_{B C D}^{γ α α} (- ω, ω) = \frac{1}{2} \frac{e^{3}}{ℏ^{2}} \frac{Γ}{ω^{2} + Γ^{2}} \int \frac{d k}{(2 π)^{3}} \times

\times \sum_{m} f_{m} \frac{\partial}{\partial k^{β}} Ω_{m}^{α γ}

\begin{aligned} σ_{I}^{γ α α} ( & - ω, ω) = \frac{e^{3}}{h^{2}} \int \frac{d k}{(2 π)^{3}} \times \\ \times \sum_{n m} \frac{(f_{m} - f_{n}) {\hat{A}}_{n m}^{α} {\hat{A}}_{m n}^{α} (\frac{\partial}{\partial k^{γ}} ϵ_{n} - \frac{\partial}{\partial k^{k}} ϵ_{m})}{{(ω - ϵ_{n} + ϵ_{m})}^{2} + Γ^{2}} \end{aligned}

\begin{aligned} σ_{S}^{γ α α} (- ω, ω) = \frac{e^{3}}{ℏ^{2}} \int \frac{d k}{(2 π)^{3}} \sum_{n m} (f_{m} - f_{n}) \times \\ \times \frac{Re [Q_{n m}^{γ α α}] (ω - ϵ_{n} + ϵ_{m}) + Γ Im [Q_{n m}^{γ α α}]}{{(ω - ϵ_{n} + ϵ_{m})}^{2} + Γ^{2}} \end{aligned}

Time reversal symmetry (TRS) implies:

\begin{matrix} T A_{m n} (k) T^{- 1} = A_{n m} (- k), \\ T r_{m n} (k) T^{- 1} = r_{n m} (- k) \\ T v_{m n} (k) T^{- 1} = - v_{n m} (- k), \end{matrix}

TRS makes the Jerk, Injection and the part proportional to

Re [Q^{γ α α}]

of the Shift current vanish after momentum integration. Thus the contribution to the total rectification conductivity is given by BCD and resonant part of the shift current:

\begin{aligned} σ_{B C D}^{γ α α} (- ω, ω) & + σ_{S}^{γ α α} (- ω, ω) = \frac{e^{3}}{ℏ^{2}} \int \frac{d k}{(2 π)^{3}} \times \\ {\frac{Γ}{ω^{2} + Γ^{2}} \sum_{m} f_{m} \frac{\partial}{\partial k^{α}} Ω_{m}^{α γ} + \\ + \sum_{n m} (f_{m} - f_{n}) \frac{Γ Im [Q_{n m}^{γ α α}]}{{(ω - ϵ_{n} + ϵ_{m})}^{2} + Γ^{2}}} \end{aligned}

Γ \to 0

limit, we can simplify the expression using the following relation:

lim_{a \to 0^{+}} \frac{a}{d_{z}^{2} + a^{2}} = π δ (d_{z})

Now we can write a general version of the sum rule by including the spectral weight, which we assume to be a smooth function of

ω

, to obtain:

\begin{aligned} \int d ω F (ω) σ^{γ α α} (- ω, ω) = π \frac{e^{3}}{ℏ^{2}} \sum_{α, n m} f_{n m} F (| ϵ_{n m} |) \times \\ \times {A_{n m}^{α} i \partial^{α} A_{m n}^{γ} + \sum_{l} A_{n m}^{α} [A_{m l}^{α} A_{l n}^{μ} - A_{m l}^{γ} A_{l n}^{α}]} \end{aligned}

where

f_{n m} = f_{n} - f_{m}

and all zero-frequency contributions vanish under the assumption of absence of spectral weight

(F (ω \to 0) = 0)

Now in the case of the operational definition of the BPRT, we consider the net change of the polarization coursed by pulse-like Electric field:

\begin{matrix} E (t) = E_{0} τ_{0} δ (t), E (ω) = \frac{E_{0} τ_{0}}{2 π} \\ Δ P = \int_{0}^{\infty} j (t) d t \end{matrix}

In this set-up, without assuming time reversal symmetry the resulting change of the polarization is:

\begin{aligned} Δ P^{γ} = & \frac{E_{0}^{α} E_{0}^{β}}{4} τ_{0}^{2} \frac{e^{3}}{ℏ^{2}} \times \\ \times {\frac{1}{Γ} (⟨ \partial^{α} \partial^{β} v^{γ} ⟩ + \sum_{a b} f_{a b} A_{a b}^{β} A_{b a}^{α} Δ_{b a}^{γ}) + \\ + ⟨ [A^{β}, i \partial^{α} A^{γ}] ⟩ + ⟨ \partial^{β} Ω^{α μ} ⟩ + \\ + ⟨ [A^{β}, [A^{α}, {\bar{A}}^{γ}]] ⟩ + (α \leftrightarrow β)} \end{aligned}

where

\sum_{a b} = \sum_{a b} \int d k / (2 π)^{d}

. As we see this result contains two extra terms proportional to

1 / Γ

which arise from the Jerk and injection current terms, but which vanish under the assumption of time-reversal symmetry.

13. APPENDIX C: COMPUTATION DETAILS

The exchange and correlation energies were considered in the generalized gradient approximation with PerdewBurke-Ernzerhof functional [36]. We projected the Bloch wave functions into high-symmetry atomic-orbital-like Wannier functions with diagonal position operator and considered all the electrons, which guarantees the full band overlap from ab-initio and Wannier functions in the energy window from

- 20

20 e V

. The calculated band structures of TaAs and LiAsSe2 are shown in Fig[5]. Based on the highly symmetric Wannier functions, we constructed tight-binding (TB) model Hamiltonians and calculated the Berry curvature dipole and second order optical conductivity tensor. The

γ

(with

α, β, γ = x, y, z

) component of Berry curvature for the

n -

th band at

k

FIG. 5: Band structures of (a) TaAs and (b)

{L i A s S e}_{2}

. Rectification conductivities of (c) TaAs and (d) LiAsSe .

point is computed by following:

Ω_{n}^{γ} (k) = \frac{ε^{γ α β}}{2} (\frac{\partial {\hat{A}}_{n}^{β}}{\partial k^{α}} - \frac{\partial {\hat{A}}_{n}^{α}}{\partial k^{β}})

where,

{\hat{A}}_{n}^{β} = ⟨ n, k | i \partial^{β} | n, k ⟩

is the Berry connection of

n

th band of Hamiltonian

H (k)

and

ε^{γ α β}

is the Levi-Civita tensor. The Berry curvature dipole was calculated from the Berry curvature in a differential manner and with integral in the whole k-space [17]:

D_{α β} = \int \frac{d k}{(2 π)^{3}} \sum_{n} f_{n} \frac{\partial Ω^{β} (k)}{\partial k^{α}} = ⟨ \frac{\partial Ω^{β} (k)}{\partial k^{α}})

where

Ω_{n}^{δ} = ε^{δ α β} \partial A^{β} / \partial k^{α}

is the Berry curvature vector and

f_{n}

is a Fermi-Dirac distribution. The relation between two and three index BCD tensors is:

ε^{α γ δ} D^{β δ} = \partial^{β} Ω^{α γ}

The second order optical conductivity tensor is calcu- lated by the following expression:

\begin{matrix} σ_{(2)}^{γ α β} (- ω, ω) = \frac{π e^{3}}{2 ℏ^{2}} \int \frac{d k}{(2 π)^{3}} \sum_{n m} (f_{n} - f_{m}) \times \\ \times {A_{n m}^{β} i \partial^{α} A_{m n}^{γ} + A_{n m}^{β} {[A^{α}, {\bar{A}}^{γ}]}_{m n} + \\ + (α \leftrightarrow β)} δ (ω - \frac{ϵ_{n} - ϵ_{m}}{ℏ}) \end{matrix}

here

ϵ_{n}

is n-th band energy. The expression above is equivalent to expressions discussed in

[5, 10]

(see Appendix D). The numerical results are consistent with the symmetry analysis. The plots showing the frequency dependence of the second order rectification conductivities of TaAs and

{L i A s S e}_{2}

are shown in Fig.[5]. These plots show only the interband shift currents, since the BCD is simply a delta function peak at zero frequency.

After a k-grid convergency test, it was found that the change of BCD for TaAs is less than

5 %

with k-grid increasing from

360 \times 360 \times 360

480 \times 480 \times 480

. A similar

k

-grid convergency was obtained with a

k

-grid of increasing from

240 \times 240 \times 240

360 \times 360 \times 360

for the second order optical conductivity in TaAs and

{L i A s S e}_{2}

14. APPENDIX D: EQUIVALENCE WITH OTHER FORMALISMS IN THE LITERATURE

Here we discuss the equivalence of our formulae for the rectification conductivity from Eq.

(60)

and others employed in the literature

[1, 2, 4, 5, 9, 32, 39, 40]

. Importantly, we show that shift currents derived in classic papers Refs.

[1, 2]

are only valid for Galilean systems, where spin orbit coupling is neglected, but not in general.

15. A. Comparison with Refs. $[5, 39]$

Here we will demonstrate that the interband rectification current in time reversal invariant materials arising from

σ_{I} + σ_{S}

from

Eqs (61 - 64)

, is identical to that from Refs.

[5, 39]

. This is the part of the rectification conductivity that controls the inter-band contribution to the QRSR and BPRT, but it is missing the intraband BCD contribution. This interband part however is the only one contributing to the SRV.

We define the following notation:

\begin{aligned} r_{n m}^{a} = \frac{v_{n m}^{a}}{i ϵ_{n m}} = {\bar{A}}_{n m}^{a}, (n \neq m) \\ r_{n m}^{a; b} = \partial^{b} r_{n m}^{a} - i r_{n m}^{a} (A_{n n}^{b} - A_{m m}^{b}) \end{aligned}

where

ϵ_{n m} = ϵ_{n} - ϵ_{m}

and

(a, b, c)

are space indices.

A l -

ternatively, this can be written as:

\begin{array}{r} r_{n m}^{a; b} = \frac{i}{ϵ_{n m}} {\frac{v_{n m}^{a} Δ_{n m}^{b}}{ϵ_{n m}} + \frac{v_{n m}^{b} Δ_{n m}^{a}}{ϵ_{n m}} - ϑ_{n m}^{a b} + \\ + \sum_{p \neq n, m} (\frac{v_{n p}^{a} v_{p m}^{b}}{ϵ_{p m}} - \frac{v_{p m}^{a} v_{n p}^{b}}{ϵ_{n p}})} \end{array}

where

Δ_{n m}^{a} = v_{n n}^{a} - v_{m m}^{a}

, and

\begin{aligned} ϑ_{n m}^{a b} = ⟨ n | \partial^{a} \partial^{b} & H | m ⟩ = \\ = \partial^{b} v_{n m}^{a} + i \sum_{p} (v_{n p}^{a} r_{p m}^{b} - r_{n p}^{b} v_{p m}^{a}) \end{aligned}

By using the identity identity

[r^{a}, r^{b}] = 0

, one obtains:

\begin{aligned} \partial^{a} A_{n m}^{b} - i {[A^{a}, {\bar{A}}^{b}]}_{n m} & = \\ = \partial^{b} A_{n m}^{a} - i A_{n m}^{a} (A_{n n}^{b} - A_{m m}^{b}) \end{aligned}

which leads to an alternative form to denote the generalized derivative:

\begin{aligned} r_{n m}^{a; b} = \partial^{b} A_{n m}^{a} - i A_{n m}^{a} (A_{n n}^{b} - A_{m m}^{b}) = \\ = \partial^{a} A_{n m}^{b} - i \sum_{p} {A_{n p}^{c} {\bar{A}}_{p m}^{a} - {\bar{A}}_{n p}^{a} A_{p m}^{c}} \end{aligned}

for

n \neq m

. By using Eq. (83) and Eq.(88) one can obtain that the conductivity from

[5, 39]

can be written as:

\begin{aligned} σ^{a b c} (0, ω, - ω) & = - \frac{i π e^{3}}{2 ℏ^{2}} \int \frac{d k}{(2 π)^{3}} \sum_{n m} f_{n m} \times \\ \times (r_{m n}^{b} r_{n m}^{c; a} + (b \leftrightarrow c)) δ (\frac{ϵ_{m n}}{ℏ} - ω) \end{aligned}

this is the same as the sum of the real part of injection and shift currents

σ_{S} + σ_{I}

assuming time-reversal invariant conditions from Eq.

(60)

, discussed in section (see Appendix B), which were the ones employed in the main text. In the main text, we are also omitting the " 0 " in the conductivity arguments, but it is implicit every time we refer to a "rectification" conductivity.

16. B. Comparison with Refs. $[1, 2, 9]$

We will here compare our formula from the current paper and Ref.[10] to the formulae from the classic Ref. [1, 2] (see eg. Ref[38, 40] for an application of this) and those of more recent work of Ref.[9] (see also Ref.[41-43] for closely related formulae to Ref.[9]). We will show that Refs.

[1, 2]

are missing crucial corrections that can only be neglected if one assumes the system to have an underlying parabolic dispersion (namely ignores spin-orbit coupling). On the other hand we will show that, upon proper regularization of infinitesimal imaginary parts in the frequency denominators, the formulae of Ref.

[9]

are equivalent to those we employ in the current paper and in our previous work Ref.[10].

Following Ref.[9] we write the second order perturbation to the Hamiltonian, as:

\hat{H} = {\hat{H}}_{0} + \hat{V} (t),

\begin{aligned} \hat{V} (t) = \frac{e}{ℏ} A^{α} {\hat{h}}^{α} + \frac{1}{2} {(\frac{e}{ℏ})}^{2} A^{α} (t) A^{β} (t) {\hat{h}}^{α β} + \dots, \end{aligned}

where

\hat{H}

is total Hamiltonian,

{\hat{H}}_{0}

is unperturbed band Hamiltonian,

E (ω) = i ω A (ω)

is electric field and

h_{a b}^{μ} = v_{a b}^{μ} = \partial^{μ} ϵ_{a} δ_{a b} + i ϵ_{a b} A_{a b}^{μ},

\begin{aligned} h_{a b}^{μ α} = \partial^{α} v_{a b}^{μ} + i v_{a b}^{μ} (A_{b b}^{α} - A_{a a}^{α}) + \\ + \sum_{c}^{'} {\frac{v_{a c}^{α} v_{c b}^{μ}}{ϵ_{c a}} - \frac{v_{a c}^{μ} v_{c b}^{α}}{ϵ_{b c}}}, \\ h_{a b}^{μ α β} = \partial^{β} h_{a b}^{μ α} + i h_{a b}^{μ α} (A_{b b}^{β} - A_{a a}^{β}) + \\ + \sum_{c}^{'} {\frac{v_{a c}^{β} h_{c b}^{μ α}}{ϵ_{c a}} - \frac{h_{a c}^{μ α} v_{c b}^{β}}{ϵ_{b c}}} \end{aligned}

where

\sum^{'}

means that one omits terms when denominator vanishes and latin sub-indices,

(a, b, c)

denote band labels, while upper greek indices

(α, β, γ)

denote space components.

From the above, the second order rectification conductivity is found to be:

σ^{μ α β} (0, ω, - ω) = \frac{e^{3}}{ω^{2}} \sum_{a b c} {f_{a} h_{a a}^{μ α β} +

\begin{aligned} + f_{a b} \frac{h_{a b}^{α} h_{b a}^{μ β}}{ω - ϵ_{a b}} + f_{a b} \frac{h_{a b}^{β} h_{b a}^{μ α}}{- ω - ϵ_{a b}} + f_{a b} \frac{h_{a b}^{α β} h_{b a}^{μ}}{ϵ_{b a}} + \\ + \frac{h_{a b}^{α} h_{b c}^{β} h_{c a}^{μ}}{ϵ_{a c}} (\frac{f_{a b}}{ω - ϵ_{b a}} - \frac{f_{c b}}{ω - ϵ_{b c}}) + \\ + \frac{h_{a b}^{β} h_{b c}^{α} h_{c a}^{μ}}{ϵ_{a c}} (\frac{f_{a b}}{- ω - ϵ_{b a}} - \frac{f_{c b}}{- ω - ϵ_{b c}})} . \end{aligned}

Here

f_{a b} = f_{a} - f_{b}

, with

f_{a}

the Fermi-Dirac occupation of band

a

. We will compare the expression above with the expression from Refs.

[1, 2]

. References

[1, 2]

employed the same gauge of Ref.[9], but assumed that the electrons have an underlying parabolic kinetic energy in addition to the periodic potential, and thus entirely neglects spin-orbit coupling effects. Therefore, the second order perturbed Hamiltonian for a single electron in Refs. [1, 2] has the form:

\begin{aligned} \hat{H} = {\hat{H}}_{0} + \hat{V} (t), A^{α} = \frac{E_{0}^{α}}{ω} \cos ω t \\ \hat{V} (t) = - \frac{e}{m_{0}} {\hat{p}}^{α} A^{α} + \frac{e^{2}}{2 m_{0}} A^{α} A^{α} \hat{1} \end{aligned}

where

m_{0}

is the bare electron mass. The second order rectification conductivity of Refs.

[1, 2]

is therefore found to be:

\begin{aligned} σ_{α β}^{μ} (ω) & = \frac{| e |^{3}}{8 π^{2} ω^{2}} \times \\ \times Re {\sum_{Ω = \pm ω} \sum_{a b c} f_{b a} \frac{v_{a b}^{α} v_{b c}^{β} v_{c a}^{μ}}{(ϵ_{a c} - i δ) (ϵ_{a b} + Ω - i δ))}} \end{aligned}

For purposes of comparing both expressions it is convenient to split them into pieces. The conductivity from Refs.

[1, 2]

we split as follows:

\begin{aligned} {K B}_{1} & = \frac{e^{3}}{ω^{2}} Re [\sum_{a b c} f_{a b} \frac{v_{a b}^{α} v_{b c}^{β} v_{c a}^{μ}}{(ϵ_{a c} - i δ) (ϵ_{a b} + ω - i δ))}] \\ {K B}_{2} & = \frac{e^{3}}{ω^{2}} Re [\sum_{a b c} f_{a b} \frac{v_{a b}^{α} v_{b c}^{β} v_{c a}^{μ}}{(ϵ_{a c} - i δ) (ϵ_{a b} - ω - i δ))}] \end{aligned}

where KB stands for Kraut-Baltz from Refs.

[1, 2]

. We also split the conductivity from Eq.(95) (Ref.[9]) into the following terms:

{P M}_{0} = \frac{e^{3}}{ω^{2}} \sum_{a b} {f_{a} h_{a a}^{μ α β} + f_{a b} \frac{h_{a b}^{α β} h_{b a}^{μ}}{ϵ_{b a}}}

\begin{aligned} {P M}_{1} = \frac{e^{3}}{ω^{2}} \sum_{a b} {f_{a b} \frac{h_{a b}^{α} h_{b a}^{μ β}}{ω - ϵ_{a b}} + f_{a b} \frac{h_{a b}^{β} h_{b a}^{μ α}}{- ω - ϵ_{a b}}}, \\ {P M}_{2} = \frac{e^{3}}{ω^{2}} \sum_{a b c} {\frac{h_{a b}^{α} h_{b c}^{β} h_{c a}^{μ}}{ϵ_{a c}} (\frac{f_{a b}}{ω - ϵ_{b a}} - \frac{f_{c b}}{ω - ϵ_{b c}})}, \\ {P M}_{3} = \frac{e^{3}}{ω^{2}} \sum_{a b c} {\frac{h_{a b}^{β} h_{b c}^{α} h_{c a}^{μ}}{ϵ_{a c}} (\frac{f_{a b}}{- ω - ϵ_{b a}} - \frac{f_{c b}}{- ω - ϵ_{b c}})}, \end{aligned}

where PM stands for Parker-Morimoto-Orenstein-Moore. Also we want to use one more index for each sub-term in these expressions. Namely, we will denote by

{P M}_{1, 1}

the firs term from expression

{P M}_{1}

, and analogously for other terms.

By comparing Eqs.

(90, 91)

and Eqs.

(96 - 97)

, we see that for systems considered in Refs.

[1, 2]

certain additional constrains are satisfied for the matrix elements as a consequence of the underlying parabolic dispersion. More specifically, the interband rectification conductivity from Refs.

[1, 2]

, can be obtained from the formula of Ref.[9], upon enforcing the following relations:

\begin{matrix} {\hat{h}}^{α} = - \frac{ℏ}{m_{0}} {\hat{p}}^{α}, \\ {\hat{h}}^{α β} = \frac{ℏ^{2}}{e^{2}} \frac{e^{2}}{2 m_{0}} δ^{α β} \hat{1} \\ {\hat{h}}^{α β γ} = {\hat{h}}^{α β γ δ} = \dots = 0 . \end{matrix}

From the expressions above one can see that the terms we have labeled

{P M}_{0}

and

{P M}_{1}

would vanish, because the tensor

{\hat{h}}_{a b}^{α β} \sim δ_{a b} δ^{α β}

, is diagonal in the

(a, b)

band indices, and therefore its contribution to the rectification conductivity vanishes after multiplying by the difference of Fermi occupation functions,

f_{a b} = f_{a} - f_{b}

, which are only non-zero for

a \neq b

. Now one can verify that the remaining terms agree between the two theories, specifically as follows:

\begin{aligned} Re [{P M}_{21}] & = \frac{e^{3}}{ω^{2}} Re [\sum_{a b c} \frac{h_{a b}^{α} h_{b c}^{β} h_{c a}^{μ}}{ϵ_{a c}} \frac{f_{a b}}{ω - ϵ_{b a}}] = \\ = \frac{e^{3}}{ω^{2}} Re [\sum_{a b c} \frac{v_{a b}^{α} v_{b c}^{β} v_{c a}^{μ}}{ϵ_{a c}} \frac{f_{a b}}{ϵ_{a b} + ω}] = {K B}_{1} \end{aligned}

Re [{P M}_{32}] = - \frac{e^{3}}{ω^{2}} Re [\sum_{a b c} \frac{h_{a b}^{β} h_{b c}^{α} h_{c a}^{μ}}{ϵ_{a c}} \frac{f_{c b}}{- ω - ϵ_{b c}}] =

= | \begin{matrix} complex \\ conjugation \end{matrix} | = - \frac{e^{3}}{ω^{2}} Re [\sum_{a b c} \frac{h_{b a}^{β} h_{c b}^{α} h_{a c}^{μ}}{ϵ_{a c}} \frac{f_{c b}}{- ω - ϵ_{b c}}] =

= | a \leftrightarrow c | = - \frac{e^{3}}{ω^{2}} Re [\sum_{a b c} \frac{h_{b c}^{β} h_{a b}^{α} h_{c a}^{μ}}{ϵ_{c a}} \frac{f_{a b}}{- ω - ϵ_{b a}}] =

\begin{aligned} = \frac{e^{3}}{ω^{2}} Re [\sum_{a b c} \frac{h_{b c}^{β} h_{a b}^{α} h_{c a}^{μ}}{ϵ_{a c}} \frac{f_{a b}}{ϵ_{a b} - ω}] = \\ = \frac{e^{3}}{ω^{2}} Re [\sum_{a b c} \frac{v_{b c}^{β} v_{a b}^{α} v_{c a}^{μ}}{ϵ_{a c}} \frac{f_{a b}}{ϵ_{a b} - ω}] = {K B}_{2} . \end{aligned}

Expressions for Re[PM

]_{22}

and

Re [P M]_{31}

can be obtained from

Re [P M]_{21}

and

Re [P M]_{32}

respectively by swapping indices

α \leftrightarrow β

. Such a symmetrization automatically implied by conductivity once it is contracted with electric field indices, as is used for current expression in Ref.[2]. Thus we conclude the conductivity from Ref.[2] is equivalent to the sum of terms

Re [P M]_{2}

and

Re [P M]_{3}

Note that in general (for example in the presence of spin-orbit coupling)

Re [P M]_{0}

and

Re [P M]_{1}

do contribute to the conductivity. These terms include not only intra-band terms, but also some contribution to interband shift currents terms that are missing in Refs. [1, 2], and are given by:

\begin{aligned} Re & [\frac{e^{3}}{ω^{2}} \sum_{a b} f_{a b} \frac{i v_{a b}^{α} v_{b a}^{μ} (A_{a a}^{β} - A_{b b}^{β})}{ω - ϵ_{a b}}] + \\ + Re [\frac{e^{3}}{ω^{2}} \sum_{a b} f_{a b} \frac{v_{a b}^{α}}{ω - ϵ_{a b}} \sum_{c}^{'} {\frac{v_{b c}^{β} v_{c a}^{μ}}{ϵ_{c b}} - \frac{v_{b c}^{μ} v_{c a}^{β}}{ϵ_{a c}}}] \\ + & + (\begin{array}{c} α \leftrightarrow β \\ ω \leftrightarrow - ω \end{array}) \end{aligned}

Or in notations used in Refs.

[1, 2]

\begin{aligned} Re [\frac{e^{3}}{2 ω^{2}} \sum_{Ω = \pm ω} \sum_{a b} f_{a b} \frac{i v_{a b}^{α} v_{b a}^{μ} (A_{a a}^{β} - A_{b b}^{β})}{Ω + ϵ_{b a} - i δ}] + \\ + & Re [\frac{e^{3}}{ω^{2}} \sum_{Ω = \pm ω} \sum_{a b} \sum_{c}^{'} f_{a b} \frac{v_{a b}^{α}}{Ω + ϵ_{b a} - i δ} \frac{v_{b c}^{β} v_{c a}^{μ}}{ϵ_{c b}}] + \\ + (α \leftrightarrow β) \end{aligned}

[1] W. Kraut and R. von Baltz, Phys. Rev. B 19, 1548

(1979)

[2] R. von Baltz and W. Kraut, Phys. Rev. B 23, 5590 (1981).

[3] B. I. Sturman and V. M. Fridkin, The Photovoltaic and Photorefractive Effects in Noncentrosymmetric Materials, edited by G. W. Taylor, Ferroelectricity and Related Phenomena (Gordon and Breach Science Publish- To show that the extra shift current contribution, which we demonstrate above, is crucial - we will compute our Quantum Rectification Sum Rule for shift currents described in Refs. [1, 2] and Ref.[9]. With the proper regularization of the infinitesimal imaginary parts of the frequency denominators of both of these references, the frequency integration of interband conductivities are:

\begin{aligned} 2 \frac{ℏ^{2}}{e^{3}} \frac{1}{π} \int_{0}^{\infty} Re [{PM}^{μ β α} (0, - ω, ω)] d ω \\ = & \sum_{a} f_{a} (i {[A^{β} \partial^{α} A^{μ}]}_{a a} + [A^{β}, {[A^{α}, {\bar{A}}^{μ}]}_{a a}) + \\ + (α \leftrightarrow β) \end{aligned}

\begin{aligned} 2 \frac{ℏ^{2}}{e^{3}} \frac{1}{π} \int_{0}^{\infty} Re [{K B}^{γ β α} (0, - ω, ω)] d ω = \\ = \sum_{a b c} f_{a b} \frac{ϵ_{b c}}{ϵ_{a b}} Re [A_{a b}^{α} A_{b c}^{β} A_{c a}^{μ}] . \end{aligned}

Therefore we see that while the expression for interband contribution to the sum rule obtained from Ref.[9] (Eq.(112)) is the same as the one derived in our previous paper Ref.

[10]

, the one that would be derived from Refs.[1, 2] (Eq.(113)) differs from them. As previously discussed, this is because Refs. [1, 2] are missing some terms which only vanish if one assumes a parabolic dispersion which completely neglects the spin-orbit coupling.

ers, 1992), vol. 8 ed.

[4] C. Aversa and J. E. Sipe, Phys. Rev. B

52, 14636

(1995).

J. E. Sipe and A. I. Shkrebtii, Phys. Rev. B

61, 5337

(2000).

[6] J. A. Brehm, S. M. Young, F. Zheng, and A. M. Rappe, The Journal of Chemical Physics 141, 204704 (2014).

[7] T. Morimoto and N. Nagaosa, Science Advances 2 (2016).

[8] N. Nagaosa and T. Morimoto, Advanced Materials 29, 1603345 (2017).

[9] D. E. Parker, T. Morimoto, J. Orenstein, and J. E. Moore, Phys. Rev. B 99, 045121 (2019).

[10] O. Matsyshyn and I. Sodemann, Phys. Rev. Lett.

123

, 246602 (2019).

[11] S. M. Young and A. M. Rappe, Phys. Rev. Lett.

109

, 116601 (2012).

[12] T. Rangel, B. M. Fregoso, B. S. Mendoza, T. Morimoto, J. E. Moore, and J. B. Neaton, Phys. Rev. Lett.

119

, 067402 (2017).

[13] A. M. Cook, B. M. Fregoso, F. de Juan, S. Coh, and J. E. Moore, Nature Communications 8, 14176 EP (2017), article.

[14] T. Morimoto, M. Nakamura, M. Kawasaki, and N. Nagaosa, Phys. Rev. Lett. 121, 267401 (2018).

[15] E. Deyo, L. E. Golub, E. L. Ivchenko, and B. Spivak, arXiv e-prints p. arXiv:0904.1917 (2009), 0904.1917.

[16] J. E. Moore and J. Orenstein, Phys. Rev. Lett.

105

, 026805 (2010).

[17] I. Sodemann and L. Fu, Phys. Rev. Lett.

115, 216806

(2015).

[18] K. Kang, T. Li, E. Sohn, J. Shan, and K. F. Mak, arXiv e-prints arXiv:

1809.08744

(2018), 1809.08744.

[19] Q. Ma, S. Xu, H. Shen, D. MacNeill, V. Fatemi, T. Chang, A. M. Mier Valdivia, S. Wu, Z. Du, C. Hsu, et al., Nature 565, 337 (2019), ISSN 1476-4687.

[20] O. O. Shvetsov, V. D. Esin, A. V. Timonina, N. N. Kolesnikov, and E. V. Deviatov, JETP Letters 109, 715721 (2019), ISSN 1090-6487.

[21] J. Son, K.-H. Kim, Y. H. Ahn, H.-W. Lee, and J. Lee, Phys. Rev. Lett. 123, 036806 (2019).

[22] S.-C. Ho, C.-H. Chang, Y.-C. Hsieh, S.-T. Lo, B. Huang, T.-H.-Y. Vu, C. Ortix, and T.-M. Chen, arXiv e-prints (2019), 1910.07509.

[23] M. Huang, Z. Wu, J. Hu, X. Cai, E. Li, L. An, X. Feng, Z. Ye, N. Lin, K. Tuen Law, et al., arXiv e-prints (2020),

2006.05615

[24] H. Weng, C. Fang, Z. Fang, B. A. Bernevig, and X. Dai, Phys. Rev. X 5, 011029 (2015).

[25] B. Q. Lv, N. Xu, H. M. Weng, J. Z. Ma, P. Richard, X. C. Huang, L. X. Zhao, G. F. Chen, C. E. Matt, F. Bisti, et al., Nature Physics 11, 724 (2015), ISSN 1745-2481.

[26] S. Huang, S. Xu, I. Belopolski, C. Lee, G. Chang, B. Wang, N. Alidoust, G. Bian, M. Neupane, C. Zhang, et al., Nature Communications

6, 7373 E P

(2015), article.

[27] S.-Y. Xu, I. Belopolski, N. Alidoust, M. Neupane, G. Bian, C. Zhang, R. Sankar, G. Chang, Z. Yuan, C.-C. Lee, et al., Science 349, 613-617 (2015).

[28] Y. Zhang, Y. Sun, and B. Yan, Phys. Rev. B 97, 041101 (2018).

[29] W. Shockley and H. J. Queisser, Journal of Applied Physics 32, 510-519 (1961).

[30] K. K. Kohli, J. Mertens, M. Bieler, and S. Chatterjee, J. Opt. Soc. Am. B 28, 470 (2011).

[31] Z. Addison and E. J. Mele, arXiv e-prints (2020), 2005.11550.

[32] S. Patankar, L. Wu, B. Lu, M. Rai, J. D. Tran, T. Morimoto, D. E. Parker, A. G. Grushin, N. L. Nair, J. G. Analytis, et al., Phys. Rev. B 98, 165113 (2018).

[33] O. Coddington, J. L. Lean, P. Pilewskie, M. Snow, and D. Lindholm, Bulletin of the American Meteorological Society 97, 1265 (2016).

[34] A. M. Glass, D. von der Linde, and T. J. Negran, Applied Physics Letters 25, 233-235 (1974).

[35] A. Zangwill, Modern Electrodynamics (Cambridge University Press, 2012).

[36] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77,3865 (1996).

[37] K. Koepernik and H. Eschrig, Phys. Rev. B

59, 1743

(1999)

[38] Y. Zhang, H. Ishizuka, J. van den Brink, C. Felser, B. Yan, and N. Nagaosa, Phys. Rev. B 97, 241118 (2018).

[39] J. Ibañez Azpiroz, S. S. Tsirkin, and I. Souza, Phys. Rev. B 97, 245143 (2018).

[40] Y. Zhang, F. de Juan, A. G. Grushin, C. Felser, and Y. Sun, Phys. Rev. B 100, 245206 (2019).

[41] G. B. Ventura, D. J. Passos, J. M. B. Lopes dos Santos, J. M. Viana Parente Lopes, and N. M. R. Peres, Phys. Rev. B 96, 035431 (2017).

[42] T. Holder, D. Kaplan, and B. Yan, Phys. Rev. Research 2, 033100 (2020).

[43] F. de Juan, Y. Zhang, T. Morimoto, Y. Sun, J. E. Moore, and A. G. Grushin, Phys. Rev. Research 2, 012017

(2020)

The Berry Phase Rectification Tensor and The Solar Rectification Vector

1. INTRODUCTION

2. THE BERRY PHASE RECTIFICATION TENSOR

3. L i A s S e 2 L i A s S e 2 LiAsSe_(2)\mathrm{LiAsSe}_{2}LiAsSe2

4. THE SOLAR RECTIFICATION VECTOR

5. A. Autocorrelations for black-body radiation

6. B. Definition of Solar Rectification Vector

7. MATERIAL PENETRATION

8. BULK RECTIFICATION OF TaAs AND L i A s S e 2 L i A s S e 2 LiAsSe_(2)\mathrm{LiAsSe}_{2}LiAsSe2

9. DISCUSSION

10. ACKNOWLEDGEMENTS

11. APPENDIX A: SECOND ORDER RESPONSE TO ELECTRIC FIELDS IN TIME AND FREQUENCY DOMAINS

12. APPENDIX B: SUM RULE COMPUTATIONS

13. APPENDIX C: COMPUTATION DETAILS

14. APPENDIX D: EQUIVALENCE WITH OTHER FORMALISMS IN THE LITERATURE

15. A. Comparison with Refs. [ 5 , 39 ] [ 5 , 39 ] [5,39][5,39][5,39]

16. B. Comparison with Refs. [ 1 , 2 , 9 ] [ 1 , 2 , 9 ] [1,2,9][1,2,9][1,2,9]

Recommended for you

Report Article

3. ${L i A s S e}_{2}$

8. BULK RECTIFICATION OF TaAs AND ${L i A s S e}_{2}$

15. A. Comparison with Refs. $[5, 39]$

16. B. Comparison with Refs. $[1, 2, 9]$