the central idea behind the field of spintronics is to use the electron s spin degree of freedom , not its charge , for information storage and processing.@xcite since the energy required to generate magnetic fields scales unfavorably at small length scales , electrical mechanisms for the manipulation and detection of magnetic moments are crucial for successful spintronics applications .
the spin transfer torque@xcite has been established as a reliable effect to manipulate the magnetization of a thin ferromagnetic layer with a spin - polarized current.@xcite whereas injection from a ferromagnet was used as a spin - polarized current source in the original realization , devices based on spin - orbit coupling have also been realized.@xcite    an extreme form of spin - orbit coupling exists in the surface states of topological insulators ( tis):@xcite these states have complete spin - momentum locking , _
i.e. _ , the electron spin and its direction of motion are perfectly correlated .
one consequence of the spin - momentum locking is that a surface charge current in a ti is automatically spin - polarized .
recent proposals have concentrated on exploiting this effect to control the dynamics of thin magnetic layers in the proximity of the surface of a three - dimensional ti@xcite or magnets coupled to the edge of a two - dimensional ti.@xcite experiments on metallic magnetic layers in contact to a ti have reported a spin - transfer torque exceeding the values found in non - topological spin - orbit materials.@xcite    in this paper , we investigate the possibility to use the spin - polarized surface currents of a topological insulator to control the magnetic moment of a molecular magnet adsorbed on the surface of a three - dimensional topological insulator or at the edge of a two - dimensional topological insulator . like a ferromagnet ,
a molecular magnet has degenerate magnetic ground states , separated by an energy barrier,@xcite although in the case of a molecular magnet the barrier is microscopic , not macroscopic , which leads to a finite relaxation time of the molecule s magnetic moment .
of particular interest are the so - called `` single - molecule magnets'',@xcite which consist of a magnetic core with a spin @xmath0 , shielded from the environment by a ( typically ) organic ligand .
for single - molecule magnets magnetic lifetimes of several years have been reported for temperatures below the `` blocking temperature '' set by the anisotropy barrier between the two magnetic ground states.@xcite single - molecule magnets have been shown to preserve their magnetic properties , including long magnetic lifetimes , when adsorbed on conducting surfaces.@xcite    the mechanism by which an applied electrical current at the surface of a ti can be used to switch the orientation of a molecular magnet is best illustrated using the example of a localized spin @xmath1 coupled to the edge of a two - dimenstional topological insulator.@xcite the electronic state at the edge of a two - dimensional ti is helical , so that electrons moving in opposite directions have opposite spin.@xcite hence , a particle current to the right is polarized as `` spin up '' , see fig .
[ fig : sketch ] . by spin conservation , backscattering of a right - moving electron into a left - moving state requires the localized spin to flip from down to up , so that a single backscattering event is sufficient to polarize the spin @xmath1 in the `` up '' state .
if the applied current is sufficiently large , the pauli principle forbids backscattering of left - moving electrons  which would be accompanied by a spin flip in the opposite direction  , so that the spin @xmath1 remains in the `` up '' state as long as the current is being applied .
weakly coupled to the helical edge of a two - dimensional topological insulator .
backscattering of a right - moving electron is accompanied by a flip to the spin `` up '' state of the localized spin , leaving it in a fully polarized state after a single backscattering event . for a sufficiently large applied current
the reverse scattering process , which would return the localized spin to the `` down '' state , is strongly suppressed.,scaledwidth=45.0% ]    this simple picture needs to be refined for molecules with a higher spin and for molecular magnets adsorbed to the two - dimensional surface of a three - dimensional ti . unlike the spin @xmath1 of the example of fig .
[ fig : sketch ] , a higher - spin molecule comes with its own anisotropy axis , and the argument based on spin conservation no longer applies in this simple form if the molecule s anisotropy axis is not aligned with the spin quantization axis of the ti edge state .
simple spin conservation arguments can not be applied to a molecular magnet on the surface of a three - dimensional topological insulator either , because in that case there is no unique spin quantization axis , as electrons can propagate in an arbitrary direction along the surface . as we will show below , in this generic situation current - induced switching of molecular magnet on a ti surface
is not perfect , although appreciable polarizations can be achieved even for randomly oriented molecules .
an interesting observation is that , unlike in the example above , for a molecular magnet on the surface of a three - dimensional ti the current - induced polarization is not a monotonously increasing function of the applied current , but has a maximum at intermediate current densities .
this article is organized as follows : in sec .
[ sec:2dti ] we consider a molecular magnet at the edge of a two - dimensional ti , expanding the example of fig .
[ fig : sketch ] to higher - spin molecules .
we calculate the probability with which current - induced switching takes place , the current - induced switching rate , and the zero - current relaxation rate that arises from the exchange coupling of the molecule s spin to the ti edge . in sec .
[ sec:3dti ] we consider the same questions for a molecular magnet on the surface of a three - dimensional topological insulator , which is the more likely candidate for an experimental realization .
we conclude with a brief outlook in section [ sec : conclusion ] .
the main text mainly focuses on molecules with integer spin @xmath2 .
the case of half - integer spin , which requires a technically more demanding analysis , is discussed in the appendix .
to introduce our notation and to provide a reference for further calculations , we start by considering a localized moment of spin @xmath3 exchange - coupled to the helical edge of a two - dimensional topological insulator .
we assume that there is only a single localized spin coupled to the edge and that the temperature is large enough that the kondo effect can be neglected . in that case transitions between different spin states
can be described using rate equations .
the coupling between the edge and the magnetic moment is described by the exchange hamiltonian @xmath4 where @xmath5 labels the coordinate along the helical edge , @xmath6 is the velocity of the helical edge state , and @xmath7 is the dimensionless exchange coupling . because of spin - momentum locking of the edge states , backscattering of edge electrons from the molecular magnet involves simultaneous flips of the spin of the edge electrons and of the localized spin of the molecular magnet .
we fix the spin @xmath5 axis along the direction of the spin quantization axis for the edge - state electrons , @xmath8 so that right - moving electrons at the edge have spin up and left - moving electrons have spin down .
we note that the spin @xmath5 axis need not be in the same direction as the laboratory @xmath5 axis .
( in fact , for the edge of a two - dimensional topological insulator , the spin quantization axis is commonly taken to be perpendicular to the plane of the ti.@xcite ) the states @xmath9 in the helical edge are labeled by their energy @xmath10 and the propagation direction @xmath11 , where we take @xmath12 for right - moving electrons ( positive @xmath5 direction ) and @xmath13 for left - moving electrons .
( note that in this notation @xmath14 merely parameterizes the energy ; its magnitude @xmath15 equals the magnitude of the momentum , but the sign of @xmath14 is that of the energy @xmath16 , not of the momentum . )    the transition rates @xmath17 between the spin states are calculated from fermi s golden rule .
the rate @xmath18 for transitions from the spin up state @xmath19 to the spin down state @xmath20 is @xmath21    \nonumber \\ & & \mbox { } \times    |\langle + , \downarrow|h_{\rm ex}| -,\uparrow \rangle|^2     \delta(\varepsilon_k - \varepsilon_{k ' } )    \nonumber \\ & = &    \frac{j^2}{2 \pi \hbar } \int\limits_{-\infty}^{+\infty } d \varepsilon    n_{-}(\varepsilon ) ( 1 - n_{+}(\varepsilon ) ) .
\label{eq : gamma1}\end{aligned}\ ] ] ( the matrix elements of @xmath22 do not depend on @xmath14 and @xmath23 , which is why we have suppressed @xmath14 , @xmath23 in our notation . ) similarly , the rate @xmath24 at which transitions from the spin down state to the spin up state take place reads @xmath25    the distribution function at temperature @xmath26 and chemical potential @xmath27 in the presence of an edge current @xmath28 is @xmath29 as one can easily verify from the relation @xmath30 $ ] .
for definiteness , we will assume throughout that the particle current is to the right , @xmath31 .
performing the integrations in eqs .
( [ eq : gamma1 ] ) and ( [ eq : gamma2 ] ) then gives @xmath32 where we abbreviated @xmath33 for high applied currents , @xmath34 , scattering processes in which a left - moving electron is scattered into a right - moving one and , correspondingly , the rate @xmath18 are exponentially suppressed .
on the other hand , for low applied currents @xmath35 , scattering of thermally excited charge carriers dominates the spin switching rates , and the application of the current @xmath36 only gives a slight asymmetry .    in the absence of coherences between
the spin states @xmath37 and @xmath38 the probabilities @xmath39 and @xmath40 to find the spin a state of spin @xmath41 can be solved from the stationary solutions of the rate equation @xmath42 ( in principle , one needs a density matrix to describe possible coherent superpositions of the degenerate states @xmath37 and @xmath38 , see , _
e.g. _ , refs .   and the appendix .
however , such coherences do not occur in the present case , in which there is no external magnetic field and a single spin quantization axis for the conduction electrons at the ti edge and the localized spin . ) substituting the rates ( [ eq : gamma12 ] ) one immediately finds @xmath43 the same result also follows from the observation that the coupling to the conduction electrons at the ti edge effectively amounts to a zeeman shift @xmath44 between the spin up and spin down states of the molecule .
we conclude that the application of a current @xmath45 leads to a complete polarization of the molecule s spin .
although the application of a current through the helical edge of a two - dimensional topological insulator causes a polarization of a spin @xmath1 exchange - coupled to the edge state , the induced magnetic moment quickly disappears as soon as the current is switched off . in our model ,
in which the coupling to the helical edge state is the only source of relaxation , the corresponding relaxation rate @xmath46 is given by the low - current limit of eq .
( [ eq : gamma12 ] ) , @xmath47     ( a ) for a spin-@xmath48 molecule and ( b ) for a spin-@xmath49 molecule.,scaledwidth=45.0% ]    longer relaxation times in the absence of an applied current requires molecules with a higher spin @xmath2 . in this case , magnetic anisotropy , @xmath50 imposes an energy barrier between states with maximal and minimal spin @xmath51 ( measured along the anisotropy axis @xmath52 ) , see fig .
[ fig : sketch_spin1 ] , which , if @xmath53 is lower than the barrier energy , leads to strongly enhanced lifetimes .
( see below for a quantitative estimate for our model . ) in our analysis , we have adopted a simple easy - axis anisotropy strength @xmath54 ( which we take to have the units of frequency ) with @xmath55 , although our considerations also carry over to other forms for the anisotropy energy in which there are two different minima .
( the case @xmath56 is for easy - plane anisotropy , which does not have two separated energy minima and will not be considered here . )
we label the spin states with the quantum number @xmath57 , @xmath58 ,  , @xmath2 for the spin component @xmath59 along the @xmath52-direction .
( for @xmath3 , @xmath60 takes the values @xmath61 , corresponding to the notation `` @xmath62 '' and `` @xmath63 '' in the previous subsection . )    if the anisotropy axis is aligned with the spin quantization axis for the helical edge state , the analysis of the previous subsection immediately carries over , and one finds @xmath64 .
\label{eq : gammass}\ ] ] for the probabilities @xmath65 one then finds @xmath66 where the prefactor is fixed by the normalization condition @xmath67 .    if the temperature is much smaller than the anisotropy energy , which is a condition that we will assume throughout this article , the rate equations can be solved directly , using detailed balance , to give the zero - current switching rate @xmath68 between the two magnetic ground states @xmath69 and @xmath70 , @xmath71 if @xmath2 is integer and @xmath72 if @xmath2 is half integer .
the application of a current @xmath36 slightly larger than @xmath73 leads to a quick and complete spin polarization of the molecule .
the switching rate @xmath74 can be estimated as the inverse of the sum of inverse transition rates @xmath75 , which gives @xmath76 for large @xmath2 .
a smaller current @xmath77 also polarizes the molecule , but since thermal activation is still needed in the switching process , the time required to reach the polarized state is still large ( though much shorter than the zero - current lifetime ) , @xmath78 if @xmath2 is integer and @xmath79 if @xmath2 is half integer .    if @xmath52 is not aligned with the @xmath5 axis , the expression for the rates @xmath80 involves the matrix elements    @xmath81    where @xmath82 and @xmath83 are the polar angles corresponding to the anisotropy axis @xmath84 .
this gives @xmath85 \sin^2 \theta    + ( 1 \mp \tau \cos \theta)^2     g[\tau { i}\mp d ( 2 s \pm 1 ) ]    \right\}. \label{eq : gammaanis}\end{aligned}\ ] ] one verifies that the zero - current rates are the same as those of eqs .
( [ eq : gammass ] ) , independent of the direction @xmath52 of the anisotropy axis , so that eqs .
( [ eq : gammaswitch0even ] ) and ( [ eq : gammaswitch0odd ] ) still apply .    because of
the degeneracy of the states @xmath37 and @xmath38 the rate equation ( [ eq : master_time - dep ] ) is not sufficient to describe transitions between these two states and , instead , one has to use a master equation for the reduced @xmath86 density matrix that is able to account for coherent superpositions of the states @xmath37 and @xmath38 .
for that reason we here restrict ourselves to the case of integer @xmath2 , for which this complication does not occur . a complete analysis and final expressions for half - integer @xmath2
are given in the appendix .
we first discuss the case of an `` intermediate '' current @xmath87 , for which the transition rates ( [ eq : gammaanis ] ) can be approximated as @xmath88    the large difference between `` upstream '' transitions to higher - energy spin states and `` downstream '' transitions to lower - energy spin states ensures that the molecule will predominantly be in one of its two magnetic ground states @xmath69 and @xmath70 . the presence of the applied current breaks the symmetry between these two ground states and leads to a preferred population of the state @xmath69 if @xmath89 and @xmath70 if @xmath90 , although there is no longer a perfect polarization for generic @xmath82 , @xmath91 switching between the magnetic ground states involves thermal activation , which is why the switching rate between the magnetic ground states is exponentially long in @xmath92 , although the rate is parametricaly larger than the spontaneous relaxation rate at zero current . defining @xmath74 as the switching rate from @xmath70 to @xmath69 we have
@xmath93 where @xmath68 is the zero - current relaxation rate of eq .
( [ eq : gammaswitch0even ] ) .
the switching rate for the inverse process is obtained upon replacing @xmath94 .    ) for @xmath95 and for angles @xmath96 , @xmath97 , @xmath98 , @xmath99 , @xmath100 , @xmath101 , @xmath102 , @xmath103 and @xmath104.,scaledwidth=45.0% ]    the difference between the case @xmath96 and generic @xmath82 is most pronounced in the `` high - current '' limit @xmath105 , in which the rates @xmath80 in eq .
( [ eq : gammaanis ] ) are dominated by the terms proportional to @xmath106 $ ] . approximating @xmath107 \approx { i}$ ] for @xmath108 , we find @xmath109 for the probabilities @xmath65 this gives immediately @xmath110 for @xmath82 close to zero this still gives a narrow distribution around @xmath111 , but perfect polarization is not reached if @xmath112 independent of the magnitude of the current @xmath36 .
examples of the distribution ( [ eq : ps2d ] ) are shown in fig .
[ fig : ps2d ] for @xmath95 and for several representative angles @xmath82 .    figure [ fig : mean2d ] shows the mean values @xmath113 of the spin component along the anisotropy axis and @xmath114 of the spin component along the current direction as a function of the angle @xmath82 for two different values of the total spin @xmath2 .
except for @xmath100 , the application of a current always results in a net spin in the current direction . for a randomly oriented molecule ,
the average moment in the current direction approaches @xmath115 in the limit of large @xmath2 . without a perfect current - induced polarization
, we should define the current - induced `` switching rate '' @xmath74 as the rate at which the steady - state probability @xmath65 is approached , starting from @xmath70 ( if @xmath89 ) .
the estimate ( [ eq : gammaswitch2dhigh ] ) for the switching rate @xmath74 obtained above for @xmath52 parallel to the @xmath5 axis remains valid as a good order - of - magnitude estimate for the case of a general orientation of the anisotropy axis in the high - current regime .
along the direction of the anisotropy axis ( solid curves ) and the mean spin component @xmath116 in the current direction ( dashed ) , as a function of the angle @xmath82 , for @xmath117 ( left panel ) and @xmath118 ( right panel).,scaledwidth=45.0% ]    the main result of the previous section is that the current - induced polarization of a molecular magnet coupled to the edge of a two - dimensional ti is complete only if the molecule s anisotropy axis is aligned with the spin quantization axis for the ti edge state , whereas the polarization is incomplete  but generally nonzero  for arbitrary orientations of the molecule .
the analysis was simplified by the fact that the edge of a two - dimensional ti has a unique propagation direction and , hence , a unique spin quantization axis .
this is the main difference with the case of a molecule coupled to the surface of a three - dimensional ti : electrons at the surface of a three - dimensional ti can propagate in all directions in the plane of the surface , even in the presence of a large applied current . hence , there is not a unique spin quantization axis , and there will always be surface - state electrons with a spin that is not aligned with the anisotropy axis of the molecular magnet .
nevertheless , as we show below , even for a molecule adsorbed on the surface of a three - dimensional ti the application of an electrical current can lead to a complete spin polarization , provided the molecule s anisotropy axis is aligned with the current direction and the magnitude of the applied current is appropriately chosen . as in the case of a two - dimensional ti ,
the spin polarization of the molecular magnet is not complete ( but generally nonzero ) if the molecule s magnetic anisotropy axis is not aligned with the current direction
.     plane .
a current is applied the positive @xmath5 direction.,scaledwidth=45.0% ]    to keep the notation close to that of the previous section , we take the ti surface to be the @xmath119 plane and take the direction of ( particle ) current flow to be the positive @xmath5 direction , see fig .
[ fig : sketch3d ] .
the electrons at the ti surface are then described by the hamiltonian @xmath120 we write the current density @xmath121 as @xmath122 where @xmath123 is the fermi wavenumber and @xmath36 has the dimension of current . for the exchange hamiltonian
we take @xmath124 where , as in the previous section , a prefactor has been included to make the strength @xmath7 of the exchange interaction dimensionless .
the hamiltonian @xmath125 for the magnetic anisotropy of the magnetic molecule is the same as in the previous subsection , see eq .
( [ eq : hs ] ) .
the eigenstates of @xmath126 are labeled by the energy @xmath10 and by the angle @xmath127 of the propagation direction in the @xmath119 plane , such that @xmath128 corresponds to the positive @xmath5 direction .
we assume that the chemical potential @xmath129 , which ensures that thermally excited carriers as well as carriers contributing to the current flow remain well away from the dirac point . in that case
the distribution function in the presence of a particle current density @xmath130 is @xmath131 equation ( [ eq : n3d ] ) both follows from a solution of the boltzmann equation and maximizes the entropy under the constraint of a fixed current density @xmath132 .    for the calculation of the transition rates between different spin states we need the matrix elements    @xmath133 ,     \label{eq : matrixelements3d}\end{aligned}\ ] ]    where we abbreviated @xmath134 .
equation ( [ eq : matrixelements3d ] ) generalizes eq .
( [ eq : matrixelements2d ] ) to the case of a surface state . as in eq .
( [ eq : matrixelements2d ] ) , @xmath82 and @xmath83 are the polar angles marking the direction of the anisotropy axis . for the transition rates
we then find , again assuming @xmath135 , @xmath136 ,    \label{eq : gammageneral}\end{aligned}\ ] ] where the function @xmath137 was introduced in eq .
( [ eq : g ] ) .
we will now analyze these rates and the resulting probabilities @xmath65 in the regimes of low current densities , @xmath138 , intermediate current densities , @xmath139 , and high current densities , @xmath140 .
( for @xmath3 no intermediate regime exists and the high - current regime is defined as @xmath34 . )    for low current densities @xmath138 the expressions ( [ eq : gammageneral ] ) for the transition rates simplify to @xmath141 the resulting zero - current switching rate @xmath68 is given by the same expression as in the two - dimensional case , see eqs .
( [ eq : gammaswitch0even ] ) and ( [ eq : gammaswitch0odd ] ) .    for intermediate current densities , @xmath142 , the exponential suppression of the transition rates to higher - energy spin states
is reduced , whereas the transition rates into the lower - energy states remains approximately as in the low - current regime . for a quantitive analysis
we restrict ourselves to the case of integer spin @xmath2 , referring to the appendix for a discussion of half - integer @xmath2 . in the intermediate - current regime
the transition rates can be approximated as @xmath143    as in the two - dimensional case , the applied current breaks the symmetry between states with positive and negative spin @xmath60 , leading to an imbalance between the populations of the ground states @xmath69 and @xmath70 , @xmath144 the population of the excited states @xmath145 remains exponentially small in @xmath146 .
the current - induced switching rate from the state @xmath70 into @xmath69 is @xmath147^s    \nonumber \\ & & \mbox { } \times    e^{2 \pi \hbar s { i}/k_{\rm b } t}.    \label{eq : gammaswitchj3deven}\end{aligned}\ ] ] the inverse rate is given by the same expression , but with @xmath148 replaced by @xmath149 .
the case @xmath150 that the molecule s anisotropy axis is aligned with the current direction requires to be considered separately . in this case
the leading approximation of eq .
( [ eq : gamma3dintermediate ] ) vanishes for the `` upstream '' rate @xmath151 with @xmath152 and subleading contributions in the small parameter @xmath153 must be taken into account .
one finds @xmath154 if @xmath152 .
asymptotically , for @xmath34 [ but still @xmath155 , this leads to a complete spin polarization of the molecule , @xmath156 [ the probability @xmath157 vanishes @xmath158 ; all other probabilities are exponentially small in @xmath146 . ]
the perfect polarization can be understood if one observes that the surface - state electrons responsible for the transitions to higher - energy spin states are predominantly electrons moving in the positive @xmath5 direction , which are then backscattered to the negative @xmath5 direction .
these electrons have a unique spin quantization axis , and their scattering leads to a well - defined change of the molecule s spin state @xmath159 .
the switching rate from the state @xmath70 to @xmath69 is given by eq .
( [ eq : gammaswitchj3deven ] ) with @xmath150 . the rate for the inverse process @xmath160 vanishes in the limit @xmath161 .    , @xmath162 , and @xmath163 for a spin-@xmath48 molecule on the surface of a three - dimensional topological insulator to be in the corresponding spin state as a function of @xmath164
we have set @xmath165 .
the inset shows the maximum value of the probability @xmath166 for a molecule with @xmath117 as a function of @xmath167.,scaledwidth=45.0% ]    remarkably , the complete polarization is lost again when the current is increased further . if @xmath105 , not only electrons moving in the positive @xmath5 direction but essentially all electrons that contribute to the current @xmath121 can scatter off the localized magnetic moment and change its spin . for high current densities
the transition rates are ( for arbitrary direction of the anisotropy axis ) @xmath168
where we abbreviated @xmath169 and one quickly obtains the steady - state distribution from here , @xmath170
f_+(\theta)^{s+s } f_-(\theta)^{s - s } }    { f_+(\theta)^{2 s + 1 } - f_-(\theta)^{2 s + 1}}.    \label{eq : ps3dhigh}\ ] ] since @xmath171 , even for a perfectly aligned molecule no perfect polarization occurs in the limit of large currents .
the probability @xmath172 to find the molecule in the maximal - spin state saturates at a value slightly above @xmath173 for molecules with large spin @xmath2 .
although no perfect polarization results for current densities @xmath105 , the rate at which the steady - state distribution is approached is considerably enhanced in comparison to the intermediate current - density regime , see eq .
( [ eq : gammaswitch2dhigh ] ) .    as an illustration
, we show the probabilities @xmath65 for @xmath174 for a molecule with @xmath117 as a function of the applied current density for @xmath175 , see fig .
[ fig : spin1 ] .
notice the ( faint ) maximum of @xmath176 for @xmath177 .
the inset shows that the maximum value of @xmath176 approaches unity in the limit when the anisotropy energy is much larger than the temperature , @xmath178 .
figure [ fig : mean3d ] shows the mean spin components @xmath179 along the molecule s anisotropy axis and @xmath116 along the current direction as a function of the angle @xmath82 for @xmath95 , for the intermediate - current regime and for the high - current regime . as in the case of a two - dimensional ti , except for @xmath100 , the application of a current always results in a net spin in the current direction .
for a randomly oriented molecule , the average moment in the current direction approaches @xmath115 in the limit of large @xmath2 .     along the direction of the anisotropy axis ( solid curve ) and
the mean spin component @xmath116 in the current direction ( dashed ) , as a function of the angle @xmath82 , for total spin @xmath95 . left panel : intermediate - current regime @xmath142 ; right panel : high - current regime @xmath180.,scaledwidth=45.0% ]    we can apply the results obtained here to a topological insulator surface with a dilute covering of molecular magnets . for a dilute covering our single - molecule analysis can be used to obtain the net polarization of the molecular layer . assuming that the molecules have randomly orientated anisotropy axis
 a valid assumption if the magnetic core is shielded by an approximately spherical shell , see , _
e.g. _ , ref .
 the net polarization is found by averaging a single molecule s average moment over the directons of the anisotropy axis @xmath52 .
the calculation is straightforward in principle , starting from the steady - state distributions ( [ eq : p3deven ] ) and ( [ eq : ps3dhigh ] ) .
[ see eqs .
( [ eq : p3dodd ] ) and ( [ eq : p3dhighodd ] ) for steady - state distributions for half - integer spin @xmath2 . ] instead of reporting the resulting expressions , which are rather cumbersome , we refer to fig .
[ fig : average ] , which shows the resulting net polarization per molecule as a function of @xmath2 .
per molecule as a function of the total spin @xmath2 , for an ensemble of randomly oriented molecular magnets on the surface of a three - dimensional topological insulator .
open circles are for the intermediate - current regime @xmath181 ; filled circles are for the high - current regime @xmath182.,scaledwidth=45.0% ]
the strict spin - momentum locking at the surface of a topological insulator provides an appealing mechanism to use electric currents to switch the spin state of a molecular magnet weakly coupled to the surface . whereas based on simple spin conservation arguments one expects a complete response for a molecular magnet at the edge of a two - dimensional topological insulator if the molecule s magnetic anisotropy axis is aligned appropriately with the spin quantization axis of the surface state , the situation is more complicated for molecular magnets on the surface of a three - dimensional topological insulator and/or for an arbitrary orientation of the magnetic anisotropy axis . for a molecule on the surface of a three - dimensional topological insulator
, a full current - induced polarization is achieved only if the molecule s magnetic anisotropy axis is aligned with the current direction and the magnitude @xmath183 of the current density is in an `` intermediate '' range @xmath184 , with @xmath185 the molecule s anisotropy energy and @xmath2 the magnitude of its spin . yet , even for a surface covered by molecules with randomly oriented anisotropy axes the application of a surface current should result in a net magnetic moment per molecule for all current strengths .
+   + the minimal current density @xmath183 required to polarize the molecule and the resulting steady - state distribution of the molecule s magnetic states depends more on properties of the topological insulator surface than properties of the molecular magnet itself .
the condition for an appreciable current - induced polarization is @xmath186 , _
i.e. _ , the disortion of the fermi surface by the applied current should be large in comparison to the blurring of the fermi surface by temperature .
this is a condition that is easier to achieve if the chemical potential is in the vicinity of the dirac point , so that @xmath123 is small , and a relatively small current density implies a large shift of the distribution function .
the role of the anisotropy energy is to set a crossover scale @xmath187 , above which current - induced transitions allow the molecule to be in all spin states , whereas for lower current densities only the magnetic ground states @xmath69 and @xmath70 are accessible .
again , the crossover between these regimes happens at lower current densities if the chemical potential is closer to the dirac point .
+   + details of the coupling between the molecule and the surface play a crucial role when it comes to setting the switching rate .
this applies both to the zero - current spontaneous switching between the degenerate magnetic ground states @xmath69 and @xmath70 and to the current - induced changes between these states .
if the temperature is below the blocking temperature , appreciable switching rates can be obtained only in the `` high - current regime '' @xmath188 , where we found the switching rate to be of order @xmath189 , where @xmath7 is the dimensionless exchange coupling . an upper limit for this rate
is found by setting @xmath190 , although we expect that the presence of the organic shell will typically lead to much smaller exchange couplings for single - molecule magnets on a topological insulator surface .
taking a current density @xmath183 at the lower end of the high - current regime , @xmath187 , we find an upper bound for switching times in the range of @xmath191 s for a molecule with a blocking temperature of a few k , but most likely significantly longer if @xmath7 is smaller .
+   + our calculations have shown that in - plane currents can be used as an effective tool to switch the orientation of the spin of molecular magnets adsorbed on the surface of a three - dimensional topological insulator , even if no complete polarization is achieved for molecules without a specially aligned magnetic anisotropy axis .
we hope that our observations , together with other theoretical studies highlighting the intriguing physics of molecular magnets on ti surfaces,@xcite will motivate further experimental and theoretical research in this direction .
we gratefully aknowledge support from the german research foundation ( dfg ) in the framework of the collaborative research center 658 .
in the case of a molecular magnet with half - integer spin @xmath2 the states @xmath37 and @xmath38 are degenerate , which means that transitions between these states can not be described by means of a rate equation , but that they need to be described with a master equation for the ( reduced ) density matrix @xmath192 the diagonal elements @xmath193 are the probabilities to find the molecule in spin state @xmath159 , @xmath194 ; the off - diagonal elements describe coherences beween these two spin states .    to allow for a unified treatment of a molecular magnet coupled to the edge of a two - dimensional ti and a molecular magnet at the surface of a three - dimensional ti
, we write matrix elements of the exchange hamiltonian ( [ eq : hex2d ] ) for the two - dimensional case as @xmath195 where we use the propagation angle @xmath128 for a right - moving electron ( labeled with @xmath196 in sec .
[ sec:2dti ] ) and @xmath197 for a left - moving electron ( labeled with @xmath198 ) .
matrix elements of the exchange hamiltonian ( [ eq : hex3d ] ) for the three - dimensional case are written @xmath199 the rates @xmath200 , see eqs .
( [ eq : gammaanis ] ) and ( [ eq : gammageneral ] ) , are then expressed as      where the double integration @xmath202 should be replaced by a double sum @xmath203 for the case of a two - dimensional ti .
the explicit expressions for the @xmath204 read @xmath205 where @xmath134 , see the text below eq .
( [ eq : matrixelements3d ] ) .
because of the large energy differences between spin states with different values of @xmath206 and the lack of direct transitions between the degenerate spin states @xmath159 and @xmath207 for @xmath208 , the rates ( [ eq : gammassapp ] ) are sufficient to calculate the ( rate of change of the ) probabilities @xmath65 that the molecule is in spin state @xmath159 for @xmath208 , see eq .
( [ eq : master_time - dep ] ) .
on the other hand , for the degenerate spin states @xmath37 and @xmath38 we must use the full density matrix @xmath209 , see eq .
( [ eq : rhodef ] ) .
the equation of motion for this density matrix has the well - known lindblad form@xcite @xmath210    \label{eq : lindblad }    \\ & & \mbox { }    + j^2 \int \frac{d\varphi}{2 \pi } \frac{d\varphi'}{2 \pi }    g[{i}(\cos \varphi'-\cos\varphi)/2 ]    \left\ { m(\varphi',\varphi ) \rho m(\varphi',\varphi)^{\dagger }    - \frac{1}{2 } \left [    m(\varphi',\varphi)^{\dagger } m(\varphi',\varphi ) \rho    +    \rho m(\varphi',\varphi)^{\dagger } m(\varphi',\varphi ) \right ]    \right\ } , \nonumber\end{aligned}\ ] ]    where the projection matrices @xmath211 , and @xmath212 are defined as @xmath213 one verifies that the lindblad equation ( [ eq : lindblad ] ) simplifies to the standard rate equation ( [ eq : master_time - dep ] ) if the density matrix @xmath209 is a diagonal matrix , @xmath214 .    the first term in eq .
( [ eq : lindblad ] ) contains the expectation value of the exchange hamiltonian projected onto the spin states @xmath215 .
it acts as an effective magnetic field driving coherent transitions between the states @xmath37 and @xmath38 .
this term plays no role for states @xmath159 with @xmath208 because it is much smaller than the anisotropy field if the exchange coupling is weak , @xmath216 .
however , for transitions between the states @xmath37 and @xmath38 there is no anisotropy field , and the exchange field competes with the similarly small scattering - induced transitions between the states @xmath37 and @xmath38 . with the same notation as above we have @xmath217 where , for a molecule coupled to the edge of a two - dimensional topological insulator , the integration @xmath218 should be replaced by a summation @xmath219 . performing the angular integration
gives @xmath220 which is indeed what one expects for an exchange field in the @xmath5 direction , recalling that the spin states @xmath159 are defined with respect to the anisotropy axis @xmath84 .
a solution of the lindblad equation proceeds by parameterizing @xmath221 where @xmath222 is the @xmath86 identity matrix and the @xmath223 are the pauli matrices in the space spanned by @xmath37 and @xmath38 .
we give the equation of motion for @xmath209 for the intermediate and high current regimes , in which we may approximate @xmath224 \approx { i}(\cos \varphi-\cos\varphi')/2 $ ] if @xmath225 and @xmath224 \approx 0 $ ] otherwise ( including the case @xmath226 ) .
@xmath227 ,    \nonumber \\
\frac{d\rho_x}{dt } & = & - \frac{j { i}}{8 } \rho_y \cos \theta     + \frac{j^2 { i}}{8 } \left\{-\rho_x [ 4 \sin^2 \theta + ( 2 s+1)^2 ]    + 2 ( \rho_z \cos \theta - 2 \rho_0 ) ( 2 s+1 ) \sin \theta \right\ }     - \frac{1}{2 } \rho_x \sum_{\pm } \gamma_{\pm 3/2,\pm 1/2 } , \nonumber \\    \frac{d \rho_y}{dt } & = &     \frac{j { i}}{16}\left [ 2 \rho_x \cos \theta + \rho_z ( 2 s+1 ) \sin
\theta    \right ]     - \frac{j^2 { i}}{8 } \rho_y    \left [ 4 \sin^2 \theta + ( 2 s+1)^2 \cos^2 \theta \right ]    - \frac{1}{2 } \rho_y \sum_{\pm } \gamma_{\pm 3/2,\pm 1/2 } , \nonumber \\    \frac{d \rho_z}{dt } & = &     - \frac{j { i}}{16 } ( 2 s+1 ) \rho_y
\theta    + \frac{j^2 { i}}{8 } ( 2 s+1 )    \left\ { ( 2 s+1)[2 \rho_0 \cos \theta - \rho_z ( \cos^2 \theta + 1 ) ]    + 2 \rho_x \sin \theta \cos \theta \right\ } \nonumber \\ & & \mbox { }    + \sum_{\pm } \left [ \pm \gamma_{\pm 1/2,\pm 3/2 } p_{\pm 3/2 }    - \frac{1}{2 } \gamma_{\pm 3/2,\pm 1/2 } ( \rho_z \pm \rho_0)\right ] ,     \label{eq : lindblad2d}\end{aligned}\ ] ]    in the three - dimensional case the equations of motion are @xmath227 ,    \nonumber \\
\frac{d\rho_x}{dt } & = & - \frac{j { i}}{8 } \rho_y \cos \theta    - \frac{j^2 { i}}{3 \pi^2 }     \left\ {     \rho_x \left [ 4 + \sin^2 \theta + 2 ( 2 s+1)^2 \right ]     + 4 \sin \theta \cos \theta \rho_z \right\ }    - \frac{j^2 { i}}{4 } \rho_0 ( 2 s+1 ) \sin
\theta     \nonumber \\ & & \mbox { }    - \frac{1}{2 } \rho_x \sum_{\pm } \gamma_{\pm 3/2,\pm 1/2 } , \nonumber \\
\frac{d \rho_y}{dt } & = &    \frac{j { i}}{16 }    \left [ 2 \rho_x \cos \theta + \rho_z ( 2 s+1 ) \sin \theta \right ]     - \frac{j^2 { i}}{3 \pi^2 } \rho_y    \left[4 + 4 \sin^2 \theta + ( 2 s+1)^2(1 + \cos^2 \theta ) \right ]    \nonumber \\ & & \mbox { }    - \frac{1}{2 } \rho_y \sum_{\pm } \gamma_{\pm 3/2,\pm 1/2 } , \nonumber \\
\frac{d \rho_z}{dt } & = &     - \frac{j { i}}{16 } ( 2 s+1 ) \rho_y \sin
\theta    + \frac{j^2 { i}}{3 \pi^2 } ( 2s+1 )    \left[2 \rho_x \sin \theta \cos \theta - \rho_z ( 2 s+1 ) ( 4 - \sin^2 \theta ) \right ]    + \frac{j^2 { i}}{8 } ( 2s+1)^2 \rho_0 \cos \theta    \nonumber \\ & & \mbox { }    + \sum_{\pm } \left [ \pm \gamma_{\pm 1/2,\pm 3/2 } p_{\pm 3/2 }    - \frac{1}{2 } \gamma_{\pm 3/2,\pm 1/2 } ( \rho_z \pm \rho_0)\right ] .
\label{eq : lindblad3d}\end{aligned}\ ] ]    in the intermediate - current regime we may set @xmath228 see eqs .
( [ eq : gammaanis ] ) and ( [ eq : gamma3dintermediate ] ) . for the high - current regime
we may approximate @xmath229 by eqs .
( [ eq : gammahigh2d ] ) and ( [ eq : gammahigh3d ] ) for the two - dimensional and three - dimensional cases , respectively .
when searching for a stationary solution , the contributions from transitions to / from the spin states @xmath230 vanish in the equations of motion for @xmath231 and @xmath232 because of detailed balance .
they remain , however , in the equations of motion for @xmath233 and @xmath234 .
once the ratio @xmath235 has been determined from the stationary solution of the lindblad equations ( [ eq : lindblad2d ] ) or ( [ eq : lindblad3d ] ) , the remaining probabilities @xmath65 for @xmath208 can be found from detailed balance , @xmath236 and from the normalization condition @xmath237 .
we now report expressions for the ratio @xmath235 for the intermediate and high current regimes @xmath238 and @xmath105 for a spin coupled to the edge of a two - dimensinal ti and to the surface of a three - dimensional ti .
the expressions for the intermediate - current regime are simplified using the inequalities @xmath239 and @xmath240 , but we make no assumptions regarding the relative magnitude of these two quantities .
( note that the intermediate current regime does not exist for @xmath3 . )
the expressions for the high - current limit are simplified using the inequality @xmath240 .    in the intermediate - current regime
we find for the case of a molecular magnet coupled to the edge of a two - dimensional ti @xmath241 with @xmath242 this result agrees with what one obtains from naive application of the rate equations if @xmath243 , but differs otherwise because of the effect of the current - induced exchange field ( [ eq : hexexpr ] ) .
note that @xmath244 .
the resulting steady - state probability distribution is @xmath245 the current - induced switching rate from @xmath70 to @xmath69 can be calculated as the transition rate from @xmath38 to @xmath37 starting with the initial condition @xmath246which is the reduced density matrix after current - induced equilibration of all states with @xmath247 , including the state @xmath38 , but without including transitions from @xmath38 to @xmath37 . within a time @xmath248 , which is short in comparison to the time required to transition from to @xmath38 to @xmath37 , a quasi - steady state
is reached for the reduced density matrix @xmath209 , in which the off - diagonal component @xmath234 becomes nonzero , @xmath249 the current - induced switching rate @xmath74 is then found as the rate of change of @xmath250 in this quasi - steady state , which gives @xmath251 where @xmath252 was defined in eq .
( [ eq : gplus2dintermed ] ) .
the switching rate for the opposite process is given the same expression , but with the replacement @xmath94 .
for the high - current limit and a two - dimensional ti we find @xmath253 with the function @xmath254 defined as @xmath255 which is different from the result one obtains from naive application of the rate equations .
for @xmath3 the ratio @xmath235 is consistent with a spin polarized in the @xmath5 direction .
the resulting probability distribution @xmath65 reads @xmath256 equation ( [ eq : gammaswitch2dhigh ] ) of the main text is a good order - of - magnitude estimate for the rate @xmath74 at which this steady - state distribution is approached .    for a molecule on the surface of a three - dimensional ti with an applied current in the intermediate - current regime
@xmath257 the stationary solution of the lindblad equation ( [ eq : lindblad3d ] ) gives the ratio @xmath253 with @xmath258 where the function @xmath259 was defined in eq .
( [ eq : fdef ] ) of the main text .
this expression simplifies to the ratio @xmath260 one obtains by considering rates only in the limit @xmath261 .
the corresonding steady - state distribution reads @xmath262 the switching rate from @xmath70 to @xmath69 is calculated in the same way as for the case of a two - dimensional ti , and one finds @xmath263 the switching rate for the inverse transition is again given by the replacement @xmath94 .
finally , for the high - current regime in three dimensions we find that the stationary solution of the lindblad equation ( [ eq : lindblad3d ] ) has @xmath253 with @xmath264 the corresponding steady - state distribution reads @xmath265 this distribution is approached at the rate given by eq .
( [ eq : gammaswitch2dhigh ] ) of the main text .
if the magnetic anisotropy axis is aligned with the current direction , the rate equation approach can be used throughout . in this case
the transition rates @xmath266 and @xmath267 for the intermediate and high - current regimes read @xmath268.\end{aligned}\ ] ] in the intermediate - current regime , the steady - state spin polarization of the molecule is complete ; the probability @xmath157 vanishes @xmath269 if @xmath2 is half integer , whereas all other probabilities are exponentially small in @xmath146 .
the switching rate out of the state @xmath70 is , in the intermediate - current regime , @xmath270 in the high - current regime the order - of - magnitude estimate ( [ eq : gammaswitch2dhigh ] ) for the switching rate also applies to the half - integer spin case .
30ifxundefined [ 1 ] ifx#1 ifnum [ 1 ] # 1firstoftwo secondoftwo ifx [ 1 ] # 1firstoftwo secondoftwo `` `` # 1''''@noop [ 0]secondoftwosanitize@url [ 0 ]
+ 12$12  & 12#1212_12%12@startlink[1]@endlink[0]@bib@innerbibempty    http://link.aps.org/doi/10.1103/revmodphys.76.323 [ * * ,   ( ) ] http://dx.doi.org/10.1146/annurev-conmatphys-070909-104123 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevapplied.4.047001 [ * * ,   ( ) ] http://www.sciencedirect.com/science/article/pii/0304885396000625 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevb.54.9353 [ * * ,   ( ) ] http://www.sciencedirect.com/science/article/pii/s0304885307010116 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.101.036601 [ * * ,   ( ) ] http://scitation.aip.org/content/aip/journal/apl/97/16/10.1063/1.3502596 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.106.036601 [ * * ,   ( ) ] http://www.nature.com/nature/journal/v476/n7359/abs/nature10309.html [ * * ,   ( ) ] http://science.sciencemag.org/content/336/6081/555 [ * * ,   ( ) ] http://www.nature.com/nmat/journal/v12/n3/abs/nmat3522.html [ * * ,   ( ) ] http://www.nature.com/articles/ncomms2709 [ * * ,   ( ) ] http://www.nature.com/nnano/journal/v8/n8/abs/nnano.2013.145.html [ * * ,   ( ) ] http://www.nature.com/nmat/journal/v12/n7/abs/nmat3675.html [ * * ,   ( ) ] http://www.nature.com/nnano/journal/v8/n7/abs/nnano.2013.102.html [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.98.106803 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/revmodphys.82.3045 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.108.187201 [ * * ,   ( ) ] http://dx.doi.org/10.7566/jpsj.83.064710 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.109.166602 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevb.93.115419 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevb.93.125303 [ * * ,   ( ) ] https://arxiv.org/abs/1509.06929v1 [ * * ] http://ieeexplore.ieee.org/document/7406697/ [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevb.90.205403 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevb.93.205130 [ * * ,   ( ) ] http://dx.doi.org/10.1038/nature13534 [ * * ,   ( ) ] http://www.nature.com/nmat/journal/v13/n7/abs/nmat3973.html [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.114.257202 [ * * ,   ( ) ] http://www.nature.com/nmat/journal/v7/n3/abs/nmat2133.html [ * * ,   ( ) ] , http://www.annualreviews.org/doi/abs/10.1146/annurev-conmatphys-070909-104053 [ * * ,   ( ) ] https://www.cambridge.org/core/journals/mrs-bulletin/article/single-molecule-magnets/65e2461c810e8de1a2ed2bf9620d3f96 [ * * ,   ( ) ] http://pubs.acs.org/doi/full/10.1021/nl2006142 [ * * ,   ( ) ] http://www.nature.com/nmat/journal/v8/n3/abs/nmat2374.html [ * * ,   ( ) ] http://dx.doi.org/10.1021/nl204141z [ * * ,   ( ) ] http://dx.doi.org/10.1021/nl503925h [ * * ,   ( ) ] http://dx.doi.org/10.1002/smll.200801594 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.102.256803 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.106.236402 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevb.92.045430 [ * * ,   ( ) ] http://www.nature.com/nphys/journal/v8/n6/full/nphys2322.html [ * * ,   ( ) ] http://www.sciencedirect.com/science/article/pii/092145269390146w [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevb.57.6602 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.90.166602 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevb.70.195345 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevb.71.195324 [ * * ,   ( ) ] http://link.aps.org/doi/10.1103/physrevlett.111.167203 [ * * ,   ( ) ] https://arxiv.org/abs/1501.02643 [ * * ] http://www.sciencedirect.com/science/article/pii/0034487772900109 [ * * ,   ] http://link.springer.com/article/10.1007/bf01608499
[ * * ,   ]