the ab initio band structure calculations of strongly - interacting low - dimensional ( quasi - one - dimensional ( q1d ) or two - dimensional ( 2d ) ) systems ( high-@xmath1 superconductor yba@xmath2cu@xmath3o@xmath4 as well as charge - density - wave ( cdw ) systems k@xmath5moo@xmath3 and bavs@xmath3 being examples @xcite ) usually reveal several bands in the vicinity of the fermi level that are , in the leading approximation , decoupled from the rest of the band structure . in the tight - binding approximation ,
the bare dispersions of such a multiband model are well understood in terms of the strong hybridization among several properly chosen orbitals per unit cell , while the hybridization with the other electronic levels is neglected .
@xcite the bare single - particle tight - binding hamiltonian is thus taken to include commensurate / incommensurate cdw or sdw ( spin - density wave ) average fields and several first - neighbour bond energies .
in the two - particle hamiltonian one usually retains the long - range coulomb interactions , the largest local ( hubbard ) interactions together with the most important short - range interactions . associated with different terms in the single - particle hamiltonian
are fluctuations in the charge and spin densities on different orbitals in the unit cell and on the corresponding bonds .
all these ingredients make strongly - interacting multiband electronic problems extremely complicated .
to explain most of experimental results on typical single - band systems , it usually suffices to describe properly the coupling between @xmath6 monopole charge fluctuations and external electromagnetic ( em ) fields and combine this coupling hamiltonian with the ordinary transport equations .
@xcite however , in multiband electronic systems , even in the cases when only one band intersects the fermi level , the analysis is considerably more complicated . not only the interband physics but also the intraband physics depends drastically on both single - particle and two - particle parameters of the hamiltonian and are affected differently by various types of fluctuations .
therefore , to understand the results of transport measurements @xcite and reflectivity measurements @xcite of multiband electronic systems in detail , we describe the @xmath6 multipolar fluctuations in the dynamical conductivity tensor in a way consistent with the results of experimental methods which probe charge and spin fluctuations which presumably occur at finite wave vectors @xcite .
the @xmath6 charge - charge correlation function is thus just one of several equally important correlation functions defined in terms of symmetrized intraband and interband charge and spin fluctuations .
common to all these correlation functions is a complicated structure of the intraband and interband electron - hole propagators . therefore , to understand the contribution of different types of fluctuations to these correlation functions , in particular to the dynamical conductivity
, we must first find the solution to quite general bethe
salpeter equations for the electron - hole propagators and then combine the result with the very definition of the response function under consideration .
this question , which is of great importance in studies of strongly - interacting low - dimensional systems , is in the focus of the present analysis .    in this article
, we thus present a general response theory which uses both an essentially exact semiclassical description of the coupling between valence electrons and external em fields @xcite and a quite general description of the intraband and interband charge fluctuations @xcite .
the result is the current - dipole kubo formula for the dynamical conductivity , @xcite with the structure of the intraband electron - hole propagators determined by means of the quantum transport equation ( i.e. , the intraband bethe
salpeter equation ) . to emphasize the connection between the present theory and the ordinary transport theory , @xcite we determine the explicit form of the intraband conductivity in the memory - function approximation .
@xcite we also find the structure of the interband conductivity in both cases , with real and with purely imaginary interband dipole vertices , to show in which way the lorentz local field corrections @xcite affect the dynamical conductivity of interacting low - dimensional systems .
the article is organized as follows . to make the reading of the article easier ,
we give in section  2 an overview of the results .
then in section  3 the notation is introduced to describe electrons and boson modes in a general multiband electronic model with boson - mediated electron - electron interactions . a complementary multiband electronic model with nonretarded electron - electron interactions will be studied in a separate article @xcite which is focused on the commensurability effects in the single - particle conductivity tensor .    starting with the quantum electrodynamical description of the interaction between valence electrons and external em fields ,
we derive the exact semiclassical form of the electron - em field coupling hamiltonian in section  4 . in section  5
we define the dynamical conductivity tensor and briefly discuss the relaxation - time approximation .
section  6 gives the structure of the self - consistent bethe
salpeter equations for the dipole vertex functions and for the related auxiliary electron - hole propagators ( to be defined latter ) in a general multiband case .
we show the equivalence of the bethe
salpeter equation for the auxiliary intraband electron - hole propagators and the quantum transport equation .
the generalized drude formula is derived from the quantum transport equation by using the memory - function approximation .
the bare interband conductivity is considered in section  7 and the related local field effects in section  8 .
section 9 contains the concluding remarks .
before discussing the details of the dynamical conductivity model proposed in this article , we present an overview of the results . in this way principal physical messages of the present work are separated from the technical details .
we are interested here in a general low - dimensional model with multiple electronic bands in the presence of boson - mediated electron - electron interactions .
the single - particle part of the total hamiltonian is assumed to be exactly solvable .
the bloch energies @xmath7 in the band labeled by the band index @xmath8 and all relevant bare vertex functions ( the dipole vertices @xmath9 , for example ) are thus taken as known functions . to study the influence of the boson - mediated electron - electron interactions on single - electron propagators @xmath10 and , in addition , on different electron - hole propagators encountered in the response functions of interest , we use high - order diagrammatic perturbation theory of the metallic state , as usual in investigations of low - dimensional electronic systems .
@xcite    the electron - hole propagators are associated with the dispersive , undamped bosons which can either correspond to the external degrees of freedom coupled to electrons / holes or the internal degrees of freedom built of the electron - hole excitations at finite values of * q*. in sections  68 , the emphasis is on the acoustic and raman - active and infrared - active optical phonons . the scattering by phase and amplitude phonon modes of the common peierls cdw model will be studied in detail within the microscopic lee  rice  anderson model in the accompanying article @xcite ( to be referred to as article ii ) .    ) for electron green s functions .
@xmath8 is the electron band index and @xmath11 is the boson branch index . the black circle and
the big gray circle represent , respectively , the bare and the renormalized electron - boson vertex . here and
herafter , the tiny and bold solid ( dashed ) lines are the bare and renormalized electron ( boson ) propagators . ]    in high - order perturbation theory , the electron green s function @xmath10 satisfies the dyson equation @xcite
@xmath12   { \cal g}_l ( { \bf k } , { \rm i } \omega_n ) = \hbar , \label{eq1}\end{aligned}\ ] ] where @xmath13 is the self - energy function of the electron ( see figure  1 ) .
in this case @xmath13 is a function of @xmath10 , resulting in a system of two coupled equations .
a strong reduction of the quasi - particle pole weight in @xmath10 at the fermi level is predicted by different strictly 1d theories @xcite and observed in angle - resolved photoemission spectroscopy ( arpes ) experiments on q1d systems @xcite . in this article
we present a systematic method of calculating the dynamical conductivity tensor in interacting low - dimensional systems which is consistent with the described treatment of the single - electron green s functions @xmath10 .      in semiclassical electrodynamics applied to a general exactly solvable single - particle problem with multiple electronic bands ,
the electric component of the em field parallel to the cartesian axis @xmath14 , @xmath15 , couples to the currents of free charges and to the dipole moments of bound charges in a way described by the coupling hamiltonian @xmath16 here , @xmath17 is the total dipole density operator .
it is the sum of the intraband and interband contributions @xmath18 with @xmath19 representing the corresponding intraband ( @xmath20 ) and interband ( @xmath21 ) dipole vertices , and @xmath22 .
it is not hard to verify that these vertices are directly related to the vertices in the current and charge density operators , @xmath23 and @xmath24 .
the relation is given by the expression @xmath25 ( @xmath26 in the longitudinal case ) .
the imaginary currents and the imaginary dipole moments do not have a straightforward physical meaning , but they make the notation used throughout this article as simple as possible . for example , they can be used to express the gauge invariance of the expression ( [ eq2 ] ) simply by rewritting @xmath27 in terms of @xmath24 and @xmath23 and using the standard relations between @xmath28 and the scalar potential @xmath29 and the vector potential @xmath30 .      in order to examine how electrons respond to applied em fields
, we introduce the dynamical conductivity tensor @xmath31 , @xcite which is the response function relating the induced total current @xmath32 to @xmath33 ; @xmath34 .
we use here the finite temperature formalism .
the dynamical conductivity @xmath31 is in this case given by analytical continuation ( @xmath35 ) of @xmath36 , which is the matsubara fourier transform of @xmath37 .
the compact form of the semiclassical coupling hamiltonian ( [ eq2 ] ) allows us to express @xmath37 in the form of the current - dipole kubo formula @xcite @xmath38 \rangle_{\rm irred}. \label{eq5}\end{aligned}\ ] ] here , @xmath39 is the total current density operator defined in a way analogous to ( [ eq3 ] ) ( see ( [ eq23 ] ) ) . in order to extract the most important physical consequences that follow from this expression , one usually uses the relations ( [ eq4 ] ) to rewrite @xmath37 either in terms of the charge - charge or the current - current kubo formulae .
@xcite in contrast to these approaches , we use here ( [ eq5 ] ) , because this expression will be shown below to provide the direct link of the present analysis to the common transport theory @xcite and to the textbook approaches to the dynamical long - range coulomb screening @xcite .     in a general multiband case .
the dynamical conductivity @xmath36 , equation ( [ eq5 ] ) , is the rpa - irreducible part of @xmath40 .
@xcite ]    to complete the response theory , we must take into account the relaxation processes describing the way in which electrons are scattered by static disorder , by other electons or by different types of boson modes .
since our primarily interest is in the low - dimensional systems where the quasi - particle pole weight in the single - electron green s function is strongly reduced by boson - mediated electron - electron interactions , we must calculate @xmath37 beyond the common weak coupling theory ( corresponding to second - order perturbation theory for the single - electron self - energy function ) .
this means that in order to work out @xmath37 we must use the general bethe
salpeter expression for the screened conductivity tensor illustrated in figure  2 and the definition @xmath41 , where @xmath42 is the external electric field @xcite .
the bold solid lines in these two diagrams are the single - electron green s functions which satisfy the dyson equation ( [ eq1 ] ) .
the intraband and interband contributions to @xmath43 can be represented by the same expression , @xmath44 equally well for @xmath20 as for @xmath21 . here
the @xmath45 are the renormalized dipole vertex functions which satisfy the self - consistent bethe
 salpeter equations for the renormalized dipole vertex , and @xmath46 . in order to emphasize the connection between the present approach and the common transport theory , as well as to simplify the notation
, we introduce the auxiliary electron - hole propagator @xmath47 by comparing the relation @xmath48 with ( [ eq6 ] ) .
therefore , our task is to solve the self - consistent equations for @xmath47 rather than that for the renormalized dipole vertices .
we limit our discussion to the low - dimensional systems where the scattering by static disorder and the scattering processes in which electrons change the band can both be neglected .
these processes have negligible influence on the properties studied here .
however , the self - consistent equations still have a very complex structure in particular as a consequence of the following features : the nesting properties of the fermi surface , the resonant nature of the scattering by the ( quasi)static cdw or sdw potentials and the local - field effects in the interband conductivity channel .
the results of the present analysis are the following : _
( i ) _ we show the formal solution to the bethe
salpeter equations for @xmath49 and to the ( intraband ) quantum transport equation , _
( ii ) _ we compare the intraband conductivity to that obtained by using the memory - function approximation , _ ( iii ) _ we calculate the bare interband conductivity and _ ( iv ) _ we rederive the well - known lorentz  lorenz form of the optical conductivity in the presence of scattering by boson modes .    _
( i ) _ although it is tempting to completely neglect the electron - em field vertex corrections in @xmath36 , this is not justified in the intraband channel , at least in the cases with low - dimensional conductivity studied here .
the method by which the intraband conductivity tensor is calculated here is very general and is based on a systematic analysis of the self - consistent equation for the auxiliary electron - hole propagator @xmath50 in which the single - electron self - energy contributions and the related vertex corrections are treated on an equal footing .
the equation is of the form @xmath51 \phi_\alpha ( { \bf k},{\bf k}_+ , { \rm i } \omega_n , { \rm i } \omega_{n+ } )   \nonumber \\ & &   \hspace{5 mm } = \frac{1}{\hbar^2 } \big [ { \cal g } ( { \bf k } , { \rm i } \omega_n )   - { \cal g } ( { \bf k}_+ , { \rm i } \omega_{n+ } ) \big ] p_\alpha ( -{\bf q } ) , \label{eq8}\end{aligned}\ ] ] with @xmath52 and @xmath53 . for simplicity
we assume that only one band intersects the fermi level and omit explicit reference to the conduction band index .
@xmath54 is the electron - hole self - energy , the structure of which is derived here for the case in which conduction electrons are scattered only by various boson modes .
these boson modes are assumed to dissipate momentum by their own means ( impurities , umklapps , etc . ) .
pinned phasons of the commensurate cdw problem , for example , enter in the present theory in the same way as the infrared - active optical phonons from section 8.2 .
this question is discussed in more details in article ii .
the electron - hole self - energy is found to be the difference between two auxiliary single - electron self - energies @xmath55 , and the contributions of the @xmath6 forward scattering processes are found to drop out of both @xmath55 and @xmath54 , in full agreement with the charge continuity equation .
equation ( [ eq8 ] ) is valid under quite general conditions , and is therefore applicable to a rich variety of problems , including purely electronic models with nonretarded local and short - range interactions .
for example , it is not hard to verify that the normal scattering processes do not lead to the resistivity in this case , as required by the general principle of microscopic reversibility .
@xcite this question will be discussed in more details in ref .  .
_ ( ii ) _ in the ordinary 3d metallic regime , equation  ( [ eq8 ] ) leads to the ordinary transport equation , with @xcite @xmath56 being the resulting intraband electrical conductivity , and @xmath26 again . here
, @xmath57 is the unity of mass , commonly taken as the free electron mass , and @xmath58 is the intraband relaxation rate . in the drude @xmath59 limit , where @xmath60 , the @xmath61 dependent term in the denominator of ( [ eq9 ] ) can be ignored , and we obtain the intraband conductivity ( [ eq10 ] ) . for @xmath62 and @xmath63 , on the other hand , the result is the well - known thomas  fermi expression @xcite for @xmath64 and for the dielectric susceptibility @xmath65 .    _
( iii ) _ the general expression for the interband single - particle conductivity is derived and compared to the approximate expression which is frequently encountered in the literature .
@xcite        _ ( iv ) _ if the scattering of valence electrons is sufficiently weak , the electron - boson vertex renormalizations are of no importance and the electron - hole propagators of figure  2 satisfy the bethe
salpeter equations shown in figure  3 with two types of contributions , the ladder contributions and the rpa ( random phase approximation ) contributions .
@xcite these equations can be completely rewritten in terms of the renormalized single - electron green s functions and the dynamically screened electron - electron interactions .
this approximation is used here to illustrate how the conductivity tensor is affected by the short - range dipole - dipole interactions .
the result is the well - known lorentz - lorenz form of the dynamical conductivity for 3d metallic @xmath0 models @xcite which is believed to be valid in the q1d conductors , as well @xcite . from the result @xmath66 we can conclude that the contributions to the conductivity tensor of the free charges and of the bound charges are decoupled and that the lorentz dipole - dipole interactions affect only the bound ( interband ) charge contribution to the conductivity tensor .
let us first introduce the quantities characterizing a general multiband electronic model with boson - mediated electron - electron interactions .
the total hamiltonian is the sum of two contributions , @xmath67 . here , @xmath68
\label{eq12 }   \end{aligned}\ ] ] is the bare electron - boson hamiltonian , and @xmath69 is the scattering hamiltonian which describes the scattering of valence electrons by each of the boson modes . as already mentioned , the bosonic modes can be external to electronic system ( e.g. phonons ) or built from electron - hole fluctuations at finite wave vectors . in ( [ eq12 ] ) , @xmath7 is the bare electron dispersion , @xmath70 is the bare boson frequency and @xmath71 is the corresponding `` mass '' parameter .
finally , @xmath72 is the generalized scattering potential . for convenience
, it is assumed to be a product of the coupling constant @xmath73 , the dimensionless electron - boson vertex @xmath74 and the boson field @xmath75 , @xmath76 in the intraband scattering approximation , in which the electron does not change the band when it is scattered by bosons , the scattering potential @xmath72 is diagonal in the band index .      ) for boson green s functions . ]    in the finite temperature formalism , the effects of @xmath77 on the single - electron and boson properties are described in terms of two exact matsubara green s functions @xcite @xmath78 \rangle ,   \nonumber \\ & & \hspace{-11 mm } { \cal d}_\nu ( { \bf q } ' , \tau ) = - \langle t_\tau \big[q_{\nu { \bf q}'}(\tau )   q^\dagger_{\nu { \bf q}'}(0 ) \big ] \rangle .
\label{eq15 } \end{aligned}\ ] ] their fourier transforms satisfy two dyson equations , equation  ( [ eq1 ] ) and @xmath79   - \hbar \sigma_\nu ( { \bf q } ' , { \rm i } \nu_n ) \big\ } { \cal d}_\nu ( { \bf q } ' , { \rm i } \nu_{n } ) = \hbar \label{eq16}\ ] ] ( see figure  4 ) .
it is convenient to adopt an abbreviated notation by writting @xmath80 this expression can also be shown in the form @xmath81 , where @xmath82 \rangle , \label{eq18}\end{aligned}\ ] ] is by definition the force - force correlation function @xcite ( ffcf ) associated with the potential @xmath72 .
( the big gray circle ) and the completely irreducible four - point vertex function @xmath83 ( the white rectangle ) . ]    in high - order perturbation theory the self - energy function of the electron , @xmath84 , is usually represented by the sum of skeleton diagrams shown in figure  5 , leading to a simple analytical expression @xmath85 here , @xmath86 is the renormalized ffcf obtained by replacing one bare electron - boson vertex in @xmath87 , @xmath88 , by the renormalized vertex @xmath89 .
equation ( [ eq19 ] ) is nothing but the three - point vertex representation of @xmath84 . in the same representation
the boson self - energy takes the form @xmath90 we have just reemphasized that the three - point vertex representation is very useful in describing the closed system of four equations characterizing strongly - interacting electron - boson systems in the absence of external em fields ( equations  ( [ eq1 ] ) , ( [ eq16 ] ) , ( [ eq19 ] ) and ( [ eq20 ] ) ) .
however , the conductivity tensor of such an electron - boson model has a compact form if the structure of the bethe  salpeter equations of section 6 is also simple , at least from the formal point of view . for this purpose
, it is much more convenient to use the four - point vertex representation .
this representation is based upon the knowledge of the irreducible four - point interaction @xmath91 .
the equivalence of the two representations in the skeleton approximation is illustrated in figure  5 , where @xmath91 is represented by the white rectangle .
the expression for @xmath84 is now of the form @xmath92 a full discussion of the bethe
salpeter equations is postponed to section 6 .
in quantum electrodynamics of tightly bound valence electrons in solids , whether metallic or insulating , the coupling between electrons and external em fields is given by the minimal gauge invariant substitution .
the expansion to the second order in the vector potential @xmath93 gives the coupling hamiltonian @xmath94 , with @xcite @xmath95 here , @xmath96 are , respectively , the total current density operator and the bare diamagnetic density operator .    to obtain the semiclassical expression for the coupling hamiltonian , equation  ( [ eq2 ] )
, we have to take into account the fact that both the interband and intraband current - current contributions to the diamagnetic current must vanish in the normal metallic or insulating state .
the cancellation of the interband contribution is a consequence of the effective mass theorem , @xcite @xmath97 leading to @xmath98 , with @xcite @xmath99 ,   \nonumber \\ & & \hspace{-5 mm } h^{\rm ext}_2 =   \frac{e^2}{2mc^2 } \sum_{l { \bf q } { \bf q}'\alpha \beta }   a_{\alpha } ( { \bf q}-{\bf q}')a_{\beta } ( { \bf q } ' )   \hat{\gamma}^{ll}_{\alpha \beta } ( -{\bf q } ) .
\label{eq25 }   \end{aligned}\ ] ] here , @xmath100 , @xmath21 , is the interband part in dipole density operator ( [ eq3 ] ) , and @xmath101 is the diamagnetic density operator . the reciprocal effective mass tensor @xmath102 is the vertex function in @xmath101 .
the intraband current - current contribution to the diamagnetic current cancels @xmath103 , confirming the assertion made in section 2 that @xmath104 we emphasize the simple structure of this coupling hamiltonian .
the most direct way to obtain this expression for @xmath27 is to consider the coupling of the scalar potential @xmath105 to the charge density operator .
@xcite      let us derive now the relations ( [ eq4 ] ) among three types of intraband and interband vertex functions underlying the gauge invariance of the present response theory . to do this
, we use very general microscopic operator equations @xmath106 , \nonumber \\ & & \hspace{-10 mm }
{ \rm i } \hbar \hat j_\alpha({\bf q } ) = { \rm i } \hbar \frac{\partial}{\partial t }   \hat p_\alpha({\bf q } ) = \big[\hat p_\alpha({\bf q } ) , h_0^{\rm el } \big ] , \label{eq26}\end{aligned}\ ] ] where @xmath107 is the electronic part in the bare hamiltonian ( [ eq12 ] ) .
furthermore , @xmath108 is the total charge density operator , with @xmath109 being the charge vertex function , while @xmath39 and @xmath110 are given by ( [ eq23 ] ) and ( [ eq3 ] ) .
the two commutators in ( [ eq26 ] ) are easily worked out , leading to the required relations @xmath111 the explicit form of each of these vertices is often obtained by calculating the coefficients in the linear term in @xmath27 in ( [ eq22 ] ) .
the long - wavelength intraband vertices are of particular interest .
they are given by model independent expressions @xmath112 here , @xmath113 is the bare electron group velocity and @xmath52 is the bare dipole vertex .
notice that the thermodynamic averages of the density operators in ( [ eq26 ] ) are the usual induced macroscopic densities .
for example , @xmath114 , @xmath115 and @xmath116 represent the common notation for the induced macroscopic current , charge and dipole densities .
they are related by two equations @xcite @xmath117 the first one is the usual charge continuity equation .
we now define the total dynamical conductivity tensor @xmath31 as the response function to the macroscopic electric field @xcite @xmath118 that is @xmath119 here , @xmath114 is the induced total current density from ( [ eq29 ] ) , @xmath42 is the external field and @xmath120 is the depolarization field ( equal to @xmath121 for the longitudinal polarization of the field in the thin slab ) .
the conductivity tensor is given by analytical continuation to the real axis of @xmath36 ( see figure  6 ) , where @xmath36 is the matsubara fourier transform of @xmath38 \rangle_{\rm irred}. \nonumber   \end{aligned}\ ] ] therefore , to find the dynamical conductivity tensor of an interacting electron - boson system , we must determine the structure of the current - dipole kubo formula ( [ eq5 ] ) which agrees with the definition of the macroscopic electric field ( [ eq30 ] ) .     from figure  2 , shown in the three - point vertex representation ( left diagram ) and in the four - point vertex representation ( right diagrams )
. the collection of diagrams in the renormalized four - point interaction must be consistent with the definition of @xmath122 in ( [ eq30 ] ) ; i.e. , the renormalized four - point interaction is irreducible in the long - range rpa contribution . ]
it is instructive first to show the structure of the bare dynamical conductivity ( [ eq5 ] ) in the simplest multiband model , namely , a two - band model with the lower band partially occupied by electrons and the upper band empty , in the limiting case where the scattering is sufficiently weak . as explained in more details in article ii , this model is of importance in accounting for the properties of the ordered state of some simple cdw systems .
@xcite    to do this , we first calculate the ideal bare conductivity tensor @xmath123 which is given by the first term in figure  6 , with @xmath124 replaced by the bare green s function @xmath125 .
summation over @xmath126 can be performed easily to give @xmath127 .
after analytical continuation , we obtain the ideal bare intraband ( @xmath20 ) and interband ( @xmath128 ) contributions described by @xmath129 } { \hbar \omega + { \rm i } \hbar \eta   + \varepsilon_{l}({\bf k})- \varepsilon_{l'}({\bf k}_+ ) } .
\label{eq32}\end{aligned}\ ] ] weak scattering processes can be inserted into this expression by replacing the adiabatic term @xmath130 by the relaxation rate @xmath131 independent of frequency and wave vector .
this is the generalization of the common relaxation - time approximation .
as required , the associated transport relaxation time @xmath132 differs from that found in the single - electron propagators .
the bare intraband contribution is of the drude form @xmath133 with @xmath134 ^ 2 \bigg(-\frac{\partial f_l({\bf k})}{\partial \varepsilon_l ( { \bf k})}\bigg ) \label{eq34}\end{aligned}\ ] ] representing the effective number of conduction electrons , @xcite and @xmath135 by assumption .
@xmath136 includes the conduction holes from the lower band and the thermally activated electrons from the upper band .
the bare interband contribution is the sum of two off - diagonal terms in ( [ eq32 ] ) , with @xmath130 replaced by @xmath131 : @xmath137 .
the main disadvantage of this simple procedure of calculating @xmath138 is that it does not apply to low - dimensional systems with sizeable interactions .
the reason for that is the following .
first , comparison with experimental data shows that the relaxation rates @xmath131 are frequency dependent .
second , the causality principle requires that @xmath139 is just the imaginary part of a complex quantity which will be called here the electron - hole self - energy function .
finally , it turns out that in the strong coupling limit the electron - hole self - energy is not a simple function of frequency @xmath140 .
instead , it depends quite drastically on two internal electron frequencies too .
we shall therefore present now a systematic treatment of scattering processes in interacting electron - boson systems with multiple electronic bands which treats all elements in the electron - hole self - energy consistently with the charge continuity equation ( [ eq29 ] ) , with the causality principle and with the definition of @xmath122 .
the resulting self - consistent equations turn out to be the generalization of the equations used by vollhardt and wlfle @xcite to study the scattering by static disorder in a single - band case .
the term bethe  salpeter equations @xcite refers to the self - consistent equations for the exact electron - hole propagators @xmath141 shown in figure  7 .
a widely - used simplification in these equations is illustrated in figure  3 . in this case , we have both the ladder contributions and the rpa - like contributions to the irreducible four - point interaction , but the three - point vertex renormalizations are neglected both in these equations and in the related dyson equations .
the interaction lines and the electron propagators are renormalized , so that these equations must be treated consistently with the dyson equations , as already pointed out in section 2 . however , in strictly 1d metallic systems as well as in strongly - correlated low - dimensional systems , it is necessary to use the bethe
salpeter equations beyond this approximation by taking into account the three - point vertex renormalizations as well .    to simplify the analysis , it is convenient to multiply the left - hand side of the diagrams in figure  7 by the corresponding bare vertex function and then sum over @xmath142 and @xmath143 . for the electron - hole propagators in the dynamical conductivity ( [ eq5 ] ) ,
the bare vertex in question is the bare dipole vertex function , and the bethe
salpeter equations from figure  7 turn into the bethe  salpeter equations for the auxiliary electron - hole propagators
@xmath144 the equations are illustrated in figure  8 .    with little loss of generality , we restrict the analysis to the case where the irreducible four - point interaction @xmath145 can be approximated by @xmath146 . in this case , the explicit form of the intraband ( @xmath20 ) and interband ( @xmath21 ) equations is the following @xmath147   \nonumber \\ & &   \hspace{0 mm } \times   \bigg\ { -p_{\alpha}^{l'l}({\bf k}_+,{\bf k } ) - \frac{1}{\beta }   \sum_{l_1l_1 ' } \sum_{{\bf q } ' { \rm i } \nu_m } u^{l'l_1',l_1l } ( { \bf q } ' , { \rm i } \nu_{m } ) \nonumber \\ & &   \hspace{10 mm } \times \phi^{l_1l_1'}_\alpha ( { \bf k}+{\bf q}',{\bf k}_++{\bf q } ' , { \rm i } \omega_m , { \rm i } \omega_{m+ } ) \bigg\ } , \label{eq36}\end{aligned}\ ] ] with @xmath148 and @xmath149 .
the @xmath150 satisfy the dyson equations for electrons ( [ eq1 ] ) and @xmath146 is a simple generalization of the completely irreducible four - point interaction @xmath151 from ( [ eq21 ] ) .    ) for auxiliary electron - hole propagators . ]    in the electronic multiband model we are interested in here , there are @xmath152 bands in the vicinity of the fermi level .
the bands are built of @xmath152 orbitals per unit cell .
associated with each of these orbitals are the intracell charge and spin fluctuations , and these intracell fluctuations are possible sources of scattering .
the monopole charge fluctuations underlying the present analysis of the dynamical conductivity tensor might be also strongly affected by the presence of the intracell fluctuations . as a result
, @xmath31 will have considerably more structure than that characterizing the ordinary single - band model , even in the case where only one band out of @xmath152 valence bands intersects the fermi level .      in the remainder of this section we consider the intraband term in ( [ eq36 ] ) from this point of view .
we assume that there is only one conduction band and again drop reference to the band index .
we use the intraband scattering approximation where the scattering processes in which the electron changes the band are neglected .
as discussed in article ii , this approximation is easily seen not to be adequate , for example , to study fluctuations in the order parameter in the pseudogap state of the ordinary cdw systems . however , it is convenient to illustrate the main ideas regarding the quantum transport equations . finally , the three - point vertex renormalizations in the renormalized interactions are also neglected , as usual in the common non - crossing approximation for the electron - boson systems . @xcite the resulting intraband bethe
salpeter equation is of the form @xmath153
\bigg\ { p_{\alpha}(-{\bf q } )   \nonumber \\ & &   \hspace{-5 mm } - \frac{1}{\beta } \sum_{{\bf q } ' { \rm i } \nu_m } { \cal f } ( { \bf q } ' , { \rm i } \nu_{m } ) \phi_\alpha ( { \bf k}+{\bf q}',{\bf k}_++{\bf q } ' , { \rm i } \omega_m , { \rm i } \omega_{m+ } ) \bigg\}. \nonumber \\ \label{eq38}\end{aligned}\ ] ] the structure of ( [ eq38 ] ) is simplified by defining the auxiliary single - electron self - energy function @xmath154 and the electron - hole self - energy @xmath155 equation ( [ eq38 ] ) now reduces to the quantum transport equation @xmath156 \phi_\alpha ( { \bf k},{\bf k}_+ , { \rm i } \omega_n , { \rm i } \omega_{n+ } )   \nonumber \\ & &   \hspace{5 mm } - \frac{1}{\hbar^2 } \big [ { \cal g } ( { \bf k } , { \rm i } \omega_n )   - { \cal g } ( { \bf k}_+ , { \rm i } \omega_{n+ } ) \big ] p_\alpha ( -{\bf q } )   \nonumber \\ & &   \hspace{-5 mm } =   - { \rm i } \pi^i ( { \bf k},{\bf k}_+ , { \rm i } \omega_n , { \rm i } \omega_{n+ } ) \phi_\alpha ( { \bf k},{\bf k}_+ , { \rm i } \omega_n , { \rm i } \omega_{n+ } ) .
\label{eq41}\end{aligned}\ ] ] equation  ( [ eq41 ] ) is the self - consistent equation for the auxiliary electron - hole propagator , in contrast to the ordinary transport equation , equation  ( [ eq47 ] ) , which is formulated in terms of the induced density @xmath157 defined by @xmath158 the conditions under which ( [ eq41 ] ) reduces to the ordinary transport equation are discussed below .    to illustrate the role of the electron - em field vertex renomalizations in ( [ eq41 ] ) ,
it is useful to define the auxiliary single - electron propagator @xmath159 as the single - particle propagator which satisfies the following dyson equation @xmath160 \widetilde { \cal g } ( { \bf k } , { \rm i } \omega_n ) = \hbar . \label{eq43}\end{aligned}\ ] ] the forward scattering processes cancel identically out in @xmath161 as a result of the charge continuity equation .
the fact that the auxiliary propagator @xmath162 has the quasi - particle pole even for relatively strong scattering processes is consistent with the fact that the residue associated with the quasi - particle pole in the exact propagator @xmath163 is proportional to @xmath164 and vanishes in the strictly 1d limit @xmath165 .
@xcite this interesting observation is in the background of a variety of approximate treatments of various pseudogapped states . in these approaches ,
the forward scattering contributions to the ffcf @xmath166 are neglected from the outset .
the expression for the intraband conductivity which correctly describes the @xmath59 drude limit , but which is not necessarily correct in the @xmath167 thomas  fermi limit , is the following @xmath168 this is nothing but ( [ eq32 ] ) for the ideal intraband conductivity in which @xmath169 is replaced by @xmath54 .      in the weak coupling limit ,
the self - energy @xmath54 is a non - singular function of two fermion frequencies and in the simplest approximation can be written in the form @xmath170 .
@xmath171 is called the memory function .
@xcite the momentum distribution function @xmath172 appearing in ( [ eq44 ] ) in the memory - function approximation reduces in the weak coupling limit to the fermi
dirac distribution function @xmath173 .
we now define the electron - boson coupling function @xmath174 , the mass enhancement factor @xmath175 , the intraband relaxation rate @xmath176 and the collision integral @xmath177 by @xmath178 summation of ( [ eq41 ] ) over @xmath126 , together with analytical continuation @xmath35 and multiplication by @xmath179 , gives the ordinary transport equation @xcite @xmath180 \delta n ( { \bf k } , { \bf q } , \omega ) + \frac{m}{m_{\alpha \alpha}^*({\bf k } , \omega ) } e \hbar v_\alpha^0 ( { \bf k } )   \nonumber \\ & &   \hspace{20 mm } \times \frac{\partial f({\bf k})}{\partial \varepsilon ( { \bf k } ) } e_\alpha ( { \bf q } , \omega ) =   \hbar i({\bf k } , \omega ) , \label{eq47}\end{aligned}\ ] ] where @xmath181 .    in solving ( [ eq47 ] ) it
is convenient to follow the usual textbook procedure .
@xcite this equation can be solved either for @xmath182 or for @xmath183 , where @xmath184 are , respectively , the induced charge and current densities of conduction electrons , and the @xmath185 , @xmath186 , are the contributions to @xmath187 which are proportional to @xmath188^i$ ] . for further considerations of the random - phase approximation in section 8.1 ,
the second route is more appropriate choice . in this case
, we consider the thin slab with the longitudinal polarization of the electric fields , introduce the scalar fields @xmath29 and @xmath189 , by the relations @xmath190 and @xmath191 , and use the definition of the dielectric susceptibility @xmath192 : @xmath193 , or @xmath194 \rho^{\rm c } ( { \bf q},\omega )   = v^{\rm ext } ( { \bf q},\omega ) .
\label{eq49}\ ] ] it is readily seen that the intraband conductivity is of the form @xmath195 ^ 2 \big ( -\partial f({\bf k})/\partial \varepsilon ( { \bf k})\big ) } { \omega + { \rm i } \gamma_{1\alpha}({\bf k},\omega ) -q_\alpha^2 v^2_\alpha ( { \bf k})/\omega}. \label{eq50}\end{aligned}\ ] ] in the drude @xmath59 limit , we set @xmath61 in the denominator equal to zero , with @xmath196 and @xmath197 , to obtain @xmath198 an expression known as the generalized drude formula . here , @xmath199 means the average over the fermi surface .
the expression ( [ eq51 ] ) remains valid when more than one band intersects the fermi level , provided that @xmath200 .
the generalization to the case with the relaxation rates @xmath201 dependent on the band index @xmath8 is straightforward .
although we have arrived at ( [ eq41 ] ) and ( [ eq47 ] ) by considering in the auxiliary self - energy @xmath55 only the scattering by boson modes , these expressions are of a more general validity , namely all relevant scattering mechanisms have to be included in @xmath202 .
therefore , we can use ( [ eq41 ] ) to study the effects of local and short - range electron - electron interactions on electrodynamic properties of strongly - correlated low - dimensional systems or to study similar effects associated with the scattering by different periodic cdw / sdw potentials .
in analogy to ( [ eq44 ] ) , we can write the bare interband conductivity in the form @xmath203 with the electron - em field vertex corrections neglected , resulting in @xmath204 . the electron - hole self - energy @xmath205 is replaced by @xmath206 , in the memory - function approximation , and by @xmath207 , in the relaxation - time approximation .    in those multiband systems in which the interband dipole vertex @xmath208 is purely imaginary , the local dipole moment vanishes and there are no local field effects .
@xcite the interband conductivity is given by @xmath209 .
this result can be equally well understood to represent the solution to the rpa equation @xmath210 \rho ( { \bf q},\omega ) = v^{\rm ext } ( { \bf q},\omega ) \label{eq53}\ ] ] for the total dielectric susceptibility @xmath211 , where @xmath212 .    in the opposite case ,
the interband dipole vertex is real and its dispersionless part is directly related to the intracell dipole transitions . particularly interesting
are the examples in which @xmath213 is completely local . in this case
we can use the procedure from the following section to determine the effects of the lorentz local fields .    at this point caution
is in order regarding the relaxation - time approximation .
it is well known that in those systems in which the damping energy @xmath214 is comparable to the interband threshold energy ( the pseudogapped state of the ordinary cdw systems is an example ) the relaxation - time approximation of ( [ eq52 ] ) overestimates the in - gap part of the interband conductivity . in this regime , @xmath215 is found to be better described by the approximate textbook expression @xcite @xmath216 } { \hbar \omega + { \rm i } \hbar \gamma_{2\alpha}({\bf k } ) - \varepsilon_{l'}({\bf k}_+ )   + \varepsilon_{l}({\bf k})}. \label{eq54}\end{aligned}\ ] ] in order to briefly illustrate our results , the typical mean - field single - particle conductivity of the peierls cdw model at temperatures well below the mean - field transition temperature @xmath217 is shown in figure 9 .
the dashed and the dot - dashed lines represent , respectively , the interband contributions calculated by using ( [ eq54 ] ) and the relaxation - time version of ( [ eq52 ] ) . notice the difference between the two curves at @xmath59 .
the solid line is the total single - particle conductivity which is the sum of the contribution of thermally activated charge carriers in ( [ eq51 ] ) and the interband contribution ( [ eq54 ] ) .
this result agrees reasonably well with the experimental spectra @xcite measured at temperatures well below the cdw transition temperature .
the detailed discussion of the mean - field approximation as well as of the effects of fluctuations in the order parameter is given in article ii .
, at temperature @xmath218 k. the bare electron dispersion is @xmath219 , with @xmath220 ev and @xmath221 mev .
the mean - field transition temperature is assumed to be @xmath222 k and @xmath223 mev is the cdw order parameter at temperature @xmath224 k. the damping energies are @xmath225 mev at @xmath226 , and @xmath227 mev at @xmath228 . the dot - dot - dashed line shows the inverse square - root singularity in @xmath229 associated with the limit @xmath230 . ]    in the next section , we shall examine the lorentz local field corrections in the relaxation - time approximation , starting with @xmath231 given by ( [ eq54 ] ) .
the damping energy @xmath232 is assumed to be small in comparison with the interband threshold energy .
most q1d systems are characterized by a large number of atoms in the unit cell ( for example , in the blue bronze k@xmath5moo@xmath3 we have 86 atoms per unit cell ) .
the bands in the vicinity of the fermi level are thus built of different molecular orbitals of the well - defined symmetry with the intermolecular distance large when compared to the first - neighbor atom - atom distance .
therefore , the bands are narrow and the typical interband threshold energy is in the infrared part of the spectrum . in most cases of interest ,
the low - energy interband electron - hole excitations are local and correspond to dipolar transitions between one molecular orbital of the @xmath233 symmetry and another of the @xmath234 symmetry .
consequently , to determine the proper structure of the total conductivity tensor @xmath235 in this case , it is important to take into account the local field effects too .
this is not the only q1d case in which the local field corrections are important . in the ordered state of the strong bond - cdw systems
the interband dipole excitations are also local .
finally , the ordinary @xmath0 models with the interband threshold energy in the visible part of the spectrum also have @xmath236 of the form ( [ eq54 ] ) , with the dipole vertex @xmath237 real .
@xcite      to demonstrate the usual procedure for dealing with the local field corrections , @xcite we consider the model with the intraband currents and with one type of bound charges .
we assume that the contribution of the bound charges is related to the intra - atomic or intramolecular dipolar transitions in @xmath236 . in section 8.2
, we shall briefly discuss the case of infrared - active optical phonons , as well . in this section
, we make a simple modification of the notation . for example
, the label intra is replaced by c and the label inter by b ( for conduction and bound charges ) .
the bare intraband dielectric susceptibility is assumed to be of the form @xmath238 with @xmath239 given by ( [ eq50 ] ) , and the bare interband dielectric susceptibility is @xmath240 with @xmath241 given by ( [ eq54 ] ) .
( a ) and @xmath242 ( b ) . ]
the electrons are assumed to be subjected to an external longitudinal electric field @xmath243 represented by the scalar potential @xmath244 . @xmath245 and @xmath246 are the fourier transforms of two types of macroscopic charge densities . @xmath247 describes the monopole - charge density of conduction electrons and @xmath248 the dipole - charge density of bound charges .
the two - component random - phase approximation describes the long - range screening in such a case . according to figure  10 ,
the rpa equations can be shown in the form @xmath249 \rho_1 ( { \bf q } , \omega )   - v_{12 } \ , \rho_2 ( { \bf q } , \omega ) = v^{\rm ext } ( { \bf q},\omega ) , \nonumber \\ & & \hspace{-7 mm } - v_{21 } \ , \rho_1 ( { \bf q } , \omega ) + \big[1/\chi_{22}^0 ( { \bf q } , \omega )   - v_{22}\big ] \rho_2 ( { \bf q } , \omega ) = v^{\rm ext } ( { \bf q},\omega ) .
\nonumber \\ \label{eq57}\end{aligned}\ ] ] here , the @xmath250 stand for two bare dielectric susceptibilities in ( [ eq55 ] ) and ( [ eq56 ] ) , and @xmath251 are the fourier transforms of different long - range interactions . the numbers @xmath252 , @xmath253 and @xmath254 correspond to the monopole - monopole , monopole - dipole and dipole - dipole interactions , respectively .
@xcite @xmath255 are easily shown to represent the common depolarization factors for the longitudinal @xmath6 fields .
@xcite it is not hard to verify that , after neglecting the lorentz term @xmath256 in @xmath257 ( i.e. , when @xmath258 ) , equations  ( [ eq57 ] ) reduce to ( [ eq53 ] ) .
equations ( [ eq57 ] ) may be rewritten now in the alternative form , in terms of two dipole densities @xmath259 , to give @xmath260 with the implicit dependence on @xmath261 and @xmath140 . here
@xmath262 is the common notation for two bare polarizabilities .
it must be emphasized that the expressions on the right - hand side of ( [ eq58 ] ) represent the macroscopic electric field @xmath122 , equation  ( [ eq30 ] ) , which is the sum of the external electric field and two depolarization fields originating from the conduction electrons and from the induced dipole moments of bound charges .
@xcite this is a useful observation which helps us to rewrite ( [ eq58 ] ) in the final form @xmath263 j^{\rm c}_{\alpha } ( { \bf q } , \omega ) = e_{\alpha}({\bf q } , \omega ) , \label{eq60 } \\
\big [ 1/\alpha^{\rm b}_{\alpha \alpha } ( { \bf q } , \omega ) \big ] p^{\rm b}_{\alpha } ( { \bf q } , \omega ) = e_{\alpha}({\bf q } , \omega ) , \label{eq61 } \end{aligned}\ ] ] with @xmath264 .
\label{eq62}\ ] ] for the thin slab with the longitudinal polarization of the fields , the field @xmath243 is identical to the dielectric displacement field @xmath265 and the relation ( [ eq62 ] ) leads to the usual form of the longitudinal macroscopic dielectric function . if other high - energy on - site dipolar excitations are also included , the total dielectric function becomes @xmath266 here , @xmath267 and @xmath268 are two renormalized polarizabilities .
in the notation of the previous sections , this result can be expressed by ( [ eq10 ] ) and ( [ eq11 ] ) .
we can therefore conclude that the contribution of the monopole charges to the conductivity tensor is decoupled from the contribution of the dipole charges .
the latter contribution is affected by the lorentz field in the usual way , giving also the possibility for the collective contribution in @xmath269 when @xmath270 .
it was implied by our discussion of the quantum transport equation that ( [ eq5 ] ) , together with the density operators ( [ eq3 ] ) and ( [ eq23 ] ) , is able to capture all relevant aspects of the electron - phonon coupling in the dynamical conductivity tensor of multiband electronic systems , at least for the scattering by acoustic and raman - active optical phonons . in order to treat carefully infrared - active optical phonons
, we must also include the direct coupling between phonons and em fields in @xmath27 and take into account the corresponding local field corrections .
the generalization of ( [ eq63 ] ) to such a case with multiple dipolar bound - charge contributions is very simple . for our present purposes it is simplest to replace @xmath271 by the sum of all bare dipolar bound - charge polarizabilities .
in the two - band case with one infrared - active phonon branch , we obtain @xmath272 with @xmath273^{-1}$ ] being the bare phonon polarizability .
in this article , we have derived the dynamical conductivity tensor in a general multiband electronic model with sizeable boson - mediated electron - electron interactions . starting with the gauge - invariant form of the coupling between valence electrons and external em fields , we have derived the general form of the self - consistent bethe
salpeter equations for the auxiliary intraband and interband electron - hole propagators .
it is demonstrated that the intraband bethe
salpeter equation is identical to the usual quantum transport equation and that it reduces to the ordinary transport equation in the weak coupling limit . in this way ,
the present response theory of multiband electronic systems is connected to the generalized drude formula for the intraband conductivity widely used in experimental analyses as well as to different fermi liquid theories based on the landau  silin transport equations . in order to illustrate the formalism presented here
, we shall give in the accompanying article @xcite the detailed description of the cdw mean - field approximation and the effects of fluctuations in the order parameter on the dynamical conductivity of the peierls cdw systems .
this research was supported by the croatian ministry of science and technology under project 119 - 1191458 - 0512 .
massidda s , yu j and freeman a j 1987 _ phys .
lett . _ a * 122 * 198 whangbo m - h and schneemeyer l f 1986 _ inor .
chem . _ * 25 * 2424 mattheiss l f 1985 _ solid state commun . _ * 93 * 791 emery v j 1987 _ phys .
lett . _ * 58 * 2794 pines d and nozires p 1989 _ the theory of quantum liquids i _
( new york : addison - wesley ) ziman j m 1972 _ electrons and phonons _ ( london : oxford university press ) ando y , komiya s , segawa k and kurita y 2004 _ phys .
lett . _ * 93 * 267001 forr l , gal r , berger h , fazekas p , penc k , kzsmrki i and mihly g 2000 _ phys .
lett . _ * 85 * 1937 uchida s , ido t , takagi h , arima t , tokura y and tajima s 1991 _ phys .
7942 degiorgi l , alavi b , mihly g and grner g 1991 _ phys .
_ b * 44 * 7808 kzsmrki i , mihly g , gal r , barii n , akrap a , forr l , homes c c and mihly l 2006 _ phys .
lett . _ * 96 * 186402 sugai s , suzuki h , takayanagi y , hosokawa t and hayamizu n 2003 _ phys . rev .
_ b * 68 * 184504 takigawa m , reyes a p , hammel p c , thompson j d , heffner r h , fisk z and ott k c 1991 _ phys .
_ b * 43 * 247 landau l d , lifshitz e m and pitaevski l p 1995 _ electrodynamics of continous media _ ( oxford : butterworth  heinemann ) kupi i 2009 _ j. raman spectrosc .
_ * 40 * 442 barii s and zelenko j 1990 _ solid state commun . _ * 74 * 367 kupi i and barii s 2007 _ phys .
_ b * 75 * 094508 kubo r , toda m and hashitsume n 1995 _ statistical physics ii _
( berlin : springer - verlag ) eliashberg g m 1961 _ zh . eksp .
fiz . _ * 41 * 410 ( 1962 _ sov .
phys . jetp _
* 14 * 886 ) abrikosov a a 1988 _ fundamentals of the theory of metals _
( amsterdam : nort - holland ) vollhardt d and wlfle p 1980 _ phys .
_ b * 22 * 4666 gtze w and wlfle p 1972 _ phys .
_ b * 6 * 1226 giamarchi t 1991 _ phys .
_ b * 44 * 2905 kupi i 2004 _ physica _ b * 344 * 27 adler s l 1962 _ phys .
rev . _ * 126 * 413 wieser n 1962 _ phys .
rev . _ * 129 * 63 upanovi p , bjeli a and barii s 1996 _ z. phys .
_ b * 101 * 387 kupi i et al unpublished abrikosov a a , gorkov l p and dzyaloshinski i e 1975 _ methods of quantum field theory in statistical physics _
( new york : dover publications ) dzyaloshinskii i e and larkin a i 1973 _ zh . eksp .
fiz . _ * 65 * 411 ( 1974 _ sov .
. jetp _ * 38 * 202 ) kupi i , rukelj z and barii s to be published ( in the main text referred to as article ii ) fetter a l and walecka j d 1971 _ quantum theory of many - particle systems _
( london : mcgraw  hill ) mahan g d 1990 _ many - particle physics _
( new york : plenum press ) solyom j 1979 _ adv .
* 28 * 201 giamarchi t 2004 _ quantum physics in 1d _ ( oxford : clarendon press ) vescoli v , zwick f , henderson w , degiorgi l , grioni m , gruner g and montgomery l k 2000 _ eur .
phys . j. _
b * 13 * 503 mitrovic s , fazekas p , sndergaard c , ariosa d , barii n , berger h , clotta d , forr l , hchst h , kupi i , pavuna d and margaritondo g 2007 _ phys .
_ b * 75 * 153103 ziman j m 1979 _ principles of the theory of solids _
( london : cambridge university press ) wooten f 1972 _ optical properties of solids _
( new york : academic press ) negele j w and orland h 1988 _ quantum many - particle systems _
( new york : addison  wesley ) kupi i 2009 _ phys .
_ b * 79 * 235104 kupi i and barii s 2005 _ fizika _ a * 14 * 47 abrikosov a a and genkin v m 1973 _ zh . eksp .
* 65 * 842 ( 1974 _ sov . phys
. jetp _ * 38 * 417 ) kupi i 2003 _ physica _ c * 391 * 251 lee p a , rice t m and anderson p w 1974 _ solid state commun . _
* 14 * 703 kim k , mckenzie r and wilkins j w 1993 _ phys .
lett . _ * 71 * 4015 ziman j m 1988 _ elements of advanced quantum theory _ ( london : cambridge university press ) millis a j and lee p a 1987 _ phys .
_ b * 35 * 3394 niki h , tuti e and barii s 1995 _ physica _ c * 241 * 247 rozenberg m j , kotliar g and kajueter h 1996 _ phys . rev .
_ b * 64 * 8452 kupi i 2002 _ physica _ b * 322 * 154