the single - impurity kondo model  @xcite is unquestionably one of the most important models studied in condensed matter physics over the past decades .
the investigation of its equilibrium properties has caused the development of important theoretical tools such as renormalization group methods or the bethe ansatz for impurity systems  @xcite .
more than 20 years ago it was also realized  @xcite that the kondo model can be used to describe transport experiments through quantum dots .
the developments in the ability to engineer devices on the nanoscale has led to the experimental realization of kondo physics in such systems  @xcite .
one particular advantage of these quantum dots is the almost full control over system parameters like temperature , bias and gate voltages , magnetic field , and exchange couplings .
these possibilities have triggered a great interest in the theoretical study of quantum dots out of equilibrium .
a wide range of theoretical methods has been applied in the past , including non - equilibrium perturbation theory  @xcite , the flow - equation method  @xcite , coulomb gas representations  @xcite , real - time and functional renormalization group methods  @xcite , a non - equilibrium extension of the numerical renormalization group ( nrg ) method  @xcite , and time - dependent density matrix renormalization group ( dmrg ) techniques  @xcite .
these studies established the importance of relaxation and decoherence effects for the understanding of non - equilibrium physics . from an experimental point of view
various measurable quantities like the steady - state current , the magnetization , and the static susceptibility have been calculated .    beside the application of these techniques
there were also attempts to employ the known integrability of certain impurity models in equilibrium , notably the anderson impurity model and the interacting resonant level model , to investigate their non - equilibrium properties .
konik et al .
@xcite calculated the differential conductance in the anderson impurity model by combining the well - known scattering states of the equilibrium system  @xcite with a landauer - bttiker formalism . in this work the chemical potentials in the leads
were coupled to dressed excitations rather than free electrons and the calculation was restricted to a subset of the excitations .
in contrast , a different approach was recently put forward by metha and andrei  @xcite to treat the interacting resonant level model .
they constructed a new set of scattering states of bethe - ansatz form which share the quantum numbers of free electrons in the incoming channel , hence allowing the application of a finite voltage in the usual manner .
however , questions concerning the existence of these scattering states and issues related to the used regularization scheme of the theory remain open .
nevertheless , the interacting resonant level model has become one of the benchmark systems in the study of non - equilibrium physics  @xcite .    despite the large number of studies of impurity models out of equilibrium only few results are known for the dynamical correlation functions .
the spin dynamics of a non - equilibrium quantum dot has been studied by using a majorana fermion representation  @xcite , which yields the qualitative low - frequency properties of the correlation functions . the transverse susceptibility in a kondo model
was studied by paaske et al .
@xcite using non - equilibrium perturbation theory together with a pseudo - fermion representation of the kondo spin .
they showed that the fourier transform of the transverse susceptibility possesses a peak if its frequency equals the value of the applied magnetic field , @xmath6 , and that the width of this peak is given by the transverse spin relaxation rate @xmath7 .
their derivation was , however , restricted to either the regime @xmath8 or @xmath9 , where @xmath0 denotes the applied voltage .
very recently , fritsch and kehrein  @xcite applied the flow - equation method to study the longitudinal correlation function as well as the magnetization and t - matrix in a kondo model in and out of equilibrium .
the numerical solution of the two - loop scaling equations allowed them to study the correlation function for all combinations of the parameters voltage , temperature , and magnetic field , provided the weak - coupling condition ( i.e. the presence of a large enough infrared cutoff ) was satisfied . in particular , these numerical solutions were used to compare the effects of an applied voltage and a finite temperature , revealing qualitative differences such as the appearance of kondo splitting in the non - equilibrium situation . in general , however , this method can not provide analytic expressions for the line shape .
in this article we will generalize the real - time renormalization group method in frequency space  @xcite to allow the calculation of dynamical correlation functions of arbitrary dot operators in systems describing spin and/or orbital fluctuations . in this setting the quantum dot is coupled to non - interacting leads which are held at different chemical potentials .
the derived two - loop rg equations can be solved analytically in the weak - coupling regime .
explicitly , we calculate the longitudinal and transverse spin - spin correlation and response functions in a two - lead kondo model in a magnetic field @xmath1 up to order @xmath10 . here
@xmath11 denotes the effective coupling at the energy scale @xmath12 where the flow of the coupling constant is cut - off . in order to satisfy the weak - coupling condition @xmath13 either the applied voltage or the magnetic field have to be sufficiently large compared to the kondo temperature , @xmath14 , where @xmath3 is the scale where the system enters the strong - coupling regime .
we note that the applied formalism does not rely on a fermionic representation of the kondo spin but rather deals with its matrix representation in liouville space directly .
the longitudinal response function possesses a peak at the spin relaxation rate @xmath15 , which gets suppressed in a finite magnetic field .
interestingly , in the case of a strong magnetic field , @xmath16 where @xmath17 denotes the renormalized magnetic field ( see eq .   for the precise value ) , the longitudinal correlation and response function show `` kink - like '' structures at the frequencies @xmath18 , which were also observed using the flow - equation method  @xcite .
here we additionally provide the line shape close to these `` kinks '' and show that the real part of the response functions shows characteristic logarithmic features at @xmath18 .
furthermore we study the longitudinal and transverse fluctuation - dissipation ratios .
as expected these ratios show a revival of the fluctuation - dissipation theorem  @xcite provided the applied voltage is small compared to the magnetic field or the considered oscillation frequency .
this article is organized as follows . in the next two sections
we will define the general set up we want to study and define the used notations
. this will include the kondo model , the notion of liouville operators as well as the definition of the symmetrized correlation function and susceptibility . in sec .
[ sec : pe ] we will then derive perturbative expansions for the kernels needed to calculate these correlation functions . this will be performed in liouville space ; the expansion is done in powers of the exchange coupling between the dot and the reservoirs .
these perturbative expansions can be applied to any model describing spin and/or orbital fluctuations as well as correlation functions of arbitrary operators . in sec .
[ sec : genrg ] we will use these results to derive the rg equations for the kernels of pure dot operators . in the following section we will further specialize to the two - lead kondo model , where we use the explicit expressions of the liouville operator and the coupling between the dot and the reservoirs to derive analytic results for the effective kernels appearing in the correlation functions .
finally , these expressions for the kernels are used in secs .
[ sec : cf ] and  [ sec : tcf ] to calculate the longitudinal and transverse correlation and response functions .
the real - time renormalization group in frequency space was applied in ref .   to calculate various quantities including the spin relaxation and dephasing rates , the renormalized magnetic field , the magnetization and the current in the anisotropic kondo model in a finite magnetic field out of equilibrium . in this reference
all notations which we will use in the following were originally set up . in order to increase the readability of the present manuscript we will briefly recall the basic formulas and notations .
we consider a quantum dot with fixed charge in the coulomb blockade regime coupled to external reservoirs . as shown in detail in ref .  , a standard schrieffer - wolff transformation leads to a hamiltonian of the form @xmath19 where @xmath20 is the reservoir part
, @xmath21 characterizes the isolated quantum dot , and @xmath0 describes the coupling between reservoirs and quantum dot .
they are given explicitly by @xmath22 here , @xmath23 are the fermionic creation ( @xmath24 ) and annihilation ( @xmath25 ) operators in the reservoirs and @xmath26 is an index characterizing all quantum numbers of the reservoir states , which contains the reservoir index @xmath27 and the spin quantum number @xmath28 .
we measure the energy @xmath29 of the reservoir states relative to the chemical potential @xmath30 of reservoir @xmath31 .
the eigenstates and eigenenergies of the isolated quantum dot are denoted by @xmath32 and @xmath33 .
the interaction @xmath0 is quadratic in the reservoir field operators , which arises from second order processes of one electron hopping on and off the quantum dot coherently .
this keeps the charge fixed and allows only spin / orbital fluctuations .
the coupling vertex @xmath34 is an arbitrary operator acting on the dot states .
it is written in its most general form , depending on the quantum numbers and energies of the reservoir states in an arbitrary way . as explained in ref .
, the rg approach can be set up in its most convenient form if one assumes that the frequency dependence of the initial vertices is rather weak and varies on the scale of the band width @xmath35 of the reservoirs . for the model we have in mind , the isotropic spin-@xmath36 kondo model in a magnetic field ,
this is certainly satisfied .
therefore , we will assume this in the following and introduce below [ see eq .
( [ contraction ] ) ] a convenient cutoff function into the free reservoir green s functions .    to achieve a more compact notation for all indices
, we write @xmath37 and sum ( integrate ) implicitly over all indices ( frequencies ) . the interaction
is then written in the compact form @xmath38 @xmath39 denotes normal - ordering of the reservoir field operators , meaning that no contraction is allowed between reservoir field operators within the normal - ordering . within the normal - ordering of eq .
( [ spin_orbital ] ) , the field operators can be arranged in an arbitrary way ( up to a fermionic sign ) , therefore the coupling vertex can always be chosen such that antisymmetry holds : @xmath40 furthermore , due to the hermiticity of @xmath0 , the vertex has the property @xmath41 where @xmath42 .     and
@xmath43 involve exchange between the electron spins of the left / right reservoir and the local spin , @xmath44 transfers an electron from one reservoir to the other during the exchange process .
we assume that the kondo model was derived from an anderson impurity model via a schrieffer - wolff transformation , which implies the relation @xmath45 . ]
the specific model we want to study is the isotropic kondo model in an external magnetic field @xmath46 ( see fig .  [
fig : kondomodel ] ) . in this case
the above relations read explicitly @xmath47 -(j_{\alpha'\alpha})_0 s^i\sigma^i_{\sigma'\sigma } & \text{for $ \eta=-\eta'=-$ } \end{array } \right.\,,\end{aligned}\ ] ] where @xmath48 , @xmath49 is the @xmath50-component of the spin-@xmath36 operator of the quantum dot , @xmath51 is a pauli matrix , and @xmath52 are the initial exchange couplings
. we will be interested in the antiferromagnetic model here , i.e. , we assume @xmath53 initially .
if one derives the kondo model via a schrieffer - wolff transformation from an anderson impurity model ( see , e.g. , ref .  )
, one further finds @xmath54 although the general formalism and many of the following formulas are also valid for an arbitrary number of reservoirs , we will consider the case of two reservoirs only with chemical potentials given by @xmath55 where @xmath0 is the applied voltage which we assume to be positive , @xmath56 .
a contraction is defined with respect to a grand - canonical distribution of the reservoirs , given by @xmath57 @xmath58 is the fermi distribution function corresponding to temperature @xmath59 ( note that the chemical potential does not enter this formula since @xmath29 is measured relative to @xmath30 ) .
furthermore , @xmath60 is the @xmath61-function in compact notation .
furthermore , we have introduced the cutoff by the band width @xmath35 into the reservoir contraction via the density of states @xmath62    in order to calculate the dynamical spin - spin correlation functions we have to know the time evolution of the density matrix @xmath63 .
formally , this follows from the solution of the von neumann equation @xmath64 where @xmath65}_-,\ ] ] is the liouvillian acting on usual operators in hilbert space via the commutator .
the form of the hamiltonian yields a similar decomposition of the liouvillian , @xmath66 with @xmath67}_-$ ] , @xmath68}_-$ ] , and @xmath69}_-$ ] .
we would like to note that the concept of liouville space and superoperators have been used in various contexts , for example in quantum statistical mechanics  @xcite .
initially , we assume that the density matrix is a product of an arbitrary dot part @xmath70 and a grandcanonical distribution @xmath71 for the reservoirs , @xmath72 furthermore , we introduce the laplace transform @xmath73 where we will frequently use the notation @xmath74 .
the stationary density matrix is defined as @xmath75 which is understood in the sense @xmath76 for any local operator @xmath77 , and can be calculated using @xmath78 the existence of a stationary state was proven in ref .   using non - equilibrium pertubation theory to all orders as well as in ref .
using the rtrg - fs , which in particular clarified the generation of the relaxation and dephasing rates under the rg flow .
the reduced density matrix of the dot is obtained by tracing out the reservoir degrees of freedom @xmath79 where @xmath80 denotes the effective liouvillian of the quantum dot formally defined in below .
the stationary reduced density matrix can then be obtained similar to , @xmath81
the quantities of interest in this article are the two - point correlation function of two operators @xmath82 and @xmath83 as well as their dynamical susceptibility with respect to the steady state , @xmath84}_+\bigr\rangle_{st } ,    \qquad   \label{eq : defphi}\\    \chi_{ab}(t)\!&=&\!{\text{i}}\,\theta(t)\bigl\langle    { \left[a(t)_\mathrm{h},b(0)_\mathrm{h}\right]}_-\bigr\rangle_{st}.    \label{eq : defchi}\end{aligned}\ ] ] where @xmath85 here the trace is taken over the dot states as well as the reservoir degrees of freedom , @xmath86 .
the time - evolution of the operators in the heisenberg picture is given by @xmath87 instead of calculating and in real time we will study their respective fourier transforms @xmath88 where @xmath89 for @xmath90 ( @xmath91 ) .
the susceptibility admits the standard decomposition @xmath92 .    in order to calculate @xmath93 and @xmath94
we introduce the auxiliary correlation functions @xmath95}_\pm\bigr\rangle_{st }    \label{eq : defbarcomega}\ ] ] with @xmath96 .
its relations to the correlation functions are given by ( see app .  [
sec : crel ] ) @xmath97 the static susceptibility is related to the dynamical susceptibility via @xmath98 where @xmath99 denotes the magnetization .
some general properties of the correlation functions can be obtained by considering their spectral representations .
let @xmath100 be a complete set of basis states of the full hamiltonian @xmath101 , i.e. @xmath102 .
furthermore , the stationary density matrix @xmath103 satisfies @xmath104}_-=0 $ ] , thus the basis states @xmath105 can be chosen such that @xmath106 using this one easily verifies the spectral representations @xmath107 these relations imply @xmath108 as well as @xmath109 , @xmath110 , and @xmath111 . in equilibrium
the matrix elements of the density matrix are given by @xmath112 with the partition sum @xmath113 , which implies the well - known fluctuation - dissipation theorem  @xcite @xmath114 as well as @xmath115 .
in this section we derive a perturbative expansion in liouville space for the auxiliary correlation functions @xmath116 , which will serve as the starting point for the derivation of the rg equations below .
a similar perturbative expansion for the effective liouvillian of the quantum dot @xmath117 has been derived in refs .  .
we will generalize these results to @xmath116 while closely following the presentation of ref .  .
as starting point to set up the formalism we assume that the operators @xmath82 and @xmath118 admit a representation similar to , @xmath119 where we recall the short - hand notation @xmath120 and sum ( integrate ) implicitly over all indices ( frequencies ) .
we further assume the operators @xmath82 and @xmath118 to be bosonic which implies @xmath121 and @xmath122 to be even .
eventually we will be concerned with the correlation functions of the spin operators on the dot , i.e. @xmath123 . in this case
the operators do not couple dot and reservoir degrees of freedom and hence only the terms with @xmath124 are non - vanishing .
however , we will keep the general forms throughout this section , which for example include the case of current operators  @xcite where @xmath125 .    in order to set up the perturbative expansions in liouville space we define the operators @xmath126}_+,\quad    l_b^\pm={\text{i}}{\left[b,.\right]}_\pm.\ ] ] then using together with we obtain after some algebra @xmath127 in the next step we use and and furthermore perform the laplace transform @xmath128 in to obtain @xmath129 here
we have used @xmath96 to ensure convergence of the integral .
the limit @xmath130 has to be taken before @xmath131 in order to reach the stationary state .
the next step is to expand the expression in the interacting part @xmath132 of the liouvillian and to integrate out the reservoir part .
this procedure was outlined for the reduced density matrix @xmath133 in detail in ref .  ; we will generalize this to the case of the auxiliary correlation functions here .
first , using @xmath134 with @xmath135 , can be formally expanded in @xmath132 , @xmath136 .
\label{eq : cab3 }    \end{split}\ ] ] second , in order to integrate out the reservoir degrees of freedom we write @xmath132 in the form @xmath137 where we implicitly sum ( integrate ) over @xmath138 as well as @xmath139 .
@xmath140 is a quantum field superoperator in liouville space for the reservoirs , defined by ( @xmath141 is an arbitrary reservoir operator ) @xmath142 here @xmath143 serves as an auxiliary index which is similar to the keldysh index indicating whether the field operator is acting on the upper or the lower part of the keldysh contour .
@xmath144 is a superoperator acting in liouville space of the quantum dot , and is defined by ( @xmath141 is an arbitrary operator on the quantum dot ) @xmath145 in the same way we define @xmath146 where the dot superoperators @xmath147 and @xmath148 act on arbitrary dot operators @xmath141 as @xmath149 for @xmath150 or @xmath151 we define @xmath152}_+$ ] and @xmath153}_\pm$ ] , respectively . as we consider only bosonic operators @xmath82 and @xmath118 which change the number of fermions by an even integer , the sign - superoperator is given by @xmath154 this operator was introduced to compensate additional signs due to interchanges of fermionic reservoir field operators as explained in detail in ref .  .
inserting the representations ( [ coupling_product ] ) as well as and into and shifting all reservoir field superoperators @xmath155 to the right using @xmath156 where we have introduced the short - hand notation @xmath157 , one can show@xcite that each term of perturbation theory can be written as a product of a dot part and an average over a sequence of field superoperators of the reservoirs with respect to @xmath71 . evaluating the latter with the help of wick s theorem
, one can represent each term of the wick decomposition by a diagram ( see fig .
[ fig : diagram_example ] for an example ) describing a certain process contributing to the auxiliary correlation function @xmath116 .
each process consists of a sequence of interaction vertices @xmath158 between the dot and the reservoirs , and a free time propagation of the dot in between ( leading to resolvents in laplace space ) .
since the reservoirs have been integrated out , the vertices are connected by reservoir contractions ( the green lines in fig .  [ fig : diagram_example ] ) .
this means that the various diagrams represent terms for the effective time evolution of the dot in the presence of dissipative reservoirs .    .
the time direction is to the left .
each vertex @xmath159 is represented by two adjacent black dots indicating the two reservoir field operators associated with each vertex .
the vertices @xmath160 and @xmath161 are represented in the same way by red dots .
for this example we have chosen @xmath125 in as it would be the case for a current - current correlation function .
the black horizontal lines connecting the vertices denote the free time propagation of the quantum system , leading to the resolvents @xmath162 in laplace space .
the green lines are the reservoir contractions arising from the application of wick s theorem .
the vertical blue lines between the vertices are auxiliary lines to determine the energy argument @xmath163 of the resolvents . ]
each diagram for the auxiliary correlation function has the form @xmath164 where @xmath68}_-$ ] , @xmath165 indicates an interaction vertex , and @xmath166 is a contraction between the reservoir field superoperators , defined by @xmath167 we stress that only the initial reduced density matrix of the dot @xmath70 defined in appears in as we have already performed the trace over the reservoir degrees of freedom ( and hence @xmath71 ) to obtain the reservoir contractions @xmath168 . to factorize the wick decomposition
, a fermionic sign has to be assigned to each permutation of reservoir field superoperators , indicated by the sign factor @xmath169 in ( [ value_diagram ] ) . for each pair of vertices connected by two reservoir lines , a combinatorial factor @xmath170 occurs , leading to the prefactor @xmath171 in ( [ value_diagram ] ) .
the value of the frequencies @xmath163 in the resolvents between the interaction vertices is determined by the sum over all variables @xmath172 of those indices belonging to the reservoir lines which are crossed by a vertical line at the position of the resolvent ( the blue lines in fig .
[ fig : diagram_example ] ) .
thereby , the index of the left vertex has to be taken of the corresponding reservoir line .
for example , the diagram shown in fig .
[ fig : diagram_example ] is given by ( the obvious dependence on the keldysh indices has been omitted for simplicity , i.e. @xmath173 and @xmath174 ) @xmath175 where the resolvents are defined by @xmath176 with @xmath177 as well as @xmath178 as can be seen from the example ( [ diagram_example ] ) , each diagram consists of a sequence of irreducible blocks ( where a vertical line always cuts at least one reservoir line ) and free resolvents @xmath179 or @xmath180 in between .
now there are two possibilities : ( i ) the vertices @xmath160 and @xmath161 do not belong to the same block ( see fig .
[ fig : block1example ] for an example ) .
( ii ) the vertices @xmath160 and @xmath161 belong to the same block ( see fig .  [
fig : block2example ] for an example ) .            in the first case ( i )
one can formally resum those terms between the vertices @xmath160 and @xmath161 which are not connected to them similar to dyson equations with the result @xmath181 where @xmath182 here the kernel @xmath183 contains the sum over all irreducible diagrams , @xmath184 where the subindex @xmath185 indicates that only irreducible diagrams are allowed where any vertical line between the vertices cuts through at least one reservoir contraction .
we further introduce irreducible blocks @xmath186 as well as @xmath187 which are given as the sum over all irreducible diagrams containing the vertices @xmath160 and @xmath188 , @xmath189 obviously , in the case of spin operators , @xmath190 , the vertex @xmath160 has no external legs and hence @xmath191 .
in contrast , although the vertex @xmath188 does not possess external legs either , there exist irreducible diagrams containing @xmath188 and at least one vertex @xmath159 to the left and one to the right of @xmath188 .
if we now proceed by resumming the irreducible blocks right to ( and not connected to ) the vertex @xmath188 similar to and perform the limit @xmath192 using , we deduce that all terms of type ( i ) contribute to @xmath193 .
\label{eq : barcabterm1}\ ] ]    in the second case ( ii ) we introduce a kernel similar to and which contains all irreducible diagrams containing both vertices @xmath160 and @xmath161 , @xmath194 in the case of spin operators , @xmath190 , there exist no irreducible diagrams connecting @xmath160 and @xmath188 , hence @xmath195 in this case . now using again for the sum of the irreducible blocks right to ( and not connected to ) the vertex @xmath188 we deduce that all terms of type ( ii ) contribute to @xmath196 .
\label{eq : barcabterm2}\ ] ] hence , taking together ( i ) and ( ii ) we finally arrive at the main result of this section , @xmath197\nonumber\\ *    & & \!\ !
-{\text{i}}\,\mbox{tr}_s\bigl[\sigma_{ab}^\pm(\omega,{\text{i}}0+)\,\rho_s^{st}\bigr ] ,    \label{eq : barcresult}\end{aligned}\ ] ] where the kernels are defined by , , , and , respectively .
in addition we note that the diagrammatic series can be partially resummed by taking all closed sub - diagrams between two fixed vertices together which contain only contractions connecting vertices between the two fixed ones .
this has the effect that the resolvents in ( [ value_kernel ] ) , , , and , are replaced by @xmath198 ( and similar for @xmath199 ) , i.e. , the full effective liouville operator occurs in the denominator . in this formulation
the number of diagrams is reduced , i.e. diagrams containing closed sub - diagrams between two vertices are no longer allowed .
in this section we will set up the generic rg equations for the kernels @xmath186 , @xmath187 , and @xmath200 , for spin operators on the dot in a model with spin / orbital fluctuations .
hence we will assume the form for the coupling between the reservoirs and the quantum dot , but will keep the vertex @xmath158 arbitrary at this stage .
the derivation will require some relations between the initial liouvillian @xmath201 , the vertex @xmath188 and the effective dot liouvillian @xmath117 which can be explicitly checked for the kondo model to be studied in the next section but have to be assumed here .
these relations are , , , and .
the generic rg equations for the vertex @xmath159 and the effective liouvillian @xmath117 have been derived and solved in ref .
, we will quote these results without derivation when they are needed .
furthermore , we will restrict ourselves to the calculation of the dynamical spin - spin correlations only , i.e. we will assume @xmath190 in what follows .
this implies in particular , that the initial values of the vertices @xmath160 and @xmath188 defined in and do not possess any external lines , i.e. @xmath124 . explicitly , @xmath202}_+,\label{eq : avertexinitial}\\    \mathcal{b}_\pm&=&{\text{i}}{\left[b,.\right]}_\pm.\label{eq : bvertexinitial}\end{aligned}\ ] ]
we will see below that this form of the vertex @xmath160 is conserved under the rg flow .
in contrast , a new effective @xmath118-type vertex @xmath203 with two external lines will be generated .
the fact that the initial vertex @xmath160 has no external lines directly implies the final results for the kernels and , namely @xmath204 hence , in the following we have to consider the kernel only .
we note that the results of this section remain valid for any pure dot operators @xmath82 and @xmath118 , i.e. any operators with @xmath124 .
the rg procedure is divided into two steps . in the first step we will integrate out the symmetric part of the reservoir contractions @xmath205 .
the reason for this is as follows : the effective dot liouvillian @xmath80 can be diagonalized as @xmath206 where @xmath207 and @xmath208 denote the eigenvalues and corresponding projectors , respectively .
this diagonalization implies for the resolvents @xmath209 now there exists a zero eigenvalue @xmath210 whose eigenstate for @xmath211 corresponds to the stationary state .
the appearance of this zero eigenvalue can lead to infrared divergencies of the frequency integrations in the perturbative expansions for the vertex @xmath159 and the effective liouvillian @xmath117 as is elaborated on in detail in ref .  .
however , after the discrete rg step we can trivially sum over the keldysh indices @xmath212 and @xmath213 by introducing @xmath214 in the resulting rg equations only the symmetric vertex @xmath215 will appear .
this vertex has the important property @xmath216 which is independent of the model specifics .
hence , after the discrete rg step the zero eigenvalue can no longer appear in any resolvent standing left to @xmath215 ( i.e. in no resolvent except the one standing left to the vertex @xmath188 ) .
this resolves the problem of infrared divergent internal frequency integrations , as the remaining eigenvalues @xmath207 have a strictly negative imaginary part ( see and ) .
we would like to refer to ref .
for a general discussion of this topic .    in the second step
we introduce a cut - off @xmath217 into the reservoir contractions via the fermi function .
we then integrate out the reservoirs by sending @xmath218 , which results in a description of the system in terms of effective dot quantities like @xmath117 .
this second , continuous , rg step is further divided into two substeps , first we integrate out the reservoir degrees of freedom down to an energy scale @xmath219 , and second we complete the flow from @xmath219 down to @xmath220 .      in the first discrete rg step we integrate out the symmetric part @xmath221 of the fermi function in the contraction ( [ liouville_contraction ] ) .
the discrete rg step for the kernel ( [ value_kernel ] ) and the vertex @xmath159 has been performed in ref .  .
here we will derive the analog results for the kernel and the vertex @xmath203 .
this is achieved by decomposing the contraction ( [ liouville_contraction ] ) according to @xmath222,\end{aligned}\ ] ] with @xmath223 . using this decomposition in
, one finds that each diagram decomposes into a series of blocks which are irreducible with respect to the symmetric part @xmath224 ( i.e. , any vertical line hits at least one symmetric contraction ) and connected to each other by antisymmetric contractions @xmath225 .
the blocks which are irreducible with respect to @xmath224 can be formally resummed into an effective kernel @xmath226 and a newly generated effective vertex @xmath227 .
the lowest order diagrams are shown in fig .
[ fig : discrete_lg ] .
( left and middle diagram ) and the effective vertex @xmath227 ( right diagram ) when the symmetric part of the contraction is integrated out .
@xmath228 ( @xmath229 ) denotes the symmetric ( antisymmetric ) contraction @xmath224 ( @xmath225 ) . ]    using the diagrammatic rules together with ( [ contraction_decomposition ] ) and the convention ( [ e_definition ] ) , we obtain for the first two diagrams @xmath230 and for the third one ( including the interchange @xmath231 ) @xmath232 we use here the original perturbation series ( [ value_kernelb ] ) so that the unperturbed liouvillian @xmath201 occurs in the resolvents .
performing the frequency integrations and assuming the band - width to be large , we obtain @xmath233 where we have performed the sum over the keldysh indices and used .
we stress at this point that the kernel @xmath234 and the renormalized vertex @xmath235 without external lines are identical , @xmath236 nevertheless we will retain the distinction between the kernel and the vertex in the following , as the former appears in the final formulas for the correlation functions , whereas the latter appears in the diagrammatic expressions for the r.h.s . of the rg equations .
we note that the frequency dependence in is generated during the flow as shown below .    after integrating out the symmetric part of the fermi function in this way
, we obtain a new diagrammatic series for the kernel analog to ( [ value_kernelb ] ) .
the liouvillian and the vertices have to be replaced by the effective ones and the contractions between the effective vertices contain only the antisymmetric part @xmath225 . due to there occur no diagrams including the new vertex @xmath237 . furthermore , since the effective quantities have become energy dependent ( also the effective vertex @xmath238 becomes energy dependent in higher order perturbation theory ) , one has to replace @xmath239 ( with @xmath240 ) in ( [ value_kernelb ] ) . since the antisymmetric part of the contraction ( [ contraction_decomposition ] ) does not depend on the keldysh indices , only the effective vertex @xmath238 averaged over the keldysh indices occurs in the new perturbative series .      in the second continuous rg procedure we deal with the remaining antisymmetric part of the fermi distribution function ,
where in each infinitesimal step a small energy shell is integrated out . instead of integrating out the energies on the real axis
, it has turned out to be more efficient to integrate out the matsubara poles of the fermi distribution function on the imaginary axis @xcite . this is achieved by introducing a formal cutoff dependence into the antisymmetric part of the fermi distribution by @xmath241 where @xmath242 are the matsubara frequencies corresponding to the temperature of reservoir @xmath31 , and @xmath243 is a theta function smeared by temperature . for @xmath244 , ( [ fermi_cutoff ] )
yields the full antisymmetric part @xmath245 of the fermi distribution . in each rg step
, one reduces the cutoff @xmath217 by @xmath246 , and integrates out the infinitesimal part @xmath247 of the fermi distribution .
for example , the new effective liouvillian at scale @xmath248 @xmath249 and similarly the new effective vertices @xmath250 and @xmath251 as well as the kernel @xmath252 can be calculated technically in the same way as for the first discrete rg step .
the only difference is that an infinitesimal small part is integrated out , so that the rg diagrams contain only one contraction involving @xmath253 .
furthermore , since the diagrams have to be irreducible with respect to this part , this contraction must connect the first with the last vertex of the diagram . using this procedure the rg equations for the dot liouvillian @xmath254 and the vertex @xmath255
have been derived in ref .  . here
we will use this technique to obtain the rg equations for @xmath256 as well as @xmath257 .     in @xmath258 .
the slash indicates the contraction where the fermi function has to be replaced by @xmath259 . ]     up to @xmath260 . ]
the diagrams contributing to the rg equations for @xmath261 and @xmath262 are shown in fig .
[ fig : rgb12 ] and fig .
[ fig : rgsigmab ] , respectively . using the definition @xmath263 together with the convention
@xmath264 we obtain the following rg equations :    @xmath265 , \label{eq : bvertexcontrg}\ ] ]    and @xmath266 \nonumber\\ & & -\frac{d\gamma^\lambda_{1}}{d\lambda}\,\gamma_2^\lambda\,\gamma_3^\lambda\ , \bar{g}_{12}(\omega)\,\pi_{12}(\omega ) \bigl [ \mathcal{b}_\pm(\omega_{12}+\bar{\omega}_{12},\xi_{12}+\bar{\omega}_{12})\ , \pi_{12}(\xi)\ , \bar{g}_{\bar{2}3}(\xi_{12}+\bar{\omega}_{12})\nonumber\\ * & & \hspace{20 mm } + \bar{g}_{\bar{2}3}(\omega_{12}+\bar{\omega}_{12})\,\pi_{13}(\omega)\ , \mathcal{b}_\pm(\omega_{13}+\bar{\omega}_{13},\xi_{13}+\bar{\omega}_{13})\ , \bigr]\pi_{13}(\xi)\ , \bar{g}_{\bar{3}\bar{1}}(\xi_{13}+\bar{\omega}_{13 } ) .
\label{eq : bkernelcontrg}\end{aligned}\ ] ]    we recall here that the kernel @xmath267 and the vertex without external lines @xmath188 equal each other , see , which yields a closed set of rg equations .
we will further show in app .
[ sec : b12 ] that the two - loop diagrams for @xmath203 as well as the one - loop diagrams containing @xmath203 itself do not contribute at second order in the coupling constant and hence can be neglected on the r.h.s . of .
the initial conditions of the rg equations are given by ( [ eq : bkernela ] ) and ( [ eq : b12a ] ) .
since @xmath268 , the solution at @xmath220 provides the result for the kernel @xmath269 from which the correlation functions can be calculated via .    as the resolvents and the vertices on the r.h.s . of the rg equations are analytic functions in all frequencies @xmath270 in the upper half of the complex plane
, all frequency integrations can be calculated analytically .
the only poles contributing are the ones of the contractions and their derivatives given by with as well as @xmath271 here @xmath272 denotes the matsubara frequency @xmath273 which lies closest to the cutoff @xmath217 . after performing the integration we find@xcite that , due to the presence of the cutoff function @xmath274 , the r.h.s . of the rg equations gives a negligible contribution for @xmath275 .
therefore , we can start the rg at @xmath276 and omit the cutoff function @xmath277 ( the precise ratio between @xmath278 and @xmath35 is determined such that no linear terms in @xmath35 in the effective liouvillian are generated  @xcite ) . as a consequence , only the matsubara poles of the fermi function in the upper half plane contribute and all real frequencies are simply replaced by matsubara frequencies . from now on , we write the frequency dependence explicitly and define the analytic continuation of the liouvillian and the vertices in imaginary frequency space by @xmath279 where we keep the real and imaginary parts of the laplace variable @xmath74 and the external frequencies @xmath280 and @xmath281 separated from now on . furthermore ,
@xmath282 , @xmath283 correspond to matsubara frequencies and the compact indices @xmath284 and @xmath285 on the l.h.s .
do no longer contain the frequencies @xmath286 . with the definition @xmath287 the rg equations and in matsubara space
can be written as    @xmath288    and @xmath289\nonumber\\ * & & \hspace{-25mm}\qquad\times \,\pi(\xi_{13},\xi'+\lambda_{t_{\alpha_1}}\!\!+\omega_3)\ , \bar{g}_{\bar{3}\bar{1}}(\xi_{13},\xi'+\lambda_{t_{\alpha_1}}\!\!+\omega_3 ; -\omega_3,-\lambda_{t_{\alpha_1 } } ) .
\label{eq : bkernelmatsubara}\end{aligned}\ ] ]    in these equations we implicitly sum over all indices and matsubara frequencies on the r.h.s . of the rg equations which do not occur on the l.h.s .
only positive matsubara frequencies smaller than the cutoff @xmath217 are allowed and each sum has to be written as @xmath290 which reduces to an integral @xmath291 for zero temperature .    in the next two subsections we will solve the rg equations ( [ eq : bvertexmatsubara ] ) and ( [ eq : bkernelmatsubara ] ) analytically in the weak coupling regime up to @xmath258 .
weak coupling is defined by the condition that the renormalized vertices @xmath292 stay small compared to one throughout the rg flow , so that the expansion in powers of @xmath159 on the r.h.s . of the rg equations
is well defined .
this condition is fulfilled if the various cutoff scales occurring in the resolvents are much larger than the kondo temperature @xmath3 at which the vertices would diverge in the absence of any cutoff scales .      as is discussed in detail in refs .
there exists a characteristic energy scale @xmath294 where @xmath295 is the renormalized magnetic field . for @xmath296
the cutoff scales @xmath297 , @xmath298 and @xmath5 can be neglected in the rg equation for the vertex @xmath215 ( see below ) .
this leads to a reference solution @xmath299 which serves as the starting point for a systematic expansion in powers of the coupling constant @xmath300 , where @xmath301 ( see ( [ g_vertex_liouville ] ) together with ) .
this yields a perturbative solution of the rg equations in the regime @xmath296 .
these results serve as initial values for the flow in the second regime @xmath302 . here
the renormalization of the vertex @xmath303 is at least of order @xmath10 , where @xmath304 .
provided the weak - coupling condition @xmath305 is satisfied all quantities can be calculated perturbatively .
this fact crucially relies on the appearance of some relaxation / dephasing rate in the resolvents , which is guaranteed by ( as then in all resolvents standing left to a vertex @xmath215 the zero eigenvalue of the liouvillian can not contribute ) .
this analysis has been performed for the liouvillian and the current kernel in the anisotropic kondo model in ref .  .
we note that the zero eigenvalue may appear in the resolvent left to the vertex @xmath188 in
. we will show below that this does not lead to any problems in the calculation of the spin - spin correlation functions in the kondo model up to order @xmath10 .    regarding the appearance of the external frequency @xmath306 as one of the cut - off parameters in we
see from the perturbative expansion that @xmath307 does not appear in the resolvents and as a cut - off parameter for all vertices right to @xmath188 .
this fact , however , will only affect the results in the regime @xmath308 .
we will therefore use as unique cut - off for all vertices appearing in the derivation of the kernel @xmath309 . in sec .
[ sec : expl ] we will show that for the spin operator in the kondo model the difference yields a correction @xmath310 and can thus be neglected .
nevertheless we stress that all vertices appearing in the stationary reduced density matrix @xmath311 in do not possess @xmath307 as cut - off parameter .
we thus deduce that in order to stay in the perturbative regime we can not rely on the external frequency @xmath307 but have to require @xmath312 .
we finally note that temperature serves as a unique cutoff for all terms on the r.h.s . of the rg equations as for @xmath313 the matsubara sums
are reduced to one term and the cutoff @xmath314 becomes independent of @xmath217 .
this trivial cutoff is set to zero in the following , i.e. we will set @xmath315 .    after these preliminary remarks let us turn to the evaluation of the rg equations .
the one - loop rg equation for the vertex @xmath215 is at zero temperature given by  @xcite @xmath316 the rg equation for the reference solution @xmath317 is obtained by assuming @xmath217 to be much larger than any other term appearing in the resolvent , which gives @xmath318 .
\label{eq : poorman}\ ] ] the initial condition for @xmath319 at @xmath320 is the bare vertex @xmath321 defined in .
the leading order solution is proportional to the coupling constant @xmath322 .
we stress that the term on the r.h.s .
, which is @xmath323 , contributes to the change of the vertex at order @xmath324 .
this is a general feature of the rg above @xmath219 . in order to calculate the change of a quantity at order @xmath325
one has to analyze those terms @xmath326 , where @xmath327 is some energy scale , and @xmath328 on the r.h.s . of the corresponding rg equation .
in contrast , terms @xmath329 @xmath330 do not contribute to the change at order @xmath325 .
the rg equation for the vertex @xmath331 is given by with the replacement  @xcite @xmath332 .
using this leading order solution we can formally expand all quantities in powers of @xmath324 , i.e. @xmath333 here @xmath68}_-$ ] is the bare dot liouvillian .
we recall , which implies @xmath334 in all orders in @xmath324 .
furthermore , we have already indicated which terms will depend on the matsubara frequencies , external frequencies @xmath335 and @xmath336 and the laplace variable @xmath337 .
the vertex @xmath215 and the liouvillian were calculated in ref .  .
we will here state those results needed for the calculation of and .
the second order vertex @xmath338 is further decomposed as @xmath339 here the vertex @xmath340 is given by  @xcite @xmath341.\ ] ] for the kondo model the vertex @xmath342 turns out to have the same matrix structure as the leading order solution @xmath299 .
this implies that both can be put together by redefining @xmath343 , which amounts to a two - loop renormalization of the kondo temperature  @xcite .
furthermore , the vertex @xmath344 is generically given by @xmath345 the zeroth - order liouvillian is given by the initial condition @xmath68}_-$ ] with while the first - order liouvillian is further decomposed as @xmath346 where @xmath347 and @xmath348 do not depend on the laplace variable .
the second - order liouvillian was calculated in ref .
, however , we will not need it for the solution of ( [ eq : bvertexmatsubara ] ) and ( [ eq : bkernelmatsubara ] ) in the regime @xmath296 .
let us now turn to the calculation of and .
the zeroth - order term of the kernel is just given by the initial condition , i.e. @xmath349}_\pm.\ ] ]    for the derivation of an rg equation for @xmath350 we have to keep all terms @xmath351 on the r.h.s .
of ( [ eq : bkernelmatsubara ] ) .
taking the zero - temperature limit and keeping only the order @xmath352 in the resolvents we obtain @xmath353 we evaluate this integral by assuming @xmath354}_-=\kappa\,\mathcal{b}_\pm^{(0 ) } ,    \label{eq : defkappa}\ ] ] which has to be checked for the specific model at hand . in the kondo model
we will find @xmath355 for @xmath356 and @xmath357 for @xmath358 ( see sec .
[ sec : expl ] below ) . can be used to shift @xmath359 to the right and evaluate the remaining integral by a partial fraction expansion @xmath360\bar{g}_{\bar{2}\bar{1}}^{(1)},\quad
\label{eq : rgb1}\end{aligned}\ ] ] where we have defined @xmath361 the leading term in is extracted by treating the terms @xmath362 separately , @xmath363 where @xmath364 can be integrated by @xmath365 with @xmath366 which has the asymptotic behavior @xmath367 as @xmath368 . using in
we obtain the frequency independent result @xmath369 the initial condition is given by , i.e. @xmath370 . as mentioned above we identify this with the vertex in first order , @xmath371 .    using the result for @xmath372
we derive in app .
[ sec : b12 ] the vertex , which is given by    @xmath373 ,    \label{eq : b12result}\end{aligned}\ ] ]    where @xmath374 ( see sec .  [
sec : expl ] below ) .
we note that @xmath375 . for large @xmath217 and using the result
we can derive the rg equation for the kernel in second order ( see app .
[ sec : bkernel ] ) @xmath376\bar{g}_{\bar{2}\bar{1}}^{(1)}\nonumber\\ *    & & \hspace{-20mm}-\frac{{\text{i}}}{2\lambda}\bigl [    \bar{g}_{12}^{(2a_1)}\,\mathcal{b}_\pm^{(0)}\,\bar{g}_{\bar{2}\bar{1}}^{(1 ) }    + \bar{g}_{12}^{(1)}\,\mathcal{b}_\pm^{(0)}\,\bar{g}_{\bar{2}\bar{1}}^{(2a_1 ) }    \bigr]-\frac{1}{2\lambda}\,\bar{g}_{12}^{(1)}\,\mathcal{b}_\pm^{(1)}\ ,      \bar{g}_{\bar{2}\bar{1}}^{(1 ) }      + \frac{1}{2\lambda}\,\bar{g}_{12}^{(1)}\ ,      { \left[z^{(1)},\mathcal{b}_\pm^{(0)}\right]}_+\ ,      \bar{g}_{\bar{2}\bar{1}}^{(1 ) } ,    \label{eq : b2kernelrg}\end{aligned}\ ] ]    where @xmath377 .
the initial condition is given by .
at this point it turns out to be useful to decompose the second - order kernel as @xmath378 here @xmath379 is given by @xmath380\bar{g}_{\bar{2}\bar{1}}^{(1)},\label{eq : sigmab2aabovelc}\end{aligned}\ ] ] which satisfies the initial condition @xmath381 the remaining terms in satisfy @xmath382,\quad\label{eq : sigmab2babovelc}\\    & & \frac{d}{d\lambda}\sigma_b^{\pm,(2c)}=    -\frac{1}{2\lambda}\,\bar{g}_{12}^{(1)}\,\mathcal{b}_\pm^{(1)}\ ,      \bar{g}_{\bar{2}\bar{1}}^{(1)}\nonumber\\ *    & & \quad+\frac{1}{2\lambda}\,\bar{g}_{12}^{(1)}\ ,      { \left[z^{(1)},\mathcal{b}_\pm^{(0)}\right]}_+\ ,      \bar{g}_{\bar{2}\bar{1}}^{(1)},\label{eq : sigmab2cabovelc}\end{aligned}\ ] ] with initial condition ( according to ( [ eq : bkernela ] ) and ( [ eq : sigmab2a_initial ] ) ) @xmath383 the rg equations in the regime @xmath296 derived in this section for a model describing spin / orbital fluctuations will be specialized and solved for the isotropic kondo model in sec .
[ sec : kondolargel ] below . in the next section
we will first study the effect of the rg flow in the regime @xmath302 .      as explained in ref .
the rg above @xmath219 has resummed all leading and subleading logarithmic contributions in @xmath384 into the renormalized vertices . at @xmath385 ,
the bare coupling constant is replaced by a renormalized one @xmath11 , all logarithmic contributions are eliminated , and a simple power series in @xmath11 remains .
thus the rg equations can be solved perturbatively provided @xmath386 .
in addition , the liouvillian in the resolvents is replaced by the full effective liouvillian @xmath80 .    when calculating diagrams containing the resolvent an obvious problem is that the frequency dependence of the effective liouvillian @xmath80 is not known explicitly .
to circumnavigate this complication we use the following approximation for the resolvents , @xmath387 here the poles @xmath388 of the resolvent follow from the self - consistency equation @xmath389 the residue satisfy  @xcite @xmath390 and hence can be set to one in the following .
starting with the one - loop rg equation for the vertex @xmath215 we observe that the terms on the r.h.s .
are already of @xmath391 .
therefore , the renormalization of the leading - order vertex @xmath299 stops at @xmath385 and we have to use its value @xmath392 in all calculations from now on . note that we indicate the use of the coupling constant at @xmath219 by the additional superscript @xmath393 .
furthermore , as we are eventually interested in the kernel @xmath394 up to @xmath391 we deduce that the higher - order vertices @xmath340 and @xmath344 will not be needed below @xmath219 .    using the replacement @xmath395 in
we obtain up to @xmath391 @xmath396 where we have already used that the flow of the kernel in order @xmath324 also stops at @xmath219 .
the second - order kernel satisfies    @xmath397    where we have used the approximation and only kept terms @xmath398 . the initial condition is given by the solution at @xmath385 .
using the decomposition we obtain @xmath399 where the flow of @xmath400 below @xmath219 is governed by .
the initial value at @xmath385 is given by , which we can rewrite as @xmath401\nonumber\\ *    & & \hspace{-30mm}\times    \bar{g}^{(1)c}_{12}\,p_i(z_i)\,\mathcal{b}_\pm^{(0)}\,p_j(z_j)\ ,    \bar{g}^{(1)c}_{\bar{2}\bar{1}}.    \label{eq : sigmab2ainitiallambdac}\end{aligned}\ ] ] in doing so we have assumed @xmath402 and neglected terms of order @xmath403 .
we will show in app .
[ sec : ass ] that these assumptions are fulfilled for the kondo model in a magnetic field .
we will further show that @xmath404 for all pairs @xmath405 for which the last line in is non - zero .
when applying the results of this section to other models of spin / orbital fluctuations one has to ensure the validity of the assumptions made above . the initial value problem with
is readily solved using @xmath406 with @xmath407 the result for the kernel at @xmath220 is then obtained using @xmath408 :    @xmath409 .
\label{eq : sigmab2aresult}\end{aligned}\ ] ]    together with @xmath410 and @xmath411 determined in the rg procedure above @xmath219 this yields the final result for the kernel in second order in @xmath11 .
it is applicable to any operator @xmath118 which does not couple the dot and reservoir degrees of freedom , i.e. whose initial value satisfies @xmath151 in .
furthermore the calculations were done for a generic model describing spin or orbital fluctuations as the initial vertex @xmath159 was assumed to have two external legs .
the only assumptions we have made regarding the model specifics are the commutation relations and , eqs .   and   as well as the specific relation between the parameter @xmath412 and the poles @xmath388 and @xmath413 , which have to be determined from the self - consistency equation .
in the next section we will apply the results derived above to the spin - spin correlation functions in the isotropic kondo model .
in particular , we will show that the assumptions discussed above are justified in this model .
finally we note that similar results for the effective liouvillian have been derived in refs .  .
in this section we will specialize the generic results derived above to the case of the spin - spin correlation functions in the isotropic , antiferromagnetic kondo model in a magnetic field .
the hamiltonian was presented in sec .
[ sec : kondo ] ; in particular , the dot hamiltonian and the coupling to the leads as given in and .      the first step is to represent the initial vertex and hamiltonian in liouville space , @xmath414}_-=h_0\,(l^{+z}+l^{-z})=h_0\,l^h,\end{aligned}\ ] ] where the spin superoperators @xmath415 are defined by their action on an arbitrary operator @xmath82 on the dot hilbert space via @xmath416 an explicit matrix representation for the superoperators @xmath417 is given in app .
[ sec : appl ] , where also further superoperators are defined .
the leading - order vertex @xmath299 was derived in refs .  .
it can be parametrized for @xmath418 as @xmath419 where @xmath420 is the vector formed by the pauli matrices .
the vertex for @xmath421 is obtained using the antisymmetry @xmath422 . inserting into the rg equation and using the antisymmetry @xmath423 ,  , and
@xmath424 we obtain @xmath425 if we assume the form , i.e. @xmath426 with @xmath427 , we obtain the usual poor - man scaling equation @xmath428 with the solution @xmath429 eq .   explicitly shows that the term @xmath323 on the r.h.s .
contributes to the renormalization in order @xmath324 .
similarly , the renormalization of a quantity in order @xmath325 is determined by the terms @xmath430 as well as @xmath328 .
in contrast , a term @xmath431 with @xmath432 does not contribute at order @xmath325 , as can be seen from @xmath433 the vertex @xmath331 is given by @xmath434    beside the leading - order vertex and the zero - order liouvillian we also need explicit expressions for the vertex @xmath340 as well as the liouvillian in first order .
these are given by  @xcite @xmath435 where @xmath436 is the coupling of the leading - order vertex @xmath299 and the trace tr is taken in the reservoir indices , @xmath437 .
furthermore we have taken the scaling limit @xmath438 , @xmath439 such that the kondo temperature @xmath3 remains constant .    as the next step it is straightforward to derive the following results for the initial vertex @xmath359 : @xmath440 where @xmath228 takes the values @xmath441 for @xmath442 and we have set @xmath443 , @xmath444 .
we stress that in this way the spin operators are directly represented by their matrices in liouville space .
hence , we do not have to use a pseudo - fermion representation of the kondo spin .
the vertex @xmath160 defined in is just given by @xmath445 .
furthermore , we recall that the kernels in zeroth order are just given by @xmath191 and @xmath446 , respectively . now
the rg equation for @xmath350 reads @xmath447 where we have used and .
the additional factor of two is due to the implicit summation over @xmath448 .
hence the solution in the scaling limit is given by @xmath449 in second order we will here determine only the terms @xmath450 and @xmath451 as their flow is cut - off at @xmath385 .
the remaining term @xmath400 will be derived in the next section .
hence we have to solve the rg equation , which using and reads @xmath452 analogously we obtain for using , and @xmath453 the solutions read @xmath454 where we have already taken the scaling limit so that the contributions from the initial condition are negligible .
as we have already explained above the rg flow of the leading order solution @xmath299 stops at the scale @xmath219 which is given by the maximal one of the external parameters , @xmath455 the value of the vertex at @xmath219 is given by with @xmath456 where due to the couplings satisfy @xmath457 with @xmath458 and @xmath459 the asymmetry ratio is defined as @xmath460 . as already mentioned the flow of all vertices right to @xmath188 does not stop at @xmath219 as defined in but rather at @xmath461 .
this affects the result for the kernel @xmath309 only in the regime @xmath462 .
we will discuss the changes in this case separately at the end of this section .
we stress , however , that the flow of all vertices @xmath215 needed for the derivation of the stationary reduced density matrix @xmath311 is cut off by @xmath461 . hence in order to stay in the weak - coupling regime we need @xmath463 as the definition of the scale @xmath219 is to some extent arbitrary as long as it remains of the order of the external energy scales in the problem , it is necessary to study the effect of a redefinition @xmath464 with @xmath465 .
this will induce a redefinition of the coupling as @xmath466 as we will show below , the redefinition @xmath464 does not change the final results for the liouvillian or the kernel @xmath309 up to order @xmath10 .
the liouvillian up to @xmath391 can be parametrized as  @xcite @xmath467 the liouvillian can be diagonalized ( see ) using the eigenvalues ( @xmath74 ) @xmath468 and the projectors @xmath469 the eigenvalues  can now be used to determine the poles @xmath388 of the resolvent defined in . solving the self - consistency equation
one finds @xmath470 , @xmath471 , and @xmath472 , where the spin relaxation and dephasing rates and the renormalized magnetic field are given up to @xmath391 by    @xmath473    the higher order terms @xmath474 for the rates were obtained in ref .  .
the voltage was defined in and we always assume @xmath475 . as these rates are obtained from the liouvillian at @xmath476 and @xmath477 , respectively , the external frequency @xmath307 does not appear as a cut - off in the definition of @xmath219 .
furthermore , we have defined the auxiliary functions @xmath478 which arises from terms like @xmath479 by neglecting the imaginary part of @xmath413 , which is proportional to @xmath480 , in the prefactor and takeing the real and imaginary parts .
we note that as @xmath472 the renormalized magnetic field automatically appears in the logarithm .
we have therefore also kept @xmath5 in the linear prefactor .
the deviation of @xmath481 from @xmath482 is only important for @xmath483 and will be neglected otherwise .
furthermore , we have omitted the imaginary part of the laplace variable , @xmath29 , since for the correlation functions calculated below we only need the liouvillian on the real axis . for @xmath484 and using , the functions in the parametrization are given up to order @xmath10 by  @xcite    @xmath485    in the last line we have replaced the bare magnetic field by the renormalized field , such that the latter appears consistently in all functions . the change is of @xmath486 .
we further note that a redefinition @xmath464 yields the same result for the liouvillian with the replacement @xmath487 . naively , the linear terms in @xmath488 and @xmath489 yield additional contributions in order @xmath490 .
these are , however , exactly cancelled by terms appearing from the logarithms using @xmath491 , where @xmath492 is given by with @xmath464 .
the facts presented above allow us to justify the assumptions and made in sec .
[ sec : genericrgbelowlambdac ] above for the specific case of the kondo model ( see app .  [
sec : ass ] ) .
after the recall of the liouvillian in second order we will finally evaluate the kernel @xmath394 in the kondo model .
the flow of the kernel also stops at @xmath219 except for one term , namely @xmath400 .
hence we find @xmath493 where the first term is given by or depending on the operator studied , and the second , fourth , and fifth is obtained from , as well as  using the replacement @xmath494 .
for the evaluation of the remaining contribution @xmath495 from let us first consider the operator @xmath358 . using @xmath496 for @xmath497 and @xmath498 ( in all other cases the l.h.s . vanishes )
we find @xmath499 we note that in the zero eigenvalue of the effective liouvillian appears in the resolvent left to @xmath500 .
as we will show in the next section the resulting term does not contribute to the correlation functions .
for the evaluation of the kernel for the calculation of the susceptibility we use    @xmath501    0 , & \text{otherwise } , \end{array}\right .
\label{eq : gpb - pg}\ ] ]    which results in @xmath502\nonumber\\ *    & & \hspace{-30 mm }    -\frac{{\text{i}}}{2}\frac{j_rj_l}{\omega-\xi+{\text{i}}(\delta-\xi ' ) }    \biggl[\left(\mathcal{h}_2(\omega+v+\tilde{h } )      + \mathcal{h}_2(\omega - v+\tilde{h})-\mathcal{h}_2(\omega+v-\tilde{h } )      -\mathcal{h}_2(\omega - v-\tilde{h})\right.\nonumber\\ *      & & \left.\quad        -\mathcal{h}_2(\xi+v+\tilde{h})-\mathcal{h}_2(\xi - v+\tilde{h } )      + \mathcal{h}_2(\xi+v-\tilde{h})+\mathcal{h}_2(\xi - v-\tilde{h})\right )    \left(l^a-\tfrac{1}{2}\,l^c\right)\nonumber\\ *    & & -\frac{1}{2}\left(\mathcal{h}_2(\omega+v+\tilde{h } )      + \mathcal{h}_2(\omega - v+\tilde{h})+\mathcal{h}_2(\omega+v-\tilde{h } )      + \mathcal{h}_2(\omega - v-\tilde{h})\right.\nonumber\\ *      & & \left.\quad        -\mathcal{h}_2(\xi+v+\tilde{h})-\mathcal{h}_2(\xi - v+\tilde{h } )      -\mathcal{h}_2(\xi+v-\tilde{h})-\mathcal{h}_2(\xi - v-\tilde{h})\right)l^h    \biggr ] .
\label{eq : sigmab2aresult2}\end{aligned}\ ] ]    next we consider the case @xmath503 and start with the evalution of the kernel @xmath504 from by using @xmath505    \mp\frac{{\text{i}}}{4}\frac{\gamma^{3z}(z_1)}{\gamma^a(z_1)}l^5_\pm , &     i=1,j=-,\\[2 mm ]    0 , & \text{otherwise}. \end{array}\right.\end{aligned}\ ] ] hence the double sum in reduces to a sum over @xmath498 , where the two terms have opposite signs and otherwise equal each other up to the appearance of the rate @xmath15 in the second term .
as can be easily shown this sum vanishes in second order , i.e. @xmath506 similarly , the kernel @xmath507 is evaluated using @xmath508    \!\!-\frac{{\text{i}}}{4}l^4_\pm , & i=\pm , j=1,\\[2 mm ]    0 , & \text{otherwise } , \end{array}\right.\ ] ] which results in    @xmath509\nonumber\\ *    & & \hspace{-30 mm }    -\frac{{\text{i}}}{2}\frac{j_rj_l}{\omega-\xi-\tilde{h}+{\text{i}}(\delta-\xi ' ) }    \biggl[\left(\mathcal{h}_1(\omega+v)+\mathcal{h}_1(\omega - v )      -\mathcal{h}_2(\xi+v+\tilde{h})-\mathcal{h}_2(\xi - v+\tilde{h})\right )    l^5_\pm\nonumber\\ *    & & + \left(\mathcal{h}_2(\omega+v-\tilde{h})+\mathcal{h}_2(\omega - v-\tilde{h } )      -\mathcal{h}_1(\xi+v)-\mathcal{h}_1(\xi - v)\right )    l^4_\pm\biggr ] .
\label{eq : sigmab2aresult2pm}\end{aligned}\ ] ]    a redefinition @xmath464 yields the same result for the kernels @xmath309 with the replacement @xmath487 as can be easily shown in the same way as for the liouvillian .    finally , let us consider the kernel in first order in the regime @xmath510 .
starting from we obtain using @xmath511 ( @xmath512 ) @xmath513 here the second vertex still depends on @xmath217 . when integrating this from @xmath461 to @xmath514 we obtain a contribution @xmath310 .
hence we deduce @xmath515 a similar analysis shows that the difference between @xmath516 and @xmath517 is at least @xmath518 .
in this section we will use the results for the liouvillian and the kernel to calculate the correlation functions and .
we first calculate the auxiliary correlation functions using .
we start with the parametrization of the liouvillian , which implies the form for the resolvent in @xmath519 .
we further deduce from the previous section that the kernels @xmath520 admit the parametrizations @xmath521 where for example @xmath522 . we stress that in contrast to the parametrization of the liouvillian we have not introduced additional factors of @xmath523 here .
the stationary reduced density matrix has the form @xmath524 with @xmath525 .
the stationary reduced density matrix can be determined using . without magnetic field
one simply finds @xmath526 .
furthermore , the rates and are given by @xmath527 we now rewrite the resolvent using the projectors and use @xmath528&=&-\frac{1}{2}\,\delta_{i1 } ,    \quad\label{eq : bla1}\\    \mbox{tr}_s\bigl[\sigma_{s^z}(\omega)\,p_i(\omega)\ ,    \sigma_{s^z}^-(\omega,{\text{i}}0+)\,\rho_s^{st}\bigr]&=&\nonumber\\ *    & & \hspace{-20mm}{\text{i}}\frac{\pi}{4}(j_r+j_l)^2\,\delta_{i1},\end{aligned}\ ] ] where we have applied as well as .
we note in particular that the term @xmath529 does not contribute to .
this yields with and @xmath530 : @xmath531 since @xmath532 we can neglect @xmath533 in . on the other hand , the real part of @xmath534 in the denominator has to be kept as it becomes large compared to @xmath307 in the small - frequency limit .
hence we arrive at @xmath535 we note that the leading term of the susceptibility is of order @xmath10 .
the correlation functions are plotted in figs .
[ fig : plot1][fig : plot3 ] .
we observe very good agreement with the results obtained by fritsch and kehrein using the flow - equation method  @xcite .     in the symmetric kondo model ( @xmath536 ) for various values of the applied voltage @xmath0 . ]     in the symmetric kondo model ( @xmath536 ) for various values of the applied voltage @xmath0 . ]     for @xmath537 and various values of the asymmetry ratio @xmath460 . ]
let us further study the behavior of the correlation functions analytically
. for small values of the frequency @xmath538 has the lorentzian form @xmath539 where we have used @xmath540 .
this result for the small - frequency regime agrees with conclusions drawn from a mapping of the spin correlators to the one - particle green s function of majorana fermions  @xcite . on the other hand , in the limit of large frequencies we find @xmath541 in agreement with the flow - equation method  @xcite .
we note that the @xmath542 appearing in the correlation function have their origin in the resolvent @xmath543 , hence the external frequency @xmath307 serves as a cut - off parameter in @xmath219 .
this results in the logarithmic corrections at large frequencies .
the susceptibility in the limiting regimes reads @xmath544 the first result shows a dependence of the gradient at small @xmath307 on the asymmetry ratio @xmath460 , while the second result indicates the revival of the fluctuation - dissipation theorem for @xmath545 .
we note that the derivation of relies on the fact that the coupling constants @xmath11 appearing in the kernel @xmath546 are cut - off by the external frequency @xmath307 , as it was discussed at the end of sec .
[ sec : expl ] .
furthermore , the susceptibility @xmath547 has a maximum at @xmath548 , where it takes the value @xmath549 this behavior was also deduced using the flow - equation method  @xcite .     for various values of the asymmetry @xmath550 and applied voltage @xmath0 . in order to get a smooth behavior at @xmath551 we have kept the @xmath552 in the definition of @xmath553 in this region .
the dotted line is a guide to the eye . ]    in order to investigate the revival of the fluctuation - dissipation theorem we introduce the longitudinal fluctuation - dissipation ratio  @xcite @xmath554 which is in equilibrium simply given by @xmath555 ( @xmath556 ) . using our results and we obtain @xmath557 which is plotted in fig .
[ fig : plot4 ] .
we find @xmath558 , i.e. the equilibrium result , whereas for small frequencies we get @xmath559 in agreement with refs .  .
we note that @xmath560 increases with increasing asymmetry @xmath550 as the coupling of the voltage to the dot becomes less effective .       in the symmetric kondo model ( @xmath536 ) for @xmath537 and various values of the applied magnetic field @xmath1 . for @xmath562
we already have @xmath16 , which implies @xmath563 in order @xmath10 ( see sec .  [
sec : longlargeh ] ) . ]    in the presence of an external magnetic field the stationary reduced density matrix is given by @xmath564 , @xmath565 , @xmath566 , with the magnetization in leading order  @xcite ( see also refs .  ) @xmath567 to evaluate the correlation function we now use @xmath568\\    & \qquad\qquad=    \left\{\begin{array}{ll}m\frac{\gamma^{3z}(\omega)}{\gamma^a(\omega)},&i=0,\\        -\frac{1}{2}-m\frac{\gamma^{3z}(\omega)}{\gamma^a(\omega)},&i=1,\\        0 , & i=\pm ,      \end{array}\right .
\end{split}\ ] ] where the term @xmath529 again does not contribute . from this a straightforward calculation using yields the correlation function up to @xmath391 @xmath569 where the zero - frequency @xmath61-peak does not appear because of our definition . the suppression of the correlation function by the finite magnetic field is shown in fig .
[ fig : plot5 ] , which agrees very well with similar plots obtained using the flow - equation method  @xcite . in the zero - frequency limit
we find @xmath570 while the leading term @xmath518 in the large frequency regime is given by ( including the logarithmic corrections in the coupling constants ) .
furthermore we observe a weak feature at @xmath571 which has for @xmath572 the line shape @xmath573.\quad    \label{eq : featurev - h}\end{aligned}\ ] ] similar features appear at @xmath574 .     in the symmetric kondo model ( @xmath536 ) for @xmath537 and various values of the applied magnetic field @xmath1 .
for @xmath575 we have @xmath16 , which implies @xmath576 in order @xmath10 ( see sec .  [
sec : longlargeh ] ) . ]    for the calculation of the susceptibility we need @xmath577\\    & \qquad\qquad=\left[\frac{{\text{i}}}{2}\,\gamma_{s^z}^{-,3z}(\omega )      + { \text{i}}\,m\,\gamma_{s^z}^{-,a}(\omega)\right]\delta_{i1 } ,    \end{split}\ ] ] which directly yields @xmath578 ,      \quad\label{eq : chipresultfiniteh}\\    \chi_{s^zs^z}''(\omega)&= &    \left[\frac{\pi}{4}(j_r+j_l)^2+m\,re\,\gamma_{s^z}^{-,a}(\omega)\right ]    \nonumber\\ *    & & \qquad\times\frac{\omega}{\omega^2+re\,\gamma^a(\omega)^2}.    \label{eq : chippresultfiniteh}\end{aligned}\ ] ] in we have kept the terms in the second line , which are of @xmath579 , as due to @xmath580 ( @xmath581 ) they become dominant in the small - frequency limit . for larger frequencies
these terms have to be neglected .
thus the static susceptibility is given in leading order by @xmath582 in agreement with the literature  @xcite .
for larger frequencies the real part of the susceptibility possesses logarithmic features at @xmath583 due to the term @xmath584 .
for example , @xmath585 where the terms represented by the dots do not contain any logarithmic features at @xmath586 . the imaginary part of the susceptibility is plotted in fig .  [ fig : plot6 ] .
it has a finite gradient at @xmath587 given by @xmath588 as well as a maximum at @xmath548 , where it takes the value @xmath589 in the large - frequency limit @xmath590 coincides with the correlation function .
furthermore , the imaginary part of the susceptibility has features at @xmath583 which have their origin in the function @xmath591 contained in @xmath592 and hence have a line shape similar to .     for @xmath593 and different values of the asymmetry @xmath550 and magnetic field @xmath1 .
in order to get a smooth behavior at @xmath594 we have kept the @xmath552 in the definition of @xmath553 . for @xmath16 the fluctuation - dissipation ratio is independently of @xmath307 given by @xmath595 .
the dotted line is a guide to the eye . ]
the fluctuation - dissipation ratio @xmath596 defined in reads in the presence of a magnetic field @xmath597 which is plotted in fig .
[ fig : plot7 ] .
larger values of the magnetic field push the system closer to its equilibrium behavior , as only those lead electrons in the energy interval @xmath598 can couple to the dot and thus induce the non - equilibrium behavior .
we note , however , that the equilibrium result @xmath595 is only reached for @xmath599 .
furthermore , increasing the asymmetry @xmath550 drives the system towards the equilibrium situation as the coupling of the voltage to the dot becomes less effective .
this effect is suppressed by increasing the magnetic field as overall less electrons couple to the dot . for small frequencies we obtain @xmath600       in the symmetric kondo model ( @xmath536 ) for @xmath602 and various values of the applied voltage @xmath0 .
the correlation function is given by @xmath603 .
for @xmath604 the susceptibility vanishes in order @xmath10 . the line shape close to @xmath586
is given by . ]     in the symmetric kondo model ( @xmath536 ) for @xmath602 and various values of the applied voltage @xmath0 .
we observe characteristic logarithmic features at @xmath18 . ]    in the case of a strong magnetic field , @xmath16 , the correlation functions up to quadratic order in the coupling are still given by , and , respectively , where the magnetization is simply @xmath605 .
one can easily show using @xmath606 that @xmath607 , which implies the equilibrium result @xmath608 .
furthermore , the correlation function vanishes identically in order @xmath10 for @xmath604 .
physically the kondo spin is in its ground state @xmath609 and the energy difference to the state @xmath610 due to the external magnetic field is given by @xmath5 .
hence one has to apply at least the frequency @xmath611 to obtain any response from the spin , where the energy @xmath0 is provided by the applied voltage .
this has to be contrasted with the result for the susceptibility in the equilibrium kondo model derived by garst et al .
they used a relation between the inelastic electron scattering and the correlation function to show that the susceptibility in equilibrium has the small - frequency behavior @xmath612 , i.e. it is non - zero for @xmath613 .
this linear behavior was also observed by costi and kieffer  @xcite as well as hewson  @xcite using a numerical renormalization group calculation . in analogy
, we expect the non - equilibrium correlation functions to be nonzero for @xmath604 in higher order in @xmath11 .
the consistent calculation of terms @xmath614 in the real - time rg procedure applied here would involve , however , 5-loop diagrams and is hence beyond the scope of this work .
the correlation function in the regime @xmath16 is plotted in fig .
[ fig : plot8 ] .
we find excellent agreement with numerical results recently obtained by fritsch and kehrein using the flow - equation method  @xcite .
in particular , we observe a splitting of the sharp edge at @xmath586 due to the applied voltage , which leads to characteristic features at @xmath18 .
using our result we can derive analytic expressions for the line shape close to these frequencies .
for example , at @xmath4 we find @xmath615 the first term shows that the gradient of @xmath547 will become negative for @xmath613 if the applied voltage is large enough , i.e. @xmath616 . in the vicinity of @xmath617
the correlation function shows similar kink - like behavior .
the physical origin of these kinks lies in the fact that at each of the energies @xmath18 a new process sets in , which involves a spin - flip on the dot costing the zeeman energy @xmath5 as well as the virtual hopping of an electron on and off the dot gaining or costing the energy @xmath618 , @xmath619 , or @xmath0 , respectively .
the real part @xmath620 of the susceptibility shows logarihmic features at @xmath18 as is shown in fig .
[ fig : plot8a].we stress that the splitting of the sharp edge at @xmath586 is a true non - equilibrium effect .
finally let us discuss the transverse correlation functions in the presence of a magnetic field .
we note that by virtue of we can restrict ourselves to the @xmath621-correlations .
the corresponding kernels up to second order in @xmath11 were calculated in , , , ,  as well as .
we will first discuss the susceptibility and present the results for the correlation function afterwards .     in the symmetric kondo model ( @xmath536 ) for @xmath602 and various values of the applied voltage @xmath0 .
]     for @xmath622 , @xmath602 , and various values of the asymmetry ratio @xmath460 .
the result is invariant under @xmath623 . ]    in order to derive the susceptibility we start with the parametrization @xmath624 where for example @xmath625 .
note that we already indicated that the explicit frequency dependence in order @xmath10 appears in @xmath626 exclusively .
( there is of course an implicit frequency dependence of @xmath627 and @xmath628 through @xmath11 . ) now using @xmath629\nonumber\\*[1 mm ]    & & \qquad={\text{i}}\left(\gamma_{s^+}^{-,3}+\bigl(\gamma_{s^+}^{-,2 }      + 2\gamma_{s^+}^{-,5}(\omega)\bigr)m\right)\delta_{i+}\qquad\end{aligned}\ ] ] we obtain @xmath630 where we have introduced the short - hand notation @xmath631 . for @xmath632
we find @xmath633 .
the transverse susceptibility has a peak at the solution of @xmath634 which is up to first order solved by @xmath635 in a finite magnetic field the spin on the dot will be in its ground state @xmath609 .
the energy difference to the excited state @xmath610 is given by @xmath5 , leading to an enhanced response of the system at this frequency .
at the peak the imaginary part of the susceptibility takes the value @xmath636 as is shown in figs .
[ fig : plot9 ] and  [ fig : plot10 ] .
the peak is suppressed by increasing the voltage , since this reduces the probability for the kondo spin to be in its ground state . on the other hand , for a fixed value of the voltage
the peak increases with increasing asymmetry ratio as the coupling of the voltage to the dot becomes less effective .
the width of the peak is up to order @xmath10 given by @xmath637 with the limiting cases @xmath638 in the equilibrium limit , @xmath639 , this corresponds to the result obtained in ref .  .
the real part of the transverse susceptibility possesses logarithmic features similar to at @xmath18 .    using a pseudo - fermion representation of the kondo spin together with non - equilibrium perturbation theory paaske et al .
@xcite previously obtained the transverse susceptibility . in order to compare these results
to we first make the approximations @xmath640 and @xmath641 . in the limit @xmath642
we then obtain using @xmath643 @xmath644 in the regime @xmath645 we have @xmath646 , thus we can neglect @xmath628 as well as @xmath626 in the numerator in , which results in @xmath647 these approximations agree with the results obtained in ref .
( we have to replace @xmath648 due to a different sign in the definition of the bare dot hamiltonian @xmath21 ) .
furthermore , in the regime @xmath649 we can use @xmath650 where we have replaced @xmath651 in the real part of @xmath626 , to obtain @xmath652 this confirms a conjecture by paaske et al .
we would like to stress , however , that our result goes beyond the approximation .    in analogy to the susceptibility one finds for the correlation function @xmath653 where we have neglected all terms of order @xmath403 in the numerator .
this allows the calculation of the transverse fluctuation - dissipation ratio @xmath654 which is plotted in fig .
[ fig : plot11 ] . for negative frequencies
@xmath655 the fluctuation - dissipation ratio takes the value @xmath656 , whereas for frequencies @xmath657 we find @xmath658 , thus recovering the equilibrium situation in these limits .
as for the longitudinal fluctuation - dissipation ration we observe that increasing the magnetic field or the asymmetry ratio @xmath550 drives the system towards the equilibrium situation .     for @xmath593 and various values of the asymmetry ration @xmath550 and the applied magnetic field @xmath1 .
the plot shows good agreement with similar results obtained in ref .  .
the dotted line is a guide to the eye . ]
in this article we have generalized the real - time renormalization group method in frequency space to allow the calculation of dynamical correlation functions of arbitrary dot operators in systems describing spin and/or orbital fluctuations .
we applied this to the two - lead kondo model in a magnetic field , where we calcualted the longitudinal and transverse spin - spin correlation and response functions up to second order in the exchange coupling .
we wish to stress that within this formalism the kondo spin is directly represented by matrices in liouville space , hence there is no need to apply a pseudo - fermion representation . specifically , we derived the two - loop rg equations for the dot operators and solved them analytically up to order @xmath10 in the weak - coupling regime .
here @xmath11 denotes the effective coupling at the energy scale @xmath12 which has to satisfy @xmath14 .
our results show several features attributed to the non - equilibrium situation , e.g. the splitting of the edge at @xmath586 of the longitudinal correlation function in a strong magnetic field or the suppression of the peak in the transverse susceptibility by a finite applied voltage .
furthermore , we find very good agreement with results for the longitudinal correlation function recently obtained by fritsch and kehrein using the flow - equation method  @xcite .
a particular advantage of our approach is the possibility to obtain analytic expressions for all correlation functions in the weak coupling limit .
we have calculated the spin - spin correlation functions for the nonequilibrium kondo model in the weak coupling regime @xmath14 .
the regime of strong coupling , @xmath659 , is still an open problem . in this case , the exchange couplings @xmath11 become of order @xmath660 and a controlled truncation of the rg equations is no longer possible . within the present rtrg - fs method it was shown in ref .   that the relaxation / dephasing rates saturate to the kondo temperature in the strong coupling regime .
as an effect the coupling constants do not diverge as in poor man scaling methods but remain finite .
however , the numerical solution of the rg equations in lowest order showed an instability against an exponentionally small change in the initial condition for the relaxation / dephasing rates .
although it was possible to find excellent agreement for the temperature dependence of the linear conductance with nrg calculations , it was necessary to fine tune the initial condition for the rates . therefore , up to now , it is not yet clear whether a controlled solution of the strong coupling regime is possible by using rtrg - fs .
the nonequilibrium kondo model describes the spin fluctuation ( or coulomb blockade ) regime of the more general nonequilibrium anderson impurity model ( for a systematic derivation of the kondo model from the anderson model using a schrieffer - wolff transformation , see e.g. ref .  ) . in this model
the single - particle spectral function is of most interest , which has recently been studied within nrg in a scattering wave basis @xcite .
the present rtrg - fs method can also be applied to this model but , as explained in detail in ref .
, in the charge fluctuation regime it is not yet clear whether a well - defined weak coupling regime exists at zero temperature . at resonance
( i.e. when the renormalized single - particle level is identical to one of the chemical potentials of the leads ) there is no energy scale except the broadening of the level itself and the expansion parameter is of order @xmath660 .
nevertheless the results obtained in ref .
( using a previous version of the real - time rg method ) after an a priori uncontrolled truncation of the rg equations were in excellent agreement with the bethe ansatz solution  @xcite for the equilibrium occupation of the local level .
these and related topics are of high interest and will be the subject of forthcoming research .
we would like to thank benjamin doyon , markus garst , stefan kehrein , verena krting , mikhail pletyukhov , and frank reininghaus for valuable discussions .
this work was supported the dfg - forschergruppe 723 `` functional renormalization group in correlated fermion systems '' and the dfg - forschergruppe 912 `` coherence and relaxation properties of electron spins '' .
obviously the time - dependent correlation function can be written as @xmath661}_+\bigr\rangle_{st}- \bigl\langle a\bigr\rangle_{st}\bigl\langle b\bigr\rangle_{st}.\ ] ] now , applying the time - translational invariance @xmath662 for all @xmath663 fixed and finite as well as the relation @xmath664}_\pm\bigr\rangle_{st}^*=
\bigl\langle{\left[a(t)_\mathrm{h},b(0)_\mathrm{h}\right]}_\pm\bigr\rangle_{st},\ ] ] which can be verified by a straightforward calculation using @xmath83 , we obtain . in the same way
can be used to derive the relation @xmath665 .
in this appendix we will calculate the vertex @xmath667 . starting from ( [ eq : bvertexmatsubara ] )
we first take the zero - temperature limit . furthermore
, we can expand the resolvents up to @xmath668 as@xcite @xmath669 with @xmath670}_+.\ ] ] using this we obtain together with the expansion @xmath671 in ( [ eq : bvertexmatsubara ] )    @xmath672}_+\right ) \frac{1}{\xi_{12}+{\text{i}}\xi'+{\text{i}}\lambda+{\text{i}}\omega_1-l_s^{(0)}-\tilde{l}_s^{(1 ) } } \left(1-\frac{z^{(1)}}{2}\right)\bar{g}_{\bar{2}1 ' } -(1\leftrightarrow 1'),\end{aligned}\ ] ]    where we have omitted the arguments of the vertices @xmath215 for simplicity .
using @xmath673}_-=    \tilde{\kappa}\,\mathcal{b}_\pm^{(0)},\quad    { \left[l_s^{(0)}+\tilde{l}_s^{(1)},\mathcal{b}_\pm^{(1)}\right]}_-=    \tilde{\kappa}\,\mathcal{b}_\pm^{(1)},\quad    { \left[l_s^{(0)}+\tilde{l}_s^{(1)},{\left[z^{(1)},\mathcal{b}_\pm^{(0)}\right]}_+\right]}_-      = \tilde{\kappa}\,{\left[z^{(1)},\mathcal{b}_\pm^{(0)}\right]}_+ ,      \label{eq : commapp}\ ] ] with @xmath674 for @xmath356 and @xmath675 for @xmath358 ( see sec .  [
sec : expl ] ) we obtain after a partial fraction expansion @xmath676}_+\right)\nonumber\\ * & & \,\times \left [ \frac{1}{\omega_{12}-\tilde{\kappa}+{\text{i}}\delta+{\text{i}}\lambda+{\text{i}}\omega_1    -l_s^{(0)}-\tilde{l}_s^{(1)}}- \frac{1}{\xi_{12}+{\text{i}}\xi'+{\text{i}}\lambda+{\text{i}}\omega_1-l_s^{(0)}-\tilde{l}_s^{(1 ) } } \right]\left(1-\frac{z^{(1)}}{2}\right)\bar{g}_{\bar{2}1 ' } + ( 1\leftrightarrow 1').\qquad\end{aligned}\ ] ] we see that no terms @xmath677 can occur on the r.h.s . and , hence , in order to determine the rg equation for @xmath667 we have to replace the vertices @xmath215 by the leading order ones @xmath299 and omit all terms containing @xmath372 , @xmath348 and @xmath678 . this yields using @xmath374 as well as @xmath679 : @xmath680\bar{g}_{\bar{2}1'}^{(1)}+(1\leftrightarrow 1 ' ) .
\label{eq : b12rgapp}\end{aligned}\ ] ]    using the rg equation for the leading - order vertex we can integrate up to higher - order corrections and obtain the solution , which in particular satisfies the initial condition @xmath681 as given by from the discrete rg step .    . ]     which contain @xmath203 itself . ]
the two - loop diagrams for the vertex @xmath203 are shown in fig .
[ fig : twoloopb12 ] . as
each diagram contains three vertices @xmath682 and we are interested in the vertex up to second order , we have to extract the terms @xmath683 .
this is done by expanding the resolvents in lowest order @xmath684 where @xmath685 ( @xmath29 ) stands for either @xmath307 ( @xmath61 ) or @xmath686 ( @xmath687 ) , respectively . using this
we immediately see that the first two diagrams are proportional to @xmath688 the fourth and fifth to @xmath689 and the third and sixth to @xmath690 thus the two - loop diagrams behave as @xmath691 and hence do not contribute to the renormalization of the second - order vertex @xmath667 .
finally we have to study the one - loop diagrams which contain @xmath203 itself .
these diagrams are shown in fig .
[ fig : rgoneloopb12withb12 ] .
as we can easily see from , the leading - order result for @xmath203 behaves for large @xmath217 as @xmath692 now the one - loop diagrams contain an additional vertex @xmath215 as well as a resolvent @xmath693 . hence we deduce that the diagrams yield terms proportional to @xmath694 which do not contribute to the renormalization of the second - order vertex @xmath667 .
in this appendix we will derive the rg equation for the second - order kernel @xmath696 .
for this we have to evaluate all terms on the r.h.s . of .
we start with the first line . using ( [ eq : restildel ] ) and expanding @xmath215 in powers of @xmath324 yields    @xmath697}_+\right)\nonumber\\ *    & & \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\times    \frac{1}{\xi_{12}+{\text{i}}\xi'+{\text{i}}\lambda+{\text{i}}\omega_2-l_s^{(0)}-\tilde{l}_s^{(1 ) } }
\left(1-\frac{z^{(1)}}{2}\right)\bar{g}_{\bar{2}\bar{1}}^{(1 ) }    \label{eq : appc1}\\    & & + \int_0^\lambda d\omega_2\ ,    \bigl({\text{i}}\,\bar{g}_{12}^{(2a_1)}+    \bar{g}_{12}^{(2b)}(\omega,\delta;\lambda,\omega_2)\bigr)\ ,    \frac{1}{\omega_{12}+{\text{i}}\delta+{\text{i}}\lambda+{\text{i}}\omega_2-l_s^{(0)}}\ ,    \mathcal{b}_\pm^{(0)}\ ,    \frac{1}{\xi_{12}+{\text{i}}\xi'+{\text{i}}\lambda+{\text{i}}\omega_2-l_s^{(0)}}\ ,    \bar{g}_{\bar{2}\bar{1}}^{(1)}\label{eq : appc2}\\    & & + \int_0^\lambda d\omega_2\,\bar{g}_{12}^{(1)}\ ,    \frac{1}{\omega_{12}+{\text{i}}\delta+{\text{i}}\lambda+{\text{i}}\omega_2-l_s^{(0)}}\ ,    \mathcal{b}_\pm^{(0)}\ ,    \frac{1}{\xi_{12}+{\text{i}}\xi'+{\text{i}}\lambda+{\text{i}}\omega_2-l_s^{(0)}}\nonumber\\ *    & & \qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad\times    \bigl({\text{i}}\,\bar{g}_{\bar{2}\bar{1}}^{(2a_1)}+    \bar{g}_{\bar{2}\bar{1}}^{(2b)}(\xi_{12},\xi'+\lambda+\omega_2 ;    -\omega_2,-\lambda)\bigr),\label{eq : appc3}\end{aligned}\ ] ]    where we have already neglected the terms @xmath348 and @xmath678 in the second and third line as they lead only to higher - order corrections . using the commutators
the first line can be treated similarly to @xmath350 derived in sec .
[ sec : wkabove ] . in the term
@xmath698 we only keep the term containing @xmath364 ( the term @xmath699 was already used to calculate @xmath350 ) , while in the other two terms we have to extract the term @xmath677 .
thus we arrive at @xmath700\bar{g}_{\bar{2}\bar{1}}^{(1)}\\    & \qquad-\frac{1}{2\lambda}\,\bar{g}_{12}^{(1)}\,\mathcal{b}_\pm^{(1)}\ ,    \bar{g}_{\bar{2}\bar{1}}^{(1 ) }    + \frac{1}{2\lambda}\,\bar{g}_{12}^{(1)}\ ,    { \left[z^{(1)},\mathcal{b}_\pm^{(0)}\right]}_+\,\bar{g}_{\bar{2}\bar{1}}^{(1)}.    \end{split }    \label{eq : appbkernelresult1a}\ ] ] using the same steps for the terms containing @xmath340 in and one finds @xmath701 .
\label{eq : appbkernelresult1b}\ ] ] for the terms of and containing @xmath344 we use in the form @xmath702 when inserted into and the first terms do not depend on the integration variable @xmath703 .
the remaining integral can be done as usual by a partial fraction expansion .
this yields ( @xmath704 for @xmath356 and @xmath357 for @xmath358 ) @xmath705\ ,    \bar{g}_{\bar{2}\bar{1}}^{(1)}\nonumber\\ *    & & \hspace{-8mm}+   \frac{{\text{i}}}{\omega\!-\!\xi\!-\!\kappa\!+\!{\text{i}}(\delta\!-\!\xi')}\ ,   \bar{g}_{12}^{(1)}\,\left [      \mathcal{k}_\lambda(\omega_{12}+{\text{i}}\delta - l_s^{(0 ) } )      \mathcal{b}_\pm^{(0)}-\mathcal{b}_\pm^{(0)}\ ,      \mathcal{k}_\lambda(\xi_{12}+{\text{i}}\xi'-l_s^{(0)})\right]\ ,    \bar{g}^{(1)}_{\bar{2}3 }      \ln\frac{2\lambda-{\text{i}}(\xi_{13}+{\text{i}}\xi'-l_s^{(0)})}{\lambda }    \,\bar{g}^{(1)}_{\bar{3}\bar{1}}.\nonumber\end{aligned}\ ] ]    if we now expand @xmath706 and the logarithm for large @xmath217 , @xmath707 and @xmath708 , and keep only the terms proportional to @xmath709 we arrive at @xmath710    \bar{g}_{\bar{3}\bar{1}}^{(1)}.    \label{eq : bg2b1stterm}\ ] ] in contrast , the second terms of and do depend on the integration variable @xmath703 .
the evaluations is , however , straightforward .
we use a partial fraction expansion for the resolvents left and right to the vertex @xmath359 as well as @xmath711 this leads to integrals of the form @xmath712 where we are only interested in the @xmath683 terms . now using the asymmetry @xmath423 we find @xmath713    \bar{g}_{\bar{3}\bar{1}}^{(1)}.    \label{eq : bg2b2ndterm}\ ] ] which has to be combined with for the full result from the terms containing @xmath344 .
the second and third line of containing the vertex @xmath667 can be treated using the same steps as were used to evaluate the @xmath344-dependent parts of and .
the result reads after some tedious but straightforward algebra @xmath714    \bar{g}_{\bar{3}\bar{1}}^{(1)}. \label{eq : appbkernelresult2}\ ] ]    finally , to extract the leading term @xmath715 of the fourth line of one can simply replace the resolvents @xmath716 by @xmath717 .
this yields @xmath718    \bar{g}_{\bar{3}\bar{1}}^{(1)}.    \label{eq : appbkernelresult3}\ ] ]    hence , the result for the rg equation of the kernel @xmath696 is obtained by summing , , , , and , using @xmath365 with in as well as @xmath719    \bar{g}_{\bar{3}\bar{1}}^{(1)}=    \frac{d}{d\lambda}\,\bar{g}_{12}^{(1)}\,\mathcal{b}_\pm^{(0)}\ ,    \bar{g}_{\bar{2}\bar{1}}^{(1)}.\ ] ] this yields .
consider an operator @xmath82 acting on the dot hilbert space having matrix elements @xmath720 with respect to the basis @xmath721 .
if @xmath722 denotes the superoperator acting on dot operators via @xmath723}_\pm$ ] then for an arbitrary dot operator @xmath118 we have @xmath724 furthermore , we represent superoperators in the matrix representation @xmath725 if @xmath726 is the product of two superoperators , then @xmath727 and the matrix of @xmath77 is simply given by the matrix product of the matrices of @xmath728 and @xmath729 .
a basis for the operators in the liouville space of the kondo dot can be built up by the spin superoperators @xmath730 and @xmath731 defined in .
an explicit representation in the basis is provided by @xmath732 @xmath733 @xmath734 furthermore we define the operators @xmath735 as well as @xmath736 we note that @xmath737 as well as @xmath738 .    in the spin sector
we will use frequently @xmath739 as well as @xmath740 the liouville operators satisfy ( the sums are over @xmath741 while the index @xmath212 takes the values @xmath143 ) @xmath742}_+\,l^{2i}&= &    \delta_{p-}\,\frac{{\text{i}}}{2}\,\text{tr}\,j\,l^h.\label{eq : appe5}\end{aligned}\ ] ] +
the results for the liouvillian presented in sec .  [ sec : resultsbelowlc ] allow us to obtain @xmath743 where we have used @xmath470 , @xmath744 and @xmath745 .
furthermore , we can expand the zero - eigenvalue projector as @xmath746 where @xmath747 .
this directly yields . to prove the final statement , @xmath404 for all pairs @xmath405 for which @xmath748 is non - vanishing
, we note that for @xmath358 these pairs are given by ( see and ) @xmath749 for @xmath500 and @xmath750 for @xmath751 , respectively . on the other hand we have @xmath357 in both cases , which implies @xmath752 for @xmath358 .
the same analysis can be performed for @xmath356 where the relevant pairs are given by @xmath753 for @xmath500 and @xmath754 for @xmath751 .