a self - consistent description of electrons , magnons , phonons and other fermionic and bosonic degrees of freedom in solids is still a challenging problem for strongly correlated systems where the simplest approximations are usually not sufficient .
one of the classical examples is the theory of itinerant electron magnetism , where the d - electrons form both the correlated energy bands as well as the quasi - localized magnetic moments and the magnon spectrum @xcite . in this problem
the conventional random phase approximation ( rpa ) gives a wrong description of the finite temperature effects whereas the proper corrections to rpa can be successfully calculated only in a very close vicinity of the stoner instability point @xcite .
the main difficulty is the coexistence of localized ( atomic - like ) and itinerant features in the behavior of d - electrons in solids @xcite .
recent progress in the theory of correlated electron systems is related to the development of dynamical mean field theory ( dmft ) @xcite , which maps the problem of fermions with local interactions on a crystal lattice onto a self - consistent effective impurity problem .
the latter can be solved very efficiently within continuous time quantum monte carlo ( ct - qmc ) @xcite or exact diagonalization ( ed)@xcite schemes . as a result ,
both atomic - like features such as multiplet formation and band dispersion can be taken into account simultaneously @xcite .
the combination of the dmft approach with realistic electronic structure calculations @xcite , in particular , for transition metal systems @xcite opens a new way for the description of electronic degrees of freedom in solids .
a generalization of the dmft approach to non - local interactions and consistent description of collective bosonic excitations @xcite boost first realistic applications of a so - called gw+dmft approach @xcite .    here
we construct a theory , that is capable to describe systems with strongly correlated electrons interacting with collective bosonic modes .
these modes can appear explicitly in an original hamiltonian ( for instance , phonons ) or arise from electron - electron interactions ( like plasmons or magnons ) .
we assume that the bose excitations can not be considered as an ideal bose gas and boson - boson correlations are important . among numerous examples of such systems
are the cases of weak itinerant ferromagnets @xcite and , more generally , any materials near quantum phase transition points @xcite .
the feynmann diagrammatic technique became the most popular tool in the quantum many body theory as a compact and elegant representation of a perturbation expansion in powers of interaction @xcite .
the main advantage of the diagrammatic approach is an opportunity to sum up infinite series of relevant diagrams going beyond the formal perturbation expansion .
the mainstream paradigm is that summing up more sophisticated series expands the area of applicability . in some cases leading series of diagrams
can be separated and summed up with the use of small parameters which are not directly related to the interaction strength .
important examples are the gell - mann bruekner theory of the high - density electron gas @xcite , the gaseous approximations of low density fermi @xcite and bose @xcite systems , the @xmath0-expansion in theory of critical phenomena @xcite , and the @xmath1 expansion in the theory of magnetic impurities @xcite .
historically , the dmft was introduced as a leading order approximation in a formal parameter @xmath2 , where @xmath3 is the dimensionality of the space .
dmft is exact in the limit of @xmath4 , at the same time attempts to construct a regular expansion in @xmath2 turn out to be not practical @xcite .
nevertheless , dmft appears to be quite a good approximation for many three- and even two - dimensional systems @xcite , since it provides reasonable interpolations between non - interacting and atomic limits for any dimensions . in terms of the diagrammatic approach the dmft corresponds to a summation of all local in space contributions to the self - energy @xcite .
recently we developed the dual fermion scheme @xcite which introduces new variables in a path integral generating an alternative diagram technique such that the bare propagators ( free dual fermions ) correspond to dmft .
this means that the main local physics is taken into account by the change of variables while the new diagrams describe only non - local corrections to the self - energy .
interaction terms in the dual variables look more complicated than the original ones : the interactions become retarded and contain not only pair terms but also multi - fermion contributions . however , an advantage of dual fermion transformation is a faster convergence of the diagram series .
for example , a simple ladder summation gives quite accurate description for the single band hubbard model in strongly correlated regime @xcite .
thus , the dual fermion approach is an example of an alternative strategy in the quantum many - body theory , when a proper choice of variables allows one to consider a simpler series of diagrams .
another example of such strategy is the use of auxiliary fermions and bosons in the hubbard model which allows one to describe the subtle effect of the growth of the effective mass near metal - insulator transitions , already on the level of the mean - field @xcite .
the consideration of collective ( bosonic ) degrees of freedom in systems of interacting fermions within the conventional technique necessarily involves summations of sophisticated diagrams .
the dispersion law of bosonic excitations ( e.g. , magnons or plasmons ) is determined by the poles of the corresponding susceptibilities ( the two particle fermionic green s functions ) and the minimal approximations which give rise to these poles , even in the case of a weak interaction corresponding to the ladder summation ( random phase approximation , rpa @xcite ) . in order to take into account magnon - magnon interactions one needs to consider more complicated , parquet - like diagrams @xcite . for strongly correlated systems
one can consider the same procedure within the dual fermion approach . up to now this is only done at the level of ladder summations @xcite . as we pointed out , the above considerations of complicated diagrams can be replaced by a proper choice of variables in the path integral .
following this strategy we propose an extension of dual fermion scheme explicitly adding bosonic variables .
it is important to note that the present theory is not a direct extension of the dual - fermion approach @xcite .
we will see that , even for the hubbard - like models with on - site interactions only , it gives a different result because of the account of collective modes ( magnons ) .
a simultaneous description of fermionic and bosonic correlations on the equal footing was started within the extended dynamical mean field ( edmft ) framework @xcite .
the strong - coupling generalization of edmft @xcite gives a useful extension of the hubbard - i like perturbation expansion for lattice models @xcite .
edmft properly includes the physics of _ local _ fermionic and bosonic correlations . at the same time , a dispersion of collective excitations within edmft is determined by the non - local part of the interactions _ only_. in particular , to describe the dispersive magnons within edmft one needs to explicitly add direct heisenberg exchange interactions to the hamiltonian .
obviously , this is not sufficient since the _ indirect _ exchange interactions ( superexchange , double - exchange , rkky , etc . )
contribute as well to the the magnon dispersion .
moreover , typically , they give the main contribution .    in the present paper
we construct a formally exact diagram technique where the bare fermionic and bosonic propagators correspond to edmft , whereas the indirect exchange is given by diagrammatic corrections .
therefore , the direct and indirect interactions enter the magnon dispersion law on equal footing .
we consider a lattice fermionic problem with the intersite hopping and generic non - local interactions . for simplicity , we present only formulae for the single - band periodic lattice , but a generalization of all expressions to a multiorbital case is straightforward . the action is @xmath5+\sum_{r , r\neq 0 , \omega , \sigma } \varepsilon_r c^\dag_{r \omega \sigma } c_{r+r \omega \sigma } + \sum_{r , r\neq 0 , \omega } v_{r \omega } \rho_{r \omega}^ * \rho_{r+r \omega}\ ] ] here @xmath6 are grassmann variables , the indices @xmath7 run over the lattice sites and translations , @xmath8 stands for the wave - vector , @xmath9 is the imaginary time , @xmath10 and @xmath11 are , respectively , the fermionic and bosonic matsubara frequencies ( @xmath12 is integer and @xmath13 is the inverse temperature ) , and @xmath14 is the spin index .
prefactors @xmath15 and @xmath16 are implied for sums over frequencies and wavevectors , respectively ( @xmath17 is the number of lattice sites ) .
the `` bosonic variable '' @xmath18 in general is a vector quantity , whose components for single band models can be the charge - density or the spin - density operators : @xmath19 , where @xmath20 with @xmath21 are pauli matrices ( @xmath22 ) , and the quantity @xmath23 is introduced to satisfy the condition @xmath24 .
to simplify the notations , we skip the index @xmath25 throughout the paper , but one should bear in mind that relevant components of @xmath18 are to be accounted in a practical calculation . for example , if the physics is determined by charge fluctuations only ( as may be the case for the coulomb interaction or an electron - electron attraction mediated by phonons ) , only the density - density term @xmath26 is important in the interaction part @xmath27 of eq.([lattice ] ) .
the heisenberg isotropic exchange interaction corresponds to @xmath28 , etc .
an important case of the hubbard model corresponds to @xmath29 , where the spin degrees of freedom @xmath30 do not show up in the bare expression part of eq.([lattice ] ) , but will appear in the following formulation , since the low - energy behaviour of system is governed by the spin fluctuations .    within the path integral representation the atomic action @xmath31
is given by : @xmath32 where @xmath33 is the local coulomb interaction .
we are interested in the single - electron green s function @xmath34 and the bosonic correlator @xmath35 .
it is an important point that an account of the bosonic quantity is fundamental : we suppose that collective modes are relevant for the physics of the system .
the general idea of effective - medium approximations is to reduce the lattice problem ( [ lattice ] ) to the so - called impurity model , describing a single atom ( or small cluster ) in a gaussian environment . for a system with bosonic and fermionic degrees of freedom ,
a single - site impurity problem is described by the action @xmath36    obtaining a numerical or analytical solution for the impurity problem is supposed to be much simpler .
we later assume that its green s functions : @xmath37 and @xmath38 are known , as well as higher - order correlators . in practice
, they can be calculated by a ct - qmc scheme @xcite or by exact diagonalization if the dimension of the hilbert space is not too large .
normally @xmath39 is called the full polarizability
. we will call it the bosonic green s function , to emphasize the formal analogy between the equations for @xmath40 and @xmath39 in our formalism .    for the so - called local approximations ,
only the green s functions of the impurity model @xmath41 are taken into account . in a nutshell ,
the initial problem ( [ lattice ] ) is replaced by an auxiliary lattice with gaussian nodes ; those nodes are characterized by the green s functions @xmath41 .
such a consideration gives us the well - known expressions for lattice green s functions @xcite , which are denoted with calligraphic letters in this paper : @xmath42 where @xmath8 is the wave vector and fourier transforms of the hopping and interaction parameters are introduced in the obvious way .
one can see that the fermionic and bosonic self - energies of these expressions are @xmath8-independent , that is , local in space .
the hybridization functions @xmath43 and @xmath44 should be defined in such a way , that the impurity problem provides us the on - site ( local ) properties of the solution for the lattice problem .
such a requirement is physically quite reasonable .
also , it is closely related to feynman s variational approach in dual fermion variables @xcite .
most common are so - called static and dynamic mean - field theories . in the static approach ,
the hybridization functions @xmath45 and @xmath46 are chosen to be frequency independent .
their values are defined from the equalities for the _
static parts _ of green s functions : @xmath47    in the dynamical edmft approach @xmath48 and @xmath49 are adjusted so that local parts of the lattice green s functions coincide with the corresponding quantities of the impurity problem _ for all frequency ( or time ) arguments _ : @xmath50    it is worthwhile to note that the standard definition of the edmft @xcite introduces an auxiliary bosonic field called @xmath51 , and the second line of eq.([dmftcondition ] ) is re - expressed as a condition for the green s function @xmath52 of that bosonic field . at the edmft level
, the conditions for @xmath53 and @xmath52 are equivalent .
we discuss this issue in detail in appendix and show that similar ( although not the same ) bosonic field can be introduced in non - local approximations beyond the edmft .
some `` mixed '' local approximations can also make sense .
for instance , if the hopping is absent ( or negligible ) , it is natural to consider a scheme with @xmath54 and to adjust bosonic hybridization only .
contrary , @xmath55 corresponds to the `` purely fermionic '' dmft scheme @xcite .
it is easy to make a connection between the lattice and impurity actions :    @xmath56+\sum_{\omega , k , \sigma } \left ( \varepsilon_k - \delta_{\omega\sigma } \right ) c^\dag_{\omega k \sigma } c_{\omega k \sigma } + \sum_{\omega , k } \left ( v_{\omega k}-\lambda_{\omega } \right ) \rho_{\omega k } ^ * \rho_{\omega k}\ ] ]    further on , we use a general matrix form of the hubbard - stratonovich transformation to decouple the term @xmath57 with @xmath58 into the pair of dual fermionic grassmann - fields @xmath59 and the term @xmath60 with @xmath61 into the pair of dual bosonic complex - fields @xmath62 ( a matrix notation is used to shorten the expressions ) : @xmath63e^{c^\dag e c } = \int d\left[c^\dag , c \right ] \det\left(\alpha^{-1}_f e \alpha^{-1}_f \right ) \int d\left[f^\dag , f \right ] e^{-f^\dag \alpha_f e^{-1 } \alpha_f f + f^\dag \alpha_f c + c^\dag \alpha_f f } \\ \\
\int d\left[\rho^ * , \rho \right]e^{\rho^ * w \rho } = \int d\left[\rho^ * , \rho \right ] \det\left(\alpha_b w^{-1 } \alpha_b
\right ) \int d\left[\eta^ * , \eta \right ] e^{-\eta^ * \alpha_b w^{-1 } \alpha_b \eta + \eta^ * \alpha_b \rho + \rho^ * \alpha_b \eta } , \end{array}\end{aligned}\ ] ] where @xmath64 are arbitrary , but local in space scaling matrices .
note that whereas the decoupling of grassmanian variables is mathematically straightforward , more comments are needed on the second line of eq.([hs ] ) .
the first comment is that @xmath18 is hermitian , so that @xmath65 , and the integral @xmath66 $ ] is taken in fact over the variables from a half - space .
the dual variables inherit this property , @xmath67 .
the second comment is that the matrix @xmath68 is not positively defined , and integrals over @xmath69 $ ] should be re - defined to ensure convergence .
here we describe the corresponding procedure , following , in particular , refs .
@xcite .    for illustration purposes ,
we suppose @xmath70 to be just a real constant factor ; the generalization to a complex matrix @xmath70 is straightforward
. let us introduce a basis where @xmath71 is diagonal .
let us number the states in this basis with an index @xmath72 , and consider some matrix element @xmath73 .
for a positive real part of @xmath74 , the corresponding part of the action can be decoupled as follows : @xmath75 where we explicitly introduced the hermitian and anti - hermitian parts of @xmath76 .
one can see that for @xmath77 and @xmath78 being interpreted as real and imaginary parts of a complex number @xmath79 , this decoupling becomes an explicit form of the second line of eq .
( [ hs ] ) , written for a particular state @xmath72
. the integral limits for positive @xmath74 can be chosen along the real axis , from @xmath80 to @xmath81 for both @xmath77 and @xmath82 .
the key observation is that for @xmath74 having a _ negative _ real part , the decoupling ( [ hsplus ] ) still holds , up to the integration path along imaginary axis , from @xmath83 to @xmath84 for @xmath77 and in the opposite direction , that is @xmath84 to @xmath83 , for @xmath82 .
thus , the integration over @xmath85 $ ] can be considered as a symbolic notation for the integration over @xmath86 $ ] with @xmath87 substituted in the integrand .
the integration path over each component of vectors @xmath88 is explicitly defined for the basis where @xmath89 is diagonal : the integration goes along the real or imaginary axis , corresponding to the sign of the real part of the diagonal matrix element @xmath74 .
one can check that such a definition of @xmath90 $ ] still allows one to use the standard algebraic manipulations of path integrals .
in particular , to obtain the green s function of dual bosonic variables @xmath91 one uses the standsrd procedure of differentiation with respect to @xmath92 , as can be explicitly shown from eq .
( [ hsplus ] ) . on the other hand
, the statistical properties of dual variables can be unusual because of the imaginary axis integration path .
for instance , similarly to the fermionic dual green s function @xcite , @xmath93 is not necessary positive - defined and consequently the imaginary part of its fourier transform is not always negative .
let us continue with the derivation of the dual formalism .
the hubbard - stratonovich decoupling in the partition function with the action ( [ lattice ] ) gives ( in matrix form ) @xmath94e^{-\sum_{r}s_{imp}[c_{r}^{\dag } , c_{r}]+c^{\dag } ( \delta -\varepsilon ) c+\rho ^{\ast } ( \lambda -v)\rho } = \det ( \alpha _ { f}^{-1}(\delta -\varepsilon ) \alpha _ { f}^{-1})\det ( \alpha _ { b}(\lambda -v)^{-1}\alpha _
{ b})\times \\ \\ \times \int d[c^{\dag } c]\int d[f^{\dag } f]\int d[\eta ^{\ast } \eta ] e^{-\sum_{r}s_{imp}[c_{r}^{\dag } , c_{r}]+f^{\dag } \alpha _ { f}c+c^{\dag } \alpha _ { f}f+\rho ^{\ast } \alpha _ { b}\eta + \eta ^{\ast } \alpha _ { b}\rho -f^{\dag } \alpha _ { f}(\delta -\varepsilon ) ^{-1}\alpha _ { f}f-\eta ^{\ast } \alpha _ { b}(\lambda -v)^{-1}\alpha _ { b}\eta }
. \end{array}\end{aligned}\ ] ] integrating out the `` localized '' fermionic grassmann variables @xmath95 with @xmath96 for the fermionic and @xmath97 for the bosonic part , we arrive at the following action in dual variables @xmath98 \\ \\
k}=(g_{\omega}^{-1}+\delta _ { \omega }    -\varepsilon _ { k})^{-1}-g_{\omega } \\ \\
\tilde{\mathcal{x}}_{\omega k}=(\chi _ { \omega}^{-1}+\lambda _ { \omega } -v_{k } ) ^{-1}-\chi _ { \omega } \\ \\
\tilde{u}[\eta , f , f^{\dag } ] = \sum_{\omega \omega } \left ( \lambda _ { \omega \omega } \eta _ { \omega } ^{\ast } f_{\omega + \omega } ^{\dag } f_{\omega } + \lambda _ { \omega \omega } ^{\ast } \eta _ { \omega } f_{\omega } ^{\dag } \ f_{\omega + \omega } \right ) + \frac{1}{4}\sum_{\omega \omega ^{\prime } \omega } \gamma _ { \omega \omega ^{\prime } \omega } f_{\omega + \omega } ^{\dag } f_{\omega ^{\prime } -\omega } ^{\dag } f_{\omega   ^{\prime } } f_{\omega}+ ... , \end{array } \label{db_action}\ ] ] where @xmath99 is a full four - point fermionic vertex of the impurity model @xmath100 ( @xmath101 is the kronecker symbol ) , and @xmath102 is a `` mixed '' quantity @xmath103 the second term in the nominator typically equals zero as @xmath104 vanishes .
the values of @xmath99 and @xmath102 are related , since @xmath105 , @xmath106 it is important to note that for the gaussian impurity problem @xmath99 vanishes , but @xmath102 takes a finite value of @xmath107 .
the following exact relations between fermionic and bosonic green s functions of original and dual variables can be proven : @xmath108 the relation for the fermionic green s functions is derived in ref .
main idea is to substitute the right - hand side of eq.([smix ] ) in the formula @xmath109 .
the last term @xmath110 appears from the differentiation of the fermionic determinant .
the bosonic relation is derived quite similarly , but the signs appearing during the derivation should be discussed .
first , the expression for the bosonic green s function has a different sign , @xmath111 .
this gives factors of @xmath112 at @xmath113 and @xmath114 .
second , since the bosonic hubbard - stratonovich transformation contains an inverse of the determinant , the same different sign appears for the last term @xmath115 , so that the signs of all the terms are changed , compared to the fermionic relation between the real and dual fermionic green s functions .
therefore the entire resulting relationship between the real and dual bosonic green s functions obeys exactly the same form , similar to the results of the supersymmetric hubbard - stratonovich transformation @xcite .
it is useful to introduce the self - energy ( for fermions ) and polarization - operator ( for bosons ) corrections to the effective - medium approximation , defined as @xmath116 we remind that @xmath117 are given by the expressions ( [ mfg ] ) .    the self - energy and polarization operator of the dual ensemble
are given by the following expressions : @xmath118 being re - expressed in these quantities , the exact relations ( [ exact ] ) have a particularly simple form : @xmath119 which results in the following expression for the fermionic and bosonic green s functions : @xmath120 one can immediately see from these formulas that the effective - medium expressions ( [ mfg ] ) correspond to the gaussian approximation for the dual variables : @xmath121 ( and consequently @xmath122 ) vanish .
following the effective medium paradigm we introduce a generalized self - consistent condition for @xmath123 and @xmath124 similar to the eq.([dmftcondition ] ) , but for `` exact '' fermionic and bosonic green s functions : @xmath125
the simplest non - local physics can already be described in the edmft . in our formalism
this corresponds to a zero - th order approximation ( `` no diagrams '' ) @xmath126 and @xmath127 , and condition ( [ dmftcondition ] ) is satisfied .
we intend to preserve the advantages of the edmft , but also include the formation of collective modes and their interaction with the electrons . for a `` purely fermionic '' hubbard - like theory
this means an account of the ladder series in dual space @xcite .    for the present theory with the dual action of eq.([db_action ] ) , there are two basic green s function lines ( fig .
[ basic_diagram ] ) : a fermionic ( full - line with arrow ) and a bosonic ( wavy - line ) as well as two interaction vertices : a triangle ( corresponding to the dual boson - fermion interaction - @xmath128 ) and a square ( describing the dual - fermion interaction - @xmath129 ) . in this paper , we restrict ourselves to a polarization operator , given by the electron ladder ended by two triangles ( fig . [ fig_ladder ] ) .
equally , one can consider this diagram as just a two - lines bubble , with a renormalized vertex ( fig .
[ fig_pi ] ) .
@xmath130 here a symbolic notation is used , so that @xmath131 is a tensor , and @xmath132 is a vector with components corresponding to different @xmath133 , and the dots in the first line denote a scalar product vector - tensor - vector . the bare dual polarization operator ( an empty loop of the dual green s function )
@xmath134 is a tensor having only diagonal components .
later we also use bare polarization operators for the initial variables and for the impurity model , so that @xmath135 scalar polarizabilities , respectively denoted by @xmath136 , @xmath137 , and @xmath138 , can be found by summating over fermionic frequencies , so that , for example , @xmath139 .    the simplest approximation one can make is just to neglect the vertex part @xmath99 . in this case @xmath140 .
simple manipulations with the second lines of eqs .
( [ mfg ] ) , ( [ selfenergy ] ) and ( [ exact2 ] ) then give @xmath141 a similar non - local contribution to the fermionic self - energy can be introduced to describe the renormalization of the fermionic green s function @xmath142 .
the corresponding diagram contains one fermionic and one bosonic line ( fig .
[ fig_sig ] ) .
we will however concentrate on the analysis of bosonic quantities and for simplicity use the edmft result @xmath143 instead of the renormalized @xmath144 propagator .
this allows for a further simplification .
indeed , in the local approximation @xmath145 and @xmath146 differ in their local part only @xcite : @xmath147 further , since we suppose that @xmath148 and since @xmath24 , the bosonic green s function for the impurity problem becomes equal @xmath149 . substituting the expression ( [ ggprime2 ] ) into eq .
( [ drpa ] ) , we obtain @xmath150 for @xmath151 , this expression is clearly related to the random phase approximation ( rpa ) . to get the `` textbook '' rpa , one should completely neglect the interaction terms of the impurity problem , so that @xmath152 and @xmath153 . the scheme with an exactly calculated @xmath154
can be called rpa+dmft . in this case
the rpa is a simple case of fig .
[ fig_pi ] . to go beyond the rpa
, one should be careful . as we will discuss in sec .
[ conserv ] in general it is wrong to work with the renormalized green s functions but to neglect the vertex @xmath99 . because this may lead to a violation of the conservation laws .
a theory using the single - bubble approximation ( [ buble ] ) with @xmath155 approximated by @xmath102 ( that is , without putting @xmath102 to @xmath156 ) is essentially the gw+edmft approach @xcite .
there is only a slight difference in the self - consistency condition , as we explain in appendix .
it is also important that for the diagram calculation of the gw+edmft approach that the local part of the bosonic green s function is excluded _
ad hoc_. as one can see from eq .
( [ edmftggprime ] ) , for the edmft choice of the bose variables , it actually does not vanish exactly .
so , our consideration proposes a certain a modification of gw+edmft . anyhow
, such modification should not strongly affect the results of the gw+edmft .
the real advantage of our consideration is the possibility to go beyond the domain of applicability of the gw+edmft , i.e. in is the situation when the single - bubble approximation becomes incorrect .
it is considered in the next sections .
the diagrams shown in figs .
[ bold3 ] , [ fig_pi ] , [ fig_sig ] form the minimal consistent set beyond the rpa to describe the formation of collective modes with electron - boson coupling taken into account explicitly .
so far , we did not consider whether the constructed approximations conserve particle number , total spin , etc .
however , the conservation laws are particularly important for bosonic excitations . for example
, the charge ( particle number ) conservation resulting in the equation @xmath157 immediately leads to the requirement @xmath158 at any finite frequency ( provided that the current - current correlator @xmath159 is finite ) . for 3d systems with coulomb interaction the long - wavelength asymptotic behaviour for the density - density correlator should be @xmath160 , where @xmath161 is the plasma frequency .
the rpa for free electrons is proven to obey this property , since @xmath162 vanishes with @xmath163 , but this is not true for the renormalized green s function @xcite . indeed , following the standard proof @xcite one writes @xmath139 as @xmath164 for free electrons of dispersionless @xmath165 ( that is , in static mean - field ) the summation over @xmath133 produces a difference of occupation numbers , @xmath166 , in the nominator , that vanishes at @xmath167 .
however this reasoning clearly breaks down for the frequency - dependent @xmath168 in dmft approximation .
it was known for a long time that the neglection of vertex corrections together with the use of _ renormalized _ green s functions , violate charge and spin conservations @xcite .
we found that @xmath169 with the dmft green s function does not vanish at small @xmath170 and therefore it follows form eq .
( [ rpa2 ] ) that _ edmft+gw is not a valid description of plasmons_. all the more , _
plasmons do not appear in simple edmft _ , as one can easily observe from the second line of eq .
( [ mfg ] ) : the local quantity @xmath171 can not behave like @xmath172 .
this is the reason why the gw approach is normally used only for calculating the single - particle green s function , while any two - particle properties are investigated within the bethe - salpeter equations @xcite .
nevertheless , such a scheme can not describe self - consistent effects of collective excitations on the single - particle spectral function .    to construct a conserving approximation
, it is useful to consider the following reformulation of the theory
. let us return to the initial lattice action ( [ lattice ] ) and do the hubbard - stratonovich decoupling for the bosonic variables only .
it will give an analog of eq .
( [ smix ] ) with the action @xmath173+\rho^ * \alpha_b \eta + \eta^ * \alpha_b \rho-\eta^ * \alpha_b ( \lambda - v)^{-1 } \alpha_b \eta,\\ s^f=\sum_r s_{at}[c^\dag_r , c_r]+\sum_{r , r\neq 0 , \omega , \sigma } \varepsilon_r c^\dag_{r \omega \sigma }
c_{r+r \omega \sigma } + \sum_{r , \omega } \lambda_{\omega } \rho_{r \omega}^*\rho_{r \omega}. \end{array}\ ] ] physically , @xmath174 describes a lattice with the interaction being local in space but retarded .
now we formally integrate over the fermionic variables .
the gaussian approximation in dual bosons for such a theory corresponds to the expression @xmath175 where @xmath176 is an exact @xmath170-dependent two - particle correlator for the lattice problem @xmath174 .
formally , the calculation of @xmath177 corresponds to a summation of all parts of the diagrams with fermionic lines in series for eq.([db_action ] ) . comparing eqs .
( [ greenf ] ) and ( [ exact ] ) we obtain @xmath178    this way , we need a conserving approximation for @xmath176 .
such an approximation is actually known within the framework of dmft . because it is a conserving theory in the baym - kadanoff sense @xcite , with the self energy being just the sum of all local diagrams .
in such a theory , the conserving approximation for two - particle quantities can be obtained as a green s function variational derivative with respect to the dispersion law , @xmath179 an important property of the baym s criterion of consistency is that the variation should be considered self - consistently , so that the @xmath43 and @xmath40 of dmft vary with @xmath0 to preserve the self - consistency condition ( [ dmftcondition ] ) .
one can also explicitly see that such a self - consistent variation preserves the charge conservation .
indeed , the charge conservation follows from gauge invariance under the transformation @xmath180 with an arbitrary dependence of the infinitesimal phase shift @xmath181 on time and space arguments . in an exact theory
, such a transformation just results in similar phase rotations for averages , for example @xmath182 .
a conservative approximation must inherit this property .
dmft is entirely defined by the self - consistency condition ( [ dmftcondition ] ) and impurity action ( [ impurity ] ) . if all the quantities entering the green s function @xmath183 ( such as @xmath0 , @xmath43 , @xmath40 ) evolve simultaneously under the gauge transformation , the self - consistency ( [ dmftcondition ] ) clearly stays fulfilled .
the entire green s function @xmath183 obeys the gauge phase rotation , which proves that the theory is conserving .
we can prove that the ladder approximation for dual bosons results in such a conserving theory .
we introduce a variation of the local part of the dispersion law @xmath184 , and corresponding variations of local quantities @xmath185 , and @xmath186 . in the tensorial notation used previously , the variations @xmath187 are vectors with a single index @xmath188 ; their components are defined as fourier - transforms : @xmath189 etc .
the straightforward variation of the self - consistency condition with the dmft green s function ( [ dmftcondition ] ) gives : @xmath190 it is worth to remind that @xmath191 and @xmath192 are tensors with double indices @xmath193 , as well as the quantities @xmath194 and @xmath195 to appear in the following formulas .
now , let us consider the full lattice susceptibility @xmath194 as a fourier - transform of eq.([susc4 ] ) with pairwise - equal space arguments : @xmath196 the self - consistency requires @xmath197 , and therefore the above - introduced variations correspond to @xmath198 similarly , for the susceptibility of the impurity model @xmath199 although the right - hand side formally includes the wavevector @xmath170 , @xmath200 is of course @xmath170-independent .    with the above two formulas formulas
, it immediately follows from ( [ suscdraft ] ) that @xmath201 such a self - consistent condition ( [ susc ] ) is equivalent to locality of the irreducibly vertex and it is used in the dmft scheme for calculations of non - local susceptibilities @xcite .
it is worthwhile to remind that the susceptibilities with a single frequency argument , which we defined previously , can be obtained by summation over the fermionic frequencies , so that @xmath202 and @xmath203 to establish a connection with the dual formalism , we substitute this expression into eq .
( [ xpi ] ) and express @xmath204 using eq .
( [ susc ] ) .
simple algebraic manipulations give @xmath205 one can recognize that @xmath206 and @xmath207 are , respectively , the previously introduced quantities @xmath102 and @xmath99 with proper indices .
furthermore , it follows from the self - consistency condition ( [ dmftcondition ] ) that @xmath208 so we conclude that the obtained expression for @xmath209 is exactly the result of the dual ladder summation ( [ buble ] ) .
finally , let us consider the plasmon dispersion given by the equation @xmath210 in the long - wavelength limit , one considers an expansion in @xmath211 , in particular @xmath212 .
the zeroth - order term in the left - hand side vanishes due to charge conservation .
the second - order term results in the dispersion equation for plasmons in a 3d correlated lattice : @xmath213 for the 2d case , the right - hand side equals @xmath214 .    while deriving the last equation , we supposed that @xmath49 stays finite and can be therefore neglected in comparison with @xmath215 at small @xmath170
. note however that this quantity still enters @xmath174 and therefore affects the left - hand side of the dispersion relation , so that the physics of local plasmonic fluctuations is taken into account .
the single - bubble approximation becomes invalid when vertex @xmath99 is large .
there is a general situation of this kind : the four - point vertex of an impurity problem with _ magnetic moment formed _ appears to have a part displaying a slow dynamics that give a large contribution at low frequencies .
this simply means that for isolated atoms the magnetic moments are free and , thus , for _ almost _ isolated atoms ( when the energy of interatomic effective interactions is much smaller than the intraatomic one ) they can be considered as _ almost _ integrals of motion .    to illustrate this statement , let us consider the atomic limit of the single - band hubbard model .
the single - electron green s function of the hubbard atom at half - filling @xmath216 has a time - scale @xmath217 , which we call fast .
contrary , the spin - spin correlator @xmath218 is independent , as one can easily check , of time arguments .
so , a time - scale for @xmath218 is the inverse temperature @xmath13 .
it means that , apart from a small time domain @xmath219 , the value of @xmath220 is determined by its non - gaussian part , which does not fall as @xmath221 increases .
the much slower dynamics of @xmath222 in comparison with @xmath40 just reflects the fact that a rotation of the spin does not change the energy of the atom , whereas the change of the particle number does .
direct calculation of @xmath99 for the hubbard atom does support this observation .
the finite - temperature expressions @xcite for @xmath223 contain a `` singular '' term proportional to @xmath224 .
so , the first two time arguments of @xmath99 in time - domain must coincide , as well as the last two : @xmath225 , but the difference @xmath226 can be large .
the properties of the impurity problem are different from those of the isolated atom . however , for the strong - coupling limit the hybridization is small . in this case
the fast dynamics roughly stay unchanged , so that the atomic green s function ( [ gatomic ] ) can be used , and the four - point vertex obeys the form @xmath227 .
the hybridization can introduce its own slow timescale @xmath228 , so that @xmath99 does depend only on @xmath226 . in the frequency representation , this means @xmath229 and at small frequencies @xmath230 . since @xmath222 includes a large factor , we can keep only the irreducible term in the expressions for @xmath102 and @xmath39 .
note that this is completely opposite to the case considered in section [ barediagsec ] .
the terms proportional to @xmath231 in eq .
( [ gamma ] ) are small in @xmath232 and therefore can be neglected .
then we obtain the following relation from eq .
( [ gamma ] ) : @xmath233 in the above equation we introduced a `` bare susceptibility '' of the impurity problem @xmath234 next , we use the relation ( [ lambdagamma ] ) and keep only the large factor @xmath99 : @xmath235 note that in the present approximation @xmath102 , like @xmath99 , only has single frequency argument .
let us consider the half - filled hubbard lattice with the nearest - neighbour hopping @xmath236 , at large @xmath33 .
the low - frequency behaviour of this system is essentially the dynamics of the heisenberg model .
the effective - medium description of spin waves in this case should result in the expression @xmath237 with @xmath238 for nearest neighbors .
although this expression is rather simple ( in our notation , it is just @xmath239 ) , it was not reproduced by any dmft - like approximations so far .
the superexchange interaction containing @xmath240 can easily be obtained from a single loop of the fermionic lines , but the rpa and gw+edmft denominators contain an inverse of this quantity ( see section [ barediagsec ] ) , so the structure of the formulas is drastically different .    in our formalism , the desired expression appears from the ladder summation , that is the calculation for the eq .
( [ buble ] ) and corresponding diagram shown in fig .
( [ fig_ladder ] ) . in this summation
, we take into account the `` slow '' vertex ( [ gammaslow ] ) only .
since @xmath99 depends only on a single time argument , the calculation is quite simple . with the exact relation ( [ exact2 ] ) , we obtain @xmath241^{-1 } .\ ] ] finally , we substitute the formulas ( [ irreducc ] ) , ( [ irreduc ] ) for @xmath39 and @xmath102 .
this results in the cancellation of @xmath99 out of the expressions : @xmath242 where @xmath132 from eq .
( [ irreduc ] ) has a meaning of an effective local stoner parameter .
both @xmath243 and @xmath244 obey the `` fast '' dynamics only . to determine the low - energy properties , in which we are interested , it is enough to calculate it at @xmath245 , and to replace the summation over frequency by an integration .
a straightforward calculation of the expression ( [ chi0 ] ) with @xmath246 gives @xmath247 . to calculate @xmath248
, we switch to real - space and use the strong - coupling limit .
it gives the green s function @xmath249 for the nearest neighbors and @xmath250 elsewhere .
now , the integration gives @xmath251 for the nearest - neighbors and we end up with the desired expression ( [ heizenberg ] ) .
this derivation can be considered as a generalization of the effective exchange formula for the magnetically ordered case @xcite to the paramagnetic state .
the standard theory of magnons in itinerant - electron systems @xcite is based on the rpa which is supposed to be exact at zero temperature ; at finite temperatures , magnon - magnon and electron - magnon interactions should be taken into account @xcite .
electron - hole excitations are primary objects in such an approach , spin waves are coherent superposition of such excitations and , in general , are only well - defined in a restricted part of the brillouin zone , due to stoner damping . in the case of strong - coupling hubbard model , we have an `` anti - rpa '' situation when local magnetic moments are primary objects . therefore , completely different approaches should be used , e.g. those based on the hubbard x - operators @xcite .
it is very non trivial that the dual boson approach describes these two opposite cases within the same formalism , thus providing a practical way to interpolate between weak and strong coupling limits .
in summary , the general scheme for non - local correlation effects within the dual boson approach is presented .
this method gives a direct means to obtain the higher - order diagrammatic contributions beyond the edmft scheme free from the double - counting contributions and it can be useful to describe the magnon - spectrum renormalization in strongly correlated systems .
the method combines the numerically exact solution of an effective single - cite impurity model with frequency - dependent electron - electron interactions as well as an analytical diagrammatic expansion of non - local fermionic and bosonic self - energies .
this approach is designed to describe of correlated systems with strong local fermion correlations and fluctuating non - local bosonic modes on equal footing .    in this paper we focused on the applications to the description of bosonic degrees of freedom .
additional contributions to the electron self energy also arise , in comparison with the dual fermion approach @xcite , e.g. , the diagram shown in fig . [ fig_sig ] .
similar diagrams are important in the cases of coexistence of disorder with electron - electron @xcite or electron - phonon @xcite interactions .
it is important to calculate these additional contributions within the hubbard model and make a detailed comparison with dual fermion theory
. this will be the subject of our further studies .
this work was supported by dfg grant no . 436 113/938/0-r and rffi 11 - 02 - 01443 .
mik acknowledges a support from the eu - india fp-7 collaboration under monami .
ail acknowledges a support from sfb 668 and for-1346 .
the formalism presented here has an important difference compare to the common introduction of the edmft : whereas we introduce dual bosons @xmath79 by the hubbard - stratonovich decoupling of the quadratic form @xmath252 , the term decoupling @xmath253 is decoupled in the standard edmft scheme @xcite .
as one can see , this does not result in any physical consequences , as the two approaches yield the identical results for the averages @xmath254 over original variables .
however , properties of the _ dual _ bosons , of course , do depend on their definition . in particular
, the difference is about local properties of the dual propagator : for the present formalism , its local part vanishes in the edmft approximation , @xmath255 whereas the `` effective - medium '' condition @xmath256 holds for the decoupling ( [ vdecoupl ] ) . here
@xmath257 is the green s function of the dual variables @xmath51 arising from the decoupling of @xmath253 : @xmath258 and @xmath259 is defined for the impurity problem with decoupled bosonic hybridization , so that @xmath260 to maintain a similarity with our formalism , we introduced in the above expressions a scaling factor @xmath70 , that is absent in the standard edmft scheme @xcite .
all physical quantities do not depend on a particular value of this factor .
we will now consider a general case of some approximation utilizing a decoupling of @xmath263 with an arbitrary @xmath264 .
the approximation can be non - gaussian , so that there is a correction @xmath265 to the edmft polarization operator : @xmath266 the green s function for the dual bosons can be found from the exact relation @xmath267 one can realize that there are two special cases .
first , let us take @xmath268 , and use @xmath269 , as in the main body of the paper .
for this choice , using ( [ w1 ] , [ exactw ] ) one can straightforwardly show that @xmath270 summing over @xmath8 and taking eq.([dbcondition ] ) into account we conclude that the condition ( [ dual0 ] ) is fulfilled @xmath271 it turns out that _ the dual propagator vanishes for variables obtained from the hubbard - stratonovich decoupling of the renormalized interaction , with the hybridization part * excluded * _ , if the condition ( [ dbcondition ] ) is fulfilled .
for the edmft approximation , @xmath265 is absent and the condition ( [ dual0 ] ) holds for the decoupling with @xmath272 used in the present paper .
now , let us consider the case @xmath273 .
it is suitable to take @xmath274 in this case .
we obtain @xmath275 again taking a summation over @xmath8 and substituting eq.([dbcondition ] ) , we obtain @xmath276 .
the second line of ( [ exactv ] ) with @xmath274 gives the same answer for @xmath261 , therefore we obtain @xmath277 _ the `` effective - medium '' condition holds for variables obtained from the hubbard - stratonovich decoupling of the renormalized interaction , with the hybridization part * included * _ , if condition ( [ dbcondition ] ) is fulfilled . for edmft
, this corresponds to the decoupling of purely @xmath253 part used in ref .
@xcite .
we stress that , strictly speaking , apart from the simple edmft case the conditions ( [ dual0 ] ) , ( [ dualem ] ) are not fulfilled simultaneously with ( [ dual0 ] ) .
this does not cause serious problems for the gw+edmft calculations , because one typically does self - consistent calculations of the edmft problem and then calculates diagram corrections for a pure edmft hybridization . on the other hand , for higher order diagrammatic approximations
the difference can be essential .
finally , we note that a specific choice of the decoupling procedure is physically irrelevant .
as we have seen , edmft can be constructed using @xmath278 or @xmath279 , and the result for initial ensemble is the same .
it can be checked that any other decoupling will also work ; what matters is the gaussian approximation for dual variables and the value of the hybridization .
the latter is defined by the condition ( [ dual0 ] ) , that does not include any dual quantities .
we use the theory based on @xmath280 decoupling because it allows a simple construction of the dual diagrammatic series without a double counting of the local contributions .
this important property is based on the fact that the edmft a dual green s function has a very simple physical meaning : according to ( [ ggprime ] ) , it is just the green s function of initial variables with local part subtracted .
e.  gull , a.  j. millis , a.  i. lichtenstein , a.  n. rubtsov , m.  troyer , p.  werner , rev .
83 ( 2011 ) 349 .
e.  gorelov , j.  koloren   , t.  wehling , h.  hafermann , a.  b. shick , a.  n. rubtsov , a.  landa , a.  k. mcmahan , v.  i. anisimov , m.  i. katsnelson , a.  i. lichtenstein , phys .
b 82 ( 2010 ) 085117 .
e.  gorelov , t.  o. wehling , a.  n. rubtsov , m.  i. katsnelson , a.  i. lichtenstein , phys . rev .
b 80 ( 2009 ) 155132
. g.  kotliar , s.  y. savrasov , k.  haule , v.  s. oudovenko , o.  parcollet , c.  a. marianetti , rev .
78 ( 2006 ) 865 .
a.  i. lichtenstein , m.  i. katsnelson , phys .
b 57 ( 1998 ) 6884
. m.  i. katsnelson , v.  y. irkhin , l.  chioncel , a.  i. lichtenstein , r.  a. de  groot , rev .
80 ( 2008 ) 315 .
r.  chitra , g.  kotliar , phys .
b 63 ( 2001 ) 115110 .