the unconventional temperature scaling of the in plane resistivity in high-@xmath3 cuprates has been of considerable interest for many years .
it is now widely accepted that overdoped cuprates exhibit a @xmath4-linear resistivity at high temperature which crosses over to a @xmath5 dependence at lower temperature [ ] .
the onset of the linear temperature dependence appears to coincide with the onset of superconductivity and the crossover temperature to a conventional @xmath5 dependence drops monotonically with the doping level to zero , roughly around optimal doping [ ] .
hussey and collaborators measured the resistivity on a number of cuprates particularly the single layer tl2201 in the overdoped region [ ] and lsco into the underdoped region [ ] . very recently the barii et al .
[ ] reported results for the resistivity for more cuprates , e.g. hg1201 and ybco over a wider density range .
they found that when the resistivity was normalized to a give a resistance per cu@xmath2o@xmath2 plaquette , @xmath6 , universal results for all cuprates followed with identical coefficients .
further they argued for a universal trend in both the linear @xmath7 and quadratic @xmath8 regimes , with the coefficients @xmath9 and @xmath10 proportional to @xmath11 , the inverse hole density .
as barii et al .
[ ] remarked , if one uses a drude formula for resistance , this behavior is indicative of a carrier density that is proportional to the density of holes rather than electrons .
such behavior is consistent with a doped mott insulator scenario .
anderson [ ] has argued that the cuprates as doped mott insulators should be in the strong coupling regime of the hubbard model at all hole densities .
he has put forward the hidden fermi liquid [ hfl ] ansatz to describe the connection between the fermi surface before gutzwiller projection and the electronic state that results from the removal of doubly occupied states by gutzwiller projection .
this hfl ansatz gives a resistivity with the observed temperature dependence , crossing over from a linear to a quadratic form as the temperature is lowered .
an earlier viewpoint interprets the cuprate phase diagram in terms of an underlying quantum critical point separating the overdoped and underdoped regions . in these theories
the effective interactions between electrons are mediated and controlled by the order parameter fluctuations associated with the quantum critical point [ ] .
another approach recently put forward by kokalj and coworkers [ ] is based on an anisotropic marginal fermi liquid self energy .
their phenomenological theory gives good fits to the resistivity and also to the angular dependent magnetoresistance [ admr ] results [ ] and a number of other properties .        in an earlier study , which employed a functional renormalization group [ frg ] treatment of the two - dimensional hubbard model to describe the overdoped region , ossadnik et al .
[ ] found good qualitative agreement with the experimental results on overdoped tl2201 for the temperature and angular dependence of the relaxation rate determined from the magnetoresistance [ ] . since overdoped tl2201 shows a conventional full fermi surface in angle - resolved photoemission spectroscopy ( arpes ) at zero magnetic field [ ] and in quantum oscillation experiments at finite magnetic fields [ ] a perturbative treatment is of interest .
note that the frg calculations are evaluated only to one - loop order and therefore are quantitatively reliable only for weak to moderate values of the hubbard onsite interaction , @xmath12 . as a result ,
a quantitative comparison to the experimental data shows a substantial underestimate of the magnitude of the calculated temperature dependent resistance , e.g. see the comparison in kokalj and mckenzie [ ] .    in their calculation of unusual transport behavior , ossadnik and
collaborators analyzed the quasiparticle scattering vertex for overdoped cuprates within a functional renormalization group ( frg ) approach [ ] .
they found an anisotropic term in the quasiparticle scattering rate with an angular form similar to the experiment and with an unconventional @xmath4-linear dependence which increased when approaching optimal doping starting from the overdoped side .
they argued that this unusual term appearing in the imaginary part in the quasiparticle self - energy could lead also to an equivalent temperature dependence for the in - plane resistivity .
the goal of this study is to expand this idea and move beyond a self energy analysis of the transport properties to a full calculation of transport properties using solutions of the boltzmann equation .
this is possible within a frg approach since the method computes the strongly renormalized quasiparticle interactions .
we determine the matrix elements for the scattering events within the frg framework presented in [ ] and use numerical techniques to solve the linearized boltzmann equation .
a detailed numerical treatment is required in view of the strong anisotropy in the transport life times [ ] due to both the anisotropic quasiparticle interactions and the strongly angle dependent fermi velocities introduced by the anisotropy of the band structure .
the importance of full transport calculations was stressed by hlubina and rice in an early paper [ ] .    as will be discussed later ,
our calculations agree qualitatively with the unconventional temperature dependence but the magnitude of the resistance is substantially lower .
this disagreement in the absolute values is in line with the restrictions on @xmath12 required to justify the one - loop approximation in the frg .
another discrepancy between our numerical results and the experiment is the lack of scaling of the resistance coefficients , @xmath13 and @xmath14 , with @xmath11 the inverse hole density as reported by barii et al .
note however , our results are for the overdoped samples . the data in fig . 5 of ref .
[ ] deviate from the inverse doping scaling in the overdoping density range . in the case of the linear coefficient ,
@xmath13 the data show a linear decrease and vanish at the onset of superconductivity in line with the earlier results of hussey and coworkers [ ] .
a good example of such an onset of a non - fermi - liquid - like @xmath4-linear contribution to the resistivity can be found in the resistivity measurements by hussey and coworkers [ ] upon reducing the hole concentration from the overdoped regime toward optimal doping in lsco  ( cf .
[ fig : experimental resistivity ] ) .
they observed a strong growth of the temperature dependent contribution to the resistivity with a decrease in the hole doping together with a remarkable trend toward linear resistivity that goes down to lowest temperatures .
they used strong magnetic fields to explore the normal state transport properties within the superconducting dome .
the analysis of the temperature dependence of the resistivity was based on a second - order polynomial fit to the experimental data over a wide temperature range : @xmath15 .
this allows to distinguish doping regimes with dominant linear versus dominant quadratic temperature dependence .
we also solve the boltzmann equation to obtain the seebeck coefficient which describes the charge transport in response to a thermal gradient .
standard fermi liquid theory predicts a metallic seebeck coefficient with a linear @xmath4-dependence .
the experiments by lalibert and coworkers [ ] showed substantial deviations from the standard form with a complex density dependence .
again our calculations reproduce the qualitative trends in the temperature and density dependence very well .    in sec . [
sec : numerical results ] we show our numerical results and compare with the experimental data of refs .
we also examine the validity of matthiessen s rule which relies on a decoupling of impurity and two - particle scattering in sec . [
sec : momentum relaxation mechanism in the cuprates and matthiessens rule ] . finally in sec .
[ sec : the seebeck effect in overdoped lsco ] , we study the thermoelectric effect of lsco , use the seebeck coefficient as a sign of non fermi liquid transport and examine the evidence for quantum critical behavior .    the numerical method that we have developed to study the transport properties of correlated materials such as lsco  is explained in detail in the appendix , app .
[ sec : model for the quasiparticle dispersion of lsco]-[sec : effective model for the normal state charge transport of the high - temperature superconductor lsco ] .
we solve the semiclassical boltzmann equation taking the full angular and energy dependence of the distribution function into account .
the two important inputs of the model are the explicit form of the band structure and the temperature dependent calculation of the quasiparticle scattering rates .
we begin by reviewing our numerical results .
we use a semiclassical model to describe the normal state charge transport of overdoped lsco . for the quasiparticle dispersion we used a two - dimensional tight binding model and a doping dependent hopping parameter renormalization scheme based on the values determined in ref . [ ] from the fermi surface in apres data on lsco .
the use of phenomenological input for the quasiparticle dispersion is important because the renormalization of the hopping integrals is of two - loop order in the frg analysis that we use to determine the scattering vertex .
the calculation of two - loop diagrams for the rg flow is , however , beyond the scope of this study . for details about the dispersion , cf .
[ sec : model for the quasiparticle dispersion of lsco ] .
the frg technique that we use to compute the renormalized quasiparticle scattering rates is outlined in app .
[ sec : strongly renormalized quasiparticle interactions ] together with brief overview of the doping and temperature dependence of the scattering vertex and a discussion of the parameters that enter the rg calculation .
taking the quasiparticle dispersion and the scattering vertex as input , the conductivity is computed by solving the linearized boltzmann equation with its full angular and energy dependence .
due to the strong anisotropy in the scattering rates and the quasiparticle velocities , a full solution of the boltzmann equation rather than a single relaxation time approximation is required .
the collision integral is computed by introducing an efficient discretization of the brillouin zone and the boltzmann equation is then solved numerically .
this method of solving the boltzmann equation and discretizing the brillouin zone is explained in app .
[ sec : effective model for the normal state charge transport of the high - temperature superconductor lsco ] .
the in - plane resistivity in our model of the normal state of overdoped lsco  is shown in fig .
[ fig : rhovstemperaturelsco ] ( data for @xmath17 ) .
we also plot the resistivity for different initial values for the unrenormalized onsite interaction parameter @xmath12 .
the data corresponding to the different values of @xmath18 are derived from renormalized scattering vertices while the data for @xmath19 are computed for particles interacting via the bare onsite repulsion . for a comparison of these calculations , cf .
[ subsec : discussion of the parameters ] .
the values for the hole - doping @xmath20 are chosen to match the experimental measurements of fig .
[ fig : experimental resistivity ] .
our results compare well to many qualitative aspects of the experimental data of ref .
first , there is the strong growth of the resistivity when reducing the hole doping from strongly overdoped values of @xmath21 to optimal doping around @xmath22 .
second , there is clear evidence for a growing linear term in the temperature dependence of @xmath23 around optimal doping . for large values of the hole doping
, @xmath16 displays a rather quadratic temperature dependence . on the other hand , for doping concentrations below @xmath24
a linear @xmath25-dependence becomes obvious .
we analyze our data in the same scheme as in the experimental paper [ ] using fits with the form @xmath26 ( sec .
[ subsubsec : linear plus quadratic fit ] ) or also by introducing a `` local scaling exponent '' @xmath27 , e.g. @xmath28 ( sec .
[ subsubsec : local scaling exponent ] ) .    for a quantitative comparison to the experiment
we convert our two - dimensional data to the experimental units of a three - dimensional sample ( cf . app .
[ sec : conversion of computation to experimental units ] ) . from fig .
[ fig : rhovstemperaturelsco ] we find that our calculated values for @xmath16 at @xmath29 are roughly by a factor four or five smaller than in experiment ( cf .
[ fig : experimental resistivity ] ) .
an underestimation of the resistivity is consistent with the fact that our unrenormalized interaction strength @xmath30 is considerably lower than expected for lsco .    in spite of the weak coupling approximation in our calculation of the scattering rates ,
the resistivity is in good qualitative agreement with the observed temperature dependence .
next we analyze the numerical resistivity data using the polynomial fit up to second order @xmath26 , in analogy to cooper et al , cf .
they fit their resistivity data over the entire temperature range and obtain the fit coefficients @xmath32 and @xmath33 as a function of the hole doping , cf . fig .
[ fig : coefficientslsco ] ( left panel ) .
they found the linear coefficient @xmath32 grows strongly with a decrease in doping while there is a constant or decreasing value of @xmath33 .
the linear @xmath32 and quadratic @xmath33 coefficients from fitting our numerical results are shown in fig .
[ fig : coefficientslsco ] ( right panel ) for the entire range of hole doping from @xmath34 to @xmath21 .
the fit was performed for temperatures ranging from @xmath35 up to @xmath36 . for the comparison of the two coefficients we introduce an average temperature scale @xmath37 in the conversion to resistivity units .
the coefficient representing the residual resistivity , @xmath38 , is discussed later .
the direct comparison of our fit coefficients to the experiment yields very good qualitative agreement .
the linear term shows a strong increase with reduced hole concentration toward optimal doping while above @xmath39 it becomes indistinguishable from zero . on the other hand ,
the quadratic contribution does not grow and even decreases almost down to zero toward optimal doping for larger @xmath12 .
a vanishing @xmath33 indicates a within our model a purely linear temperature dependence .
whether or not the quadratic coefficient of the fit to the experimental data vanishes is not really clear from ref .
[ ] , although it appears that the @xmath33 determined with the conventional @xmath31 fit indeed decreases and finally vanishes around @xmath40 , while the quadratic coefficient determined from a parallel resistor fit remains constant .
the coefficients in the overdoped regime disagree with the reported @xmath11 dependence by barii et al .
as already mentioned in the introduction , in the overdoped regime , their data suggest a deviation from the remarkable @xmath11 dependence of the linear term in the underdoped regime .
the experimental data of both ref .
[ ] and [ ] indicate that the linear term eventually goes to zero for sufficiently large hole density , @xmath41 .
an alternative analysis of the experiments was recently put forward by hussey et al .
taking the derivative @xmath42 they distinguished two regimes which are not immediately obvious in the resistivity data , but in which @xmath16 scales differently with temperature .
hussey and coworkers identified a high - temperature regime where @xmath43 is constant , i.e. the resistivity depends linearly on @xmath25 . toward lower temperatures @xmath44 deviates downward from this constant value . with this rather abrupt change
one can associate a temperature @xmath45 , called  coherence temperature  in ref .
[ ] which decreases with lowering the hole concentration and disappears probably around the same value of @xmath20 where the characteristic temperature for the pseudogap vanishes extrapolated from the underdoped side . in fig .
[ fig : phase diagram and coherence temperature ] the doping dependence of @xmath45 measured in ref .
[ ] is shown .
it was experimentally shown that @xmath45 coincides with the loss of the quasiparticle coherence peak [ ] .
plotting @xmath46 for different doping levels in our calculation , we can analogously identify distinct regimes in the temperature scaling of the resistivity , cf .
[ fig : coherence_temperature ] ( left panel ) .
there is a low - temperature regime with finite slope of @xmath46 , an intermediate temperature regime with finite , but smaller slope , and a high - temperature regime where @xmath46 is constant .
the temperatures that separate these three regimes will be referred to as @xmath47 and @xmath48 .
the temperature scales @xmath47 and @xmath48 are displayed in fig .
[ fig : exponentphasediagram ] .
the doping dependence of @xmath48 shows a very similar trend as @xmath45 .
this is interesting because there is no real loss of quasiparticle coherence in our semiclassical model based on the boltzmann transport equation .
note that we have not considered the values for @xmath47 and @xmath48 for @xmath49 as it becomes difficult to distinguish the regimes for @xmath40    it is interesting to relate the onset of the unconventional scaling , @xmath47 , to the proximity of the fermi energy to the van hove level , @xmath50 .
we introduce the temperature scale @xmath51 , which is shown as the red squares in fig .
[ fig : exponentphasediagram ] .
the factor of 4 is introduced for convenience , as we consider only those low - energy states relevant for the transport , that live within an energy range of @xmath52 from the fermi level .
the doping dependence of @xmath53 and @xmath47 is very similar and we conclude that proximity to a van hove singularity plays an important role in the development of non - fermi liquid behavior .    to compare to the analysis in ref .
[ ] we fit the high - temperature linear part of @xmath43 with a constant @xmath54 and the low - temperature regime by a linear function @xmath55 . in the high - temperature regime
, @xmath54 represents the coefficient of the dominating linear term in the resistivity , while @xmath56 corresponds to the coefficient of the dominating linear term in the limit @xmath57 .
the results for these coefficients are shown in fig .
[ fig : coherence_temperature ] ( right panel ) .
( @xmath58 ) as a function of the hole doping of lsco .
the coherence temperature is where the temperature derivative of the conductivity deviates from its constant value at high temperature ( data adopted from ref .
[ ] ) . , scaledwidth=33.0% ]          a different way to illustrate the regime change on the resistivity data is to define an exponent of @xmath16 for each temperature formulated as @xmath59 the exponent @xmath60 is obtained by taking the logarithmic derivative of the @xmath16 , @xmath61 where @xmath38 denotes the residual resistivity .
the results for @xmath60 over a wide range in doping are shown in fig .
[ fig : exponentphasediagram ] .     defined in eqs . and for our model of lsco  in color - code .
the temperature scales @xmath47 and @xmath48 , cf .
[ subsubsec : scaling regimes ] , are shown as black squares . the temperature scale @xmath53 corresponds to the energy difference between the chemical potential and the energy level of the van hove singularity , @xmath62 ( we consider low - energy states up to an energy of @xmath52 around the fermi energy ) .
note that the unconventional scaling @xmath63 goes down to the lowest temperatures at optimal doping .
the lifshitz transition is indicated as the horizontal dashed line . ]    even though such a local scaling law has large uncertainties , it is intriguing to find a very good qualitative match between the function of @xmath27 and the generic scaling behavior of the cuprates , cf . and ref .
in the overdoped low - temperature region of the phase diagram , the local exponent rises to the value of @xmath64 corresponding to a standard fermi liquid .
around optimal doping , @xmath27 is always well below @xmath65 down to the lowest temperature of our study . in the high - temperature regime ,
the onset of unconventional scaling ( @xmath66 ) is shifted to higher temperature with increasing doping .      using the non - equilibrium distribution function
, we can analyze the angular dependence of this unconventional scaling and compare to the difference in the temperature dependence between isotropic and anisotropic scattering rates in the frg study [ ] .
the quasiparticle scattering rates are strongly anisotropic and especially the quasiparticles close to the saddle points ( anti - nodal - direction ) in the brillouin zone have very short transport life times .
the scattering rates that cause these short life times are strongly enhanced as the doping level decreases toward optimal doping . to examine the role of strong scattering of the anti - nodal quasiparticles for the unconventional temperature dependence of the overall resistivity
, we interpret the angular patches of our brillouin zone discretization ( cf .
[ fig : discretization ] in the appendix ) as individual resistors and parametrize the low - energy quasiparticle states as a parallel resistor network .
each resistor is characterized by the conductivity @xmath67 and resistivity @xmath68 where the index @xmath69 labels the the angular patch . the total conductivity and resistivity
are determined by @xmath70 in fig .
[ fig : angularresolutioncoefficients ] we depict the the angle - resolved scaling by plotting the linear and quadratic contributions to the patch resistivity as functions of the angular patch index @xmath69 and the hole concentration .
the temperature dependence of each resistor was analyzed by a linear plus quadratic fit , @xmath71 . in order to compare linear and quadratic term ,
we give them a weight factor that represents their contribution to the total resistivity and normalize their sum , the total resistivity @xmath23 , to unity , @xmath72 thus @xmath73 represents the normalized contribution of the angular patch @xmath69 to the term in the total resistivity that scales with temperature exponent @xmath74 .    ) . on top
, we have plotted the weighted linear term @xmath75 , and the lower plot shows the weighted quadratic term @xmath76 , defined in eq . .
these @xmath73 s are weighted according to the contribution of the angular patch @xmath69 to the total charge transport .
the patches range in angular direction from @xmath77 ( patch 1 ) to @xmath78 ( patch 10 ) .
the @xmath73 s are evaluated at the temperature @xmath37 , the average temperature of the fit region . ]    from fig .
[ fig : angularresolutioncoefficients ] one can clearly see that in the strongly overdoped regime only the quadratic term is present and is essentially independent of the angle . on the other hand , below @xmath79 ,
the linear term grows showing strong anisotropy .
the linear term appears first close to the nodal direction ( patches 4 to 7 ) , and constantly spreads out .
below @xmath80 the linear term starts to dominate the transport .
the quadratic term , however , remains finite in the nodal direction down to the lowest doping levels .
it is interesting to see that around optimal doping , where the fermi surface is very close to the van hove singularity , the charge transport from the angular patches 1 and 10 , that contain the saddle points , essentially vanishes .
this is consistent with the strong growth of the scattering rates from the frg analysis .
these angular patches are short - circuited as a consequence of strong scattering , cf .
our discussion here is fully compatible with the behavior of the angle - dependent scattering rates found by the frg study by ossadnik and collaborators [ ] .
we can also analyze the current distribution as a function of the temperature . in fig .
[ fig : angular distribution of the current ] we display the angular distribution of the charge current for two different doping levels . at the lowest temperatures
the distribution function changes its angular dependence strongly and the weight of the current is shifted from small angles to larger angles , while above a certain temperature the angular distribution remains invariant .    these regimes can be associated with two different kinds of scattering mechanisms . while at low temperature impurity scattering dominates , for sufficiently high temperature the transport is mainly determined by particle - particle scattering . since only umklapp
scattering contributes noticeably to the resistivity , the two - particle scattering contribution is anisotropic , while the @xmath81-potential impurities give isotropic scattering .
this difference accounts for a qualitative change in the angular dependence of the distribution function .
we will return to the interplay of the two scattering mechanisms in the context of a discussion of matthiessen s rule , sec .
[ sec : momentum relaxation mechanism in the cuprates and matthiessens rule ] .
we turn now to the role of the renormalization of the scattering vertex in the unconventional nature of the transport properties .    as explained in app .
[ sec : choice of parameters ] we used a value of @xmath82 for lsco .
when we extend our calculation to a series of lower initial values of @xmath83 the temperature dependence of the resistivity and the coefficients of the corresponding linear plus quadratic fit are shown fig .
[ fig : rhovstemperaturelsco ] and in fig .
[ fig : coefficientslsco ] .
in addition we show the resistivity for the unrenormalized interaction @xmath84 which corresponds to an frg treatment with infinite cutoff . in this case
the scattering rates are determined only by the scattering phase space which is specified by the fermi function and the geometry of the fermi surface .    comparing the resistivity curves for @xmath17 and @xmath85
, we observe very similar overall trends as function of doping .
the high - temperature resistivity grows upon reducing the hole - concentration and the @xmath4-linear part increases at the same time .
it is obvious that the scattering phase space plays an important role in the change of the behavior of resistivity for different hole - concentrations .
overall the behavior is , however , less pronounced for the unrenormalized than renormalized vertex as is clear in fig .
[ fig : coefficientslsco ] .
the renormalization of the vertex causes a sizable increase of the resistivity as one can see in the comparison of the curves for @xmath86 and @xmath85 , and it also strongly enhances the trend toward linear resistivity .
this feature is also easily identified from the @xmath12-dependence of the quadratic term in fig .
[ fig : coefficientslsco ]
. for stronger initial interaction , and therefore stronger renormalization effects , the quadratic term is substantially reduced , even as the total resistivity , the sum of linear and quadratic term is increased .
thus , it is evident that correlation effects included in the vertex renormalization by the frg scheme , make a substantial contribution of the unconventional transport properties , in particular , in the regime close to optimal doping .
the vertex renormalization leads to a non - trivial momentum space structure of the effective scattering vertex which can be translated to a longer ranged real space interaction compared to the local hubbard-@xmath12 . as a second effect of the frg analysis , the momentum space structure but also the absolute values of the scattering vertex become temperature dependent . in order to discriminate between these two effects we modified our calculations in a way as to include only the momentum space renormalization of the vertex , but neglect the temperature dependence
this was implemented by stopping the rg flow with a large cutoff @xmath87 . in this case
_ large _ refers to the cutoff exceeding the largest temperature of our study .
while this completely suppresses the temperature dependence of the scattering vertex , it still introduces momentum space anisotropy .
the result resembles the behavior of the simple unrenormalized onsite repulsion @xmath88 of fig .
[ fig : rhovstemperaturelsco ] .
the introduction of momentum anisotropy through the rg flow has almost no qualitative effect on the temperature dependence of the resistivity . from this comparison
we can conclude that the pronounced linear resistivity in the limit of @xmath57 is a consequence of the strongly temperature dependent quasiparticle scattering rates .
this supports the idea of ref .
[ ] where the unconventional temperature dependence of the scattering rates through the rg flow were described .    at high temperature , however , the growth of the linear - temperature term when the doping is decreased toward optimal doping , is mainly due to the fermi surface geometry and the temperature dependence of the scattering volume , while the momentum anisotropy of the scattering interaction is responsible for the enhanced value of the resistivity .
in the cuprates the resistivity is caused by electron - electron interactions which raises the question , how the total momentum of the entire distribution function is relaxed . in the following
we will discuss which contributions are combined in the resistivity calculations of the previous section .      in impurity scattering processes ,
the momentum difference between initial and final state is transferred to the lattice and determines the residual resistivity as the zero - temperature limit .
a deeper analysis shows that the resistivity from impurity scattering can have a weak temperature dependence due to thermal broadening of the distribution function and a possible bias toward higher or lower quasiparticle velocity . in our calculation the impurities are modeled as @xmath81-potentials .    a second and more complex contribution to the electrical resistivity originates from electron - electron scattering with finite momentum transfer to the lattice , so - called umklapp scattering , where a momentum corresponding to a reciprocal lattice vector is transferred to the lattice .
note that the fermi surface has to be sufficiently large in order to allow for umklapp assisted scattering between different momenta on the fermi surface .    for a square lattice we can distinguish between two kinds of umklapp scatterings : with momentum transfers of the type @xmath89 and @xmath90 to the lattice .
the first kind of umklapp scattering has a large phase space for a sufficiently large fermi surface with regions where the @xmath91-component of the fermi vector satisfies the condition @xmath92 and other regions with @xmath93 .
the second type of umklapp process is only possible , if the fermi surface intersects with the so - called umklapp surface , which is the rotated square ( diamond ) that connects the four saddle points at @xmath94 and @xmath95 .
this umklapp scattering has relatively small phase space due to the stringent momentum constraints . within our model for lsco
the @xmath89 umklapp process largely dominates over the @xmath90 umklapp process for all doping levels as far as the charge transport resistance is concerned . on the strong overdoped side ,
@xmath90 umklapp scattering is even completely absent in the low - temperature limit .    due to the geometrical constraints for umklapp scattering
the momentum space structure of the scattering vertex can play an important role .
the temperature dependence of the resistivity is controlled by umklapp scattering .
if the scattering rates become especially high on certain parts of the fermi surface , their contributions to transport will be suppressed yielding the possibility for unconventional temperature scaling .
often different contributions to the resistivity are assumed to be decoupled and one adds them as serial resistors .
this decoupling leads to matthiessen s rule , which is only justified , if each individual contribution can be well described by a single relaxation time . in this case
the collision integral in the boltzmann equation directly splits into the scattering channels . in the general case ,
when such a single relaxation time approximation is insufficient , the collision integral does not decouple .    in fig .
[ fig : angular distribution of the current ] we have seen that in the low - temperature regime the angular distribution of the current changes its form .
this was attributed to a crossover from an impurity scattering dominated regime to a two - particle scattering regime .
these different scattering mechanisms lead to different current distributions due to their difference in the scattering geometry . while impurity scattering is isotropic , the umklapp processes responsible for the momentum relaxation in the two - particle scattering channel are highly anisotropic leading to a non - trivial interplay of impurity scattering with two - particle interactions and a potential breakdown of matthiessen s rule .    to examine the deviations from matthiessen s rule we compare our calculation with the sum of the  individual  contributions to the resistivity .
first we remove the impurities to obtain a  bare
electron - electron resistivity @xmath96 , second we estimate @xmath97 by considering free ( non - interacting ) electrons scattered by impurities .    in a first step
we compare the linear and the quadratic coefficients of a @xmath98-fit for the clean limit , @xmath99 , to the coefficients in the original calculation including impurities of fig .
[ fig : coefficientslsco ] .
the quadratic coefficient grows weakly down to @xmath24 and remains constant at @xmath100 , while for the full calculation the quadratic coefficient decreased below @xmath79 and even went to zero for sufficiently strong coupling .
furthermore , the linear term is somewhat smaller without impurities .
this difference in the behavior of the coefficients is a violation of matthiessen s rule .    in order to make this more evident ,
we plot the deviation from matthiessen s rule , @xmath101 , quantitatively defined as @xmath102/\rho_{\text{full } } , \label{eq : deviation from matthiessens rule}\end{aligned}\ ] ] as a function of the temperature for different doping levels .
the result is shown in a color plot in fig .
[ fig : matthiessenruledeviation ] .
we immediately observe that the deviation from matthiessen s rule is strongest at the optimally doped regime at low temperature .
the difference in the scattering geometry for impurity scattering and two - particle scattering is most pronounced in this part of the diagram . here
the deviation @xmath101 reaches values as large as 16% .    it should be noted that the one - loop frg approximation is less reliable at low temperature and low doping , as pointed out earlier .
hence the magnitude of @xmath101 in our calculation carries considerable uncertainties in this part of the diagram .
nevertheless , the violation of matthiessen s rule is an generic property for systems featuring a crossover from isotropic impurity scattering to strongly anisotropic inter - particle scattering dominated transport .
[ t ]     and normalize it to @xmath103 .
it is clearly visible that the deviation @xmath101 ( cf .
is most pronounced in the critical region of the phase diagram . ]
the seebeck effect is a further example which allows us to test fermi liquid properties in transport .
this effect as the electric response to a temperature gradient is based on thermo - diffusion of charge carriers and can be analyzed with boltzmann transport theory too .
lalibert et al .
measured the seebeck coefficient @xmath104 of @xmath105 [ eu - lsco ] , cf . ref .
while their study aims at probing the fermi surface reconstruction scenarios in the underdoped regime of the phase diagram , they also present data up to a hole - doping level of @xmath106 .
data of the seebeck coefficient are displayed in fig .
[ fig : experimental seebeck eu - lsco and numerical seebeck ] ( left panel ) for the doping regime @xmath107 .
europium doping of lsco  is believed to have a stronger influence on the properties of underdoped samples , but little effect is expected in the overdoped regime .        for underdoped samples @xmath108 shows a sign change , starting negative at low temperature it turns positive at higher temperature . on the other hand , for overdoped samples
@xmath108 remains positive for all temperatures measured and shows a rather strong increase toward low temperature near optimal doping .
the inset of fig .
[ fig : experimental seebeck eu - lsco and numerical seebeck ] ( left panel ) demonstrates the @xmath109 dependence of @xmath110 at low temperature , which has been interpreted as a signature of quantum criticality .
we have numerically calculated the seebeck coefficient for lsco  within our semiclassical approach , as explained in app .
[ subsec : calculation of the seebeck coefficient ] .
the results are shown in figs .
[ fig : experimental seebeck eu - lsco and numerical seebeck ] ( right panel ) and [ fig : seebeck phase diagram ] . a good overview on the behavior of the seebeck coefficient as a function of doping and temperature is obtained from fig . [ fig : seebeck phase diagram ] ( left panel ) . at the doping
@xmath111 we observe a sign change as we cross the lifshitz transition at low temperature .
interestingly , for @xmath112 increasing the temperature leads to a sign change from negative at low to positive at higher temperature . turning to the other limit of strong overdoping
the seebeck coefficient turns negative again as expected for electrons is a band with small filling .
we also observe a doping dependent temperature scale where @xmath104 displays a rapid increase .
interestingly , this temperature compares well with @xmath47 .    in order to compare our numerical with the experimental results ,
we have plotted @xmath110 in fig .
[ fig : experimental seebeck eu - lsco and numerical seebeck ] for some selected values of hole doping .
the inset shows the low - temperature regime for the doping level @xmath113 on a logarithmic temperature scale .
the general structure of the calculated seebeck coefficient is in good qualitative agreement with the experimental data measured for eu - lsco  fig .
[ fig : experimental seebeck eu - lsco and numerical seebeck ] . in both cases
, there is a critical doping level below which the seebeck coefficient changes sign with increasing temperature , starting from negative in the @xmath57 limit . above this critical doping level , @xmath108 becomes very large for low temperature and strongly decreases with temperature .
the low - temperature rise of @xmath108 becomes less pronounced with increasing hole doping . in the low - temperature regime
, @xmath108 acquires a @xmath109 dependence which is clearly visible in both the experimental data and the numerical data below @xmath114 .
a small discrepancy between our model and the experiment is the value of the critical doping level .
we have designed the band structure that we use to describe lsco  such that the lifschitz transition is located at @xmath115 , cf .
[ sec : model for the quasiparticle dispersion of lsco ] .
the van hove singularity at the lifschitz transition is the origin for the critical behavior in our model and thus constrains the critical doping level to @xmath116 . in the experimental study , however , a critical doping level of @xmath117 is reported .      in the last section
, we have discussed the low - temperature @xmath109 dependence of @xmath108 and the sign changes with doping or temperature .
these properties of the seebeck coefficient are often taken as indication for the presence of a quantum critical point [ ] .
these features are found in both the experimental data and the numerical simulation .
we now want to locate the areas of critical behavior in the phase diagram .    for conventional metals ,
the seebeck coefficient is a linear function of the temperature manifesting fermi liquid behavior .
deviations from this linear dependence may be used as a signature for non - fermi liquid physics .
thus , we analyze the the temperature derivative of @xmath110 which should vanish for fermi liquids . fig .
[ fig : seebeck phase diagram ] ( right panel ) displays this function .
this value is largest at low temperature around optimal doping and a sign change appears at the lifshitz transition .
this is the region where the strongest differences with the fermi liquid picture emerge . from this finding
we may identify the vicinity of the van hove singularity and the presence of a lifshitz transition as the origin of non - fermi liquid behavior and the signatures of quantum criticality .
interestingly , fig .
[ fig : seebeck phase diagram ] ( right panel ) allows us also to see the temperature of the sharp increase of @xmath108 , cf .
[ fig : seebeck phase diagram ] ( left panel ) , appearing as a maximum of @xmath118 .
this characteristic temperature @xmath119 separates the high - doping low - temperature regime coming closest to fermi liquid behavior ( @xmath120-dependence of resistivity ) from the lower - doping high - temperature regime which is characterized by the non - fermi liquid @xmath4-linear resistivity , as discussed in sec .
[ subsubsec : scaling regimes ] .
we investigate the temperature and doping dependence of charge transport coefficients using a combination of an experimentally determined band structure and a renormalized scattering matrix calculated in a one - loop renormalization group approximation .
the band dispersion matches the quasiparticle spectrum in arpes for lsco [ ] , and as such includes the renormalization of the dispersion with hole density due to interactions .
the electron - electron scattering vertex develops momentum - space and temperature dependence as the hole density is reduced .
we have chosen a moderate to weak value for the bare onsite coulomb repulsion consistent with the one - loop approximation .
we also refrain from including the pseudogap region which would require a solution of the strong coupling regime that results from divergences in the frg flow .
the boltzmann transport equation is solved numerically with collision integrals for electron - electron and impurity scattering determined by fermi s golden rule .
this semiclassical approach determines the quasiparticle distribution functions taking into account the full angular and energy dependence of the collision integral .
we find good qualitative agreement with the main trends in the experiments on the in - plane resistivity of lsco   [ ] .
this includes the increase of the temperature dependent term in the resistivity upon decreasing the hole - concentration and the trend toward a linear temperature dependence of the resistivity @xmath121 around optimal doping .
a parameterization of our results using a polynomial fit in @xmath4 up to second order , inspired by the experimental analyses [ ] , gives consistency between the trends in experiment and our model .
while the @xmath4-linear behavior of @xmath121 is often attributed to the existence of a quantum critical point around optimal doping [ ] , our discussion shows that this feature may also be understood as a consequence of the unusual temperature dependence of the renormalized scattering vertices and the increased scattering phase space around this doping regime , where a lifshitz transition of the band structure leads the chemical potential very close to a van hove singularity .
the magnitude of the calculated resistivity is substantially lower than in the experiment , consistent with the moderate to weak interaction values and one - loop rg approximation    motivated by the experimental data analysis [ ] we also identified crossover regimes in @xmath122 and find rather sharp changes in the temperature dependence .
there are two characteristic temperatures @xmath123 and @xmath48 where the latter denotes the lower bound of the @xmath124-linear regime of @xmath121 and @xmath125 the upper bound for a purely @xmath126-dependence ( fermi liquid regime ) .
the intermediate temperature range may be considered as combination of both .
the resulting phase diagram gives a very good account of the trends of @xmath121 observed in the experiment .
in particular , we can show that there is a clear correlation of the unconventional temperature scaling with the van hove singularity , quantified through the effective distance of the chemical potential from the van hove level , @xmath127 , cf .
[ fig : exponentphasediagram ] .    by comparing the temperature dependence of @xmath16 for renormalized and bare scattering rates
, we could demonstrate the importance of the scattering vertex renormalization .
while in the strongly overdoped regime , the difference in the numerical data determined from renormalized and unrenormalized scattering rates is not very pronounced , around optimal doping the phase space driven growth of a linear term is strongly enhanced due to renormalization of the quasiparticle interactions .
it also turns out that the impurity scattering is important in this respect as well , contributing to the unconventional temperature dependence .
this is surprising as impurities contribute usually only a weakly temperature dependent offset to resistivity .
we found deviations from matthiessen s rule in large parts of the considered phase diagram .
this deviation is especially pronounced as we approach optimal doping at intermediate to low temperature due to the different momentum dependence of impurity and electron - electron scattering .
matthiessen s rule , however , is recovered in the more conventional , strongly overdoped regime at low temperature .
we have also studied the seebeck coefficient within our model .
we find qualitative agreement between our numerical data and the main trends in the measurements performed on a slightly different compound , eu - lsco .
we believe that the eu - doping does not strongly influence the seebeck coefficient in the optimal to overdoped doping regime which we are interested in .
the calculated seebeck coefficient , also displays a critical region in the phase diagram with pronounced deviations from fermi liquid behavior .
we have related the critical behavior to proximity of the van hove singularity to the fermi energy which is linked to the lifschitz transition .
we are grateful to carsten honerkamp , nigel hussey , lszl forr , and stevan arsenijevi for helpful discussions .
this study was supported by the swiss nationalfonds , the nccr manep , the hittec project of the competence center energy & mobility and the synergia teo .
we use a band structure parametrized within a tight - binding model on a square lattice , including nearest , next - nearest , and next - to - next - nearest neighbor hopping @xmath128 the hopping parameters @xmath129 , @xmath130 , @xmath131 , and the chemical potential @xmath132 define the geometry of the fermi surface and the filling of the band .    based on their arpes study , yoshida et al .
[ ] have fitted the low - energy quasiparticle spectrum of lsco  and the fermi surface in the doping range from @xmath133 to @xmath134 using the model . fixing he magnitude of @xmath129 to @xmath135 and the ratio @xmath136 for all dopings they used @xmath132 and @xmath130 as fitting parameters .
their results are shown in fig .
[ fig : discretization ] .
the relation between @xmath132 and @xmath130 as a function of the doping is linear to a very good approximation ,    @xmath137    yoshida et al .
identified a lifshitz transition at @xmath138 .
our hopping parameter renormalization scheme , eq .
, locates the lifshitz transition at @xmath139 .
changing @xmath140 to @xmath141 shifts the lifshitz transition within our model to the experimental doping level .
note , that a second arpes study by ino et al .
[ ] also confirms the presence of a lifschitz transition at a doping level of about @xmath142 .
the calculation of the renormalized vertex in ref .
[ ] is based on the frg equation for the one - particle irreducible generating functional derived in refs . [ ] which gives hierarchically coupled flow equations for the one - particle irreducible vertices .
the rg - flow follows from a wilsonian flow scheme with a sharp momentum cutoff .
the rg equations are solved using the standard truncation of all vertices with more than four legs , so that the self energy and the four - point vertex are the only remaining quantities in the calculation . for the vertex ,
the momentum dependence is discretized and constrained to momenta lying on the fermi surface .
the frequency dependence of the vertex is suppressed . for the rg - flow ,
the flow parameter @xmath87 is taken from @xmath143 to zero .
in general , the truncated flow tends to diverge for low temperature at finite energy scales @xmath87 , indicating the appearance of an instability .
the energy scale at which the divergence occurs is related to the critical temperature @xmath3 of the corresponding instability , in our case superconductivity .
note , however , that due to the truncation of the rg equations the @xmath3 is largely overestimated .    for a study of normal state transport , the divergence of the vertex should be avoided .
while in the experiment superconductivity is suppressed by a magnetic field which is difficult to incorporate in the rg analysis , we keep thge cutoff energy scale @xmath144 finite and so suppress the leading @xmath145-wave pairing instability in the rg flow .
the choice for @xmath146 is discussed in app .
[ sec : choice of parameters ] .
the quasiparticle scattering vertex @xmath147 is peaked for momenta @xmath148 lying close to the saddle - points @xmath149 and @xmath150 and has especially strong contributions for momentum transfer of @xmath151 , cf .
[ fig : vertex ] .
further discussions of the rg flow and interpretations of vertex diagrams as in fig .
[ fig : vertex ] are given in refs .
the frg method is a weak coupling analysis . on the underdoped side , below a temperature scale @xmath152
a pseudogap opens as an indication that the system is driven into a strong coupling phase .
the strong - correlation physics are not included in our analysis .
hence , we restrict our study to dopings from optimal to strongly overdoped .      the onsite repulsion @xmath12 for the cuprates is a large energy scale of the order of the band width as at half - filling a mott insulating phase is realized . using such large values of @xmath12 , however , leads to a divergence of the rg flow within the normal state temperature range , since the instability temperatures are overestimated .
this urges us to start with moderate value of the onsite repulsion of @xmath153 which then is renormalized to much larger values for the low - energy quasiparticles after integrating out the high - energy states .
moreover , we suppress the @xmath145-wave pairing instability by keeping the cutoff @xmath144 finite in the rg flow , as mentioned above . for this purpose
we choose @xmath154 corresponding to the lowest temperature in our study .
the boltzmann equation relates the substantial derivative of the quasiparticle distribution function to the collision integral , taking all scattering events into account which yield a change of the distribution function , @xmath155\cdot{\boldsymbol{\nabla}}_{p}f({\boldsymbol{r}},{\boldsymbol{p}},t)\phantom{\int}\nonumber\\      & = \left.\frac{\partial f({\boldsymbol{r}},{\boldsymbol{p}},t)}{\partial t}\right|_{\text{coll}}. \label{eq : the boltzmann equation}\end{aligned}\ ] ] we consider two contributions in the collision integral , two - particle collisions labeled by `` ee '' and impurity scattering labeled by `` imp '' : @xmath156    the first contribution is expressed as , @xmath157with@xmath158 where the transition rates @xmath159 are generated through fermi s golden rule from the matrix elements of the quasiparticle interaction vertex @xmath160 , cf .
[ sec : strongly renormalized quasiparticle interactions ] .
the second contribution to the collision integral is simpler due to its single particle form , @xmath161with@xmath162    we ignore electron - phonon contributions here because , in the cuprates , the momentum relaxation of the quasiparticles is dominated by electron - electron interactions as was concluded e.g. from the temperature dependence of the hall coefficient [ ] and from the observation of a quadratic resistivity for overdoped lsco  [ ] .
it is important to notice that two - particle scattering can only yield a finite resistivity , if momentum is transferred to the lattice , which only can be realized through umklapp scattering ( see sec .
[ sec : momentum relaxation mechanism in the cuprates and matthiessens rule ] ) .
umklapp scattering has special geometrical constraints and , therefore , introduces strong anisotropy in the brillouin zone and requires a detailed analysis of the collision integral beyond a single - relaxation - time approach .
we linearize the boltzmann equation in terms of the deviation from equilibrium @xmath163 .
the collision integral is interpreted as a linear integral operator acting on @xmath163 .
we evaluate the integral kernel on discrete patches of a finite mesh ( see fig.[fig : discretization ] ( a ) ) .
note in this context that that the energy conservation appearing as delta function in the scattering rates has to be treated carefully .
in discretized momentum space , the linearized boltzmann equation is reduced to a set of linear equations which are solved numerically .    in order to match the fermi surface ,
we use polar angle and the energy to represent each momentum vector . for our curved - patch discretization
we can adjust the grid spacing in radial ( energy ) direction to the temperature scale and always cover the relevant scattering phase space .
a sketch of this discretization is given in fig .
[ fig : discretization ] ( a ) . for the calculations that we present here we have chosen a discretization of 24 patches in radial direction and 120 patches in angular direction .
the latter are not distributed equally but are much finer in anti - nodal than in nodal direction to take care of the flat dispersion in the proximity of the saddle points .      for a study of the seebeck effect
, the boltzmann equation is extended by a new drift term that accounts for the response of the distribution function to the external thermal gradient , @xmath164 with the electro - chemical potential defined as @xmath165 .
the collision integral remains unchanged .
the distribution function that solves this equation carries a heat current and an electrical current .
these two currents are defend as the linear response to the external fields as @xmath166 considering an open circuit geometry with @xmath167 , the typical setup for thermopower measurement , the seebeck coefficient is defined as the proportionality factor between the applied thermal gradient and the induced electrochemical potential , @xmath168 . for simplicity
we assume all forces aligned along one direction and , thus , we can use scalar quantities .
it is important to note that we do not make use of analytic simplifications as mott s formula or sommerfeld expansion .
according to the arpes study [ ] used above to find the effective tight - binding parameters of lsco , a value of @xmath135 for the nearest neighbor hopping parameter @xmath129 fits the arpes data best .
this defines an energy scale which translates the energy scales of temperatures in units of @xmath129 to @xmath169 .
if we want to study the temperature range up to room temperature of @xmath36 the corresponding energy scale is @xmath170 . in order to cover at least one order of magnitude in temperatures we choose @xmath171 and @xmath172 .
our model is based on a two - dimensional interacting electron gas on a lattice .
the experimental setup is based on the measurement of the resistivity of a three - dimensional but layered sample . for a quantitative comparison of the experimental and numerical data , a conversion of the two - dimensional resistivity to a three - dimensional resistivity
is required . for a layered system
this is very easily achieved by multiplication of the two - dimensional resistivity with the interlayer distance .
the c - axis lattice constant of lsco  is given by @xmath173 [ ] and the interlayer distance by half the lattice constant , because the unit cell of the 2 - 1 - 4-compounds contains 2 layers .
f. lalibert , j. chang , n. doiron - leyraud , e. hassinger , r. daou , m. rondeau , b.j .
ramshaw , r. liang , d.a .
bonn , w.n .
hardy , s. pyon , t. takayama , h. takagi , i. sheikin , l. malone , c. proust , k. behnia , and l. taillefer : nature comm .
* 2 * ( 2011 ) 432    t. yoshida , x.j .
zhou , k. tanaka , w.l .
yang , z. hussain , z .- x .
shen , a. fujimori , s. sahrakorpi , m. lindroos , r.s .
markiewicz , a. bansil , seiki komiya , yoichi ando , h. eisaki , t. kakeshita , and s. uchida : phys .
b * 74 * ( 2006 ) 224510