in recent years , it has become widely accepted that electron - phonon ( ep ) interaction plays a crucial role in systems with strong electronic correlations , such as high - temperature superconducting cuprates@xcite or colossal magnetoresistive manganites.@xcite a large amount of the available data on optical and transport properties has been interpreted in terms of polaronic states . in view of the complex underlying physics , however , the application of single - polaron theories seems to be problematic . on the other hand , a reliable theory of strongly coupled ep systems with finite carrier densities
is not available yet . while limiting cases , such as weak coupling or large phonon frequencies  compared to the electronic bandwidth  can to some degree
be described theoretically , especially the adiabatic regime of small phonon frequencies at intermediate - to - strong ep coupling constitutes a long - standing , open problem . since analytical results
can only be obtained at the expense of uncontrolled approximations , progress toward an understanding of many - polaron systems requires unbiased numerical methods .
owing to the availability of high performance computers , the latter can yield insight into the properties of basic microscopic models such as the spinless holstein model considered here .
nevertheless , apart from work on the half - filled band case , where the physics is dominated by a cross over from a luttinger liquid to a peierls state , the many - polaron problem has received little attention in the past .    motivated by this situation , in the present work
, we address the important issue of the effect of carrier density on the character of the quasiparticles of the system . while for very strong couplings
no significant changes are expected due to the existence of rather independent small ( self - trapped ) polarons with negligible residual interaction , a density - driven cross over from a state with large polarons to a metal with weakly dressed electrons may occur in the intermediate coupling regime .
this problem has recently been investigated experimentally in terms of optical measurements on @xmath0 films.@xcite while small polarons have been identified as the charge carriers in the low - temperature metallic phase of the ca doped system , the weaker electron - lattice coupling in the sr doped system gives rise to large - polaron - like behavior of the optical conductivity .
the effect of the finite carrier density on such extended polaronic quasiparticles has been analyzed@xcite using a recently developed weak - coupling theory.@xcite    for the case of the holstein model , the abovementioned density - driven transition from small to large polarons and finally to weakly ep - dressed electrons is expected to be possible only in one dimension ( 1d ) , where large polarons extending over more than one lattice site exist at weak coupling .
existing work indicates that in higher dimensions , the electrons remain quasifree at weak coupling , and become self - trapped above a critical value of the ep coupling strength.@xcite by contrast , approximate variational calculations suggest that large polarons may also exist in @xmath1.@xcite the situation is different for frhlich - type models@xcite with long - range ep interaction , in which large polaron states exist even for strong coupling and in higher dimensions .    here
, we make use of quantum monte carlo ( qmc ) and exact diagonalization ( ed ) methods , and calculate the photoemission spectra for finite clusters at zero and finite temperature . due to the limitations of existing algorithms ,
we develop a new grand - canonical qmc approach which is free of autocorrelations and therefore enables us to study small phonon frequencies and low temperatures .
the organization of this paper is as follows . in sec .
[ sec : model ] , we present the model and discuss previous work , while sec .  [ sec : methods ] contains details about the methods we employed .
section  [ sec : results ] gives a discussion of our numerical data , and finally , we conclude in sec .  [
sec : conclusions ] .
the new qmc technique is described in the appendix .
the spinless holstein model is defined by the hamiltonian @xmath2 where @xmath3 ( @xmath4 ) creates ( annihilates ) a spinless fermion at lattice site @xmath5 , @xmath6 ( @xmath7 ) creates ( annihilates ) a phonon of energy @xmath8 ( @xmath9 ) at site @xmath5 , and @xmath10 . while the first term describes hopping processes between neighboring lattice sites @xmath11 , the second term corresponds to the elastic and kinetic energy of the phonons . finally , the last term in eq .
( [ eq : model : h ] ) represents a local coupling of the lattice displacement to the electron occupation @xmath12
.    the parameters of the model are the hopping integral @xmath13 , the phonon energy @xmath8 , and the coupling constant @xmath14 .
the latter is related to the polaron binding energy @xmath15 via @xmath16 . by defining the dimensionless quantities @xmath17 , and @xmath18 ,
we are left with two independent parameters . in the sequel
, we express all energies in units of @xmath13 .
the lattice constant is taken to be unity , and periodic boundary conditions are applied .
the holstein hamiltonian  ( [ eq : model : h ] ) represents a generic model for a many - polaron system . due to the restriction to spinless fermions
, on - site bipolaron formation is suppressed.@xcite    quantum fluctuations of the lattice enable the electrons to hop even in the strong - coupling regime , while in models with classical phonons translation symmetry is broken , so that no itinerant polaron states exist .
along this line , the half - filled spinless holstein model has been studied in one dimension using qmc,@xcite ed,@xcite dmrg@xcite and variational methods.@xcite similar to the spinfull model at half filling , the spinless model with one electron per two lattice sites displays a quantum phase transition from a metallic state to a peierls insulating state with charge - density - wave ( cdw ) order at @xmath19 . while the insulating cdw state of the spinless model is stable for any @xmath20 in the adiabatic limit @xmath21 , it can be destroyed by quantum phonon fluctuations for @xmath22.@xcite capone _ et al._@xcite used ed to study the cross over from a single polaron to a many - electron system away from half filling . although their work is actually for the holstein - hubbard model , results are very similar to the spinless holstein model in the limit of large @xmath23 considered in ref .  .
the authors found that the critical coupling for polaron formation is unaffected by the electron density in the antiadiabatic regime .
finally , we would like to point out previous analytical efforts to take into account squeezing and correlation effects in many - polaron systems.@xcite
in this work we use qmc , and ed in combination with the kernel polynomial method ( kpm ) , to calculate , in particular , single - particle spectra .
these two approaches are characterized by different advantages and shortcomings .
ed in the form used here yields exact spectra at @xmath19 with extremely high energy resolution .
however , even for relatively small clusters of ten sites  the size used here  this method requires parallel supercomputers to handle the enormous hilbert space ( matrix dimension @xmath24@xmath25 ) .
consequently , the momentum resolution , determined by the cluster size , is rather crude .
in contrast , the grand - canonical qmc simulations can be performed on personal computers , and for larger clusters .
the drawback is that ( a ) the method is limited to finite temperatures , and ( b ) the calculation of spectral properties requires the analytic continuation from imaginary time to real frequencies , which is an ill - conditioned problem .
therefore , regularization techniques such as the maximum entropy method ( mem ) used here have to be applied , which offer only limited energy resolution and may introduce uncontrolled errors .
for this reason , we will also present a direct comparison of ed and qmc results .      in applying existing qmc methods for the holstein model to the parameter regime considered here ,
one faces extremely long autocorrelations times , as well as large statistical errors.@xcite this makes it very difficult to obtain accurate results for spectral properties , due to the use of the mem . to overcome these problems , the algorithm presented in the appendix is based on the canonical lang - firsov transformation.@xcite the latter explicitly contains the shift of the lattice equilibrium position in the presence of an electron , which is favorable for the qmc simulations .
moreover , the absence of a direct coupling between the electron density and the lattice displacements in the transformed hamiltonian  ( [ eq : app : hspinless ] ) permits exact , uncorrelated sampling of the phonon degrees of freedom ( momenta ) in terms of principal components .
here we have mainly used a value of @xmath26 for the inverse temperature @xmath27 . for
the critical parameters studied below , accurate simulations at even lower temperatures become very demanding due to a minus - sign problem ( sec .
[ sec : app : sign_results ] ) . nevertheless , we shall present results for the density of states at @xmath28 and half filling but for a reduced cluster size .
note that already for @xmath26 we have @xmath29 for the values
@xmath30 used in the sequel , so that direct thermal excitation of phonons is expected to be rather unimportant .
for the trotter decomposition , we chose @xmath31 ( see appendix ) for static observables . owing to the larger numerical effort associated with the calculation of dynamical properties , for the latter we use @xmath32the larger error being small compared to the uncertainties introduced by the mem inversion .
the trotter error in the holstein model becomes important mainly for large phonon frequencies and very strong ep interaction.@xcite since we do not consider this regime here , the abovementioned values of @xmath33 yield satisfactorily small systematic errors .
away from half filling , the electron density has to be adjusted to the values of interest by varying the chemical potential @xmath34 . due to the computational effort associated with this trial - and - error procedure ,
the actual band filling @xmath35 , say , is usually only very close to a desired value , e.g. , @xmath36 , but not exactly the same .
furthermore , within the qmc simulations , the filling can only be determined up to a statistical error @xmath37 , similar to other observables . for the results presented here ,
the relative deviation @xmath38 of the actual value in the qmc simulations from the value @xmath39 reported is always smaller than @xmath40 .
the same is true for the relative statistical error @xmath41 .    to perform the continuation of the one - electron green function from imaginary time to real frequencies
, we use the mem ( for a review see , e.g. , ref .  ) . the algorithm employed here
incorporates an improved evidence approximation,@xcite which is less susceptible to overfitting  resulting in peaks whose position , weights and even existence are not clear  than the conventional evidence approximation.@xcite on the other hand , groups of narrowly spaced peaks are represented by an envelope if the accuracy of the qmc data is not sufficient .
the proposal of an alternative evidence approximation  leading in general to stronger regularization and , consequently , lower energy resolution than the conventional approximation  may not seem advantageous at first .
however , considering the great variety of different spectra compatible with the qmc data within the error bars due to the extremely ill - conditioned analytical continuation problem , we argue that it rather corresponds to choosing the `` most uninformative '' solution in the spirit of the maximum entropy principle .
in fact , we found that the alternative evidence approximation yields much more consistent results than the `` classic '' evidence approximation@xcite when we subjected artificial test data to different noise while keeping the signal - to - noise ratio fixed .
a detailed comparison of different regularization schemes available will be presented elsewhere .    in order to scrutinize the resulting spectra , the following
test has been performed .
following the idea of silver _
et al._,@xcite the results from the mem applied to the original qmc data @xmath42 have been taken as an `` exact '' spectral function @xmath43 , from which mock data @xmath44 for the imaginary - time green function has been created by using eq .
( [ eq : app : maxent ] ) and adding gaussian random noise of absolute size @xmath45@xmath46 , comparable to the statistical errors of the qmc results . again using the mem to invert these data , the thereby obtained spectral function @xmath47 allows to check the stability of the mem result @xmath43 with respect to statistical errors .
we find that the features in our results for @xmath43 remain stable if @xmath48 varies within the statistical errors , so that the data presented below are satisfactorily reliable apart from the aforementioned broadening and the limited energy resolution .
while our method is completely free of autocorrelations between subsequent phonon configurations , results for the imaginary - time green function at different times are usually not statistically independent.@xcite however , in the present case , owing to the exact sampling and the large number of measurements , these correlations are very small and have been neglected since they have negligible influence on the results .
an additional , stringent test of the qmc results consists of a comparison to ed data which will be presented in sec .
[ sec : results : intermediate ] .
finally , we note that all qmc simulations have been run on personal computers ( e.g. , intel xeon 2600 mhz ) , with cpu times varying between several hours and a couple of days .
the combination of ed and kpm has been described in detail in ref .  .
in this work , we are mainly interested in the single - particle spectral functions @xmath49      \ , ,    \end{split }    \\    & & \begin{split }      a^-(k,{\omega } )      =       \sum_l      |      &      { { \langle}\psi_{l , k}^{(n_\text{e}-1)}| }      c^{{\phantom{\dag}}}_{q - k }      { |\psi_{0,q}^{(n_\text{e})}{\rangle } }      |^2\ ,      \\      &      \times      \delta[{\omega}+ ( e_{l , k}^{(n_\text{e}-1 ) } - e_{0,q}^{(n_\text{e } ) } ) ]      \ , ,    \end{split}\end{aligned}\ ] ] where @xmath50 denotes the @xmath51th eigenstate with @xmath52 electrons , momentum @xmath53 and energy @xmath54 . the cluster size has been fixed at @xmath55 , with as many as 15 dynamical phonons and resulting truncation errors @xmath56 for the ground - state energy.@xcite additionally , the convergence with the number of phonons has been monitored by means of the phonon distribution function.@xcite it is important to point out that such a small number of dynamical phonons only yields accurate results for the parameters considered here if the symmetric phonon mode with @xmath57 is separated from the hamiltonian analytically before the diagonalization.@xcite in the present case , this phonon mode  corresponding to a homogeneous shift of the oscillator equilibrium positions  together with that for @xmath58 represents the dominant contribution to phonon excitations , so that the separation is very efficient at reducing the computational effort .
the largest matrix dimension in the present calculations was about @xmath59 , running on six compute nodes of a hitachi sr8000-f1 .
the energy resolution of the kpm spectral function with 512 moments has been enhanced by a factor of 8 by the mem,@xcite resulting in a significantly better resolution compared to , e.g. , the spectral decoding method .
we would like to point out that in contrast to the case of qmc , the ed data to which the mem is applied does not have a statistical error .
it is long known that within the holstein model  ( [ eq : model : h ] ) , a single electron undergoes a cross over to a small polaron at a critical value of the ep coupling.@xcite the important question to be addressed here is what happens if the density of electrons ( or polarons ) is increased .
we do not expect significant changes in the strong - coupling regime due to the rather
localized " nature of the small polaron state . on the contrary , for intermediate ep interaction and @xmath60 ,
each electron is surrounded by a phonon cloud extending over several lattice sites .
now , at larger densities , a substantial overlap of the single - particle wavefunctions occurs , leading to a dissociation of the individual polarons and finally to a restructuring of the whole many - particle ground state .
furthermore , it has been shown that squeezing effects reduce the polaronic band - narrowing in the many - electron case,@xcite thereby favoring such a scenario . since within the holstein model , a large polaron state
is expected to be stable in 1d , we will only consider this case in the sequel .    in view of the situation in most polaronic materials
, we mainly consider the adiabatic regime , taking @xmath61 .
to illustrate the important differences between @xmath62 and @xmath63 , however , some results for the antiadiabatic value @xmath64 will also be shown .
as we shall see below , the spectra in the intermediate coupling region turn out to have a fairly complex structure .
it is therefore helpful to begin with the limiting cases of weak and strong ep coupling .      for weak coupling @xmath65 ,
the sign problem is not severe ( sec .
[ sec : app : sign_results ] ) and the qmc simulations can easily be performed for large lattices with @xmath66 , thereby making the dispersion of quasiparticle features well visible . as indicated in eq .
( [ eq : app : maxent ] ) , the mem inversion yields @xmath67 , where @xmath34 denotes the chemical potential .
one - electron spectral function @xmath67 from qmc for different band fillings @xmath39 , @xmath66 , @xmath26 , @xmath61 , and weak coupling @xmath65.,title="fig:",scaledwidth=23.85% ]   one - electron spectral function @xmath67 from qmc for different band fillings @xmath39 , @xmath66 , @xmath26 , @xmath61 , and weak coupling @xmath65.,title="fig:",scaledwidth=23.85% ] +   one - electron spectral function @xmath67 from qmc for different band fillings @xmath39 , @xmath66 , @xmath26 , @xmath61 , and weak coupling @xmath65.,title="fig:",scaledwidth=23.85% ]   one - electron spectral function @xmath67 from qmc for different band fillings @xmath39 , @xmath66 , @xmath26 , @xmath61 , and weak coupling @xmath65.,title="fig:",scaledwidth=23.85% ]    figure  [ fig : results : qmc_lambda0.1 ] shows the evolution of the one - electron spectral function @xmath67 with increasing electron density @xmath39 . at first sight , we see that the spectra bear a close resemblance to the free - electron case , i.e. , there is a strongly dispersive band running from @xmath68 to @xmath69 .
the latter can be attributed to weakly dressed electrons with an effective mass similar to the noninteracting case . as expected , the height ( width ) of the peaks increases ( decreases ) significantly in the vicinity of the fermi momentum @xmath70 , which is determined by the crossing of the band with the chemical potential .
one - electron density of states @xmath71 from qmc for different band fillings @xmath39 ,
@xmath66 , @xmath26 , @xmath61 and @xmath65.,scaledwidth=47.5% ]    since the peaks in fig .  [ fig : results : qmc_lambda0.1 ] are generally very sharp , it is important to check the results from the inversion with the mem . by computing the approximate number of electrons corresponding to the fillings @xmath720.4 of the 32-site system , we can roughly determine @xmath70 .
although the nonzero ep interaction will affect the band dispersion @xmath73 , the weak coupling considered here justifies the use of a rigid band approximation , as reflected in the free - particle character of the spectra in fig .
[ fig : results : qmc_lambda0.1 ] .    for @xmath72 ,
i.e. , about three electrons for @xmath66 , the fermi momentum is expected to lie between @xmath74 and @xmath75 , corresponding to the second and third curve from the top in fig .  [
fig : results : qmc_lambda0.1](a ) .
in fact , we see that the largest peaks are those for @xmath74 and @xmath76 , having almost the same height . above and below these momenta ,
the peaks are noticeably smaller . a similar analysis carried out for the other fillings in figs .
[ fig : results : qmc_lambda0.1](b)(d ) also shows a very good agreement between the location of the highest peaks and the estimated values of @xmath70 .    finally , we would like to point out that the apparent absence of any phonon signatures in fig .
[ fig : results : qmc_lambda0.1 ] is not a defect of the mem , but results form the large scale of the @xmath77-axis chosen . as a consequence , the peaks running close to the bare band
dominate the spectra and suppress any small phonon peaks present .
if we enlarge the resolution , for all densities @xmath720.4 , we observe the band flattening@xcite at large wavevectors which originates from the intersection of the approximately free - electron dispersion with the optical phonon energy occurring at @xmath78 .    to complete our discussion of the weak - coupling regime ,
we show in fig .
[ fig : results : dos_wc ] the one - electron density of states ( dos ) @xmath71 given by eq .
( [ eq : app : dosgen ] ) .
clearly , for small @xmath39 , there is a peak with large spectral weight at the fermi level . in contrast , for large @xmath39 , the tendency toward formation of a peierls ( band ) insulating state at @xmath79 suppresses the dos at the fermi level , although we are well below the critical value of @xmath80 at which the cross over to the insulating state takes place at @xmath19.@xcite the additional small features separated from @xmath34 by the bare phonon energy @xmath8 will be discussed in the intermediate coupling case below .
one - electron spectral function @xmath67 from qmc for different band fillings @xmath39 , @xmath66 , @xmath26 , @xmath61 , and strong coupling @xmath81.,title="fig:",scaledwidth=23.85% ]   one - electron spectral function @xmath67 from qmc for different band fillings @xmath39 ,
@xmath66 , @xmath26 , @xmath61 , and strong coupling @xmath81.,title="fig:",scaledwidth=23.85% ] +   one - electron spectral function @xmath67 from qmc for different band fillings @xmath39 , @xmath66 , @xmath26 , @xmath61 , and strong coupling @xmath81.,title="fig:",scaledwidth=23.85% ]   one - electron spectral function @xmath67 from qmc for different band fillings @xmath39 , @xmath66 , @xmath26 , @xmath61 , and strong coupling @xmath81.,title="fig:",scaledwidth=23.85% ]    we now turn our attention to the opposite limit of the spinless holstein model , namely the case of strong ep coupling @xmath82 .
while calculations with many other numerical methods such as , e.g. , standard grand - canonical qmc or ed , become very difficult in this regime , the qmc approach used here yields quite accurate results even for large systems . the improved statistics in the strong - coupling regime is also reflected in the dependence of the average sign on @xmath80 ( cf .
the discussion of fig .
[ fig : app : sign_beta_omega ] in the appendix ) .
as before , we study the dependence of the spectral function on the band filling , again for @xmath66 .
the results are depicted in fig .
[ fig : results : qmc_lambda2.0 ] .    in the low - density regime @xmath72 [ fig .
[ fig : results : qmc_lambda2.0](a ) ] , we expect the well - known almost flat polaron band having exponentially reduced spectral weight ( given by @xmath83 in the single - electron , strong - coupling limit ) which , nevertheless , can give rise to coherent transport at @xmath19 .
unfortunately , the qmc / mem has difficulties resolving such extremely weak signatures , especially at larger @xmath39 , where the statistical errors are more noticeable , and for temperatures such that the thermal energy is comparable to the polaron band width .
obtained by integration over the spectrum is not sensitive enough to yield the correct value of @xmath34 . ]
therefore the polaron band can hardly be resolved in figs .
[ fig : results : qmc_lambda2.0](a ) or ( b ) , and is not even detectable in figs .  [ fig : results : qmc_lambda2.0](c ) and ( d )
. however , if we consider a density @xmath84 , corresponding to a single electron on a cluster with @xmath85 , at a lower temperature @xmath28 , the qmc results do indeed give a polaron band around @xmath34 ( fig .
[ fig : results : qmc_lambda2.0_polaron ] ) . due to the sign problem
, similar calculations for @xmath86 would be very demanding .
besides , the spectrum consists of two incoherent features located above and below the chemical potential , which reflect the phonon - mediated transitions to high - energy electron states . here , the maximum of the photoemission spectra ( @xmath87 ) follows a tight - binding cosine dispersion .
the incoherent part of the spectra is broadened according to the phonon distributions .
we would like to mention that the spectrum of a single polaron at strong ep coupling has been studied analytically by alexandrov and ranninger.@xcite     dispersion of the polaron band in @xmath67 for @xmath84 , @xmath85 and @xmath28 ( drawn to a larger scale as compared to fig .
all other parameters as in fig .
[ fig : results : qmc_lambda2.0].,scaledwidth=23.85% ]    with increasing band filling [ see panels ( b ) and ( c ) ] , the chemical potential is shifted to lower energies @xmath88 , but is expected to be still located in a narrow band with little spectral weight . regardless of the aforementioned shortcomings , the qmc results provide a clear picture of the physics .
there exists a finite gap to the photoemission ( inverse photoemission ) parts of the spectrum , so that the system typifies as a polaronic metal .
we shall see below that a completely different behavior is observed at intermediate coupling .
notice that the incoherent inverse photoemission ( photoemission ) signatures are more pronounced at small ( large ) wavevectors .    finally , for @xmath36 [ fig .
[ fig : results : qmc_lambda2.0](d ) ] , the incoherent features lie rather close to the fermi level , thus being accessible by low - energy excitations .
now , the photoemission spectrum for @xmath89 becomes almost symmetric to the inverse photoemission spectrum for @xmath90 and already reveals the gapped structure which occurs at @xmath79 due to the cdw - formation accompanied by a peierls distortion .    as in the weak - coupling case
discussed above , the properties of the system also manifest itself in the dos , shown in fig .
[ fig : results : dos_sc ] . owing to the strong ep interaction ,
the spectral weight at the chemical potential is exponentially small for all fillings @xmath39 .    at half filling , the dos exhibits particle - hole symmetry , and the system can be described as a peierls insulator , consisting of a polaronic superlattice .
in contrast to the weak - coupling case , the ground state is characterized as a polaronic insulator rather than a band insulator .
moreover , here the gap at the fermi level persists even at the finite temperature used in the simulation , while for @xmath65 and @xmath79 the insulating state found at @xmath19 is destroyed by thermal fluctuations ( fig .
[ fig : results : dos_wc ] ) . note that , owing to the sign problem , the result for @xmath79 in fig .
[ fig : results : dos_sc ] is for a reduced cluster size @xmath85 .
however , in the strong - coupling regime considered here , finite - size effects are expected to be very small .
one - electron density of states @xmath71 from qmc for different band fillings @xmath39 and cluster sizes @xmath91 , @xmath26 , @xmath61 and @xmath81.,scaledwidth=47.5% ]      +          as discussed in the introduction , a cross over from a polaronic state to a system with weakly dressed electrons can be expected in the intermediate coupling regime . here
we choose @xmath92 , which corresponds to the critical value for the small - polaron cross over in the one - electron problem ( cf . , e.g. , fig .  5 in ref .  ) .
we first discuss the qmc results .
owing to the sign problem , which is particularly noticeable for @xmath92 ( see appendix ) , we have to decrease the system size as we increase the electron density @xmath39 .
figure  [ fig : results : qmc_lambda1.0 ] shows the spectral function for @xmath92 and increasing band filling . owing to the overlap of large polarons in the intermediate coupling regime , we start with the very low density case @xmath84 [ fig .  [
fig : results : qmc_lambda1.0](a ) ] . compared to the behavior for @xmath82 [ fig .
[ fig : results : qmc_lambda2.0](a ) ] , we notice that the polaron band now lies much closer to the incoherent features , and that there is a mixing of these two parts of the spectrum at small values of @xmath53
. nevertheless , the almost flat polaron band is well visible for large @xmath53 .    with increasing density , the polaron band merges with the incoherent peaks at higher energies .
this is the above - anticipated density - driven cross over from a polaronic to a ( diffusive ) metallic state , with the main band crossing the fermi level . at @xmath36
, we find a band with large effective width ranging from @xmath93 to @xmath94 .
a more detailed discussion of the spectral function will be given below when we present ed results .
( color online ) one - electron density of states @xmath71 from qmc for different band fillings @xmath39 , cluster sizes @xmath91 and inverse temperatures @xmath95 . here
@xmath61 and @xmath96.,scaledwidth=47.5% ]    further information about the density dependence can be obtained from the one - electron dos .
the latter is presented in fig .
[ fig : results : dos_ic ] for different fillings @xmath840.5 . as in fig .  [ fig : results : qmc_lambda1.0 ] , the cluster size is reduced with increasing @xmath39 in order to cope with the sign problem . to illustrate the rather small influence of finite - size effects , fig .
[ fig : results : dos_ic ] also contains results for @xmath55 .    for low density @xmath84 , the dos in fig .
[ fig : results : dos_ic ] lies in between the results for weak and strong coupling discussed above .
although the spectral weight at the chemical potential is strongly reduced compared to @xmath65 , @xmath97 is still significantly larger than for @xmath81 .    when the density is increased to @xmath98 , the dos at the chemical potential increases , as a result of the dissociation of the polarons .
as we increase @xmath39 further , a pseudogap begins to form at @xmath34 , which is a precursor of the cdw gap at half filling and zero temperature .
finally , in the case of half filling @xmath79 , the dos has become symmetric with respect to @xmath34 .
there are broad features located either side of the chemical potential , which take on maxima close to @xmath99 .
however , apart from the strong - coupling case , where the single - polaron binding energy is still a relevant energy scale , the position of these peaks is rather determined by the energy of the upper and lower bands , split by the formation of a peierls state .
the gap of size @xmath100 expected for the insulating charge - ordered state at @xmath19 is partially filled in due to the finite temperature considered here .
furthermore , we find additional , much smaller features roughly separated from @xmath34 by the bare phonon frequency @xmath8 , whose height decreases with decreasing temperature , as revealed by the results for @xmath28 ( fig .  [
fig : results : dos_ic ] ) .
these peaks  not present at @xmath19 ( see ref .  )
arise from thermally activated transitions to states with additional phonons excited , and are also visible in figs .
[ fig : results : dos_wc ] and  [ fig : results : dos_sc ] . while for weak coupling [ @xmath101 , @xmath65 , fig .
[ fig : results : dos_wc ] ] , the maximum of these features is almost exactly located at @xmath102 , it moves to @xmath103 for intermediate coupling [ @xmath104 , @xmath92 , fig .
[ fig : results : dos_ic ] ] , and finally to @xmath105 for strong coupling [ @xmath106 , @xmath82 , fig .
[ fig : results : dos_sc ] ] .
although the exact positions of the peaks are subject to uncertainties due to the mem , this reflects the shift of the maximum in the phonon distribution function with increasing coupling .
the mem yields an envelope of the multiple peaks separated by @xmath8 .
+        at this point , we would like to mention the exact relation@xcite @xmath107 for the first moment of the one - particle spectral function of the spinless holstein model , with @xmath108 . while the zeroth moment is identical to the normalization of @xmath43 for each @xmath53 , @xmath109 depends in a nontrivial way on the parameters of the system . in principle , eq .
( [ eq : results : moment ] ) may be used as a boundary condition in the mem.@xcite however , in the present case , we have found that while the normalization is virtually exact for all results shown , the first moment deviates from the exact values .
moreover , the use of the exact values for @xmath110 as a condition in the inversion causes the mem not to converge for some parameters .
additional calculations on small systems reveal that this disagreement results from the finite trotter error , despite the relatively small value @xmath32 used . away from parameters
such that @xmath111 , @xmath110 as determined from the mem spectra deviates from the exact values by maximally 1020 percent .
nevertheless , the peak positions and weights in the spectra shown are fairly accurate , as illustrated by the comparison with ed below .
finally , we would like to point out that calculations of @xmath43 for @xmath112 become very time - consuming due to the increasing requirement in disk space ( see appendix ) and longer simulation times .
despite the wealth of information contained in the qmc results presented so far , it is obvious that an interpretation of the various excitations is far from being easy .
in addition to the finite temperature , the qmc method used here does not allow to separately calculate the photoemission and inverse photoemission parts @xmath113 and @xmath114 , respectively .
moreover , the use of the mem limits the energy resolution and introduces uncertainties concerning the positions and weights of structures in the spectra .
these circumstances make it difficult , e.g. , to distinguish between phononic and electronic contributions , especially if @xmath115 .
to gain further insight , we therefore supplement the qmc data with ed results .      in order to overcome the abovementioned limitations of qmc , here
we use ed in combination with the kpm , as described in sec .
[ sec : methods : ed ] , to calculate the one - particle spectral functions defined by eq .
( [ eq : methods : apm ] ) . to test the quality of the qmc spectra
, we directly compare the two methods , keeping in mind that the qmc calculations have been performed at a finite temperature , whereas ed yields ground - state results .
( color online ) comparison of the integrated spectral weight obtained from the qmc and ed data of fig .
[ fig : results : ed_qmc_lambda1.0](a ) for three different momenta @xmath53.,scaledwidth=45.0% ]    figure  [ fig : results : ed_qmc_lambda1.0 ] displays the ed results for the same parameters as in fig .
[ fig : results : qmc_lambda1.0 ] , i.e. , @xmath61 , @xmath92 , and @xmath720.4 . additionally shown
are the integrated spectral weights of @xmath113 and @xmath114 ( dashed lines ) , as well as qmc results for the same system size and @xmath116 ( solid thick lines ) , which have been shifted by @xmath34 .
the fermi level lies between the highest - energy peak in @xmath113 and the lowest - energy peak in @xmath114 .
we see from fig .
[ fig : results : ed_qmc_lambda1.0 ] that the qmc results are compatible with the ed data . because of the finite temperature and the use of the mem , qmc , of course , can not reproduce the sharp ed peaks . to further test the agreement between the two methods , we compare in fig .
[ fig : results : ed_qmc_weight ] the integrated spectral weight of the data shown in fig .
[ fig : results : ed_qmc_lambda1.0](a ) , and find a very satisfactory accord .
moreover , this agreement becomes even better for larger densities @xmath39 , for which the sharp peaks in the ed data evolve into bands [ see fig .  [
fig : results : ed_qmc_lambda1.0](b)(d ) ] .
in addition to the possibility of distinguishing between @xmath114 and @xmath113crucial for an identification of polaron bands at intermediate coupling  the higher energy resolution of ed also allows to resolve signatures of phonon excitations .    starting with the case @xmath72 [ fig .
[ fig : results : ed_qmc_lambda1.0](a ) ] , we notice that there is a polaron band at the fermi level , excitations to / from which are given by the highest ( lowest ) peak of @xmath113 ( @xmath114 ) .
the photoemission part @xmath113 reflects the poisson - like phonon distribution of the one - electron ground state ( @xmath72 and @xmath55 ) .
the integrated spectral weight gives a measure for the electronic weight of the various poles in the one - electron spectrum .
for example , it reveals that the phonon peaks in @xmath113 have very little weight .
in fact , the integrated weight jumps to a finite value at the first peaks near the fermi level for @xmath117 and @xmath118 , while it changes very little as one moves further down in energy . for @xmath119 ,
the small spectral weight contained in @xmath113 is continuously distributed among the phonon peaks .
we also observe the well - known flattening of the polaron band at large values of @xmath53 .
similar to the one - electron case,@xcite the low - energy states have mostly electronic character at small @xmath53 , and become mostly phononic at large values of @xmath53 .
while this effect is expected to occur in the low - density regime , we find that is persists even for @xmath36 [ fig .
[ fig : results : ed_qmc_lambda1.0](d ) ] . moreover , fig .  [
fig : results : ed_qmc_lambda1.0 ] reveals that the maximum in the incoherent contribution follows closely the free - electron dispersion for all densities @xmath720.4 .    with increasing electron density , the equally spaced peaks in @xmath113 broaden significantly , until at @xmath36 they have merged to form a broad band .
the polaron band is still visible at @xmath98 [ fig .
[ fig : results : ed_qmc_lambda1.0](b ) ] , but becomes indistinguishable from the incoherent excitations at even larger @xmath39 . eventually , at @xmath36 [ fig .
[ fig : results : ed_qmc_lambda1.0](d ) ] , we find an almost symmetric behavior of @xmath114 and @xmath113 with respect to the fermi level at @xmath120 . compared to the strong - coupling case [ fig .
[ fig : results : qmc_lambda2.0](d ) ] , there exist incoherent excitations with almost zero energy .
this indicates the metallic character of the system expected for large polaron densities and intermediate ep coupling .
( color online ) one - electron spectral functions @xmath121 and @xmath122 from ed for @xmath55 , @xmath123 , and @xmath81 . here
@xmath72 ( @xmath113 : red dot - dashed lines , @xmath114 : black solid lines ) and @xmath124 ( @xmath113 : green dashed two - dotted lines , @xmath114 : blue dotted lines ) .
dashed lines correspond to the integrated spectral weight of @xmath121 and @xmath122 , respectively.,scaledwidth=45.0% ]    until here , we have only presented results for the adiabatic case @xmath61 .
although the latter is most relevant for strongly correlated materials , it is important to compare the above findings to the antiadiabatic strong - coupling regime . to this end , we chose @xmath64 as well as @xmath82 . the condition for the formation of a single small polaron in this case , however , is @xmath125 , because the phonon frequency is the important energy scale .
in contrast , the criterion @xmath126 is purely based on the balance of kinetic and potential energy.@xcite    the spectral function for @xmath72 and @xmath124 is displayed in fig .
[ fig : results : w4 ] for @xmath127 , i.e. , at the transition point where self - trapping occurs for a single particle .
contrary to the case @xmath92 and @xmath61 considered above , the absorption features are now separated by the phonon energy @xmath8 , with their width being determined by electronic excitations .    a comparison of the spectral functions for @xmath72 and @xmath124 in fig .
[ fig : results : w4 ] indicates that there is no density - driven cross over of the system as observed in the adiabatic case . in particular ,
owing to the large phonon energy , there are no low - energy excitations close to the polaron band , so that the latter remains well separated from the incoherent features even for @xmath124 .
furthermore , the spectral weight of the polaron band also remains almost unchanged as we increase the density from @xmath72 to @xmath124 .
consequently , almost independent small polarons are formed also at finite electron densities , in accordance with previous findings of capone _
et al._@xcite for small systems .
normalized kinetic energy as a function of ep coupling @xmath80 for different band fillings @xmath39 and adiabatic ratios @xmath128 . also shown
are the results for a single polaron , obtained with the method of ref .  .
errorbars are smaller than the symbols used and @xmath31.,scaledwidth=45.0% ]    previously , the electronic kinetic energy has been studied extensively to monitor the polaron and bipolaron cross over in the holstein model with one and two electrons ( see , e.g. , refs .   and  ) . here ,
at least in the adiabatic regime , polaron formation manifests itself in a drop of the kinetic energy around the critical coupling . in the following , however ,
we illustrate that the kinetic energy [ eq .
( [ eq : app : ek ] ) ] does not reflect the significant changes in the one - particle spectrum we found with increasing density in the intermediate coupling regime ( see sec .
[ sec : results : intermediate ] ) .    to compare different band fillings
, we consider the normalized quantity @xmath129 , shown in fig .
[ fig : results : ek ] as a function of ep coupling for different fillings and phonon frequencies , for a fixed system size @xmath130 and @xmath26 .    in the adiabatic regime @xmath61 ,
owing to the overlap of the phonon clouds surrounding the individual quasiparticles , the carriers clearly become more mobile with increasing density , but the change in the kinetic energy is nevertheless surprisingly small .
in addition to these two fillings , we have also included in fig .
[ fig : results : ek ] the result for a single polaron with the same parameters , calculated with the qmc method of ref .  , which is indistinguishable from the result for @xmath72 .
as expected from the almost identical spectra for @xmath72 and 0.3 presented in sec .
[ sec : results : nonadiabatic ] for the antiadiabatic regime , @xmath131 is almost independent of @xmath39 in this case .
strongly correlated electron - phonon systems in dependence on the carrier density have been studied by means of unbiased numerical methods . using exact diagonalization together with the kernel polynomial method , as well as a novel quantum monte carlo approach
, we were able to present accurate results for the one - dimensional spinless holstein model in the most difficult regime of small phonon frequencies and intermediate to strong electron - phonon interaction .
we have calculated the one - particle spectral function , the density of states and the electronic kinetic energy at different band fillings .
while exact diagonalization is restricted to relatively small clusters , thereby limiting the momentum resolution of the spectra , quantum monte carlo can be used to study larger systems , but yields a lower energy resolution for dynamical quantities .
the reliability of the results from the maximum entropy inversion has been scrutinized , and we find a very good agreement between the two methods .    in the adiabatic case ,
characteristic for , e.g. , the manganites , we observe a carrier density - driven cross over from a polaronic state to a metallic system at weak and intermediate electron - phonon coupling .
note that the nature of the quasiparticles is changed in the course of this transition .
this points toward an insufficiency of simple single - polaron theories to explain experimental results for polaronic materials .
in fact , very similar effects due to carrier interaction have recently been observed experimentally for the manganites.@xcite    on the contrary , for large phonon frequencies and strong coupling , the individual polarons remain virtually unaffected upon increasing the electron density .
consequently , the system remains polaronic even at large band fillings .
finally , the kinetic energy is shown to be too insensitive in order to reveal the above cross over .    for weak electron - phonon coupling ,
the system becomes a peierls or band insulator at half filling , whereas for strong coupling a polaronic insulating state with an exponentially small polaron band at the fermi level , and a finite gap to the incoherent excitations , prevails for all fillings .
finally , reliable calculations of the optical conductivity are highly desirable in order to explain the large amount of available experimental data for polaronic materials .
work along this line is in progress .
this work was supported by the austrian science fund ( fwf ) , project no .
p15834 , the deutsche forschungsgemeinschaft through spp1073 , and the bavarian competence network for high performance computing ( konwihr ) .
m.  h. is grateful to doc ( doctoral scholarship program of the austrian academy of sciences ) , and to hpc - europa .
h.  f. acknowledges the hospitality at graz university of technology . we would like to thank a. alvermann , h. g. evertz , t. lang and a. weie for useful discussion .
furthermore , we acknowledge generous computer time granted by the lrz munich .
in this appendix we extend the one - electron qmc algorithm developed in ref .  , henceforth also referred to as i , to the spinless holstein model with many electrons .
we begin by applying the lang - firsov transformation@xcite @xmath132 with @xmath133 to the hamiltonian  ( [ eq : model : h ] ) in first quantization and with dimensionless phonon variables @xmath134 and @xmath135 , which takes the form@xcite @xmath136 the coupling constant @xmath137 is related to that in eq .
( [ eq : model : h ] ) by @xmath138 . using the result @xmath139 ( see i ) , we define the grand - canonical hamiltonian @xmath140 where @xmath34 denotes the chemical potential and @xmath15 is the polaron binding energy
.      we would like to point out that in the spinfull case , the transformed hamiltonian would contain an attractive hubbard term .
for the determinant qmc method to be applicable , the latter has to be decoupled using auxiliary fields.@xcite in contrast , for the spinless model considered here , the algorithm is almost identical to the one - electron case .
the following derivation assumes a hypercubic lattice in d dimensions , consisting of @xmath142 sites .
we use the suzuki - trotter decomposition@xcite @xmath143 with @xmath144 .
the trace appearing in the partition function @xmath145 can be split up into a bosonic and a fermionic component leading to @xmath146 where @xmath147 , and with periodic boundary conditions @xmath148 .
we have chosen the time ordering , which is arbitrary , in increasing order .
the phonon coordinates @xmath134 in @xmath149 can be integrated out analytically in the same manner as in i. moreover , the momenta @xmath135 can be replaced by their eigenvalues on each time slice , and the partition function takes the form @xmath150 with @xmath151 , @xmath152^{nl/2}$ ] and @xmath153 the bosonic action @xmath154 may be expressed in terms of principal components ( see i ) .
the fermion degrees of freedom can be integrated out exactly leading to@xcite @xmath155 where the @xmath156 matrix @xmath157 is given by @xmath158 with @xmath159 here @xmath160 is the usual tight - binding hopping matrix.@xcite to save computer time , we employ the checkerboard breakup@xcite @xmath161 using eq .
( [ eq : app : checker ] ) , the numerical effort scales as @xmath162 instead of @xmath163 , while the error due to this additional approximation is of the same order @xmath164 as the trotter error in eq .
( [ eq : app : suzuki - trotter ] ) .    defining the bosonic and fermionic weights @xmath165 and @xmath166 respectively , the partition function finally becomes @xmath167 one of the advantages of blankenbecler _
et al._s@xcite formalism as well as the current approach is the close relation to the one - electron green function @xmath168 working in real space and imaginary time , we have@xcite @xmath169 and @xmath170 we would like to mention that despite the formal similarity to the method of blankenbecler _
et al._,@xcite the numerical realization of the present approach is quite different . while the grand - canonical method of ref .
benefits enormously with respect to performance from a local updating scheme for the phonons , here we use a global updating in terms of the principal components together with the reweighting method described in i. although this requires us to recalculate the full matrix @xmath171 in each mc step , the resulting statistically independent configurations clearly outweigh the loss in performance , especially for small @xmath8 for which autocorrelation times can exceed @xmath172 sweeps when using the original grand - canonical qmc method.@xcite the approach proposed here allows one to perform uncorrelated measurements after each update .
moreover , our calculations show that numerical stabilization by means of a time - consuming singular value decomposition is not necessary for all parameters considered in sec .
[ sec : results ] . finally , owing to the phase factors in the transformed hopping term , all the matrices become complex - valued . an important practical test in simulations
therefore is the reality of expectation values of observables , i.e. , the disappearance of the imaginary part within statistical errors .    the phase factors in the hopping term [ eq .
( [ eq : app : hspinless ] ) ] give rise to a minus - sign problem , i.e. , @xmath173 can take on negative values , whereas @xmath174 is always positive .
while for a fixed number of electrons the sign problem diminishes with increasing system size,@xcite the opposite is true if we fix the electron density . in that case
, we find an exponential decrease of the average sign with increasing cluster size , and a strong dependence on the band filling ( sec .
[ sec : app : sign_results ] ) .
the calculation of observables within the formalism presented here is similar to the standard determinant qmc method.@xcite an important difference is that we have to use the lang - firsov - transformed operators , i.e. , @xmath175 as mentioned above , the mc sampling is based purely on the bosonic weight @xmath174 .
this corresponds to a reweighting of the probability distribution so that the fermionic weight @xmath173 is treated as part of the observables.@xcite as for a single electron,@xcite the overlap of the probability distributions involved is large enough to permit this procedure .
hence , we have @xmath176 where the expectation value @xmath177 for an equal - time observable @xmath178 is defined as @xmath179      we begin with the electron density @xmath180 the expectation value in eq .
( [ eq : app : n ] ) can be calculated from the diagonal elements of the green function @xmath181 [ eq .  ( [ eq : app : ga ] ) ] , i.e. , @xmath182 .
similarly , the absolute value of the electronic kinetic energy per site is given by @xmath183 equal - time two - particle correlation functions such as @xmath184 may be calculated with the help of wick s theorem.@xcite      an observable of great interest is the time - dependent one - particle green function @xmath185 which is related to the spectral function @xmath186 through @xmath187 the inversion of the above relation is ill - conditioned and requires the use of the mem.@xcite fourier transformation leads to @xmath188 the allowed imaginary times are @xmath189 , with nonnegative integers @xmath190 . within the qmc approach ,
we have@xcite @xmath191 another interesting quantity is the one - electron density of states @xmath192 where @xmath193 .
it may be obtained numerically via @xmath194 and subsequent use of the mem .
average sign @xmath195 [ eq .  ( [ eq : app : sign ] ) ] of the fermionic weight @xmath173 as a function of electron - phonon coupling @xmath80 in one dimension ( a ) for different inverse temperatures @xmath95 , and ( b ) for different values of the adiabatic ratio @xmath128 .
here and in subsequent figures @xmath31 , lines are guides to the eye only , and errorbars are suppressed if smaller than the symbols shown.,title="fig:",scaledwidth=35.0% ] +   average sign @xmath195 [ eq .  ( [ eq : app : sign ] ) ] of the fermionic weight @xmath173 as a function of electron - phonon coupling @xmath80 in one dimension ( a ) for different inverse temperatures @xmath95 , and ( b ) for different values of the adiabatic ratio @xmath128 .
here and in subsequent figures @xmath31 , lines are guides to the eye only , and errorbars are suppressed if smaller than the symbols shown.,title="fig:",scaledwidth=35.0% ]        while hamiltonian  ( [ eq : app : hspinless ] ) is symmetric with respect to a particle - hole transformation for @xmath141 ( half filling ) , this symmetry is broken if we use the checkerboard approximation  ( [ eq : app : checker ] ) , so that @xmath196 for the above choice of the chemical potential . to simplify calculations ,
the results for the average sign at @xmath79 presented below have therefore been obtained using the exact hopping term . for general band fillings , we find that the checkerboard breakup requires different values of @xmath34 to perform simulations at the same electron densities , but results are identical within statistical errors if @xmath34 is adjusted accordingly .    the derivation of the qmc algorithm in sec .
[ sec : app : qmc_partition ] is valid for any dimension d of the lattice under consideration .
nevertheless , here we only report results for the sign problem in @xmath197 , the case considered in sec .
[ sec : results ] , and make some remarks about the influence of the dimensionality at the end .    the average sign is defined as @xmath198 with the fermionic weight in the present case given by @xmath166 [ eq .  ( [ eq : app : omega ] ) ] .
we begin with the dependence of @xmath195 on the electron - phonon coupling strength . for simplicity , we show results for @xmath79 , while the effect of band filling will be discussed later .
the choice @xmath79 is convenient since we know the chemical potential , and we shall see below that the sign problem is most pronounced for a half - filled band .
moreover , all existing qmc results for the spinless holstein model are for half filling , so that it is interesting to see how the sign problem affects simulations within the current approach .
average sign as a function of ( a ) band filling @xmath39 , and ( b ) system size @xmath91.,title="fig:",scaledwidth=35.0% ] +   average sign as a function of ( a ) band filling @xmath39 , and ( b ) system size @xmath91.,title="fig:",scaledwidth=35.0% ]    figure  [ fig : app : sign_beta_omega](a ) reveals that the average sign takes on a minimum near the point where  depending on the filling
the cross over to small polarons or a cdw state occurs . in the adiabatic regime @xmath60 ,
the cross over condition is @xmath126 [ fig .
[ fig : app : sign_beta_omega](a ) ] , whereas for @xmath199 we have @xmath200 [ fig .
[ fig : app : sign_beta_omega](b ) ] . at weak and strong coupling ,
the average sign is close to 1 , so that accurate simulations can be carried out . as shown in fig .
[ fig : app : sign_beta_omega](b ) , we find that @xmath195 becomes very small for @xmath60 and intermediate @xmath80 , while it increases noticeably in the nonadiabatic regime @xmath199 .
owing to the absence of any autocorrelations in our approach , the present method therefore represents a significant improvement of existing algorithms , as the latter face severe autocorrelations even for @xmath201 , thereby limiting the temperature range and cluster size .
as illustrated in fig .
[ fig : app : sign_n_n](a ) , the average sign depends strongly on the band filling @xmath39 . while @xmath202 in the vicinity of @xmath203 or @xmath204 , a significant reduction is visible near half filling @xmath79 .
the minimum occurs at @xmath79 , and the results display particle - hole symmetry as expected .
here we have chosen @xmath26 , @xmath61 and @xmath92 , for which the sign problem is most noticeable according to fig .
[ fig : app : sign_beta_omega ] . note that there is no sign problem in simulations for , e.g. , the hubbard model at half filling , where the ( fermionic ) weights of the @xmath205 and @xmath206 electrons are equal and their product  entering the partition function  is positive.@xcite this
can also be expected within the current approach if we consider the half  filled spinfull holstein or holstein - hubbard model , and work along this line is in progress .    finally , in fig .
[ fig : app : sign_n_n](b ) , we report the average sign as a function of system size , again for @xmath79 .
the behavior is strikingly different from the one - electron case considered in i. while in the latter @xmath207 as @xmath208 , here the average sign decreases nearly exponentially with increasing system size .
obviously , this limits the applicability of our method .
however , as illustrated in sec .
[ sec : results ] , we can nevertheless obtain accurate results at low temperatures , small phonon frequencies , and over a large range of the electron - phonon interaction .
moreover , we would like to point out that for such parameters , existing methods suffer strongly from autocorrelations , rendering simulations extremely difficult .
the dependence of the sign problem on the dimension of the system is similar to the single - electron case.@xcite for the same parameters @xmath91 , @xmath128 , @xmath209 and @xmath80 , the minimum in @xmath195 at intermediate @xmath80 becomes more pronounced as one increases the dimension of the cluster .    to conclude with
, we would like to point out that , in principle , the sign problem can be compensated by performing sufficiently long qmc runs , although we have to keep in mind that the statistical errors increase proportional to @xmath210 ( ref .  ) .
owing to the purely phononic updates , the present algorithm is very fast and we have made up to about @xmath211 single measurements to obtain the results presented in sec .  [
sec : results ] .
the use of the reweighting method requires us to store the measured values of observables in every mc step , and to perform the jackknife error analysis@xcite at the end of the run .
for example , in the case of the green function  ( [ eq : app : gb_tau ] ) , @xmath212 complex numbers have to be stored for each momentum @xmath53 in every step .
depending on the observables of interest , this sets a practical limit for the maximal number of measurements due to restrictions in available disk space .