emergence and manipulation of spins are a major research topic in spintronics .
applying a controlled light is one of the promising techniques to manipulate spins .
recently , the spin angular momentum of a circularly polarized light has been observed to induce the magnetization in solid - state materials through spin - orbit interactions@xcite .
this technique has further been applied to the ultrafast magnetization switching , whose time is much shorter than that by an applied magnetic field@xcite .    when a light is irradiated on a surface of a three - dimensional topological insulator ( ti ) ,
spin is predicted to emerge in the perpendicular direction to the electric field of the light@xcite . here
, a ti is an anomalous material with strong spin - orbit interactions .
electrons insulate in the bulk , while they conduct on the surface of the ti , where the exotic surface state of the ti is caused by both the spin - orbit interaction and the topological electric structure@xcite . on the surface of the ti , the direction of the spin and that of the momentum are perfectly locked to be perpendicular to each other , which is dubbed spin - momentum locking . because of this spin - momentum locking ,
the charge current generated along the direction of the electric field causes the spin density in the perpendicular direction@xcite .
such a manipulation of spin and charge current using a light may make it possible to develop magneto - optical devices based on tis@xcite .
recently , magneto - optical effects and optical excitation using a twisted light beam , whose phase is twisted around the direction of the propagation of light , have been theoretically predicted@xcite and experimentally carried out@xcite .
a twisted light has the following two intriguing properties distinct from a plane wave@xcite .
first , the phase of the light is twisted around the center of the beam , and hence , has a singularity at the center . as a result ,
strength of the light becomes zero at the center of the beam .
second , because of the twisted phase , the strength of the light strongly depends on the space , whose distributions are manipulated by the angular momentum of the light .
the above properties can be well - understood by writing down the electric field of the light .
the electric field of the twisted light beam traveling along the @xmath0 axis at @xmath1 , @xmath2 , can be described by@xcite @xmath3,\end{aligned}\ ] ] where @xmath4 is the two - dimensional polar coordinates at @xmath1 , @xmath5 is the time , and @xmath6 and @xmath7 are the momentum and frequency of the twisted light beam , respectively .
here , @xmath8 denotes the magnitude of the electric field , which depends on the space and becomes zero at the center @xmath9 for a nonzero @xmath10 due to the phase singularity .
@xmath11 and @xmath12 represent the @xmath0 components of the spin and orbital angular momentum of the light , respectively .
the former corresponds to the direction of the circular polarization , i.e. , @xmath13 represents a right - handed ( left - handed ) circularly polarized wave , while the latter describes the winding of the electric field in the @xmath1 plane .
in fact , the electric field of a twisted light has the topological quantity .
we will see later that the winding number of a twisted light given by eq .
( [ eq:1 - 1 ] ) is proportional to @xmath14 and @xmath15[see the discussion below eq .
( [ eq : winding - number ] ) ] .
so far , it has been theoretically predicted that in the presence of the spin - orbit interaction , unconventional photo - induced spin excitation and current emerge due to the spatial dependence of the strength of the electric field of the twisted light .
however , the orbital angular momentum of the light is not transfered to the spin polarization since the beam waist is much larger than the width of the electron wave function , which is in the order of the lattice constant , and hence each electrons feels locally uniform electric field .
actually the experimental investigation of the photo - induced spin polarization@xcite could not detect the orbital angular momentum dependence in the semiconductor with the rashba and dresselhaus type spin orbit interaction as an exception there is a theoretical prediction for cylindrical quantum disks@xcite .    in this paper
, we theoretically study spin and charge generation due to the electric field of the twisted light beam on a disordered surface of a doped ti by using the green s function technique .
we analytically calculate the linear response function of the spin density to a space - time dependent external electric field .
we find that the local and nonlocal spin densities are induced by the electric field and the gradient of the electric field , respectively , via the spin - momentum locking .
here , the local spin density comes from the charge current that flows along the electric field , whereas the nonlocal one couples to the diffusive charge current due to the impurity scatterings on the disordered surface of the ti .
in addition , the gradient of the electric field also induces the charge density and the spin current . applying the obtained results to the electric field of a twisted light beam
, we find that various spin distributions appear depending on the orbital as well as spin angular momentum of the light .
moreover the spin distributions have topological structures i.e. , magnetic vortex - like textures , characterized with winding numbers , which dependent on both @xmath16 and @xmath10 .
the induced spin structure evolves in time but its winding number remains a constant .
the manipulation of such a topological spin structure could be applicable for the spintronics related to magnetic vortices and skyrmions .
this paper is organized as follows . in sec .
[ sec : ii ] , we introduce the model hamiltonian for the disordered surface of the ti in the presence of a space - time dependent electromagnetic field .
we also present the green s functions on the disordered surface of the ti . in sec .
[ sec : iii ] , we calculate the induced spin density on the surface of the ti within the linear response to the applied electric field .
section iv discusses the general properties of the charge density , spin density , charge current , and spin current induced by the electric field on the surface of the ti .
section [ sec:4 ] discusses the properties of the twisted - light - induced spin and charge distributions .
section [ sec:5 ] summarizes the paper .
appendices [ sec : a]-[sec : g ] give the detailed calculations used in sec . iii .
we consider a disordered surface of a three - dimensional ti and apply a space - time dependent electromagnetic field to it .
we assume that impurities on the disordered surface are nonmagnetic .
the hamiltonian we consider is given by @xmath17\psi ,       \\ \label{eq:2 - 3 } % \ref{eq:2 - 3 } \mathcal{h}_{\textrm{em } } & = - ev_{\textrm{f } } \int d\bm{x } \psi^\dagger   ( \hat { \bm{\sigma } } \times \bm{a}_{\textrm{em}})_z \psi ,      \\ \label{eq:2 - 4 } % \ref{eq:2 - 4 } \mathcal{v}_{\textrm{imp } } & =   \int d\bm{x } \
u_\textrm{i } \psi^\dagger \psi ,    \end{aligned}\ ] ] where @xmath18 and @xmath19 are the creation and annihilation operators of conduction electrons on the surface of the ti , @xmath20 are the pauli matrices , and @xmath21 is the elementary charge of electrons . here
, we assume a doped ti , and @xmath22 and @xmath23 are the fermi energy and the fermi velocity , respectively , on the surface of the doped ti .
@xmath24 is the gauge coupling between conduction electrons and the electromagnetic field .
the vector potential of the electromagnetic field @xmath25 generally depends on the space and time , and the electric field and the magnetic field are respectively given by @xmath26 and @xmath27 .
@xmath28 describes the potential due to the nonmagnetic impurity scatterings , where @xmath29 is the potential energy density with @xmath30 being the number of the impurities , @xmath31 a constant , @xmath32 the position of the @xmath33-th impurity on the surface , and @xmath34 the area of the surface . here
, the contribution from the impurity potential is treated as the impurity average @xmath35 , where @xmath36 is the fourier transform of @xmath37 .    to calculate the spin density and the charge density
, we use the green s function method . in the absence of the electromagnetic field ,
the retarded green s function is given by @xcite @xmath38^{-1}.\end{aligned}\ ] ] here , a variable with a hat denotes a two - by - two matrix . by calculating the self - energy @xmath39 within the self - consistent born approximation @xcite and expanding it with respect @xmath40 up to the linear terms , @xcite eq .
( [ eq:2 - 5 ] ) is rewritten as @xmath41^{-1},\end{aligned}\ ] ] where @xmath42 is the modified fermi velocity due to nonmagnetic impurity scatterings with @xmath43 being a constant depending on the properties of the ti , and the imaginary part of the self - energy @xmath44 defines the transport relaxation time @xmath45 . here , @xmath46 is the concentration of the impurities on the surface and @xmath47 is the density of states of electrons on the surface . since we are considering a metallic state , @xmath48 satisfies @xmath49 . by comparing eqs .
( [ eq:2 - 2 ] ) and ( [ eq:2 - 6 ] ) , we see that the effective hamiltonian for the surface electrons affected by impurities is given by the right - hand side of eq .
( [ eq:2 - 2 ] ) with replacing @xmath50 with @xmath51 .
accordingly , @xmath50 in eq .
( [ eq:2 - 3 ] ) is replaced by @xmath51 .
[ sec:3 ] we calculate the spin density induced by an applied electric field on a disordered surface of a doped ti by using the keldysh green s function method within the linear response to @xmath52 . the spin density @xmath53 is described by using the lesser component of the keldysh - green s function in the same position and time @xmath54 as @xmath55 \hspace{5mm}(i = x , y , z),\end{aligned}\ ] ] where @xmath56 denotes the trace over the spin indices .
then , from the dyson s equation for @xmath57 the induced spin density within the linear response to @xmath52 is given by @xmath58 a_{{\textrm{em } } , u } ( { \bm{q},\omega } ) , \end{aligned}\ ] ] where @xmath59 is the levi - civita symbol , @xmath60 is the spin - spin response function , and @xmath61 and @xmath7 are the momentum and frequency of @xmath62 , respectively .
the response function @xmath60 can be decomposed as @xmath63 and @xmath64 is represented by @xmath65^ <      & = \sum_{\bm{k } , \omega } [ \hat{g}_{\bm{k}-\frac{q}{2},\omega-\frac{\omega}{2 } } \hat{\lambda}_\nu ( \omega,\bm{q},\omega ) \hat{g}_{\bm{k}+\frac{q}{2},\omega+\frac{\omega}{2 } } ] ^{<}.\end{aligned}\ ] ] here , @xmath66 is the vertex function due to @xmath28 , whose diagram is shown in the fig .
[ fig : diagram ] , and is given by    . the dashed and wavy lines mean the potential due to the nonmagnetic impurity scatterings and gauge coupling between conduction electrons and the electromagnetic field , respectively . ]    @xmath67^n}_{\nu \mu}\hat{\sigma}_\mu,\end{aligned}\ ] ]    where @xmath68 is the two - by - two identity matrix and @xmath69 is a @xmath70 matrix defined from the following equation @xmath71    expanding eq .
( [ eq:3 - 3 ] ) with respect to the retarded and advanced green s functions , @xmath72 and @xmath73 , based on the formula@xcite @xmath74 , where @xmath75 is the fermi distribution function , the spin - spin response function can be divided into three terms : @xmath76 here , @xmath77 ( = ) is defined by @xmath78^n}_{\nu \mu}\hat{\sigma}_{\mu } ,      \\ \label{eq:3 - 5 - 2 } % \ref{eq:3 - 5 - 2 }      \hat { \gamma }
^{\textrm{ab}}_\nu(\omega,\bm{q},\omega )          & \equiv n_{\textrm{i } } u_{\textrm{i}}^2   \sum_{\bm{k } }   \hat{g}^{\textrm{a}}_{\bm{k}-\frac{q}{2},\omega-\frac{\omega}{2 } } \hat{\sigma}_\nu \hat{g}^{\textrm{b}}_{\bm{k}+\frac{q}{2},\omega+\frac{\omega}{2 } } ,          \\          & = \sum_{\mu=0,x , y , z}\gamma^{\textrm{ab}}_{\nu \mu}(\omega,\bm{q},\omega)\hat{\sigma}_\mu.\end{aligned}\ ] ] by expanding @xmath79 and @xmath80 with respect to @xmath81 and @xmath7 , we find that they are in the order of @xmath82 and @xmath83 ( @xmath84 ) in @xmath85 ( @xmath86 ) can be approximated with @xmath87 ( see appendices [ sec : a]-[sec : a - c ] for the detailed calculation ) .
then , by expanding the fermi distribution function with respect to @xmath7 , the dominant contributions of eqs .
( [ eq:3 - 7 ] ) , ( [ eq:3 - 8 ] ) and ( [ eq:3 - 9 ] ) are written by @xmath88 in addition , @xmath89 and @xmath90 are shown to be much smaller than @xmath91 ( the detailed calculation is given in appendix [ sec : d ] ) .
thus , @xmath92 is approximately given by @xmath93 in the low - temperature limit , we approximate the derivative of the fermi distribution function as @xmath94 .
then , the integral over @xmath95 in @xmath96 reduces to @xmath97 we further expand @xmath98 as @xmath99 , and rewrite eq .
( [ eq:3 - 13 ] ) as @xmath100 here , we define @xmath101 , which is calculated up to the quadratic terms in @xmath102 and the primary terms in @xmath7 as ( see appendices [ sec : a ] and [ sec : b ] for the detailed derivation ) @xmath103 ,              \\
\label{eq:3 - 16 } % \ref{eq:3 - 16 } \hat{i}_{\zeta = x , y }      & = \frac{\pi\nu_e}{2\eta } \biggl [ \sum_{\nu = x , y } \biggl\ { \frac{1}{2}\biggl(1- i\omega \tau - \frac{3}{4}\ell^2 q^2 \biggr)\delta_{\zeta \nu } + \frac{1}{4}\ell^2 q_\zeta q_\nu   \biggr\ } \hat{\sigma}_\nu        + \sum_{\alpha = x , y } \frac{i}{2}\ell \hat{\sigma}_0 q_\alpha \epsilon_{\zeta \alpha z } \biggr ] ,       \\
\label{eq:3 - 16 - 2 } % \ref{eq:3 - 16 - 2 } \hat{i}_{\zeta = z }      & = o \left ( \frac{\hbar}{\epsilon_\textrm{f } \tau } \right ) , \end{aligned}\ ] ] where @xmath104 is the mean free path .
since @xmath105 is negligibly small as compared with @xmath106 and @xmath107 , we consider only the contributions from @xmath108 .
since eqs .
( [ eq:3 - 15 ] ) and ( [ eq:3 - 16 ] ) do not include @xmath109 , they are represented by using the pauli matrices as @xmath110 where @xmath111 is the 3 @xmath112 3 symmetric matrix given by    @xmath113    on the other hand , from eq .
( [ eq:3 - 4 - 2 ] ) , the vertex function @xmath114 can be described by @xmath115_{\nu \alpha}\hat{\sigma}_\alpha + \cdots       \\      & = \sum_{\alpha=0,x , y}[(1-\gamma^{\textrm{ra}})^{-1}]_{\nu\alpha } \hat{\sigma}_\alpha,\end{aligned}\ ] ] where the second equality holds when @xmath116 is satisfied . by using @xmath117 , @xmath118 , and eq .
( [ eq:3 - 18 ] ) , one can see @xmath119 indeed satisfies @xmath116 .
then , the matrix @xmath120 is calculated as @xmath121 from which we obtain the spin - spin response function as @xmath122_{\zeta ' \nu } \hat{\sigma}_{\zeta'},\end{aligned}\ ] ] where we have used the fact that @xmath123 and @xmath124 are symmetric matrices .
thus , from eqs .
( [ eq:3 - 2 ] ) and ( [ eq:3 - 20 ] ) , the @xmath125 components of the spin density are given by @xmath126 } %    \textrm{tr } [ \hat{\pi}_{\mu\nu } ] a_{\textrm{em},u }      & = - \frac{e\tilde{v}_\textrm{f}\nu_e   \tau}{2l^2 } \epsilon_{z \nu u } \partial_t\sum_{\bm{q},\omega } e^{i(\omega t -\bm{q}\cdot\bm{x } ) }      [ \gamma^{\textrm{ra}}\lambda^{\textrm{ra}}]_{\mu\nu }   a_{{\textrm{em } } , u}.\end{aligned}\ ] ] substituting eq .
( [ eq:3 - 19 ] ) in eq .
( [ eq:3 - 20 - 2 ] ) , we obtain @xmath127 }      % \ { [ \hat{\gamma}\hat{\lambda}^{\textrm{ra}}]_{xx }   a_{{\textrm{em } } , y }      %    -[\hat{\gamma}\hat{\lambda}^{\textrm{ra}}]_{xy }   a_{{\textrm{em } } , x } \ }      % \\
\notag       & = \frac{1}{2}e\tilde{v}_\textrm{f}\nu_e   \tau e_y        + \frac{e\tilde{v}_\textrm{f}\nu_e   \tau}{2l^2 }   \partial_t \sum_{\bm{q},\omega } e^{i(\omega t -\bm{q}\cdot\bm{x } ) }      \frac {   \ell^2 ( q_y^2   a_{{\textrm{em } }
, y } + q_y q_x   a_{{\textrm{em } } , x})}{q^2 \ell^2 + i\omega \tau } ,      \\ \label{eq:3 - 22 } % \ref{eq:3 - 22 } s_y & = -\frac{1}{2}e\tilde{v}_\textrm{f}\nu_e   \tau e_x       - \frac{e\tilde{v}_\textrm{f}\nu_e   \tau}{2l^2 }   \partial_t \sum_{\bm{q},\omega } e^{i(\omega t -\bm{q}\cdot\bm{x } ) }      \frac {   \ell^2 ( q_x^2   a_{{\textrm{em } } , x } + q_y q_x   a_{{\textrm{em } } , y})}{q^2 \ell^2 + i\omega \tau}.\end{aligned}\ ] ]    the second terms of eqs .
( [ eq:3 - 21 ] ) and ( [ eq:3 - 22 ] ) can be described by using the charge density @xmath128 on the surface . here , @xmath129
a_{{\textrm{em } } , u } ( \bm{q},\omega)$ ] , is obtained from the charge - spin response function @xmath130 as @xmath131_{0\nu }   a_{{\textrm{em } } , u } \biggr\ }      \\       & =      \frac { e^2 \tilde{v}_{\textrm{f}}\nu_e \tau}{l^2 } \ell \partial_t \nabla_\nu \sum_{\bm{q},\omega}e^{i(\omega t-\bm{q}\cdot\bm{x } ) } \frac { 1 } { q^2 \ell^2 + i\omega \tau }     a_{{\textrm{em } } , \nu }      \label{eq:3 - 23 } % \ref{eq:3 - 23 }      \\      & =   - 2 e^2 \nu_e d \tau   \bm{\nabla } \cdot \langle \bm{e } \rangle_{\textrm{d } } ,       \label{eq:3 - 24 } % \ref{eq:3 - 24}\end{aligned}\ ] ] where @xmath132 is the diffusion constant . here
, @xmath133 is defined by the convolution of @xmath134 and the diffusive propagation function @xmath135 as @xmath136.\end{aligned}\ ] ] the diffusive propagation function @xmath135 is also the green s function satisfying the following differential equation @xmath137 equations ( [ eq:3 - 24 ] ) and ( [ eq:3 - 25 ] ) show that due to the impurities the effect of the applied electric field on the surface electrons is not instantaneous but diffusively propagates on the surface of the ti . equation ( [ eq:3 - 25 ] ) gives the definition of such a nonlocal electric field .
suppose that the surface of the ti is isotropic , the gradient of the charge density is given by @xmath138 which is related to the spin density [ eqs .
( [ eq:3 - 21 ] ) and ( [ eq:3 - 22 ] ) ] as @xmath139 on the other hand , the spin density on the surface of the ti is related to the charge current via @xmath140 where we have replaced @xmath23 in eq .
( [ eq:2 - 3 ] ) with @xmath141 so as to take into account the effects of impurities . by substituting eq .
( [ eq:3 - 31 ] ) in eq .
( [ eq:3 - 33 ] ) , we obtain @xmath142 we have confirmed that eqs .
( [ eq:3 - 24 ] ) and ( [ eq:3 - 32 ] ) satisfy the charge conservation law : @xmath143 ( see appendix [ sec : f ] for the detailed calculation ) .
using the results in sec .
[ sec:3 ] , we discuss the property of the charge , spin , charge current , and spin current densities induced by the electric field applied on the disordered surface of the doped ti .      we find that as shown in eq .
( [ eq:3 - 24 ] ) the charge density @xmath128 is induced by the divergence of the nonlocal electric field:@xmath144 .
therefore , when we apply a uniform electric field , no charge density is induced . from eqs .
( [ eq:3 - 25])-([eq:3 - 27 ] ) we obtain the diffusion equation for the charge transport : @xmath145 which indicates that the divergence of the applied electric field works as a source of the diffusive propagation of the charge density .
we find that from the left side of the equation above , @xmath146 , the diffusion constant is @xmath147 , a twice of that on the surface of a metal@xcite .
the factor @xmath148 comes from the difference in the self - energy due to impurity scattering : the self - energy on the surface of an isotropic metal is given by @xmath149 , where @xmath150 is the density of states in the metal , whereas that on the surface of the ti is @xmath151 ; the factor @xmath152 , which originates from the linear dispersion of the surface of the ti , leads to the coefficient @xmath147 . here ,
the diffusive equation of motion qualitatively agrees with the previous works@xcite .
we turn to the discussion on the spin density given by eq .
( [ eq:3 - 31 ] ) .
we find that the spin density can be divided into the one induced by the local electric field @xmath134 and that by the nonlocal electric field @xmath133 .
we define the local spin density @xmath153 as the first term of eq .
( [ eq:3 - 31 ] ) : @xmath154 the local spin density lies in the @xmath155 plane and is perpendicular to the local electric field @xmath134@xcite , which is the consequence of the spin - momentum locking on the surface of the ti .
the spin due to the local electric field is a kind of the edelstein effect@xcite . on the other hand , the nonlocal spin density @xmath156 , the second term of eq .
( [ eq:3 - 31 ] ) , is generated from the spatial gradient of the charge density .
the direction of @xmath157 is in the @xmath155 plane and perpendicular to the gradient of the charge density @xmath158 . here
, @xmath157 is also written in terms of @xmath159 as @xmath160 ,   %    \\ \label{eq : spin - y10 } % \ref{eq : spin - y10 } % s^{(\textrm{nl})}_y %    & =   - \frac{e\tilde{v}_\textrm{f } \nu_e \ell^2 } { 2 } [ \nabla_x^2   \langle e_x \rangle_\textrm{d}+ \nabla_x \nabla_y \langle e_y % \rangle_\textrm{d } ] .
\label{eq : spin - x10 } % \ref{eq : spin - x10 } \bm{s}^{(\textrm{nl } ) }       & =   -\frac{e\tilde{v}_\textrm{f } \nu_e\tau \ell^2 } { 2 } \left ( \bm{z } \times \bm{\nabla } \right)(\bm{\nabla } \cdot \langle \bm{e } \rangle_{\textrm{d}}),\end{aligned}\ ] ] which means that the spin density is generated by the second spatial derivative of the nonlocal electric field @xmath159 .
the nonlocal spin density diffusively propagates through the impurity scatterings on the surface of the ti .
the diffusion propagation of the spin is described by @xmath161 we find that the diffusion propagation of the spin is triggered by an inhomogeneous electric field , @xmath162 .
( this property is also predicted in ref .
@xmath163 ) . hence ,
when we apply a uniform electric field on the surface , the nonlocal spin density is not generated .
as in the case of charge density , the diffusion constant for the spin density is @xmath147 .
we note that both the local and nonlocal spin densities are proportional to @xmath51 . since @xmath51 s on the top and
the bottom sides of the ti have opposite signs , the direction of the induced spin on the top surface of the ti is perfectly opposite to that on the bottom surface of the ti , when we apply the same electric field on both the top and bottom side of the ti .
the charge current on the surface of the ti is proportional to the spin density as shown in eq . ( [ eq:3 - 33])@xcite .
the origin lies on the spin - momentum locking on the surface of the ti . from eq .
( [ eq:3 - 32 ] ) , we find that the charge current @xmath164 can also be divided into the local and nonlocal parts .
the first term of eq .
( [ eq:3 - 32 ] ) gives the local charge current @xmath165 .
@xmath166 is the electric current directly induced by the applied electric field .
we find from the above result that the longitudinal conductivity is given by @xmath167 , which agrees with the existing work@xcite .
the second term of eq .
( [ eq:3 - 32 ] ) corresponds to the nonlocal charge current density @xmath168 .
@xmath169 is the diffusion current and is generated by the spatial gradient of the charge density on the disordered surface of the doped ti .
the charge current can be rewritten by using the nonlocal electric field @xmath170 as @xmath171 since both the local and nonlocal charge currents are proportional to @xmath172 , their directions are the same for top and bottom surface of the ti .
next , we calculate the spin current due to the applied electric field on the disordered surface of the doped ti .
the spin current @xmath173 is defined by @xmath174 where the subscript and superscript of @xmath173 denote the direction of the flow and spin , respectively , and @xmath175 represents the spin relaxation . using the hamiltonian in eq .
( [ eq:2 - 1 ] ) and the heisenberg equation for @xmath176 , we obtain @xmath177 note that the direction of the flow and spin is perpendicular each to other . this is the consequence of the spin - momentum locking on the surface of the ti .
we also note that the spin current is proportional to the charge density , and from eq .
( [ eq:3 - 24 ] ) , proportional to the divergence of the nonlocal electric field : @xmath178 hence , when we apply a spatially uniform electric field on the surface , no spin current is induced . besides , we find that the spin current is an odd - function of @xmath51 , which means the spin current depends on the chirality on the surface of the ti .
namely , the relative direction between flow and spin of @xmath173 on the top side of the ti is opposite to that on the bottom side of the ti .    in the conventional spin - orbit coupled systems ,
the spin current is generated by an applied electric field@xcite , which called the spin hall effect .
besides , the generated spin current can be converted into the charge current via the spin - orbit interaction.@xcite these effect can be understand from the coupling between the spin current and the charge current : @xmath179@xcite . on the surface of the ti , on the other hand , we find from eqs .
( [ eq:4 - 3 - 2 ] ) and ( [ eq:4 - 4 - 3 ] ) that the nonlocal charge current is proportional to the gradient of the spin current@xcite : @xmath180 again , this is the consequence of the spin - momentum locking and eq .
( [ eq:4 - 4 - 4 ] ) generally holds for the system of electrons on a surface of tis@xcite .
this property in the ti is distinct from that in a conventional metal .
the direction of the charge current is parallel to the spatial gradient of @xmath173 .
this relation is plausible due to the following reasons .
first , the charge density @xmath128 is proportional to the spin current .
second , a diffusive particle current generally proportional to a spatial gradient of particles .
we note that there is no relation between the spin current and the local charge current @xmath166 .    finally , we comment on the property of the spin relaxation torque .
the relaxation torque @xmath175 defined in eq .
( [ eq:4 - 4 - 1 ] ) is obtained within the linear response to the electric field as @xmath181            + o [ ( \bm{z } \times \bm{\nabla})_\alpha ( \bm{\nabla } \cdot \langle \dot{\bm{e } } \rangle_{\textrm{d } } ) ] .\end{aligned}\ ] ] here , @xmath182 can be divided into the local and nonlocal terms , which correspond to the first and second terms , respectively , in the first square bracket in the most right - hand side of eq .
( [ eq:4 - 4 - 5 ] ) .
the local one is given by the time derivative of the applied electric field , and its direction is perpendicular to both @xmath183 and @xmath184 .
the nonlocal one is proportional to the second derivative of the nonlocal electric field @xmath185 .
these above results and properties are the same as the spin density on the surface of the ti with magnetism@xcite .
using the results obtained in sec .
[ sec:4 ] , we discuss the properties of the spin and charge densities due to the electric field of a twisted light beam with various orbital angular momentum .
first , we explain the property of the electric field of the twisted light beam with the laguerre - gaussian modes@xcite .
we assume that the twisted light beam propagates along the @xmath0 axis and the electric field of the twisted light beam lies in the @xmath155 plane at the top surface of the ti ( @xmath1 ) .
the schematic of the system is illustrated in fig .
[ fig : fig1 ] .
the twisted light beam satisfies the wave equation , @xmath186 , where @xmath187 is the velocity of light in a vacuum .
then , the electric field @xmath188 on the surface is written by@xcite @xmath189 where @xmath190 is the distance from the center of the light on the top surface ( @xmath1 ) and @xmath191 is the azimuthal angle .
the helicity @xmath192 denotes the right - hand ( left - hand ) circularly polarized light and corresponds to the spin angular momentum of light @xmath193 .
the orbital angular momentum of light , @xmath194 , determines the whirling pattern of the electric field on the plane at @xmath1 , which can be manipulated in experiments@xcite .
the phase @xmath195 depends on @xmath196 and the distance from the light source , @xmath0 , and is given by@xcite @xmath197 }   - \frac{q_z r^2 } { 2r(z_0)},\end{aligned}\ ] ] where the first term denotes guoy phase , @xmath198 $ ] is the radius of the beam curvature with taking the origin of the @xmath0 axis at the beam waist , @xmath199 is the rayleigh range , @xmath200 is the waist size of the @xmath201 mode , @xmath202 is the wavelength , and @xmath6 is the wave vector of the light .
the integer @xmath203 is one of the indices that specify the laguerre - gaussian mode and denotes the number of oscillations of the electric field @xmath204 in the radial direction . here
, @xmath205 is given by @xmath206 } } ( \sqrt{2}u)^{|m^z_{\textrm{l}}| } l_{|m^z_{\textrm{l}}|}^{p } ( 2u^2 ) \exp { ( - u^2)},\end{aligned}\ ] ] where @xmath207^{\frac{1}{2 } } ] $ ] , @xmath208 is a constant , and @xmath209 is the laguerre polynomials defined by @xmath210 in the following discussion , we consider only the @xmath211 modes .
we also assume that the twisted light beam is focused at the surface of the ti , i.e. , @xmath212 .
then , the phase @xmath213 becomes zero , and the magnitude @xmath204 defined in eq .
( [ eq : mathcal(e ) ] ) reduces to @xmath214    figures [ fig : fig2](a)-[fig : fig2](d ) show the snapshots of the electric field for @xmath215 and @xmath216 and @xmath217 . in both cases of @xmath218 and @xmath219 ,
the amplitude of the electric field exponentially decays with @xmath220 .
in addition to this , for the cases of nonzero @xmath15 , the magnitude of the electric field vanishes at @xmath9 because of the phase singularity .
this is a characteristic property of the twisted light beam . besides , the direction of the electric field depends on the polar angle around the center of the incident light : while the direction of the electric field is uniform for @xmath221 [ fig .
[ fig : fig2](a ) ] , the direction of the electric field at @xmath222 rotates by @xmath223 in the counter - clockwise direction as one goes around the beam center from @xmath224 to @xmath223 [ fig .
[ fig : fig2](b ) ] . similarly for the cases of @xmath225 and @xmath226 , the direction of the electric field changes by @xmath227 and @xmath228 , respectively [ figs .
[ fig : fig2](c ) and [ fig : fig2](d ) ] .
note that the configurations shown in figs .
3(a)3(d ) are snapshots and they evolve in time depending on @xmath14 : @xmath215 means that the electric field at a fixed point rotates in the clockwise direction as time evolves [ fig .
[ fig : fig2](e ) ]
.    the topological properties of the twisted light beam discussed above can be understood by introducing the winding number of the electric field . in general , the winding number of a 2d vector field @xmath229 on a closed loop @xmath230 is defined by @xmath231 \equiv \frac{1}{2\pi}\oint_{c } \delta_{ij}\epsilon^{\mu\nu}dx_i \frac{n_\mu}{|\bm{n}| } \frac{\partial}{\partial x_j } \left(\frac{n_\nu}{|\bm{n}|}\right),\end{aligned}\ ] ] where @xmath232 is the 2d levi - civita symbol , and @xmath233 is supposed to be nonzero on @xmath230 .
the winding number corresponds to the number of times the 2d unit vector @xmath234 rotates about the @xmath0 axis as one traces the contour @xmath230 . for the electric field given by eq .
( [ eq : ef - of - opticalvotex ] ) , the winding number @xmath235 is defined on a contour that encloses @xmath9 . substituting eq .
( [ eq : ef - of - opticalvotex ] ) in eq .
( [ eq : winding - number ] ) , we obtain @xmath236 .
for example , for the electric field with @xmath237 , and @xmath226 , we have @xmath238 , and @xmath217 , respectively , which are consistent with the configurations shown in fig .
[ fig : fig2 ] .
the result @xmath239 is also consistent with the fact that the direction of the electric field with @xmath240 is spatially uniform even for @xmath219 and that the whirling direction for @xmath241 is opposite to that for @xmath215 .    ,
@xmath242 , @xmath243 , and @xmath244 , where @xmath14 and @xmath15 denote the spin and orbital angular momentum of light , respectively . shown
are density plots of the magnitude of the electric field , and the black arrows show the direction of the electric field .
@xmath200 is the waist size for the @xmath245 mode .
( e ) time evolution of the electric field with @xmath246 ( top ) and @xmath226 ( bottom ) , where @xmath247 with @xmath7 being the angular frequency of the light . ]
we consider the charge density due to the electric field of the twisted light beam on the disordered surface of the ti .
the setup we consider is schematically described in fig .
[ fig : fig1 ] . the induced charge density @xmath128 is given by eq . ( [ eq:3 - 24 ] ) , which can be rewritten as @xmath248 by using eq .
( [ eq : ef - of - opticalvotex ] ) , @xmath249 and @xmath250 the divergence of the electric field for @xmath251 is given by @xmath252}.      \label{eq : nabla - e}%\ref{eq : nabla - e}\end{aligned}\ ] ] then , @xmath253 becomes @xmath254 where @xmath255 denotes the total angular momentum of light .
we further simplify eqs .
( [ eq : rho_e_laguerre ] ) and ( [ eq : rho_e ] ) .
first of all , these equations are valid only for @xmath256 .
this condition means that the period of the oscillation of @xmath253 , which is the same as that of the electric field , @xmath247 , is much slower than the electron relaxation time @xmath48 . on the other hand ,
the length scale of the spatial variation of @xmath253 is in the order of the beam waist @xmath200 [ see eq . ( [ eq : mathcal(e ) ] ) and fig . [
fig : fig2 ] ] , which is comparable to the wavelength @xmath202 of the light .
since @xmath202 satisfies @xmath257 , the spatial variation of @xmath253 is much slower than the mean - free path @xmath258 , where @xmath187 is the speed of light and we have used @xmath259 for realistic tis@xcite .
then , since @xmath48 and @xmath258 determines the decay time and decay length of the diffusion propagator @xmath260 , respectively , @xmath261 in the integrand of eq .
( [ eq : rho_e_laguerre ] ) can be approximated as @xmath262 , and the convolution can be approximately described by @xmath263 where @xmath264 is a constant coefficient and is estimated by @xmath265    , @xmath266 , @xmath243 , and @xmath267 .
the dashed lines in ( a ) , ( b ) , and ( c ) indicate axes of the inversion symmetry .
( e ) time evolution of the charge density with @xmath246 ( top ) and @xmath226 ( bottom ) . in all figures ,
we use @xmath268 mm , @xmath269 mev , @xmath270 m / s , and @xmath271 s. ]    figure [ fig : fig3 ] represents the nonlocal charge density due to the twisted light beam , @xmath128 , for a @xmath272 and @xmath273 and @xmath217 .
we find that the distribution of @xmath128 depends on the @xmath0-component of the total angular momentum of light , @xmath274 . for @xmath275 [ fig . [
fig : fig3](b ) ] , the charge density is isotropically induced from the center and the sign of the induced charge changes at @xmath276 . on the other hand ,
the charge density for nonzero @xmath277 distributes anisotropically .
the symmetry of the distribution of @xmath128 with @xmath278 [ figs .
[ fig : fig3](a ) and [ fig : fig3](c ) ] and @xmath279 [ fig .
[ fig : fig3](d ) ] are the same as that for the electron wave functions with the @xmath280 and @xmath281 orbital , respectively . the dashed lines in figs .
[ fig : fig3](a ) , [ fig : fig3](c ) and [ fig : fig3](d ) indicate the axes of the inversion symmetry of the charge density .
as time evolves , the distribution of the induced charge density rotates around the beam center ( @xmath282 ) or oscillates around the zero value @xmath283 with the frequency of light @xmath7 [ see fig . [
fig : fig3](e ) ] .
the time evolution of the charge density also depends on the total angular momentum of light : when the sign of the total angular momentum is minus ( plus ) , the distribution of the charge density rotates in the clockwise ( counterclockwise ) direction around the phase singularity during the irradiation .
when we turn off the incident light , the charge density diffusively propagates on the disordered surface of the ti with obeying the diffusive equation of motion represented in eq . (
[ eq:3 - 3 - 1 ] ) .
eventually , the induced charge vanishes .
we turn to discuss the spin density induced by the electric field of the twisted light beam in the same setup as that considered in sec .
[ sec:5 - 2 ] . as discussed in sec .
[ sec:4 - 2 ] , the induced spin density can be divided into the local and nonlocal ones as @xmath284 . the local spin density @xmath153 is described in eq .
( [ eq : local - spin ] ) and its snapshots for several @xmath15 are shown in figs .
[ fig : fig4](a1)[fig : fig4](d1 ) .
the direction of @xmath153 is perfectly perpendicular to the electric field .
we find that the dynamical vortex - like spin structure is generated by the twisted light beam and the winding number of the local spin density , which is defined by eq .
( [ eq : winding - number ] ) with @xmath285 , is identical to that of the electric field , @xmath286=\omega_v[\bm{e}]=-\sigma_{\textrm{l}}^z m_{\textrm{l}}^z\end{aligned}\ ] ] [ see figs .
[ fig : fig4](a1)[fig : fig4](d1 ) and figs .
[ fig : fig2](a)[fig : fig2](d ) ] . on the other hand ,
the nonlocal spin density @xmath157 is proportional to the spatial gradient of the charge density [ see the second term of eq .
( [ eq:3 - 31 ] ) ] , and can be estimated by using eq .
( [ eq:5 - 2 - 3 ] ) as @xmath287 .
the snapshots of @xmath157 are shown in figs .
[ fig : fig4](a2)[fig : fig4](d2 ) .
we find that dynamical vortex - like spin structures appear and the spin density becomes zero at the center of the vortex . here
, we note that @xmath157 and @xmath288 share the same winding number as they are perfectly perpendicular to each other .
since the winding number @xmath289 ) is @xmath290 around the maxima and minima ( the saddle points ) of @xmath128 , the centers of the spin vortices locate at the extrema ( minima , maxima , and saddle points ) of @xmath128 , and therefore , they align on the symmetry axis of the distribution of @xmath128 [ the dotted lines in fig .
[ fig : fig3 ] ] . for the cases shown in figs .
[ fig : fig4](a2)[fig : fig4](d2 ) , all spin vortices have the winding number @xmath291 except for the one at the center of fig .
[ fig : fig4](d2 ) , which corresponds to the saddle point of @xmath128 and has the winding number @xmath217 .
since @xmath157 is related to @xmath292 rather than @xmath134 , the configuration of the spin vortices depends on the total angular momentum @xmath293 .
note however that with the parameters for a realistic system , @xmath294 is in the order of @xmath295 and @xmath296 is negligibly small as compared with @xmath297 .
when we turn off the beam , @xmath296 becomes prominent and diffusively propagates .
we expect that the photo - induced spin texture can be observed by pump prove technique with the twisted light beam and circularly polarized light beam@xcite .
the profile of the charge current is similar to that of the spin because of the spin momentum locking on the surface of the ti . in fig .
[ fig : fig5 ] , we show the snapshots of the charge current for various angular momentum of light , where left ( right ) four panels depict the local ( nonlocal ) components . reflecting the relation @xmath298 , the magnitude of @xmath299 has the same profile as that of @xmath300 , while the direction of @xmath299 is obtained by rotating @xmath301 by @xmath302 about the @xmath0 axis . as in the case of the spin density ,
the local ( nonlocal ) part of the charge current is related to @xmath303 ( @xmath304 ) and hence its configuration is mainly determined by @xmath305 ( @xmath293 ) .
figure [ fig : fig6 ] shows the light - induced spin current . as one can see from eq .
( 55 ) , the magnitude of the spin current is proportional to @xmath306 and the direction of the spin and its current perfectly perpendicular to each other .
these properties also come from the spin - momentum locking on the surface of the ti .
we find that the spatial profile of the spin current is different from that of the ( spin - polarized ) charge current , which are shown in figs .
[ fig : fig6 ] and [ fig : fig5 ] respectively .
in fact , they are related to each other via eq .
( [ eq:4 - 4 - 4 ] ) and only the nonlocal charge current couples to the spin current .
we study the charge density , the spin density , the charge current density , and the spin current density induced by a twisted light beam shined on a disordered surface of a doped ti by using the keldysh - green s function method .
we have discussed the responses of charge and spin to the space - time dependent electric field .
the obtained results are summarized in tab .
i. the effect of the electric field on the electric charge is twofold .
first , it induces a charge current along the direction of the electric field .
second , the inhomogeneity of the electric field causes a gradient of the charge density , which then leads to a diffusive charge current .
we call the former the local charge current and the latter the nonlocal charge current , based on whether the charge current depends only on the local electric field or is affected by the nonlocal one . since the spin and momentum of electrons on a surface of a ti are locked to be perpendicular to each other , the emergence of the charge current implies that the spin density is induced in the perpendicular direction to the charge current . our calculation based on the linear response theory gives the analytic description for the local and nonlocal spin densities as well as the local and nonlocal charge current densities .
we also find that the induced charge density also gives rise to a spin current , which is related to the nonlocal part of the charge current via eq .
( [ eq:4 - 4 - 4 ] ) .    by taking into account the spatial and temporal configuration of the electric field associated with a twisted light beam
, we have shown that various types of spin vortices arise . since the local spin density is perpendicular to the electric field , their winding numbers are identical and determined by the product of the spin and orbital angular momentum of the twisted light beam . in this paper , we have assumed that the time and length scales of the diffusion of electrons is much faster than those of incident light .
in such a situation , we can approximate the nonlocal electric field with the bare electric field , and the nonlocal quantities are described using the divergence of the electric field .
thus , the configurations of the nonlocal densities and currents are determined by the total angular momentum of the twisted light beam .    [ cols="<,^,^,^,^,^",options="header " , ]
the authors would like to thank a. dutt for valuable discussions .
this work was supported by grant - in - aid for young scientists ( b ) ( no .
15k17726 ) , by grants - in - aid for scientific research on innovative areas `` topological materials science '' ( no . 15h05853 and no .
15h05855 ) from the ministry of education , culture , sports , science , and technology , japan ( mext ) , by grant - in - aid for challenging exploratory research ( no .
15k13498 ) , and by the core research for evolutional science and technology ( crest ) of the japan science . k.t . acknowledges support from a grant - in - aid for japan society for the promotion of science ( jsps ) fellows .
we estimate @xmath307 in eq .
( [ eq:3 - 14 ] ) by expanding with respect to @xmath81 and @xmath7 and by using the green s function @xmath308             + \frac{1}{8 } \sum_{\xi,\xi'=x , y } q_\xi q_{\xi ' } \frac{\partial^2 \hat{g}^\textrm{r}_{\bm{k } } } { \partial k_\xi \partial k_{\xi ' } }           + o(q^3 , q\omega , \omega^2 ) , \end{aligned}\ ] ] where we use the short hand notation @xmath309 .
@xmath307 is decomposed into four terms as @xmath310 in appendix [ sec : e ] , we list useful formulas for the integral of the green s functions , which are used in the following calculations .      first , to calculate @xmath312 in eq .
( [ eq : coefficient pi-1 ] ) , we divide @xmath313 into the even and odd functions of @xmath40 : @xmath314 where @xmath315 , @xmath316^*$ ] , @xmath317 , @xmath318 , @xmath319 , and @xmath320 are defined by @xmath321 @xmath315 and @xmath322 are the even functions of @xmath40 .
@xmath323 is represented by @xmath324           \\ \label{eq : app-12 } % \ref{eq : app-12 }             & \hspace{10 mm }          + \hbar \tilde{v}_{\textrm{f } } \sum_{\ell = x , y}\left [ h k_\ell \hat{\sigma}_\mu ( \bm{z } \times \hat{\bm{\sigma}})_\ell   + h^ * k_\ell   ( \bm{z } \times \hat{\bm{\sigma}})_\ell \hat{\sigma}_\mu \right].\end{aligned}\ ] ] the first and second terms are corresponding to the even and odd functions of @xmath40 , respectively . in the following calculation
, we simply assume that the surface of the ti is isotropic as a function of @xmath40 : @xmath325 . by taking an average over the direction of @xmath40 , @xmath326 , in eq .
( [ eq : app-12 ] ) and using @xmath327 and @xmath328 , we obtain @xmath329 the second term in the above equation becomes @xmath330 thus , we obtain @xmath331 as @xmath332 the integral of eq .
( [ eq : app-15 ] ) is obtained by @xmath333 [ ( h^*)^2 - \hbar^2 \tilde{v}_{\textrm{f}}^2 k^2 ] }       %    &   = \frac{1}{2\pi } \int_0^\infty    \frac{kdk } { [ h^2 - \hbar^2 \tilde{v}_{\textrm{f}}^2 k^2 ] [ ( h^*)^2 - \hbar^2 \tilde{v}_{\textrm{f}}^2 k^2 ] }   %    \\      = \frac{\nu_e } { 2\epsilon_{\textrm{f } } } \int_0^\infty \frac{dx } { [ h^2
-x ] [ ( h^*)^2 -x ] }           \\       = & \frac{\nu_e } { 2\epsilon_{\textrm{f } } [ h^2 - ( h^*)^2 ] } \int^\infty_0 dx \biggl [ \frac{1 } { x -h^2 } - \frac{1 } { x - ( h^*)^2 }   \biggr],\end{aligned}\ ] ] where @xmath334 is the density of states on the surface . here ,
the above integral is given by @xmath335       =   \log { \biggl| \frac{x - h^2}{x-(h^*)^2 } \biggr|_{x\to0}^{x\to\infty } }       =    i \biggl [ \arg{(x -h^2 ) } - \arg{(x -(h^*)^2 ) } \biggr]_{x\to0}^{x\to\infty}.\end{aligned}\ ] ] in the above equation , we have used @xmath336 , where @xmath337 , @xmath338 . here , @xmath339 is @xmath340 thus , @xmath341 , @xmath342 } |_{x\to0}^{x\to\infty } = - \pi + o(\hbar/\epsilon_{\textrm{f}}\tau ) ) $ ] , and @xmath343 are satisfied .
we have @xmath344 and @xmath345 eq .
( [ eq : app-16 ] ) can be estimated around the fermi energy @xmath346 as @xmath347      to calculate @xmath349 in eq .
( [ eq : coefficient pi-2 ] ) , we divide @xmath350 into the even and odd functions of @xmath40 . here
, @xmath351 and @xmath352 becomes @xmath353 + 2 h^ *   \hbar \tilde{v}_{\textrm{f } } k_\ell   ( \bm{z } \times \hat{\bm{\sigma}})_\ell ,    \\
\notag % % % % % % % \hat{\mathcal{q } } \hat{\sigma}_\mu(\hat{\mathcal{q}}^\dag)^2       = & \biggl [ h \{(h^*)^2 + \hbar^2 \tilde{v}_{\textrm{f}}^2 k^2   \ } \hat{\sigma}_\mu          + 2 h^ * \hbar^2 \tilde{v}_{\textrm{f}}^2 k_\ell k_{\ell ' } ( \bm{z}\times \hat{\bm{\sigma}})_\ell   \hat{\sigma}_\mu ( \bm{z}\times \hat{\bm{\sigma}})_{\ell ' }    \biggr ]      \\           & + \biggl [ 2 |h|^2 \hbar \tilde{v}_{\textrm{f } }    k_\ell ( \bm{z}\times \hat{\bm{\sigma}})_\ell \hat{\sigma}_\mu \hat{\sigma}_\mu             + \ { ( h^*)^2 + \hbar^2 \tilde{v}_{\textrm{f}}^2 k^2 \ } \hbar \tilde{v}_{\textrm{f } } k_\ell ( \bm{z}\times \hat{\bm{\sigma}})_\ell    \hat{\sigma}_\mu               \biggr].\end{aligned}\ ] ] the first and second terms of the above equations are the even and odd functions of @xmath40
. then , we have @xmath354 from which @xmath355 $ ] is obtained by @xmath356 \hat{\sigma}_\mu -\textrm{h.c . }
\biggr\ }      \\
\notag       & = -\frac{\hbar}{2 }    \biggl [    \frac{\nu_e}{8\eta^2 }            \left ( i\pi - \frac{2\eta^2}{\epsilon_{\textrm{f}}^2 }   \right )   \left(1+i\frac{\eta}{\epsilon_{\textrm{f } } } \right )             -\textrm{c.c . } \biggr ] \hat{\sigma}_\mu + o\left(\frac{\hbar}{\epsilon_{\textrm{f}}\tau}\right )      \\ \label{eq : coefficient c-2 no2 } % \ref{eq : coefficient c-2 no2 }       &   =    \frac{-i\hbar   \pi \nu_e}{8\eta^2 }   \hat{\sigma}_\mu + o\left(\frac{\hbar}{\epsilon_{\textrm{f}}\tau}\right ) ,      \\\notag \hat{i}^{(1)}_{z }        & = -\frac{\hbar}{2 }    \biggl\ { \sum_{k }    d^\textrm{r } ( d^\textrm{a})^2           \bigl [   h(h^*)^2 + ( h - 2h^*)\hbar^2 \tilde{v}_{\textrm{f}}^2 k^2    \bigr ] \hat{\sigma}_\mu -\textrm{h.c . } \biggr\ }      \\
\notag       & = -\frac{\hbar}{2}\left[\frac{-i\pi\nu_e}{4\epsilon_\textrm{f}^2}\left(1+i\frac{\eta}{\epsilon_{\textrm{f } } } \right ) - \textrm{c.c.}\right]\hat{\sigma}_z + o\left(\frac{\hbar}{\epsilon_{\textrm{f}}\tau}\right )      \\ \label{eq : coefficient c-2 no3 } % \ref{eq : coefficient c-2 no3 }       &   =    \frac{i\hbar   \pi \nu_e}{4\epsilon_\textrm{f}^2 }   \hat{\sigma}_\mu + o\left(\frac{\hbar}{\epsilon_{\textrm{f}}\tau}\right).\end{aligned}\ ] ]      @xmath358 is given by @xmath359 + 2\hbar^2 \tilde{v}_{\textrm{f}}^2 k_\xi [ h + \hbar \tilde{v}_{\textrm{f } } \bm{k}\cdot ( \bm{z } \times \bm{\sigma } ) ] } { [ h^2 -\hbar^2 \tilde{v}_{\textrm{f}}^2 k^2 ] ^2 }          & =   \hbar \tilde{v}_{\textrm{f } } ( \bm{z}\times \hat{\bm{\sigma}})_\xi d^\textrm{r } + 2\hbar^2\tilde{v}_{\textrm{f}}^2 k_\xi \hat{\mathcal{q } } ( d^\textrm{r})^2 .
\label{eq : differential1 gr-0 } % \ref{eq : differential1 gr-0}\end{aligned}\ ] ] by using eq .
( [ eq : differential1 gr-0 ] ) , @xmath360 becomes @xmath361              \hat{\sigma}_\mu d^\textrm{a } \hat{\mathcal{q}}^\dag       \\
= &    \hbar \tilde{v}_{\textrm{f } } ( \bm{z}\times \hat{\bm{\sigma}})_\xi \hat{\sigma}_\mu   \langle   \hat{\mathcal{q}}^\dag\rangle_k |d^\textrm{r}|^2                      + 2\hbar^2\tilde{v}_{\textrm{f}}^2 \langle   k_\xi \hat{\mathcal{q } } \hat{\sigma}_\mu   \hat{\mathcal{q}}^\dag \rangle_k ( d^\textrm{r})^2 d^\textrm{a}. \end{aligned}\ ] ] here , @xmath362 and @xmath363 are given by @xmath364 @xmath365 in eq .
( [ eq : coefficient pi-3 ] ) is given by @xmath366.\end{aligned}\ ] ] here , @xmath367 can be transformed as @xmath368 as a result , @xmath365 is obtained by @xmath369 + o\left ( \frac{\hbar}{\epsilon_\textrm{f } \tau}\right )      \\
\label{eq : a29}%\ref{eq : a29 }      & = \frac { i\pi \nu_e } { 8\eta^2 }   \hbar \tilde{v}_{\textrm{f } } \epsilon_{\mu\xi z } + o\left ( \frac{\hbar}{\epsilon_\textrm{f } \tau}\right ) ,      \\
\notag \hat{i}^{(2)}_{z \xi }       & =   \frac{\hbar \tilde{v}_{\textrm{f } } i\pi\nu_e } { 8 \eta^2 } \epsilon_{\xi z u }          \biggl [    \left ( 1 + 2\frac{\eta^2 } { \epsilon_{\textrm{f}}^2 } \right ) \delta_{z u }               + \frac{\eta}{\epsilon_{\textrm{f } } }   \epsilon_{z u\nu } \hat{\sigma}_\nu             \biggr ] + o\left ( \frac{\hbar}{\epsilon_\textrm{f } \tau}\right )      \\   \label{eq : a29-z}%\ref{eq : a29-z }      & = \frac { i\pi \nu_e } { 8\epsilon_\textrm{f}\eta }   \hbar \tilde{\xi}_{\textrm{f } } + o\left ( \frac{\hbar}{\epsilon_\textrm{f } \tau}\right ) .\end{aligned}\ ] ]      @xmath371 in eq .
( [ eq : coefficient pi-4 ] ) is represented by using the partial integral as @xmath372.\end{aligned}\ ] ] in order to consider @xmath373 , we use the following equations : @xmath374   %        \\          & =    2\hbar^2 \tilde{v}_{\textrm{f}}^2               \biggl\ {              \delta_{\xi\xi ' } \hat{\mathcal{q } } ( d^\textrm{r})^2                   + 4\hbar^2 \tilde{v}_{\textrm{f}}^2 k_\xi k_{\xi ' } \hat{\mathcal{q } } ( d^\textrm{r})^3                  +   \hbar \tilde{v}_{\textrm{f } } [   k_\xi ( \bm{z } \times \hat{\bm{\sigma}})_{\xi ' }                                        + k_{\xi ' } ( \bm{z } \times \hat{\bm{\sigma}})_{\xi } ]    ( d^\textrm{r})^2              \biggr\ } ,      \\
\notag \frac{\partial^2 \hat{g}^\textrm{r}_{\bm{k}}}{\partial k_\xi \partial k_{\xi ' } } \hat{\sigma}_\mu \hat{g}^\textrm{a}_{\bm{k } }       & =    2\hbar^2 \tilde{v}_{\textrm{f}}^2               \biggl\ {              \delta_{\xi\xi ' } \hat{\mathcal{q } } ( d^\textrm{r})^2\hat{\sigma}_\mu \hat{\mathcal{q}}^\dag d^\textrm{a }                       + 4\hbar^2 \tilde{v}_{\textrm{f}}^2 k_\xi k_{\xi ' } \hat{\mathcal{q } } ( d^\textrm{r})^3\hat{\sigma}_\mu \hat{\mathcal{q}}^\dag d^\textrm{a }          \\&\hspace{20 mm }                   +   \hbar \tilde{v}_{\textrm{f } } [   k_\xi ( \bm{z } \times \hat{\bm{\sigma}})_{\xi ' }                                        + k_{\xi ' } ( \bm{z } \times \hat{\bm{\sigma}})_{\xi } ]    ( d^\textrm{r})^2\hat{\sigma}_\mu \hat{\mathcal{q}}^\dag d^\textrm{a }              \biggr\}.\end{aligned}\ ] ] the average of @xmath40 in @xmath323 , @xmath375 and @xmath376 become @xmath377\end{aligned}\ ] ] here , we have used @xmath378 .
then , @xmath379 is given by @xmath380             \\
\notag                & +                     2 |h|^2 \hbar^2 \tilde{v}_{\textrm{f}}^2 k^2   \delta_{\xi \xi ' } \hat{\sigma}_\mu ( d^\textrm{r})^3 d^\textrm{a }           \\ \notag                 & +   \frac{1}{2 } \hbar^4 \tilde{v}_{\textrm{f}}^4 k^4    ( d^\textrm{r})^3 d^\textrm{a }           [   ( \bm{z } \times \hat{\bm{\sigma}})_\xi \hat{\sigma}_\mu ( \bm{z } \times \hat{\bm{\sigma}})_{\xi ' }              +   ( \xi \leftrightarrow \xi ' )          ]              \biggr ]           \\               & +   \textrm{h.c . }
\biggl\ }   % second       \\
\notag        = \frac{\hbar^2 \tilde{v}_{\textrm{f}}^2 } { 2 }           \sum_{\bm{k } }           & \biggl [               |h|^2 \delta_{\xi\xi ' } \hat{\sigma}_\mu   \biggl (   |d^\textrm{r}|^2 ( d^\textrm{r } + d^\textrm{a } )                   + 2 \hbar^2\tilde{v}_{\textrm{f}}^2 k^2 |d^\textrm{r}|^2 \ { ( d^\textrm{r})^2 + ( d^\textrm{a})^2 \ } \biggr )       \\
\notag                  & + \biggl [ \hbar^2\tilde{v}_{\textrm{f}}^2 k^2 |d^\textrm{r}|^2 ( d^\textrm{r } + d^\textrm{a } )                       + \hbar^4 \tilde{v}_{\textrm{f}}^4 k^4 |d^\textrm{r}|^2 \ { ( d^\textrm{r})^2 + ( d^\textrm{a})^2 \ } \biggr ] \\ \label{eq : differential-4 grgrgrga 4 } % \ref{eq : differential-4 grgrgrga 4 }                               & \   \times \bigl [   \epsilon_{\xi z \mu } ( \bm{z}\times \hat{\bm{\sigma}})_{\xi ' }               + \epsilon_{\xi ' z \mu } ( \bm{z}\times \hat{\bm{\sigma}})_\xi - \delta_{\xi \xi ' } \hat{\sigma}_\mu \bigr ]                 \biggr ] .\end{aligned}\ ] ] here , we have used @xmath381.\end{aligned}\ ] ] @xmath382 is given by using equations in appendix [ sec : e ] as @xmath383    \hat{\sigma}_\nu + o\left ( \frac{\hbar}{\epsilon_\textrm{f } \tau } \right ) ,      \\\label{eq : a38-z}%\ref{eq : a38-z } \hat{i}^{(3)}_{z \xi\xi ' }          = & \hbar^2 \tilde{v}_{\textrm{f}}^2   \frac{\pi \nu_e } { 8\epsilon_\textrm{f}^2\eta } + o\left ( \frac{\hbar}{\epsilon_\textrm{f } \tau } \right).\end{aligned}\ ] ]    using @xmath384          =         q_\alpha q_\beta \bigl [    3    \delta_{\alpha \beta } \delta_{\nu \mu } - 2 \delta_{\alpha \nu } \delta_{\beta \mu } \bigr]$ ] , we obtain @xmath385    \hat{\sigma}_\nu.\end{aligned}\ ] ]    thus , @xmath386 is obtained by @xmath387    \sigma_\nu %    \\ \hat{i}_{\mu = x , y }          & \simeq    \frac{\pi   \nu_e } { 4\eta }             \biggl [ ( 1 - i
\omega \tau   - \frac{3}{2}d\tau q^2 )   \hat{\sigma}_\mu          +   i q_\alpha \ell    \epsilon_{\mu \alpha z }              + d\tau    q_{\mu }   q_\nu \hat{\sigma}_\nu          \biggr ] ,      \\ \hat{i}_{\mu = z }          &   \simeq o\left ( \frac{\hbar}{\epsilon_\textrm{f } \tau } \right),\end{aligned}\ ] ] where @xmath388 and @xmath104 are the diffusion constant and the mean free path of the surface electrons , respectively .
we will calculate @xmath389 using the same formalism in the appendix [ sec : a ] . here , @xmath390 can be expanded with respect to @xmath81 and @xmath7 within @xmath391 and @xmath392 as @xmath393      we will calculate @xmath395 in eq .
( [ eq : coefficient pi-1 ] ) . by using @xmath396 and
@xmath397 $ ] , @xmath395 becomes @xmath398.\end{aligned}\ ] ] the above equation can be estimated around the fermi energy , @xmath346 as @xmath399   \sum_{k }    |d^\textrm{r}|^2          =   \frac{\pi\nu_e } { 2\eta } .\end{aligned}\ ] ]      @xmath400 in eq .
( [ eq : coefficient pi-2 ] ) is calculated by using @xmath401 @xmath402 , and @xmath403 as @xmath404           \\           & =   -\frac{i\pi\nu_e } { 4\eta^2}\hbar   + o\left(\frac{\hbar}{\epsilon_{\textrm{f}}\tau}\right ) .
%     \\ % \omega \hat{i}^{(1)}_{0 }    & = \frac{\pi\nu_e}{\eta}(-i\omega \tau)\end{aligned}\ ] ]      @xmath405 in eq .
( [ eq : coefficient pi-3 ] ) is obtained by @xmath406|d^\textrm{r}|^2       % = & \frac{\hbar \tilde{v}_{\textrm{f } } \pi\nu_e } { 2\eta } \epsilon_{\xi z \ell } \hat{\sigma}_\ell \biggl [   \frac{i\eta}{2\epsilon_{\textrm{f}}^2 } + \frac{i } { 4\eta }   \biggr ]      % \\
\simeq    \frac{\pi\nu_e } { 2\eta } \frac{i \tilde{v}_{\textrm{f } } \tau}{2 }   \epsilon_{\xi z \ell } \hat{\sigma}_\ell .
\end{aligned}\ ] ] here , we have used the following equations @xmath407|d^\textrm{r}|^2      % \\      & = \hbar \tilde{v}_{\textrm{f } } \epsilon_{\xi z \ell } \hat{\sigma}_\ell \bigl [   h^ *   + 2\hbar^2 \tilde{v}_{\textrm{f}}^2 k^2 \epsilon_{\textrm{f } }    d^\textrm{r } \bigr]|d^\textrm{r}|^2 .
%    \\ % \langle \hat{g}^\textrm{r}_{\bm{k } } \frac{\partial \hat{g}^\textrm{a}_{\bm{k}}}{\partial k_\xi }     \rangle_k %    & = \hbar \tilde{v}_{\textrm{f } } \epsilon_{\xi z \ell } \sigma_\ell \bigl [ h + 2\hbar^2 \tilde{v}_{\textrm{f}}^2 k^2 \epsilon_{\textrm{f } } d^\textrm{a } \bigr]|d^\textrm{r}|^2 .
\end{aligned}\ ] ]      @xmath409 becomes @xmath410 the average @xmath411 is reduced to be @xmath412 after we integrate the above equation as a function of @xmath40 , we obtain @xmath413 |d^\textrm{r}|^4 ]            - ( \hbar \tilde{v}_{\textrm{f } } k)^2 |d^\textrm{r}|^2 ( d^\textrm{a } + d^\textrm{r } ) \biggr\}q_\alpha q_\beta      \\ %    & = \delta_{\alpha\beta } \biggl\ { \frac{\pi\nu_e}{4\eta \epsilon_{\textrm{f}}^2 } %        + \frac{\pi\nu_e}{8\eta^3 }   %         + \frac{\pi\nu_e}{4 \eta \epsilon_\textrm{f}^2 } \biggr\ }      & = \frac{\pi\nu_e   } { 8\eta^3 } q^2    + o\left(\frac{\hbar}{\epsilon_{\textrm{f}}\tau}\right).\end{aligned}\ ] ] here , we have used @xmath414 % \\ & = - q_\alpha q_\beta       \delta_{\alpha \beta } ( d^\textrm{a } + d^\textrm{r}),\end{aligned}\ ] ] thus , @xmath415 becomes @xmath416
we estimate @xmath418 by using the same formalism in appendices [ sec : a ] and [ sec : b ] . here , @xmath419 can be expanded as @xmath420 where coefficients in the above equation are given by @xmath421 = \left\ { \begin{array}{ll } \displaystyle-\frac{n_\textrm{i}u_\textrm{i}^2\nu_e}{2\epsilon_f}\hat{\sigma}_\mu & ( \mu = x , y)\\ \displaystyle- \frac{in_\textrm{i}u_\textrm{i}^2\pi\nu_e}{2\epsilon_f}\arg(-\hbar\omega-\epsilon_\textrm{f}+i\eta)\hat{\sigma}_z & ( \mu = z)\\ \displaystyle-\frac{n_\textrm{i}u_\textrm{i}^2\nu_e}{\epsilon_f } + \frac{in_\textrm{i}u_\textrm{i}^2\pi\nu_e}{2\epsilon_f}\arg(-\hbar\omega-\epsilon_\textrm{f}+i\eta ) & ( \mu=0 ) \end{array}\right . ,       \\ \hat{c}^{\textrm{rr}(1)}_{\mu}\label{eq : crr1 }      = & \frac{n_\textrm{i}u_\textrm{i}^2}{2}\sum _ { k}\left [ \hat{g}^\textrm{r}\hat{\sigma}_\mu\frac{\partial \hat{g}^\textrm{r}}{\partial \omega } -\frac{\partial \hat{g}^\textrm{r}}{\partial \omega}\hat{\sigma}_\mu \hat{g}^\textrm{r } \right ]          = 0 ,       \\ \hat{c}^{\textrm{rr}(2)}_{\mu\xi}\label{eq : crr2 }      = & \frac{n_\textrm{i}u_\textrm{i}^2}{2}\sum_{k}\left[\hat{g}^\textrm{r}\hat{\sigma}_\mu\frac{\partial \hat{g}^\textrm{r}}{\partial k_\xi } -\frac{\partial \hat{g}^\textrm{r}}{\partial k_\xi}\hat{\sigma}_\mu \hat{g}^\textrm{r } \right ]      = \left\ {          \begin{array}{ll }          \displaystyle\frac{n_\textrm{i}u_\textrm{i}^2\hbar v_f\nu_e}{2\epsilon_f[(\hbar\omega+\epsilon_f)^2+\eta^2]}[i(\hbar\omega+\epsilon_f)+\eta]\delta_{\xi \mu }     \hat{\sigma}_z & ( \mu = x , y )              \\
\displaystyle-\frac{n_\textrm{i}u_\textrm{i}^2\hbar v_f\nu_e}{2\epsilon_f[(\hbar\omega+\epsilon_f)^2+\eta^2]}[i(\hbar\omega+\epsilon_f)+\eta]\hat{\sigma}_\xi & ( \mu = z )              \\          0 & ( \mu=0 )          \end{array}\right . ,       \\ \notag   \hat{c}^{\textrm{rr}(3)}_{\mu\xi\xi ' }      = & \frac{n_\textrm{i}u_\textrm{i}^2}{4}\sum_{k}\left[\hat{g}^\textrm{r}\hat{\sigma}_\mu\frac{\partial^2 \hat{g}^\textrm{r}}{\partial k_\xi\partial k_{\xi ' } } + \frac{\partial^2 \hat{g}^\textrm{r } } { \partial k_\xi\partial k_{\xi'}}\hat{\sigma}_\mu \hat{g}^\textrm{r}\right ]      \\\label{eq : crr3 }      = & \left\ { \begin{array}{ll } \displaystyle\frac{n_\textrm{i}u_\textrm{i}^2\hbar^2v_f^2\nu_e}{12\epsilon_f[(\hbar\omega+\epsilon_f)^2+\eta^2]^2}\left[(\hbar\omega+\epsilon_{\textrm{f}})^2-\eta^2 - 2i\eta(\hbar\omega+\epsilon_{\textrm{f}})\right]&\\      \displaystyle\hspace{20mm}\times[\epsilon_{\xi z \mu}(\bm{z}\times\hat{\bm{\sigma}})_{\xi'}+\epsilon_{\xi ' z \mu}(\bm{z}\times\hat{\bm{\sigma}})_\xi -2\delta_{\xi\xi'}\hat{\sigma}_\mu]&(\mu = x , y)\\ \displaystyle-\frac{n_\textrm{i}u_\textrm{i}^2\hbar^2v_f^2\nu_e}{4\epsilon_f[(\hbar\omega+\epsilon_f)^2+\eta^2]^2}\left[(\hbar\omega+\epsilon_{\textrm{f}})^2-\eta^2 - 2i\eta(\hbar\omega+\epsilon_{\textrm{f}})\right]\delta_{\xi\xi'}\hat{\sigma}_z&(\mu = z)\\ \displaystyle\frac{n_\textrm{i}u_\textrm{i}^2\hbar^2v_f^2\nu_e}{12\epsilon_f[(\hbar\omega+\epsilon_f)^2+\eta^2]^2}\left[(\hbar\omega+\epsilon_{\textrm{f}})^2-\eta^2 - 2i\eta(\hbar\omega+\epsilon_{\textrm{f}})\right]\delta_{\xi\xi'}\hat{\sigma}_0&(\mu=0 ) \end{array}\right . .\end{aligned}\ ] ] from the eqs . (
[ eq : crr0])-([eq : crr3 ] ) and @xmath422 , the elements of @xmath417 is negligibly small as compared with the ones of @xmath423 , since @xmath82 is satisfied .
we estimate the response function composed of only the retarded ( advanced ) green s functions @xmath425 ( @xmath426 ) using the same formalism in appendices [ sec : a ] , [ sec : b ] and [ sec : a - c ] . from eqs .
( [ eq:3 - 11 ] ) and ( [ eq:3 - 11 - 2 ] ) , @xmath427 are written as @xmath428       \\&\hspace{20 mm }                  + \frac{1}{2 } \omega f'_{\omega } \biggl [    \hat{g}^\textrm{r}_{\bm{k}-\frac{\bm{q}}{2},\omega-\frac{\omega}{2 } } \hat{\sigma}_\nu   \hat{g}^\textrm{r}_{\bm{k}+\frac{\bm{q}}{2},\omega+\frac{\omega}{2 } }           +   \left(\hat{g}^\textrm{r}_{\bm{k}+\frac{\bm{q}}{2},\omega+\frac{\omega}{2 } } \hat{\sigma}_\nu \hat{g}^\textrm{r}_{\bm{k}-\frac{\bm{q}}{2},\omega-\frac{\omega}{2}}\right)^\dag \biggr]\biggr\}. \label{eq : pirr+aa}%\ref{eq : pirr+aa}\end{aligned}\ ] ] the magnitude of second term of the above equation is smaller than that of @xmath91 [ see eqs .
( [ eq:3 - 10 ] ) and appendices [ sec : a ] , [ sec : b ] and [ sec : a - c ] ] .
the first term is expanded as @xmath429       \\       = & \hat{d}^{(0)}_{\nu}+\omega \hat{d}^{(1)}_{\nu}+\sum_{\xi = x , y}q_\xi \hat{d}^{(2)}_{\nu\xi}+o(q^2 , q\omega , \omega^2),\end{aligned}\ ] ] where @xmath430 , @xmath431 and @xmath432 in the above equation are given by @xmath433 = 0 ,      \\\label{eq : d1nu}%\ref{eq : d1nu } \hat{d}^{(1)}_{\nu }      = & -\frac{1}{2}\sum_{\bm{k } , \omega}\left [ \left(\hat{g}_{\bm{k},\omega}^\textrm{r}\hat{\sigma}_\mu\frac{\partial \hat{g}_{\bm{k},\omega}^\textrm{r}}{\partial \omega } -\frac{\partial \hat{g}_{\bm{k},\omega}^\textrm{r}}{\partial \omega}\hat{\sigma}_\mu \hat{g}_{\bm{k},\omega}^\textrm{r}\right ) + \textrm{h.c . }
\right ]          = 0 ,      \\\label{eq : d2nu}%\ref{eq : d2nu } \hat{d}^{(2)}_{\nu\xi }      = & -\frac{1}{2}\sum_{\bm{k } , \omega}f_{\omega}\left[\left(\hat{g}_{\bm{k},\omega}^\textrm{r}\hat{\sigma}_\mu\frac{\partial \hat{g}_{\bm{k},\omega}^\textrm{r}}{\partial k_\xi } -\frac{\partial \hat{g}_{\bm{k},\omega}^\textrm{r}}{\partial k_\xi}\hat{\sigma}_\mu \hat{g}_{\bm{k},\omega}^\textrm{r}\right ) + \textrm{h.c . } \right ]      = -\frac{\pi \tilde{v}_{\textrm{f}}\nu_{\textrm{e}}}{2\epsilon_{\textrm{f}}}\hat{\sigma}_z + o\left(\frac{\hbar}{\epsilon_{\textrm{f}}\tau}\right).\end{aligned}\ ] ] here , we have used @xmath434 in the above equation , where @xmath435 is a step function . from eqs . ( [ eq:3 - 2 ] ) , ( [ eq:3 - 6 ] ) , ( [ eq:3 - 33 ] ) and ( [ eq : d2nu ] ) , there are spin density and charge current induced by the magnetic field @xmath436 . since we consider only the electric field , we ignore the densities induced by the magnetic field .    we find that the order of @xmath437 are @xmath438 and @xmath439 are negligibly small as compared with @xmath91 .
we will show the following integrals as a functions of @xmath40 .
the integrals are obtained by @xmath440    d^\textrm{r } ( d^\textrm{a})^2       & \simeq   \frac{\nu_e}{8\eta^2 \epsilon_{\textrm{f } } }           \biggl [ i\pi - \frac{2\eta^2}{\epsilon_{\textrm{f}}^2 }   \biggr ] ,       \\\label{eq : d5}%\ref{eq : d5 } \sum_{k }   [   h(h^*)^2 + ( h - 2h^*)\hbar^2 \tilde{v}_{\textrm{f}}^2 k^2 ] d^\textrm{r } ( d^\textrm{a})^2       & \simeq \frac{-i\pi\nu_e}{4\epsilon_\textrm{f}^2}\left(1+i\frac{\eta}{\epsilon_{\textrm{f } } } \right ) ,      \\
\sum_k \hbar^2\tilde{v}_{\textrm{f}}^2 k^2   |d^\textrm{r}|^2 ( d^\textrm{a } - d^\textrm{r } )      & \simeq \frac { i\pi \nu_e } { 8\eta^2 \epsilon_{\textrm{f } } } ,       \\ \sum_k   |d^\textrm{r}|^2       & \simeq \frac{\pi\nu_e}{4\eta \epsilon_{\textrm{f}}^2 } ,      \\ % 1 \sum_k \bigl [ ( d^\textrm{r})^2 d^\textrm{a } + ( d^\textrm{a})^2 d^\textrm{r } \bigr ]       & \simeq   -\frac{\nu_e}{2 \epsilon_{\textrm{f}}^5 } ,
\\ \sum_k ( \hbar \tilde{v}_{\textrm{f } } k ) ^2 \bigl [ ( d^\textrm{r})^2 d^\textrm{a } + ( d^\textrm{a})^2 d^\textrm{r } \bigr ]       & \simeq - \frac { \pi \nu_e } { 4 \eta \epsilon_{\textrm{f}}^2 } ,       \\
\sum_k ( \hbar \tilde{v}_{\textrm{f } } k ) ^2 ( d^\textrm{r})^3 d^\textrm{a }       & \simeq   -\frac{\pi \nu_e}{64 \epsilon_{\textrm{f}}^2 \eta^3 } \left ( 1- i\frac{2\eta}{\epsilon_{\textrm{f } } } \right ) ,        \\ % 4 \sum_k ( \hbar \tilde{v}_{\textrm{f } } k ) ^2 \ { ( d^\textrm{r})^3 d^\textrm{a } + ( d^\textrm{a})^3 d^\textrm{r}\ }      & \simeq -\frac{\pi \nu_e}{32 \epsilon_{\textrm{f}}^2 \eta^3 } ,      \\ \sum_k ( \hbar \tilde{v}_{\textrm{f } } k ) ^4 ( d^\textrm{r})^3 d^\textrm{a }      & \simeq -\frac{\pi \nu_e } { 32 \eta^3 } ,       \\ \sum_k(\hbar v_fk)^2(d^\textrm{r})^{n }      & \simeq-\frac{1}{n-1}\sum_k(d^\textrm{r})^{n-1 } \ \ \ \ \ ( n\geq3 ) ,      \\ \sum_k(\hbar v_fk)^4(d^\textrm{r})^{n } & = \frac{2}{(n-1)(n-2)}\sum_k(d^\textrm{r})^{n-2 } \ \ \ \ \ ( n\geq4 ) ,      \\ \sum_k(\hbar v_fk)^2(d^\textrm{r})^{3 }      & \simeq-\frac{1}{2}\sum_k(d^\textrm{r})^{2 } ,      \\ \sum_k(\hbar v_fk)^2(d^\textrm{r})^{4 }      & \simeq -\frac{1}{3}\sum_k(d^\textrm{r})^{3}= \frac{\pi\nu_e } { 32\epsilon_{\textrm{f}}^2 \eta^3 } ,      \\ \sum_k(\hbar v_fk)^4(d^\textrm{r})^{4 }      & \simeq \frac{1}{3}\sum_k(d^\textrm{r})^{2 } = \frac{\pi\nu_e } { 32 \eta^3}.\end{aligned}\ ] ] where @xmath441 is defined by @xmath442 here , we have used @xmath443 in the above equation .
to check validity of our results , we substitute the charge current and charge density in eqs .
( [ eq:3 - 32 ] ) and ( [ eq:3 - 24 ] ) into the charge conservation law @xmath143 . from eq .
( [ eq:3 - 23 ] ) , @xmath444 becomes @xmath445 }           \frac { i \omega^2 q_\nu   } { q^2 \ell^2 + i\omega \tau } a_{\textrm{em},\nu } .
%        \\ %    & =   \frac{\hbar e \nu_e j_{ex } } { \eta v } \sum_{\bm{q},\omega}e^{i[\omega t-\bm{q}\cdot\bm{x } ] }   %        \biggl [ %        \frac { i \ell \omega^2}{q^2 \ell^2 + i\omega \tau } \epsilon_{\zeta \nu z } q_\zeta   s^\nu_{\bm{q},\omega } %        \biggr]\end{aligned}\ ] ] from eq .
( [ eq:3 - 21 ] ) , ( [ eq:3 - 22 ] ) and @xmath446 , @xmath447 becomes @xmath448 }             \biggl [   - \omega q_x a_{{\textrm{em } } , x } + \biggl\ { \frac { \omega q_x^2 q_y \ell^2 } { q^2 \ell^2 + i\omega \tau } a_{{\textrm{em } } , y } + \frac { \omega q_x^3 \ell^2 } { q^2 \ell^2 + i\omega \tau } a_{{\textrm{em } } , x }   \biggr\ } \biggr ] ,             \\
\nabla_y j_{y }   & = \frac{e^2 \tilde{v}_{\textrm{f}}^2 \nu_e \tau } { l^2 }             \sum_{\bm{q},\omega}e^{i[\omega t-\bm{q}\cdot\bm{x } ] }             \biggl [   -\omega q_y a_{{\textrm{em } } , y } + \biggl\ { \frac { \omega q_x q_y^2 \ell^2 } { q^2 \ell^2 + i\omega \tau } a_{{\textrm{em } } , x } +   \frac{\omega q_y^3 \ell^2 } { q^2 \ell^2 + i\omega \tau } a_{{\textrm{em } } , y }   \biggr\ } \biggr ] ,       \\   % % % % % % % \nabla_x j_x + \nabla_y j_y         & =    -\frac { e^2 \tilde{v}_{\textrm{f}}^2 \ell \nu_e } { l^2 }             \sum_{\bm{q},\omega}e^{i[\omega t-\bm{q}\cdot\bm{x } ] }              \frac { i \omega^2   q_\nu        } { q^2 \ell^2 + i\omega \tau }     a_{{\textrm{em } } , \nu}. %        \\ %        & =   \frac { \hbar e   \nu_e j_{ex } } { \eta   v }   \sum_{\bm{q},\omega}e^{i[\omega t-\bm{q}\cdot\bm{x } ] } \biggl [ \frac { i\omega^2 \tilde{v}_{\textrm{f } }   \tau } { q^2 \ell^2 + i\omega \tau }   \epsilon_{\zeta \nu z } q_\zeta s^\nu_{\bm{q},\omega } \biggr].\end{aligned}\ ] ] therefore , @xmath128 and @xmath164 follow the charge conservation law , @xmath449 .
diffusive green s function on the disordered surface of the ti can be integrated as follows : @xmath450 } ,       \\
\notag   \sum_{q_x } e^{-(2d t q_x^2 + i q_x x ) }       & \sim \frac{1}{2\pi}\int_{-\infty}^\infty dq_x e^{-(2d t q_x^2 + i q_x x ) } %    \\
\notag   %    & = \frac{e^{- \frac{2x^2}{dt}}}{2\pi}\int_{-\infty}^\infty dq_x e^{-2d t \ { [ q_x + i \frac{x}{dt}]^2 \ } }      \\      & = \frac{\sqrt{\pi}}{2\pi \sqrt{2dt } } \exp{\left(- \frac{2x^2}{dt}\right)}.\end{aligned}\ ] ] thus , from the above equations , @xmath135 in the coordinates space is obtained by @xmath451}.\end{aligned}\ ] ]    a. v. kimel , a. kirilyuk , p. a. usachev , r. v. pisarev , a. m. balbashov , and th .
rasing : nature ( london ) * 435 * , 655 ( 2005 ) .
f. hansteen , a. kimel , a. kirilyuk , and th .
rasing : phys .
lett . * 95 * , 047402 ( 2005 ) : phys .
b * 73 * , 014421 ( 2006 ) .
r. iida , t. satoh , t. shimura , k. kuroda , b. a. ivanov , y. tokunaga , and y. tokura : phys .
b * 84 * , 064402 ( 2011 ) .
a. kirilyuk , a. v. kimel , and t. rasing : rev .
phys . * 82 * , 2731 ( 2010 ) . c. d. stanciu , f. hansteen , a. v. kimel , a. kirilyuk , a. tsukamoto , a. itoh , and t. rasing , phys .
99 * , 047601 ( 2007 ) .
k. vahaplar _
et al_. , phys .
lett . * 103 * , 117201(2009 ) .
s. raghu , s. b. chung , x .-
qi , and s .- c .
zhang , phys .
lett . * 104 * , 116401 ( 2010 ) .
t. misawa , t. yokoyama , and s. murakami : phys .
b * 84 * , ( 2011 ) 165407 .
x. liu and j. sinova , phys .
lett . * 111 * , 166801 ( 2013 ) .
g. f. quinteiro and p. i. tamborenea , phys .
b * 79 * , 155450 ( 2009 ) . g. f. quinteiro and p. i. tamborenea , phys . rev .
b * 82 * , 125207 ( 2010 )
. g. f. quinteiro , p. i. tamborenea , and j. berakdar , opt .
express * 19 * 26733 ( 2011 ) .
j. wtzel , a. s. moskalenko , and j. berakdar , opt .
express * 20 * 27792 ( 2012 ) .
m. b. faras , g. f. quinteiro , and p. i. tamborenea , eur .
j. b * 86 * , 432 ( 2013 ) .
m. cygorek , p. i. tamborenea , and v. m. axt , phys .
b * 92 * , 115301 ( 2015 ) .
y. ueno , y. toda , s. adachi , r. morita , and t. tawara , opt .
express * 17 * , 20567 ( 2009 ) .
n. b. clayburn , j. l. mccarter , j. m. dreiling , m. poelker , and t. j. gay , phys .
b * 87 * , 035204 ( 2013 ) .
l. marrucci , c. manzo , and d. paparo , phys .
lett . * 96 * , 163905 ( 2006 ) .
l. allen , m. w. beijersbergen , r. j. c. spreeuw , and j. p. woerdman , phys .
a * 45 * , 8185 ( 1992 ) .
g. molina - terriza , j. p. torres , and l. torner , nature phys . * 3 * , 305 - 310 ( 2007 ) .
m. padgett , j. courtial , and l. allen , phys . today 57(5 ) , 35 ( 2004 ) .
j. t. mendonca , b. thide , and h. then , phys .
* 102 * , 185005 ( 2009 ) .
v. e. lembessis , m. babiker , and d. l. andrews , phys .
a * 79 * , 011806(r ) ( 2009 ) .
g. d. mahan , _ many - particle physics _
( springer , 2000 ) , p.514 .
j. fujimoto , a. sakai , h. kohno , phys .
b * 87 * , 085437 ( 2013 ) .
a. sakai and h. kohno , phys .
b * 89 * , 165307 ( 2014 ) .
k. taguchi , k. shintani , y. tanaka , phys .
b * 92 * , 035425 ( 2015 ) .      a. a. burkov and d. g. hawthorn , phys
* 105 * , 066802 ( 2010 ) .
p. schwab , r. raimondi , and c. gorini , epl ( europhysics letters ) , * 93 * , 67004 ( 2011 ) .
k. taguchi , t. yokoyama and y. tanaka , phys .
b * 89 * , 085407 ( 2014 ) .
s. murakami , n. nagaosa , and s. c. zhang , science 301 , 1348 - 1351 ( 2003 ) .
j. sinova , d. culcer , q. niu , n. a. sinitsyn , phys .
lett . * 92 * , 126603 ( 2004 ) .
e. l. rashuba , phys .
b * 68 * , 241315 ( 2003 ) . t. kimura , y. otani , k. tsukagoshi , and y. aoyagi , j. magn .
272 - 276 , e1333 ( 2004 ) .
y. k. kato , r. c. myers , a. c. gossard , and d. d. awschalom , science 306 , 1910 , e1333 ( 2004 ) .
e. saitoh , m. ueda , h. miyajima , and g. tatara , appl .
. lett . * 88 * , 182509 ( 2006 ) .
a. a. taskin , s. sasaki , k. segawa , and y. ando , phys .
lett . * 109 * , 066803 ( 2012 ) .
using @xmath452 , @xmath453 , and @xmath454 , we can obtain the relation : @xmath455 .
because of @xmath456 , @xmath457 becomes negligible small as @xmath458 .
streakily speaking , there is no problem to take @xmath459 .