quantum dots and single - molecule junctions have been considered as possible building blocks for nano - electronics and quantum - information processing.@xcite recent technological progress has made it possible to manufacture and study electron transport trough ultra - small quantum - dot devices , nanotubes or single molecules.@xcite these devices are designed as a small central region comprising of a quantum dot or a single molecule which is coupled to at least two leads where the finite bias voltage is applied to .
the investigation of such devices is of fundamental importance for our understanding of open quantum systems out of equilibrium .    due to the quantization of the charge ,
the physical properties of such junctions are dominated by many - body effects at temperatures below the charging energy @xmath2 , where @xmath3 is the capacitance of the device .
the experimental devices are often fully controllable by external gate electrodes or elongation of the scanning tunneling microscope tip.@xcite this gives the opportunity to directly study true many - body correlation effects , such as the kondo effect ( see for example ref .   and ) , under the influence of an external bias voltage . however , the theoretical understanding of the interplay between coherent transport favored by many - body correlations and current - driven dephasing at finite bias is still at its infancy and further investigations are needed .
the present work has two main objectives . on the one hand
, we establish the reliability of the recently introduced scattering - states numerical renormalization group ( snrg ) approach@xcite to quantum transport by comparing results for small values of @xmath0 to the diagrammatic kadanoff - baym - keldysh approach , which becomes exact in the limit @xmath1 . on the other hand
, we will discuss discrepancies and reveal shortcomings of those diagrammatic approaches at intermediate values of the charging energy .
we investigate quantum transport through a quantum - dot device using a minimal model@xcite where the complex interacting region is replaced by a single spinful orbital which is coupled to two noninteracting leads . a single coulomb matrix element @xmath0 accounts for the charging energy of the device .
we calculate nonequilibrium spectral functions@xcite and current - voltage ( iv ) characteristics using the snrg as well as different approximations@xcite within the diagrammatic kadanoff - baym - keldysh expansion in the local coulomb interaction @xmath0 .    over the past 40 years
, the keldysh technique  @xcite has proven to be the most successful approach to nonequilibrium dynamics . in the context of quantum transport through nano - junctions direct expansions in the interaction@xcite as well as self - consistent re - summation schemes have been employed.@xcite however , such diagrammatic expansions rely on a small expansion parameter , and are , therefore , confined to weak coupling .
but quantum - impurity models@xcite commonly used in the theory of quantum transport on the molecular level often exhibit infra - red divergences in perturbation theory@xcite which also restrict the diagrammatic keldysh approaches to certain parameter regimes usually to high temperature or to large bias .
in contrast to equilibrium conditions , where complete and accurate solutions can be obtained using a variety of nonperturbative techniques such as the bethe _
ansatz_,@xcite conformal field theory,@xcite or wilson s numerical renormalization group ( nrg ) approach,@xcite techniques for calculating quantum - transport out of equilibrium remain largely at the development stage .
recent advancements on the analytical side @xcite include suitable adaptations of the wegner s flow - equation@xcite and the real - time renormalization - group method.@xcite these methods can successfully access large voltages , but are generally confined to the weak - coupling regime . based on the scattering - states approach to quantum transport@xcite the bethe ansatz was extended to quantum - impurity models out of equilibrium,@xcite but remains limited to a certain class of models .    on the numerical side ,
progress has been made in several directions .
currents have been extracted from time - dependent density matrix renormalization group@xcite calculations using finite 1d wires , and the results agree well with bethe ansatz results for certain models.@xcite quantum monte carlo approaches based on scattering states@xcite can access the intermediate coupling regime@xcite at finite bias .
recent real - time formulations of continuous - time quantum monte carlo@xcite and an iterative real - time path integral approach@xcite to quantum transport offer the appealing advantage of working directly in the continuum limit , but are confined to relatively short time scales .
access to low temperatures and long times is hampered in the former case by a severe sign problem , and by the extrapolation to long memory times in the latter case .
hence neither approach can presently be applied to nonequilibrium dynamics of correlated systems with a small underlying energy scale , as is the case with ultra - small quantum dots when tuned to the kondo regime .
the usage of lippmann - schwinger scattering states has been well established in quantum - field theory@xcite for over 50 years and also successfully adapted to the description of quantum transport through strongly interacting nano - devices coupled to ballistic leads.@xcite these states fulfill the correct boundary condition of the open quantum system : ( i ) they break time - reversal symmetry and , therefore , are ( ii ) complex and current - carrying and ( iii ) describe ballistic transport in the leads combined with scattering events in the small interacting quantum - dot region .
this time - reversal symmetry breaking is required for current carrying systems and reveals itself naturally in all diagrammatic approaches by the occurrence of retarded and advanced green s functions .
it is a consequence of any regularization when performing the limit to an infinitely large system .
in particular , the work of hershfield @xcite and doyon and andrei @xcite has rigorously shown that these boundary conditions remain unaltered when a local interaction is switched on .
the noninteracting current - carrying system evolves into the new steady - state of the _ interacting _ system , and the steady - state density operator retains a boltzmannian form.@xcite the explicit construction of those scattering states allows to exactly solve the dc and ac kondo model at the toulouse point@xcite as well as the interacting resonant level model.@xcite    recently , an extension@xcite to wilson s numerical renormalization group has been developed for steady - state quantum transport through nano - devices which is able to deal with the crossover from weak to strong coupling for arbitrary bias voltages .
it is based on oguri s idea@xcite of discretizing the single - particle scattering states which are the solutions of the lippman - schwinger equation@xcite for the noninteracting problem and , therefore , fulfill the correct boundary condition of an open quantum system .
this scattering - states numerical renormalization group approach@xcite ( snrg ) evolves the analytically known density matrix of a noninteracting system to the density matrix of the fully interacting problem by employing the time - dependent nrg ( td - nrg).@xcite the nrg is ideally suited to the problem , being known to provide accurate solutions of quantum - impurity models on all relevant interaction strengths at zero bias.@xcite since the td - nrg can access exponentially long time scales,@xcite dwell times on the order of the inverse kondo - temperature are easily accessible .
this paper is organized as follows .
after the model used is defined , we provide the details of the different theoretical approaches in sec .
[ sec : theory ] .
we summarize the basic ideas of the snrg method introduced in ref .   in sec .
[ sec : snrg ] and state all necessary equations of the diagrammatic nonequilibrium techniques in sec .
[ sec : kbk - approach ] .
the main body of the paper is in sec .
[ sec : results ] , where we present and discuss the results obtained for the various methods .
in order to set the stage for a detailed comparison between the snrg and diagrammatic approaches at finite bias , we begin with a discussion of the magnetic and charge fluctuation scales at zero bias in sec .  [
sec : equiscale ] .
since the nrg provides an accurate solution in this regime for arbitrary coupling strengths and temperatures , this reveals the validity range of the diagrammatic approaches .
we show that  in contrast to the nrg  some of the diagrammatic expansions fail to produce a single low - energy scale for intermediate and large values of @xmath0 .
however , in the weak correlation regime all these approaches agree excellently for arbitrary voltages at small @xmath0 and yield the same nonequilibrium green functions as well as iv characteristics which are presented in sec .  [ sec : weak - coupling ] .
discrepancies between the different approaches at intermediate values of the coulomb interaction are discussed in sec .
[ sec : inter - coupling ] , where the spectral functions and iv curves of a symmetric and an asymmetric junction are considered .
we conclude with summary and a short outlook in sec .
[ sec : summary ] .
quantum impurity models are used to describe quantum transport on the molecular level .
their hamiltonian @xmath4 @xmath5 consists of three parts : an impurity part @xmath6 modeling the interacting device with a finite number of degrees of freedom , one or several bosonic or fermionic baths represented by @xmath7 , and the coupling of these subsystems by @xmath8 .    throughout this paper , we restrict ourselves to junctions modeled by the single impurity anderson model with one spinful orbital coupled to a left ( l ) and a right ( r ) lead and an on - site coulomb repulsion @xmath0 @xmath9 here , @xmath10 is the single - particle energy of the quantum dot , @xmath11 measures its orbital occupancy and @xmath12 represent the elementary hybridization - matrix elements coupling the dot to the two leads .
the different chemical potentials @xmath13 in both leads appear as a shift of the band centers and are functions of the external voltage @xmath14 .    for simplicity , we assume that both leads have the same density of states , @xmath15 , characterized by the same band width @xmath16 but different band centers .
this hamiltonian is commonly used to model a single coulomb - blockade resonance in ultra - small quantum dots.@xcite        in the absence of the local coulomb repulsion @xmath17 , the single - particle problem is diagonalized exactly in the continuum limit@xcite by the following scattering - states creation operators @xmath18   \punkt\end{aligned}\ ] ] @xmath19 labels left ( right ) moving scattering states created by @xmath20 .
the local retarded resonant level green s function @xmath21^{-1}\end{aligned}\ ] ] enters as an expansion coefficient . defining @xmath22 , we will use @xmath23 and @xmath24 - i\gamma(\w)\end{aligned}\ ] ] in the following .    in the limit of infinitely large leads
 volume @xmath25  the single - particle spectrum remains unaltered , and these scattering states diagonalize the hamiltonian ( [ eqn : siam ] ) for @xmath26 : @xmath27    the scattering states are solutions of the lippmann - schwinger equation@xcite and therefore break time - reversal symmetry , which constitutes a necessary boundary condition to describe a current carrying open quantum system .
this is encoded in the small imaginary part @xmath28 entering eq .
( [ eq : scattering - states - operators ] )  ( [ eq : deltadef ] ) required for convergence when performing the continuum limit @xmath25 in the leads .
the complex expansion coefficients in ( [ eq : scattering - states - operators ] ) are given by _
retarded _ functions , e.g.  @xmath29 , which causes the scattering states to be complex and current carrying . for zero bias voltage ,
time - reversal symmetry manifests itself in the identical spectrum for left and right movers which are time - reversal pairs in that limit .    to avoid any contribution from bound states
, we will implicitly assume a wide band limit : @xmath30 , where @xmath31 and @xmath32 .
hershfield has shown that the density operator for such a noninteracting current carrying quantum system retains its boltzmannian form@xcite @xmath33 even for finite bias .
the @xmath34 operator accounts for the different occupation of the left- and right - moving scattering states , and @xmath13 for the different chemical potentials of the leads .
therefore , all steady - state expectation values of operators can be calculated using @xmath35 which includes the finite bias . in the absence of a coulomb repulsion
@xmath0 , this is a trivial and well - understood problem .
it was shown@xcite that the current expectation value using this density - operator @xmath36 reproduced the standard result@xcite for noninteracting devices .
the knowledge of the analytical form of @xmath35 , however , makes this steady - state model accessible to a nrg approach.@xcite    the expansion coefficients of @xmath37 in eq .
( [ eq : scattering - states - operators ] ) contain the complex single - particle green function @xmath38 which we separate in modulus and phase @xmath39 this phase is absorbed into the new scattering states @xmath40 by a local gauge transformation .
the impurity operator @xmath41 is expanded into left- and right - mover contributions @xmath42 using the inversion of eq .
( [ eq : scattering - states - operators ] ) .
these two new operators @xmath43 are then defined as @xmath44 and obey the anti - commutator relation @xmath45 .
the scattering - states numerical renormalization group approach@xcite ( snrg ) starts from a logarithmic discretization of the scattering - states continuum @xmath46 in intervals @xmath47 $ ] and @xmath48 $ ] @xmath49 , controlled by the parameters@xcite @xmath50 and @xmath51 $ ] .
the intervals for @xmath52 are defined as @xmath53 $ ] and @xmath54 $ ] .
an average over various @xmath55-values@xcite is used to mimic the conduction band continuum .
then , the discretized version of the noninteracting hamiltonian ( [ eq : h0-u0 ] ) is mapped onto a semi - infinite wilson chain    @xmath56    whose tight - binding matrix elements @xmath57 decay exponentially @xmath58 for large @xmath59 .
in contrast to the standard nrg,@xcite the impurity degree of freedom has been included into @xmath60 since not the leads but the full scattering states have been discretized .
any complex phase in the tight - binding parameters @xmath57 can be absorbed into the creation ( anihilation ) operators @xmath61 @xmath62 of an electron on the chain link @xmath59 with spin @xmath63 and mover @xmath64 by a local gauge transformation .
we use @xmath65 defined in eq .
( [ eqn : d - alpha - cont ] ) as starting vector @xmath66 for the householder transformation@xcite and obtain the tight - binding coefficients of the wilson chain ( [ eqn : h0-lambda ] ) by the usual procedure.@xcite it is straight forward to shown that the energy of the first chain link corresponds to the energy of the original quantum - dot orbital : @xmath67 .      in order to include the local coulomb interaction ,
the density operator @xmath68 must be expanded in the new orbitals @xmath65 .
it consist of two contributions : a density term and a backscattering term @xmath69 , where @xmath70 and the backscattering @xmath71 term is defined as @xmath72 the local coulomb interaction @xmath73 @xmath74 leads to a mixing of left and right movers since @xmath71 does not commute with @xmath75 .
however , the term @xmath76 , @xmath77 commutes with @xmath34 and can be absorbed into the steady - state density operator @xmath78/z$ ] with @xmath79 using the arguments given in ref .  .
starting from an equilibrated system for times @xmath80 , the initial hamiltonian @xmath81 is changed to @xmath82 by a sudden quench at @xmath83 .
then , the density operator @xmath84 evolves from its initial value @xmath35 at @xmath83 as @xmath85 if @xmath86 describes a quantum impurity problem and @xmath87 is an impurity operator , it was recently shown that the real - time dynamics of the expectation value of @xmath88 can be calculated@xcite by evaluating @xmath89 where @xmath90 denotes the matrix elements of the operator @xmath87 and @xmath91 the nrg eigenenergy of the eigenstate @xmath92 to @xmath93 at nrg iteration @xmath94 .
the sum restriction @xmath95 indicates that at least one of the states @xmath96 must be a discard state at iteration @xmath94 .
excitations between two retained states will be refined in the following iterations @xmath97 and , therefore , will contribute at a later iteration .
@xmath98 labels the environment degrees of freedom of the wilson chain links to be incorporated in subsequent iterations @xmath97 and the @xmath99 are just tensor - product states of the eigenstates of the @xmath94th iteration and the yet uncoupled rest chain .
the reduced density matrix @xmath100 traces out all environment degrees of freedom @xmath98 .
the initial conditions are encoded into the density operator @xmath35 calculated with the initial hamiltonian @xmath81 .
the calculation of the overlap matrix between the nrg eigenstates of @xmath81 and @xmath82 allows for the basis set transformation of @xmath101 into the basis of the final hamiltonian provided that @xmath35 remains restricted to the last wilson shell.@xcite this transformed @xmath101 enters eq .
( [ eqn : o - vs - t - td - nrg ] ) .
the discarded states form a complete basis set@xcite for the fock - space of the _ entire _ wilson chain of length @xmath102 , i.e.  @xmath103 where @xmath104 labels all discarded states at iteration @xmath94 .
the iterative diagonalization thus procures the set of ( approximate ) eigenstates for the _ whole _ energy range from high energies on the order of the bandwidth down to very low energies such as the kondo scale .
this is indispensable because nonequilibrium processes usually involve all energy scales and can not be confined to a finite low energy window set by the last wilson shell as in the usual equilibrium nrg .      in sec .
[ sec : scattering - state ] we have argued that the analytic form the steady - state nonequilibrium density operator is known for the noninteracting case .
this allows for applying the nrg approach to construct a faithful representation of @xmath105 .
we assume that when switching on the coulomb interaction @xmath73 for infinitely large leads ( i ) a steady state is reached after some characteristic but finite time and ( ii ) it is unique and independent of the initial condition .
as described earlier , the boundary condition of time - reversal symmetry breaking is imposed on the scattering states and the nonequilibrium density operator at @xmath83 for @xmath26 .
the interaction quench at @xmath83 , i.e.  switching on a local scattering potential , and the subsequent unitary time evolution do not affect this boundary condition , and the time - evolved operators characterize the interacting current carrying open quantum system .
the time average of the density operator @xmath106 projects out the steady - state contributions to the time - evolved density operator @xmath107 even in a finite - size system : only the energy diagonal terms contribute in accordance with the steady - state condition @xmath108 = 0 $ ] .
even though @xmath109 remains unknown analytically , we can construct it systematically using the time - dependent nrg@xcite described above .
the steady - state retarded green s function is defined as @xmath110_s } \theta(t ) , \label{eqn : steady - state - gf }    \end{aligned}\ ] ] where @xmath111 , @xmath112_s$ ] denotes the commutator ( @xmath113 ) for bosonic , and the anti - commutator ( @xmath114 ) for fermionic correlation functions .
this green s function can be calculated using the time - dependent nrg@xcite and extending ideas developed for equilibrium green s functions.@xcite the completeness relation for the basis of discarded states as introduced above is given by @xmath115 where @xmath116 denotes the first iteration at which the nrg truncation is employed , @xmath102 is the total number of iterations ( i.e.  the length of the wilson chain ) , and @xmath104 only runs over the states which are discarded at iteration @xmath94 . for each iteration @xmath94
, we can partition the completeness relation ( [ eqn : completness - nrg ] ) into two parts , @xmath117 , where the first part incorporates the iterations @xmath116 to @xmath94 and the second the iterations @xmath118 to @xmath102 . since @xmath119 spans the part of the fock - space which contains all kept states @xmath120 after iteration @xmath94 , the identity @xmath121 must hold .
the different contributions to the green s function are calculated for each energy scale @xmath122 by expanding the ( anti-)commutator in eq .
( [ eq : neqgf ] ) and inserting the completeness relations eqs .
( [ eqn : completness - nrg ] ) and ( [ eqn : kept - states ] ) repeatedly . by making use of the fact that local operators @xmath123 and @xmath124 are diagonal in the environment degree of freedom @xmath98 , reduced density matrices @xmath101 occur naturally when tracing out the environment @xmath98 here as well .
although the excitation energies remain confined to the same energy scale , terms connecting different energy scales @xmath125 and @xmath126 are implicitly included through the reduced density matrices such as defined in ( [ eqn : red - dm ] ) . similar to the real - time dynamics
, the summation over all @xmath94 then ensures that all energy scales @xmath125 contribute to the green s functions.@xcite a detailed derivation is given in ref .  .
it was shown that the algorithm is identical to the equilibrium algorithm@xcite if @xmath127 .
laplace transforming @xmath128 yields the steady - state spectral function for the retarded green s function which is used to calculate the current ( see eq .
( [ eq : ss - current ] ) below ) .
we employ the nonequilibrium perturbation theory as formulated by kadanoff and baym@xcite and keldysh @xcite on the usual keldysh time contour , for example , see refs .  .
since we are only interested in the steady - state properties , the information and correlations of the initial conditions are assumed to be lost .
this is archived by sending the initial time @xmath129 and dropping all correlation functions which involve the initial state .
it is again assumed , that the system reaches a steady state which is translational invariant in time .
therefore , the single - particle green s function does only depend on the difference between the two formerly independent times of particle creation and annihilation .
the laplace transform of the time difference then leads to the formulation in frequency space for all green s functions of the steady state .    in the nonequilibrium steady - state formulation two independent components of the contour ordered green s function survive which are chosen to be the retarded and lesser green s functions , @xmath130 and @xmath131 respectively .
the advanced and greater functions are related via @xmath132    the two relevant green s functions can be expressed as @xcite @xmath133    \\
\feff(\w ) & = &   f_l(\w)\gamma_l(\w)+f_r(\w)\gamma_r(\w ) \label{eqn : f - eff }    \quad,\end{aligned}\ ] ] where , again , @xmath134 are the hybridization functions of the leads , @xmath135 their imaginary parts ( see eq .
( [ eq : deltadef ] ) ) and @xmath136 are the fermi functions of the corresponding leads .
the retarded and lesser self - energies , @xmath137 and @xmath138 respectively , include all correlation effects induced by the coulomb interaction @xmath0 .
@xmath139 accounts for the frequency independent hartree energy shift .
diagrammatic expansions in the coulomb interaction of the self - energy @xcite have been investigated for systems in equilibrium @xcite as well as in nonequilibrium.@xcite the self - energies can be evaluated either non self - consistently , where bare propagators are used as inner lines , or in terms of skeleton diagrams , where fully - dressed propagators are taken into account .    in this study
we focus on three different approximations for the self - energy : ( a ) the bare expansion up to second order in @xmath0 , where hartree - fock ( hf ) propagators are used as internal lines .
the latter are just the noninteracting propagators , but with a shifted level position @xmath140 . this approximation is labeled @xmath141 and its diagrammatic representation is schematically shown in fig .
[ fig : perturb](a ) and ( b ) .
( b ) the self - consistent evaluation of the second - order skeleton diagram of fig .
[ fig : perturb](a ) and ( b ) .
this approximation is called second born approximation ( 2ba ) but in contrast to the usual 2ba , no exchange contribution exists for the single impurity anderson model  ( [ eqn : siam ] ) with only one spinful orbital .
( c ) in the gw approximation@xcite ( gwa ) the bare coulomb interaction @xmath0 is screened by an infinite series of particle - hole excitations , which can be summed as indicated in fig .
[ fig : perturb](c ) .
no contributions with odd orders in the interaction occur in this series due to the definition of the matrix elements of the coulomb interaction in our model  ( [ eqn : siam ] ) , where we set matrix elements between electrons with the same spin explicitly to zero .
the 2ba and the gwa are both evaluated self - consistently , and thus the self - energies can be derived from a luttinger - ward functional.@xcite therefore , both constitute conserving approximations in the sense of kadanoff and baym.@xcite it can be shown that elementary sum rules such as charge and current conservation are obeyed.@xcite in contrast , the non self - consistent @xmath142 approximation is not conserving which can lead to the violation of current conservation , as it will be demonstrated later .
the hartree shift is produced by the average occupation of the quantum dot @xmath143 and analytic expressions for the self - energies read @xmath144 where the effective interactions are given by @xmath145 and the particle - hole bubbles are @xmath146 in the above expressions the advanced and greater green s functions can be determined via eq .
( [ eq : advgf ] ) and ( [ eq : grtgf ] ) and @xmath147 denotes the opposite spin of @xmath148 .
equations ( [ eq : grdef])-([eq : pless ] ) form a closed set , which is solved self - consistently for the 2ba and gwa . for the @xmath142 approximation
all particle - hole propagators ( [ eq : pr])-([eq : pless ] ) are evaluated only once with bare green s functions @xmath149 and eq .
( [ eq : wrba ] ) is used as the effective interaction . the hartree shift is included in order to determine the desired filling .
the effective fermi function @xmath150 was defined in eq .
( [ eqn : f - eff ] ) .
the gwa @xcite has been successfully applied to overcome some shortcomings of local - density calculations and estimate the screening of the coulomb interactions in solid state physics .
recently , it has been employed to calculate quantum transport through nanoscale devices.@xcite in the context of the single impurity anderson model it was shown to accurately describe the equilibrium properties in the weakly - interacting regime and in asymmetric situations with a nearly empty or nearly full impurity orbital.@xcite in the strongly interacting kondo regime , i.e.  @xmath151 , the gwa produces a narrow peak in the spectral function at the fermi level , which could be interpreted as remnants of the expected many - body resonance.@xcite however , the line shape of this low - energy resonance as well as the high - energy hubbard peaks at @xmath152 and @xmath153 are not correctly reproduced by this approximation.@xcite additionally , for very large interactions strength @xmath154 , all three perturbative approaches favor an unphysical magnetic ground state , which is actually forbidden by the mermin - wagner theorem.@xcite in the nonequilibrium situation , the proximity to bifurcation points of these sets of equations leads to unphysical hysteretic response.@xcite      the current flowing from lead @xmath64 onto the impurity region can be expressed as @xcite @xmath155\end{aligned}\ ] ] where @xmath156/\pi$ ] is the frequency and voltage dependent spectral function of the retarded impurity green s function .
since the steady - state current onto the interacting region from the left must be equal to the current leaving to the right lead , i.e.  @xmath157 , we can symmetrize the left and the right currents with a linear combination@xcite and write it as @xmath158 in the wide band limit , @xmath159 and @xmath160 holds such that the term proportional to @xmath161 drops out of ( [ eq : currentsym - def ] ) and we obtain @xmath162   \pi\gamma \rho^r_\s(\w , v )    \ : .
\end{aligned}\ ] ] where we have defined @xmath163 @xmath164 @xmath163 reaches the universal conductance quantum @xmath165 for a symmetric point - contact junction , @xmath166 , and is strongly suppressed in the tunneling regime @xmath167 .    for the voltage drop across the two contacts of the impurity to the leads
we employ a serial resistor model where the chemical potentials in the leads are given by @xmath168 and @xmath169 .    at zero temperature ,
the zero bias conductance @xmath170 is proportional to the spectral function at the fermi level . in the zero temperature fermi liquid and for a symmetric junction @xmath171 .
the conductance is given by its universal value @xmath172 which shows in the slope at zero bias of the iv characteristics , i.e.  @xmath173 .
we also define a leakage current @xmath174 i g^<_\s(\w )          \\
\nonumber   & & + \frac{4e}{h}\sum_\s \int \!d\w\ : \feff(\w ) \pi \rho^r_\s(\w)\end{aligned}\ ] ] which must vanish due to current conservation , @xmath175 , in a physical junction .
therefore , deviations from @xmath176 measures shortcomings of an approximation .
in this section we compare and discuss the results obtained from the different diagrammatic keldysh approaches with the snrg . for simplicity , we used symmetric structureless leads characterized by a constant density of states with a half - bandwidth @xmath177 , i.e.   @xmath178 .
the total @xmath32 is used as the energy scale : all energies , voltages and temperatures are measured in units of @xmath179 throughout the paper .    for the snrg
a rather large @xmath180 was chosen , and @xmath181 states were retained in each nrg - iteration step .
@xmath55-averaging@xcite with either @xmath182 or @xmath183 different @xmath55-values was performed , and the broadening parameter for the spectral function@xcite was chosen @xmath184 . for large @xmath0
some unphysical wiggles may emerge in the spectral function , as it is explained below . in principle , these wiggles can be minimized by choosing a smaller @xmath185 , incorporating more states or performing the @xmath55-averaging with a larger number of @xmath55-values .
we did not include an external magnetic field , and no magnetic solutions are encountered for the parameter values used in this paper .
therefore , we will drop the spin - index from now .
the two spin components of the spectral functions and self - energies are identical , e.g.  @xmath186 and @xmath187 respectively .
before we apply finite bias voltages , we will compare the different equilibrium low - energy scales obtained with the diagrammatic approaches to nrg results .
while the diagrammatic approach becomes exact only in the weak - coupling limit @xmath1 , the nrg produces the correct scales for all interaction strengths .
we will identify the validity range of the diagrammatic expansion . in that regime
the diagrammatic approach produces correct results even in nonequilibrium , and we will therefore use it to benchmark the snrg for finite voltages .
the single impurity anderson model in equilibrium for @xmath188 always forms a local fermi liquid.@xcite the spectral function for a symmetric junction approaches the zero temperature limiting value @xmath189 in accordance with the friedel sum rule.@xcite the fermi - liquid formation is associated with a characteristic low - energy scale , which is identified with the kondo temperature @xmath190 at large coulomb repulsions and near half - filling .
the snrg coincides with the usual nrg@xcite in equilibrium , which accurately describes the crossover from high to low temperatures and provides the correct low - energy scale @xmath190 depending exponentially on @xmath0.@xcite the @xmath142 approximation , however , predicts a low - energy scale which is perturbative in @xmath0 and too large.@xcite the gwa does produce a narrow many - body resonance in the spectral function at the fermi level . extracting a low - energy scale from the full width at half maximum ( fwhm ) for an asymmetric junction ( @xmath191 ) , as shown in fig .  5 of ref .
, suggests an exponential variation with the ionic level position @xmath10 .
however , the exponent has the wrong prefactor as compared to the exact analytic form.@xcite    in order to extract the low - energy scale from our model calculations we employ two different methods : we calculate the temperature dependent zero bias conductance @xmath192 and fit it to a phenomenological form.@xcite since @xmath193 is directly determined by the spectral function , it is sensitive to the amount of spectral weight in the temperature window @xmath194 .
the scale @xmath195 extracted in this way constitutes the energy scale relevant for the zero - bias charge transport in the system .
this procedure yields the same result as the aforementioned extraction from the fwhm of the resonance at the fermi level .
the second way utilizes the screening of the effective local magnetic moment , @xmath196 , where @xmath197 is the magnetic susceptibility , @xmath198 the magnetization and @xmath199 an external magnetic field .
we calculate @xmath198 for a finite but small external magnetic field @xmath200 , and extract the susceptibility via the difference quotient . in the fermi - liquid regime
the effective magnetic moment follows an universal curve as function of temperature from which the low - energy scale is determined by defining @xmath201.@xcite the resulting @xmath202 sets the scale for magnetic excitations in the system and is directly linked to the kondo - screening of the local magnetic moment .    for large values of the coulomb interaction , the scales
@xmath202 and @xmath203 should coincide ( apart from a constant of order one ) and vary as @xmath204 for a symmetric quantum dot . for very small values of the coulomb repulsion @xmath205 both should approach @xmath206 .
the charge scale @xmath195 is expected to be roughly constant and on the order of @xmath207 for @xmath208 since for such small interactions charge fluctuations to and from the leads dominate the physics , and the spectral function stays very close to its hf form . on the other hand ,
the magnetic scale is known to decrease exponentially for all @xmath0 .
@xcite     extracted from the various approximations as described in the text .
a fit to the magnetic scale of the nrg to show the exponential decay @xmath209 is also included in the plot . ]
figure  [ fig : gfu4-ii ] shows the two scales extracted from nrg and gwa calculations for a symmetric junction in equilibrium .
the nrg results show the expected @xmath0-dependencies : the charge scale @xmath195 is on the order of @xmath206 for small @xmath210 and decreases exponentially for large @xmath211 .
the magnetic scale @xmath202 decreases exponentially for all @xmath0 as it is evident from the comparison with a fit function @xmath212 also included in the plot .
furthermore , there exists only one universality scale for large @xmath0 which manifests itself by @xmath213 ( not shown ) .    on the other hand ,
the scales obtained from the gwa agree with the nrg only for small @xmath0 .
the charge scale @xmath195 perfectly agrees with the nrg curve for @xmath214 .
significant deviations are observed for larger @xmath0 , where the gwa-@xmath195 decreases faster than the nrg . for @xmath0
significantly larger than the ones shown in the plot , no scales could be extracted due to the artificial symmetry breaking already reported in the literature.@xcite    we added a second gwa charge scale @xmath215 to the graph which is obtained from the width of the low - energy feature at @xmath216 of @xmath217 ( and not at the fwhm as for @xmath195 ) .
the correspondingly extracted scale should coincide with @xmath195 , apart from a prefactor .
but it is found that both scales follow the same trend only for small @xmath0 and already for @xmath218 a much stronger decrease than the expected @xmath204 is observed in @xmath215 .    )
zero - temperature spectral functions for @xmath219 and @xmath220 and a quantum - point contact @xmath221 .
the inset shows the temperature evolution of the spectral function right at the fermi level , @xmath222 .
the nrg - parameters are @xmath223 , @xmath224 , @xmath183 , @xmath225 and @xmath226 nrg iterations were performed . ]
therfore , the extraction of the charge scale within the gwa at intermediate @xmath0 is somewhat ambiguous .
a comparison of the zero - temperature equilibrium spectral function of the gwa and nrg for @xmath227 is depicted in fig .
[ fig : gfu4-iii ] .
the low - energy feature of the gwa spectral function is too narrow and exhibits a rather spiky line - shape which suggest at too low charge scale .
this is supported by the evolution of @xmath217 as a function of temperature which is shown in the inset of fig .
[ fig : gfu4-iii ] .
the logarithmic increase of @xmath217 which occurs at temperatures on the order of the relevant charge scale also reveals that the charge scale is predicted as too low in gwa compared to the nrg .
however , a considerable broadening occurs away from the fermi level which leads to the same fwhm for the gwa as in the nrg and consequently the larger @xmath195 emerges in thermodynamic quantities like @xmath228 .
the magnetic scale @xmath202 extracted from the gwa exhibits some peculiar @xmath0-dependence . for small @xmath229
the scale agrees with the nrg .
however , it develops a minimum at @xmath230 and then increases again for increasing @xmath0 ! this clearly indicates a failure of the gwa to describe magnetic properties for intermediate and large interactions . since the gwa effective moments @xmath231 show universality as functions of the dimensionless temperature @xmath232 for low temperatures ( not shown ) , the increase in @xmath202 implies a too strong screening of magnetic moments .
the effective coulomb interaction is over - screened by @xmath233 ( see eq .
( [ eq : wrgw ] ) ) .
the electrons remain itinerant even at rather large @xmath0 , and the gwa fails to capture the atomic limit .
therefore , the magnetic screening scale @xmath202 remains large in the gwa and actually increases with @xmath0 .
the scales extracted from the other two diagrammatic approximations all coincide with the gwa for small @xmath234 . for larger @xmath0
the @xmath142 approximation produces the same difference as the gwa between the charge and the magnetic scale , whereas within the 2ba both scales decrease with increasing @xmath0 , but in a polynomial rather than exponential fashion .
we have established that the diagrammatic approach produces reliable results for interactions up to the order of the hybridization strength @xmath234 which we will use in the following section to benchmark the snrg in that regime .
we study the nonequilibrium properties of a symmetric and an asymmetric junction in the weakly correlated regime @xmath236 .
we use a very low temperature @xmath237 , which is sufficiently small compared to all other scales in the problem so it can be considered as @xmath238 with impunity .    for such small interactions the diagrammatic approaches and the snrg yield identical results for all voltages .
figure  [ fig : weakcoupall](a ) shows the nonequilibrium spectral function of a symmetric junction in the quantum - point contact regime , i.e.  @xmath239 and @xmath240 . for small voltages
@xmath241 the spectra are even indistinguishable from the hartree - fock ( hf ) result . only at larger @xmath242 small deviations around the fermi level as can be observed .     and @xmath243 for different voltages .
results for the spectral functions are shown for hf , gwa and snrg , while the self - energy is shown for the gwa only .
( b ) current as a function of voltage for a symmetric ( @xmath244 and @xmath243 ) and an asymmetric junction ( @xmath245 , @xmath246 and @xmath247 ) .
@xmath248 is normalized to @xmath249 ( symmetric ) and @xmath250 ( asymmetric ) and measured in units of @xmath179 .
a small temperature of @xmath251 was used for all calculations .
parameters for the snrg calculations are @xmath180 , @xmath252 , @xmath182 , @xmath253 and 12 nrg - iterations were performed .
, title="fig : " ]   and @xmath243 for different voltages .
results for the spectral functions are shown for hf , gwa and snrg , while the self - energy is shown for the gwa only .
( b ) current as a function of voltage for a symmetric ( @xmath244 and @xmath243 ) and an asymmetric junction ( @xmath245 , @xmath246 and @xmath247 ) .
@xmath248 is normalized to @xmath249 ( symmetric ) and @xmath250 ( asymmetric ) and measured in units of @xmath179 .
a small temperature of @xmath251 was used for all calculations .
parameters for the snrg calculations are @xmath180 , @xmath252 , @xmath182 , @xmath253 and 12 nrg - iterations were performed .
, title="fig : " ]    the imaginary part of the retarded self - energy @xmath254 $ ] for that junction obtained with the gwa is shown in the inset of fig .
[ fig : weakcoupall](a ) for various voltages . the overall scale of @xmath254 $ ] is much smaller than @xmath206 and , therefore , the total self - energy @xmath255 is dominated by the charge - fluctuation scale @xmath206 .
but the general influence of a finite bias voltage can already be observed here : the quasiparticle scattering amplitude is increased by the interplay between the voltage - induced fluctuations and the interaction .
the characteristic fermi - liquid quadratic minimum in @xmath254 $ ] at the fermi level is destroyed with increasing voltage and the local fermi liquid prevailing in equilibrium ( @xmath256 ) is suppressed at large enough bias .
the evolution of the minimum in @xmath254 $ ] with voltage bears some resemblance with a temperature evolution .
the quasi - particle coherence is destroyed by a finite voltage in a similar fashion as with increasing temperature .
the resulting iv characteristics of a symmetric and an asymmetric junction are shown in fig .
[ fig : weakcoupall](b ) for @xmath257 .
the current is normalized to @xmath258 and measured in units of @xmath206 .
the rescaled current always saturates at @xmath259 for large voltages independent of @xmath0 ( not shown ) , as required by eq .
( [ eq : currentsym ] ) .
the initial slope at zero voltage of the iv curve remains unaltered for all values of @xmath0 , in accordance with the fermi liquid nature of the model at small bias and zero temperature .    while a symmetric junction with symmetric coupling to the leads always results in symmetric spectral functions , @xmath260 , and antisymmetric iv characteristics , @xmath261 ( see eq .
( [ eq : currentsym ] ) ) , an asymmetric junction in combination with asymmetric coupling yields a non - antisymmetric iv characteristics with @xmath262 .
this is clearly visible in fig .
[ fig : weakcoupall](b ) .
the bias window ranges from @xmath263 to @xmath264 and is not symmetric around the fermi level . in combination with the shift of spectral weight to higher energies in @xmath265
the center of the spectral function is at @xmath266  this leads to a smaller contribution to the current for negative voltages .
in contrast to the spectral functions , the iv characteristics of the snrg and gwa agree perfectly with the hf results for all voltages . the current is rather insensitive to the detailed distribution of spectral weight and measures only the total amount in the bias window @xmath267 $ ] .
the snrg produces the correct results for small values of the interaction @xmath229 , and has thus no principal limitations .
therefore , the expectation that it is reliable at arbitrary interaction strengths as well is warranted .      as demonstrated in the previous section
the interaction plays a minor role for small @xmath269 . on the other hand , with an odd number of electrons on the quantum dot and at very large @xmath269 and @xmath270 , the system develops a kondo effect as @xmath188
( see for example ref .   and ) .
the snrg was shown @xcite to correctly describe the strongly correlated kondo regime out of equilibrium .
the enhancement of the conductance in the coulomb - blockade region was reproduced for small bias , and the destruction of the many - body resonance at the fermi level with increasing voltage has been studied .    in this section
we will focus on the intermediate interaction regime , where correlations become increasingly important . since the diagrammatic keldysh approximations already show deficiencies in equilibrium
 see , for example , in sec .
[ sec : equiscale ] or ref .
 discrepancies will extend to finite voltages .
and @xmath271 at @xmath272 .
the bias voltages are @xmath273 ( a ) , @xmath274 ( b ) and @xmath275 ( c ) .
snrg - parameters are @xmath180 , @xmath181 , @xmath183 , @xmath225 and 8 nrg - iterations were performed .
, title="fig : " ] +   and @xmath271 at @xmath272 .
the bias voltages are @xmath273 ( a ) , @xmath274 ( b ) and @xmath275 ( c ) .
snrg - parameters are @xmath180 , @xmath181 , @xmath183 , @xmath225 and 8 nrg - iterations were performed .
, title="fig : " ] +   and @xmath271 at @xmath272 .
the bias voltages are @xmath273 ( a ) , @xmath274 ( b ) and @xmath275 ( c ) .
snrg - parameters are @xmath180 , @xmath181 , @xmath183 , @xmath225 and 8 nrg - iterations were performed .
, title="fig : " ]    the nonequilibrium spectral function for various voltages and intermediate interaction @xmath276 is shown in fig .
[ fig : u4-i](a)-(c ) for a symmetric junction ( @xmath277 ) with symmetric coupling to the leads ( @xmath240 ) at @xmath272 .
now all diagrammatic approximations yield different results .    at low voltages ,
the snrg and @xmath142 approximation reproduce the slight humps at energies @xmath278 which are the first indicators of upper and lower hubbard satellites forming at large @xmath0 .
the gwa only produces the broad high - energy tails , without the indication of forming separate peaks and the 2ba completely fails to produce the enlarged spectral weight at high energies .
as the voltage is raised , the coulomb interaction causes additional dephasing , leading to increasingly broadened spectra .
however , the @xmath142 approximation produces systematically too broad high - energy tails and an unphysical plateau around @xmath279 , which even develops a slight dip as seen in fig .
[ fig : u4-i](c ) .
we attribute this to a tendency to overestimate the coulomb repulsion .
this might already be guessed from the equilibrium spectral functions , where the @xmath142 approximation unexpectedly produces the high - energy hubbard satellites for arbitrary large coulomb repulsion .
these are connected to the ionic many - body states of the isolated atom which are not expected to be described by a second - order perturbation theory .
however , the analytic structure of the retarded self - energy , eq .  ( [ eq : sigr ] ) , ( [ eq : wrba ] ) , ( [ eq : pr ] ) and ( [ eq : gfu0 ] ) , has two direct consequences : ( i ) for small coupling to the leads or large @xmath0 it favors a @xmath280 behavior .
this results in a two - peak structure with the peak - positions and widths roughly given by @xmath281 and @xmath206 , respectively .
( incorporating a screened and dynamic coulomb interaction , as it is done in the 2ba and gwa , leads to a prefactor smaller than @xmath282 and additional imaginary parts enter in the frequency dependence of @xmath283 .
the hubbard satellites are then moved to lower energies and broadened . )
( ii ) the fermi functions entering eq .
( [ eq : pr ] ) through @xmath284 lead to a narrowing of the integration interval for decreasing temperature . at zero temperature
this always produces a vanishing imaginary part of the self - energy at the fermi level , @xmath285=0$],@xcite given that the noninteracting propagators are non - singular at @xmath279 .
( this reasoning also holds for the 2ba and gwa ) .
the combination of ( i ) and ( ii ) gives rise to the two - peak structure in the spectral function for large @xmath0 and the emergence of an additional peak at the fermi level for low temperatures which is usually interpreted as the kondo resonance .
but in principle there is no justification why the @xmath142 approximation should be reliable for large values of @xmath0 under arbitrary conditions .
already for the asymmetric model in equilibrium the phase - space argument ( ii ) does not guarantee the correct description of the low - temperature fermi liquid anymore , and it is well - known that the @xmath142 approximation produces unphysical results.@xcite therefore , the large differences to all other methods at finite voltages , as it is observed here for @xmath219 , is not surprising .
the snrg tends to produce additional features in the spectral function for large bias voltages at the positions of the chemical potentials of the leads , @xmath286 .
these are the humps visible for larger voltages in the curves of fig .
[ fig : u4-i ] .
they are artifacts of the nrg discretization and dependent on the broadening procedure of the nrg spectral functions.@xcite in equilibrium , the nrg only resolves spectral information above a cutoff frequency @xmath287 where @xmath288 is on the order of the temperature @xmath289 .
the nrg broadening parameters @xcite are usually adjusted such that artefacts are minimized .
additionally the spectral function is interpolated between @xmath290 .
this translates itself to the present implementation of the snrg which does not provide spectral information in the intervals @xmath291 $ ] centered around the two chemical potentials . here , @xmath292 and @xmath293 .
furthermore , the time - dependent nrg introduces additional discretization errors@xcite which increase with increasing value of @xmath0 .
@xmath55-averaging over different discretizations@xcite improves the spectral functions and these artifacts could be removed by adjusting the broadening parameter depending on the voltage . in this paper , however , we keep the broadening@xcite parameter fixed at @xmath184 independent of the bias and performed @xmath55-averaging with @xmath182 and @xmath294 .
let us focus on the different behavior of the spectral functions around @xmath279 depicted in fig .
[ fig : u4-i](a)-(c ) .
the height @xmath217 is reduced for increasing @xmath242 , and the spectral functions of the snrg , 2ba and gwa approach each other and eventually coincide . for @xmath275 ,
the 2ba curve is still a little higher compared to snrg and gwa , but at even larger voltages ( not shown ) it also falls on top of the snrg and gwa . in contrast , the @xmath142 spectral function does not approach this large voltage limit , and the zero frequency value of the spectral function is considerably reduced compared to the other approaches .
this is in accord with an enhanced scattering amplitude at the fermi level , visible in the imaginary part of the self - energy depicted in the insets .    upon increasing the voltage ,
the @xmath142 approximation does not follow a systematic trend since @xmath217 is larger than the snrg at low voltages and smaller at high @xmath242 .
the other approximations show systematic deviations as the 2ba is always larger than the snrg , while the gwa is always smaller .    in the present calculations ,
the temperature @xmath272 is only about on fifth of the equilibrium kondo temperature for these parameter values , i.e.  @xmath295 . as discussed in sec .
[ sec : equiscale ] , the gwa produces a too small charge scale , which results in an even higher effective temperature .
this leads to a reduction of the spectral function at the fermi level in addition to the effect of the small bias .
the imaginary part of the self - energy is correspondingly too large compared to the snrg , as can be seen in the inset of fig .
[ fig : u4-i](a ) .
the 2ba , on the other hand , overestimates the low - energy scale , which explains the trends in @xmath217 and @xmath296 $ ] .
increasing the current through the junction by applying a larger bias enhances the charge fluctuations on the local orbital . as already mentioned in section [ sec : weak - coupling ] , theses additional fluctuations introduce dephasing@xcite and destroy the coherent quasiparticles which constitute the low - temperature fermi liquid .
the accompanying destruction of the characteristic quadratic minimum in @xmath297 $ ] around @xmath298 is observed in the insets .
the system is driven away from the equilibrium fermi - liquid fixed point , and the spectral functions at the fermi level decreases .
at very large voltage @xmath299 the coherent quasiparticles are completely suppressed , as it can be seen from the large imaginary part of the self - energy around @xmath300 .
the spiky features in the snrg self - energy for @xmath275 are due to the aforementioned discretization errors and have no physical meaning
.       obtained from the spectral functions presented in fig .
[ fig : u4-i ] for a quantum - point contact with @xmath243 .
the current is normalized to @xmath249 and measured in units of @xmath179 . ]
a coulomb interaction @xmath219 leads to a reduction of the current compared to its hf value as depicted in fig .
[ fig : u4-iv ] .
this is characteristic for the onset of the coulomb blockade .
all approaches predict this reduction but slight differences can be noticed .
the @xmath142 approximation overestimates the coulomb blockade resulting in a current which is systematically smaller than the snrg result . even though the @xmath142 spectral function differs strongly from all other approaches for large @xmath242 , this failure to describe the correct single - particle dynamics
is concealed in the current as all approaches yield identical results .
it again shows the insensitivity of the current to the detailed distribution of spectral weight in @xmath265 .
the 2ba slightly overestimates the current for intermediate voltages , which is again explained by the too large low - energy scale @xmath195 and the accompanying underestimation of correlation effects .
the gwa current merges with the snrg result for @xmath301 , which  together with the satisfactory spectral function for these voltages  indicates a good description of the nonequilibrium properties for intermediate to large @xmath242 .      the spectral function for a quantum dot with a level position @xmath302 , coulomb interaction @xmath276 and asymmetric coupling @xmath247 is shown in fig .
[ fig : u4asymgf ] .
the asymmetry between positive and negative voltages is directly visible in the spectral functions .    , @xmath302 , with asymmetric coupling @xmath247 for small ( upper ) and large ( lower ) bias , @xmath303 and @xmath304 , respectively .
the vertical dashed lines indicate the location of the left and right chemical potentials .
the inset shows the total occupation @xmath305 of the impurity as a function of the bias voltage obtained with the diagrammatic approaches .
snrg - parameters are the same as for fig .
[ fig : u4-i ] . , title="fig : " ] , @xmath302 , with asymmetric coupling @xmath247 for small ( upper ) and large ( lower ) bias , @xmath303 and @xmath304 , respectively .
the vertical dashed lines indicate the location of the left and right chemical potentials .
the inset shows the total occupation @xmath305 of the impurity as a function of the bias voltage obtained with the diagrammatic approaches .
snrg - parameters are the same as for fig .
[ fig : u4-i ] . , title="fig : " ]    apart from voltage - induced broadening which was already discussed in previous sections , an additional shift of spectral weight in @xmath265 is observed with increasing @xmath242 . while the snrg moves spectral weight to higher energies for positive bias and towards @xmath279 for negative @xmath242 ,
the diagrammatic approaches produce the opposite trends .    due to the stronger coupling to the left lead , @xmath306 ,
the left - moving scattering states dominate the excitations and the spectral function of the impurity orbital , as can be seen from eq .
( [ eqn : d - lr - expand ] ) and  ( [ eqn : d - alpha - cont ] ) . the effective noninteracting single - particle excitation energy of an @xmath64-mover is given by @xmath307 , which implies almost symmetric parameters for the left - movers at the negative voltage @xmath308 since then @xmath309 .
therefore , the asymmetry of the spectral function is expected to be reduced and @xmath265 to be closer to that of a symmetric junction , which is indeed observed in the snrg .
the diagrammatic approaches underestimate correlations in the ionic many - body states , and the occupancy of the impurity is overestimated for negative voltages .
it increases almost linearly with negative voltage as can be seen from the inset of fig .
[ fig : u4asymgf](b )
. therefore , the hartree shift  ( [ eq : sighartree ] ) also increases , and spectral weight is moved towards higher energies opposite to what would be expected from the physical argument presented above . as a consequence
the spectral functions are strongly attracted to @xmath310 .
the effective single - particle excitation energy of a left - moving scattering state for a positive voltage @xmath275 is greater than zero , @xmath311 .
this produces an intermediate valence situation for the left - movers , where correlations renormalize the effective excitation energies to even larger frequencies@xcite and a shift of spectral weight to higher energies results .
an additional drag of spectral weight towards the chemical potential of the weaker coupled right lead , @xmath312 , is expected .
the snrg produces such a shift as can be seen in fig .
[ fig : u4asymgf](b ) . the diagrammatic approaches , however , underestimate the level - renormalization in the presence of strong valence fluctuations , a tendency already observable in equilibrium ( not shown ) .
additionally , the reduced occupancy ( inset ) diminishes the hartree energy which again leads to a shift towards the stronger coupled chemical potential @xmath313 .    ,
@xmath276 and @xmath247 , calculated with the spectral functions depicted in fig .
[ fig : u4asymgf ] .
the result for the symmetric junction already shown in fig .
[ fig : u4-iv ] is included for comparison .
( b ) the leakage current @xmath314 for the keldysh approaches obtained from eq .
( [ eq : currentleak ] ) .
the currents are normalized to @xmath250 and measured in units of @xmath179 .
, title="fig : " ] , @xmath276 and @xmath247 , calculated with the spectral functions depicted in fig .
[ fig : u4asymgf ] .
the result for the symmetric junction already shown in fig .
[ fig : u4-iv ] is included for comparison .
( b ) the leakage current @xmath314 for the keldysh approaches obtained from eq .
( [ eq : currentleak ] ) .
the currents are normalized to @xmath250 and measured in units of @xmath179 .
, title="fig : " ]    figure  [ fig : u4asymiv ] shows the iv characteristics of this junction .
the asymmetry of @xmath315 is clearly visible when compared to the result from the symmetric junction ( also included in the plot ) . for @xmath316 ,
the rescaled current is very close to its values from the symmetric junction and the discrepancies between the 2ba , gwa and snrg follow the already discussed characteristics : the 2ba underestimates correlations and yields a slightly too large current .
the gwa has the correct distribution of spectral weight in the bias window and produces a rather good estimate for the current , despite its deficiencies in the description of the single - particle spectra .
the failure to produce the correct shifts of spectral weights in @xmath265 causes the current in the diagrammatic approaches to be smaller than the snrg for negative voltages .
the @xmath142 approximation , however , reveals its non - conserving nature in the violation of current conservation for this asymmetric junction .
this is illustrated in fig .
[ fig : u4asymiv](b ) displaying the leakage current @xmath317 of eq .
( [ eq : currentleak ] ) .
in contrast to the conserving 2ba and gwa methods and the snrg , @xmath317 does not vanish for the @xmath142 approximation !
thus , left and right current of eq .
( [ eq : currentsym - def ] ) do not have the same magnitude , i.e.  @xmath318 , and the current calculated from eq .
( [ eq : ss - current ] ) does not make sense , since different linear combinations @xmath319 ( @xmath320 ) yield different results .
therefore , we did not include the calculated iv curves in fig .
[ fig : u4asymiv](a ) .
increasing the asymmetry further , i.e.  @xmath321 , recovers the equilibrium spectral functions of a quantum dot coupled to a single lead in all approaches ( not shown ) . in the snrg , the backscattering term @xmath71 , eq .
( [ eq : back - scattering ] ) , is suppressed , and the model approaches an equilibrium single - channel problem . in the diagrammatic approaches , the nonequilibrium conditions enter only through the effective fermi function @xmath322 , eq .
( [ eqn : f - eff ] ) , which approaches its equilibrium value for @xmath323 . in this regime ,
all differences in the spectral functions of the presented approaches are given by the known discrepancies already present in equilibrium .
in the recently developed snrg approach to open quantum systems the scattering states of a noninteracting quantum impurity model are used to construct the nonequilibrium green s functions for the steady state at finite bias voltage .
we have established the reliability of the snrg by benchmarking it against the diagrammatic kadanoff - baym - keldysh approach , which becomes exact in the limit @xmath1 .
it has been shown that the spectra and the current - voltage characteristics agree excellently for small coulomb interactions for symmetric and asymmetric junctions at arbitrary bias voltage .    for intermediate values of @xmath0
we have compared the snrg to three different approximations obtained from the keldysh approach , namely the second - order perturbation theory ( @xmath142 ) , the fully self - consistent second - order ( 2ba ) and the gw approximation ( gwa ) . as correlation effects
play an increasingly important role discrepancies occur between the different methods .
these were explained by the insufficient treatment of the coulomb interaction within the diagrammatic approaches .
the fermi liquid at zero bias voltages is characterized by a single low - energy scale which is captured accurately by the snrg , but is not properly reproduced by the diagrammatic approaches .
no single low - energy scale can be extracted from the @xmath142 approximation and the gwa at intermediate and large @xmath0 .
while the scale associated with charge fluctuations decreases with increasing @xmath0 , the magnetic scale exhibits a qualitatively different @xmath0-dependency , since it develops a minimum for intermediate @xmath0 and increases again towards larger @xmath0 ! the gwa shows a tendency to over - screen magnetic moments for increasing values of @xmath0 and fails to reproduce the atomic limit .
these deficiencies translate themselves to finite bias and explain the discrepancies at small to intermediate voltages .    at large bias
voltages the self - consistent diagrammatic approaches ( 2ba and gwa ) reproduce the snrg spectral functions for a symmetric junction , while the second - order perturbation theory yields an unphysical plateau around the fermi level .
all diagrammatic approximations and the snrg capture the onset of the coulomb blockade in the iv characteristics of the symmetric junction .
the small discrepancies are explained by the deficiencies in the treatment of the interaction .
however , the failure of the @xmath142 approximation to correctly describe the single - particle dynamics at large bias is masked in the current , since there only the total spectral weight in the bias window enters .
in contrast to the other methods , the @xmath142 approximation reveals its non - conserving nature by producing a finite leakage current for an asymmetric junction , which is unphysical .
this raises the question about the reliability of the results obtained within that method or extensions of it,@xcite even for a symmetric junction .
they are only well - justified for cases where @xmath324 for all frequencies .
the voltage dependent redistribution of spectral weight for an asymmetric junction is not well - reproduced by the diagrammatic approaches .
this has been attributed to too large hartree shifts due to the wrong occupation number of the impurity and the inaccurate renormalization of the single - particle level in intermediate - valence situations .
it leads to the underestimation of the current for large negative voltages .
the snrg provides access to the description of nonequilibrium steady - state properties of nanoscale junctions for arbitrary coulomb interaction and voltages .
it opens promising perspectives for future investigations , such as the influence of charge fluctuations when approaching the strongly - correlated regime or the effects of an applied magnetic field .
we are grateful to j.  e.  han , a.  millis , a.  schiller , p.  schmitteckert , g.  schn , h.  schoeller , p.  werner , m.  wegewijs and g.  zarand for helpful discussions .
we acknowledge financial support from the deutsche forschungsgemeinschaft under an 275/6 - 1 and supercomputer support by the nic , forschungszentrum jlich under project no .  hhb000 .