after almost two decade of intense research about the cuprates , there is still plenty of open questions in this problem .
however , it is recognized that the electrons which move on the @xmath2 planes are the most relevant to describe their physical properties @xcite . in the undoped regime , these compounds are insulators and exhibit antiferromagnetic order at sufficient low temperatures @xcite .
the physical properties of the insulating phase can be well described by the heisenberg model @xcite . upon doping , these systems
suppress the antiferromagnetic order and become superconductors . in this scenario
there is no doubt that the @xmath3-@xmath3 electron correlations play a fundamental role .
the study of the electronic structure near the fermi level @xmath4 in such strongly correlated systems is very important to understand their physical properties @xcite .
earlier angle - resolved photoemission experiments ( arpes ) have showed the presence of flat bands close to @xmath4 in a region centered around the point @xmath5 in the @xmath6-type cuprates like @xmath7@xmath8@xmath9 and @xmath10@xmath11@xmath12 @xcite . due to the presence of strong correlations , to study some physical properties of these cuprate compounds ,
the one - band hubbard model @xcite can be used .
_ @xcite have done monte carlo calculations in the one - band hubbard model .
their results show bands with a flat region near @xmath5 point for a given doping which agreed with the previously mentioned arpes results @xcite .
beenen and edwards @xcite , using the roth s two - pole approximation @xcite in the one - band hubbard model , have studied the normal state of the model obtaining flat quasi - particle bands , which agree well with those found with monte carlo simulations @xcite .
the roth s two - pole approximation has been proposed to improve the hubbard - i approximation @xcite by considering a decoupling scheme which produces an additional energy shift ( the roth s band shift ) in the peaks of the spectral function .
that result is in agreement with those obtained by harris and lange @xcite .
they looked at the moments of individual peaks in the spectral function .
the presence of the exchange term @xmath13 in the roth s band shift exhibits in it a spin dependence . as consequence ,
the roth s method raises the possibility of magnetic solution in the hubbard model while this feature is not present in the hubbard - i approximation .
recently , due to the good agreement between the roth s and the monte carlo data , beenen and edwards have extended the roth s two - pole approximation in order to investigate the superconducting properties of the one - band hubbard model .
their main achievement has been to show the emergence of the pairing with @xmath14 symmetry in a given amount of doping . in that approach , the gap equation for @xmath3-wave
symmetry depends on a particular four operator correlation function which , in principle , can be found extending the roth s formalism to obtain two particle green s functions .
however , the authors have introduced two decoupling schemes to calculate the gap .
the first one ( the factorization procedure ) has been formulated to treat the problem for intermediated values of @xmath15 ( the coulomb interaction ) and it provides an upper estimate for the gap and @xmath16 ( the critical temperature ) .
the second one is adapted for very large @xmath15 scenario which preserves the proper limit for @xmath17 where the gap function vanishes .
nevertheless , the one - band models neglect the presence of the oxygen sites . due to the strong correlations at the @xmath18-sites
, the oxygen sites may be occupied by holes when the system is doped @xcite . for instance
, the hubbard one - band model suffers some limitations to describe the low - energy physical properties of the cuprate superconductors @xcite . in the doped regime
, the one - band hubbard model gives a wrong description of various properties , like , for example , the asymmetric magnetic doping - temperature phase diagram @xcite .
therefore , a model which take into account also the oxygen can be more adequate to treat the cuprate systems in the doped regime @xcite .
this raises the question whether it is possible to extend the beenen and edwards analysis to investigate the @xmath3-wave symmetry superconductivity when the hybridization is present .
recently , the present authors have used the extended hubbard model @xcite with the roth s method to study the role of hybridization in the superconductivity following closely the approach introduced by beenen and edwards @xcite . as discussed in the references @xcite and @xcite , the flattening of the bands is directly related to the band shift .
the presence of flat bands at fermi level @xmath4 in the @xmath6-type cuprates @xcite suggests a high density of states at the fermi level , which can favor pair formation .
therefore , considering that the main responsible elements for the density of states at the fermi level are the @xmath3-electrons , it has been assumed in ref.@xcite that the @xmath19 pairs are the most relevant ones for superconductivity @xcite .
the band shift plays an important role in the study of the superconducting properties of the model using the roth s or some similar procedures . in reference
@xcite , the factorization procedure @xcite has been used to investigate the effects of the hybridization on the superconductivity .
it has been shown @xcite that the hybridization has strong effects in the shift and , therefore , in some superconducting physical properties such as the critical temperature @xmath16 . however , as a first approach , in reference @xcite , the band shift has been evaluated taking into account the hybridization effects , but disregarding temperature effects , superconducting properties and , most important , it has been considered in the limit @xmath20 . as a consequence of this limit , many correlation functions , which appear in the shift , are vanished .
the important point is that these correlation functions are very relevant in the sense of to include correctly the hybridization effects .
therefore , it would be necessary to calculate the shift with the @xmath15 finite in order to include the hybridization effects in a more complete way .    in this work
, the superconductivity problem has been studied using the roth s method , following closely reference @xcite , but adapted to the @xmath3-@xmath6 extended hubbard model . here , special attention is devoted to the effects of the hybridization and superconductivity in the band shift . in order to have the effects of the hybridization included properly in the superconductivity , the gap function
is obtained using the factorization procedure @xcite and the shift is evaluated with finite @xmath15 .
this procedure is justified because it preserves some correlation functions present in the band shift , which are non vanishing for finite @xmath15 . as consequence
, it captures the effects of the hybridization properly .
some preliminary results of this approach have been given in ref .
@xcite .
there are some shortcomings in the beenen and edwards approach @xcite . for instance , the @xmath14 pairing is quite dependent on the choice of the decoupling scheme for the correlation functions related to the gap .
however , in the present work , the main goal is to study the effects of hybridization .
therefore , as discussed in the previous paragraph , the natural choice is the decoupling scheme for intermediated @xmath15 , which is also the simplest one .
one is allowed to find in that procedure , at least , a better estimate for the gap ( and therefore for @xmath16 ) as a function of hybridization within the same decoupling procedure .
the paper presents the following organization . in section [ sec:2 ]
it is introduced the model and given a short introduction of the roth s method @xcite .
also , some analytic expressions for quasi - particle bands and the green s functions are derived . in section [ sec:3 ] , the factorization procedure proposed by beenen and edwards @xcite
is applied for the present case . in section [ sec:4 ] ,
the band shift is discussed in detail .
the numerical results are showed and discussed in section [ sec:5 ] .
finally , in section [ sec:6 ] , a short summary and some concluding remarks are given .
the model considered here assumes overlapping bands .
it is characterized by a narrow @xmath3-like band with a large density of states and a wide @xmath6-like band with low density of states .
the extended hubbard model is defined as : @xmath21 where @xmath22 is the chemical potential .
the @xmath23 , @xmath24 and @xmath25 , @xmath26 are the creation and annihilation operators of the @xmath3- and @xmath6-electrons , respectively , with spin @xmath27 at a lattice site @xmath28 .
the @xmath29 and @xmath30 are the centers of the on site energies of the occupied orbitals of the copper and oxygen respectively .
the second term of the hamiltonian given in eq .
( [ eq2.0 ] ) describes a narrow @xmath3-band with a hopping amplitude @xmath31 .
the hamiltonian ( [ eq2.0 ] ) considers also a @xmath6-band which is wider than the @xmath3-band .
the following relation between @xmath31 and @xmath32 can be established @xmath33 with @xmath34 .
also , @xmath35 to coincide the bottom of the @xmath3- and @xmath6-bands with the @xmath36 point @xmath37 as suggested by experimental results @xcite .
the third term corresponds to the coulomb interaction @xmath15 that represents the repulsion between two holes in the same @xmath3-orbital .
the last term of the hamiltonian ( [ eq2.0 ] ) is the @xmath0 hybridization and describes the nearest neighbor hopping process between the @xmath3-orbital of the cu - atom and the @xmath6-orbital of the o - atom .
considering a rectangular two dimensional lattice , the unperturbed @xmath3- and @xmath6-energy bands are given by @xmath38 and @xmath39 where @xmath40 is the lattice constant .    in this work , the hamiltonian given in eq .
( [ eq2.0 ] ) has been investigated using the roth s two - pole approximation @xcite to obtain the green s function in the zubarev s formalism . in the roth s procedure , a set of operators @xmath41
is introduced in order to describe the relevant one particle excitations of the system .
these operators satisfy in some approximation the following relation : @xmath42 _ { \left ( -\right ) } = \sum_{m}k_{nm}a_{m } \label{eq2.1}.\ ] ] anticommuting both sides of eq .
( [ eq2.1 ] ) with each operator of the set @xmath41 and taking the thermal average , the equation ( [ eq2.1 ] ) becomes : @xmath43 where @xmath44 and @xmath45 are the energy and normalization matrices , given by @xmath46 _ { \left ( -\right ) } , a_{m}^{\dag}\right ] _ { \left ( + \right ) } \right\rangle \label{eq2.3}\ ] ] and @xmath47_{\left ( + \right ) } \rangle \label{eq2.4}.\ ] ] in matrix notation , eq .
( [ eq2.2 ] ) is written as @xmath48 , where , if @xmath49 is nonsingular , then the @xmath50 matrix can be obtained .
with the equation of motion ( in the zubarev s formalism ) of the green s function @xmath51 and the eqs .
( [ eq2.1])-([eq2.4 ] ) , it is possible to obtain the following general green s functions @xmath52_{\left ( + \right ) } \rangle \label{eq2.6}.\ ] ] in the particular case , where @xmath53 , the elements of the green s function matrix @xmath54 are given by eq .
( [ eq2.5 ] ) .
thus , using the matrices @xmath55 and @xmath49 , the matrix @xmath54 is given by @xmath56 where @xmath57    considering the fact that the operators of the set @xmath58 describe the particle excitations of the system , the choice of these operators is very relevant to study the physical properties of the system . in order to discuss superconductivity , beenen and edwards , in their approach with the one - band hubbard model ,
mixed electron and hole operators and evaluated anomalous correlation functions @xcite .
therefore , using a set of four operators @xmath59 @xmath60 , it has been obtained a four - pole approximation to the green s functions . however , in order to discuss the role of the hybridization , it is necessary to adapt the formalism to include a @xmath6-operator in the original set of operators used by beenen and edwards .
thus , the new set of operators is given by @xmath61    in the present work , only the singlet pairing is considered , and particularly the _ d_-wave symmetry . in this particular case , @xmath62 and @xmath63 , where @xmath64 are the nearest neighbors of @xmath28 . using the set of operators given in eq .
( [ eq2.9 ] ) , and introducing the symmetries discussed above , the elements of the energy matrix defined in eq .
( [ eq2.3 ] ) can be obtained as : @xmath65\nonumber \\ \label{eq2.16}\end{aligned}\ ] ] where @xmath66 and @xmath67 are the fourier transform of @xmath68 and @xmath69 respectively .
the matrix @xmath70 present in the energy matrix @xmath71 is given by : @xmath72 % \nonumber \\ \label{eq2.10}\end{aligned}\ ] ] where @xmath73 .
it is assumed that the system considered here is translationally invariant , then @xmath74 .
the quantity @xmath75 is the fourier transform of @xmath76 in eq .
( [ eq2.11 ] ) , the band shift @xmath77 is defined as : @xmath78 where @xmath79 can be separated into two explicit contributions @xmath80 the terms @xmath81 and @xmath82 are associated with the hopping @xmath83 and the hybridization @xmath69 , respectively .
the hybridized term of @xmath79 may be written as @xmath84\delta _ { ij } \label{eq2.14},\ ] ] and the part associated to the hopping @xmath83 is given by @xmath85 the calculation of the correlation functions presented in eqs .
( [ eq2.14 ] ) and ( [ eq2.15 ] ) will be discussed in detail in section [ sec:4 ] .
one of the most important elements of the matrix @xmath71 is @xmath86 , where @xmath87 and @xmath88 the correlation function @xmath89 gives the gap of the superconductor state in the @xmath90wave case .
the elements of the normalization matrix @xmath91 are given from eq .
( [ eq2.4 ] ) as : @xmath92 and @xmath93 the remaining elements of the normalization matrix @xmath91 , due to the _ d_-wave symmetry and the anticommution rules , have been found to be zero .    using the energy and the normalization matrices @xmath71 and @xmath91 , respectively , the matrix green s function @xmath94 defined in eq .
( [ eq2.7 ] ) can be obtained . for simplicity ,
only the _ most relevant _ elements ( for the purposes of this work ) of this @xmath95 @xmath96 matrix are shown . following the roth s notation @xcite ,
the correlation function @xmath97 is related to the green s function @xmath98 as : @xmath99 where @xmath100 is the fermi function .
the chemical potential @xmath22 is obtained in the standard way , using the element @xmath101 of the matrix @xmath96 and the relation given in eq .
( [ eq2.21 ] ) .
the matrix element @xmath101 is given by @xmath102 } { \overline{d}\left ( \omega \right ) } , \label{eq2.22}\ ] ] where @xmath103 and @xmath104 are elements of the energy matrix @xmath105 , defined in eq .
( [ eq2.16 ] ) . in eq .
( [ eq2.22 ] ) , it is also necessary to introduce the following definitions : @xmath106 with @xmath107 @xmath108 and @xmath109 the quantities @xmath110 and @xmath111 are defined as @xmath112 and , @xmath113 in the particular case , when @xmath114 , @xmath115 and @xmath69 are zero , @xmath110 and @xmath111 represent the quasi - particle bands in the paramagnetic normal state of the one - band hubbard model . the term @xmath116 is given by : @xmath117 where @xmath118 is the fourier transform of @xmath77 given in eq .
( [ eq2.12 ] ) .
the denominator of the green s function @xmath119 given in eq .
( [ eq2.22 ] ) is defined as : @xmath120\nonumber \\ \label{eq2.30}\end{aligned}\ ] ] where @xmath121 with @xmath122\nonumber\\ & & \times[(\omega+e_{11})(\omega n_{-\sigma}^d+e_{22})-(\omega n_{-\sigma}^d+e_{12})^2].\nonumber\\ \label{eq2.32}\end{aligned}\ ] ] in eq .
( [ eq2.32 ] ) , @xmath123 and @xmath124 are elements of the energy matrix @xmath125 given in eq .
( [ eq2.10 ] ) .
the use of a set of five operators @xmath126 results in a five - pole approximation to the green s functions .
then , the @xmath127 defined in eq .
( [ eq2.30 ] ) may be also written as : @xmath128 where the quasi - particle bands @xmath129 ( with @xmath130 ) satisfy @xmath131 .
therefore , the resulting green s function can be written as a sum of five terms : @xmath132 where @xmath133 express the spectral weights which satisfy @xmath134
in the case of @xmath3-wave symmetry , the traditional correlation function @xmath135 is always zero .
therefore , this correlation function can not be used to determinate the pairing gap in the @xmath3-wave channel @xcite . in the factorization procedure proposed by beenen and edwards in ref .
@xcite , the correlation function given by eq .
( [ eq2.18 ] ) is rewritten as @xmath136\langle d_{i-\sigma}d_{l\sigma}\rangle \label{eq3.0}\ ] ] where the symmetry @xmath137 is conserved , and the products @xmath138 and @xmath139 are split up .
it is also introduced @xmath140 which allows to rewrite eq .
( [ eq3.0 ] ) as @xmath141 where @xmath142 can be calculated from @xmath143 .    considering the _ d_-wave symmetry ,
the fourier transform of @xmath144 given by eq .
( [ eq2.17 ] ) becomes @xmath145 \label{eq3.3}\ ] ] where @xmath146 is the gap - function amplitude .
due to the @xmath3-wave symmetry , @xmath147 for @xmath148 in the @xmath149 direction and @xmath150 when @xmath148 is in the @xmath151 direction .
the fourier transform of the correlation function @xmath152 is given by @xmath153 where @xmath154 is the number of sites in the system .
the correlation function @xmath155 can be evaluated using the green s function @xmath156 and the relation given by eq .
( [ eq2.21 ] ) .
the green s function @xmath156 can be rewritten as : @xmath157 where @xmath158 and @xmath127 is defined in eq .
( [ eq2.30 ] ) .    combining the equation ( [ eq2.17 ] ) with the eqs .
( [ eq3.2 ] ) to ( [ eq3.6 ] ) , the gap equation can be written as : @xmath159 where @xmath160 with @xmath161 ^ 2 f_{q\sigma } ^{13}(\omega ) \label{eq3.9}.\ ] ]
using the definition ( [ eq2.13 ] ) in eq .
( [ eq2.12 ] ) , the band shift @xmath162 can be written as : @xmath163 where @xmath164 + \lambda _ { ij\sigma } ^d}{n_{-\sigma } ^{d } ( 1-n_{-\sigma } ^{d } ) } \label{eq4.1}\ ] ] and @xmath165 the quantity @xmath166 is given by eq .
( [ eq2.14 ] ) .
the correlation function @xmath167 present in @xmath166 can be obtained from de green s function @xmath168v_k^{pd } } { \overline{d}\left ( \omega \right ) } \label{eq4.3}.\ ] ] the remaining correlation function @xmath169 present in @xmath166 is calculated from the green s function @xmath170v_k^{dp } } { \overline{d}\left ( \omega \right ) } \label{eq4.4},\ ] ] where @xmath171 with @xmath172 defined in eq .
( [ eq2.23 ] ) .
the quantity @xmath173 , in terms of the elements of the energy matrix ( [ eq2.16 ] ) , is given by : @xmath174 the green s function @xmath175 tends to zero as @xmath20 , consequently , the correlation function @xmath176 also vanishes recovering the result of ref .
@xcite for @xmath177 .
the quantity @xmath178 present in eq .
( [ eq4.1 ] ) is given by eq .
( [ eq2.15 ] ) .
the fourier transform of @xmath179 is given by : @xmath180 substituting eq .
( [ eq2.15 ] ) into eq .
( [ eq4.1 ] ) and then putting the result into eq .
( [ eq4.7 ] ) , the fourier transform of @xmath181 can be written as : @xmath182 the correlation functions present in @xmath183 are evaluated following the original roth s procedure @xcite . introducing extra operators @xmath184 ,
the correlation functions of the form @xmath185 can be calculated by using eqs .
( [ eq2.6 ] ) and ( [ eq2.21 ] ) . in refs .
@xcite , the sum present in eq .
( [ eq2.6 ] ) has been considered only over the operators which describe the normal state of the system . in the present work ,
the sum includes also the hole operators which describe the superconducting properties of the system .
thus , @xmath183 is given by : @xmath186 where the term @xmath187 is directly related to the gap function @xmath188 through the green s functions @xmath189 and @xmath190 ( see appendix a ) . the quantities @xmath191 , @xmath192 and @xmath187
are given in appendix a.
in this section , the numerical results obtained in this work are presented .
one of the most important parameters of the model given in eq .
( [ eq2.0 ] ) is the @xmath0 hybridization @xcite , which is defined as @xmath193 \label{eq5.0}.\ ] ] in this work , as in ref .
@xcite , the hybridization has been assumed @xmath194-independent @xmath195 , where @xmath196 is the average over the brillouin zone . in ref .
@xcite sengupta and ghatak have also used a @xmath194-independent hybridization due to the fact that the pairs occur within a small energy interval around the fermi level , therefore the dispersion of the hybridization can be neglected .
the total occupation number is given by @xmath197 , where @xmath198 is obtained combining @xmath143 ( eq . ( [ eq2.22 ] ) ) and the relation given in eq .
( [ eq2.21 ] ) .
the charge transfer energy @xmath199 is positive .
this means that the first hole added to the system will energetically prefer to occupy the @xmath3-orbital of the copper ions @xcite .
all results presented in this section are obtained with @xmath200 and @xmath201 .
consequently , @xmath202 , as estimated in ref .
@xcite .
llllll & @xmath15 & @xmath188 & @xmath203 & @xmath204 & @xmath187 + beenen and + edwards @xcite & @xmath205 & 0 & 0 & 0 & 0 +   + ref.@xcite & @xmath206 & 0 & 0 & @xmath205 & 0 +   + present work & @xmath205 & @xmath205 & @xmath205 & @xmath205 & @xmath205 +    as discussed in refs .
@xcite , the band shift @xmath207 ( see eq .
( [ eq4.0 ] ) ) can be evaluated considering different approximations . in the limit @xmath20 , some terms of the band shift vanish ( see ref .
@xcite ) . in ref .
@xcite , the present authors estimate @xmath207 in the limit of @xmath20 , but with finite @xmath15 in other parts of the problem . in ref .
@xcite , beenen and edwards evaluated @xmath207 in the normal state ( where @xmath208 ) and considering @xmath203 equal to zero and finite @xmath15 using the one - band hubbard model ( hybridization null ) . in the present work , the correlation functions present in @xmath207 given in eq .
( [ eq4.8 ] ) are evaluated following closely the procedure used by roth in ref .
@xcite . nevertheless , here , also the hole operators given in the set of eq .
( [ eq2.9 ] ) are used to evaluate the correlation functions . as consequence , a new term ( @xmath187 ) appears in @xmath207
( [ eq4.9 ] ) ) .
the approximations used to evaluate @xmath207 are shown in table [ tab:1 ] . in figure [ fig1](a ) , the quasi - particle bands @xmath129 , with @xmath209 ( see eq .
( [ eq2.33 ] ) ) , are plotted along the symmetry lines @xmath210 , in the two - dimensional brillouin zone .
the quasi - particle energies @xmath129 , in the superconducting state , are relative to the chemical potential @xmath22 .
the circles bf show the @xmath211 band , where the center of @xmath211 is shifted by @xmath212 relative to the zero of energy .
all results shown in this paper are obtained with @xmath213 .
the dashed line corresponds to the noninteracting @xmath214 band @xmath215 relative to the noninteracting chemical potential .
the figure [ fig1](b ) shows the superconducting gap between the electron and hole bands in the neighborhood of the @xmath5 point , while on the @xmath216 diagonal ( fig . [ fig1](a ) ) the gap is zero .
this fact reflects the @xmath3-wave symmetry proposed in this work .
the dashed lines show the absence of the gap in the normal state . in figure [ fig1](c ) ,
the region near to the @xmath217 point shows the gaps produced by the @xmath0 hybridization @xmath204 .
the dashed lines show the result for @xmath218 . in figure
[ fig1.1 ] , the electron and hole quasi - particle bands are shown for two different values of hybridization . as can be observed , the hybridization shifts the quasi - particle bands to lower energy by breaking the symmetry in relation to the @xmath194 axis .
the figure [ fig2 ] shows the spectral weights @xmath133 for two different hybridizations .
the dashed line corresponds to the sum of the five spectral weights which is equal to one ( see eq .
( [ eq10.4 ] ) ) . in figure [ fig2](b ) , the effects of the hybridization on the spectral weights are shown .
such effects cause a small change in the chemical potential and consequently in the superconductivity .
the figure [ fig3 ] shows the behavior of gap function amplitude @xmath219 as a function of the hybridization @xmath204 .
it is clear that there is a decreasing of @xmath219 with increasing @xmath204 .
the analysis of the function @xmath220 introduced in eq .
( [ eq3.81 ] ) and defined in eq .
( [ eq3.9 ] ) is important to understand the behavior of the gap function amplitude showed in figure [ fig3 ] .
the figure [ fig4 ] shows the function @xmath220 for @xmath221 and two different values of hybridization .
as can be seen in the dashed line , the magnitude of the function @xmath220 decreases when the hybridization is enhanced .
moreover , the function is shifted to lower energy , breaking the symmetry respect to @xmath222 . the symmetry break has been also observed in figure [ fig1.1 ] for the electron and hole quasi - particle bands . for @xmath221 , the product @xmath223 given in eq .
( [ eq3.81 ] ) vanishes when @xmath224 .
that is because the fermi function @xmath100 is zero for that range of @xmath225 . as consequence of the shift and the suppression of @xmath220 , the value of @xmath226 , which is given by the integral in eq .
( [ eq3.81 ] ) , decreases when the hybridization increase .
however , from eq .
( [ eq3.8 ] ) , it is necessary a minimum value for @xmath226 to obtain a nonzero solution for @xmath227 .
but , for very strong values of hybridization , the minimum value for @xmath226 is not reached and only the zero solution exists .    according to this analysis , there is a critical value of hybridization ( @xmath228 ) , above which , the superconductivity is suppressed .
similar results , which show a critical value for the hybridization , were also obtained in ref .
@xcite , for a @xmath194-dependent hybridization and using the hartree - fock approximation for the electron - electron interaction . in refs .
@xcite , although the high @xmath229 was not considered , the hybridization effects play an important role for resonant states .
the discussion above is also valid if the values of the temperature @xmath203 are raised with @xmath230 constant .
the only difference is that in this case the fermi function becomes sloping smoothly , changing the product @xmath223 .
the effect of the temperature in the fermi function causes a decreasing of @xmath226 and consequently of @xmath229 .    since the hybridization is directly related to the applied pressure @xcite , the transition temperature @xmath229 may have a dependence on pressure through the hybridization
however , the pressure dependence of @xmath229 is very complicated in high temperature superconductors . as discussed in ref .
@xcite , at least , in conventional superconductivity where the electron pairing is mediated by phonons , two effects are responsible for the pressure dependence of @xmath229 . the first one is related to the lattice vibrations , while the second one comes from the electronic contribution .
as long the pressure is increased , the lattice vibrations tend to increase @xmath229 , whereas the effects of the electronic contribution associated with the hybridization cause a decreases of @xmath229 .
the figure [ fig5 ] shows the function amplitude @xmath219 , as a function of temperature @xmath203 for several values of hybridization .
if the hybridization is enhanced , the critical temperature @xmath229 decreases .
therefore , the present results agree with the discussion above in the scenario where the electronic effects dominate .
this behavior for @xmath229 is also shown in figure [ fig6 ] , in phase diagrams displaying @xmath229 versus the total occupation number @xmath231 .
the numerical results obtained show that there is a critical value of hybridization where the superconducting phase vanishes .
these results agree with the ones obtained by the authors of ref .
@xcite for heavy - fermion superconductivity with an @xmath232-boson treatment .
the solid lines in figures [ fig6](a)-(b ) show the present result for @xmath229 , where the effects of the temperature , the superconductivity and the coulomb interaction have been included in the calculation of the band shift @xmath207 .
the dashed lines correspond to the results obtained in ref .
@xcite , where the band shift has been evaluated considering @xmath221 , @xmath208 and @xmath20 ( see table [ tab:1 ] ) .
the difference between the results can be explained by the analysis of the eq .
( [ eq2.14 ] ) where some of the correlation functions present in the band shift vanish in the @xmath17 limit . in eq .
( [ eq2.14 ] ) , the correlation function @xmath233 , which is directly related to the hybridization effects in the band shift , vanishes for @xmath17 .
it is important to highlight that the correlation functions in eq .
( [ eq2.14 ] ) are both negative .
therefore , for large @xmath15 , the correlation function @xmath234 decreases and the hybridized shift @xmath235 is enhanced .
however , for intermediate values of @xmath15 , both correlation functions remain finite . as consequence , the hybridization effects in the band shift and , therefore , in the superconductivity
, are weakened .
the figures [ fig6](c)-(d ) show the present results when the value of @xmath15 is increased .
the main consequence is , within the factorization procedure , to shift the window of doping where superconductivity is found , as in ref @xcite .    in figure [ fig7](a ) , the chemical potential is show as a function of the total occupation number @xmath231 for @xmath236 and two different hybridizations . in ref .
@xcite , the authors criticized the roth s method because the compressibility @xmath237 is negative in the vicinity of half - filling in the beenen and edwards result . in ref .
@xcite , by using a composite operator approach and imposing the pauli principle , the authors have showed that the compressibility remains negative .
however , they also showed that the pairing decreases the strength of the negative compressibility .    in the present work ,
a careful study about the nature of the negative compressibility and the effect of the hybridization near half - filling in roth s approximation has been carried out .
it has been verified that the most important contribution to provide negative compressibility comes from the spin - term @xmath238 present in the @xmath3-part @xmath239 of the roth s band shift @xmath207 ( see eqs .
[ eq4.8 ] and [ apx26 ] ) . in reference
@xcite , it has been showed that the correlation function @xmath238 plays an important role on the flattening of the quasi - particle bands . the correlation function @xmath238 increases with occupation and its effect is pronounced near half - filling . nevertheless , when the hybridization is present , the numerical results show that it acts in the sense of suppressing the negative compressibility near half - filling . because the hybridization considered here is @xmath194-independent @xcite , the hybridization term @xmath235 of the band shift is constant within the brillouin zone .
its main effect is to shift the poles of the green s functions and consequently to change the value of the chemical potential suppressing the negative compressibility . in figure [ fig7](a ) , it is clear that the effect of the hybridization in the chemical potential decreases the negative compressibility .
the figure [ fig7](b ) shows the gap function amplitude @xmath219 as a function of the total occupation number .
this result agrees with those obtained in figure [ fig3 ] , where @xmath219 decreases with increasing of @xmath230 .
in this work , the roth s two - pole approximation is extended to study the superconducting properties of the extended hubbard model given in eq .
( [ eq2.0 ] ) .
the quality of the roth s two - pole approximation had been investigated in a previous work by beenen and edwards @xcite . in their work
, they showed the remarkable agreement between the roth s and the monte carlo results @xcite for the one - band hubbard model in the paramagnetic normal state .
moreover , the flat bands obtained with roth s procedure show a qualitative agreement with the arpes experiment data @xcite in cuprates .
it is important to point out that the flattening observed in the quasi - particle bands which produces a peak in the density of states , can be connected with the van hove scenario . in cuprate systems
the van hove singularity is present in the vicinity of the fermi energy .
therefore , it is believed that the van hove scenario play a fundamental role in order to clarify the mechanism which drives the transition to superconductivity in these interesting materials @xcite .
the accuracy of the roth s results is very related to the adequate evaluation of the band shift .
therefore , the focus of the present work has been to evaluate the roth s band shift taking into account relevant effects as coulomb interaction , temperature , superconductivity and hybridization .
also , the effect of the hybridization in the superconducting of the model has been studied .
this work has been carried out following the factorization procedure proposed by beenen and edwards @xcite . in order to study superconductivity , beenen and edwards proposed to include hole operators in the original set of operators that describes the normal state of the system .
these operators can introduce the pairing formation in the @xmath3-band . the factorization procedure proposed by beenen and edwards @xcite and the @xmath3-wave symmetry are considered to obtain the gap function amplitude .
the hybridization effects are considered by also including a @xmath6-operator .
thus , the set of operators is enlarged to five , which results in a five - pole approximation to the green s functions .
the hybridization effects present in the band shift come from some correlation functions .
the important point is that part of them vanish when @xmath17 , as it have been done in ref .
@xcite . in order to consider properly the hybridization effects
, the band shift should be obtained for finite @xmath15 .
in fact , the obtained phase diagrams show that the presence of superconducting order exists in a larger range of doping when compared with the @xmath17 limit @xcite , for the same hybridization .
therefore , this result suggests that , in the @xmath17 limit , the hybridization effects are overestimated .
that is the ultimate justification for the use of the factorization procedure @xcite , which is valid for intermediated values of @xmath15 for the gap function .
the beenen and edwards s @xcite results are recovered taking @xmath218 in the present work .
the hybridization @xmath204 breaks the symmetry between the electron and hole quasi - particle bands , respect to @xmath194 axis .
also , the gap amplitude function @xmath219 and the critical temperature @xmath229 are suppressed with increasing the hybridization @xmath204 .
the results show that the chemical potential does not change significantly away the half - filling .
however , near half - filling , it is showed that the negative compressibility decreases with increasing @xmath204 .
the correlation functions present in the @xmath3-part of the band shift @xmath207 were discussed in detail . when the hole operators are also considered to obtain this correlation functions
, a new term appears in the @xmath3-part of the band shift @xmath207 .
the new term is directly associated with the superconducting properties of the system .
nevertheless , this term is quite small and therefore may be disregarding in the calculation of the band shift .
the authors are grateful to the grupo de fsica estatstica - ifm , universidade federal de pelotas where part of numerical calculations were performed .
this work was partially supported by the brazilian agencies capes ( coordenao de aperfeioamento de pessoal de nvel superior ) and fapergs ( fundao de amparo  pesquisa do rio grande do sul ) .
the correlation functions present in the band shift @xmath239 can be evaluated by introducing extra @xmath240 operators , as in the original roth s procedure . combining the eq .
( [ eq2.6 ] ) and the relation given in eq .
( [ eq2.21 ] ) , it is possible to write @xmath241_{(+)}\rangle \label{apx1},\ ] ] where @xmath126 and @xmath242 are members of the set of operators given in eq .
( [ eq2.9 ] ) . for evaluate @xmath243
, it has been necessary to introduce the following @xmath240 operators : @xmath244 and @xmath245    by considering the operator given in the eq .
( [ apx2 ] ) , the correlation function @xmath246 can be written as : @xmath247 where the right side of eq .
( [ apx4 ] ) may be obtained using the relation given by eq .
( [ apx1 ] ) .
therefore , it is necessary to evaluate the anticommutators @xmath248_{(+)}$ ] for the set of operators @xmath249 given in eq .
( [ eq2.9 ] ) . for @xmath250 ,
the @xmath251 operators are given by : @xmath252 @xmath253 @xmath254 @xmath255 and @xmath256 thus , the following results have been obtained @xmath257_{(+)}\rangle = n_{0\sigma}^{d } -e^{i\vec{k}\cdot \vec{r}_j}\langle d_{0\sigma}^{\dagger}d_{j\sigma}\rangle \label{apx9.1},\ ] ] @xmath258_{(+ ) }
\rangle=\langle n_{0-\sigma } ^{d}n_{j\sigma}^{d}\rangle -e^{i\vec{k}\cdot \vec{r}_j}\langle d_{0\sigma}^{\dagger}n_{j-\sigma}^{d}d_{j\sigma}\rangle \label{apx10},\ ] ] @xmath259_{(+)}\rangle=0 \label{apx11},\ ] ] @xmath260_{(+)}\rangle= -\langle n_{j\sigma } ^{d}d_{0\sigma}^{\dagger } d_{0-\sigma}^{\dagger}\rangle \label{apx12},\ ] ] @xmath261_{(+)}\rangle=0 \label{apx13}\ ] ] where , it has been assumed that the brackets are real and unchanged when the indices 0 and @xmath262 are interchanged . also ,
due to translational invariance of the system , @xmath263 .
considering the relations given by eqs .
( [ apx1 ] ) and ( [ apx4 ] ) with the results from eq .
( [ apx9.1 ] ) to eq .
( [ apx13 ] ) , the correlation function @xmath246 can be written as : @xmath264 where @xmath265 . in eq .
( [ apx14 ] ) , it has been introduced the following definitions : @xmath266 @xmath267 @xmath268 @xmath269 and @xmath270 where @xmath119 is given in eq .
( [ eq2.22 ] ) .
the remaining green s functions @xmath271 and @xmath190 are given respectively by @xmath272 } { \overline{d}\left ( \omega \right ) } \label{apx19.1}\ ] ] and @xmath273 where @xmath274 is defined in eq .
( [ eq4.5 ] ) and @xmath275 in eq .
( [ eq2.33 ] ) .
it is also necessary to define @xmath276 @xmath277 and @xmath278 where @xmath189 is given in eq .
( [ eq3.6 ] ) .    the correlation function @xmath279 present in eq .
( [ apx14 ] ) , can be obtained by repeating the procedure above using the operator @xmath280 ( given by eq .
( [ apx3 ] ) ) .
thus , @xmath281 where @xmath282 @xmath283 and @xmath284 with @xmath285 reversing the spin labels i.e. , @xmath286 in eq .
( [ apx19.6 ] ) and substituting the result into eq .
( [ apx14 ] ) , then @xmath287 \label{apx23}.\end{aligned}\ ] ]    for evaluate the two last correlation functions present in eq .
( [ eq4.8 ] ) , the following operators have been introduced @xmath288 and @xmath289 using @xmath290 , and following the procedure outlined above , the correlation function @xmath291 is given by @xmath292 .
\nonumber \\ \label{apx26}\end{aligned}\ ] ] similarly , using @xmath293 @xmath294 where the @xmath3-wave symmetry has been considered , therefore , @xmath295 .
the four @xmath296 operators introduced up to now are exactly the same operators used by roth in ref .
@xcite to obtain the band shift @xmath297 in the normal state and without hybridization .
however , in the present work , due to the presence of the hole operators ( see eq .
( [ eq2.9 ] ) ) , a new @xmath240 operator , which is given by @xmath298 has been introduced . with this operator , the correlation function @xmath299 present in eq .
( [ apx27 ] ) may be evaluated .
thus ,        the result given in eq .
( [ apx29 ] ) can be used in eq .
( [ apx23 ] ) to obtain @xmath303 \label{apx31}.\end{aligned}\ ] ] finally , with the results ( [ apx29 ] ) , ( [ apx30 ] ) and ( [ apx31 ] ) into eq .
( [ eq4.9 ] ) , the following result has been obtained @xmath304 where @xmath305 @xmath306 @xmath307}{1-\beta_{\sigma}\beta_{-\sigma}}\nonumber\\ & & -\frac{\alpha_{j\sigma}^{(1)}n_{0j-\sigma}^{(1 ) } + \beta_{j\sigma}^{(1 ) }   m_{j-\sigma}^{(1)}}{1+\beta_{\sigma } } \label{apx35}\end{aligned}\ ] ] with @xmath308