1. Introduction
In the strong underdoping low-temperature region of the phase diagram of the cuprates, ARPES experiments prove the existence of gapless Fermi arcs (see e.g., [
1] and references therein). The Fermi liquid paradigm then would suggest a metallic resistivity with a temperature dependence
, with
at least greater than 1. In the same region of the phase diagram, however, the experiments show that the in-plane resistivity has a minimum at finite temperature, i.e., a metal-insulator crossover (see e.g., [
2]). Furthermore if the resistivity is suitably normalized it exhibits a universal behaviour [
3] which appears hardly compatible with an explanation of that crossover in terms of disorder localization.
We propose that the coexistence of the two phenomena can be explained in term of the FL* fractionalized Fermi liquid concept, introduced in [
4,
5], within an approach to the low-energy physics of the cuprates in terms of the
t-
-
J model.
There have been several proposals of a FL* nature of the above quoted region of the phase diagram of the cuprates, named “pseudogap phase” [
6,
7] (and references therein), but to our knowledge not a careful discussion of the consistency between Fermi arcs and metal-insulator crossover (MIC).
An FL* is an exotic fractionalized Fermi liquid which has hole-like (in the case relevant for the cuprates) excitations near a Fermi surface satisfying a generalization of Luttinger theorem, together with fractionalized excitations emerging from a topological order, in particular holons charged spinless and spinons neutral of spin , interacting via a so-called slave-particle gauge field . The hole is a bound state or a resonance arising from holon-spinon binding generated by the gauge attraction.
The present paper is part of a project aimed at analyzing the consequences of the proposed FL* nature of the cuprates. The general idea underlying this project is that one can explain the Fermi liquid (FL) versus non-Fermi liquid (NFL) dichotomies encountered in the cuprates, in the present case the FL nature of Fermi arcs and the NFL nature of the MIC, in the following way: If some physical response is dominated by the hole excitations (spinon-holon interaction), then its behaviour is close to that of a Fermi liquid, if instead it is dominated by spinon-spinon (and possibly also holon-holon) interaction then it has a clear Non-Fermi liquid character.
Of course it is known that the
t-
-
J model is insufficient to explain some phenomena appearing in the cuprates such as charge density waves (see e.g., [
8]), but apparently at least these phenomena are qualitatively irrelevant for the issue of the present paper at strong underdoping [
9]. More generally, although clearly such a simple model cannot reproduce real experimental data, we show that even with an approximate treatment many of their doping and temperature dependencies can be indeed qualitatively understood. Furthermore with the introduction of experimentally derived but doping and temperature independent scales also a semi-quantitative agreement is obtained in many cases.
2. General Outline of the Approach
In the description of the low-energy physics of the cuprates in terms of a two-dimensional (2D)
t-
-
J model the sites correspond to the Cu sites of a CuO plane in the cuprates and the Zhang-Rice singlets to the empty sites. The Hamiltonian is given by:
where
denotes a link between lattice sites
i and
j,
denotes (next-)nearest neighbor sites and
the spin index which is assumed summed up also in the following if repeated.
denotes the Gutzwiller projection implementing the no-double occupation and describing the Mott physics of the system. It can be tackled with a spin-charge decomposition [
10,
11] of the hole field:
where the holon
h is a spinless fermion, so that by Pauli principle the no double occupation is implemented, and the spinon
s is a boson, satisfying in each site the constraint
[
12]. This decomposition involves an unphysical degree of freedom, because one can multiply at each site
s and
h by arbitrary opposite phase factors leaving
c unchanged. The corresponding local
gauge invariance is implemented by introducing a slave-particle gauge field
which in turn produces an attraction between holon and spinon. To show this effect of gauge fluctuations beyond perturbation theory, as a rough approximation, we apply in [
13] a kind of eikonal resummation.
However, a comparison with the corresponding one-dimensional (1D) model [
14,
15] suggests to modify the statistics of spinons and holons at lattice level, turning them into semions, i.e., their fields under equal time oriented exchange acquire a phase
(see e.g., [
16]). These phase factors are opposite for spinons and holons, so that
c retains its fermionic nature. This modified statistics is crucial for our discussion (as it was in 1 dimension), so let’s give some details.
The change of statistics is implemented in the Lagrangian formalism by coupling the holon to a U(1) gauge field
, related to the charge, with a Chern-Simons action with coefficient -1/2 and coupling the spinon to an SU(2) gauge field
, related to the spin, with Chern-Simons action with coefficient 1/2 [
12,
17]. This procedure corresponds to dressing the holon with a -1/2 charge flux and the spinon with a 1/2 spin flux. Due to the topological properties of Chern-Simons theory, if no approximations are made the result is an exact rewriting of the original model. However, if suitable approximations are introduced this rewriting provides a treatment of the model that with dimensional reduction was successful in deriving the correct large-scale behaviour of the model in 1D [
14]. Therefore we suggested that this provides a good approach also for the 2D model, relevant for the low-energy physics of the cuprates. The parameter region corresponding to what we call in this paper the pseudogap (PG) “phase” of the cuprates has an upper boundary in the (in-plane) doping (
)-temperature (
T) phase diagram identified with the inflection point of in-plane resistivity [
2]. Therefore the PG “phase” appears at low
and
T.
In PG the holon dressing by the charge gauge field
induces a uniform charge-flux
per plaquette [
12] and via Hofstadter mechanism [
18] it converts the spinless femionic holons
h into lattice Dirac fermions with dispersion defined in the Magnetic Brillouin Zone (MBZ). More precisely, the results of [
19] prove that in the
t-
model for elementary circuits, triangles and squares, the optimal flux at half-filling is
for triangles and
for squares. Therefore we assume that the “uniform” component of
B, denoted as in [
20] by
at sufficiently small doping and temperature is given as follows: let
i be a site in the even Néel sublattice and
be its nearest-neighbour sites in the 1 and 2 directions respectively and
its next-nearest-neighbour sites. Then
The fluctuations of
around its “uniform” component turn the holons into semions. Since semionic spinless excitations have the same Fermi surface (FS) of spin 1/2 fermions [
21,
22] the above
flux converts the fermionic holon with tight-binding dispersion
into a pair of lattice “Dirac fields”,
, with pseudospin index
a related to the two Néel sublattices and with dispersion:
restricted to the magnetic Brillouin zone. This produces small holon Fermi surfaces near
with Fermi momenta
in the MBZ, if identified with the union of two square regions, denoted as R(ight) and L(eft) centered at
and
, respectively. In the hole spinon-holon resonance the Dirac structure of such holons induces [
13] a momentum-dependent factor in the wave-function renormalization constant,
which reduces the spectral weight in the outer boundary of the MBZ, one key ingredient for a phenomenology of Fermi arcs at strong underdoping.
Let us turn to the gauge field associated to spin. The leading
component of the
field is along the magnetization axis of the undoped model and it introduces an interaction with the spinons, whose dynamics at large scales without
is described by a CP1 model, leading to an additional term
where, as result of a mean-field approximation,
. Neglecting a staggered term with mean zero, the field
describes quantum vortices centered on the positions of lattice holons, with opposite chirality for the two Néel sublattices, so that we call them “antiferromagnetic (AF) spin” vortices. In a quenched treatment of such vortices we derive the mean-field expectation value
, with
, which opens a mass gap for the spinon, consistent with th AF correlation length at small
extracted from the neutron experiments [
23]. Thus, propagating in the gas of slowly moving AF spin vortices, the AF spinons, originally gapless in the undoped Heisenberg model, acquire a finite gap, leading to a short range AF order.
This spinon gap is crucial because the spinon contribution to the slave-particle gauge action is then of Maxwell type, hence subleading with respect to the contribution of gapless holons. The holon contribution is characterized by the Reizer singularity in the dominant transverse component, killing the logarithmic confinement induced by the spinon contribution. Without the spinon gap the CP1 model would be in the broken symmetry “Higgs” phase and the spinon contribution to the gauge action would be gapped by the Anderson-Higgs mechanism, thus becoming the dominant one. As shown below the Reizer singularity of the slave-particle gauge field is a key factor for the discussion of resistivity and Fermi arcs. In particular the gauge interaction turns the spinon mass into a “complex temperature and frequency dependent mass” where c is a constant of and the holon diamagnetic susceptibility proportional to . Furthermore the spinon gap modifies the chemical potential of the holes with respect to that of the holons so that the Fermi momentum is of and the binding with the spinons produces a hole FS given by 4 small components centered at . The -flux introduces a topological order and the relevant version of the topological Luttinger theorem with the relation “fraction of the area of the BZ enclosed by the hole FS” = is satisfied. This is possible due to the modification of the holon chemical potential and actually, somewhat surprisingly, even a simple doping independent combination of the holon Fermi momenta and the spinon momenta corresponding to the spinon gap produces an area of the hole FS quite well compatible with that of the Luttinger theorem, for all the dopings in the physical range.
The above results show that the system behaves like a FL* liquid. Furthermore with the choice of
made above the FS of the hole near the diagonals of the BZ, where most of the spectral weight is concentrated, almost overlaps with the large tight-binding FS, see
Figure 1, as necessary for a consistency with ARPES data, although, admittedly, the overlap is not enough.
The last ingredient crucial in our discussion is a long-range attraction between the AF spin vortices centered on holons in two different Néel sublattices. If we treat in mean-field the spinons (
in Equation (
4) we see that the interaction mediated by AF spin vortices on holons is of 2D Coulomb type. From the study of planar Coulomb systems we know that below a crossover temperature
, a finite density of incoherent holon pairs appears. We already remarked that the low-energy holon modes in PG have two small FS in the MBZ; the 2D Coulomb interaction acts on these modes coupling opposites sites of each FS with an attractive potential with a range given by the screening length of the Coulomb gas. A BCS treatment provides p-wave pairings on the two FS, that glued in the entire BZ for the hole resonance produce a d-wave order parameter, as first suggested in a different setting in [
24]. However, we do not have condensation of holon pairs because the fluctuations of the phase of the pairing field are too strong [
25], restoring the slave-particle gauge symmetry broken by the BCS approximation.
3. FL* Approach to Fermi Arcs
Let us summarize in precise terms the output of the above considerations on the self-energy of the hole resonance. The binding induced by the gauge fluctuations dominated by the Reizer singularity computed with the above cited eikonal approximation produces an imaginary part of the self-energy of the form:
Notice that
is linear in
for sufficiently small frequencies, therefore the hole-resonances in PG behave like a Marginal Fermi liquid [
26]. Indeed, such behaviour has been found first in optimally [
27] and then also in underdoped [
28] cuprates.
The contribution to the real part of the self-energy,
corresponding to Equation (
5) is estimated by a Kramers-Kroenig procedure, with an UV cutoff at the charge-transfer gap
, as suggested by experimental data.
An additional term can arise due to standard Landau hole-hole interaction.
Let us turn to the effect of the attraction induced by the AF spin vortices.
In our approximation we keep constant (up to its
d-structure) the modulus of the holon-pairing parameter
, obtained as solution of a BCS-like equation [
25], near the Fermi surface. However, since the charge pairs are not condensed its phase is strongly fluctuating. Therefore, the field describing these phase fluctuations has a doping and temperature dependent gap, denoted by
, decreasing with
T, that modifies the standard BCS form of self-energy near the hole Fermi surface. For
we adopt a form suggested in [
29]:
where
is (an estimate of) the temperature of condensation of holon pairs, which however never occur in the parameter range considered (in [
30] it was set to 0). A reasonable approximation for the contribution to the self-energy turns out to be [
31]:
where
is the hole dispersion modified by
(sub-leading modification not taken into account in [
31]) and
the scattering rate of the hole without charge pairing. To understand its implications, let us consider its behaviour in two regions of energy: For
expanding the self-energy up to the second order in terms of
we get
Thus for low frequencies the effect of holon pairing appears through a contribution to the wavefunction renormalization constant not modifying the Marginal FL nature, but with a strong dependence on the direction of
that, together with the Dirac structure, heavily suppresses the spectral weight along the hole FS as we move towards the MBZ boundary from point in the diagonal inside the MBZ, see
Figure 1. The “shadow band” outside of the MBZ perhaps can then be invisible because it disappears in the background incoherent continuum.
For
we expand the self-energy in powers of
and we get
which is the self-energy of a d-wave superconductor. The spectral weight is symmetric in frequencies and non-vanishing on the FS, but it exhibits two symmetric SC-like peaks roughly at
. It turns out that approximately for
such SC-like peaks cover the FL one in the spectral weight. Due to the
k-dependence of
near the nodal region the behaviour is always FL-like, with an effective FS, while as we move towards the antinodal directions most of the spectral weight is concentrated on the SC-like modes, and even if there is a FS, it has negligible effects. We identify the Fermi arcs as the region of the FS where the FL peak is the dominant one, see
Figure 2. As a consequence of the temperature dependence of
the Fermi arcs shrink as we lower
T, see
Figure 3, in agreement with experiments [
32,
33].
The above computations clearly show that in PG on the hole resonance the effect of the gauge fluctuations is to produce an isotropic Marginal FL behaviour at small energy. Instead, the effect of the vortex attraction is an anisotropic d-wave like contribution on the Fermi arcs and an effective gap outside of them. This suppression of the spectral weight away from the diagonals of the BZ induced by AF vortices, combined with the suppression of the wave-function renormalization constant due to the Dirac structure (see Equation (
3)) produces a strong suppression of the spectral weight outside of the MBZ. Perhaps in this region the corresponding ARPES signal disappears in the thermal continuum background so that the Fermi arc might appear as consistent with the tight-binding Fermi surface. Small deviations near the boundary of the MBZ are however expected, as in [
34].
Actually [
31], above the continuum background, besides the coherent hole peak there is an additional contribution to the hole spectral weight of the incoherent part describing the high energy—momentum contribution of the unbinded spinon—holon system, see
Figure 2. Its leading contribution on the hole FS at energy
can be approximately estimated as
if
, where
is the density of states of the holon. Due to the linear dispersion of the holon in the PG,
is approximately linear in its argument and the incoherent component of the hole gives a contribution to the spectral weight at positive
growing approximately linearly in
starting from
. This additional linear contribution agrees with ARPES data in strongly underdoped cuprates (see e.g., [
35]).
4. FL* Approach to Metal-Insulator Crossover
In strongly underdoped cuprates, or more precisely in PG, the in-plane resistivity exhibits a quite peculiar behaviour: it increases more than linearly at intermediate temperatures, with a behaviour in the same range, but at lower temperatures, unless it drops going to zero at the critical temperature, it acquires a insulating behaviour. For sufficiently underdoped samples, therefore, one can identify a minimum in resistivity at a temperature, denoted by marking a metal—insulator crossover. As previously mentioned, at higher temperatures, above the behaviour, an inflection point in the resistivity as a function of temperature is frequently used in defining the lower pseudogap temperature .
Often the MIC is attributed to disorder-induced localization. However, this interpretation is at odds at least with the following facts: (i) The role of true disorder has been analyzed carefully in [
36] and although it can induce a MIC, the slope of the metallic behaviour in this case is independent of the degree of disorder, whereas such slope in the standard situation depends on the doping concentration. Perhaps instead the disorder is responsible for the logarithmic behaviour found below the MIC. (ii) If the resistivity is suitably normalized it exhibits a universal behaviour [
3], being both doping and even material essentially independent, in terms of
. Due to this universality one can see a relation between
and
; the last one has been shown to be independent of disorder [
36], reinforcing the idea that the MIC is not caused by disorder.
If the MIC is intrinsic it would be incompatible with the FL nature of Fermi arcs. In fact, following the analysis of the previous section, on the basis of the existence of gapless Fermi arcs with marginal FL behaviour, using Drude formula one would expect the in-plane resistivity to exhibit a metallic, linear in
T behaviour in PG. However, in a FL* theory since spinons and holons interact via the slave-particle gauge field, one should see if it is predominant the contribution of the hole-hole interaction, due to spinon-holon coupling, or the contribution due to spinon-spinon and holon-holon interactions. It turns out that the second alternative gives the leading effect for in-plane resistivity. To estimate the spinon-spinon interaction we employ the Kubo formula for the spinon current and adopt the eikonal approximation discussed above. For the holons, since the minimal coupling to Chern-Simons gauge fields doesn’t change local gauge-invariant quantities, we can use the Drude formula. Then the spinon and holon contributions are summed up following the Ioffe-Larkin rule [
37] i.e., the resistivity of the two subsystems are added:
. Due to the slave-particle gauge interaction the holon contribution is given by the standard expression [
38]
, where
is the holon Fermi energy. As discussed in detail in [
13] the Reizer singularity of the slave-particle gauge field modify the spinon current correlator giving rise to a contribution to the in-plane resistivity given by
with
is a kind of “thermal de Broglie wavelength”. Notice that the constant
c here is the same used in the discussion of Fermi arc and
. For T
100K up to almost
one finds a metallic behaviour,
, but in the limit
the
exhibits an insulating behaviour,
for low
T. Such spinon contribution is dominating w.r.t. the metallic holon contribution, so that the resistivity emerging from the Ioffe-Larkin rule at low
T is still insulating and the overall behavior is semi quantitatively consistent with the experimental data (except at lower
T where, as remarked above, perhaps standard localization occurs). This metal-insulator crossover is determined by the interplay between the AF correlation length
and the thermal de Broglie wave length
. Hence in the limit
we find a “peculiar” localization due to short-range AF order, producing an insulating behaviour due to the gauge interaction in the same range of parameters where Fermi arcs exist. Since the current considered in this computation is the same considered in the computation made with the hole resonances, the spinon behaviour of the resistivity is the dominating one thus providing a solution to the puzzle of the FL arcs vs NFL resistivity dichotomy. Furthermore the spinon dominance explain the universality noticed in [
3]. In fact [
39], let’s normalize the resistivity as follows:
Then, as follows from Equation (
10), neglecting the holon contribution
would be a universal function of the normalized temperature
and the characteristic temperature
itself is essentially determined by spinons. In the Ioffe-Larkin rule used a relative coefficient between the holon and the spinon contribution appears as a free parameter, independent of both
T and
, due to the use of a scale renormalization in the continuum limit in the computation of the spinon contribution. An optimized coefficient, while adding a slight degree of non-universality, allows for a reasonably correct fit of experimental data, over a wide range of temperatures and dopings from
to
. Let us notice as a marginal note [
40] that the peculiar dependence on
T of the “complex mass”
produces for many decades of K the linear in
behavior found in [
41,
42], without attributing it to the quasi-particle structure there advocated. A final remark: the inflection point point in temperature of the in-plane resistivity identifying
in our formalism depends on a delicate cancellation between the doping dependence of the spinon gap
and the Fermi momenta of the holons,
and the result well compares with the experimental data. For this cancellation to occur it is crucial that
, it would not occur for the hole Fermi momenta
.