in a small metallic island where the charging energy @xmath5 exceeds the temperature @xmath6 ( we use the unit @xmath7 ) , the coulomb interaction affects transport properties through the island .
the resulting phenomenon is called the coulomb blockade ( cb ) and such system is named the single - electron transistor ( set ) .
cb has been attracted much attention in the last decade @xcite and the nature of the transport properties of set has been clarified .
set is interesting because it is regarded as one of the most simple examples of strongly correlated system which can be brought into the non - equilibrium state by applied bias voltage .
early investigations consider the case where the tunneling conductance is so small that the higher order quantum fluctuation of island charge is negligible .
recently , the quantum fluctuation in set has attracted much attention as one of the basic problems in this field .
the quantum fluctuation is quantitatively characterized by the dimensionless parallel conductance : @xmath8 where @xmath9 is the quantum resistance and @xmath10 is the parallel tunneling resistance of the source and the drain junctions .
there has been much development on theoretical investigation in the whole range of @xmath11 . especially in the _ weak tunneling regime _
( @xmath12 ) , the life - time broadening of a charge state level is much smaller than the typical level spacing of charge states , and thus the effective two - state model , which is equivalent to the multichannel anisotropic kondo model in the equilibrium state@xcite , well describes the low - energy physics . with this model
, it is predicted that the quantum fluctuation of charge causes the renormalization of the conductance and the charging energy below the kondo temperature " @xmath13 @xcite .
the renormalization of the conductance is confirmed experimentally as @xmath14 dependence of the conductance peak at low temperature@xcite .
it is also predicted that in the non - equilibrium state , the dissipative charge fluctuation causes the life - time broadening of a charge state level and smears the structures of @xmath15-@xmath16 characteristic@xcite .
though investigations have revealed much about the quantum fluctuation , most of them have been limited to averaged quantities . in order to understand the nature of the quantum fluctuation ,
investigations on the higher - order correlation function of fluctuation operators are required .
a good starting point may be the investigation on the second moment of fluctuation operators , i.e. , the noise@xcite .
the charge noise and the current noise in the weak tunneling regime is also important for practical applications , because it determines the performance of set electrometers @xcite .
the current noise is defined by the auto - correlation function of the current fluctuation operator @xmath17 as @xmath18 where @xmath19 means the statistical average . until recently
, investigations on the noise have been done using the framework of the orthodox `` theory@xcite , which takes account of the lowest order quantum fluctuation , namely the sequential tunneling ( st ) process .
recently , several authors@xcite discussed the higher order quantum fluctuation in cb regime within the co - tunneling theory''@xcite . however , there is no approximation covering both st and cb regime .
the aim of the present work is to construct a theoretical framework , which covers both of these regimes for arbitrary @xmath11 and clarify how the quantum fluctuation affects the noise .
the keldysh formalism@xcite has been one of the most powerful methods to study the non - equilibrium properties of mesoscopic systems . however to apply this method to set in the two - state limit , one must overcome a technical difficulty : the spin-1/2 operator , which is introduced to restrict charge number states by the strong coulomb interaction , prevents one from utilizing wick s theorem .
the most successful treatment to overcome this problem is given in ref .
@xcite , in which a formulation of perturbative expansion for the reduced density matrix in the real time domain is developed and the inelastic resonant tunneling processes is treated .
the method of ref .
@xcite enables one to classify various tunneling processes using diagrammatic techniques , and can be also applied to other systems with the strong local correlation , such as quantum dot@xcite . in spite of these successes , it seems to be still difficult to apply this method for the calculation of higher order correlation functions , since this method requires to solve a special integro - differential equation even for the calculation of the average in the presence of large quantum fluctuation@xcite .    in this paper
, we investigate the current noise and the charge noise in the regime of large quantum fluctuation of charge out of equilibrium .
we adopt the schwinger - keldysh approach and the drone - fermion representation of the effective spin-1/2 operator@xcite .
schwinger - keldysh approach enables us to calculate any order moment systematically by the functional derivative technique@xcite satisfying the charge conservation@xcite , and it helps us in manipulating many complicated terms .
the drone - fermion representation allows us to utilize the fermionic wick s theorem and to take effects of the strong correlation into account . with the help of this technique
, we can extensively take account of the higher order processes of tunneling .
we will show that our approximation reproduces the resonant tunneling approximation ( rta)@xcite as for the average current and the average charge .
the outline of this paper is as follows . in sec .
[ sec : keldysh ] , we briefly summarize the keldysh formalism and introduce an approximate generating functional .
we also show that the average and the noise expressions can be derived using the functional derivation . in sec .
[ sec : main ] , we actually calculate the average current , the average charge , the current noise and the charge noise . in sec .
[ sec : results ] we show numerical results for the noise and give some discussions on the non - equilibrium fluctuation , the thermal fluctuation and the renormalization effect .
section [ sec : summary ] summarizes our results .
in this section , we give preliminary definitions of the schwinger - keldysh approach and we summarize three useful representations : the _ closed time - path _ , the _ single time _ and the _ physical _ representations ( sec . 2 of ref .
@xcite ) . for simplicity ,
we consider the following action of a free fermion with a linear source term to explain the basic formalism : @xmath20 where the closed time path @xmath21 consists of the forward branch @xmath22 , the backward branch @xmath23 and the imaginary time - path @xmath24 as shown in fig .
[ fig : closedtimepath]@xcite .
@xmath25 is a grassmann variable satisfying the anti - periodic boundary condition @xmath26 .
the complex variable @xmath27 is defined only on the forward and backward branch @xmath28 .
the generating functional for the connected closed time - path green function ( gf ) is defined as , @xmath29 \exp \left ( -s/{i}\hbar \right ) .
\nonumber\ ] ] gf in the closed time - path representation is obtained by the second derivative of the generating functional with respect to @xmath27 ( @xmath30 ) : @xmath31 hereafter , we use arguments @xmath32 instead of @xmath33 for short .    though the closed time - path representation makes the formulation compact , in order to obtain the physical quantities
, we sometimes need the single time representation in which the time on @xmath21 is projected onto the real axis . in this representation
the degrees of freedoms of fields are doubled which we denote as , @xmath34 etc . here
@xmath35 is defined on @xmath36 and @xmath37 is the real time . in the same way as @xmath38 , gf is transformed into @xmath39 matrix in the keldysh space : @xmath40 here , arguments @xmath41 and @xmath42 are the real time and each component is defined with the statistical average in the path integral representation @xmath43 a \exp ( -s/{i}\hbar)/z \right|_{j=0}$ ] as @xmath44 .
diagonal components @xmath45 and @xmath46 are the causal and the anti - causal gf , respectively .
off - diagonal components are correlation functions , which are written in the operator representation as @xmath47 and @xmath48 @xcite . here
the statistical average is defined as @xmath49/z_0 $ ] where the hamiltonian operator @xmath50 is @xmath51 and the partition function @xmath52 is given as @xmath53 $ ] .
it is known that four components of eq .
( [ eqn : ctpgfsingle ] ) are not independent .
this redundancy is removed by the keldysh rotation@xcite .
after the rotation , we obtain the so - called physical representation : @xmath54 gfs denoted by superscripts , @xmath55 , @xmath56 and @xmath57 are advanced , retarded and keldysh components , respectively . in the practical calculations , instead of @xmath58 and @xmath59 ,
the center - of - mass coordinate @xmath60 and the relative coordinate @xmath61 , @xmath62 are used in most cases .
there are the following relations between components of gf : @xmath63 \rangle = g^c(1,2 ) .
\label{eqn : comute}\end{aligned}\ ] ] here we introduce a notation @xmath64 whose physical meaning is the spectral density in the energy space .
equation ( [ eqn : anti - comute ] ) is derived from the normalization of the step function ( see eq .
( 2.67 ) in ref .
@xcite ) .
the normalization condition of the density matrix results in an important equation @xmath65 ( see sec .
2.4 of ref .
@xcite)@xcite which is equivalent to the following equations : @xmath66
in this section , we introduce our model hamiltonian and derive the generating functional for set , based on which we construct the perturbation theory .
figure [ fig : system ] shows the equivalent circuit of a set . a metallic island exchanges electrons with the left ( right ) lead via a tunnel junction characterized by the tunneling matrix element @xmath67 .
the island is coupled to leads and a gate via capacitors @xmath68 , @xmath69 and @xmath70 .
we consider in the weak tunneling regime @xmath71 and use the two - state model . in this paper , we limit ourselves to the symmetric case : @xmath72 and @xmath73 .
the total hamiltonian consists of the unperturbed part @xmath50 and the tunneling part @xmath74 . omitting a trivial constant ,
the unperturbed part is given as @xmath75 where @xmath76 is the annihilation operator of an electron with wave vector @xmath77 in the left ( right ) lead ( @xmath78 ) or in the island ( @xmath79 ) .
the subscript @xmath80 numbers the transverse channels including spin degree of freedom . the density of states is considered as constant in each region : @xmath81 ( @xmath82 ) .
the second term is the charging energy and the effective spin-1/2 operator @xmath83 acts on the lowest two charge states .
the energy difference between two charge states is given by @xmath84 where @xmath85 is the gate charge .
the tunneling hamiltonian @xmath86 describes the electron tunneling across the junctions and simultaneous change of the charge state of the island .
@xmath87 is the phase difference between the left and the right leads and parameters @xmath88 characterizes the voltage drop between the left ( right ) lead and the island .
the tunneling hamiltonian is adiabatically turned on in the remote past and off in the distant future .
it is the widely adopted procedure , which ensures the time translational invariance to describe a stationary state .    in order to utilize wick s theorem for fermions
, we employ the mapping of the effective spin-1/2 operator onto two fermion operators @xmath89 and @xmath90 @xcite : @xmath91 , where @xmath92 is a majorana fermion operator ( @xmath93 ) .
this representation is called drone - fermion representation @xcite , because @xmath94 is a drone " whose only job is to make spin-1/2 operators of different spins commute , rather than anti - commute@xcite .    employing the hamiltonian operator in the drone - fermion representation and following the standard manner to introduce a path integral @xcite , we obtain the generating functional in the path integral representation : @xmath95 \exp \left ( -\frac{s}{{i}\hbar } \right ) ,   \label{eqn : z}\end{aligned}\ ] ] where all field variables are grassmann variables satisfying the anti - periodic boundary condition .
the action is given by @xmath96 where , @xmath97 .
( [ eqn : action0 ] ) we introduced auxiliary source fields , @xmath98 and @xmath99 , in order to calculate the average and the noise by the functional derivation .
it is noticed that the degrees of freedoms are doubled , as shown in eq .
( [ eqn : fieldsingletime ] ) . after the derivation
, these variables are put as @xmath100 and @xmath101 to be related with the parameters of the actual system .    by introducing a linear source term @xmath102 , where @xmath38 is a grassmann variable
, all fields can be traced out@xcite . in the limit of large transverse - channel number
, @xmath103 is expressed as@xcite @xmath104 \right ) \nonumber \\ & \times & \left . \exp \left ( -\frac{1}{2 { i}\hbar } \int_c { { \rm d}}1 { { \rm d}}2 \ , j(1 ) g_{\phi}(1,2 ) j(2 ) \right ) \right|_{j=0 } \nonumber \\ & \times & 2 { \rm e}^{{\rm tr } \left [ \ln g_c^{-1 } \right ] } ,   \label{eqn : startz}\end{aligned}\ ] ] where we omitted the partition function of noninteracting electrons . the trace @xmath105 and the products represent the integration along @xmath21 as follows , @xmath106 = \int_c { { \rm d}}1 { { \rm d}}2 \ , g_c(1,2 ) j(2 ) g_{\phi}(2,1 ) .
\nonumber\ ] ] the particle - hole gf @xmath107 , in the closed time - path form is written as @xmath108 where @xmath109 is the number of the transverse channels .
the gf for free electron in the lead @xmath110 and the island @xmath111 , @xmath112 , is given as @xmath113 which satisfies the anti - periodic boundary condition : @xmath114 .
the @xmath115-field and the @xmath116-field gfs defined by @xmath117 also satisfy the anti - periodic boundary condition .
( [ eqn : startz ] ) , we construct a systematic perturbation expansion in terms of @xmath118 .
for example , we show the diagrammatic representation of the zeroth ( @xmath119 ) , the first ( @xmath120 ) and the second order contribution ( @xmath121 ) to the generating functional @xmath122 in fig . [
fig : diagram1 ] . here , solid lines , dotted lines and wavy lines represent gfs for @xmath115-field , @xmath116-field and particle - hole , respectively .
practical forms are given as follows : @xmath123 \label{eqn : w0 } , \\
w^{(1 ) } & = & { i}\hbar \ , { \rm tr } \left [ g_c \sigma_c \right ] \label{eqn : w1 } , \\ w^{(2)}_{\text{$c$-field } } & = & \frac{{i}\hbar}{2 } { \rm tr } \left[(g_c\sigma_c)^2 \right ] \label{eqn : w2c},\end{aligned}\ ] ] where @xmath124 is the term corresponding to the first diagram in fig .
[ fig : diagram1 ] ( c ) , which we call the _
@xmath115-field correction_. the other four diagrams can be written in the same way . here
, the self - energy of the @xmath115-field @xmath125 is defined as @xmath126    in a previous paper@xcite , we have shown that the first order contribution causes the divergence at the degeneracy point @xmath127 for average charge . in order to regularize the divergence
, we proposed an approximate generating functional obtained by summing up @xmath115-field corrections @xmath128 to infinite order : @xmath129 \label{eqn : wrta } \\ & = & - { i}\hbar \left ( { \rm tr}\left[g_c^{-1}\right ] - \sum_{n=1}^{\infty } \frac{1}{n } { \rm tr } \left[(g_c \sigma_c)^n \right ] \right ) .
\nonumber\end{aligned}\ ] ] figure [ fig : diagram2 ] ( a ) shows the diagrammatic representation of @xmath130 .
the circle represents the self - energy of @xmath115-field and the thick line represents the full @xmath115-field gf defined by the dyson equation @xmath131 whose diagrammatic representation is shown in fig .
[ fig : diagram2 ] ( b ) . in sec .
[ sec : reformulation ] we will show that @xmath130 reproduces results of rta for the normal metal island . in sec .
[ sec : mainbody ] , we calculate the charge noise and the current noise based on @xmath130 .      in this section , we summarize expressions for the average and the noise on the basis of functional derivative .
relations between physical quantities in the generating functional representation and those in the operator representation are demonstrated in appendix [ appendix : expressions ] .
we also show that the gauge invariance of generating functional leads to the charge conservation law .
the exact average current expression is given by the functional derivative of exact generating functional @xmath122 with respect to the phase difference@xcite : @xmath132 the center - of - mass coordinate of the phase difference is determined by the josephson relation@xcite .
the relative coordinates are put to zero because they are the fictitious variables . from now on , we suppress the equations in the subscript after the vertical bar for short .
the average charge is calculated by the functional derivation in terms of the scalar potential of @xmath115-field as @xmath133 the current noise and the charge noise defined as eq .
( [ eqn : siidef ] ) are given by the second derivative with respect to @xmath134 and @xmath135 @xcite : @xmath136 where @xmath137 is obtained from the first term by replacing subscripts @xmath138 with @xmath115 .
it is convenient to rewrite eq .
( [ eqn : i - exact ] ) in the form of the linear combination of the tunneling current at the left junction @xmath139 and that at the right junction @xmath140 as @xmath141 . in the same way ,
( [ eqn : sii ] ) is written in terms of the correlation function of @xmath142 and @xmath143 ( @xmath144 ) , which we denote by @xmath145 , as @xmath146 here , @xmath147 and @xmath145 are written as @xmath148 by regarding @xmath149 and @xmath150 as formally independent variables .
the generating functional eq .
( [ eqn : z ] ) is invariant under the gauge transformation , i.e. , the phase transformation of the @xmath115-field and the change of the @xmath115-field scalar potential : @xmath151 where @xmath152 is defined on @xmath28 .
the relation between the gauge invariance and the charge conservation in the non - equilibrium state has been analyzed in ref .
for our system , the following expressions of the current continuity and the charge conservation for correlation functions @xmath153 can be proved ( appendix .
[ appendix : conservationlaw ] ) .
in this section , we derive approximate current noise and charge noise expressions , which is main purpose of this paper .
we summarize gfs in the physical representation in sec . [
sec : freegreenfunctions ] which are needed for practical calculations . in sec .
[ sec : reformulation ] we calculate the average charge and the average current and show that our approximation completely reproduces the results of rta .
we also give some notes on the diagrammatic rule suitable for the functional derivative technique . in sec .
[ sec : mainbody ] , we calculate the current noise and the charge noise based on the diagrammatic rule . we will check that our results satisfy the charge conservation law and the fluctuation - dissipation theorem .
we summarize the fourier transformation of the retarded and the keldysh component of gfs . hereafter in this subsection , we put the auxiliary source fields as shown in the subscript of eq .
( [ eqn : i - exact ] ) . the solutions of differential equations ( [ eqn : gfc ] ) and ( [ eqn : gfd ] ) imposing anti - periodic boundary condition are@xcite @xmath154 @xmath155 where @xmath156 is a positive infinitesimal number and the @xmath157-function in the energy space is defined as @xmath158 the advanced component is the complex conjugate of the retarded component : @xmath159 .
the two components of particle - hole gf are given by , @xmath160 where @xmath161 ( appendix .
[ appendix : loop ] ) .
@xmath162 is the dimensionless conductance for tunnel junction @xmath163 written in terms of the tunnel resistance @xmath164 as @xmath165 .
the spectral density of the particle - hole propagator is given by @xmath166 where the lorentzian cut - off function is introduced@xcite .
the keldysh component and the component defined by eq .
( [ eqn : comute ] ) for the @xmath115-field self - energy eq .
( [ eqn : partselfenergy ] ) and those for the particle - hole gf are related to each other : @xmath167 the retarded component of the self - energy is given as @xmath168 where @xmath169 is the digamma function .
the full @xmath115-field gf is obtained by solving the dyson equation in the closed time - path representation eq .
( [ eqn : fullgc])@xcite : @xmath170 where we used the definition of @xmath157-function eq .
( [ eqn : deltae ] ) . here
we remark the following : equation ( [ eqn : cg ] ) shows that in the limit of @xmath171 , the charge states are independent of the initial equilibrium distribution , because the keldysh component of @xmath115-field gf represents the distribution of charge states .
this fact suggests that @xmath130 describes a physically reasonable non - equilibrium stationary state , which should not depend on any initial state .
we transform above expressions into single - time representation .
employing eqs .
( [ eqn : anti - comute ] ) and ( [ eqn : comute ] ) we obtain @xmath172 here @xmath173 and @xmath174 are written with the fermi function @xmath175 and the bose distribution function @xmath176 .
functions @xmath177 and @xmath178 are given by @xmath179 and @xmath180 , respectively .
the reason why we adopt the generating functional approach is that once an approximate generating functional is obtained , one can calculate any order moment systematically by the functional derivation . in the following sections we perform practical calculations , employing @xmath130 introduced in sec .
[ subsec : w ] .
detailed discussions on our approximation are retained in appendix [ appendix : perturbartion ] .
the average current is calculated by substituting @xmath130 into eq .
( [ eqn : ir])@xcite : @xmath181 \right| \label{eqn : trace1}.\end{aligned}\ ] ] here , we used the fact that the self - energy eq .
( [ eqn : partselfenergy ] ) includes the phase factor through the particle - hole gf eq .
( [ eqn : particlehole ] ) . in the language of feynmann diagrams , eq .
( [ eqn : trace1 ] ) can be rewritten in a compact form : @xmath182{diagram11a.eps } \end{array } - \begin{array}{c } \includegraphics[width=7ex]{diagram11b.eps } \end{array } \right ] \right| \label{eqn : i1}.\end{aligned}\ ] ] solid dots with @xmath37 represents @xmath183 .
the circle with @xmath163 is the partial self - energy defined by eq .
( [ eqn : partselfenergy ] ) . here
we obtain the diagrammatic rule similar to that of refs .
@xcite .
( i ) : :    the diagrams corresponding to the functional derivative with respect    to @xmath184 is obtained by a    series of the following operations ; to put a solid dot onto all    possible positions of vertices in a closed diagram , to assign    @xmath163 on a circle connected to the solid dot , to    multiple @xmath185 , and to assign minus sign if a solid line    comes into the solid dot .
next we project the fictitious time on @xmath21 to the real axis .
as the tunneling hamiltonian is zero on @xmath24 , @xmath186 is zero for @xmath187 or @xmath188 .
hence , the diagrams in eq .
( [ eqn : i1 ] ) are rewritten as , @xmath189 \right| , \label{eqn : trace2}\end{aligned}\ ] ] where we performed the keldysh rotation .
@xmath190 is the pauli matrix in the keldysh space @xcite .
the trace is carried out over the keldysh space and products represents the integration along the real - time as @xmath191 = \int_{-\infty}^{\infty } { { \rm d}}1 { { \rm d}}2 \ , { \rm tr } \left [ \tilde{g}_{\phi}(1,2 ) \tilde{\varphi}(2 ) \tilde{g}_{\phi}(2,1 ) \right ] .
\nonumber\ ] ] the phase in the physical representation @xmath192 is written as @xmath193 which leads to a useful relation for the functional derivative technique : @xmath194 by using the property of the gf in the physical representation @xmath195 , and that of a pauli matrix @xmath196 , we can see that the second term of eq .
( [ eqn : trace2 ] ) is minus the complex conjugate of the first term . to generalize this property
, we obtain a helpful rule to reduce the number of diagrams : a diagram which is complex conjugate of a certain diagram is obtained by changing the direction of all lines and putting minus sign if the diagram includes odd numbers of vertices without solid dot .    by performing the fourier transformation , we rewrite eq .
( [ eqn : i1 ] ) as : @xmath197{diagram11a.eps } \end{array } \right| = 2 e{\rm re } \int \frac{{{\rm d}}\varepsilon}{h } { \rm tr } \left [ \tilde{g}_c(\varepsilon ) \frac{{\mbox{\boldmath$\tau$}}^0}{2 }   \tilde{\sigma}_{\rm r}(\varepsilon ) { \mbox{\boldmath$\tau$}}^1 \right ] \nonumber \\ \label{eqn : diagram - trace } \\ & = & e \int \frac{{{\rm d}}\varepsilon}{h } \frac { \sigma_{\rm r}^c(\varepsilon ) g_c^k(\varepsilon ) } { 2 } - ( c \leftrightarrow k ) .
\label{eqn : iphysical}\end{aligned}\ ] ] the second term @xmath198 is obtained from the first term by swapping superscripts @xmath57 and @xmath21 . by generalizing eq .
( [ eqn : diagram - trace ] ) , we obtain the rule to calculate the diagram :    ( ii ) : :    put @xmath39 matrix gfs in the physical representation    to the corresponding lines and circles , put    @xmath199 or    @xmath200 to a vertex with or without    solid dot .
subsequently , carry out the trace over keldysh space and    the integration over the frequency @xmath201 .    employing eqs .
( [ eqn : anti - comute ] ) and ( [ eqn : comute ] ) , eq .
( [ eqn : iphysical ] ) is transformed into the single time representation : @xmath202 by using eqs .
( [ eqn : correlationalpha ] ) , ( [ eqn : correlationself ] ) and ( [ eqn : correlatonfullgc1 ] ) , and noting @xmath203 , we obtain the final form of the current expression which has the same form as the landauer formula : @xmath204 where , @xmath205 is the effective transmission probability of lead electrons thorough the island : @xmath206 from the second form , we can see our result is equivalent to that of rta@xcite .    in the same way as the average current , the average charge is evaluated by substitute @xmath130 into eq .
( [ eqn : q - exact ] ) @xcite , @xmath207{diagram12.eps } \end{array } \right| , \nonumber\end{aligned}\ ] ] where the solid dot with @xmath37 corresponds to @xmath208 whose practical expression is equal to the right - hand side of eq .
( [ eqn : deltaphysical1r ] ) in the physical representation .
the diagrammatic rule of functional derivation in terms of the scalar potential is given as follows :    ( i ) : :    the diagrams corresponding to the functional derivative with respect    to @xmath98 is obtained by inserting a solid dot into all    possible positions of @xmath115-field gf .    following the rule * ( ii ) * , the diagram for average charge can be also calculated : @xmath209 \nonumber \\ & = & \frac{1}{2}+\int \frac{{{\rm d}}\varepsilon}{\pi } \frac{\alpha^r(\varepsilon)}{\alpha^k(\varepsilon ) } { \rm i m } g_c^r ( \varepsilon ) \label{eqn : qrta}.\end{aligned}\ ] ] the imaginary part of @xmath210 , which has a peak at @xmath211 , describes the excitation property of the charge state .
when the broadening of the peak is sufficiently small and @xmath212 , @xmath210 is approximately given by @xmath213 here @xmath214 is the renormalization factor @xmath215 .
the low energy cut - off @xmath216 is @xmath217 where one of three parameters must be much larger than the other two@xcite . the imaginary part of self - energy @xmath218 represents the life - time broadening effect .
when the renormalization effect is negligible @xmath219 , it is written as @xmath220 , where @xmath221 and @xmath222 @xmath223 ( @xmath224 ) is equal to the tunneling rate into ( out of ) the island through the junction r estimated by fermi s golden rule .
it is noticed that the condition @xmath225 is enough to neglect the broadening of the peak . by using the approximation eq .
( [ eqn : approxg ] ) , eq .
( [ eqn : qrta ] ) for equilibrium state reproduces the result of rta@xcite , @xmath226 . according to ref .
@xcite , eq .
( [ eqn : irta ] ) is consistent with the co - tunneling theory and orthodox theory for two - state limit@xcite .
we can also confirm that eq .
( [ eqn : qrta ] ) is consistent with the orthodox theory in the following way : in the limit of @xmath227 with keeping @xmath228 or @xmath6 finite , the renormalization factor approaches unity and the imaginary part of eq .
( [ eqn : approxg ] ) reduces to the @xmath157-function @xmath229 .
thus , eq . ( [ eqn : qrta ] )
reproduces the result of the orthodox theory @xmath230 where @xmath231 .      here , we calculate the current noise and the charge noise based on the reformulated rta . in the following discussions we limit ourselves to zero frequency component @xmath232 where we used a fact that at a stationary state , the time translational invariance is satisfied and the correlation function depends only on the difference of @xmath37 and @xmath233 . from the definition eq .
( [ eqn : sii ] ) , one obtain the current noise by applying the rule * ( i ) * twice . reducing the number of diagrams and using the rule * ( ii ) * , which can be also applied to the calculation of zero frequency noise diagram ( appendix .
[ appendix : noiserule ] ) , we obtain    @xmath234{diagram13a.eps } \end{array } \delta_{\rm r , r ' } + \begin{array}{c } \includegraphics[width=8ex]{diagram13b.eps } \end{array } - \begin{array}{c } \includegraphics[width=8ex]{diagram13c.eps } \end{array } \right ] \right| - \frac{(\delta \rightarrow c)}{4 } \nonumber \\ & = & 2 \int \frac{{{\rm d}}\varepsilon}{h } { \rm tr } \left [ \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{\sigma}_{\rm r}(\varepsilon )   \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{g}_c(\varepsilon )   \ , \delta_{\rm r , r ' } + \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{\sigma}_{\rm r}(\varepsilon ) { \mbox{\boldmath$\tau$}}^1   \tilde{g}_c(\varepsilon ) { \mbox{\boldmath$\tau$}}^1 \tilde{\sigma}_{\rm r'}(\varepsilon ) \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{g}_c(\varepsilon ) - \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{\sigma}_{\rm r}(\varepsilon ) { \mbox{\boldmath$\tau$}}^1 \tilde{g}_c(\varepsilon ) \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{\sigma}_{\rm r'}(\varepsilon ) { \mbox{\boldmath$\tau$}}^1 \tilde{g}_c(\varepsilon ) \right ] \nonumber \\ & - & ( { \mbox{\boldmath$\tau$}}^0 \rightarrow { \mbox{\boldmath$\tau$}}^1 ) ,   \nonumber\end{aligned}\ ] ]    where we used a fact that the second term of eq . ( [ eqn : sii ] ) is obtained from the first term by changing @xmath235 to @xmath200 in the physical representation ( appendix .
[ appendix : noiserule ] ) .
after some straightforward calculations , we obtain the following expression : @xmath236 \label{eqn : siirta } \\ & = & \frac{2}{r_{\rm k } } \int { { \rm d}}\varepsilon [ t^f(\varepsilon ) \ { f^{-}_{\rm l}(\varepsilon ) f^{+}_{\rm r}(\varepsilon ) + f^{+}_{\rm l}(\varepsilon ) f^{-}_{\rm r}(\varepsilon ) \ } - t^f(\varepsilon)^2 \ { f^{-}_{\rm l}(\varepsilon ) - f^{-}_{\rm r}(\varepsilon ) \}^2 ] \label{eqn : siirta_2},\end{aligned}\ ] ] where @xmath237 is the other side of @xmath163 ( for example , when the index @xmath238 is @xmath239 the index of the other side @xmath240 is @xmath241 ) .
the charge noise is evaluated by adopting the definition eq .
( [ eqn : sqq ] ) and by applying rules * ( i ) * and * ( ii ) * : @xmath242{diagram14.eps } \end{array } \right| -\frac{(\delta \leftrightarrow c)}{4 } = e^2 \hbar^2 \int \frac{{{\rm d}}\varepsilon}{h } { \rm tr } \left [ \tilde{g}_c(\varepsilon ) \frac{{\mbox{\boldmath$\tau$}}^0}{2 }   \tilde{g}_c(\varepsilon ) \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \right ] - ( { \mbox{\boldmath$\tau$}}^0 \rightarrow { \mbox{\boldmath$\tau$}}^1 ) \nonumber \\ & = & e^2 \hbar^2 \int \frac{{{\rm d}}\varepsilon}{h } g_c^{-+}(\varepsilon ) g_c^{+-}(\varepsilon ) = \sum_{\rm r , r'=l , r } \frac{e^4 r_{\rm k}}{2 \pi^2 } \int { { \rm d}}\varepsilon \left| g_c^r(\varepsilon ) \right|^4 \alpha^k_{\rm r}(\varepsilon ) \alpha^k_{\rm r'}(\varepsilon ) f^{+}_{\rm r}(\varepsilon ) f^{-}_{\rm r'}(\varepsilon ) , \label{eqn : sqqrta}\end{aligned}\ ] ]    where we used eqs .
( [ eqn : correlationalpha ] ) and ( [ eqn : correlatonfullgc1 ] ) .    equation ( [ eqn : siirta_2 ] ) has the same form as the current noise expression of a point contact without coulomb interaction@xcite .
this result is anticipated , because the tunneling current is expressed in the same form as the landauer formula .
however there is an important difference as mentioned in ref .
@xcite : the effective transmission probability includes the coulomb correlation and the inelastic relaxation effect .    here we note the following points .
first , our approximation satisfies the charge conservation law : as the approximate generating functional @xmath130 is invariant under the gauge transformation eq .
( [ eqn : gauge_t ] ) , we can show the conservation law , eqs .
( [ eqn : iconserve ] ) and ( [ eqn : sconserve ] ) , for the approximate expressions . especially for zero frequency component
, they reduce to equations for the current conservation law , @xmath243 and @xmath244 .
it is worth noticing that the gauge invariance is automatically satisfied when @xmath130 consists of closed diagrams .
secondly , we can show that our result satisfies the fluctuation - dissipation theorem : at @xmath245 , the current noise expression is reduced to the johnson - nyquist formula @xmath246 where @xmath247 is the conductance expressed as @xmath248 .
we discuss the current noise expression in the regime where @xmath225 and the renormalization effect is negligible , @xmath219 . using eqs .
( [ eqn : transmission ] ) and ( [ eqn : approxg ] ) , the first term of eq .
( [ eqn : siirta ] ) is approximately given by @xmath249 where @xmath250 for @xmath251 , @xmath252 is the co - tunneling current from lead r/@xmath253 to lead @xmath253/r , and eq .
( [ eqn : siicotwidth ] ) reproduces the co - tunneling theory @xcite in the two - state limit .
equation ( [ eqn : siicotwidth ] ) , which we call the _ co - tunneling theory with life - time broadening _ , is equivalent to the previously proposed equation in ref .
@xcite in the limit of @xmath254 .
the validity of this approximation is discussed in the next section from the point of view of the numerical results .    by taking the limit @xmath227 , with paying attention to the relation @xmath255
\}^2 = \delta(\varepsilon-\delta_0 ) $ ] , we obtain the result of the orthodox theory in the two - state limit : @xmath256 = \alpha_0^{-1 } 2 e \left ( i_{+}- \frac{2 i_{-}^2}{e \gamma } \right ) \label{eqn : siiorth},\ ] ] where @xmath257 .
@xmath258 is the average current of the orthodox theory .
the second term is related to the part of eq .
( [ eqn : siirta ] ) proportional to @xmath259 , which represents the coulomb correlation and reduces the current noise from the poissonian value .
it is noticed that even at small tunneling conductance , the second term is important in st regime .    in the same way as above discussions , for @xmath260 and @xmath219
, we obtain the co - tunneling theory@xcite with life - time broadening from eq .
( [ eqn : sqqrta ] ) : @xmath261 in the limit of @xmath227 , we can confirm that our result reproduces the orthodox theory@xcite : @xmath262 = \alpha_0 4 e^2 \gamma_{+ } \gamma_{-}/ \gamma^3 \label{eqn : sqqorth}.\ ] ]
in this section we present results obtained by the numerical calculation of eqs .
( [ eqn : siirta ] ) and ( [ eqn : sqqrta])@xcite , based on which we discuss the non - equilibrium fluctuation in sec .
[ sec : discussion1 ] and that of thermal fluctuation in sec .
[ sec : discussion2 ] .
we also discuss the effect of renormalization on the noise at low temperature caused by quantum fluctuation in sec .
[ sec : discussion3 ] .      in order to clarify the nature of non - equilibrium current fluctuation
, we consider the case of zero temperature , where there is no thermal fluctuation .
furthermore , we limit ourselves to the condition of the high bias voltage , @xmath263 , where the renormalization effect is negligible . figure [ fig : s ] shows the current noise ( the top panel ) and the charge noise ( the bottom panel ) as a function of the excitation energy @xmath264 for small @xmath11 ( ( a-1 ) and ( b-1 ) ) and those for large @xmath11 ( ( a-2 ) and ( b-2 ) ) .
plots are normalized by values of the orthodox theory at @xmath127 .
solid lines , dashed lines and dotted lines show results of our approximation ( eqs .
( [ eqn : siirta ] ) and ( [ eqn : sqqrta ] ) ) , the orthodox theory ( eqs .
( [ eqn : siiorth ] ) and ( [ eqn : sqqorth ] ) ) and the co - tunneling theory ( eqs .
( [ eqn : siicotwidth ] ) and ( [ eqn : sqqcotwidth ] ) with @xmath251 ) , respectively .
when @xmath11 is small , our results well reproduce the orthodox theory ( ( a-1 ) and ( b-1 ) ) . for large @xmath11 ( ( a-2 ) and ( b-2 ) ) and in cb regime ( @xmath265 ) , they agree well with the co - tunneling theory .
figure [ fig : in ] shows the average current ( a ) and the average charge ( b ) estimated by our approximation ( solid lines ) and the orthodox theory ( dashed lines ) . in cb regime ,
we can see both the average value and the noise are enhanced by the quantum fluctuation .
around @xmath127 , the average current and the current noise are strongly suppressed due to the life - time broadening effect @xcite .
next , we discuss the validity of the co - tunneling theory with life - time broadening ( see eqs .
( [ eqn : siicotwidth ] ) and ( [ eqn : sqqcotwidth ] ) ; dot - dashed lines in fig . [
fig : s ] ) . as for the charge noise in the limit of @xmath227
, it reproduces the result of the orthodox theory as well as our approximation ( fig .
[ fig : s ] ( b-1 ) : the solid line and the dot - dashed line almost overlap each other ) . in st regime
it overestimates the current noise ( fig .
[ fig : s ] ( a-1 ) ) , because it does not take the coulomb correlation effect into account as mentioned before .
we want to stress again that our result reproduces the orthodox theory in the limit of @xmath227 , which can not be achieved by the co - tunneling theory with life - time broadening .    the physical picture of the non - equilibrium current fluctuation
is understood more clearly with the help of the fano factor defined by @xmath266 .
the fano factor is unity when the tunneling event of electrons is poissonian process and is suppressed when tunneling events are correlated .
the orthodox theory predicts the sub - poissonian behavior in st regime because of the coulomb correlation : suppose one electron tunnels into the island through one junction , the next tunneling event must be the out going process of another electron through the other junction .
in cb regime , the co - tunneling theory predict @xmath267 .
it means that co - tunneling events , viz .
the simultaneous tunneling events of two electrons through the two junctions , occur randomly .
the figure [ fig : fano ] shows the excitation energy dependence of the fano factor obtained by our approximation . in the small @xmath11 limit
, our approximation reproduces the orthodox theory in st regime and the co - tunneling theory in cb regime and smoothly interpolate two theories ( the dotted line actually almost coincides with the result of orthodox theory in st regime ) . for larger @xmath11 the fano factor
is further suppressed ( the dashed and the solid line ) .
especially , our result predicts the value smaller than @xmath268 at the degeneracy point .
it is a distinctive result because the orthodox theory predicts the inequality @xmath269 ( p. 137
@xcite ) .
next we consider the physical meaning of our result . in cb regime and near the threshold voltage ,
the origin of the suppression of fano factor is considered to be the enhancement of the effective transmission probability @xmath205 , because the increase in the current noise is much larger then that in the charge noise ( see figs .
[ fig : s ] ( a-2 ) and [ fig : s ] ( b-2 ) ) . as the fano factor of the shot noise
is approximately given by @xmath270@xcite , the enhancement of the transmission probability results in the suppression of the fano factor .
around the degeneracy point , the origin of the suppression of fano factor is different , because the normalized current and the current noise are suppressed as shown in fig . [
fig : in ] ( a ) and fig .
[ fig : s ] ( a-2 ) .
we consider that the fano factor is suppressed because of the dissipation , i.e. the life - time broadening effect : rta takes account of the dissipation process which is the leak of an electron from the island while another electron tunnels into the island and relaxes to the local equilibrium state of the island .
the suppression of the fano factor by the dissipation was previously predicted for the 1d electron channel coupled with a boson bath@xcite .
next we discuss the effect of the thermal fluctuation .
figure [ fig : nyquist ] shows the temperature dependence of the current noise ( a ) , the average current ( b ) and the fano factor ( c ) at a threshold for various @xmath11 .
they are normalized by the value of the orthodox theory at @xmath271 .    as the temperature increases , the average current and the current noise increase because of the thermal fluctuation . at sufficiently high temperature @xmath272
, orthodox theory predicts that the average current saturates at @xmath273 and the thermal fluctuation dominates the current noise , i.e. @xmath274 which is the similar form as the johnson - nyquist noise for the ohmic resistance@xcite ( the plot for @xmath275 almost coincides with the orthodox theory ) .
our result further shows that the average current and the current noise are suppressed as @xmath11 increases ( panels ( a ) and ( b ) ) .
it is considered to be attributed to the higher order tunneling effect : the lifetime broadening caused by the thermal fluctuation is enhanced for the large tunnel conductance .
the panel ( c ) is the fano factor versus temperature plot .
the fano factor is independent of @xmath11 , which means that the correlation between tunneling events does not depend on @xmath11 .
next we consider the renormalization effect at low bias voltage and temperature : @xmath276 .
figures [ fig : sqq - low ] ( a ) and ( b ) show the charge noise normalized by @xmath277 at @xmath127 as a function of the temperature and the bias voltage , respectively .
the charge noise is suppressed for large @xmath11 , which is attributed to the renormalization of the system parameters . in the regime @xmath278 , where the life - time broadening effect is negligible ( @xmath225 )
, we can approximate eq .
( [ eqn : sqqrta ] ) as @xmath279 instead of eq .
( [ eqn : sqqorth ] ) . here @xmath280 and @xmath281 are tunneling rate eq .
( [ eqn : tunnelingrate ] ) written by using renormalized parameters such as @xmath282 where the renormalization factor @xmath214 is @xmath283 .
the lower cut - off energy @xmath216 is @xmath284 for the panel ( a ) and @xmath285 for the panel ( b ) .
it is natural to interpret eq .
( [ eqn : chargernormal ] ) that the charge of a carrier is modified as @xmath286 by the renormalization effect .
the interpretation is similar to that of the doubling of shot noise at the normal - metal(n)-superconductor(s ) interface . since at ns interface the carrier is @xmath287-charged particle , viz .
a cooper pair , the shot noise is twice as large as that at nn interface @xcite .
though the normalized charge noise is suppressed with increasing @xmath11 , one see that the charge noise always diverges at @xmath288 for arbitrary @xmath11 in the weak tunneling regime .
since the charge noise is related to the charge susceptibility `` for excitation energy @xmath264 , the divergence means that the number of charge changes by one '' at the degeneracy point when we sweep the excitation energy .
it is confirmed by the fact that the slope of excitation energy dependence of the average charge ( an inset of fig .
[ fig : sqq - low ] ( b ) ) diverges at the degeneracy point .
next we discuss the renormalization effect on the fano factor .
figure [ fig : fano3 ] shows the excitation energy dependence of the fano factor for various bias voltage .
we can see at small bias voltage where the charging energy renormalization is pronounced @xmath289 , the valley structures of curves are widened .
the same behavior can be seen in the differential conductance shown in refs .
we also see that the fano factor is suppressed with increasing bias voltage at @xmath127 .
this suppression is also attributed to the dissipation as discussed in sec .
[ sec : discussion1 ] : as the bias voltage increases , the dissipative charge fluctuation is enhanced and thus the fano factor is suppressed .
by using the drone - fermion representation and the schwinger - keldysh approach , we have calculated the current noise and the charge noise in the regime of large quantum fluctuations of charge out of equilibrium .
we have reformulated and extended rta in a charge conserving way .
our approximation interpolates previous theories , the orthodox theory and the co - tunneling theory : our result coincides with the orthodox theory in the limit of @xmath227 and is consistent with the co - tunneling theory in cb regime .
the approximation is verified from the fact that the result satisfies the fluctuation - dissipation theorem . in previous papers
, we also checked numerically that the energy sensitivity does not exceed the quantum limit@xcite .
we showed that at zero temperature and @xmath263 , the life - time broadening caused by non - equilibrium dissipative charge fluctuation suppresses the current noise in st regime .
it also suppresses the fano factor more than the coulomb correlation does .
especially the fano factor is suppressed below the minimum value predicted by the orthodox theory @xmath268 around @xmath127 .
the origin of the suppression is attributed to the charge fluctuation which appears as the enhancement of the transmission probability in cb regime and the dissipation in st regime .    at @xmath290
, we showed that the average current and the current noise deviate from the predictions of the orthodox theory with increasing @xmath11 .
however , the fano factor is independent of @xmath11 and is proportional to the temperature .
it means that the current noise is dominated by the thermal fluctuation and the correlation between the tunneling events does not depends on @xmath11 .    at small bias voltage and temperature @xmath276 ,
the charge noise is suppressed as compared with the prediction of the orthodox theory .
we showed that it can be interpreted as the renormalization for the unit of island charge .
although the charge is renormalized , the charge noise diverges at @xmath288 for arbitrary @xmath11 in the weak tunneling regime .
it means that the quantum fluctuation does not wash out the charge quantization .    in this paper , we have limited ourselves to the discussions on the second moment and the zero frequency component , because we think them primitive .
the investigation on the frequency dependence of noise will be important to estimate the performance of high - speed set electrometer completely@xcite .
the investigation on the higher order moment and the full counting statistics@xcite will help us to understand carriers of strongly correlated system out of equilibrium .
we would like to thank y. isawa , j. martinek , yu .
v. nazarov and g. johansson for variable discussions and comments .
this work was supported by a grant - in - aid for scientific research ( c ) , no .
14540321 from mext .
was supported by mext , grant - in - aid for encouragement of young scientists , no .
in this appendix , we show the relation between expressions for average and noise in the operator representation and those in the generating functional representation .
the variation of the exact action eq .
( [ eqn : action0 ] ) accompanied by the infinitesimal variation of @xmath291 is given as the twisted " combination ( sec .
9.3.2 in ref .
@xcite ) of @xmath292 and @xmath291 , @xmath293 where the center - of - mass coordinate of charge is @xmath294 . employing this source term , we can show that the eq .
( [ eqn : q - exact ] ) is equivalent to @xmath295 where @xmath296 : @xmath297 where we used eq .
( [ eqn : anti - comute ] ) to obtain the final form .
the second derivative @xmath298 is calculated as @xmath299 here the first term includes the correlation function of 4 field variables on a same branch , @xmath300 , which is not well defined at @xmath301 .
usually , an additional operation to determine the order of field variables is required to remove the uncertainty .
alternatively , we subtract a term @xmath302 from eq .
( [ eqn : sqq_k ] ) as shown in our definition eq .
( [ eqn : sqq ] ) .
here we used the normalization eq .
( [ eqn : normalization ] ) , @xmath303 . as a result , the first term in eq .
( [ eqn : sqq_k ] ) is replaced by @xmath304 , which does not include the uncertainty .
( [ eqn : anti - comute ] ) , we show that our definition is equal to the standard charge noise expression @xmath305 we should stress that our definition does not change the final result because eq . ( [ eqn : sqq2nd ] ) is 0 from the normalization .
the exact current expression is obtained in the same way . the source term corresponding to eq .
( [ eqn : ds ] ) is given by @xmath306 which leads to relations @xmath307 where the current operator at junction @xmath163 is defined as @xmath308
in this appendix , we demonstrate the charge conservation law .
as @xmath122 is invariant under the transformation eq .
( [ eqn : gauge_t ] ) , we obtain an identity @xmath309 we can derive the other equation , which is obtained from above equation by replacing @xmath138 with @xmath115 . however , the latter equation is not important in the following discussions . by putting the auxiliary source fields as the values given in the subscripts of eq .
( [ eqn : i - exact ] ) , and employing eqs .
( [ eqn : i - exact ] ) and ( [ eqn : q - exact ] ) , we obtain the current continuity equation , eq .
( [ eqn : iconserve ] ) .
next we demonstrate the charge conservation law for correlation functions . by acting the operator @xmath310 or @xmath311 on eq .
( [ eqn : identityw1 ] ) , we obtain following two equations : @xmath312 by comparing the left - hand side of the former equation and the right - hand side of the latter equation , by setting the auxiliary source fields as the values given in the subscripts of eq .
( [ eqn : i - exact ] ) , and by using eqs .
( [ eqn : sii ] ) and ( [ eqn : sqq ] ) , we obtain the charge conservation law for correlation functions , eq .
( [ eqn : sconserve ] ) .
in this appendix , we calculate the particle - hole gf . we begin with the tunneling action for the large transverse channel obtained from eq .
( [ eqn : action0 ] ) by tracing out the electron degrees of freedom@xcite : @xmath313 . in the physical representation , it is rewritten as @xmath314 , where the vector field @xmath315 is defined in the same way as eq .
( [ eqn : fieldphysical ] ) .
each component of @xmath316 can be calculated by utilizing the functional derivation .
for example , ( 1,2)-component is @xmath317 .
\nonumber \\ \label{eqn : tracea1}\end{aligned}\ ] ] as the functional derivative is @xmath318 in the physical representation , the trace yields @xmath319/2 $ ] . here
we put @xmath320 and @xmath321 .
the other components are evaluated in the same way .
then four components are given as @xmath322 where we omitted arguments and coefficients . to obtain this form
, we used the normalization eq .
( [ eqn : normalization ] ) or the relation @xmath323 ( see eqs . ( 2.64 ) and ( 2.65 ) in ref .
@xcite ) . the following calculations are same as those in ref .
@xcite . employing the fourier transforms of the gf defined in eq . ( [ eqn : gfl ] ) , @xmath324 and @xmath325 , we obtain eq .
( [ eqn : alphark ] ) .
another loop diagram , the self - energy eq .
( [ eqn : partselfenergy ] ) , can be calculated in the same way .
four components are given in the same form as those of eq .
( [ eqn : ph ] ) . by using eqs .
( [ eqn : drk ] ) and ( [ eqn : alphark ] ) , we obtain @xmath326 from these equations , eqs .
( [ eqn : selfenergy ] ) and ( [ eqn : selfenergyr ] ) can be derived .
in this appendix , we describe some results of finite order perturbation theory and explain why we introduced eq .
( [ eqn : wrta ] ) which takes account of all orders for @xmath115-field corrections .
first , we calculate the first order contribution of average charge by employing @xmath120 , eq .
( [ eqn : q - exact ] ) and rules * ( i ) * and * ( ii)*.    @xmath327{diagram15.eps } \end{array } \right| = -\frac{{i}}{2 \pi } \int { { \rm d}}\varepsilon { \rm tr } \left [ \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{g}_c(\varepsilon ) { \mbox{\boldmath$\tau$}}^1   \tilde{\sigma}_c(\varepsilon ) { \mbox{\boldmath$\tau$}}^1   \tilde{g}_c(\varepsilon )   \right ] \nonumber \\ & = & -{i}\int \frac{{{\rm d}}\varepsilon}{8 \pi } ( g_c^k(\varepsilon ) g_c^p(\varepsilon ) \sigma_c^p(\varepsilon ) + g_c^k(\varepsilon ) g_c^c(\varepsilon ) \sigma_c^c(\varepsilon)+ 2 |g_c^r(\varepsilon)|^2 \sigma_c^k(\varepsilon ) ) , \nonumber\end{aligned}\ ] ]    where gf denoted with superscript @xmath328 is given as @xmath329 etc .
in the equilibrium state , the second and the third terms of the second line are negligibly small @xmath330 and the first term is simplified to @xmath331 where we utilized the relation @xmath332 . in the limit of zero temperature , eq .
( [ eqn : q2nd ] ) leads to the log - divergence@xcite as @xmath333 .    above result
suggests that the @xmath115-field correction is responsible for the divergence .
it is further confirmed by calculating second order contribution of average current generated from @xmath115-field correction @xmath334 . employing rules * ( i ) * and * ( ii ) * , we obtain @xmath335{diagram16.eps } \end{array } \right| \nonumber \\ & = & 2 e { \rm re } \int \frac{{{\rm d}}\varepsilon}{h } { \rm tr } \left [ \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{\sigma}_{\rm r}(\varepsilon )   { \mbox{\boldmath$\tau$}}^1 \tilde{g}_c(\varepsilon ) { \mbox{\boldmath$\tau$}}^1   \tilde{\sigma}_c(\varepsilon ) { \mbox{\boldmath$\tau$}}^1   \tilde{g}_c(\varepsilon ) \right ] \nonumber \\ & = & e \int \frac{{{\rm d}}\varepsilon}{2 h }    \sigma^k_{\rm \bar{r}}(\varepsilon ) \sigma^c_{\rm
r}(\varepsilon)-(c \leftrightarrow k ) + \delta i^{(2)}_{\rm r } \nonumber \\ & = & \gamma^{+}_{\rm r}(\eta ) - \gamma^{-}_{\rm r}(\eta ) + \delta i^{(2)}_{\rm r}.   \label{eqn : i2}\end{aligned}\ ] ] the first and the second terms of eq .
( [ eqn : i2 ] ) , which are consistent of the expression for co - tunneling current @xcite , also diverge at the degeneracy point . from above discussions , we can deduce that the most simple way to regulate the divergence is to sum up @xmath115-field corrections @xmath128 up to infinite @xmath80 as shown in eq .
( [ eqn : wrta ] ) .
it should be noted that the correction term @xmath336 diverges in the limit of @xmath171 .
this divergence disappears when we consider eq .
( [ eqn : wrta ] ) as discussed in sec .
[ sec : reformulation ] and sec .
[ sec : mainbody ] .
in this appendix , we demonstrate that the rule * ( ii ) * can be also applied to the calculation of zero - frequency noise
. we also demonstrate the rule to calculate the second term of the noise expression eq .
( [ eqn : sii ] ) or eq .
( [ eqn : sqq ] ) .
for example , we consider the current noise related to @xmath120 . though we consider the simple case , the following discussions can be generalized . from the definition eq .
( [ eqn : siipart ] ) , the noise diagrams is obtained as : @xmath337{diagram17.eps } \end{array } \right| \delta_{\rm r , r ' } - \frac{(\delta \rightarrow c)}{4}. \label{eqn : sii_1}\ ] ] the integration in terms of @xmath37 in the first term is calculated as @xmath338{diagram17.eps } \end{array } \right| = \int { { \rm d}}t { \rm tr } \left . \left
[ \frac{\delta \varphi_{\rm r}}{\delta \varphi_{{\rm r } \ , \delta}(t ) } \sigma_{\rm r } \frac{\delta \varphi_{\rm r'}}{\delta \varphi_{{\rm r ' } \ , \delta}(t ' ) } g_c \right ] \right| \nonumber \\ & = & \int \frac{{{\rm d}}\varepsilon}{h } { \rm tr } \left [ \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{\sigma}_{\rm r}(\varepsilon ) \frac{{\mbox{\boldmath$\tau$}}^0}{2 } \tilde{g}_c(\varepsilon ) \right ] .
\label{eqn : trace4}\end{aligned}\ ] ] and thus we can see that the rule * ( ii ) * can be applied for the calculation of the zero - frequency current noise . as for the second term
, the derivation @xmath339 appears instead of eq .
( [ eqn : deltaphysical1r ] ) .
thus we can derive the rule that the second term is obtained from the first term by replacing @xmath235 with @xmath200 as @xmath340/h $ ] @xcite . as a result ,
( [ eqn : sii_1 ] ) is expressed as @xmath341 \delta_{\rm r , r'}/(2 h ) - ( { \mbox{\boldmath$\tau$}}^0 \rightarrow { \mbox{\boldmath$\tau$}}^1 ) .
\nonumber\ ] ] as for the charge noise , we can repeat the same discussions as above .