the landau criterion  @xcite concerns the threshold velocity @xmath1 of an obstacle moving through stationary superfluids at zero temperature , beyond which excitations are created and superfluidity is lost . by galilean invariance , the criterion also applies to the threshold velocity @xmath1 of a uniformly flowing superfluid ; see , e.g. , ref .
the value of the threshold velocity @xmath1 is determined by the following formula @xmath2 where @xmath3 is the excitation spectrum . for superconductors
the landau criterion gives @xmath4 , where @xmath5 is the pairing gap at zero temperature and zero flow , and @xmath6 is the fermi momentum of the entities that are paired ( with @xmath7 being the associated angular wavenumber ) .
however , it was found by rogers  @xcite that in three spatial dimensions the critical velocity @xmath8the uniform velocity of the superfluid condensate at which the superconducting state becomes unstable  of a clean superconductor exceeds @xmath1 .
furthermore , as discussed , e.g. , in ref  @xcite , for superflow velocities @xmath9 in the range @xmath10 , gapless excitations occur in clean superconductors .
moreover , the ratio @xmath11 has been found to be @xmath12 ( @xmath13 ) .
in contrast , it was found  @xcite that in two dimensions @xmath14 . as for the case of one dimension ,
bagwell reported  @xcite that @xmath15 , as found in two dimensions .
in contrast to the case of three dimensions , in neither one nor two dimensions is gapless superconductivity predicted to occur in the presence of flow .    in this paper , we re - analyze the critical velocity of a clean superconductor in one spatial dimension via a mean - field treatment , and obtain the dependence of this velocity on the temperature @xmath16 , along with the associated temperature - velocity phase diagram .
in particular , we find that at zero temperature the critical velocity is smaller by a factor of @xmath0 than the landau critical velocity @xmath1 .
even though for @xmath17 a gap in the quasiparticle excitation spectrum remains , the superconducting state first becomes unstable there .
this is in contrast to a previous report  @xcite , which held that the critical velocity is @xmath1 , i.e. , the velocity at which the gap in the quasiparticle spectrum closes and excitations proliferate so as to destroy superconductivity .
the smaller value of the critical velocity , which we obtain here , is the result of a pre - emptive clogston - chandrasekhar  like discontinuous phase transition  @xcite , and is analogous to the threshold uniform exchange - field in a superconductor , previously investigated by sarma  @xcite and by maki and tsuneto  @xcite . at low temperatures (
i.e. , for @xmath18 , where @xmath19 is the critical temperature in the absence of flow ) , the transition from the superconducting to the normal state , occurring due to the presence of flow , remains discontinuous .
in contrast , for @xmath20 the transition is continuous , i.e. , it is associated with a continuously vanishing order parameter , e.g. , as @xmath16 approaches the flow - velocity  dependent transition temperature @xmath21 from below .    from the experimental standpoint , it is usually more convenient to control the condensate _
current _ rather than its _ velocity_. ( an exception is provided by a closed loop threaded by magnetic flux . ) the more widely appropriate physical observable is thus the critical current , i.e. , the maximum equilibrium current that a superconductor can sustain , and this is what experiments on superconductors frequently measure .
it should , however , be remarked that recent experiments on trapped condensates of atomic gases make it possible to probe the critical velocity directly  @xcite .
langer and fisher  @xcite introduced a fresh perspective on the issue of the critical velocity of a sueprfluid .
they explained that  as a matter of principle  superflow is inherently unstable , even at velocities below those that allow quasiparticle excitations to proliferate .
this instability results from the possible occurrence of intrinsic , topologically allowed , collective excitation events , in which , e.g. , a vortex ring nucleates locally , traverses an arrhenius energy barrier , and grows without bound , thus eradicating a quantum of flow velocity . in the case of narrow channels , such as the one - dimensional superconductors discussed here
, the fluctuations take the form of phase slip events , which can be either thermal  @xcite or quantum  @xcite .
but , before one takes into account the effects of fluctuations , it is important to understand and clarify the behavior at the mean - field level , which is mainly what we consider in the present paper .
in other words , the work reported here assumes that the rate of such toplogical fluctuations is negligibly small , so that our critical - velocity results , established on the basis of thermodynamic stability and the possibility of quasiparticle production , retain a use . in the present paper we focus on issues of phases of thermodynamic equilibrium ,
competitions between them , and metastability ; we do not attempt to address issues of kinetics , such as the rates at which phase transitions proceed and metastable states evolve into stable ones .
the remainder of the paper is organized as follows . in sec .
[ sec : hami ] we diagonalize the hamiltonian of a superconductor in the presence of flow , and in sec .  [ sec : freeen ] we calculate the free energy of this system . in sec .
[ sec : orderpara ] we use the self - consistency of the order parameter to calculate the dependence of the order parameter on the flow velocity at various temperatures . in sec .
[ sec : transitiont ] we calculate the dependence of the transition temperature on the flow velocity , assuming the transition to be continuous , i.e. , in the limit of linear instability .
double solutions for the transition temperature turn out to exist for larger values of the flow velocity .
we resolve this issue in sec .
[ sec : phase ] by identifying the globally stable solution , i.e. , the one that corresponds to the lower free energy . in sec .
[ sec : super ] we address the existence of metastable solutions , and discuss the extent to which these provide a superflow - based analog of the phenomena of supercooling and superheating .
we also obtain the analog of the superheating limit , together with the equilibrium phase boundary and the supercooling limit .
having ascertained the true , equilibrium , velocity - dependent order parameter , we calculate in sec .
[ sec : current ] the dependence of the supercurrent on the flow velocity , and hence determine the critical current and superfluid density . in sec .
[ sec : disorder ] we briefly discuss the effect of disorder on the nature of the transition . in sec .
[ sec : higherd ] we contrast the critical - velocity results for one dimension with those for two and three dimensions .
we conclude in sec .
[ sec : conclude ] .
we shall use a microscopic approach to discuss an effectively one - dimensional superconducting system , in which there is spatially uniform flow at velocity @xmath22 . here , @xmath23 is the electron mass and @xmath24 is the corresponding momentum associated with the flow ( or equivalently , the angular wavenumber , as we shall set @xmath25 ) .
( what we mean by `` effectively one - dimensional '' is that the transverse dimensions are much smaller than the zero - temperature coherence length @xmath26 . )
we begin by writing down a hamiltonian that is equivalent to eq .
( iii.4 ) of rogers  @xcite and eq .  ( 1 ) of nozires and schmitt - rink  @xcite : @xmath27 where @xmath28 and @xmath29 are respectively the annihilation and creation operators of electrons of spin - projection @xmath30 at ( one - dimensional ) momentum @xmath31 , @xmath32 is the bcs pairing potential  @xcite , and @xmath33 is the chemical potential . following rogers  @xcite , we assume that @xmath32 depends on the difference between @xmath31 and @xmath34 , and hence is independent of @xmath24 . to justify our ignoring any @xmath24 dependence
, we note that the critical momentum @xmath24 is no bigger than roughly @xmath35 , and this is roughly @xmath36 , i.e. , several orders of magnitude smaller than @xmath7 , where @xmath37 is the atomic spacing and @xmath26 is the zero - temperature superconducting coherence length .
this form of hamiltonian anticipates that any cooper pairs are formed via the pairing of electrons in states @xmath38 and @xmath39 , and that the resulting pairs have center - of - mass momentum @xmath40 .
next , we make the mean - field approximation and adopt the bcs form for @xmath32  @xcite , thus arriving at the hamiltonian @xmath41 in which @xmath42 is the magnitude of @xmath32 in the momentum range for which the pairing potential is nonzero , @xmath43 , and @xmath44 is the self - consistency condition on the order parameter ( and similarly for @xmath45 ) .
this mean - field hamiltonian can be solved via the bogoliubov - valatin transformation @xmath46 with @xmath47 where , for convenience , we have made the definitions @xmath48 and @xmath49 .
we choose @xmath50 to be real and non - negative , and thus may drop the absolute value on @xmath50 in the definition of @xmath51 . with this procedure ,
the hamiltonian becomes diagonal : @xmath52 this expression also holds for higher dimensions , provided we replace @xmath31 by @xmath53 , @xmath9 by @xmath54 , and @xmath55 by @xmath56 .
as the hamiltonian has been diagonalized , we can readily evaluate the helmholtz free energy @xmath57 of the system : @xmath58,\ ] ] where @xmath59 and @xmath60 are fermi distribution functions , defined via @xmath61 @xmath62 indicates an average , weighted by the equilibrium density matrix , and @xmath63 is the inverse temperature ( with boltzmann s constant set to unity ) . we therefore arrive at the result : @xmath64 one may equally well calculate the partition function @xmath65 , so as to obtain the free energy via @xmath66 .
the order parameter is to be determined self - consistently , via @xmath67 we note that because the momentum sum runs over both positive and negative values one , can safely replace @xmath68 by @xmath69 ( or vice versa ) in this equation , and thus arrives at the result : @xmath70 which is identical to eq .
( 21 ) of ref .
this equation can also be derived by demanding stationarity of the free energy , i.e. , @xmath71 .
the above results , eqs .
( [ eqn : freeen0 ] ) to  ( [ eqn : self_con2 ] ) hold for two- and three - dimensional systems as well , provided we replace @xmath31 by @xmath53 , @xmath9 by @xmath54 , and @xmath55 by @xmath56 , but in the following we shall focus mainly on the case of one dimension .
from eq .  ( [ eqn : self_con1 ] ) , or equivalently eq .  ( [ eqn : self_con2 ] ) , we can determine the values of the order parameter at any temperature and superflow velocity .
of course , there is always the trivial solution @xmath72 , which corresponds to the normal state .
however , in this section our focus is on nontrivial solutions .
let us first examine the limit of zero temperature . in this limit , the self - consistency condition  ( [ eqn : self_con2 ] ) becomes @xmath73 where @xmath74 is the heaviside step function . by linearizing the @xmath31-dependent spectrum around @xmath75 , and exchanging @xmath76 for the @xmath77 value of the order parameter @xmath78
, we arrive at the following result for @xmath79 : @xmath80 this condition yields two branches of solutions for the order parameter : ( 1 ) if @xmath81 , we obtain @xmath82 ; on the other hand , ( 2 ) if @xmath83 , the condition becomes @xmath84 which leads to @xmath85 .
these two zero - temperature solutions were first obtained by sarma  @xcite in connection with the exchange - field effect in superconductors . in fig .
[ fig:1 ] we show these two zero - temperature solutions [ i.e. , the horizontal and parabolic curves emanating from the upper right point @xmath86 , together with the corresponding solutions at nonzero temperatures . by continuity ,
it is not surprising that there is a low - temperature regime in which there is multivaluedness in the solutions for @xmath50 , as we shall see in the following subsection and in fig .
[ fig:1 ] .      for these , we can solve eq .
( [ eqn : self_con1 ] ) , or equivalently eq .  ( [ eqn : self_con2 ] ) , numerically , and obtain the velocity - dependent order parameter @xmath87 at arbitrary temperatures , as illustrated in fig .  [ fig:1 ] for several values of the temperature .
we remark that in solving for the order parameter we have linearized the spectrum @xmath88 about fermi momentum @xmath7 ; this linearization is valid as long as @xmath89 . from fig .
[ fig:1 ] it is evident that over a certain ( higher ) range of velocities @xmath9 there are two solutions for @xmath50 , these solutions differing from those obtained by bagwell  @xcite , by whom multiple solutions were not found .
this type of multiplicity feature was first observed by sarma  @xcite in the context of the exchange field in superconductors .
in fact , the self - consistency conditions holding in the exchange - field case are identical to the ones holding here , provided one makes an identification between @xmath90 and the exchange - field energy @xmath91 .
the two situations appear to be essentially equivalent .    the feature of double solutions that we have found for the order parameter , as shown in fig
.  [ fig:1 ] , is in contrast with previous findings  @xcite ( see , in particular , fig .  3 therein ) .
we suspect that this discrepancy results from the use of an iterative scheme for solving the self - consistency equation numerically ; when multiple solutions exist , such scheme yields results that depend sensitively on the initial conditions .
we have instead used the bisection method , which locates all possible solutions .
if one assumes that the transition between the superconducting and normal states is continuous , one can determine the transition temperature via the self - consistency condition taken in the limit of vanishing @xmath50 .
we note that , owing to the relation of the fermi function to the matsubara sum ( see , e.g. , ref .
@xcite ) , eq .  ( [ eqn : self_con1 ] ) is equivalent to @xmath92 where @xmath93 with @xmath94 are the matsubara frequencies .
we denote by @xmath95 the value of the temperature that solves eq .
( [ eqn : selfcon3 ] ) in the limit @xmath96 . in this limit , imposing a cutoff @xmath97 on the mastubara frequencies , and integrating out the momentum @xmath31 , we arrive at the condition @xmath98 in which we have redefined @xmath99 to mean @xmath100 . by exchanging @xmath76 for the zero - flow critical temperature @xmath19 and performing the summation , eq .
( [ eqn : tc ] ) becomes @xmath101 where @xmath102 is the di - gamma function  @xcite .
the solutions of eq .
( [ eqn : tc0 ] ) for @xmath95 are shown in fig .
[ fig : tc0_v ] , which , like fig
.  [ fig:1 ] , exhibits a double - solution feature for a range of velocities .
the solution at @xmath103 and @xmath104 corresponds to the branch @xmath105 .
as we shall see below , this branch corresponds to an _
unstable _ solution ; a correct description of the @xmath95 vs.  @xmath9 phase diagram requires the consideration of free energies .
in the light of the results obtained in sec .  [ sec : transitiont ] , we ask the question : what is the true transition temperature ? to determine this , one has to take into account the multiple solutions of the order - parameter self - consistency condition and compare their free energies ( e.g. , see ref .
@xcite ) . with this prescription
we ask : will our transition temperature @xmath21 have the same functional form as the @xmath106 of sarma  @xcite and of maki and tsuneto  @xcite ( smt ) , in which @xmath107 is the magnetic exchange field ?    in the case considered by smt , a comparison is made between the free energies of a spin - unpolarized superconducting state and a normal state that is partially spin - polarized , both states being subject to an exchange field . in the present case
the comparison is between a flowing superconducting state and a stationary normal state .
viewed from a reference frame that flows with the superconducting state , the normal state flows and can be regarded as the analog of smt s polarized state .
furthermore , our superconducting state is the analog of smt s superconducting state .
this perspective suggests that the correct procedure for our purposes is to compare the free energies of the flowing superconducting state and the stationary normal state , as was first done by bagwell  @xcite .
the physical reason for the equivalence between the exchange - field case considered by sarma and the superflow case considered here is as follows . in the case of sarma , due to the exchange field , the up and down electrons in a cooper pair have an energy difference of @xmath108 . in the case of flow ,
the electron pairing is between states @xmath109 and @xmath110 , and which have an energy difference of @xmath111 near the fermi surface . in the former case
, the normal state can be polarized , and thus can reduce its energy by an amount proportional to @xmath112 . in the latter case
, the stationary normal state has an energy that is lower than the flowing normal state by an amount proportional to @xmath113 .
hence , the two scenarios appear to be equivalent .    to understand the root of our state - selection procedure , we begin with an analogy .
consider the liquid and crystalline states of a particular material .
following callen  @xcite , we observe that the information sufficient to specify uniquely an equilibrium state is greater for the crystalline state than for the liquid state , owing to the spontaneously broken translational symmetry of the crystal and the resulting low - energy goldstone field , i.e. , the displacement field .
thus , when assessing the relative thermodynamic stability of a liquid and a crystalline state , one must specify a displacement field for the crystal but one must not ( and indeed can not ) do so for the liquid state . it is meaningful to speak of an equilibrium crystalline state that supports a specified static shear stress , but not so for a liquid , as stress would induce a nonequilibrium dissipative steady state of shear flow , in which entropy would constantly be being produced . in summary
, there is a unique equilibrium state for the liquid , but there is a family of equilibrium states for the crystal , differing in their displacement field ( or strain field , or its thermodynamic conjugate , the stress field ) .
one is entitled to consider the relative thermodynamic stability of a liquid and any one of this family of crystalline states .
returning now to the case of superconductivity , the precise analogy is on the one hand between the crystal and the superconductor , and on the other hand between the liquid and the normal metal .
one is thus entitled to consider the relative thermodynamic stability of a stationary normal state and any one of the family of the flowing superconducting states . that there is a family of flowing superconducting states originates in the spontaneously broken gauge symmetry of the superconducting state , which amounts to a controllable thermodynamic field , ( i.e. , the phase field ) , which is the analog of the displacement field of the crystal .
( the analogy is : phase @xmath114 displacement ; velocity @xmath114 strain ; current density @xmath114 stress.)by contrast , the normal state is unique .    to obtain the phase diagram
, one should then compare the free energy @xmath115 of a flowing superconducting state to the free energy @xmath116 of a stationary ( i.e. , not flowing ) normal state  @xcite .
the energy difference between the flowing and stationary normal states plays the role of the paramagnetic energy .
we shall see that the true equilibrium transition temperature then has exactly the functional form obtained by sarma  @xcite and by maki and tsuneto  @xcite .
the free energy  ( [ eqn : freeen ] ) of the superconducting state at a fixed chemical potential @xmath33 is given explicitly by @xmath117.\ ] ] the stationary , normal state at chemical potential @xmath118 thus has free energy @xmath119.\ ] ] our goal is to find the difference between two energies at the same chemical potential , which is given by @xmath120-\left[f_n(0,\mu)-f_n\big(0,\mu-\frac{q^2}{2m}\big)\right],\ ] ] where the term in the second pair of square brackets equals @xmath121 ( i.e. , the analog of the paramagnetic energy ) , with @xmath122 being the total number of electrons .
hence , from eqs .
( [ eqn : freeenn0 ] ) and  ( [ eqn : difff ] ) we arrive at the result that @xmath123\nonumber\\ & = & n\frac{q^2}{2 m } + \frac{\delta^2}{g}-t\sum_k \ln\left[\frac{\cosh\beta e_k+\cosh\beta{kv}}{\cosh\beta \overline{\veps}_k+1}\right ] .
\label{eqn : fsfn}\end{aligned}\ ] ]    what we should do to obtain the phase diagram is to compare the possible nonzero solution(s ) for the order parameter to determine which solution possesses the lower free energy @xmath115 and whether this free energy is lower than that of a stationary ( i.e. , nonmoving ) normal state , i.e. , whether @xmath124 . for convenience ,
let us make the definition @xmath125.\ ] ] the solution for @xmath50 that has the largest value of @xmath126 , provided this difference is greater than zero , is the true equilibrium solution ; otherwise the normal state will be the true equilibrium state .
the kinetic flow term @xmath127 can be re - expressed ( by relating the total electron number to the density of state ) as @xmath128 , where @xmath129 is the density of states per spin , and we have , for convenience , set the total `` volume '' of the system to be unity .
note that from eq .
( [ eqn : dfq2 ] ) we have at zero temperature @xmath130 this expression also holds for higher - dimensional systems , provided we replace @xmath31 by @xmath53 , @xmath9 by @xmath54 , and @xmath55 by @xmath56 .
the effective condensation energy @xmath131 for the solution @xmath82 [ see the discussion below eq .
( [ eqn : orderparameter1 ] ) ] is @xmath132 . balancing @xmath131 against the additional energy @xmath133 due to the flow , we obtain @xmath134 .
this corresponds to the clogston - chandrasekhar limit in the exchange - field case  @xcite .    in fig .
[ fig : freeen ] , we plot @xmath131 versus the momentum @xmath24 for the two superconducting solutions at several temperatures , as well as @xmath135 .
the crossing of the curve representing @xmath135 by the curve representing the superconducting state signifies a discontinuous transition . by finding such transitions at various values of @xmath9 ,
we arrive at the equilibrium phase boundary @xmath21 , shown as the curve abd in fig .
[ fig : sfphase ] .
as discussed in sec .  [ sec : phase ] , we have discussed the true equilibrium phase transition curve abd ; see fig .
[ fig : sfphase ] .
but what is the physical meaning of the branch bc ?    it would be a limit of metastability if a normal state with flow were a true equilibrium state .
suppose that the system were in an equilibrium normal state with @xmath136 ( i.e. , lying to the right of the equilibrium phase boundary abd ) , and were then ` quenched ' by the rapid change of @xmath16 and/or @xmath9 into the metastable region bcd .
the system would remain in the normal state , but metastably so , until a fluctuation were to occur that would nucleate a droplet of the superconducting phase large enough to grow and complete the conversion of the state to the true equilibrium state , i.e. , the superconducting state . the metastability limit ( which is also known as a spinodal line )
is formally obtained by solving for the temperature obeying the equation [ c.f .
( [ eqn : self_con1 ] ) ,  ( [ eqn : selfcon3 ] ) and  ( [ eqn : tc ] ) ] @xmath137 however , as we argue in sec .
[ sec : phase ] , a normal state with a constant flow is not a true equilibrium state ; therefore , strictly speaking , the curve bc does not represent a meaningful supercooling curve .
on the other hand , the branch be is a true limit of metastability for the superconducting state
.    the system can be ` quenched ' into a metastable superconducting state from an equilibrium superconducting state [ i.e. , from @xmath138 to the left of the equilibrium phase boundary ] by rapidly changing @xmath138 to a value in the region bed .
the system will remain superconducting , but metastably so , until a normal - phase nucleation event carries it to the equilibrium ( i.e. , normal ) state . in this case , the metastability limit is determined by simultaneously solving the equations @xmath139 and @xmath140 ( but @xmath141 ) numerically ; see fig .  [ fig : sfphase ] .
we note that a corresponding diagram in the case of exchange - field effect in superconductors was first obtained by maki and tsuneto  @xcite . to answer the question of how long it would take for system in a metastable state to find
its equilibrium state is a kinetic one that would require further investigation .
in this section we examine the dependence of the supercurrent on the flow velocity , and thereby determine the critical current , i.e. , the maximum equilibrium supercurrent that the system can sustain at various temperatures . in general , the charge current @xmath142 carried by the system is defined via @xmath143 as @xmath144 is the total number of electrons @xmath122 , we can re - write eq .
( [ eqn : iq ] ) as @xmath145 where @xmath146 . by using eqs .
( [ eqn : bogovale ] ) and  ( [ eqn : f12 ] ) , the second term can be simplified to @xmath147 , and thus the expression for @xmath142 becomes @xmath148 where the electron charge is @xmath149 and the temperature @xmath16 is implicit in the fermi function @xmath68 ; see eq .
( [ eqn : f12 ] ) .
we shall compute the current in the superconducting state as a function of @xmath16 and @xmath9 , and then , by choosing the value @xmath150 of @xmath9 in the range @xmath151 $ ] that maximizes @xmath152 at fixed @xmath16 , we shall obtain the critical current , denoted by @xmath153 $ ] , at that temperature .
equation  ( [ eqn : iq2 ] ) , applied at @xmath154 , tells us that the zero - temperature current is given by @xmath155,\\ & = & i_q^0\left[\frac{v}{v_l } - \sqrt{\frac{v^2}{v_l^2}-1}\,\,\theta(v - v_l)\right],\end{aligned}\ ] ] where @xmath156 .
although , formally , @xmath157 achieves its maximum value ( viz . , @xmath158 ) at @xmath159 , this maximum value is unattainable , because the superconducting state gives way to the normal state via a pre - emptive transition at @xmath160 ( unless the system falls out of equilibrium and remains metabstably in the superconducting state ) .
therefore , the maximum equilibrium critical current at zero temperature occurs at @xmath17 , and has the value @xmath161 . at @xmath17
the superconducting state first becomes unstable , so the critical current that is attainable in practice is slightly less than @xmath161 .
this is unlike the behavior in three dimensions , as well as the higher of the nonzero temperatures ; for these situations the maximum current is achieved at a value of superfluid velocity for which the superconducting state is still globally stable .    in fig .
[ fig : ivsv ] we show the equilibrium supercurrent @xmath142 vs. the superfluid velocity  @xmath9 for various values of the temperature . the value of @xmath162 that corresponds to a maximum in the current ( i.e. , the critical current @xmath163 ) is not necessarily @xmath164 .
( it _ is _
@xmath8 for lower temperatures , but not for higher ones . ) the temperature dependence of the critical current @xmath165 is shown in fig .  [
fig : icvst ] .
a related quantity of interest is the superfluid density @xmath166 , defined as @xmath167 i.e. , the response of the current density @xmath168 ( whose expression is identical to @xmath142 because we have set the system ` volume ' to unity ) to an infinitesimal change in the ` driving ' velocity @xmath9 .
the temperature dependence of @xmath166 is shown in fig .
[ fig : nsvst ] . the vanishing of @xmath166 near @xmath95 is ,
as expected , linear , and the departure of @xmath166 from its zero - temperature value has , at low temperatures , a thermally activated form .
so far , we have discussed clean systems , i.e. , systems not disordered by any scattering electrons by impurities
. what effect will disorder have ?
in particular , would the disorder change the nature of the superconducting - to - normal transition ? in this section
we briefly discuss the latter issue . in the case of disorder ,
the green function technique is better suited for dealing with disorder than is trying to diagonalize the hamiltonian with inclusion of an arbitrary configuration of the impurities .
an equivalent calculation in the context of the exchange field was done by maki and tsuneto  @xcite and , more recently , by wei and goldbart  @xcite in the context of little - parks effect in small rings . here , in contrast with ref .
@xcite , we are not concerned with the finiteness of the system size , and we simply quote the relevant results for the equation obeyed by the critical temperature , i.e. , eq .
( b39 ) of ref .
@xcite : @xmath169,\end{aligned}\ ] ] where , in the present case , @xmath170 and @xmath171 , the parameter @xmath172 is the elastic mean - free time , and @xmath102 is the di - gamma function  @xcite .    in the clean limit ( i.e. , @xmath173 ) , eq .  ( [ eqn : arbitrarydisorder ] ) reduces to eq .
( [ eqn : tc0 ] ) , for which multiple solutions for @xmath21 exist . for strong disorder ( @xmath174 ) ,
( [ eqn : arbitrarydisorder ] ) reduces to @xmath175 for which no multiple solutions for @xmath21 exist .
this suggests that the assumption of a vanishing order parameter in the search for the critical temperature is valid and , hence , that the transition is continuous .
moreover , from eq .
( [ eqn : tcdisorder ] ) , in the strong - disorder limit the critical velocity is determined via @xmath176 numerically , we have found that for @xmath177 multiple solutions do not occur .
this , then , gives the threshold for the disorder strength that divides the discontinuous and continuous transition regimes . in fig .
[ fig : tvsvdisorder ] we show the solutions of eq .
( [ eqn : arbitrarydisorder ] ) for various values of @xmath178 . in the strong - disorder limit
we obtain from eq .
( [ eqn : vcdisorder ] ) that @xmath179 ( where @xmath180 is the elastic mean - free path and @xmath26 is zero - temperature coherence length ) , in particular that the critical velocity exceeds @xmath1 .
but will the critical _ current _ also exceed its clean - limit value ? according to bardeen  @xcite ,
the superfluid density in the strong - disorder limit is reduced by a factor of roughly @xmath181 .
an order - of - magnitude estimate gives for the critical current the product of the superfluid density @xmath182 ( where @xmath183 denotes the superfluid density in the clean limit ) and the critical velocity @xmath184 .
this , in turn , gives that the critical supercurrent is reduced by roughly a factor @xmath185 , which is much less than unity in the strong - disorder limit .
therefore , even though the critical velocity is increased by disorder , the critical current is reduced by it .
we now digress to discuss the scenario of a closed superconducting loop threaded by magnetic flux , which is pertinent to the experiment performed by liu et al .
@xcite and the situation considered theoretically in ref .
@xcite . there ,
the boosted velocity is induced by the magnetic flux .
the measurement of @xmath95 vs. flux @xmath186 ( with @xmath187 ) presented in fig .  4 of ref .
@xcite shows non-@xmath188-periodic behavior , with the first dome ( near @xmath189 ) being higher than the second dome ( near @xmath190 ) .
naturally , the data explore only a limited range of flux , so it is not possible to determine with certainty whether higher flux values would yield dome heights that oscillate with flux or decay .
oscillations would be a consequence of flux through the hole ; decay would be a consequence of flux through the sample , which would cause pair breaking . to ascertain which of these two effects ( oscillation or decay ) is likely to dominate in the range of fluxes probed by the experiment , it is useful to compare several lengthscales .
first , the sample radius @xmath191 is smaller than the coherence length @xmath26 , by a factor of roughly @xmath192 . following ref .
@xcite , this suggests the smallness of the radius as the origin of the dome - height mismatch . on the other hand ,
a more complete characterization of the system would take into account the elastic mean - free path @xmath180 , which for the sample studied in ref .
@xcite is shorter than the radius , by a factor of approximately @xmath193 .
this puts the sample in the dirty regime , and as a result , the impact of the smallness of the radius is suppressed and tends to restore the @xmath188 periodicity of @xmath95 discussed in ref .
@xcite and given by @xmath194 where @xmath195 and @xmath196 . what else could cause the lack of @xmath188 periodicity ? because the magnetic flux is not confined inside the hole in the loop , but is also present throughout the sample , owing to its finite thickness , the orbital pair - breaking effect comes into play
. to account for this effect , we use the scenario of the orbital pair - breaking effect of a parallel magnetic field applied to a film of finite thickness @xmath197 ( which is about @xmath198 in the experiment ) , which leads us to add to the argument of the second term on the r.h.s .  of eq .
( [ eqn : tcdisorder ] ) the term @xmath199 see , e.g. , ref .  @xcite .
( the term should be correct up to a numerical factor due to geometry . )  as we see from fig .
[ fig : liu ] , the resulting dependence of the critical temperature on the flux @xmath200 seems to reproduce , at least qualitatively , the features reported by liu et al .  @xcite .
we have seen that in the clean limit the discontinuous transition at flow velocity @xmath201 occurs due to a competition between the condensation energy @xmath202 and the flow kinetic energy @xmath203 ( for one dimension ; see footnote  @xcite ) , where @xmath204 . the principle that we employed in sec .  [
sec : phase ] is to compare the free energies of the flowing superconducting states and of the stationary normal state . as further confirmation of the validity of this principle
, we now apply it to higher - dimensional settings , and compare the results with the known one mentioned in the introduction .    in two dimensions ,
the solution of the self - consistency equation  ( [ eqn : self_con2 ] ) gives @xmath82 for @xmath205 , and @xmath72 for @xmath206 .
this indicates a discontinuous transition at @xmath207 .
indeed , this result is consistent with free - energy considerations : by using the condensation energy @xmath208 and the flow kinetic energy @xmath209 ( the difference from the one - dimensional case arising from the density of states ; see footnote  @xcite ) , we conclude that @xmath210 .    in three dimensions
the solution of the self - consistency equation  ( [ eqn : self_con2 ] ) gives @xmath82 for @xmath211 , and for @xmath212 ( see ref .
@xcite ) gives the following implicit equation for @xmath50 : @xmath213=\ln \left(\frac { v}{v_l}\right),\end{aligned}\ ] ] where we have ignored terms of order @xmath214 or smaller ( with @xmath97 denoting the debye frequency ) .
equation  ( [ eqn:3ddv ] ) gives that @xmath72 at @xmath215 .
is this linear - instability result consistent with free - energy considerations ( i.e. , the balancing of the two energies ) ? if one were to navely use the condensation energy @xmath216 to balance the flow energy @xmath217
@xcite , one would obtain @xmath218 !
does this mean that the superconducting state really becomes unstable at @xmath219 ?
such a result would be in conflict with the result , mentioned in sec .
[ sec : intro ] , that @xmath220 ?
how is this apparent conflict resolved ? for @xmath205 , the effective condensation energy @xmath131 is indeed @xmath221 .
but for @xmath206 , @xmath222 , so it would be incorrect to equate @xmath221 and @xmath223 in order to deduce @xmath8 . instead , one should use the expression for @xmath131 appropriate to the range @xmath206 .
in fact , an evaluation of eq .
( [ eqn : fzerot ] ) in three dimensions gives ( for @xmath206 ) @xmath224 ^ 3 \nonumber\\ & & + n_0\ , \delta^2 \sinh^{-1}\frac{\sqrt{(k_f v)^2-{\delta}^2}}{\delta}-n_0\,\delta^2\ln\frac{\sqrt{(k_f\ , v)^2-{\delta}^2}+k_f\ , v}{\delta},\end{aligned}\ ] ] where the dependence of @xmath50 on @xmath9 is given by eq .
( [ eqn:3ddv ] ) . from eqs .
( [ eqn : fzerot ] ) and  ( [ eqn : df3d ] ) we have that for @xmath225 , @xmath226 and @xmath227 . only when @xmath9 increases to @xmath228 ( i.e. , when @xmath72 ) do we have @xmath229 .
thus , the transition at @xmath230 is continuous , as the order parameter vanishes continuously there , and , hence , the conflict is resolved .
the dependence of the quantity @xmath231 on flow velocity @xmath9 is plotted in fig .  [
fig : dfq2vsv ] , which clearly shows that the superconducting state is stable for flow velocities in the range @xmath232 .
we have revisited the issue of the critical velocity of an effectively one - dimensional superconductor using mean - field theory .
as we have discussed , this issue is equivalent to that of the threshold uniform exchange - field in a superconductor , first investigated by sarma  @xcite and by maki and tsuneto  @xcite .
this equivalence is due to the correspondence of ( i )  the zeeman frustration energy @xmath108 between the two electrons of a cooper pair ( in the case of exchange field ) , and ( ii )  the kinetic frustration energy @xmath233 ( in the case of flow ) . in particular
, we have found that at zero temperature the critical velocity is a factor of @xmath0 smaller than the landau critical velocity , in contrast to a previous finding  @xcite .
our result originates in a clogston - chandrasekhar  like discontinuous transition between the superconducting and the normal state .
the physical reason for the discontinuity of the transition is that it results from a balance between the condensation and flow energies .
( recall that the bulk superconducting phase transition in magnetic field results from a balance between the condensation ( @xmath234 ) and magnetic field ( @xmath235 ) energies , and that the transition is discontinuous . )
the transition remains discontinuous for temperatures below approximately @xmath236 ; above that , it is continuous .
we have also studied the issue of critical supercurrents , determined the current - temperature phase diagram , and examined metastability ( and its limits ) in the temperature versus flow - velocity phase diagram . as a test of our underlying principle , namely the comparison of the free energies of flowing superconducting and stationary normal states , we have also examined the two- and three - dimensional cases , and have thus obtained results that are in agreement with previous findings .    in addition , we have commented on the effects of electron scattering by impurities and , in particular , we have argued that strong disorder renders continuous the aforementioned discontinuous phase transition . even though disorder can increase the value of the critical _ velocity _ , the physically measurable quantity , i.e. , the critical _ current _ ,
is still reduced by the disorder .    throughout this paper , we have adopted a mean - field approach , thus ignoring the effects of fluctuations such as phase slips .
furthermore , we have limited ourselves to the analysis of equilibrium properties as well as issues of metastability
. the intriguing issue of how a flowing superconducting state undergoes the transition to the stationary normal state , which essential dynamical processes take place , and on what timescale all remain as research directions for the future .
_ acknowledgments_. the authors acknowledge useful discussions with anton burkov , tony leggett and frank wilhelm .
this work was supported by iqc , nserc , and orf ( tcw ) and by doe grant no .
defg02 - 07er46453 through the frederick seitz materials research laboratory at the university of illinois ( pmg ) .
99 l. d. landau , j. phys .
5 * , 71 ( 1941 ) .
e. varoquaux , c. r. physique * 7 * , 1101 ( 2006 ) .
k. t. rogers , ph.d .
thesis , university of illinois , 1960 .
j. bardeen , rev .
* 34 * , 667 ( 1962 ) .
for example , see a. m. zagoskin , _ quantum theory of many - body systems _ ( springer - verlag , new york , 1998 ) .
p. f. bagwell , phys .
b * 49 * , 6841 ( 1994 ) .
a. m. clogston , phys .
* 9 * , 266 ( 1962 ) ; b. s. chandrasekhar , appl .
. lett . * 1 * , 7 ( 1962 ) .
g. sarma , j. phys .
solids , * 24 * , 1029 ( 1963 ) .
k. maki and t. tsuneto , prog .
* 31 * , 945 ( 1964 ) . c. raman , m. khl , r. onofrio , d. s. durfee , c. e. kuklewicz , z. hadzibabic , and w. ketterle , phys .
lett . * 83 * , 2502 ( 1999 ) ; d. e. miller , j. k. chin , c. a. stan , y. liu , w. setiawan , c. sanner , and w. ketterle , phys .
lett . * 99 * , 070402 ( 2007 ) .
j. s. langer and m. e. fisher , phys .
* 19 * , 560 ( 1967 ) . j. s. langer and v. ambegaokar , phys
. rev . * 164 * , 498 ( 1967 ) .
a. d. zaikin , d. s. golubev , a. van otterlo , and g. t. zimanyi , phys .
* 78 * , 1552 ( 1997 ) .
p. nozires and s. schmitt - rink , j. low temp .
phys . * 59 * , 195 ( 1985 ) .
j. bardeen , l. n. cooper , and j. r. schrieffer , phys . rev .
* 106 * , 162 ( 1957 ) .
m. abramowitz and i. a. stegun , _ handbook of mathematical functions _ , ( dover , new york , 1972 ) , pp .
259 and pp . 556 .
wei and p. m. goldbart , phys .
b * 77 * , 224512 ( 2008 ) . because the expression for the density of states depends on spatial dimension , we have that : ( a ) in one dimension , @xmath237 ; ( b ) in two dimensions , @xmath238 ; and ( c ) in three dimensions , @xmath239 .


understanding strongly correlated electron systems represents a major challenge in solid state physics .
dynamical mean - field theory ( for a review see ref . ) has emerged as a powerful tool to address this challenge . within this approach
, the lattice problem for an interacting electron system is mapped onto a quantum impurity problem with a bath  determined self - consistently@xcite . in its simple form ,
the construction becomes exact in the limit of infinite lattice coordination as shown in ref . .
in addition , several generalizations of the approach , such as the extended dmft@xcite and the dmft treatment of correlated hopping@xcite , requiring different scalings of the parameters in the hamiltonian , have been presented .
more generally , dmft is a tool with strong computational power , and provides valuable insight into the physics of strongly - correlated electrons . in the case of simple model hamiltonians with purely local interactions , such as the hubbard model
, it was able to provide a consistent description for several non - perturbative properties such as , for example , the mott - hubbard metal insulator transition@xcite .
currently , many efforts are made to extend the dmft approach in order to make it suitable in addressing more realistic problems .
the purpose of this work is to develop a strong coupling approach for interacting lattice problems . within this approach ,
we construct a systematic generalization of the dmft technique capable of describing any non - local term of a generic lattice hamiltonian .
the present formulation has several advantages .
first , it unifies various dmftschemes .
for example , the treatment of correlated hopping@xcite , the single site dmft  of the hubbard model ( see section iii ) , and the extended dmft@xcite are all limiting cases of this unified approach .
second , besides unifying different dmft schemes , this formulation is sufficiently general to be a starting point for expansions around dmft .
one possibility is described in the appendix .
in addition , the formulation in terms of hubbard operators makes it an ideal framework for a dmft  based renormalization group approach along the lines of ref . .
third , being based on an expansion around the atomic limit , it provides new and valuable insight into the diagrammatics underlying cluster dmft schemes , which will be a subject of a subsequent publication .
the central result of this work is the formulation of a perturbation theory , for a general hamiltonian of the form    @xmath0    containing a local term , @xmath1 , and a non - local one , @xmath2 . in terms of hubbard operators , @xmath3
, we can write the two contributions as @xmath4and @xmath5 in these equations @xmath6 and @xmath7 represent single - site states . for simplicity , we consider the case of spin @xmath8 fermions , so that for each site @xmath9 there will be four states @xmath10 : @xmath11 ( empty site ) , @xmath12 ( single occupied site with spin up ) , @xmath13 ( single occupied site , spin down ) , and @xmath14 ( double occupied site ) .
the hubbard operators @xmath3 describe transitions between these states .
these operators have a fermionic character if the occupation numbers of the states @xmath15 and @xmath16 differ by one , and bosonic character otherwise .
the algebra of the hubbard operators is defined by the multiplication rule @xmath17together with the conserving condition @xmath18and the commutation relation @xmath19_{\pm } = \delta _ { ij}(\delta _ { \beta \beta ^{\prime } } x_{i}^{\alpha \alpha ^{\prime } } \pm \delta _ { \alpha \alpha ^{\prime } } x_{i}^{\beta ^{\prime } \beta } ) .
\label{oper3}\]]in eq .
( [ oper3 ] ) the anticommutator is used if both operators are fermionic and the commutator otherwise .
the canonic fermion operators @xmath20 can be expressed in terms of hubbard operators as @xmath21the parameters @xmath22 in eq .
( [ ham2 ] ) represent the single - site energies and will be determined , in general , by the chemical potential , @xmath23 , the on - site hubbard interaction , @xmath24 , and the external fields .
the coupling constants @xmath25 may include contributions coming , for example , from hopping , @xmath26 , spin - spin interaction , @xmath27 , or non - local coulomb interaction , @xmath28 . within the present approach ,
all these contributions are treated on equal footing , and we can regard @xmath29 as generalized hopping  matrix elements .
our main object of interest is the green s function of the hubbard operators .
however , the theory is expressed naturally in terms of a generalized irreducible two - point cumulant , m , and a dressed hopping , @xmath30 . using the language of a standard field theory ( see the appendix ) , we interpret m as a generalized self energy and @xmath30 as the corresponding green s function associated with a set of auxiliary canonical fermionic and bosonic fields , and show that they have a clear diagrammatic interpretation in the locator expansion .
in addition , our formulation contains the average @xmath31 of the hubbard operator , which can be viewed as a generalized magnetization , together with a generalized effective magnetic  field , h , represented , in the language of the standard field theory , by the mean value of the auxiliary field .
formally , the relationship between these quantities and the central object of the theory ,  g , is given by the equations : @xmath32 and @xmath33 where e represents the bare coupling constants .
next , we introduce the functional @xmath34 = -\mbox{tr}~\ln[1-me ] - \mbox{tr}[m{\cal e } ] -\frac{1}{2}qeq + \psi[{\cal e } , q ] ,                 \nonumber\end{aligned}\ ] ] where @xmath35 represents a generalized baym - kadanoff - type functional that can be obtained as a sum of all vacuum - to - vacuum skeleton diagrams . using this functional perspective
, we can view @xmath30 and m , as well as q and h , as pairs of conjugate quantities and we have @xmath36 and @xmath37 , where the effective magnetic  field can be expressed as @xmath38 .
the single site dynamical mean field theory is a local approximation for  @xmath39 , which can be obtained as a stationarity condition for the functional @xmath40 .
this condition translates into a self - consistent impurity problem defined by the statistical operator @xmath41 where @xmath42 represents an external field and @xmath43 the hybridization .
these quantities are subjected to the self - consistency conditions : @xmath44^{-1 } =   [ m_{loc}^{-1}(i\omega_n ) - \delta(i\omega_n)]^{-1 } ,   \nonumber\end{aligned}\ ] ] and @xmath45    the strategy that we use in presenting our g - dmft scheme consists in constructing three formally distinct formulations . on the one hand , this allows us to establish the equivalence of these approaches to the problem of correlated electrons and to extract a generalized unitary picture . on the other hand
, these various angles clearly reveal the aspects that are relevant for the possible extensions of the theory .
the natural starting point for our construction is given by a generalized expansion around the atomic limit .
this expansion , in terms of the coupling constants of the non - local contributions to the hamiltonian , is derived in sec .
the technique is characterized by a close formal analogy with the canonical  perturbation expansion in terms of inter - site hopping .
consequently , we will focus on the new aspects generated by the description of the problem in terms of hubbard operators , the main goal being to obtain the relation between the green s function for the hubbard operators and the irreducible two - point cumulant . using this relation ,
we derive in sec .
iii the generalized dmft equations .
a simple illustration of the implementation of our scheme is given for the hubbard model . in sec .
iv we formulate the theory in terms of a functional of the renormalized hopping .
the generalized dmft equations are shown to be the result of a simple local approximation on a generalized baym - kadanoff - type functional . an alternative derivation is described in the appendix . by decoupling the non - local term via a hubbard - stratonovich transformation ,
we formulate the theory in terms of a set of canonical fermionic and bosonic fields . in this language ,
the standard treatment based on wick s theorem applies . a possibility of expanding around the dmft solution
is also presented .
expansions around the atomic limit have been powerful tools for studying both models with localized spins@xcite , and models of interacting itinerant fermions . following the pioneering work of hubbard@xcite
, metzner@xcite developed a renormalized series expansion for the single - band hubbard model .
generalizing the `` linked - cluster expansion '' ideas@xcite , this approach involves only connected diagrams and unrestricted lattice sums , and enables one to construct self - consistent approximations . by analyzing this renormalized expansion , it was shownreview that , in the limit of infinite spatial dimensions , the dynamical mean - field theory equations for the hubbard model are recovered .
an extension of this approach to the case of correlated hopping was developed by shvaika@xcite .
the purpose of this section is to obtain a generalization of the dmft equations that describe the physics of a generic lattice model for interacting fermions , starting from a renormalized strong - coupling expansion around the atomic limit@xcite . within this scheme ,
all the non - local contributions to the hamiltonian , for example hopping terms , spin - spin interaction terms , or non - local coulomb interaction contributions , are treated on equal footing .
let us consider a system described by the hamiltonian ( [ ham1 ] ) .
the first step in our derivation is to write an expansion for the grand - canonical potential and the greens functions for the @xmath46-operators in terms of `` hopping '' matrix elements @xmath47 and bare cumulants . starting with the statistical operator @xmath48 we can write the grand - canonical potential as @xmath49 where @xmath50 is the grand - canonical potential in the atomic limit ,
@xmath51 is the inverse temperature , and the unperturbed ensemble average is given by @xmath52 expanding the exponential in eq .
( [ sigm ] ) we obtain the nth order contribution to @xmath53 : @xmath54 where @xmath55 represents the imaginary - time ordering operator , and @xmath56 are the hubbard operators in the interaction representation , @xmath57 using the multiplication properties ( [ oper1 ] ) of the hubbard operators , it is straightforward to obtain the explicit imaginary time dependence , @xmath58 because @xmath59 is a sum of local operators , the ensemble average in ( expan ) factorizes into independent local averages that can be evaluated using the algebra of the x - operators .
however , the explicit calculation of @xmath60 is cumbersome , as the summation over the site variables @xmath61 is restricted . to overcome this problem , following metzner@xcite
, we introduce the bare cumulant defined by @xmath62 where the generating functional is @xmath63 = \ln\bigg\langle \hat{t}\exp\left\ { - \sum_{\alpha \beta } \int_0^{\beta } d\tau ~\eta_{\alpha \beta}(\tau ) x^{\alpha \beta}(\tau)\right\}\bigg\rangle_0 .
\label{wfunc}\ ] ] the fields @xmath64 are either complex numbers or grassmann variables , depending on the nature ( bosonic or fermionic ) of the corresponding operators @xmath65 . in the eqs .
( [ cumul0 ] ) and ( [ wfunc ] ) the site index , @xmath9 , for the @xmath66 field and the hubbard operators , was omitted for simplicity .
using the bare cumulants allows one to express the grand - canonical potential as an expansion containing unrestricted sums . each term will be a product of `` hopping '' matrix elements , @xmath47 , and local cumulants , and
can be represented diagrammatically .
further , using linked - cluster type arguments@xcite one can show that only the connected diagrams contribute .
the basic diagrammatic rules are similar to those given by metzner@xcite for the hubbard model .
however , the internal lines corresponding to `` hopping '' matrix elements will carry two extra pairs of indices @xmath67 and @xmath68 ( one for each end of the line ) representing the single - site states .
also , the nth order cumulant , @xmath69 , represented by a n - valent point vertex , will be attached to n lines , with the corresponding labels @xmath70 . if we split the states associated with a cumulant in two subsets , @xmath71 and @xmath72 , the multiplication rule ( [ oper1 ] )
requires that each single - particle state ( @xmath73 , @xmath74 , or @xmath75 ) be found the same number of times in each of the subsets .
also , in contrast with the simple case of an expansion in the hopping matrix , @xmath26 , in general , cumulants of order n with n odd may be non - zero .
this is the case if @xmath76 contains terms with bosonic x - operators having a non - zero average value .
such terms may occur , for example , if a non - local coulomb interaction @xmath77 is considered . to summarize these observations ,
we present in fig .
[ fig1 ] the diagrams yielding the leading contributions to @xmath60 .
contains terms with bosonic x - operators having non - zero average values ( q - operators).,scaledwidth=50.0% ]    our next task is to write an expansion for the green s functions of the hubbard operators . we define @xmath78\ } , \label{green1}\ ] ] where the statistical operator @xmath79 is given by eq .
( [ sigm ] ) , and the hubbard operators are in the heisenberg representation , @xmath80 the perturbation expansion for the green s functions is similar to the expansion for the grand - canonical potential .
the consequence of the presence of two extra x - operators in ( [ green1 ] ) is that the corresponding diagrams will be rooted , i.e. one ( for @xmath81 ) or two of the vertices will have fixed site indices and some fixed single - site state indices .
the general rules are analogous to those for the hubbard model@xcite , with the difference that extra indices ( associated with the single - site states ) are carried by the `` hopping '' matrix elements and the bare cumulants .
a special case occurs when bosonic x - operators with non - zero average values , @xmath82 , which we will call q - operators , are present in @xmath76 . in this case , as mentioned before , cumulants of order n with n odd may be non - zero , leading to disconnected rooted diagrams in the expansion of @xmath83 . by simply regrouping the contributions from disconnected diagrams we can write : @xmath84 where @xmath85 represents the average of the operator @xmath86 with respect to the full hamiltonian , and @xmath87 has a diagrammatic representation containing only connected diagrams . in the rest of this work we will always use the `` connected '' green s functions defined by ( [ green2 ] ) . obviously , if @xmath76 does not contain q - operators , @xmath88 and @xmath89 coincide , as the averages in ( [ green2 ] ) vanish .
to summarize , we present in fig .
[ fig2 ] some leading contributions to @xmath88 .     only diagrams containing cumulants with an even number of legs are shown.,scaledwidth=50.0% ]    at this point , a brief discussion about the fourier transformation of the green s functions and of the cumulants is required .
due to the fact that the hubbard operators have either bosonic - like or fermionic - like characters , the fourier transformed green s functions , @xmath90 , will have odd matsubara frequency , if both x - operators are fermionic , or even matsubara frequency otherwise . on the other hand
, the cumulants may depend on both even and odd matsubara frequencies .
recalling the definition , eqs .
( [ cumul0]-wfunc ) , a pair of indices @xmath91 associated with a fermionic - like hubbard operator will induce a dependence on an odd frequency @xmath92 , while an even frequency will be associated with a pair with bosonic character . a line connecting two vertices , and yielding the factor @xmath93 , will carry the frequency corresponding to the nature of the two pairs @xmath94 and @xmath68 .
this nature is always the same as a result of the fact that the hamiltonian conserves the number of particles .
in the rest of this work we will not use different notations for odd and even frequencies , but we should keep in mind this aspect of the problem .    .
the simplified notation @xmath95 , @xmath96 was used .
the diagrams containing one - particle reducible decorations of bare cumulants can be summed and replaced by diagrams with decorations equal to @xmath97 .
the average @xmath98 is represented diagrammatically by a square.,scaledwidth=70.0% ]    the expansion that we derived is rather complicated and it may be of little help if used directly . however , we know that in order to obtain self - consistent approximations and , in particular , dmft - type equations , a renormalized expansion is necessary .
the key in obtaining a renormalized expansion is to introduce the notion of irreducibility with respect to one line ( representing @xmath93 ) .
we define the irreducible two - point cumulant , @xmath99 , as the sum of all connected diagrams having two external legs ( labeled by @xmath100 and @xmath101 , respectively ) , with the property that the two legs can not be separated by cutting one line .
this definition differs slightly from the standard definition@xcite , which requires that the diagram be not disconnected by cutting a single line .
the difference comes from the existence of cumulants with odd number of legs .
if @xmath76 does not contain q - operators , the standard diagrammatics , along with the usual definition of irreducibility are valid .
some typical , low order contributions to the irreducible cumulant @xmath99 are shown in fig.[fig3 ] .
notice that we can classify the diagrams into classes , each class containing one `` basic '' diagram that is one - particle irreducible with respect to all `` hopping '' lines plus diagrams derived from it by decorating the bare cumulants with reducible contributions .
it is useful to perform partial summations of the diagrams in each class .
the resulting diagrams will have all the reducible decorations replaced by factors equal to @xmath102 .
we can now write a renormalized expansion for the green s functions in terms of irreducible cumulants and `` hopping '' matrix elements . the diagrammatic representation of this expansion
is shown in fig.[fig3 ] .
analytically , the corresponding equation will be : @xmath103 where , to reduce the number of indices , we introduced a simplified notation for a pair of single - site states , @xmath104 .
it is also convenient to regard a quantity @xmath105 as an element of a matrix @xmath106 . using this matrix notation , and fourier transforming the spatial dependence
, one can write the relation between the irreducible cumulant and the x - operator green s function , @xmath107^{-1 } ,   \label{lark2}\ ] ] where @xmath108 is the fourier transform of the coupling constant matrix . what exactly are the green s function included in the matrix @xmath109 ?
obviously , those that are coupled by the generalized `` hopping '' @xmath110 .
specifically , if @xmath111 and @xmath112 are the sets of single - site pairs that label the coupling constants in @xmath76 , each pair appearing only once in the corresponding set , the matrices will have dimension p and the matrix elements will be labeled @xmath113 , @xmath114 , and @xmath115 , respectively .
a particular example of this construction will be presented in the next section for the hubbard model .
equations ( [ lark1 ] ) and ( [ lark2 ] ) represent the main result of this section .
they are valid in arbitrary dimensions and for any hamiltonian that has the form given by eqs .
( [ ham1]-[ham3 ] ) when written in terms of hubbard operators .
in the standard dynamical mean - field theory@xcite , the single - particle self energy for the electron green s function is local in the limit of infinite spatial dimensions .
however , for a general lattice hamiltonian containing fermionic x - operators , this property is violated , as shown by schiller@xcite in the case of correlated hopping .
nevertheless , one can overcome this problem by adopting a strong - coupling perspective@xcite and regarding the two - point irreducible cumulant , rather than the self - energy , as the the quantity in which the theory should be formulated .
if @xmath76 contains only fermionic x - operators , and assuming that the coupling constants in @xmath76 are properly scaled , to keep the energy ( per site ) finite , one can show@xcite that the irreducible cumulants become local , @xmath116 consequently , the natural quantity in which such a theory should be formulated is the irreducible cumulant , rather than the self - energy .
this becomes immediately transparent using eqs .
( [ lark2 ] ) and ( [ mm0 ] ) , and the fact that the electron green s function is given by the sum @xmath117 .
\label{elgreen}\ ] ] although @xmath118 is k - independent , @xmath119 will have , in general , a complicated dependence on the wave - number , and the corresponding self - energy will be non - local , breaking the standard dmft scheme . formally , a self - energy matrix can be introduced by defining @xmath120 with @xmath121 being the identity matrix . however , in general , this quantity does not have the desired analytic properties@xcite and may diverge as @xmath122 .
if we go further and generalize the model to include both fermionic and bosonic degrees of freedom , the picture presented above fails , as one can not scale the coupling constants in such a way that one has both finite energy ( per site ) and local irreducible cumulants in the limit of infinite dimensions@xcite .
in contrast , this work assumes a different perspective : rather than viewing dmft as an approximation that becomes exact in the limit of infinite coordination , we see it as the first order approximation in a hierarchy of ( increasing size ) cluster approximations . in this picture , the exact limit is obtained for an infinite - large cluster .
consequently , our g - dmft scheme will be defined by the condition that the two - point irreducible cumulant be local . unlike the case of extended dmft
, we have no small parameter ( other than the inverse of the cluster size ) justifying this approach . obviously , if the model contains only fermionic degrees of freedom , the two perspectives described above are equivalent .    based on the locality of @xmath123
, we can approach the central step in a dmft - type scheme , namely the self - consistent mapping of the lattice problem onto an effective impurity problem .
let us consider the single - impurity problem with the statistical operator @xmath124\right\ } ,   \label{impur}\ ] ] where @xmath125 is an external field that couples to the q - operators , @xmath126 is the set of indices corresponding to the q - operators , @xmath127 is the hybridization matrix and @xmath128 contains one term , for example @xmath129 , from the local hamiltonian ( [ ham2 ] ) . using the expansion procedures described in the previous section , with the hybridization replacing the `` hopping '' matrix @xmath130
, we obtain for the impurity green s function a relation similar to ( [ lark2 ] ) , @xmath131^{-1},\ ] ] where all the quantities are local , and @xmath132 is the irreducible cumulant with respect to the hybridization . for practical calculations , as well as to gain some physical insight , it is useful to express the impurity problem in a hamiltonian form . taking into account that in our matrix formulation the indices
@xmath67 corresponding to fermionic and bosonic hubbard operators do not mix , all the matrices , including @xmath133 , are block - diagonal with one block corresponding to the fermionic indices and one to the bosonic sector . consequently , we will introduce two types of auxiliary degrees of freedom , one associated with a fermionic bath , described by operators @xmath134 , and one associated with a bosonic bath , described by @xmath135
. the impurity hamiltonian will be @xmath136 + \sum_{k}\sum_{(\alpha \beta)\in\mathcal{b } } [ v_b^{\alpha \beta}(k ) ~b_{k\sigma_{\alpha \beta}}^{\dagger } x^{\alpha \beta } + h.c .
\label{himp}\end{aligned}\ ] ] the first term in the right hand side describes an isolated impurity , the second term represents the coupling to the external field and the summation is restricted to labels corresponding to bosonic operators having a non - zero mean value , the q - operators ( @xmath126 ) , the third and forth terms describe the fermionic and bosonic baths , respectively , and the last two terms represent the coupling of the impurity with the two baths . in the last two terms ,
the summations over the x - operator labels @xmath137 are restricted to the fermionic ( @xmath138 ) and bosonic sectors ( @xmath139 ) , respectively .
the expression of the hybridization in terms of coupling constants @xmath140 can be obtained by integrating out the auxiliary degrees of freedom .
we have @xmath141 where @xmath142 if the element is in the fermionic sector of the @xmath43 matrix and @xmath143 for the bosonic case .
returning to our main problem , we obtain the self - consistency conditions by requiring that the lattice irreducible cumulant be equal to the effective impurity problem irreducible cumulant , @xmath144 . by comparing the diagrammatic expansion of @xmath132 in terms of bare cumulants ( dressed with external field lines ) and hybridization functions with the expansion of @xmath145 , we observe that the equality is satisfied if we identify a line representing @xmath146 with the sum of all the paths starting with @xmath147 and ending with @xmath148 , as well as the sum betweeen an external field line @xmath149 and the one - particle reducible decoration @xmath150 with the lattice one - particle reducible decoration @xmath151 . explicitly , the expansion for the impurity irreducible cumulant is identical with the expansion of the lattice irreducible cumulant provided @xmath152 where @xmath153 represents the matrix green s function for a lattice with the site `` 0 '' removed , together with a condition that determines self - consistently the external field @xmath154 , @xmath155 notice that equation ( [ self1 ] ) represents the generalization of the well known dmft self - consistency condition@xcite .
in addition , if the non - local hamiltonian @xmath76 contains hubbard operators that acquire a non - zero mean value , our lattice model is mapped into an impurity model with an external field determined self - consistently by equation ( [ selfh ] ) .
our final step consists in finding a self - consistent equation that determines the hybridization @xmath156 in terms of the lattice green s function . comparing eq .
( [ lark1 ] ) for @xmath157 ( after making @xmath158 ) with the equivalent expression for @xmath159 , @xmath160 we obtain for @xmath161 : @xmath162^{-1}.\end{aligned}\ ] ] using equation ( [ self1 ] ) for @xmath163 and the fourier transform for the local green s function , we obtain the self - consistent condition that determines the hybridization , @xmath164^{-1}=\left[\mathbf{m}_0^{-1}(i\omega_n ) - \mathbf{\delta}(i\omega_n)\right]^{-1}.   \label{self0}\ ] ]    formally , eq .
( [ self0 ] ) closely resembles the standard form of the dmft self - consistency condition@xcite .
we should , however , keep in mind the fact that this is a matrix equation involving green s functions for the hubbard operators . to solve a problem described by the hamiltonian ( ham1-[ham3 ] ) within this scheme , implies to solve the impurity problem , ( [ impur ] ) , find an irreducible cumulant matrix @xmath118 , then use the self - consistency condition ( [ self0 ] ) together with equation ( [ selfh ] ) to determine the new hybridization , @xmath165 , and the new field @xmath166 .
the process is repeated until convergence is reached .      in order to get a feeling of how this approach works , as well as to check its validity
, we can do a simple exercise by applying the procedure to the hubbard model , for which we should be able to recover the known dmft results . in this particular case ,
the hamiltonian will be @xmath167 + \sum_{i , j}\sum_{\sigma}~t_{ij}[x_i^{\sigma 0}x_j^{0 \sigma } + x_i^{2 \bar{\sigma}}x_j^{\bar{\sigma } 2 } + \sigma x_i^{\sigma 0}x_j^{\bar{\sigma } 2 } + \sigma x_i^{2 \bar{\sigma}}x_j^{0 \sigma } ] , \label{hbham}\ ] ] where @xmath168 represents the chemical potential , @xmath24 the on - site coulomb interaction and @xmath26 the hopping matrix elements , which are equal to @xmath169 if i and j are nearest neighbors and zero otherwise .
we can introduce the hopping matrix , @xmath170 where @xmath171 is the fourier transform of @xmath26 .
next , we introduce the green s function matrix and the irreducible cumulant matrix : @xmath172 our goal is to show that the matrix equation ( [ self0 ] ) reduces in this particular case to the well known dmft self - consistency condition , @xmath173^{-1 } + \sigma(i\omega_n),\ ] ] where @xmath174 represents the self - energy of the single - particle electron green s function and @xmath175 the hybridization of the effective single - impurity problem .
the electron single - particle green s function , @xmath119 , is given by a sum of all the elements of @xmath109 , as shown in eq .
( [ elgreen ] ) .
similarly , we define the quantity @xmath176.\ ] ] it is straightforward to show that the matrix green s function has the structure : @xmath177^{-1 } = \frac{1}{1-\epsilon_k \mathcal{m}(i\omega_n)}\mathbf{m}(i\omega_n ) + \delta \frac{\epsilon_k}{1-\epsilon_k \mathcal{m}(i\omega_n ) } \left(\begin{array}{cc } -1 & \sigma \\ \sigma & -1 \\ & \end{array}\right ) ,   \label{struct}\ ] ] where @xmath178 . using eq .
( [ struct ] ) and the expression ( [ elgreen ] ) for the electron green s function , together with the definition of @xmath179 , we have @xmath180 which shows that @xmath181 is the single particle irreducible cumulant , @xmath182 , where @xmath174 is the single particle self - energy . on the other hand , making the summation over @xmath183 in eq .
( [ struct ] ) and introducing the result in the self - consistency relation ( [ self0 ] ) gives : @xmath184 with @xmath185 taking into account the structure of @xmath186 given by eq .
( eldelta ) , the impurity problem described by the statistical operator ( impur ) reduces to an anderson single impurity problem with a hybridization @xmath187 .
finally , we have @xmath188^{-1},\ ] ] which represents the standard dmft self - consistency equation .
functionals of green s functions provide an elegant and powerful tool for formulating dynamical mean - field theories . within this framework ,
the dmft equations are obtained as a direct result of a simple approximation on the baym - kadanoff functional@xcite .
the extended dynamical mean - field equations have also been rigorously derived using a functional techniqueedmft5 .
moreover , the functional perspective seems to offer the natural framework for cluster generalizations of the dynamical mean - field theory that include effects of short - range correlations .
the purpose of this section is to introduce a strong - coupling generalization of the functional dmft formulation@xcite .
the construction of a generalized baym - kadanoff - type functional will allow us to present an alternative derivation of the generalized dmft equations , and , more importantly , will set the basis for a future cluster extension of the g - dmft approach .
the main idea of the functional interpretation of dmft is to formulate a local approximation for the baym - kadanoff functional @xmath189 $ ] , which is the sum of all vacuum - to - vacuum skeleton ( two - particle irreducible ) diagrams constructed with the full propagator @xmath88 and the interaction vertices .
this functional has the property that @xmath190 where @xmath191 is the self - energy .
as the coordination number of the lattice goes to infinity , the self - energy becomes local and @xmath192 will depend only on the local green s functions , @xmath193 : @xmath194,\ ] ] where @xmath195 $ ] is a functional of the local green s function at site @xmath196 only .
this scheme can not be transfered directly to approach our strong - coupling problem . as we mentioned above , within our strong - coupling perspective , the natural quantity to describe the system is the two - point irreducible cumulant , @xmath197 , rather than the self - energy .
our task is to identify a conjugate quantity ( which can not be the green s , function but rather something with dimensions of energy ) , and to construct a functional of this quantity that will determine @xmath198 by a relation similar to ( phideriv ) .    following refs . and , we start by introducing two time dependent external sources , one , @xmath199 , coupled to the x - operator green s functions , and the other , @xmath200 , coupled to the average @xmath201 of the bosonic hubbard operators that occur in the non - local hamiltonian @xmath76 .
the corresponding generating functional will be @xmath202 = -\frac{1}{\beta } \ln tr \hat{\rho}[j,\lambda],\ ] ] where the statistical operator @xmath203 is given by @xmath204 = e^{-\beta h_0 } \hat{t } \exp\left\{-\int_0^{\beta}d\tau\left[\int_0^{\beta}d\tau^{\prime}\left(h_1(\tau , \tau^{\prime } ) + \sum_{i , j}\sum_{\mu,\nu}j_{ji}^{\nu\mu}(\tau-\tau^{\prime } ) x_i^{\mu}(\tau)x_j^{\nu}(\tau^{\prime})\right ) + \sum_{i , \mu } \lambda_i^{\mu}(\tau ) x_i^{\mu}(\tau ) \right ] \right\}.\ ] ] the field @xmath199 is equivalent with introducing frequency dependent coupling constants @xmath205 .
the green s functions can be expressed as functional derivatives of the generating functional , @xmath206 $ ] , with respect to the frequency dependent field , @xmath207}{\delta e_{ji}^{\nu\mu}(i\omega_n)}\bigg|_{j=0 , \lambda=0 } ,   \label{fderivg}\ ] ] where the average value @xmath208 is non - zero only if @xmath209 is a q - operator .
the generating functional @xmath206 $ ] reduces to the grand - canonical potential , @xmath60 , in the absence of the external sources .
the next step is to integrate eq .
( [ fderivg ] ) .
as we are interested in a strong - coupling equivalent of the baym - kadanoff functional , we will express the green s function in terms of the two - point irreducible cumulant , rather than the self - energy@xcite . using the notation @xmath210 for a tensor with elements
@xmath211 , we can formally write equation ( [ lark1 ] ) as @xmath212^{-1}\hat{\mathbf{m } } ,   \label{hatg}\ ] ] where @xmath213 is the identity tensor , @xmath214 . by integrating ( [ fderivg ] ) with @xmath215 given by ( [ hatg ] )
we obtain @xmath202 = \omega_0 + \frac{1}{\beta } \sum_{\omega_n}tr~\ln[\hat{\mathbf{i } } - \hat{\mathbf{m}}\hat{\mathbf{e } } ] + \frac{1}{\beta } \sum_{\omega_n}tr[\hat{\mathbf{m}}\hat{\mathbf{{\mathcal{e } } } } ] + \xi[j,\lambda ] ,   \label{gama1}\ ] ] with @xmath216 the grand - canonical potential in the atomic limit .
the quantity @xmath217 is the renormalized coupling constant ( or renormalized generalized hopping ) and can be represented diagrammatically as a sum of chains of `` hopping '' lines and irreducible cumulants .
analytically we have @xmath218_{ij}^{\mu\nu}.   \label{cale}\ ] ] the last term in eq .
( [ gama1 ] ) is a functional with the property @xmath219}{\delta e_{ji}^{\nu\mu}(i\omega_n ) } = tr\left[\hat{\mathbf{m}}\frac{\delta \hat{\mathbf{{\mathcal{e}}}}}{\delta e_{ji}^{\nu\mu}(i\omega_n)}\right ] -\beta\langle x_i^{\mu } \rangle\langle x_j^{\nu } \rangle .
\label{xideriv1}\ ] ] also , because the functional derivative of @xmath206 $ ] with respect to @xmath220 gives the average value of @xmath209 , we will have @xmath221 + \langle x_i^{\mu } \rangle .
\label{xideriv2}\ ] ]      for the case with no q - operators in the non - local hamiltonian our construction is complete . in this particular case
the only relevant dependence in the functionals is on the renormalized coupling @xmath234 .
the functional @xmath235 $ ] is formally analogous to the functional @xmath236 $ ] introduced by chitra and kotliar@xcite from a weak - coupling perspective . in our picture the strong coupling correspondent of the self - energy is , as noted before , the irreducible cumulant , while the renormalized and bare coupling constants correspond to the full and non - interacting green s function , respectively .
also , the strong - coupling analog of the luttinger - ward functional @xmath189 $ ] is @xmath237 $ ] .
similarly to its weak - coupling counterpart , @xmath237 $ ] can be expressed as a sum of all vacuum - to - vacuum skeleton diagrams containing bare cumulants and renormalized `` hopping '' lines .
all the skeleton diagrams are two - particle irreducible , i.e. they can not be separated into disconnected parts by cutting two @xmath234 lines .    .
the contribution of the external field @xmath238 is represented only in diagram ( i ) by dotted lines .
the numbers in parentheses represent the symmetry factors ( if different from 1 ) of the corresponding diagrams ( the number of ways in which a graph can be made isomorphic with itself ) .
the renormalized coupling @xmath234 is represented by double lines.,scaledwidth=50.0% ]    in general , however , the diagrammatic expression of @xmath230 is more complicated and involves , in addition to the bare cumulants , both renormalized coupling constants , @xmath239 , and bare coupling constants , @xmath110 .
also , lines corresponding to the external field @xmath220 have to be added .
the generalized skeleton diagrams have to satisfy two conditions : i ) each diagram is two - particle irreducible with respect to the renormalized `` hopping '' ; ii ) by cutting any bare `` hopping '' line the diagram is divided into two disconnected parts .
these conditions determine a tree - like structure for the generalized skeleton diagrams with skeleton blocks coupled by bare @xmath227 lines .
this structure is illustrated in fig .
[ fig4 ] for a few low - order contributions to @xmath230 . due to the tree - like structure of the diagrams ,
it is rather difficult to find a local approximation for the functional @xmath230 . in order to overcome this problem
, we proceed to the last step of our construction and define the legendre transform of @xmath240 with respect to the average of bosonic hubbard operators , @xmath241 , @xmath242 = \omega^*[{\mathcal{e } } , \lambda ; e ] -    \sum_i\sum_{\mu\in \mathcal{b}_q } \lambda_i^{\mu}q_i^{\mu } ,   \label{gamaeq1}\ ] ] where the bare coupling @xmath110 is treated as a parameter and is not written explicitly as an argument of the functional @xmath243 . using the condition @xmath244 to eliminate @xmath220
, we obtain @xmath242 = \omega_0 + \frac{1}{\beta } \sum_{\omega_n}tr~\ln[{\hat{\mathbf{i } } - \hat{\mathbf{m}}[{\mathcal{e } } , q]\hat{\mathbf{e } } } ] +   \frac{1}{\beta } \sum_{\omega_n}tr[\hat{\mathbf{m}}[{\mathcal{e } } , q]\hat{\mathbf{{\mathcal{e } } } } ] - \frac{1}{2}\sum_{i , j}\sum_{\mu,\nu\in \mathcal{b}_q } q_i^{\mu }   e_{ij}^{\mu\nu } q_j^{\nu } + \psi[{\mathcal{e } } , q ] .
\label{gama2}\ ] ] the properties with respect to @xmath234 of the new functionals are similar to those of @xmath229 and @xmath230 , namely the stationarity condition @xmath245 yields equation ( [ cale ] ) , and we have @xmath246}{\delta { \mathcal{e}}_{ji}^{\nu\mu}(i\omega_n ) } = m_{ij}^{\mu\nu}(i\omega_n ) .
\label{fderivm1}\ ] ] in addition , stationarity with respect to @xmath247 , @xmath248 , yields @xmath249}{\delta q_i^{\mu } } = \sum_{j , ~\beta\in \mathcal{b}_q } e_{ij}^{\mu\nu } ~q_j^{\nu}.   \label{fderivq}\ ] ] the diagrammatic representation of the new luttinger - ward - type functional @xmath35 contains skeleton diagrams that are two - particle irreducible with respect to the @xmath234 lines and have cumulants dressed with @xmath250 contributions . the structure of the diagrams
is illustrated in fig [ fig5 ] . using the diagrammatic representation of the functionals @xmath35 and @xmath230 , as well as that of @xmath247 , and taking into account the symmetry factors associated with each diagram , it is straightforward to check that @xmath251 = \xi^*[{\mathcal{e}},\lambda=0 ; e ] + \frac{1}{2}\sum_{i , j}\sum_{\mu,\nu\in \mathcal{b}_q } q_i^{\mu } e_{ij}^{\mu\nu } q_j^{\nu},\ ] ] where @xmath252 .
the average of a diagonal operator , @xmath253 , is represented by a square .
the symmetry factors are written in parentheses . if @xmath254 , the remaining basic skeleton diagrams ( a , e , g , ... ) are identical with the remaining diagrams in fig .
[ fig4 ] ( a , b , c , ... ) , and @xmath255 . in general
, we have the basic diagrams plus diagrams obtained by decorating the bare cumulants with lines and squares representing @xmath256.,scaledwidth=40.0% ]    we are prepared now to introduce the generalized dmft from a functional perspective . by analogy with the standard approach@xcite
, we construct an approximate theory by making the ansatz that the luttinger - ward functional @xmath35 depends only on local renormalized `` hopping '' matrices , i.e. , @xmath251 = \sum_i \psi[\mathbf{{\mathcal{e}}}_{ii},q],\ ] ] where @xmath257 is a functional that depends only on the local renormalized coupling constants at site i. in diagrammatic terms , it means that all the bare cumulants contained in a skeleton diagram contributing to @xmath257 are defined on the same site .
a direct consequence of this ansatz is that the two - point irreducible cumulant of the theory is local , @xmath258 with @xmath259 being the kronecker @xmath260 symbol .
the stationarity condition for @xmath40 yields @xmath261_{ii } = \frac{1}{n}\sum_k \mathbf{e}(k)(\mathbf{i } - \mathbf{m}\mathbf{e}(k))^{-1 } ,   \label{caleloc}\ ] ] where @xmath198 is given by ( [ mdr ] ) and , for the second equality , we assumed that the system is homogeneous .
relation ( [ caleloc ] ) represents our generalized dmft equation and it should be viewed as a functional equation for the local renormalized coupling matrix .
to cast it in a more familiar form , let us observe that using the expressions for the green s function , eq .
( [ hatg ] ) , and for the renormalized coupling , eq .
( [ cale ] ) , we have @xmath262 taking into account the local character of the irreducible cumulant , we obtain , for a homogeneous system , @xmath263^{-1 } , \label{caleloc1}\ ] ] where @xmath264 has to be understood as a functional of @xmath234 .    to be able to solve eq .
( [ caleloc ] ) or eq .
( [ caleloc1 ] ) we need to know the functional dependence of @xmath35 , at least approximately .
a direct approach is rather hard in practice , therefore it is more convenient to generate the functionals from a purely local theory .
a simple inspection of the diagrammatic expression for @xmath35 ( see figure [ fig5 ] ) , shows that , in the case when all the bare cumulants are defined on the same site @xmath129 , the same diagrams will correspond to an impurity problem defined by the statistical operator given by eq .
( [ impur ] ) .
in fact , taking into account the previous construction of the strong - coupling functionals , we have @xmath265 = \frac{\omega_0}{n }   & + & \frac{1}{\beta } \sum_{\omega_n}tr~\ln[{\mathbf{i } - \mathbf{m}_{imp}[{\mathcal{d } } , q]\mathbf{\delta } } ] + \frac{1}{\beta } \sum_{\omega_n}tr[\mathbf{m}_{imp}[{\mathcal{d } } , q]\mathbf{{\mathcal{d } } } ]   \\ \nonumber   & - & \frac{1}{2}\sum_{\mu , \nu\in \mathcal{b}_q } q^{\mu } \delta^{\mu\nu}(0)q^{\nu }   -\sum_{\mu \in \mathcal{b}_q } h_{imp}^{\mu}q^{\mu }   + \psi_{imp}[{\mathcal{d } } , q ] ,   \label{gimp } \end{aligned}\ ] ] where @xmath266^{-1}$ ] is the renormalized hybridization , @xmath267 and @xmath268 is the hybridization at zero frequency .
in addition , we have to keep in mind that the bare impurity cumulants are dressed by external field lines , @xmath269 . the diagrammatic expressions for @xmath270 $ ] and @xmath271 $ ] are identical and we have @xmath272 provided @xmath273 equations ( scsce ) represent the self - consistency condition for our strong - coupling generalized dynamical mean - field theory .
the first of these equations can be rewritten in the more familiar form , eq .
( [ self0 ] ) , by simply using the expressions of the renormalized quantities @xmath234 and @xmath274 .
a direct consequence of the first self - consistency condition is the equality of the irreducible cumulants , @xmath275 .
finally , let us mention that this correspondence with the impurity problem is also useful in establishing a relation between the grand - canonical potential of the lattice and the grand - canonical potential of the impurity . at stationarity
both @xmath276 and @xmath40 reduce to the corresponding grand - canonical potential , while @xmath272 .
therefore we have @xmath277 - \frac{1}{\beta}\sum_{\omega_n}\left[tr \ln[\mathbf{i } - \mathbf{m}\mathbf{\delta } ]   - \frac{1}{n } \sum_k tr \ln [ \mathbf{i } - \mathbf{m}\mathbf{e}(k ) ] \right].\ ] ]    this work was supported by the nsf under grant dmr-0096462 and by the rutgers university center for materials theory .
we would like to thank c. henley and a. schiller for useful discussions and s. pankov and a. shvaika for helpful critical observations .
the purpose of this appendix is twofold : 1 ) to present an alternative formulation of our strong - coupling theory that is closely related to the more familiar perturbation expansions based on wick s theorem and translate the basic quantities involved in this approach in the corresponding language .
2 ) to propose a scheme that allows one to expand around the dmft solution .
the central idea is to decouple the nonlocal term using a generalized hubbard - stratonovich transformation and to trade in the x - operators for a set of canonical fermionic and bosonic fields . in terms of these new auxiliary fields ,
the standard techniques of the perturbation theory can be applied .
this formulation represents a generalization of the strong - coupling expansion used by pairault , snchal , and tremblay@xcite for the study of the hubbard model .    starting with the partition function @xmath278\right\ } ,   \label{partfunc}\ ] ]
we perform a hubbard - stratonovich transformation and express the non - local part as a gaussian integral over the auxiliary field @xmath279 , @xmath280\right\ } ,   \label{zdezeta}\ ] ] with @xmath281 = - \ln\bigg\langle \hat{t } \exp\left\{-\int_0^{\beta } d\tau \sum_{i , \mu } \left[\zeta_i^{\mu}(\tau)x_i^{\mu}(\tau ) + x_i^{\mu}(\tau)\zeta_i^{\mu}(\tau)\right]\right\}\bigg\rangle_0 . \label{sdezeta}\ ] ] the field @xmath279 contains both fermionic ( grassmann variables ) and bosonic ( complex ) components , as determined by the corresponding x - operators . the free propagator appearing in the quadratic term of eq .
( zdezeta ) is @xmath282 .
notice that the auxiliary fields have canonical fermionic or bosonic statistics and all the complications related to the unconventional algebra of the hubbard operators are buried in the average @xmath283 with respect to the unperturbed ensemble . due to the nature of the local hamiltonian @xmath59
, the interaction term @xmath284 $ ] can be written as a sum of purely local contributions , each having an expression similar to eq .
( [ sdezeta ] ) but without the summation over the site index @xmath9 .
expanding in the @xmath285-fields , we have explicitly @xmath281 = \sum_i\sum_{n=0}^{\infty } \frac{1}{n!}\sum_{\mu_1 , ... , \mu_n } c_i^n(\mu_1 , ... , \mu_n;\tau_1 , ... ,
\tau_n)~ \zeta_i^{\mu_n}(\tau_n) ...
\zeta_i^{\mu_1}(\tau_1 ) ,   \label{intvert}\ ] ] where @xmath286 is the n - point bare cumulant defined by eq . ( [ cumul0 ] ) .
we can now develop a standard perturbation theory for the @xmath285-field , having @xmath287 as the free propagator , and vertices given by @xmath286 .
the relation between the green s function of the x - operators and that of the auxiliary fields can be easily obtained if we notice that , due to the form of the interaction term ( [ sdezeta ] ) , we have @xmath288}{\delta \zeta_j^{\nu}(\tau^{\prime})\delta \zeta_i^{\mu}(\tau ) } \exp\{-s_0[\zeta]\},\ ] ] where @xmath289 $ ] represents the quadratic term in eq .
( [ zdezeta ] ) . after performing two integrations by part
, we obtain @xmath290 where we used the tensor notation , and we introduced the full propagator , @xmath291 , for the auxiliary fields , defined as @xmath292 where @xmath293 is the ( imaginary ) time ordering operator and @xmath294 represents the average over the statistical ensemble . in terms of the self - energy for the @xmath285 variables , @xmath295 , one has @xmath296 which shows that @xmath297 corresponds to the irreducible cumulant @xmath198 of the x - operator formulation . also , the full green s function @xmath298 is the renormalized coupling costant @xmath299 of the original theory .
in addition , in the presence of bosonic hubbard operators with non - zero mean value , we have @xmath300 showing that the average of the auxiliary field represents the effective magnetic   field of the original formulation .
this equation completes the correspondence between the the formulation in terms of x - operators and that in terms of auxiliary @xmath285-fields .
it is now straightforward to re - derive our generalized mean - field equations using the standard techniques of a perturbation theory for the auxiliary fields .
for example , the baym - kadanoff functional , @xmath301 $ ] , can be defined in the usual way as the sum of all vacuum - to - vacuum skeleton diagrams .
the corresponding free - energy functional can be written as @xmath302 = \frac{1}{\beta}\sum_{\omega_n } tr\ln\left[\hat{\mathbf{\chi}}_0^{-1 } - \hat{\mathbf{\sigma}}_{\zeta}[\chi , q]\right ] + \frac{1}{\beta}\sum_{\omega_n } tr[\hat{\mathbf{\sigma}}_{\zeta}[\chi , q ] \hat{\mathbf{\chi } } ] -\frac{1}{2}\sum_{i , j}\sum_{\mu,\nu}q_i^{\mu}(\chi_0)_{ij}^{\mu\nu}q_j^{\nu } + \phi[\chi , q ] ,   \label{fren}\ ] ] where the possibility of a non - zero average , @xmath247 , for the x - operators has been explicitly considered .
taking into account the correspondence between the auxiliary field formulation and the original x - operator formulation , we can see that @xmath303 $ ] is the equivalent of the functional @xmath304 $ ] , while the baym - kadanoff functional @xmath305 $ ] corresponds to @xmath306 $ ] , up to an @xmath307-independent term equal to the atomic grand - canonical potential @xmath216 .
we can construct now a strong - coupling expansion to a certain order in @xmath308 , as in the work of pairault , snchal , and tremblay@xcite .
however , any finite order expansion will intrinsically contain the problem of the exponentially large degeneracy of the atomic ground state , leading to the break - down of the solution in the low temperature limit@xcite .
next , we can recover our strong coupling dmft scheme as a local approximation for the baym - kadanoff functional .
this renormalized perturbation approach avoids the problem of the low temperature barrier .
further , we can exploit the idea of introducing auxiliary degrees of freedom and construct a formulation which , in principle , allows one to determine corrections to the dmft approximation .
the basic idea is to introduce the auxiliary fields in such a way as to obtain the dmft result as the lowest order approximation .
starting again from the partition function ( partfunc ) , instead of performing directly a hubbard - stratonovich transformation , we first add and subtract in the exponent a term equal to @xmath309 where @xmath310 is the hybridization function for the impurity problem associated with the dmft approximation .
notice that this is a purely local quantity .
next , we introduce the auxiliary field @xmath311 that decouples the quadratic term containing @xmath312 . the effective theory for
the new variables will have a free term @xmath313 = - \int_0^{\beta } d\tau \int_0^{\beta } d\tau^{\prime}\sum_{<i , j > } \sum_{\mu , \nu } \xi_i^{\mu}(\tau ) ( \chi_0^{-1})_{ij}^{\mu\nu}(\tau-\tau^{\prime } ) \xi_j^{\nu}(\tau^{\prime}),\ ] ] and an interaction term @xmath314 = -\ln\bigg\langle \hat{t } \exp\left\
{ -\int_0^{\beta } d\tau   d\tau^{\prime}~x_i^{\mu}\delta_{ii}^{\mu\nu } x_i^{\nu } + \xi_i^{\mu}x_i^{\mu }   + x_i^{\mu}\xi_i^{\mu}\right\}\bigg\rangle_0 ,   \label{sintdmft}\ ] ] where , for simplicity , we omitted the summations and the imaginary time dependence .
again , the interaction part can be written as a sum of purely local terms . however , the vertices of the new theory are no longer the bare cumulants of the hubbard operators . instead ,
let us notice that each term @xmath315 $ ] of the interaction part represents the generating functional @xmath316 $ ] for the impurity problem defined by eq .
( [ impur ] ) . in the present formulation ,
the external fields @xmath317 are represented by the average values @xmath318 .
a direct consequence of this correspondence is that the second order cumulant generated by @xmath319 is equal to the impurity green s function @xmath320 associated with the generalized dmft solution of our problem .
further , the relation between the propagator for the @xmath311-field and the x - operator green s function can be derived as before and we have @xmath321 with @xmath322 the self - energy for the @xmath311-field .
let us construct now the lowest order approximations for the pertubation expansions formulated in terms of auxiliaryy fields @xmath285 and @xmath311 , respectively . in the first case ,
the simplest diagram contributing to @xmath297 is equal to the atomic green s function .
the corresponding result represents the hubbard - i approximation@xcite . in the second case ,
the zero - order approximation for the self - energy is given by the dmft impurity green s function , @xmath323 .
using the correspondence relation ( [ gdechi ] ) with @xmath324 , as well as the equation @xmath325 , we have @xmath326^{-1}.\ ] ] this shows that , in the formulation using auxiliary @xmath311 fields , the generalized dmft result is recoverd as the lowest order approximation .
obviously , in order to determine the hybridization function @xmath43 , one has to first find the solution of the problem within the dmft scheme .
however , the strong coupling expansion constructed here allows us to determine in a systematic way the corrections to the dmft result by computing higher order contributions to the self - energy @xmath322 .