the local coulomb interaction @xmath6 of electrons occupying quantum dots leads to a variety of effects .
many of them can conveniently be studied in transport through the dot within the linear regime .
theoretically as well as experimentally well - investigated examples are the coulomb blockade ( cb ) peaks of the transmission ( conductance ) @xcite as well as the plateaus of width @xmath2 of the transmission ( conductance ) induced by the kondo effect @xcite .
additional features of interacting multi - level dots that have recently attracted considerable theoretical attention are the population inversions of the dot levels @xcite , the phase lapses of the transmission phase @xmath1 or , equivalently , the zeros of the transmission amplitude @xmath7 ( transmission zeros ) @xcite and additional correlation - induced resonances of @xmath8 @xcite .
they appear in certain parts of the parameter space when the level occupancies and the transmission amplitude are investigated as functions of the level positions , which can be tuned via a nearby plunger gate voltage .
such effects were hitherto mostly studied in a minimal model involving only two levels . a very important step towards a unified understanding of population inversions , phase lapses and correlation - induced resonancess in spin - polarized two - level dots
was recently taken by a multi - stage mapping of the problem on a generalized kondo model and a subsequent renormalization group and bethe ansatz analysis of the effective hamiltonian @xcite .
theoretical studies of phase lapses ( transmission zeros ) are of primary interest in connection with a series of linear response transmission measurements by the weizman group @xcite on aharonov - bohm rings containing a quantum dot in one arm . under suitable conditions
both the phase @xmath1 and magnitude @xmath8 of the transmission amplitude @xmath9 of the dot can be extracted from the aharonov - bohm oscillations of the current through the ring @xcite . when this is done as function of a plunger gate voltage @xmath10 that linearly shifts the dot s single - particle energy levels downward , @xmath11 ( @xmath12 is a level index ) , a series of well - separated cb peaks of rather similar width and height was observed in @xmath13 , across which @xmath14 continuously increased by @xmath0 , as expected for breit - wigner - like resonances . in each cb valley between any two successive peaks , @xmath1 always jumped sharply downward by @xmath0 .
this phase lapse behavior was found to be `` universal '' , occurring in a large succession of valleys for every _ many_-electron dot studied in refs .
this universality is puzzling , since naively the behavior of @xmath15 is expected to be `` mesoscopic '' , i.e.  to show a phase lapse in some cb valleys and none in others , depending on the dot s shape , the parity of its orbital wavefunctions , etc . only recently @xcite , also
the _ few_-electron regime was probed experimentally : as @xmath10 was increased to successively fill up the dot with electrons , starting from electron number @xmath16 , @xmath14 was observed to behave mesoscopically in the _ few_-electron regime , whereas the above - mentioned universal phase lapse behavior emerged only in the _ many_-electron regime @xmath17 .    it was suggested in ref .
@xcite that a generic difference between the few- and many - electron dots may be that for the latter , transport might simultaneously occur through several partially filled single - particle levels in parallel .
a possible reason could be that the mean ( non - interacting ) level spacing @xmath4 of the topmost filled levels decreases as the number of electrons increases , while the charging energy @xmath2 still implies well separated transmission resonances @xcite .
this scenario forms the basis of a recent systematic study by us of the interplay of level spacing , level width and charging energy on the phase lapses for up to four interacting levels and spin - polarized electrons @xcite .
we showed that the universal phase lapse and transmission zero behavior appearing at small @xmath4 can be understood as resulting from a fano - type interference effect @xcite involving transport through two or more effective dot levels , whose position and width have been renormalized by the coulomb interaction and coupling to the leads .
the importance of several overlapping levels for phase lapses had earlier been pointed out by silvestrov and imry in a rather specific model of a single wide and several narrow levels with strong interaction @xcite .    here
we supplement our earlier study @xcite by discussing the relation between phase lapses and population inversions and by investigating the role of finite temperatures @xmath18 as well as spin , focusing on @xmath19 levels .
when spin is included , the kondo effect plays a role for an odd average occupation of the dot , but the phase lapse scenario is unaffected by this .
experimentally the behavior of the phase in the presence of the kondo effect was investigated in refs .
@xcite and @xcite . as in ref .
@xcite we are concerned with the _ generic _ behavior and thus investigate the entire parameter space , going beyond subspaces of higher symmetry [ such as left - right ( l - r ) or 1 - 2 symmetry of the coupling between the left and right leads and the two levels ] . for low temperatures and sufficiently small single - particle spacings @xmath4 of the dot levels
, we find @xmath0 phase lapses between two transmission peaks in an overwhelmingly large part of the parameter space of the level - lead couplings .
we point out that the two level case is special compared to models with @xmath20 , as for generic level - lead couplings a transmission zero and phase lapse occurs between the two transmission peaks _ even at _ u=0 . the effect of the interaction is merely to increase the separation of the transmission peaks . for large @xmath4 the appearance or not of phase lapses between transmission peaks
depends on the relative sign of the level - lead couplings in analogy to the noninteracting case @xcite .    for spin - polarized dots
we in addition compare our @xmath21 results with the ones of recent mean - field studies @xcite . in these works level - lead couplings beyond the subspaces with increased symmetry were studied , and a remarkably more complex behavior was found once the symmetries were broken .
the importance of considering such generic parameter sets was independently pointed out in ref .
we here elucidate how the phase lapse behavior is affected by quantum fluctuations , which are expected to be strong in low - dimensional systems .
we find that upon including quantum fluctuations , the part of parameter space exhibiting universal @xmath0 phase lapses between well separated cb peaks becomes larger than suggested by the mean - field study .
in particular , we do not recover certain peculiar features of the mean - field results of refs .
@xcite and @xcite , namely the occurrence , in certain regimes of parameter space , of a phase lapse of less than @xmath0 ( instead of precisely @xmath0 ) , accompanied by the disappearance of the corresponding transmission zero @xcite .
these features thus turn out to be artifacts of the mean - field approximation , which misses the rather simple scenario for the phase lapse behavior of a two - level dot at small @xmath4 : for generic level - lead couplings a phase lapse and transmission zero between two transmission peaks is already present at @xmath5 ; increasing the coulomb interaction the peaks become well separated while the phase lapse and transmission zero remain in the valley between them .    in the model of a single wide and several narrow levels
@xcite a relation between phase lapses and population inversions was discussed .
therefore , in phase lapse studies quite often also the level occupancies @xmath22 , @xmath23 , are investigated .
we emphasize that the generic appearance of a phase lapse and transmission zero even at @xmath5 renders the two - level model unsuitable for establishing a general relation between phase lapses and population inversions , as the latter only appear at sufficiently large @xmath2 .
furthermore , we show that discontinuities of the @xmath22 as a function of @xmath10 are an artifact of the mean - field solution ( see refs .
@xcite and @xcite ) . within our approaches
discontinuities are only found for l - r symmetric level - lead couplings with a relative plus sign of the underlying hopping matrix elements and degenerate levels , a case which was earlier identified as being nongeneric @xcite , because the transmission shows only a single peak .
this paper is organized as follows . in sec .
[ model ] we introduce our model for the spin - polarized and spin - degenerate two - level dot .
we discuss the relation between the measured magnetic flux dependence of the interferometer s linear conductance and the magnitude and phase of the dot s transmission amplitude .
the latter can be computed from the one - particle green function of the dot .
we present a brief account of our techniques to obtain the latter , the numerical ( nrg ) @xcite and functional renormalization group ( frg ) methods . for an introduction to the use of the frg to quantum dots
see refs .
@xcite and @xcite .
in contrast to the implementation of the nrg used in ref .
@xcite , where we were restricted to l - r symmetric dots at @xmath21 , we have now adopted the full density matrix ( fdm ) nrg method of @xcite , which enables us to investigate dots with arbitrary level - lead state overlap matrix elements @xmath24 ( with @xmath25 ) as well as to study finite temperatures . in secs .
[ resultsnon]-[resultsspin ] we present our results of the @xmath10 dependence of @xmath8 and @xmath1 .
first we briefly discuss the noninteracting two - level dot with generic level - lead couplings and point out that the phase lapse scenario differs from the one for more than two levels .
we then investigate interacting , spin - polarized dots , study the relation between phase lapses and population inversions and compare to the mean - field results for the phase lapses .
the issue of continuous versus discontinuous @xmath10 dependence of the level occupancies @xmath26 and @xmath27 is commented on .
next we study the role of finite temperatures . finally , we consider spin - degenerate levels at small @xmath3 which implies the appearance of kondo physics at odd average dot filling . using nrg and frg we show that the spin does not alter the universal phase lapse scenario .
our findings are summarized in sec .
[ summary ] .
in this section we introduce our model for the two - level dot .
we argue that it is the energy dependent ( effective ) transmission amplitude @xmath28 which one has to compute if one is interested in comparing to the measurements of refs .
@xcite of the magnitude of the transmission amplitude and its phase .
the amplitude @xmath29 can be determined from the matrix elements of the dot s interacting one - particle green function .
we furthermore discuss aspects of the nrg and the frg specific to our problem .
our hamiltonian consists of three parts @xmath30 the two semi - infinite leads are modeled as noninteracting one - dimensional tight - binding chains and for simplicity are assumed to be equal @xmath31 the hopping strength in the leads is @xmath32 .
we use standard second quantized notation with @xmath25 indicating the left and right leads , where the quantum numbers @xmath33 and @xmath34 label wannier states and spin , respectively .
the dot is described by @xmath35 where the term with @xmath36 and @xmath37 is excluded from the sum in the interacting part .
we define @xmath38 . in experimental systems the inter- and intra - level coulomb repulsion can be expected to be comparable in size and to avoid a proliferation of parameters we assumed them to be equal .
this assumption is not essential ; by relaxing it we have checked that our results are robust against inter - level variations of the interaction strengths .
finally , the coupling between dot and lead states is given by @xmath39 with real overlap matrix elements @xmath24 .
for later purposes we define @xmath40 and @xmath41 .    for simplicity , part of our studies will be performed on a model of spinless electrons , for which the spin index will be dropped .
the resulting model may be regarded as a spin - polarized version of the spinful model obtained if the latter is put in a very large magnetic field .
the experimental two - path interferometer is characterized by wide base regions and very narrow point contacts towards the emitter and collector .
electrons can thus only pass once through either the upper or lower arm ( one of them containing the dot ) before reaching the collector .
it is furthermore reasonable to assume that the voltage drops at the point contacts and the dot remains in equilibrium @xcite . under these conditions
the magnetic flux @xmath42 dependent part of the conductance of the aharonov - bohm interferometer @xmath43 is given by @xcite @xmath44 } \end{aligned}\ ] ] with the flux quantum @xmath45 , the derivative @xmath46 of the fermi function as well as the transmission amplitudes @xmath47 of the arm containing the dot and @xmath48 of the reference arm .
the transmission @xmath49 is the product of the transmission @xmath50 through the dot and the transmission @xmath51 through the rest of the interferometer arm containing the dot .
it is reasonable to assume that @xmath51 as well as @xmath52 are only weakly energy and gate voltage dependent and thus the @xmath10-dependence of the aharonov - bohm oscillations of the measured linear conductance of the interferometer is dominated by the @xmath10 dependence of the magnitude and phase of the transmission amplitude through the dot . as usual @xcite , we compute the energy - averaged transmission phase @xmath1 and magnitude @xmath8 of the dot for a fixed spin direction as the phase and absolute value of @xmath53 where the frequency integral represents an energy average , weighted by the derivative of the fermi function , as appropriate for a finite electron temperature in the leads . in the limit @xmath54 , @xmath55 reduces to a @xmath4-function and @xmath56 is equal to @xmath57 .
we here take the chemical potential @xmath58 .    using scattering theory @xmath28 ( for fixed spin direction )
can be related to the spin - independent matrix elements ( in the @xmath23 indices of the wannier states ) of the dot s one - particle retarded green function @xmath59 @xmath60 with ( after taking the wide band limit ; see below ) @xmath61 where @xmath62 denotes the local density of states at the end of each semi - infinite lead . without loss of generality
we have assumed that @xmath63 , @xmath64 and @xmath65 .
the spin - independent dot occupancies @xmath22 ( per spin direction ) , that we will also investigate , follow from the green function @xmath66 by integrating over frequency ( or can be computed directly when using nrg ) . for l - r symmetry of the level - lead couplings and at temperature @xmath21 , @xmath67 can also be expressed in terms of the spin independent occupancies @xcite @xmath68\pi ) } e^{i(n_e+n_o)/\pi } \ ; ,    \label{t_lr}\ ] ] where @xmath69 for @xmath70 and @xmath71 for @xmath72 , respectively . here
we will compute @xmath59 in two ways , using both a truncated , that is approximate , frg scheme , and a numerically exact method , the nrg .
the truncated frg is an approximation scheme to obtain the self - energy @xmath73 ( and thus the one - particle green function ) and higher order vertex functions for many - body problems @xcite . as a first step in the application of this approach to quantum dots
one integrates out the noninteracting leads within the functional integral representation of our many - body problem @xcite .
the leads provide a frequency dependent one - particle potential on the dot levels . on the imaginary frequency axis it is given by @xmath74 where @xmath75 denotes the spin - independent green function of the isolated semi - infinite leads taken at the last lattice site @xmath76
as we are not interested in band effects we take the wide band limit .
the potential then reduces to @xmath77 after this step , instead of dealing with an infinite system we only have to consider the dot of two interacting levels .    in the computation of the interacting one - particle green function projected onto the dot system
the sum of the dot hamiltonian with @xmath78 and @xmath79 can be interpreted as a frequency dependent `` single - particle hamiltonian '' and in the following will be denoted by @xmath80 . for the spin - polarized case
it is a @xmath81 matrix in the quantum number @xmath23 . including spin , because of the additional quantum number @xmath82 , @xmath83 is a @xmath84 matrix which is block - diagonal in @xmath34 ( spin conservation ) .
as we are here not interested in the role of a magnetic field lifting the spin - degeneracy of each level the @xmath85 and @xmath86 blocks are equivalent .
the resolvent @xmath87^{-1}$ ] obtained from @xmath88 is equivalent to the noninteracting propagator of our two - level many - body problem projected on the dot levels . in the generating functional of the one - particle irreducible vertex functions we replace @xmath89 by @xmath90^{-1}\ ] ] with @xmath91 being an infrared cutoff running from @xmath92 down to @xmath93 .
taking the derivative with respect to @xmath91 one can derive an exact , infinite hierarchy of coupled differential equations for vertex functions , such as the self - energy and the one - particle irreducible two - particle interaction . in particular , the flow of the self - energy @xmath94 ( one - particle vertex ) is determined by @xmath94 and the two - particle vertex @xmath95 , while the flow of @xmath95 is determined by @xmath94 , @xmath95 , and the flowing three - particle vertex . the latter could be computed from a flow equation involving the four - particle vertex , and so on . at the end of the frg flow
@xmath96 is the self - energy @xmath73 of the original , cutoff - free problem we are interested in @xcite from which the green function @xmath59 can be computed using the dyson equation .
a detailed derivation of the frg flow equations for a general quantum many - body problem that only requires a basic knowledge of the functional integral approach to many - particle physics @xcite and the application of the method for a simple toy problem is presented in ref .
@xcite . for an overview of the application to quantum dots see refs .
@xcite and @xcite .
we here truncate the infinite hierarchy of flow equations by only keeping the self - energy and the frequency - independent part of the two - particle vertex .
higher order terms can be neglected if the bare two - particle interaction is not too large . by comparison to nrg data this approximation scheme was earlier shown to provide excellent results for a variety of dot systems @xcite . for further comparison
see fig .
[ fig7 ] below .
the present scheme leads to a frequency - independent self - energy ( see below ) . as finite frequency effects ( inelastic processes ) become important at temperatures @xmath18 , but these are not accurately treated by the level of approximation used here , in the present paper we shall show frg results only for @xmath21 .
the truncation leads to the coupled differential flow equations @xmath97 and @xmath98 \big\ } \label{finalflowgamma}\end{aligned}\ ] ] where @xmath99 , @xmath100 , etc .
are multi - indices representing the quantum numbers @xmath101 and @xmath102^{-1 } \ ; .\ ] ] in the model with spin - degenerate levels each index @xmath99 , @xmath100 etc
.  can take four different values @xmath23 and @xmath103 which gives 16 equations for @xmath94 and 256 for the two - particle vertex . for a spin - polarized two - level dot
the multi - indices take two values and one obtains 4 equations for @xmath94 and 16 for the two - particle vertex .
the number of independent equations can be significantly reduced ( see below ) taking into account the antisymmetry of the two - particle vertex and the spin symmetry ( for spin - degenerate levels ) both being preserved by eqs .
( [ finalflowsigma ] ) and ( [ finalflowgamma ] ) .
the initial conditions at @xmath104 are given by @xmath105 while @xmath106 is given by the bare antisymmetrized two - body interaction .
in the spin - polarized case the only nonzero components of the two - particle vertex at @xmath107 are @xmath108 in the model including spin the initial conditions take the form @xmath109 all other components which do not arise out of these by permutations ( @xmath110 and @xmath111 ) are zero . the self - energy matrix and thus the one - particle green function is completely independent of the spin direction and in the following we suppress the spin indices .
as already mentioned the present approximation leads to a frequency - independent self - energy .
this allows for a simple single - particle interpretation of its matrix elements .
the sum of the @xmath112 and the bare level position correspond to the flowing effective level positions , @xmath113 , while @xmath114 is a hopping between the levels 1 and 2 generated in the frg flow .
the frg formalism then reduces to a set of coupled differential flow equations for @xmath115 , @xmath116 and a few ( one in the spin - polarized case and seven for spin - degenerate levels ) independent components of the two - particle vertex .
these flow equations can easily be integrated numerically using standard routines .
it is important to note that although we start out with intra- and inter - level coulomb interactions of equal strengths they generically become different during the frg flow ( because of the different @xmath117 ) .
furthermore , additional interaction terms which are initially zero will be generated in the flow .
the set of equations significantly simplifies if the flow of the vertex is neglected while the results remain qualitatively the same . within this additional approximation and for a spin - polarized dot the flow equations for @xmath115 and @xmath116 are explicitly given in ref .  @xcite . in certain limiting cases
it is even possible to analytically solve the differential equations @xcite .
however , in the present work , the flow of the vertex is retained which clearly improves the quality of approximation @xcite .    at the end of the frg flow , the full green function takes the form @xmath118^{-1}_{j , j ' } = i \omega \delta_{j , j ' } - h_{j , j'}(i\omega)$ ] with an effective , _ noninteracting _ ( but @xmath10- , @xmath2- and @xmath119-dependent ) `` hamiltonian '' @xmath120 in a last step we have to perform the analytic continuation to the real frequency axis @xmath121 . this is straightforward , as the only frequency dependence of @xmath122 is the trivial one of the lead contribution eq .
( [ leadpotdefwideband ] ) .
then @xmath28 can be computed using eq .
( [ transformula ] ) .
the numerical renormalization group was invented by k.  g.  wilson in 1974 as a nonperturbative renormalization scheme for the kondo model @xcite .
it was later extended to the fermionic @xcite anderson model which describes a localized electronic state coupled to a fermionic bath .
the nrg allows thermodynamic and dynamic properties of such strongly correlated systems to be calculated at zero and finite temperature @xcite .
the key idea of nrg is to discretize the conduction band of the bath logarithmically , leading to a tight - binding chain for which the hopping matrix elements between the successive sites fall off exponentially with @xmath123 , where @xmath124 is the discretization parameter , typically @xmath125 , and @xmath126 is the site index .
this energy scale separation ensures that the problem can be solved iteratively by adding one site at a time and diagonalizing the enlarged system at each step , thereby resolving successively smaller and smaller energy scales .
thus , by choosing the length @xmath127 of the chain so large that the corresponding energy scale @xmath128 is smaller than all other energies in the problem , all relevant energy scales can be resolved and treated properly .
since the dimension of the hilbert space of the chain increases exponentially with the length of the chain , a truncation scheme has to be adopted , according to which only the lowest @xmath129 eigenstates of the chain are retained at each iteration .
recently , it was shown that by also keeping track of discarded states a complete , but approximate , basis of states can be constructed @xcite .
this can be used to calculate spectral functions which rigorously satisfy relevant sum rules @xcite .    in order to obtain the transmission through the dot @xmath130 eq .
( [ transformula ] ) we follow @xcite and @xcite to compute the imaginary part of the local green functions at temperature @xmath3 , using the lehmann representation @xmath131 )    \nonumber    \\ *    & -\pi &    \frac{e^{-\omega_n / t}}{z}\ ,      \langle n| d_{j\prime,\sigma}^\dag |m\rangle    \ ,    \langle m| d_{j,\sigma } |n\rangle    \ ,    \mathcal{\delta}(\omega+[\omega_m-\omega_n]),\end{aligned}\ ] ] with @xmath132 , the many - body eigenstates @xmath133 and eigenenergies @xmath134 .
since these are causal functions , the real part can be accessed by performing a kramers - kronig transformation @xcite .    in the next three sections we present our results . in sec .
[ resultsnon ] for the @xmath5 case . in sec .
[ resultsgeneric ] we present the generic phase lapse scenario for interacting spin - polarized dots , compare to the mean - field results and investigate the role of finite temperatures . finally , in sec .
[ resultsspin ] we study the spinful two - level dot .
the large number of parameters makes it essential to analyze the transmission for the noninteracting case before considering the effect of two - particle interactions .
we focus on @xmath21 .
a closed expression for @xmath13 and @xmath135 ( for a fixed spin direction ) at @xmath5 can be obtained from eqs .
( [ transamp ] ) and ( [ transformula ] ) by replacing @xmath136 by @xmath137 , @xmath138^{1/2}}{\left [ \left ( \gamma_1^l \gamma_2^r + \gamma_2^l \gamma_1^r   - 2 s \sqrt{\gamma_1^l \gamma_1^r   \gamma_2^l \gamma_2^r}- \varepsilon_{1 } \varepsilon_{2}\right)^2 + \left (   \varepsilon_{1 } \gamma_2 + \varepsilon_{2 } \gamma_1 \right)^2\right]^{1/2 } } \ ; , \\
\label{u0phase } \mbox { } \hspace{-1.5 cm } \theta(v_g ) & = & \arctan { \left [ \frac{\varepsilon_1\gamma_2+\varepsilon_2\gamma_1}{\varepsilon_1\varepsilon _ 2-\left(\sqrt{\gamma_1^l\gamma_2^r}-s\sqrt{\gamma_1^r\gamma_2^l}\right)^2 } \right ] } \ , \mbox{mod } \ , \pi \ ; , \end{aligned}\ ] ] with @xmath139 . for later use
we define @xmath140 for a fixed set of @xmath117 the @xmath4 dependence of @xmath13 and @xmath135 is shown in the first columns of figs.[fig1 ] and [ fig2 ] for @xmath141 .
the results are qualitatively the same as those obtained for @xmath5 . for generic level - lead couplings
@xmath117 the gate voltage dependence of eq .
( [ u0conduct ] ) in the limit of small and large @xmath142 is dominated by two peaks ( of height @xmath143 ) and a transmission zero .
associated with the transmission zero is a @xmath0 phase lapse at the same gate voltage .
the transmission zero ( and phase lapse ) follows from perfect destructive interference at a particular @xmath10 and in the limit of a strong asymmetry in the coupling of the two levels to the leads , @xmath144 or vice versa , can be understood as resulting from a fano anti - resonance @xcite .
apparently the fano anti - resonance with vanishing transmission is robust if one goes away from this limit towards more symmetric level - lead couplings . across each of the transmission resonances
@xmath1 increases roughly by @xmath0 as expected for a breit - wigner resonance .
further details of @xmath13 and @xmath135 depend on @xmath145 .
for @xmath70 the transmission zero ( and phase lapse ) is located between the two conductance peaks for all @xmath4 .
for @xmath146 the resonance peak positions depend on the asymmetry of the @xmath117 and the separation of the peaks is small if the @xmath117 are close to l - r symmetry , that is close to @xmath147 . for l - r symmetric dots and @xmath148 the transmission zero ( and phase lapse )
disappears ( not shown in the figures ) .
this is an example of a submanifold in parameter space with nongeneric behavior .
a complete account of such cases ( which also remain nongeneric for @xmath6 ) is given in refs .
@xcite and @xcite .
as they require fine tuning these parameter sets are presumably irrelevant in connection with the experiments and we will here only briefly mention results obtained in such cases .    for @xmath72 and fixed @xmath117 the position of the transmission zeros and phase lapses with respect to the cb peaks is different for small or large @xmath142 ( see fig .  [ fig2 ] ) . at small @xmath142
it is located between the two conductance peaks , whereas for large @xmath142 it lies on one of the outer sides of these peaks @xcite . in the crossover regime between these limiting cases
the height of one of the peaks decreases , while the other becomes broader and splits up into two resonances separated by a minimum with non - vanishing conductance [ see fig .  [ fig2 ] ( g ) and ( j ) ] .
the crossover scale @xmath149 depends on the choice of @xmath117 . for large @xmath142 , @xmath8 has three local maxima , although the height of one of the maxima is significantly smaller than the height of the other two ( not shown in fig .
[ fig2 ] ) . for fixed , asymmetric @xmath117 and @xmath146
the separation of the two conductance peaks for @xmath72 is significantly larger than for @xmath70 [ compare figs .
[ fig1 ] ( a ) and [ fig2 ] ( a ) ] .
it is important to note that for small @xmath142 essential features of the universal phase lapse regime established in the experiments are already found at @xmath5 : regardless of the sign @xmath145 for generic @xmath117 ( that is with the exception of a few cases with increased symmetry ) two transmission resonances are separated by a transmission zero and @xmath0 phase lapse . at @xmath5
the peak separation is too small and the shape of the @xmath10 dependence of the transmission and phase close to the peaks is qualitatively different from those observed experimentally ( namely lorentzian - like for the magnitude of the transmission , s - shaped for the phase ) .
as we show next the latter problems do not arise for sufficiently large interaction @xmath2 , which in particular leads to an increased separation @xmath150 of the transmission peaks .
in ref .  @xcite it was shown that frg and nrg results for @xmath13 and @xmath135 agree quantitatively up to fairly large @xmath2 . for a generic set of couplings @xmath151 we present frg data for @xmath135 and @xmath13 together with the occupancies of the levels @xmath22 for different @xmath2 and @xmath4 in figs .
[ fig1 ] ( @xmath70 ) and [ fig2 ] ( @xmath72 ) . increasing @xmath152 the separation of the transmission peaks in the limit of small and large @xmath142 increases and
is eventually given by @xmath150 .
even though this charging effect appears to be straightforward it is important to note that in particular the groundstate at small @xmath142 is highly correlated .
this becomes explicit from the mapping of the present problem on a generalized single impurity anderson and kondo model as discussed in ref .
an indication of strong correlation effects are the correlation - induced resonances of the transmission found in ref .
@xcite , which we briefly mention below . with increasing @xmath2 , even at small
@xmath142 the gate voltage dependence of @xmath135 across the transmission resonances becomes s - shaped and the resonances more lorentzian - like [ see third columns of figs .  [ fig1 ] and [ fig2 ] ] .
obviously , for @xmath70 the transmission zero and phase lapse remain between the two transmission peaks for all @xmath4 and @xmath2 ( see fig .
[ fig1 ] ) . for @xmath72
this only holds for sufficiently small level spacings as , similar to the @xmath5 case , with increasing @xmath142 a crossover sets in to a regime in which the transmission zero and phase lapse are no longer between the peaks .
analogously to the @xmath5 case , the crossover scale @xmath149 depends on the particular choice of @xmath117 . as can be seen from the second row of fig .
[ fig2 ] ( the cb peaks at large @xmath2 have still almost equal height ) , with increasing @xmath152 , @xmath149 is pushed towards larger values .
the coulomb interaction thus stabilizes the parameter regime of universal phase lapses .
this shows that the effect of the coulomb interaction leading to universal @xmath0 phase lapses between separated cb peaks in a two - level dot is rather straightforward : for small @xmath142 the phase lapse and transmission zero are already present at @xmath5 , and the effect of finite @xmath2 is simply that the cb peaks become well separated because of charging effects .
they also lead to a lorentzian - like lineshape of the peaks and an s - like variation of @xmath1 across them.the present scenario has to be contrasted to the one obtained for @xmath20 levels discussed in ref .  @xcite .
the generic appearance of @xmath153 transmission zeros and phase lapses separating the transmission peaks at small @xmath142 and @xmath5 is specific to the case with @xmath19 levels . for @xmath20
the number of transmission zeros and phase lapses at @xmath5 strongly depends on the parameters and the mechanism leading to universal @xmath0 phase lapses at sufficiently large @xmath2 ( at small @xmath142 ) is much more involved @xcite .
this shows that although important insights can be gained from studying two - levels to achieve a _
complete _ understanding of the phase lapse scenario it is essential to study dots with coulomb interaction and more than two levels @xcite .     and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 .
the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=35.2%,height=141 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=29.2%,height=141 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=33.0%,height=141 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=35.2%,height=117 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=29.2%,height=117 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=33.0%,height=117 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=35.2%,height=117 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=29.2%,height=117 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 .
the parameters are @xmath154 and @xmath70 . for better visibility @xmath155 were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=33.3%,height=117 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 .
the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=35.2%,height=117 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=29.2%,height=117 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=33.0%,height=117 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=35.2%,height=137 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 .
the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=29.2%,height=137 ]   and @xmath142 that govern the gate voltage @xmath10 dependence of the magnitude of the transmission @xmath8 ( black ) , the transmission phase @xmath1 ( red ) and the level occupancies ( green and blue ) of a spin - polarized two - level dot at @xmath21 . the parameters are @xmath154 and @xmath70 . for better visibility @xmath155
were shifted by @xmath156 .
the depicted behavior is the generic one and in particular qualitatively independent of the actual choice of @xmath151 ( up to certain cases of increased symmetry ; for examples see the text ) .
the behavior at @xmath141 is qualitatively the same as the one at @xmath5 .
the results were obtained using the truncated frg.,title="fig:",scaledwidth=33.0%,height=137 ]    , but for @xmath72.,title="fig:",scaledwidth=35.2%,height=141 ] , but for @xmath72.,title="fig:",scaledwidth=29.2%,height=141 ] , but for @xmath72.,title="fig:",scaledwidth=33.0%,height=141 ] , but for @xmath72.,title="fig:",scaledwidth=35.2%,height=117 ] , but for @xmath72.,title="fig:",scaledwidth=29.2%,height=117 ] , but for @xmath72.,title="fig:",scaledwidth=33.0%,height=117 ] , but for @xmath72.,title="fig:",scaledwidth=35.2%,height=117 ] , but for @xmath72.,title="fig:",scaledwidth=29.2%,height=117 ] , but for @xmath72.,title="fig:",scaledwidth=33.3%,height=117 ] , but for @xmath72.,title="fig:",scaledwidth=35.2%,height=117 ] , but for @xmath72.,title="fig:",scaledwidth=29.2%,height=117 ] , but for @xmath72.,title="fig:",scaledwidth=33.0%,height=117 ] , but for @xmath72.,title="fig:",scaledwidth=35.2%,height=137 ] , but for @xmath72.,title="fig:",scaledwidth=29.2%,height=137 ] , but for @xmath72.,title="fig:",scaledwidth=33.0%,height=137 ]    the lineshape of @xmath8 shows characteristic differences in the limits of small and large @xmath142 . in the universal regime at small @xmath142 and for sufficiently large
@xmath152 the two cb peaks have equal width of order @xmath157 ( not @xmath158 ) and equal height which is consistent with the expectation that at each peak the transport occurs through both bare levels simultaneously .
a similar behavior is observed in the experiments . in the mesoscopic regime ( @xmath159 )
the width of the @xmath160th peak is given by @xmath158 and the relative height @xmath161 by @xmath162 , independent of the value of @xmath2 .
this effectively noninteracting lineshape can be understood from the gate voltage dependence of the effective level positions @xmath163 at the end of the frg flow .
when one level is charged the effective level position of the other level is pushed downwards by @xmath2 . besides this , the gate voltage dependence of the @xmath163 remains linear , leading to two transmission peaks at gate voltages @xmath164 with separation @xmath150 , but with the same width and height as for @xmath5 .
the hopping between the two effective levels generated in the frg flow is small and can be neglected .
apart from the cb peaks , for sufficiently large @xmath152 the transmission shows additional features at small @xmath165 [ see figs .
[ fig1 ] ( b ) , ( c ) , ( e ) and [ fig2 ] ( c ) ] .
these are the correlation - induced resonances mentioned in the introduction , which have been found to be most pronounced at @xmath148 and in this case occur for interactions larger than a critical @xmath166 which depends on the @xmath117 and @xmath145 @xcite .
their appearance indicates that the groundstate at small @xmath142 is strongly correlated ( as mentioned in ref .
@xcite the correlation - induced resonances are not captured by a mean - field analysis ; see below ) . associated with the correlation - induced resonances
is a sharp increase of @xmath1 [ see figs .
[ fig1 ] ( b ) and ( e ) ] . at large @xmath152
the correlation - induced resonances are exponentially ( in @xmath152 ) sharp features that vanish quickly with increasing @xmath3 ( see below ) , which might be one of the reasons why up to now they have not been observed in experiments .
the correlation - induced resonances are not directly linked to the universal phase lapse scenario .    for increasing @xmath152 at fixed @xmath142 and decreasing @xmath142 at fixed @xmath152
we observe an increased tendency towards population inversion of the @xmath22 .
we define that a population inversion occurs if ( 1 ) @xmath167 at a certain @xmath168 and ( 2 ) one @xmath22 has positive and one negative slope at @xmath168 so that the filling of one level causes a tendency for the other to empty . for large @xmath152 and small @xmath142
it is mainly the more strongly coupled level ( in figs .
[ fig1 ] and [ fig2 ] , this is the level 2 shown in blue ) whose population increases across both cb peaks while it is depopulated in between .
this behavior is reminiscent of the one discussed in the model with a broad and several narrow levels @xcite , where a relation between population inversion and phase lapse behavior was proposed .
remarkably , for sufficiently large @xmath152 , we find population inversion even for small asymmetries @xmath169 ( which is only 1.5 in the example of figs.[fig1 ] and [ fig2 ] ) .
we emphasize that despite this resemblance to the observation of ref .
@xcite , the @xmath19 model is not appropriate to establish a general relation between the appearance of population inversions and @xmath0 phase lapses at small @xmath142 @xcite . while the latter are already present at @xmath5 , the former only develop with increasing @xmath2 [ compare figs .  [ fig1 ] ( a ) , ( d ) or [ fig2 ] ( a ) , ( d ) to figs
.  [ fig1 ] ( c ) , ( f ) or [ fig2 ] ( c),(f ) ] .
note that the gate voltage @xmath168 at which the population inversion occurs is generically _ not _ identical to the position of the phase lapse and transmission zero [ see figs .  [ fig1 ] ( b ) , ( c ) , ( f ) and [ fig2 ] ( c ) , ( f ) ] @xcite .
however , for l - r symmetric @xmath117 eq .  ( [ t_lr ] ) ensures that if a population inversion occurs its position is identical to the one of the phase lapse and transmission zero .    as can be seen in figs .
[ fig1 ] ( c ) , ( f ) and [ fig2 ] ( c ) , ( f ) for small @xmath142 and large @xmath152 the @xmath22 show a rather strong gate voltage dependence between the cb peaks .
nevertheless the total dot occupancy @xmath170 is only weakly @xmath10 dependent and close to 1 within the entire cb valley .
this is reminiscent of the plateau - like occupancy in the local moment regime of the single impurity anderson model showing the kondo effect . as discussed in ref .
@xcite a relation to this model can indeed be established .
we note in passing that with the exception of the nongeneric case of l - r symmetric @xmath117 , @xmath70 and @xmath148 , the @xmath22 are continuous functions of @xmath10 .      in refs .
@xcite and @xcite golosov and gefen ( gg ) analyze the phase lapse scenario of the spin - polarized interacting two - level dot within the mean - field approximation .
however , they anticipated themselves that quantum fluctuations not captured in the mean - field approach could be important .
examples of this had been pointed out already in refs .
@xcite and @xcite .
thus , gg emphasized that the effects of such fluctuations on their results need to be studied in subsequent work .
the present subsection is devoted to this task .     and different @xmath171 moving from ggs `` phase '' 2 into `` phase '' 1 . figures ( e)-(h ) are for constant @xmath171 and different @xmath172 moving from ggs `` phase 3 '' into `` phase 1 '' .
the inset in fig .
( e ) shows a zoom in of the gate voltage region around the phase lapse . for a detailed comparison to the mean - field results
see the text .
the results were obtained at @xmath21 using the truncated frg.,title="fig:",scaledwidth=45.5%,height=143 ]   and different @xmath171 moving from ggs `` phase '' 2 into `` phase '' 1 .
figures ( e)-(h ) are for constant @xmath171 and different @xmath172 moving from ggs `` phase 3 '' into `` phase 1 '' .
the inset in fig .
( e ) shows a zoom in of the gate voltage region around the phase lapse . for a detailed comparison to the mean - field results
see the text .
the results were obtained at @xmath21 using the truncated frg.,title="fig:",scaledwidth=43.5%,height=143 ]   and different @xmath171 moving from ggs `` phase '' 2 into `` phase '' 1 .
figures ( e)-(h ) are for constant @xmath171 and different @xmath172 moving from ggs `` phase 3 '' into `` phase 1 '' .
the inset in fig .
( e ) shows a zoom in of the gate voltage region around the phase lapse . for a detailed comparison to the mean - field results
see the text .
the results were obtained at @xmath21 using the truncated frg.,title="fig:",scaledwidth=45.5%,height=173 ]   and different @xmath171 moving from ggs `` phase '' 2 into `` phase '' 1 .
figures ( e)-(h ) are for constant @xmath171 and different @xmath172 moving from ggs `` phase 3 '' into `` phase 1 '' .
the inset in fig .
( e ) shows a zoom in of the gate voltage region around the phase lapse . for a detailed comparison to the mean - field results
see the text .
the results were obtained at @xmath21 using the truncated frg.,title="fig:",scaledwidth=43.5%,height=173 ]   and different @xmath171 moving from ggs `` phase '' 2 into `` phase '' 1 .
figures ( e)-(h ) are for constant @xmath171 and different @xmath172 moving from ggs `` phase 3 '' into `` phase 1 '' .
the inset in fig .
( e ) shows a zoom in of the gate voltage region around the phase lapse . for a detailed comparison to the mean - field results
see the text .
the results were obtained at @xmath21 using the truncated frg.,title="fig:",scaledwidth=45.5%,height=143 ]   and different @xmath171 moving from ggs `` phase '' 2 into `` phase '' 1 .
figures ( e)-(h ) are for constant @xmath171 and different @xmath172 moving from ggs `` phase 3 '' into `` phase 1 '' .
the inset in fig .
( e ) shows a zoom in of the gate voltage region around the phase lapse . for a detailed comparison to the mean - field results
see the text .
the results were obtained at @xmath21 using the truncated frg.,title="fig:",scaledwidth=43.5%,height=143 ]   and different @xmath171 moving from ggs `` phase '' 2 into `` phase '' 1 .
figures ( e)-(h ) are for constant @xmath171 and different @xmath172 moving from ggs `` phase 3 '' into `` phase 1 '' .
the inset in fig .
( e ) shows a zoom in of the gate voltage region around the phase lapse . for a detailed comparison to the mean - field results
see the text .
the results were obtained at @xmath21 using the truncated frg.,title="fig:",scaledwidth=45.5%,height=173 ]   and different @xmath171 moving from ggs `` phase '' 2 into `` phase '' 1 .
figures ( e)-(h ) are for constant @xmath171 and different @xmath172 moving from ggs `` phase 3 '' into `` phase 1 '' .
the inset in fig .
( e ) shows a zoom in of the gate voltage region around the phase lapse . for a detailed comparison to the mean - field results
see the text .
the results were obtained at @xmath21 using the truncated frg.,title="fig:",scaledwidth=43.5%,height=173 ]    gg consider the subspace of level - lead couplings defined by @xmath173 . performing a unitary transformation on the dot states the part of the hamiltonian eq .
( [ generalham ] ) containing dot operators can be transformed to ( see refs .
@xcite and @xcite ) @xmath174 \;,\end{aligned}\ ] ] with the transformed operators and parameters indicated by a hat .
the change of basis leads to a direct hopping @xmath175 between the transformed levels .
we here focus on the relative sign @xmath176 .
gg then introduce the two new dimensionless parameters @xmath172 and @xmath171 as @xmath177 and @xmath178 .
varying @xmath172 and @xmath171 gg investigate the phase lapse behavior for fixed level spacing ( in the new basis ) @xmath179 and interaction @xmath180 .
as the parameters in the new basis are rather complicated combinations of the original ones , the variation of @xmath171 and @xmath172 corresponds to the variation of the @xmath117 ( within the above specified subspace ) , @xmath4 and even @xmath145 .
our simple picture that increasing the level spacing @xmath4 of the untransformed model leads from the universal to the mesoscopic regime can not easily made explicit using the parameters of gg . to make the comparison of our results to the mean - field study definite we nevertheless follow the steps of gg .
varying @xmath181 at @xmath176 we move around in the right part of what is called the `` phase diagram '' by gg ( fig .  4 of ref .
@xcite ) . in fig .
[ fig3 ] we show the behavior of @xmath8 and @xmath1 varying @xmath171 at fixed @xmath172 [ figs .  [ fig3 ] ( a)-(d ) ] and @xmath172 at fixed @xmath171 [ figs.[fig3 ] ( e)-(h ) ] , respectively . in the first case , upon increasing @xmath171 at constant @xmath182
, we move from ggs `` phase 2 '' , ( red in fig .  4 of ref .
@xcite ) , with the transmission zero and phase lapse outside the two cb peaks , into `` phase 1 '' ( blue in fig .  4 of ref .
@xcite ) , with the transmission zero and phase lapse between the peaks .
the mean - field approximation correctly captures the presence of these two regimes .
we find a smooth crossover between them ( which is why we prefer the notion of different `` regimes '' rather than `` phases '' ) : the @xmath10 value at which the transmission zero and phase lapse occur _ smoothly _ crosses from lying outside the right cb peak to lying between the two cb peaks .
this is similar to the smooth crossover we observe in the @xmath72 case of the untransformed model when moving from the mesoscopic to the universal regime ( see columns of fig .
[ fig2 ] ) .    for fixed @xmath183 and increasing @xmath172 we move from `` phase 3 '' ( green in fig .  4 of ref .
@xcite ) into `` phase 1 '' .
in contrast to the mean - field approximation where an abrupt transition from `` phase 3 '' to `` phase 1 '' occurs , figs .
[ fig3 ] ( e)-(h ) show a rather smooth evolution .
the mean - field `` phase 3 '' is characterized by discontinuous population switching and a phase lapse between the cb peaks which , surprisingly , is smaller than @xmath0 @xcite .
furthermore , in this parameter regime the mean - field results show no transmission zero , as discussed in ref .
however , from figs .
[ fig3 ] ( e)-(h ) it is apparent that `` phase 3 '' , is an artifact of the mean - field approximation : upon taking fluctuations into account via frg , the @xmath0 phase lapse and the transmission zero are found to remain in the cb valley .
the evolution with increasing @xmath172 is similar to the @xmath70 case of the untransformed model when the level spacing is increased at fixed @xmath117 ( see the second and third columns of fig .
[ fig1 ] ) .
`` phase 3 '' then corresponds to the parameter regime with small level spacing and a sizable @xmath152 in which correlations are of particular importance leading e.g.  to the correlation - induced resonances [ see figs .
[ fig3 ] ( e)-(h ) ]
. that the mean - field approximation fails to properly describe this strongly correlated regime is not surprising and has been recognized earlier @xcite . in ref .
@xcite a connection between the small @xmath4 regime of the present model and the local moment ( kondo ) regime of the single - impurity anderson model was established .
thus , the artifacts of the mean - field approximation for the spinless two - level dot are reincarnations of the the well - known artifacts it produces when applied to the anderson model in the local moment regime .    upon taking fluctuations into account ,
the discontinuities of the @xmath155 in `` phase 3 '' are washed out .
the only choice of parameters for which we find discontinuous behavior is the one with l - r symmetric @xmath117 , @xmath70 and @xmath148 , a case which was already identified as being nongeneric @xcite .
any arbitrarily small deviation from these conditions leads to a continuous gate voltage dependence of @xmath155 . for parameters close to the nongeneric point
the change in @xmath22 at first sight appears to be rather sharp and an extremely high resolution in @xmath10 is required to identify the behavior as continuous .
the fact that the mean - field treatment at @xmath6 and small level spacings incorrectly causes the transmission zero to disappear and the corresponding phase lapse to become smaller than @xmath0 is its most consequential problem . by generating such features ,
the mean - field treatment masks the very simple scenario that emerges upon properly including fluctuations : for increasing @xmath2 , the @xmath0 phase lapse and transmission zero found for small @xmath142 at @xmath5 remain in the cb valley , while the transmission peaks become well separated , lorentzian - like and the phase acquires an s - shape across the resonances .
we next investigate how the phase and magnitude of the transmission are affected by finite temperatures .
to this end , we use the new fdm - nrg algorithm recently proposed in ref .  @xcite . in fig .
[ fig4 ] we show nrg data for @xmath13 and @xmath135 at different @xmath3 .
we consider generic level - lead couplings @xmath184 , @xmath185 , @xmath186 and @xmath187 ( universal regime ) as well as @xmath188 ( mesoscopic regime ) .
as expected , with increasing temperature the sharp @xmath0 phase lapses are gradually smeared out , the transmission zero vanishes and the change of @xmath1 at the phase lapse becomes smaller than @xmath0 .
furthermore , the cb peaks decrease and broaden . for @xmath188 and @xmath72 [ fig .
[ fig4 ] ( d ) ] the phase lapse lies outside the cb peaks and outside the window of gate voltages shown .     and
@xmath135 obtained by nrg for @xmath189 , @xmath185 for @xmath186 and in the universal ( @xmath187 ) and the mesoscopic ( @xmath188 ) regime .
we used the nrg parameters @xmath190 and @xmath191.,title="fig:",scaledwidth=49.5%,height=166 ]   and @xmath135 obtained by nrg for @xmath189 , @xmath185 for @xmath186 and in the universal ( @xmath187 ) and the mesoscopic ( @xmath188 ) regime .
we used the nrg parameters @xmath190 and @xmath191.,title="fig:",scaledwidth=47.5%,height=166 ]   and @xmath135 obtained by nrg for @xmath189 , @xmath185 for @xmath186 and in the universal ( @xmath187 ) and the mesoscopic ( @xmath188 ) regime .
we used the nrg parameters @xmath190 and @xmath191.,title="fig:",scaledwidth=49.5%,height=196 ]   and @xmath135 obtained by nrg for @xmath189 , @xmath185 for @xmath186 and in the universal ( @xmath187 ) and the mesoscopic ( @xmath188 ) regime .
we used the nrg parameters @xmath190 and @xmath191.,title="fig:",scaledwidth=47.5%,height=196 ]    in the mesoscopic regime [ figs .
[ fig4 ] ( b ) and ( d ) ] the explicit temperature dependence of the green function entering eq .
( [ transamp ] ) via eq.([transformula ] ) is rather weak and the temperature dependence of the cb peaks and the phase lapse can be understood from the behavior in the noninteracting model , but with level spacing @xmath150 . for small
@xmath3 the height @xmath192 of the @xmath160th cb peak scales as @xmath193 and the width @xmath194 of the phase lapse as @xmath195 @xcite .
the relevant scale for sizable temperature effects in the peak height is thus @xmath158 while it is @xmath150 in the smearing of the phase lapses .
since we have chosen @xmath196 a reduction of @xmath161 is visible for temperatures at which the phase lapse is still fairly sharp [ see fig .
[ fig4 ] ( b ) ] .    due to the importance of correlation effects at small @xmath142 ,
the @xmath3 dependence of @xmath13 and @xmath135 in the universal regime is different from the noninteracting case . here
the explicit temperature dependence of @xmath197 is much stronger and can not be neglected .
the resulting @xmath3 dependence of @xmath13 and @xmath135 is shown in figs .
[ fig4 ] ( a ) and ( c ) .
a comparison to fig.[fig4 ] ( b ) shows that in the universal regime the smearing of the phase lapse sets in at a lower energy scale than in the mesoscopic regime .
this scale depends on the relative sign @xmath145 of the level - lead hopping matrix elements [ compare figs .
[ fig4 ] ( a ) and ( c ) ] . furthermore , in contrast to the mesoscopic regime the scales on which the cb peaks and the phase lapse are affected by temperature are comparable
. a more detailed investigation of the temperature dependence in the universal regime , which also discusses the fate of the correlation - induced resonances , is beyond the scope of the present work and is left as subject for future studies .    , but for @xmath198 and @xmath72.,title="fig:",scaledwidth=49.5%,height=166 ] , but for @xmath198 and @xmath72.,title="fig:",scaledwidth=47.5%,height=166 ] , but for @xmath198 and @xmath72.,title="fig:",scaledwidth=49.5%,height=196 ] , but for @xmath198 and @xmath72.,title="fig:",scaledwidth=47.5%,height=196 ]    , but for @xmath198 and @xmath72.,title="fig:",scaledwidth=49.5%,height=166 ] , but for @xmath198 and @xmath72.,title="fig:",scaledwidth=47.5%,height=166 ] , but for @xmath198 and @xmath72.,title="fig:",scaledwidth=49.5%,height=196 ] , but for @xmath198 and @xmath72.,title="fig:",scaledwidth=47.5%,height=196 ]
we finally investigate the effect of the spin degree of freedom on the discussed phase lapse scenario , at @xmath21 . in figs .
[ fig5 ] ( @xmath70 , @xmath199 ) and [ fig6 ] ( @xmath72 , @xmath198 ) we show frg data for the evolution of @xmath13 , @xmath135 and @xmath200 ( for a fixed spin direction ) with increasing @xmath4 for a generic @xmath201 . the overall dependence of @xmath135 on @xmath4 is similar to the one observed in the spinless case ( compare to figs .
[ fig1 ] and [ fig2 ] ) .
in particular , the behavior at small @xmath142 appears to be almost unaffected by the presence of the spin degree of freedom [ figs .
[ fig5 ] ( a ) , ( b ) and [ fig6 ] ( a ) , ( b ) ] . for large @xmath142 [ figs .
[ fig5 ] ( d ) and [ fig6 ] ( d ) ] the transmission resonances are located at odd average total filling of the two - level dot indicated by shoulders in the @xmath22 [ figs .
[ fig5 ] ( d ) and [ fig6 ] ( d ) ] . at these fillings and for
sufficiently large @xmath152 the kondo effect is active and the resonances can not be regarded as lorentzian - like cb peaks . instead they show a plateau - like shape known from the spinful single - level dot ( see ref .
@xcite and references therein ) . across the kondo plateaus of @xmath8
the s - shaped increase of the phase is interrupted by a shoulder at @xmath202 as expected for the kondo effect @xcite .
it would be very interesting to study how each of the kondo plateaus of @xmath8 with increasing temperature crosses over to two cb peaks and the phase behaves in the generated cb valley .
this question is left for future investigations .
the behavior of the phase in the presence of the kondo effect was experimentally investigated at temperatures comparable to the kondo temperature @xcite , and much below the kondo temperature @xcite .
as we study the zero temperature case it is proper to compare our calculations to the measurements at low temperatures . indeed fig .
[ fig6 ] ( c ) for intermediate @xmath142 qualitatively reproduces the experimental results at low temperature as shown in fig .
3 c of ref.@xcite .
in particular the increase of the phase by more then @xmath0 and the absence of clearly developed kondo plateaus are reproduced .    in fig .
[ fig5 ] ( a ) we left out the frg data around @xmath203 as for these gate voltages some of the components of the flowing two - particle vertex become large .
this indicates the breakdown of our present truncation scheme @xcite and the results for @xmath8 , @xmath1 and @xmath22 become unreliable . for an explicit comparison to nrg data of @xmath8
see fig .
[ fig7 ] ( a ) .
we note in passing that for @xmath70 correlation - induced resonances occur also in the model with spin [ fig .
[ fig5 ] ( a ) ] @xcite .    in fig .
[ fig7 ] we compare frg and nrg results for l - r symmetric level - lead couplings .
the computational resources required to obtain nrg data away from l - r symmetry become large and such data are not required for the aim of the present paper .
we can then use eq .
( [ t_lr ] ) and must only compute the occupancies @xmath22 , which is numerically less demanding . for @xmath204 ( as exclusively shown ) it is only the absence of the correlation - induced resonances for @xmath70 [ compare figs .
[ fig5 ] ( a ) and [ fig7 ] ( a ) ] which is different from the results of the generic @xmath151 shown in figs .
[ fig5 ] and [ fig6 ] . with the exception of the @xmath205 regime in the case of small @xmath4 and
@xmath70 the frg and nrg data compare quite well . in this case
we only show the frg data for @xmath8 as the results for the phase and occupancies become rather erratic .
the reason for the breakdown of the currently used truncated frg is explained in refs .
@xcite and @xcite and is related to the fact that at small @xmath4 and small @xmath10 the correlations in effect become extremely large .    ,
@xmath135 and @xmath200 ( per spin direction ) in the universal ( small @xmath142 ) and mesoscopic regime ( large @xmath142 ) of a spinful two - level dot at @xmath21 .
the parameters are @xmath199 and @xmath206 . in ( a ) the frg data for @xmath135 and @xmath200 around @xmath203
are not shown .
already the unsatisfactory comparison between frg and nrg for @xmath8 indicate that the frg becomes unreliable in this regime . for more details on this , see the text . for the nrg parameters we used @xmath207 , @xmath208 for @xmath70 and @xmath209 for @xmath72.,title="fig:",scaledwidth=47.5%,height=166 ] , @xmath135 and @xmath200 ( per spin direction ) in the universal ( small @xmath142 ) and mesoscopic regime ( large @xmath142 ) of a spinful two - level dot at @xmath21 .
the parameters are @xmath199 and @xmath206 . in ( a ) the frg data for @xmath135 and
@xmath200 around @xmath203 are not shown .
already the unsatisfactory comparison between frg and nrg for @xmath8 indicate that the frg becomes unreliable in this regime . for more details on this , see the text . for the nrg parameters we used @xmath207 , @xmath208 for @xmath70 and @xmath209 for @xmath72.,title="fig:",scaledwidth=47.5%,height=166 ] , @xmath135 and @xmath200 ( per spin direction ) in the universal ( small @xmath142 ) and mesoscopic regime ( large @xmath142 ) of a spinful two - level dot at @xmath21 .
the parameters are @xmath199 and @xmath206 . in ( a )
the frg data for @xmath135 and @xmath200 around @xmath203 are not shown .
already the unsatisfactory comparison between frg and nrg for @xmath8 indicate that the frg becomes unreliable in this regime . for more details on this , see the text .
for the nrg parameters we used @xmath207 , @xmath208 for @xmath70 and @xmath209 for @xmath72.,title="fig:",scaledwidth=47.5%,height=196 ] , @xmath135 and @xmath200 ( per spin direction ) in the universal ( small @xmath142 ) and mesoscopic regime ( large @xmath142 ) of a spinful two - level dot at @xmath21 .
the parameters are @xmath199 and @xmath206 . in ( a ) the frg data for @xmath135 and @xmath200 around @xmath203
are not shown .
already the unsatisfactory comparison between frg and nrg for @xmath8 indicate that the frg becomes unreliable in this regime . for more details on this , see the text . for the nrg parameters we used @xmath207 , @xmath208 for @xmath70 and @xmath209 for @xmath72.,title="fig:",scaledwidth=47.5%,height=196 ]
in the present paper we studied the appearance of phase lapses in an interacting two - level quantum dot considering the entire parameter space using nrg and a truncated frg scheme .
as a starting point we briefly discussed the noninteracting case at temperature @xmath21 and pointed out that for generic level - lead couplings , that is up to cases with increased symmetry , essential features of the universal phase lapse scenario are already established at @xmath5 .
for single - particle level spacings @xmath4 small compared to the level broadenings @xmath158 the transmission is characterized by two transmission peaks of equal width with a transmission zero and an associated @xmath0 phase lapse between them ( universal regime ) . for large @xmath142 at @xmath5 the appearance or not of a transmission zero and phase lapse between the two transmission peaks depends on the relative sign @xmath145 of the level - lead couplings ( mesoscopic regime ) . within a spinless model we have shown that the separation of the two transmission peaks increases linearly with the interaction @xmath2 while the @xmath0 phase lapse and transmission zero remain in the valley between them . furthermore , with increasing @xmath2 the increase of the phase across the peaks takes an s - shape and the peaks become lorentzian - like , thus assuming shapes resembling those observed experimentally @xcite . for @xmath72 and increasing @xmath4 , a crossover occurs to a regime in which the @xmath0 phase lapse and transmission zero lies outside the two cb peaks . the crossover scale @xmath149 increases with increasing interaction and thus the coulomb repulsion stabilizes the universal phase lapse behavior .
we have investigated the relation between phase lapses and population inversions of the level occupancies @xmath22 .
we have shown that a mean - field treatment of the present problem correctly reproduces certain features of the behavior discussed above , but is not able to produce the universal phase lapse scenario at small @xmath142 due to artifacts of the approximation such as a phase lapse by less than @xmath0 @xcite and a vanishing of the transmission zero @xcite .
furthermore , the discontinuous gate voltage dependence of the @xmath22 found in the mean - field approximation turned out to be an artifact .
next , we studied how the phase lapse behavior is affect by temperatures @xmath18 . the universal phase lapse ( at small @xmath142 )
is smeared out but remains visible for not too large @xmath3 . in the mesoscopic regime with @xmath210 the smearing of the phase lapse and the decrease of the cb peaks can be understood in detail in analogy to the noninteracting case . for @xmath211 correlations
are more important and a detailed understanding of the temperature dependence of the transmission @xmath56 requires further studies .    the phase lapse behavior in both the universal and mesoscopic regimes is also stable if the spin degree of freedom is included . for sufficiently large @xmath152 in this case
the kondo effect is active at odd average dot filling , leading to minor modifications of the scenario discussed above .
in particular , at large @xmath142 the @xmath21 transmission peaks are kondo plateaus rather than lorentzian - like cb peaks . across these kondo
plateaus the phase shows a shoulder at @xmath202 .
in contrast , the behavior at small @xmath142 appears to be almost unaffected by the spin degree .
a study of the combined effect of finite temperature and spin is left for future work .
we thank r.  berkovits , p.  brouwer , y.  gefen , l. glazman , d.  golosov , m.  heiblum , j.  imry , v.  kashcheyevs , j.  knig , f.  marquardt , m.  pustilnik , h.  schoeller , k.  schnhammer , and a.  silva for valuable discussions . this work was supported in part , for vm , by the dfg through sfb602 ; for th , aw and jvd by the dfg through sfb631 and de730/3 - 2 , and also by spintronics rtn ( hprn - ct-2002 - 00302 ) , nsf ( phy99 - 07949 ) and dip - h.2.1 ; and for yo by dip - h.2.1 , bsf and the humboldt foundation .