since last few decades one can observe a common pursuit towards miniaturization of electronic systems . according to the empirical moore s law
, the number of transistors per microchip is doubled every two years .
it is however obvious that the possibility of further miniaturization will be stopped in the near future due to the loss of chips stability when the device components achieve critical dimensions in the nanometer range .
thus , the most challenging task facing contemporary science and technology is to implement structures , alternative to the silicon - based devices , whose size could be reduced further .    the same tendency can also be observed in the development of hard discs , whose memory cells become smaller and smaller every year .
this generates new challenges related to information reading / writing .
an important step towards further miniaturization was the discovery of the giant magnetoresistance ( gmr ) effect .
the gmr was discovered in 1988 in artificially layered metallic structures consisting of ferromagnetic 3d films separated by nonmagnetic metallic layers @xcite .
it turned out that electrical resistance of magnetic metallic multilayers depends on their magnetic state , and usually drops when magnetic configuration varies from antiparallel alignment to the parallel one @xcite .
the gmr effect occurs for current flowing in the film plane as well as perpendicularly to it .
this phenomenon turned out to be very useful for applications in highly sensitive read heads , and allowed reading information from smaller memory cells using their weak magnetic field .
the discovery of gmr initiated broad interest in spin polarized electronic transport in nanoscopic systems .
it turned out that electron spin provides an additional degree of freedom , which considerably broadens the range of applications of mesoscopic systems in novel electronic devices . the spin - based nanoelectronics
called now spin electronics or shortly spintronics  is a relatively new area of mesoscopic physics dealing with the interplay of charge and spin degrees of freedom @xcite .
although a lot of theoretical and experimental works on the spintronic properties of mesoscopic systems have been carried out , this field is still in an early stage of development .
an effect similar to the aforementioned gmr also occurs when nonmagnetic metallic layer in a trilayer structure is replaced with a nonmagnetic insulating barrier , and the current flows owing to the phenomenon of quantum - mechanical tunnelling through the barrier .
this effect was discovered in 1975 in ferromagnetic planar junctions by jullire @xcite , and is of current interest due to applications in magnetic storage technology ( magnetic random access memories ) and in spin - electronics devices @xcite . as in the case of gmr ,
the tunnel magnetoresistance ( tmr ) consists in a decrease ( increase is also possible ) in the junction resistance when magnetic configuration of the junction changes from the antiparallel to parallel one .
tunneling in complex junctions , particularly in mesoscopic ones , where charging effects become important , is still not fully explored . a specific kind of such systems are double - barrier junctions with a small central electrode ( called an island ) .
such systems are known as single - electron transistors ( sets ) , mainly because electrons in a biased device flow one by one through the system and the transfer of single electrons can be controlled by a gate voltage .
electronic transport in such devices was extensively studied in the past decade , but mainly in the nonmagnetic limit @xcite .
recent experiments on magnetic nano - structured materials revealed new phenomena associated with the interplay of ferromagnetism and discrete charging effects .
a typical example is a ferromagnetic single - electron transistor , i.e. , a small metallic nanoparticle ( semiconducting quantum dots or molecules are also possible ) coupled by tunnel barriers to ferromagnetic electrodes @xcite .
first ferromagnetic single - electron transistors were fabricated by ono et al .
@xcite and later by brckl et al .
transport in ferromagnetic single - electron transistors with nonmagnetic metallic islands
both normal and superconducting  was also measured @xcite .
one should bear in mind , that the interplay of spin and charge effects was already studied long time ago in granular systems , in which magnetic nanoparticles were randomly dispersed in a nonmagnetic matrix @xcite .
recently granular films were investigated again by several groups @xcite , but both the size of the grains and also their location were strictly controlled .    in the case of sufficiently large metallic islands ( but still in the nanometer range ) , discrete structure of the electronic states in the grain is not resolved and is irrelevant . to observe the discrete electronic states in transport characteristics one should either diminish size of the metallic nanoparticles @xcite , or use semiconducting quantum dots based on two - dimensional electron gas @xcite .
an alternative strategy is to use ferromagnetic semiconducting materials instead of metallic ones as the electrodes @xcite .
magnetic impurities in the middle of the tunnel barrier of a ferromagnetic tunnel junction @xcite also can be considered as quantum dots with a very strong coulomb interaction .
another group of single - electron devices are molecular ferromagnetic transistors @xcite and , especially , ferromagnetically contacted carbon nanotubes @xcite .    in the following
we will review basic transport characteristics of ferromagnetic single - electron transistors in the sequential tunnelling , cotunneling , and strong coupling ( kondo ) regimes .
in particular , we will consider such properties of the device like charge and spin currents , tunnel magnetoresistance , spin accumulation , shot noise , kondo effect , and others . in section 2
we review basic principles of single - electron transport in the case of a large metallic island attached to ferromagnetic leads .
the size of the island is however small enough so that the charging energy is the dominant energy scale in the system .
the limits of fast and slow spin relaxation in the island are also discussed .
then we consider shot noise and the role of discrete electronic structure in such devices . in section 3
we consider devices based on double metallic islands . in turn , in section 4 we deal with electronic transport through single - level quantum dots ( qds ) in the sequential tunnelling and cotunneling regimes .
transport through multi - level qds , including transport through two - level quantum dots and carbon nanotubes attached to ferromagnetic leads is discussed in section 5 .
transport through quantum dots in the kondo regime is briefly addressed in section 6 .
final conclusions are in section 7 .      the systems considered in this review consist of a mesoscopic central part ( island ) coupled by tunnel barriers to external ferromagnetic leads .
the central part is characterized by an addition energy , which corresponds to the energy needed for adding a single electron , and includes contributions from the electrostatic charging energy and the discrete single - particle level separation . if the charging energy is the most relevant energy scale , the systems exhibit the single - electron charging effects @xcite .
an electron can tunnel to the central part only when the energy provided by the transport voltage surpasses the corresponding charging energy ; otherwise , the current is exponentially suppressed and the system is in the coulomb blockade regime @xcite .
once the bias voltage is larger than the threshold voltage , the electrons can tunnel one by one through the system leading to the step - like current - voltage characteristics
so - called coulomb staircase .
the blockade can also be overcome by applying a gate voltage @xmath0 that leads to sawtooth - like variation of electric current with @xmath0  so - called coulomb oscillations @xcite .
when the electrodes are made of a ferromagnetic material , the system exhibits further interesting phenomena resulting from the interplay of charge and spin degrees of freedom @xcite . in particular , the tunnelling current flowing through the system depends on the relative alignment of the magnetic moments of ferromagnetic electrodes , giving rise to the tmr effect @xcite , which is described quantitatively by the ratio @xmath1 where @xmath2 and @xmath3 denote the total system resistance in the parallel and antiparallel magnetic configurations , respectively , and @xmath4 , @xmath5 are the corresponding currents .
a simple theoretical model of tmr was introduced by jullire @xcite , who considered a single planar ferromagnetic tunnel junction and showed that tmr in such a device is given by @xmath6 , where the spin polarization @xmath7 of the lead @xmath8 ( @xmath9 for the left / right lead ) is defined as @xmath10 with @xmath11 being the spin - dependent density of states of lead @xmath8 for the spin - majority @xmath12 and spin - minority @xmath13 electrons .     a sketch illustrating ( a ) a single first - order ( sequential tunnelling ) process , and ( b ) a second - order ( cotunneling ) transport process . ]      in the following considerations we will distinguish between the three different transport regimes :    _ sequential tunnelling _  in the regime of weak coupling between the island and leads , and out of the coulomb blockade , electron transport is dominated by processes of the first order in the coupling parameter .
electrons flow then consecutively one by one due to the tunnelling events , and the transition rate from an initial state @xmath14 to a final state @xmath15 can be determined from the fermi golden - rule , @xmath16 where @xmath17 is the relevant tunnelling hamiltonian , and @xmath18 ( @xmath19 ) is the energy of the initial ( final ) state .
this transport regime is known as the sequential tunnelling regime @xcite .
an example of a single first - order process is sketched in fig .
[ sec1:fig1]a .
_ cotunneling _  although the sequential tunnelling in the coulomb blockade regime is exponentially suppressed , the current still flows due to higher - order tunnelling processes involving tunnelling of a single , two or more electrons _ via _ intermediate virtual states @xcite .
these processes are known as cotunneling .
an exemplary cotunneling process is illustrated in fig .
[ sec1:fig1]b .
the second order perturbation theory gives the cotunneling rate from an initial state @xmath14 to a final state @xmath15 @xcite @xmath20 where the summation is over all virtual states @xmath21 , and @xmath22 is the energy od the state @xmath23 ( @xmath24 ) .
in the coulomb blockade regime this tunnelling rate is only algebraically suppressed , contrary to the sequential transport processes which are then suppressed exponentially .
because of that even at low temperatures and in the strong coulomb blockade regime the rates of cotunneling processes do not vanish .
the second - order corrections become also important on resonance for intermediate coupling strengths .    _
strong coupling _  for strong coupling of the metallic island to electrodes , the tunnelling processes lead to logarithmic corrections to conductance , and the perturbation theory fails at the degeneracy points of two consecutive charge states . in the case of quantum dots additionally the kondo effect appears at low temperatures ( below the kondo temperature @xmath25 , @xmath26 ) leading to an enhanced conductance in the linear response regime @xcite .
in this section we review spin - polarized transport in a metallic ferromagnetic single - electron transistor ( fm set ) .
the device consists of a metallic nanoparticle as the central electrode ( island ) , which is coupled through tunnel barriers to external reservoirs of spin polarized electrons .
a gate voltage is attached capacitively to the island , which allows to control position of the corresponding fermi level .
electronic transport in nonmagnetic sets was already extensively studied in the past two decades @xcite .
recently , the attention was also drawn to electron tunnelling in magnetic systems @xcite , which was stimulated by recent progress in nanotechnology .
it has been shown theoretically that some qualitatively new effects may arise from the interplay of charging effects and spin degrees of freedom .
these include , for example , oscillations of tmr with increasing bias voltage , spin accumulation , enhancement of tmr in the coulomb blockade regime , etc . @xcite .
the enhancement of tmr in the cotunneling regime and the oscillations of tmr as a function of the transport voltage have also been observed experimentally @xcite .
the capacitance @xmath27 of few - nanometer - size particles is of the order of @xmath28 f @xcite .
consequently , the corresponding charging energy , @xmath29 , establishes a new relevant energy scale .
if the charging energy is larger than the thermal energy , @xmath30 , where @xmath31 denotes temperature and @xmath32 is the boltzmann constant , the effects due to discreteness of charge become observable in transport characteristics @xcite .
the single - electron transistors considered in this section are illustrated schematically in fig .
[ sec2:fig1 ] , where part ( a ) shows a device whose all three electrodes are ferromagnetic , whereas part ( b ) shows a system with ferromagnetic source and drain electrodes and a nonmagnetic island .
generally , magnetic moments of the leads may form an arbitrary magnetic configuration .
however , we focuss on systems whose moments form either parallel ( p ) or antiparallel ( ap ) magnetic configuration , as shown in fig .  [ sec2:fig1 ] . in the system shown in part ( a ) the magnetic moments of external electrodes point in the same direction , while magnetic moment of the island is either parallel or antiparallel to them . in the system shown in part ( b ) ,
magnetic moments of external electrodes in the parallel configuration are aligned , while in the antiparallel configuration they are anti - aligned .
the two magnetic alignments can be easily achieved by sweeping magnetic field through the hysteresis loop , provided the respective ferromagnetic components have different coercive fields .
one may also make use of exchange anisotropy to fix magnetic moment of a particular layer and rotate magnetic moment of the second layer with a weak magnetic field .
generally , there is a chance that electron tunnelling through a barrier will change its spin orientation .
however , we consider only spin - conserving tunnelling processes through the two barriers .
first , we assume that the islands are relatively large , so the effects due to quantization of the corresponding energy levels can be neglected . for smaller islands , however , the discrete energy spectrum may modify transport characteristics @xcite , and this will be considered later ( see section [ discreteset ] ) .
schematics of ferromagnetic single - electron transistors .
the parallel and antiparallel magnetic configurations of the system are also specified .
the island is separated from external electrodes by tunnel barriers .
each barrier is characterized by its spin - dependent resistance @xmath33 and capacitance @xmath34 ( @xmath35 ) .
the system is symmetrically biased , @xmath36 , @xmath37 , and there is also a gate voltage @xmath38 applied to the island . ]    to describe charge and spin transport we need to specify a model hamiltonian of the system .
first , the electrostatic energy required to add @xmath39 excess electrons to the island , while keeping constant voltages @xmath40 and @xmath41 in the left and right electrodes and the gate voltage @xmath42 , is given by @xcite @xmath43 where a constant term ( independent of @xmath39 ) is irrelevant and has been dropped .
here , the total island s capacitance @xmath27 is the sum of the capacitances of the left and right junctions and of the gate , @xmath44 .
the external charge @xmath45 accounts for the effect of applied voltages , and can be continuously tuned .
the total hamiltonian of the device takes then the general form @xmath46 here @xmath47 describes the decoupled leads and island , while @xmath48 takes into account the lead - island coupling .
the ferromagnetic leads and the island are described by @xmath49 for @xmath50 , where @xmath51 are the fermi operators for electrons with a wavevector @xmath52 and spin @xmath53 in the electrodes and island ( @xmath50 ) , and @xmath54 is the corresponding single - particle energy .
the last part of the hamiltonian , @xmath55 describes tunnelling processes between the leads and island , with @xmath56 being the relevant matrix elements .
when the resistances of both tunnel barriers are much larger than the quantum resistance , @xmath57 , and the system is not in the coulomb blockade regime , transport is dominated by sequential tunnelling processes . as a consequence , the charge is well localized in the island and the _ orthodox _ tunnelling theory is applicable @xcite .    in order to calculate electric current in the sequential transport regime
, one may use the method based on the master equation which is a detailed balance of electrons tunnelling to and off the island . in the stationary case
the master equation reads    @xmath58 p(n , v)\nonumber \\    & & + \sum\limits_{\sigma }    \left[\gamma _ { l\sigma}^+(n-1 )    + \gamma _ { r\sigma}^+(n-1)\right]p(n-1,v )    \nonumber\\    & & + \sum\limits_{\sigma }    \left[\gamma _ { l\sigma}^-(n+1 )    + \gamma _ { r\sigma}^-(n+1)\right]p(n+1,v),\end{aligned}\ ] ]    where @xmath59 is the probability to find the island in a state with @xmath39 additional electrons when a bias voltage @xmath60 ( @xmath61 ) is applied , and @xmath62 is the spin - dependent rate for tunnelling of electrons with spin @xmath53 from the lead @xmath8 to island ( upper sign ) and from the island to lead @xmath8 ( lower sign ) , when the island is occupied by @xmath39 excess electrons .
these tunnelling rates depend on the bias voltage ( not indicated explicitly ) , and can be expressed by means of the fermi golden rule as @xmath63 -1}\,,\ ] ] where @xmath64 describes a change in the electrostatic energy of the system caused by the corresponding tunnelling event , when in the initial state there were @xmath39 additional electrons on the island . in the above equation @xmath33 denotes the spin - dependent tunnel resistance of the @xmath8-th barrier , given by @xmath65 where @xmath66 has been assumed for simplicity . the spin dependence of the resistance is a consequence of the spin - dependent density of electron states at the fermi level in the respective electrodes and the corresponding tunnelling matrix elements . in particular , in eq .
( [ eq : setbarrierr ] ) it is due to the spin - dependent density of states of the @xmath8-th lead , @xmath67 , spin - dependent density of states of the island , @xmath68 , as well as due to the spin - dependent tunnelling matrix elements @xmath69 .    in a stationary state ,
the net transition rate between the charge states with @xmath39 and @xmath70 excess electrons on the island vanishes .
the probability @xmath59 can be then determined from the following recursion relation @xcite @xmath71 where @xmath72 and @xmath73 , corresponding to transition rates for tunnelling to and off the island , respectively .
generally , energy of an electron after tunnelling event is relaxed to the relevant fermi level in a short time scale .
one can assume that the energy relaxation time is the shortest time scale , shorter than the time between two successive tunnelling events .
however , such a restriction can not be imposed on the spin relaxation time which can be relatively long . in a general case
, a nonequilibrium magnetic moment may accumulate on the island due to the spin dependence of tunnelling processes , which leads to spin splitting of the corresponding fermi level . for arbitrary spin relaxation times
, the splitting of the fermi level can be determined from the spin balance @xcite @xmath74 where @xmath75 is the island s volume , @xmath76 is the electron charge ( @xmath77 ) , @xmath78 denotes the spin relaxation time on the island , while @xmath79 and @xmath80 are the currents flowing through the left and right junctions in the spin channel @xmath53 .
the last term in eq .
( [ eq : setspinconser ] ) takes into account intrinsic spin - flip processes on the island .
> from this condition it is possible to calculate self - consistently the shifts of the fermi levels due to spin accumulation for both spin orientations .
the electric current flowing through the left junction can be calculated from the following formula : @xmath81 p(n , v )    \,.\ ] ] similar formula also holds for @xmath82 . in the stationary state the currents flowing through both junctions are equal , @xmath83 .    in the following
we discuss two limiting cases : the limit of fast spin relaxation and the limit of slow spin relaxation on the island . in the former case the spin of an electron tunnelling to
the island relaxes before a next tunnelling event takes place . in the latter case , on the other hand
, the electron spin is conserved for a time much longer than the time between successive tunnelling events .      in the limit of fast spin relaxation ,
sets with nonmagnetic islands behave like nonmagnetic junctions , exhibiting no tmr effect . therefore
, the following discussion in this subsection is limited to fm sets with ferromagnetic islands , see fig .
[ sec2:fig1]a .
generally , one may assume that the spin - dependent resistances for the parallel ( @xmath84 ) and antiparallel ( @xmath85 ) configurations fulfill the condition @xmath86 , where @xmath87 .
this formula follows from eq .
( [ eq : setbarrierr ] ) when assuming that all spin effects are included into the spin dependent density of states ( transfer matrix element are independent of spin and magnetic configuration ) . moreover ,
since there is no spin accumulation on the island in the limit of fast spin relaxation , there is no associated spin splitting of the fermi level , @xmath88 .
current in the parallel and antiparallel magnetic configurations ( a ) and the resulting tmr ( b ) as a function of the bias voltage for a fm set with a ferromagnetic island in the limit of fast spin relaxation .
the parameters are : @xmath89 k , @xmath90 m@xmath91 , @xmath92 m@xmath91 , @xmath93 m@xmath91 , @xmath94 m@xmath91 , @xmath95 ( for @xmath35 ) , @xmath96 af , @xmath97 af , @xmath98 af , @xmath99 , and the offset charge @xmath100 . the dots in part ( b ) present the experimental data taken from ref .  @xcite . ]    the current flowing through the system in the parallel and antiparallel configurations is displayed in fig .
[ sec2:fig2]a . for both magnetic configurations
, the @xmath101 curves reveal the well - known coulomb steps . moreover ,
these two curves are different ; the current flowing in the parallel configuration is generally larger than the current flowing in the antiparallel configuration , see fig .
[ sec2:fig2]a .
this difference leads in turn to nonzero tmr effect , as shown in fig .
[ sec2:fig2]b .
the tmr effect has a component that oscillates as a function of the bias voltage .
the amplitude of these oscillations , however , decreases as the transport voltage increases . for the parameters assumed here , tmr reaches local maxima at the voltages corresponding to the positions of coulomb steps .
however , this is not a general rule , and for other parameters tmr can have local minima at the coulomb steps . the global maximum value of tmr in fig .
[ sec2:fig2]b appears at the first step , i.e. , at the threshold voltage .
when the temperature increases , the effects due to discrete charging , i.e. the coulomb steps and enhancement of tmr at the coulomb steps become diminished and disappear at @xmath102 @xcite . the oscillatory behavior of the tmr effect with increasing transport voltage was observed experimentally , for example by ernult _
et al . _
@xcite , see the dots in fig .
[ sec2:fig2]b .
the curves presented in fig .  [ sec2:fig2 ] are calculated for the parameters corresponding to those in ref .
@xcite . to get a good agreement with experimental observations ,
a nonzero offset charge @xmath103 on the island ( due to external charges ) has been assumed . as can be seen in fig .
[ sec2:fig2]b , there is a satisfactory agreement between the theoretical curve and experimental data for voltages above the threshold , @xmath104mv , while in the coulomb blockade regime transport calculated using the sequential tunnelling approximation is not properly described .      in the case of a fm set with a nonmagnetic island ,
see fig .
[ sec2:fig1]b , a nonzero tmr in the sequential tunnelling regime can exist only when the spin relaxation time is sufficiently long , i.e. significantly longer than the time between successive tunnelling events .
we note that the longest spin relaxation times were measured for aluminium and copper @xcite ( for example , the relaxation time for copper was estimated to be of the order of @xmath105 s ) .
if this is the case , a nonequilibrium magnetic moment builds up on the island due to spin accumulation .
this moment leads a nonvanishing tmr .
in other words , the island becomes magnetized in a nonequilibrium situation , and the created moment depends on the bias and gate voltages .
basic characteristics of a fm set with a nonmagnetic island as a function of the bias voltage in the limit of slow spin relaxation on the island : ( a ) spin splitting of the fermi level , ( b ) @xmath101 characteristics , and ( c ) tunnel magnetoresistance .
the parameters are : @xmath106 , @xmath107 m@xmath108 , @xmath109 m@xmath91 , @xmath110 m@xmath108 , @xmath111 m@xmath91 , @xmath112 af , and @xmath99 . ]    since the density of states at the fermi level in a nonmagnetic island is independent of the spin orientation , one finds @xmath113 .
apart from this , the resistances in the antiparallel configuration are @xmath114 and @xmath115 , see also fig .
[ sec2:fig1]b .
typical transport characteristics of a fm set with nonmagnetic island in the limit of slow spin relaxation are shown in fig .
[ sec2:fig3 ] as a function of the bias voltage .
the splitting of the fermi level is displayed in fig .
[ sec2:fig3]a for both magnetic configurations .
as one could expect , the splitting takes place only in the antiparallel alignment , while in the parallel configuration there is no spin accumulation . in the former case
the ratio of tunnelling rates for electrons with opposite spin orientations becomes the same for electrons tunnelling to and off the island only when a nonequilibrium magnetic moment is built on the island . in the parallel configuration ,
however , this condition is already fulfilled without any spin accumulation . however , this is true only in the case when both junctions are characterized by equal spin asymmetries .
when the spin asymmetries of both electrodes are different , spin accumulation also occurs in the parallel configuration .    as illustrated in fig .
[ sec2:fig3]a , behavior of the fermi level splitting with increasing transport voltage can be decomposed into two components .
one component monotonously increases , while the second one oscillates with increasing bias voltage .
this oscillatory behavior can be accounted for in the following way .
let us assume that the voltage is slightly above that corresponding to a certain coulomb step and begins to increase .
then , the spin accumulation also increases , until a local maximum value is reached .
the local maximum occurs at a voltage , at which the chemical potential of the depleted spin channel approaches the value which allows the next charge state in the island .
this , in turn , enhances tunnelling rate ( onto the island ) of electrons corresponding to the depleted spin channel , and consequently reduces the spin splitting of the fermi level .
when the voltage increases further , a local minimum in the spin accumulation is then reached at a voltage , where the chemical potential of the second ( accumulated ) spin channel approaches the value which allows the next charge state on the island .
the same scenario repeats at each coulomb step leading to the oscillatory component in the spin accumulation . in the sequential tunnelling regime , where only the first - order tunnelling processes are taken into account , spin accumulation is exponentially small in the coulomb blockade regime , as can be seen in fig .
[ sec2:fig3]a .
the current as a function of the bias voltage for the parallel and antiparallel configurations is shown in fig .
[ sec2:fig3]b . as before ,
characteristic coulomb steps are clearly visible . moreover , owing to different stationary spin accumulations in the parallel and antiparallel configurations , the current in the parallel configuration is also different from that in the antiparallel configuration .
this gives rise to the tmr effect which is presented in fig .
[ sec2:fig3]c .
now , the bias dependence is more complex than it was in the case of fm sets with magnetic islands in the absence of spin accumulation .
it is interesting to note , that tmr can change sign in the vicinity of the coulomb steps in current - voltage curves , which is a consequence of different spin accumulations in the two magnetic configurations .    generally , it is more difficult to obey the slow spin relaxation limit in ferromagnetic islands than in nonmagnetic ones .
anyway , some experimental data show that this is achievable . in such a case
, the spin accumulation builds up in the magnetic island and has significant influence on transport characteristics . in refs .
@xcite it has been shown that spin accumulation modifies the staircase-like variation of the electric current with the bias voltage .
since the spin accumulation depends on the magnetic configuration of the junction , this can also lead to an enhanced tmr effect .
moreover , in some voltage regions tmr can change sign . in addition , the coulomb steps for the two magnetic configurations become slightly shifted , as was also observed experimentally @xcite .
the difference between current - voltage curves in the fast and slow spin relaxation limits varies continuously with spin relaxation rates .
this difference was used very recently to evaluate spin relaxation time in co nanoparticles @xcite .      the electric current and spin accumulation in the limit of sequential tunnelling ( see the preceding subsection ) are exponentially small in the coulomb blockade regime . however , charge transport may occur not only due to spin - dependent sequential tunnelling , but also due to spin - dependent cotunneling processes @xcite , which give dominant contribution in the coulomb blockade regime and also lead to spin accumulation @xcite .
close to resonance ( in the vicinity of the threshold voltage ) both sequential and cotunneling currents may be comparable .
the current @xmath116 is then equal to the sum of first , @xmath117 , and second , @xmath118 , order contributions , @xmath119 .
to calculate cotunneling current far from the resonance one could use eq .
( [ sec1:eq : cotrate ] ) .
this formula , however , can not be used when voltage approaches the threshold voltage ( at resonance ) . to calculate current and associated spin accumulation in the whole voltage range , the real - time diagrammatic formalism @xcite has been used .
the spin accumulation on the island ( or equivalently spin splitting of the electrochemical potential ) is then determined from the spin balance equation , similarly as in the sequential tunnelling regime , see eq .
( [ eq : setspinconser ] ) .
differential conductance in fm sets with nonmagnetic islands as a function of the bias voltage @xmath60 in the ( a ) parallel and ( b ) antiparallel configurations , calculated in the whole ( sequential and cotunneling ) transport regime in the limit of slow spin relaxation in the island .
the curves have been calculated for @xmath120 and symmetric junctions with @xmath121 .
the spin polarization is @xmath122 and @xmath123 for ni and fe electrodes , respectively .
the dotted line corresponds to the sequential tunnelling limit in the set with fe electrodes .
( after ref .
@xcite ) ]    the differential conductance in the whole ( sequential and cotunneling ) transport regime is shown in fig .  [ sec2:fig4 ] for sets with ni and fe electrodes and for nonmagnetic islands ( in the limit of slow spin relaxation ) .
upper part corresponds to the parallel configuration while the lower one to the antiparallel one .
first , we find well resolved splitting of the conductance peaks in the antiparallel alignment , while no splitting can be seen in the parallel configuration .
this splitting in conductance peaks is a direct consequence of the spin splitting of the corresponding electrochemical potential of the island , and therefore can be used to detect and measure spin accumulation .
the absence of conductance splitting in the parallel configuration is simply a consequence of the absence of spin accumulation in this configuration .
such a splitting of the conductance peaks in the cotunneling regime was recently observed experimentally @xcite .
generally , there are several experimental techniques by which the spin accumulation can be detected indirectly @xcite .
the peculiarities of transport characteristics of fm sets offer new possibilities .
( color online ) the influence of discrete electronic structure on transport characteristics in a fm set with nonmagnetic island in the limit of slow spin relaxation .
voltage dependence of the tunnel current @xmath116 ( a ) , the differential conductance @xmath124 ( b ) , spin accumulation @xmath125 ( c ) , standard deviation @xmath126^{1/2 } $ ] ( d ) , and tunnel magnetoresistance @xmath127 ( e ) , calculated at @xmath128 . the solid and dashed curves in ( a ) , ( c ) and ( d ) correspond to the antiparallel and parallel configurations , respectively .
the other parameters are : @xmath129 , @xmath130 , @xmath131 , @xmath132 , @xmath133 and @xmath134 for the parallel alignment ( @xmath135 for the antiparallel alignment ) .
( after ref .
@xcite ) ]    discussion up to now was limited to the case when addition of a single electron costed only an electrostatic energy .
when , however , the size of metallic island becomes reduced further , electron spectrum of the island can not be considered as a continuous one , and discrete structure of the energy levels plays an important role .
addition of an electron costs then not only the electrostatic charging energy , but also the energy equal to the level spacing @xmath136 . in this subsection
we describe briefly the effect of discreteness on the transport characteristics .
the relevant approach was developed in refs . @xcite .
the electron transport in the stationary state is then governed by the solution of the generalized master equation @xcite @xmath137 where @xmath138 denotes the probability to find @xmath139 and @xmath140 excess electrons on the island ( @xmath141 is the total number of excess electrons ) . the first term in eq .
( [ seq2eq:1 ] ) describes how the probability of a given charge and spin state decays due to electron tunnelling to and from the island , whereas other terms describe the rate at which this probability increases .
the @xmath91-terms account for spin - flip relaxation processes .
the coefficients entering eq .
( [ seq2eq:1 ] ) are defined as @xmath142 and @xmath143 $ ] , where @xmath144 are the tunnelling rates for electrons with spin @xmath53 , tunnelling to ( @xmath145 ) the island from the lead @xmath146 or backward ( @xmath147 ) .
these coefficients are given by @xcite @xmath148 here , @xmath149 is the fermi function ( @xmath150 ) , whereas @xmath151 ( @xmath152 ) describes the probability that the energy level @xmath153 is occupied by an electron with spin @xmath53 for a particular configuration @xmath154 .
the parameter @xmath155 is the bare tunnelling rate of electrons between the lead @xmath156 and the energy level @xmath153 of the island , and @xmath157 is the transition probability from the state @xmath158 to @xmath159 of the island due to the spin - flip processes .
the energies @xmath160 and @xmath161 are given by @xmath162 and @xmath163 where @xmath164 $ ] and @xmath165 .
> from the solution @xmath138 of the master equation ( eq .  [ seq2eq:1 ] ) , one can determine current flowing through the island , @xmath166 for further discussion we assume that the discrete energy levels @xmath153 are spin degenerate ( nonmagnetic situation ) and equally separated with the level spacing @xmath136 .    to emphasize the role of spin accumulation
let us assume that the intrinsic spin relaxation time on the island is long enough to neglect all intrinsic spin - flip processes .
the corresponding @xmath101 characteristics for the parallel and antiparallel alignments are shown in fig .
[ sec2:fig6]a . in both cases
the electric current is blocked below a threshold voltage , and a typical coulomb staircase appears above it , with additional small steps due to the discrete levels of the island .
the effects due to discrete charging and discrete electronic structure are more clearly seen in the differential conductance shown in fig .
[ sec2:fig6]b , where the small peaks correspond to new discrete levels taking part in transport .
the difference between the @xmath116-@xmath60 characteristics for the parallel and antiparallel configurations is due to a different spin accumulation in both geometries . in fig .
[ sec2:fig6]c we see the average value of the difference between the numbers of spin - up and spin - down excess electrons on the island , @xmath167 , i.e. , the spin accumulation .
there is no significant spin accumulation in the parallel configuration .
the number @xmath168 of spins accumulated on the island fluctuates in time around its average value @xmath169 , as shown in fig .
[ sec2:fig6]d , where the standard deviation @xmath170^{1/2}$ ] is plotted against the voltage @xmath60 .
it is worth to note that although there is almost no spin accumulation in the parallel configuration , the corresponding fluctuations are relatively large .
the difference between the @xmath101 curves in the parallel and antiparallel configurations leads to the tunnel magnetoresistance , shown in fig .
[ sec2:fig6]e .
the broad peaks correspond to the coulomb steps , while the fine structure originates from the discrete structure of the density of states of the island .
one can extend the presented formalism also to the case of a fm set with a ferromagnetic island @xcite and calculate nonequilibrium spin accumulation on the island ( when spin relaxation time is sufficiently long ) . using similar
approach inoue and brataas @xcite analyzed current - induced magnetization reversal induced by spin accumulation rather than by spin torque .
they found that the magnetization reversal is possible when a free energy change due to nonequilibrium spin accumulation is comparable to the anisotropy energy .
> from the application point of view , an important transport characteristics of the system is the corresponding current noise .
functionality of a system depends on the relevant noise to signal ratio , which should be as small as possible . however , the noise is also a source of additional information on the system properties like quantum and coulomb correlations  @xcite .
the shot noise in ferromagnetic single - electron transistors was studied theoretically in ref .
@xcite , where the method developed for spinless electrons in single - electron transistors  @xcite was extended to magnetic systems .
the time correlation function of any two quantities @xmath171 and @xmath172 can be expressed as  @xcite @xmath173 where , @xmath174 is the conditional probability to find the system in the final state with @xmath175 and @xmath176 excess electrons at time @xmath177 , if there was @xmath139 and @xmath140 excess electrons in the initial time t=0 .
the probability @xmath178 can be determined from eq .
( [ seq2eq:1 ] ) .
voltage dependence of the current shot noise at @xmath179 ( a ) in a fm set with nonmagnetic island and in the slow spin relaxation limit for @xmath180 , @xmath181 for the antiparallel alignment ( @xmath182 , for the parallel configuration ) ( other parameters as in fig .
[ sec2:fig6 ] ) . in part
( b ) @xmath183 is split into two components : @xmath184 ( upper curves ) and @xmath185 ( lower curves ) .
( after ref .
@xcite ) ]    following this procedure , the current - current correlation function in a fm set considered in the preceding subsection has been calculated in ref .
the corresponding fourier transform can be presented as @xmath186 where @xmath187 is the schottky value ( the frequency independent part ) , while the second term in eq .
( [ sec2eq:5 ] ) is the frequency dependent component .
figure  [ sec2:fig7]a shows the bias dependence of the zero - frequency current noise @xmath183 .
the current noise is smaller in the antiparallel configuration than in the parallel one .
this is because in the presence of spin accumulation ( which is significant only in the antiparallel alignment ) the amplitude of fluctuations is smaller . in fig .
[ sec2:fig7]b @xmath183 is split into two parts ; the frequency independent component @xmath184 and the contribution @xmath185 arising from the frequency dependent part of the current noise , see eq .
( [ sec2eq:5 ] ) .
the component @xmath184 is almost constant @xmath188 at the plateaux of the @xmath101 curve and increases with opening of new channels .
dynamical correlations between the currents are described by @xmath189 .
its value in the limit @xmath190 can be positive between the @xmath101 steps and negative when new channels become open .
this is evident for the antiparallel alignment at @xmath191 , when opening a tunnelling channel for electrons with @xmath192 leads to negative dynamical correlations .
this effect is almost compensated by an increase in @xmath184 , and therefore one gets only a small reduction of the current noise @xmath183 .    in the power spectrum of the current
@xmath189 one can distinguish two distinct relaxation times , one in the high and another one in the low frequency regions @xcite . in a wide voltage range
the corresponding relaxation times are very close to the effective relaxation times for the charge and spin noise .
the asymmetry between the tunnelling channels for electrons with the opposite spins leads to activation of the spin component in the current noise . in ref .
@xcite the components @xmath193 and @xmath194 , corresponding to the charge and spin noise , respectively , have been extracted from @xmath195 .
it has been shown that the charge component is almost constant whereas the spin component increases with spin polarization @xmath196 of the leads , and for @xmath197 can be much larger than the charge component .
the analysis showed that both charge and spin fluctuations are relevant for the shot noise in fm sets .
it has been also pointed out that a super - poissonian shot noise , the fano factor @xmath198 , can occur due to lifting of the spin degeneracy .
some new features of transport characteristics appear in spin - polarized electronic transport through double - island structures @xcite . a typical double - island device , shown schematically in fig .  [ sec3:fig1 ] , consists of two metallic islands separated from each other and from external reservoirs by tunnel barriers .
a gate voltage is additionally attached capacitively to each island .
generally , both islands and external electrodes can be either magnetic or nonmagnetic @xcite . in the following , we consider the situations when at least two electrodes ( external or central ) are ferromagnetic , and their magnetic moments form only collinear ( parallel and antiparallel ) configurations . apart from this , we limit discussion to continuous density of states in the islands .
schematic of a double - island device .
the first ( 1 ) and second ( 2 ) islands are capacitively coupled to two gate voltages , @xmath199 and @xmath200 , and separated from each other and from the left and right electrodes by tunnel barriers . ]
hamiltonian of the double - island systems is similar to that used for single - island devices , but now the electrostatic energy is given by @xcite @xmath201 where @xmath202 is the number of excess electrons on the first ( second ) island and @xmath203 describes the respective charge induced by applied voltages , @xmath204 .
apart from this , @xmath205 and @xmath206 denote the charging energies of the corresponding islands , whereas @xmath207 describes electrostatic coupling of the islands , @xmath208 with @xmath209 being the total capacitance of the first ( second ) island , @xmath210 , and @xmath211 denoting the capacitance of the middle junction .    to find the stationary current flowing through the system one needs then to know the probabilities @xmath212 of finding the system in the charge state with @xmath213 and @xmath214 additional electrons on the first and second islands , respectively , when a bias voltage @xmath60 is applied .
these probabilities can be calculated in a recursive way from the following steady - state master equation @xmath215 p(n_1,n_2,v)\nonumber \\    & & + \sum\limits_{\sigma } \gamma _ { l    1}^\sigma(n_{1}-1,n_{2})p(n_{1}-1,n_{2},v)\nonumber\\    & & + \sum\limits_{\sigma } \gamma_{1 l }    ^\sigma(n_{1}+1,n_{2})p(n_{1}+1,n_{2},v ) \nonumber\\    & & + \sum\limits_{\sigma}\gamma _ { 1    2}^\sigma(n_{1}+1,n_{2}-1)p(n_{1}+1,n_{2}-1,v)\nonumber\\    & & + \sum\limits_{\sigma}\gamma    _ { 2 1}^\sigma(n_{1}-1,n_{2}+1)p(n_{1}-1,n_{2}+1,v ) \nonumber\\    & & + \sum\limits_{\sigma}\gamma _ { r    2}^\sigma(n_{1},n_{2}-1)p(n_{1},n_{2}-1,v)\nonumber\\    & & + \sum\limits_{\sigma}\gamma _ { 2    r}^\sigma(n_{1},n_{2}+1)p(n_{1},n_{2}+1,v ) \,,\end{aligned}\ ] ] with the normalization condition @xmath216 .
the corresponding transition rates are given by eq .
( [ eq : setrate ] ) , with the respective changes in the total system electrostatic energy .
for example for tunnelling from the first island to the second one , the change in the electrostatic energy can be written as @xmath217 here , @xmath218 and @xmath219 denote the corresponding shifts of the chemical potentials for spin @xmath53 in the first and second islands , respectively .
the spin asymmetry of tunnelling through the barrier @xmath8 will be characterized in this section by the asymmetry factor @xmath220 , defined as @xmath221 .
this asymmetry factor corresponds to the ratio of the respective spin - dependent densities of states . in particular , for fm / nm junctions it is given by @xmath222 , whereas for fm / fm junctions @xmath223 , provided that the two ferromagnetic electrodes are built of the same material .
we note that the relation between spin dependent barrier resistances in the parallel and antiparallel configurations of magnetic moments on the opposite sides of the barrier is the same as that described in sec.2 . when a bias voltage is applied to the system
, a nonequilibrium spin accumulation may appear on the islands . generally , the shifts of the fermi level for both spin orientations are different .
however , one can assume that the ratio of the fermi level shifts for the spin - up and spin - down electrons , defined as @xmath224 , fulfills the relation @xmath225 , for the first ( @xmath226 ) and second ( @xmath227 ) island , respectively , with @xmath228 being the spin - dependent density of states of the @xmath229-th island .
as a consequence , for nonmagnetic islands one directly gets , @xmath230 .
electric current flowing through the system can be determined from the following formula : @xmath231    p(n_{1},n_{2},v)\,,\end{aligned}\ ] ] which corresponds to the current flowing through the left junction , but @xmath232 in the stationary state .       the average electron number on the islands ( a ) , currents ( b ) in the parallel and antiparallel configurations and the resulting tmr ( c ) as a function of the bias voltage .
the parameters are : @xmath233 af , @xmath234 .
the spin asymmetries of resistances in the parallel configuration are @xmath235 , @xmath236 , whereas the total junction resistances are @xmath237 m@xmath91 . in the antiparallel configuration , @xmath238 , @xmath239 , and @xmath240 . ]    consider first transport characteristics of a system built of two ferromagnetic islands and nonmagnetic external electrodes , shown in fig .
[ sec3:fig2 ] .
different magnetic configurations of the system are specified in the inset of fig .
[ sec3:fig2]c . when the spin relaxation time in the islands is much shorter than the time between two successive tunnelling events , no spin accumulation builds up on the islands . since there is an asymmetry between the two barrier resistances , some charge accumulates on the islands , as displayed in fig .
[ sec3:fig2]a .
first of all , the magnitude of excess charge on the islands increases with increasing voltage . for the barrier asymmetry assumed in fig .
[ sec3:fig2 ] , the electrons easier tunnel to the second island from the right lead than out of the second island to the first one . as a consequence ,
the electrons accumulate on the second island . on the other hand , the electrons easier tunnel out of the first island to the left lead than from the second island to the first one .
thus , the number of excess electrons on the first island decreases with increasing the bias voltage  there are holes accumulated on the first island .
this occurs in both magnetic configurations .
the electric current flowing through the system in the parallel and antiparallel configurations is shown in fig .
[ sec3:fig2]b .
the coulomb steps due to discrete charging are clearly evident , and the difference in currents flowing in the parallel and antiparallel configurations gives rise to the tunnel magnetoresistance , shown in fig .
[ sec3:fig2]c .
the tmr effect oscillates as a function of the bias voltage and the amplitude of these oscillations decreases as the voltage is increased in a similar way as in the case of single - island fm sets , see fig .
. moreover , some dips occur now in tmr at the voltages corresponding to the steps in the current - voltage characteristics .     the shifts of the fermi levels for spin - up electrons ( a ) , currents in the parallel and antiparallel configurations ( b ) and tmr ( c ) as a function of the bias voltage .
the parameters are : @xmath233 af , @xmath241 , @xmath230 , @xmath242 , @xmath243 , @xmath235 , @xmath244 , whereas @xmath245 m@xmath91 . in the antiparallel configuration , @xmath246 , @xmath247 . ]
when the spin relaxation time is longer than the time between two successive tunnelling events , a nonequilibrium magnetic moment appears on each island .
the corresponding shifts of the fermi level due to spin accumulation can be calculated from spin current conservation , see eq .
( [ eq : setspinconser ] ) , written for each island @xcite .    figure  [ sec3:fig3 ] shows the shifts of the fermi levels and currents in the parallel and antiparallel configurations , as well as the resulting tmr calculated in the limit of long spin relaxation time for the system built of ferromagnetic electrodes and nonmagnetic islands .
the magnetic moments of external electrodes can form either parallel or antiparallel configurations , as illustrated in the inset of fig .
[ sec3:fig3]c .
first of all , the nonequilibrium spin accumulation , shown in fig .
[ sec3:fig3]a , exists not only in the antiparallel configuration but also in the parallel one . in the antiparallel configuration
the shifts of the fermi level for a given spin orientation are equal on both islands , whereas in the parallel configuration they are opposite .
the effects due to discrete charging lead to an oscillatory behavior of the fermi level shift , in a similar way as in the case of single - island fm sets with nonmagnetic islands discussed above .
the currents flowing through the system in both magnetic configurations are shown in fig .
[ sec3:fig3]b . due to
nonequilibrium spin accumulation induced on the islands , these currents are different , which leads to nonzero tunnel magnetoresistance , as displayed in fig .
[ sec3:fig3]c .
it is interesting to note , that now tmr changes sign in certain transport voltage regions .
these effects are clearly due to magnetic moments accumulated on the islands . if the spin relaxation time becomes shorter than the time between successive tunnelling events , spin accumulation disappears and , consequently , tmr also vanishes .
the shifts of the fermi levels for majority - spin electrons ( a ) , currents ( b ) and differential conductance ( c ) in the parallel and antiparallel configurations , and the resulting tmr ( d ) as a function of the bias voltage .
the parameters are : @xmath248 k , @xmath249 af , @xmath250 af , @xmath251 af , @xmath252 , @xmath253 , and @xmath254 , whereas @xmath255 .
the total junction resistances are @xmath256 m@xmath91 , @xmath257 m@xmath91 . in the antiparallel configuration , @xmath258 , and @xmath259 , for @xmath260 .
( after ref .
@xcite ) ]    recently , several experiments on spin - polarized transport through granular systems were reported @xcite .
for example , in ref .
@xcite , the tunnelling current was driven from a tip of scanning tunnelling microscope through ferromagnetic grains to ferromagnetic electrode .
transport measurements of such devices showed pronounced coulomb steps in the current - voltage characteristics .
moreover , negative differential conductance was observed .
it was further proposed that the negative differential conductance could be a consequence of a nonequilibrium spin accumulation .
such systems can be modelled theoretically by double - island devices whose two islands and the right electrode are ferromagnetic , whereas the left electrode is nonmagnetic , corresponding to the nonmagnetic tip of scanning tunnelling microscope , see the inset of fig .
[ sec3:fig4]c .
transport characteristics of such a device are displayed in fig .
[ sec3:fig4 ] for the parameters taken from ref .
the shifts of the fermi level for spin - up electrons are shown in fig .
[ sec3:fig4]a .
these shifts are different in both magnetic configurations .
generally , spin accumulation on the first island is larger than the accumulation on the second island .
the reason for this is the fact that the rate for electron tunnelling from the first island to the left lead is smaller than the rate for tunnelling of electrons to or from the second island , which is due to asymmetry of barriers .
the currents flowing through the system in the parallel and antiparallel configurations are illustrated in fig .
[ sec3:fig4]b .
moreover , negative differential conductance occurs in both magnetic configurations , however , it is more pronounced in the antiparallel configuration , as shown in fig .
[ sec3:fig4]c .
it can be seen that negative differential conductance increases with increasing the bias voltage .
the resulting tmr effect is displayed in fig .
[ sec3:fig4]d .
it is interesting to note that tmr oscillates between negative and positive values .
( color online ) differential conductance in the antiparallel configuration as a function of the bias voltage calculated for different values of the spin relaxation time @xmath261 and for @xmath262/ev .
the other parameters as well as magnetic configuration of the system are the same as in fig .
[ sec3:fig4 ] .
( after ref .
@xcite ) ]     ( color online ) tmr as a function of the bias voltage calculated for different values of spin relaxation time @xmath261 and for @xmath262/ev .
the other parameters as well as magnetic configuration of the system are the same as in fig .
[ sec3:fig4 ] .
( after ref .
@xcite ) ]    oscillations of the sign of tmr and differential conductance result from spin accumulations in both islands and are absent in the limit of fast spin relaxation ( no spin accumulation ) in the islands .
the results presented in fig .  [ sec3:fig4 ] correspond to the long spin relaxation limit .
in such a limit , some spin accumulation may occur even for a very small current flowing through the system , giving rise to ndc and oscillations of the tmr sign for small bias voltages .
however , both effects disappear when the spin relaxation time is shorter than the time between successive tunnelling events .
thus , for a finite relaxation time , one may expect absence of ndc and tmr oscillations for small voltages and the onset of these effects at larger voltages . this is because at some voltage there is a crossover from the fast to slow spin relaxation limits .
in fact such a behavior is consistent with experimental data of ref .
@xcite .
the disappearance of ndc with decreasing spin relaxation time @xmath263 is shown explicitly in fig .
[ sec3fig : g_ap ] , where the bias dependence of differential conductance is presented for different values of the spin relaxation time .
this figure clearly shows that ndc disappears when spin relaxation time decreases , in agreement with the above discussion . in the limit of fast spin relaxation time , the differential conductance is positive , although its periodic modulation still remains .    similarly , periodic oscillations of the sign of tmr also disappear with decreasing spin relaxation time @xmath264 .
this behavior is shown in fig .
[ sec3fig : tmr ] , where the bias dependence of tmr is shown for several values of @xmath264 .
first , the transitions to negative tmr disappear with decreasing @xmath264 .
the tmr becomes then positive , although some periodic modulations survive .
second , the phase of the modulations shifts by about @xmath265 when the spin relaxation varies from fast to slow limits .    in the case of analyzed systems , the negative differential conductance occurs due to nonequilibrium spin accumulation on the islands .
it is however worth noting that the negative differential conductance may also exist in single - electron devices built of nonmagnetic materials @xcite .
in the previous two sections we considered spin - dependent transport through single - electron devices , whose central electrodes ( islands ) were described by continuous ( or discrete with small level separation ) energy spectrum , and the most relevant and dominant energy scale was the electrostatic charging energy . in ultra - small metallic islands ( nanoparticles ) or in semiconducting quantum dots , the level spacing is comparable with the charging energy or even larger .
this also happens in the case of a molecule attached to metallic leads . in this limit one arrives at slightly more sophisticated single - electron devices @xcite .    in this section
we consider a fm set based on a semiconductor quantum dot coupled to ferromagnetic leads .
this model also applies to molecules attached to ferromagnetic electrodes .
single - electron transistors based on quantum dots are of current interest not only because of new and interesting physics emerging in those systems , but , more importantly , due to possible future applications and due to the possibility of manipulation of a single electron charge and a single electron spin @xcite . furthermore , quantum dots are also interesting for future applications in quantum computing @xcite .
transport properties of quantum dots coupled to nonmagnetic leads have already been extensively studied both theoretically and experimentally @xcite .
however , further interesting effects occur in the case of quantum dots coupled to ferromagnetic leads , e.g. spin accumulation , parity effect on tunnel magnetoresistance , zero - bias anomaly in the coulomb blockade regime , exchange field , splitting of the kondo anomaly , and others .
most of the works concerned theoretical description of spin - polarized transport in the weak coupling regime , as well as in the strong coupling regime , where the kondo physics emerges @xcite .
sequential transport through a single - level quantum dot coupled to ferromagnetic leads was studied for both collinear @xcite and non - collinear @xcite configurations of the electrodes magnetic moments .
spin - polarized transport in the cotunneling regime has also been addressed for collinear systems @xcite , as well as for systems magnetized non - collinearly @xcite .
furthermore , the resonant tunnelling was also considered @xcite .
quantum dots coupled to ferromagnetic leads may be realized experimentally in various ways , including ultrasmall metallic ( e.g. aluminum ) nanoparticles @xcite , single molecules @xcite , granular structures @xcite , self - assembled dots in ferromagnetic semiconductors @xcite , carbon nanotubes @xcite , and magnetic tunnel junctions @xcite .
quite recently , semiconductor quantum dots based on two - dimensional electron gas were successfully attached to ferromagnetic leads @xcite .
let us begin our discussion with the case when the dot is weakly coupled to the leads .
one can then use a perturbative approach in which the coupling is considered as a small perturbation to the unperturbed quantum dot and the leads .
transport is then dominated either by the first order or by the second order ( cotunneling ) contributions . later on we consider a more general situation when the coupling may be strong and a new kondo physics emerges at low temperatures .
schematic of a quantum dot coupled to ferromagnetic leads .
the magnetic moments of external electrodes can be aligned either in parallel or antiparallel .
the system is symmetrically biased and there is a gate voltage attached to the dot . ]
a schematic of a quantum dot coupled to ferromagnetic leads is presented in fig .
[ sec4:fig1 ] .
the magnetizations of the leads can be either parallel or antiparallel .
there is also a gate voltage attached to the dot .
the system is modelled by an anderson - like hamiltonian of the general form @xcite @xmath266 where the first and second terms describe the left and right reservoirs of noninteracting electrons , @xmath267 for @xmath268 [ see eq .
( [ eq : hamiltonian2 ] ) ] .
the third term of the hamiltonian , @xmath269 , represents the dot and includes two components : one describes noninteracting electrons in the dot level and the other represents the coulomb interaction of two electrons residing in this level , @xmath270 where @xmath271 is the correlation energy , @xmath272 is the energy of an electron in the dot with spin @xmath53 , and @xmath273 ( @xmath274 ) is the corresponding creation ( annihilation ) operator .
the position of the dot level can be tuned by the gate voltage , but is independent of the symmetrically applied transport voltage . generally , the dot level may be spin - split , for example due to a stray field of the electrodes or due to an external magnetic field .
the corresponding level splitting is denoted by @xmath275 , whereas @xmath276 is the energy of the spin - degenerate dot level .
in general , four different states of the dot are possible : empty dot ( @xmath277 ) , singly occupied dot with a spin - up ( @xmath278 ) or spin - down ( @xmath279 ) electron , and doubly occupied dot ( @xmath280 ) , where @xmath281 are the corresponding eigenfunctions .    interaction between the leads and quantum dot is incorporated in the tunnelling hamiltonian , @xmath282 , given by @xmath283 where @xmath284 are the tunnel matrix elements .
tunnelling gives rise to an intrinsic broadening @xmath285 of the dot levels , @xmath286 .
the parameters @xmath287 and @xmath288 describe contributions to the level widths due to coupling of the dot to the lead @xmath8 .
the respective contribution @xmath289 can be expressed in terms of the fermi golden rule as @xmath290 .
assuming the tunnel matrix elements @xmath291 to be independent of the wave vector @xmath292 , one can write @xmath293 the coupling parameters are usually expressed in terms of the spin polarization @xmath7 of the lead @xmath8 , defined by eq .
( [ eq : polarization ] ) , as @xmath294 , where @xmath295 . as reported in @xcite , typical values of the dot - lead coupling strength @xmath296 in the weak coupling regime are of the order of tens of @xmath297ev .    in order to investigate transport properties of the system in the whole range of parameters ,
one has to use a method which is more sophisticated than that based on the fermi golden rule .
two such techniques are commonly used : the method based on the equation of motion for the electron green functions , and the real - time diagrammatic technique .
discussion in this section is based mainly on the latter one @xcite .
this technique is based on a systematic perturbation expansion of expectation value of the current operator and the density matrix order by order in the dot - lead coupling strength @xmath296 .
even , if we limit further considerations based on this technique to the second order processes , this allows us to describe transport in the resonance regime , where the usual second order perturbation term diverges .
one also could go beyond the second order theory and describe more subtle effects .    within the real - time diagrammatic technique ,
the density matrix elements of the quantum dot , @xmath298 , are given by a kinetic equation in the liouville space @xmath299 where @xmath300 is the irreducible self - energy corresponding to transition forward in time from state @xmath301 to @xmath302 and then backward in time from state @xmath303 to @xmath304 .
when the tunnelling processes are spin - conserving and magnetic moments of the leads are collinear , the density matrix is diagonal . after performing the perturbation expansion
, one gets the following first - order ( sequential tunnelling ) and second - order ( cotunneling ) master equations : @xmath305 respectively , where the probabilities obey the normalization condition @xmath306 .    systematic perturbation expansion of electric current with respect to the coupling strength @xmath296 can be performed in a similar way , and the first - order and second - order contributions to current are given by the expressions @xcite @xmath307 where the coefficients @xmath308 and @xmath309 are the first- and second - order self - energies , modified as compared to @xmath310 , to account for the number of electrons transferred through the barriers @xcite .
transport characteristics of quantum dots are generally studied in both the linear and nonlinear response regimes @xcite .
it is thus important to distinguish between different transport regimes of quantum dots , which are sketched in fig .
[ sec4fig : regimes ] and labelled by the corresponding capital letters .     a sketch illustrating different transport regimes .
the respective regimes are separated by solid lines and labelled correspondingly . ]
first of all , by changing position of the dot level ( by the gate voltage for instance ) or applying the bias voltage , one can cross over from one regime to another .
the three regions ( a , b and a ) around zero bias correspond to the regime where sequential tunnelling is exponentially suppressed and current flows mainly due to cotunneling .
the charge state of the dot is then fixed ( strictly at zero temperature ) to zero electrons in the regime a , one electron in the regime b , and two electrons in the regime a. the first - order tunnelling processes are possible once the bias voltage is increased above the threshold voltage , allowing for finite occupation of two adjacent charge states ( zero and one for regime c , and one and two for regime c ) . in the regime d all four dot states @xmath311 are possible . by performing a particle - hole transformation , the behavior in regime a and c can be mapped to that in regime a and c , respectively .      to discuss the role of second - order processes we first briefly describe the case of a quantum dot coupled to nonmagnetic leads ( @xmath312 ) .
linear conductance as a function of the level position .
the dashed line corresponds to the first - order contribution @xmath313 , the dotted line presents the second - order conductance @xmath314 and the solid line shows the sum @xmath315 .
the different transport regimes are also specified .
the parameters are : @xmath316 , @xmath317 , and @xmath318 .
( after ref .
@xcite ) ]    figure  [ sec4fig : nonmagncond ] shows the first - order ( dashed line ) and second - order ( dotted line ) contributions to the linear conductance as well as the total conductance ( solid line ) .
the conductance is shown there as a function of the dot s level energy .
when the dot level crosses the fermi level of electrodes , there is a resonance peak in the linear conductance .
another resonance appears when @xmath319 crosses the fermi level .
the resonance peaks acquire a certain width as a result of the level broadening due to coupling to the leads ( thermal fluctuations also contribute ) .
it is interesting to note , that the second - order contribution becomes negative at resonances , which indicates that the second - order processes renormalize the first - order ( sequential ) contribution . except for resonances ,
the dot is either in the empty ( regime a ) or doubly occupied ( regime a ) state , or in the coulomb blockade regime ( regime b ) . in all these three cases
the cotunneling contribution to electric current becomes dominant .
it is also worth noting , that the second - order processes lead to renormalization of the dot level energy @xcite .      when the external electrodes are ferromagnetic , qualitatively new features appear in the transport characteristics .
these are particularly pronounced in the tunnel magnetoresistance , defined by eq .
( [ eq : tmr ] ) .
below we describe some of them .
( color online ) the first - order ( a ) and first- plus second - order ( b ) tunnel magnetoresistance as a function of the bias and gate voltages for the parameters : @xmath320 , @xmath321 , and @xmath322 .
( after ref .   @xcite),title="fig : " ]   ( color online ) the first - order ( a ) and first- plus second - order ( b ) tunnel magnetoresistance as a function of the bias and gate voltages for the parameters : @xmath320 , @xmath321 , and @xmath322 .
( after ref .
@xcite),title="fig : " ]    tunnel magnetoresistance calculated in the first order of the dot - lead coupling strength @xmath296 ( sequential transport limit ) is shown in fig .
[ sec4fig : tmr3d]a as a function of the bias and gate voltages .
it is evident that the tmr generally acquires then two different values , depending on the transport region . for the regions a ( and a ) , b and d , the tmr value is @xmath323 while for the region c ( and c ) it is @xmath324 this behavior can be accounted for as follows . for the regions a ( a ) and b the first - order linear conductance in the parallel and antiparallel configurations is @xmath325 which leads to the magnetoresistance equal to @xmath326 . to account for
the behavior of tmr in the regions c and d let us consider the zero temperature limit .
in the region c there are then three dot states taking part in transport : @xmath327 ( because of the particle - hole symmetry , the results are also applicable to regime c ) . on finds then the first - order currents in the parallel and antiparallel configurations to be @xmath328 and @xmath329 , which leads to @xmath330 .
similar analysis for the region d gives @xmath331 .
as we already know from previous sections , first - order ( sequential ) transport does not describe properly the blockade regions ( the regions a , a and b ) .
transport in these regions is expected to be strongly modified by the second - order ( cotunneling ) processes .
accordingly , the above picture is also strongly modified .
the corresponding tmr is displayed in fig .
[ sec4fig : tmr3d]b .
since the sequential tunnelling dominates transport above the threshold voltage ( outside the blockade regions , i.e. in the regions c and d ) , the total tmr is only slightly modified there as compared to that in the first - order . in the regions a ,
a and b , however , the first - order tunnelling processes are exponentially suppressed and the current is dominated by the second - order ( cotunneling ) processes . as follows from fig .
[ sec4fig : tmr3d]b , cotunneling has then a significant influence on tmr .
apart from this , tmr in the regions a and a behaves differently from that in the region b.     the total linear tunnel magnetoresistance ( solid line ) as a function of the level position .
the dashed line represents the first - order tunnel magnetoresistance , whereas the dotted line corresponds to the jullire s value .
the parameters are @xmath332 , @xmath321 , and @xmath322 .
( after ref .
@xcite ) ]    let us look now in more details at the vertical cross - section of fig .
[ sec4fig : tmr3d]b along the zero bias ( linear response ) .
the corresponding linear tmr as a function of the level position ( gate voltage ) is shown in fig .
[ sec4fig : regimeaba ] . by changing position of the dot level
, one crosses over from the regime a through the regime b to the regime a. for comparison , the jullire s value of tmr , as well as the first - order contribution to tmr ( denoted as @xmath333 ) are shown there .
note , that the first - order tmr is independent of the gate voltage and is equal to a half of the jullire s value .
first of all , one can note that the second - order processes modify the tmr substantially . unlike the first - order term , the total tmr ( first - order plus second - order contributions ) does depend on the level position . the interesting feature of tmr shown in fig .  [ sec4fig :
regimeaba]b is a strong parity effect .
the tmr reaches maximum when there is an even number of electrons on the dot ( zero for the regime a or two for the regime a ) , and has minimum for an odd ( one in the regime b ) number of electrons .
a universal feature is that for even electron number , tmr exactly coincides with the jullire s value .
the system behaves then like a single ferromagnetic tunnel junction .
such a situation can take place when tunnelling processes are coherent .
the only second - order processes that contribute to conductance in these regions are the non - spin - flip cotunneling processes , in which the electron spin is conserved .
such processes indeed are fully coherent .
the corresponding cotunneling rates are proportional to the product of the density of states of the left and right leads , thus , one can express the second - order linear conductance as @xmath334 for the parallel and antiparallel configuration , respectively . as a consequence ,
the tmr is then equal to that of a single planar ferromagnetic tunnel junction @xmath335 because of the particle - hole symmetry , the same result can be obtained for regime a.    the situation becomes , however , more complex for an odd number of electrons on the dot .
apart from the non - spin - flip processes , there are also spin - flip cotunneling processes that change spin of the dot .
such second - order processes give rise to spin relaxation on the dot and lead to reduction of tmr in the regime b. the dependence of tmr on the position of the energy level in the regime b reflects the relative relation of the spin - flip to non - spin - flip cotunneling processes .
the minimum value of tmr appears for @xmath336 , see fig .
[ sec4fig : regimeaba]b , where one finds @xmath337    the total current ( a ) in the parallel ( solid line ) and antiparallel ( dashed line ) magnetic configurations as a function of the bias voltage in units of @xmath338 .
part ( b ) shows the first - order contribution to the tmr ( dashed line ) and the total tmr ( solid line ) .
the parameters are : @xmath332 , @xmath339 , @xmath321 , and @xmath322 .
( after ref .
@xcite ) ]    electric current flowing through the system in the nonlinear response regime is shown in fig .
[ sec4fig : regimedbd]a for both magnetic configurations and for a symmetric anderson model ( @xmath336 ) .
the corresponding tmr is also shown there , see fig .
[ sec4fig : regimedbd]b .
the contribution due to the second - order tunnelling processes is significant ( or even dominant ) in the coulomb blockade regime .
the sequential tmr is shown by a dashed line in fig .
[ sec4fig : regimedbd]b .
one can see that sequential tmr is constant as a function of the bias voltage .
however , this is not a universal behavior and occurs only for symmetric anderson model .
generally , first - order tmr depends on the applied bias voltage , as illustrated in fig  [ sec4fig : tmr3d]a .
the second - order processes lead to a strong and nontrivial dependence of the total tmr on the bias voltage . for large values of the transport voltage ( regime d ) ,
the first - order tunnelling processes dominate , and therefore the total tmr is only slightly modified as compared to that in the sequential transport regime
. however , for voltages below the threshold voltage ( regime b ) , the cotunneling processes lead to a strong bias dependence of tmr in the coulomb blockade regime , and to a deep minimum of tmr in the zero - bias limit . in the case of metallic islands ,
cotunneling processes usually lead to an enhancement of the tmr effect in the coulomb blockade regime @xcite . here
, we have the opposite situation , i.e. , suppression of the effect . at low bias voltage , @xmath340 ,
the single - barrier spin - flip processes reduce the tmr .
this is however no longer the case for nonlinear response regime , @xmath341 , where the spin accumulation diminishes the amount of spin - flip processes and the tmr increases .
consequently , the tmr effect in regime b increases with rising the bias voltage within the limits @xmath342 .
the minimal value is reached at @xmath343 and @xmath336 , whereas the maximal value is approached for bias voltages large as compared to thermal energy but still far away from the onset of sequential tunnelling .
( color online ) the total tunnel magnetoresistance as a function of the bias and gate voltages in the presence of external magnetic field for ( a ) @xmath344 and ( b ) @xmath345 .
the parameters are : @xmath332 , @xmath321 , and @xmath322.,title="fig : " ]   ( color online ) the total tunnel magnetoresistance as a function of the bias and gate voltages in the presence of external magnetic field for ( a ) @xmath344 and ( b ) @xmath345 .
the parameters are : @xmath332 , @xmath321 , and @xmath322.,title="fig : " ]    new features of transport characteristics appear when the spin degeneration of the dot level is lifted , @xmath346 , e.g. due to an external magnetic field applied to the system .
it turns out that a finite zeeman splitting , @xmath347 , changes the transport characteristics substantially .
figure  [ sec4fig : tmr2dfield]a illustrates the gate and bias voltage dependence of tunnel magnetoresistance for @xmath348 .
first of all , a finite zeeman splitting affects the tmr mainly in the regimes a ( a ) and b , where the second - order processes dominate . furthermore , tmr exhibits a distinctively different behavior in the regions a and a. it is enhanced in the case of empty dot , and reduced in the case of doubly occupied dot ( it may even become negative in the latter case ) .
this effect depends on the orientation of applied magnetic field ; when the field is applied in the opposite direction , @xmath349 , there is an enhancement of tmr in regime a and suppression in regime a , as illustrated in fig .
[ sec4fig : tmr2dfield]b .
thus , by changing the sign of the zeeman splitting ( magnetic field orientation ) or applying a gate voltage , which changes the dot occupation , it is possible to reduce or enhance the tmr effect considerably .
another interesting feature of tmr displayed in fig .
[ sec4fig : tmr2dfield]a appears in the region b. there is a strong asymmetry of tmr with respect to the bias reversal .
the tmr is decreased for positive and increased for negative bias voltage .
this can also be seen in fig .
[ sec4fig : tmr2dfield]b , but for positive transport voltage tmr is then larger than for negative bias voltage .
the crossover between those two values of tmr takes place roughly at the zero bias . as a consequence , by changing the bias voltage in the small range one
can substantially tune magnitude of the tmr effect .
as we know from previous subsection , the most interesting features due to ferromagnetism of the electrodes , e.g. the zero - bias anomaly , occur in the blockade regime , where sequential tunnelling processes play a minor role .
this anomaly appears only in the case when the dot is singly occupied at equilibrium ( @xmath350 , @xmath351 ) and the system is in a deep coulomb blockade regime ; latexmath:[$\gamma , k_{\rm b}t \ll    sequential tunnelling is exponentially suppressed , and cotunneling gives the dominant contribution to electric current @xcite .
therefore , the sequential tunnelling can be completely ignored . in this subsection
we look more carefully at these features taking into account only second - order contribution .    in order to calculate the cotunneling current in the deep blockade regime
one can employ a simplified second - order perturbation theory and master equation for the occupation probabilities .
the rate of a cotunneling processes from lead @xmath8 to lead @xmath353 , which change the dot state from @xmath281 to @xmath354 , can be written as @xmath355 with @xmath356 and @xmath357 denoting the energies of initial and final states , @xmath358 being the state of the system with an electron in the lead @xmath8 and the dot in state @xmath281 , whereas @xmath359 is a virtual state with @xmath360 denoting the corresponding energy . among different cotunneling processes one can distinguish the single - barrier ( @xmath361 ) and double - barrier ( @xmath362 ) cotunneling as well as spin - flip ( @xmath363 ) and non - spin - flip ( @xmath364 ) cotunneling .
the spin - flip processes change the spin state of the dot , whereas the non - spin - flip processes do not change the dot state .
the current flows through the system due to double - barrier cotunneling processes . on the other hand
, the single - barrier processes do not contribute directly to electric current , however , they can change the dot occupations , and this way also influence the current .    the cotunneling current flowing through the system from the left to right lead
is given by @xmath365,\ ] ] where @xmath366 denotes the corresponding occupation probability .
the probabilities @xmath366 can be found from the master equation , @xmath367 \,,\ ] ] together with the normalization condition @xmath368 .      in order to elucidate and understand the anomalous behavior of tmr in the coulomb blockade regime , we show in fig .
[ sec4fig : regimeb]a the differential conductance in the small bias regime for both parallel and antiparallel configurations and for several temperatures .
figure  [ sec4fig : regimeb]b displays the corresponding tmr effect .
first of all , the tmr effect in the regime b for @xmath369 increases with lowering temperature and approaches the jullire s value , whereas the minimum at zero bias does not depend on temperature .
the differential conductance in the parallel alignment has characteristics typical of the cotunneling regime , with a smooth parabolic dependence on the bias voltage . for antiparallel configuration ,
on the other hand , differential conductance has a local maximum at zero bias , followed by local minimum with increasing bias , as illustrated in fig .
[ sec4fig : regimeb]a .
this zero - bias anomaly stems from the interplay of the spin - flip and non - spin - flip single - barrier and double - barrier cotunneling processes @xcite .
the minimum in the tmr effect is a direct consequence of this anomalous behavior of differential conductance in the antiparallel configuration .
the differential conductance ( a ) for the parallel and antiparallel configurations and the tunnel magnetoresistance ( b ) as a function of the bias voltage for different values of temperature .
the maximum in differential conductance for antiparallel configuration at zero bias is clearly demonstrated .
the other parameters are the same as in fig .
[ sec4fig : regimedbd ] .
figure was generated using the scheme for the perturbation expansion in the coulomb blockade regime .
( after ref .
@xcite ) ]    the zero - bias anomaly in the cotunneling regime is qualitatively similar to the anomaly due to the kondo effect , which occurs in the strong coupling limit @xcite .
there are , however , some distinct differences .
first of all , processes responsible for the zero - bias anomaly in the cotunneling regime are of the second order in tunnelling processes , while these leading to the kondo effect are of higher order .
the conductance in the cotunneling regime is much smaller than in the kondo regime , where almost perfect transmission ( @xmath370 ) through the dot is possible owing to the kondo peak in the density of states at the fermi level .
furthermore , the kondo peak occurs at temperatures lower than the so - called kondo temperature , @xmath371 , and exists also in the parallel configuration @xcite .    to understand the mechanism of the zero - bias anomaly it is crucial to distinguish between different types of cotunneling processes : the single - barrier cotunneling processes shown in fig .
[ sec5fig : zba2]a and double - barrier cotunneling processes illustrated in fig .
[ sec5fig : zba2]b . both single - barrier and double - barrier processes
can be either spin - flip or non - spin - flip ones . the current flows due to double - barrier cotunneling , whereas the single - barrier cotunneling can influence the current in an indirect way , by changing the spin state of the dot . in the antiparallel configuration
, there is a finite spin accumulation on the dot , as presented in fig .
[ sec5fig : zba2]c .
single - barrier ( a ) and double - barrier ( b ) cotunneling processes , and the occupation probabilities for spin - up and spin - down electrons in the antiparallel configuration ( c ) . the parameters are @xmath372 , @xmath373 , @xmath374 , and @xmath322 .
( after ref .
@xcite ) ]    the different occupation probabilities for spin - up and spin - down electrons appear due to spin asymmetry in tunnelling processes . in equilibrium ,
both rates are equal and there is no spin accumulation , @xmath375 .
when a bias voltage is applied , and the system is in the antiparallel configuration , the relative amount of single - barrier cotunneling is diminished as compared to the double - barrier cotunneling .
this is because the rate of single - barrier cotunneling is proportional to thermal energy , whereas that of double - barrier cotunneling is proportional to the bias voltage . in the nonlinear response regime , the magnetic state of the dot
is mainly determined by the spin - flip processes that transfer an electron from the left to the right leads . the one shown in fig .
[ sec5fig : zba2]b changes the dot spin from @xmath376 to @xmath377 . because the rate of this process is proportional to a product of densities of states for majority electrons , the corresponding rate is larger than that of the other process that changes the dot spin from @xmath377 to @xmath376 , where only the minority spins are involved .
this results in a nonequilibrium spin accumulation , @xmath378 , that increases with increasing voltage , as shown in fig .
[ sec5fig : zba2]c .
the initial state for the dominant spin - flip cotunneling process that contributes to current is @xmath376 , as sketched in fig .
[ sec5fig : zba2]b .
thus , the conductance is diminished by spin accumulation .
this is the mechanism by which spin accumulation gives rise to nonzero tunnel magnetoresistance effect , @xmath379 , see fig .
[ sec4fig : regimeb]b .
since the spin accumulation reduces electronic transport , any spin - flip process that reduces the spin accumulation should enhance the conductance . in particular , single - barrier spin - flip cotunneling is an example of such a process .
as pointed above , the rate of single - barrier processes scales with @xmath380 while that of double - barrier cotunneling is proportional to @xmath381 , which explains the mechanism of the zero - bias anomaly . at low bias voltage , @xmath382
, single - barrier spin - flip processes play a significant role
they decrease the spin accumulation opening this way the system for the fastest double - barrier cotunneling . as a consequence ,
the current increases relatively fast with applied bias , leading to a maximum in differential conductance . for @xmath383 ,
on the other hand , the relative role of single - barrier processes is negligible as compared to double - barrier cotunneling , and the conductance is reduced .
thus , the interplay between the rates of double - barrier and single - barrier cotunneling processes leads to the maximum in the differential conductance at the zero bias .
intrinsic spin relaxation in the dot can result , for instance , from spin - orbit interaction , coupling of the electron spin to nuclear spins , etc .
however , we will not consider a particular microscopic mechanism of the intrinsic spin - flip processes , but simply assume that the spin - flip relaxation is described by a spin - relaxation time @xmath263 , and is taken into account _ via _ a relaxation term in the appropriate master equation for the occupation probabilities @xcite , @xmath384 where @xmath385 , @xmath386 denotes the probability that the dot is occupied by a spin-@xmath53 electron , and @xmath387 is the cotunneling rate from lead @xmath388 to lead @xmath389 with a change of the dot spin from @xmath53 to @xmath390 ( @xmath391 ) .
the last term describes the spin relaxation processes , which in the case of spin - degenerate dot level reduces to @xmath392 .    as before , ( see section 2 ) one can distinguish between the fast and slow spin relaxation limits .
the former ( latter ) limit corresponds to the situation when the time between successive cotunneling events , @xmath393 , is significantly longer ( shorter ) than the intrinsic spin relaxation time @xmath263 .
a typical spin relaxation time for quantum dots can be relatively long , up to @xmath297s @xcite .
on the other hand , the time between successive cotunneling events can be estimated taking into account the fact that the rate of spin - flip cotunneling is generally larger than that of non - spin - flip cotunneling ( for a finite parameter @xmath271 ) .
assuming @xmath394 , one then finds @xmath395 with @xmath396 and @xmath397 . assuming typical parameters @xcite , one can roughly estimate @xmath393 to range from @xmath398 ns to @xmath399 ns .     the differential conductance in the antiparallel configurations ( a ) and tunnel magnetoresistance ( b ) as a function of the bias voltage for different spin relaxation @xmath400 .
the parameters are : @xmath401 , @xmath402 , @xmath403 , and @xmath404 .
( after ref .
@xcite ) ]    the anomalous behavior of differential conductance in the antiparallel configuration results from spin asymmetry of the tunnelling processes , which gives rise to spin accumulation , @xmath405 .
the intrinsic spin - flip processes in the dot reduce the spin accumulation and this way suppress the anomaly .
this behavior is displayed in fig .
[ sec5fig1](a ) , where different curves correspond to different values of the parameter @xmath8 defined as @xmath400 .
thus , @xmath406 describes the case with no intrinsic spin relaxation , whereas the curves corresponding to nonzero @xmath8 describe the influence of intrinsic relaxation processes .
first of all , one can note that small amount of intrinsic spin - flip processes enhances the zero - bias anomaly [ see the curve for @xmath407 in fig .
[ sec5fig1](a ) ] .
this is because such processes play then a role similar to that of single - barrier spin - flip cotunneling . in the case of fast spin relaxation ,
on the other hand , the spin accumulation is suppressed and the anomaly disappears , as can be seen in fig .  [ sec5fig1](a ) for @xmath408 .
it is also worth noting that spin - flip processes in the dot enhance the overall conductance in the antiparallel configuration . in the parallel configuration ,
however , the differential conductance does not depend on intrinsic relaxation .    modifying the conductance in the antiparallel configuration ,
spin relaxation processes suppress the tunnel magnetoresistance , as shown in fig .
[ sec5fig1](b ) . since the zero - bias maximum in conductance is suppressed in the fast spin relaxation limit , the corresponding dip in tmr at small voltages is also suppressed by the intrinsic relaxation processes .
more specifically , the dip in tmr broadens with increasing @xmath8 and disappears in the limit of fast relaxation ( see the curve for @xmath408 ) .
an interesting feature of tmr in the presence of spin - flip scattering in the dot is the crossover from positive to negative values when @xmath8 increases , as illustrated in fig .
[ sec5fig1](b ) .
thus , the difference between conductances in the parallel and antiparallel magnetic configurations persists even for fast spin relaxation in the dot , contrary to the sequential tunnelling regime , where such a difference disappears @xcite .
this seemingly counterintuitive behavior can be understood by taking into account the following two facts : ( i ) absence of spin accumulation in the dot for fast spin relaxation ( @xmath409 ) , and ( ii ) difference in the fastest cotunneling processes contributing to the current in the two magnetic configurations .
the fastest double - barrier cotunneling processes involve only the majority - spin electrons of the two leads  thus , in the parallel configuration the fastest cotunneling processes are the non - spin - flip ones .
they take place either _ via _ the empty - dot virtual state ( for one orientation of the dot spin ) or _ via _ the doubly occupied dot virtual state ( for the second orientation of the dot spin ) .
the dominant contribution to the current is then proportional to @xmath410 . on the other hand , in the antiparallel magnetic configuration
the fastest cotunneling processes are the spin - flip ones , which can occur only for one particular orientation of the dot spin . however , for this spin orientation cotunneling can take place _ via _ both empty and doubly occupied dot virtual states .
the corresponding dominant contribution to electric current is then proportional to @xmath411 ^ 2=1/\varepsilon^2 + 1/(\varepsilon + u)^2 -2/[\varepsilon ( \varepsilon + u)]$ ] .
it is thus clear that the difference in currents flowing through the system in the antiparallel and parallel configurations is equal to @xmath412 $ ] , which results from the interference term . since @xmath413 and @xmath414
, this interference contribution is positive . as a result ,
the current in the antiparallel configuration is larger than the current in the parallel configuration .
the minimum in tmr at zero bias in the case of a symmetric anderson model can be expressed as @xmath415 whereas for @xmath416 and @xmath417 one finds @xmath418 }    { 2(1-p^2)(3-p^2)k_{\rm b}t\gamma^2+(3+p^2)\varepsilon^2    h/\tau_{\rm sf}}.\ ] ] the latter formula approximates the value of tmr corresponding to the bias voltage at which the differential conductance has a local minimum . in the slow spin relaxation limit
one finds @xmath419 , and @xmath420 .
however , in the limit of fast spin relaxation tmr becomes negative and is given by @xmath421      the discussion up to now was limited to the case of degenerate dot level .
the situation changes when @xmath422 , e.g. , due to an external magnetic field .
the level splitting is described by the parameter @xmath423 , where the magnetic field is assumed to be along the magnetic moment of the left electrode . in fig .
[ sec5fig4 ] we show the bias voltage dependence of the differential conductance in the parallel and antiparallel configurations for different values of parameter @xmath8 .
the differential conductance in the nonlinear response regime for different spin relaxation in the parallel ( a ) and antiparallel ( b ) configurations .
the parameters are : @xmath401 , @xmath424 , @xmath425 , @xmath403 and @xmath322 .
( after ref .
@xcite ) ]    in the limit of no intrinsic spin relaxation in the dot ( solid line in fig .  [ sec5fig4 ] ) and at low bias voltage , the dot is occupied by a spin - up electron and the current flows mainly due to non - spin - flip cotunneling .
the spin - flip cotunneling processes are suppressed for @xmath426 , which results in the steps in differential conductance at @xmath427 .
the suppression of spin - flip inelastic cotunneling was recently used as a tool to determine the spectroscopic @xmath428-factor @xcite .
when @xmath429 becomes larger than @xmath430 , spin - flip cotunneling is allowed , consequently the conductance increases .
however , there is a large asymmetry of differential conductance in the antiparallel configuration with respect to the bias reversal . to understand this asymmetry , it is crucial to realize that when the splitting @xmath431 is larger than @xmath432 , the single - barrier spin - flip cotunneling processes can occur only when the dot is occupied by a spin - down electron
thus , the single - barrier processes can assist the fastest double - barrier cotunneling processes , but only for positive bias .
this is because the fastest processes can occur when the dot is occupied by a spin - down electron for negative bias and by a spin - up electron for positive bias , leading to larger conductance for positive than for negative bias voltage .
no such asymmetry occurs in the parallel configuration [ see fig .  [ sec5fig4](a ) ] , as now the system is fully symmetric with respect to bias reversal .
the situation changes when intrinsic spin - flip relaxation processes occur in the dot . for the parameters assumed in fig .
[ sec5fig4 ] , the spin relaxation in the dot affects the conductance only for @xmath433 , while for @xmath434 the conductance is basically independent of @xmath8 ( see fig .  [ sec5fig4 ] ) .
this is because for @xmath434 and @xmath435 , the dot is predominantly occupied by a spin - up electron and the transitions to the spin - down state due to relaxation processes are energetically forbidden . as a consequence ,
the current flows mainly due to non - spin - flip cotunneling , irrespective of spin relaxation time .
this scenario holds for both magnetic configurations of the system .
when @xmath433 , the spin - flip cotunneling processes can take place and the dot can be either in the spin - up or spin - down state . in the parallel configuration , fig .
[ sec5fig4](a ) , the conductance is slightly reduced by the spin - flip relaxation processes . this can be understood by realizing the fact that the fastest non - spin - flip cotunneling processes in the parallel configuration are more probable when the dot is occupied by a spin - down electron than by a spin - up one [ due to smaller energy denominator , see eq .
( [ eq : cotunnelingrate ] ) ] .
since the spin - down state ( as that of larger energy ) relaxes relatively fast to the spin - up state ( which has definitely smaller energy ) , this leads to a reduction in the conductance . on the other hand ,
the differential conductance in the antiparallel configuration is enhanced by the relaxation processes for positive bias and diminished for negative bias voltages .
consider first the situation for positive bias .
as already discussed above for @xmath406 , an important role in that transport regime is played by the single - barrier spin - flip cotunneling processes , which open the system for the fast double - barrier cotunneling by reversing spin of the dot from the spin - down to the spin - up state .
the relaxation processes play a role similar to that of the single - barrier cotunneling , and lead to a certain increase in the conductance . for negative bias voltage , in turn , the fast double - barrier cotunneling processes occur when the dot is occupied by a spin - down electron .
the probability of such events is decreased by spin relaxation , leading to a reduced conductance .
an interesting consequence of the enhancement ( reduction ) of the differential conductance for positive ( negative ) bias voltage is an increase of the asymmetry with respect to the bias reversal  see the curve for @xmath408 in fig .
[ sec5fig4](b ) .
an interesting situation occurs when the quantum dot is coupled asymmetrically to the left and right leads @xmath436 .
the differential conductance for a system with one electrode nonmagnetic and the other one made of a ferromagnet with large spin polarization ( in the following referred to as strong ferromagnet ) is shown in fig .
[ sec5fig5 ] for the case when the dot is described by an asymmetric anderson model ( @xmath437 ) .
the differential conductance in the nonlinear response regime for the asymmetric anderson model , @xmath438 , for different spin relaxation .
the parameters are : @xmath401 , @xmath439 , @xmath440 and @xmath441 .
( after ref .
@xcite ) ]    consider first the situation in the absence of intrinsic spin relaxation in the dot ( solid curves in fig .
[ sec5fig5 ] ) . when @xmath442 , the influence of single - barrier cotunneling can be neglected @xcite .
the cotunneling processes which transfer charge from one lead to another take place _ via _ two possible virtual states
empty dot ( an electron residing in the dot tunnels to one of the leads and another electron from the second lead enters the dot ) and doubly occupied dot ( an electron of spin opposite to that in the dot enters the dot and then one of the two electrons leaves the dot ) .
consider first positive bias ( @xmath443 , electrons flow from right to left , i.e. , from normal metal to strong ferromagnet ) , and assume for clarity of discussion that the strong ferromagnet is a half - metallic one with full spin polarization ( only spin - up electrons can then tunnel to the left lead ) .
when a spin - down electron enters the dot , it has no possibility to leave the dot for a long time .
the allowed cotunneling processes occur then _ via _ doubly occupied dot virtual states . in the absence of intrinsic spin relaxation in the dot , the only processes which can reverse the dot spin are the single - barrier cotunneling ones , which however play a minor role when @xmath444 .
thus , the current flows due to non - spin - flip cotunneling _ via _ doubly - occupied dot virtual states , whereas cotunneling through empty - dot virtual states is suppressed .
the situation is changed for negative bias ( electrons flow from strong ferromagnet to normal metal ) .
now , the dot is mostly occupied by a spin - up electron , which suppresses cotunneling _ via _ doubly - occupied dot virtual state and the only contribution comes from cotunneling _ via _ empty - dot virtual state .
the ratio of cotunneling rates through the empty dot and doubly occupied dot virtual states is approximately equal to @xmath445^{-2}$ ] . in the situation presented in fig .
[ sec5fig5 ] one finds @xmath446 .
accordingly , the conductance for negative bias is much larger than for positive bias voltage .    when @xmath447 becomes of the order of @xmath432 or smaller , the rate of single - barrier cotunneling is of the order of the rate of double - barrier cotunneling .
therefore , the single - barrier processes can play an important role in transport .
more precisely , single - barrier cotunneling processes can reverse spin of an electron in the dot and thus can open the system for the fast cotunneling processes
.    intrinsic spin - flip processes in the dot have similar influence on electronic transport as in the case discussed in the previous section .
as before , relaxation processes remove the asymmetry with respect to the bias reversal and suppress the zero - bias anomaly .
thus , the diode - like behavior can appear only in the limit of slow spin relaxation , and is suppressed in the limit of fast spin relaxation , as shown in fig .
[ sec5fig5 ] by the curves corresponding to @xmath408 .
new features of transport characteristics occur when the leads magnetizations form an arbitrary noncollinear magnetic configuration @xcite .
such a configuration can be controlled by a weak external magnetic field , weak enough to neglect the corresponding zeeman splitting of the dot level and weaker than the effective exchange field exerted on the dot by ferromagnetic leads .
in particular , the differential conductance is significantly modified by the exchange field in noncollinear configurations .
this also applies to the tmr effect , which for arbitrarily aligned leads magnetizations can be defined as , @xmath448/i(\varphi)$ ] , where @xmath449 is an angle between the leads magnetic moments .
the system , whose spin moments @xmath450 and @xmath451 of the left and right lead , respectively , form an arbitrary configuration is shown in fig .  [ sec6fig:1 ] .
( color online ) schematic of a quantum dot coupled to ferromagnetic leads with non - collinearly aligned magnetizations .
the net spin moments of the left @xmath450 and right @xmath451 lead form an angle @xmath449 .
there is a symmetric bias voltage applied to the system . ]
the corresponding hamiltonian , written for the quantization axis of the dot equivalent to that of the left lead , is given by eq .
( [ sec4eq : qdhamiltonian ] ) , where the tunnel hamiltonian can be decomposed into two terms @xmath452 .
the term @xmath453 describes tunnelling processes between the left electrode and the dot and takes the same form as for collinear configurations , while @xmath454 describes tunnelling processes between the dot and right lead , and acquires the form @xmath455 \,.\end{aligned}\ ] ]    the real - time diagrammatic method described at the beginning of this section for collinear magnetic configurations can be extended to include the effects due to noncollinearity of the leads magnetic moments .
the density matrix of the quantum dot for an arbitrary magnetic configuration is given by @xmath456 the diagonal elements of the density matrix correspond to the respective occupation probabilities , while the off - diagonal elements @xmath457 and @xmath458 describe the dot spin @xmath459 , with @xmath460 , @xmath461 , and @xmath462 .
the density matrix elements can be determined from the corresponding kinetic equation , which in the steady state and for spin - degenerate dot level can be written as @xcite @xmath463 by expanding the self - energies and density matrix elements , one can calculate the sequential and cotunneling current .    in the following we focus on transport in the coulomb blockade regime , @xmath464 , and
when the dot is singly occupied , @xmath465 .
> from the discussion above we know that there is an anomalous behavior of the differential conductance in the small bias regime for antiparallel magnetic configuration .
now , we will consider how this anomaly changes when magnetic configuration varies continuously from antiparallel to parallel alignment .
we consider first the symmetric anderson model , @xmath336 .
the exchange field vanishes then for an arbitrary magnetic configuration , and both differential conductance and tmr change monotonically when going from antiparallel to parallel magnetic configurations .
however , for @xmath466 , the effects of exchange field become important and lead to nontrivial behavior of transport characteristics , as will be discussed later .
( color online ) differential conductance @xmath467 ( a ) and tunnel magnetoresistance tmr ( b ) as a function of the bias voltage @xmath60 for indicated values of the angle @xmath449 between magnetic moments , calculated for symmetric anderson model .
the parameters are : @xmath468 , @xmath469 , @xmath403 , and @xmath470 .
( after ref .
@xcite ) , title="fig : " ]   ( color online ) differential conductance @xmath467 ( a ) and tunnel magnetoresistance tmr ( b ) as a function of the bias voltage @xmath60 for indicated values of the angle @xmath449 between magnetic moments , calculated for symmetric anderson model .
the parameters are : @xmath468 , @xmath469 , @xmath403 , and @xmath470 .
( after ref .
@xcite ) , title="fig : " ]    the differential conductance @xmath471 as a function of the bias voltage is shown in fig .
[ sec6fig:2](a ) for several values of the angle @xmath449 .
the conductance varies monotonically with the angle between leads magnetic moments .
when the angle increases from zero to @xmath265 , the anomaly emerges at small values of @xmath449 and its relative height increases with increasing angle , reaching a maximum value at @xmath472 .
the exchange field for a symmetric anderson model is negligible , and the average dot spin tends to zero in the linear response regime .
the angular dependence of the linear conductance can be then expressed as @xmath473 \,.\ ] ] the variation of @xmath471 with @xmath449 at low bias is thus characterized by the factor @xmath474 , which leads to maximum ( minimum ) conductance in the parallel ( antiparallel ) magnetic configuration .
such behavior is typical of a normal spin - valve effect .
the bias dependence of the associated tmr is shown in fig .
[ sec6fig:2](b ) for several values of the angle @xmath449 .
the zero - bias anomaly in the differential conductance [ see fig .
[ sec6fig:2](a ) ] leads to the corresponding anomaly ( dip ) in tmr at small bias voltages .
the dip in tmr decreases when magnetic configuration departs from the antiparallel alignment , and eventually disappears in the parallel configuration .
the variation of tmr with the angle is monotonic , similarly as the angular variation of the differential conductance .
the dependence of tmr on the angle @xmath449 at zero bias is given approximately by the formula @xmath475 now , the angular dependence of tmr is governed by @xmath476 , which gives maximum tmr in the antiparallel configuration and zero tmr in the parallel one .
( color online ) the differential conductance @xmath467 ( a ) and tunnel magnetoresistance tmr ( b ) as a function of the bias voltage @xmath60 for indicated values of the angle @xmath449 for asymmetric anderson model @xmath477 and @xmath403 .
the other parameters are the same as in fig .
[ sec6fig:2 ] .
( after ref .
@xcite ) , title="fig : " ]   ( color online ) the differential conductance @xmath467 ( a ) and tunnel magnetoresistance tmr ( b ) as a function of the bias voltage @xmath60 for indicated values of the angle @xmath449 for asymmetric anderson model @xmath477 and @xmath403 .
the other parameters are the same as in fig .
[ sec6fig:2 ] .
( after ref .
@xcite ) , title="fig : " ]    the transport characteristics , described above for the symmetric anderson model , are significantly modified when the model becomes asymmetric , i.e. when @xmath466 .
this can be realized by shifting the dot level position by a gate voltage applied to the dot .
different contributions to the exchange field do not cancel then , and the resulting effective exchange field becomes nonzero . and this exchange field has a significant influence on transport properties .
the strength of effective exchange field grows with deviation of the anderson model from the symmetric one , described quantitatively by @xmath478 ( with @xmath479 for the symmetric model ) .
the leading contribution to the exchange field comes from the first - order diagrams of the perturbation expansion in terms of the real - time diagrammatic technique , whereas the current flows due to the second - order tunnelling processes .
this leads to two different time scales which determine transport characteristics .
transport properties are thus a result of the interplay of the first- and second - order processes .
although the first - order processes do not contribute directly to current , they influence transport _ via _ modification of the dot spin .
the differential conductance is shown in fig .
[ sec6fig:3](a ) as a function of bias voltage for an asymmetric anderson model , @xmath480 , and for different values of @xmath449 .
the low - bias differential conductance becomes enhanced in a certain range of the angle @xmath449 , leading to a nonmonotonic dependence of the conductance @xmath471 on the angle between the leads magnetizations .
this nonmonotonic behavior is also visible in tmr , see fig .
[ sec6fig:3](b ) .
there is also a range of the angle @xmath449 , where tmr changes sign and becomes negative in the small bias region , i.e. the corresponding conductance is larger than that in the parallel configuration .
the most characteristic features of transport characteristics in the presence of exchange field are the enhanced differential conductance at a non - collinear alignment and its rapid drop when the system approaches the antiparallel configuration . the key role in this behavior is played by the first - order processes giving rise to the exchange field .
these processes lead to the precession of spin in the dot , which facilitates tunnelling processes and leads to an increase in the conductance as compared to that in the parallel configuration .
when the configuration becomes close to the antiparallel one , the first - order processes become suppressed and the conductance drops to that for antiparallel alignment .
the nonmonotonic behavior of the differential conductance with @xmath449 leads to a nonmonotonic dependence of tmr .
the effects due to exchange field give rise to a local minimum in tmr at a non - collinear magnetic configuration .
moreover , in this transport regime tmr changes sign and becomes negative .
when the magnetic configuration is close to the antiparallel one , tmr starts to increase rapidly reaching maximum for @xmath481 .
the negative tmr and its sudden increase when the configuration tends to the antiparallel one are a consequence of the processes leading to nonmonotonic behavior of the differential conductance , as described above .
to understand more intuitively the above presented behavior of the differential conductance and tmr at low bias voltage and close to the antiparallel configuration , one should consider two different time scales .
one time scale , @xmath482 , is established by the virtual first - order processes responsible for the spin precession due to exchange field , @xmath483 the second time scale , @xmath393 , is associated with second - order processes which drive the current through the system . at low temperature and low bias voltage ,
the cotunneling rate can be expressed as @xmath484 the rate @xmath393 depends linearly on the applied voltage , whereas @xmath482 is rather independent of @xmath60 . as a consequence , at low bias and for non - collinear configuration
, the exchange field plays an important role leading to a nonmonotonic dependence of differential conductance on the angle between the leads magnetizations . when magnetic configuration is close to the antiparallel one , the spin precession rate is deceased ( @xmath485 ) and , at certain angle , the rate of spin precession becomes comparable to the cotunneling rate .
this gives rise to a sudden drop ( increase ) in differential conductance ( tmr ) .
we note that the nonmonotonic dependence of differential conductance and magnetoresistance has also been observed in quantum dots in the strong coupling limit @xcite .
up to now we discussed theoretical aspects of electronic transport through the simplest quantum dots , i.e. the dots with only one orbital level . real dots
are usually more complex and their electronic spectrum includes many orbital levels taking part in electronic transport .
this may significantly change transport characteristics and also may lead to qualitatively new features @xcite . in this section
we will consider some of the new effects in transport through multi - dot systems .
the schematic of a two - level quantum dot coupled to ferromagnetic leads with collinear magnetizations is shown in fig .
[ sec7fig:1 ] .
( color online ) energy diagram of a two - level quantum dot coupled to ferromagnetic leads . for clarity reasons the energy diagram is shown here for @xmath486 .
the leads magnetizations can form either parallel or antiparallel configurations .
the arrows indicate the net spin of the leads . ]
the quantum dot is described by the following hamiltonian @xmath487 where @xmath488 is the particle number operator , @xmath489 , @xmath490 ( @xmath491 ) is the creation ( annihilation ) operator of an electron with spin @xmath53 on the @xmath229th level ( @xmath492 ) , and @xmath493 is the corresponding single - particle energy .
the on - level coulomb repulsion between two electrons of opposite spins is described by @xmath271 , whereas the inter - level repulsion energy is denoted by @xmath494 .
the forth term in eq .
( [ sec7eq : hamiltonian ] ) describes the zeeman energy , with @xmath495 being the zeeman splitting of the energy levels ( @xmath496 is an external magnetic field along the magnetic moment of the left electrode ) . to present the main features of transport characteristics ,
it is convenient to introduce the level spacing @xmath497 and also define @xmath498 .    in the following
we analyze the current @xmath116 , differential conductance @xmath471 , and the fano factor @xmath499 in the parallel and antiparallel magnetic configurations , as well as the corresponding tmr for two - level quantum dots .
the fano factor , @xmath500 , describes the deviation of the zero - frequency shot noise @xmath501 from the poissonian shot noise @xmath502 .
the presented results have been obtained within the real - time diagrammatic technique @xcite .
( color online ) the differential conductance @xmath467 as a function of the bias voltage and level position in the parallel magnetic configuration ( a ) and tmr ( b ) for the parameters : @xmath503 , @xmath504 , @xmath505 , @xmath506 , @xmath507 , and @xmath508 ( @xmath509 , @xmath492 ) .
( after ref .
@xcite ) ]    typical variation of the differential conductance with the bias voltage @xmath60 and level position @xmath276 ( gate voltage ) is shown in fig .
[ sec7fig:2](a ) for the parallel magnetic configuration . conductance in the antiparallel configuration is qualitatively similar to that in the parallel configuration , but generally smaller , which leads to nonzero tmr effect .
the diamonds in fig .
[ sec7fig:2](a ) around @xmath343 correspond to the coulomb blockade regions .
when lowering position of the dot levels , the charge of the dot changes successively .
the dot is empty for @xmath510 , occupied by one electron for @xmath511 , doubly occupied for @xmath512 , occupied by three electrons for @xmath513 , and the two orbital levels of the dot are fully occupied for @xmath514 . in all these transport regions
the dot is in a well - defined charge state , and the sequential tunnelling is exponentially suppressed . if the bias voltage is increased above a certain threshold voltage , the current flows due to first - order ( sequential ) tunnelling processes .
when the thermal energy is low enough , one observes then a well - resolved step in the current as a function of the bias voltage . in the density plots shown in fig .
[ sec7fig:2]a , this can be seen in the form of lines that clearly separate the coulomb blockade regions from transport regions associated with consecutive charge states taking part in transport .
when the bias voltage increases further , additional steps [ and consequently lines in fig .
[ sec7fig:2](a ) ] arise at voltages where new states becomes active in transport .
the corresponding tmr as a function of the bias and gate voltages is shown in fig .
[ sec7fig:2](b ) .
it is worth noting that tmr takes now several well - defined values
. such behavior of tmr is significantly different from that for a single - level quantum dot , where tmr in the sequential tunnelling regime acquires only two values @xcite . as in single - level dots , tmr in the linear response regime
is independent of the gate voltage and is given by @xmath515 , i.e. half of the tmr in the julliere model @xcite .
figure  [ sec7fig:2 ] also shows that , when increasing the bias voltage @xmath60 and keeping constant position of the dot levels , the current and tmr acquire some specific and well - defined values in different transport regions . as shown in ref .
@xcite , the current and tmr at these plateaus can be approximated by simple analytical formulas .
( color online ) current ( a ) in the units of @xmath516 and fano factor ( b ) in the parallel ( solid line ) and antiparallel ( dashed line ) configurations as well as tunnel magnetoresistance ( a ) as a function of the bias voltage for @xmath517 .
the other parameters are the same as in fig .
[ sec7fig:2 ] .
( after ref .
@xcite),height=309 ]    this behavior is presented in more details in fig .
[ sec7fig:4](a ) , where the bias voltage dependence of the current and tmr is shown explicitly for the case when the dot level is above ( @xmath518 ) the fermi level of the leads at equilibrium .
the current in both magnetic configurations and the associated tmr exhibit characteristic plateaus which correspond to different transport regions .
the corresponding fano factors @xmath519 and @xmath520 in the parallel and antiparallel magnetic configurations are shown in fig .
[ sec7fig:4]b .
similarly as current and tmr , the fano factor acquires roughly constant values , different in different transport regions . due to the spin asymmetry in the coupling of the dot to external leads , the bias dependence of the fano factor is significantly different from that in the corresponding nonmagnetic situations @xcite . for both magnetic configurations of the system ,
the fano factor depends on the polarization factor @xmath196 ( differently in the two configurations , in general ) . if @xmath521 , the fano factor becomes divergent due to the thermal noise , which dominates in this transport regime ; in the case of @xmath343 , the noise is given by @xmath522 , with @xmath523 being the linear conductance , leading to a divergency of the fano factor @xcite .
apart from this , in some transport regions we find @xmath524 , while in the other ones @xmath525 .
the ratio @xmath526 depends generally on the spin polarization of the leads @xmath196 @xcite .
furthermore , we note that if the leads are half - metallic , the fano factor in the parallel configuration diverges as @xmath527 .
this increase of the fano factor is due to the enhanced spin asymmetry in transport processes through the dot @xcite . on the other hand , in the antiparallel configuration the fano factor tends to unity for @xmath527 , except for the coulomb blockade regime with two electrons trapped in the dot , where the shot noise is super - poissonian .
when one of the leads is half - metallic ( @xmath528 ) and the other one is nonmagnetic ( @xmath318 ) , transport characteristics become asymmetric with respect to the bias reversal . furthermore , the current can be suppressed in certain bias regions , and this suppression is accompanied by the occurrence of ndc .
this basically happens when the electrons residing in the dot have spin opposite to that of electrons in the half - metallic drain electrode . in fig .
[ sec7fig:5 ] we show the current and fano factor for a quantum dot coupled to half - metallic ( left ) and nonmagnetic ( right ) lead as a function of the bias voltage .
for the situation shown in fig .
[ sec7fig:5 ] , current is suppressed in certain regions of positive bias voltage . as follows from fig .
[ sec7fig:5](a ) , there are three such blockade regions , labelled with the consecutive numbers . on the other hand , for negative bias voltage , the current changes monotonically with the transport voltage , as shown in the inset of fig .
[ sec7fig:5](a ) .
the current ( a ) in units of @xmath516 and fano factor ( b ) as a function of the bias voltage for @xmath529 , @xmath441 , while the other parameters are the same as in fig .
[ sec7fig:2 ] .
the inset in part ( a ) shows the current in the whole range of the bias voltage .
( after ref .
@xcite ) ]    the mechanism of the blockade can be described as follows . in the blockade region
( 1 ) , @xmath530 , the dot is in the state @xmath531 and the spin - down electron residing in the dot has no possibility to tunnel further to the left lead , which leads to suppression of the current . here
, the first ( second ) _ ket _ corresponds to the first ( second ) orbital level of the dot .
the blockade region ( 2 ) , @xmath532 , is associated with the occupation of the states @xmath531 and @xmath533 .
the current is then prohibited due to the full occupation of the single - particle spin - down states .
when increasing the bias voltage further , @xmath534 , the blockade of the current becomes suppressed [ see the plateau between regions ( 2 ) and ( 3 ) in fig .  [ sec7fig:5](a ) ] , which is due to a finite occupation of state @xmath535 .
although tunnelling of spin - down electrons is then blocked , the current is still carried by spin - up electrons . in turn , the blockade region ( 3 ) occurs for @xmath536 , where the dot is in the triplet state @xmath537 and tunnelling is also suppressed .
thus , the current is blocked when the total dot spin @xmath538 is either @xmath539 or @xmath540 , i.e. , the spin of electrons on the dot is opposite to that of electrons in the half - metallic lead .
there is no suppression of current for negative voltage .
this is because the electrons residing in the dot can always tunnel to the nonmagnetic drain electrode .
it is also interesting to note that a pure triplet state is formed in the region ( 3 ) @xcite .    the current blockade in the regions ( 1 ) to ( 3 ) in fig .
[ sec7fig:5](a ) for positive bias is not due to the charging effects as in the coulomb blockade regime , but due to a particular occupation of the dot spin state .
such blockade is frequently referred to as the pauli spin blockade , and has already been found in single - dot and double - dot systems @xcite .
it is worth noting that the current blockade in fig .
[ sec7fig:5](a ) is not complete  there is a small leakage current in each blockade region , which results from the fact that the assumed spin polarization of the half - metallic lead is not exactly equal to unity .
when spin polarization is increased , the current in the blockade regions decreases ( compare the curves for @xmath541 and @xmath542 ) .
the spin blockade of charge current leads to the super - poissonian shot noise , i.e. , the corresponding fano factor is larger than unity . on the other hand ,
the fano factor outside the spin blockade regions is sub - poissonian ( smaller than unity ) , as shown in fig .
[ sec7fig:5](b ) .
the enhancement of the shot noise in the pauli blockade regions is a consequence of large spin asymmetry in the tunnelling processes .
the occurrence of a spin - down electron on the dot prevents further tunnelling processes for a longer time , while spin - up electrons on the dot can escape much faster , allowing further tunnelling processes .
this bunching of fast tunnelling processes gives rise to large current fluctuations , and consequently also to fano factors much larger than unity .
recently an interacting three - terminal quantum dot with ferromagnetic leads was considered in ref .
the dot operated as a beam splitter - one contact was a source and the other two acted as drains .
the authors found a dynamical spin blockade ( spin - dependent bunching of tunnelling events ) and positive zero - frequency cross - correlations of the current in the drain electrodes .      in all transport characteristics discussed above only sequential tunnelling processes
were taken into account , while higher order processes , in particular the cotunneling ones , have been neglected . from our discussion on single - level dots we already know that higher order contributions are particularly pronounced in the blockade regions , where they have a significant influence on tmr and fano factor .
the role of cotunneling processes in two - level dots was studied in a recent paper @xcite .
it has been shown there that the cotunneling processes lead to a significant enhancement ( or reduction ) of tmr in the blockade regime . outside the blockade regions ,
however , tmr is determined mainly by sequential transport , so the contribution from cotunneling processes is rather minor .
similarly , the cotunneling processes significantly modify the fano factor in the blockade regions , while outside these regions the fano factor is only weakly sensitive to the cotunneling .
more specifically , the fano factor in the blockade regions is reduced in comparison to that obtained in the sequential transport limit . however , it is still larger than one , indicating super - poissonian shot noise in the blockade regime due to bunching of the fastest cotunneling processes .
there is also another interesting phenomenon in the blockade regime , which resembles the zero bias anomaly observed in single - level dots in the antiparallel magnetic configuration .
this takes place when each orbital level of the dot at equilibrium is occupied by a single electron , and the system is in the coulomb blockade regime @xcite .
the hamiltonian @xmath543 of the system is given by eq .
( [ sec4eq : qdhamiltonian ] ) , where @xmath544 additionally includes the exchange term and can be expressed by the formula @xcite @xmath545 the last term in @xmath546 corresponds to the exchange energy according to the hund s rule , with @xmath547 being the respective exchange coupling and @xmath548 denoting a vector of pauli spin matrices .
six different two - particle states of the dot are possible , these are : three singlets @xmath549 , @xmath550 , @xmath551 , and three triplets @xmath552 , @xmath553 and @xmath554 . in the case of finite level spacing , @xmath555 , and @xmath556 , the lowest singlet state is @xmath557 .
transport characteristics strongly depend on the ground state of the system .
it is therefore useful to introduce the difference between the energy of the lowest lying singlet @xmath558 and triplet @xmath559 states , @xmath560 .
for @xmath561 , the ground state of the dot is a singlet , @xmath557 , whereas for @xmath562 , the ground state is a triplet , which is three - fold degenerate , @xmath563 , @xmath564 , @xmath565 . on the other hand , for @xmath566 ,
the dot is in a mixed state and the occupation of singlet and each triplet is equal at equilibrium and given by @xmath567 .
( color online ) the differential conductance in the parallel ( solid line ) and antiparallel ( dashed line ) magnetic configuration of the system as a function of the bias voltage for several values of @xmath568 as indicated : ( a ) @xmath569 , ( b ) @xmath570 , ( c ) @xmath571 , ( d ) @xmath572 .
the other parameters are : @xmath573 , @xmath574 , @xmath575 , @xmath576 , and @xmath322 .
( after ref .
@xcite ) ]    figure  [ sec7fig:8 ] presents the bias voltage dependence of differential conductance @xmath471 in the parallel and antiparallel magnetic configuration and for several values of @xmath577 .
let us consider the situation in fig  [ sec7fig:8](c ) , which corresponds to degenerate singlet and triplet states , @xmath571 .
the differential conductance in the parallel configuration exhibits a smooth parabolic dependence on the bias voltage , whereas in the antiparallel configuration there is a maximum at the zero bias voltage .
this effect bears a resemblance to that found in the case of singly occupied one - level quantum dots discussed in previous sections @xcite . here
, however , the mechanism leading to the maximum is different  the zero - bias peak appears due to cotunneling through the singlet and triplet states of the dot . in the case of @xmath571 and at low bias voltages , all the four dot states , i.e. @xmath557 , @xmath563 , @xmath565 , @xmath564 ,
participate in transport on an equal footing .
consequently , the current flows due to both spin - flip and non - spin - flip cotunneling processes .
to understand the mechanism leading to the zero - bias peak , one should bear in mind that in the antiparallel configuration the spin - majority electrons of one lead tunnel to the spin - minority electron band of the other lead .
for example , for positive bias voltage ( electrons tunnel then from the right to left lead ) , the spin-@xmath578 electrons can easily tunnel to the left lead ( the spin-@xmath578 electrons are the majority ones ) , while this is more difficult for the spin-@xmath579 electrons ( they tunnel to the minority electron band ) .
thus , the occupation of state @xmath565 ( @xmath537 ) is increased with increasing bias voltage , while the occupation of state @xmath564 ( @xmath580 ) decreases .
the unequal occupations of these triplet states lead to a nonequilibrium spin accumulation in the dot , @xmath581 .
it is further interesting to note that in the antiparallel configuration the possible non - spin - flip cotunneling processes are proportional to @xmath582 and @xmath583 , whereas the spin - flip cotunneling is proportional to @xmath584 and @xmath585 .
it is clear that the fastest cotunneling processes are the ones involving only the majority spins , i.e. @xmath584 .
however , because of nonequilibrium spin accumulation , with increasing the bias ( @xmath586 ) , the dot becomes dominantly occupied by majority electrons of the right lead , @xmath587 , and the processes proportional to @xmath584 are suppressed . as a consequence ,
the differential conductance drops with the bias voltage , leading to the zero - bias peak .
this is thus the nonequilibrium accumulation of spin @xmath588 which is responsible for the maximum in @xmath471 at low bias voltage .
because the dot is coupled symmetrically to the leads , there is no spin accumulation in the parallel configuration and the differential conductance exhibits a smooth parabolic dependence .
the maximum in the conductance at zero bias disappears for @xmath561 as well as for @xmath589 . instead of the narrow maximum , a relatively broad minimum develops then in the conductance for both magnetic configurations , as shown in fig  [ sec7fig:8](a , b ) or for one magnetic configuration as in fig  [ sec7fig:8](d ) .
consider for instance the case @xmath561 , shown in fig .
[ sec7fig:8](a , b ) . at low bias voltage
the dot is occupied by two electrons on the lowest energy level and the ground state is singlet , @xmath590 .
current is mediated then only by non - spin - flip cotunneling processes .
once @xmath591 , the triplet states start participating in transport leading to an increase in the differential conductance at @xmath592 .
thus , the suppression of cotunneling through @xmath588 states gives rise to a broad minimum in the differential conductance .
the width of this minimum is determined by the splitting between the singlet and triplet states , @xmath593 .
( color online ) tunnel magnetoresistance as a function of the bias voltage for @xmath594 , as indicated in the figure .
the other parameters are the same as in fig .
[ sec7fig:8 ] .
( after ref .
@xcite ) ]    the bias voltage dependence of tmr is shown in fig .
[ sec7fig:10 ] . generally , tmr displays a nontrivial dependence on the ground state of the dot .
when the dot is occupied by a singlet , @xmath561 , tmr displays a maximum plateau at low bias , which can be approximated by , @xmath595 @xcite .
this is the julliere value of tmr . here , it results from the fact that in this transport regime current flows due to non - spin - flip cotunneling . on the other hand , for @xmath596 ,
the tmr exhibits a minimum at zero bias , which can be seen in fig .
[ sec7fig:10 ] for @xmath571 and @xmath597 .
when the number of discrete levels in the dot ( relevant for electronic transport ) becomes larger than two , transport characteristics of the dots attached to ferromagnetic leads become more complex and reveal further interesting features .
such features were observed experimentally mainly in molecule - based quantum dots . when a large natural molecule is weakly attached to metallic leads ,
it can be treated simply as a multi - level qd @xcite . as a specific case
, we discuss transport properties of a single - wall metallic carbon nanotube ( cnt ) weakly coupled to ferromagnetic leads @xcite .
sequential transport in cnts weakly coupled to nonmagnetic and ferromagnetic leads was considered in a recent paper @xcite , where tmr and the fano factor have been calculated within the real time diagrammatic technique .
the system was modelled by the hamiltonian in the form introduced by oreg _
@xcite , @xmath598 ^ 2 \nonumber \\
+ \delta u \sum_{\mu j } n_{\mu j\uparrow } n_{\mu j\downarrow }     + j \sum_{\mu j , \mu^\prime j^\prime }     n_{\mu j\uparrow } n_{\mu^\prime j^\prime\downarrow }     \,,\end{aligned}\ ] ] where @xmath599 , and @xmath600 ( @xmath601 ) is the creation ( annihilation ) operator of an electron with spin @xmath53 on the @xmath229th level in the subband @xmath297 ( @xmath602 ) .
the corresponding energy @xmath603 of the @xmath229th discrete level in the subband @xmath297 is given by @xmath604 , where @xmath275 is the mean level spacing and @xmath605 describes the energy mismatch between the discrete levels corresponding to the two subbands .
the second term in eq .
( [ eq : hnt ] ) stands for the electrostatic energy of a charged cnt , with @xmath271 denoting the charging energy and @xmath606 being the charge on the nanotube induced by gate voltages .
the next term corresponds to the on - level coulomb interaction with @xmath607 being the relevant on - site coulomb parameter . finally , the last term in eq .
( [ eq : hnt ] ) describes the exchange energy , with @xmath547 being the relevant exchange parameter .
( color online ) the differential conductance in the parallel configuration ( a ) and tunnel magnetoresistance ( b ) as a function of bias and gate voltages .
the parameters are : @xmath608 mev , @xmath609 , @xmath610 , @xmath611 , @xmath612 , @xmath613 , @xmath614 , @xmath615 , and @xmath616 mev . ]    in fig .
[ sec7fig:11 ] we show the differential conductance in the parallel configuration and associated tmr , calculated in the first order approximation .
the results are for a specific choice of parameters , namely @xmath617 .
the latter condition indicates that the following sequence of the ground states is realized when filling the cnt with electrons while sweeping the gate voltage ( for @xmath343 ) : @xmath618 , where @xmath501 is the total spin of the nanotube .
this means that in a certain transport regions , the systems is in the triplet state at equilibrium .
[ sec7fig:11](a ) clearly reveals the fourfold pattern of the conductance spectra , associated with filling of the consecutive levels in the two electron subbands of the cnt .
the corresponding tmr is shown in fig .
[ sec7fig:11](b ) , and reveals a complex variation with both gate and bias voltages .    as shown in ref .
@xcite , shot noise calculated in the first order approximation is super - poissonian in the blockade regions , and sub - poissonian in the transport regions were sequential tunnelling dominates .
the corresponding fano factor in the latter case is close to 1/2 , while in the former one is above 1 .
detailed calculations taking into account second order ( cotunneling ) processes show that the fano factor in the blockade regions is generally larger than 1 , although it is reduced in comparison to that derived in the first order approximation , similarly as in the case of two - level dot .
when a quantum dot or a large molecule is strongly coupled to external leads , low - temperature transport characteristics reveal features that are typical of the kondo phenomenon @xcite . to analyze the kondo effect @xcite
, one may use various theoretical techniques .
for instance , one can use the perturbative methods , like the real - time diagrammatic technique discussed in previous sections .
however , one has to go beyond the approximations used in the discussion of the sequential and cotunneling transport .
another approximate method which is frequently used to study electronic transport through quantum dots , including the kondo regime , is the equation of motion method for the relevant green function .
alternatively , one may employ nonperturbative methods , such as for example the numerical renormalization group method @xcite .
the kondo effect @xcite in electronic transport through a quantum dot has been first predicted theoretically in refs .
since then , it has become well documented experimentally @xcite .
the dot with an odd number of electrons has a local spin which at low temperatures , @xmath619 , and in the presence of strong coupling to the electrodes behaves effectively like a magnetic kondo impurity .
exchange coupling with the leads electrons gives rise to spin fluctuations in the dot and screening of the dot s spin .
this , in turn , leads to the formation of a singlet state and kondo resonance peak at the fermi level in the dot s density of states .
this kondo - correlated state leads to an increased transmission through the dot , and also gives rise to a sharp zero - bias anomaly in the current - voltage characteristics .
the successful observation of the kondo effect in semiconductor - based quantum dots , as well as in molecular quantum dots based for instance on carbon nanotubes @xcite or other single molecules @xcite attached to metallic electrodes opened new possibility to study the influence of ferromagnetism on the kondo phenomenon @xcite .    in this section
we briefly discuss how the kondo effect is modified by ferromagnetism of the electrodes .
first , we note that in the extreme case of half - metallic leads , i.e. in the absence of minority - spin electrons , the screening of the dot spin is not possible in the parallel magnetic configuration .
consequently , no kondo - correlated state can be then formed in this particular magnetic configuration .
the possibility of the kondo effect in a quantum dot attached to ferromagnetic electrodes was discussed in a number of papers @xcite .
it was shown , that the kondo resonance in the parallel configuration is split and suppressed in the presence of ferromagnetic leads @xcite .
however , it was also demonstrated that this splitting can be compensated by an appropriate external magnetic field which restores the kondo effect @xcite .
let us begin with the collinear ( parallel and antiparallel ) magnetic configurations .
for the antiparallel configuration and vanishing magnetic field and bias voltage , the model is equivalent to a quantum dot coupled to a single lead with density of states @xmath620 @xcite .
thus , the usual kondo resonance forms in such a case , which is the same as for nonmagnetic electrodes @xcite .
the situation changes in the parallel configuration , where there is an overall asymmetry for up and down spins , say @xmath621 .
this significantly reduces ( or even suppresses ) the kondo effect .
some basic relations concerning the kondo effect , and particularly those concerning the kondo temperature , can be derived from a simple poor man s scaling @xcite performed in two stages @xcite .
first one reduces the energy scale of the effective electron bandwidth @xmath622 starting from @xmath623 .
this leads to a renormalization of the level position @xmath624 according to the scaling equations @xcite @xmath625 this leads to a spin splitting of the level , which in the presence of a magnetic field simply adds to the initial zeeman splitting @xmath626 . as a result
, one finds @xmath627 the scaling of eq .
( [ eq : hal_scaling ] ) is terminated at @xmath628 @xcite .
the strong - coupling limit can be reached by tuning the external magnetic field @xmath496 in such a way that the total effective zeeman splitting vanishes , @xmath629 . in the second stage @xcite ,
spin fluctuations are integrated out .
the schrieffer - wolff transformation @xcite allows to map the anderson model , with renormalized parameters @xmath630 and @xmath631 , to the effective kondo hamiltonian @xcite @xmath632 with @xmath633 in the large-@xmath271 limit . although initially identical , the three coupling constants @xmath634 , @xmath635 , and @xmath636 are renormalized differently during the second stage of scaling .
the relevant scaling equations are @xmath637 all the coupling constants reach the stable strong - coupling fixed point @xmath638 at the kondo energy scale , @xmath639 . for the parallel configuration ,
the kondo temperature in leading order depends on the polarization @xmath196 in the leads following the formula @xmath640 the kondo temperature is maximal for nonmagnetic leads , @xmath641 , and vanishes for @xmath642 .
the unitary limit for the parallel configuration can be reached by tuning the magnetic field , as discussed above . in this case
, the maximum conductance through the quantum dot is the same as for nonmagnetic leads , i.e. @xmath643 per spin .
the scaling procedure was extended in ref .
@xcite to noncollinear configurations , where instead of eq .
( [ eq : kondo_temperature ] ) one finds @xmath644 this formula shows explicitly how the kondo temperature varies with the angle @xmath645 between magnetic moments of the leads .
numerical renormalization - group ( nrg ) technique @xcite is one of the most powerful and accurate methods available currently to study strongly - correlated systems in the kondo regime .
this technique has been adapted recently to the case of a quantum dot coupled to ferromagnetic leads @xcite .
a simple way to model the ferromagnetic leads in the standard nrg procedure is to take the density of states in the leads to be constant and independent of electron spin , @xmath646 , the bandwidths to be equal @xmath647 , and lump all the spin - dependence into the spin - dependent coupling parameter , @xmath648 , which can be taken as independent of energy .
the nrg method , with recent improvements related to high - energy features and finite magnetic field @xcite , is a well - established method to study the kondo impurity ( quantum dot ) physics . using this technique
one can calculate the level occupation @xmath649 , and the spin - dependent single - particle spectral density @xmath650 for arbitrary temperature @xmath31 , magnetic field @xmath496 and polarization @xmath196 . using
the spectral function one can find the spin - resolved linear conductance @xmath651 with @xmath652 denoting the fermi function .
the kondo effect is usually observed as a zero bias anomaly in the differential conductance .
this requires theoretical methods applicable to nonequilibrium situations .
one of such techniques is based on the nonequilibrium green functions .
the latter are usually determined from the relevant equation of motion and the appropriate decoupling scheme . to describe the main features of the kondo effect in transport through quantum dots connected to ferromagnetic leads ( at least qualitatively ) , one may use the decoupling scheme introduced in ref .
@xcite , but generalized by a self - consistent determination of the level energy to account for the exchange field @xcite .
alternatively , one may calculate the exchange field first and then include it in the formalism by hand @xcite .
differential conductance in the p ( a , d ) and ap ( b , e ) configurations , and the corresponding tmr ( c , f ) .
contributions to the conductance from spin - up ( dashed lines ) and spin - down ( dotted lines ) channels are also indicated .
left column ( a , b , c ) corresponds to zero magnetic field , whereas the right one ( d , e , f ) to the case when the compensating magnetic field @xmath653 t is applied ( for the lande factor @xmath654 ) .
the numerical results are for large @xmath271 limit , while the other parameters are : @xmath655 , @xmath656 , @xmath657 .
( after ref .
@xcite ) ]    figure  [ sec8fig:3 ] shows the differential conductance as a function of the transport voltage .
for zero magnetic field there is a pronounced splitting of the peak in the parallel configuration [ fig .
[ sec8fig:3](a ) ] , which can be tuned away by an appropriate external magnetic field [ fig .
[ sec8fig:3](b ) ] . in the antiparallel configuration , in turn , the opposite happens , i.e. there is no splitting at @xmath658 [ fig .
[ sec8fig:3](d ) ] and a finite splitting at @xmath659 [ fig .
[ sec8fig:3](e ) ] with an additional asymmetry in the peak amplitudes as a function of the bias voltage .
figure  [ sec8fig:3](c , f ) shows the corresponding tmr , which becomes negative without magnetic field and changes sign to positive when a compensating magnetic field is applied .    within the eom approach the effect of spin - dependent quantum charge fluctuations
is accounted in the self - consistent but intuitive manner .
the real - time diagrammatic technique @xcite enables one to construct a systematic approach , where the effect of ferromagnetic electrodes can be analyzed without any additional assumptions .
recently the resonant tunnelling approximation ( rta ) was extended by utsumi _
_ @xcite to account for influence of the electrodes ferromagnetism on the kondo phenomenon .
this technique gives more reasonable results and allows for further systematic insight into the physics of the transport through quantum dots in the kondo regime .
the kondo phenomenon in noncollinear magnetic configuration was considered in refs @xcite .
the results show that in symmetric systems , the kondo anomaly gradually disappears when the magnetic configuration varies from antiparallel to parallel one .
however , the main drop of the kondo peak appears already at small deviations from the antiparallel alignment .
recently pasupathy et al .
@xcite studied electronic transport through a single @xmath660 molecule attached to ferromagnetic nickel electrodes in the kondo regime .
it was shown that the kondo correlations appear even in the presence of ferromagnetic leads .
the zero - bias anomaly in the nonequilibrium conductance , however , was split for the parallel alignment of the leads magnetization in agreement with theoretical predictions . for the antiparallel alignment , on the other hand , no splitting of the zero - bias anomaly was observed .
some residual splitting observed in this geometry for some samples can be interpreted as an effect of asymmetric coupling @xmath661 .
similar behavior has also been observed in electronic transport through carbon nanotubes coupled to ferromagnetic leads @xcite .
in this review we addressed certain aspects of electronic transport in mesoscopic tunnel junctions consisting of magnetic metallic nanoparticles , semiconducting quantum dots , or molecules attached to ferromagnetic electrodes .
we are aware that the review does not cover all the problems considered in such systems and many references have not been included .
there are some other aspects of electronic and spin transport which have been omitted .
one of such aspects is spin torque acting on the island and/or electrodes in ferromagnetic single electron transistors based on magnetic metallic nanoparticles @xcite .
another topic omitted in this review is the role of electron coupling to a phonon bath .
such a coupling leads to new interesting phenomena in transport characteristics @xcite .
for instance , the electron - phonon coupling leads to phonon satellite peaks in the differential conductance .
however , the most important spin effects generated by spin - dependent tunnelling through the barriers separating the central part from the ferromagnetic electrodes have been addressed .
we note that spin polarized transport through nanostructures , where the coulomb effects become important , is rapidly developing in recent years .
this applies not only to theoretical part of the subject , but also to the experimental side , where recent progress in nanotechnology allows to attach external leads to individual nanoparticles , quantum dots , and single molecules .
an interesting and new field within this topic is the spintronics based on molecular magnets .
such molecules are good candidates for magnetic memory cells and also are considered as potential candidates for information processing devices .
there is no doubt , that the rapid progress is the physics and technology of spin polarized electronic transport through nanostructures is stimulated by application possibilities in magnetoelectronics , spintronics , and information technology .
we acknowledge discussions with j. knig , j. martinek , w. rudziski , r. wirkowicz , and m. wilczyski .
this work was supported by funds of the ministry of science and higher education as a research project in years 2006 - 2009 .
i. w. also acknowledges support from the foundation for polish science .
100                                      l. p. kouwenhoven , c. m. marcus , p. l. mceuen , s. tarucha , r. m. westervelt , and n. s. wingreen , in proceedings of the nato advanced study institute on mesoscopic electron transport , ed . by l. l. sohn , l. p. kouwenhoven , and g. schn ( kluwer series e345 , 1997 ) .
a. jensen , j. nygrd and j. borggreen in _ proceedings of the international symposium on mesoscopic superconductivity and spintronics _ , edited by h. takayanagi and j. nitta , pp .
33 - 37 , ( world scientific 2003 ) .