1. Introduction
Semimetals are materials which can support gapless quasiparticle excitations in two or three dimensions, in the vicinity of isolated band touching points in the Brillouin zone, thus possessing discrete Fermi points (rather than Fermi surfaces). They come in different varieties, for example, the Fermi points may appear at linear band crossings (e.g., graphene, Weyl semimetals), or at quadratic band crossings [
1,
2] (e.g., Luttinger semimetals). A more non-trivial example of such semimetal is the double-Weyl semimetal, which consists of two bands touching each other linearly along one momentum direction, but quadratically along the remaining directions [
3,
4,
5,
6,
7]. Some of these three-dimensional (3d) semimetals (e.g., Weyl and double-Weyl semimetals) possess a non-zero Berry curvature at the Fermi nodes. In this paper, we focus on the 3d double-Weyl semimetals [
2,
8,
9], which, in the momentum space, have double the monopole charge of Weyl semimetals.
A double-Weyl semimetal can be realized by applying a Zeeman field to an isotropic Luttinger semimetal [
2]. They are also predicted to appear [
9,
10] in SrSi
, and in the ferromagnetic phase of HgCr
Se
. Our aim is to study the circular photogalvanic effect (CPGE), also known as chiral photocurrent. The CPGE refers to the DC current that is generated as a result of shining circularly polarized light on the surface of an optically active metal [
11,
12,
13,
14]. In fact, the CPGE refers to the part of the photocurrent that switches sign with the sign of the helicity of the incident polarized light. This is a non-linear response, as it is second order in the applied AC electric field, and at low frequencies, it depends on the orbital Berry phase of the Bloch electrons. Hence, CPGE is a measure of the topological charge at a Fermi node possessing a non-trivial Berry curvature.
The quantization of the CPGE has been demonstrated in earlier works for the topological Weyl nodes [
15,
16]. In this paper, we will consider the issue of quantization of CPGE for the double-Weyl nodes. First, we will show that in the absence of interactions, the CPGE is indeed proportional to the topological charge of the node at low-enough frequencies. Secondly, we will examine the effect of Hubbard interactions on this quantized value.
2. The Continuum Hamiltonian for a Double-Weyl Semimetal
The Hamiltonians describing a pair of double-Weyl nodes can be written in the form [
2,
8,
9]
with
Here,
are the three Pauli matrices, and the “±” sign reflects the two opposite chiralities of the two nodes. The energy eigenvalues are:
For each the given two-band Hamiltonians, we can define an
Berry curvature, which is analogous to a magnetic field in momentum space. This Berry curvature is given by:
where
. It is easy to check that this magnetic field is divergenceless
as long as it is computed in regions away from the points of singularity where
. The band touching point is such a singularity, where we have:
Thus, each double-Weyl node is a source of two Berry flux quanta. These nodes come in pairs, sourcing equal and opposite flux quanta, such that the sum of Berry flux quanta from both the double-Weyl nodes vanishes, which is the desired physical scenario as the Brillouin zone is a closed manifold without any boundary through which no net flux can emanate.
3. Quantization of CPGE in the Absence of Interactions
The CPGE tensor is defined as [
15,
16]:
where
is the electric charge. To perform the integrals, we change variables as follows:
All non-diagonal components
evaluate to zero. Clearly, we see that
where is
the monopole charge of the corresponding double-Weyl node. The time derivative of the injection current is defined as the second-order response
to an electric field
. Therefore, the CPGE is also quantized.
Now let us compute the second-order photocurrent from the field-theoretic definition, using Feynman diagrams. First, we need the three components of the paramagnetic current operator (using
), which are given by:
From now on we will drop the “±” subscript/superscript and concentrate only on the double-Weyl node with charge , unless stated otherwise. This is justified when the DC contribution to the photocurrent can be calculated separately for each node, such as when the nodes are well separated in the momentum space.
The expression for the second-order photocurrent is given by:
where
and the contributions
and
are given by Feynman diagrams of the type shown in
Figure 1. In the second line, we have used the relation between the electric field and the vector potential, which is:
.
We compute the analytical expressions for
in the Matsubara formalism, such that
where
T is the temperature,
n is an integer, and
. In the zero-temperature limit, we can use
Furthermore, from the expression for
, we can obtain
by using the relation:
In the absence of interactions, we can calculate the contributions from each node separately. The Green’s function for the first double-Weyl node is given by:
where we have introduced the projectors
onto the conduction (“+”) and the valence (“−”) bands, and have chosen the chemical potential
to be negative for definiteness (i.e.,
). Similarly, the Green’s function for the second double-Weyl node is given by:
where we have chosen
for definiteness.
Performing all the integrals, we finally get:
for
One can check that
, and hence the computation of
is sufficient to know all the non-zero components of
.
We need to find the physical response through the analytical continuation of the above expressions from Matsubara frequencies to real frequencies. This is a subtle procedure which should be carried out carefully. Choosing
for definiteness, the analytical continuation is performed by taking [
17,
18]
The logarithms then transform according to
We then need to set
with
. After the analytical continuation, we find that in this limit,
An identical contribution comes from
(on using Equation (
14)). Adding these together, we find that the current expression in Equation (
12) reduces to:
In the time domain, this corresponds to
This agrees with Equation (
10).
This result from the non-interacting case has been obtained for the first double-Weyl node with the chemical potential
. Analogously, for the second node, we would obtain:
Consequently, in the frequency range , only the first node contributes to the CPGE, while the contribution from the second node is zero due to Pauli blocking.
4. Corrections to the Quantized CPGE Due to Short-Ranged Hubbard Interactions
In this section, we consider the first-order perturbative corrections originating from four-fermion interactions. The interaction Hamiltonian for short-ranged Hubbard interactions is given by:
where
is the Hubbard interaction strength (positive
corresponds to the attractive interaction), and
denotes the fermion field with nodal index
and pseudospin index
s. The first and the second terms describe the intranodal and internodal scattering processes, respectively. These are shown diagrammatically in
Figure 2. In the diagrams, we have used a solid line to represent the Green’s function for the first double-Weyl node, and a dashed line to depict the Green’s function for the second double-Weyl node. In the following subsections, we will compute the first-order self-energy and vertex corrections due to the Hubbard interactions.
4.1. First-Order Self-Energy Corrections
The contributions to the first-order self-energy correction are given by the Feynman diagrams shown in
Figure 3. For the short-ranged Hubbard interaction, scatterings between double-Weyl nodes of opposite chiralities must be taken into account, which are given by the second term of Equation (
24). The analytic expression for
Figure 3a reads as:
where
is the number of holes below the double-Weyl point in the first node. In a similar fashion, the contribution from
Figure 3b evaluates to:
with
denoting the number of electrons above the double-Weyl point in the second node.
Finally, the contributions from
Figure 3c,d evaluate to:
resulting in the total self-energy
The effect of this self-energy is to simply shift the chemical potential by an amount
Clearly, this does not change the CPGE current, as it only modifies the frequency range where the quantized value of the CPGE is valid.
4.2. First-Order Vertex Corrections
The Feynman diagrams contributing to first-order vertex corrections are shown in
Figure 4. When the vertex
, with the external Matsubara frequency set to
for definiteness,
Figure 4a contributes as:
The only non-vanishing contribution from
Figure 4a comes for
, which gives:
where
is the UV momentum cutoff.
The contribution from the diagram in
Figure 4b is analogous, but has an overall opposite sign due to the opposite chirality of the second node, and with
:
Adding these two contributions together, we find that for the first node, the vertex with
(and external frequency
) is renormalized according to:
which is finite and does not contain the UV cutoff anymore. This gives the correction:
Performing the analytical continuation
, and setting
, we find that this contributes as:
which leads to the correction
for the current in the
direction. Here, we have neglected the corrections to the chemical potentials, since they only change the frequency range within which CPGE for the non-interacting case is non-zero.
In a similar fashion, we get:
By symmetry in the
plane, we infer that
Due to the intrinsic anisotropy of the problem, it is not surprising that the corrections for the current in the direction is different from that in the plane.
5. Summary and Outlook
We have computed the CPGE for the double-Weyl semimetal, first in the absence of interactions and then in the presence of short-ranged Hubbard interactions. In the non-interacting case, for low-enough frequency ranges of the applied electric field, the CPGE gets contribution only from one double-Weyl node and has a quantized value proportional to the topological charge of the corresponding node. However, switching on Hubbard interactions affects this result, destroying the quantization. This is similar to the results found for the case of CPGE currents in Weyl semimetals [
18]. The only difference is that the corrections for the current in the
direction is different from that in the
plane, due to the anisotropic dispersion of the starting Hamiltonian. These results imply that unlike the quantum Hall effect in gapped phases or the chiral anomaly in field theories, the quantization of the CPGE in topological semimetals is not protected.
In future, it will be interesting to look at the corrections coming from the Coulomb interactions. The computations will be cumbersome for this case compared to the Weyl semimetal, due to the anisotropic dispersion of the double-Weyl Hamiltonian. It will also be interesting to see the effect of short-ranged correlated disorder on the CPGE, using the well-known techniques [
19,
20,
21].