Extremely Correlated Fermi Liquid theory of the $t$-$J$ model in 2 dimensions: Low Energy properties

Low energy properties of the metallic state of the 2-dimensional tJ model are presented at various densities and temperatures for second neighbor hopping t', with signs that are negative or positive corresponding to hole or electron doping. The calculation employs a closed set of equations for the Greens functions obtained from the extremely correlated Fermi liquid theory. These equations, when used in $d=\infty$ reproduce most of the known low energies features of the $U=\infty$ Hubbard model. In 2-dimensions we are able to study the variations due to the superexchange J. The resulting Dyson self energy is found to be momentum dependent as expected. The density and temperature dependent quasiparticle weight, decay rate and the peak spectral heights over the Brillouin zone are calculated. We also calculate the resistivity, Hall conductivity and cotangent of the Hall angle in experimentally relevant units. These display significant thermal sensitivity for density n>~ 0.8, signifying an effective Fermi-liquid temperature scale which is two or three orders of magnitude below the bare bandwidth. Flipping the sign of the hopping t', i.e. studying hole versus electron doping, is found to induce a change in curvature of the temperature dependent resistivity from convex to concave at low temperatures. Our results provide a natural route for understanding the observed difference in the temperature dependent resistivity of strongly correlated electron-doped and hole-doped matter.


I. INTRODUCTION
The t-J model in 2-dimensions (2-d) has been argued to be of fundamental importance for understanding strongly correlated matter, including the high Tc superconductors 1,2 . Due to the difficulties inherent in a strong coupling problem, very few techniques are available for extracting its low temperature physics. Towards this end we have recently developed the extremely correlated Fermi liquid (ECFL) theory 3,4 . It is an analytical method for treating very strong correlations of lattice Fermions, employing Schwinger's technique of functional differential equations togather with several important added ingredients. While further details can be found in Ref. (3,4), a brief summary of the main idea behind the ECFL theory seems appropriate. We consider the Hubbard model with a large interaction U → ∞, and hence the name of the theory. A well known expansion in the inverse powers of U leads to the t-J model (defined below 2 ). Taking the infinite U limit forces one to abandon the conventional Feynman diagram based perturbation theory in U , and to make a fresh start.
The ECFL theory starts with the graded Lie-algebra of the Gutzwiller projected, i.e. infinite-U limit Fermi operators Eqs. (2,3). This leads to an exact functional differential equation for the Greens functions, known as the Schwinger equation of motion Eq. (18 or 22). In this equation, a parameter λ is introduced; λ is bounded in the range ∈ [0, 1] and represents the evolution from the free Fermi limit. We then use a systematic expansion in the parameter λ, for solving the Schwinger equations perturbatively in λ. In this scheme we start with the uncorrelated Fermi gas at λ = 0 and end up at the fully correlated projected Fermion problem at λ = 1. The scheme thus represents a generalization of the usual perturbation theory for canonical Fermionic models, in order to handle a non-canonical Fermionic problem such as the t-J model. The context of interacting Bosons provides a useful parallel. In the well known problem of representing spin S variables in terms of canonical Bosons, one uses the expansion parameter 1 2S with a similar range 1 2S ∈ [0, 1]. We may think of λ as being analogous to the parameter 1 2S as shown in Ref. (4). The introduction of the parameter λ and the λ-expansion scheme thus enabled are among the main technical advances introduced in the ECFL theory.
This approach has been recently benchmarked [5][6][7] against the numerically exact results from the single impurity Anderson model, and the d = ∞ Hubbard model from Dynamical Mean Field Theory (DMFT) 9,10 . These tests provide quantitative support to our general scheme described below, especially for low energy response. Our scheme has no specific limitation to d = ∞, and is expected to be reasonably accurate in any dimension d > 1, including 2 dimensions, a case of great experimental importance due to the High Tc cuprate materials. It is applied here to probe the metallic state of the t-J model in 2-d. We present results for the electron self energy, the spectral functions, the resistivity the Hall constant and the Hall angle at various temperatures and electron density n = N/N s (number of electrons per site). We also use the notation of hole density δ = (1 − n) instead of n, following the convention used in several experimental studies of doped Mott systems..
We explore various values of the parameters of the t-J model, including the second neighbor hopping, which turns out to play a very important role in determining the effective Fermi Liquid (FL) temperature scale. We investigate the resistivity due to mutual collisions of electrons at low temperatures, and its dependence on the parameters of the model. We pay special attention to the resistivity since this easily measured -but notoriously hard to calculate object, reveals the lowest energy scale physics of charge excitations in metallic systems, and therefore is of central importance.

II. METHODS
In this section we summarize the equations used in the present calculation, together with the arguments leading to them-further details may be found in earlier papers on this The t-J model is a two component Fermi system on a lattice, defined on the restricted subspace of three local states, obtained by excluding all doubly occupied configurations. The allowed states at a single site are |a with a = 0, ↑, ↓, and the double occupancy state | ↑↓ is removed by the (Gutzwiller) projection operator P G = Π i (1 − n i↑ n i↓ ). We use the Hubbard The λ expansion and the auxiliary Greens function The main task is to compute solutions of the Schwinger-Dyson equation, i.e. the functional differential equation Eq. (18) or Eq. (22). If symmetry-breaking, such as magnetism or superconductivity is ignored, then a liquid state ensues, where we would like the solution to connect continuously with the Fermi gas. For this purpose we seek guidance from standard Feynman-Dyson perturbation theory for canonical models. The repulsive Hubbard model is an ideal example, where the corresponding Schwinger-Dyson equation can be schematically written as: Comparing with Eq. (22), we see that the left-hand sides are of the same form, but the righthand sides differ, in Eq. (22) the local Greens function γ multiplies the delta function. In turn this extra term originates from the second (non canonical) term in the anti-commutator in Eq. (2), and is therefore the signature term of extremely strong correlations.
Within the Schwinger viewpoint of Eq. (23), we can view the skeleton graph perturbation theory (Feynman-Dyson) as an iterative scheme in U , i.e. using the n th order results to generate the (n + 1) th order terms by functional differentiation. In the ECFL theory the iterative scheme used is defined by generalizing Eq. (22) to The explicit solutions in the ECFL theory start from this basic equation. More explicitly, in Eq. (24) the exact Eq. (18) is generalized to include the λ parameter 38 by scalinĝ We found in Ref.
(3) that an efficient method for proceeding with the iterative scheme is to first perform a factorization of the Greens function into two parts. The first is an auxiliary Greens function g satisfying a canonical equation, thus admitting a Dysonian expansion with its attendant advantage of summing a geometric series with every added term of the denominator. There remain some terms that cannot be pushed into the denominator, these are collected together as the caparison function µ. In the matrix notation used above we first decompose the Greens function as: this implies a product in the k, ω domain as written below in Eq. (32). The differential operator X in equation Eq. (24) is distributed over the two factors of Eq. (25) using the Leibniz product rule, as X.g. µ = X.g. µ + X.g. µ (26) where the contraction symbol Xa indicates the term being differentiated by the functional derivative terms in X, while the matrix indices follow the dots. Using g −1 .g = 1 Eq. (24) is now written as This equation factors exactly into two equations upon insisting that g has a canonical structure: and We can then use g.g −1 = 1 to simplify the term X.g.g −1 = −X.g.g −1 , giving rise to a Dyson self-energy expressed in terms of a Dyson vertex functions as usual. The idea then is to iterate the pair of Equations (28,29) jointly in λ. Details of the skeleton expansion nature can be found in Ref. (3,5,6). The main point to note is that while g −1 , µ in Eq. (28) and Eq. (29) are expanded in powers of λ, the function g is kept unexpanded as a basis term (or "atom") of the skeleton expansion, temporarily ignoring its relationship as the inverse of g −1 . The equal time value of the variable γ in Eq. (11) is taken from the exact sum-rule for G in Eq. (9). The initial values at λ = 0 are g = g 0 and µ = 1, and we must remember to use the product form Eq. (25) to determine the local Greens function γ in Eq. (11). We should note that when the source is turned off V → 0 we recover space and time translation invariance so that Eq. (25) is simply G( k, iω j ) = g( k, iω j ). µ( k, iω j ), with the Matsubara frequency ω j = (2j + 1)πβ. At low T, the leading singularities of G are co-located with those of g provided the caparison function µ is sufficiently smooth-this situation is realized in all studies done so far.

C.
The shift identities and second chemical potential u 0 Before proceeding with the iterative scheme, it is important to discuss a simple but crucial symmetry of the t-J model-the shift invariance, first noted in Ref. (11). In an exact treatment shifting t ij → t ij + c t δ ij with c t arbitrary, is easily seen to be innocuous, it merely adds to Eq. (4) a term −c t σN σ whereby the center of gravity of the band is displaced.
(HereN σ is the number operator for electrons with spin σ.) However in situations such as the λ expansion, the Gutzwiller constraint is released at intermediate values, here it has the effect of adding terms derivable from a local (i.e. Hubbard type) interaction term. 39 To see this consider the fundamental commutator term [H tJ , X 0σ i i ] in Eq. (13), here under the shift t ij → t ij + c t δ ij , the third term gives rise to an extra term c t Xσ iσi i X 0σ i i . This term vanishes only in a Gutzwiller projected state, the equations of motion by themselves do not eliminate it. Its appearance is tantamount to adding a Hubbard like term ct An efficient method to do so is to explicitly introduce an extra Lagrange multiplier u 0 through a term λu 0 i N i↑ N i↓ to the Hamiltonian Eq. (4). This amounts to replacing in all terms other than in the bare propagator g 0 . The u 0 term makes no difference when λ is set at unity in the exact series, since double occupancy is excluded.
In practice, we set λ = 1 in equations that are truncated at various orders of λ, and the magnitude of u 0 is fixed through a second constraint. We thus have two variables to fix, namely u 0 and µ. We also have two constraints, the number sum-rules n σ = g σσ (iτ − , iτ ), and n σ = G σσ (iτ − , iτ )(Eq. (9)). In the absence of a magnetic field the number densities n σ reduce as n σ → n 2 , where n is the number of particles per site. After turning off the sources, in the momentum-frequency space we can further introduc- Here ε k and J k are the Fourier transforms of −t ij and J ij . In the right hand side of Eq. (31), the second and third terms arise respectively from the equal-time limit of λY 1 and λX.g.g −1 in Eq. (28) respectively. The two self energies Φ, Ψ are explicitly λ dependent, they vanish at infinite frequency for any λ. Thus we write The auxiliary Greens function satisfies a second sum-rule that is identical to Eq. (9), both may written in the Fourier domain: Eq. (25) can now be written explicitly in the non-Dysonian form proposed in Ref. (3) We observe that simple Fermi liquid type self energies Ψ and Φ can, in the combination above, lead to highly asymmetric (in frequency) Dyson self energies 3,5,6,11 . Finally we note that our calculations are performed in terms of spectral function obtainable from analytic continuation of the Matsubara frequencies into the upper complex half plane of frequencies: and similarly defined spectral functions for variables g, Φ, Ψ etc.

D. Summary of equations to second order in λ
In the following, we use the minimal second order equations 5-7 obtained by expanding Eq. (28) and Eq. (29) to second order in λ. The calculation is straightforward and a systematic notation is detailed in Ref. (6), which is followed here. We use the abbreviation 8 k ≡ ( k, iω n ), and also redefine Φ(k) = χ(k) + ε k Ψ(k), keeping in mind that one set of terms in Φ have an external common factor of ε k multiplied by all terms in Ψ. We next collect the answers below in terms of the two self energies χ, Ψ We now expand Ψ and χ from Eq. (28) and Eq. (29) in powers of λ. To the lowest two orders we find Ψ = λΨ [1] . The next step is to introduce u 0 explicitly: we write ε k → ε k = ε k − u 0 2 in every occurrence of ε k , except in the bare propagation term (the term with an underbrace) in Eq. (36).
Note that the shift with u 0 also applies to the term χ [0] , it now reads . We note the expressions for χ [1] , Ψ [1] from Ref. (6) Eq. (65-67): We now set λ = 1 and record the final equations: where we used a nearest neighbor J ij and cubic symmetry in the simplifications. We can verify that the above expressions obey the shift invariance: if we shift ε k → ε k + c 0 , the arbitrary constant c 0 can be absorbed by shifting µ → µ + c 0 and u 0 → u 0 + 2c 0 , and is thus immaterial. The band energy is given explicitly as , where t and t are the first and second neighbor hopping amplitudes.

E. High energy cutoff scheme
The self consistent solution of the second order equations of Eqs. (38,39,40,41) plus the number sum-rules, can be found numerically by discretizing the momentum and frequency variables on a suitable grid. This procedure can be carried out in a straightforward way for low T < ∼ t and high hole densities δ > ∼ 0.3 (low particle densities n < ∼ 0.7). At lower hole densities or at high temperature T t, the equations run into convergence problems. The origin of this problem is the formation of weak and featureless tails of the spectral functions extending to quite high energies. These tails are known to be artificial, since they do not occur in the exact numerical solutions where available. Thus the second order theory seems insufficient in the regime of low hole densities δ < ∼ 0.2, where much of the current interest lies. A technically rigorous resolution of the problem of weak tails seems possible. However it requires the non-trivial calculation of higher order terms in the λ expansion. Such higher order terms oscillate in sign and hence cancellations at high energies are expected.
In view of the substantial magnitude of the program of summing the λ series to high orders, it seems worthwhile to investigate simpler and physically motivated approximations for improving the lowest order scheme. It turns out that there are a few interesting alternatives in this direction. In Ref. (5) we showed one convenient way to handle the high energy tail problem practically, through the introduction of a high energy cutoff. The choice of an objective cutoff was rationalized by considering two physically different limits, that of high particle density n → 1 and the simpler high temperature limit, where related tails are found.
The cutoff is chosen using the analytically available high T limit results and then applied to all densities and T.
The cutoff scheme of Ref. (5) is not rigorous, but enables us to extract meaningful results for low energy excitations from the second order λ equations, out to fairly low hole densities where W T is a smooth window (even) function shown in Fig. (3) Ref. (5). It is centered on the bare band energy and has width 4D, where 2D is the bandwidth (∼ 8t in this case).
The constant N k is fixed by the normalization condition ρ g ( k, ω)dω = 1. In the present case of 2-d we can impose this cutoff window at each k individually, so that only k states very far from the chemical potential are affected.
The two chemical potentials µ and u 0 are determined through the number sum rules written in terms of the Fermi function f (ω) = (1 + e βω ) −1 and the spectral functions:  Following 2 we assume t ∼ 0.45 eV, giving a bandwidth ∼ 3.6 eV. §Single-particle spectrum: The quasiparticle energy E( k) and quasiparticle weight Z( k) are found from G as usual 5 . In Fig. (1) we display the hole density δ and t dependence of the low temperature Z(k F ), along the nodal (i.e. 11 ) direction. The typical magnitudes of Z are comparable or lower than those reported in d = ∞ 5 . A new and important feature is the strong sensitivity of Z(k F ) to the sign and magnitude of t /t. Both decreasing t (at fixed δ) and decreasing δ (at fixed t ) reduce Z. This feature is basic to understanding our main results. We next study the decay rate of the electrons found as the half-width at half-maximum of the spectral function ρ G ( k, ω) at fixed k. We display the T variation of Γ and − m Σ at the Fermi surface for three representative values of t /t in Fig. (2). Both variables display considerable variation with modest change of T .
The case of t > 0 shows a distinct quadratic T dependence, but for t ≤ 0 we note the strong reduction, or absence, of such a quadratic dependence. Below we note a closely parallel T and t dependence of the resistivity. In Fig. (3) we display the photoemission accessible peak heights of the spectral function {t ρ G ( k, ω)} max over the BZ at three representative values of t /t, at three temperatures T = 63, 210, 334K. The peaks track the non-interacting FS, changing from hole-like in Panels (a,b,c) to strongly electron-like in Panels (g,h,i). Several features are noteworthy. The peaks are higher in the nodal relative to the anti-nodal direction at low T. We observe the high sensitivity to warming, in going from T = 63K to T = 334K a small (∼ 0.7%) change in T relative to the bandwidth causes a five to fifteen-fold drop in the spectral peaks at the Fermi surface. This is correlated to the thermal variation of − m Σ at the same set of t , shown in the inset of Fig. (2), since the intensity at k F is essentially the inverse of this object.
Meanwhile the background spectral weight rises rapidly in all cases, to a roughly similar magnitude. The figure shows that at low T the curve with t > 0 has much higher peaks than t ≤ 0, giving the impression of weaker correlations. However the drop on warming is the largest in this case, which signifies another facet of strong correlations. §Resistivity: We now study the behavior of the resistivity from electron-electron scattering. We use the popular bubble approximation, factoring the current correlator as , where the bare current vertex is the velocity v α k = ∂ε k ∂kα . In tight binding theory v α k oscillates in sign, resulting in a vanishing average over the Brillouin zone. This oscillation is expected to reduce magnitude of the vertex corrections 16  Recalling that our bandwidth is ∼ 3.6 eV, we observe that the tiny 0.35% variation of temperature relative to the bandwidth, in warming from 63K to 210K drops the peak height by a factor between 5 and 10, followed by a more gentle fall to 334K.
This extraordinary thermal sensitivity is characteristic of our solution of the t-J model, it is also reflected in other variables discussed here such as the resistivity. each sheet represented by the 2-d t-J model, the DC resistivity ρ xx can be written in terms of the electron spectral function as follows. We define a dimensionless resistivityρ xx whose inverse is given byσ where A k ≡ 1 Ns k A( k), while the momentum resolved relaxation scale is: and f ≡ 1/(1 + exp βω) is the Fermi function. This object resembles the spectral peaks in In Ref. (36,37) the authors note that in a metallic system this parameter is expected to be greater than unity, and its least value is k F λ m = 1 for the case of unitary (impurity) scattering. Thus we expect that ρ xx ≤ ρ 0 , i.e.ρ xx ≤ 1 in a good metal. The Ioffe-Regel-Mott resistivity scale used in Ref. (5,7,9) provides a similar measure for quantifying the magnitudes of the resistivity found in strongly correlated metals. However we should keep in mind that both estimates suffer from ambiguities in defining a precise threshold value of resistivity, since factors of 2 (or of 2π) cannot be ruled out in Eq. (47). In all curves we see that the curvature changes from positive (for t ≥ 0) to negative (for t < 0) at say 150K. To understand the role of t /t we note that the resistivity in Eq. (45) depends on t /t through the velocity v x k , in addition to a dependence through the self energies Eqs. (39,38). To gauge their relative importance it is useful to examine a local approximation of Eq. (45) where the two functions are averaged separately over momentum: The velocity squared average is independent of the sign of t , therefore the local approximation, shown in the inset of Panel(b), probes only the dependence through Eqs. (39,38).
Comparing the inset and main figure in Panel(b), we see that at t = 0.4t both resistivity curves display a positive curvature. At t = −0.4t we see that ρ Local is essentially linear in T, while ρ xx shows a negative curvature. The behavior of −Σ in the inset of Fig. (2) qualitatively resembles the resistivity. The difference is actually related to the velocity factors, which are very different effect between t < 0 and t > 0. These cause the integrals to have very different thermal variation. §Hall response: Within the bubble scheme, we may also calculate the Hall conductivity 10,23-25 as σ xy = −2π 2 /ρ 0 × ( Φ Φ 0 ) ×σ xy , the dimensionless conductivity: and Here Φ = Ba 2 0 is the flux 17 and Φ 0 = hc/(2|e|) is the flux quantum. In terms of these we can compute the Hall constant R H and Hall angle with v 0 = (a 2 0 c 0 ) 17,26 . In Fig. (5) we display the computed Hall variables. In Panel (a) tan Θ H is shown for two values of t /t displaying hole-like and electron-like behavior. A decrease in hole density reduces the magnitude in either case. In Panel (b) we display the computed cot(Θ H ) versus T 2 with three values of t giving an electron-like FS. We note that cot Θ H is approximately linear with T 2 27-29 and is strongly affected by the magnitude of t . The two distinct cot(Θ H ) ∝ T 2 regimes seen in Fig. (5-b) are also seen in many experiments, the crossover to a strange metal corresponding to the bending temperature 29 . Our results for the hole-like case t < 0 at low T are less robust due to the small magnitude of σ xy and errors from the oscillating sign of η in Eq. (49). In Panel (c) we show the Hall constant R H at three densities for representative values of t /t. Its sign is electron-like for t > 0 and hole-like for t ≤ 0, tracking the change in topology of the Fermi surface in Fig. (2). The magnitude of R H is substantially affected by changing t . This is a strong correlation effect, and discourages envisaging any simple relationship between the Hall number and hole density. By varying t', the second neighbor hopping at a fixed t and J, we found in Fig. (1) a remarkable variation of the quasiparticle weight Z(k F ) that is characteristic of the 2-d square lattice, with no simple analog in d = ∞. We found t < 0 leads to a considerable reduction in its magnitude, while t > 0 leads to a larger value and thus a more robust Fermi liquid. A direct calculation of the single particle spectral width Γ = −ZΣ confirms this observation in Fig. (2), and when studied as a function of the temperature, shows a much larger magnitude, and hence broader spectral lines.
Our two striking results concern the spectral heights over the Brillouin zone, and the resistivity as a function of T at various densities and t'. The spectral height is the peak value of ρ G ( k, ω) scanned over ω, and equals the inverse of the least magnitude of mΣ( k, ω).
In Fig. (3) we present both the T evolution (going horizontally) and the t' evolution (going vertically) of this important object visible in ARPES. We note that t < 0 model with a very small Γ also displays a rapid loss of coherence on warming. The quasiparticle peaks drop rapidly, while the valleys, representing the background spectral weight in photoemission, catch up with the peaks in magnitude. A similar variation happens for t = 0 but the drop of the peak heights is more pronounced. The case of t > 0 has the largest drop of peak heights, while its effective Fermi temperature is the largest of the three cases. It follows that the electron doped case has a more robust Fermi liquid appearance for T lower than its Fermi scale. Our study provides absolute scale values for these observable heights, and it should be interesting to study these experimentally for comparison. Towards that objective we note that t > 0 maps to the electron doped High Tc superconductors, while t ≤ 0 maps to the hole doped cases, as we may also deduce from the shapes of the Fermi surfaces in the above figure.
The other striking results concerns the resistivity. We are able to calculate the longitudinal resistivity ρ xx on a doubly absolute scale, both the magnitude of ρ xx and that of T are given in physical units by using reasonable values for the basic parameters of the t-J model and the lattice constants Fig. (4). We find essentially the experimentally observed scales for both axes, and there is room for further adjustments of bare scales if needed.
The main finding is that as δ is varied towards half filling, the regime of linear resistivity increases in the hole-like cases (t ≤ 0) and the quadratic dependence regime shrinks to very low T scales, falling below the known superconducting transition temperatures. The other important finding is that the concavity (convexity) of resistance versus T, usually taken to denote a (non) Fermi liquid behavior, requires an enlarged viewpoint; we find that the sign of t' flips the two cases. As an example, the case t ≤ 0 has a pronounced convex regime at low T. This could be naively ascribed to a non-Fermi liquid behavior, but in reality is a crossover range to the strange metal regime.
We also present results in these results are encouraging, we believe that further work is needed to unravel what we appears to be a knee in the cot(Θ H ) versus T 2 curve, and also to better estimate the density dependence of the Hall constant.
In this work, we used a scheme from the ECFL theory where the second order λ expansion terms are supplemented with a high energy cutoff. This scheme has been benchmarked in d = ∞ against DMFT 5,7 for computing transport and other low energy excitations, giving good agreement with exact numerical results. As detailed in Ref. Hence we might expect to find further refinements of the absolute values of the quasiparticle weight and resistivity to emerge from further studies of the λ expansion. However it seems likely that the crucial variation of resistivity and Hall constant with the magniude and sign of t found here will persist in more exact future results. Hence it seems that we can draw some useful conclusions already regarding the difference between hole and electron doping.
We have shown a range of results for the 2-d t-J model, obtained by varying different parameters within our scheme. It is interesting that the magnitudes of various transport variables, presented here in physical units 30 , are roughly on the scale of reported measurements [18][19][20][21][22][31][32][33] . Although it is not our primary aim here to produce exact fits, we note that the agreement can be improved in many cases with suitable changes of the bare (band) parameters.
In the range of parameters considered here, a metallic state has been posited, and therefore the role of the exchange J is limited; we find very little variation of the transport quantities with a change in J. The transport parameter variation with density seems very similar to that found in d = ∞ in Ref. (5,7,9) where a large variety of Gutzwiller correlated metallic states were shown to arise 7 , with their origin in the U = ∞ or Gutzwiller correlation rather than with J. The added feature in d = 2 is the important role played by t', as stressed here. We expect magnetic, superconducting and possibly other broken symmetry states at the lowest T and δ to arise, largely due to the effect of J. Further work is necessary to find reliable calculational schemes for these broken symmetry states.
A few broad conclusions suggest themselves. The parameter t /t plays a key role in determining the low-energy scales. In Fig. (1) we see that the quasiparticle weight Z has a large variation with t . The origin of this sensitivity lies in the self energies in Eqs. (39,38), where combinations of the band energies ε k play the role of an effective interaction. Varying t /t therefore changes the self-energies strongly, in contrast to the usual weak change via the band parameters in Eq. (37).
Our main findings are as follows. (I) The spectral functions are highly sensitive to thermal variation; in Fig. (3) we observe a five to fifteen fold drop in intensity with a variation of k B T about 1/100 th the bandwidth ∼ 3.6 eV. This is in severe conflict with expectations from conventional theories of metals. (II) We note from Fig. (4) that with t ≤ 0, a Fermi liquid (FL) resistivity ρ ∝ T 2 is seen only at very low T. The very low T (FL) regime is followed by a "strange metal" regime, also at low T , where we find a ρ vs T curve, with zero or negative curvature. This regime parallels the Gutzwiller-correlated strange metal regime reported in