electronic transport through dna is crucial for its biological functions such as the repair mechanism after radiation damage and biosynthesis . as proposed in @xcite electron transport through dna proceeds along a one - dimensional pathway constituted by the overlap between @xmath1orbitals in neighboring base pairs .    with the perspective to possible applications in molecular electronics biomaterials
are supposed to play an important role because they may offer the ultimate way to miniaturization and their attributed electrical transport properties are expected to be markedly different from those of traditional macroscopic conductors @xcite . in this context enormous effort has been directed towards investigations of conductivity of dna viewed as a promising conduit for charge transport on nanoscales @xcite .
so far the findings concerning conductivity of dna have been controversial .
there are the results of recent measurements indicating that dna behaves as a well conducting one - dimensional molecular wire @xcite-@xcite .
in contrast it was reported that dna is insulating @xcite and for short oligomers built up from base pairs of the same type semiconductivity was observed @xcite .
this qualitative discrepancy is shared also by theoretical findings so that it remains unclear whether dna is conductive or not .
for example a theoretical investigation based on first principle calculations suggests the absence of dc - conductivity in dna @xcite .
on the other hand there stand the studies using various assumptions for modeling the dna structure which focus on different aspects such as the influence of aperiodicity , temperature driven fluctuations and aggregation effects on the energetic control of charge migration which have all come to the conclusion that dna is a conductor @xcite-@xcite .
the detailed understanding of the conduction mechanism of dna remains still an unsolved task .
different attempts to model the charge transport of dna were based on transport via coherent tunneling @xcite , classical diffusion under the conditions of temperature - driven fluctuations @xcite , incoherent phonon - assisted hopping @xcite,@xcite , variable range hopping between localized states @xcite and charge carriages mediated by polarons @xcite,@xcite and solitons @xcite .    in the present work we use for the theoretical description of charge transport in dna a nonlinear approach facilitating the concept of polaron respectively breather solutions .
the structure of the bent double helix of @xmath0dna is modeled by a sterical network of oscillators in the frame of the base pair picture @xcite,@xcite,@xcite taking into account deformations of the hydrogen bonds within a base pair and twist motions between adjacent base pairs .
the electron motion is described by a tight - binding system .
it is assumed that the electron motion is predominantly influenced by vibrational modes of the double helix .
the nonlinear interaction between the electron and the vibrational modes causes the formation of polarons respectively electron - vibron breathers supporting charge localization and transport .
the aim of the current investigation is to elucidate the role played by static parametrical disorder and complex structural effects due to irregular deviations from the ordered helical shape of the double helix ( reflected in random arrangements of the positions of the bases ) with regard to the charge transport in dna .
there arises at least one interesting question , namely : is long - lived and stable charge transport in dna achievable at all under the impact of disorder ? to tackle this problem in our paper we exploit analytical and numerical methods from the theory of ( nonlinear ) dynamical systems in combination with the theory of ( linear ) disordered systems .
the paper is organized as follows : in the first section we describe the model for the charge transport along the steric structure of the bent dna double helix .
the second section is concerned with the construction of polaron respectively static electron - vibron breather solutions of the nonlinear lattice system when parametrical and structural disorder effects are taken into account . in section [ section : moving ] we discuss the mobility properties of the polarons respectively electron - vibron breathers and their transport efficiency in dependence on the degree of disorder contained in the on - site energies and/or in structural arrangements of the helix . in a subsequent section
we invoke the floquet analysis to explain also the mobility of the polarons and breathers in terms of their spectral properties
. finally , we present a summary and discussion of the results .
dna chains are quite flexible and exhibit structural fluctuations influencing the transport properties .
our model is designed to cover the basic features of the dna double helix structure needed for a proper description of the charge transport dynamics under the impact of fluctuating base motions . despite the fact that the charge transport proceeds across a one - dimensional channel of base stacks spanning along the sugar - phosphate backbone so that dna can effectively be viewed as a one - dimensional molecular wire @xcite some significant details of the three - dimensional dna structure have to be incorporated in the model . for the physical description of dna molecules
their detailed chemical structure can be neglected .
they can be simply considered as bent double - stranded systems with an adequate number of parameters sufficient for modeling the considered physical process . in the frame of simple dna models such as the base pair picture the bases
have been treated as single nondeformable entities .
the helicoidal structure of dna is then conveniently described in a cylindrical reference system where each base pair possesses two degrees of freedom , namely a radial variable measuring the transversal displacements of the base pair ( that is deformations of the h - bond ) and its angle with a reference axis in a plane perpendicular to the helix backbone which defines the twist of the helix @xcite . in fig .
[ fig : sketch ] we depict schematically the structure of the dna model .    ,
@xmath2 , @xmath3 and the radial and angular variables @xmath4 and @xmath5 , respectively are indicated . ]    with relevance to charge migration the most influential fluctuating motions of dna originate from tranverse vibrations of the bases relative to each other within each base plane , i.e. the base pair distance stretchings / compressions and the torsional variations of the helicoidal twist dominating longitudinal ( acoustic strand ) phonons @xcite because the sugar - phosphate chains are rather stretch proof @xcite .
therefore it is justified to discard the impact of longitudinal vibrations and the considered motion is restricted to the base pair planes @xcite .
furthermore , the angular twist and radial vibrational motions evolve independently on two different time scales and can be regarded as decoupled degrees of freedom in the harmonic treatment of the normal mode vibrations of dna @xcite .
the hamiltonian for the electron migration along a strand in dna consists of three parts @xmath6 with @xmath7 representing the part corresponding to the particle charge transport over the base pairs , @xmath8 describes the dynamics of the h - bond vibrations and @xmath9 models the dynamics of the relative twist angle between two consecutive base pairs .
the electronic part is given by a tight - binding system @xmath10 the index @xmath11 denotes the site of the @xmath12th base on one of the two strands and @xmath13 determines the probability that the electron occupies this site .
@xmath14 is the local electronic energy and @xmath15 stands for the transfer matrix element which is responsible for the transport of the electron along the stacked base pairs .
the vibronic part @xmath8 takes into account radial displacements of the base units from their equilibrium positions within the base pair plane connected with deformations of the hydrogen bonds linking two bases .
the corresponding part of the hamiltonian treating the bond vibrations classically and harmonically is given by @xmath16\,.\ ] ] with the momenta @xmath17 conjugate to the radial coordinates @xmath18 measuring the deviation from the equilibrium length of a hydrogen bond along the line bridging two bases of a base pair , @xmath19 is the frequency and @xmath20 is the reduced mass of a base pair .
the twist motion part reads as @xmath21\,,\label{eq : htwist}\ ] ] where @xmath5 is the relative angle between two adjacent base pairs quantifying displacements from its equilibrium twist angle @xmath2 and @xmath22 is the reduced moment of inertia .
we emphasize that in the context of charge transport in dna the nonlinear dynamics of large amplitude structural transitions are not of interest because charge transport processes are thought to be mainly facilitated by soft fluctuational vibrational modes allowing the harmonic approximation of their dynamics .
the interaction between the electronic variable and the structure variables @xmath18 and @xmath5 stems from the parameter dependence of the electronic parameters @xmath14 and @xmath23 . the diagonal term expressing the most efficient coupling is of the form @xmath24 and takes into account the modulation of the on - site electronic energy @xmath25 by the radial vibrations of the base pairs .
in turn the actual charge occupation has its impact on the local radial distortion of the helix . as the transfer matrix elements @xmath23 are concerned we assume that they depend on the three - dimensional distance between two consecutive bases along a strand in the following fashion @xmath26 the quantity @xmath27 regulates how strong @xmath23 is influenced by the distance and the latter is determined by @xmath28 with @xmath29 expanding the expression ( [ eq : dist ] ) up to first order around the equilibrium positions
gives @xmath30\,.\label{eq : distance}\ ] ]    realistic parameters for dna molecules are given by @xcite,@xcite : @xmath31 , @xmath32 , @xmath33 , @xmath34 , @xmath35\,\times10^{12}\,s^{-1}$ ] , @xmath36 , @xmath37 and @xmath38 .
we scale the time according to @xmath39 and introduce the dimensionless quantities : @xmath40 subsequently we omit the tildes .    with the use of the expression ( [ eq : distance ] ) the equations of motion derived from the hamiltonian ( [ eq : hel])-([eq : htwist ] ) read as @xmath41+[c_{n}^*c_{n-1}+c_{n}c_{n-1}^*]\,\right\}\label{eq :
dotr}\\ \ddot{\theta}_{n\,n-1}&=&-\omega^2\,\theta_{n\,n-1}\,-\ , \alpha\,v\,\frac{r_0 ^ 2}{l_0}\ , \sin \theta_0\,[c_{n}^*\,c_{n-1}\,+\,c_{n}\,c_{n-1}^*]\ , , \label{eq : dottheta}\end{aligned}\ ] ] and the quantity @xmath42 appearing in eq .
( [ eq : dotc ] ) determines the time scale separation between the fast electron motion and the slow bond vibrations .
we remark that in the limit case of @xmath43 and constant @xmath44 the set of coupled equations represents the holstein system widely used in studies of polaron dynamics in one - dimensional lattices @xcite . furthermore in the linear limit case emerging for @xmath45 the anderson model
is obtained for random @xmath46 @xcite,@xcite .    regarding the dynamical energy exchange between the various degrees of freedom we state that due to the large differences between the electronic and vibronic frequencies we are able to prove rigorously that the energy exchange is suppressed on time scales growing exponentially with the ratio of the frequencies .
corresponding bounds on the energy exchange suppression are derived with the help of nekhoroshev - like arguments @xcite,@xcite . furthermore applying the concept of the anti - integrable limit
@xcite we are able to prove the existence of breather solutions for the coupled system ( [ eq : dotc])-([eq : dottheta ] ) .
the mathematical details are delegated to a forthcoming article @xcite .
the values of the scaled parameters are given by @xmath47 , @xmath48\,\times 10^{-2}$ ] , @xmath49 , @xmath50 and @xmath51 . in our model study
we treat the electron - mode coupling strengths @xmath52 and @xmath27 for which no reliable data are available as adjustable parameters . in the forthcoming investigations
we adopt the values of these free parameters such that not too strong deformations of the helix result which is essential for our used harmonic treatment of the dynamics of the structural coordinates .
being interested in the simulation of a nonlinear charge transport mechanism in dna we embark as a first step on the construction of localized stationary solutions of the coupled system ( [ eq : dotc])-([eq : dottheta ] ) . exploiting the fact that the adiabaticity parameter @xmath53 is small , i.e. the fast charge transport along the base stacks and the slow bond vibrations and the even slower rotational twist motions evolve on distinct time scales , the inertia in eqs .
( [ eq : dotr ] ) and ( [ eq : dottheta ] ) are negligible .
therefore , one can solve in the adiabatic limit the resulting static equations which results in the instantaneous displacements @xmath54 + [ c_{n}^*c_{n-1}+c_{n}c_{n-1}^*]\,\right\}\ , , \label{eq : rinst}\\ \theta_{n\,n-1}&=&-\frac{\alpha\,v}{\omega^2}\ , \frac{r_0 ^ 2}{l_0}\,\sin \theta_0\ , [ c_n^*c_{n-1}+c_nc_{n-1}^ * ] \label{eq : thetainst}\,.\end{aligned}\ ] ] upon insertion of ( [ eq : rinst ] ) and ( [ eq : thetainst ] ) into eq .
( [ eq : dotc ] ) one obtains a nonlinear discrete schrdinger equation for the electronic amplitude @xmath55 + [ c_n^*c_{n-1}+c_n c_{n-1}^*]\,\right\}\,\right]\,c_n \nonumber\\ & -&(\,1+b_{n+1\,n}\,)\,c_{n+1}\ , -\,(\,1+b_{n\,n-1}\,)\,c_{n-1 } \label{eq : dnls}\,,\end{aligned}\ ] ] with @xmath56\right.\right.\nonumber\\ & + & \left.\alpha\,(1-\cos \theta_0)\ , \left(\,[c_{n+1}^*c_{n}+c_{n+1}c_n^*]\right.\right.\nonumber\\ & + & \left.\left.2\,[c_n^*c_{n-1}+c_n c_{n-1}^ * ] + [ c_{n-1}^*c_{n-2}+c_{n-1 } c_{n-2}^*]\,\right)\right]\nonumber\\ & + & \left .
\frac{\alpha\,v}{\omega^2}\,r_{0}^2\,\sin^2\theta_0\ , [ c_n^*c_{n-1}+c_n c_{n}c_{n-1}^*]\right\}\end{aligned}\ ] ]    in order to search for stationary localized solutions ( polaron states ) we substitute @xmath57 $ ] in eq .
( [ eq : dnls ] ) and obtain the difference system @xmath58 + [ \phi_n^*\phi_{n-1}+\phi_n \phi_{n-1}^*]\,\right\}\,\right ] \,\phi_n \nonumber\\ & -&(\,1+b_{n+1\,n}\,)\,\phi_{n+1}\ , -\,(\,1+b_{n\,n-1}\,)\,\phi_{n-1 } \
, , \label{eq : statthree}\end{aligned}\ ] ] where @xmath59 is given by @xmath60\right.\right.\nonumber\\ & + & \left.\alpha\,(1-\cos \theta_0)\ , \left(\,[\phi_{n+1}^*\phi_{n}+\phi_{n+1}\phi_n^ * ] \right.\right.\nonumber\\ & + & \left.\left.2\,[\phi_n^*\phi_{n-1}+\phi_n \phi_{n-1}^ * ] + [ \phi_{n-1}^*\phi_{n-2}+\phi_{n-1 } \phi_{n-2}^*]\,\right ) \right ] \nonumber\\ & + & \left .
\frac{\alpha\,v}{\omega^2}\,r_{0}^2\,\sin^2\theta_0\ , [ \phi_n^*\phi_{n-1}+\phi_n \phi_{n-1}^*]\right\}\end{aligned}\ ] ]    the ground state of the system ( [ eq : statthree ] ) is computed with the numerical map method outlined in @xcite,@xcite .    .
( b ) the static radial displacements @xmath4 .
( c ) the static angular twists @xmath61 . ]    in fig .
[ fig : profile ] we depict typical patterns of polaron states .
first we discuss the ordered ( periodic ) case arising e.g. for synthetically produced dna molecules consisting of a single type of base pairs ( e.g. poly(g)-poly(c ) dna polymers ) surrounded by vacuum which provides realistic ideal conditions for the consideration of dna conductivity in experimental and theoretical studies @xcite .
the electronic contribution is shown in fig .
[ fig : profile ]  ( a ) . in the ordered case
the electronic wave function is localized at the central lattice site and the amplitudes decay monotonically and exponentially with growing distance from this central site .
we remark that the spatial extension of the polaron decreases with increased coupling strength @xmath52 and there is a smooth transition from large to small polarons .
analogously , enlarging the off - diagonal coupling strength @xmath27 leads also to enhanced degree of localization .
the size of the polarons plays a crucial role for their mobility in the sense that for large up to medium polarons coherent motion can be activated whereas small polarons are immobile due to their pinning to the discrete lattice @xcite,@xcite .
the associated patterns of the static radial and angular displacements of the bases are shown in fig .
[ fig : profile ]  ( b ) and ( c ) , respectively . like the electronic wave function
the radial and angular displacements are exponentially localized at the central lattice site ( base pair ) .
notice that due to the overall minus sign of the radial and angular excitation patterns the h - bridges get compressed while the helix experiences a local unwinding around the occupation peak of the localized electron .     directed along the corresponding pinning mode component .
the kick amplitude is @xmath62 .
( a ) uniformly moving electronic breather .
( b ) the associated moving breather in the radial displacements .
( c ) the angular breather contribution consisting of a small amplitude breather pinned at the starting site and a relatively large amplitude moving breather component propagating in alliance with the electron and radial breathers .
, title="fig : " ]   directed along the corresponding pinning mode component .
the kick amplitude is @xmath62 .
( a ) uniformly moving electronic breather .
( b ) the associated moving breather in the radial displacements .
( c ) the angular breather contribution consisting of a small amplitude breather pinned at the starting site and a relatively large amplitude moving breather component propagating in alliance with the electron and radial breathers .
, title="fig : " ]   directed along the corresponding pinning mode component .
the kick amplitude is @xmath62 .
( a ) uniformly moving electronic breather .
( b ) the associated moving breather in the radial displacements .
( c ) the angular breather contribution consisting of a small amplitude breather pinned at the starting site and a relatively large amplitude moving breather component propagating in alliance with the electron and radial breathers .
, title="fig : " ]    for a more realistic study we take static diagonal disorder in the on - site electronic energy @xmath46 into account .
the random values caused for example by the inhomogeneous broadening of the sites of distinct ion pairs with different energies were simulated by random potentials @xmath63 $ ] with mean value @xmath64 and different mean standard deviations @xmath65 .
( note that any @xmath66 term on the r.h.s .
( [ eq : dotc ] ) with constant @xmath67 can be eliminated by a gauge transformation @xmath68 . )
generally the localized excitation patterns do not change qualitatively when diagonal disorder is taken into account .
compared with the ordered case one merely observes a shift of the excitation peak away from the central lattice site depending on the respective realization of the random potential .
however , since in the presence of diagonal disorder the translational invariance of the lattice system is broken the localized amplitude pattern is no longer reflection symmetric with respect to a distinguished lattice site as opposed to the corresponding symmetric localized state of the ordered system .
moreover we find that the combined effect of the two localization mechanisms , viz .
nonlinear polaron and linear anderson - mode formation respectively , leads to enhancement of the degree of localization compared to the case when only one of the two mechanisms acts .
furthermore , real dna molecules exhibit random structural imperfections of their double helix caused e.g. by the random base sequence , the deforming impact of the chemical surroundings when dna gets buffeted by water molecules and/or the varying hydrophobic potential of the base pair interactions depending on the ambient aqueous solvent that may leave the involved helix structures in irregularly distorted shapes .
accordingly , the structural disorder is incorporated in our model by randomly distributed equilibrium positions of the bases so that the resulting irregular helical dna matrix deviates from the perfectly regular helix structure .    to be precise , we consider randomly distributed torsional coordinates @xmath5 and radii @xmath4 .
since the positions ( actually , the distances ) govern the value of the corresponding transfer matrix elements @xmath69 nondiagonal disorder is induced in the latter at the same time .
random transfer integrals along the helix were used to represent random base pair sequences for which the charge transport along dna was discussed recently in a model relying on localized electron states @xcite .
in this section we study the charge transport supported by mobile polarons respectively breathers propagating along dna . to activate polaron motion an established method is provided by proper stimulation of localized internal normal modes of certain shape and frequency which are also called pinning modes @xcite .
the identification of possibly existing pinning modes is related to the linear stability analysis of the polarons . in order to investigate the stability of the polarons with respect to linear perturbations we derive the linearized system for which we impose small perturbations on the stationary polaron solution @xmath70 , @xmath71 in the form : @xmath72\,\exp(-i\,e t/\tau ) \label{eq : cperturb}\\ r_{n}(t)&=&r_{n}^{(0)}+ \delta r_{n}(t)\\ \theta_{n}(t)&=&\theta_{n}^{(0)}+ \delta \theta_{n}(t)\,,\label{eq : thetaperturb}\end{aligned}\ ] ] where @xmath73 , @xmath74 and @xmath75 are small time - dependent variables .
substitution of ( [ eq : cperturb])-([eq : thetaperturb ] ) into eqs .
( [ eq : dotc])-([eq : dottheta ] ) and retaining only the terms linear in the perturbations results in the system of tangent equations .
the stability analysis is then performed using floquet theory .    to this end
we integrate the linear system of the tangent equations over one period @xmath76 yielding the floquet map @xmath77 with the floquet matrix f determining the evolution of the initial set of variables @xmath78 after times @xmath79 . to ensure linear stability all eigenvalues @xmath80 of the floquet matrix @xmath81 have to lie on the unit circle , i.e. the floquet arguments @xmath82 have to be real - valued .
intensive numerical investigations in a wide parameter range have proved that all the polarons gained from the map method are linearly stable .
moreover the floquet analysis supplies the frequencies of the normal modes via the relation @xmath83 $ ] with @xmath84 expressed in @xmath85 and @xmath86 is an arbitrary integer number .
we identified the pinning modes for given sets of parameters .
according to @xcite we initiate the motion of polarons through suitable perturbations of the velocity variables @xmath87 and/or @xmath88 targeted in the direction of the pinning mode , that is we use for the numerical integration of the system ( [ eq : dotc])-([eq : dottheta ] ) the following initial conditions @xmath89 with the normalized momentum parts @xmath90 of the pinning mode and @xmath91 are the perturbation strengths ( for details see @xcite ) .    in fig .
[ fig : fig3 ] we present results for the activation of breather motion along dna in the ordered case .
we integrated the set of coupled nonlinear equations ( [ eq : dotc])-([eq : dottheta ] ) with a fourth - order runge - kutta method and the accuracy of the computation was checked through monitoring the conservation of the total energy as well as the norm @xmath92 .
the applied kicking strength is @xmath93 corresponding to an amount of radial kinetic energy of the order of @xmath94 .
interestingly , we found that there exists no angular pinning mode component @xmath95 ( see further below ) so that alone the excitation of the radial component @xmath96 produces long - lived stable polaron and breather motions .
hence , our findings suggest that higher frequency radial fluctuational modes of dna instigate the propagation of localized structures along dna rather than low frequency twist fluctuational modes . the fig .  [
fig : fig3 ]  ( a ) demonstrates that the electronic component of the polaron moves with uniform velocity along the lattice maintaining its localized shape throughout the journey . in the same manner the radial polaron part travels as a localized wave along the strand so that the electron propagation is accompanied by a local and temporal contraction of the radii .
however , the angular static polaron component splits up into two parts when the radial pinning mode is initially imposed . at the starting position there
remains a standing breather of comparatively large amplitudes which reverse periodically their sign in the course of time .
thus the helix experiences alternately a local winding and unwinding around its central base pair .
in addition a second breather of relatively small and exclusively negative amplitude propagates in unison with the electronic and radial polaron components .
the result is that the region of the helix which is traversed by the electron gets locally untwisted . as the quantitative assessment of the charge transport
is concerned we note that the breather travels across seventy sites ( base pairs ) in a thousand ( dimensionless ) time units meaning that physically the charge conduction across the base pairs of dna over a distance of @xmath97 takes @xmath98 leading to a higher transport rate than a typical one found for protein conductivity @xcite .
this underlines the status ascribed to dna molecules as promising candidates establishing fast and efficient charge transport in ( albeit only theorized yet ) molecular electronical engineering utilizing them as molecular wires .    in the presence of diagonal static disorder , that is when
the on - site electronic energies @xmath46 are represented by random numbers distributed in the interval @xmath99 with width @xmath65 , we find that mobile breathers exist up to a critical degree of disorder @xmath100 and conductivity persists .
( let us stress that it was indeed experimentally found that the random sequence @xmath0dna is an electrical conductor @xcite . ) beyond this critical value continuous propagation of the breathers along the lattice is prevented .
this behavior is illustrated in fig .
[ fig : fig4 ] displaying the spatio - temporal evolution of the electronic amplitude @xmath101 . while for a relatively small amount of disorder @xmath102
the electronic breather ( together with its radial and torsional counterparts ) is able to travel uniformly along the lattice for a ten times increased @xmath103 the breather motion is restricted to a neighborhood of the initially excited position around which it wanders in an oscillatory manner .     of the random distribution of the on - site electronic energy @xmath104 .
( a ) width @xmath102 .
the electron breather moves with uniform velocity along a strand .
( b ) width @xmath103 .
confined electron breather with positions oscillating itinerantly around the initial lattice position .
, title="fig : " ]   of the random distribution of the on - site electronic energy @xmath104 .
( a ) width @xmath102 .
the electron breather moves with uniform velocity along a strand .
( b ) width @xmath103 .
confined electron breather with positions oscillating itinerantly around the initial lattice position .
, title="fig : " ]    regarding mobility there act two competing mechanisms in the combined disordered nonlinear system , namely a linear anderson - type one related with the effects of disorder which attempts to pin the modes ( bear in mind that the linear anderson - modes are supposed to be immobile ) and opposed to this there operates the nonlinear mechanism utilizing the pinning mode to activate motion .
the relative strength of these two mechanisms @xmath105 determines then breather mobility . in order to summarize the results for the breather mobility in the presence of diagonal disorder we display in [ fig : fig5 ] the temporal behavior of the first momentum of the electronic occupation probability defined as @xmath106    ) , dashed - dotted line ( @xmath107 ) , dotted line ( @xmath108 ) and long dashed line ( @xmath109 ) . ]    which determines the position of the breather center as a function of time .
the breathers are mobile and propagate unidirectionally with ( virtually ) constant velocity along the lattice up to an amount of disorder @xmath110 .
generally , the larger the width @xmath65 the lower becomes this velocity .
finally , for an overcritical amount of disorder the breather performs oscillatory motions around its starting lattice position .
the frequency of these oscillations diminishes gradually with more enlarged @xmath65 and simultaneously the excursions from the starting site get continuously stronger confined . despite the non - directed character of its motion even for a moderate amount of disorder @xmath102 the breather is able to travel as far as twenty lattice sites away from its starting position .
eventually linear anderson - type disorder effects overwhelm more and more the nonlinear activation mechanism for breather motion through the excitation of the nonlinear pinning mode .    concerning the impact of structural disorder we considered random distributions of the base pair spacings @xmath111 $ ] and the twist angles @xmath112 $ ] with different mean standard deviations @xmath113 .    for representative results we averaged over several realizations of structural disorder .
the results for the temporal evolution of the electron breather center are depicted in fig .
[ fig : fig6 ] .
we conclude that structural disorder does not significantly affect the mobility of the breathers . for growing degree of disorder
merely the propagation velocity gets reduced . even in the case of fairly strong disorder @xmath114 ( dashed line )
the breather is still able to propagate steadily .
and angles @xmath61 .
assignment of the line types to the increasing mean standard deviations ranging up to @xmath115 of the respective mean value : full line ( @xmath116 ) , dashed dotted line ( @xmath117 ) and dashed line ( @xmath115 ) . ]
we studied also the polaron mobility under the impact of combined static diagonal and structural disorder and found that the strong pinning effect generated by static diagonal disorder prevails over the more harmless and merely velocity - reducing influence of structural irregularities .
particularly the last result is in compliance with the findings in @xcite where it was reported that hopping conductivity in proteins and dna is influenced first of all by the molecular conformations and not by their sequences .
it is instructive to consider also the linear limit case of the coupled nonlinear electron - vibration system even though the results obtained throughout this paper bear essentially nonlinear character .
the study of the linear limit system itself and its floquet analysis can be applied as a reference and as a tool for obtaining useful information for the interpretation of the floquet analysis of the fully nonlinear system .
the linear limit system is derived from the eqns .
( [ eq : dotc])-([eq : dottheta ] ) for @xmath45 giving @xmath118 in the first equation the terms @xmath119 , which in a homogeneous system reduce to @xmath120 , have been eliminated through the gauge transformation @xmath121 . in the linear limit the electronic and vibrational degrees of freedom are decoupled .    to construct the electronic state we substitute in eq .
( [ eq : dotlinealc ] ) , @xmath122 with time - independent @xmath123 yielding the linear stationary discrete schrdinger equation : @xmath124 . with the bloch mode solutions
@xmath125 we get then the dispersion relation : @xmath126 .
the linear electronic energy spectrum runs from @xmath127 to @xmath128 .
the polaron in the full system displayed in fig .
[ fig : profile ] possesses an energy value @xmath129 which is slightly below the lower band edge @xmath130 assigned to a wave number @xmath131 and in the linear regime the frequencies are in the range @xmath132 $ ] . conclusively , the time scale for the evolution of the angular ( radial ) variables @xmath5 ( @xmath4 ) represented by independent harmonic oscillators with frequencies @xmath133 and @xmath134 , respectively , is two orders ( one order ) of magnitude larger than the one of the electron .
the linear analysis is also useful to understand the origin of many of the numerically calculated floquet arguments for the fully nonlinear system .
the values of the variables @xmath4 and @xmath5 after a period @xmath135 are : @xmath136 and @xmath137 .
therefore , their corresponding floquet arguments are given by @xmath138 and @xmath139 , for @xmath4 and @xmath5 , respectively , which are attributed to angular positions in the complex plane with values around @xmath140 and @xmath141 for polarons with energies slightly below @xmath130 and with increased @xmath142 these angular positions shift to lower values on the unit circle .
the numerical integration of the system of tangent equations over one period of the polaron s oscillation @xmath143 allows us to obtain the actual floquet matrix .
its eigenvalues can be written as @xmath144 with @xmath82 being the floquet arguments .
the system is linearly stable provided all the floquet arguments are real . as noted earlier in this paper within the range of parameters studied in this paper linear stability is assured , but this is not the only information that can be gained from the floquet analysis .    in order to relate the positions of the numerous eigenvalues of the floquet matrix in the complex plane to the various degrees of freedom contained in the system it is elucidating to study first an anti - coupling limit in the electronic subsystem . introducing an auxiliary tunable coupling parameter @xmath145 , which corresponds to the transfer integral in the scaled system , the anti - coupling limit is obtained with @xmath146 and
the actually studied coupled system with @xmath147 .    for the decoupled electron system ( @xmath146 ) the value of the energy @xmath148 can be easily derived with the help of the eqs .
( [ eq : rinst ] ) and ( [ eq : statthree ] ) , giving @xmath149 and @xmath150 , which leads to @xmath151 for the localized state @xmath152 and @xmath153 as a typical value . therefore , the @xmath4-multipliers are in the vicinity of @xmath154 , the @xmath5 s lie at small angles at both sides of @xmath155 and the @xmath13 s are located at @xmath155 on the unit circle , in accordance with the analysis performed in [ sub : linear ] .
as we switch on the coupling @xmath156 and increase it up to the actual value @xmath147 , we can follow the evolution of the multipliers .
eventually , for @xmath147 the @xmath4 s and @xmath5 s have moved to positions with about half their initial angle values while spreading a little on the unit circle ( the @xmath5 s to lesser extent though ) .
the @xmath13 s separate from @xmath155 and invade the whole circle , as a consequence of the nonlinear coupling between the electron amplitudes and the vibrations .
the exceptions are the pair at @xmath155 , corresponding to the phase mode , and another isolated pair close to this point .
the latter pair is associated with the spatially antisymmetric pinning mode , responsible for the mobility of the polaron .
note that the linear approximation described in [ sub : linear ] is an approximation for the radial and angular variables .
the floquet multipliers for the decoupled and full systems , respectively are shown in fig .
[ fig : floquet1 ] .    )
-([eq : dotlinealth ] ) at the anti - coupling limit @xmath146 .
right figure : the fully coupled nonlinear system . ]    in order to illuminate further the role played by the vibrational modes with respect to the initiation of mobility of the polaron we cast the floquet matrix @xmath81 in the form of a supermatrix formed by submatrices : @xmath157 and the indices of the submatrices refer to the origin of their entries . from the arguments of the eigenvalues of the submatrices an / or combinations of them , it is possible to deduce additional information about the spectral features attributed to the various degrees of freedom contained in the coupled nonlinear system .
the matrix @xmath158 appears in the born - oppenheimer approximation . to recover the pinning mode of the complete system it suffices to consider only the subsystem consisting of the electronic part coupled to the radial motions for which the floquet matrix reduces to the submatrix @xmath159 , which means that the mobility of the polaron is mainly a consequence of the coupling between the electron amplitudes and the radial variables and that the electron will not move within the born - oppenheimer approximation .
the figs .
[ fig : floquet2 ] and [ fig : floquet3 ]    illustrate these statements .
the left ( right ) part of the fig .
[ fig : floquet2 ] shows the multipliers of @xmath160 ( @xmath161 ) . in comparison with the linear case studied in the preceding section
we note that while the radial components of the multipliers shifted their positions to lower angles and exhibit some spread the angular multipliers hardly alter their position .
hence , perturbations of the angular twist motions in the context of the coupled system resemble the behavior of the linear limit system . the left part of fig .
[ fig : floquet3 ] depicts the eigenvalues of @xmath162 in the vicinity of @xmath155 .
the pair of pinning mode eigenvalues is recognizable . in the right part of the figure
the normalized pinning mode is portrayed .
( left figure ) and @xmath161 ( right figure ) . ]     in the vicinity of @xmath155 .
the pinning mode eigenvalues are recognizable .
the right part of the figure shows the normalized pinning mode . ]
moreover , it should be pointed out that for larger elongations of @xmath5 , beyond the range where the linear expansion of @xmath163 is valid , there exists also an angular pinning mode as a consequence of the coupling between the electronic amplitudes and the angular variables .
the fact that the disorder eventually prevents the polaron from being moved can also be explained with the help of the floquet multipliers . when the degree of disorder is increased the pinning mode eigenvalues move progressively away from @xmath155 towards larger arguments @xmath84 and approach the rest of the @xmath164 eigenvalues .
[ fig : fdisorder ] displays the growth of the arguments of the pinning mode as a function of the disorder strength @xmath65 .
eventually the system is too far away from the instability due to the marginal mode , which appears when the pair of pinning mode eigenvalues collapses at @xmath155 .
this is easy to understand , as the disorder brings about a decoupling of the variables , favoring localized , stationary , isolated modes .     demonstrating the growth of the argument of the pinning mode with increased amount of disorder so that the pinning mode eventually merges in the rest of the modes . ]
in the present work we have studied an electron transfer mechanism in dna relying on the coupling of the charge carrying unit to the vibrational modes of dna responsible for the formation of polarons and electron - vibron breathers . the steric structure of the double helix of @xmath0dna has been described in the context of the base - pair picture considering angular and radial vibrational motions of the base pairs .
we have focused interest on the initiation of long - range and stable polaron and breather motion along the dna structure when static and/or structural disorder inherent to any real dna molecule is taken into account .
we remark that the proposed polaron - like transport mechanism applies to electron hole conduction too . in the first part of the paper we have constructed exact stationary localized solutions ( polarons and static electron - vibron breathers ) with the help of a nonlinear map approach .
we have found that the interplay of disorder and nonlinearity amplifies the degree of localization in comparison with cases when either disorder or nonlinearity alone is present . with concern to mobile electron - vibron breather solutions
, whose motion is activated through a kick mechanism utilizing the pinning mode , with the current work we have demonstrated that they are fairly robust when the electron - vibration system is subjected to randomness .
generally , even for a moderate amount of disorder the breathers retain their localized shapes and support directed , long - ranged coherent electron transfer along strands of the bent double helix .
in particular the last result suggests that dna seems suitable for the design of ( one - dimensional ) functional nanostructures as ingredients of molecular nanoelectronical devices .
furthermore , we have observed that the electron breather is accompanied by vibrational breathers in the radial as well as angular components , respectively causing local and temporal deformations of the traversed region of the double helix . conclusively ,
efficient charge transport in dna proceeds on the condition that the double helix undergoes structural changes which is exemplary for the interplay of structure and function in flexible and adaptive biomolecules .
one of the authors ( d.h . )
acknowledges support by the deutsche forschungsgemeinschaft via a heisenberg fellowship ( he 3049/1 - 1 ) .
j.f.r.a . acknowledges d.h . and the institut fr theoretische physik for their warm hospitality .
j.a . is grateful for a scholarship supplied by the studienstiftung des deutschen volkes .
the authors are also grateful to the support under the locnet eu network hprn - ct-1999 - 00163 .    99 d.d .
eley and d.i .
spivey , trans .
faraday soc .
* 58 * , 411 ( 1962 ) .
mirkin , r.l .
letsinger r.c .
mucic and j.j .
storhoff , nature * 382 * , 607 ( 1996 ) ; a.a .
voityuk , n. rsch , m. bixon , and j. jortner , j. phys . chem . *
104 * , 9740 ( 2000 ) ; f.c .
grozema , y.a .
berlin and l.d.a .
siebbels , j. am .
soc . * 122 * 10903 ( 2000 ) ; y.a .
berlin , a.l .
burin and m.a .
ratner , j. am .
soc . * 123 * , 260 ( 2001 ) .
m. ratner , nature * 397 * , 480 ( 1999 ) .
. hall , r.e .
holmkin and j.k .
barton , nature * 382 * , 731 ( 1996 ) .
d.b . hall and j.k .
barton , j. am .
* 119 * , 5045 ( 1997 ) .
arkin , e.d.a .
stemp , s.c .
pulver and j.k .
barton , chem .
biol . * 4 * , 369 ( 1997 ) .
y. okahata , t. kobayashi , k. tanaka , and m. shimomura , j. am .
120 * , 6165 ( 1998 ) .
e. meggers , m.e .
michel - beyerle and b. giese , j. am .
soc . * 120 * , 12950 ( 1998 ) .
fink and c. schnenberger , nature * 398 * , 407 ( 1999 ) .
p. tran , b. alavi and g. gruner , phys .
* 85 * , 1564 ( 2000 ) .
e. braun , y. eichen , u. sivan and g. ben - yoseph , nature * 391 * , 775 ( 1998 ) .
d. porath , a. bezryadin , s. de vries and c. dekker , nature * 403 * , 635 ( 2000 ) .
de pablo , f. moreno - herrero , j. colchero , j. g@xmath165mez herrero , p. herrero , a.m. bar@xmath166 , p. ordej@xmath166n , j.m .
soler and e. artacho , phys .
* 85 * , 4992 ( 2000 ) .
j. jortner , m. bixon , t. langenbacher and m.e .
michel - beyerle , proc .
usa * 95 * , 12759 ( 1998 ) .
y .- j . ye , r .- s .
chen , a. martinez , p. otto and j. ladik , sol .
. comm . * 112 * , 139 ( 1999 ) .
r. bruinsma , g. grner , m.r .
dorsogna and j. rudnick , phys .
lett . * 85 * , 4393 ( 2000 ) .
yu and xueyu song , phys .
lett . * 86 * , 6018 ( 2001 ) .
m. hjort and s. stafstrm , phys .
lett . * 87 * , 228101 - 1 ( 2001 ) .
d. ly et al , j. am .
* 118 * , 8747 ( 1996 ) .
j. jortner , proc .
usa , * 95 * , 12759 ( 1998 ) .
e. conwell and s.v .
rakhmanova , proc .
usa , * 97 * , 4556 ( 2000 ) .
d. ly , l. sanii and g.b .
schuster , j. am .
soc . , * 121 * , 9400 ( 1999 ) .
z. hermon , s. caspi and e. ben - jacob , europhys .
lett . * 43 * , 482 ( 1998 ) .
m. peyrard and a.r .
bishop , phys .
lett . * 62 * , 2755 ( 1989 ) .
yakushevich , quart .
* 26 * , 201 ( 1993 ) .
m. barbi , s. cocco and m. peyrard , phys .
lett . a * 253 * , 358 ( 1999 ) .
l. stryer _ biochemistry _ , ( freeman , new york , 1995 ) . s. cocco and r. monasson , j. chem . phys .
* 112 * , 10017 ( 2000 ) .
holstein , ann .
ny * 8 * , 325 , 343 ( 1959 ) . p.w .
anderson , phys .
* 109 * , 1492 ( 1958 ) .
mott , adv . phys .
* 10 * , 107 ( 1961 ) .
d. bambusi and a. giorgilli , j. stat . phys . *
71 * , 569 ( 1993 ) .
d. hennig , j. math . phys .
* 39 * , 3568 ( 1998 ) .
mackay and s. aubry , nonlinearity * 7 * , 1623 ( 1994 ) .
d. hennig , in preparation .
g. kalosakas , s. aubry and g.p .
tsironis , phys .
b * 58 * , 3094 ( 1998 )
. n.k . voulgarakis and g.p .
tsironis , phys . rev .
b * 63 * , 14302 ( 2001 ) .
m. peyrard and m.d .
kruskal , physica d * 14 * , 88 ( 1984 ) .
s. flach and c.r .
willis , phys .
295 * , 181 ( 1998 ) .
d. chen , s. aubry , and g.p .
tsironis , phys .
lett . * 77 * , 4776 ( 1996 ) .
y .- j . ye and j. ladik , int .
. chem . * 52 * , 491 ( 1994 ) .


understanding strong correlation effect away from equilibrium has been one of the most interesting yet challenging problems in condensed matter physics .
a prominent realization of such phenomena is embodied in quantum transport through a nanostructure under finite bias voltage . to understand the interplay of the correlation effect and nonequilibrium nature , nonequilibrium version of the single impurity anderson model ( siam ) and its extensions
have been serving and continue to do so as a central theoretical model .
the siam is indeed considered to be one of the best studied strongly correlated model , and despite its apparent simplicity , it exhibits rich physics already in equilibrium , such as the coulomb blockade and the kondo physics that have been observed in experiments .
equilibrium properties of the siam have been well understood thanks to concerted efforts of several theoretical approaches over years  by perturbative treatment , fermi liquid description , as well as exact results by the bethe ansatz method and numerical renormalization group ( nrg ) calculations  ( see for instance , . ) in contrast , the situation of the nonequilibrium siam is not so satisfactory .
each of the above approaches has met some difficulty in treating nonequilibrium phenomena . a theoretical approach that can deal with the strong correlation effect in nonequilibrium is still called for .    notwithstanding , a number of analytical and numerical methods have been devised to investigate nonequilibrium stationary phenomena : nonequilibrium perturbation approaches  @xcite and its modifications  @xcite , the non - crossing approximation  @xcite , the functional renormalization group treatment  @xcite , quantum monte carlo calculations on the keldysh contour  @xcite , the iterative real - time path integral method  @xcite and so on .
unfortunately those approaches fail to give a consistent picture concerning the finite bias effect on the dot spectral function , particularly regarding a possible splitting of the kondo resonance .    as for the equilibrium siam ,
the second - order perturbation regarding the coulomb interaction @xmath0 on the dot  @xcite is known to capture essential features of kondo physics and agrees qualitatively well with exact results obtained by the bethe ansatz and nrg approaches  @xcite .
such good agreement seems to persist in nonequilibrium stationary state at finite bias voltage .
for the two - terminal particle - hole ( ph ) symmetric siam , a recent study by @xcite showed that the nonequilibrium second order perturbation calculation of the spectral function agrees with that calculated by the diagrammatic quantum monte carlo simulation , excellently up to interaction strength @xmath1 ( where @xmath2 is the total relaxation rate due to leads ) ; pretty well even for @xmath3 at bias voltage @xmath4 .
a typical magnitude of @xmath5 of a semiconductor quantum dot is roughly @xmath6 depending on the size and the configuration of the dot .
therefore , there is a good chance of describing a realistic system within the validity of nonequilibrium perturbation approach .
the great advantage of semiconductor dot systems is to allow us to control several physical parameters .
those include changing gate voltage as well as configuring a more involved structure such as a multi - terminal dot  @xcite or an interferometer embedding a quantum dot .
theoretical treatments often limit themselves to a system with the ph symmetry where the dot occupation number is fixed to be one half per spin .
although assuming the ph symmetry makes sense and comes in handy in extracting the essence of the kondo resonance , we should bear in mind that such symmetry is not intrinsic and can be broken easily in realistic systems , by gate voltage , asymmetry of the coupling with the leads , or asymmetric drops of bias voltage  @xcite .
the ph asymmetry commonly appears in a multi - terminal dot or in an interferometer embedding a quantum dot .
it is also argued that the effect of the ph asymmetry might be responsible for the deviation observed in nonequilibrium transport experiments from the `` universal '' behavior of the ph symmetric siam  @xcite . to work on realistic systems ,
it is imperative to understand how the ph asymmetry affects nonequilibrium transport .    in this paper
, we examine the second - order nonequilibrium perturbation regarding the coulomb interaction @xmath0 of the multi - terminal siam .
the ph symmetry is not assumed , and miscellaneous types of asymmetry of couplings to the leads and/or voltage drops are incorporated as a generic multi - terminal configuration . our main focus is to provide solid analytical results of the behavior of the nonequilibrium self - energy and hence the dot spectral function for the full range of frequency and bias voltage , within the validity of the second - order perturbation theory of @xmath0 .
the result encompasses fermi liquid behavior as well as incoherent non - fermi liquid contribution , showing analytically that increasing finite bias voltage leads to a crossover from the kondo resonance to the coulomb blockade behaviors .
the present work contrasts preceding perturbative studies  @xcite whose evaluations relied on either numerical means or a small - parameter expansion of bias voltage and frequency .
the only notable exception , to the author s knowledge , is a recent work by @xcite , which succeeded in evaluating analytically the second - order self - energy for the two - terminal ph symmetric dot . intending to apply such analysis to a wider range of realistic systems and examine the effect that the two - terminal ph symmetric siam can not capture , we extend their approach to a generic multi - terminal dot where the ph symmetry may not necessarily be present .
an embarrassing drawback of using the nonequilibrium perturbation theory is that when one has it naively apply to the ph asymmetric siam , it may disrespect the preservation of the steady current .
@xcite as a result , one then needs some current - preserving prescription , and different self - consistent schemes have been proposed and adopted .  @xcite
as will be seen , the current preserving condition involves all the frequency range , not only of low frequency region that validates fermi liquid description but also of the incoherent non - fermi liquid part ( see eq .   below ) .
therefore , an approximation based on the low energy physics , particularly the fermi liquid picture , should be used with care .
the self - energy we will construct analytically is checked to satisfy the spectral sum rule at finite bias voltage , so that we regard it as giving a consistent description for the full range of frequency in nonequilibrium . by taking its advantage
, we also demonstrate a self - consistent , current - preserving calculation of the nonequilibrium spectral function for a system where the ph symmetry is not present .
the paper is organized as follows . in sec .
[ sec : model ] , we introduce the multi - terminal siam in nonequilibrium .
we review briefly how to obtain the exact current formula by clarifying the role of the current conservation at finite bias voltage .
section  [ sec : self - energy ] presents analytical expression of the retarded self - energy for a general multi - terminal dot up to the second - order of the interaction strength .
subsequently in sec .
[ sec : discussion ] , we examine and discuss its various analytical behaviors including high- and low - bias voltage limits . sec .
[ sec : spectral - fn ] is devoted to constructing a non - equilibrium spectral function using the self - energy obtained in the previous section .
we focus our attention on two particular situations : ( 1 )  self - consistent , current - preserving evaluation of the nonequilibrium spectral function for the two - terminal ph asymmetric siam , and ( 2 )  multi - terminal effect of the ph symmetric siam .
finally we conclude in sec .
[ sec : conclusion ] .
mathematical details leading to our main analytical result eq .   as well as other necessary material regarding dilogarithm is summarized in appendices .
the model we consider is the single impurity anderson model connecting with multiple leads @xmath7 whose chemical potentials are sustained by @xmath8 .
the total hamiltonian of the system consists of @xmath9 , where @xmath10 , @xmath11 and @xmath12 represent the dot hamiltonian with the coulomb interaction , the hopping term between the dot and the leads , and the hamiltonian of a noninteracting lead @xmath13 , respectively .
they are specified by @xmath14 where @xmath15 is the dot electron number operator with spin @xmath16 and @xmath17 are electron operators at the lead @xmath13 . in the following ,
we consider the spin - independent transport case , but an extension to the spin - dependence situation such as in the presence of magnetic field or ferromagnetic leads is straightforward . when applying the wide - band limit ,
all the effect of the lead @xmath13 is encoded in terms of its chemical potential @xmath8 and relaxation rate @xmath18 , where @xmath19 is the density of states of the lead @xmath13 . the dot level @xmath20 controls the average occupation number on the dot ; it corresponds roughly to @xmath21 , @xmath22 , @xmath23 for @xmath24 , @xmath25 , and @xmath26 , respectively .
the ph symmetry is realized when @xmath27 and @xmath28 ( see eqs .   and below ) .
we here briefly summarize how the current through the dot is determined in a multi - terminal setting .
special attention is paid to the role of the current conservation because it has been known that nonequilibrium perturbation calculation does not respect it in general .
@xcite we illustrate how to ensure the current conservation by requiring minimally .
the argument below is valid regardless of a specific approximation scheme chosen , whether nonequilibrium perturbation or any other approach .
following the standard protocol of the keldysh formulation  @xcite , we start with writing the current @xmath29 flowing from the lead @xmath13 to the dot in terms of the dot s lesser green function @xmath30 and the retarded one @xmath31 : @xmath32,\end{aligned}\ ] ] where @xmath33 is the fermi distribution function at the lead @xmath13 .
as the present model preserves the total spin as well as charge , the net spin current flowing to the dot should vanish in the steady state , which imposes the integral relation between @xmath30 and @xmath31 , @xmath34 = 0 .
\label{eq : gless - gr}\end{aligned}\ ] ] here we have introduced the total relaxation rates @xmath35 and the effective fermi distribution @xmath36 weighted by the leads , @xmath37 when we ignore the energy dependence of the relaxation rates @xmath38 , we can recast eq .   into a more familiar form @xmath39 because @xmath40 is the definition of the exact dot occupation number . note the quantity @xmath41 is nothing but the exact dot spectral function out of equilibrium . we emphasize that eq .
( or equivalently eq .  ) is the minimum , exact requirement that ensures the current preservation .
it constrains the exact @xmath42 and @xmath43 that depend on the interaction as well as bias voltage in a nontrivial way .
one can accordingly eliminate @xmath44 in @xmath29 , to reach the landauer - buttiker type current formula at the lead @xmath13 , @xmath45 or , the current conservation allows us to write it as @xmath46,\end{aligned}\ ] ] where @xmath47 is the exact number of states with spin @xmath16 at finite temperature in general , defined by @xmath48 it tells us that differential conductance @xmath49 with fixing all other @xmath50 s is proportional to the nonequilibrium dot spectral function , provided changing @xmath8 does not affect the occupation number  @xcite .
such situation is realized , for instance , when a probe lead couples weakly to the dot .
the case of a noninteracting dot always satisfies the current preserving condition eq .   as @xmath51 holds for any @xmath52 ;
the distribution function of dot electrons @xmath53 is equal to @xmath54 .
this is not the case for an interacting dot , however .
as for the interacting case , not so much can be said .
we only see the special case with the two - terminal ph symmetric dot satisfy eq .   by choosing @xmath55 irrespective of interaction strength . except for this ph symmetric case ,
a general connection between @xmath42 and @xmath43 is not known so far .
it is remarked that , based on the quasiparticle picture , a noninteracting relation @xmath56 is sometimes used to deduce an approximate form of @xmath42 out of @xmath43 for an interacting dot .
such approximation is called the ng s ansatz  @xcite .
although it might be be simple and handy , its validity is far from clear .
we will not rely on such additional approximation below .
it is also important to distinguish in eq .
the electron occupation number @xmath57 from the quasiparticle occupation number @xmath58 , as the two quantities are different at finite bias voltage since the luttinger relation holds only in equilibrium  @xcite .
contribution to the dot occupation number comes from all range of frequency , including the incoherent part .
one sees fulfilling the spectral weight sum rule @xmath59 crucial to have the dot occupation number @xmath57 normalized correctly . in general
, one needs to determine @xmath57 appropriately to satisfy eq .   as a function of interaction and chemical potentials of the leads .
the applicability of quasi - particle approaches that ignores the incoherent part is unclear .
in this section , we evaluate analytically the nonequilibrium retarded self - energy up to the second order of interaction strength @xmath0 for the multi - terminal siam .
we first examine the contribution at the first - order and the role of current preservation
. then we present the analytical result of the second - order self - energy in terms of dilogarithm .    following the standard treatment of the keldysh formulation ,  @xcite the nonequilibrium green function and the self - energy
take a matrix structure @xmath60 satisfying symmetry relations @xmath61 and @xmath62 . the retarded green function is defined by @xmath63 ; the retarded self - energy , by @xmath64 .    to proceed evaluation , it is convenient to classify self - energy diagrams into the two types : the hartree - type diagram ( fig .
[ fig : hartree - type ] ) that can be disconnected by cutting a single interaction line , and the rest which we call the correlation part and reassign the symbol @xmath65 to .
the latter starts at the second - order .
the resulting green function ( matrix ) takes a form of @xmath66^{-1},\end{aligned}\ ] ] where @xmath67 represents a pauli matrix of the keldysh structure , and @xmath68 refers to the exact occupation number of the dot electron with the opposite spin .
accordingly , the corresponding retarded green function becomes @xmath69 where @xmath70 is the hartree level of the dot .    .
the double line refers to the exact green function . ]      before starting evaluating the correlation part @xmath71 that starts contributing at the second order , it is worthwhile to examine the current preserving condition eq .   at the first order . at this order
, it reduces to the self - consistent hartree - fock equation for the dot occupation number @xmath72 : @xmath73 .
\label{eq : hf - self - consistency}\end{aligned}\ ] ] it shows how the two - terminal ph symmetric siam is special by choosing @xmath74 , @xmath75 and @xmath76 ; the second term of the right - hand side vanishes by having the solution @xmath77 even at finite bias voltage .
it also indicates that the current preservation necessarily has the occupation number depend on asymmetry of the lead couplings as well as interaction strength for the ph asymmetric siam .
indeed , for a small deviation from the ph symmetry and bias voltage , we see the hartree - fock occupation number behave as @xmath78 where @xmath79 is the average chemical potential weighted by leads , @xmath80 note @xmath79 vanishes when no bias voltage applies , as we incorporate the overall net offset by leads into @xmath20 .              following the standard perturbation treatment of the keldysh formulation , we see there is only one diagram contributing to @xmath81 at the second order ( fig .  [
fig : sigma2 ] ) after eliminating the hartree - type contribution .
the contribution is written as @xmath82 where @xmath83 refer to to the unperturbed green functions ( including the hartree term ) , whose concrete expressions are found in appendix a. the polarization matrix @xmath84 is defined by @xmath85 ( fig .
[ fig : pi ] ) .. ]    as was shown by the current formula in the previous section , we need only the dot spectral function to study quantum transport , hence @xmath71 suffices .
therefore it is more advantageous to work on the rak - representation , where the polarization parts become    @xmath86 , \\      & \pi^{a}_{12 } = \frac{i\hbar}{2 } \left [ g^{a}_{12 } g^{k}_{21 } +        g^{k}_{12 } g^{r}_{21 }      \right ] , \\      & \pi^{k}_{12 } = \frac{i\hbar}{2 } \left [ g^{k}_{12 } g^{k}_{21 } +        g^{r}_{12 } g^{a}_{21 } + g^{a}_{12 } g^{r}_{21 } \right ] ,    \end{aligned}\ ] ]    [ eq : def - pi ]    and their fourier transformations are given in appendix b. accordingly , we can express the retarded self - energy @xmath71 as @xmath87,\end{aligned}\ ] ] where @xmath88 the above second - order expression of @xmath71 is standard , but it has so far been mainly used for numerical evaluation , quite often restricted for the two - terminal ph symmetric siam .
we intend to evaluate eqs .
( [ eq : i1],[eq : i2 ] ) analytically for the generic multi - terminal siam , along the line employed in ref .  .
delegating all the mathematical details to appendices c and d , we summarize our result of the analytical evaluation of @xmath71 as follows :    @xmath89 , \label{eq : sigmar}\end{aligned}\ ] ]    here , functions @xmath90 ( @xmath91 ) are found to be ( using @xmath92 ) ,    [ eq : xi123 ] @xmath93 \nonumber \\
& \qquad \qquad   + \sum_{a , b,\beta }        \tfrac{4\gamma_{a }          \gamma_{b}}{\gamma^{2 } } \big [ { \mathop{\mathrm{li}_{2}}\nolimits}(\tfrac{-\varepsilon+ \beta \zeta_{a\bar{\sigma } } } { \beta          \zeta_{b \bar{\sigma } } + i\gamma } ) + \tfrac{1}{4 } { \mathop{\mathrm{log}}\nolimits}^{2 }        ( \tfrac{-\zeta_{a\bar{\sigma}}+i\gamma}{\zeta_{b\bar{\sigma } } +          i\gamma } )        \big ] , \\
\xi_{2}(\varepsilon ) & =   6\pi^{2 }        - \sum_{a , b,\beta } \tfrac{4\gamma_{a } \gamma_{b}}{\gamma^{2 } } \big [        \lambda(\tfrac{\varepsilon -\zeta_{a\sigma } +          2i\gamma}{\beta \zeta_{b\bar{\sigma } } + i\gamma } ) +        \lambda(\tfrac{\varepsilon - \beta \zeta_{b\bar{\sigma } }          + 2i\gamma}{\zeta_{a\sigma } + i\gamma } ) + \tfrac{1}{2 } { \mathop{\mathrm{log}}\nolimits}^{2 }        ( \tfrac{\zeta_{a\sigma } + i\gamma}{\beta \zeta_{b \bar{\sigma } } +          i\gamma } ) \big ] \nonumber \\ & \qquad   \qquad - \sum_{a , b,\beta }        \tfrac{4\gamma_{a } \gamma_{b}}{\gamma^{2 } } \big [ \lambda        ( \tfrac{\varepsilon+ \beta \zeta_{a\bar{\sigma } } +          2i\gamma } { \beta \zeta_{b\bar{\sigma } } + i\gamma } ) +        \tfrac{1}{4 } { \mathop{\mathrm{log}}\nolimits}^{2 } ( \tfrac{-\zeta_{a\bar{\sigma } } +          i\gamma}{\zeta_{b\bar{\sigma } } + i\gamma } )        \big ] , \\
\xi_{3 } & = \left[\sum_{a }        \tfrac{2\gamma_{a}}{\gamma }        { \mathop{\mathrm{log}}\nolimits}\left ( \tfrac{-\zeta_{a\bar{\sigma } } +            i\gamma}{\zeta_{a\bar{\sigma } } + i\gamma}\right ) \right]^{2 } ,      \end{aligned}\ ] ] [ eq : xi-123 ]    where the summations over @xmath94 as well as terminals @xmath95 are understood .
function @xmath96 is dilogarithm , whose definition as well as various useful properties are summarized in appendix c ; @xmath97 is defined by   is equivalent to that given in ref .  , but we prefer writing it in this form because its analyticity is more transparent . ]
@xmath98 { \mathop{\mathrm{log}}\nolimits}z .
\label{def : lambda}\end{aligned}\ ] ] analytical formula @xmath71 given in eqs .   and
constitutes the main result of the present paper .
consequently , the nonequilibrium spectral function of the multi - terminal siam is given analytically for full range of frequency and bias voltage , once one chooses @xmath57 to satisfy eq .  .
the result also applies to a more involved structured system , such as an interferometer embedding a quantum dot , by simply replacing @xmath20 and @xmath38 to take account of those geometric effect .
having obtained an explicit analytical form of the second - order self - energy @xmath99 at arbitrary frequency and bias voltage , we now examine its various limiting behaviors .
most of those limiting behaviors have been known for the two - terminal ph symmetric siam , so it s assuring to reproduce those expressions in such a case . simultaneously , our results below provide multi - terminal , ph asymmetric extensions of those asymptotic results .
we can reproduce the equilibrium result by setting all the chemical potentials equal , @xmath100 . then we immediately see @xmath101 and @xmath102 , \\ & \xi_{2 } = 8\cdot \big[\frac{3\pi^{2}}{4 } - 3\lambda(\tfrac{\varepsilon_{\sigma } + 2i\gamma}{i\gamma } )    \big],\end{aligned}\ ] ] where @xmath103 . as a result ,
the correlation part of the self - energy in equilibrium becomes @xmath104 .
\label{eq : sigma - equilibrium}\end{aligned}\ ] ] the ph symmetric case in particular corresponds to @xmath105 .
it reproduces the perturbation results by yamada and yoshida  @xcite up to the second - order of @xmath0 , when we expand the above for small @xmath52 .
the ph symmetric result is indeed identical with the one obtained in ref .   for arbitrary frequency ( see also eq .   below ) .
@xcite have evaluated analytically the self - energy and the spectral function for the two - terminal ph symmetric siam . in our notation
, it corresponds to the case @xmath106 , and @xmath107 .
when we parametrize the two chemical potentials by @xmath108 with @xmath109 in eqs .  , the self - energy can be written as @xmath110 , \label{eq : ph - symmetric - sigmar}\end{aligned}\ ] ] where @xmath111 , \\ & \xi_{2 } =   6\pi^{2 }   - 6\sum_{a , b } \big [      \lambda(\tfrac{\omega -a ev/2 + 2i\gamma}{b ev/2 +   i\gamma } )   + \frac{1}{4 } { \mathop{\mathrm{log}}\nolimits}^{2 } ( \tfrac{a ev/2 + i\gamma}{b ev/2 + i\gamma } ) \big].\end{aligned}\ ] ] the above result are identical with what ref .   obtained .
we now employ the small - parameter expansion of @xmath71 around the half - filled equilibrium system . here
we assume parameters @xmath112 and @xmath113 are much smaller than the total relaxation rate @xmath2 .
the expansion of @xmath114 is found to contain the first- and second - order terms regarding @xmath115 and @xmath116 , while @xmath117 do only the second - order terms .
therefore the result of the expansion up to the second - order of these small parameters is presented as @xmath118.\end{aligned}\ ] ] functions @xmath90 can be expanded straightforwardly by using the taylor expansion of dilogarithm in appendix c. they are found to behave @xmath119 , \\ & \xi_{2}(\varepsilon ) \approx   - 8 \ ,   \big [   \frac{3 ( -\varepsilon^{2 } + 2    \varepsilon\bar{\mu } + 2    \bar{\mu}^{2 } -     \overline{\mu^{2 } } - 2    \overline{\mu^{2}})}{4(i\gamma)^{2 } }    \big ] \\ & \xi_{3 } \approx \frac{16 \bar{\mu}^{2}}{(i\gamma)^{2}},\end{aligned}\ ] ] where @xmath79 is defined in eq .
and we have introduced @xmath120 combining all of these , we reach the small bias / frequency behavior of the self - energy @xmath81 as @xmath121 \nonumber \\ & \qquad \qquad -   \frac{iu^{2}}{2\pi^{2}\gamma^{3 } } \big [   ( \omega - \bar{\mu})^{2 } + 3\ , ( \delta \mu)^{2 } \big ] .
\label{eq : small - parameter - expansion}\end{aligned}\ ] ]    small bias expansion of @xmath122 for the two - terminal system has been discussed and determined by the argument based on the ward identity .
@xcite the dependence appearing in eq
.   fully conforms to it ( except for the presence of the bare interaction instead of the renormalized one ) .
indeed , correspondence is made clear by noting the parameters @xmath79 and @xmath123 for the two - terminal case become @xmath124 the presence of linear term in @xmath52 and @xmath125 for the two - terminal ph asymmetric siam was also emphasized recently .
@xcite      one expects naively that the limit of large bias voltage @xmath126 corresponds to the high temperature limit @xmath127 in equilibrium ; it was shown to be so for the two - terminal ph symmetric siam .
@xcite we now show that the same applies to the multi - terminal siam where bias voltages of the leads are pairwisely large , _
i.e. _ , half of them are positively large , and the others are negatively large .    in the large bias voltage limit , all the arguments of dilogarithm functions appearing in eqs .
become @xmath128 , where the values of dilogarithm are known ( see appendix c ) .
accordingly , the pairwisely large bias limit of @xmath90 is found to be @xmath129 correspondingly , the retarded self - energy becomes @xmath130 it shows that the result of the multi - terminal siam is the same with that of the two - terminal ph symmetric siam except for a frequency shift .
accordingly , the retarded green function @xmath131 in this limit is given by @xmath132 the form indicates that for sufficiently strong interaction @xmath133 , the dot spectral function has two peaks at @xmath134 with broadening @xmath135 , so the system is driven into the the coulomb blockade regime . on the other hand , for weak interaction @xmath136
, it has only one peak with two different values of broadening that reduce to @xmath2 and @xmath137 in the @xmath138 limit .    what is the role of the current preservation condition eq .   in this limit ?
it just determines the dot occupation number explicitly .
in fact , the condition becomes @xmath139 with @xmath140 , and @xmath57 is independent of the interaction strength because bias voltage sets the largest scale .
one can evaluate the above integral exactly to have @xmath141 where @xmath142 . as a result , expanding it up to the second order of @xmath143 leads to @xmath144 the first term corresponds to the occupation number that one expects naturally from the effective distribution @xmath36 ; it corresponds , for instance , @xmath145 for the two - terminal dot with @xmath146 with @xmath147 .
the second term is the deviation from it , which is proportional to the average of the inverse chemical potential weighted by the leads .
we now turn our attention to the behavior of the nonequilibrium dot spectral function , using our analytical expression of the self - energy eqs .
( [ eq : sigmar],[eq : xi123 ] ) .
below we particularly focus our attention on the two cases : the two - terminal ph asymmetric siam where current preservation has been an issue , and the multi - terminal ph symmetric siam where the role of multiple leads has been raising questions . in all of the calculations below ,
we have checked numerically the validity of the spectral weight sum rule at each configuration of bias voltages .      as was emphasized in sec .
[ sec : current - formula ] , when a dot system does not retain the ph symmetry , the stationary current is not automatically conserved and one must impose the current preservation condition eq .   explicitly . as the right - hand side of eq .   also depends on the dot occupation number @xmath57 , this requires us to determine @xmath57 self - consistently by using the retarded green function in a certain approximation  the second - order perturbation theory in the present case .    ) at finite bias voltage @xmath148 .
the interaction strength is chosen as @xmath149 .
the dotted line represents the result of @xmath150 and @xmath151 . ]
figure  [ fig : ph - symmetric ] shows the result of nonequilibrium dot spectral function of the two - terminal ph symmetric siam at bias voltage @xmath152 , and @xmath153 , which is essentially the same result with ref .
( of a different set of parameters ) .
the occupation number is fixed to be @xmath154 in this case , so its self - consistent determination is unnecessary .
the results were compared favorably with those obtained by diagrammatic quantum monte carlo calculations  @xcite ; a relatively good quantitative agreement was observed up to @xmath155 ( where the bethe ansatz kondo temperature @xmath156  @xcite while the estimated half width of the kondo resonance @xmath157 ) and bias voltage @xmath158 . applying bias voltage
gradually suppresses the kondo resonance without splitting it , and the two peaks at @xmath159 are developed at larger bias voltages , which corresponds to the discussion in the previous section .
figure  [ fig : ph - asymmetric ] shows the result of our self - consistent calculation of the nonequilibrium spectral function for the two - terminal ph asymmetric siam at ( a ) @xmath160 and ( b ) @xmath161 .
a paramagnetic - type solution is assumed in determining @xmath57 . as in the ph symmetric siam
, one sees increasing bias voltage not split but suppress the kondo resonance while it develops a peak around @xmath162 .
the kondo resonance peak is suppressed more significantly at @xmath160 than at @xmath163 , because the kondo temperature of the former ( @xmath164 ; @xmath165 ) is smaller than that of the latter ( @xmath166 ; @xmath167 ) .
an interesting feature of the ph asymmetric siam is that spectral weight of the kondo resonance seems shifting gradually toward @xmath168 with increasing bias voltage , without exhibiting a three - peak structure in the ph symmetric case .
this suggests a strong spectral mixing between the kondo resonance and a coulomb peak at finite bias voltage .
because of it , the interval of the two peaks at finite bias is observed as roughly @xmath169 and gets widened up to @xmath0 for larger @xmath170 .
the bias dependence somehow looks similar to what was obtained by assuming equilibrium noninteracting effective distribution for @xmath57  @xcite ( which is hard to justify in our opinion ) , though we emphasize our present calculation only relies on the current - preservation condition without using any further assumption .
it is remarked that the effect shown by bias voltage is quite reminiscent of finite temperature effect that was observed in the ph asymmetric siam in equilibrium  @xcite .    , and ( b ) @xmath171 .
all the other parameters are the same with fig .
[ fig : ph - symmetric ] . as an eye guide
, the ph symmetric result of @xmath150 and @xmath151 is shown as a dotted line.,title="fig : " ] , and ( b ) @xmath171 .
all the other parameters are the same with fig .
[ fig : ph - symmetric ] . as an eye guide
, the ph symmetric result of @xmath150 and @xmath151 is shown as a dotted line.,title="fig : " ]    .
results of the interaction strengths @xmath172 , @xmath173 , and @xmath174 are shown , while dotted lines refer to the noninteracting case as an eye guide .
( a ) spectral function of the ph symmetric siam at @xmath175 .
( b ) spectral function of the ph asymmetric siam at @xmath176.,title="fig : " ] .
results of the interaction strengths @xmath172 , @xmath173 , and @xmath174 are shown , while dotted lines refer to the noninteracting case as an eye guide .
( a ) spectral function of the ph symmetric siam at @xmath175 .
( b ) spectral function of the ph asymmetric siam at @xmath176.,title="fig : " ]    more insight can be gained by examining how the spectral structure depends on the interaction strength at finite bias voltage .
figures  [ fig : ldos - u ] ( a , b ) show a structural crossover from a noninteracting resonant peak ( the dotted line ) to correlation peaks , for ( a ) the ph symmetric siam @xmath175 , and ( b ) the ph asymmetric siam @xmath177 .
the ph symmetric siam shows introducing @xmath0 leads to developing the correlation two peaks as well as the kondo peak that is suppressed by finite bias voltage .
in contrast , the bias voltage effect on the ph asymmetric siam is more involved , because the kondo resonance is apparently shifted and mixed with one of the correlation peaks , eventually showing the two - peak structure at @xmath178 in large bias voltage limit .      , and @xmath179 ) .
other parameters are chosen as the same as in fig .
[ fig : ph - symmetric ] .
the dotted line correspond to the two - terminal noninteracting unbiased case , while the dashed line , to the two - terminal interacting unbiased case . ]    to examine finite bias affects further and see particularly how the presence of multi - terminals affects the nonequilibrium spectral function , we configure a special setup of the multi - terminal siam that preserves the ph symmetry : the dot is connected with @xmath180 identical terminals , with bias levels being distributed equidistantly between @xmath181 and @xmath182 , and each of relaxation rates is set to be @xmath183 .
the latter ensures that the unbiased spectral function is the same , hence the kondo temperature .
results of the nonequilibrium spectral function are shown in fig .
[ fig : nldos - by - v ] .
again , we confirm that no splitting of the kondo resonance is observed in this multi - terminal setting .
one sees further that the increasing the number of terminals _ enhances _ the kondo resonance .
this can be understood by weakening the bias suppression effect on the kondo resonance for a larger @xmath180 .
more precisely , one may estimate the suppressing effect from small bias behavior eq .  .
hence @xmath184 is a control parameter . in the present multi - terminal ph symmetric
setting , the quantity @xmath185 is found to be @xmath186 therefore @xmath185 decreases with increasing @xmath180 , which results in weakening the suppression and enhancing the kondo resonance for a larger @xmath180 .
( @xmath187 ) .
other parameters are the same as in fig .
[ fig : nldos - by - v ] .
( b )  the effective fermi distribution @xmath188 at zero temperature for the ph symmetric multi - terminal dot ( @xmath189 ) .
the inset shows relative locations of bias levels with fixed @xmath190 as a function of the number of terminals.,title="fig : " ]   ( @xmath187 ) .
other parameters are the same as in fig .
[ fig : nldos - by - v ] .
( b )  the effective fermi distribution @xmath188 at zero temperature for the ph symmetric multi - terminal dot ( @xmath189 ) .
the inset shows relative locations of bias levels with fixed @xmath190 as a function of the number of terminals.,title="fig : " ]    the preceding argument also tells us that if the spectral function bears any multi - terminal signatures at all , they would be more conspicuous by examining it with fixing @xmath185 rather than @xmath125 .
this is done in fig .
[ fig : nldos - by - dmu ] ( a ) ; figure  [ fig : nldos - by - dmu ] ( b ) shows how the effective dot distribution @xmath188 and the relative locations of bias levels ( in the inset ) evolve for a fixed @xmath185 when @xmath180 increases .
no multi - terminal signature in the nonequilibrium spectral function is seen in fig .
[ fig : nldos - by - dmu ] ( a ) ; results of different @xmath180 actually collapse , not only around zero frequency but in the entire frequency range .
it suggests that the suppression of the kondo resonance deeply correlates the development of coulomb peaks , and a mixing between those spectral weight is important .
the parameter @xmath185 controls a crossover from the kondo resonance to the coulomb blockade structure .
we may also understand the similarity between bias effect and temperature effect by the connection through the large @xmath180 limit of the effective fermi distribution @xmath188 , as shown in fig .
[ fig : nldos - by - dmu ] ( b ) .
although there is a consensus that bias voltage starts suppressing the kondo peak , and eventually destroys it with developing the two coulomb peaks when bias voltage is much larger than the kondo temperature , there is a controversy as to whether the kondo resonance peak will be split or not in the intermediate range of bias voltage .
all the results obtained by the second - order perturbation consistently indicate that there is no split of the kondo resonance ; finite bias voltage starts to suppress the kondo resonance , and develops the two coulomb blockade peaks by shifting the spectral weight from the kondo resonance .
we should mention that some other approximations draw a different conclusion . here
we make a few remarks on apparent discrepancy seen in various theoretical approaches as well as experiments , as well as some speculation based on our results .
typically , several approaches that rely on the infinite @xmath0 limit , notably non - crossing approximation , equation of motion method , and other approaches investigating the kondo hamiltonian , observed the splitting of the kondo resonance under finite bias voltage  @xcite . those results , however , have to be interpreted with great care , in our view .
generally speaking , the spectral function obtained by those approaches does not obey the spectral weight sum rule : ignoring the doubly occupied state typically leads to the sum rule @xmath191  @xcite , rather than the correct value .
therefore only half of the spectral weight can be accounted for in those methods .
simultaneously , such ( false ) sum rule in conjugation with the bias - suppression of the kondo resonance can not hep but lead to a two - peak structure of the spectrum within the range of attention .
splitting of the kondo resonance might be an artifact of the approximation .
not fulfilling the correct sum rule , those approaches may not be able to distinguish whether finite bias will split the kondo resonance or simply suppresses it with developing the coulomb peaks .
as for the two - terminal ph symmetric siam , fourth - order contribution regarding the coulomb interaction @xmath0 has been evaluated numerically  @xcite .
the results seem unsettled , though .
while @xcite suggested the fourth - order term may yield the splitting of the kondo resonance in the intermediate bias region @xmath192 for sufficiently large interaction @xmath193 , which the second - order calculation fails to report , another numerical study indicates that the spectral function remains qualitatively the same with the second - order result  @xcite .
experimental situation is not so transparent , either .
while the splitting of the kondo resonance was reported in a three - terminal conductance measurement in a quantum ring system  @xcite , a similar spitting observed in the differential conductivity was attributed to being caused by a spontaneous formation of ferromagnetic contacts , not purely to bias effect  @xcite .
it is also pointed out that it has been recently recognized that the rashba spin - orbit coupling induces spin polarization nonmagnetically in a quantum ring system with a dot when applying finite bias voltage  @xcite ; hence such spin magnetization might possibly lead to the splitting of the kondo resonance .
the kondo resonance is a manifestation of singlet formation between the dot and the lead electrons .
one may naively think that when several chemical potentials are connected with the dot , such singlet formation would take place at each lead _ separately _ , causing multiple kondo resonances . the results of the multi - terminal ph symmetric siam presented in the previous section tempt us to speculate a different picture . let us suppose that ( almost ) the same dot distribution function @xmath194 is realized for a fixed @xmath185 with a different terminal number @xmath180 , as fig .
[ fig : nldos - by - dmu ] ( a ) suggests . note the assumption is fully consistent to the ward identity for low bias , but it invalidates a quasiparticle ansatz @xmath195 that explicitly depends on @xmath180 . in the large @xmath180 limit with a fixed @xmath185 ,
the effective fermi distribution @xmath188 resembles the fermi distribution at finite temperature @xmath196 .
accordingly bias voltage may well give effect similar to finite temperature .
it is seen in low and large bias limit for a dot with or without the ph symmetry .
it implies that a dot electron can not separately form a singlet with the lead at each chemical potential , because it needs to implicates states at different chemical potentials through coupling with other leads .
our second - order perturbation results seem to support this view .
in summary , we have evaluated analytically the second - order self - energy and green function for a generic multi - terminal single impurity anderson model in nonequilibrium .
various limiting behaviors have been examined analytically .
nonequilibrium spectral function that preserves the current is constructed and is checked to satisfy the spectral weight sum rule .
the multi - terminal effect is examined for the ph - symmetric siam , particularly . within the validity of the present approach
, it is shown that the kondo peak is not split due to bias voltage .
it is found that most of the finite bias effect is similar to that of finite temperature in low- and high - bias limits with and without the ph symmetry .
such nature could be understood by help of the ward identity and the connection through the @xmath197 terminal limit .
the present analysis serves as a viable tool that can cover a wide range of experimental situations .
although there is still a chance that high - order contributions might generate a new effect such as split kondo resonances in a limited range of parameters , it is believed that the second - order perturbation theory can capture the essence of the kondo physics in most of realistic situations . moreover , having a concrete analytical form that satisfies both the current conservation and the sum rule , the present work provides a good , solid , workable result that more sophisticated future treatment can base on .
the author gratefully acknowledges a.  sunou for fruitful collaboration that delivered some preliminary results in this work .
the author also appreciates r.  sakano and a.  oguri for helpful discussion at the early stage of the work . the work was partially supported by grant - in - aid for scientific research ( c ) no .
22540324 and 26400382 from mext , japan .
[ sec : appendix ]
we start with the nonequilibrium green function @xmath198 without the coulomb interaction on the dot . its keldysh structure is specified by @xmath199 where @xmath36 is the effective dot distribution defined in eq .  .
we incorporate the hartree - type diagram into the unperturbed part by replacing @xmath20 to @xmath200 .
note @xmath68 is the exact dot occupation that needs to be determined consistently later . after mapping to the rak representation ,
those components are given by @xmath201    \left [ g^{r}_{\sigma } ( \omega ) - g^{a}_{\sigma } ( \omega ) \right ] .
\label{eq : gk0}\end{aligned}\ ] ] the function @xmath202 reduces to @xmath203 at zero temperature .
taking the fourier transformation of eqs .  , using eq .  , and making further manipulations , we can rewrite @xmath204 and @xmath205 as @xmath206^ { * } ,   \label{eq : pir}\\   \pi^{k}(\varepsilon ) & =   2i\sum_{a , b } \frac{\gamma_{a } \gamma_{b}}{\gamma^{2 } }   \coth      \tfrac{\beta(\varepsilon-\mu_{ab})}{2 }    { \mathop{\mathrm{im}}\nolimits}\left [ \frac{\gamma      b_{ab}(\varepsilon)}{\pi \varepsilon ( \varepsilon+2i\gamma ) } \right ] , \label{eq : pik}\end{aligned}\ ] ] where @xmath207 , @xmath208 is the inverse temperature , and @xmath209 is given by @xmath210 \!\left [    g^{a}(\varepsilon ' ) - g^{r}(\varepsilon'+\varepsilon ) \right].\end{aligned}\ ] ] in the present work , we are interested in the zero temperature limit , for which @xmath211 becomes @xmath212 .
the function @xmath213 in this limit is evaluated as ( with @xmath112 ) @xmath214 this corresponds to a multi - terminal extension of the result obtained for the two - terminal ph symmetric siam .
to complete evaluating the remaining integral over @xmath215 of eqs .
( [ eq : i1],[eq : i2 ] ) , we take full advantage of various properties of dilogarithm function @xmath96 . a concrete integral formula we have utilized
will be given in appendix d. for the sake of completeness , we here collect its definition and properties necessary to complete our evaluation .
[ [ definition ] ] definition + + + + + + + + + +    dilogarithm @xmath96 with a complex argument @xmath216 is not so commonly found in literature .
as it is a multi - valued function , we need to specify its branch structure properly .
one way to define dilogarithm @xmath96 all over the complex plane consistently is to use the integral representation @xmath217 the multivaluedness of dilogarithm @xmath218 originates from logarithm in the integrand .
here we designate the principal value of logarithm as @xmath219 , which is defined by @xmath220 according to eq .  , @xmath96 has a branch - cut just above the real axis of @xmath221 .
accordingly , its values just above and below the real axis are different for @xmath221 : @xmath222 but @xmath223 .
some special values are known analytically .
we need @xmath224 , @xmath225 , @xmath226 and @xmath227 for evaluation later .
[ [ sec : functional - relations ] ] functional relations + + + + + + + + + + + + + + + + + + + +    dilogarithm @xmath96 has interesting symmetric properties regarding its argument @xmath228 ; values at @xmath228 , @xmath229 , @xmath230 , @xmath231 , @xmath232 and @xmath233 are all connected with one another . those points are ones generated by symmetric operations @xmath234 and @xmath235 defined by @xmath236 and @xmath237 forms a group .
other operations correspond to @xmath238 applying a series of integral transformations in eq .
, one can connect the values of dilogarithm at these values with one another  @xcite .
note , those functional relations are usually presented only for real arguments . extending them for complex variables
needs examining its branch - cut structure carefully . by following and extending the derivations in ref   for complex @xmath216
, we prove that the following functional relations are valid for any complex variable @xmath228 .
@xmath239^{2 } ,   \label{eq : dilog - func - rel - start}\\ & { \mathop{\mathrm{li}_{2}}\nolimits}(tz ) =      -{\mathop{\mathrm{li}_{2}}\nolimits}(z ) + \tfrac{\pi^{2}}{6 } - { \mathop{\mathrm{log}}\nolimits}(tz ) { \mathop{\mathrm{log}}\nolimits}z , \\ & { \mathop{\mathrm{li}_{2}}\nolimits}(tstz ) = - { \mathop{\mathrm{li}_{2}}\nolimits}(z ) -\tfrac{1}{2 }      { \mathop{\mathrm{log}}\nolimits}^{2}(stz )   \nonumber \\ & \qquad \qquad \qquad - [ { \mathop{\mathrm{log}}\nolimits}(tz ) + { \mathop{\mathrm{log}}\nolimits}(stz ) ] { \mathop{\mathrm{log}}\nolimits}z , \\ & { \mathop{\mathrm{li}_{2}}\nolimits}\big ( tsz \big ) = { \mathop{\mathrm{li}_{2}}\nolimits}(z ) -      \tfrac{\pi^{2}}{6 }   \nonumber \\ & \qquad \qquad \qquad -\tfrac{1}{2 } { \mathop{\mathrm{log}}\nolimits}^{2}(sz )   - { \mathop{\mathrm{log}}\nolimits}(sz ) { \mathop{\mathrm{log}}\nolimits}(tz ) , \\ & { \mathop{\mathrm{li}_{2}}\nolimits}\big(stz \big ) = { \mathop{\mathrm{li}_{2}}\nolimits}(z ) +      \tfrac{\pi^{2}}{6 }   \nonumber \\ & \qquad \qquad \qquad + \tfrac{1}{2 } { \mathop{\mathrm{log}}\nolimits}^{2}(tz ) + { \mathop{\mathrm{log}}\nolimits}(tz ) { \mathop{\mathrm{log}}\nolimits}(tstz ) .
\label{eq : dilog - func - rel - end}\end{aligned}\ ] ] to our knowledge , the above form of extension of functional relations of dilogarithm has not been found in literature .
[ [ the - taylor - expansion ] ] the taylor expansion + + + + + + + + + + + + + + + + + + + +    to examine various limiting behaviors , we need the taylor expansion of dilogarithm , which is derived straightforwardly from eq .  , @xmath240 the presence of @xmath241 reflects the branch - cut structure of @xmath96 . in particular , we utilize the following expansion in our analysis . @xmath242
here we derive and present the central integral formula for evaluating eqs .
( [ eq : i1],[eq : i2 ] ) . by performing a simple integral transformation in eq .
, we have the integration @xmath243 where all the parameters @xmath244 as well as @xmath228 may be taken as complex numbers .
combined with fractional decomposition , we see the following integral be evaluated in terms of dilogarithm : @xmath245
the remaining task to complete calculating @xmath71 in the form of eqs .
( [ eq : sigmar][eq : xi123 ] ) is to collect all the relevant formulas and organize them in a form that conforms to eq .  .
to write concisely , we introduce the following notations @xmath247 where the hartree level @xmath248 is defined as before .
we express the terms @xmath249 and @xmath250 defined in eqs .
( [ eq : i1],[eq : i2 ] ) as @xmath251 here the leads @xmath95 in @xmath249 as well as @xmath252 in @xmath250 carry spin @xmath253 , while @xmath13 in @xmath250 does spin @xmath16 .
singularity on energy integration is prescribed by the principal values .
equation enables us to perform and express the above integrals in terms of dilogarithm .
the resulting expressions are still complicated , but we can simplify them further using functional relations of dilogarithm eqs  ( [ eq : dilog - func - rel - start][eq : dilog - func - rel - end ] ) .
these require straightforward but rather laborious manipulations . in this way
, we reach the final expression of @xmath81 of eq .  .