phonons are important in metallic systems .
they affect the electronic behaviour in diverse ways causing the dominant temperature dependent contribution to the resistivity , an enhancement of the specific heat , and an attractive effective electron - electron interaction , which may induce superconductivity .
one can expect strong electron - phonon interactions in strongly correlated systems , such as heavy fermion compounds@xcite , where the radius of the rare earth ions is very sensitive to the @xmath0-electron occupation .
the coupling of electronic states to the lattice also plays an important role in the anomalous electronic behaviour of the manganites@xcite .
the superconductivity in the fullerides@xcite appears to be induced by electron - phonon interactions within a strongly correlated electron band .
however , the fully quantum mechanical treatment of lattice effects in these compounds has so far received comparatively little attention .
the reason for that lies in the difficulty in handling strong electron - phonon coupling together with strong electron - electron repulsion .
recent advances in non - perturbative techniques , most notably the development of the dynamical mean field theory ( dmft)@xcite , enable one to investigate the interplay of electron and phonon interactions in these systems .
one of the most prominent models in the field of strongly correlated electron systems is the hubbard model@xcite which been used extensively to study the effects of strong local electron - electron interactions . especially within the framework of dmft@xcite the nature of the mott transition
could be clarified@xcite .
the holstein model@xcite has been used to study polaronic effects in the absence of electron - electron repulsion , mainly in the limit of low electron density .
systems with finite electron density have received less attention , but have been studied within the dmft using monte carlo@xcite , exact diagonalization ( ed)@xcite , perturbative approaches@xcite and the numerical renormalization group ( nrg)@xcite .
even less studied has been the interplay of electron - electron interactions and the electron - phonon coupling@xcite .
the most natural model incorporating both these terms is the holstein - hubbard model defined by the hamiltonian @xmath1 here @xmath2 describes the electron
electron interaction within a band of dispersion @xmath3 .
the electron density @xmath4 at site @xmath5 couples linearly to the local displacement operator @xmath6 with an electron - phonon coupling @xmath7 .
the phonons are assumed to be dispersionless ( local einstein phonons ) with energy @xmath8 and @xmath9 is the mass of the vibrating ions .     and a semielliptic band of width @xmath10 as calculated in ref .  .
the dashed line is the locus of @xmath11 .
points on the dot - dashed line fulfil @xmath12.,scaledwidth=46.0% ]    recently , the @xmath13 phase diagram of the half filled holstein - hubbard model has been calculated@xcite using various local approximations , and is shown in fig .
[ phases.eps ] .
the results are for a semi - elliptic band of width @xmath10 and phonon frequency @xmath14 ( @xmath15 ) .
all types of long - range order are excluded .
we can distinguish three different phases : metallic , bipolaronic and mott insulating phase .
the metallic phase is always found to be a fermi liquid .
the latter two phases display a gap in the one - electron spectra and will be referred to as gapped phases .
the parameters of fig .
[ phases.eps ] , namely @xmath16 and @xmath17 are taken throughout this paper , unless mentioned otherwise .
one way of simplifying the discussion of the holstein - hubbard model would be to integrate out the phonons in order to derive an effective model for the electrons .
the effective action is then governed by the @xmath18-dependent potential@xcite @xmath19 this reduces to the static quantity @xmath20 in the limits @xmath21 and @xmath22 .
therefore , on small energy scales , we expect an effective attractive or repulsive interaction depending on the value of @xmath23 . for small phonon energies
@xmath8 , however , there is no reason to expect the static quantity @xmath23 to be sufficient to characterize the behaviour of the holstein - hubbard model .
the properties of this model are expected to depend on both @xmath2 and @xmath7 individually .
this is also apparent from the phase diagram in fig .
[ phases.eps ] .
the dot - dashed line , where the double occupancy takes the value of a free system @xmath12 is quite distinct from the line @xmath24 .
the aim of this study is to clarify the interplay of these different types of interactions .
an interesting question about this model is whether the suppression of local charge fluctuations by a significant hubbard @xmath2 will effectively decouple electrons and phonons on all energy scales .
it is possible , however , that for the behaviour on low energy scales the effective onsite-@xmath2 is strongly renormalized and hence the electron phonon coupling will have a comparatively stronger effect on energy scales @xmath25 ( see ref  ) .    in this paper
, we present results on dynamical electronic and phonon response functions for the holstein - hubbard model as calculated within the dynamical mean field theory@xcite ( dmft ) .
the effective impurity problem is solved using wilson s numerical renormalization group ( nrg ) method@xcite extended to treat the electron
phonon coupling@xcite .
the combination of dmft and nrg has proven to be a reliable and effective method for calculations of both metallic and insulating phases at @xmath13 .
it has been applied to the mott transition in the pure hubbard model@xcite and to the pure holstein model@xcite . in our calculations
we use the discretization parameter @xmath26 and retain up to @xmath27 states .
the value of @xmath26 gives rise to a bandwidth correction factor @xmath28 ( see eq . ( 5.42 ) in ) which we take into account in the nrg procedure .
for the calculation of dynamical response functions we use the method described in ref .  .
apart from the usual local one - electron green s function @xmath29 we will be interested in two different phonon green s functions , @xmath30 since we will treat only homogeneous phases here and exclude all types of long - range order , we have dropped the site indices in all of the definitions .
the spectral weight of @xmath31 will be denoted by @xmath32 , and similar for all other green s functions .
the dynamical charge and spin susceptibility are given by @xmath33 respectively .
the operator @xmath34 represents the electronic term that couples to the phonons , @xmath35 it should be noted that our definition of @xmath36 differs from the usual one by a factor of four .
we also calculate the quasiparticle weight @xmath37 given by @xmath38 within the local approximation of the dmft .
the different regimes of the phase diagram of fig .
[ phases.eps ] give rise to the following structure of the paper : in sec .
[ sec : g_dependence ] we choose a fixed value of @xmath2 and study dynamical properties as a function of the phonon coupling @xmath7 .
we find continuous transitions from a metallic to a bipolaronic state for small values of @xmath39 and @xmath40 ( sec .
[ ssec : weaku ] ) . for a larger value of @xmath41 , as discussed in sec .
[ ssec : strongu ] , the transition to the bipolaronic state becomes discontinuous .
complementarily , we look at the @xmath2-dependence of dynamical properties in sec .  [
sec : u_dependence ] keeping @xmath7 fixed to @xmath42 .
for this value of @xmath7 the system can be in all three phases , depending on @xmath2 . in sec  [ sec : polaronic ] we study the properties of the model on the two polaronic lines as given by @xmath24 and @xmath43 ( dashed and dot - dashed lines in fig .
[ phases.eps ] ) .
we also show spectra for @xmath24 for different values of @xmath8 . in sec .
[ sec : overview ] we give an overview of our results , and we present our conclusions in sec .
[ sec : conclusions ] .
in this section , we fix the value of the hubbard interaction  @xmath2 and study dynamical properties for a varying electron - phonon coupling @xmath7 .
we distinguish between the cases of weak and strong @xmath2 .       and various values of @xmath7 .
the inset shows the variation of the quasiparticle weight @xmath37 for @xmath40 ( solid line ) and @xmath39 ( dashed line).,scaledwidth=46.0% ]    first of all , we look at the one - electron spectra for a fixed @xmath44 as shown in fig .
[ g_spec_u10_gx.eps ] .
there is a strong similarity of these results with those of the pure holstein model ( @xmath39 ) calculated previously@xcite . at weak coupling
a narrow peak emerges at the fermi energy .
the top of this feature is rather flat and the imaginary part of the self energy has a shallow @xmath45-dependence .
this one would expect from lowest order perturbation theory in @xmath7 , where the imaginary part of the self energy vanishes in the range @xmath46 for @xmath39 ( see eq .
( 19 ) in ref .  ) . upon increasing @xmath7
further , there is a rapid narrowing of this feature until it disappears at a critical coupling of @xmath47 and a gap opens .
similar to the pure holstein model , there is no preformed gap in this case in contrast to the mott transition in the hubbard model .
the metallic regime corresponds to a renormalized fermi liquid .
the quasiparticle weight @xmath37 is shown as a function of @xmath7 in the inset together with that for the pure holstein model ( @xmath39 ) . for @xmath40 ,
the quasiparticles are already slightly renormalized at @xmath48 .
however , initially the renormalization with @xmath7 is weaker than in the case @xmath39 but more rapid on the approach to @xmath49@xcite . as with the pure holstein model
, there is no evidence of multiple solutions near the phase transition .
this is in contrast to the situation near the mott transition for the hubbard model , where there is a finite coexistence region .    , top ) and @xmath40 ( bottom ) .
the values of @xmath7 are @xmath50 ( bottom ) , @xmath51 ( top ) in the upper panel and @xmath50 to @xmath52 in steps of @xmath53 in the lower panel .
the dotted lines are scaled up by a factor of @xmath54.,title="fig:",scaledwidth=23.0% ] , top ) and @xmath40 ( bottom ) .
the values of @xmath7 are @xmath50 ( bottom ) , @xmath51 ( top ) in the upper panel and @xmath50 to @xmath52 in steps of @xmath53 in the lower panel .
the dotted lines are scaled up by a factor of @xmath54.,title="fig:",scaledwidth=23.0% ] + , top ) and @xmath40 ( bottom ) .
the values of @xmath7 are @xmath50 ( bottom ) , @xmath51 ( top ) in the upper panel and @xmath50 to @xmath52 in steps of @xmath53 in the lower panel .
the dotted lines are scaled up by a factor of @xmath54.,title="fig:",scaledwidth=23.0% ] , top ) and @xmath40 ( bottom ) .
the values of @xmath7 are @xmath50 ( bottom ) , @xmath51 ( top ) in the upper panel and @xmath50 to @xmath52 in steps of @xmath53 in the lower panel .
the dotted lines are scaled up by a factor of @xmath54.,title="fig:",scaledwidth=23.0% ]    figure  [ csss_spec_u10_gx.eps ] shows the spectra of the corresponding charge- and spin susceptibilities at @xmath40 and @xmath39 for comparison . in the free system , @xmath55 as defined in eq .
( [ eq : chi_csss ] ) , which can be seen in the lowest curves of the two upper panels . looking at the case @xmath39 first ( upper panel ) , we see the emergence and buildup of a low - energy peak in @xmath56 and the suppression of weight in @xmath57 for @xmath58 . in the bipolaronic state ( @xmath59 ) , both charge and spin susceptibilities have the same gap , roughly twice that seen in the electronic spectrum .
they have an identical peak above this gap .
these can be seen in the dotted curves of fig .
[ csss_spec_u10_gx.eps ] , which have been scaled up by a factor of @xmath54 to make these features visible .
the coincidence of these peaks can be explained by the fact that spin fluctuations necessarily require the breakup of a bipolaronic pair .
their weight is very small due to the strong bipolaronic binding .
there is an additional very small peak in the charge susceptibility in the gap close to @xmath8 from the residual couplings to the phonon mode . in the absence of charge
ordering , there must be charge fluctuations at @xmath60 which , however , can not be observed numerically .    for @xmath40 ( lower panel )
we observe the same overall trends .
however , starting from a reduced charge susceptibility due to the finite @xmath2 , the peak in @xmath56 develops more slowly .
the spin susceptibility , which is enhanced due to the finite @xmath2 , changes very little in the metallic phase .
the low - energy spin fluctuations are almost completely suppressed with the emergence of the gap in the bipolaronic state .
the charge fluctuations can be directly related to the phonon propagator via the equation ( see appendix [ app : d_chic ] ) @xmath61 where @xmath62 is the free phonon propagator .
however , use of this equation with an _ approximate _ form of @xmath63 can lead to severe numerical errors , as discussed in appendix  [ app : d_chic ] . in fig .
[ bb_spec_u10_gx.eps ] , we give the results for the phonon spectra for @xmath39 and @xmath40 calculated via the self energy@xcite , similar to the procedure we use for the electronic spectra@xcite . the results for @xmath39 differ significantly from those published earlier@xcite based on eq .
[ eq : d_chic ] .
in contrast to the earlier results , the main phonon peak is found initially to soften significantly in a similar way to the migdal elisahberg result .
the two - peak structure develops only upon approaching the transition to the bipolaronic state . in the gapped state ,
the high - frequency peak narrows and tends towards @xmath8 .
the same trend can be seen for @xmath40 ( right - hand plot ) .
the slower initital softening correlates with the suppression of charge fluctuations seen in @xmath63 .     in the pure holstein model ( @xmath39 ) . in the metallic range ( @xmath64 ) , the direct calculation ( @xmath65 ) agrees well with the integration over the phonon green s function ( @xmath66 ) and charge susceptibility ( @xmath67 ) ) . for details
see appendices  [ app : d_chic ] and  [ app : n_vs_x ] .
the dot - dashed line shows the result from the migdal eliashberg calculation.,scaledwidth=46.0% ]    on the approach to the transition , significant negative spectral weight for @xmath68 can be seen in the spectra shown in fig .
[ bb_spec_u10_gx.eps ] . as a consequence of the spectral theorem , this relates directly to the average number of excited phonons in the ground state .
the spectra of the green s function @xmath69 show similar features as @xmath31 , and integrated up to @xmath60 yield the average lattice fluctuations as shown in fig .
[ x2_av_u00_gx.eps ] .
these results agree well with the direct evaluation of the expectation value @xmath70 in the metallic phase . for low to intermediate values of @xmath7 ,
the nrg results do not deviate significantly from the midgal eliashberg calculation and are significantly smaller than reported in ref .   based on eq .
[ eq : d_chic ] . on the approach to the transition
, there is a large increase in the lattice fluctuation . in the gapped phase ,
we observe a large difference between the direct calculation and the integration over the negative part of the phonon spectral density .
one way of looking at it is that due to the limited numerical resolution we miss the contribution from the peak at @xmath71 in the spectrum of @xmath69 .
an other way of explaining this difference is that there is no dynamics connecting the two degenerate ground states with @xmath72 and @xmath73 in the gapped phase .
these two ground states have nonvanishing average displacements @xmath74 , see appendix  [ app : n_vs_x ] , and integration over the phonon spectrum yields the fluctuations about @xmath74 , i.e. , @xmath75 .
the calculation of @xmath76 in one of the ground states quantitatively explains the different values of the two methods in the gapped phase .       and various values of @xmath7 .
the central peak broadens with increasing @xmath7 and at @xmath77 disappears in a first order phase transition to the bipolaronic state .
the inset shows the quasiparticle weight @xmath37 ( @xmath66 ) and the double occupancy @xmath78 ( @xmath79 ) as functions of @xmath7.,scaledwidth=46.0% ]    next we look at the transition to the bipolaronic state for a larger fixed value of @xmath41 .
figure  [ g_spec_u50_gx.eps ] shows the corresponding one - particle electron spectral functions for various values of @xmath7 .
their @xmath7-dependence is in sharp contrast to the small @xmath2 case shown in fig .
[ g_spec_u10_gx.eps ] .
initially for @xmath48 , we have the three - peak structure of a strongly correlated hubbard model . with increasing  @xmath7 ,
the central resonance broadens slightly as is reflected in the increasing quasiparticle weight @xmath37 ( see inset ) .
this can be ascribed to the partial cancellation of the hubbard repulsion by phonon mediated electron - electron attraction , see eq .
( [ eq : ueff_def ] ) .
the effect , however , is much weaker than one would deduce from @xmath23 .
the position of the high - energy hubbard bands is virtually unaffected by the phonons , as expected from eq .
( [ eq : ueff_def ] ) evaluated at @xmath80 .
when approaching the critical coupling @xmath81 , the quasiparticle peak disappears abruptly and @xmath37 jumps to zero , indicating a discontinuous transition to the bipolaronic phase .
the double occupancy , shown in the inset to fig .
[ g_spec_u50_gx.eps ] , is much less than a quarter in the metallic state but makes a sudden jump to @xmath82 on the point of the phase transition .
a coexistence of metallic and gapped solutions is found in the narrow range @xmath83 .    .
the values of @xmath7 correspond to those in fig .
[ bb_spec_u50_gx.eps ] .
note the different scales of the @xmath84-axes .
dotted lines have been scaled up by a factor of @xmath54.,title="fig:",scaledwidth=23.0% ] .
the values of @xmath7 correspond to those in fig .
[ bb_spec_u50_gx.eps ] .
note the different scales of the @xmath84-axes .
dotted lines have been scaled up by a factor of @xmath54.,title="fig:",scaledwidth=23.0% ]    the spectra of the charge and spin susceptibilities are shown in fig .
[ csss_spec_u50_gx.eps ] .
they display similar trends , but more pronounced , than the corresponding results for @xmath40 ( lower panels of fig .
[ csss_spec_u10_gx.eps ] ) .
the low - energy peak in @xmath63 only becomes clearly visible when @xmath7 reaches the value of @xmath85 . in the bipolaronic phase , the peaks above the gap in the charge and spin susceptibilities are not identical as in the @xmath39 case , but very similar
however , the gaps are too large to be visible on the @xmath18-scale of fig .
[ csss_spec_u50_gx.eps ] .
the only peak visible , when enhanced by a factor of @xmath86 , is the peak in the charge susceptibility within the gap near @xmath8 .     and various values of @xmath7.,scaledwidth=46.0% ]    in fig .
[ bb_spec_u50_gx.eps ] , the spectra of the corresponding phonon propagators are plotted . over a large initial range of @xmath87
there is very little softening , corresponding to a complete suppression of charge fluctuations .
when approaching @xmath49 , the phonon mode softens and develops a low - energy shoulder .
after the discontinuous transition to the bipolaronic phase , the mode hardens back to @xmath8 and narrows .
in this section we will fix the electron - phonon coupling at the value @xmath42 and study the dependence of dynamic response functions on the hubbard repulsion  @xmath2 .
the value of @xmath42 is chosen because all three phases can be found for it , depending on @xmath2 ( see fig .
[ phases.eps ] ) .
figure  [ g_spec_ux_g45.eps ] shows the one - electron spectral function in the vicinity of both phase transitions .
the right - hand plot for large values of @xmath2 is very similar to what is observed in the mott transition of the pure hubbard model .
the lower and the upper hubbard peak at @xmath88 are well developed and move to higher energies as @xmath2 increases . at the same time
, the central peak narrows and vanishes at the mott transition . as for the pure hubbard model
, there is a preformed gap in the one - electron spectrum and a broad region of numerical coexistance of metallic and insulating solutions .
the electronic spectra are very similar to those of the pure hubbard model near the mott transition .
the phonons do not seem to alter much the picture of the phase transition .
the transition from the bipolaronic to the metallic state is completely different , as can be seen in the left - hand plot of fig .
[ g_spec_ux_g45.eps ] .
the phase transition occurs close to @xmath40 and the spectra show a similar behaviour as discussed in sec .
[ ssec : weaku ] , where @xmath7 is varied . in the bipolaronic state , with the initial increase of @xmath2 , the two polaron bands move towards the fermi level . on the metallic side , just after the transition , we see a very sharp resonance which broadens when @xmath2 is increased further . there is no signature of a preformed gap once the system is metallic .
again , there is no evidence of a significant coexistence region .
the spectra of the charge and spin susceptibilities for @xmath42 are shown in fig .
[ csss_spec_ux_g45.eps ] .
we clearly see that the low - energy features of the two susceptibilities signal the appearance of two different instabilities of the system , depending on @xmath2 . in the bipolaronic phase
, both spectra are gapped and we see only the peak near @xmath8 in the charge susceptibility . upon entering the metallic phase , the zero frequency peak in the charge susceptibility moves to finite @xmath18 .
this peak loses weight and slightly disperses to higher energies as the system becomes more metallic . in this small @xmath2 regime ,
spin fluctuations are suppressed .
the converse occurs when approaching the mott transition at @xmath89 .
the charge fluctuations are suppressed and peak in the spin susceptibility builds up , signalling instability towards antiferromagnetic ordering . in the mott insulating state
, both susceptibilities are gapped .
the spin fluctuations above the gap ( beyond the scale of fig .
[ csss_spec_ux_g45.eps ] ) must largely be due to the charge fluctuations , as in the pure hubbard case .
however these peaks are not identical , in contrast to the bipolaronic gapped state .
the difference in the peaks can be interpreted as a result of the spin correlations induced by @xmath2 . however , due to the phonon coupling the spin fluctuation peak is reduced compared to the pure hubbard case .     and various values of @xmath2 .
the transition from the bipolaronic to the metallic state is clearly visible by the two - peak structure whereas the transition to the mott insulator does not obviously affect the phonon propagator.,scaledwidth=46.0% ]    figure  [ bb_spec_ux_g45.eps ] displays the phonon propagator for this scan . in the bipolaronic phase at @xmath39 , electrons and phonons are almost decoupled , and the phonon spectrum shows only one peak just below @xmath14 .
when increasing @xmath2 , this mode softens and splits into two peaks close to the phase transition . on the metallic side , the low - energy peak survives and continuously hardens back to @xmath8 , as the system approaches the mott transition .
there is no signature of the mott transition in the phonon spectrum .
the effect of an increasing @xmath2 is just to suppress continuously the charge fluctuations which results in a decoupling of electrons and phonons which drives the hardening of the phonon peak .
another feature of the suppression of charge fluctuations is the narrowing of the phonon peak with increase of @xmath2 , as also observed in ref .  .
in our case , however , the shift of the peak appears to be more marked .
the narrowing can also be seen in comparing the curves with the same ( low ) value of @xmath7 in figs .
[ bb_spec_u10_gx.eps ] and [ bb_spec_u50_gx.eps ] .    .
the weight for @xmath90 indicates the degree of double occupancy ; @xmath42 and various values of @xmath2 .
the inset shows the double occupancy calculated directly ( @xmath67 ) and via the spectral theorem from @xmath91 ( @xmath79 ) .
circles indicate the total spectral weight of @xmath91 which should evaluate to @xmath92 at half filling .
the discrepancies are most pronounced near the phase transitions at @xmath93 and @xmath94.,scaledwidth=46.0% ]    further insight into the nature of these transitions can be gained by looking at the behaviour of the higher order electron green s function @xmath95 this green s function is needed for the calculation of the self energy in the nrg procedure@xcite .
the integrated spectral weight for @xmath96 equals the double occupancy @xmath97@xcite .
moreover , the total spectral weight of @xmath98 yields the average density @xmath99 , which should be @xmath100 in the particle - hole symmetric case treated here .
figure  [ f_spec_ux_g45.eps ] shows the spectral function of @xmath91 for @xmath42 . in the bipolaronic phase ,
almost all of the weight is located below the fermi energy signalling a high value of the double occupancy ( see inset ) .
as we enter the metallic phase , a large part of this weight is transferred to higher energies . in the mott insulator ,
virtually no weight is left at @xmath68 and the double occupancy is completely suppressed .
the inset to fig .
[ f_spec_ux_g45.eps ] also shows the total spectral weight of @xmath91 , which should evaluate to @xmath101 .
we see that this is indeed more or less the case , except for values of @xmath2 close to the phase transitions @xmath93 and @xmath94 . for these values
, the truncation of the hilbert space in the nrg introduces the most significant errors .
there we also find the strongest discrepancies between the direct calculation of the double occupancy from the ground state expectation value ( marked with @xmath67 ) and the indirect , less accurate calculation via @xmath91 ( marked with @xmath79 ) .
in this section we discuss the dynamical properties along two particular lines in the @xmath102 plane ( see fig .  [ phases.eps ] ) .
the first of these lines is given by @xmath103 where one might naively expect free quasiparticles close to the fermi surface .
the other line is the location of points where the double occupancy @xmath43 , as in the free system .    only for very large values of @xmath8
can the hubbard term be completely cancelled by the phonon - mediated attractive term ( see eq .
[ eq : ueff_def ] ) , in which case both polaronic lines would coincide .
we investigate these lines to understand the degree of compensation of these competing interactions for small phonon frequency @xmath104 .       at @xmath14 (
dashed line in fig .
[ phases.eps ] ) .
the inset shows the variation of the quasiparticle weight  @xmath37 and double occupancy  @xmath78.,scaledwidth=46.0% ]    figure  [ g_spec_ueff0_ux.eps ] shows the one - electron spectral functions for various values of @xmath7 and @xmath2 such that @xmath24 . for small values of @xmath2 ( and @xmath7 ) , we see the sharp feature in the electronic spectrum as discussed in sec .
[ sec : g_dependence ] .
however , upon increasing @xmath2 further , hubbard bands develop and the central feature becomes the central quasiparticle peak similar to the pure hubbard model for large @xmath2(@xmath105 ) .
this indicates that phonons play their assumed role of compensating the hubbard repulsion only to a very limited extent .
the only remarkable difference to the pure hubbard model is that on the line @xmath24 the phase transition is shifted towards a higher value of @xmath106 , as compared to @xmath107 for @xmath48 .
the inset to fig .
[ g_spec_ueff0_ux.eps ] shows that the quasiparticle weight decreases steadily from @xmath108 to zero and the double occupancy from @xmath43 to a small but finite value .    .
the inset shows the expectation value of the lattice fluctuations as a function of @xmath2.,scaledwidth=46.0% ]    the phonon spectra , as displayed in fig .
[ bb_spec_ueff0_ux.eps ] , are always made up of only one single peak that initially moves towards lower energies as more phonons are excited and lattice fluctuations increase ( see inset ) . at @xmath109
this trend reverses and the phonon peak hardens back when approaching the phase transition to the gapped state .
simultaneously , the number of excited phonons decreases to zero and the lattice fluctuations reach the value of the quantum mechanical zero point fluctuations .
the reason for this non - monotonic behaviour of the lattice fluctuations as function of @xmath7 is not obvious .
the maximum in the lattice fluctuations as shown in the inset correlates with the value of @xmath2 where the transition to the bipolaronic state changes from second order ( for @xmath110 ) to first order ( for @xmath111 ) .
this could be a possible explanation as one would expect to see stronger lattice fluctuations as one aproaches a second - order transition to the bipolaronic state .      .
( dot - dashed line in fig .
[ phases.eps ] ) .
the inset shows the variation of the quasiparticle weight  @xmath37 along this line.,scaledwidth=46.0% ]    the line in the phase diagram where @xmath113 starts at the uncorrelated system @xmath114 and ends at @xmath115 where it merges with the phase boundary to the bipolaronic state . as a consequence ,
for larger values of @xmath2 , the double occupancy has a jump at the transition from the metallic state ( @xmath116 ) to the bipolaronic state ( @xmath117 ) ( see fig .  3 in ref .   and inset of fig .  [ g_spec_u50_gx.eps ] ) .
therefore this transition has to be first order in this range .
figure  [ g_spec_d25_ux.eps ] shows the one - electron spectral functions along this line .
we observe the appearance of the narrow feature at the fermi energy that is a typical for the electron - phonon coupling .
the band develops broad shoulders , but they never separate into distinct subbands .    .
the inset shows the expectation value of the lattice fluctuations as a function of @xmath2 along this line.,scaledwidth=46.0% ]    the phonon propagator , as displayed in fig .
[ bb_spec_d25_ux.eps ] shows a remarkable softening of the original peak and the build - up of negative spectral weight . as a consequence ,
the lattice fluctuations increase almost linearly with @xmath2 and are quite pronounced for @xmath118 .       and various values of @xmath8.,scaledwidth=46.0% ]
so far we have worked with a fixed value of @xmath14 . however ,
the degree of compensation of the competing interactions is dependent on @xmath8 .
here we return to the case @xmath24 and study the variation of @xmath8 with fixed @xmath41 and @xmath119 .    in fig .
[ g_spec_ueff0_u50_wx.eps ] we plot the one - electron spectral density for a range of values of @xmath8 . for the largest value of @xmath120
, the spectrum is virtually the same as that of the free system , because we have almost complete compensation .
we now examine the behaviour as we progressively reduce @xmath8 . first the high - energy phonon subbands gain weight and move towards the fermi energy .
their position is roughly given by @xmath8 .
there is a commensurate reduction of the width of the central peak . for @xmath121
the phonon subbands are visible as shoulders of the central peak . upon further reduction of @xmath8
, these shoulders move back towards higher energies again and , for @xmath122 , become the hubbard bands located at @xmath123 . for the smallest value of @xmath14 , the spectrum is identical to the one discussed in sec .
[ ssec : ueff0 ] .     and @xmath41 for various values of @xmath8.,scaledwidth=46.0% ]    in fig .
[ bb_spec_ueff0_u50_wx.eps ] we plot the spectra of the phonon propagator as a function of the relative energy scale @xmath124 . for large @xmath125 , where the compensation is almost complete ,
we see the narrow phonon peak at @xmath8 .
in addition we see a small feature near @xmath60 . as @xmath8
is decreased , the low - energy feature gains weight and the mode at @xmath8 broadens and moves to slightly higher energies , as one would expect when @xmath8 approaches but lies above the upper band edge . at @xmath126 ,
the two features have approximately equal weight . for smaller values of @xmath8
the upper peak has almost zero weight .
the remaining peak is the soft mode discussed in the previous sections .
( o ) , the rescaled double occupancy @xmath127 ( @xmath79 ) , and lattice fluctuations @xmath128 ( @xmath67 ) as a function of @xmath8 for @xmath24 .
the dot - dashed line represents the number of excited phonons @xmath129.,scaledwidth=46.0% ]    in fig .  [ zdp_ueff0_u50_wx.eps ]
we illustrate how different static quantities reflect the trend of the holstein - hubbard model towards a free system in the limit of large @xmath8 .
we plot @xmath37 and @xmath127 as a function of @xmath8 .
the double occupation value converges relatively rapidly to that of the free state .
the quasiparticle weight , however , converges rather more slowly .
in fact , only for @xmath130 is it approaching the uncorrelated limit @xmath108 .    in the same fig .  [ zdp_ueff0_u50_wx.eps ]
we also plot the lattice fluctuations @xmath131 .
this curve goes monotonically to zero as @xmath22 asymptotically behaving as @xmath132 for large @xmath8 .
however , if we plot the rescaled values @xmath133 , we find a a maximum which coincides with the region of maximum softening of the phonon mode .
for larger values of @xmath8 this quantity decreases to the quantum limit of @xmath134 .
one would expect the phonons to decouple from the electrons in the limit of large @xmath8 .
it is interesting to observe , however , that @xmath135 ( dot - dashed line in fig .
[ zdp_ueff0_u50_wx.eps ] ) fits well the high-@xmath8 behaviour of the lattice fluctuations which indicates that the ground state of the phonons in this limit is a coherent state rather than an eigenstate of the free phonon hamiltonian @xmath136 .
for an overview we will find it useful to classify the metallic state in the phase diagram into three regions , separated by the two polaronic lines , in each of these we find qualitatively different behaviour .    we consider first of all the region below the line @xmath24 . here
one interesting question is how the coupling to the phonons affects the transition to the mott insulating state .
our results indicate that this transition is always very similar to that found in the pure hubbard model : the appearance of a preformed gap , a rather large coexistence region between insulating and metallic solutions and the quasiparticle weight continuously going to zero .
there is only a small increase in the critical value of @xmath2 which is found empirically to depend upon @xmath7 to good approximation as @xmath138 the small influence of the phonon coupling can be readily explained by suppression of charge fluctuations for large @xmath2 . as the electron - phonon coupling in the holstein - hubbard model is to the charge fluctuations ,
the suppression of these renders the coupling ineffective .
in fact , the whole region @xmath137 appears to be dominated by the hubbard term . for @xmath139 the second term in the retarded interaction as defined in eq .
[ eq : ueff_def ] is negligible and the effective interaction is essentially given by @xmath2 . also on the lowest energy scales
@xmath140 , @xmath141 .
next we look at the complementary region defined by @xmath142 which is dominated by the coupling to phonons . here
the similarity is with the pure holstein model .
the metal to bipolaronic transition takes place without a preformed gap and with the characteristic softening of the phonon mode .
the transition in this region appears also to be continuous with no significant coexistence region .
the main effect of the hubbard  @xmath2 is to push the phase transition and its precursors to somewhat larger values of @xmath7 .
finally , in the region bounded by the two curves @xmath24 and @xmath43 there is more complex interplay of the two interactions , as can be seen in our results on the response functions . beyond the point where the line @xmath43 merges into the phase boundary the transition to the bipolaronic state
is first order and thus qualitatively different from that discussed in the last paragraph . along the lower boundary line of this region ,
i.e. , for @xmath24 , the effective interaction is still repulsive for @xmath139 .
the question arises whether on the lowest energy scale there is a complete cancellation of the two competing interactions
. this question can be investigated by examining the lowest - lying one- and two - particle excitations from the interacting ground state .
the calculation of the effective interaction @xmath143 between the renormalized quasiparticles has been considered in detail for impurity models in ref .  .
a similar analysis is possible using the effective impurity model within the dmft . from ref .
, @xmath144 where @xmath145 ( @xmath146 ) is the energy of the lowest - lying single - particle ( particle - particle ) excitation in the @xmath147-th iteration of the nrg procedure and @xmath148 is the nrg discretization parameter . as for the impurity model , we find that , within the numerical accuracy of the nrg , @xmath143 vanishes on the line @xmath24 and is positive ( negative ) below ( above ) this line
. therefore , along this line we have a system of renormalized ( @xmath149 ) but non - interacting quasiparticles . on the high - energy scale @xmath150 ,
the system is still dominated by @xmath2 .
this can be seen by the formation of the hubbard bands as shown in fig .
[ g_spec_ueff0_ux.eps ] and the suppression of the double occupancy .
in this study we have investigated the particle - hole symmetric holstein - hubbard model .
we have used the dynamical mean field theory in combination with the numerical renormalization group to calculate dynamical correlation functions for the full hamiltonian with quantum phonons .
this non - perturbative approach allows us to access all parameter regimes of the model at zero temperature .
a non - perturbative technique is indeed necessary because perturbative methods such as midgal - eliashberg approach for the holstein model are known to break down@xcite in the strong - coupling limit .
we find three regions with qualitatively different behaviour . in the region of the mott metal insulator transition
we find the phonon effects are largely suppressed by the on - site repulsion @xmath2
. this conclusion may be specific for the non - degenerate holstein - hubbard model where the coupling is purely to the local charge density . a system with orbital degeneracy and jahn - teller phonons shows@xcite strong phonon effects even for strong @xmath2 . in our case of the
non - degenerate holstein - hubbard model , the electronic effects on phonons yield only a modest softening of the phonon mode in the metallic regime .
this softening disappears in the mott insulator .
weak electron - electron interactions seem to have little effect other than to delay the onset of the transition to the bipolaronic phase . for larger values of @xmath2
the transition changes from continuous to discontinuous . for small values of @xmath2
we see a complete softening of the phonon peak when approaching the transition .
the softening is much less pronounced for larger @xmath2 .
the line @xmath24 in the metallic phase does seem to have some significance in that for @xmath137 the spin susceptibility has a dominant low - energy peak .
for @xmath151 the dominant low - energy peak is in the charge susceptibility .
the line @xmath24 also appears to be significant for the low - energy fermi liquid behaviour as the on - site quasiparticle interaction changes sign on or close to this line .
this is not obvious as one might have expected that on the low energy scale of the electron - phonon interaction , the electron - electron interaction is no longer given by the bare @xmath2 but some renormalized value @xmath152 .
this would imply that one should use @xmath152 in eq .
[ eq : ueff_def ] on the energy scale @xmath140 , so that the change of sign of the quasiparticle interaction @xmath143 , as defined in eq .
( [ eq : tildeueff ] ) , would occur when @xmath153 .
however , our results indicate that this is not the case and that @xmath143 vanishes on or close to the line @xmath154 .
a similar situation was found in the case the anderson - holstein impurity model@xcite .
the holstein - hubbard model is very rich and shows also diverse forms of behaviour that could not be addressed in this paper .
apart from the results presented here , this model is known to exhibit various types of symmetry - breaking phases@xcite .
we intend to extend our calculations to investigate the competition between charge order , antiferromagnetism , superconductivity and the effects of doping@xcite .
this should enable us to make direct contact with experimental results on such compounds as the fullerides , where this model with the assumption of the coupling of the electron density to the local phonon mode should be directly applicable . for applications to jahn - teller systems , such as in the manganites , the model would have to be extended to include degenerate orbitals .
we wish to thank the epsrc ( grant gr / s18571/01 ) for financial support .
this work was partially supported by sunnynames llp .
we acknowledge fruitful discussions with a. gogolin , d. edwards and y. = ono .
here we reconsider the relation between the phonon propagator and the charge susceptibility ( eq .  ( [ eq : d_chic ] ) ) for the holstein - hubbard model and comment on pitfalls in its exploitation .
[ [ equations - of - motion ] ] equations of motion : + + + + + + + + + + + + + + + + + + + +    for the derivation of eq .
( [ eq : d_chic ] ) , we can use the standard equation of motion for the fourier transform of the double - time green s function , @xmath156_\eta \rangle+ { { \big<\!\big < { [ a , h]_{-}\,;b } \big>\!\big > } } _ { \omega}\ ] ] for the boson green s function @xmath31 we take @xmath157 , @xmath158 , and @xmath159 . for @xmath69 we take @xmath160 and @xmath159 . the following procedure works for both @xmath31 and @xmath69 , but we will use the former in what follows .
the equation of motion gives @xmath161 where @xmath162 .
taking the equation of motion for the right - hand operator , we have @xmath163 hence the result ( @xmath164 is the free phonon propagator ) , @xmath165 which is eq .
( [ eq : d_chic ] ) linking the phonon propagator @xmath31 with the charge susceptibility @xmath166 .    as a green s function
, the phonon propagator has a series of poles of first order , but none of second order .
for this to be the case , eq .
( [ eq : d_chic ] ) tells us that the charge susceptibility has to vanish as @xmath167 in order to ensure that the pole of @xmath31 at @xmath168 is first order .
this is also evident if we express @xmath63 in terms of the irreducible particle - hole bubble @xmath169 .
it takes the form @xmath170 with the non - interacting phonon propagator @xmath171 .
since @xmath172 diverges at @xmath168 , we have @xmath173 .    [ [ spectral - density ] ] spectral density : + + + + + + + + + + + + + + + + +    we will derive an expression for the spectral density of the phonon propagator , i.e. , for @xmath174 that implicitly includes the fact that @xmath173 .    consider the system with a finite basis such that , @xmath175 for @xmath176 , where @xmath177 are the discrete set of excitations in the charge density response function . then @xmath178 the contributions to the spectral density @xmath179 arise solely from the poles at @xmath168 and @xmath180 .
these do not coincide because @xmath173 .
the contribution to @xmath179 from the poles at @xmath180 are straightforward to evaluate and give @xmath181 the contribution from the pole at @xmath168 is a little trickier to evaluate .
the term in @xmath182 contributes @xmath183 where @xmath184 denotes the derivative of the delta function .
however , we have @xmath185 where @xmath186 is the derivative of @xmath63 and is given by @xmath187 using @xmath173 and collecting all the terms together we obtain @xmath188 this expression can be used to calculate the spectral density of the phonon propagator @xmath179 from the peaks @xmath189 of the charge susceptibility @xmath63 .
in fact , the numerical evaluation of eq .
( [ eq : rho_b ] ) yields the correct result for the phonon propagator and the lattice fluctuations , as shown in fig .
[ x2_av_u00_gx.eps ] for the case of the pure holstein model .
the condition @xmath173 has , of course , been used explicitly in the derivation of eq .
( [ eq : rho_b ] ) .    the problem of the direct evaluation of eq .
( [ eq : d_chic ] ) for the phonon propagator is the numerical error in @xmath63 .
in fact , the truncation of the hilbert space in the nrg procedure entails that the condition ( [ eq : chic_w0 ] ) is not exactly met .
therefore @xmath63 has no root exactly at @xmath168 .
this leads to a double pole of @xmath31 at @xmath168 and to a distorted result for the phonon propagator at energies close to @xmath8 , when derived from eq .
( [ eq : d_chic ] ) .
lattice fluctuation , as calculated by integrating @xmath69 from eq .
( [ eq : d_chic ] ) are strongly overestimated for the same reason .
in this appendix we derive an expression linking the average local displacement @xmath191 to the local electron density . for this derivation , we modify the hamiltonian by coupling an extra term linear in @xmath192 at each site .
the extra term in the hamiltonian reads @xmath193 with the coupling @xmath194 . upon differentiating the expression for the partition function @xmath195 with respect to @xmath196
, we obtain @xmath197 we then apply a canonical transformation to the hamiltonian , @xmath198 , with @xmath199 given by @xmath200 this is a displaced oscillator transformation at each lattice site . the phonon and electron operators are transformed as @xmath201 after some algebra we find that the terms depending on @xmath202 in the transformed hamiltonian @xmath203 read @xmath204 the partition function @xmath195 is , of course , not changed by the canonical transformation of the hamiltonian . upon expressing @xmath195 in terms of @xmath205 , differentiating with respect to @xmath196 and then putting @xmath206 we find the result @xmath207 which relates the average local displacement to the local average electron density