Brief summary of the charge current, spin current, spin, and charge density induced by localized spin dynamics on the disordered surface of a TI. The charge current and spin density have both a local and nonlocal contribution. The spin current and charge density are described by the nonlocal contribution.
We theoretically show spin and charge transport on the disordered surface of a three‐dimensional topological insulator with a magnetic insulator when localized spin of the magnetic insulator depends on time and space. To ascertain the transports, we use a low‐energy effective Hamiltonian on the surface of a topological insulator using the exchange interaction and calculate analytically using Green's function techniques within the linear response to the exchange interaction. As a result, the time‐dependent localized spin induces the charge and spin current. These currents are detected from change in the half‐width value of the ferromagnetic resonance of the localized spin when the magnetic resonance of the localized spin is realized in the attached magnetic insulator. We also show spin and charge current generation in a three‐dimensional Weyl–Dirac semimetal, which has massless Dirac fermions with helicity degrees of freedoms. The time‐dependent localized spin drives the charge and spin current in the system. The charge current as well as the spin current in the Weyl–Dirac system are slightly different from those on the surface of the topological insulator.
- Spin pumping
- Spin–momentum locking
- Surface of topological insulator
- Weyl–Dirac semimetal
- Massless Dirac fermions
A crucial issue in spintronics is the generation and manipulation of a charge and spin current by magnetism, since these mechanisms can be applicable to magnetic devices. One way to generate charge and spin flow is called “spin pumping,” which pumps from the angular momentum of a magnetization's localized spin into that of electrons through the dynamics of magnetization as well as spin-orbit interactions [1, 2]. No other way of doing this has so far been discussed in the field of metallic spintronics.
Ever since the discovery of a topological insulator (TI) [3–6], spintronics using topology has been studied. A TI has a gapless surface, its bulk is insulating, but its surface is metallic as a result of two‐dimensional massless Dirac fermions on the surface [4–6]. Because of spin–orbit interactions on the surface the spin and momentum of Dirac fermions are perfectly linked to each other. The relation between a TI's spin and momentum is dubbed “spin–momentum locking,” and the direction in which they travel is perpendicular to each other. Because of spin–momentum locking, unconventional spin‐related phenomena—such as magnetoresistance [7–15], the magnetoelectric effect [16–21], diffusive charge–spin transport [22–25], and the spin pumping effect [26–31]—have so far been the only phenomena theoretically and experimentally studied.
Of the unconventional phenomena on the surface of a TI, spin pumping is one of the most interesting when it comes to spintronics. Here spin pumping on the surface is different from that in metals. As a result of spin–momentum locking on the surface the localized spin plays the role of a vector potential, whereas time‐differential localized spin effectively plays the role of an electric field acting on electrons on the surface of the TI [16–19]. As a result, even in the absence of an applied electric field, the charge current is generated by time‐dependent localized spin as shown in Figure 1. The induced charge current flows along z×∂tS [28–31]. Actually, this spin pumping effect on the surface has been experimentally demonstrated in the junction of the TI by attaching magnetic materials [29–31] by changing the half‐width value when ferromagnetic resonance is realized in the attached magnetic material [32, 33].
Recently, it has been reported that localized spin on the surface of a TI subject to magnetism depends on space and its spin texture seems to be a magnetic domain wall . It is predicted that in the presence of a spatial‐dependent spin structure the charge current, which reflects the spin structure, is induced . Moreover, the spin current as well as the charge density are induced when an inhomogeneous spin structure exists on the surface. Detailed results are shown in Section 2. This study may help the study of spin pumping on the surface of a TI with inhomogeneous spin textures [35, 36].
Recently, the next generation of spintronics has been theoretically and experimentally studied in Weyl–Dirac semimetals. A Weyl–Dirac semimetal possessing three‐dimensional massless Dirac fermions has attracted much attention in condensed matter physics [37–39]. Such a semimetal has been experimentally demonstrated [40–45]. In addition, Weyl–Dirac semimetals have been theoretically predicted in a superlattice heterostructure based on the TI. Such a heterostructure has been realized in the GeTe/Sb2Te3 superlattice .
Spin‐momentum locking occurs in a Weyl–Dirac semimetal, but the locking is slightly different from that on the surface of the TI. As a result of spin–momentum locking, the spin polarization (charge density) and the charge current (spin current) are linked to each other. Moreover, Dirac fermions have helicity degrees of freedom, which are decomposed into left‐ and right‐handed fermions. Note that the total charge flow of Dirac fermions of left‐ and right‐handed Weyl fermions is preserved. In a Weyl–Dirac semimetal the anomaly‐related effect, which is discussed in the field of relativistic high‐energy physics, has also been discussed in condensed matter physics. Studies up to the moment have asserted that the charge current is generated by magnetic properties with helicity degrees of freedoms [47–54]. Our goal is to introduce the helicity‐dependent spin pumping effect, one of the characteristic properties of Dirac fermions (as shown in Section 3).
2. Spin and charge transport due to spin pumping on the surface of a topological insulator
We will calculate the charge and spin current due to spin pumping on the surface of a TI with an attached magnetic insulator (MI) (as illustrated in Figure 1). To do so, we consider the following low‐energy effective Hamiltonian, which describes the surface of the TI with localized spin of the MI :
where HTI‐surface, Hex, and Himp are the low‐energy effective Hamiltonian on the surface of the TI, the spin–exchange interaction between the localized spin of the attached MI and the spin of conduction electrons, and the nonmagnetic impurity scattering on the surface of the TI, respectively. HTI‐surface is given by
where vF,B, EF, k, and σ are the velocity of bare electrons, the Fermi energy, momentum on the surface of the TI, and the Pauli matrix in the spin space, respectively; and Ψk = t(Ψ↑, Ψ↓) and Ψ† are the annihilation and creation operators of the electrons with up‐spin and down‐spin index (↑, ↓). Hex is represented by
where Jex is the coefficient of the exchange interaction; and S is the localized spin of the MI on the surface. Here, we take into account the space‐ and time‐dependence of S. Himp is given by
Nonmagnetic impurity scattering is taken into account for a delta function type [19, 21–24] and is considered within the Born approximation. Because of impurity scattering the Fermi velocity of bare Dirac fermions in Eq. (2) is modified by vF,B → vF [23, 24, 28]. We use vF in what follows.
2.2. Charge and spin current due to localized spin dynamics
We will calculate the charge and spin current as well as the charge and spin density due to magnetization dynamics on the surface of the TI. They are given using the Keldysh–Green function and the lesser component of Green's function  as
As a result of spin–momentum locking the charge current <ji> and spin current <jis,α> are proportional to the spin density <si> and charge density <ρ>, respectively, where < > = denotes the expectation value in H and e(<0) is the charge of Dirac fermions. These relations are derived from the Heisenberg equation. The superscript and subscripts in <jis,α> show the direction of flow and spin of the spin current, respectively. In what follows the charge and spin current are considered within the linear response to Hex. This assumption can be a good approximation because the energy scale of the exchange interaction is smaller than that of the bandwidth of the low‐energy effective Hamiltonian on the surface of the TI [28, 31]. We also assume the TI has a metallic surface of and that a finite Fermi surface exists on the surface; that is, h/(EFτ) ≪ 1 where h is the Planck constant and τ is the relaxation time of nonmagnetic impurity scatterings on the surface.
The charge and spin current are represented by
where q and Ω are the momentum and frequency of the localized spin Sjq,Ω, respectively; L2 is the area of the surface; and Πj is the response function. These currents are obtained from Πj:
where fω, gω,kR = ([gω,kA]†) = [hw/(2π) + EF – hvF/(2π)(k × σ)z + ih/(4πτ)]‐1, and Λj are the Fermi distribution function, the retarded (advanced) Green function, and the vertex correction of the ladder diagram, respectively. The above functions are estimated in the regime Ωτ ≪ 1 and qvFτ ≪ 1, which are satisfied when the dynamics of the localized spin is lower than the terahertz regime and the spatial gradient of the localized spin is slow compared with the electron mean free path. Then, by expanding Ω and q within Ωτ ≪ 1 and qvFτ ≪ 1, the vertex function Λj is represented by Λj = [Γ + Γ2 + Γ3 + …]jnσn = [ (1 – Γ)‐1]jnσn with
where Γjn (||Γjn|| < 1) is the 3 × 3 matrix (j, n = 0, x, y). As a result, Πj is given by
where νTIe is the density of state on the surface of the TI; <S||>D is the convolution between S|| and the diffusion propagator D on the surface; and S|| is the in‐plane localized spin S||= S – Szz as
Eq. (10) shows that time‐dependent localized spin induces the charge and spin current and that they can be decomposed into local and nonlocal contributions. The first term in Eq. (10) is the charge current due to spin dynamics at that position on the surface; its direction is along the z × ∂tS|| direction [28–31]. On the other hand, the second term in Eq. (10) is the charge current due to both time‐ and spatial‐dependent localized spin; its direction is along the [z × (z × ∂)] [∂ × ∂t<S||>D]z direction, where <S||<D indicates the nonlocal contribution from the localized spin. Moreover, this term is zero when the spin texture is spatially uniform. Therefore, this charge current is caused by diffusion with the dynamics of the spatial inhomogeneous spin texture. It is noted that this second term is also described by the spatial gradient of the charge density.
Because of spin–momentum locking on the surface of the TI, spin polarization is given by <j> = 2evF(z × <s>). Thus, the properties of spin polarization are similar to those of the charge current.
Eq. (11) shows that the spin current is induced by the time‐ and spatial‐dependent S||. Then, the charge density is also generated because of spin–momentum locking. In addition, the second term of the charge current <j>nonlocal = 2evFνTIeJexτl2[z × (z × ∂)] [∂ × (∂t<S||>D) ]z is given by the spatial gradient of the charge density, whereas <j>nonlocal is proportional to the spatial gradient of the spin current  as
This result shows that the relation between the charge and spin current is different from that in the metallic spintronics system . The spin current is proportional to the charge current and spin current flow is perpendicular to the charge flow and its spin polarization.
Note that no out‐of‐plane localized spin Szz contributes in a dominant way to charge and spin current generation. The reason is σz does not couple with momentum p in Eq. (2). Therefore, we believe a warping effect [1–4] on the surface of the TI is likely. Spin polarization along the z‐direction is also generated by spin dynamics.
2.3. Spin torque
Based on these results, we look at localized spin dynamics after generation of the charge current and spin polarization on the surface . We assume there is an external static and AC magnetic field on the surface—as shown in Figure 2(a). The static magnetic field arranges the localized spin texture and the AC magnetic field triggers its spin dynamics. The propagation direction of the microwave is parallel to the static magnetic field, which is along the y‐axis. The texture is called the “longitudinal conical spin order” .
After applying the microwave the dynamics of the localized spin of the MI is induced by the in‐plane AC magnetic field of the microwave. The dynamics of the spin induces the charge and spin current. Then, from <j> = 2evF(z × <s>), Eq. (10), and spin–momentum locking, spin polarization is also generated. Note that spin polarization induced in this way acts on the localized spin as an effective exchange field for localized spin. This contribution is given from the Landau–Lifshitz–Gilbert (LLG) equation of motion [32, 33] on the surface of the TI:
where M = -M(S/S) is magnetization of the MI; γ is the gyromagnetic ratio; μ is permeability; αG is a Gilbert damping constant; H = H0 + HAC are the applied magnetic fields on the surface; and H0 and Hac denote the static and AC magnetic fields, respectively. The first term of Eq. (13) denotes the field‐like torque that drives the dynamics of M. The second term is the damping torque that prevents its dynamics—these terms are schematically illustrated in Figure 2(b). The third term Te is the spin torque due to spin polarization of Dirac fermions on the surface of the TI (<s>). It is given by Te =2Jexa2 (M × <s>), where a is the lattice constant on the surface.
Since Te is proportional to <s>, this torque is decomposed into local (TeL) and nonlocal spin torque (TeD) terms as
where M|| = -M(S||/S) is in‐plane magnetization; and κ is the dimensionless coefficient.
As a result of spin–momentum locking the torque is Te ∼ M × <j>. This relation is useful as a way of detecting nonzero <j> [29–31]. Moreover, nonlocal spin torque can be described by T<eD ∼ M × <j>nonlocal . Eq. (12) shows this torque may well be affected by spin current contribution, whereas the relation between spin torque and spin current can be used as a way to detecting <jis,α> on the surface:
We now consider magnetization when it is given by M = ( Mx, My, Mz) = [m cos(q.x – Ωt), My, m sin(q.x – Ωt)] with q = (qx, qy), m2 ≪ My, My ∼ M, and ∂tMy ∼ 0. Such a time‐dependent magnetic structure is assumed to be realized when the magnitude of the static magnetic field is larger than that of the AC magnetic field HAC=(hx, 0, hz); that is, hx, hz ≪ H0. Then, substituting this M into Eqs. (13) and (14), the LLG equation becomes a linear approximation:
where ωH = γμH and ωM = γμM are the angular frequencies of H0 and M, respectively; and κqy2l2∂t<M||,x>D indicates the contribution from TeD. Substituting <M||,x>D into Eq. (13), we obtain magnetic permeability around the surface:
where χxx and χzx denote the longitudinal and transverse magnetic permeability as
The coefficients κ'q,Ω and ζq,Ω are given by
Hence, κ'q,Ω and ζq,Ω depend on q and Ω. If q = 0, κ'q,Ω and ζq,Ω go to κ'q,Ω → κ and ζq,Ω → 1, respectively. Then, the magnetic permeability is isotropic. In q = 0, χxx and χzz go to the same value and χxz = –χxx is satisfied. On the other hand, in nonzero q the magnetic permeability is anisotropic: χxx ≠ χzz.
Figure 3(a) shows the dependence of Ω on the imaginary part of the longitudinal magnetic permeability for several q. The resonance frequency Ωr and the half‐width value ΔΩ are slightly changed by the nonzero q. Figure 3(b) shows change in the resonance frequency for several fH in greater detail—fH is the frequency described by the static magnetic field H0. The resonance frequency decreases with increasing ql from ql = 0 into ql = q0l, where q0 satisfies q02l2 ∼ Ωτ. In addition, when ql takes a large value the resonance frequency returns to that frequency at ql=0. Figure 3(c) shows the momentum dependence of the half‐width value on renormalization by the angular frequency of the applied magnetic field for several ωH =2πfH frequencies when realistic parameters on the surface of the TI are taken into account . Change in the normalized half‐width value ΔΩ/ωH(ql) decreases from ΔΩ/ωH(ql = 0) = 0.038 to ΔΩ/ωH(ql ∼ 0.1) = 0.030 with increasing ql.
Any change in the half‐width value ΔΩ indicates an induced charge current on the surface of the TI because TeL + TeD is proportional to the induced charge current [29–31]. In addition, the half‐width value ΔΩ of TeD represents a contribution from the spatial gradient of the spin current. These relations can be useful for the detection of spin current. For example, let us suppose that in a finite ql, ΔΩ includes the contribution from the whole of the spin torque Te and the contribution from the nonlocal term TeD is reduced by the large magnetic field. In a higher magnetic field the inhomogeneous spin texture would be expected to be aligned along the magnetic field direction and the spin texture to be spatially uniform; hence, the value of ql decreases and reaches zero. Then, the half‐width value ΔΩ(TeD = 0) has no contribution from TeD. The difference ΔΩ(TeD ≠ 0) – ΔΩ(TeD = 0) includes a pure contribution from TeD.
Any contribution from the applied magnetic field is of course a concern. Note that an in‐plane static magnetic field contributes no finite charge current generation or spin polarization , whereas the contribution from the applied magnetic field is negligible.
3. Spin pumping in a Weyl–Dirac semimetal
We now consider spin pumping in a Weyl–Dirac semimetal subjected to magnetism. To do this, a Weyl–Dirac semimetal subjected to spin and momentum locking, such as a supperlattice hetrod structure constructed from a TI/normal insulator/TI [39, 46 – 48], is considered. The low‐energy effective Hamiltonian describing the Weyl–Dirac semimetals that have a spin–exchange interaction is given by [50–52]
The first term takes the form
where Ψk = t(Ψ↑, +Ψ↓, +Ψ↑, –Ψ↓,–) is the annihilation operator of an electron with spin (↑, ↓) and helicity τ(+, –) degrees of freedom; τα=0,x,y,z and σα=0,x,y,z are Pauli matrices of the helicity and spin, respectively; vF is the Fermi velocity; μ is the chemical potential of the Weyl–Dirac semimetal; and b and μ5 denote the difference in position between each Weyl cone in momentum and energy space, respectively. In Dirac semimetals, we set b = 0 and μ5= 0.
The second term of Eq. (20) indicates the spin–exchange interaction:
where Jex is the coefficient of the exchange interaction; and S is the localized spin that depends on space and time in the Weyl semimetal. This exchange interaction is independent of the helicity index. Himp shows the nonmagnetic impurity scattering in the Weyl–Dirac semimetal:
This Hamiltonian is similar to that on the surface of the TI—see Eq. (4). In the following calculations, Hex and Himp are treated as perturbations within the same formalism as laid out in Section 2.2. This treatment is allowed when the energy scale of the exchange interaction is smaller than the bandwidth of Eq. (21). Then, the low‐energy effective Hamiltonian in Eq. (21) gives a good approximation.
3.2. Response function within the linear response to the exchange interaction
To calculate the charge and spin current within the linear response to Hex, we use Green's functions. From Eq. (20) the charge current <ji> (i = x, y, z) can be defined by
where < > denotes the expectation value in Eq. (20). Such a charge current can be decomposed as <j> = <j+> + <j–>, where <j±> = ±evF<Ψ†±σΨ±> and Ψ±=t(Ψ↑, ±Ψ↓,±) is the annihilation operator around each helicity τ = ±1. Since there is no mixing between Ψ†↓,+ and Ψ↑,– in the Hamiltonian in Eq. (20) the charge current around each helicity can be calculated separately. Because of spin–momentum locking in the Weyl–Dirac semimetal the charge current links to the spin polarization. The charge current of each helicity is proportional to the spin polarization of each helicity as
The spin current in the Weyl‐Dirac semimetal can be defined from the Heisenberg equation for the spin operator:
where js,αi is the spin current operator; and Tr,α is the spin relaxation term. Spin current density can also be decomposed into the spin current density of each helicity and can be calculated separately:
The superscript and subscripts of js,αi denote the direction of flow and spin of the spin current of each helicity, respectively. The direction of spin is perfectly parallel to that of the flow, and the spin current density is proportional to the charge density of each helicity.
The relaxation term can also be decomposed by <Tr>= <Tr,+> + <Tr,–> and can be given by <Tr,+>:
The charge and spin current can be obtained by calculating the response function (Figure 4):
where q and Ω are the momentum and frequency of the localized spin, respectively; L3 is the volume of the system; and Πj± is the response function within the helicity index:
where fω, gω,k,±R= ([gω,k,±A]†) = [hw/(2π) + μ± – hvFk.σ/(2π) + ih/(4πτ±)]–1, and Λj± are the Fermi distribution function, the retarded (advanced) Green's function, and the vertex correction of the ladder diagram at each helicity ±, respectively; and τ± is the relaxation time of nonmagnetic impurities in the material, which depends on the helicity, because it depends on the chemical potential μ±. When Ωτe,+ ≪ 1 and qvFτe,+ ≪ 1 are satisfied the vertex function Λj± is given by Λj±= [Γ± + Γ±2 + Γ±3 + …]jnσn = [(1 – Γ±)–1]jnσn with
where Γjn± (||Γjn±|| < 1) is a 4 × 4 matrix (j, n = 0, x, y, z). As a result, Πj± can be represented by
where the above matrix component (m, m'= x, y, z) is
As a result, the charge current and spin current can be obtained by
where νWS± and <ρ±> are the charge density of state and the induced charge density of each helicity, respectively:
where <S>D,± is the convolution between the localized spin and diffusive propagation of each helicity. <S>D,± can be defined by
Note that the above results are obtained when b = 0. These results are easily generalized when b = 0. Because b behaves like a static Zeeman field acting on the whole of the band of HWeyl, it can shift as a result of Pauli paramagnetism. However, b cannot drive a net current because b is static . On the other hand, the dynamics of the localized spin are only effective near the Fermi surface, the structure of which does not depend on b. Therefore, we obtain the same charge and spin current in Eqs. (33)–(35) even when b ≠ 0.
3.3. Charge and spin current due to spin pumping effects
From the above results the total charge current can be given by
The charge current is triggered by the dynamics of localized spin. The first term indicates the local term of the dynamics of localized spin. Its direction is parallel to ∂tS. The second term shows the nonlocal term, which is generated by the convolution between localized spin and the diffusion propagator. A nonzero nonlocal term is given when the spin texture of localized spin depends on space. Thus, the driving force needed to induce the charge current is the same as localized spin dynamics, which plays the role of driving the charge current on the surface of the TI. Note that the charge current in a Weyl–Dirac semimetal depends strongly on the valley index. When the difference between τ+ and τ– as well as between νWS+ and νWS– is realized, there is a population imbalance between two of the bands of each helicity. Then, nonzero <j> is induced when νWS+τe,+ ≠ νWS–τe,– or μ5 = 0.
Note that the property of the charge current at each helicity links to that of spin polarization because of spin–momentum locking. However, after summation of the indices of helicity the relation between total charge current and total spin polarization is changed. As a result, even in the absence of population imbalance, nonzero spin polarization can be given by Eq. (25) as
Local spin polarization is along ∂tS and nonlocal spin polarization is along the spatial gradient of ∂t(∂iSi).
Total spin current can be generated by spatial divergence of localized spin dynamics, where localized spin is the convolution with diffusion. As a result, the spin current can be regarded as a nonlocal spin current. Note that a nonzero spin current is generated when localized spin depends on space. Such a spin current becomes nonzero even in the absence of population imbalance.
The nonlocal spin current in this case is obtained from the diffusive motion of spin density, which is driven by the dynamics of localized spin, where localized spin depends on space. The spin diffusive motion of each helicity can be given from Eqs. (25), (34), (37) as
As a result, the diffusive motion of total spin becomes:
Time‐dependent and spatial‐dependent localized spin, ∂[∂t(∂iSi)], triggers the diffusive motion of total spin density. When there is no population imbalance between each helicity, diffusive motion can accompany total spin density without any charge flow. Hence, a pure spin current can be generated.
Our results on the charge and spin current due to spin pumping on the surface of a three‐dimensional TI (Section 2)  and in the bulk of a three‐dimensional Weyl–Dirac semimetal (Section 3)  are summarized in this chapter.
Section 2 summarizes our results on spin pumping on the surface of a TI attached to an MI. The results are calculated using the standard Keldysh–Green function method within the linear response to the exchange interaction between the conduction spin and localized spin of an MI. The purpose of this work is to derive charge and spin current generation due to localized spin dynamics on the disordered surface of a TI; in particular, when the localized spin depends on space on the surface. The main results on the surface of a TI are summarized in Table 1. Time‐dependent localized spin on the surface is a prerequisite to obtaining nonzero charge and spin current generation. Moreover, Table 1 shows that when the spin texture is spatially inhomogeneous, not only the local charge but also the nonlocal charge and spin current are generated by time‐dependent localized spin. The flow and spin polarization of the spin current are perfectly perpendicular to each other because of spin–momentum locking. The magnitude of the spin current is proportional to the charge density, which is induced by divergence between time‐dependent localized spin and the diffusive propagator on the surface—see Eqs. (10) and (11). Such pumping effects are caused by time‐dependent localized spin, which plays a role in driving the charge current and can be regarded as an effective electric field EsTI.
Recently, it has been reported that the localized spin texture at the junction of the TI/MI is spatially inhomogeneous . We suppose that the spin current we have obtained at the junction is generated when the spin texture moves temporally.
On the basis of these results, in Section 2.3 we discussed a way of detecting the charge current and spin current induced on the surface of a TI attached to an MI by using ferromagnetic resonance. We assume that the dynamics of localized spin is triggered by the applied static and AC magnetic field of the microwave. The dynamics of the localized spin induced both the charge current and the spin current. Such induced currents are related to spin density, and spin polarization acts on the localized spin in much the same way as spin torque—see Eq. (14). Hence, the half‐width of ferromagnetic resonance changes as shown in Figure 3.
Spin pumping in a Weyl–Dirac semimetal hosting massless Weyl–Dirac fermions is summarized in section 3. The results are obtained within the same formalism as laid out in Sections 2.1 and 2.2. The charge and spin current as well as the charge and spin density are given in Table 2. Semimetals are subject to spin–momentum locking. The spin direction of Weyl–Dirac fermions brought about by spin–momentum locking is perfectly parallel/antiparallel to its momentum and its locking is determined by the helicity degrees of freedom of Weyl–Dirac fermions. As a result, the charge current and spin polarization induced depend on the helicity indices. Eqs. (38) and (39) show that localized spin dynamics induces the charge current and spin polarization, respectively; hence, localized spin plays the role of an effective electric field Es,±WDS acting on electrons . Moreover, Es,±WDS depends on the helicity index, whereas total charge current is proportional to μ5. As a result, a nonzero charge current is generated when there is population imbalance between each helicity. On the other hand, localized spin dynamics also drives the spin current when localized spin depends on time and position in a Weyl–Dirac semimetal. The spin current is finite even in the absence of population imbalance. Then, the spin current does not accompany charge flow. These results may be of use to next‐generation spintronics devices based on Weyl–Dirac semimetals.
This work was supported by a Grant‐in‐Aid for JSPS Fellows (Grant No. 13J03141), by a Grant‐in‐Aid for Challenging Exploratory Research (Grant No. 15K13498), and by the Core Research for Evolutional Science and Technology (CREST) of Japan Science.