the problem of a single hole in the @xmath1-@xmath2 model interacting also with the lattice degrees of freedom has recently attracted a notable interest in connection with the physical properties of the high-@xmath4 cuprates .
@xcite in particular , in parent and strongly underdoped compounds , angle resolved photoemission spectroscopy ( arpes ) reveals a low energy peak whose dispersion is well described by the @xmath1-@xmath2 model , while its anomalously large broadness has been ascribed to incoherent multi - phonon shake - off processes.@xcite a similar interplay between electron - electron and electron - phonon interactions should in principle be reflected in the optical conductivity spectra . as a matter of fact ,
the most remarkable features observed in the underdoped region are an ubiquitous mid - infrared ( mir ) peak at @xmath5 ev , and a weaker peak around @xmath6 ev , the latter being more strongly doping and temperature dependent.@xcite several interpretations have been proposed for the origin of these features , including midgap or impurity states , @xcite charge / spin stripes,@xcite polaronic excitations,@xcite and the interaction with the antiferromagnetic background.@xcite    the first proposal of the @xmath1-@xmath2 model as a suitable basis to discuss the optical spectra of the cuprates was advanced by zhang and rice @xcite who observed that the @xmath7 behavior@xcite of the optical absorption above @xmath5 ev could be naturally associated to the incoherent motion in an antiferromagnetic background .
successively , the optical conductivity of the @xmath1-@xmath2 model has been investigated in detail using several techniques , such as exact - diagonalization,@xcite analytical approximations@xcite and dynamical mean - field theory ( dmft).@xcite however , and in spite of the above discussed relevance of the electron - phonon coupling , the optical conductivity of the @xmath1-@xmath2 model _ in the presence _ of the lattice degrees of freedom has not been thoroughly investigated .
numerical calculations based on exact - diagonalization of finite clusters were employed for instance in ref . to evaluate @xmath0 .
alternatively , the optical conductivity was calculated analytically in ref .
based on a non - crossing born approximation , which is however unable to describe the polaron formation .    in this paper
we present results for the optical conductivity of a single hole in the holstein-@xmath1-@xmath2 model obtained in the framework of dynamical mean - field theory .
one - hole spectral properties at zero temperature were discussed in a previous publication where antiferromagnetic correlations were shown to enhance the effects of the electron - phonon coupling.@xcite a similar result was found also in the antiferromagnetic phase of the holstein - hubbard model using dmft techniques.@xcite a serious drawback of dmft , which is obtained as the exact solution of the lattice problem in the limit of infinite dimensions , is that the magnetic background is treated in a classical way .
this , together with the fact that trugman loops are negligible in infinite dimensions , prevents the possibility to account for coherent hole - propagation , which is related to the spin - flip fluctuations . as a consequence
, no drude peak can be observed in the optical conductivity . on the other hand ,
the incoherent contributions of @xmath0 are mainly dominated by _
properties , such as the local electron - phonon scattering and spin - string excitations within the magnetic polaron , which are well captured by this approach.@xcite    bearing the above limitations in mind , the aim of the present work is thus to focus on the incoherent part of the optical conductivity of a single - hole and to investigate in detail its features in the different physical regimes of the holstein-@xmath1-@xmath2 model .
the dependence of the optical spectra on the microscopic parameters is analyzed with special attention to the intermediate coupling region , where the interplay between magnetic and lattice degrees of freedom is strongest .
we show that a crucial role is played by the formation of the lattice polaron , which drives a transfer of spectral weight towards higher frequencies , opening a pseudogap at low frequencies .
conversely , starting from the polaronic phase , two different mechanisms can be clearly identified as being responsible for the disappearance of the pseudogap : @xmath8 ) reducing the effective exchange energy scale suppresses the positive feedback of magnetism on polaron formation and can lead to a _ closing _ of the pseudogap as the system crosses back to the non - polaronic regime ; @xmath9 ) increasing the temperature , which does not alter the lattice / magnetic interplay , leads to a _ filling _ of the pseudogap more similar to what is expected in purely polaronic models . in the immediate vicinity of the polaron crossover ,
the spectra are characterized by a coexistence of a magnon peak at low frequency and a broad polaronic band at higher frequency , which closely resembles the experimental situation observed in the cuprates .    on theoretical grounds ,
the definition of the optical conductivity of a single hole is a delicate matter which needs particular care .
we provide an analytical derivation which generalizes the results of refs . to the holstein-@xmath1-@xmath2 model .
this approach permits us to identify the role of the different one - particle properties on the optical conductivity .
comparison with numerical data is also discussed , showing a good agreement between our findings and exact diagonalization results .
the paper is organized as follows . in sec .
[ s - model ] we discuss the exact solution in infinite dimensions for the one - hole green s function of the holstein-@xmath1-@xmath2 model at finite temperature . results for the one - particle spectral features are discussed in sec .
[ s - spectral ] . an analytical expression for the optical conductivity @xmath0
is derived in sec .
[ s - opt ] where we also investigate the different polaronic features and their dependence on the microscopic parameters . in sec .
[ s - sc ] we present a further simplified expression for @xmath0 valid in the lattice polaron regime .
the main results are briefly summarized in sec .
[ s - summary ] where also the consequences of including spin fluctuations ( here neglected ) is also discussed .
finally , a detailed derivation of the analytical expression for the one - particle spectral function and the optical conductivity is reported in the appendices .
in the following we consider the case of a single hole in an antiferromagnetic ( af ) background interacting with local holstein phonons . using the linear spin - wave approximation@xcite and neglecting terms that vanish in the limit of large coordination number @xmath10 ,
we can write the hamiltonian as:@xcite @xmath11 here @xmath12 is the creation operator for boson spin defects , @xmath13 is the single spinless hole operator and @xmath14 represents the density of spin defects which is finite at nonzero temperature .
note that in the thermodynamical limit the presence of a single hole does not affect the magnetic state , which can be thus evaluated ( in the @xmath10 limit ) in the absence of spin dynamics .
the density of spin defects can be obtained from the magnetization @xmath15 via the curie - weiss equation : @xmath16 which defines a nel temperature @xmath17 .
concerning the electron - lattice interaction , we shall mainly focus on the adiabatic regime @xmath18 which is relevant to the experimental systems of interest . in this regime ,
a dimensionless electron - phonon coupling constant can be defined as @xmath19 , the polaron energy in units of the hopping integral .
an exact solution for the thermodynamical and the one - particle spectral properties of eq .
( [ ham ] ) at @xmath20 was obtained in ref . in terms of a continued fraction . in order to derive the optical conductivity ,
the one - particle green s function must be generalized to finite temperature , which involves the following steps : @xmath8 ) one has to allow for thermally excited phonons ; @xmath9 ) the presence of thermally excited spin defects requires the introduction of a `` spin - resolved '' green s function , to distinguish hole excitations created on sites with / without spin defects ; @xmath21 ) finally , the reduced magnetization introduces an effective exchange coupling @xmath22 .
the details of a formal derivation of the one - particle propagator at finite temperature are reported in appendix [ app - green - t ] ; we summarize here the main results .    following ref .
, we define @xmath23 as the green s function for one hole created on a site _ in the absence _ of spin defects .
a careful analysis ( see appendix [ app - green - t ] ) shows that the green s function @xmath24 for a hole created on a site with a spin defect is simply @xmath25 .
in addition , at finite temperature the green s function @xmath26 itself is defined as a thermal average over the phonons : @xmath27 where @xmath28 represents the propagation of one hole created on a site with @xmath29 excited phonons and @xmath30 is the single - site phonon partition function @xmath31 .
following ref . , we can derive a self - consistent expression for @xmath32 in terms of a continued fraction .
we can write : @xmath33 where    @xmath34    carrier motion and exchange interactions are taken into account by the bath propagator @xmath35 where the `` hopping '' self - energy is defined as @xmath36 .
\label{sigma_t}\ ] ] finally , the spin - defect averaged green s function , which is the physically probed quantity in photoemission experiments , is obtained as : @xmath37 note that the factors @xmath38 and @xmath39 in front of @xmath40 and @xmath41 account for the probability of a site to be , respectively , free or populated by a spin defect .
note also that the green s function @xmath42 appearing in eqs .
( [ bathg],[sigma_t ] ) is a phonon averaged quantity , so that the solution of eqs .
( [ g - gnn])-([gtot ] ) involves the simultaneous self - consistency of all of the @xmath43 , which is much more computationally expensive than the solution at zero temperature .
it is easy to check that in the absence of electron - phonon interaction eq .
( [ gtot ] ) recovers the thermal green s function for the pure @xmath1-@xmath2 model as defined in eq .
( 2.12 ) by stumpf and logan.@xcite on the other hand , it should be stressed that the present solution , although it formally recovers the results for the pure holstein model in the limit @xmath44 ( paramagnetic case ) [ eqs .
( 40)-(42 ) of ref .
] , is still described by a purely local green s function @xmath45 due to the assumption of a classical ( although disordered ) spin background .
the main drawback of this assumption is thus that no coherent dispersive peak is obtained in this framework even in the @xmath46 limit of the @xmath1-@xmath2 holstein model , and , consequently , no drude peak appears in the optical conductivity .
notwithstanding this limitation the present approach , as we will show below , is still quite able to reproduce in more than qualitative agreement the incoherent features of both the one - particle green s function and of the optical conductivity .
before discussing the optical conductivity , let us briefly present our results for the one - particle spectral function in the holstein-@xmath1-@xmath2 model at finite temperature .
previous approaches have focused separately on thermal effects either in the pure holstein or in the pure @xmath1-@xmath2 model . the temperature evolution of the hole spectral function @xmath47 for the @xmath1-@xmath2 model in the infinite dimensional limit has been analyzed in ref . .
at @xmath20 , it consists of a series of @xmath48-function magnon peaks whose distribution reflects the strength of the magnetic polaron : for sufficiently large @xmath49 the profile is rapidly decaying with energy ( reminiscent of the string picture of small magnetic polarons ) whereas for small @xmath49 it acquires a more symmetric shape , reducing to a semi - circular function in the limit @xmath50 .
the effect of a non - zero temperature within this context is to broaden each @xmath48-function with a bandwidth @xmath51 which is ruled by the thermal spin defect probability @xmath39 .
we shall term this effect the _ intrinsic magnetic broadening_. such intrinsic magnetic broadening @xmath51 is however exponentially small for @xmath52 , and it becomes significant only close to @xmath53 .
a semi - circular shape is recovered also in the paramagnetic phase @xmath54 , where @xmath55 .
as the electron - lattice coupling is turned on , each magnon peak splits into several sub - peaks spaced by @xmath56 , reflecting the dressing of the hole by phononic excitations . in this context
the strength of the electron - phonon coupling rules not only the number of phonon satellites but also its spectral weight profile .
just as in the pure holstein model , while the number of phonon peaks is quite small in the weak coupling regime , in the polaronic state a large number of phonon satellites appear with a characteristic gaussian profile .
the envelope of the phononic fine structure has a spread which is governed by the energy associated with the lattice fluctuations : it is given by @xmath57 in the quantum limit , and increases @xmath58 as @xmath59 .
@xcite    in fig .
we show the temperature evolution of a typical polaronic spectral function , at @xmath60 , @xmath61 and @xmath62 . throughout the paper ,
when not specified , we shall take the hopping matrix element @xmath1 as the energy unit .
the temperatures considered here are @xmath63 , corresponding to @xmath64 . at the lowest temperature ( @xmath65 )
a phononic fine structure can be clearly seen , superimposed on the magnon peaks .
the width of each phononic peak is due to the intrinsic magnetic broadening @xmath51 described in ref . .
it is exponentially small at this temperature , so that a small lorentzian broadening @xmath66 has been introduced for clarity . on the other hand ,
the spread of the multiphonon structure gives rise to an overall width to the magnon peaks that is in good agreement with the expected value @xmath67 .    upon increasing the temperature to @xmath68 ,
two different effects are visible .
first , the width @xmath51 of each of the fine peaks increases due to the intrinsic magnetic broadening , leading to a much smoother curve ( this effect actually overcomes the small lorentzian broadening @xmath69 introduced previously ) . also , the overall spread of the multiphonon profiles increases due to the thermal phonon fluctuations , as expected for @xmath70
finally , at @xmath71 , the system is so close to the paramagnetic phase that neither the phonon peaks nor the magnon structure can be resolved .    at this point
, it is useful to comment about the effects of coherent hole propagation , that are implicitly neglected in our approach .
these would induce a finite dispersion of order @xmath72 to the lowest energy magnon peak.@xcite it is clear that such dispersion would be visible only at sufficiently low temperatures and at moderate electron - phonon coupling strengths , when the energy scale @xmath2 is smaller than both the intrinsic magnetic broadening and the gaussian phonon spread , whereas in the opposite case it will presumably be hidden below a featureless background .
these considerations give further support to the present dmft approach in the polaronic and/or high temperature regime , where neglecting the coherent hole propagation would not affect significantly the spectral properties .
as we shall see below , this is even more true for what regards the finite - frequency optical conductivity , where any dispersive peak would be convolved in any case with high - energy featureless structures .
the evaluation of the optical conductivity for a single hole is a delicate matter , which is only partly simplified in the context of dmft due to the absence of vertex corrections.@xcite a controlled procedure is derived in terms of an expansion of the inverse fugacity at finite temperature , performing the limit @xmath73 to enforce the thermodynamically vanishing particle density .
we can thus define the optical conductivity per hole , @xmath74 , which is a finite quantity and which presents the same features as in the dilute ( but finite ) hole density limit . applying this formalism
one derives a similar expression as obtained in refs .
, here adapted to take into account the electron - phonon interaction .
leaving once more the technical details in appendix [ app - sigma ] , we report here the main results .
the optical conductivity per hole is expressed as @xmath75 ,   \label{sigma } \end{aligned}\ ] ] where @xmath76 and @xmath77 the last line defines a `` weighted spectral function '' , @xmath78 , which represents thermally excited states and plays an important role in the determination of the optical conductivity .
let us remark that , although eqs .
( [ sigma])-([wdos ] ) are formally analogous to those for the pure @xmath1-@xmath2 model , the electron - phonon interaction appears implicitly in them in the evaluation of the local green s function @xmath42 .
note also that , in the paramagnetic limit @xmath79 , eqs .
( [ sigma])-([wdos ] ) do not recover the results of the holstein model,@xcite because eq .
( [ sigma ] ) involves the convolution of two _ local _ rather than * k*-_dependent _ propagators . as discussed above , this is due to the classical treatment of the spin degrees of freedom which does not allow for coherent transport , so that no drude peak is recovered in the present analysis .
this , however , has only a minor influence on the finite - frequency optical conductivity , which is dominated by local incoherent excitations .    in order to assess the validity of the present treatment , we compare in fig .
[ f - feshke ] the optical conductivity of the holstein-@xmath1-@xmath2 model in infinite dimensions , as described by eqs .
( [ sigma])-([wdos ] ) , with numerical calculations using lanczos diagonalization for a single hole in the 2d holstein-@xmath1-@xmath2 model on a @xmath80 cluster.@xcite for technical reasons ( see discussion below ) , dmft data are averaged with a gaussian filter of amplitude @xmath81 such that phonon resonances are still well separated .
microscopic values are @xmath82 , @xmath62 , @xmath83 , and @xmath20 ( for ref . ) and @xmath84 for the present results .
these values correspond to a case where the lattice / magnetic polaron is formed and incoherent contributions to the optical conductivity are indeed dominant .
the good agreement of the overall shape confirms the feasibility of our analysis to investigate the finite frequency optical conductivity .      as discussed in the previous section ,
the one - hole spectral function consists of a set of narrow bands whose width @xmath51 is controlled by the intrinsic magnetic broadening .
this intrinsic width vanishes at low temperature due to the absence of spin - wave dispersion , making the direct evaluation of eqs .
( [ sigma])-([wdos ] ) quite hard to perform .
special care is thus needed in order to calculate numerically the optical conductivity .
we shall be mainly interested in the features of the optical conductivity that are related to the electron - lattice coupling .
we shall therefore retain the details of the spectra on the scale of the phonon frequency @xmath56 . from the practical point of view
, the primary object of our analysis will be a gaussian average of the optical conductivity @xmath0 with a gaussian filter of amplitude @xmath81.@xcite note that , even though the gaussian average preserves the total spectral weight of the original set of data , the accuracy of the final result will be limited by the finite frequency sampling @xmath85 . in explicit terms , if the spacing @xmath85 is larger than the intrinsic peak - width @xmath51 , the data sampling will probe the spectral features in a random way , yielding a highly inaccurate result for both the shape and spectral weight of the optical conductivity .
this is shown in the left panels of fig .
[ filtro ] where we plot the dependence of the ( gaussian averaged ) optical conductivity on the sampling spacing @xmath85 for different temperatures .
at high temperature @xmath86 the intrinsic magnetic broadening @xmath51 is large enough so that both the shape and the total spectral weight of the optical conductivity are well captured even with a relatively large mesh ( @xmath87 ) . at lower temperature @xmath68
however @xmath51 becomes so small that a much finer mesh ( @xmath88 ) is needed in order to get accurate results . at @xmath89 ,
finally , no convergence is achieved even for the finest sampling mesh considered in this paper , @xmath90 ( in fig .
[ filtro ] , for graphical reasons , we plot curves only up to @xmath88 ) .
clearly , since the intrinsic peak - width @xmath51 vanishes exponentially at low @xmath91 , the problem evidenced here can not be solved by merely reducing the sampling mesh .    to overcome this difficulty we add to the system a small uncorrelated disorder with semielliptic distribution of amplitude @xmath92,@xcite which is able to yield a finite peak - width even in the zero temperature limit .
we choose @xmath93 to assure a sufficiently dense mesh for an accurate sampling .
we then scale @xmath94 keeping fixed the ratio @xmath95 in order to approach the correct physical limit in the absence of disorder .
results obtained with a finite @xmath92 are shown in the right panels of fig .
[ filtro ] .
no appreciable difference is visible at high temperature @xmath86 where the thermally driven magnetic broadening @xmath51 is large enough to guarantee the convergence even in the absence of disorder . on the other hand
our procedure provides a clear improvement already at @xmath68 where the convergence as function of the sampling spacing is more easily achieved in the presence of disorder ( note that the converged results coincide with the results for the most dense mesh in the _ absence _ of disorder , showing that no spurious structures appear due to our scaling procedure ) .
finally , for @xmath89 the disorder scaling procedure is the only way to guarantee convergence of both the shape and the total spectral weight of the optical conductivity .      using the above described procedure
, we shall now concentrate on the evolution of the optical conductivity in the adiabatic regime , fixing the ratio @xmath61 .
this value is qualitatively representative of the cuprates , where the half - bandwidth @xmath96 ev and typical optical phonon frequencies @xmath97 mev .
this regime is also the most interesting one from the theoretical point of view , since in this case the interplay between lattice and spin degrees of freedom has its most dramatic effects .
[ fig : scanl ] shows the evolution of the optical absorption at fixed @xmath62 and low temperature @xmath98 , upon varying the electron - lattice coupling . at @xmath99 ,
the result is very reminiscent of the spectra calculated for the pure @xmath1-@xmath2 model in ref . .
it consists of a series of magnetic peaks , dominated by the sharp single - magnon peak located at @xmath100 and rapidly decaying at higher frequency . in this weak - coupling regime @xmath99 ,
the electron - phonon interaction simply gives rise to a multi - phonon fine - structure with a gaussian profile .
each magnon peak acquires thus a phonon - driven width without modifying however the overall distribution of spectral weight .    the main effect of increasing the electron - lattice coupling is a progressive shift of the spectral weight towards higher frequencies .
this is an evidence of the formation of the lattice polaron , which occurs through a gradual crossover in the presence of a finite @xmath56 .
increasing @xmath3 also modifies the shape of the low - energy absorption edge , converting the sharp magnon - peak at @xmath101 into a smoother gaussian lineshape , typical of polaronic absorption . in the strong coupling regime ( @xmath102 ) ,
characteristic of a small lattice polaron , the position of the maximum in the optical conductivity is expected to scale linearly as @xmath103 .
this is in good qualitative agreement with our data reported in fig .
[ fig : scanl ] where however the effects of a finite hopping integral @xmath1 ( ref . ) and of the @xmath7 behavior at high frequency ( see ref . )
result in a slight redshift of the maximum of the polaronic structure .
the evolution of the optical conductivity with @xmath3 can be understood by analyzing the building blocks of eqs .
( [ rozz],[wdos ] ) . in fig .
[ fig : dostot](a ) , we report both the spectral function @xmath104 and the weighted spectral function @xmath105 for two typical values @xmath99 and @xmath60 . in the first case ( dark red curve ) , there is no energy separation between the single - hole excitations in @xmath104 and the thermally excited states in @xmath105 .
the low - frequency gap in the optical absorption seen in fig .
[ fig : scanl ] arises due to the explicit shift of the spectral function by the quantity @xmath106 in eq .
( [ sigma ] ) , representing the energy cost to create one spin - defect as the hole hops in the af background .
it is now interesting to compare these features with the results for @xmath60 [ light green curve in fig .
[ fig : dostot](a ) ]
. as we can see , the spectral functions @xmath107 for @xmath60 and @xmath99 are qualitatively similar , the only major difference being the increased number of phonon satellites involved in each magnetic peak , reflecting the increased number of phonons in the polaron cloud . in order to understand the modification of the optical conductivity we thus focus on the weighted spectral function @xmath105 .
the latter undergoes a much more drastic change reflecting in an explicit way the signature of polaron formation .
of particular relevance is the shift of @xmath105 to much more negative energies which characterizes the lattice trapping .
note also that the magnetic peaks that are clearly visible at @xmath99 are completely washed out at @xmath60 , merging into a broad polaronic - like spectrum centered at higher `` binding '' energies .
this change in the nature of the thermally excited states is at the origin of both the opening of a polaronic pseudogap and of the smoothening of the features observed in @xmath0 .
a qualitatively similar evolution is observed upon varying @xmath2 , which gives a clear illustration of the positive interplay between magnetic and lattice polaron effects .
this is shown in fig .
[ fig : scanj ] where we report the optical conductivity for different values of @xmath2 at constant @xmath60 and at the same temperature @xmath108 as in fig .
[ fig : scanl ] ( note that this corresponds to different @xmath109 as @xmath53 scales with @xmath2 ) .
remarkably , even though @xmath3 is kept constant , reducing the magnetic exchange from the initial value @xmath110 leads to a gradual loss of the lattice polaronic features and to the undressing of the hole from multi - phonon excitations .
this results in a shift of spectral weight from the polaronic peak to a magnon peak at lower frequency .
the magnon peak is already visible as a shoulder in the data of fig .
[ fig : scanj ] at @xmath111 and @xmath112 , and it emerges more clearly at lower values of @xmath2 . at @xmath113 , in particular ,
both the magnetic peak and the broad polaronic band are visible in the absorption spectrum .
the featureless nature of the absorption curves at low @xmath49 can be ascribed to the increasing disorder of the af environment as @xmath2 diminishes and @xmath114 increases.@xcite    once again , the comparison of @xmath105 and @xmath104 for @xmath62 [ light green curve in [ fig : dostot](a ) ] and for @xmath115 [ dark blue curve in [ fig : dostot](b ) ] allows to visualize in a simple way the loss of the polaronic features . even in this case
the most relevant quantity is the weighted spectral function @xmath105 whose main excitations , for @xmath116 , are shifted to much less negative energies than for @xmath62 , closing the gap between @xmath105 and @xmath104 .
note in addition that , although less evident than at @xmath62 , a magnetic peak in the spectral function @xmath104 is still visible even for @xmath115 .
the convolution of @xmath104 with @xmath105 gives rise to the small magnetic peak at @xmath117 observed in the optical conductivity .    in fig .
[ fig : scant ] we report the evolution of the optical spectra versus temperature in the polaronic regime at @xmath60 , @xmath61 and @xmath118 . although fig .
[ fig : scant ] can look at a first glance quite similar to fig .
[ fig : scanl ] , no shift of the peak of @xmath0 is observed here . instead , there is a progressive _ filling _ of the low frequency gap with temperature that is quite similar to what is observed in the pure holstein model.@xcite this can be pointed out again by the comparative study of @xmath104 with @xmath105 . for @xmath71 [ light brown curve in [ fig : dostot](b ) ] the weighted spectral function @xmath78 still presents an extended peak at negative energy @xmath119 , signalizing that the lattice polaron is not completely destroyed .
however , the broadening of @xmath78 is now significantly enhanced as @xmath120 . this feature , along with a similar broadening of the spectral function @xmath121 due to the approaching of the paramagnetic limit @xmath122 ( see sec .
[ s - spectral ] ) , leads to a significant overlap of the two spectral functions and to a continuous filling of the gap in the optical conductivity as the temperature is increased .      in the previous appendix
we have derived an exact expression for the green s function of a single hole in the holstein-@xmath1-@xmath2 model in infinite dimensions .
here we investigate within the same framework the optical conductivity per hole , @xmath0 , in the zero density limit .
to this aim we provide an alternative derivation with respect to ref . .
the main advantage of the present approach is to deal at the same level with both the spin and phonon degrees of freedom , allowing thus for an immediate generalization of the @xmath1-@xmath2 model to the holstein-@xmath1-@xmath2 model . as a result
we obtain a final expression of the optical conductivity as a functional of the one - hole green s function which is formally similar to the one of ref .
but where the local green s function takes into account the electron - phonon interaction .    the formal way to derive the optical conductivity @xmath0 for a single charge , is to consider the limit @xmath203 , where @xmath204 , @xmath205 are quantities defined in the grand - canonical ensemble @xmath206 , and where @xmath207 it the total number of charges .
the limit @xmath208 is enforced by expanding to lowest order in terms of the fugacity @xmath209 ( @xmath73 ) . in this limit only the subsector @xmath210 survives both in @xmath204 and in @xmath205 .
let us first consider the hole density , which we can write as : @xmath211 where @xmath212 is a complete set of eigenstates with eigenvalues @xmath213 of one single hole ( subspace @xmath210 ) . here
, since the state @xmath212 must contain one hole at site @xmath8 , the hole - number operator is simply defined as @xmath214 without any spin defect .
inserting now into eq .
( [ ndef ] ) a complete set of eigenstates @xmath215 in the subspace @xmath216 ( no hole ) and a @xmath48-function , we obtain @xmath217 .
\label{n1 } \end{aligned}\ ] ] note that only the @xmath215 which do not contain any spin defect at site @xmath8 would contribute to eq .
( [ n1 ] ) , otherwise , after the hole @xmath13 creation , we would end up with a state with one hole and one spin defect present , which is forbidden in the hilbert space .
we have thus : @xmath218 ,   \label{n2 } \end{aligned}\ ] ] where @xmath219 is the statistical probability to have a site with no spin defect and where now only the @xmath215 with no spin defect at the site @xmath8 are selected . the square bracket in eq .
( [ n2 ] ) is just the local spectral function @xmath220 , namely @xmath221 [ remind that @xmath222 , so that we have @xmath223    let us turn now to the optical conductivity or , more precisely , to the current - current response function @xmath224 which is related to @xmath204 through the relation @xmath225 .
according to the previous argumentation , in the limit @xmath226 we can limit our analysis to the @xmath210 subspace and write : @xmath227 where @xmath228 is the imaginary time in the matsubara space and @xmath229 is the current operator to be defined below .
the current density operator can be written as @xmath233 , where we are summing explicitly on all possible directions ( this prescription compensates for the well - known vanishing of the current - current response in infinite dimensions@xcite ) . as discussed in refs . and as appearing in eq .
( [ ham ] ) , @xmath234 if there is a spin defect on the site @xmath180 ( @xmath235 ) , while @xmath236 if there is _ no _ spin defect on the site @xmath180 ( @xmath237 ) . due to the classical nature of the magnetic background , it is easy to realize that a retraceable path approximation is enforced also in the current - current response function just as in the one - particle green s function .
we obtain thus : @xmath238        let us summarize the physical meaning of eq .
( [ jj5 ] ) .
the eigenstate @xmath241 contains one hole at the site @xmath180 , and eq .
( [ jj5 ] ) describes the hopping of the hole to site @xmath8 .
if no spin defect is initially present at @xmath8 ( @xmath145 ) , this process involves the creation of a spin defect at @xmath180 , and hence @xmath149 , @xmath242 . on the other hand , in the alternative case where a spin defect is initially present at @xmath8 ( @xmath146 ) ,
the hopping process destroys the spin defect at the site @xmath8 , so that @xmath243 and @xmath244 .
let us now ask ourselves the following question : how much does the one - hole state @xmath241 differ from the `` free - like '' state @xmath246 in the absence of holes ? it is clear that in the @xmath182 case only the phonon - spin configuration on the site @xmath180 is affected by the presence of the hole , whereas all the other sites would be unaffected . in the presence of hole dynamics ( @xmath247 ) , however , all the other sites are in principle affected .
we remind also that , because of the classical magnetic background , the hole dynamics obeys a retraceable path approximation just as in a bethe lattice , as depicted in fig .
[ f - sketch ] .
let us consider now a given specific link of two nearest neighbor sites , @xmath248 .
in a true bethe lattice , the whole system can be divided in two subspaces , ( 1 ) and ( 2 ) , connected by the hopping term @xmath249 .
this means that , to leading order , the presence of one hole on the site @xmath180 would not affect the subspace ( 1 ) .
similar considerations hold true for a generic lattice under the retraceable path conditions : additional links between the two subspaces ( dashed lines in fig .
[ f - sketch ] ) , will contribute only to @xmath250 and they can be neglected in the @xmath251 limit . for practical purposes we can thus split the total hamiltonian as @xmath252 , where @xmath253 accounts for the phonon - spin degrees of freedom of the subspace ( 1 ) ( not including @xmath180 ) , while @xmath254 contains the phonon - spin degrees of freedom of the subspace ( 2 ) ( not including @xmath8 ) . in a similar way the eigenstates @xmath241 ,
@xmath246 can be written ( in the leading order of a @xmath255 expansion ) as @xmath256 and @xmath257 . employing these results , and noting that @xmath258 , @xmath259 , @xmath260 , we have thus :      where @xmath262 , @xmath263 are the partition functions of the corresponding subspaces and @xmath264 denotes a sum over the @xmath265 nearest neighbors of the site @xmath8 .
performing the analytical continuation @xmath266 , and introducing once more appropriate @xmath48-functions , we end up thus with : @xmath267 }   { 4 z \omega }   \int d\omega   \mbox{e}^{-\beta \omega }   \left\ {   \sum_{\langle j \rangle_i,\alpha_1,\gamma_1 }   \frac{\mbox{e}^{-\beta e_{\gamma_1}}}{z^{\rm t}_1 }   \left|   \left\langle \alpha_1 \left| h_j^\dagger a_j \right| \gamma_1 \right\rangle   \right|^2   \delta(\omega - e_{\alpha_1}+e_{\gamma_1 } )   \right\ }   \nonumber\\   & & \times   \left\ {   \sum_{i,\alpha_2,\lambda_2 }   \frac{\mbox{e}^{-\beta e_{\alpha_2}}}{z^{\rm t}_1 }   \left|\left\langle \alpha_2 \left| h_i   \right|\lambda_2\right\rangle\right|^2   \delta(\omega+\omega - e_{\lambda_2}+e_{\alpha_2 } )   \right\ }   \nonumber\\   & = &   \frac{t^2\pi n \mbox{e}^{\beta \mu } \left[1-\mbox{e}^{-\beta\omega}\right ] }   { 4\omega }   p_0 p_1\int d\omega   \mbox{e}^{-\beta \omega }   \bar{\rho}_1(\omega)\bar{\rho}_0(\omega+\omega ) ,   \label{pi0 } \end{aligned}\ ] ]      the same derivation can be now employed for the case when a spin defect is present on the site @xmath8 .
after few straightforward calculations , we obtain : @xmath270 }   { 4\omega }   p_0 p_1   \nonumber\\   & &   \times   \int d\omega   \mbox{e}^{-\beta \omega }   \bar{\rho}_0(\omega)\bar{\rho}_1(\omega+\omega ) .
\label{pi1 } \end{aligned}\ ] ] summing the two contributions ( [ pi0 ] ) , ( [ pi1 ] ) , and reminding @xmath271 , @xmath272 , after a change of variable we obtain : @xmath273p_0 }   { 4\omega }   \int d\omega   \mbox{e}^{-\beta \omega}\bar{\rho}(\omega )   \nonumber\\   & &   \times   % \int dz   % \mbox{e}^{-\beta z}\bar{\rho}(z )   \left[p_1\bar{\rho}(\omega+\omega+\tilde{j}/2 )   + p_0 \bar{\rho}(\omega+\omega-\tilde{j}/2)\right ] ,   \label{pit } \end{aligned}\ ] ] where we made use also of the relation @xmath274 .
finally , dividing eq .
( [ pit ] ) by ( [ nfin ] ) we obtain the dimensionless quantity eq .
( [ sigma ] ) .
note that , dividing eq .
( [ pit ] ) by ( [ nfin ] ) , the common factors @xmath275 cancel out , so that the limit @xmath276 is well defined .
shen , f. ronning , d.h .
lee , n.j.c .
ingle , w. meevasana , f. baumberger , a. damascelli , n.p .
armitage , l.l .
miller , y. kohsaka , m. azuma , m. takano , h. takagi , and z .- x .
shen , phys .
lett . * 93 * , 267002 ( 2004 ) .                                for a review of early works
see e. dagotto , rev .
phys . * 66 * , 763 ( 1994 ) .
see also later papers as : j. jakli and p. prelovek , phys .
b * 52 * , 6903 ( 1995 ) ; m.m .
zemlji and p. prelovek , phys .
b * 72 * , 075108 ( 20055 ) .
in the previous section we have discussed the technical difficulties involved in the calculation of the optical conductivity . even adopting the proposed scaling procedure ( cf .
[ tech ] ) , the computational cost can still be quite demanding , especially in the strong coupling lattice polaron regime where the number of phonons to be taken into account scales as @xmath123 . in this section
we present a simple analytical formulation of the optical conductivity which is valid precisely in the adiabatic ( @xmath124 ) small polaron regime and which involves almost no numerical effort .
the basis of this simplified formulation is the polaronic nature of the weighted spectral function @xmath105 pointed out in fig .
[ fig : dostot ] .
we have seen indeed that at moderate values of the electron - lattice coupling , just above the polaron crossover , the function @xmath105 becomes essentially independent of the exchange term @xmath2 .
a closer look shows that in the lattice polaron regime the spectral function @xmath105 can be well approximated with the result of the atomic limit @xmath125,@xcite @xmath126,\ ] ] where @xmath127 .
this function for the case @xmath60 , @xmath62 , @xmath61 , @xmath89 is shown in fig .
[ fig : dostot]a compared with the full numerical solution for the same parameters .
the good agreement can be ascribed to the atomic nature of the thermally induced excitations , which do not hybridize with the hopping continuum described by the self - energy eq .
( [ sigma_t ] ) ( an equivalent result was already pointed out in ref . ) .
( [ rhowsc ] ) provides thus a simple analytical expression for @xmath78 which can be employed in eq .
( [ sigma ] ) for the evaluation of the optical conductivity .
a simple analytical approximation can also be derived for the one - hole spectral function @xmath107 involved in eq .
( [ sigma ] ) . at this level , we are mainly interested in the overall shape of the optical conductivity @xmath0 , disregarding the multi - phonon fine structure that will be anyway washed out at large @xmath3 as shown in fig . [
fig : scanl ] . in this perspective
we can consider formally the limit @xmath128 and use the approach devised in refs . as well as in ref .
to treat the adiabatic limit .
in particular , the sum over the discrete energy levels in eq .
( [ g - gnn ] ) can be replaced by an integration over gaussianly distributed local random energies @xmath129 which account for the thermal / quantum fluctuations of the phonon field .
the local propagator thus becomes@xcite @xmath130 where @xmath131 is a gaussian function with the same variance as in eq .
( [ rhowsc ] ) .
self - consistency is achieved by employing eqs .
( [ bathg])-([sigma_t ] ) . the comparison between the approximate formula for @xmath107 valid in the polaronic regime and the full numerical solution is also shown in fig .
[ fig : dostot]a for @xmath60 , @xmath62 , @xmath61 , @xmath89 .
note that , although the bath propagator @xmath132 does not depend explicitly on the local random energies @xmath129 , it still depends parametrically on @xmath3 and @xmath2 through the self - consistency conditions eqs .
( [ bathg])-([sigma_t ] ) , so that it still retains the relevant lattice and magnetic spectral structures .
once the approximate analytical expressions for the spectral functions @xmath104 and @xmath105 are provided , we can now simply evaluate the optical conductivity by using eq .
( [ sigma ] ) .
let us stress that the numerical solution of the approximate expressions for @xmath104 and @xmath105 , defined in eqs .
( [ rhowsc])-([eq : avprop ] ) requires a remarkably smaller computational cost than the full numerical solution , especially in the lattice polaronic regime where eqs .
( [ rhowsc])-([eq : avprop ] ) are valid and where the computational cost of the full numerical solution is highest .
the comparison of the approximate formula with the full numerical solution as function of the electron - phonon coupling constant @xmath3 is illustrated in fig .
[ fig : sc ] .
the agreement is excellent in the small polaron regime @xmath133 and it is still quite satisfactory even for moderate electron - phonon coupling @xmath134 .
in this work , we have provided an analytical treatment for the optical conductivity @xmath0 of one hole in the holstein-@xmath1-@xmath2 model at finite temperature , in the limit of infinite dimensions . in this context
a dynamical mean - field solution can be derived where the local self - consistent problem is solved exactly . due to intrinsic limitations enforced in infinite dimensions , we are not able to account for the coherent propagation of a single hole in the antiferromagnetic background .
nevertheless , we have shown how the incoherent optical processes , related to local excitations ( einstein phonons and local spin - flips ) are well described in our approach , as confirmed by the comparison with numerical results obtained by lanczos diagonalization .
our main aim has been to investigate the incoherent features of @xmath0 in an intermediate coupling region where the positive interplay between the magnetic and lattice degrees of freedom is more relevant , sustaining the formation of a spin - lattice polaron . in this context
we have studied the evolution of polaronic features in the optical conductivity as function of the different microscopic parameters , as the electron - phonon coupling @xmath3 , the temperature @xmath91 and the effective exchange energy @xmath135 .
we remind that @xmath135 does not represent the bare exchange energy but rather a mean - field - like weiss exchange coupling which depends on the local magnetization @xmath136 . in the cuprates , for example ,
this parameter can be tuned by varying the hole doping starting from the parent af phase .
we have shown that the role of the electron - phonon coupling is twofold .
on one hand it rules the formation of the lattice polaron , changing the incoherent part of @xmath0 from a typical spin - polaron spectrum , characterized by magnetic peaks and by an optical gap at @xmath137 , to a broad lattice - polaron - like shape located at higher frequency , from which the magnetic peaks are essentially washed out . on the other hand
it also tunes the amount of quantum lattice fluctuations , reflected in the emergence of multi - phonon satellite peaks .
this fine structure survives also when the lattice polaron is destroyed at small @xmath3 and it gives rise to an effective broadening of the magnetic peaks which can be much larger than the intrinsic broadening driven by the thermal magnetic fluctuations .
the present approach allows us to distinguish between two different mechanisms leading to a suppression of the polaronic pseudogap : @xmath8 ) at intermediate values of the electron - phonon coupling , the reduction of the effective exchange energy @xmath135 leads to a shift of spectral weight from high energy lattice polaronic features to low energy magnetic excitations .
this results in a _ closing _ of the pseudogap as the hole undresses from its lattice polaron cloud ; @xmath9 ) conversely , increasing the temperature within the polaronic regime , gives rise to a _ filling _ of the pseudogap .
both mechanisms are actually observed in the optical spectra of the underdoped cuprates.@xcite    as a final remark , we briefly comment on the robustness of our results in physical systems where quantum spin fluctuations are present , allowing for coherent motion of the holes . as discussed above ,
one of the main effects is the emergence of a dispersive pole with small spectral weight in the one - particle spectral function . in the optical spectra , this gives rise to a drude - like low - frequency response , but it does not affect the high - frequency incoherent part ( moreover the coherent spectral weight is strongly reduced when the electron - lattice coupling is turned on ) .
regarding the high - frequency part , the intrinsic dispersion of the spin fluctuations is itself expected to smear the magnetic peaks .
this would affect our results only in the weak electron - phonon coupling regime , where no lattice polaron is formed .
in this case an intrinsic broadening of the magnetic peaks in the optical spectra due to the spin - fluctuation dispersion should be considered for a quantitative analysis . on the other hand ,
the effects of the dispersion of spin - fluctuations are expected to be barely visible in the lattice polaron regime , where the smearing due to the multi - phonon satellite structure around each magnetic peak dominates .
in this appendix we provide a detailed derivation of the green s function of one hole in the holstein-@xmath1-@xmath2 model at finite temperature in the infinite dimensional limit .
a formal derivation for the pure @xmath1-@xmath2 model was discussed in ref . , while the derivation for the full holstein-@xmath1-@xmath2 model at @xmath20 was provided in ref . .
on this ground , here we limit ourselves to the derivation of the finite temperature self - consistent equations in terms of continued fractions .
let us start by writing the hamiltonian in eq .
( [ ham ] ) as @xmath138 , where @xmath139 represents the non - local hopping term while @xmath140 contains all the other , purely local , contributions .
we also define the green s function as :      where @xmath142 denotes the state with @xmath143 phonons and @xmath144 spin defects ( @xmath145 : no spin defect , @xmath146 : spin defect ) on the site @xmath8 and with a generic set @xmath147 of phonons and @xmath148 of spin defects on all the other sites . with these notations @xmath149
when @xmath145 and @xmath150 when @xmath146 .
in addition , @xmath151 is the total partition function of the system in the absence of holes , @xmath152 being the number of sites ; in this case @xmath139 does not contribute and @xmath151 factorizes into a phonon and a spin part as @xmath153 .
each of them can be in addition factorized with respect to the site index , e.g. @xmath154 .
similar considerations hold true for the exponential terms @xmath155 which apply on the states with no holes on the right side of eq .
( [ gdef ] ) . reminding @xmath156 ,
we can write , after few straightforward steps , in the fourier space :    @xmath157 } }   { z^{\rm ph}_i z^{\rm spin}_i }   \sum_{\{n\}_{j \neq i } , \{s\}_{j \neq i } }   \frac{\mbox{e}^{-\beta \sum_{j \neq i}[n_j \omega_0 + s_j \tilde{j}/2 ] } }   { z^{\rm ph}_{j \neq i } z^{\rm spin}_{j \neq i } }   \nonumber\\   & &   \times   \left\langle \{s\}_{j \neq i},\{n\}_{j \neq i};s_i , n_i \right|   c_i^\dagger   \frac{1}{\omega+n_i\omega_0+s_i\tilde{j}/2-h }   c_i   \left|n_i , s_i;\{n\}_{j \neq i},\{s\}_{j \neq i } \right\rangle .   \label{gdefw } \end{aligned}\ ] ]    this can be rewritten as @xmath158 where @xmath159 } }   { z^{\rm ph}_{j \neq i } z^{\rm spin}_{j \neq i } }   \nonumber\\   & &   \times   \left\langle \{s\}_{j \neq i},\{n\}_{j \neq i};s_i , n_i \right|   c_i^\dagger   \frac{1}{\omega+s_i\tilde{j}/2-h }   c_i   \left|n_i , s_i;\{n\}_{j \neq i},\{s\}_{j \neq i } \right\rangle .
\label{gnn - def } \end{aligned}\ ] ] similarly , the factor @xmath160 defines the local population of spin defects which can be evaluated within the mean - field theory enforced by the infinite dimensional limit , namely @xmath161 , @xmath162 .
we have thus : @xmath163 with @xmath164 } }   { z^{\rm ph}_{j \neq i }
z^{\rm spin}_{j \neq i } }   \left\langle \{s\}_{j \neq i},\{n\}_{j \neq i};s_i , n_i \right|   c_i^\dagger   \frac{1}{\omega+s_i\tilde{j}/2-h }   c_i   \left|n_i , s_i;\{n\}_{j \neq i},\{s\}_{j \neq i } \right\rangle . \end{aligned}\ ] ]    finally , using the definition of @xmath165 , one can see that @xmath166 , and we can write @xmath167 where @xmath168 } }   { z^{\rm ph}_{j \neq i } z^{\rm spin}_{j \neq i } }   \left\langle \{s\}_{j \neq i},\{n\}_{j \neq i};n_i \right|   \frac{1}{\omega - h }   \left|n_i;\{n\}_{j \neq i},\{s\}_{j \neq i } \right\rangle .   \label{bargnn - def } \end{aligned}\ ] ]      we can write eq . ( [ bargnn - def ] ) using the short - hand notation : @xmath172 where @xmath173 , and where the operator @xmath174 } }   { z^{\rm ph}_{j \neq i } z^{\rm spin}_{j \neq i } }   \nonumber\\   & &   \hspace{-2cm}\times   \left\langle \{s\}_{j \neq i},\{n\}_{j \neq i}\right|   \ldots   \left|\{n\}_{j \neq i},\{s\}_{j \neq i } \right\rangle , \end{aligned}\ ] ] denotes the average over the spin and phonon configurations on all the sites but @xmath8 .
the operator @xmath175 represents a direct generalization of the quantity @xmath176 introduced by stumpf and logan in ref . to include the phonon degrees of freedom .
it is also convenient to note that ( for @xmath177 ) : @xmath178 } }   { z^{\rm ph}_j z^{\rm spin}_j }   \nonumber\\   & &   \hspace{-2cm}\times   \left\langle s_j , n_j , n_i\right|   \ldots   \left|n_i , n_j , s_j\right\rangle ,   \label{pop } \end{aligned}\ ] ] where @xmath179 is defined in similar way as @xmath175 as the average over the spin and phonon configurations on all the sites except @xmath8 and @xmath180 .
having introduced the necessary definitions , from now on we can follow the derivation in ref .
properly adapted to the finite temperature case .
in particular we can introduce the local green s function @xmath181 defined as the atomic @xmath182 limit of eq .
( [ bargnn ] ) .
note that in the atomic limit @xmath183 so that @xmath184 we can also generalize eq .
( [ bargnn ] ) for off - diagonal local phonon matrix elements : @xmath185 whose analytical expression will be provided in eq .
( [ glocnm ] ) .
we can now employ the standard relation @xmath186 which , on a classical spin background , gives rise to the retraceable path approximation .
we have thus @xmath187 where @xmath188 represents the dynamics of the hole after hopping on the neighboring site @xmath180 .
since , in the leading term of a @xmath189 expansion , the hopping process and the further dynamics of the hole do not involve the phonon degrees of freedom on site @xmath8 , such evolution occurs _ in the presence _ of @xmath190 phonons on the site @xmath8 , which are reflected in an shift of the frequency argument , @xmath191 .
note that @xmath192 still contains the full thermal average on the phonons at site @xmath180 , so that it is related to the thermally averaged green s function defined in eq .
( [ g - gnn - app ] ) .
in addition , @xmath192 will depend also on the initial spin configuration @xmath193 at site @xmath180 : hopping to a site @xmath180 free of spin defects will _ create _ a spin defect on the site @xmath8 , while hopping to a site @xmath8 with a spin defect will restore the initial magnetic background at site @xmath8 _ destroying _ a spin defect at @xmath8 . using these considerations
we have thus : @xmath194   \nonumber\\   & & \hspace{-14 mm }   =   \frac{t^2}{4 }   \left [   ( 1-x)\bar{g}_j(\omega-\tilde{j}/2 )   + x
\bar{g}_j(\omega+\tilde{j}/2 )   \right ] ,   \label{sigma_t - app } \end{aligned}\ ] ] and we can write eq .
( [ recursion ] ) in the compact form : @xmath195 where we have dropped the unnecessary site indices and we have added the index @xmath196 to denote a self - energy term arising from the hopping processes .
we can indeed employ now the usual procedure to obtain the electron - phonon self - energy in terms of a continued fraction .
in particular , in addition to the diagonal elements of eq .
( [ glocnn ] ) , we specify also the non - diagonal ones which read : @xmath197^{nm}=   \left[\omega - n\omega_0\right]\delta_{n , m}+g x^{n , m}.   \label{glocnm } \end{aligned}\ ] ]