the theory of high - temperature superconductivity remains one of the most challenging problems in solid - state physics . in many metallic compounds , over a broad range of compositions and temperatures ,
only two phases are encountered , i.e.  the normal ( fermi liquid ) and a magnetic or superconducting phase . in sharp contrast , the high - temperature superconductors ( htsc ) and other strongly correlated electron systems , such as heavy - fermion and a variety of transition - metal oxide systems , exhibit many ordered phases , which appear to compete and sometimes coexist.@xcite in the htsc , besides ordered antiferromagnetic ( af ) and superconducting ( sc ) phases , compelling evidence exists for charge- and spin-``stripe '' phases and phases with coexisting sc and af order.@xcite the common denominator and the underlying reason for these competing orders is certainly the presence of strong electronic correlations and mott - hubbard physics : the interplay between kinetic - energy and coulomb - correlation effects induces an extreme sensitivity to external parameters ( doping , temperature , pressure , etc . ) and a rather difficult to predict `` outcome '' , i.e.  the characteristic low - energy excitations and the phase diagram at low temperatures .
the central challenge in the field of high-@xmath1 superconductivity is , therefore , the connection of the ( known ) microscopic interactions at the level of electrons and ions , which are at high energy ( @xmath2ev ) and temperature @xmath3 , with the `` emerging phenomena '' at @xmath4 , i.e.  competing and nearly degenerate orders .
ideally , one should employ a systematic renormalization - group approach to integrate out the irrelevant degrees of freedom and , thereby , correctly bridge high to low energies and eventually go to @xmath4 .
it is , however , by no means obvious how to do this when strong correlations are present , such as in the htsc .    in this context cluster techniques , which systematically approach the infinite - size ( low - energy ) limit ,
provide a powerful alternative .
@xcite in this paper , we will discuss and apply a variational - cluster approach ( vca ) which was proposed recently .
@xcite it is based on the self - energy - functional theory ( sft ) , @xcite which provides a general variational scheme to use dynamical information from an exactly solvable `` reference system '' ( in our case , an isolated cluster ) to go to the infinite - size lattice fermion problem at low temperatures and at @xmath4 , in particular . in our earlier work @xcite this scheme
was formulated to study phases with spontaneously broken symmetry . for the cluster sizes used
, it was shown that the vca correctly reproduces long - range af order for the two - dimensional ( @xmath5 ) hubbard model and the absence of this order in one dimension .
this non - trivial `` test '' implies that the vca goes well beyond ordinary mean - field theory .
another crucial test is provided by the dynamical information contained in the one - particle green s function @xmath6 .
compared to variational schemes based on wave functions , @xcite an important advantage of the vca consists in the fact that it quite naturally gives the one - electron green s function @xmath6 . for the @xmath5 hubbard model ,
it was recently demonstrated @xcite that the vca , with the lattice tiled by ( @xmath7 ) clusters , correctly reproduces low - temperature quantum monte - carlo ( qmc ) data , in particular , the coherent and incoherent `` bands '' experimentally known from arpes data .
@xcite    these tests provide also the foundation for attacking the question whether the `` minimal '' microscopic model , namely the @xmath5 one - band hubbard model , reproduces the essential features of the electron- and hole - doped htsc phase diagram .
we will not go into a lengthy discussion of what interactions should be retained at the electron - ion level .
but , when choosing the @xmath5 one - band hubbard model , @xcite i.e.@xmath8 where @xmath9 , @xmath10 are the usual annihilation and creation operators , @xmath11 denote the hopping - matrix elements , @xmath12 the density at site @xmath13 with spin @xmath14 and @xmath15 the local coulomb repulsion , one has introduced gross simplifications , leaving out other orbital ( e.g.  @xmath16 ) degrees of freedom , long - range coulomb interaction , electron - phonon coupling , etc . nevertheless , this choice appears to be legitimate , last not least in view of the amazing agreement achieved between numerical simulations and experimental results for the normal - state properties of the cuprates ( see , for example , refs .  ) .
the ground - state phase diagram of the model was recently investigated using the vca by snchal et al .
@xcite and by two of us .
@xcite there are important technical differences , but the `` upshot '' of the two works is as follows : for the cluster sizes used in the vca , the @xmath4 phase diagram of the hubbard model ( including hopping terms up to second or third - nearest neighbors ) turns out to be qualitatively similar to that of the electron- and hole - doped cuprates .
the model correctly describes the overall phase diagram , such as the occurrence of the af and sc phases and predicts the corresponding doping ranges in qualitative agreement with the experiments for the cuprate materials .
the present paper has several purposes : first , we would like to stress that for an application of the vca to the high-@xmath1 problem it is of crucial importance to treat the on - site energies in the reference system as variational parameters .
we will show that this ensures a thermodynamically consistent determination of the average particle number .
compared to the study of ref .
this represents an important methodical extension . on the other hand , without the variational optimization of the on - site energy , one has to tolerate an inconsistency in the determination of the average particle number .
the effects of this error shall be demonstrated by model calculations .
the issue of thermodynamic consistency is also discussed for the off - diagonal elements of the one - particle density matrix and , in case of spontaneous symmetry breaking , for the respective order parameter .
it is interesting to note that there are no such problems in the dynamical mean - field theory ( dmft ) @xcite and its cluster extensions .
@xcite here the on - site energies are kept at their `` physical '' values from the very beginning rather than being determined from the self - consistency condition .
the ( cellular ) dmft , however , can be considered as a special approximation within the general sft framework .
@xcite hence , the question arises why a fixed on - site energy does _ not _ spoil thermodynamic consistency in the case of ( cellular ) dmft .
it is interesting to note that , in case of interacting bosons , the issue of consistency has also been shown to be very important , recently .
@xcite    second , an accurate analysis of the behavior of the chemical potential as a function of the particle density close to the transition to a non - magnetic state , as well as a corresponding maxwell construction , indicate the presence of an inhomogeneous ground state with macroscopically large regions of low- and high - particle density .
using @xmath17-site clusters , we can get a rough estimate of the ground - state phase diagram and investigate the instability of the homogeneous ( af , sc ) phases against charge inhomogeneities .
this represents an important complement to the work of ref .  , where larger cluster sizes up to @xmath18 sites have been considered but without an appropriate analysis via a maxwell construction .
finally , our numerical results give valuable insights into different questions of the high-@xmath1 problem , as they provide direct access to the single - particle excitation spectrum in the strong - coupling regime at zero temperature .
this has to be contrasted with the work of maier et al .
@xcite who have been able to treat clusters with @xmath19 sites within the dynamical cluster approximation ( dca ) , but are restricted to finite temperatures , intermediate coupling @xmath20 and imaginary time .
a review comparing the application of different cluster ( sft , dca , and cellular dmft ) as well as weak - coupling methods for the hubbard model has appeared recently .
@xcite    our paper is organized as follows : we start with a summary of the central ideas of the sft in general ( sec .  [
sec : sft ] ) .
thermodynamic consistency with respect to the average particle number is discussed in sec .
[ sec : n ] .
consistency with respect to the off - diagonal elements of the density matrix is addressed in appendix [ sec : cicj ] and the case of the c - dmft is discussed in appendix [ sec : cdmft ] .
our theoretical considerations are completed by describing some computational details of the vca calculations in sec .
[ sec : vca ] .
section [ spe ] then presents the results for the @xmath4 spectral function @xmath21 for electron- and hole - doping .
we discuss how the characteristically different doping dependencies of the spectra give rise to different fermi - surface evolutions upon doping .
these fermi - surface evolutions can then be tied up with the characteristic differences in the electron- and hole - doped phase diagrams , such as the enhanced robustness of the af order in the electron - doped case .
the ground - state phase diagram will be presented and discussed in detail in sec .
finally , sec .
[ summ ] contains our main conclusions and a summary .
the central idea of the self - energy - functional theory ( sft ) @xcite is to make use of the universality of the luttinger - ward functional @xmath22 $ ] @xcite or of its legendre transform @xmath23 $ ] : for a system with hamiltonian @xmath24 , where @xmath25 are the one - particle and @xmath26 the interaction parameters , the functional dependence @xmath27 $ ] or @xmath28 $ ] is independent of @xmath25 .
this universality is obvious as the luttinger - ward functional is defined via a skeleton - diagram expansion involving dressed propagators and vertices only .
@xcite    concentrating on the self - energy @xmath29 instead of the single - particle green s function @xmath6 , the grand potential of the system at temperature @xmath3 and chemical potential @xmath30 can be written as a functional of @xmath29 : @xmath31 =   { \mbox{tr}}\ln ( { { \bm g}}_{0 , { { \bm t}}}^{-1 } -                  { { \bm \sigma}})^{-1 } + f_{{{\bm u } } } [ { { \bm \sigma } } ] \ : , \label{eq : sf}\ ] ] where @xmath32 is the free green s function and @xmath33 with the usual trace @xmath34 and the matsubara frequencies @xmath35 for integer @xmath36 . at the physical self - energy @xmath37 ,
the grand potential is stationary : @xmath38 = 0 $ ] .
why is it more advantageous to express @xmath39 as a functional of the self - energy @xmath29 rather than @xmath6 ?
this has to do with the `` short - range '' character of @xmath29 as a function of its real - space coordinates , @xcite which in general is due to the fact that @xmath29 is qualitatively related to a dynamically screened particle - particle interaction . for the hubbard model , in particular
, one believes that important effects are sufficiently accounted for by a local @xcite or short - ranged @xcite self - energy , at least for high lattice dimensions .
a local or short - ranged self - energy , however , can be well generated by a cluster of finite size , and for the subsequent optimization of the cluster trial self - energy , the self - energy functional ( [ eq : sf ] ) can be used .
this concept allows for the construction of a class of conceptually clear and thermodynamically consistent approximations , including the dynamical mean - field theory ( dmft ) @xcite and a cluster extension of the dmft @xcite ( see ref .   for a detailed discussion ) .    due to the universality of @xmath40
$ ] , we have @xmath41 = f_{{{\bm u } } } [ { { \bm \sigma } } ] + { \mbox{tr}}\ln          ( { { \bm g}}_{0,{{\bm t}}'}^{-1}-{{\bm \sigma}})^{-1 } \label{eq : sfp}\ ] ] for the self - energy functional of a so - called `` reference system '' , which is given by a hamiltonian with the same interaction part @xmath26 but modified one - particle parameters @xmath42 : @xmath43 .
although it has different microscopic parameters , the reference system is assumed to be in the same macroscopic state as the original system , so it has the same temperature @xmath3 and the same chemical potential @xmath30 . by a proper choice of its one - particle part ,
the problem posed by the reference system @xmath44 can be much simpler than the original problem posed by @xmath45 , such that the self - energy of the reference system , @xmath46 , can be computed exactly within a certain subspace of parameters @xmath42 .
combining eqs .
( [ eq : sf ] ) and ( [ eq : sfp ] ) , we can eliminate the functional @xmath40 $ ] . inserting as a trial self - energy the self - energy of the reference system ,
then yields : @xmath47           =           \omega_{{{\bm t } } ' , { { \bm u } } }          +          { \mbox{tr}}\ln ( { { \bm g}}_{0,{{\bm
t}}}^{-1}-{{\bm \sigma}}_{{{\bm t}}',{{\bm u}}})^{-1 }          -           { \mbox{tr}}\ln { { \bm g}}_{{{\bm t}}',{{\bm u } } } \ : , \label{eq : ocalc}\ ] ] where @xmath48 and @xmath49 are the grand potential and the green s function of the reference system .
stationary points are obtained , and this is the approximation , by restricting the variation to the subspace of trial self - energies @xmath46 : @xmath50 }    { \partial { { \bm t } } ' } = 0 \qquad \mbox{for } \quad { { \bm t } } ' = { { \bm t}}'_{\rm s } \ : .
\label{eq : stat}\ ] ] varying the trial self - energy means to vary the one - particle parameters @xmath42 of the reference system .
for further details of the approach see ref .  .    here
, we will focus on the hubbard model , eq .
( [ eq : hub ] ) , as the original model given by @xmath45 .
different possible choices for @xmath44 and the corresponding systematics of approximations generated in this way are discussed in ref .  .
in the following , we will concentrate on two cluster approaches : ( i ) the variational cluster approach ( vca ) @xcite and ( ii ) the cellular dmft ( c - dmft ) .
@xcite the vca can be seen as a variational generalization of the cluster - perturbation theory .
@xcite it is obtained by partitioning the infinite lattice into disconnected ( identical ) clusters of @xmath51 sites each and choosing @xmath44 to consist of the intra - cluster parts only , i.e.  the inter - cluster hopping is switched off in @xmath44 .
the c - dmft is a cluster variant of the dmft . in the context of the sft
, it is obtained in the same way as the vca but with an additional coupling of each of the @xmath51 correlated cluster sites to a continuous bath , i.e.  to an infinite number of uncorrelated additional bath sites .
the on - site energies of the bath sites as well as their coupling to the original sites are treated as variational parameters .
once a reference system is specified , one should , in principle , vary _ all _ one - particle parameters of @xmath44 .
this procedure would give the optimal result but requires a search for a stationary point in a high - dimensional parameter space .
from a pragmatic point of view it is thus advisable to concentrate on a few parameters only which have to be selected by physical arguments . here
, we argue that the variation of the on - site energies is important to achieve thermodynamic consistency with respect to the average particle number . in case of the hubbard model , this means to consider the site - independent energy @xmath52 as ( one of the ) variational parameter(s ) .
the average particle number @xmath53 can be calculated in two different ways : on the `` zero - particle level '' by differentiation of the grand potential with respect to the chemical potential @xmath30 : @xmath54 and on the `` one - particle level '' by frequency integration of the one - particle excitation spectrum : @xmath55 here , @xmath56 is the local ( possibly spin - dependent ) spectral density with @xmath57 for short , and @xmath58 is the fermi function . for simplicity ,
we exclude non - collinear magnetic states and assume all expectation values to be diagonal in the spin index .
thermodynamic consistency means that both ways of calculating @xmath53 yield the same result .
since @xmath59 $ ] is the _ approximate _ sft grand potential at the stationary point @xmath60 , and since the spectral density or , equivalently , the one - particle green s function is the _ approximate _ green s function given by @xmath61 , the equivalence of ( [ eq : n1 ] ) and ( [ eq : n2 ] ) is by no means understood _ a priori_.    to prove thermodynamic consistency , we start from eq .
( [ eq : n1 ] ) . according to eq .
( [ eq : ocalc ] ) , there is a twofold @xmath30 dependence of @xmath62 $ ] : ( i ) an _ explicit _ @xmath30 dependence due to the chemical potential in the free green s function of the original model , @xmath63 , and ( ii ) an _ implicit _ @xmath30 dependence due to the @xmath30 dependence of the self - energy @xmath64 , the green s function @xmath65 and the grand potential @xmath66 of the reference system : @xmath67 note that the implicit @xmath30 dependence is due to the chemical potential of the reference system which , by construction , is in the same macroscopic state as the original system ( with the same temperature @xmath3 and the same chemical potential @xmath30 ) _ as well as _ due to the @xmath30 dependence of the stationary point @xmath68 itself .
this is a subtlety which , however , can be ignored since @xmath69 for @xmath60 because of stationarity condition ( [ eq : stat ] ) .
( actually , only those elements of @xmath68 show up a @xmath30 dependence that are treated as variational parameters . according to the chain rule
, however , the derivative of @xmath39 has to be performed just with respect to those elements , with a vanishing result due to the stationarity condition . )
the self - energy , the green s function and the grand potential of the reference system are defined as grand - canonical averages .
hence , their @xmath30 dependence due to the grand - canonical hamiltonian @xmath70 is ( apart from the sign ) the same as their dependence on @xmath71 : @xmath72 , etc .
consequently , we have : @xmath73 the first derivative is readily calculated : @xmath74 here , the contour @xmath75 encloses the first - order poles of the fermi function at @xmath76 in counter - clockwise direction .
using cauchy s theorem , we can proceed to an integration over real frequencies . inserting into eq .
( [ eq : exim ] ) , we get : @xmath77 for @xmath60 .
the first term on the right - hand side is just the expression for the average particle number given by ( [ eq : n2 ] ) .
the second term on the right - hand side vanishes _ provided that _ the variational condition ( [ eq : stat ] ) is satisfied , i.e.  _ provided that _
@xmath71 is included in the set of variational parameters . in this case one
has thermodynamic consistency .
if @xmath71 was not treated as a variational parameter but kept at the value given by the original system , @xmath78 , one would have a finite @xmath79 in ( [ eq : nn ] ) , and the two expressions ( [ eq : n1 ] ) and ( [ eq : n2 ] ) for the average particle number would yield different results .
this completes the proof .
( [ eq : n1 ] ) for the average particle number of the lattice model can be compared with @xmath80 which gives the average particle number of the reference system .
again , there are two ways to get @xmath81 : either as the derivative of the reference system s grand potential @xmath82 or by frequency integration of the reference system s spectral density @xmath83 .
as the reference system is solved exactly , both ways must yield the same result .
note , however , that @xmath84 in general .
the above reasoning can straightforwardly be generalized to the off - diagonal ( @xmath85 ) elements of the one - particle density matrix @xmath86 .
this is discussed in appendix [ sec : cicj ] .    the effect of thermodynamic ( in)consistency is illustrated for the single - band hubbard model of eq .
( [ eq : hub ] ) in fig .
[ checkn ] .
the figure shows @xmath87 as a function of @xmath30 as obtained from vca calculations described in sec .
[ sec : vca ] below .
solid lines display the result obtained by frequency integration of the spectral density , eq .
( [ eq : n2 ] ) , while dashed lines show the result of the numerical @xmath30 derivative of @xmath39 , eq .
( [ eq : n1 ] ) .
we compare the results obtained by considering @xmath71 as a variational parameter ( a ) , with the ones obtained by setting @xmath88 ( b ) . in case ( a )
the two curves coincide within numerical accuracy , as expected , while in case ( b ) a considerable discrepancy is observed .
this discrepancy increases upon approaching the transition to the non - magnetic state , i.e.  precisely in the interesting region , where it reaches about @xmath89 doping .     filling @xmath87 as a function of chemical potential @xmath30 obtained by integration of the spectral density ( eq .  ( [ eq : n2 ] ) , solid lines ) and via the derivative of @xmath39 ( eq .  ( [ eq : n1 ] ) , dashed lines ) .
results are obtained , respectively , by considering @xmath71 as a variational parameter ( a ) , and by setting @xmath88 ( b ) .
calculations for the electron - doped case with @xmath90 , nearest - neighbor hopping @xmath91 and next - nearest - neighbor hopping @xmath92 . ]    while thermodynamic consistency with respect to the particle number is an issue for most approximations within the sft , there is one exception : the cellular dmft .
this is discussed in appendix [ sec : cdmft ] .      in detail
, the calculation proceeds as follows : we consider the @xmath5 single - band hubbard model of eq .
( [ eq : hub ] ) with nearest- ( @xmath93 ) and next - nearest ( @xmath94 ) neighbor hoppings and the hubbard @xmath15 . for the purpose of our qualitative description of the htsc , it is sufficient to take parameters which are typical for both hole- and electron - doped high-@xmath1 cuprates , i.e.  @xmath95 and @xmath96 .
the energy scale is set by choosing @xmath91 .
different values of the parameters ( e.g. a change of @xmath15 or @xmath94 within @xmath97 ) or a third - nearest neighbor hopping have been also incorporated , checked and found not to qualitatively change our conclusions ( for example , concerning the weak phase separation detected in the electron - doped case ) .
the hamiltonian of the reference system @xmath44 is given by a set of decoupled clusters of finite size .
for an individual cluster , the hamiltonian reads : @xmath98 it consists of the hubbard hamiltonian confined to the finite cluster @xmath99 plus two symmetry - breaking terms ( weiss fields ) @xmath100 and @xmath101 with @xmath102 and @xmath103 where @xmath104 is the strength of the staggered and @xmath105 the strength of the nearest - neighbor @xmath0-wave pairing field .
@xmath106 is the af wave vector , and @xmath107 denotes the @xmath0-wave form factor which is non - vanishing for nearest - neighbor lattice sites only and is equal to @xmath108 ( @xmath109 ) for @xmath110 in x ( y ) direction .
the sum in eq .
( [ hsc ] ) is restricted to sites @xmath13 and @xmath111 belonging to the same cluster . according to the discussion in sec .
[ sec : n ] , the site - independent energy @xmath112 shall be treated as a variational parameter .
it shows up in the local term @xmath113 which is already included in @xmath99 .
the optimization of @xmath71 has to be done simultaneously with the optimization of the parameters @xmath104 and @xmath105 .    due to the optimization of the weiss field strengths @xmath104 and @xmath105
, one can account for spontaneous af and @xmath0-wave sc symmetry breaking .
the respective af and sc order parameters , @xmath114 and @xmath115 , are defined as @xmath116 in the limit @xmath117 where @xmath118 are the strengths of external _ physical _ staggered and pairing fields .
these physical fields should not be confused with the ficticious weiss fields with strengths @xmath104 and @xmath105 .
adding the respective physical field terms to the hamiltonian @xmath45 and performing the derivative with respect to @xmath119 of the sft grand potential ( at the respective optimal ficticious field strengths @xmath120 ) , yields @xmath114 and @xmath115 _ consistently _ with their representations ( on the `` one - particle level '' ) as frequency integrals over the usual and anomalous one - particle spectral density .
this consistency is shown in the same way as in sec .
[ sec : n ] for the average particle number and is a consequence of treating the ficticious fields as variational parameters .
the quality of the approximation is decisively influenced by the cluster size used . on the one hand , for an appropriate characterization of the phase transition ( as considered in sec .
iv , below ) , one needs a sufficient accuracy in the grand potential @xmath39 .
this accuracy is , first of all , determined by the requirement that the clusters chosen must be large enough to fully account for the `` short - range '' spatial dependence of the self - energy as has been discussed below eq .
( [ eq : sf ] ) . on the other hand , for a given cluster size , an as accurate as possible numerical evaluation has to be employed .    with respect to the latter , we found it both convenient and accurate to evaluate the frequency integrals contained in the trace in eq .
( [ eq : ocalc ] ) not by numerical integration ( where the required accuracy is difficult to achieve ) , but by converting these integrals to a sum over the poles of the green s function ( see ref .   for details ) .
this , however , requires the computation of all many - body eigenstates in a given sector of the cluster hamiltonian @xmath121 ( eq .
( [ eq : ham ] ) ) including symmetry - breaking fields .
this technically limits the cluster sizes to be considered .
therefore , in the present work , we have chosen an infinite lattice tiled with @xmath122 clusters .
cluster consisting of @xmath122 sites , as well as larger cluster sizes have recently been systematically studied by kyung et al .
@xcite in their influence on various physical quantities of the @xmath5 hubbard model by means of the c - dmft .
already the smallest , i.e.  @xmath122 , cluster has been found to account for more than 95% of the `` correlation effect '' in the single - particle spectrum ( for a precise definition see eq .
( 16 ) of ref .  ) .
this suggests that at least some of the relevant questions in a strongly correlated electron system , modelled by a hubbard - type hamiltonian , may be studied , even rather accurately , with a such a small reference cluster .
nevertheless , for a more accurate determination of the sc state and of the phase transition , larger clusters have to be considered eventually . only in this way
, phase - fluctuation effects between different d - wave cooper pairs ( which determine the spatial dependence of the self energy in the sc state ) can be accounted for .
we expect , for example , that calculations for larger clusters would display a smaller sc gap than the one seen in fig .
[ spectrelhol]b around @xmath123 in the electron - doped case . on the other hand
, we expect our results on the global phase diagram to be quite robust .
this is corroborated by the only weak dependence on ( moderate ) model changes , such as the next - nearest hopping or the hubbard interaction.@xcite
and @xmath124 ) case at @xmath90 ( @xmath91 ) . for comparison with the qmc data ,
the next - nearest - neighbor hopping is set to @xmath125 .
, title="fig : " ] +   and @xmath124 ) case at @xmath90 ( @xmath91 ) . for comparison with the qmc data , the next - nearest - neighbor hopping is set to @xmath125 .
, title="fig : " ] +    the evolution of the single - particle excitation spectrum @xmath126 as a function of hole doping , and , in particular , the transition from a `` small '' fermi surface ( hole pockets around ( @xmath127 ) ) to an lda - like fermi surface closed around @xmath128 are key observations in unlocking the mystery of the cuprates .
earlier qmc calculations for the @xmath5 hubbard model @xcite found that the single - particle spectral weight @xmath21 semi - quantitatively reproduces both the momentum ( @xmath129-symmetry ) and , in particular , doping dependence of the `` high - energy '' pseudogap of the order of the exchange energy @xmath130 as found in photoemission experiments around ( @xmath131 ) .
the corresponding qmc data are reproduced in fig .
[ qmc ] for comparison with the results of the vca shown in fig .
[ spectrtp0 ] for hole dopings @xmath132 and @xmath124 .
we first discuss the qmc data . in the under - doped regime ,
the `` pseudogap '' feature near ( @xmath131 ) moves to lower binding energy as doping is increased . at about @xmath133 ( full circle in fig .
[ qmc ] ) the pseudogap vanishes in overall accordance with the experimental findings .
@xcite in the experiments , the `` high - energy '' pseudogap is identified with the centroids of spectral weight near ( @xmath131 )
.    however , the experimental arpes spectra also display a `` low - energy '' pseudogap in the normal ( superconducting ) state above ( below ) @xmath1 with energy @xmath134 , inferred from the leading edge in the spectral density , which also opens up in the under - doped regime and vanishes in the over - doped regime .
@xcite this empirical correlation between the disappearance of the `` high - energy '' pseudogap and the decrease of the sc gap and , therefore , pairing strength suggested already several years ago that the high - energy features at ( @xmath131 ) ( which are correctly described by the `` high - temperature '' qmc simulations ) are closely related to the pairing interaction . @xcite
it is clear , however , that there is a strong need to perform calculations at much lower temperatures , i.e.  temperatures below @xmath135 . only then the `` low - energy '' sc or normal state pseudogap can be detected and
only then the question `` where do holes enter first ? '' can be answered correctly . according to the results of the `` high - temperature '' ( @xmath136 ) qmc simulations shown in fig .
[ qmc ] for @xmath132 , holes enter first around ( @xmath137 ) .
this can be understood by referring to the idea @xcite that a higher temperature @xmath3 effectively acts as an increased doping which destroys the magnetic brillouin zone .
this allows holes to first enter into the `` arc '' of single - particle excitations spanned around @xmath106 .     and
@xmath92 ( @xmath91 ) in the hole - doped ( a ) and in the electron - doped ( b ) case .
results are shown for dopings in the mixed af+sc state ( see sec .  [ gspd ] ) , i.e.  for @xmath138 ( a ) and @xmath139 ( b ) , respectively .
, title="fig : " ]   and @xmath92 ( @xmath91 ) in the hole - doped ( a ) and in the electron - doped ( b ) case .
results are shown for dopings in the mixed af+sc state ( see sec .  [ gspd ] ) , i.e.  for @xmath138 ( a ) and @xmath139 ( b ) , respectively .
, title="fig : " ]    however , a portion of the corresponding `` large '' fermi surface seems to disappear already at the lower temperature @xmath140 .
this is indicated in fig .
[ qmc](a ) by the downturn of the quasi - particle - like band between ( @xmath131 ) and ( @xmath137 ) . in the experiment ( fig .  [ qmc](b ) ) , this behavior in the under - doped regime has been interpreted as the opening of a pseudogap in the underlying fermi surface near the ( @xmath141 ) to ( @xmath137 ) line .
@xcite this opening is obviously intimately related to the above question , namely where doped holes first enter , i.e.  to the possibility of hole pockets developing at very small temperatures and low dopings around ( @xmath127 ) in the hole - doped case .
in contrast to the finite-@xmath3 qmc result , the corresponding vca calculations for @xmath4 ( see fig .  [ spectrtp0 ] ) show that holes first go into the coherent band around ( @xmath142 ) forming `` hole pockets '' consistent with experiments .
@xcite note that for the dopings considered in fig .
[ spectrtp0 ] , the system is in a mixed af+sc state with non - vanishing af order parameter ( see also fig .  [ phas - hole ] ) .
therefore , although the low - lying spectral weight is found around @xmath143 , the spectrum still shows an af gap at this wave vector .    related to this , it appears that even for @xmath124 the quasi - particle - like dispersion does not cross the fermi energy along the nodal direction ( between ( @xmath144 ) and ( @xmath145 ) in contrast to the experimental situation in fig .
although some ( weak ) accumulation of spectral weight close and above @xmath146 can be seen in fig .
[ spectrtp0 ] , we believe that this non - crossing is ( partly ) a shortcoming of our @xmath147 calculation and that the self - energy in the nodal direction requires larger clusters .    to address this issue for more realistic model parameters , i.e.  including a next - nearest - neighbor hopping @xmath148 , as well as to study the corresponding doping evolution in the electron - doped case , we discuss the vca results for @xmath126 shown in fig .
[ spectrelhol ] .
the calculations have been performed for the hole- ( a ) and for the electron - doped system ( b ) in the mixed af+sc phase ( see sec .
[ gspd ] ) , i.e.  for @xmath149 and @xmath139 , respectively .
let us concentrate on the very small ( @xmath150 ) hole doping first .
in agreement with the corresponding experiments @xcite in hole - doped cuprates , holes indeed first enter at ( @xmath127 ) .
although the system is in a mixed af+sc phase ( see sec .
[ gspd ] ) , the sc gap is zero at this nodal point , so that doping into nodal states , i.e.  normal metallic screening , apparently destabilizes the af solution already for rather low doping values ( @xmath151 ) .
this explains the phase diagram of fig .
[ phas - hole ] .
evolution of the low - energy spectrum upon hole doping .
parameters are as in fig .
[ spectrelhol ] .
the weight is obtained by integrating the low - energy @xmath152-spectrum down to @xmath153 below the fermi energy . ]     the same as in fig .
[ fermis - holdop ] but for the electron - doped case . ]    a similar picture can be inferred by looking at the evolution of the fermi surface as a function of doping .
this can be extracted from fig .
[ fermis - holdop ] where the low - energy spectral weight is plotted in the brillouin zone . in qualitative agreement with experiments , hole pockets start forming around @xmath143 for low doping , while a large fermi surface centered about @xmath128 starts building up at higher doping .    in contrast , in electron - doped systems , doped electrons initially form pockets around ( @xmath131 ) ( see figs .  [ spectrelhol](b ) and [ fermis - eldop ] , see also refs .  ) , in agreement with experiments .
@xcite fig .
[ fermis - eldop ] shows the spectral weight obtained by integrating @xmath154 down to @xmath153 below the fermi energy ( @xmath146 ) .
one has to be careful when interpreting fig .
[ fermis - eldop ] , since for @xmath155 there is still an af gap at the fermi surface near @xmath143 besides the sc one near @xmath123 .
therefore , the excitation spectrum is completely gapped , and one has to go away from @xmath146 in order to find some weight .
since we integrate only in a small energy window around the fermi level , the scale in fig .
[ fermis - eldop ] is an order of magnitude smaller than in the hole - doped case , fig .
[ fermis - holdop ] , for which the af gap near @xmath143 is shifted away from the fermi surface .
nevertheless , one finds that the lowest - lying states in the electron - doped case are around @xmath123 . here
, the density of states is large and provides a large `` reservoir '' for electron doping .
this , in combination with the fact that the chemical potential lies in the sc gap ( of the af+sc phase ) , stabilizes the af solution for a larger doping range than in the hole - doped case , allowing for the af gap to decrease more gradually ( see sec .
[ gspd ] ) . also in this case , the breakdown of the magnetic solution occurs as soon as the chemical potential reaches the bottom of the band at ( @xmath127 ) .
the observation that doping into ( @xmath127 ) generically makes the af phase unstable , suggests that also the occurrence of phase separation , to be discussed in sec .
[ gspd ] , is rather generic and quite independent of model details .    the low - energy spectral weight , as plotted for different dopings in fig .
[ fermis - eldop ] , agrees quite well with the arpes measurements of armitage et al .
@xcite ; doped electrons away from half - filling first form pockets around @xmath123 . at higher
doping , fermi - surface segments start developing along the @xmath156 line . finally , these segments connect and form a large fermi surface around @xmath128 near optimal doping .
in this section the ground - state phase diagram and the af to sc transition are presented and discussed ( see also refs .  ) . in particular ,
we will focus on the similarities and differences between the hole- and electron - doped system as well as their relation with the evolution of the single - particle spectrum .    in order to determine the @xmath4 phase diagram
, we proceed as described in sec .
[ sec : vca ] , i.e.  two symmetry - breaking terms ( weiss fields ) @xmath100 and @xmath101 are included and their respective strengths , @xmath104 and @xmath105 are treated as variational parameters in addition to @xmath71 . as discussed in sec .
[ sec : n ] , the use of the additional variational parameter @xmath71 is required in order to have a consistent determination of the particle density .    the phase diagram for the hubbard model with @xmath90 and next - nearest - neighbor hopping @xmath92 ( we set @xmath91 ) as obtained from our calculations , is shown in fig .
[ phas - hole ] for the hole - doped and in fig.[phas - elec ] for the electron - doped case . in the upper part of each figure ,
the chemical potential @xmath30 is plotted as a function of @xmath157 . in the corresponding lower parts
, we display the af @xmath158 and sc @xmath159 order parameter as a function of doping @xmath157 .
note that the order parameters @xmath114 and @xmath115 are different from the weiss fields @xmath104 and @xmath105 , respectively . quite generally , however , a nonvanishing stationary value for a weiss field produces a nonvanishing order parameter , respectively , although the latter can be much smaller .
antiferromagnetic and superconducting order parameters , @xmath114 and @xmath115 , and chemical potential @xmath30 as functions of hole doping @xmath157 .
@xmath115 and @xmath30 are plotted for the af+sc ( green , @xmath160 , @xmath161 ) as well as for the pure sc homogeneous solutions ( blue , @xmath162 , @xmath163 ) .
note that @xmath115 is scaled by a factor 5 for convenience . for @xmath164
the system exhibits a coexistence of af and d - wave sc order .
phase separation occurs between the doping levels @xmath165 and @xmath166 . for @xmath167 pure d - wave sc is realized . in the phase separation region @xmath168 ,
the homogeneous solution becomes unstable , and the system prefers to separate into a mixture of two densities corresponding to @xmath165 and @xmath166 .
the chemical potential @xmath169 is determined by the maxwell construction shown in the upper figure . at @xmath170 the slope of the af+sc solution changes sign . ]
same as fig .
[ phas - hole ] but for electron doping .
note the enhanced robustness of the af state and the strongly reduced scale @xmath171 as compared to hole doping . ]
let us discuss hole doping first ( see fig .
[ phas - hole ] ) . for low dopings @xmath157
we find a homogeneous symmetry - broken phase in which both the af as well as the sc order parameter @xmath114 and @xmath115 are non - zero .
this corresponds to a phase where af and sc order microscopically and coherently coexist .
a homogeneous phase with pure sc ( @xmath172 and @xmath173 ) is obtained for larger dopings .
the behavior of @xmath114 vs.  @xmath157 in fig .
[ phas - hole ] seems to suggest that the transition to the non - magnetic state is continuous ( second order ) as a function of doping .
however , a glance at the non - monotonous behavior of the chemical potential @xmath30 plotted as a function of @xmath157 in the upper part of the figure , indicates the occurrence of a charge instability .
the system tends to separate into a hole - poor ( @xmath165 ) and a hole - rich ( @xmath166 ) phase .
the two dopings @xmath165 and @xmath166 , as well as the chemical potential @xmath169 in the phase - separated region , are identified by the maxwell construction shown in the upper part of fig .
[ phas - hole ] .
@xmath170 is the point where the slope of @xmath174 changes sign .
of course , this phase separation is obtained within a treatment that , in the present work , just considers for homogeneous order parameters only and that neglects surface effects .
furthermore , the vca is a mean - field - type approach on length scales beyond the size of the individual cluster .
therefore , the above result has to be interpreted with care .
we expect it to signal a _
tendency _ towards the formation of microscopic inhomogeneities , such as stripes , checkerboard patterns , etc .
@xcite this tendency will be further studied in future work by considering larger clusters and/or by allowing for a more general variational solution which explicitly describes stripe inhomogeneities ( for example , by considering coupled clusters with different dopings , as in ref .  ) . in this case , one might expect the phase transition to become more continuous and phase separation to disappear eventually .
let us now discuss the electron - doped case and the similarities and differences as compared to hole - doping .
the first observation is that , while the phase diagrams in figs .
[ phas - hole ] and [ phas - elec ] are _ qualitatively _ similar , the phase in which long - range af order is realized is spreading to significantly larger doping values in the electron - doped case , in overall agreement with the experimental situation .
[ phas - elec ] shows that phase separation occurs in the electron - doped case as well , although the associated energy scale @xmath171 is smaller with respect to the hole - doped case by about an order of magnitude . in ref .
it is argued that this can be related to the different pseudogap and sc transition scales in hole- and electron - doped materials . this may give support to theories @xcite which are based on the notion that fluctuations of competing phases , or of the related order parameters , are responsible for the pseudogap phenomenon .
vs.  @xmath104 ( with the two other parameters @xmath105 and @xmath71 fixed at their stationary - point values ) for different @xmath30 in the hole - doped case .
parameters : @xmath90 , @xmath92 ( @xmath91 ) . ]    depending on the value of @xmath30 , there may be two solutions , an af+sc and a pure sc one , corresponding to two stationary points of @xmath39 . in order to gain insight into the first - order transition between these two solutions , it is instructive to observe the behavior of the grand potential @xmath39 as a function of the variational parameter @xmath104 associated with this transition for different @xmath30 . for simplicity , we consider here the hole - doped case only .
results are plotted in fig .
[ omegamin ] .
note that @xmath104 is varied , while the other two variational parameters , @xmath105 and @xmath71 , are fixed at their values at the respective stationary point . at low doping
, there is a single minimum of the grand potential @xmath39 at a finite value @xmath104 only ( fig .
[ omegamin ] , @xmath175 ) .
in contrast , the paramagnetic ( @xmath176 ) solution is given by a local maximum .
upon further doping , @xmath39 _ additionally _ develops a local minimum at vanishing af variational parameter @xmath176 ( fig .
[ omegamin ] , @xmath177 ) . for @xmath178 , the minimum at @xmath179 becomes lower than the one at finite @xmath104 indicating a first - order phase transition to a non - magnetic state ( fig .
[ omegamin ] , @xmath180 ) . eventually , the local maximum lying between the two minima merges with the minimum at finite @xmath104 and _ the af solution disappears _ ( fig .
[ omegamin ] , @xmath181 ) .
this disappearance just occurs when the chemical potential @xmath30 enters the quasi - particle band around @xmath143 ( @xmath182 ) .
although fig .
[ omegamin ] displays the results for the hole - doped case , this behavior occurs qualitatively in the electron - doped case as well , although with a much smaller characteristic energy scale as discussed above .
again , these results reflect a qualitatively similar behavior for electron and hole doping .
by means of a recently developed quantum - cluster approach , we have carried out a detailed analysis of the phase transition from the antiferromagnetic to the superconducting phase in the hubbard model at zero temperature .
the main results concerning the nature of the af to sc transition are summarized in figs .
[ phas - hole ] and [ phas - elec ] . at low dopings
the af phase actually mixes with a weak @xmath0-wave sc component .
a similar coexistence phase is observed , for example , in prcecuo , @xcite and was also obtained in previous mean - field and cluster calculations . @xcite upon further increasing doping , we find a transition to a pure @xmath0-wave sc phase .    the phase - separation scenario , which is found for hole but also for electron doping within our approximation ,
should carefully be interpreted as a general _ tendency _ of the system to form microscopically inhomogeneous phases .
indeed , by allowing for a more general spatial dependence of the order parameter , the macroscopic phase separation will probably be replaced by other microscopically inhomogeneous phases , such as stripes , checkerboard order , etc . in particular , this might be expected to be the case if long - range coulomb interaction is taken into account additionally .
@xcite the situation is subtle in the electron - doped case . here ,
in contrast to previous theoretical calculations , our results also suggest phase separation , although the corresponding energy scale ( see fig .  [ phas - elec ] ) is one order of magnitude smaller than in hole - doped compounds .    our results for the single - particle excitations presented in sec .
[ spe ] support the idea that it is , in fact , the single - particle spectrum which holds the key for understanding the qualitative differences seen in the ground - state phase diagram between the hole- and electron - doped cuprates .
this concerns , for example , the robustness of the af state in the electron - doped case . here
, electrons initially form pockets around ( @xmath131 ) in accordance with experiments . @xcite the fact that the density of states is large there , as well as the presence of a sc `` gap '' at this nodal point , stabilizes the af state for a larger doping regime as compared to the hole - doped case .
the breakdown of the magnetic solution appears as soon as `` normal '' metallic screening sets in , i.e.  when @xmath30 touches the band at ( @xmath142 ) .
@xcite    one should note that the variational cluster approach ( vca ) is able to treat the fluctuations correctly up to the range of the cluster size only .
therefore , the question arises whether the sc solution we ( and also others @xcite ) obtain within the af phase is a true long - range sc phase or whether it is only a signal of strong pairing fluctuations within the af phase leading to a sc pseudogap
. the latter hypothesis could be supported by the fact that results obtained with different cluster sizes @xcite seem to indicate a size dependence of the sc order parameter , and by the fact that the sc order parameter is about a factor three smaller in the af+sc phase than in the pure sc one . the presence or not of such a microscopic coexistence phase
may depend on material details .
certainly , our results suggest that the sc gap @xmath115 ( or pseudogap in the case of fluctuations ) is important in order to stabilize the af phase in electron - doped materials , as discussed in sec .
[ spe ] .
similar vca calculations have recently been carried out by snchal et al .
@xcite using clusters up to 10 sites but without the variation of the on - site energies .
these authors show that the single - band hubbard model is sufficient to explain the different overall shapes of the phase diagrams for hole- and electron - doped cuprates . however , their results seem to suggest that in the electron - doped case the af to sc transition is continuous and associated with a quantum - critical point , in contrast to our results .
this clearly shows the importance for a consistent determination of the average particle number .
the recent substantial progress in relating the `` high - energy '' physics of the hubbard model and its variants to the low - energy physics of competing phases is to a large extent due to the development of different quantum - cluster theories .
apart from the vca , these are the cluster extensions of the dmft , such as the cellular dmft ( c - dmft ) and the dynamical cluster approximation ( dca ) .
kyung et al.@xcite have shown that some of the important problems in strongly correlated electron systems may be studied highly accurately using comparatively small clusters .
maier et al .
@xcite performed a systematic cluster - size study of the @xmath5 hubbard model using rather large clusters ( up to @xmath183 sites ) .
converged results suggest a finite-@xmath3 instability to @xmath0-wave sc state .
because of the qmc minus - sign problem , however , results were limited to @xmath184 where the typical correlation energy @xmath15 and the magnetic energy scale of the htsc is not yet achieved . on the other hand ,
the present vca studies are clearly not yet converged with respect to the cluster size , as one can read off from , for example , the relatively large sc gap in the single - particle excitations , displayd in fig .
[ fermis - holdop ] .
an extension to larger clusters is , at least in principle , also possible within the vca .
this , however , necessarily implies the use of stochastic ( qmc ) methods as solvers for the cluster reference system .
we would like to thank d.  j.  scalapino for many useful discussions .
one of us ( wh ) would additionally like to acknowledge the hospitality of the kavli institute for theoretical physics in santa barbara ( supported by nsf grant no .
phy99 - 0794 ) .
the work was supported by the dfg forschergruppe 538 , by the konwihr supercomputing network in bavaria and by the fwf project n.  p18505-n16 .
the reasoning in sec .
[ sec : n ] can straightforwardly be generalized to the off - diagonal ( @xmath85 ) elements of the one - particle density matrix @xmath86 . in this case
, thermodynamic consistency means that , for a selected pair of sites @xmath185 , the derivative @xmath186 is equivalent with the integral @xmath187 where , for convenience , the hopping is formally assumed to be spin - dependent .    starting with eq .
( [ eq : hop1 ] ) , we note that the @xmath188 dependence of @xmath59 $ ] is due to the explicit @xmath188 dependence of the free green s function @xmath189 in eq .
( [ eq : ocalc ] ) , and due to the implicit @xmath188 dependence of @xmath68 : @xmath190}{\partial t_{ij\sigma } }     \nonumber \\    & = &    t \sum_{\omega_n } e^{i\omega_n 0^+ } \left (    \frac{1}{{{\bm g}}^{-1}_{0,{{\bm t}}}(i\omega_n ) - { { \bm \sigma}}_{{{\bm t}}'_{\rm s } , { { \bm u}}}(i\omega_n ) }     \right)_{ji\sigma }    \nonumber \\    & + &    \frac{\partial \omega_{{{\bm t } } , { { \bm u } } } [ { { \bm \sigma}}_{{{\bm t}}',{{\bm u } } } ] }    { \partial { { \bm t } } ' } \bigg|_{{{\bm t } } ' = { { \bm t}}'_{\rm s } } \cdot    \frac{\partial { { \bm t}}'_{\rm s}}{\partial t_{ij\sigma } } \ : .\end{aligned}\ ] ] because of the stationarity of the grand potential , the second term on the r.h.s .  can be ignored .
we are thus left with the first term which , after transforming the matsubara sum into an integration over real frequencies , exactly yields the desired expression ( [ eq : hop2 ] ) .    in case that not all elements of @xmath42 are treated as variational parameters , the derivative of @xmath39 in the second term on the r.h.s .
is a derivative with respect to those elements only .
if , for the given pair of sites @xmath13 and @xmath111 , @xmath191 is selected as a variational parameter , the contribution due to the implicit @xmath188 dependence vanishes again . if , on the other hand , @xmath191 is not treated as a variational parameter but fixed at @xmath192 from the beginning , the second term will give a finite contribution which spoils the thermodynamic consistency .
the c - dmft is obtained by choosing as a reference system disconnected clusters with a continuous non - interacting bath coupled to each cluster site . carrying out the @xmath42 partial derivatives in eq .
( [ eq : stat ] ) and using eq .  ( [ eq : ocalc ] ) for the sft grand potential , we get the euler equation in the form : @xmath193 here @xmath13 and @xmath111 run over the sites of the original lattice , excluding bath sites .
the one - particle bath parameters , namely the on - site energies of the bath sites and the hybridization of the bath sites with the correlated ( physical ) sites , have to be treated as ( a continuous set of ) variational parameters . in the c
- dmft it is assumed that bath parameters can be found such that the first factor in ( [ eq : euler ] ) vanishes : @xmath194 for arbitrary @xmath195 and for sites @xmath13 and @xmath111 belonging to the same cluster . note that , by construction of the reference system , @xmath196 and also @xmath197 vanish if @xmath13 and @xmath111 belong to different clusters . consequently ,
if bath parameters can be found such that eq .
( [ eq : cdmft ] ) holds , the euler equation ( [ eq : euler ] ) will be satisfied , too : the self - energy functional is stationary at the c - dmft self - energy .
( [ eq : cdmft ] ) is just the self - consistency equation of the c - dmft ( see refs .  ) .
it is easy to see that @xmath198 within the c - dmft : this simply follows by comparing eqs .
( [ eq : nn ] ) and ( [ eq : nnp ] ) for @xmath53 and @xmath81 and by using the c - dmft self - consistency equation ( [ eq : cdmft ] ) .
consider now the high - frequency expansions of the green s function and of the self - energy of the reference system .
using @xmath199 , @xmath200 , and @xmath201 for short , we have ( see ref .   for high - frequency expansions in the hubbard model ) :
@xmath202 here , @xmath203 is the average occupation in the reference system . using @xmath204 , inserting the expansions into the c - dmft self - consistency equation ( [ eq : cdmft ] ) , and expanding in powers of @xmath205 once more
, we immediately find : @xmath206 for ( correlated ) sites @xmath13 and @xmath111 within the same cluster , and , in particular , @xmath78 where @xmath207 and @xmath208 .
this means that , within the c - dmft , consistency with respect to the particle number is assured by setting the on - site energies within a cluster of the reference system to their `` physical '' values .
it goes without saying that a consistent determination of @xmath53 requires an _ exact _ solution of the reference system which , in case of the c - dmft , is by no means trivial .
an obvious idea is to replace the continuous bath by a few uncorrelated sites only to allow for an application of the lanczos method ( see ref .
, for an example ) .
this , however , immediately implies that the self - consistency equation ( [ eq : cdmft ] ) can not be fulfilled exactly any longer and , strictly speaking , the determination of the particle number becomes inconsistent _ unless _ @xmath71 is treated as a variational parameter within the sft framework .                  c. dahnken , m. aichhorn , w. hanke , e. arrigoni and m. potthoff , phys .
b * 70 * , 245110 ( 2004 ) ; c. dahnken , m. potthoff , e. arrigoni , and w. hanke , to appear in fizika nizkikh temperatur ( cond - mat/0504618 ) .                                                    n.  p. armitage , f.  ronning , d.  h. lu , c.  kim , a.  damascelli , k.  m. shen , d.  l. feng , h.  eisaki , z .- x .
shen , p.  k. mang , n.  kaneko , m.  greven , y.  onose , y.  taguchi , and y.  tokura , phys . rev .
* 88 * , 257001 2002 .
n.  p. armitage , d.  h. lu , c.  kim , a.  damascelli , k.  m. shen , f.  ronning , d.  l. feng , p.  bogdanov , z.  x. shen , y.  onose , y.  taguchi , y.  tokura , p.  k. mang , n.  kaneko , and m.  greven , phys .
* 87 * , 147003 2001 .