Universal linear-temperature resistivity: possible quantum diffusion transport in strongly correlated superconductors

The strongly correlated electron fluids in high temperature cuprate superconductors demonstrate an anomalous linear temperature (T) dependent resistivity behavior, which persists to a wide temperature range without exhibiting saturation. As cooling down, those electron fluids lose the resistivity and condense into the superfluid. However, the origin of the linear-T resistivity behavior and its relationship to the strongly correlated superconductivity remain a mystery. Here we report a universal relation dρ/dT=(μ0kB/ℏ)λL2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d{\boldsymbol{\rho }}/{\boldsymbol{dT}}=({{\boldsymbol{\mu }}}_{{\bf{0}}}{{\boldsymbol{k}}}_{{\boldsymbol{B}}}/{\boldsymbol{\hslash }})\,{{\boldsymbol{\lambda }}}_{{\boldsymbol{L}}}^{{\bf{2}}}$$\end{document}, which bridges the slope of the linear-T-dependent resistivity (dρ/dT) to the London penetration depth λL at zero temperature among cuprate superconductor Bi2Sr2CaCu2O8+δ and heavy fermion superconductors CeCoIn5, where μ0 is vacuum permeability, kB is the Boltzmann constant and ħ is the reduced Planck constant. We extend this scaling relation to different systems and found that it holds for other cuprate, pnictide and heavy fermion superconductors as well, regardless of the significant differences in the strength of electronic correlations, transport directions, and doping levels. Our analysis suggests that the scaling relation in strongly correlated superconductors could be described as a hydrodynamic diffusive transport, with the diffusion coefficient (D) approaching the quantum limit D ~ ħ/m*, where m* is the quasi-particle effective mass.

In quantum mechanics, the uncertainty principle gives rise to quantum fluctuations of the system that may impose some universal bound on its physical properties. Calculations based on the AdS/CFT (Anti de-Sitter/ Conformal Field Theory) have suggested a lower bound for the liquid viscosity, η/s ≥ ħ/4πk B 1 , where η is the shear viscosity and s the entropy. Recent experiments also revealed a quantum bound D s ≥ ħ/m for the spin diffusivity D s in a strongly interacting Fermi gas 2,3 . Here ħ is the reduced Planck constant and m is the mass of particles. It is therefore interesting to ask if such a lower bound may be realized in the electron transport of strongly correlated quantum critical systems.
One of the distinguished features of strongly correlated cuprate superconductors is the linear-temperature (T) dependent resistivity 4 , which could extend to very high temperature 4 and violate the Mott-Ioffe-Regel (MIR) limit 5 . The linear relationship has also been observed in some heavy fermion superconductors 6,7 , starting at around the superconducting transition temperature T c and extending to high temperatures of about 10-20 times of T c . Many different mechanisms have been proposed to explain the microscopic origin of the linear-T resistivity behavior including quantum critical theories and the more exotic AdS/CFT calculations. On the other hand, recent experiment suggested that the linear-T resistivity in different materials may share a similar scattering rate 8 .
In the present work, we investigated the linear-T resistivity in a number of strongly correlated superconductors and demonstrate a connection between its coefficient and the superfluid density responsible for the charge carrying in the superconducting state. We show that this can be understood by a diffusion transport of heavy

Results
We start with the heavy fermion superconductor CeCoIn 5 . Among all strongly correlated superconductors, CeCoIn 5 is remarkably similar to the high T c cuprate superconductors in several aspects 9 . For example, it has also a two-dimensional Fermi surface 10,11 , its superconducting phase is near to an antiferromagnetic phase [12][13][14] , and its superconducting gap has d-wave symmetry [15][16][17] . Besides, CeCoIn 5 is one of the purest strongly correlated superconductors 6,7 , with a tunable linear-T resistivity under modest applied pressure [18][19][20][21] . To examine its transport properties, we have therefore grown high quality CeCoIn 5 single crystal samples by an indium self-flux method 6 and performed detailed transport measurements under pressure to avoid disorder related effects. Figure 1(a) demonstrates the T-dependent resistivity curve of CeCoIn 5 under pressure from 0 GPa to 1.0 GPa. All of them exhibit a perfect linear-in-T resistivity from around T c to about 20 K as indicated by the dashed lines. The inset of Fig. 1(a) shows the T-dependent resistivity of CeCoIn 5 up to 300 K. For comparison, we also plot in Fig. 1(b) the resistivity of Bi 2 Sr 2 CaCu 2 O 8+x from underdoped to overdoped regime with the oxygen contents from x = 0.2135 to 0.27 22 . Figure 1(c) demonstrates the dρ/dT versus λ L 2 for both compounds, using the experimental results for the penetration depth measured previously by muon spin spectroscopy 23 and ac susceptibility 24 . We see remarkably that all the investigated samples fall on the same straight line described by , with a coefficient that is determined entirely by the fundamental constants (μ 0 : the vacuum permeability; k B : the Boltzmann constant; ħ: the reduced Planck constant). This indicates a universal origin for the charge transport in both compounds.
The above relation between dρ/dT and λ L 2 can be extended to various other strongly correlated superconductors with linear-T resistivity. The data are summarized in Fig. 2 on a log-log scale. Most resistivity data were taken from experimental results on high-quality single crystal samples in order to obtain the intrinsic linear-in-T coefficient. The values of the penetration depth were obtained by muon spin spectroscopy 25 , optical conductivity measurement 26,27 and some other techniques. Note that for superconducting thin films, the experimental magnetic penetration depth generally deviates from the London penetration depth λ L due to structural disorders in the films 28,29 . Even in high quality ultrathin films, there is a large difference in superfluid density between the film and the bulk materials with same T c 30,31 . Consequently all the data of the London penetration depth shown in Fig. 2 were taken only from bulk materials. It is worth noting that Fig. 2 also includes the transport data for cuprate superconductors along different transport directions, e.g., YBa 2 Cu 3 O 6.93 along the a, b and c-axis. Cuprates generally exhibit a metallic in-plane resistivity but an insulating-like resistivity along the c-axis below certain temperature, which reflects the two-dimensional nature of the system. Correspondingly, the penetration depth along the c-axis is determined by a Josephson-coupling between superconducting layers [32][33][34] , which is different from the in-plane one 26 . Thus it is amazing to observe that the same scaling relation holds true for both directions. Combining the data for all the strongly correlated superconductors summarized here, we see that the scaling, ρ µ , spans over several orders of magnitude. Note that the in-plane LSCO data in the extremely underdoped regime 0.07 ≤ p ≤ 0.12 demonstrates a systematic deviation from the scaling relationship as shown in Fig. 2. The deviation could be understood in terms of the complex competing phase, like charge density wave and pseudogap, which become significant in the underdoped regime.
, where ρ is in the unit of μΩcm, λ L is in μm, and T is in K. See Table 1  The above scaling relation is consistent with several well-known experimental facts. First, considering the special case at T = T c and neglecting the residual resistivity, the scaling relation where σ c is the dc conductivity at T c 26 . Second, the Drude formula 35 is often used to describe the resistivity of conventional metals, ρ = m*/n n e 2 τ, where m* is the effective mass of the quasi-particles, n n is the carrier density of quasi-particles, e is the charge of electrons, and τ is the relaxation time. If we naively match the Drude formula with the above scaling relation for a non-quasiparticle system and assume that the normal fluid and the superfluid are composed of the same charge carriers, λ L = (m*/μ 0 n n e 2 ) 1/2 , we obtain immediately a material-independent scattering rate τ −1 = k B T/ħ for all these strongly correlated superconductors. This is consistent with the universal scattering rate recently observed in the linear-in-T resistivity region among good and "bad" metals 8 . However, one can not take it for granted that the normal fluid in Drude model and superfluid in London equation are always the same. Actually, experiments showed that only part of normal carriers condensate into superfluid 36 . In addition, the measurements of the London moment already revealed the mass of Cooper pairs are undressed and have twice of the electron's bare mass, regardless the conventional metal superconductors 37 , heavy fermion superconductors 38 or cuprates 39 , which is different from the effective mass in the Drude formula. These results suggest that the mass and carrier density of the superfluid (n s ) and the normal fluid (n n ) are different in strongly correlated superconductors. So one can not directly obtain the universal scaling relation simply by replacing ρ with Drude model and λ L with London equation. The universal scaling relation 2 has much deeper physics, which directly links the superfluid at zero temperature to the normal fluid responsible for the linear-in-T resistivity in strongly correlated superconductors. It reveals an underlying relation between the superfluid and normal carriers: n s /m e = n n /m*. And indeed experimental evidence shows that about one of fourth normal carriers 36 condensates into superfluid in optimal doped cuprates while the effective mass of optimum cuprates is about 3-4 times of the electron free mass 36,40 , which validate n s /m e = n n /m*.
The above result provides important information on the nature of the electron transport in the quantum critical regime. Recently, several experiments have shown that electrons in solid can exhibit hydrodynamic flows similar to a classical viscous liquid, if the electron fluid equilibrates by the electron-electron collisions [41][42][43] . Thus the electron transport in strongly correlated superconductors, where electron-electron interactions play a major role in the scattering processes, might in principle have a hydrodynamic description. Consequently, its linear-in-T resistivity could be described by the well-known Einstein's relation 44 , an important law for the hydrodynamic transport, which states that the mobility (μ) of a particle in a fluid is related to its diffusion coefficient (D), namely, D = μk B T. Hence we have ρ = k B T/n n e 2 D and in the linear-in-T regime, the diffusion coefficient D must be a temperature-independent constant. Combining this and the scaling relation immediately yields D = ħ/m*, which is the quantum limit of the charge diffusion coefficient for the quasi-particles with an effective mass, m*. This is one of the most important consequence of our observations. Actually, the quantum limit of the diffusion coefficient was recently observed in cold fermionic atomic gases in the unitary limit of scattering 2, 3 . It implies that quantum diffusion transport might be a universal property of strongly correlated fermionic systems where the electron scatterings are so strong that the transport becomes highly incoherent. In fact, it was proposed recently that the transport in an incoherent metal is controlled by the collective diffusion of energy and charge 45 , supporting the proposed scenario of quantum diffusion transport in the present work. Thus, the obtained scaling relation suggests the superfluid could also be governed by the quantum diffusion, since it connects the ground state with the normal state in the strongly correlated superconductors.
Our results also provide some insights on the nature of strongly correlated superconductivity, which is often born out of strongly correlated normal fluid in the quantum critical regime. Since the latter already approaches  Table 1 for details, including errors.
www.nature.com/scientificreports/ 5 SCIentIFIC RepORTS | 7: 9469 | DOI:10.1038/s41598-017-09792-z the quantum diffusion limit before it transits into the superfluid state, it implies a zero-point motion of the superfluid. Some people considered the quantum diffusion as a necessary condition for the presence of superfluid 46,47 . In fact, the quantum diffusion might explain the Uemura results for superconducting transition temperatures. Y. Uemura et al. observed that the underdoped cuprate superconductors exhibit a Bose-Einstein-condensation (BEC)-like superconducting transition but with a reduced transition temperature 25,48 . Actually, the BEC generally occurs when the thermal de Broglie wavelength λ dB is comparable to the distance between bosons, where λ dB characterizes a length scale within which the bosons can be regarded as quantum mechanical wave-packets. However, the quantum diffusion gives a new length scale ξ τ = D Th with τ = ħ/k B T, which characterizes the length scale that carriers can travel before losing their quantum coherence. Since the diffusion length of electron pairs under certain temperature, it makes ξ Th a new dephasing length to determine the BEC temperature. Thus the BEC temperature (T B ) is reduced to T c /T B = (ξ Th /λ dB ) 2 = 1/π as observed in the Uemura plot 25,48 .

Conclusion
In summary, we observed a universal scaling relation 2 , which connects linear-T-dependent resistivity to superconducting superfluid density at zero temperature in strongly correlated superconductors. Our analysis suggests that the quantum diffusion might be the origin of this scaling relation. In this case, the charge transport is viewed as a diffusion process of quasi-particles with a diffusion coefficient that approaches the quantum limit, D ~ ħ/m*.

Method
The high quality CeCoIn 5 single crystal samples are grown by an indium self-flux method 6 . High quality crystals were chosen to perform the transport measurements. Four leads were attached to the single crystal, with the current applied parallel to the crystallographic a axis. The resistivity was measured both in ambient pressure as well as under hydrostatic pressure P.