topological transport phenomena have been attracting a great deal of interest , because they exhibit universal properties that are insensitive to perturbations and independent of material details .
a classical example of such a transport phenomenon is the integer quantum hall ( iqh ) effect in two - dimensional ( 2d ) electron systems , first discovered in 1980 ,  @xcite which is characterized by an integer quantization of the hall conductivity in unit of @xmath2 .
the iqh effect has been observed in a large variety of materials , ranging from traditional semiconductors , to oxides ,  @xcite graphene ,  @xcite and topological insulators ( tis ) .
@xcite laughlin  @xcite interpreted the iqh effect in terms of an adiabatic charge pump .
thouless , kohmoto , nightingale , and nijs  @xcite established a relation between the quantized hall conductivity of the iqh system and a topological invariant , the first chern number .
thouless and niu  @xcite also related the amount of charge pumped in a 1d charge pump to the chern number .
a variant of the iqh effect , the quantum spin hall ( qsh ) effect , was proposed recently ,  @xcite which has been experimentally realized in hgte quantum wells  @xcite and inas / gasb bilayers .
@xcite extension of the idea of the qsh effect has led to the discovery of 3d tis .
@xcite a qsh system , which is also called a 2d ti , has an insulating band gap in the bulk and a pair of gapless helical edge states at the sample boundary .
when the electron spin is conserved , a qsh system can be viewed as two independent iqh systems without landau levels .
@xcite different from the charge , the spin does not obey a fundamental conservation law .
in general , when the spin conservation is absent , unconventional topological invariants , either the @xmath3 index  @xcite or the spin chern numbers ,  @xcite are needed to describe the qsh systems . the time - reversal ( tr ) symmetry
is considered to be a prerequisite for the qsh effect , which protects both the @xmath3 index and gapless nature of the edge states .
however , based upon the spin chern numbers , it was shown that the bulk topological properties remain intact even when the tr symmetry is broken .
@xcite this finding evokes interest to pursue direct investigation and possibly utilization of the robust topological properties of the tis , besides using their symmetry - protected gapless edge states which are more fragile in realistic environments .    unlike the first chern number underlying the iqh systems , which is embedded into the hall conductivity , up to now the topological invariants in the tis
have not been directly observable .
several experimental methods were proposed , but have not been realized .
one was to measure the topological magnetoelectric effect ,  @xcite for which experimental complexities exist .
@xcite fu and kane  @xcite put forward an abstract 1d model , in which the spin pumping was related to the @xmath4 index in the limit of weak coupling .
however , how this fictitious model could be implemented is still unknown .
furthermore , from the viewpoint of application , generalization of the idea of the @xmath4 pump to higher dimension is meaningless , because according to the @xmath4 theory ,  @xcite only the states at the tr - invariant point of the brillouin zone can contribute to the spin pumping , and so the pumping rate can not be enhanced by increment of dimension . in a recent work ,  @xcite the more general case of finite coupling between the pump and electrode is investigated by using the scattering matrix method .
it was found that the spin pumping in the model of fu and kane can survive finite scattering of magnetic impurities , and so may be attributed to the spin chern numbers rather than the @xmath3 index .
some other authors  @xcite proposed to pump quantized charge through the helical edge states by precessing a magnet covering the edge of a 2d ti , so that the number of gapless edge channels can be counted through electrical measurement .
this method is indirect , in the sense that the topological invariants are intrinsic properties of the bulk electron wavefunctions , which do not immediately determine the charge pumping in the edge channels .    here
we predict an intriguing bulk topological pumping effect , directly driven by nonzero spin chern numbers , in a qsh system electrically modulated in adiabatic cycles . as a consequence of the topological spectral flows of the spin - polarized wannier functions ( spwfs ) in the bulk of the system
, spin can be pumped into a nonmagnetic electrode continuously without net charge transfer .
the total amount of spin pumped per cycle is proportional to the ( cross - section ) width of the sample , and insensitive to the material parameters and spin - mixing effect due to the rashba spin - orbit coupling .
this electrical spin pump establishes a basis , on which spintronic applications taking advantage of the robust topological properties of the tis can be developed .
especially , if a half - metallic electrode with spin polarization parallel ( or antiparallel ) to the @xmath5-axis is used , a quantized charge pumping conductivity , @xmath6 ( or @xmath7 ) , can be measured by electrical means , demonstrating a way to observe the spin chern numbers @xmath0 directly .
let us consider a 2d model hamiltonian @xmath8 with @xmath9-m(t)\hat{\sigma}_{z}\ .
\label{hamil}\ ] ] here @xmath10 is the electron charge , @xmath11 is the 2d momentum , @xmath12 is the vector potential of an @xmath13 electric field @xmath14 applied along the @xmath15 direction with @xmath16 and frequency @xmath17 being designated , and @xmath18 .
this model can describe both the qsh materials , the hgte quantum wells ,  @xcite and inas / gasb bilayers ,  @xcite in the linear order in momentum .
for definity of discussion , we confine ourselves to the hgte quantum wells , for which @xmath19 with @xmath20 are the pauli matrices for spin , and @xmath21 for the electron and hole bands .
as will be discussed below , the time - dependent mass term @xmath22 can be induced by varying the voltages of the dual gates .
@xmath23 represents the rashba spin - orbit coupling  @xcite @xmath24\ .
\label{rashba}\ ] ] to the linear order in momentum , the rashba spin - orbit coupling is nonvanishing only in the electron band .
@xcite    within the adiabatic approximation , for a bulk sample there exists a finite energy gap between the conduction and valence bands for @xmath25 or @xmath26 . at @xmath27 and @xmath26 ,
the conduction and valence bands touch at @xmath28 and @xmath29 or @xmath30 with @xmath31 . to clarify the topological properties underlying the spin / charge pumping , we consider @xmath32 as a parameter , and calculate the spin chern numbers @xmath0 in the standard way ,  @xcite on the torus of the two variables @xmath33 and @xmath34 with @xmath35 as the period .
the spin chern numbers are obtained as @xmath36 for @xmath37 , and vanish elsewhere .
not surprisingly , the band touching points @xmath38 serve as the critical points .
( vertical axis ) .
the parameters are taken to be @xmath39 , @xmath40 , @xmath41 , and @xmath42 , with @xmath43 as the lattice constant and @xmath44 ( @xmath45 ) as the hopping integral of the tight - binding hamiltonian . ,
width=288 ]    we now consider a system consisting of a pump for @xmath46 , an electrode for @xmath47 , and a potential barrier in between .
the total hamiltonian of the system reads @xmath48 where @xmath49 has been given above , and @xmath50 is the hamiltonian of the electrode .
a possible experimental setup for realizing this hamiltonian is explained in appendix a in more details . in the barrier region ,
the term @xmath51 opens an insulating gap of size @xmath52 , which accounts for contact deficiencies between the pump and electrode .
the crucial role of the nonzero spin chern numbers in the spin / charge pumping process can be visualized by using the spwfs , which were first introduced in ref .
we construct a tight - binding hamiltonian for the effective 1d system at any given @xmath53 according to eq .
( [ htot ] ) , and diagonalize the total hamiltonian of the pump and electrode numerically . following the same procedure as calculating the spin chern numbers ,  @xcite the space occupied by electrons is partitioned into two spin sectors after diagonalizing the spin operator @xmath54 in the occupied space . by definition ,
the states in the two spin sectors are essentially the maximally spin - polarized states
. then we construct the wannier functions  @xcite for the spin - up and spin - down sectors , respectively , which are called the spwfs .    the evolution of the centers of mass of the spwfs for @xmath55 and @xmath56 is shown in fig .
[ fig_1 ] .
we see that the wannier centers for the spin - up sector move right and those for the spin - down sector move left , each center shifting on average a lattice constant per cycle . within the adiabatic approximation , time @xmath57 plays the same role as the momentum of an additional dimension ,  @xcite namely , @xmath58 .
therefore , when @xmath53 is considered as a parameter , the evolution of the wannier functions of the effectively 1d system related to various @xmath59 with time @xmath60 can be understood from the static properties of a 2d system associated with various @xmath59 and @xmath61 . in the general theory ,  @xcite the relationship between the chern number and the spectral flows of the wannier functions in a 2d system has been established . according to this theory , the average displacement of each of the centers of the spwfs in the spin - up ( spin - down ) sector with changing @xmath61 ( or @xmath60 ) from @xmath62 to @xmath63 , in units of the lattice constant , must equal to the spin chern number @xmath64 ( @xmath65 ) . therefore , the nontrivial transfer of the spwfs observed in fig .
1 is a direct manifestation of the nonzero spin chern numbers @xmath66 in the pump ( for @xmath67 ) .
more interestingly , we see that such spectral flows can go across the finite barrier ( @xmath68 ) , and extend into the electrode , even though the barrier and electrode are topologically trivial .
physically , because the system needs to recover its original eigenstates when each cycle ends , the nontrivial spectral flows of the spwfs in the ti need to constitute closed loops through formation of edge states at the boundary ,  @xcite or extend into the electrode .
however , localized edge states can not exist at the finite barrier due to quantum tunneling effect , so the transfer of the spectral flows of the spwfs into the electrode occurs .
this result will be further confirmed by direct calculation based upon the scattering matrix theory in the next section .
the spwfs are just another equivalent representation of the occupied space , and so the counter spectral flows of the wannier centers in the two spin sectors represent the true movements of the electrons .
if the rashba spin - orbit coupling were neglected , the wannier functions would be the eigenstates of @xmath69 .
the nontrivial spectral flows indicate that at the given @xmath53 , in each cycle a spin - up electron goes from the pump into the electrode , and a spin - down electron moves oppositely .
therefore , no net charge transfer occurs but a quantized spin of @xmath70 is pumped into the electrode . when the small rashba spin - orbit coupling is turned on , while the topological spectral flows remain intact , as seen from fig .
[ fig_1 ] , the spin polarizations of the wannier functions are no longer fully parallel to the @xmath5-axis , and may also vary with time . as a consequence ,
the amount of spin pumped per cycle will deviate from the quantized value .
in general , the amount of the spin pumped can be conveniently calculated by using the scattering matrix formula .
@xcite the @xmath5-component of the spin pumped per cycle is given by  @xcite    @xmath71    where @xmath72 ( @xmath73 ) is the reflection amplitude for an electron at the fermi energy incident from the spin-@xmath74 channel of the electrode and reflecting back into the spin-@xmath75 channel . in the following calculations , the fermi energy is set to be zero ( @xmath76 ) , and the rashba spin - orbit coupling is treated as a perturbation . as shown in appendix b , to the linear order in @xmath77 , we obtain @xmath78}{\mbox{ch}(2\gamma_{0}d)-\sin(2\theta ) \mbox{sh}(2\gamma_{0}d)}+{\cal o}(\epsilon^2)\ , \label{rupup0}\ ] ]    @xmath79 } { \mbox{ch}(2\gamma_{0}d)-\sin(2\theta)\mbox{sh}(2\gamma_{0}d ) } + { \cal o}(\epsilon^2)\ , \label{rdnup0}\ ] ]    and @xmath80 and @xmath81 , where @xmath82 and @xmath83 $ ] .
we note that the dimensionless quantity @xmath84 appears as the small expansion parameter . since @xmath85 and @xmath86 are always real , the contributions from the third and fourth terms in eq .
( [ ds1 ] ) vanish .
consequently , @xmath87 this expression has a geometric explanation : the amount of spin pumped per cycle equals to the difference between the areas enclosed by the directional trajectories of @xmath88 and @xmath89 on the complex plane , multiplied by @xmath90 . due to the relation @xmath91 , yielding @xmath92 , the two terms in eq .
( [ ds2 ] ) make an equal contribution , so that we can focus on the first term . while the expression ( [ rupup0 ] ) for @xmath88 is independent of the rashba spin - orbit coupling , as will be shown soon , a combination of eqs .
( [ rupup0 ] ) and ( [ rdnup0 ] ) allows us to evaluate the amount of spin pumped up to the second order in @xmath93 .
, as a function of @xmath94 for four sets of @xmath95 .
the other parameters are taken to be @xmath96 and @xmath97 with @xmath98 .
inset : trajectories of @xmath88 in a cycle on the complex plane.,width=230 ]    we first consider the case of @xmath96 . from eq .
( [ rupup0 ] ) , it is easy to show @xmath99 . in fig .
[ fig_2 ] , we plot the argument @xmath100 of @xmath101 as a function of @xmath94 for several parameter sets . for either @xmath102 ( ideal contact ) or @xmath103 ( strong potential barrier ) , @xmath100 always increments @xmath104 in a cycle as long as @xmath105 . in this case
, the trajectories of @xmath101 always form a unit circle on the complex plane , oriented counterclockwise , as shown in the inset of fig .
[ fig_2 ] , suggesting @xmath106 ( for @xmath67 ) . for @xmath107 , however , the situation is quite different .
@xmath100 does not change after going through a cycle , and the trajectory of @xmath101 does not enclose a finite area , so that @xmath108 .
apparently , the present result conforms to the spin chern numbers given by eq .
( [ schern ] ) and the spectral flows of the spwfs .
( in unit of @xmath109 ) as a function of @xmath110 for @xmath39 and three different values of @xmath77 .
the other parameters are taken to be the same as in fig .
[ fig_2 ] .
inset : the trajectory of @xmath88 in a cycle for @xmath56 and @xmath111 , with the unit circle indicated by the dotted line.,width=249 ]    next we study the correction to @xmath112 due to nonzero rashba spin - orbit coupling . by expressing @xmath113 in the polar coordinate system , where @xmath114 is the argument at @xmath115 , and @xmath116 and @xmath117 stand for
the second - order corrections to @xmath118 and @xmath119 , respectively , due to the rashba spin - orbit coupling , eq .
( [ ds2 ] ) becomes @xmath120 + { \cal o}(\epsilon^3)$ ] .
we notice that @xmath116 is a small quantity fluctuating around @xmath62 and periodic in time , @xmath121 , so that @xmath122 . using the identity @xmath123
, we then obtain @xmath124 , where @xmath125 can be calculated from eq .
( [ rupup0 ] ) and @xmath85 has been given by eq .
( [ rdnup ] ) .
this is an expression for @xmath112 accurate to the second order in @xmath93 . at @xmath102
, we obtain @xmath112 as @xmath126(c_{+}-c_{-})\ , \label{ds_r}\ ] ] for @xmath37 , and @xmath108 elsewhere .
as expected , nonzero rashba spin - orbit coupling causes @xmath112 to deviate from its quantized value , i.e. , @xmath127 . in real materials
, @xmath77 is usually much smaller ( by an order of magnitude or more ) than @xmath128 ,  @xcite so that the deviation is less than a percent for an ideal connection between the pump and electrode .
for @xmath129 , @xmath112 can be evaluated numerically and its calculated result is plotted in fig .
[ fig_3 ] as a function of @xmath110 for three different strengths of the rashba spin - orbit coupling . for @xmath96 ,
@xmath112 is quantized to @xmath130 , independent of @xmath110 . for @xmath131 and @xmath132 , weak potential barrier ( @xmath133 )
has little effect on @xmath112 .
this is reasonable as the leading - order correction of small @xmath110 must be @xmath134 .
appreciable deviations from the quantized value occur for strong potential barrier ( e.g. , @xmath135 ) .
we note that @xmath116 does not affect the orbit of @xmath88 , as it represents a variation in the tangent direction of the orbit , and the orbit can be determined by @xmath136 .
inset shows the trajectory of @xmath88 on the complex plane for @xmath56 and @xmath111 .
for such a strong potential barrier , the orbit of @xmath88 deviates from the unit circle visibly .
the above result suggests that improving the contact quality between the pump and electrode is helpful for obtaining a nearly integer - quantized value of the pumped spin .    by summing over @xmath53 between @xmath30 and @xmath137
, we obtain for the total spin pumped per cycle @xmath138 with @xmath139 ( c_{+}-c_{-})\ , \ ] ] for a good contact ( @xmath133 ) .
@xmath140 is in scale with width @xmath141 of the pump . by noting that @xmath142 is proportional to the number of cycles per unit time
, @xmath143 can be considered as the spin pumping conductivity .
now we discuss a possible way to experimentally observe the spin chern numbers , by using a half - metallic electrode , in which conducting channels for electron spin antiparallel to the spin polarization are absent .
we first consider the case , where the spin polarization of the electrode is parallel to the @xmath5 axis .
the hamiltonian of the electrode is taken as @xmath144 . in this case , as shown in appendix c , @xmath88 is still given by eq .
( [ rupup0 ] ) but @xmath145 .
it follows that for any @xmath32 between @xmath146 and @xmath147 , the charge pumped per cycle is integer - quantized and equal to @xmath148 .
similarly , for the spin polarization of the electrode antiparallel to the @xmath5 axis , the charge pumped is equal to @xmath149 .
therefore , the total charge pumped per cycle is given by @xmath150 with @xmath151 where the spin chern number @xmath152 ( @xmath153 ) is taken for the spin polarization of the electrode parallel ( antiparallel ) to the @xmath5 axis .
we emphasize that eqs .
( [ totcharge ] ) and ( [ ccond ] ) obtained above are valid for finite rashba spin - orbit coupling and finite potential barrier between the pump and electrode , indicating that the quantized charge pumping is robust against small perturbations .
experimentally , @xmath154 can be obtained by measuring the electrical current in the electrode , and from eqs .
( [ totcharge ] ) and ( 13 ) , @xmath0 can be evaluated , yielding an experimental method to measure the spin chern numbers directly .
the sign inversion of @xmath154 with reversing the spin polarization of the half - metallic electrode , as indicated by eqs .
( [ totcharge ] ) and ( [ ccond ] ) , is a hallmark of the present spin chern charge pump , which can be used to distinguish it from the conventional thouless charge pump  @xcite .
we have used the single - electron approximation , where the electron interaction is not taken into account . in particular ,
in the half - metallic electrode case , one of the spin channel is blocked at the boundary , which naturally induces some charge and spin accumulations , and consequently changes the potential profile . however , owing to the screening effect , the change in the potential profile is expected to be localized at the boundary , which in effect modifies the potential barrier between the pump body and the electrode . as has been shown above , the charge pumping effect is independent of the existence and details of the potential barrier , and so we believe that the pumping effect will survive the charge and spin accumulations .
up to now , all the results obtained from the scattering matrix formula are apparently in complete agreement with the spin chern numbers given by eq .
( [ schern ] ) .
these results can not be explained within the framework of the @xmath3 theory .
@xcite while one can define a @xmath3 index at the tr - invariant point @xmath155 , the effective 1d hamiltonian given by eqs .
( [ hamil ] ) and ( [ rashba ] ) for any given nonzero @xmath32 does not preserve the tr symmetry , as its tr partner is at @xmath156 , making the @xmath4 index invalid .
the @xmath4 theory predicted that the tr symmetry is crucial for the topological spin pumping ,  @xcite suggesting that only the states at the tr - invariant point @xmath155 can contribute to the spin pumping .
this clearly contradicts the present result that all the states with @xmath157 contribute equally , which is obtained directly from the scattering matrix formula .
this point is also evidenced by the fact that the total amount of charge or spin pumped per cycle is in proportion to the sample width @xmath141 .
for the same reason , the pumping effect found in this work is also essentially different from that via edge states in refs .
@xcite , where the amount of spin or charge pumped per cycle is proportional to the number of the gapless edge channels .    in conclusion
, our work uncovers a bulk topological pumping effect due to direct transfer of the spwfs between the pump and electrode , without the participation of edge states .
this measurable effect reveals the bulk topological properties of the system that are neither captured by the @xmath3 index nor reflected by the number of gapless edge channels .
it can be accurately described by the spin chern numbers .
this spin chern pump may lay the foundation for direct experimental study and possibly utilization of the robust topological properties of the tis .    the previous experimental work  @xcite evidenced the difficulty of modulating in time the properties of an open quantum dot without generating undesired bias voltages due to stray capacitances .
this problem might not be significant in our pumping setup , where a much larger bulk sample of the ti can be used and the stray capacitances can be greatly reduced .
moreover , a possible way around the obstacle is to use the @xmath13 josephson effect to induce periodically time - dependent andreev reflection amplitudes in a hybrid normal - superconducting system .
@xcite concrete design of a spin chern pump based upon the josephson effect will await future work .
while the proposed spin pumping scheme may have the advantage of low noises , its practical application in spintronic devices still relies on the discovery of new tis with bulk band gaps much greater than room temperature , which determine the temperature range where the spin chern pumping effect can survive .
currently , precessing magnetization is a feasible method to generate robust spin currents in spintronic devices at room temperature .
@xcite    this work was supported by the state key program for basic researches of china under grants numbers 2015cb921202 , 2014cb921103 ( ls ) , 2011cb922103 and 2010cb923400 ( dyx ) , the national natural science foundation of china under grant numbers 11225420 ( ls ) , 11174125 , 91021003 ( dyx ) and a project funded by the papd of jiangsu higher education institutions .
we also thank the us nsf grants numbers dmr-0906816 and dmr-1205734 ( dns ) .
in what follows we expand on the model setup and possible experimental realization of the spin chern pump in more details .      a possible experimental realization of the hamiltonian eq .
( [ hamil ] ) for the pump is illustrated in fig .
a hgte / cdte quantum - well heterostructure with dual gates ( top and bottom ) is placed between two conductive plates .
it is known that when the width of the quantum well ( thickness of the hgte film ) is above a critical size @xmath158 ,  @xcite the band structure is inverted , characterized a negative mass term @xmath159 in the hamiltonian , corresponding to the qsh state .
if the width of the quantum well falls below @xmath160 , the band structure will be aligned in a `` normal '' way with a positive mass term @xmath161 , corresponding to a normal insulator . as has been discussed in refs .
@xcite and @xcite , the topological phase transition between the qsh phase and normal insulator can also be tuned by applying a gate voltage , which effectively reduces the width of the quantum well .
it is assumed that the quantum well under consideration has a width somewhat greater than @xmath160 , and so has a negative mass term @xmath159 initially . with increasing the gate voltage , the electron mass increases , and
can invert its sign .
usually , increasing the gate voltage may also adjust the carrier density .
nevertheless , it has been shown  @xcite that , by using dual gates and properly tuning their voltages @xmath162 and @xmath163 , it is generally possible to change the electron mass in the desired manner to be @xmath164 , while keeping the electron fermi energy still in the band gap .    the effect of the conductive plates is easily understood . when a voltage drop @xmath165 is applied across the plates , a uniform electric field @xmath166
will be generated in the space between the two plates .
the electrons in the quantum well experience a vector potential @xmath167 with @xmath168 defined as @xmath169 .
if the electric field is chosen to be @xmath170 , one gets @xmath12 with @xmath16 , as desired .
we point out that the exact time dependencies of @xmath22 and @xmath168 are not essential for realizing the spin chern pump , provided that they have the same periodicity and a constant relative phase shift .
+     direction.,width=163 ]      the pumping effect is insensitive to material details of the electrode .
the electrode is taken to be a normal metal with a 2d parabolic hamiltonian @xmath171 .
when @xmath172 is sufficiently large , for a given @xmath173 , we can linearize the effective 1d hamiltonian @xmath174 at the right and left fermi points @xmath175 with @xmath176 .
a pauli matrix @xmath177 is introduced to describe the two branches . to be consistent with the form of the hamiltonian in the pump , we use @xmath178 and @xmath179 , respectively , to represent the right - moving and left - moving branches for @xmath180 and oppositely for @xmath181 . as a result ,
the hamiltonian of the electrode becomes @xmath182 at @xmath183 , where @xmath184 and @xmath185 .
the spin pumping effect is usually dominated by small @xmath53 , so that we can further approximate @xmath186 , with purpose to minimize the number of adjustable parameters in the model .
for the present dirac - like hamiltonian , an ordinary potential barrier has a very weak effect on the electron transmission due to the klein paradox .
therefore , we take the hamiltonian for the barrier to be @xmath187 .
the inclusion of potential @xmath188 opens up an insulating energy gap of size @xmath52 around the fermi level , which presumably is more efficient for describing the contact deficiencies and structural mismatch between the pump and electrode .
the half metal , e.g. , cro@xmath189 , la@xmath190sr@xmath191mno@xmath192 , etc . , is a substance that acts as a conductor to electrons of one spin orientation , but as an insulator to those of the other spin orientation . from the viewpoint of the electronic structure , one of the spin subbands is metallic , whereas the fermi level falls into an energy gap of the other spin subband . to simulate the half - metallic electrode , @xmath174 is taken to be @xmath193 , where @xmath194 stands for the spin polarization of the electrode parallel and antiparallel to the @xmath5-axis , respectively .
the second term opens an energy gap of size @xmath195 around the fermi level for electron spin antiparallel to the spin polarization of the electrode , without affecting the other spin subband . as a result ,
the electron density of states is fully spin - polarized at the fermi energy .
@xmath196 is set to be infinity in the final result .
we now solve the scattering problem for an electron at the fermi energy incident from the electrode .
the fermi energy will be taken to be @xmath183 , which is in the band gap of the pump .
therefore , the incident electron will be fully reflected back into the electrode .
the rashba spin - orbit coupling is treated as a perturbation , and the result will be calculated to the linear order in the small quantity @xmath197 .
the wavefunctions of the pump ( @xmath198 ) , barrier ( @xmath199 ) and electrode ( @xmath200 ) are denoted by @xmath201 , @xmath202 , and @xmath203 , respectively .
we have two boundary conditions : @xmath204 and @xmath205 .
we use @xmath206 and @xmath207 to represent the eigenstates of @xmath69 , and @xmath208 and @xmath179 to represent those of @xmath209 . on the basis @xmath210 , @xmath211 , @xmath212 , and @xmath213 , the hamiltonian of pump ( the eqs .
( 1 ) and ( 2 ) in the manuscript ) can be expanded as a @xmath214 matrix    @xmath215    where @xmath216 . for energy @xmath217
, the eigen - equation is obtained from eq .
( [ hp_mat ] ) @xmath218 ^ 2-m^2(t)r^2_{0}\tilde{k}^2=0\ , \label{eigen_mat}\ ] ] with @xmath219 , and up to a normalization factor , the eigenfunctions are @xmath220 where @xmath221 and @xmath222/r_{0}$ ] .
we need to solve @xmath59 from the eigen - equation eq .
( [ eigen_mat ] ) .
we notice that the equation is a 4th - degree polynomial of @xmath59 with real coefficients , so complex conjugate roots must appear in pairs .
moreover , eq .  ( [ eigen_mat ] ) is even in @xmath59 , so positive and negative roots appear in pairs . in combination ,
( [ eigen_mat ] ) must have four roots of the form @xmath223 , @xmath224 , @xmath225 , and @xmath226 . by substitution of the four roots into eq .
( [ eigenfunc ] ) , we can in principle obtain four different eigenfunctions .
for the present scattering problem , we only need the two eigenfunctions that are decaying into the pump , which correspond to the two roots with negative imaginary parts .    for @xmath96 ,
it is easy to obtain for the roots for the two decaying modes : @xmath227 with @xmath228 , which are two - fold degenerate .
the corresponding two decaying eigenfunctions are given by @xmath229 @xmath230 where @xmath231 $ ] .
for @xmath232 , we write the roots of @xmath233 as @xmath234 , and also write @xmath235 as @xmath236 . to the second order in @xmath237 , we can solve for @xmath238 from the eigen - equation eq .
( [ eigen_mat ] ) @xmath239 noticing that the expression for @xmath240 given below eq .
( [ eigenfunc ] ) has a factor @xmath77 in the denominator , we keep @xmath238 to the second order , for the purpose to calculate @xmath240 to the linear order . by using the relation @xmath241
, we derive from eq .
( [ dkx2 ] ) @xmath242 with these relations , we obtain for @xmath243 and @xmath240 @xmath244\\ & \pm&i\frac{r_{0}}{2}(\hbar\eta)^2\sin^3(2\theta)+{\cal o}(\epsilon^2)\ , \end{aligned}\ ] ] and @xmath245 we can always eliminate any common factor that appears in all the four components of eq .
( [ eigenfunc ] ) , whenever possible . by eliminating a common factor @xmath246
, we rewrite @xmath243 and @xmath240 as @xmath247 and @xmath248 then the two decaying wavefunctions can be derived to be @xmath249 \end{array } \right)e^{\eta_{\mp}x}\ , \label{eigenfuncr0}\ ] ] where @xmath250 $ ] .
some remarks are in order . with respect to the wavefunctions at @xmath96 ,
namely , eqs .
( [ eigenfuncwithoutr00 ] ) and ( [ eigenfuncwithoutr01 ] ) , nonzero @xmath77 leads to nonperturbative change in the wavefunctions eq .
( [ eigenfuncr0 ] ) , in the sense that eq .
( [ eigenfuncr0 ] ) will not recover eqs .
( [ eigenfuncwithoutr00 ] ) and ( [ eigenfuncwithoutr01 ] ) in the limit @xmath251 .
this is reasonable , just as what always happens in the degenerate perturbation theory of the quantum mechanics . for the present problem ,
the wavefunction @xmath201 in the pump is always expressed as an arbitrary linear superposition of the two decaying modes : @xmath252 . by defining @xmath253 $ ] and @xmath254 $ ]
, we can rewrite @xmath201 as @xmath255 .
the final result for the reflection amplitudes depends only on @xmath256 , so we explicitly write out the expression for @xmath256 as follows @xmath257 where @xmath258 @xmath259 now we see that in the limit @xmath251 , the total wavefunction @xmath260 will go back to the form of a superposition of the two decaying wavefunctions given in eqs .
( [ eigenfuncwithoutr00 ] ) and ( [ eigenfuncwithoutr01 ] ) . in conclusion , while small rashba spin - orbit coupling may cause a nonperturbative change of the individual decaying wavefunctions , it modifies the `` space '' spanned by the two decaying modes in a perturbative manner .
it is this `` space '' which determines the final result of the reflection amplitudes .
this is the physical reason why in the final result , the rashba spin - orbit coupling modifies the reflection amplitudes in a perturbative manner .
the wavefunction in the potential barrier can be written as    @xmath261    where @xmath82 .
for an electron incident from the spin - up channel , the wavefunction in the electrode is given by @xmath262 first , matching the wavefunctions eqs .
( [ wavebarrier ] ) and ( [ waveelectrode ] ) at @xmath263 , one obtain @xmath264e^{-\gamma_{0}d}\ , \label{c_1}\\ c_{2}&=&\frac{1}{2}[(1+i)+r_{\uparrow\uparrow}(1-i)]e^{\gamma_{0}d}\ , \label{c_2}\\ c_{3}&=&\frac{r_{\downarrow\uparrow}}{2}(1+i)e^{-\gamma_{0}d}\ , \label{c_3}\\ c_{4}&=&\frac{r_{\downarrow\uparrow}}{2}(1-i)e^{\gamma_{0}d}\ .
\label{c_4}\end{aligned}\ ] ] in the next step , we will match wavefunctions at @xmath265 .
substituting eqs .
( [ c_1]-[c_4 ] ) into eq .
( [ wavebarrier ] ) , we can write eq .
( [ wavebarrier ] ) at @xmath266 as @xmath267 where @xmath268 . now equating eq .
( [ psim ] ) with eq .
( [ psip0 ] ) , we obtain @xmath269 it follows from eq .
( [ abcd0 ] ) @xmath78}{\mbox{ch}(2\gamma_{0}d)-\sin(2\theta ) \mbox{sh}(2\gamma_{0}d)}+{\cal o}(\epsilon^2)\ , \label{rupup}\ ] ] and @xmath79 } { \mbox{ch}(2\gamma_{0}d)-\sin(2\theta)\mbox{sh}(2\gamma_{0}d ) } + { \cal o}(\epsilon^2)\ . \label{rdnup}\ ] ]      the reflection amplitudes for an electron incident from the spin - down channel can be solved similarly .
now the wavefunction in the electrode is given by @xmath270 the forms of the wavefunctions in the pump and barrier remain to be the same . by some algebra , we arrive at @xmath271 equating eq .
( [ psim ] ) with eq .
( [ psip1 ] ) , we obtain @xmath272 it follows from eq .
( [ abcd1 ] ) @xmath273}{\mbox{ch}(2\gamma_{0}d)-\sin(2\theta ) \mbox{sh}(2\gamma_{0}d)}+{\cal o}(\epsilon^2)\ , \label{rdndn}\ ] ] and @xmath274 } { \mbox{ch}(2\gamma_{0}d)-\sin(2\theta)\mbox{sh}(2\gamma_{0}d ) } \left(\frac{r_{0}}{v_{\mbox{\tiny f}}}\right)+{\cal o}(\epsilon^2)\ . \label{rupdn}\ ] ]      the total hamiltonian of the system is invariant under the transformation @xmath275 so the corresponding transformation of eq .
( [ waveelectrode ] ) @xmath276 must also be an eigenstate of @xmath277 .
this result tells that an electron incident from the spin - down channel will be reflected into the spin - down channel with amplitude @xmath278 , and also into the spin - up channel with amplitude @xmath279 . comparing it with eq .
( [ waveelectrode1 ] ) and noticing that @xmath280 is equivalent to @xmath281 , we find immediately the following relations @xmath282 and @xmath283 the reflection amplitudes given in eqs .
( [ rupup ] ) , ( [ rdnup ] ) , ( [ rdndn ] ) , and ( [ rupdn ] ) apparently satisfy these relations .
to simulate the half - metallic electrode , the hamiltonian of the electrode is taken to be @xmath193 .
consider first the case , where the spin polarization of the electrode is parallel to the @xmath5-axis . now
the hamiltonian of the potential barrier becomes @xmath284 for simplicity , we will take @xmath285 limit in the final result . in this case , an electron incident from the spin - up channel has some probability to be reflected into the spin - down channel from the pump , but the reflected wave will decay quickly to @xmath62 within the barrier , and has no chance to reach the electrode .
therefore , @xmath145 and the wavefunction in the electrode becomes @xmath286 the wavefunction in the potential barrier can be written as      where @xmath82 and @xmath288 .
the wavefunction in the pump remains to be same .
repeating the same calculation as in sec .
ii , it is straightforward to obtain @xmath78}{\mbox{ch}(2\gamma_{0}d)-\sin(2\theta ) \mbox{sh}(2\gamma_{0}d)}+{\cal o}(\epsilon^2)\ . \label{rupup_down}\ ] ] this expression is identical in form to eq .
( [ rupup ] ) obtained in sec .
ii for a nonmagnetic electrode .
however , an important difference is @xmath145 . as a result ,
@xmath289 up to any order in @xmath237 .
this means that @xmath290 in eq .
( [ rupup_down ] ) must be a correction only to the argument of @xmath88 , which as discussed in the manuscript , will not modify the orbit of @xmath88 on the complex plane .
the orbit is always a unit circle for @xmath291 , and the amount of charge pumped per cycle by the @xmath53 state is quantized to @xmath292 .
similarly , for the spin polarization of the electrode antiparallel to the @xmath5-axis , the charge pumped per cycle is @xmath293 .
99 k. klitzing , g. dorda , and m. pepper , phys .
* 45 * , 494 ( 1980 ) .
a. tsukazaki , a. ohtomo , t. kita , y. ohno , h. ohno , and m. kawasaki , science * 315 * , 1388 ( 2007 ) .
y. zhang , y. w. tan , h. l. stormer , and p. kim , nature * 438 * , 201 - 204 ( 2005 ) . c. brne , c. x. liu , e. g. novik , e. m. hankiewicz , h. buhmann , y. l. chen , x. l. qi , z. x. shen , s. c. zhang , and l. w. molenkamp , phys .
. lett . * 106 * 126803 ( 2011 ) .
r. b. laughlin , phys .
b * 23 * , 5632 ( 1981 ) .
d. j. thouless , m. kohmoto , m. p. nightingale , and m. den nijs , phys . rev .
lett . * 49 * , 405 ( 1982 ) .
q. niu , and d. j. thouless , j. phys . a : math
* 17 * , 2453 ( 1984 ) .
d. j. thouless , phys . rev .
b * 27 * , 6083 ( 1983 ) .    c. l. kane , and e. j. mele , phys .
lett . * 95 * , 226801 ( 2005 ) .
b. a. bernevig , and s. c. zhang , phys .
lett . * 96 * , 106802 ( 2006 ) .
m. knig , s. wiedmann , c. brne , a. roth , h. buhmann , l. w. molenkamp , x. l. qi , and s. c. zhang , science * 318 * , 766 ( 2007 ) .
i. knez and r. r. du , frontiers of phys .
* 7 * , 200 ( 2012 ) .
j. e. moore and l. balents , phys .
b * 75 * , 121306 ( r ) ( 2007 ) .
l. fu and c. l. kane , phys .
b * 76 * , 045302 ( 2007 ) .
m. z. hasan , and c. l. kane , rev .
* 82 * , 3045 ( 2010 ) .
x. l. qi and s. c. zhang , physics today * 63 * , 33 ( 2010 ) .
f. d. m. haldane , phys .
* 61 * , 2015 ( 1988 ) . c. l. kane and e. j. mele , phys .
lett . * 95 * , 146802 ( 2005 ) .
d. n. sheng , z. y. weng , l. sheng , and f. d. m. haldane , phys .
lett . * 97 * , 036808 ( 2006 ) .
e. prodan , phys .
b * 80 * , 125327 ( 2009 ) .
y. yang , z. xu , l. sheng , b. g. wang , d. y. xing , and d. n. sheng , phys .
* 107 * , 066602 ( 2011 ) .    x. l. qi , t. l. hughes , and s. c. zhang , phys
b * 78 * , 195424 ( 2008 ) .
a. m. essin , j. e. moore , and d. vanderbilt , phys .
* 102 * , 146805 ( 2009 ) .
l. fu , and c. l. kane , phys .
b * 74 * , 195312 ( 2006 ) . c. q. zhou , y. f. zhang , l. sheng , r. shen , d. n. sheng , and d. y. xing , phys . rev .
b * 90 * , 085133 ( 2014 ) .
x. l. qi , t. l. hughes , and s. c. zhang , nature physics * 4 * , 273 ( 2008 ) .
f. mahfouzi , b. k. nikoli , s. h. chen , and c. r. chang , phys .
b * 82 * , 195440 ( 2010 ) .
b. a. bernevig , t. l. hughes , and s. c. zhang , science * 314 * , 1757 - 1761 ( 2006 ) .
m. knig , h. buhmann , l. w. molenkamp , t. l. hughes , c. x. liu , x. l. qi , and s. c. zhang , j. phys .
* 77 * , 031007 ( 2008 ) .
w. yang , k. chang , and s. c. zhang , phys .
* 100 * , 056602 ( 2008 ) . c. x. liu ,
t. l. hughes , x. l. qi , k. wang , and s. c. zhang , phys .
lett . * 100 * , 236601 ( 2008 ) .
d. g. rothe , r. w. reinthaler , c. x. liu , l. w. molenkamp , s. c. zhang , and e. m. hankiewicz , new j. phys .
* 12 * , 065012 ( 2010 ) .
g. h. wannier , rev .
mod . phys . * 34 * , 645 ( 1962 ) .
x. l. qi , phys .
lett . * 107 * , 126803 ( 2011 ) .
h. li , l. sheng , and d. y. xing , phys . rev .
lett . * 108 * , 196806 ( 2013 ) .
m. bttiker , h. thomas , and a. prtre , z. phys .
b * 94 * , 133 ( 1994 ) .
p. w. brouwer , phys .
b * 58 * , r10135 ( 1998 ) .
m. switkes , c. m. marcus , k. campman , and a. c. gossard , science * 283 * , 1905 ( 1999 ) .
f. giazotto , p. spathis , s. roddaro , s. biswas , f. taddei , m. governale , and l. sorba , nature physics * 7 * , 857 ( 2011 ) . g. e. w. bauer , and y. tserkovnyak , physics * 4 * , 40 ( 2011 ) .